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

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

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

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

Mitochondria in cardiomyocyte Ca signaling

Ca²⁺ signaling is of vital importance to cardiac cell function and plays an important role in heart failure. It is based on sarcolemmal, sarcoplasmic reticulum and mitochondrial Ca²⁺ cycling. While the first two are well characterized, the latter remains unclear, controversial and technically challenging.In mammalian cardiac myocytes, Ca²⁺ influx through L-type calcium channels in the sarcolemmal membrane triggers Ca²⁺ release from the nearby junctional sarcoplasmic reticulum to produce Ca²⁺ sparks. When this triggering is synchronized by the cardiac action potential, a global [Ca²⁺]_i transient arises from coordinated Ca²⁺ release events. The ends of intermyofibrillar mitochondria are located within 20 nm of the junctional sarcoplasmic reticulum and thereby experience a high local [Ca²⁺] during the Ca²⁺ release process. Both local and global Ca²⁺ signals may thus influence calcium signaling in mitochondria and, reciprocally, mitochondria may contribute to the local control of calcium signaling. In addition to the intermyofibrillar mitochondria, morphologically distinct mitochondria are also located in the perinuclear and subsarcolemmal regions of the cardiomyocyte and thus experience a different local [Ca²⁺].Here we review the literature in regard to several issues of broad interest: (1) the ultrastructural basis for mitochondrion – sarcoplasmic reticulum cross-signaling; (2) mechanisms of sarcoplasmic reticulum signaling; (3) mitochondrial calcium signaling; and (4) the possible interplay of calcium signaling between the sarcoplasmic reticulum and adjacent mitochondria.Finally, this review discusses experimental findings and mathematical models of cardiac calcium signaling between the sarcoplasmic reticulum and mitochondria, identifies weaknesses in these models, and suggests strategies and approaches for future investigations.     

Thermal modulation of epicardial Ca2+ dynamics uncovers molecular mechanisms of Ca2+ alternans

Ca²⁺ alternans (Ca-Alts) are alternating beat-to-beat changes in the amplitude of Ca²⁺ transients that frequently occur during tachycardia, ischemia, or hypothermia that can lead to sudden cardiac death. Ca-Alts appear to result from a variation in the amount of Ca²⁺ released from the sarcoplasmic reticulum (SR) between two consecutive heartbeats. This variable Ca²⁺ release has been attributed to the alternation of the action potential duration, delay in the recovery from inactivation of RYR Ca²⁺ release channel (RYR2), or an incomplete Ca²⁺ refilling of the SR. In all three cases, the RYR2 mobilizes less Ca²⁺ from the SR in an alternating manner, thereby generating an alternating profile of the Ca²⁺ transients. We used a new experimental approach, fluorescence local field optical mapping (FLOM), to record at the epicardial layer of an intact heart with subcellular resolution. In conjunction with a local cold finger, a series of images were recorded within an area where the local cooling induced a temperature gradient. Ca-Alts were larger in colder regions and occurred without changes in action potential duration. Analysis of the change in the enthalpy and Q₁₀ of several kinetic processes defining intracellular Ca²⁺ dynamics indicated that the effects of temperature change on the relaxation of intracellular Ca²⁺ transients involved both passive and active mechanisms. The steep temperature dependency of Ca-Alts during tachycardia suggests Ca-Alts are generated by insufficient SERCA-mediated Ca²⁺ uptake into the SR. We found that Ca-Alts are heavily dependent on intra-SR Ca²⁺ and can be promoted through partial pharmacologic inhibition of SERCA2a. Finally, the FLOM experimental approach has the potential to help us understand how arrhythmogenesis correlates with the spatial distribution of metabolically impaired myocytes along the myocardium.     

Cytoplasmic Ca2+ does not inhibit the cardiac muscle sarcoplasmic reticulum ryanodine receptor Ca2+ channel,although Ca2+-induced Ca2+ inactivation of Ca2+ release is observed in native vesicles

Single channel properties of cardiac and fast-twitch skeletal muscle sarcoplasmic reticulum (SR) release channels were compared in a planar bilayer by fusing SR membranes in a Cs⁺-conducting medium. We found that the pharmacology, Cs⁺ conductance and selectivity to monovalent and divalent cations of the two channels were similar. The cardiac SR channel exhibited multiple kinetic states. The open and closed lifetimes were not altered from a range of 10^–7 to 10^–3 M Ca²⁺, but the proportion of closed and open states shifted to shorter closings and openings, respectively.However, while the single channel activity of the skeletal SR channel was activated and inactivated by micromolar and millimolar Ca²⁺, respectively, the cardiac SR channel remained activated in the presence of high [Ca²⁺]. In correlation to these studies, [³H]ryanodine binding by the receptors of the two channel receptors was inhibited by high [Ca²⁺] in skeletal but not in cardiac membranes in the presence of adenine nucleotides. There is, however, a minor inhibition of [³H]ryanodine binding of cardiac SR at millimolar Ca²⁺ in the absence of adenine nucleotides.When Ca²⁺-induced Ca²⁺ release was examined from preloaded native SR vesicles, the release rates followed a normal biphasic curve, with Ca²⁺-induced inactivation at high [Ca²⁺] for both cardiac and skeletal SR. Our data suggest that the molecular basis of regulation of the SR Ca²⁺ release channel in cardiac and skeletal muscle is different, and that the cardiac SR channel isoform lacks a Ca²⁺-inactivated site.This work was supported by research grants from the National Institutes of Health HL13870 and AR38970, and the Texas Affiliate of the American Heart Association, 91A-188. M. Fill was the recipient of an NIH fellowship AR01834.     

Transport of Ca2+ from Sarcoplasmic Reticulum to Mitochondria in Rat Ventricular Myocytes

Studies with electron microscopy have shown that sarcoplasmic reticulum (SR) andmitochondria locate close to each other in cardiac muscle cells. We investigated the hypothesis thatthis proximity results in a transient exposure of mitochondrial Ca²⁺ uniporter (CaUP) to highconcentrations of Ca²⁺ following Ca²⁺ release from the SR and thus an influx of Ca²⁺into mitochondria. Single ventricular myocytes of rat were skinned by exposing them to aphysiological solution containing saponin (0.2 mg/ml). Cytosolic Ca²⁺ concentration ([Ca²⁺]_c)and mitochondrial Ca²⁺ concentration ([Ca²⁺]_m) were measured with fura-2 and rhod2,respectively. Application of caffeine (10 mM) induced a concomitant increase in[Ca²⁺]_c and [Ca²⁺]_m.Ruthenium red, at concentrations that block CaUP but not SR release, diminished thecaffeine-induced increase in [Ca²⁺]_m but not[Ca²⁺]_c. In the presence of 1 mM BAPTA, a Ca²⁺ chelator,the caffeine-induced increase in [Ca²⁺]_m was reduced substantially less than [Ca²⁺]_c. Moreover,inhibition of SR Ca²⁺ pump with two different concentrations of thapsigargin caused anincrease in [Ca²⁺]_m, which was related to the rate of [Ca²⁺]_c increase. Finally, electronmicroscopy showed that sites of junctions between SR and T tubules from which Ca²⁺ is released,or Ca²⁺ release units, CRUs, are preferentially located in close proximity to mitochondria.The distance between individual SR Ca²⁺ release channels (feet or ryanodine receptors) isvery short, ranging between approximately 37 and 270 nm. These results are consistent withthe idea that there is a preferential coupling of Ca²⁺ transport from SR to mitochondria incardiac muscle cells, because of their structural proximity.     

Molecular architecture of Ca2+ signaling control in muscle and heart cells

Ca²⁺ signaling in skeletal and cardiac muscles is a bi-directional process that involves cross-talk between signaling molecules in the sarcolemmal membrane and Ca²⁺ release machinery in the intracellular organelles. Maintenance of a junctional membrane structure between the sarcolemmal membrane and the sarcoplasmic reticulum (SR) provides a framework for the conversion of action potential arrived at the sarcolemma into release of Ca²⁺ from the SR, leading to activation of a variety of physiological processes. Activity-dependent changes in Ca²⁺ storage inside the SR provides a retrograde signal for the activation of store-operated Ca²⁺ channel (SOC) on the sarcolemmal membrane, which plays important roles in the maintenance of Ca²⁺ homeostasis in physiology and pathophysiology. Research progress during the last 30 years had advanced our understanding of the cellular and molecular mechanisms for the control of Ca²⁺ signaling in muscle and cardiovascular physiology. Here we summarize the functions of three key molecules that are located in the junctional membrane complex of skeletal and cardiac muscle cells: junctophilin as a “glue” that physiologically links the SR membrane to the sarcolemmal membrane for formation of the junctional membrane framework, mitsugumin29 as a muscle-specific synaptophysin family protein that contributes to maintain the coordinated Ca²⁺ signaling in skeletal muscle, and TRIC as a novel cation-selective channel located on the SR membrane that provides counter-ion current during the rapid process of Ca²⁺ release from the SR.     

Kinetics on Demand Is a Simple Mathematical Solution that Fits Recorded Caffeine-Induced Luminal SR Ca2+ Changes in Smooth Muscle Cells

The process of Ca²⁺ release from sarcoplasmic reticulum (SR) comprises 4 phases in smooth muscle cells. Phase 1 is characterized by a large increase of the intracellular Ca²⁺ concentration ([Ca²⁺]_i) with a minimal reduction of the free luminal SR [Ca²⁺] ([Ca²⁺]_FSR). Importantly, active SR Ca²⁺ ATPases (SERCA pumps) are necessary for phase 1 to occur. This situation cannot be explained by the standard kinetics that involves a fixed amount of luminal Ca²⁺ binding sites. A new mathematical model was developed that assumes an increasing SR Ca²⁺ buffering capacity in response to an increase of the luminal SR [Ca²⁺] that is called Kinetics-on-Demand (KonD) model. This approach can explain both phase 1 and the refractory period associated with a recovered [Ca²⁺]_FSR. Additionally, our data suggest that active SERCA pumps are a requisite for KonD to be functional; otherwise luminal SR Ca²⁺ binding proteins switch to standard kinetics. The importance of KonD Ca²⁺ binding properties is twofold: a more efficient Ca²⁺ release process and that [Ca²⁺]_FSR and Ca²⁺-bound to SR proteins ([Ca²⁺]_BSR) can be regulated separately allowing for Ca²⁺ release to occur (provided by Ca²⁺-bound to luminal Ca²⁺ binding proteins) without an initial reduction of the [Ca²⁺]_FSR.     

The Role of Dyadic Organization in Regulation of Sarcoplasmic Reticulum Ca Handling during Rest in Rabbit Ventricular Myocytes

The dyadic organization of ventricular myocytes ensures synchronized activation of sarcoplasmic reticulum (SR) Ca²⁺ release during systole. However, it remains obscure how the dyadic organization affects SR Ca²⁺ handling during diastole. By measuring intraluminal SR Ca²⁺ ([Ca²⁺]_SR) decline during rest in rabbit ventricular myocytes, we found that ∼76% of leaked SR Ca²⁺ is extruded from the cytosol and only ∼24% is pumped back into the SR. Thus, the majority of Ca²⁺ that leaks from the SR is removed from the cytosol before it can be sequestered back into the SR by the SR Ca²⁺-ATPase (SERCA). Detubulation decreased [Ca²⁺]_SR decline during rest, thus making the leaked SR Ca²⁺ more accessible for SERCA. These results suggest that Ca²⁺ extrusion systems are localized in T-tubules. Inhibition of Na⁺-Ca²⁺ exchanger (NCX) slowed [Ca²⁺]_SR decline during rest by threefold, however did not prevent it. Depolarization of mitochondrial membrane potential during NCX inhibition completely prevented the rest-dependent [Ca²⁺]_SR decline. Despite a significant SR Ca²⁺ leak, Ca²⁺ sparks were very rare events in control conditions. NCX inhibition or detubulation increased Ca²⁺ spark activity independent of SR Ca²⁺ load. Overall, these results indicate that during rest NCX effectively competes with SERCA for cytosolic Ca²⁺ that leaks from the SR. This can be explained if the majority of SR Ca²⁺ leak occurs through ryanodine receptors in the junctional SR that are located closely to NCX in the dyadic cleft. Such control of the dyadic [Ca²⁺] by NCX play a critical role in suppressing Ca²⁺ sparks during rest.     

The Role of Dyadic Organization in Regulation of Sarcoplasmic Reticulum Ca2+ Handling during Rest in Rabbit Ventricular Myocytes

The dyadic organization of ventricular myocytes ensures synchronized activation of sarcoplasmic reticulum (SR) Ca²⁺ release during systole. However, it remains obscure how the dyadic organization affects SR Ca²⁺ handling during diastole. By measuring intraluminal SR Ca²⁺ ([Ca²⁺]_SR) decline during rest in rabbit ventricular myocytes, we found that ∼76% of leaked SR Ca²⁺ is extruded from the cytosol and only ∼24% is pumped back into the SR. Thus, the majority of Ca²⁺ that leaks from the SR is removed from the cytosol before it can be sequestered back into the SR by the SR Ca²⁺-ATPase (SERCA). Detubulation decreased [Ca²⁺]_SR decline during rest, thus making the leaked SR Ca²⁺ more accessible for SERCA. These results suggest that Ca²⁺ extrusion systems are localized in T-tubules. Inhibition of Na⁺-Ca²⁺ exchanger (NCX) slowed [Ca²⁺]_SR decline during rest by threefold, however did not prevent it. Depolarization of mitochondrial membrane potential during NCX inhibition completely prevented the rest-dependent [Ca²⁺]_SR decline. Despite a significant SR Ca²⁺ leak, Ca²⁺ sparks were very rare events in control conditions. NCX inhibition or detubulation increased Ca²⁺ spark activity independent of SR Ca²⁺ load. Overall, these results indicate that during rest NCX effectively competes with SERCA for cytosolic Ca²⁺ that leaks from the SR. This can be explained if the majority of SR Ca²⁺ leak occurs through ryanodine receptors in the junctional SR that are located closely to NCX in the dyadic cleft. Such control of the dyadic [Ca²⁺] by NCX play a critical role in suppressing Ca²⁺ sparks during rest.     

Control of Sarcoplasmic Reticulum Ca2+ Release by Stochastic RyR Gating within a 3D Model of the Cardiac Dyad and Importance of Induction Decay for CICR Termination

The factors responsible for the regulation of regenerative calcium-induced calcium release (CICR) during Ca²⁺ spark evolution remain unclear. Cardiac ryanodine receptor (RyR) gating in rats and sheep was recorded at physiological Ca²⁺, Mg²⁺, and ATP levels and incorporated into a 3D model of the cardiac dyad, which reproduced the time course of Ca²⁺ sparks, Ca²⁺ blinks, and Ca²⁺ spark restitution. The termination of CICR by induction decay in the model principally arose from the steep Ca²⁺ dependence of RyR closed time, with the measured sarcoplasmic reticulum (SR) lumen Ca²⁺ dependence of RyR gating making almost no contribution. The start of CICR termination was strongly dependent on the extent of local depletion of junctional SR Ca²⁺, as well as the time course of local Ca²⁺ gradients within the junctional space. Reducing the dimensions of the dyad junction reduced Ca²⁺ spark amplitude by reducing the strength of regenerative feedback within CICR. A refractory period for Ca²⁺ spark initiation and subsequent Ca²⁺ spark amplitude restitution arose from 1), the extent to which the regenerative phase of CICR can be supported by the partially depleted junctional SR, and 2), the availability of releasable Ca²⁺ in the junctional SR. The physical organization of RyRs within the junctional space had minimal effects on Ca²⁺ spark amplitude when more than nine RyRs were present. Spark amplitude had a nonlinear dependence on RyR single-channel Ca²⁺ flux, and was approximately halved by reducing the flux from 0.6 to 0.2 pA. Although rat and sheep RyRs had quite different Ca²⁺ sensitivities, Ca²⁺ spark amplitude was hardly affected. This suggests that moderate changes in RyR gating by second-messenger systems will principally alter the spatiotemporal properties of SR release, with smaller effects on the amount released.     

A quantitative model for linking Na+/Ca2+ exchanger to SERCA during refilling of the sarcoplasmic reticulum to sustain [Ca2+] oscillations in vascular smooth muscle

We have developed a quantitative model for the creation of cytoplasmic Ca²⁺ gradients near the inner surface of the plasma membrane (PM). In particular we simulated the refilling of the sarcoplasmic reticulum (SR) via PM–SR junctions during asynchronous [Ca²⁺]_i oscillations in smooth muscle cells of the rabbit inferior vena cava. We have combined confocal microscopy data on the [Ca²⁺]_i oscillations, force transduction data from cell contraction studies and electron microscopic images to build a basis for computational simulations that model the transport of calcium ions from Na⁺/Ca²⁺ exchangers (NCX) on the PM to sarcoplasmic/endoplasmic reticulum Ca²⁺ ATPase (SERCA) pumps on the SR as a three-dimensional random walk through the PM–SR junctional cytoplasmic spaces. Electron microscopic ultrastructural images of the smooth muscle cells were elaborated with software algorithms to produce a very clear and dimensionally accurate picture of the PM–SR junctions. From this study, we conclude that it is plausible and possible for enough Ca²⁺ to pass through the PM–SR junctions to replete the SR during the regenerative Ca²⁺ release, which underlies agonist induced asynchronous Ca²⁺ oscillations in vascular smooth muscle.     

Assessment of Sarcoplasmic Reticulum Ca Depletion During Spontaneous Ca Waves in Isolated Permeabilized Rabbit Ventricular Cardiomyocytes

In this study, Ca²⁺ release due to spontaneous Ca²⁺ waves was measured both from inside the sarcoplasmic reticulum (SR) and from the cytosol of rabbit cardiomyocytes. These measurements utilized Fluo5N-AM for intra-SR Ca²⁺ from intact cells and Fluo5F in the cytosol of permeabilized cells. Restricted subcellular volumes were resolved with the use of laser-scanning confocal microscopy. Local Ca²⁺ signals during spontaneous Ca²⁺ release were compared with those induced by rapid caffeine application. The free cytoplasmic [Ca²⁺] increase during a Ca²⁺ wave was 98.1% ± 0.3% of that observed during caffeine application. Conversion to total Ca²⁺ release suggested that Ca²⁺ release from a Ca²⁺ wave was not significantly different from that released during caffeine application (104% ± 6%). In contrast, the maximum decrease in intra-SR Fluo-5N fluorescence during a Ca²⁺ wave was 82.5% ± 2.6% of that observed during caffeine application. Assuming a maximum free [Ca²⁺] of 1.1 mM, this translates to a 96.2% ± 0.8% change in intra-SR free [Ca²⁺] and a 91.7% ± 1.6% depletion of the total Ca²⁺. This equates to a minimum intra-SR free Ca²⁺ of 46 ± 7 μM during a Ca²⁺ wave. Reduction of RyR2 Ca²⁺ sensitivity by tetracaine (50 μM) reduced the spontaneous Ca²⁺ release frequency while increasing the Ca²⁺ wave amplitude. This did not significantly change the total depletion of the SR (94.5% ± 1.1%). The calculated minimum [Ca²⁺] during these Ca²⁺ waves (87 ± 19 μM) was significantly higher than control (p < 0.05). A computational model incorporating this level of Ca²⁺ depletion during a Ca²⁺ wave mimicked the transient and sustained effects of tetracaine on spontaneous Ca²⁺ release. In conclusion, spontaneous Ca²⁺ release results in substantial but not complete local Ca²⁺ depletion of the SR. Furthermore, measurements suggest that Ca²⁺ release terminates when luminal [Ca²⁺] reaches ∼50 μM.     

Mitochondrial Ca2+-Handling in Fast Skeletal Muscle Fibers from Wild Type and Calsequestrin-Null Mice

Mitochondrial calcium handling and its relation with calcium released from sarcoplasmic reticulum (SR) in muscle tissue are subject of lively debate. In this study we aimed to clarify how the SR determines mitochondrial calcium handling using dCASQ-null mice which lack both isoforms of the major Ca²⁺-binding protein inside SR, calsequestrin. Mitochondrial free Ca²⁺-concentration ([Ca²⁺]_mito) was determined by means of a genetically targeted ratiometric FRET-based probe. Electron microscopy revealed a highly significant increase in intermyofibrillar mitochondria (+55%) and augmented coupling (+12%) between Ca²⁺ release units of the SR and mitochondria in dCASQ-null vs. WT fibers. Significant differences in the baseline [Ca²⁺]_mito were observed between quiescent WT and dCASQ-null fibers, but not in the resting cytosolic Ca²⁺ concentration. The rise in [Ca²⁺]_mito during electrical stimulation occurred in 20−30 ms, while the decline during and after stimulation was governed by 4 rate constants of approximately 40, 1.6, 0.2 and 0.03 s⁻¹. Accordingly, frequency-dependent increase in [Ca²⁺]_mito occurred during sustained contractions. In dCASQ-null fibers the increases in [Ca²⁺]_mito were less pronounced than in WT fibers and even lower when extracellular calcium was removed. The amplitude and duration of [Ca²⁺]_mito transients were increased by inhibition of mitochondrial Na⁺/Ca²⁺ exchanger (mNCX). These results provide direct evidence for fast Ca²⁺ accumulation inside the mitochondria, involvement of the mNCX in mitochondrial Ca²⁺-handling and a dependence of mitochondrial Ca²⁺-handling on intracellular (SR) and external Ca²⁺ stores in fast skeletal muscle fibers. dCASQ-null mice represent a model for malignant hyperthermia. The differences in structure and in mitochondrial function observed relative to WT may represent compensatory mechanisms for the disease-related reduction of calcium storage capacity of the SR and/or SR Ca²⁺-leakage.     

Ca2+ 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.     

Mechanisms by which cytoplasmic calcium wave propagation and alternans are generated in cardiac atrial myocytes lacking T-tubules-insights from a simulation study

This study investigated the mechanisms underlying the propagation of cytoplasmic calcium waves and the genesis of systolic Ca²⁺ alternans in cardiac myocytes lacking transverse tubules (t-tubules). These correspond to atrial cells of either small mammals or large mammals that have lost their t-tubules due to disease-induced structural remodeling (e.g., atrial fibrillation). A mathematical model was developed for a cluster of ryanodine receptors distributed on the cross section of a cell that was divided into 13 elements with a spatial resolution of 2 μm. Due to the absence of t-tubules, L-type Ca²⁺ channels were only located in the peripheral elements close to the cell-membrane surface and produced Ca²⁺ signals that propagated toward central elements by triggering successive Ca²⁺-induced Ca²⁺ release (CICR) via Ca²⁺ diffusion between adjacent elements. Under control conditions, the Ca²⁺ signals did not fully propagate to the central region of the cell. However, with modulation of several factors responsible for Ca²⁺ handling, such as the L-type Ca²⁺ channels (Ca²⁺ influx), SERCA pumps (sarcoplasmic reticulum (SR) Ca²⁺ uptake), and ryanodine receptors (SR Ca²⁺ release), Ca²⁺ wave propagation to the center of the cell could occur. These simulation results are consistent with previous experimental data from atrial cells of small mammals. The model further reveals that spatially functional heterogeneity in Ca²⁺ diffusion within the cell produced a steep relationship between the SR Ca²⁺ content and the cytoplasmic Ca²⁺ concentration. This played an important role in the genesis of Ca²⁺ alternans that were more obvious in central than in peripheral elements. Possible association between the occurrence of Ca²⁺ alternans and the model parameters of Ca²⁺ handling was comprehensively explored in a wide range of one- and two-parameter spaces. In addition, the model revealed a spontaneous second Ca²⁺ release in response to a single voltage stimulus pulse with SR Ca²⁺ overloading and augmented Ca²⁺ influx. This study provides what to our knowledge are new insights into the genesis of Ca²⁺ alternans and spontaneous second Ca²⁺ release in cardiac myocytes that lack t-tubules.     

Axial stretch enhances sarcoplasmic reticulum Ca leak and cellular Ca reuptake in guinea pig ventricular myocytes: Experiments and models

Cardiac cellular calcium (Ca²⁺) handling is the well-investigated mediator of excitation–contraction coupling, the process that translates cardiac electrical activation into mechanical events. The reverse—effects of mechanical stimulation on cardiomyocyte Ca²⁺ handling—are much less well understood, in particular during the inter-beat period, called ‘diastole’. We have investigated the effects of diastolic length changes, applied axially using a pair of carbon fibres attached to opposite ends of Guinea pig isolated ventricular myocytes, on the availability of Ca²⁺ in the main cellular stores (the sarcoplasmic reticulum; SR), by studying the rest-decay of SR Ca²⁺ content [Ca²⁺]_SR, and the reloading of the SR after prior depletion of Ca²⁺ from the cell.Cells were loaded with Fura-2 AM (an indicator of the cytosolic ‘free’ Ca²⁺ concentration, [Ca²⁺]_i), and pre-conditioned by field-stimulation (2 Hz) at 37 °C, while [Ca²⁺]_i transients and sarcomere length (SL) were recorded simultaneously. After reaching a steady state in the behaviour of observed parameters, stimulation was interrupted for between 5 and 60 s, while cells were either held at resting length, or stretched (controlled to cause a 10% increase in SL, to aid inter-individual comparison). Thereafter, each cell was returned to its original resting length, followed by swift administration of 10 mM of caffeine (in Na⁺/Ca²⁺-free solution), which causes the release of Ca²⁺ from the SR (caffeine), but largely prevents extrusion of Ca²⁺ from the cytosol to the cell exterior (Na⁺/Ca²⁺-free solution). By comparing the [Ca²⁺]_i in cells exposed/not exposed to diastolic stretch of different duration, we assessed the rest-decay dynamics of [Ca²⁺]_SR. To assess SR reloading after initial Ca²⁺ depletion, the same stretch protocol was implemented after prior emptying of the cell by application of 10 mM of caffeine in normal Tyrode solution (which causes Ca²⁺ to be released from the SR and extruded from the cell via the Na⁺/Ca²⁺ exchanger; NCX).Axial stretch enhanced the rate of both rest-decay and reloading of [Ca²⁺]_SR. Application of 40 μM streptomycin, a blocker of stretch-activated ion channels, did not affect the stretch-induced increase in SR reloading. This behaviour was reproduced in a computer simulation study, using a modified version of the 2006 Iribe–Kohl–Noble model of single cardiac myocyte Ca²⁺ handling, suggesting that stretch increases both Ca²⁺ leak from the SR and Ca²⁺ influx via the sarcolemma. This may have important implications for the mobilisation of Ca²⁺ in stretched cells, and could contribute to the regional ‘matching’ of individual cardiomyocyte contractility to dynamic, and regionally varying, changes in mechanical loads, such as diastolic pre-load, of cardiac tissue.     

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

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

Calcium buffering properties of sarcoplasmic reticulum and calcium-induced Ca2+ release during the quasi-steady level of release in twitch fibers from frog skeletal muscle

Experiments were performed to characterize the properties of the intrinsic Ca²⁺ buffers in the sarcoplasmic reticulum (SR) of cut fibers from frog twitch muscle. The concentrations of total and free calcium ions within the SR ([Ca_T]_SR and [Ca²⁺]_SR) were measured, respectively, with the EGTA/phenol red method and tetramethylmurexide (a low affinity Ca²⁺ indicator). Results indicate SR Ca²⁺ buffering was consistent with a single cooperative-binding component or a combination of a cooperative-binding component and a linear binding component accounting for 20% or less of the bound Ca²⁺. Under the assumption of a single cooperative-binding component, the most likely resting values of [Ca²⁺]_SR and [Ca_T]_SR are 0.67 and 17.1 mM, respectively, and the dissociation constant, Hill coefficient, and concentration of the Ca-binding sites are 0.78 mM, 3.0, and 44 mM, respectively. This information can be used to calculate a variable proportional to the Ca²⁺ permeability of the SR, namely d[Ca_T]_SR/dt ÷ [Ca²⁺]_SR (denoted release permeability), in experiments in which only [Ca_T]_SR or [Ca²⁺]_SR is measured. In response to a voltage-clamp step to −20 mV at 15°C, the release permeability reaches an early peak followed by a rapid decline to a quasi-steady level that lasts ∼50 ms, followed by a slower decline during which the release permeability decreases by at least threefold. During the quasi-steady level of release, the release amplitude is 3.3-fold greater than expected from voltage activation alone, a result consistent with the recruitment by Ca-induced Ca²⁺ release of 2.3 SR Ca²⁺ release channels neighboring each channel activated by its associated voltage sensor. Release permeability at −60 mV increases as [Ca_T]_SR decreases from its resting physiological level to ∼0.1 of this level. This result argues against a release termination mechanism proposed in mammalian muscle fibers in which a luminal sensor of [Ca²⁺]_SR inhibits release when [Ca_T]_SR declines to a low level.     

Examination of the Effects of Heterogeneous Organization of RyR Clusters,Myofibrils and Mitochondria on Ca2+ Release Patterns in Cardiomyocytes

Spatio-temporal dynamics of intracellular calcium, [Ca²⁺]_i, regulate the contractile function of cardiac muscle cells. Measuring [Ca²⁺]_i flux is central to the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease. However, current imaging techniques are limited in the spatial resolution to which changes in [Ca²⁺]_i can be detected. Using spatial point process statistics techniques we developed a novel method to simulate the spatial distribution of RyR clusters, which act as the major mediators of contractile Ca²⁺ release, upon a physiologically-realistic cellular landscape composed of tightly-packed mitochondria and myofibrils. We applied this method to computationally combine confocal-scale (~ 200 nm) data of RyR clusters with 3D electron microscopy data (~ 30 nm) of myofibrils and mitochondria, both collected from adult rat left ventricular myocytes. Using this hybrid-scale spatial model, we simulated reaction-diffusion of [Ca²⁺]_i during the rising phase of the transient (first 30 ms after initiation). At 30 ms, the average peak of the simulated [Ca²⁺]_i transient and of the simulated fluorescence intensity signal, F/F₀, reached values similar to that found in the literature ([Ca²⁺]_i ≈1 μM; F/F₀≈5.5). However, our model predicted the variation in [Ca²⁺]_i to be between 0.3 and 12.7 μM (~3 to 100 fold from resting value of 0.1 μM) and the corresponding F/F₀ signal ranging from 3 to 9.5. We demonstrate in this study that: (i) heterogeneities in the [Ca²⁺]_i transient are due not only to heterogeneous distribution and clustering of mitochondria; (ii) but also to heterogeneous local densities of RyR clusters. Further, we show that: (iii) these structure-induced heterogeneities in [Ca²⁺]_i can appear in line scan data. Finally, using our unique method for generating RyR cluster distributions, we demonstrate the robustness in the [Ca²⁺]_i transient to differences in RyR cluster distributions measured between rat and human cardiomyocytes.