A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation

Figures

• 

  Figure 1. A schematic diagram of ventricular myocytes modeled as a rectangular solid.

• 

  Figure 2. The placement of calcium release sites (A) at one Z-depth, and (B) at one Z-disc. The inset in (A) shows the CRUs on two T-tubules at two adjacent Z-discs. The distribution of inter-CRU distance is derived based on the experimental data.

• 

  Figure 3. The schematic diagram shows a grid point as a cube of size 200 nm on each dimension in the 3-dimensional space. In this cube, we put a single dyad. Assuming that the width of the dyad is 300 nm, we model the efflux of calcium flowing into 4 adjacent grid points, and it is equally split into 4 parts. The grid location is given by the coordinate (i,j,k).

• 

  Figure 4. Using the simulation data, a multitude of information can be extracted. Apart from a snapshot of a pseudo-line scan calcium transient during Ca2+ transient along the longitudinal direction, as shown in (A), we are able to show (B) the dynamics of free calcium content, (C) the dynamics of calcium-bound fluorescent, (D) the dynamics of total calcium content.

• 

  Figure 5. The probability of one CRU triggering the neighboring one at different distance and the delay. (A) The probability of triggering a Ca2+ spark under normal diastolic conditions ([Ca2+]myo = 0.096 µM, and [Ca2+]nsr = 1.02 mM). (B) The Ca2+ spark triggering delay under normal conditions. (C) The probability of triggering a Ca2+ spark under high cytosolic calcium ([Ca2+]myo = 0.4 µM, and [Ca2+]nsr = 1.02 mM). (D) The Ca2+ spark triggering delay under high cytosolic calcium.

• 

  Figure 6. The probability of one CRU triggering the neighboring one at different distances and the delay. (A) The probability of triggering a spark at the second CRU at different distances from an activated CRU under the high overload condition ([Ca2+] = 0.156 µM, and [Ca2+]nsr = 1.70 mM). (B) The delay in triggering the second CRU under that high overload condition. (C) The probability of triggering a spark at the second CRU at different distances from an activated CRU under the low overload condition ([Ca2+]myo = 0.156 µM, and [Ca2+]nsr = 1.30 mM). (D) The delay in triggering the second CRU under that low overload condition.

• 

  Figure 7. (A) A simulated calcium spark. (B) Free calcium shows the underlying structure of the release site (the delayed activation of the two satellite clusters are invisible under fluorescence profile. (C) The profile of a calcium spark giving FWHM = 1.85 um (each color represents the snapshot at different time points after the peak (e.g., bk = 0 means the black line at 0 ms is delayed). (D) The free calcium profile using back-calculation method agrees with experimental estimates, however, it underestimates the real free myoplasmic calcium amplitude.

• 

  Figure 8. The probability of one CRU triggering the neighboring CRU at (A) different distance and (B) the delay. The overload conditions ([Ca2+]myo = 0.156 µM, and [Ca2+]nsr = 1.30 mM) where each CRU has 3 satellite clusters of 10 RyR2s, each at distance 0.2 µm were used.

• 

  Figure 9. (A) Simulation geometry to study spark-induced waves (X = CRU location). Black X = the activated CRU, Blue X = the CRU to be activated by diffusing Ca2+. (B) The proposed simulation geometry with an intermediate cluster (Green X = intermediate RyR2 cluster).

• 

  Figure 10. (A) The Po,trigger of Ca2+ release from 9 activated CRUs on one Z-line on the CRU at different distances. (B) The time delay for the activation at different distances.

• 

  Figure 11. (A) The Po,trigger of Ca2+ release from 9 activated CRUs on one Z-line, with 1 intermediate RyR2 cluster in the middle, on the CRU at the next Z-line of different distance. (B) The time delay for the activation at different distances.

• 

  Figure 12. (A) Calcium overload ([Ca]nsr = 1.7 mM, [Ca]i = 0.15 μM), this computational model of the rat ventricular myocyte can reproduce a repetitive sustained calcium wave which typically initiates at one end of the cell. The initiation site typically occurs where release sites are closer together or at a boundary. (B) Calcium waves in the rat ventricular myocyte under [Ca]o = 5 mM overload experimental conditions, the repetitive waves occur at a particular site for each cell. This suggests that there are more density of release sites surrounding the region that allows mass calcium release high enough to trigger the wave. Some waves can sustain to the next end, while some decay and stop in between, which suggests a stochastic nature of the waves (personal communication from Brian Hagen).