Ca2+ Release through Ryanodine Receptors in the Sarcoplasmic Reticulum and Ca2+ Sequestration by the Mitochondria in Smooth Muscle Cells

The contribution of the Ca²⁺-induced Ca²⁺ release (CICR) mechanism in excitation-contraction (E-C) coupling and the tightness of the coupling between Ca²⁺ influx and Ca²⁺ release are still controversial in smooth muscle cells (SMC). In SMC isolated from the guinea-pig vas deferens or urinary bladder, a depolarizing stimulus initially induced spot-like increases in the intracellular Ca²⁺ concentration ([Ca²⁺] _i ), called “Ca²⁺ hot spots,” at several superficial areas in the cell. When a weak stimulus (a small or a short depolarizing step) was applied, only a few Ca²⁺ hot spots appeared transiently in the superficial area but did not spread into other regions, to trigger global [Ca²⁺] _i rise. Such depolarization-evoked local Ca2+ transients were distinctive from spontaneous Ca²⁺ sparks, since the former were susceptible to Ca²⁺ blockers, ryanodine, and inhibitors of the Ca²⁺ pump in the sarcoplasmic reticulum (SR), suggesting pivotal roles of Ca²⁺ influx through voltage-dependent Ca²⁺ channels (VDCC) and Ca²⁺ release from the SR through ryanodine receptors (RyR) for the activation of Ca²⁺ spots. Frequently discharging Ca²⁺ spark sites (FDS) under resting conditions were located exactly in the same areas as Ca²⁺ hot spots evoked by depolarization, indicating the existence of distinct local junction sites for tight coupling between VDCC in the plasmalemma and RyR in the SR. Co-localization of clusters of RyR and large-conductance Ca²⁺-activated K⁺ (BK) channels was also suggested. The fast and tight coupling for CICR in these junctional sites was triggered also by an action potential, whereas a slower spread of Ca²⁺ wave to the whole-cell areas suggests the loose coupling in propagating CICR to other cell areas. It can therefore be postulated that CICR may occur in two steps upon depolarization; the initial CICR in distinct junctional sites shows tight coupling between Ca²⁺ influx and release, and the following CICR may propagate slow Ca²⁺ waves to other areas. Ryanodine receptors form a multiprotein complex with molecules such as calsequestrin, junctin, triadin, junctophilins, and FK506-binding proteins, which directly or indirectly regulate the RyR activity and the tight coupling. Moreover, an evoked Ca²⁺ spot may enhance Ca²⁺ uptake by neighboring mitochondria and their ATP production to increase energy supply to the Ca²⁺ pump of the SR in the microdomain.     

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.     

Impaired S-Nitrosylation of the Ryanodine Receptor Caused by Xanthine Oxidase Activity Contributes to Calcium Leak in Heart Failure

S-Nitrosylation is a ubiquitous post-translational modification that regulates diverse biologic processes. In skeletal muscle, hypernitrosylation of the ryanodine receptor (RyR) causes sarcoplasmic reticulum (SR) calcium leak, but whether abnormalities of cardiac RyR nitrosylation contribute to dysfunction of cardiac excitation-contraction coupling remains controversial. In this study, we tested the hypothesis that cardiac RyR2 is hyponitrosylated in heart failure, because of nitroso-redox imbalance. We evaluated excitation-contraction coupling and nitroso-redox balance in spontaneously hypertensive heart failure rats with dilated cardiomyopathy and age-matched Wistar-Kyoto rats. Spontaneously hypertensive heart failure myocytes were characterized by depressed contractility, increased diastolic Ca²⁺ leak, hyponitrosylation of RyR2, and enhanced xanthine oxidase derived superoxide. Global S-nitrosylation was decreased in failing hearts compared with nonfailing. Xanthine oxidase inhibition restored global and RyR2 nitrosylation and reversed the diastolic SR Ca²⁺ leak, improving Ca²⁺ handling and contractility. Together these findings demonstrate that nitroso-redox imbalance causes RyR2 oxidation, hyponitrosylation, and SR Ca²⁺ leak, a hallmark of cardiac dysfunction. The reversal of this phenotype by inhibition of xanthine oxidase has important pathophysiologic and therapeutic implications.     

CRISPR/Cas9 Gene editing of RyR2 in human stem cell-derived cardiomyocytes provides a novel approach in investigating dysfunctional Ca2+ signaling

Type-2 ryanodine receptors (RyR2s) play a pivotal role in cardiac excitation-contraction coupling by releasing Ca²⁺ from sarcoplasmic reticulum (SR) via a Ca²⁺ -induced Ca²⁺ release (CICR) mechanism. Two strategies have been used to study the structure-function characteristics of RyR2 and its disease associated mutations: (1) heterologous cell expression of the recombinant mutant RyR2s, and (2) knock-in mouse models harboring RyR2 point mutations. Here, we establish an alternative approach where Ca²⁺ signaling aberrancy caused by the RyR2 mutation is studied in human cardiomyocytes with robust CICR mechanism. Specifically, we introduce point mutations in wild-type RYR2 of human induced pluripotent stem cells (hiPSCs) by CRISPR/Cas9 gene editing, and then differentiate them into cardiomyocytes. To verify the reliability of this approach, we introduced the same disease-associated RyR2 mutation, F2483I, which was studied by us in hiPSC-derived cardiomyocytes (hiPSC-CMs) from a patient biopsy. The gene-edited F2483I hiPSC-CMs exhibited longer and wandering Ca²⁺ sparks, elevated diastolic Ca²⁺ leaks, and smaller SR Ca²⁺ stores, like those of patient-derived cells. Our CRISPR/Cas9 gene editing approach validated the feasibility of creating myocytes expressing the various RyR2 mutants, making comparative mechanistic analysis and pharmacotherapeutic approaches for RyR2 pathologies possible.     

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.     

Oxidative Stress and Ca2+ Release Events in Mouse Cardiomyocytes

Cellular oxidative stress, associated with a variety of common cardiac diseases, is well recognized to affect the function of several key proteins involved in Ca²⁺ signaling and excitation-contraction coupling, which are known to be exquisitely sensitive to reactive oxygen species. These include the Ca²⁺ release channels of the sarcoplasmic reticulum (ryanodine receptors or RyR2s) and the Ca²⁺/calmodulin-dependent protein kinase II (CaMKII). Oxidation of RyR2s was found to increase the open probability of the channel, whereas CaMKII can be activated independent of Ca²⁺ through oxidation. Here, we investigated how oxidative stress affects RyR2 function and SR Ca²⁺ signaling in situ, by analyzing Ca²⁺ sparks in permeabilized mouse cardiomyocytes under a broad range of oxidative conditions. The results show that with increasing oxidative stress Ca²⁺ spark duration is prolonged. In addition, long and very long-lasting (up to hundreds of milliseconds) localized Ca²⁺ release events started to appear, eventually leading to sarcoplasmic reticulum (SR) Ca²⁺ depletion. These changes of release duration could be prevented by the CaMKII inhibitor KN93 and did not occur in mice lacking the CaMKII-specific S2814 phosphorylation site on RyR2. The appearance of long-lasting Ca²⁺ release events was paralleled by an increase of RyR2 oxidation, but also by RyR-S2814 phosphorylation, and by CaMKII oxidation. Our results suggest that in a strongly oxidative environment oxidation-dependent activation of CaMKII leads to RyR2 phosphorylation and thereby contributes to the massive prolongation of SR Ca²⁺ release events.     

Modeling Local X-ROS and Calcium Signaling in the Heart

Stretching single ventricular cardiac myocytes has been shown experimentally to activate transmembrane nicotinamide adenine dinucleotide phosphate oxidase type 2 to produce reactive oxygen species (ROS) and increase the Ca²⁺ spark rate in a process called X-ROS signaling. The increase in Ca²⁺ spark rate is thought to be due to an increase in ryanodine receptor type 2 (RyR2) open probability by direct oxidation of the RyR2 protein complex. In this article, a computational model is used to examine the regulation of ROS and calcium homeostasis by local, subcellular X-ROS signaling and its role in cardiac excitation-contraction coupling. To this end, a four-state RyR2 model was developed that includes an X-ROS-dependent RyR2 mode switch. When activated, [Ca²⁺]_i-sensitive RyR2 open probability increases, and the Ca²⁺ spark rate changes in a manner consistent with experimental observations. This, to our knowledge, new model is used to study the transient effects of diastolic stretching and subsequent ROS production on RyR2 open probability, Ca²⁺ sparks, and the myoplasmic calcium concentration ([Ca²⁺]_i) during excitation-contraction coupling. The model yields several predictions: 1) [ROS] is produced locally near the RyR2 complex during X-ROS signaling and increases by an order of magnitude more than the global ROS signal during myocyte stretching; 2) X-ROS activation just before the action potential, corresponding to ventricular filling during diastole, increases the magnitude of the Ca²⁺ transient; 3) during prolonged stretching, the X-ROS-induced increase in Ca²⁺ spark rate is transient, so that long-sustained stretching does not significantly increase sarcoplasmic reticulum Ca²⁺ leak; and 4) when the chemical reducing capacity of the cell is decreased, activation of X-ROS signaling increases sarcoplasmic reticulum Ca²⁺ leak and contributes to global oxidative stress, thereby increases the possibility of arrhythmia. The model provides quantitative information not currently obtainable through experimental means and thus provides a framework for future X-ROS signaling experiments.     

Mg2+ activates the ryanodine receptor type 2 (RyR2) at intermediate Ca2+ concentrations

To clarify whether activity of the ryanodine receptor type 2 (RyR2) is reduced in the sarcoplasmic reticulum (SR) of cardiac muscle, as is the case with the ryanodine receptor type 1 (RyR1), Ca²⁺-dependent [³H]ryanodine binding, a biochemical measure of Ca²⁺-induced Ca²⁺ release (CICR), was determined using SR vesicle fractions isolated from rabbit and rat cardiac muscles. In the absence of an adenine nucleotide or caffeine, the rat SR showed a complicated Ca²⁺ dependence, instead of the well-documented biphasic dependence of the rabbit SR. In the rat SR, [³H]ryanodine binding initially increased as [Ca²⁺] increased, with a plateau in the range of 10–100 µM Ca²⁺, and thereafter further increased to an apparent peak around 1 mM Ca²⁺, followed by a decrease. In the presence of these modulators, this complicated dependence prevailed, irrespective of the source. Addition of 0.3–1 mM Mg²⁺ unexpectedly increased the binding two- to threefold and enhanced the affinity for [³H]ryanodine at 10–100 µM Ca²⁺, resulting in the well-known biphasic dependence. In other words, the partial suppression of RyR2 is relieved by Mg²⁺. Ca²⁺ could be a substitute for Mg²⁺. Mg²⁺ also amplifies the responses of RyR2 to inhibitory and stimulatory modulators. This stimulating effect of Mg²⁺ on RyR2 is entirely new, and is referred to as the third effect, in addition to the well-known dual inhibitory effects. This effect is critical to describe the role of RyR2 in excitation-contraction coupling of cardiac muscle, in view of the intracellular Mg²⁺ concentration. [³H]ryanodine binding; CICR; stimulation by physiological Mg²⁺, excitation-contraction coupling in the heart     

Two ryanodine receptor isoforms in nonmammalian vertebrate skeletal muscle: possible roles in excitation-contraction coupling and other processes

The ryanodine receptor (RyR) is a Ca²⁺ release channel in the sarcoplasmic reticulum in vertebrate skeletal muscle and plays an important role in excitation–contraction (E–C) coupling. Whereas mammalian skeletal muscle predominantly expresses a single RyR isoform, RyR1, skeletal muscle of many nonmammalian vertebrates expresses equal amounts of two distinct isoforms, α-RyR and β-RyR, which are homologues of mammalian RyR1 and RyR3, respectively. In this review we describe our current understanding of the functions of these two RyR isoforms in nonmammalian vertebrate skeletal muscle. The Ca²⁺ release via the RyR channel can be gated by two distinct modes: depolarization-induced Ca²⁺ release (DICR) and Ca²⁺-induced Ca²⁺ release (CICR). In frog muscle, α-RyR acts as the DICR channel, whereas β-RyR as the CICR channel. However, several lines of evidence suggest that CICR by β-RyR may make only a minor contribution to Ca²⁺ release during E–C coupling. Comparison of frog and mammalian RyR isoforms highlights the marked differences in the patterns of Ca²⁺ release mediated by RyR1 and RyR3 homologues. Interestingly, common features in the Ca²⁺ release patterns are noticed between β-RyR and RyR1. We will discuss possible roles and significance of the two RyR isoforms in E–C coupling and other processes in nonmammalian vertebrate skeletal muscle.     

Dynamic Interreceptor Coupling Contributes to the Consistent Open Duration of Ryanodine Receptors

Ca²⁺ spark is the elementary Ca²⁺ signaling event in muscle excitation-contraction coupling. The rise time of Ca²⁺ spark is rather stable under different conditions, suggesting consistent open duration of ryanodine receptors (RyRs) in vivo. It has been proposed that the array-based behavior of RyRs plays an important role in shaping Ca²⁺ spark dynamics, particularly in controlling the open duration of RyR clusters. Therefore, we investigated the possible role of inter-RyR coupling in stabilization of the open time of arrayed RyRs under several potential perturbations, for instance, array size, inter-RyR coupling noise, and up-regulation or down-regulation of the activity of partial RyRs in the array. We found that RyR arrays with dynamic coupling showed consistent open duration against the perturbations, whereas the RyR array with constant coupling did not. On the other hand, the open probability and amplitude of RyR arrays with dynamic interreceptor coupling were sensitive to the perturbations. These two points were consistent with experimental observations, indicating that the RyR array with dynamic coupling could recapitulate in vivo open properties of RyRs. Our findings support the idea that dynamic coupling is a feasible in vivo working mechanism of RyR arrays.     

A Human Ventricular Myocyte Model with a Refined Representation of Excitation-Contraction Coupling

Cardiac Ca²⁺-induced Ca²⁺ release (CICR) occurs by a regenerative activation of ryanodine receptors (RyRs) within each Ca²⁺-releasing unit, triggered by the activation of L-type Ca²⁺ channels (LCCs). CICR is then terminated, most probably by depletion of Ca²⁺ in the junctional sarcoplasmic reticulum (SR). Hinch et al. previously developed a tightly coupled LCC-RyR mathematical model, known as the Hinch model, that enables simulations to deal with a variety of functional states of whole-cell populations of a Ca²⁺-releasing unit using a personal computer. In this study, we developed a membrane excitation-contraction model of the human ventricular myocyte, which we call the human ventricular cell (HuVEC) model. This model is a hybrid of the most recent HuVEC models and the Hinch model. We modified the Hinch model to reproduce the regenerative activation and termination of CICR. In particular, we removed the inactivated RyR state and separated the single step of RyR activation by LCCs into triggering and regenerative steps. More importantly, we included the experimental measurement of a transient rise in Ca²⁺ concentrations ([Ca²⁺], 10–15 μM) during CICR in the vicinity of Ca²⁺-releasing sites, and thereby calculated the effects of the local Ca²⁺ gradient on CICR as well as membrane excitation. This HuVEC model successfully reconstructed both membrane excitation and key properties of CICR. The time course of CICR evoked by an action potential was accounted for by autonomous changes in an instantaneous equilibrium open probability of couplons. This autonomous time course was driven by a core feedback loop including the pivotal local [Ca²⁺], influenced by a time-dependent decay in the SR Ca²⁺ content during CICR.     

On the Adjacency Matrix of RyR2 Cluster Structures

In the heart, electrical stimulation of cardiac myocytes increases the open probability of sarcolemmal voltage-sensitive Ca²⁺ channels and flux of Ca²⁺ into the cells. This increases Ca²⁺ binding to ligand-gated channels known as ryanodine receptors (RyR2). Their openings cause cell-wide release of Ca²⁺, which in turn causes muscle contraction and the generation of the mechanical force required to pump blood. In resting myocytes, RyR2s can also open spontaneously giving rise to spatially-confined Ca²⁺ release events known as “sparks.” RyR2s are organized in a lattice to form clusters in the junctional sarcoplasmic reticulum membrane. Our recent work has shown that the spatial arrangement of RyR2s within clusters strongly influences the frequency of Ca²⁺ sparks. We showed that the probability of a Ca²⁺ spark occurring when a single RyR2 in the cluster opens spontaneously can be predicted from the precise spatial arrangements of the RyR2s. Thus, “function” follows from “structure.” This probability is related to the maximum eigenvalue (λ ₁) of the adjacency matrix of the RyR2 cluster lattice. In this work, we develop a theoretical framework for understanding this relationship. We present a stochastic contact network model of the Ca²⁺ spark initiation process. We show that λ ₁ determines a stability threshold for the formation of Ca²⁺ sparks in terms of the RyR2 gating transition rates. We recapitulate these results by applying the model to realistic RyR2 cluster structures informed by super-resolution stimulated emission depletion (STED) microscopy. Eigendecomposition of the linearized mean-field contact network model reveals functional subdomains within RyR2 clusters with distinct sensitivities to Ca²⁺. This work provides novel perspectives on the cardiac Ca²⁺ release process and a general method for inferring the functional properties of transmembrane receptor clusters from their structure.     

Distribution and Function of Cardiac Ryanodine Receptor Clusters in Live Ventricular Myocytes

The cardiac Ca²⁺ release channel (ryanodine receptor, RyR2) plays an essential role in excitation-contraction coupling in cardiac muscle cells. Effective and stable excitation-contraction coupling critically depends not only on the expression of RyR2, but also on its distribution. Despite its importance, little is known about the distribution and organization of RyR2 in living cells. To study the distribution of RyR2 in living cardiomyocytes, we generated a knock-in mouse model expressing a GFP-tagged RyR2 (GFP-RyR2). Confocal imaging of live ventricular myocytes isolated from the GFP-RyR2 mouse heart revealed clusters of GFP-RyR2 organized in rows with a striated pattern. Similar organization of GFP-RyR2 clusters was observed in fixed ventricular myocytes. Immunofluorescence staining with the anti-α-actinin antibody (a z-line marker) showed that nearly all GFP-RyR2 clusters were localized in the z-line zone. There were small regions with dislocated GFP-RyR2 clusters. Interestingly, these same regions also displayed dislocated z-lines. Staining with di-8-ANEPPS revealed that nearly all GFP-RyR2 clusters were co-localized with transverse but not longitudinal tubules, whereas staining with MitoTracker Red showed that GFP-RyR2 clusters were not co-localized with mitochondria in live ventricular myocytes. We also found GFP-RyR2 clusters interspersed between z-lines only at the periphery of live ventricular myocytes. Simultaneous detection of GFP-RyR2 clusters and Ca²⁺ sparks showed that Ca²⁺ sparks originated exclusively from RyR2 clusters. Ca²⁺ sparks from RyR2 clusters induced no detectable changes in mitochondrial Ca²⁺ level. These results reveal, for the first time, the distribution of RyR2 clusters and its functional correlation in living ventricular myocytes.     

Dihydropyridine receptor-ryanodine receptor interactions in skeletal muscle excitation-contraction coupling

Much recent progress has been made in our understanding of the mechanism of sarcoplasmic reticulum Ca²⁺ release in skeletal muscle. Vertebrate skeletal muscle excitation-contraction (E-C) coupling is thought to occur by a mechanical coupling mechanism involving protein-protein interactions that lead to activation of the sarcoplasmic reticulum (SR) ryanodine receptor (RyR)/Ca²⁺ release channel by the voltage-sensing transverse (T–) tubule dihydropyridine receptor (DHPR)/Ca²⁺ channel. In a subsequent step, the released Ca²⁺ amplify SR Ca²⁺ release by activating release channels that are not linked to the DHPR. Experiments with mutant muscle cells have indicated that skeletal muscle specific DHPR and RyR isoforms are required for skeletal muscle E-C coupling. A direct functional and structural interaction between a DHPR-derived peptide and the RyR has been described. The interaction between the DHPR and RyR may be stabilized by other proteins such as triadin (a SR junctional protein) and modulated by phosphorylation of the DHPR.     

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.     

Moment closure for local control models of calcium-induced calcium release in cardiac myocytes

In prior work, we introduced a probability density approach to modeling local control of Ca²⁺-induced Ca²⁺ release in cardiac myocytes, where we derived coupled advection-reaction equations for the time-dependent bivariate probability density of subsarcolemmal subspace and junctional sarcoplasmic reticulum (SR) [Ca²⁺] conditioned on Ca²⁺ release unit (CaRU) state. When coupled to ordinary differential equations (ODEs) for the bulk myoplasmic and network SR [Ca²⁺], a realistic but minimal model of cardiac excitation-contraction coupling was produced that avoids the computationally demanding task of resolving spatial aspects of global Ca²⁺ signaling, while accurately representing heterogeneous local Ca²⁺ signals in a population of diadic subspaces and junctional SR depletion domains. Here we introduce a computationally efficient method for simulating such whole cell models when the dynamics of subspace [Ca²⁺] are much faster than those of junctional SR [Ca²⁺]. The method begins with the derivation of a system of ODEs describing the time-evolution of the moments of the univariate probability density functions for junctional SR [Ca²⁺] jointly distributed with CaRU state. This open system of ODEs is then closed using an algebraic relationship that expresses the third moment of junctional SR [Ca²⁺] in terms of the first and second moments. In simulated voltage-clamp protocols using 12-state CaRUs that respond to the dynamics of both subspace and junctional SR [Ca²⁺], this moment-closure approach to simulating local control of excitation-contraction coupling produces high-gain Ca²⁺ release that is graded with changes in membrane potential, a phenomenon not exhibited by common pool models. Benchmark simulations indicate that the moment-closure approach is nearly 10,000-times more computationally efficient than corresponding Monte Carlo simulations while leading to nearly identical results. We conclude by applying the moment-closure approach to study the restitution of Ca²⁺-induced Ca²⁺ release during simulated two-pulse voltage-clamp protocols.     

Size Matters: Ryanodine Receptor Cluster Size Heterogeneity Potentiates Calcium Waves

Ryanodine receptors (RyRs) mediate calcium (Ca)-induced Ca release and intracellular Ca homeostasis. In a cardiac myocyte, RyRs group into clusters of variable size from a few to several hundred RyRs, creating a spatially nonuniform intracellular distribution. It is unclear how heterogeneity of RyR cluster size alters spontaneous sarcoplasmic reticulum (SR) Ca releases (Ca sparks) and arrhythmogenic Ca waves. Here, we tested the impact of heterogeneous RyR cluster size on the initiation of Ca waves. Experimentally, we measured RyR cluster sizes at Ca spark sites in rat ventricular myocytes and further tested functional impacts using a physiologically detailed computational model with spatial and stochastic intracellular Ca dynamics. We found that the spark frequency and amplitude increase nonlinearly with the size of RyR clusters. Larger RyR clusters have lower SR Ca release threshold for local Ca spark initiation and exhibit steeper SR Ca release versus SR Ca load relationship. However, larger RyR clusters tend to lower SR Ca load because of the higher Ca leak rate. Conversely, smaller clusters have a higher threshold and a lower leak, which tends to increase SR Ca load. At the myocyte level, homogeneously large or small RyR clusters limit Ca waves (because of low load for large clusters but low excitability for small clusters). Mixtures of large and small RyR clusters potentiates Ca waves because the enhanced SR Ca load driven by smaller clusters enables Ca wave initiation and propagation from larger RyR clusters. Our study suggests that a spatially heterogeneous distribution of RyR cluster size under pathological conditions may potentiate Ca waves and thus afterdepolarizations and triggered arrhythmias.     

Regulation of the RyR channel gating by Ca<Superscript>2+</Superscript> and Mg<Superscript>2+</Superscript>

Ryanodine receptors (RyRs) are the Ca²⁺ release channels in the sarcoplasmic reticulum in striated muscle which play an important role in excitation-contraction coupling and cardiac pacemaking. Single channel recordings have revealed a wealth of information about ligand regulation of RyRs from mammalian skeletal and cardiac muscle (RyR1 and RyR2, respectively). RyR subunit has a Ca²⁺ activation site located in the luminal and cytoplasmic domains of the RyR. These sites synergistically feed into a common gating mechanism for channel activation by luminal and cytoplasmic Ca²⁺. RyRs also possess two inhibitory sites in their cytoplasmic domains with Ca²⁺ affinities of the order of 1 μM and 1 mM. Magnesium competes with Ca²⁺ at these sites to inhibit RyRs and this plays an important role in modulating their Ca²⁺-dependent activity in muscle. This review focuses on how these sites lead to RyR modulation by Ca²⁺ and Mg²⁺ and how these mechanisms control Ca²⁺ release in excitation-contraction coupling and cardiac pacemaking.     

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.     

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.