A model of propagating calcium-induced calcium release mediated by calcium diffusion

The effect of sudden local fluctuations of the free sarcoplasmic [Ca++]i in cardiac cells on calcium release and calcium uptake by the sarcoplasmic reticulum (SR) was calculated with the aid of a simplified model of SR calcium handling. The model was used to evaluate whether propagation of calcium transients and the range of propagation velocities observed experimentally (0.05-15 mm s(-1)) could be predicted. Calcium fluctuations propagate by virtue of focal calcium release from the SR, diffusion through the cytosol (which is modulated by binding to troponin and calmodulin and sequestration by the SR), and subsequently induce calcium release from adjacent release sites of the SR. The minimal and maximal velocities derived from the simulation were 0.09 and 15 mm s(-1) respectively. The method of solution involved writing the diffusion equation as a difference equation in the spatial coordinates. Thus, coupled ordinary differential equations in time with banded coefficients were generated. The coupled equations were solved using Gear's sixth order predictor-corrector algorithm for stiff equations with reflective boundaries. The most important determinants of the velocity of propagation of the calcium waves were the diastolic [Ca++]i, the rate of rise of the release, and the amount of calcium released from the SR. The results are consistent with the assumptions that calcium loading causes an increase in intracellular calcium and calcium in the SR, and an increase in the amount and rate of calcium released. These two effects combine to increase the propagation velocity at higher levels of calcium loading.