Nanoscale organization of ryanodine receptor distribution and phosphorylation pattern determines the dynamics of calcium sparks

Super-resolution imaging techniques have provided a better understanding of the relationship between the nanoscale organization and function of ryanodine receptors (RyRs) in cardiomyocytes. Recent data have indicated that this relationship is disrupted in heart failure (HF), as RyRs are dispersed into smaller and more numerous clusters. However, RyRs are also hyperphosphorylated in this condition, and this is reported to occur preferentially within the cluster centre. Thus, the combined impact of RyR relocalization and sensitization on Ca²⁺ spark generation in failing cardiomyocytes is likely complex and these observations suggest that both the nanoscale organization of RyRs and the pattern of phosphorylated RyRs within clusters could be critical determinants of Ca²⁺ spark dynamics. To test this hypothesis, we used computational modeling to quantify the relationships between RyR cluster geometry, phosphorylation patterns, and sarcoplasmic reticulum (SR) Ca²⁺ release. We found that RyR cluster disruption results in a decrease in spark fidelity and longer sparks with a lower amplitude. Phosphorylation of some RyRs within the cluster can play a compensatory role, recovering healthy spark dynamics. Interestingly, our model predicts that such compensation is critically dependent on the phosphorylation pattern, as phosphorylation localized within the cluster center resulted in longer Ca²⁺ sparks and higher spark fidelity compared to a uniformly distributed phosphorylation pattern. Our results strongly suggest that both the phosphorylation pattern and nanoscale RyR reorganization are critical determinants of Ca²⁺ dynamics in HF.     

Ca2+ entry-independent effects of L-type Ca2+ channel modulators on Ca2+ sparks in ventricular myocytes

During the cardiac action potential, Ca²⁺ entry through dyhidropyridine receptor L-type Ca²⁺ channels (DHPRs) activates ryanodine receptors (RyRs) Ca²⁺-release channels, resulting in massive Ca²⁺ mobilization from the sarcoplasmic reticulum (SR). This global Ca²⁺ release arises from spatiotemporal summation of many localized elementary Ca²⁺-release events, Ca²⁺ sparks. We tested whether DHPRs modulate Ca²⁺sparks in a Ca²⁺ entry-independent manner. Negative modulation by DHPR of RyRs via physical interactions is accepted in resting skeletal muscle but remains controversial in the heart. Ca²⁺ sparks were studied in cat cardiac myocytes permeabilized with saponin or internally perfused via a patch pipette. Bathing and pipette solutions contained low Ca²⁺ (100 nM). Under these conditions, Ca²⁺ sparks were detected with a stable frequency of 3–5 sparks·s^–1·100 µm^–1. The DHPR blockers nifedipine, nimodipine, FS-2, and calciseptine decreased spark frequency, whereas the DHPR agonists Bay-K8644 and FPL-64176 increased it. None of these agents altered the spatiotemporal characteristics of Ca²⁺ sparks. The DHPR modulators were also without effect on SR Ca²⁺ load (caffeine-induced Ca²⁺ transients) or sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA) activity (Ca²⁺ loading rates of isolated SR microsomes) and did not change cardiac RyR channel gating (planar lipid bilayer experiments). In summary, DHPR modulators affected spark frequency in the absence of DHPR-mediated Ca²⁺ entry. This action could not be attributed to a direct action of DHPR modulators on SERCA or RyRs. Our results suggest that the activity of RyR Ca²⁺-release units in ventricular myocytes is modulated by Ca²⁺ entry-independent conformational changes in neighboring DHPRs. exitation-contraction coupling; ryanodine receptor; sarco(endo)plasmic reticulum Ca²⁺-ATPase; dihydropyridine receptor; sarcoplasmic reticulum     

Single Ryanodine Receptor Channel Basis of Caffeine's Action on Ca Sparks

Caffeine (1, 3, 7-trimethylxanthine) is a widely used pharmacological agonist of the cardiac ryanodine receptor (RyR2) Ca²⁺ release channel. It is also a well-known stimulant that can produce adverse side effects, including arrhythmias. Here, the action of caffeine on single RyR2 channels in bilayers and Ca²⁺ sparks in permeabilized ventricular cardiomyocytes is defined. Single RyR2 caffeine activation depended on the free Ca²⁺ level on both sides of the channel. Cytosolic Ca²⁺ enhanced RyR2 caffeine affinity, whereas luminal Ca²⁺ essentially scaled maximal caffeine activation. Caffeine activated single RyR2 channels in diastolic quasi-cell-like solutions (cytosolic MgATP, pCa 7) with an EC₅₀ of 9.0 ± 0.4 mM. Low-dose caffeine (0.15 mM) increased Ca²⁺ spark frequency ∼75% and single RyR2 opening frequency ∼150%. This implies that not all spontaneous RyR2 openings during diastole are associated with Ca²⁺ sparks. Assuming that only the longest openings evoke sparks, our data suggest that a spark may result only when a spontaneous single RyR2 opening lasts >6 ms.     

Ryanoids and imperatoxin affect the modulation of cardiac ryanodine receptors by dihydropyridine receptor Peptide A

Ca²⁺-entry via L-type Ca²⁺ channels (DHPR) is known to trigger ryanodine receptor (RyR)-mediated Ca²⁺-release from sarcoplasmic reticulum (SR). The mechanism that terminates SR Ca²⁺ release is still unknown. Previous reports showed evidence of Ca²⁺-entry independent inhibition of Ca²⁺ sparks by DHPR in cardiomyocytes. A peptide from the DHPR loop II-III (PepA) was reported to modulate isolated RyRs. We found that PepA induced voltage-dependent “flicker block” and transition to substates of fully-activated cardiac RyRs in planar bilayers. Substates had less voltage-dependence than block and did not represent occupancy of a ryanoid site. However, ryanoids stabilized PepA-induced events while PepA increased RyR2 affinity for ryanodol, which suggests cooperative interactions. Ryanodol stabilized Imperatoxin A (IpTx_A) binding but when IpTx_A bound first, it prevented ryanodol binding. Moreover, IpTx_A and PepA excluded each other from their sites. This suggests that IpTx_A generates a vestibular gate (either sterically or allosterically) that prevents access to the peptides and ryanodol binding sites. Inactivating gate moieties (“ball peptides”) from K⁺ and Na⁺ channels (ShakerB and KIFMK, respectively) induced well resolved slow block and substates, which were sensitive to ryanoids and IpTx_A and allowed, by comparison, better understanding of PepA action. The RyR2 appears to interact with PepA or ball peptides through a two-step mechanism, reminiscent of the inactivation of voltage-gated channels, which includes binding to outer (substates) and inner (block) vestibular regions in the channel conduction pathway. Our results open the possibility that “ball peptide-like” moieties in RyR2-interacting proteins could modulate SR Ca²⁺ release in cells.     

Reconciling depressed Ca2+ sparks occurrence with enhanced RyR2 activity in failing mice cardiomyocytes

Abnormalities in cardiomyocyte Ca²⁺ handling contribute to impaired contractile function in heart failure (HF). Experiments on single ryanodine receptors (RyRs) incorporated into lipid bilayers have indicated that RyRs from failing hearts are more active than those from healthy hearts. Here, we analyzed spontaneous Ca²⁺ sparks (brief, localized increased in [Ca²⁺]_i) to evaluate RyR cluster activity in situ in a mouse post-myocardial infarction (PMI) model of HF. The cardiac ejection fraction of PMI mice was reduced to ∼30% of that of sham-operated (sham) mice, and their cardiomyocytes were hypertrophied. The [Ca²⁺]_i transient amplitude and sarcoplasmic reticulum (SR) Ca²⁺ load were decreased in intact PMI cardiomyocytes compared with those from sham mice, and spontaneous Ca²⁺ sparks were less frequent, whereas the fractional release and the frequency of Ca²⁺ waves were both increased, suggesting higher RyR activity. In permeabilized cardiomyocytes, in which the internal solution can be controlled, Ca²⁺ sparks were more frequent in PMI cells (under conditions of similar SR Ca²⁺ load), confirming the enhanced RyR activity. However, in intact cells from PMI mice, the Ca²⁺ sparks frequency normalized by the SR Ca²⁺ load in that cell were reduced compared with those in sham mice, indicating that the cytosolic environment in intact cells contributes to the decrease in Ca²⁺ spark frequency. Indeed, using an internal “failing solution” with less ATP (as found in HF), we observed a dramatic decrease in Ca²⁺ spark frequency in permeabilized PMI and sham myocytes. In conclusion, our data show that, even if isolated RyR channels show more activity in HF, concomitant alterations in intracellular media composition and SR Ca²⁺ load may mask these effects at the Ca²⁺ spark level in intact cells. Nonetheless, in this scenario, the probability of arrhythmogenic Ca²⁺ waves is enhanced, and they play a potential role in the increase in arrhythmia events in HF patients.     

Control of Sarcoplasmic Reticulum Ca2+ Release by Stochastic RyR Gating within a 3D Model of the Cardiac Dyad and Importance of Induction Decay for CICR Termination

The factors responsible for the regulation of regenerative calcium-induced calcium release (CICR) during Ca²⁺ spark evolution remain unclear. Cardiac ryanodine receptor (RyR) gating in rats and sheep was recorded at physiological Ca²⁺, Mg²⁺, and ATP levels and incorporated into a 3D model of the cardiac dyad, which reproduced the time course of Ca²⁺ sparks, Ca²⁺ blinks, and Ca²⁺ spark restitution. The termination of CICR by induction decay in the model principally arose from the steep Ca²⁺ dependence of RyR closed time, with the measured sarcoplasmic reticulum (SR) lumen Ca²⁺ dependence of RyR gating making almost no contribution. The start of CICR termination was strongly dependent on the extent of local depletion of junctional SR Ca²⁺, as well as the time course of local Ca²⁺ gradients within the junctional space. Reducing the dimensions of the dyad junction reduced Ca²⁺ spark amplitude by reducing the strength of regenerative feedback within CICR. A refractory period for Ca²⁺ spark initiation and subsequent Ca²⁺ spark amplitude restitution arose from 1), the extent to which the regenerative phase of CICR can be supported by the partially depleted junctional SR, and 2), the availability of releasable Ca²⁺ in the junctional SR. The physical organization of RyRs within the junctional space had minimal effects on Ca²⁺ spark amplitude when more than nine RyRs were present. Spark amplitude had a nonlinear dependence on RyR single-channel Ca²⁺ flux, and was approximately halved by reducing the flux from 0.6 to 0.2 pA. Although rat and sheep RyRs had quite different Ca²⁺ sensitivities, Ca²⁺ spark amplitude was hardly affected. This suggests that moderate changes in RyR gating by second-messenger systems will principally alter the spatiotemporal properties of SR release, with smaller effects on the amount released.     

Spatial and temporal Ca, Mg, and ATP dynamics in cardiac dyads during calcium release

We have constructed a three-dimensional reaction-diffusion model of the mammalian cardiac calcium release unit. We analyzed effects of diffusion coefficients, single channel current amplitude, density of RyR channels, and reaction kinetics of ATP²⁻ with Ca²⁺ and Mg²⁺ ions on spatiotemporal concentration profiles of Ca²⁺, Mg²⁺, and ATP²⁻ in the dyadic cleft during Ca²⁺ release. The model revealed that Ca²⁺ concentration gradients persist near RyRs in the steady state. Even with low number of open RyRs, peak [Ca²⁺] in the dyadic space reached values similar to estimates of luminal [Ca²⁺] in ∼1 ms, suggesting that during calcium release the Ca²⁺ gradient moves from the cisternal membrane towards the boundary of the dyadic space with the cytosol. The released Ca²⁺ bound to ATP²⁻, and thus substantially decreased ATP²⁻ concentration in the dyadic space. The released Ca²⁺ could also replace Mg²⁺ in its complex with ATP²⁻ during first milliseconds of release if dissociation of MgATP was fast. The results suggest that concentration changes of Ca²⁺, Mg²⁺, and ATP²⁻ might be large and fast enough to reduce dyadic RyR activity. Thus, under physiological conditions, termination of calcium release may be facilitated by the synergic effect of the construction and chemistry of mammalian cardiac dyads.     

Enhanced ryanodine receptor-mediated calcium leak determines reduced sarcoplasmic reticulum calcium content in chronic canine heart failure

In this study, we investigated the role of elevated sarcoplasmic reticulum (SR) Ca²⁺ leak through ryanodine receptors (RyR2s) in heart failure (HF)-related abnormalities of intracellular Ca²⁺ handling, using a canine model of chronic HF. The cytosolic Ca²⁺ transients were reduced in amplitude and slowed in duration in HF myocytes compared with control, changes paralleled by a dramatic reduction in the total SR Ca²⁺ content. Direct measurements of [Ca²⁺]_SR in both intact and permeabilized cardiac myocytes demonstrated that SR luminal [Ca²⁺] is markedly lowered in HF, suggesting that alterations in Ca²⁺ transport rather than fractional SR volume reduction accounts for the diminished Ca²⁺ release capacity of SR in HF. SR Ca²⁺ ATPase (SERCA2)-mediated SR Ca²⁺ uptake rate was not significantly altered, and Na⁺/Ca²⁺ exchange activity was accelerated in HF myocytes. At the same time, SR Ca²⁺ leak, measured directly as a loss of [Ca²⁺]_SR after inhibition of SERCA2 by thapsigargin, was markedly enhanced in HF myocytes. Moreover, the reduced [Ca²⁺]_SR in HF myocytes could be nearly completely restored by the RyR2 channel blocker ruthenium red. The effects of HF on cytosolic and SR luminal Ca²⁺ signals could be reasonably well mimicked by the RyR2 channel agonist caffeine. Taken together, these results suggest that RyR2-mediated SR Ca²⁺ leak is a major factor in the abnormal intracellular Ca²⁺ handling that critically contributes to the reduced SR Ca²⁺ content of failing cardiomyocytes.     

Increased calcium leak associated with reduced calsequestrin expression in hyperthyroid cardiomyocytes

IntroductionCalcium (Ca²⁺) leak during cardiac diastole is chiefly mediated by intracellular Ca²⁺ channel/Ryanodine Receptors. Increased diastolic Ca²⁺ leak has been proposed as the mechanism underlying the appearance of hereditary arrhythmias. However, little is known about alterations in diastolic Ca²⁺ leak and the specific roles played by key intracellular Ca²⁺-handling proteins in hyperthyroidism, a known arrhythmogenic condition.AimWe sought to determine whether there were modifications in diastolic Ca²⁺ leak, based on the recording of Ca²⁺ sparks and Ca²⁺ waves; we also investigated changes in the expression and activity of key Ca²⁺ handling proteins, including ryanodine receptors, Sarco-Endoplasmic Reticulum Ca²⁺ ATPase pump and calsequestrin in isolated left-ventricular cardiomyocytes isolated from hyperthyroid rats.Materials and methodsElectrocardiography (ECG) recordings were performed in control and hyperthyroid rats. Ca²⁺ sparks, Ca²⁺ waves, and electrically-stimulated Ca²⁺ transients were recorded in Fluo-3-loaded cardiomyocytes from both experimental groups using confocal microscopy. In addition, left-ventricular homogenates and Ryanodine Receptor-enriched membrane fractions were prepared for assessing [³H]-ryanodine binding, hydrolytic ATPase activity of SERCA pump and expression levels of key proteins by Western blot, and cDNA for real-time qPCR.Results and conclusionsExtrasystoles were observed in hearts of hyperthyroid rats by ECG recordings. Arrhythmogenic activity, high incidence of Ca²⁺ waves, and de novo Ca²⁺ wavelets −in the absence of sarcoplasmic reticulum Ca²⁺ overload- were recorded in these cardiomyocytes. The exacerbated diastolic Ca²⁺ leak and arrhythmogenic activities were related to a diminished expression of calsequestrin along with increased SERCA pump activity, which, in effect, promoted a gain-of-function in RyRs without alterations in SR Ca²⁺ load, RyR expression or its Ca²⁺ sensitivity.     

A calcium-induced calcium release mechanism mediated by calsequestrin

Calcium (Ca²⁺)-induced Ca²⁺ release (CICR) is widely accepted as the principal mechanism linking electrical excitation and mechanical contraction in cardiac cells. The CICR mechanism has been understood mainly based on binding of cytosolic Ca²⁺ with ryanodine receptors (RyRs) and inducing Ca²⁺ release from the sarcoplasmic reticulum (SR). However, recent experiments suggest that SR lumenal Ca²⁺ may also participate in regulating RyR gating through calsequestrin (CSQ), the SR lumenal Ca²⁺ buffer. We investigate how SR Ca²⁺ release via RyR is regulated by Ca²⁺ and calsequestrin (CSQ). First, a mathematical model of RyR kinetics is derived based on experimental evidence. We assume that the RyR has three binding sites, two cytosolic sites for Ca²⁺ activation and inactivation, and one SR lumenal site for CSQ binding. The open probability (P_o) of the RyR is found by simulation under controlled cytosolic and SR lumenal Ca²⁺. Both peak and steady-state P_o effectively increase as SR lumenal Ca²⁺ increases. Second, we incorporate the RyR model into a CICR model that has both a diadic space and the junctional SR (jSR). At low jSR Ca²⁺ loads, CSQs are more likely to bind with the RyR and act to inhibit jSR Ca²⁺ release, while at high SR loads CSQs are more likely to detach from the RyR, thereby increasing jSR Ca²⁺ release. Furthermore, this CICR model produces a nonlinear relationship between fractional jSR Ca²⁺ release and jSR load. These findings agree with experimental observations in lipid bilayers and cardiac myocytes.     

Structural and functional properties of ryanodine receptor type 3 in zebrafish tail muscle

The ryanodine receptor (RyR)1 isoform of the sarcoplasmic reticulum (SR) Ca²⁺ release channel is an essential component of all skeletal muscle fibers. RyR1s are detectable as “junctional feet” (JF) in the gap between the SR and the plasmalemma or T-tubules, and they are required for excitation–contraction (EC) coupling and differentiation. A second isoform, RyR3, does not sustain EC coupling and differentiation in the absence of RyR1 and is expressed at highly variable levels. Anatomically, RyR3 expression correlates with the presence of parajunctional feet (PJF), which are located on the sides of the SR junctional cisternae in an arrangement found only in fibers expressing RyR3. In frog muscle fibers, the presence of RyR3 and PJF correlates with the occurrence of Ca²⁺ sparks, which are elementary SR Ca²⁺ release events of the EC coupling machinery. Here, we explored the structural and functional roles of RyR3 by injecting zebrafish (Danio rerio) one-cell stage embryos with a morpholino designed to specifically silence RyR3 expression. In zebrafish larvae at 72 h postfertilization, fast-twitch fibers from wild-type (WT) tail muscles had abundant PJF. Silencing resulted in a drop of the PJF/JF ratio, from 0.79 in WT fibers to 0.03 in the morphants. The frequency with which Ca²⁺ sparks were detected dropped correspondingly, from 0.083 to 0.001 sarcomere⁻¹ s⁻¹. The few Ca²⁺ sparks detected in morphant fibers were smaller in amplitude, duration, and spatial extent compared with those in WT fibers. Despite the almost complete disappearance of PJF and Ca²⁺ sparks in morphant fibers, these fibers looked structurally normal and the swimming behavior of the larvae was not affected. This paper provides important evidence that RyR3 is the main constituent of the PJF and is the main contributor to the SR Ca²⁺ flux underlying Ca²⁺ sparks detected in fully differentiated frog and fish fibers.     

Coupled calcium release channels and their regulation by luminal and cytosolic ions

Contraction in skeletal and cardiac muscle occurs when Ca²⁺ is released from the sarcoplasmic reticulum (SR) through ryanodine receptor (RyR) Ca²⁺ release channels. Several isoforms of the RyR exist throughout the animal kingdom, which are modulated by ATP, Ca²⁺ and Mg²⁺ in the cytoplasm and by Ca²⁺ in the lumen of the SR. This review brings to light recent findings on their mechanisms of action in the mammalian isoforms RyR-1 and RyR-2 with an emphasis on RyR-1 from skeletal muscle. Cytoplasmic Mg²⁺ is a potent RyR antagonist that binds to two classes of cytoplasmic site, identified as low-affinity, non-specific inhibition sites and high-affinity Ca²⁺ activation sites (A-sites). Mg²⁺ inhibition at the A-sites is very sensitive to the cytoplasmic and luminal milieu. Cytoplasmic Ca²⁺, Mg²⁺ and monovalent cations compete for the A-sites. In isolated RyRs, luminal Ca²⁺ alters the Mg²⁺ affinity of the A-site by an allosteric mechanism mediated by luminal sites. However, in close-packed RyR arrays luminal Ca²⁺ can also compete with cytoplasmic ions for the A-site. Activation of RyRs by luminal Ca²⁺ has been attributed to either Ca²⁺ feedthrough to A-sites or to Ca²⁺ regulatory sites on the luminal side of the RyR. As yet there is no consensus on just how luminal Ca²⁺ alters RyR activation. Recent evidence indicates that both mechanisms operate and are likely to be important. Allosteric regulation of A-site Mg²⁺ affinity could trigger Ca²⁺ release, which is reinforced by Ca²⁺ feedthrough.     

Sarcoplasmic Reticulum K (TRIC) Channel Does Not Carry Essential Countercurrent during Ca Release

The charge translocation associated with sarcoplasmic reticulum (SR) Ca²⁺ efflux is compensated for by a simultaneous SR K⁺ influx. This influx is essential because, with no countercurrent, the SR membrane potential (V_m) would quickly (<1 ms) reach the Ca²⁺ equilibrium potential and SR Ca²⁺ release would cease. The SR K⁺ trimeric intracellular cation (TRIC) channel has been proposed to carry the essential countercurrent. However, the ryanodine receptor (RyR) itself also carries a substantial K⁺ countercurrent during release. To better define the physiological role of the SR K⁺ channel, we compared SR Ca²⁺ transport in saponin-permeabilized cardiomyocytes before and after limiting SR K⁺ channel function. Specifically, we reduced SR K⁺ channel conduction 35 and 88% by replacing cytosolic K⁺ for Na⁺ or Cs⁺ (respectively), changes that have little effect on RyR function. Calcium sparks, SR Ca²⁺ reloading, and caffeine-evoked Ca²⁺ release amplitude (and rate) were unaffected by these ionic changes. Our results show that countercurrent carried by SR K⁺ (TRIC) channels is not required to support SR Ca²⁺ release (or uptake). Because K⁺ enters the SR through RyRs during release, the SR K⁺ (TRIC) channel most likely is needed to restore trans-SR K⁺ balance after RyRs close, assuring SR V_m stays near 0 mV.     

Sarcoplasmic Reticulum K+ (TRIC) Channel Does Not Carry Essential Countercurrent during Ca2+ Release

The charge translocation associated with sarcoplasmic reticulum (SR) Ca²⁺ efflux is compensated for by a simultaneous SR K⁺ influx. This influx is essential because, with no countercurrent, the SR membrane potential (V_m) would quickly (<1 ms) reach the Ca²⁺ equilibrium potential and SR Ca²⁺ release would cease. The SR K⁺ trimeric intracellular cation (TRIC) channel has been proposed to carry the essential countercurrent. However, the ryanodine receptor (RyR) itself also carries a substantial K⁺ countercurrent during release. To better define the physiological role of the SR K⁺ channel, we compared SR Ca²⁺ transport in saponin-permeabilized cardiomyocytes before and after limiting SR K⁺ channel function. Specifically, we reduced SR K⁺ channel conduction 35 and 88% by replacing cytosolic K⁺ for Na⁺ or Cs⁺ (respectively), changes that have little effect on RyR function. Calcium sparks, SR Ca²⁺ reloading, and caffeine-evoked Ca²⁺ release amplitude (and rate) were unaffected by these ionic changes. Our results show that countercurrent carried by SR K⁺ (TRIC) channels is not required to support SR Ca²⁺ release (or uptake). Because K⁺ enters the SR through RyRs during release, the SR K⁺ (TRIC) channel most likely is needed to restore trans-SR K⁺ balance after RyRs close, assuring SR V_m stays near 0 mV.     

Recruiting RyRs to Open in a Ca2+ Release Unit: Single-RyR Gating Properties Make RyR Group Dynamics

In cardiac myocytes, clusters of type-2 ryanodine receptors (RyR2s) release Ca²⁺ from the sarcoplasmic reticulum (SR) via a positive feedback mechanism in which fluxed Ca²⁺ activates nearby RyRs. Although the general principles of this are understood, less is known about how single-RyR gating properties define the RyR group dynamics in an array of many channels. Here, we examine this using simulations with three models of RyR gating that have identical open probabilities: the commonly used two-state Markov gating model, one that utilizes multiple exponentials to fit single-channel open time (OT) and closed time (CT) distributions, and an extension of this multiexponential model that also includes experimentally measured correlations between single-channel OTs and CTs. The simulations of RyR clusters that utilize the multiexponential gating model produce infrequent Ca²⁺ release events with relatively few open RyRs. Ca²⁺ release events become even smaller when OT/CT correlations are included. This occurs because the correlations produce a small but consistent bias against recruiting more RyRs to open during the middle of a Ca²⁺ release event, between the initiation and termination phases (which are unaltered compared to the uncorrelated simulations). In comparison, the two-state model produces frequent, large, and long Ca²⁺ release events because it had a recruitment bias in favor of opening more RyRs. This difference stems from the two-state model’s single-RyR OT and CT distributions being qualitatively different from the experimental ones. Thus, the details of single-RyR gating can profoundly affect SR Ca²⁺ release even if open probability and mean OTs and CTs are identical. We also show that Ca²⁺ release events can terminate spontaneously without any reduction in SR [Ca²⁺], luminal regulation, Ca²⁺-dependent inactivation, or physical coupling between RyRs when Ca²⁺ flux is below a threshold value. This supports and extends the pernicious attrition/induction decay hypothesis that SR Ca²⁺ release events terminate below a threshold Ca²⁺ flux.     

Ryanodine receptor allosteric coupling and the dynamics of calcium sparks

Puffs and sparks are localized intracellular Ca²⁺ elevations that arise from the cooperative activity of Ca²⁺-regulated inositol 1,4,5-trisphosphate receptors and ryanodine receptors clustered at Ca²⁺ release sites on the surface of the endoplasmic reticulum or the sarcoplasmic reticulum. While the synchronous gating of Ca²⁺-regulated Ca²⁺ channels can be mediated entirely though the buffered diffusion of intracellular Ca²⁺, interprotein allosteric interactions also contribute to the dynamics of ryanodine receptor (RyR) gating and Ca²⁺ sparks. In this article, Markov chain models of Ca²⁺ release sites are used to investigate how the statistics of Ca²⁺ spark generation and termination are related to the coupling of RyRs via local [Ca²⁺] changes and allosteric interactions. Allosteric interactions are included in a manner that promotes the synchronous gating of channels by stabilizing neighboring closed-closed and/or open-open channel pairs. When the strength of Ca²⁺-mediated channel coupling is systematically varied (e.g., by changing the Ca²⁺ buffer concentration), simulations that include synchronizing allosteric interactions often exhibit more robust Ca²⁺ sparks; however, for some Ca²⁺ coupling strengths the sparks are less robust. We find no evidence that the distribution of spark durations can be used to distinguish between allosteric interactions that stabilize closed channel pairs, open channel pairs, or both in a balanced fashion. On the other hand, the changes in spark duration, interspark interval, and frequency observed when allosteric interactions that stabilize closed channel pairs are gradually removed from simulations are qualitatively different than the changes observed when open or both closed and open channel pairs are stabilized. Thus, our simulations clarify how changes in spark statistics due to pharmacological washout of the accessory proteins mediating allosteric coupling may indicate the type of synchronizing allosteric interactions exhibited by physically coupled RyRs. We also investigate the validity of a mean-field reduction applicable to the dynamics of a ryanodine receptor cluster coupled via local [Ca²⁺] and allosteric interactions. In addition to facilitating parameter studies of the effect of allosteric coupling on spark statistics, the derivation of the mean-field model establishes the correct functional form for cooperativity factors representing the coupled gating of RyRs. This mean-field formulation is well suited for use in computationally efficient whole cell simulations of excitation-contraction coupling.     

The Role of Dyadic Organization in Regulation of Sarcoplasmic Reticulum Ca Handling during Rest in Rabbit Ventricular Myocytes

The dyadic organization of ventricular myocytes ensures synchronized activation of sarcoplasmic reticulum (SR) Ca²⁺ release during systole. However, it remains obscure how the dyadic organization affects SR Ca²⁺ handling during diastole. By measuring intraluminal SR Ca²⁺ ([Ca²⁺]_SR) decline during rest in rabbit ventricular myocytes, we found that ∼76% of leaked SR Ca²⁺ is extruded from the cytosol and only ∼24% is pumped back into the SR. Thus, the majority of Ca²⁺ that leaks from the SR is removed from the cytosol before it can be sequestered back into the SR by the SR Ca²⁺-ATPase (SERCA). Detubulation decreased [Ca²⁺]_SR decline during rest, thus making the leaked SR Ca²⁺ more accessible for SERCA. These results suggest that Ca²⁺ extrusion systems are localized in T-tubules. Inhibition of Na⁺-Ca²⁺ exchanger (NCX) slowed [Ca²⁺]_SR decline during rest by threefold, however did not prevent it. Depolarization of mitochondrial membrane potential during NCX inhibition completely prevented the rest-dependent [Ca²⁺]_SR decline. Despite a significant SR Ca²⁺ leak, Ca²⁺ sparks were very rare events in control conditions. NCX inhibition or detubulation increased Ca²⁺ spark activity independent of SR Ca²⁺ load. Overall, these results indicate that during rest NCX effectively competes with SERCA for cytosolic Ca²⁺ that leaks from the SR. This can be explained if the majority of SR Ca²⁺ leak occurs through ryanodine receptors in the junctional SR that are located closely to NCX in the dyadic cleft. Such control of the dyadic [Ca²⁺] by NCX play a critical role in suppressing Ca²⁺ sparks during rest.     

Intracellular calcium release channels mediate their own countercurrent: the ryanodine receptor case study

Intracellular calcium release channels like ryanodine receptors (RyRs) and inositol trisphosphate receptors (IP₃Rs) mediate large Ca²⁺ release events from Ca²⁺ storage organelles lasting >5 ms. To have such long-lasting Ca²⁺ efflux, a countercurrent of other ions is necessary to prevent the membrane potential from becoming the Ca²⁺ Nernst potential in <1 ms. A recent model of ion permeation through a single, open RyR channel is used here to show that the vast majority of this countercurrent is conducted by the RyR itself. Consequently, changes in membrane potential are minimized locally and instantly, assuring maintenance of a Ca²⁺-driving force. This RyR autocountercurrent is possible because of the poor Ca²⁺ selectivity and high conductance for both monovalent and divalent cations of these channels. The model shows that, under physiological conditions, the autocountercurrent clamps the membrane potential near 0 mV within ∼150 μs. Consistent with experiments, the model shows how RyR unit Ca²⁺ current is defined by luminal [Ca²⁺], permeable ion composition and concentration, and pore selectivity and conductance. This very likely is true of the highly homologous pore of the IP₃R channel.     

Sarcoplasmic Reticulum Structure and Functional Properties that Promote Long-Lasting Calcium Sparks

Calcium (Ca) sparks are the fundamental sarcoplasmic reticulum (SR) Ca release events in cardiac myocytes, and they have a typical duration of 20–40 ms. However, when a fraction of ryanodine receptors (RyRs) are blocked by tetracaine or ruthenium red, Ca sparks lasting hundreds of milliseconds have been observed experimentally. The fundamental mechanism underlying these extremely prolonged Ca sparks is not understood. In this study, we use a physiologically detailed mathematical model of subcellular Ca cycling to examine how Ca spark duration is influenced by the number of functional RyRs in a junctional cluster (which is reduced by tetracaine or ruthenium red) and other SR Ca handling properties. One RyR cluster contains a few to several hundred RyRs, and we use a four-state Markov RyR gating model. Each RyR opens stochastically and is regulated by cytosolic and luminal Ca. We varied the number of functional RyRs in the single cluster, diffusion within the SR network, diffusion between network and junctional SR, cytosolic Ca diffusion, SERCA uptake activity, and RyR open probability. For long-lasting Ca release events, opening events within the cluster must occur continuously because the typical open time of the RyR is only a few milliseconds. We found the following: 1) if the number of RyRs is too small, it is difficult to maintain consecutive openings and stochastic attrition terminates the release; 2) if the number of RyRs is too large, the depletion of Ca from the junctional SR terminates the release; and 3) very long release events require relatively small-sized RyR clusters (reducing flux as seen experimentally with tetracaine) and sufficiently rapid intra-SR Ca diffusion, such that local junctional intra-SR [Ca] can be maintained by intra-SR diffusion and overall SR Ca reuptake.     

Superresolution Modeling of Calcium Release in the Heart

Stable calcium-induced calcium release (CICR) is critical for maintaining normal cellular contraction during cardiac excitation-contraction coupling. The fundamental element of CICR in the heart is the calcium (Ca²⁺) spark, which arises from a cluster of ryanodine receptors (RyR). Opening of these RyR clusters is triggered to produce a local, regenerative release of Ca²⁺ from the sarcoplasmic reticulum (SR). The Ca²⁺ leak out of the SR is an important process for cellular Ca²⁺ management, and it is critically influenced by spark fidelity, i.e., the probability that a spontaneous RyR opening triggers a Ca²⁺ spark. Here, we present a detailed, three-dimensional model of a cardiac Ca²⁺ release unit that incorporates diffusion, intracellular buffering systems, and stochastically gated ion channels. The model exhibits realistic Ca²⁺ sparks and robust Ca²⁺ spark termination across a wide range of geometries and conditions. Furthermore, the model captures the details of Ca²⁺ spark and nonspark-based SR Ca²⁺ leak, and it produces normal excitation-contraction coupling gain. We show that SR luminal Ca²⁺-dependent regulation of the RyR is not critical for spark termination, but it can explain the exponential rise in the SR Ca²⁺ leak-load relationship demonstrated in previous experimental work. Perturbations to subspace dimensions, which have been observed in experimental models of disease, strongly alter Ca²⁺ spark dynamics. In addition, we find that the structure of RyR clusters also influences Ca²⁺ release properties due to variations in inter-RyR coupling via local subspace Ca²⁺ concentration ([Ca²⁺]_ss). These results are illustrated for RyR clusters based on super-resolution stimulated emission depletion microscopy. Finally, we present a believed-novel approach by which the spark fidelity of a RyR cluster can be predicted from structural information of the cluster using the maximum eigenvalue of its adjacency matrix. These results provide critical insights into CICR dynamics in heart, under normal and pathological conditions.