The structure and functioning of the couplon in the mammalian cardiomyocyte

Cardiomyocyte ryanodine receptor clusters expand and coalesce after application of isoproterenol

Earlier work has shown that ventricular ryanodine receptors (RyR2) within a cluster rearrange on phosphorylation as well as with a number of other stimuli. Using dSTORM, we investigated the effects of 300 nmol/liter isoproterenol on RyR2 clusters. In rat ventricular cardiomyocytes, there was a symmetrical enlargement of RyR2 cluster areas, a decrease in the edge-to-edge nearest neighbor distance, and distribution changes that suggested movement to increase the cluster areas by coalescence. The surface area covered by the phosphorylated clusters was significantly greater than in the control cells, as was the cluster density. This latter change was accompanied by a decreased cluster fragmentation, implying that new tetramers were brought into the sarcoplasmic reticulum. We propose a possible mechanism to explain these changes. We also visualized individual RyR2 tetramers and confirmed our earlier electron-tomographic finding that the tetramers are in a disorganized but non-random array occupying about half of the cluster area. Multiclusters, cluster groups defined by the maximum distance between their members, were analyzed for various distances. At 100 nm, the areas occupied by the multiclusters just exceeded those of the single clusters, and more than half of the multiclusters had only a single subcluster that could initiate a spark. Phosphorylation increased the size of the multiclusters, markedly so for distances >100 nm. There was no relationship between the number of subclusters in a group and the area covered by it. We conclude that isoproterenol induces rapid, significant, changes in the molecular architecture of excitation-contraction coupling.

Caveolin-3 and Arrhythmias: Insights into the Molecular Mechanisms

Caveolin-3 is a muscle-specific protein on the membrane of myocytes correlated with a variety of cardiovascular diseases. It is now clear that the caveolin-3 plays a critical role in the cardiovascular system and a significant role in cardiac protective signaling. Mutations in the gene encoding caveolin-3 cause a broad spectrum of clinical phenotypes, ranging from persistent elevations in the serum levels of creatine kinase in asymptomatic humans to cardiomyopathy. The influence of Caveolin-3(CAV-3) mutations on current density parallels the effect on channel trafficking. For example, mutations in the CAV-3 gene promote ventricular arrhythmogenesis in long QT syndrome 9 by a combined decrease in the loss of the inward rectifier current (I_K1) and gain of the late sodium current (I_Na-L). The functional significance of the caveolin-3 has proved that caveolin-3 overexpression or knockdown contributes to the occurrence and development of arrhythmias. Caveolin-3 overexpression could lead to reduced diastolic spontaneous Ca²⁺ waves, thus leading to the abnormal L-Type calcium channel current-induced ventricular arrhythmias. Moreover, CAV-3 knockdown resulted in a shift to more negative values in the hyperpolarization-activated cyclic nucleotide channel 4 current (I_HCN4) activation curve and a significant decrease in I_HCN4 whole-cell current density. Recent evidence indicates that caveolin-3 plays a significant role in adipose tissue and is related to obesity development. The role of caveolin-3 in glucose homeostasis has attracted increasing attention. This review highlights the underlining mechanisms of caveolin-3 in arrhythmia. Progress in this field may contribute to novel therapeutic approaches for patients prone to developing arrhythmia.

Triiodothyronine maintains cardiac transverse-tubule structure and function

Subclinical hypothyroidism and low T3 syndrome are commonly associated with an increased risk of cardiovascular disease (CVD) and mortality. We examined effects of T3 on T-tubule (TT) structures, Ca²⁺ mobilization and contractility, and clustering of dyadic proteins. Thyroid hormone (TH) deficiency was induced in adult female rats by propyl-thiouracil (PTU; 0.025%) treatment for 8 weeks. Rats were then randomized to continued PTU or triiodo-L-thyronine (T3; 10 μg/kg/d) treatment for 2 weeks (PTU + T3). After in vivo echocardiographic and hemodynamic recordings, cardiomyocytes (CM) were isolated to record Ca²⁺ transients and contractility. TT organization was assessed by confocal microscopy, and STORM images were captured to measure ryanodine receptor (RyR2) cluster number and size, and L-type Ca²⁺ channel (LTCC, Ca_v1.2) co-localization. Expressed genes including two integral TT proteins, junctophilin-2 (Jph-2) and bridging integrator-1 (BIN1), were analyzed in left ventricular (LV) tissues and cultured CM using qPCR and RNA sequencing. The T3 dosage used normalized serum T3, and reversed adverse effects of TH deficiency on in vivo measures of cardiac function. Recordings of isolated CM indicated that T3 increased rates of Ca²⁺ release and re-uptake, resulting in increased velocities of sarcomere shortening and re-lengthening. TT periodicity was significantly decreased, with reduced transverse tubules but increased longitudinal tubules in TH-deficient CMs and LV tissue, and these structures were normalized by T3 treatment. Analysis of STORM data of PTU myocytes showed decreased RyR2 cluster numbers and RyR localizations within each cluster without significant changes in Ca_v1.2 localizations within RyR clusters. T3 treatment normalized RyR2 cluster size and number. qPCR and RNAseq analyses of LV and cultured CM showed that Jph2 expression was T3-responsive, and its increase with treatment may explain improved TT organization and RyR-LTCC coupling.

Nano-scale morphology of cardiomyocyte t-tubule/sarcoplasmic reticulum junctions revealed by ultra-rapid high-pressure freezing and electron tomography

Detailed knowledge of the ultrastructure of intracellular compartments is a prerequisite for our understanding of how cells function. In cardiac muscle cells, close apposition of transverse (t)-tubule (TT) and sarcoplasmic reticulum (SR) membranes supports stable high-gain excitation-contraction coupling. Here, the fine structure of this key intracellular element is examined in rabbit and mouse ventricular cardiomyocytes, using ultra-rapid high-pressure freezing (HPF, omitting aldehyde fixation) and electron microscopy. 3D electron tomograms were used to quantify the dimensions of TT, terminal cisternae of the SR, and the space between SR and TT membranes (dyadic cleft). In comparison to conventional aldehyde-based chemical sample fixation, HPF-preserved samples of both species show considerably more voluminous SR terminal cisternae, both in absolute dimensions and in terms of junctional SR to TT volume ratio. In rabbit cardiomyocytes, the average dyadic cleft surface area of HPF and chemically fixed myocytes did not differ, but cleft volume was significantly smaller in HPF samples than in conventionally fixed tissue; in murine cardiomyocytes, the dyadic cleft surface area was higher in HPF samples with no difference in cleft volume. In both species, the apposition of the TT and SR membranes in the dyad was more likely to be closer than 10 nm in HPF samples compared to CFD, presumably resulting from avoidance of sample shrinkage associated with conventional fixation techniques. Overall, we provide a note of caution regarding quantitative interpretation of chemically-fixed ultrastructures, and offer novel insight into cardiac TT and SR ultrastructure with relevance for our understanding of cardiac physiology.

Physiological and Pathological Relevance of Selective and Nonselective Ca2+ Channels in Skeletal and Cardiac Muscle

Contraction of the striated muscle is fundamental for human existence. The action of voluntary skeletal muscle enables activities such as breathing, establishing body posture, and diverse body movements. Additionally, highly precise motion empowers communication, artistic expression, and other activities that define everyday human life. The involuntary contraction of striated muscle is the core function of the heart and is essential for blood flow. Several ion channels are important in the transduction of action potentials to cytosolic Ca2+ signals that enable muscle contraction; however, other ion channels are involved in the progression of muscle pathologies that can impair normal life or threaten it. This chapter describes types of selective and nonselective Ca2+ permeable ion channels expressed in the striated muscle, their participation in different aspects of muscle excitation and contraction, and their relevance to the progression of some pathological states.

Distribution of cardiac sodium channels in clusters potentiates ephaptic interactions in the intercalated disc

KEY POINTS

It has been proposed that ephaptic conduction, relying on interactions between the sodium (Na⁺ ) current and the extracellular potential in intercalated discs, might contribute to cardiac conduction when gap junctional coupling is reduced, but this mechanism is still controversial. In intercalated discs, Na⁺ channels form clusters near gap junction plaques, but the functional significance of these clusters has never been evaluated. In HEK cells expressing cardiac Na⁺ channels, we show that restricting the extracellular space modulates the Na⁺ current, as predicted by corresponding simulations accounting for ephaptic effects. In a high-resolution model of the intercalated disc, clusters of Na⁺ channels that face each other across the intercellular cleft facilitate ephaptic impulse transmission when gap junctional coupling is reduced. Thus, our simulations reveal a functional role for the clustering of Na⁺ channels in intercalated discs, and suggest that rearrangement of these clusters in disease may influence cardiac conduction.

ABSTRACT

It has been proposed that ephaptic interactions in intercalated discs, mediated by extracellular potentials, contribute to cardiac impulse propagation when gap junctional coupling is reduced. However, experiments demonstrating ephaptic effects on the cardiac Na⁺ current (I_Na ) are scarce. Furthermore, Na⁺ channels form clusters around gap junction plaques, but the electrophysiological significance of these clusters has never been investigated. In patch clamp experiments with HEK cells stably expressing human Na_v 1.5 channels, we examined how restricting the extracellular space modulates I_Na elicited by an activation protocol. In parallel, we developed a high-resolution computer model of the intercalated disc to investigate how the distribution of Na⁺ channels influences ephaptic interactions. Approaching the HEK cells to a non-conducting obstacle always increased peak I_Na at step potentials near the threshold of I_Na activation and decreased peak I_Na at step potentials far above threshold (7 cells, P = 0.0156, Wilcoxon signed rank test). These effects were consistent with corresponding control simulations with a uniform Na⁺ channel distribution. In the intercalated disc computer model, redistributing the Na⁺ channels into a central cluster of the disc potentiated ephaptic effects. Moreover, ephaptic impulse transmission from one cell to another was facilitated by clusters of Na⁺ channels facing each other across the intercellular cleft when gap junctional coupling was reduced. In conclusion, our proof-of-principle experiments demonstrate that confining the extracellular space modulates cardiac I_Na , and our simulations reveal the functional role of the aggregation of Na⁺ channels in the perinexus. These findings highlight novel concepts in the physiology of cardiac excitation.

The role of spatial organization of Ca2+ release sites in the generation of arrhythmogenic diastolic Ca2+ release in myocytes from failing hearts

In heart failure (HF), dysregulated cardiac ryanodine receptors (RyR2) contribute to the generation of diastolic Ca2+ waves (DCWs), thereby predisposing adrenergically stressed failing hearts to life-threatening arrhythmias. However, the specific cellular, subcellular, and molecular defects that account for cardiac arrhythmia in HF remain to be elucidated. Patch-clamp techniques and confocal Ca2+ imaging were applied to study spatially defined Ca2+ handling in ventricular myocytes isolated from normal (control) and failing canine hearts. Based on their activation time upon electrical stimulation, Ca2+ release sites were categorized as coupled, located in close proximity to the sarcolemmal Ca2+ channels, and uncoupled, the Ca2+ channel-free non-junctional Ca2+ release units. In control myocytes, stimulation of β-adrenergic receptors with isoproterenol (Iso) resulted in a preferential increase in Ca2+ spark rate at uncoupled sites. This site-specific effect of Iso was eliminated by the phosphatase inhibitor okadaic acid, which caused similar facilitation of Ca2+ sparks at coupled and uncoupled sites. Iso-challenged HF myocytes exhibited increased predisposition to DCWs compared to control myocytes. In addition, the overall frequency of Ca2+ sparks was increased in HF cells due to preferential stimulation of coupled sites. Furthermore, coupled sites exhibited accelerated recovery from functional refractoriness in HF myocytes compared to control myocytes. Spatially resolved subcellular Ca2+ mapping revealed that DCWs predominantly originated from coupled sites. Inhibition of CaMKII suppressed DCWs and prevented preferential stimulation of coupled sites in Iso-challenged HF myocytes. These results suggest that CaMKII- (and phosphatase)-dependent dysregulation of junctional Ca2+ release sites contributes to Ca2+-dependent arrhythmogenesis in HF.

Ultrastructural Analysis of Self-Associated RyR2s

In heart, type-2 ryanodine receptor (RyR2) forms discrete supramolecular clusters in the sarcoplasmic reticulum known as calcium release units (CRUs), which are responsible for most of the Ca(2+) released for muscle contraction. To learn about the substructure of the CRU, we sought to determine whether RyR2s have the ability to self-associate in the absence of other factors and if so, whether they do it in a specific manner. Purified RyR2 was negatively stained and imaged on the transmission electron microscope, and RyR2 particles closely associated were further analyzed using bias-free multivariate statistical analysis and classification. The resulting two-dimensional averages show that RyR2s can interact in two rigid, reproducible configurations: "adjoining", with two RyR2s alongside each other, and "oblique", with two partially overlapped RyR2s forming an angle of 12°. The two configurations are nearly identical under two extreme physiological Ca(2+) concentrations. Pseudo-atomic models for these two interactions indicate that the adjoining interaction involves contacts between the P1, SPRY1 and the helical domains. The oblique interaction is mediated by extensive contacts between the SPRY1 domains (domains 9) and P1 domains (domains 10) of both RyR2s and not through domain 6 as previously thought; in addition its asymmetric interface imposes steric constrains that inhibit the growth of RyR2 as a checkerboard, which is the configuration usually assumed, and generates new configurations, i.e., "branched" and "interlocked". This first, to our knowledge, structural detailed analysis of the inter-RyR2 interactions helps to understand important morphological and functional aspects of the CRU in the context of cardiac EC coupling.

Cardiac voltage-gated calcium channel macromolecular complexes

Over the past 20 years, a new field of research, called channelopathies, investigating diseases caused by ion channel dysfunction has emerged. Cardiac ion channels play an essential role in the generation of the cardiac action potential. Investigators have largely determined the physiological roles of different cardiac ion channels, but little is known about the molecular determinants of their regulation. The voltage-gated calcium channel Ca(v)1.2 shapes the plateau phase of the cardiac action potential and allows the influx of calcium leading to cardiomyocyte contraction. Studies suggest that the regulation of Ca(v)1.2 channels is not uniform in working cardiomyocytes. The notion of micro-domains containing Ca(v)1.2 channels and different calcium channel interacting proteins, called macro-molecular complex, has been proposed to explain these observations. The objective of this review is to summarize the currently known information on the Ca(v)1.2 macromolecular complexes in the cardiac cell and discuss their implication in cardiac function and disorder. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Nonuniform and variable arrangements of ryanodine receptors within mammalian ventricular couplons

RATIONALE

Single-tilt tomograms of the dyads in rat ventricular myocytes indicated that type 2 ryanodine receptors (RYR2s) were not positioned in a well-ordered array. Furthermore, the orientation and packing strategy of purified type 1 ryanodine receptors in lipid bilayers is determined by the free Mg2+ concentration. These observations led us to test the hypothesis that RYR2s within the mammalian dyad have multiple and complex arrangements.

OBJECTIVES

To determine the arrangement of RYR2 tetramers in the dyads of mammalian cardiomyocytes and the effects of physiologically and pathologically relevant factors on this arrangement.

METHODS AND RESULTS

We used dual-tilt electron tomography to produce en-face views of dyads, enabling a direct examination of RYR2 distribution and arrangement. Rat hearts fixed in situ; isolated rat cardiomyocytes permeabilized, incubated with 1 mmol/L Mg2+, and then fixed; and sections of human ventricle, all showed that the tetramer packing within a dyad was nonuniform containing a mix of checkerboard and side-by-side arrangements, as well as isolated tetramers. Both phosphorylation and 0.1 mmol/L Mg2+ moved the tetramers into a predominantly checkerboard configuration, whereas the 4 mmol/L Mg2+ induced a dense side-by-side arrangement. These changes occurred within 10 minutes of application of the stimuli.

CONCLUSIONS

The arrangement of RYR2 tetramers within the mammalian dyad is neither uniform nor static. We hypothesize that this is characteristic of the dyad in vivo and may provide a mechanism for modulating the open probabilities of the individual tetramers.

Three-Dimensional Reconstruction of Cardiac Sarcoplasmic Reticulum Reveals a Continuous Network Linking Transverse-Tubules

Rationale:

The organization of the transverse-tubular (t-t) system and relationship to the sarcoplasmic reticulum (SR) underpins cardiac excitation–contraction coupling. The architecture of the SR, and relationship with the t-ts, is not well characterized at the whole-cell level. Furthermore, little is known regarding changes to SR ultrastructure in heart failure.

Objective:

The aim of this study was to unravel interspecies differences and commonalities between the relationship of SR and t-t networks within cardiac myocytes, as well as the modifications that occur in heart failure, using a novel high-resolution 3-dimensional (3D) imaging technique.

Methods and Results:

Using serial block face imaging coupled with scanning electron microscopy and image analysis, we have generated 3D reconstructions of whole cardiomyocytes from sheep and rat left ventricle, revealing that the SR forms a continuous network linking t-ts throughout the cell in both species. In sheep, but not rat, the SR has an intimate relationship with the sarcolemma forming junctional domains. 3D reconstructions also reveal details of the sheep t-t system. Using a model of tachypacing-induced heart failure, we show that there are populations of swollen and collapsed t-ts, patches of SR tangling, and disorder with rearrangement of the mitochondria.

Conclusions:

We provide the first high-resolution 3D structure of the SR network showing that it forms a cell-wide communication pipeline facilitating Ca 2+ diffusion, buffering, and synchronicity. The distribution of the SR within the cell is related to interspecies differences in excitation–contraction coupling, and we report the first detailed analysis of SR remodeling as a result of heart failure.

Interaction of membrane/lipid rafts with the cytoskeleton: impact on signaling and function: membrane/lipid rafts, mediators of cytoskeletal arrangement and cell signaling

The plasma membrane in eukaryotic cells contains microdomains that are enriched in certain glycosphingolipids, gangliosides, and sterols (such as cholesterol) to form membrane/lipid rafts (MLR). These regions exist as caveolae, morphologically observable flask-like invaginations, or as a less easily detectable planar form. MLR are scaffolds for many molecular entities, including signaling receptors and ion channels that communicate extracellular stimuli to the intracellular milieu. Much evidence indicates that this organization and/or the clustering of MLR into more active signaling platforms depends upon interactions with and dynamic rearrangement of the cytoskeleton. Several cytoskeletal components and binding partners, as well as enzymes that regulate the cytoskeleton, localize to MLR and help regulate lateral diffusion of membrane proteins and lipids in response to extracellular events (e.g., receptor activation, shear stress, electrical conductance, and nutrient demand). MLR regulate cellular polarity, adherence to the extracellular matrix, signaling events (including ones that affect growth and migration), and are sites of cellular entry of certain pathogens, toxins and nanoparticles. The dynamic interaction between MLR and the underlying cytoskeleton thus regulates many facets of the function of eukaryotic cells and their adaptation to changing environments. Here, we review general features of MLR and caveolae and their role in several aspects of cellular function, including polarity of endothelial and epithelial cells, cell migration, mechanotransduction, lymphocyte activation, neuronal growth and signaling, and a variety of disease settings. This article is part of a Special Issue entitled: Reciprocal influences between cell cytoskeleton and membrane channels, receptors and transporters. Guest Editor: Jean Claude Hervé.

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.

Adding a new dimension to cardiac nano-architecture using electron microscopy: coupling membrane excitation to calcium signaling

Advances in microscopic imaging technologies and associated computational methods now allow descriptions of cellular anatomy to go beyond 2-dimensions, revealing new micro-domain dynamics at unprecedented resolutions. In cardiomyocytes, electron microscopy (EM) first described junctional membrane complexes between the sarcolemma and sarcoplasmic reticulum over a half-century ago. Since then, 3-dimensional EM technologies such as electron tomography have become successful in determining the realistic nano-geometry of membrane junctions (dyads and peripheral junctions) and associated structures such as transverse tubules (T-tubules, aka. T-system). Concomitantly, super-resolution light microscopy has gone beyond the diffraction-limit to determine the distribution of molecules, such as ryanodine receptors, with 10(-8) meter (10nm) order accuracy. This review provides the current structural perspective and functional interpretation of membrane junction complexes, which are the central machinery controlling cardiac excitation-contraction coupling via calcium signaling.

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.

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.

Modelling cardiac calcium sparks in a three-dimensional reconstruction of a calcium release unit

Triggered release of Ca2+ from an individual sarcoplasmic reticulum (SR) Ca(2+) release unit (CRU) is the fundamental event of cardiac excitation–contraction coupling, and spontaneous release events (sparks) are the major contributor to diastolic Ca(2+) leak in cardiomyocytes. Previous model studies have predicted that the duration and magnitude of the spark is determined by the local CRU geometry, as well as the localization and density of Ca(2+) handling proteins. We have created a detailed computational model of a CRU, and developed novel tools to generate the computational geometry from electron tomographic images. Ca(2+) diffusion was modelled within the SR and the cytosol to examine the effects of localization and density of the Na(+)/Ca(2+) exchanger, sarco/endoplasmic reticulum Ca(2+)-ATPase 2 (SERCA), and calsequestrin on spark dynamics. We reconcile previous model predictions of approximately 90% local Ca(2+) depletion in junctional SR, with experimental reports of about 40%. This analysis supports the hypothesis that dye kinetics and optical averaging effects can have a significant impact on measures of spark dynamics. Our model also predicts that distributing calsequestrin within non-junctional Z-disc SR compartments, in addition to the junctional compartment, prolongs spark release time as reported by Fluo5. By pumping Ca(2+) back into the SR during a release, SERCA is able to prolong a Ca(2+) spark, and this may contribute to SERCA-dependent changes in Ca(2+) wave speed. Finally, we show that including the Na(+)/Ca(2+) exchanger inside the dyadic cleft does not alter local [Ca(2+)] during a spark.