Ca Sparks Do Not Explain all Ryanodine Receptor-Mediated SR Ca Leak in Mouse Ventricular Myocytes

Diastolic Ca leak from the sarcoplasmic reticulum (SR) of ventricular myocytes reduces the SR Ca content, stabilizing the activity of the SR Ca release channel ryanodine receptor for the next beat. SR Ca leak has been visualized globally using whole-cell fluorescence, or locally using confocal microscopy, but never both ways. When using confocal microscopy, leak is imaged as “Ca sparks,” which are fluorescent objects generated by the local reaction-diffusion of released Ca and cytosolic indicator. Here, we used confocal microscopy and simultaneously measured the global ryanodine-receptor-mediated leak rate (J_leak) and Ca sparks in intact mouse ventricular myocytes. We found that spark frequency and J_leak are correlated, as expected if both are manifestations of a common phenomenon. However, we also found that sparks explain approximately half of J_leak. Our strategy unmasks the presence of a subresolution (i.e., nonspark) release of potential physiological relevance.     

Intracellular Ca alternans: coordinated regulation by sarcoplasmic reticulum release, uptake, and leak

Beat-to-beat alternation in the cardiac intracellular Ca (Ca_i) transient can drive action potential (AP) duration alternans, creating a highly arrhythmogenic substrate. Although a steep dependence of fractional sarcoplasmic reticulum (SR) Ca release on SR Ca load has been shown experimentally to promote Ca_i alternans, theoretical studies predict that other factors are also important. Here we present an iterated map analysis of the coordinated effects of SR Ca release, uptake, and leak on the onset of Ca_i alternans. Predictions were compared to numerical simulations using a physiologically realistic AP model as well as to AP clamp experiments in isolated patch-clamped rabbit ventricular myocytes exposed to 1), the Ca channel agonist BayK8644 (100 nM) to increase SR Ca load and release fraction, 2), overexpression of an adenoviral SERCA2a construct to increase SR Ca uptake, and 3), low-dose FK506 (20 μM) or ryanodine (1 μM) to increase SR Ca leak. Our findings show that SR Ca release, uptake, and leak all have independent direct effects that promote (release and leak) or suppress (uptake) Ca_i alternans. However, since each factor affects the other by altering SR Ca load, the net balance of their direct and indirect effects determines whether they promote or suppress alternans. Thus, BayK8644 promotes, whereas Ad-SERCA2a overexpression, ryanodine, and FK506 suppress, Ca_i alternans under AP clamp conditions.     

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.     

How Does Stochastic Ryanodine Receptor-Mediated Ca Leak Fail to Initiate a Ca Spark?

Spontaneous calcium (Ca) sparks are initiated by single ryanodine receptor (RyR) opening. Once one RyR channel opens, it elevates local [Ca] in the cleft space ([Ca]_Cleft), which opens other RyR channels in the same Ca release unit (CaRU) via Ca-induced Ca-release. Experiments by Zima et al. (J. Physiol. 588:4743–4757, 2010) demonstrate that spontaneous Ca sparks occur only when intrasarcoplasmic-reticulum (SR) [Ca] ([Ca]_SR) is above a threshold level, but that RyR-mediated SR Ca leak exists without Ca sparks well below this threshold [Ca]_SR. We examine here how single RyR opening at lower [Ca]_SR can fail to recruit Ca sparks at a CaRU, while still contributing to SR Ca leak. We assess this using a physiologically detailed mathematical model of junctional SR Ca release in which RyR gating is regulated by [Ca]_SR and [Ca]_Cleft. We find that several factors contribute to the failure of Ca sparks as [Ca]_SR declines: 1), lower [Ca]_SR reduces driving force and thus limits local [Ca]_Cleft achieved and the rate of rise during RyR opening; 2), low [Ca]_SR limits RyR open time (τ_O), which further reduces local [Ca]_Cleft attained; 3), low τ_O and fast [Ca]_Cleft dissipation after RyR closure shorten the opportunity for neighboring RyR activation; 4), at low [Ca]_SR, the RyR exhibits reduced [Ca]_Cleft sensitivity. We conclude that all of these factors conspire to reduce the probability of Ca sparks as [Ca]_SR declines, despite continued RyR-mediated Ca leak. In addition, these same factors explain the much lower efficacy of L-type Ca channel opening to trigger local SR Ca release at low [Ca]_SR during excitation-contraction coupling. Conversely, all of these factors are fundamentally important for increasing the propensity for pro-arrhythmic Ca sparks and waves in cardiac myocytes at high [Ca]_SR.     

Predicting Local SR Ca Dynamics during Ca Wave Propagation in Ventricular Myocytes

Of the many ongoing controversies regarding the workings of the sarcoplasmic reticulum (SR) in cardiac myocytes, two unresolved and interconnected topics are 1), mechanisms of calcium (Ca²⁺) wave propagation, and 2), speed of Ca²⁺ diffusion within the SR. Ca²⁺ waves are initiated when a spontaneous local SR Ca²⁺ release event triggers additional release from neighboring clusters of SR release channels (ryanodine receptors (RyRs)). A lack of consensus regarding the effective Ca²⁺ diffusion constant in the SR (D_Ca,SR) severely complicates our understanding of whether dynamic local changes in SR [Ca²⁺] can influence wave propagation. To address this problem, we have implemented a computational model of cytosolic and SR [Ca²⁺] during Ca²⁺ waves. Simulations have investigated how dynamic local changes in SR [Ca²⁺] are influenced by 1), D_Ca,SR; 2), the distance between RyR clusters; 3), partial inhibition or stimulation of SR Ca²⁺ pumps; 4), SR Ca²⁺ pump dependence on cytosolic [Ca²⁺]; and 5), the rate of transfer between network and junctional SR. Of these factors, D_Ca,SR is the primary determinant of how release from one RyR cluster alters SR [Ca²⁺] in nearby regions. Specifically, our results show that local increases in SR [Ca²⁺] ahead of the wave can potentially facilitate Ca²⁺ wave propagation, but only if SR diffusion is relatively slow. These simulations help to delineate what changes in [Ca²⁺] are possible during SR Ca²⁺release, and they broaden our understanding of the regulatory role played by dynamic changes in [Ca²⁺]_SR.     

Regulation of sarcoplasmic reticulum Ca2+ release by cytosolic glutathione in rabbit ventricular myocytes

Of the major cellular antioxidant defenses, glutathione (GSH) is particularly important in maintaining the cytosolic redox potential. Whereas the healthy myocardium is maintained at a highly reduced redox state, it has been proposed that oxidation of GSH can affect the dynamics of Ca²⁺-induced Ca²⁺ release. In this study, we used multiple approaches to define the effects of oxidized glutathione (GSSG) on ryanodine receptor (RyR)-mediated Ca²⁺ release in rabbit ventricular myocytes. To investigate the role of GSSG on sarcoplasmic reticulum (SR) Ca²⁺ release induced by the action potential, we used the thiol-specific oxidant diamide to increase intracellular GSSG in intact myocytes. To more directly assess the effect of GSSG on RyR activity, we introduced GSSG within the cytosol of permeabilized myocytes. RyR-mediated Ca²⁺ release from the SR was significantly enhanced in the presence of GSSG. This resulted in decreased steady-state diastolic [Ca²⁺]_SR, increased SR Ca²⁺ fractional release, and increased spark- and non-spark-mediated SR Ca²⁺ leak. Single-channel recordings from RyR’s incorporated into lipid bilayers revealed that GSSG significantly increased RyR activity. Moreover, oxidation of RyR in the form of intersubunit crosslinking was present in intact myocytes treated with diamide and permeabilized myocytes treated with GSSG. Blocking RyR crosslinking with the alkylating agent N-ethylmaleimide prevented depletion of SR Ca²⁺ load induced by diamide. These findings suggest that elevated cytosolic GSSG enhances SR Ca²⁺ leak due to redox-dependent intersubunit RyR crosslinking. This effect can contribute to abnormal SR Ca²⁺ handling during periods of oxidative stress.     

Control of Ca Release by Action Potential Configuration in Normal and Failing Murine Cardiomyocytes

Cardiomyocytes from failing hearts exhibit spatially nonuniform or dyssynchronous sarcoplasmic reticulum (SR) Ca²⁺ release. We investigated the contribution of action potential (AP) prolongation in mice with congestive heart failure (CHF) after myocardial infarction. AP recordings from CHF and control myocytes were included in a computational model of the dyad, which predicted more dyssynchronous ryanodine receptor opening during stimulation with the CHF AP. This prediction was confirmed in cardiomyocyte experiments, when cells were alternately stimulated by control and CHF AP voltage-clamp waveforms. However, when a train of like APs was used as the voltage stimulus, the control and CHF AP produced a similar Ca²⁺ release pattern. In this steady-state condition, greater integrated Ca²⁺ entry during the CHF AP lead to increased SR Ca²⁺ content. A resulting increase in ryanodine receptor sensitivity synchronized SR Ca²⁺ release in the mathematical model, thus offsetting the desynchronizing effects of reduced driving force for Ca²⁺ entry. A modest nondyssynchronous prolongation of Ca²⁺ release was nevertheless observed during the steady-state CHF AP, which contributed to increased time-to-peak measurements for Ca²⁺ transients in failing cells. Thus, dyssynchronous Ca²⁺ release in failing mouse myocytes does not result from electrical remodeling, but rather other alterations such as T-tubule reorganization.     

β1a490–508, a 19-Residue Peptide from C-Terminal Tail of Cav1.1 β1a Subunit,Potentiates Voltage-Dependent Calcium Release in Adult Skeletal Muscle Fibers

The α1 and β1a subunits of the skeletal muscle calcium channel, Cav1.1, as well as the Ca²⁺ release channel, ryanodine receptor (RyR1), are essential for excitation-contraction coupling. RyR1 channel activity is modulated by the β1a subunit and this effect can be mimicked by a peptide (β1a490–524) corresponding to the 35-residue C-terminal tail of the β1a subunit. Protein-protein interaction assays confirmed a high-affinity interaction between the C-terminal tail of the β1a and RyR1. Based on previous results using overlapping peptides tested on isolated RyR1, we hypothesized that a 19-amino-acid residue peptide (β1a490–508) is sufficient to reproduce activating effects of β1a490–524. Here we examined the effects of β1a490–508 on Ca²⁺ release and Ca²⁺ currents in adult skeletal muscle fibers subjected to voltage-clamp and on RyR1 channel activity after incorporating sarcoplasmic reticulum vesicles into lipid bilayers. β1a490–508 (25 nM) increased the peak Ca²⁺ release flux by 49% in muscle fibers. Considerably fewer activating effects were observed using 6.25, 100, and 400 nM of β1a490–508 in fibers. β1a490–508 also increased RyR1 channel activity in bilayers and Cav1.1 currents in fibers. A scrambled form of β1a490–508 peptide was used as negative control and produced negligible effects on Ca²⁺ release flux and RyR1 activity. Our results show that the β1a490–508 peptide contains molecular components sufficient to modulate excitation-contraction coupling in adult muscle fibers.     

Effects of Partial Sarcoplasmic Reticulum Calcium Depletion on Calcium Release in Frog Cut Muscle Fibers Equilibrated with 20 mM EGTA

Resting sarcoplasmic reticulum (SR) Ca content ([Ca_SR]_R) was varied in cut fibers equilibrated with an internal solution that contained 20 mM EGTA and 0–1.76 mM Ca. SR Ca release and [Ca_SR]_R were measured with the EGTA–phenol red method (Pape et al. 1995. J. Gen. Physiol. 106:259–336). After an action potential, the fractional amount of Ca released from the SR increased from 0.17 to 0.50 when [Ca_SR]_R was reduced from 1,200 to 140 μM. This increase was associated with a prolongation of release (final time constant, from 1–2 to 10–15 ms) and of the action potential (by 1–2 ms). Similar changes in release were observed with brief stimulations to −20 mV in voltage-clamped fibers, in which charge movement (Q _cm) could be measured. The peak values of Q _cm and the fractional rate of SR Ca release, as well as their ON time courses, were little affected by reducing [Ca_SR]_R from 1,200 to 140 μM. After repolarization, however, the OFF time courses of Q _cm and the rate of SR Ca release were slowed by factors of 1.5–1.7 and 6.5, respectively. These and other results suggest that, after action potential stimulation of fibers in normal physiological condition, the increase in myoplasmic free [Ca] that accompanies SR Ca release exerts three negative feedback effects that tend to reduce additional release: (a) the action potential is shortened by current through Ca-activated potassium channels in the surface and/or tubular membranes; (b) the OFF kinetics of Q _cm is accelerated; and (c) Ca inactivation of Ca release is increased. Some of these effects of Ca on an SR Ca channel or its voltage sensor appear to be regulated by the value of [Ca] within 22 nm of the mouth of the channel.     

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.     

The Role of Dyadic Organization in Regulation of Sarcoplasmic Reticulum Ca2+ 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.     

Alterations in ryanodine receptors and related proteins in heart failure

Sarcoplasmic reticulum (SR) Ca^2 + release plays an essential role in mediating cardiac myocyte contraction. Depolarization of the plasma membrane results in influx of Ca^2 + through l-type Ca^2 + channels (LTCCs) that in turn triggers efflux of Ca^2 + from the SR through ryanodine receptor type-2 channels (RyR2). This process known as Ca^2 +-induced Ca^2 +release (CICR) occurs within the dyadic region, where the adjacent transverse (T)-tubules and SR membranes allow RyR2 clusters to release SR Ca^2 + following Ca^2 + influx through adjacent LTCCs. SR Ca^2 + released during systole binds to troponin-C and initiates actin–myosin cross-bridging, leading to muscle contraction. During diastole, the cytosolic Ca^2 + concentration is restored by the resequestration of Ca^2 + into the SR by SR/ER Ca^2 +-ATPase (SERCA2a) and by the extrusion of Ca^2 + via the Na⁺/Ca^2 +-exchanger (NCX1). This whole process, entitled excitation–contraction (EC) coupling, is highly coordinated and determines the force of contraction, providing a link between the electrical and mechanical activities of cardiac muscle. In response to heart failure (HF), the heart undergoes maladaptive changes that result in depressed intracellular Ca^2 + cycling and decreased SR Ca^2 + concentrations. As a result, the amplitude of CICR is reduced resulting in less force production during EC coupling. In this review, we discuss the specific proteins that alter the regulation of Ca^2 + during HF. In particular, we will focus on defects in RyR2-mediated SR Ca^2 + release. This article is part of a Special Issue entitled: Heart failure pathogenesis and emerging diagnostic and therapeutic interventions.     

Association between plasma metal elements and platelet dysfunction in trauma-induced coagulopathy rat model

BackgroundDisorders of metal elements and platelet dysfunction are common in patients with trauma-induced coagulopathy (TIC).AimThe aim of this study was to explore the potential role of plasma metal elements in platelet dysfunction in TIC.MethodsThirty Sprague-Dawley rats were divided into control, hemorrhage shock (HS) and multiple injury (MI) groups. At timepoints of 0.5 and 3 h after trauma and being documented as _{HS 0.5 h}, HS_3 h, _{MI 0.5 h} or MI_3 h, blood samples were harvested for inductively coupled plasma mass spectrometer, conventional coagulation function and thromboelastograph.ResultsThe plasma zinc (Zn), vanadium (V) and cadmium (Ca) decreased initially in HS _0.5 h and recovered slightly in HS_3 h, whereas their plasma concentrations continued to decrease from beginning till MI_3 h (p < 0.05). In HS, plasma Ca, V and nickel were negatively correlated to the time taken to reach the initial formation (R), whereas R was positively correlated to plasms Zn, V, Ca and selenium in MI (p < 0.05). In MI, plasma Ca was positively correlated to maximum amplitude, and plasma V was positively correlated to platelet count (p < 0.05).ConclusionThe plasma concentrations of Zn, V and Ca appeared to contribute to platelet dysfunction in _{HS 0.5 h}, HS_3 h, _{MI 0.5 h} and MI_3 h, which were trauma type sensitive.     

Genetic modulation of the SERCA activity does not affect the Ca leak from the cardiac sarcoplasmic reticulum

The Ca²⁺ content in the sarcoplasmic reticulum (SR) determines the amount of Ca²⁺ released, thereby regulating the magnitude of Ca²⁺ transient and contraction in cardiac muscle. The Ca²⁺ content in the SR is known to be regulated by two factors: the activity of the Ca²⁺ pump (SERCA) and Ca²⁺ leak through the ryanodine receptor (RyR). However, the direct relationship between the SERCA activity and Ca²⁺ leak has not been fully investigated in the heart. In the present study, we evaluated the role of the SERCA activity in Ca²⁺ leak from the SR using a novel saponin-skinned method combined with transgenic mouse models in which the SERCA activity was genetically modulated. In the SERCA overexpression mice, the Ca²⁺ uptake in the SR was significantly increased and the Ca²⁺ transient was markedly increased. However, Ca²⁺ leak from the SR did not change significantly. In mice with overexpression of a negative regulator of SERCA, sarcolipin, the Ca²⁺ uptake by the SR was significantly decreased and the Ca²⁺ transient was markedly decreased. Again, Ca²⁺ leak from the SR did not change significantly. In conclusion, the selective modulation of the SERCA activity modulates Ca²⁺ uptake, although it does not change Ca²⁺ leak from the SR.     

Determination of heat transfer coefficients in plastic French straws plunged in liquid nitrogen

The knowledge of the thermodynamic process during the cooling of reproductive biological systems is important to assess and optimize the cryopreservation procedures. The time–temperature curve of a sample immersed in liquid nitrogen enables the calculation of cooling rates and helps to determine whether it is vitrified or undergoes phase change transition. When dealing with cryogenic liquids, the temperature difference between the solid and the sample is high enough to cause boiling of the liquid, and the sample can undergo different regimes such as film and/or nucleate pool boiling.In the present work, the surface heat transfer coefficients (h) for plastic French straws plunged in liquid nitrogen were determined using the measurement of time–temperature curves. When straws filled with ice were used the cooling curve showed an abrupt slope change which was attributed to the transition of film into nucleate pool boiling regime. The h value that fitted each stage of the cooling process was calculated using a numerical finite element program that solves the heat transfer partial differential equation under transient conditions. In the cooling process corresponding to film boiling regime, the h that best fitted experimental results was h = 148.12 ± 5.4 W/m² K and for nucleate-boiling h = 1355 ± 51 W/m² K. These values were further validated by predicting the time–temperature curve for French straws filled with a biological fluid system (bovine semen-extender) which undergoes freezing. Good agreement was obtained between the experimental and predicted temperature profiles, further confirming the accuracy of the h values previously determined for the ice-filled straw. These coefficients were corroborated using literature correlations.The determination of the boiling regimes that govern the cooling process when plunging straws in liquid nitrogen constitutes an important issue when trying to optimize cryopreservation procedures. Furthermore, this information can lead to improvements in the design of cooling devices in the cryobiology field.     

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.     

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.     

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.     

Ca2+ Overload and Sarcoplasmic Reticulum Instability in tric-a Null Skeletal Muscle

The sarcoplasmic reticulum (SR) of skeletal muscle contains K⁺, Cl⁻, and H⁺ channels may facilitate charge neutralization during Ca²⁺ release. Our recent studies have identified trimeric intracellular cation (TRIC) channels on SR as an essential counter-ion permeability pathway associated with rapid Ca²⁺ release from intracellular stores. Skeletal muscle contains TRIC-A and TRIC-B isoforms as predominant and minor components, respectively. Here we test the physiological function of TRIC-A in skeletal muscle. Biochemical assay revealed abundant expression of TRIC-A relative to the skeletal muscle ryanodine receptor with a molar ratio of TRIC-A/ryanodine receptor ∼5:1. Electron microscopy with the tric-a^−/− skeletal muscle showed Ca²⁺ overload inside the SR with frequent formation of Ca²⁺ deposits compared with the wild type muscle. This elevated SR Ca²⁺ pool in the tric-a^−/− muscle could be released by caffeine, whereas the elemental Ca²⁺ release events, e.g. osmotic stress-induced Ca²⁺ spark activities, were significantly reduced likely reflecting compromised counter-ion movement across the SR. Ex vivo physiological test identified the appearance of “alternan” behavior with isolated tric-a^−/− skeletal muscle, i.e. transient and drastic increase in contractile force appeared within the decreasing force profile during repetitive fatigue stimulation. Inhibition of SR/endoplasmic reticulum Ca^{2+ ATPase} function could lead to aggravation of the stress-induced alternans in the tric-a^−/− muscle. Our data suggests that absence of TRIC-A may lead to Ca²⁺ overload in SR, which in combination with the reduced counter-ion movement may lead to instability of Ca²⁺ movement across the SR membrane. The observed alternan behavior with the tric-a^−/− muscle may reflect a skeletal muscle version of store overload-induced Ca²⁺ release that has been reported in the cardiac muscle under stress conditions.     

Dynamics of calcium sparks and calcium leak in the heart

We present what we believe to be a new mathematical model of Ca²⁺ leak from the sarcoplasmic reticulum (SR) in the heart. To our knowledge, it is the first to incorporate a realistic number of Ca²⁺-release units, each containing a cluster of stochastically gating Ca²⁺ channels (RyRs), whose biophysical properties (e.g., Ca²⁺ sensitivity and allosteric interactions) are informed by the latest molecular investigations. This realistic model allows for the detailed characterization of RyR Ca²⁺-release properties, and shows how this balances reuptake by the SR Ca²⁺ pump. Simulations reveal that SR Ca²⁺ leak consists of brief but frequent single RyR openings (∼3000 cell⁻¹ s⁻¹) that are likely to be experimentally undetectable, and are, therefore, “invisible”. We also observe that these single RyR openings can recruit additional RyRs to open, due to elevated local (Ca²⁺), and occasionally lead to the generation of Ca²⁺ sparks (∼130 cell⁻¹ s⁻¹). Furthermore, this physiological formulation of “invisible” leak allows for the removal of the ad hoc, non-RyR mediated Ca²⁺ leak terms present in prior models. Finally, our model shows how Ca²⁺ sparks can be robustly triggered and terminated under both normal and pathological conditions. Together, these discoveries profoundly influence how we interpret and understand diverse experimental and clinical results from both normal and diseased hearts.