Calcium activation of ryanodine receptor channels--reconciling RyR gating models with tetrameric channel structure

Elementary intracellular Ca signals approximated as a transition of release channel system from a metastable state

Cardiac muscle contraction is initiated by an elementary Ca signal (called Ca spark) which is achieved by collective action of Ca release channels in a cluster. The mechanism of this synchronization remains uncertain. We approached Ca spark activation as an emergent phenomenon of an interactive system of release channels. We constructed a weakly lumped Markov chain that applies an Ising model formalism to such release channel clusters and probable open channel configurations and demonstrated that spark activation is described as a system transition from a metastable to an absorbing state, analogous to the pressure required to overcome surface tension in bubble formation. This yielded quantitative estimates of the spark generation probability as a function of various system parameters. We performed numerical simulations to find spark probabilities as a function of sarcoplasmic reticulum Ca concentration, obtaining similar values for spark activation threshold as our analytic model, as well as those reported in experimental studies. Our parametric sensitivity analyses also showed that the spark activation threshold decreased as Ca sensitivity of RyR activation and RyR cluster size increased.

How does flecainide impact RyR2 channel function?

Flecainide, a cardiac class 1C blocker of the surface membrane sodium channel (NaV1.5), has also been reported to reduce cardiac ryanodine receptor (RyR2)-mediated sarcoplasmic reticulum (SR) Ca2+ release. It has been introduced as a clinical antiarrhythmic agent for catecholaminergic polymorphic ventricular tachycardia (CPVT), a condition most commonly associated with gain-of-function RyR2 mutations. Current debate concerns both cellular mechanisms of its antiarrhythmic action and molecular mechanisms of its RyR2 actions. At the cellular level, it targets NaV1.5, RyR2, Na+/Ca2+ exchange (NCX), and additional proteins involved in excitation-contraction (EC) coupling and potentially contribute to the CPVT phenotype. This Viewpoint primarily addresses the various direct molecular actions of flecainide on isolated RyR2 channels in artificial lipid bilayers. Such studies demonstrate different, multifarious, flecainide binding sites on RyR2, with voltage-dependent binding in the channel pore or voltage-independent binding at distant peripheral sites. In contrast to its single NaV1.5 pore binding site, flecainide may bind to at least four separate inhibitory sites on RyR2 and one activation site. None of these binding sites have been specifically located in the linear RyR2 sequence or high-resolution structure. Furthermore, it is not clear which of the inhibitory sites contribute to flecainide's reduction of spontaneous Ca2+ release in cellular studies. A confounding observation is that flecainide binding to voltage-dependent inhibition sites reduces cation fluxes in a direction opposite to physiological Ca2+ flow from SR lumen to cytosol. This may suggest that, rather than directly blocking Ca2+ efflux, flecainide can reduce Ca2+ efflux by blocking counter currents through the pore which otherwise limit SR membrane potential change during systolic Ca2+ efflux. In summary, the antiarrhythmic effects of flecainide in CPVT seem to involve multiple components of EC coupling and multiple actions on RyR2. Their clarification may identify novel specific drug targets and facilitate flecainide's clinical utilization in CPVT.

Magnesium Ions Moderate Calcium-Induced Calcium Release in Cardiac Calcium Release Sites by Binding to Ryanodine Receptor Activation and Inhibition Sites

Ryanodine receptor channels at calcium release sites of cardiac myocytes operate on the principle of calcium-induced calcium release. In vitro experiments revealed competition of Ca2+ and Mg2+ in the activation of ryanodine receptors (RyRs) as well as inhibition of RyRs by Mg2+. The impact of RyR modulation by Mg2+ on calcium release is not well understood due to the technical limitations of in situ experiments. We turned instead to an in silico model of a calcium release site (CRS), based on a homotetrameric model of RyR gating with kinetic parameters determined from in vitro measurements. We inspected changes in the activity of the CRS model in response to a random opening of one of 20 realistically distributed RyRs, arising from Ca2+/Mg2+ interactions at RyR channels. Calcium release events (CREs) were simulated at a range of Mg2+-binding parameters at near-physiological Mg2+ and ATP concentrations. Facilitation of Mg2+ binding to the RyR activation site inhibited the formation of sparks and slowed down their activation. Impeding Mg-binding to the RyR activation site enhanced spark formation and speeded up their activation. Varying Mg2+ binding to the RyR inhibition site also dramatically affected calcium release events. Facilitation of Mg2+ binding to the RyR inhibition site reduced the amplitude, relative occurrence, and the time-to-end of sparks, and vice versa. The characteristics of CREs correlated dose-dependently with the effective coupling strength between RyRs, defined as a function of RyR vicinity, single-channel calcium current, and Mg-binding parameters of the RyR channels. These findings postulate the role of Mg2+ in calcium release as a negative modulator of the coupling strength among RyRs in a CRS, translating to damping of the positive feedback of the calcium-induced calcium-release mechanism.

In silico simulations reveal that RYR distribution affects the dynamics of calcium release in cardiac myocytes

The dyads of cardiac myocytes contain ryanodine receptors (RYRs) that generate calcium sparks upon activation. To test how geometric factors of RYR distribution contribute to the formation of calcium sparks, which cannot be addressed experimentally, we performed in silico simulations on a large set of models of calcium release sites (CRSs). Our models covered the observed range of RYR number, density, and spatial arrangement. The calcium release function of CRSs was modeled by RYR openings, with an open probability dependent on concentrations of free Ca2+ and Mg2+ ions, in a rapidly buffered system, with a constant open RYR calcium current. We found that simulations of spontaneous sparks by repeatedly opening one of the RYRs in a CRS produced three different types of calcium release events (CREs) in any of the models. Transformation of simulated CREs into fluorescence signals yielded calcium sparks with characteristics close to the observed ones. CRE occurrence varied broadly with the spatial distribution of RYRs in the CRS but did not consistently correlate with RYR number, surface density, or calcium current. However, it correlated with RYR coupling strength, defined as the weighted product of RYR vicinity and calcium current, so that CRE characteristics of all models followed the same state-response function. This finding revealed the synergy between structure and function of CRSs in shaping dyad function. Lastly, rearrangements of RYRs simulating hypothetical experiments on splitting and compaction of a dyad revealed an increased propensity to generate spontaneous sparks and an overall increase in calcium release in smaller and more compact dyads, thus underlying the importance and physiological role of RYR arrangement in cardiac myocytes.

Therapeutic Strategies to Target Calcium Dysregulation in Alzheimer's Disease

Alzheimer’s disease (AD) is the most common form of dementia, affecting millions of people worldwide. Unfortunately, none of the current treatments are effective at improving cognitive function in AD patients and, therefore, there is an urgent need for the development of new therapies that target the early cause(s) of AD. Intracellular calcium (Ca²⁺) regulation is critical for proper cellular and neuronal function. It has been suggested that Ca²⁺ dyshomeostasis is an upstream factor of many neurodegenerative diseases, including AD. For this reason, chemical agents or small molecules aimed at targeting or correcting this Ca²⁺ dysregulation might serve as therapeutic strategies to prevent the development of AD. Moreover, neurons are not alone in exhibiting Ca²⁺ dyshomeostasis, since Ca²⁺ disruption is observed in other cell types in the brain in AD. In this review, we examine the distinct Ca²⁺ channels and compartments involved in the disease mechanisms that could be potential targets in AD.

Structure and Function of the Human Ryanodine Receptors and Their Association with Myopathies-Present State, Challenges, and Perspectives

Cardiac arrhythmias are serious, life-threatening diseases associated with the dysregulation of Ca2+ influx into the cytoplasm of cardiomyocytes. This dysregulation often arises from dysfunction of ryanodine receptor 2 (RyR2), the principal Ca2+ release channel. Dysfunction of RyR1, the skeletal muscle isoform, also results in less severe, but also potentially life-threatening syndromes. The RYR2 and RYR1 genes have been found to harbor three main mutation “hot spots”, where mutations change the channel structure, its interdomain interface properties, its interactions with its binding partners, or its dynamics. In all cases, the result is a defective release of Ca2+ ions from the sarcoplasmic reticulum into the myocyte cytoplasm. Here, we provide an overview of the most frequent diseases resulting from mutations to RyR1 and RyR2, briefly review some of the recent experimental structural work on these two molecules, detail some of the computational work describing their dynamics, and summarize the known changes to the structure and function of these receptors with particular emphasis on their N-terminal, central, and channel domains.

Unbalance Between Sarcoplasmic Reticulum Ca2 + Uptake and Release: A First Step Toward Ca2 + Triggered Arrhythmias and Cardiac Damage

The present review focusses on the regulation and interplay of cardiac SR Ca²⁺ handling proteins involved in SR Ca²⁺ uptake and release, i.e., SERCa2/PLN and RyR2. Both RyR2 and SERCA2a/PLN are highly regulated by post-translational modifications and/or different partners' proteins. These control mechanisms guarantee a precise equilibrium between SR Ca²⁺ reuptake and release. The review then discusses how disruption of this balance alters SR Ca²⁺ handling and may constitute a first step toward cardiac damage and malignant arrhythmias. In the last part of the review, this concept is exemplified in different cardiac diseases, like prediabetic and diabetic cardiomyopathy, digitalis intoxication and ischemia-reperfusion injury.

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.

Roles for the Endoplasmic Reticulum in Regulation of Neuronal Calcium Homeostasis

By influencing Ca²⁺ homeostasis in spatially and architecturally distinct neuronal compartments, the endoplasmic reticulum (ER) illustrates the notion that form and function are intimately related. The contribution of ER to neuronal Ca²⁺ homeostasis is attributed to the organelle being the largest reservoir of intracellular Ca²⁺ and having a high density of Ca²⁺ channels and transporters. As such, ER Ca²⁺ has incontrovertible roles in the regulation of axodendritic growth and morphology, synaptic vesicle release, and neural activity dependent gene expression, synaptic plasticity, and mitochondrial bioenergetics. Not surprisingly, many neurological diseases arise from ER Ca²⁺ dyshomeostasis, either directly due to alterations in ER resident proteins, or indirectly via processes that are coupled to the regulators of ER Ca²⁺ dynamics. In this review, we describe the mechanisms involved in the establishment of ER Ca²⁺ homeostasis in neurons. We elaborate upon how changes in the spatiotemporal dynamics of Ca²⁺ exchange between the ER and other organelles sculpt neuronal function and provide examples that demonstrate the involvement of ER Ca²⁺ dyshomeostasis in a range of neurological and neurodegenerative diseases.

Imaging and analysis of genetically encoded calcium indicators linking neural circuits and behaviors

Confirming the direct link between neural circuit activity and animal behavior has been a principal aim of neuroscience. The genetically encoded calcium indicator (GECI), which binds to calcium ions and emits fluorescence visualizing intracellular calcium concentration, enables detection of in vivo neuronal firing activity. Various GECIs have been developed and can be chosen for diverse purposes. These GECI-based signals can be acquired by several tools including two-photon microscopy and microendoscopy for precise or wide imaging at cellular to synaptic levels. In addition, the images from GECI signals can be analyzed with open source codes including constrained non-negative matrix factorization for endoscopy data (CNMF_E) and miniscope 1-photon-based calcium imaging signal extraction pipeline (MIN1PIPE), and considering parameters of the imaged brain regions (e.g., diameter or shape of soma or the resolution of recorded images), the real-time activity of each cell can be acquired and linked with animal behaviors. As a result, GECI signal analysis can be a powerful tool for revealing the functions of neuronal circuits related to specific behaviors.

Stochastic and deterministic approaches to modelling calcium release in cardiac myocytes at different spatial arrangements of ryanodine receptors

Calcium release sites (CRSs) play a key role in excitation-contraction coupling of cardiac myocytes. Recent studies based on electron tomography and super-resolution imaging revealed that CRSs are not completely filled with ryanodine receptors (RyRs) and that the spatial arrangement of RyRs is neither uniform nor static. In this work, we studied the effect of spatial arrangement of RyRs on RyR activation using simulations based on Monte Carlo (MC) and mean-field (MF) approaches. Both approaches showed that activation of RyRs is sensitive to the arrangement of RyRs in the CRS. However, the MF simulations did not reproduce results of MC simulations for non-compact CRSs, suggesting that the approximations used in the MF approach are not suitable for simulation studies of RyRs arrangements observed experimentally. MC simulations revealed the importance of realistic spatial arrangement of RyRs for adequate modelling of calcium release in cardiac myocytes.

Drosophila CaV2 channels harboring human migraine mutations cause synapse hyperexcitability that can be suppressed by inhibition of a Ca2+ store release pathway

Gain-of-function mutations in the human CaV2.1 gene CACNA1A cause familial hemiplegic migraine type 1 (FHM1). To characterize cellular problems potentially triggered by CaV2.1 gains of function, we engineered mutations encoding FHM1 amino-acid substitutions S218L (SL) and R192Q (RQ) into transgenes of Drosophila melanogaster CaV2/cacophony. We expressed the transgenes pan-neuronally. Phenotypes were mild for RQ-expressing animals. By contrast, single mutant SL- and complex allele RQ,SL-expressing animals showed overt phenotypes, including sharply decreased viability. By electrophysiology, SL- and RQ,SL-expressing neuromuscular junctions (NMJs) exhibited enhanced evoked discharges, supernumerary discharges, and an increase in the amplitudes and frequencies of spontaneous events. Some spontaneous events were gigantic (10-40 mV), multi-quantal events. Gigantic spontaneous events were eliminated by application of TTX-or by lowered or chelated Ca2+-suggesting that gigantic events were elicited by spontaneous nerve firing. A follow-up genetic approach revealed that some neuronal hyperexcitability phenotypes were reversed after knockdown or mutation of Drosophila homologs of phospholipase Cβ (PLCβ), IP3 receptor, or ryanodine receptor (RyR)-all factors known to mediate Ca2+ release from intracellular stores. Pharmacological inhibitors of intracellular Ca2+ store release produced similar effects. Interestingly, however, the decreased viability phenotype was not reversed by genetic impairment of intracellular Ca2+ release factors. On a cellular level, our data suggest inhibition of signaling that triggers intracellular Ca2+ release could counteract hyperexcitability induced by gains of CaV2.1 function.

Regulation of the RyR channel gating by Ca2+ and Mg2

Ryanodine receptors (RyRs) are the Ca2+ 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 Ca2+ 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 Ca2+. RyRs also possess two inhibitory sites in their cytoplasmic domains with Ca2+ affinities of the order of 1 μM and 1 mM. Magnesium competes with Ca2+ at these sites to inhibit RyRs and this plays an important role in modulating their Ca2+-dependent activity in muscle. This review focuses on how these sites lead to RyR modulation by Ca2+ and Mg2+ and how these mechanisms control Ca2+ release in excitation-contraction coupling and cardiac pacemaking.

The N-Terminal Region of the Ryanodine Receptor Affects Channel Activation

Mutations in the cardiac ryanodine receptor (RyR2), the ion channel responsible for release of calcium ions from intracellular stores into cytoplasm, are the cause of several inherited cardiac arrhythmias. At the molecular level, disease symptoms can be mimicked by domain peptides from mutation-prone regions of RyR2 that bind to RyR2 and activate it. Here we show that the domain peptide DP_cpvtN2, corresponding to the central helix of the N-terminal region of RyR2, activates the RyR2 channel. Structural modeling of interaction between DP_cpvtN2 and the N-terminal region of RyR2 in the closed and open conformation provided three plausible structures of the complex. Only one of them could explain the dependence of RyR2 activity on concentration of DP_cpvtN2. The structure of the complex was at odds with the previously proposed “domain switch” mechanism of competition between domain peptides and ryanodine receptor domains. Likewise, in structural models of the N-terminal region, the conformational changes induced by DP_cpvtN2 binding were different from those induced by mutation of central helix amino acids. The activating effect of DP_cpvtN2 binding and of mutations in the central helix could be explained by their similar effect on the transition energy between the closed and open conformation of RyR2.

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

Calcium waves can form and propagate at low frequencies of spontaneous calcium sparks if the calcium dependence of spark frequency is sufficiently steep, or the number of open RyRs is sufficiently large.

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

Modification of cardiac RYR2 gating by a peptide from the central domain of the RYR2

The effect of a domain peptide DPCPVTc from the central region of the RYR2 on ryanodine receptors from rat heart has been examined in planar lipid bilayers. At a zero holding potential and at 8 mmol L−1 luminal Ca2+ concentration, DPCPVTc induced concentrationdependent activation of the ryanodine receptor that led up to 20-fold increase of PO at saturating DPCPVTc concentrations. DPCPVTc prolonged RyR2 openings and increased RyR2 opening frequency. At all peptide concentrations the channels displayed large variability in open probability, open time and frequency of openings. With increasing peptide concentration, the fraction of high open probability records increased together with their open time. The closed times of neither low- nor high-open probability records depended on peptide concentration. The concentration dependence of all gating parameters had EC50 of 20 μmol L−1 and a Hill slope of 2. Comparison of the effects of DPCPVTc with the effects of ATP and cytosolic Ca2+ suggests that activation does not involve luminal feed-through and is not caused by modulation of the cytosolic activation A-site. The data suggest that although “domain unzipping” by DPCPVTc occurs in both modes of RyR activity, it affects RyR gating only when the channel resides in the H-mode of activity.

Luminal Ca2+ controls activation of the cardiac ryanodine receptor by ATP

The synergic effect of luminal Ca²⁺, cytosolic Ca²⁺, and cytosolic adenosine triphosphate (ATP) on activation of cardiac ryanodine receptor (RYR2) channels was examined in planar lipid bilayers. The dose–response of RYR2 gating activity to ATP was characterized at a diastolic cytosolic Ca²⁺ concentration of 100 nM over a range of luminal Ca²⁺ concentrations and, vice versa, at a diastolic luminal Ca²⁺ concentration of 1 mM over a range of cytosolic Ca²⁺ concentrations. Low level of luminal Ca²⁺ (1 mM) significantly increased the affinity of the RYR2 channel for ATP but without substantial activation of the channel. Higher levels of luminal Ca²⁺ (8–53 mM) markedly amplified the effects of ATP on the RYR2 activity by selectively increasing the maximal RYR2 activation by ATP, without affecting the affinity of the channel to ATP. Near-diastolic cytosolic Ca²⁺ levels (<500 nM) greatly amplified the effects of luminal Ca²⁺. Fractional inhibition by cytosolic Mg²⁺ was not affected by luminal Ca²⁺. In models, the effects of luminal and cytosolic Ca²⁺ could be explained by modulation of the allosteric effect of ATP on the RYR2 channel. Our results suggest that luminal Ca²⁺ ions potentiate the RYR2 gating activity in the presence of ATP predominantly by binding to a luminal site with an apparent affinity in the millimolar range, over which local luminal Ca²⁺ likely varies in cardiac myocytes.

Construction of calcium release sites in cardiac myocytes

Local character of calcium release in cardiac myocytes, as defined by confocal recordings of calcium sparks, implies independent activation of individual calcium release sites based on ryanodine receptor (RyR) channel recruitment. We constructed virtual calcium release sites (vCRSs) composed of a variable number of RyR channels distributed in clusters in accordance with the experimentally observed cluster size distribution. The vCRSs consisted either of a single virtual calcium release unit (vCRU), in which all clusters shared a common dyadic space, or of multiple virtual calcium release units (CRUs) containing one cluster each and having separate dyadic spaces. We explored the stochastic behavior of vCRSs to understand the activation and recruitment of RyRs during calcium sparks. RyRs were represented by the published allosteric gating model that included regulation by cytosolic Ca(2+) and Mg(2+). The interaction of Mg(2+) with the RyR Ca(2+)-binding sites and the refractory period of vCRSs were optimized to accord with the experimentally observed calcium dependence of calcium spark frequency. The Mg(2+)-binding parameters of RyRs that provided the best description of spark frequency depended on the number of RyRs assembled in the vCRSs. Adequate inhibitory effect of Mg(2+) on the calcium dependence of RyR open probability was achieved if the vCRSs contained at least three clusters. For the distribution of the number of open RyRs in evoked calcium sparks to correspond to the experimentally observed distribution of spark calcium release fluxes, at least three clusters had to share a common virtual CRU, in which ∼3 RyRs open to form an average spark. These results reconcile the small cluster size and stochastic placement of RyRs in the release sites with the estimates of the amount of RyR protein, volume density of calcium release sites, and the size of calcium release sites in rat cardiac myocytes.

Calcium spike variability in cardiac myocytes results from activation of small cohorts of ryanodine receptor 2 channels

In mammalian cardiac myocytes, the elementary calcium releases triggered by step voltage stimuli manifest either as solitary or as twin spikes that vary widely in kinetics and amplitude for unknown reasons. Here we examined the variability of calcium spikes measured using line-scanning confocal microscopy in patch-clamped rat ventricular myocytes. Amplitude distributions of the single and of the first of twin spikes were broader than those of the second spikes. All could be best approximated by a sum of a few elementary Gaussian probability distribution functions. The latency distributions of the single and the first spikes were identical, much shorter and less variable than those of the second spikes. The multimodal distribution of spike amplitudes and the probability of occurrence of twin spikes were stochastically congruent with activation of only a few of the many RyR2 channels present in the release site cluster. The occurrence of twin release events was rare due to refractoriness of release, induced with a probability proportional to the number of RyR2s activated in the primary release event. We conclude that the variability of the elementary calcium release events supports a calcium signalling mechanism that arises from stochastics of RyR2 gating and from inactivation of local origin.

A mechanistic description of gating of the human cardiac ryanodine receptor in a regulated minimal environment

Cardiac muscle contraction, triggered by the action potential, is mediated by the release of Ca²⁺ from the sarcoplasmic reticulum through ryanodine receptor (RyR)2 channels. In situ, RyR2 gating is modulated by numerous physiological and pharmacological agents, and altered RyR2 function underlies the occurrence of arrhythmias in both inherited and acquired diseases. To understand fully the mechanisms underpinning the regulation of RyR2 in the normal heart and how these systems are altered in pathological conditions, we must first gain a detailed knowledge of the fundamental processes of RyR2 gating. In this investigation, we provide key novel mechanistic insights into the physical reality of RyR2 gating revealed by new experimental and analytical approaches. We have examined in detail the single-channel gating kinetics of the purified human RyR2 when activated by cytosolic Ca²⁺ in a stringently regulated environment where the modulatory influence of factors external to the channel were minimized. The resulting gating schemes are based on an accurate description of single-channel kinetics using hidden Markov model analysis and reveal several novel aspects of RyR2 gating behavior: (a) constitutive gating is observed as unliganded opening events; (b) binding of Ca²⁺ to the channel stabilizes it in different open states; (c) RyR2 exists in two preopening closed conformations in equilibrium, one of which binds Ca²⁺ more readily than the other; (d) the gating of RyR2 when bound to Ca²⁺ can be described by a kinetic scheme incorporating bursts; and (e) analysis of flicker closing events within bursts reveals gating activity that is not influenced by ligand binding. The gating schemes generated in this investigation provide a framework for future studies in which the mechanisms of action of key physiological regulatory factors, disease-linked mutations, and potential therapeutic compounds can be described precisely.

Calcium alternans in a couplon network model of ventricular myocytes: role of sarcoplasmic reticulum load

Intracellular calcium (Ca) alternans in cardiac myocytes have been shown in many experimental studies, and the mechanisms remain incompletely understood. We recently developed a "3R theory" that links Ca sparks to whole cell Ca alternans through three critical properties: randomness of Ca sparks; recruitment of a Ca spark by neighboring Ca sparks; and refractoriness of Ca release units. In this study, we used computer simulation of a physiologically detailed mathematical model of a ventricular myocyte couplon network to study how sarcoplasmic reticulum (SR) Ca load and other physiological parameters, such as ryanodine receptor sensitivity, SR uptake rate, Na-Ca exchange strength, and Ca buffer levels affect Ca alternans in the context of 3R theory. We developed a method to calculate the parameters used in the 3R theory (i.e., the primary spark rate and the recruitment rate) from the physiologically detailed Ca cycling model and paced the model periodically to elicit Ca alternans. We show that alternans only occurs for an intermediate range of the SR Ca load, and the underlying mechanism can be explained via its effects on the 3Rs. Furthermore, we show that altering the physiological parameters not only directly changes the 3Rs but also alters the SR Ca load, having an indirect effect on the 3Rs as well. Therefore, our present study links the SR Ca load and other physiological parameters to whole cell Ca alternans through the framework of the 3R theory, providing a general mechanistic understanding of Ca alternans in ventricular myocytes.

Control of Ca2+ release by action potential configuration in normal and failing murine cardiomyocytes

Cardiomyocytes from failing hearts exhibit spatially nonuniform or dyssynchronous sarcoplasmic reticulum (SR) Ca(2+) 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(2+) release pattern. In this steady-state condition, greater integrated Ca(2+) entry during the CHF AP lead to increased SR Ca(2+) content. A resulting increase in ryanodine receptor sensitivity synchronized SR Ca(2+) release in the mathematical model, thus offsetting the desynchronizing effects of reduced driving force for Ca(2+) entry. A modest nondyssynchronous prolongation of Ca(2+) release was nevertheless observed during the steady-state CHF AP, which contributed to increased time-to-peak measurements for Ca(2+) transients in failing cells. Thus, dyssynchronous Ca(2+) release in failing mouse myocytes does not result from electrical remodeling, but rather other alterations such as T-tubule reorganization.

Frequency and release flux of calcium sparks in rat cardiac myocytes: a relation to RYR gating

Cytosolic calcium concentration in resting cardiac myocytes locally fluctuates as a result of spontaneous microscopic Ca²⁺ releases or abruptly rises as a result of an external trigger. These processes, observed as calcium sparks, are fundamental for proper function of cardiac muscle. In this study, we analyze how the characteristics of spontaneous and triggered calcium sparks are related to cardiac ryanodine receptor (RYR) gating. We show that the frequency of spontaneous sparks and the probability distribution of calcium release flux quanta of triggered sparks correspond quantitatively to predictions of an allosteric homotetrameric model of RYR gating. This model includes competitive binding of Ca²⁺ and Mg²⁺ ions to the RYR activation sites and allosteric interaction between divalent ion binding and channel opening. It turns out that at rest, RYRs are almost fully occupied by Mg²⁺. Therefore, spontaneous sparks are most frequently evoked by random openings of the highly populated but rarely opening Mg₄RYR and CaMg₃RYR forms, whereas triggered sparks are most frequently evoked by random openings of the less populated but much more readily opening Ca₂Mg₂RYR and Ca₃MgRYR forms. In both the spontaneous and the triggered sparks, only a small fraction of RYRs in the calcium release unit manages to open during the spark because of the limited rate of Mg²⁺ unbinding. This mechanism clarifies the unexpectedly low calcium release flux during elementary release events and unifies the theory of calcium signaling in resting and contracting cardiac myocytes.

Luminal Mg2+, a key factor controlling RYR2-mediated Ca2+ release: cytoplasmic and luminal regulation modeled in a tetrameric channel

In cardiac muscle, intracellular Ca(2+) and Mg(2+) are potent regulators of calcium release from the sarcoplasmic reticulum (SR). It is well known that the free [Ca(2+)] in the SR ([Ca(2+)](L)) stimulates the Ca(2+) release channels (ryanodine receptor [RYR]2). However, little is known about the action of luminal Mg(2+), which has not been regarded as an important regulator of Ca(2+) release. The effects of luminal Ca(2+) and Mg(2+) on sheep RYR2 were measured in lipid bilayers. Cytoplasmic and luminal Ca(2+) produced a synergistic increase in the opening rate of RYRs. A novel, high affinity inhibition of RYR2 by luminal Mg(2+) was observed, pointing to an important physiological role for luminal Mg(2+) in cardiac muscle. At diastolic [Ca(2+)](C), luminal Mg(2+) inhibition was voltage independent, with K(i) = 45 microM at luminal [Ca(2+)] ([Ca(2+)](L)) = 100 microM. Luminal and cytoplasmic Mg(2+) inhibition was alleviated by increasing [Ca(2+)](L) or [Ca(2+)](C). Ca(2+) and Mg(2+) on opposite sides of the bilayer exhibited competitive effects on RYRs, indicating that they can compete via the pore for common sites. The data were accurately fitted by a model based on a tetrameric RYR structure with four Ca(2+)-sensing mechanisms on each subunit: activating luminal L-site (40-microM affinity for Mg(2+) and Ca(2+)), cytoplasmic A-site (1.2 microM for Ca(2+) and 60 microM for Mg(2+)), inactivating cytoplasmic I(1)-site (approximately 10 mM for Ca(2+) and Mg(2+)), and I(2)-site (1.2 microM for Ca(2+)). Activation of three or more subunits will cause channel opening. Mg(2+) inhibition occurs primarily by Mg(2+) displacing Ca(2+) from the L- and A-sites, and Mg(2+) fails to open the channel. The model predicts that under physiological conditions, SR load-dependent Ca(2+) release (1) is mainly determined by Ca(2+) displacement of Mg(2+) from the L-site as SR loading increases, and (2) depends on the properties of both luminal and cytoplasmic activation mechanisms.

Local calcium release activation by DHPR calcium channel openings in rat cardiac myocytes

The principal role of calcium current in the triggering of calcium release in cardiac myocytes is well recognized. The mechanism of how calcium current (I(Ca)) controls the intensity of calcium release is not clear because of the stochastic nature of voltage-dependent gating of calcium channels (DHPRs) and of calcium-dependent gating of ryanodine receptors (RyRs). To disclose the relation between DHPR openings and the probability of calcium release, local calcium release activation by I(Ca) was investigated in rat ventricular myocytes using patch-clamp and confocal microscopy. Calcium spikes were activated by temporally synchronized DHPR calcium current triggers, generated by instantaneous 'tail' I(Ca) and modulated by prepulse duration, by tail potential, and by the DHPR agonist BayK 8644. The DHPR-RyR coupling fidelity was determined from the temporal distribution of calcium spike latencies using a model based on exponentially distributed DHPR open times. The analysis provided a DHPR mean open time of approximately 0.5 ms, RyR activation time constant of approximately 0.6 ms, and RyR activation kinetics of the 4th order. The coupling fidelity was low due to the inherent prevalence of very short DHPR openings but was increased when DHPR openings were prolonged by BayK 8644. The probability of calcium release activation was high, despite low coupling fidelity, due to the activation of many DHPRs at individual release sites. We conclude that the control of calcium release intensity by physiological stimuli can be achieved by modulating the number and duration of DHPR openings at low coupling fidelity, thus avoiding the danger of inadvertently triggering calcium release events.

Stochastic binding of Ca2+ ions in the dyadic cleft; continuous versus random walk description of diffusion

Ca(2+) signaling in the dyadic cleft in ventricular myocytes is fundamentally discrete and stochastic. We study the stochastic binding of single Ca(2+) ions to receptors in the cleft using two different models of diffusion: a stochastic and discrete Random Walk (RW) model, and a deterministic continuous model. We investigate whether the latter model, together with a stochastic receptor model, can reproduce binding events registered in fully stochastic RW simulations. By evaluating the continuous model goodness-of-fit for a large range of parameters, we present evidence that it can. Further, we show that the large fluctuations in binding rate observed at the level of single time-steps are integrated and smoothed at the larger timescale of binding events, which explains the continuous model goodness-of-fit. With these results we demonstrate that the stochasticity and discreteness of the Ca(2+) signaling in the dyadic cleft, determined by single binding events, can be described using a deterministic model of Ca(2+) diffusion together with a stochastic model of the binding events, for a specific range of physiological relevant parameters. Time-consuming RW simulations can thus be avoided. We also present a new analytical model of bimolecular binding probabilities, which we use in the RW simulations and the statistical analysis.

Quantification of the effects of a ryanodine receptor channel mutation on interaction with a ryanoid

Understanding the nature of the interaction of the plant alkaloid ryanodine with its receptor channel (RyR) is important to aid interpretation of physiological studies and provide structure-function information about RyR. We present here the first quantitative description of the relative single-channel kinetic effects of a single-point mutation in RyR2. We exploit the well-characterized ryanoid 8beta-amino-9alpha-hydroxyryanodine that displays reversible kinetics with RyR2. We explicitly show that the effect of the Q4863A mutation is to increase the apparent dissociation constant by increasing the apparent dissociation rate of the ryanoid. The voltage-dependence of the interaction displays no change. We infer that Q4863 is not involved with the voltage-drop but is able to influence ryanoid-bound structural changes. We discuss structural mechanisms by which this mutation could affect ryanoid interaction.

A domain peptide of the cardiac ryanodine receptor regulates channel sensitivity to luminal Ca2+ via cytoplasmic Ca2+ sites

The clustering of cardiac RyR mutations, linked to sudden cardiac death (SCD), into several regions in the amino acid sequence underlies the hypothesis that these mutations interfere with stabilising interactions between different domains of the RyR2. SCD mutations cause increased channel sensitivity to cytoplasmic and luminal Ca(2+). A synthetic peptide corresponding to part of the central domain (DPc10:(2460)G-P(2495)) was designed to destabilise the interaction of the N-terminal and central domains of wild-type RyR2 and mimic the effects of SCD mutations. With Ca(2+) as the sole regulating ion, DPc10 caused increased channel activity which could be reversed by removal of the peptide whereas in the presence of ATP DPc10 caused no activation. In support of the domain destablising hypothesis, the corresponding peptide (DPc10-mut) containing the CPVT mutation R2474S did not affect channel activity under any circumstances. DPc10-induced activation was due to a small increase in RyR2 sensitivity to cytoplasmic Ca(2+) and a large increase in the magnitude of luminal Ca(2+) activation. The increase in the luminal Ca(2+) response appeared reliant on the luminal-to-cytoplasmic Ca(2+) flux in the channel, indicating that luminal Ca(2+) was activating the RyR2 via its cytoplasmic Ca(2+) sites. DPc10 had no significant effect on the RyR2 gating associated with luminal Ca(2+) sensing sites. The results were fitted by the luminal-triggered Ca(2+) feed-through model and the effects of DPc10 were explained entirely by perturbations in cytoplasmic Ca(2+)-activation mechanism.

A kinetic model of single and clustered IP3 receptors in the absence of Ca2+ feedback

Ca2+ liberation through inositol 1,4,5-trisphosphate receptor (IP3R) channels generates complex patterns of spatiotemporal cellular Ca2+ signals owing to the biphasic modulation of channel gating by Ca2+ itself. These processes have been extensively studied in Xenopus oocytes, where imaging studies have revealed local Ca2+ signals ("puffs") arising from clusters of IP3R, and patch-clamp studies on isolated oocyte nuclei have yielded extensive data on IP3R gating kinetics. To bridge these two levels of experimental data, we developed an IP3R model and applied stochastic simulation and transition matrix theory to predict the behavior of individual and clustered IP3R channels. The channel model consists of four identical, independent subunits, each of which has an IP3-binding site together with one activating and one inactivating Ca2+-binding site. The channel opens when at least three subunits undergo a conformational change to an "active" state after binding IP3 and Ca2+. The model successfully reproduces patch-clamp data; including the dependence of open probability, mean open duration, and mean closed duration on [IP3] and [Ca2+]. Notably, the biexponential distribution of open-time duration and the dependence of mean open time on [Ca2+] are explained by populations of openings involving either three or four active subunits. As a first step toward applying the single IP3R model to describe cellular responses, we then simulated measurements of puff latency after step increases of [IP3]. Assuming that stochastic opening of a single IP3R at basal cytosolic [Ca2+] and any given [IP3] has a high probability of rapidly triggering neighboring channels by calcium-induced calcium release to evoke a puff, optimal correspondence with experimental data of puff latencies after photorelease of IP3 was obtained when the cluster contained a total of 40-70 IP3Rs.

Ca2+ stores regulate ryanodine receptor Ca2+ release channels via luminal and cytosolic Ca2+ sites

The free [Ca2+] in endoplasmic/sarcoplasmic reticulum Ca2+ stores regulates excitability of Ca2+ release by stimulating the Ca2+ release channels. Just how the stored Ca2+ regulates activation of these channels is still disputed. One proposal attributes luminal Ca2+-activation to luminal facing regulatory sites, whereas another envisages Ca2+ permeation to cytoplasmic sites. This study develops a unified model for luminal Ca2+ activation for single cardiac ryanodine receptors (RyR2) and RyRs in coupled clusters in artificial lipid bilayers. It is shown that luminal regulation of RyR2 involves three modes of action associated with Ca2+ sensors in different parts of the molecule; a luminal activation site (L-site, 60 microM affinity), a cytoplasmic activation site (A-site, 0.9 microM affinity), and a novel cytoplasmic inactivation site (I2-site, 1.2 microM affinity). RyR activation by luminal Ca2+ is demonstrated to occur by a multistep process dubbed luminal-triggered Ca2+ feedthrough. Ca2+ binding to the L-site initiates brief openings (1 ms duration at 1-10 s(-1)) allowing luminal Ca2+ to access the A-site, producing up to 30-fold prolongation of openings. The model explains a broad data set, reconciles previous conflicting observations and provides a foundation for understanding the action of pharmacological agents, RyR-associated proteins, and RyR2 mutations on a range of Ca2+-mediated physiological and pathological processes.

Kinetic properties of the cardiac L-type Ca2+ channel and its role in myocyte electrophysiology: a theoretical investigation

The L-type Ca(2+) channel (Ca(V)1.2) plays an important role in action potential (AP) generation, morphology, and duration (APD) and is the primary source of triggering Ca(2+) for the initiation of Ca(2+)-induced Ca(2+)-release in cardiac myocytes. In this article we present: 1), a detailed kinetic model of Ca(V)1.2, which is incorporated into a model of the ventricular mycoyte where it interacts with a kinetic model of the ryanodine receptor in a restricted subcellular space; 2), evaluation of the contribution of voltage-dependent inactivation (VDI) and Ca(2+)-dependent inactivation (CDI) to total inactivation of Ca(V)1.2; and 3), description of dynamic Ca(V)1.2 and ryanodine receptor channel-state occupancy during the AP. Results are: 1), the Ca(V)1.2 model reproduces experimental single-channel and macroscopic-current data; 2), the model reproduces rate dependence of APD, [Na(+)](i), and the Ca(2+)-transient (CaT), and restitution of APD and CaT during premature stimuli; 3), CDI of Ca(V)1.2 is sensitive to Ca(2+) that enters the subspace through the channel and from SR release. The relative contributions of these Ca(2+) sources to total CDI during the AP vary with time after depolarization, switching from early SR dominance to late Ca(V)1.2 dominance. 4), The relative contribution of CDI to total inactivation of Ca(V)1.2 is greater at negative potentials, when VDI is weak; and 5), loss of VDI due to the Ca(V)1.2 mutation G406R (linked to the Timothy syndrome) results in APD prolongation and increased CaT.