Numerical analysis of ryanodine receptor activation by L-type channel activity in the cardiac muscle diad

Biological noise is a key determinant of the reproducibility and adaptability of cardiac pacemaking and EC coupling

Each heartbeat begins with the generation of an action potential in pacemaking cells in the sinoatrial node. This signal triggers contraction of cardiac muscle through a process termed excitation-contraction (EC) coupling. EC coupling is initiated in dyadic structures of cardiac myocytes, where ryanodine receptors in the junctional sarcoplasmic reticulum come into close apposition with clusters of CaV1.2 channels in invaginations of the sarcolemma. Cooperative activation of CaV1.2 channels within these clusters causes a local increase in intracellular Ca2+ that activates the juxtaposed ryanodine receptors. A salient feature of healthy cardiac function is the reliable and precise beat-to-beat pacemaking and amplitude of Ca2+ transients during EC coupling. In this review, we discuss recent discoveries suggesting that the exquisite reproducibility of this system emerges, paradoxically, from high variability at subcellular, cellular, and network levels. This variability is attributable to stochastic fluctuations in ion channel trafficking, clustering, and gating, as well as dyadic structure, which increase intracellular Ca2+ variance during EC coupling. Although the effects of these large, local fluctuations in function and organization are sometimes negligible at the macroscopic level owing to spatial-temporal summation within and across cells in the tissue, recent work suggests that the "noisiness" of these intracellular Ca2+ events may either enhance or counterintuitively reduce variability in a context-dependent manner. Indeed, these noisy events may represent distinct regulatory features in the tuning of cardiac contractility. Collectively, these observations support the importance of incorporating experimentally determined values of Ca2+ variance in all EC coupling models. The high reproducibility of cardiac contraction is a paradoxical outcome of high Ca2+ signaling variability at subcellular, cellular, and network levels caused by stochastic fluctuations in multiple processes in time and space. This underlying stochasticity, which counterintuitively manifests as reliable, consistent Ca2+ transients during EC coupling, also allows for rapid changes in cardiac rhythmicity and contractility in health and disease.

Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going?

Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation-contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.

Ageing Causes Ultrastructural Modification to Calcium Release Units and Mitochondria in Cardiomyocytes

Ageing is associated with an increase in the incidence of heart failure, even if the existence of a real age-related cardiomyopathy remains controversial. Effective contraction and relaxation of cardiomyocytes depend on efficient production of ATP (handled by mitochondria) and on proper Ca2+ supply to myofibrils during excitation-contraction (EC) coupling (handled by Ca2+ release units, CRUs). Here, we analyzed mitochondria and CRUs in hearts of adult (4 months old) and aged (≥24 months old) mice. Analysis by confocal and electron microscopy (CM and EM, respectively) revealed an age-related loss of proper organization and disposition of both mitochondria and EC coupling units: (a) mitochondria are improperly disposed and often damaged (percentage of severely damaged mitochondria: adults 3.5 ± 1.1%; aged 16.5 ± 3.5%); (b) CRUs that are often misoriented (longitudinal) and/or misplaced from the correct position at the Z line. Immunolabeling with antibodies that mark either the SR or T-tubules indicates that in aged cardiomyocytes the sarcotubular system displays an extensive disarray. This disarray could be in part caused by the decreased expression of Cav-3 and JP-2 detected by western blot (WB), two proteins involved in formation of T-tubules and in docking SR to T-tubules in dyads. By WB analysis, we also detected increased levels of 3-NT in whole hearts homogenates of aged mice, a product of nitration of protein tyrosine residues, recognized as marker of oxidative stress. Finally, a detailed EM analysis of CRUs (formed by association of SR with T-tubules) points to ultrastructural modifications, i.e., a decrease in their frequency (adult: 5.1 ± 0.5; aged: 3.9 ± 0.4 n./50 μm2) and size (adult: 362 ± 40 nm; aged: 254 ± 60 nm). The changes in morphology and disposition of mitochondria and CRUs highlighted by our results may underlie an inefficient supply of Ca2+ ions and ATP to the contractile elements, and possibly contribute to cardiac dysfunction in ageing.

Cx43 hemichannel microdomain signaling at the intercalated disc enhances cardiac excitability

Cx43, a major cardiac connexin, forms precursor hemichannels that accrue at the intercalated disc to assemble as gap junctions. While gap junctions are crucial for electrical conduction in the heart, little is known about the potential roles of hemichannels. Recent evidence suggests that inhibiting Cx43 hemichannel opening with Gap19 has antiarrhythmic effects. Here, we used multiple electrophysiology, imaging, and super-resolution techniques to understand and define the conditions underlying Cx43 hemichannel activation in ventricular cardiomyocytes, their contribution to diastolic Ca2+ release from the sarcoplasmic reticulum, and their impact on electrical stability. We showed that Cx43 hemichannels were activated during diastolic Ca2+ release in single ventricular cardiomyocytes and cardiomyocyte cell pairs from mice and pigs. This activation involved Cx43 hemichannel Ca2+ entry and coupling to Ca2+ release microdomains at the intercalated disc, resulting in enhanced Ca2+ dynamics. Hemichannel opening furthermore contributed to delayed afterdepolarizations and triggered action potentials. In single cardiomyocytes, cardiomyocyte cell pairs, and arterially perfused tissue wedges from failing human hearts, increased hemichannel activity contributed to electrical instability compared with nonfailing rejected donor hearts. We conclude that microdomain coupling between Cx43 hemichannels and Ca2+ release is a potentially novel, targetable mechanism of cardiac arrhythmogenesis in heart failure.

The central domain of cardiac ryanodine receptor governs channel activation, regulation, and stability

Structural analyses identified the central domain of ryanodine receptor (RyR) as a transducer converting conformational changes in the cytoplasmic platform to the RyR gate. The central domain is also a regulatory hub encompassing the Ca²⁺-, ATP-, and caffeine-binding sites. However, the role of the central domain in RyR activation and regulation has yet to be defined. Here, we mutated five residues that form the Ca²⁺ activation site and 10 residues with negatively charged or oxygen-containing side chains near the Ca²⁺ activation site. We also generated eight disease-associated mutations within the central domain of RyR2. We determined the effect of these mutations on Ca²⁺, ATP, and caffeine activation and Mg²⁺ inhibition of RyR2. Mutating the Ca²⁺ activation site markedly reduced the sensitivity of RyR2 to Ca²⁺ and caffeine activation. Unexpectedly, Ca²⁺ activation site mutation E3848A substantially enhanced the Ca²⁺-independent basal activity of RyR2, suggesting that E3848A may also affect the stability of the closed state of RyR2. Mutations in the Ca²⁺ activation site also abolished the effect of ATP/caffeine on the Ca²⁺-independent basal activity, suggesting that the Ca²⁺ activation site is also a critical determinant of ATP/caffeine action. Mutating residues with negatively charged or oxygen-containing side chains near the Ca²⁺ activation site significantly altered Ca²⁺ and caffeine activation and reduced Mg²⁺ inhibition. Furthermore, disease-associated RyR2 mutations within the central domain significantly enhanced Ca²⁺ and caffeine activation and reduced Mg²⁺ inhibition. Our data demonstrate that the central domain plays an important role in channel activation, channel regulation, and closed state stability.

Elementary calcium signaling in arterial smooth muscle

Vascular smooth muscle cells (VSMCs) of small peripheral arteries contribute to blood pressure control by adapting their contractile state. These adaptations depend on the VSMC cytosolic Ca²⁺ concentration, regulated by complex local elementary Ca²⁺ signaling pathways. Ca²⁺ sparks represent local, transient, rapid calcium release events from a cluster of ryanodine receptors (RyRs) in the sarcoplasmic reticulum. In arterial SMCs, Ca²⁺ sparks activate nearby calcium-dependent potassium channels, cause membrane hyperpolarization and thus decrease the global intracellular [Ca²⁺] to oppose vasoconstriction. Arterial SMC Ca_v1.2 L-type channels regulate intracellular calcium stores content, which in turn modulates calcium efflux through RyRs. Ca_v3.2 T-type channels contribute to a minor extend to Ca²⁺ spark generation in certain types of arteries. Their localization within cell membrane caveolae is essential. We summarize present data on local elementary calcium signaling (Ca²⁺ sparks) in arterial SMCs with focus on RyR isoforms, large-conductance calcium-dependent potassium (BK_Ca) channels, and cell membrane-bound calcium channels (Ca_v1.2 and Ca_v3.2), particularly in caveolar microdomains.

Junctional sarcoplasmic reticulum motility in adult mouse ventricular myocytes

Excitation-contraction (EC) coupling is the coordinated process by which an action potential triggers cardiac myocyte contraction. EC coupling is initiated in dyads where the junctional sarcoplasmic reticulum (jSR) is in tight proximity to the sarcolemma of cardiac myocytes. Existing models of EC coupling critically depend on dyad stability to ensure the fidelity and strength of EC coupling, where even small variations in ryanodine receptor channel and voltage-gated calcium channel-α 1.2 subunit separation dramatically alter EC coupling. However, dyadic motility has never been studied. Here, we developed a novel strategy to track specific jSR units in dissociated adult ventricular myocytes using photoactivatable fluorescent proteins. We found that the jSR is not static. Instead, we observed dynamic formation and dissolution of multiple dyadic junctions regulated by the microtubule-associated molecular motors kinesin-1 and dynein. Our data support a model where reproducibility of EC coupling results from the activation of a temporally averaged number of SR Ca²⁺ release units forming and dissolving SR-sarcolemmal junctions. These findings challenge the long-held view that the jSR is an immobile structure and provide insights into the mechanisms underlying its motility.

Regulation of Mitochondrial ATP Production: Ca2+ Signaling and Quality Control

Cardiac ATP production primarily depends on oxidative phosphorylation in mitochondria and is dynamically regulated by Ca²⁺ levels in the mitochondrial matrix as well as by cytosolic ADP. We discuss mitochondrial Ca²⁺ signaling and its dysfunction which has recently been linked to cardiac pathologies including arrhythmia and heart failure. Similar dysfunction in other excitable and long-lived cells including neurons is associated with neurodegenerative diseases such as Alzheimer's disease (AD), amyotrophic lateral sclerosis (ALS), and Parkinson's disease (PD). Central to this new understanding is crucial Ca²⁺ regulation of both mitochondrial quality control and ATP production. Mitochondria-associated membrane (MAM) signaling from the sarcoplasmic reticulum (SR) and the endoplasmic reticulum (ER) to mitochondria is discussed. We propose future research directions that emphasize a need to define quantitatively the physiological roles of MAMs, as well as mitochondrial quality control and ATP production.

Triclosan Affects Ca2+ Regulatory Module and Musculature Development in Skeletal Myocyte during Early Life Stages of Zebrafish (Danio rerio)

Advanced technologies for toxicity tests are designed to identify biomarkers with superior predictive power or end points of the complex web of biological pathways. However, the data obtained need to be fully characterized for dose-response, physiological systems, and relevance to a system or (sub) population before biological interpretation and decision making. In this study, the toxicity of triclosan (TCS) on zebrafish was selected as a case study to correlate the observed morphological effects with existing data and identify the critical events by receptor activity sensitivity analysis. Triclosan exhibited weak acute toxicity against zebrafish and significantly affected the development of trunk muscles at 0.52, 1.04, and 1.73 μM. Through receptor-mediated screening, we found that the adverse effects of TCS induce Ryanodine receptor 1 (RyR1) activity and distort Ca²⁺ signaling. The trunk skeletal muscle abnormalities occurred only when the dihydropyridine receptor (DHPR) was blocked, demonstrating that TCS mainly influenced the Ca²⁺ regulatory module associated with signaling between DHPRs and RyR1; DHPRs mainly regulated the orthograde and retrograde signaling in skeletal muscles. This unexpected result could integrate the mode of action of TCS and provide insight for high-throughput screening and toxicity prediction using zebrafish.

Dynamic and Irregular Distribution of RyR2 Clusters in the Periphery of Live Ventricular Myocytes

Cardiac ryanodine receptors (RyR2s) are Ca²⁺ release channels clustering in the sarcoplasmic reticulum membrane. These clusters are believed to be the elementary units of Ca²⁺ release. The distribution of these Ca²⁺ release units plays a critical role in determining the spatio-temporal profile and stability of sarcoplasmic reticulum Ca²⁺ release. RyR2 clusters located in the interior of cardiomyocytes are arranged in highly ordered arrays. However, little is known about the distribution and function of RyR2 clusters in the periphery of cardiomyocytes. Here, we used a knock-in mouse model expressing a green fluorescence protein (GFP)-tagged RyR2 to localize RyR2 clusters in live ventricular myocytes by virtue of their GFP fluorescence. Confocal imaging and total internal reflection fluorescence microscopy was employed to determine and compare the distribution of GFP-RyR2 in the interior and periphery of isolated live ventricular myocytes and in intact hearts. We found tightly ordered arrays of GFP-RyR2 clusters in the interior, as previously described. In contrast, irregular distribution of GFP-RyR2 clusters was observed in the periphery. Time-lapse total internal reflection fluorescence imaging revealed dynamic movements of GFP-RyR2 clusters in the periphery, which were affected by external Ca²⁺ and RyR2 activator (caffeine) and inhibitor (tetracaine), but little detectable movement of GFP-RyR2 clusters in the interior. Furthermore, simultaneous Ca²⁺- and GFP-imaging demonstrated that peripheral RyR2 clusters with an irregular distribution pattern are functional with a Ca²⁺ release profile similar to that in the interior. These results indicate that the distribution of RyR2 clusters in the periphery of live ventricular myocytes is irregular and dynamic, which is different from that of RyR2 clusters in the interior.

3D dSTORM imaging reveals novel detail of ryanodine receptor localization in rat cardiac myocytes

KEY POINTS

Using 3D direct stochastic optical reconstruction microscopy (dSTORM), we developed novel approaches to quantitatively describe the nanoscale, 3D organization of ryanodine receptors (RyRs) in cardiomyocytes. Complex arrangements of RyR clusters were observed in 3D space, both at the cell surface and within the cell interior, with allocation to dyadic and non-dyadic pools. 3D imaging importantly allowed discernment of clusters overlapping in the z-axis, for which detection was obscured by conventional 2D imaging techniques. Thus, RyR clusters were found to be significantly smaller than previous 2D estimates. Ca2+ release units (CRUs), i.e. functional groupings of neighbouring RyR clusters, were similarly observed to be smaller than earlier reports. Internal CRUs contained more RyRs in more clusters than CRUs on the cell surface, and yielded longer duration Ca2+ sparks.

ABSTRACT

Cardiomyocyte contraction is dependent on Ca2+ release from ryanodine receptors (RyRs). However, the precise localization of RyRs remains unknown, due to shortcomings of imaging techniques which are diffraction limited or restricted to 2D. We aimed to determine the 3D nanoscale organization of RyRs in rat cardiomyocytes by employing direct stochastic optical reconstruction microscopy (dSTORM) with phase ramp technology. Initial observations at the cell surface showed an undulating organization of RyR clusters, resulting in their frequent overlap in the z-axis and obscured detection by 2D techniques. Non-overlapping clusters were imaged to create a calibration curve for estimating RyR number based on recorded fluorescence blinks. Employing this method at the cell surface and interior revealed smaller RyR clusters than 2D estimates, as erroneous merging of axially aligned RyRs was circumvented. Functional groupings of RyR clusters (Ca2+ release units, CRUs), contained an average of 18 and 23 RyRs at the surface and interior, respectively, although half of all CRUs contained only a single 'rogue' RyR. Internal CRUs were more tightly packed along z-lines than surface CRUs, contained larger and more numerous RyR clusters, and constituted ∼75% of the roughly 1 million RyRs present in an average cardiomyocyte. This complex internal 3D geometry was underscored by correlative imaging of RyRs and t-tubules, which enabled quantification of dyadic and non-dyadic RyR populations. Mirroring differences in CRU size and complexity, Ca2+ sparks originating from internal CRUs were of longer duration than those at the surface. These data provide novel, nanoscale insight into RyR organization and function across cardiomyocytes.

Ryanodine receptor dispersion disrupts Ca2+ release in failing cardiac myocytes

Reduced cardiac contractility during heart failure (HF) is linked to impaired Ca²⁺ release from Ryanodine Receptors (RyRs). We investigated whether this deficit can be traced to nanoscale RyR reorganization. Using super-resolution imaging, we observed dispersion of RyR clusters in cardiomyocytes from post-infarction HF rats, resulting in more numerous, smaller clusters. Functional groupings of RyR clusters which produce Ca²⁺ sparks (Ca²⁺ release units, CRUs) also became less solid. An increased fraction of small CRUs in HF was linked to augmented ‘silent’ Ca²⁺ leak, not visible as sparks. Larger multi-cluster CRUs common in HF also exhibited low fidelity spark generation. When successfully triggered, sparks in failing cells displayed slow kinetics as Ca²⁺ spread across dispersed CRUs. During the action potential, these slow sparks protracted and desynchronized the overall Ca²⁺ transient. Thus, nanoscale RyR reorganization during HF augments Ca²⁺ leak and slows Ca²⁺ release kinetics, leading to weakened contraction in this disease.

The muscle cells of the heart coordinate how they contract and relax in order to produce the heartbeat. During heart failure, these cells become less able to contract. As a result the heart becomes inefficient, pumping less blood around the body.

For the cardiac muscle cells to contract, the levels of calcium ions in the cells needs to rapidly increase. In failing hearts, these increases in calcium ion levels are smaller, slower and less well coordinated. It was not known what causes these changes, making it difficult to treat heart failure.

Calcium ions are released in cardiac muscle cells through protein channels called ryanodine receptors. These receptors form clusters that allow them to synchronize when they open and close. Could the reorganization of ryanodine receptors account for the problems seen in failing hearts? To investigate, Kolstad et al. examined rat hearts using a technique called super-resolution microscopy. This showed that the clusters of ryanodine receptors break apart during heart failure to form smaller clusters. Further experiments showed that calcium ions ‘leak’ from these smaller clusters, reducing the amount of calcium that can be released into cardiac muscle cells during each heartbeat. Released calcium also spreads between the dispersed clusters, resulting in a slower rise of the calcium levels in the cells. Both changes contribute to weakened contractions of cells in failing hearts.

Therefore, heart failure can be traced back to very small rearrangements of the ryanodine receptors. This understanding will help researchers as they investigate new ways to treat heart failure.

Role of defective Ca2+ signaling in skeletal muscle weakness: Pharmacological implications

The misbehaving attitude of Ca²⁺ signaling pathways could be the probable reason in many muscular disorders such as myopathies, systemic disorders like hypoxia, sepsis, cachexia, sarcopenia, heart failure, and dystrophy. The present review throws light upon the calcium flux regulating signaling channels like ryanodine receptor complex (RyR1), SERCA (Sarco-endoplasmic Reticulum Calcium ATPase), DHPR (Dihydropyridine Receptor) or Cav1.1 and Na+/Ca²⁺ exchange pump in detail and how remodelling of these channels contribute towards disturbed calcium homeostasis. Understanding these pathways will further provide an insight for establishing new therapeutic approaches for the prevention and treatment of muscle atrophy under stress conditions, targeting calcium ion channels and associated regulatory proteins.

The intricacies of atrial calcium cycling during excitation-contraction coupling

Blatter discusses the initiation and spread of Ca release, Ca store depletion, and release termination in atrial myocytes.

MarkoLAB: A simulator to study ionic channel's stochastic behavior

Mathematical models of the cardiac cell have started to include markovian representations of the ionic channels instead of the traditional Hodgkin & Huxley formulations. There are many reasons for this: Markov models are not restricted to the idea of independent gates defining the channel, they allow more complex description with specific transitions between open, closed or inactivated states, and more importantly those states can be closely related to the underlying channel structure and conformational changes.

METHODS

We used the LabVIEW^® and MATLAB^® programs to implement the simulator MarkoLAB that allow a dynamical 3D representation of the markovian model of the channel. The Monte Carlo simulation was used to implement the stochastic transitions among states. The user can specify the voltage protocol by setting the holding potential, the step-to voltage and the duration of the stimuli.

RESULTS

The most studied feature of a channel is the current flowing through it. This happens when the channel stays in the open state, but most of the time, as revealed by the low open probability values, the channel remains on the inactive or closed states. By focusing only when the channel enters or leaves the open state we are missing most of its activity. MarkoLAB proved to be quite useful to visualize the whole behavior of the channel and not only when the channel produces a current. Such dynamic representation provides more complete information about channel kinetics and will be a powerful tool to demonstrate the effect of gene mutations or drugs on the channel function.

CONCLUSIONS

MarkoLAB provides an original way of visualizing the stochastic behavior of a channel. It clarifies concepts, such as recovery from inactivation, calcium- versus voltage-dependent inactivation, and tail currents. It is not restricted to ionic channels only but it can be extended to other transporters, such as exchangers and pumps. This program is intended as a didactical tool to illustrate the dynamical behavior of a channel. It has been implemented in two platforms MATLAB^® and LabVIEW^® to enhance the target users of this new didactical tool. The computational cost of implementing a stochastic simulation is within the range of a personal computer performance; making MarkoLAB suitable to be run during a lecture or presentation.

Transverse tubule remodelling: a cellular pathology driven by both sides of the plasmalemma?

Transverse (t)-tubules are invaginations of the plasma membrane that form a complex network of ducts, 200-400 nm in diameter depending on the animal species, that penetrates deep within the cardiac myocyte, where they facilitate a fast and synchronous contraction across the entire cell volume. There is now a large body of evidence in animal models and humans demonstrating that pathological distortion of the t-tubule structure has a causative role in the loss of myocyte contractility that underpins many forms of heart failure. Investigations into the molecular mechanisms of pathological t-tubule remodelling to date have focused on proteins residing in the intracellular aspect of t-tubule membrane that form linkages between the membrane and myocyte cytoskeleton. In this review, we shed light on the mechanisms of t-tubule remodelling which are not limited to the intracellular side. Our recent data have demonstrated that collagen is an integral part of the t-tubule network and that it increases within the tubules in heart failure, suggesting that a fibrotic mechanism could drive cardiac junctional remodelling. We examine the evidence that the linkages between the extracellular matrix, t-tubule membrane and cellular cytoskeleton should be considered as a whole when investigating the mechanisms of t-tubule pathology in the failing heart.

Studying dyadic structure-function relationships: a review of current modeling approaches and new insights into Ca2+ (mis)handling

Excitation–contraction coupling in cardiac myocytes requires calcium influx through L-type calcium channels in the sarcolemma, which gates calcium release through sarcoplasmic reticulum ryanodine receptors in a process known as calcium-induced calcium release, producing a myoplasmic calcium transient and enabling cardiomyocyte contraction. The spatio-temporal dynamics of calcium release, buffering, and reuptake into the sarcoplasmic reticulum play a central role in excitation–contraction coupling in both normal and diseased cardiac myocytes. However, further quantitative understanding of these cells’ calcium machinery and the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease requires accurate knowledge of cardiac ultrastructure, protein distribution and subcellular function. As current imaging techniques are limited in spatial resolution, limiting insight into changes in calcium handling, computational models of excitation–contraction coupling have been increasingly employed to probe these structure–function relationships. This review will focus on the development of structural models of cardiac calcium dynamics at the subcellular level, orienting the reader broadly towards the development of models of subcellular calcium handling in cardiomyocytes. Specific focus will be given to progress in recent years in terms of multi-scale modeling employing resolved spatial models of subcellular calcium machinery. A review of the state-of-the-art will be followed by a review of emergent insights into calcium-dependent etiologies in heart disease and, finally, we will offer a perspective on future directions for related computational modeling and simulation efforts.

A novel mechanism of tandem activation of ryanodine receptors by cytosolic and SR luminal Ca2+ during excitation-contraction coupling in atrial myocytes

KEY POINTS

In atrial myocytes excitation-contraction coupling is strikingly different from ventricle because atrial myocytes lack a transverse tubule membrane system: Ca²⁺ release starts in the cell periphery and propagates towards the cell centre by Ca²⁺ -induced Ca²⁺ release from the sarcoplasmic reticulum (SR) Ca²⁺ store. The cytosolic Ca²⁺ sensitivity of the ryanodine receptor (RyRs) Ca²⁺ release channel is low and it is unclear how Ca²⁺ release can be activated in the interior of atrial cells. Simultaneous confocal imaging of cytosolic and intra-SR calcium revealed a transient elevation of store Ca²⁺ that we termed 'Ca²⁺ sensitization signal'. We propose a novel paradigm of atrial ECC that is based on tandem activation of the RyRs by cytosolic and luminal Ca²⁺ through a 'fire-diffuse-uptake-fire' (or FDUF) mechanism: Ca²⁺ uptake by SR Ca²⁺ pumps at the propagation front elevates Ca²⁺ inside the SR locally, leading to luminal RyR sensitization and lowering of the cytosolic Ca²⁺ activation threshold.

ABSTRACT

In atrial myocytes Ca²⁺ release during excitation-contraction coupling (ECC) is strikingly different from ventricular myocytes. In many species atrial myocytes lack a transverse tubule system, dividing the sarcoplasmic reticulum (SR) Ca²⁺ store into the peripheral subsarcolemmnal junctional (j-SR) and the much more abundant central non-junctional (nj-SR) SR. Action potential (AP)-induced Ca²⁺ entry activates Ca²⁺ -induced Ca²⁺ release (CICR) from j-SR ryanodine receptor (RyR) Ca²⁺ release channels. Peripheral elevation of [Ca²⁺ ]_i initiates CICR from nj-SR and sustains propagation of CICR to the cell centre. Simultaneous confocal measurements of cytosolic ([Ca²⁺ ]_i ; with the fluorescent Ca²⁺ indicator rhod-2) and intra-SR ([Ca²⁺ ]_SR ; fluo-5N) Ca²⁺ in rabbit atrial myocytes revealed that Ca²⁺ release from j-SR resulted in a cytosolic Ca²⁺ transient of higher amplitude compared to release from nj-SR; however, the degree of depletion of j-SR [Ca²⁺ ]_SR was smaller than nj-SR [Ca²⁺ ]_SR . Similarly, Ca²⁺ signals from individual release sites of the j-SR showed a larger cytosolic amplitude (Ca²⁺ sparks) but smaller depletion (Ca²⁺ blinks) than release from nj-SR. During AP-induced Ca²⁺ release the rise of [Ca²⁺ ]_i detected at individual release sites of the nj-SR preceded the depletion of [Ca²⁺ ]_SR , and during this latency period a transient elevation of [Ca²⁺ ]_SR occurred. We propose that Ca²⁺ release from nj-SR is activated by cytosolic and luminal Ca²⁺ (tandem RyR activation) via a novel 'fire-diffuse-uptake-fire' (FDUF) mechanism. This novel paradigm of atrial ECC predicts that Ca²⁺ uptake by sarco-endoplasmic reticulum Ca²⁺ -ATPase (SERCA) at the propagation front elevates local [Ca²⁺ ]_SR , leading to luminal RyR sensitization and lowering of the activation threshold for cytosolic CICR.

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

The cardiac Ca(2+) 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(2+) sparks showed that Ca(2+) sparks originated exclusively from RyR2 clusters. Ca(2+) sparks from RyR2 clusters induced no detectable changes in mitochondrial Ca(2+) level. These results reveal, for the first time, the distribution of RyR2 clusters and its functional correlation in living ventricular myocytes.

Multiple modes of ryanodine receptor 2 inhibition by flecainide

Catecholaminergic polymorphic ventricular tachycardia (CPVT) causes sudden cardiac death due to mutations in cardiac ryanodine receptors (RyR2), calsequestrin, or calmodulin. Flecainide, a class I antiarrhythmic drug, inhibits Na(+) and RyR2 channels and prevents CPVT. The purpose of this study is to identify inhibitory mechanisms of flecainide on RyR2. RyR2 were isolated from sheep heart, incorporated into lipid bilayers, and investigated by single-channel recording under various activating conditions, including the presence of cytoplasmic ATP (2 mM) and a range of cytoplasmic [Ca(2+)], [Mg(2+)], pH, and [caffeine]. Flecainide applied to either the cytoplasmic or luminal sides of the membrane inhibited RyR2 by two distinct modes: 1) a fast block consisting of brief substate and closed events with a mean duration of ∼1 ms, and 2) a slow block consisting of closed events with a mean duration of ∼1 second. Both inhibition modes were alleviated by increasing cytoplasmic pH from 7.4 to 9.5 but were unaffected by luminal pH. The slow block was potentiated in RyR2 channels that had relatively low open probability, whereas the fast block was unaffected by RyR2 activation. These results show that these two modes are independent mechanisms for RyR2 inhibition, both having a cytoplasmic site of action. The slow mode is a closed-channel block, whereas the fast mode blocks RyR2 in the open state. At diastolic cytoplasmic [Ca(2+)] (100 nM), flecainide possesses an additional inhibitory mechanism that reduces RyR2 burst duration. Hence, multiple modes of action underlie RyR2 inhibition by flecainide.

Genetic interactions found between calcium channel genes modulate amyloid load measured by positron emission tomography

Late-onset Alzheimer's disease (LOAD) is known to have a complex, oligogenic etiology, with considerable genetic heterogeneity. We investigated the influence of genetic interactions between genes in the Alzheimer's disease (AD) pathway on amyloid-beta (Aβ) deposition as measured by PiB or AV-45 ligand positron emission tomography (PET) to aid in understanding LOAD's genetic etiology. Subsets of the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohorts were used for discovery and for two independent validation analyses. A significant interaction between RYR3 and CACNA1C was confirmed in all three of the independent ADNI datasets. Both genes encode calcium channels expressed in the brain. The results shown here support previous animal studies implicating interactions between these calcium channels in amyloidogenesis and suggest that the pathological cascade of this disease may be modified by interactions in the amyloid-calcium axis. Future work focusing on the mechanisms of such relationships may inform targets for clinical intervention.

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(2+) spark evolution remain unclear. Cardiac ryanodine receptor (RyR) gating in rats and sheep was recorded at physiological Ca(2+), Mg(2+), and ATP levels and incorporated into a 3D model of the cardiac dyad, which reproduced the time course of Ca(2+) sparks, Ca(2+) blinks, and Ca(2+) spark restitution. The termination of CICR by induction decay in the model principally arose from the steep Ca(2+) dependence of RyR closed time, with the measured sarcoplasmic reticulum (SR) lumen Ca(2+) 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(2+), as well as the time course of local Ca(2+) gradients within the junctional space. Reducing the dimensions of the dyad junction reduced Ca(2+) spark amplitude by reducing the strength of regenerative feedback within CICR. A refractory period for Ca(2+) spark initiation and subsequent Ca(2+) 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(2+) in the junctional SR. The physical organization of RyRs within the junctional space had minimal effects on Ca(2+) spark amplitude when more than nine RyRs were present. Spark amplitude had a nonlinear dependence on RyR single-channel Ca(2+) 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(2+) sensitivities, Ca(2+) 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.

Ca²⁺ waves in the heart

Ca(2+) waves were probably first observed in the early 1940s. Since then Ca(2+) waves have captured the attention of an eclectic mixture of mathematicians, neuroscientists, muscle physiologists, developmental biologists, and clinical cardiologists. This review discusses the current state of mathematical models of Ca(2+) waves, the normal physiological functions Ca(2+) waves might serve in cardiac cells, as well as how the spatial arrangement of Ca(2+) release channels shape Ca(2+) waves, and we introduce the idea of Ca(2+) phase waves that might provide a useful framework for understanding triggered arrhythmias.

Ca(2+) release events in cardiac myocytes up close: insights from fast confocal imaging

The spatio-temporal properties of Ca²⁺ transients during excitation-contraction coupling and elementary Ca²⁺ release events (Ca²⁺ sparks) were studied in atrial and ventricular myocytes with ultra-fast confocal microscopy using a Zeiss LSM 5 LIVE system that allows sampling rates of up to 60 kHz. Ca²⁺ sparks which originated from subsarcolemmal junctional sarcoplasmic reticulum (j-SR) release sites in atrial myocytes were anisotropic and elongated in the longitudinal direction of the cell. Ca²⁺ sparks in atrial cells originating from non-junctional SR and in ventricular myocytes were symmetrical. Ca²⁺ spark recording in line scan mode at 40,000 lines/s uncovered step-like increases of [Ca²⁺]_i. 2-D imaging of Ca²⁺ transients revealed an asynchronous activation of release sites and allowed the sequential recording of Ca²⁺ entry through surface membrane Ca²⁺ channels and subsequent activation of Ca²⁺-induced Ca²⁺ release. With a latency of 2.5 ms after application of an electrical stimulus, Ca²⁺ entry could be detected that was followed by SR Ca²⁺ release after an additional 3 ms delay. Maximum Ca²⁺ release was observed 4 ms after the beginning of release. The timing of Ca²⁺ entry and release was confirmed by simultaneous [Ca²⁺]_i and membrane current measurements using the whole cell voltage-clamp technique. In atrial cells activation of discrete individual release sites of the j-SR led to spatially restricted Ca²⁺ release events that fused into a peripheral ring of elevated [Ca²⁺]_i that subsequently propagated in a wave-like fashion towards the center of the cell. In ventricular myocytes asynchronous Ca²⁺ release signals from discrete sites with no preferential subcellular location preceded the whole-cell Ca²⁺ transient. In summary, ultra-fast confocal imaging allows investigation of Ca²⁺ signals with a time resolution similar to patch clamp technique, however in a less invasive fashion.

Facilitation of cytosolic calcium wave propagation by local calcium uptake into the sarcoplasmic reticulum in cardiac myocytes

The widely accepted paradigm for cytosolic Ca(2+) wave propagation postulates a 'fire-diffuse-fire' mechanism where local Ca(2+)-induced Ca(2+) release (CICR) from the sarcoplasmic reticulum (SR) via ryanodine receptor (RyR) Ca(2+) release channels diffuses towards and activates neighbouring release sites, resulting in a propagating Ca(2+) wave. A recent challenge to this paradigm proposed the requirement for an intra-SR 'sensitization' Ca(2+) wave that precedes the cytosolic Ca(2+) wave and primes RyRs from the luminal side to CICR. Here, we tested this hypothesis experimentally with direct simultaneous measurements of cytosolic ([Ca(2+)](i); rhod-2) and intra-SR ([Ca(2+)](SR); fluo-5N) calcium signals during wave propagation in rabbit ventricular myocytes, using high resolution fluorescence confocal imaging. The increase in [Ca(2+)](i) at the wave front preceded depletion of the SR at each point along the calcium wave front, while during this latency period a transient increase of [Ca(2+)](SR) was observed. This transient elevation of [Ca(2+)](SR) could be identified at individual release junctions and depended on the activity of the sarco-endoplasmic reticulum Ca(2+)-ATPase (SERCA). Increased SERCA activity (β-adrenergic stimulation with 1 μM isoproterenol (isoprenaline)) decreased the latency period and increased the amplitude of the transient elevation of [Ca(2+)](SR), whereas inhibition of SERCA (3 μM cyclopiazonic acid) had the opposite effect. In conclusion, the data provide experimental evidence that local Ca(2+) uptake by SERCA into the SR facilitates the propagation of cytosolic Ca(2+) waves via luminal sensitization of the RyR, and supports a novel paradigm of a 'fire-diffuse-uptake-fire' mechanism for Ca(2+) wave propagation in cardiac myocytes.

Super-resolution imaging of EC coupling protein distribution in the heart

The cardiac ryanodine receptor (RyR) plays a central role in the control of contractile function of the heart. In cardiac ventricular myocytes RyRs and associated Ca(2+) handling proteins, including membrane Ca(2+) channels, Ca(2+) pumps and other sarcolemmal and sarcoplasmic reticulum proteins interact to set the time course and amplitude of the electrically triggered cytosolic Ca(2+) transient. It has become increasingly clear that protein distribution and clustering on the nanometer scale is critical in determining the interaction of these proteins and the resulting properties of cardiac Ca(2+) handling. Such intricate near-molecular scale detail cannot be visualized with conventional fluorescence microscopy techniques (e.g. confocal microscopy) but it has recently become accessible with optical super-resolution techniques. These techniques retain the advantages of fluorescent marker technology, i.e. high specificity and excellent contrast, but have a spatial resolution approaching 10nm, i.e. objects not much further apart than 10nm can be distinguished, previously only attainable with electron microscopy. We review the use of these novel imaging techniques for the study of protein distribution in cardiac ventricular myocytes and discuss technical considerations as well as recent findings using super-resolution imaging. An emphasis is on single molecule localization based super-resolution approaches and their use to reveal the complexity of RyR cluster morphology, placement and relationship to other excitation-contraction coupling proteins. Super-resolution imaging approaches have already demonstrated their utility for the study of cardiac structure-function relationships and we anticipate that their use will rapidly increase and help improve our understanding of cardiac Ca(2+) regulation.

The proximal C-terminus of α1C subunits is necessary for junctional membrane targeting of cardiac L-type calcium channels

In cardiac myocytes, LTCCs (L-type calcium channels) form a functional signalling complex with ryanodine receptors at the JM (junctional membrane). Although the specific localization of LTCCs to the JM is critical for excitation–contraction coupling, their targeting mechanism is unclear. Transient transfection of GFP (green fluorescent protein)–α1S or GFP–α1C, but not P/Q-type calcium channel α1A, in dysgenic (α1S-null) GLT myotubes results in correct targeting of these LTCCs to the JMs and restoration of action-potential-induced Ca2+ transients. To identify the sequences of α1C responsible for JM targeting, we generated a range of α1C–α1A chimaeras, deletion mutants and alanine substitution mutants and studied their targeting properties in GLT myotubes. The results revealed that amino acids L1681QAGLRTL1688 and P1693EIRRAIS1700, predicted to form two adjacent α-helices in the proximal C-terminus, are necessary for the JM targeting of α1C. The efficiency of restoration of action-potential-induced Ca2+ transients in GLT myotubes was significantly decreased by mutations in the targeting motif. JM targeting was not disrupted by the distal C-terminus of α1C which binds to the second α-helix. Therefore we have identified a new structural motif in the C-terminus of α1C that mediates the targeting of cardiac LTCCs to JMs independently of the interaction between proximal and distal C-termini of α1C.

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.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

Nanoscale organization of junctophilin-2 and ryanodine receptors within peripheral couplings of rat ventricular cardiomyocytes

The peripheral distributions of the cardiac ryanodine receptor (RyR) and a junctional protein, junctophilin-2 (JPH2), were examined using single fluorophore localization-based super-resolution microscopy in rat ventricular myocytes. JPH2 was strongly associated with RyR clusters. Estimates of the colocalizing fraction of JPH labeling with RyR was ~90% within 30 nm of RyR clusters. This is comparable to fractions estimated from confocal data (~87%). Similarly, most RyRs were associated with JPH2 labeling in super-resolution images (~81% within 30 nm of JPH2 clusters). The shape of associated RyR clusters and JPH2 clusters were very similar, but not identical, suggesting that JPH2 is dispersed throughout RyR clusters and that the packing of JPH2 into junctions and the assembly of RyR clusters are tightly linked.

Simulation strategies for calcium microdomains and calcium-regulated calcium channels

In this article, we present an overview of simulation strategies in the context of subcellular domains where calcium-dependent signaling plays an important role. The presentation follows the spatial and temporal scales involved and represented by each algorithm. As an exemplary cell type, we will mainly cite work done on striated muscle cells, i.e. skeletal and cardiac muscle. For these cells, a wealth of ultrastructural, biophysical and electrophysiological data is at hand. Moreover, these cells also express ubiquitous signaling pathways as they are found in many other cell types and thus, the generalization of the methods and results presented here is straightforward.The models considered comprise the basic calcium signaling machinery as found in most excitable cell types including Ca(2+) ions, diffusible and stationary buffer systems, and calcium regulated calcium release channels. Simulation strategies can be differentiated in stochastic and deterministic algorithms. Historically, deterministic approaches based on the macroscopic reaction rate equations were the first models considered. As experimental methods elucidated highly localized Ca(2+) signaling events occurring in femtoliter volumes, stochastic methods were increasingly considered. However, detailed simulations of single molecule trajectories are rarely performed as the computational cost implied is too large. On the mesoscopic level, Gillespie's algorithm is extensively used in the systems biology community and with increasing frequency also in models of microdomain calcium signaling. To increase computational speed, fast approximations were derived from Gillespie's exact algorithm, most notably the chemical Langevin equation and the τ-leap algorithm. Finally, in order to integrate deterministic and stochastic effects in multiscale simulations, hybrid algorithms are increasingly used. These include stochastic models of ion channels combined with deterministic descriptions of the calcium buffering and diffusion system on the one hand, and algorithms that switch between deterministic and stochastic simulation steps in a context-dependent manner on the other. The basic assumptions of the listed methods as well as implementation schemes are given in the text. We conclude with a perspective on possible future developments of the field.

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).

Proton-sensing Ca2+ binding domains regulate the cardiac Na+/Ca2+ exchanger

The cardiac Na(+)/Ca(2+) exchanger (NCX) regulates cellular [Ca(2+)](i) and plays a central role in health and disease, but its molecular regulation is poorly understood. Here we report on how protons affect this electrogenic transporter by modulating two critically important NCX C(2) regulatory domains, Ca(2+) binding domain-1 (CBD1) and CBD2. The NCX transport rate in intact cardiac ventricular myocytes was measured as a membrane current, I(NCX), whereas [H(+)](i) was varied using an ammonium chloride "rebound" method at constant extracellular pH 7.4. At pH(i) = 7.2 and [Ca(2+)](i) < 120 nM, I(NCX) was less than 4% that of its maximally Ca(2+)-activated value. I(NCX) increases steeply at [Ca(2+)](i) between 130-150 nM with a Hill coefficient (n(H)) of 8.0 ± 0.7 and K(0.5) = 310 ± 5 nM. At pH(i) = 6.87, the threshold of Ca(2+)-dependent activation of I(NCX) was shifted to much higher [Ca(2+)](i) (600-700 nM), and the relationship was similarly steep (n(H) = 8.0±0.8) with K(0.5) = 1042 ± 15 nM. The V(max) of Ca(2+)-dependent activation of I(NCX) was not significantly altered by low pH(i). The Ca(2+) affinities for CBD1 (0.39 ± 0.06 μM) and CBD2 (K(d) = 18.4 ± 6 μM) were exquisitely sensitive to [H(+)], decreasing 1.3-2.3-fold as pH(i) decreased from 7.2 to 6.9. This work reveals for the first time that NCX can be switched off by physiologically relevant intracellular acidification and that this depends on the competitive binding of protons to its C(2) regulatory domains CBD1 and CBD2.

Na+ currents are required for efficient excitation-contraction coupling in rabbit ventricular myocytes: a possible contribution of neuronal Na+ channels

Ca2+ transients were activated in rabbit ventricular cells by a sequence of action potential shaped voltage clamps. After activating a series of control transients, Na+ currents (INa) were inactivated with a ramp from -80 to -40 mV (1.5 s) prior to the action potential clamp. The transients were detected with the calcium indicator Fluo-4 and an epifluorescence system. With zero Na+ in the pipette INa inactivation produced a decline in the SR Ca2+ release flux (measured as the maximum rate of rise of the transient) of 27 ± 4% (n = 9, P < 0.001) and a peak amplitude reduction of 10 ± 3% (n = 9, P < 0.05). With 5 mm Na+ in the pipette the reduction in release flux was greater (34 ± 4%, n = 4, P < 0.05). The ramp effectively inactivates INa without changing ICa, and there was no significant change in the transmembrane Ca2+ flux after the inactivation of INa. We next evoked action potentials under current clamp. TTX at 100 nm, which selectively blocks neuronal isoforms of Na+ channels, produced a decline in SR Ca2+ release flux of 35 ± 3% (n = 6, P < 0.001) and transient amplitude of 12 ± 2% (n = 6, P < 0.05). This effect was similar to the effect of INa inactivation on release flux. We conclude that a TTX-sensitive INa is essential for efficient triggering of SR Ca2+ release. We propose that neuronal Na+ channels residing within couplons activate sufficient reverse Na+-Ca2+ exchanger (NCX) to prime the junctional cleft with Ca2+. The results can be explained if non-linearities in excitation-contraction coupling mechanisms modify the coupling fidelity of ICa, which is known to be low at positive potentials.

Changes in the organization of excitation-contraction coupling structures in failing human heart

BACKGROUND

The cardiac myocyte t-tubular system ensures rapid, uniform cell activation and several experimental lines of evidence suggest changes in the t-tubular system and associated excitation-contraction coupling proteins may occur in heart failure.

METHODS AND RESULTS

The organization of t-tubules, L-type calcium channels (DHPRs), ryanodine receptors (RyRs) and contractile machinery were examined in fixed ventricular tissue samples from both normal and failing hearts (idiopathic (non-ischemic) dilated cardiomyopathy) using high resolution fluorescent imaging. Wheat germ agglutinin (WGA), Na-Ca exchanger, DHPR and caveolin-3 labels revealed a shift from a predominantly transverse orientation to oblique and axial directions in failing myocytes. In failure, dilation of peripheral t-tubules occurred and a change in the extent of protein glycosylation was evident. There was no change in the fractional area occupied by myofilaments (labeled with phalloidin) but there was a small reduction in the number of RyR clusters per unit area. The general relationship between DHPRs and RyR was not changed and RyR labeling overlapped with 51±3% of DHPR labeling in normal hearts. In longitudinal (but not transverse) sections there was an ∼30% reduction in the degree of colocalization between DHPRs and RyRs as measured by Pearson's correlation coefficient in failing hearts.

CONCLUSIONS

The results show that extensive remodelling of the t-tubular network and associated excitation-contraction coupling proteins occurs in failing human heart. These changes may contribute to abnormal calcium handling in heart failure. The general organization of the t-system and changes observed in failure samples have subtle differences to some animal models although the general direction of changes are generally similar.

Local Ca2+ releases enable rapid heart rates in developing cardiomyocytes

The ability to generate homogeneous intracellular Ca(2+) oscillations at high frequency is the basis of the rhythmic contractions of mammalian cardiac myocytes. While the specific mechanisms and structures enabling homogeneous high-frequency Ca(2+) signals in adult cardiomyocytes are well characterized, it is not known how these kind of Ca(2+) signals are produced in developing cardiomyocytes. Here we investigated the mechanisms reducing spatial and temporal heterogeneity of cytosolic Ca(2+) signals in mouse embryonic ventricular cardiomyocytes. We show that in developing cardiomyocytes the propagating Ca(2+) signals are amplified in cytosol by local Ca(2+) releases. Local releases are based on regular 3-D sarcoplasmic reticulum (SR) structures containing SR Ca(2+) uptake ATPases (SERCA) and Ca(2+) release channels (ryanodine receptors, RyRs) at regular intervals throughout the cytosol. By evoking [Ca(2+)](i)-induced Ca(2+) sparks, the local release sites promote a 3-fold increase in the cytosolic Ca(2+) propagation speed. We further demonstrate by mathematical modelling that without these local release sites the developing cardiomyocytes lose their ability to generate homogeneous global Ca(2+) signals at a sufficiently high frequency. The mechanism described here is robust and indispensable for normal mammalian cardiomyocyte function from the first heartbeats during the early embryonic phase till terminal differentiation after birth. These results suggest that local cytosolic Ca(2+) releases are indispensable for normal cardiomyocyte development and function of developing heart.

Optical single-channel resolution imaging of the ryanodine receptor distribution in rat cardiac myocytes

We have applied an optical super-resolution technique based on single-molecule localization to examine the peripheral distribution of a cardiac signaling protein, the ryanodine receptor (RyR), in rat ventricular myocytes. RyRs form clusters with a mean size of approximately 14 RyRs per cluster, which is almost an order of magnitude smaller than previously estimated. Clusters were typically not circular (as previously assumed) but elongated with an average aspect ratio of 1.9. Edge-to-edge distances between adjacent RyR clusters were often <50 nm, suggesting that peripheral RyR clusters may exhibit strong intercluster signaling. The wide variation of cluster size, which follows a near-exponential distribution, is compatible with a stochastic cluster assembly process. We suggest that calcium sparks may be the result of the concerted activation of several RyR clusters forming a functional "supercluster" whose gating is controlled by both cytosolic and sarcoplasmic reticulum luminal calcium levels.

Calcium-Induced Calcium Release in Skeletal Muscle

Calcium-induced calcium release (CICR) was first discovered in skeletal muscle. CICR is defined as Ca2+ release by the action of Ca2+ alone without the simultaneous action of other activating processes. CICR is biphasically dependent on Ca2+ concentration; is inhibited by Mg2+, procaine, and tetracaine; and is potentiated by ATP, other adenine compounds, and caffeine. With depolarization of the sarcoplasmic reticulum (SR), a potential change of the SR membrane in which the luminal side becomes more negative, CICR is activated for several seconds and is then inactivated. All three types of ryanodine receptors (RyRs) show CICR activity. At least one RyR, RyR1, also shows non-CICR Ca2+ release, such as that triggered by the t-tubule voltage sensor, by clofibric acid, and by SR depolarization. Maximum rates of CICR, at the optimal Ca2+ concentration in the presence of physiological levels of ATP and Mg2+ determined in skinned fibers and fragmented SR, are much lower than the rate of physiological Ca2+ release. The primary event of physiological Ca2+ release, the Ca2+ spark, is the simultaneous opening of multiple channels, the coordinating mechanism of which does not appear to be CICR because of the low probability of CICR opening under physiological conditions. The coordination may require Ca2+, but in that case, some other stimulus or stimuli must be provided simultaneously, which is not CICR by definition. Thus CICR does not appear to contribute significantly to physiological Ca2+ release. On the other hand, CICR appears to play a key role in caffeine contracture and malignant hyperthermia. The potentiation of voltage-activated Ca2+ release by caffeine, however, does not seem to occur through secondary CICR, although the site where caffeine potentiates voltage-activated Ca2+ release might be the same site where caffeine potentiates CICR.

Local control of excitation-contraction coupling in human embryonic stem cell-derived cardiomyocytes

We investigated the mechanisms of excitation-contraction (EC) coupling in human embryonic stem cell-derived cardiomyocytes (hESC-CMs) and fetal ventricular myocytes (hFVMs) using patch-clamp electrophysiology and confocal microscopy. We tested the hypothesis that Ca²⁺ influx via voltage-gated L-type Ca²⁺ channels activates Ca²⁺ release from the sarcoplasmic reticulum (SR) via a local control mechanism in hESC-CMs and hFVMs. Field-stimulated, whole-cell [Ca²⁺]_i transients in hESC-CMs required Ca²⁺ entry through L-type Ca²⁺ channels, as evidenced by the elimination of such transients by either removal of extracellular Ca²⁺ or treatment with diltiazem, an L-type channel inhibitor. Ca²⁺ release from the SR also contributes to the [Ca²⁺]_i transient in these cells, as evidenced by studies with drugs interfering with either SR Ca²⁺ release (i.e. ryanodine and caffeine) or reuptake (i.e. thapsigargin and cyclopiazonic acid). As in adult ventricular myocytes, membrane depolarization evoked large L-type Ca²⁺ currents (I_Ca) and corresponding whole-cell [Ca²⁺]_i transients in hESC-CMs and hFVMs, and the amplitude of both I_Ca and the [Ca²⁺]_i transients were finely graded by the magnitude of the depolarization. hESC-CMs exhibit a decreasing EC coupling gain with depolarization to more positive test potentials, “tail” [Ca²⁺]_i transients upon repolarization from extremely positive test potentials, and co-localized ryanodine and sarcolemmal L-type Ca²⁺ channels, all findings that are consistent with the local control hypothesis. Finally, we recorded Ca²⁺ sparks in hESC-CMs and hFVMs. Collectively, these data support a model in which tight, local control of SR Ca²⁺ release by the I_Ca during EC coupling develops early in human cardiomyocytes.

Towards computational modeling of excitation-contraction coupling in cardiac myocytes: reconstruction of structures and proteins from confocal imaging

Computational models of excitation-contraction (EC) coupling in myocytes are valuable tools for studying the signaling cascade that transduces transmembrane voltage into mechanical responses. A key component of these models is the appropriate description of structures involved in EC coupling, such as the sarcolemma and ion channels. This study aims at developing an approach for spatial reconstruction of these structures. We exemplified our approach by reconstructing clusters of ryanodine receptors (RyRs) together with the sarcolemma of rabbit ventricular myocytes. The reconstructions were based on dual labeling and three-dimensional (3D) confocal imaging of segments of fixed and permeabilized myocytes lying flat or on end. The imaging led to 3D stacks of cross-sections through myocytes. Methods of digital image processing were applied to deconvolve, filter and segment these stacks. Finally, we created point meshes representing RyR distributions together with volume and surface meshes of the sarcolemma. We suggest that these meshes are suitable for computational studies of structure-function relationships in EC coupling. We propose that this approach can be extended to reconstruct other structures and proteins involved in EC coupling.

Calcium sparks

The calcium ion (Ca(2+)) is the simplest and most versatile intracellular messenger known. The discovery of Ca(2+) sparks and a related family of elementary Ca(2+) signaling events has revealed fundamental principles of the Ca(2+) signaling system. A newly appreciated "digital" subsystem consisting of brief, high Ca(2+) concentration over short distances (nanometers to microns) comingles with an "analog" global Ca(2+) signaling subsystem. Over the past 15 years, much has been learned about the theoretical and practical aspects of spark formation and detection. The quest for the spark mechanisms [the activation, coordination, and termination of Ca(2+) release units (CRUs)] has met unexpected challenges, however, and raised vexing questions about CRU operation in situ. Ample evidence shows that Ca(2+) sparks catalyze many high-threshold Ca(2+) processes involved in cardiac and skeletal muscle excitation-contraction coupling, vascular tone regulation, membrane excitability, and neuronal secretion. Investigation of Ca(2+) sparks in diseases has also begun to provide novel insights into hypertension, cardiac arrhythmias, heart failure, and muscular dystrophy. An emerging view is that spatially and temporally patterned activation of the digital subsystem confers on intracellular Ca(2+) signaling an exquisite architecture in space, time, and intensity, which underpins signaling efficiency, stability, specificity, and diversity. These recent advances in "sparkology" thus promise to unify the simplicity and complexity of Ca(2+) signaling in biology.

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.

Electrogenic Na+/Ca2+-exchange of nerve and muscle cells

The plasma membrane Na(+)/Ca(2+)-exchanger is a bi-directional electrogenic (3Na(+):1Ca(2+)) and voltage-sensitive ion transport mechanism, which is mainly responsible for Ca(2+)-extrusion. The Na(+)-gradient, required for normal mode operation, is created by the Na(+)-pump, which is also electrogenic (3Na(+):2K(+)) and voltage-sensitive. The Na(+)/Ca(2+)-exchanger operational modes are very similar to those of the Na(+)-pump, except that the uncoupled flux (Na(+)-influx or -efflux?) is missing. The reversal potential of the exchanger is around -40 mV; therefore, during the upstroke of the AP it is probably transiently activated, leading to Ca(2+)-influx. The Na(+)/Ca(2+)-exchange is regulated by transported and non-transported external and internal cations, and shows ATP(i)-, pH- and temperature-dependence. The main problem in determining the role of Na(+)/Ca(2+)-exchange in excitation-secretion/contraction coupling is the lack of specific (mode-selective) blockers. During recent years, evidence has been accumulated for co-localisation of the Na(+)-pump, and the Na(+)/Ca(2+)-exchanger and their possible functional interaction in the "restricted" or "fuzzy space." In cardiac failure, the Na(+)-pump is down-regulated, while the exchanger is up-regulated. If the exchanger is working in normal mode (Ca(2+)-extrusion) during most of the cardiac cycle, upregulation of the exchanger may result in SR Ca(2+)-store depletion and further impairment in contractility. If so, a normal mode selective Na(+)/Ca(2+)-exchange inhibitor would be useful therapy for decompensation, and unlike CGs would not increase internal Na(+). In peripheral sympathetic nerves, pre-synaptic alpha(2)-receptors may regulate not only the VSCCs but possibly the reverse Na(+)/Ca(2+)-exchange as well.

Spatial and temporal Ca2+, Mg2+, and ATP2- dynamics in cardiac dyads during calcium release

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

Effect of changes in action potential spike configuration, junctional sarcoplasmic reticulum micro-architecture and altered t-tubule structure in human heart failure

Using a Monte-Carlo model of L-type Ca2+ channel (DHPR) gating, we have examined the effect of changes in the early time course of the action potential as seen in human heart failure on excitation contraction coupling. The time course of DHPR Ca2+ influx was coupled into a simple model of sarcoplasmic reticulum Ca2+ release. Our model shows that the loss of the initial spike in human heart failure should reduce the synchrony of Ca2+ spark production and lead to the appearance of late Ca2+ sparks and greater non-uniformity of intracellular Ca2+. Within the junctional space of the cardiac dyad, a small increase in the mean distance of a DHPR from a RyR results in a marked decrease in the ability of the DHPR-mediated increase in local [Ca2+] concentration to activate RyRs. This suggests that the efficiency of EC coupling may be reduced if changes in micro-architecture develop and such effects have been noted in experimental models of heart failure. High resolution imaging of t-tubules in tachycardia-induced heart failure show deranged t-tubule structure. While in normal human hearts t-tubules run mainly in a radial direction, t-tubules in the heart failure samples were oriented more toward the long axis of the cell. In addition, t-tubules may become dilated and bifurcated. Our data suggest that changes in the micro-architecture of the cell and membrane structures associated with excitation-contraction coupling, combined with changes in early action potential configuration can reduce the efficiency by which Ca2+ influx via DHPRs can activate SR calcium release and cardiac contraction. While the underlying cause of these effects is unclear, our data suggest that geometric factors can play an important role in the pathophysilogy of the human heart in failure.

Temperature and Ca2+dependence of [3H]ryanodine binding in the burbot (Lota lotaL.) heart

Opening and closing of the cardiac ryanodine (Ry) receptor (RyR) are coordinated by the free intracellular Ca2+concentration, thus making the Ca2+binding properties of the RyR important for excitation-contraction coupling. Unlike mammalian cardiac RyRs, which lose their normal function at low temperatures, RyRs of ectothermic vertebrates remain operative at 2–4°C, as indicated by Ry sensitivity of contractile force. To investigate the mechanisms of low temperature adaptation of ectothermic RyRs, we compared Ca2+-dependent kinetics of [3H]ryanodine binding in cardiac preparations of a fish (burbot, Lota lota) and a mammal (rat). The number of ventricular [3H]ryanodine binding sites determined at 20°C was 1.54 times higher in rat than burbot heart (0.401 ± 0.039 and 0.264 ± 0.019 pmol/mg protein, respectively) ( P < 0.02), while the binding affinity ( Kd) for [3H]ryanodine was similar (3.38 ± 0.63 and 4.38 ± 1.14 nM for rat and burbot, respectively) ( P = 0.47). The high-affinity [3H]ryanodine binding to burbot and rat cardiac preparations was tightly coordinated by the free Ca2+concentration at both 20°C and 2°C and did not differ between the two species. Half-maximal [3H]ryanodine binding occurred at 0.191 ± 0.027 μM and 0.164 ± 0.034 μM Ca2+for rat and at 0.212 ± 0.035 μM and 0.188 ± 0.039 μM Ca2+for burbot ( P = 0.65), at 2°C and 20°C, respectively. In two other fish species, rainbow trout ( Oncorhynchus mykiss) and crucian carp ( Carassius carassius), the Ca2+-binding affinity at 20°C was 4.4 and 5.9 times lower, respectively, than in the burbot. At 20°C, the rate of [3H]ryanodine binding to the high-affinity binding site was similar in rat and burbot but was drastically slowed in rat at 2°C. At 2°C, [3H]ryanodine failed to dissociate from rat cardiac RyRs, and at 10°C and 20°C, the rate of dissociation was two to three times slower in rat than burbot preparations. The latter finding is compatible with a channel gating mechanism, where the closing of the Ca2+release channel is impaired or severely retarded by low temperature in rat but less so in burbot preparations. The stronger effect of low temperature on association and dissociation rate of [3H]ryanodine binding in rat compared with burbot suggests that RyRs of the ectothermic fish, unlike those of endothermic rat, are better able to open and close at low temperatures.

A 3D Monte Carlo analysis of the role of dyadic space geometry in spark generation

In multiple biological systems, vital intracellular signaling processes occur locally in minute periplasmic subspaces often referred to as signaling microdomains. The number of signaling molecules in these microdomains is small enough to render the notion of continuous concentration changes invalid, such that signaling events are better described using stochastic rather than deterministic methods. Of particular interest is the dyadic cleft in the cardiac myocyte, where short-lived, local increases in intracellular Ca2+ known as Ca2+ sparks regulate excitation-contraction coupling. The geometry of dyadic spaces can alter in disease and development and display significant interspecies variability. We created and studied a 3D Monte Carlo model of the dyadic cleft, specifying the spatial localization of L-type Ca2+ channels and ryanodine receptors. Our analysis revealed how reaction specificity and efficiency are regulated by microdomain geometry as well as the physical separation of signaling molecules into functional complexes. The spark amplitude and rise time were found to be highly dependent on the concentration of activated channels per dyadic cleft and on the intermembrane separation, but not very sensitive to other cleft dimensions. The role of L-type Ca2+ channel and ryanodine receptor phosphorylation was also examined. We anticipate that this modeling approach may be applied to other systems (e.g., neuronal growth cones and chemotactic cells) to create a general description of stochastic events in Ca2+ signaling.