Dynamic simulation of the effect of calcium-release activated calcium channel on cytoplasmic Ca2+ oscillation

A mathematical model is proposed to illustrate the activation of STIM1 (stromal interaction molecule 1) protein, the assembly and activation of calcium-release activated calcium (CRAC) channels in T cells. In combination with De Young-Keizer-Li-Rinzel model, we successfully reproduce a sustained Ca(2+) oscillation in cytoplasm. Our results reveal that Ca(2+) oscillation dynamics in cytoplasm can be significantly affected by the way how the Orai1 CRAC channel are assembled and activated. A low sustained Ca(2+) influx is observed through the CRAC channels across the plasma membrane. In particular, our model shows that a tetrameric channel complex can effectively regulate the total quantity of the channels and the ratio of the active channels to the total channels, and a period of Ca(2+) oscillation about 29 s is in agreement with published experimental data. The bifurcation analyses illustrate the different dynamic properties between our mixed Ca(2+) feedback model and the single positive or negative feedback models.