Ryanodine receptor gating controls generation of diastolic calcium waves in cardiac myocytes.

The role of cardiac ryanodine receptor (RyR) gating in the initiation and propagation of calcium waves was investigated using a mathematical model comprising a stochastic description of RyR gating and a deterministic description of calcium diffusion and sequestration. We used a one-dimensional array of equidistantly spaced RyR clusters, representing the confocal scanning line, to simulate the formation of calcium sparks. Our model provided an excellent description of the calcium dependence of the frequency of diastolic calcium sparks and of the increased tendency for the production of calcium waves after a decrease in cytosolic calcium buffering. We developed a hypothesis relating changes in the propensity to form calcium waves to changes of RyR gating and tested it by simulation. With a realistic RyR gating model, increased ability of RyR to be activated by Ca2+ strongly increased the propensity for generation of calcium waves at low (0.05-0.1-µM) calcium concentrations but only slightly at high (0.2-0.4-µM) calcium concentrations. Changes in RyR gating altered calcium wave formation by changing the calcium sensitivity of spontaneous calcium spark activation and/or the average number of open RyRs in spontaneous calcium sparks. Gating changes that did not affect RyR activation by Ca2+ had only a weak effect on the propensity to form calcium waves, even if they strongly increased calcium spark frequency. Calcium waves induced by modulating the properties of the RyR activation site could be suppressed by inhibiting the spontaneous opening of the RyR. These data can explain the increased tendency for production of calcium waves under conditions when RyR gating is altered in cardiac diseases.

Electrodiffusion kinetics of ionic transport in a simple membrane channel.

We employ numerical techniques for solving time-dependent full Poisson-Nernst-Planck (PNP) equations in 2D to analyze transient behavior of a simple ion channel subject to a sudden electric potential jump across the membrane (voltage clamp). Calculated spatiotemporal profiles of the ionic concentrations and electric potential show that two principal exponential processes can be distinguished in the electrodiffusion kinetics, in agreement with original Planck's predictions. The initial fast process corresponds to the dielectric relaxation, while the steady state is approached in a second slower exponential process attributed to the nonlinear ionic redistribution. Effects of the model parameters such as the channel length, height of the potential step, boundary concentrations, permittivity of the channel interior, and ionic mobilities on electrodiffusion kinetics are studied. Numerical solutions are used to determine spatiotemporal profiles of the electric field, ionic fluxes, and both the conductive and displacement currents. We demonstrate that the displacement current is a significant transient component of the total electric current through the channel. The presented results provide additional information about the classical voltage-clamp problem and offer further physical insights into the mechanism of electrodiffusion. The used numerical approach can be readily extended to multi-ionic models with a more structured domain geometry in 2D or 3D, and it is directly applicable to other systems, such as synthetic nanopores, nanofluidic channels, and nanopipettes.