Ca2+-mobility in the sarcoplasmic reticulum of ventricular myocytes is low

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Adhesion to laminin-1 and collagen IV induces the formation of Ca2+ microdomains that sensitize mouse T cells for activation

During an immune response, T cells migrate from blood vessel walls into inflamed tissues by migrating across the endothelium and through extracellular matrix (ECM). Integrins facilitate T cell binding to endothelial cells and ECM proteins. Here, we report that Ca²⁺ microdomains observed in the absence of T cell receptor (TCR)/CD3 stimulation are initial signaling events triggered by adhesion to ECM proteins that increase the sensitivity of primary murine T cells to activation. Adhesion to the ECM proteins collagen IV and laminin-1 increased the number of Ca²⁺ microdomains in a manner dependent on the kinase FAK, phospholipase C (PLC), and all three inositol 1,4,5-trisphosphate receptor (IP₃R) subtypes and promoted the nuclear translocation of the transcription factor NFAT-1. Mathematical modeling predicted that the formation of adhesion-dependent Ca²⁺ microdomains required the concerted activity of two to six IP₃Rs and ORAI1 channels to achieve the increase in the Ca²⁺ concentration in the ER-plasma membrane junction that was observed experimentally and that required SOCE. Further, adhesion-dependent Ca²⁺ microdomains were important for the magnitude of the TCR-induced activation of T cells on collagen IV as assessed by the global Ca²⁺ response and NFAT-1 nuclear translocation. Thus, adhesion to collagen IV and laminin-1 sensitizes T cells through a mechanism involving the formation of Ca²⁺ microdomains, and blocking this low-level sensitization decreases T cell activation upon TCR engagement.

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.

A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation

Calcium (Ca2+) plays a central role in the excitation and contraction of cardiac myocytes. Experiments have indicated that calcium release is stochastic and regulated locally suggesting the possibility of spatially heterogeneous calcium levels in the cells. This spatial heterogeneity might be important in mediating different signaling pathways. During more than 50 years of computational cell biology, the computational models have been advanced to incorporate more ionic currents, going from deterministic models to stochastic models. While periodic increases in cytoplasmic Ca2+ concentration drive cardiac contraction, aberrant Ca2+ release can underly cardiac arrhythmia. However, the study of the spatial role of calcium ions has been limited due to the computational expense of using a three-dimensional stochastic computational model. In this paper, we introduce a three-dimensional stochastic computational model for rat ventricular myocytes at the whole-cell level that incorporate detailed calcium dynamics, with (1) non-uniform release site placement, (2) non-uniform membrane ionic currents and membrane buffers, (3) stochastic calcium-leak dynamics and (4) non-junctional or rogue ryanodine receptors. The model simulates spark-induced spark activation and spark-induced Ca2+ wave initiation and propagation that occur under conditions of calcium overload at the closed-cell condition, but not when Ca2+ levels are normal. This is considered important since the presence of Ca2+ waves contribute to the activation of arrhythmogenic currents.

Three-Dimensional Model of Sub-Plasmalemmal Ca2+ Microdomains Evoked by the Interplay Between ORAI1 and InsP3 Receptors

Ca²⁺ signaling plays an essential role in T cell activation, which is a key step to start an adaptive immune response. During the transition from a quiescent to a fully activated state, Ca²⁺ microdomains characterized by reduced spatial and temporal extents are observed in the junctions between the plasma membrane (PM) and the endoplasmic reticulum (ER). Such Ca²⁺ responses can also occur in response to T cell adhesion to other cells or extracellular matrix proteins in otherwise unstimulated T cells. These non-TCR/CD3-dependent Ca²⁺ microdomains rely on d-myo-inositol 1,4,5-trisphosphate (IP₃) signaling and subsequent store operated Ca²⁺ entry (SOCE) via the ORAI/STIM system. The detailed molecular mechanism of adhesion-dependent Ca²⁺ microdomain formation remains to be fully elucidated. We used mathematical modeling to investigate the spatiotemporal characteristics of T cell Ca²⁺ microdomains and their molecular regulators. We developed a reaction-diffusion model using COMSOL Multiphysics to describe the evolution of cytosolic and ER Ca²⁺ concentrations in a three-dimensional ER-PM junction. Equations are based on a previously proposed realistic description of the junction, which is extended to take into account IP₃ receptors (IP₃R) that are located next to the junction. The first model only considered the ORAI channels and the SERCA pumps. Taking into account the existence of preformed clusters of ORAI1 and STIM2, ORAI1 slightly opens in conditions of a full ER. These simulated Ca²⁺ microdomains are too small as compared to those observed in unstimulated T cells. When considering the opening of the IP₃Rs located near the junction, the local depletion of ER Ca²⁺ allows for larger Ca²⁺ fluxes through the ORAI1 channels and hence larger local Ca²⁺ concentrations. Computational results moreover show that Ca²⁺ diffusion in the ER has a major impact on the Ca²⁺ changes in the junction, by affecting the local Ca²⁺ gradients in the sub-PM ER. Besides pointing out the likely involvement of the spontaneous openings of IP₃Rs in the activation of SOCE in conditions of T cell adhesion prior to full activation, the model provides a tool to investigate how Ca²⁺ microdomains extent and interact in response to T cell receptor activation.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Complex patterns of subcellular cardiac alternans

Cardiac alternans, in which the membrane potential and the intracellular calcium concentration exhibit alternating durations and peak amplitudes at consecutive beats, constitute a precursor to fatal cardiac arrhythmia such as sudden cardiac death. A crucial question therefore concerns the onset of cardiac alternans. Typically, alternans are only reported when they are fully developed. Here, we present a modelling approach to explore recently discovered microscopic alternans, which represent one of the earliest manifestations of cardiac alternans. In this case, the regular periodic dynamics of the local intracellular calcium concentration is already unstable, while the whole-cell behaviour suggests a healthy cell state. In particular, we use our model to investigate the impact of calcium diffusion in both the cytosol and the sarcoplasmic reticulum on the formation of microscopic calcium alternans. We find that for dominant cytosolic coupling, calcium alternans emerge via the traditional period doubling bifurcation. In contrast, dominant luminal coupling leads to a novel route to calcium alternans through a saddle-node bifurcation at the network level. Combining semi-analytical and computational approaches, we compute areas of stability in parameter space and find that as we cross from stable to unstable regions, the emergent patterns of the intracellular calcium concentration change abruptly in a fashion that is highly dependent upon position along the stability boundary. Our results demonstrate that microscopic calcium alternans may possess a much richer dynamical repertoire than previously thought and further strengthen the role of luminal calcium in shaping cardiac calcium dynamics.

Modulation of Cardiac Alternans by Altered Sarcoplasmic Reticulum Calcium Release: A Simulation Study

Background: Cardiac alternans is an important precursor to arrhythmia, facilitating formation of conduction block, and re-entry. Diseased hearts were observed to be particularly vulnerable to alternans, mainly in heart failure or after myocardial infarction. Alternans is typically linked to oscillation of calcium cycling, particularly in the sarcoplasmic reticulum (SR). While the role of SR calcium reuptake in alternans is well established, the role of altered calcium release by ryanodine receptors has not yet been studied extensively. At the same time, there is strong evidence that calcium release is abnormal in heart failure and other heart diseases, suggesting that these changes might play a pro-alternans role.

Aims: To demonstrate how changes to intracellular calcium release dynamics and magnitude affect alternans vulnerability.

Methods: We used the state-of-the-art Heijman–Rudy and O’Hara–Rudy computer models of ventricular myocyte, given their detailed representation of calcium handling and their previous utility in alternans research. We modified the models to obtain precise control over SR release dynamics and magnitude, allowing for the evaluation of these properties in alternans formation and suppression.

Results: Shorter time to peak SR release and shorter release duration decrease alternans vulnerability by improved refilling of releasable calcium within junctional SR; conversely, slow release promotes alternans. Modulating the total amount of calcium released, we show that sufficiently increased calcium release may surprisingly prevent alternans via a mechanism linked to the functional depletion of junctional SR during release. We show that this mechanism underlies differences between “eye-type” and “fork-type” alternans, which were observed in human in vivo and in silico. We also provide a detailed explanation of alternans formation in the given computer models, termed “sarcoplasmic reticulum calcium cycling refractoriness.” The mechanism relies on the steep SR load–release relationship, combined with relatively limited rate of junctional SR refilling.

Conclusion: Both altered dynamics and magnitude of SR calcium release modulate alternans vulnerability. In particular, slow dynamics of SR release, such as those observed in heart failure, promote alternans. Therefore, acceleration of intracellular calcium release, e.g., via synchronization of calcium sparks, may inhibit alternans in failing hearts and reduce arrhythmia occurrence.

Regional acidosis locally inhibits but remotely stimulates Ca2+ waves in ventricular myocytes

Aims

Spontaneous Ca²⁺ waves in cardiomyocytes are potentially arrhythmogenic. A powerful controller of Ca²⁺ waves is the cytoplasmic H⁺ concentration ([H⁺]_i), which fluctuates spatially and temporally in conditions such as myocardial ischaemia/reperfusion. H⁺-control of Ca²⁺ waves is poorly understood. We have therefore investigated how [H⁺]_i co-ordinates their initiation and frequency.

Methods and results

Spontaneous Ca²⁺ waves were imaged (fluo-3) in rat isolated ventricular myocytes, subjected to modest Ca²⁺-overload. Whole-cell intracellular acidosis (induced by acetate-superfusion) stimulated wave frequency. Pharmacologically blocking sarcolemmal Na⁺/H⁺ exchange (NHE1) prevented this stimulation, unveiling inhibition by H⁺. Acidosis also increased Ca²⁺ wave velocity. Restricting acidosis to one end of a myocyte, using a microfluidic device, inhibited Ca²⁺ waves in the acidic zone (consistent with ryanodine receptor inhibition), but stimulated wave emergence elsewhere in the cell. This remote stimulation was absent when NHE1 was selectively inhibited in the acidic zone. Remote stimulation depended on a locally evoked, NHE1-driven rise of [Na⁺]_i that spread rapidly downstream.

Conclusion

Acidosis influences Ca²⁺ waves via inhibitory Hi+ and stimulatory Nai+ signals (the latter facilitating intracellular Ca²⁺-loading through modulation of sarcolemmal Na⁺/Ca²⁺ exchange activity). During spatial [H⁺]_i-heterogeneity, Hi+-inhibition dominates in acidic regions, while rapid Nai+ diffusion stimulates waves in downstream, non-acidic regions. Local acidosis thus simultaneously inhibits and stimulates arrhythmogenic Ca²⁺-signalling in the same myocyte. If the principle of remote H⁺-stimulation of Ca²⁺ waves also applies in multicellular myocardium, it raises the possibility of electrical disturbances being driven remotely by adjacent ischaemic areas, which are known to be intensely acidic.

Three-dimensional spatio-temporal modelling of store operated Ca2+ entry: Insights into ER refilling and the spatial signature of Ca2+ signals

The spatial organisation of Orai channels and SERCA pumps within ER-PM junctions is important for enhancing the versatility and specificity of sub-cellular Ca²⁺ signals generated during store operated Ca²⁺ entry (SOCE). In this paper, we present a novel three dimensional spatio-temporal model describing Ca²⁺ dynamics in the ER-PM junction and sub-PM ER during SOCE. We investigate the role of Orai channel and SERCA pump location to provide insights into how these components shape the Ca²⁺ signals generated and affect ER refilling. We find that the organisation of Orai channels within the ER-PM junction controls the amplitude and shape of the Ca²⁺ profile but does not enhance ER refilling. The model shows that ER refilling is only weakly affected by the location of SERCA2b pumps within the ER-PM junction and that the placement of SERCA2a pumps within the ER-PM junction has much greater control over ER refilling.

Overview of the Microenvironment of Vasculature in Vascular Tone Regulation

Hypertension is asymptomatic and a well-known “silent killer”, which can cause various concomitant diseases in human population after years of adherence. Although there are varieties of synthetic antihypertensive drugs available in current market, their relatively low efficacies and major application in only single drug therapy, as well as the undesired chronic adverse effects associated, has drawn the attention of worldwide scientists. According to the trend of antihypertensive drug evolution, the antihypertensive drugs used as primary treatment often change from time-to-time with the purpose of achieving the targeted blood pressure range. One of the major concerns that need to be accounted for here is that the signaling mechanism pathways involved in the vasculature during the vascular tone regulation should be clearly understood during the pharmacological research of antihypertensive drugs, either in vitro or in vivo. There are plenty of articles that discussed the signaling mechanism pathways mediated in vascular tone in isolated fragments instead of a whole comprehensive image. Therefore, the present review aims to summarize previous published vasculature-related studies and provide an overall depiction of each pathway including endothelium-derived relaxing factors, G-protein-coupled, enzyme-linked, and channel-linked receptors that occurred in the microenvironment of vasculature with a full schematic diagram on the ways their signals interact. Furthermore, the crucial vasodilative receptors that should be included in the mechanisms of actions study on vasodilatory effects of test compounds were suggested in the present review as well.

A model of cardiac ryanodine receptor gating predicts experimental Ca2+-dynamics and Ca2+-triggered arrhythmia in the long QT syndrome

Abnormal Ca²⁺ handling is well-established as the trigger of cardiac arrhythmia in catecholaminergic polymorphic ventricular tachycardia and digoxin toxicity, but its role remains controversial in Torsade de Pointes (TdP), the arrhythmia associated with the long QT syndrome (LQTS). Recent experimental results show that early afterdepolarizations (EADs) that initiate TdP are caused by spontaneous (non-voltage-triggered) Ca²⁺ release from Ca²⁺-overloaded sarcoplasmic reticulum (SR) rather than the activation of the L-type Ca²⁺-channel window current. In bradycardia and long QT type 2 (LQT2), a second, non-voltage triggered cytosolic Ca²⁺ elevation increases gradually in amplitude, occurs before overt voltage instability, and then precedes the rise of EADs. Here, we used a modified Shannon-Puglisi-Bers model of rabbit ventricular myocytes to reproduce experimental Ca²⁺ dynamics in bradycardia and LQT2. Abnormal systolic Ca²⁺-oscillations and EADs caused by SR Ca²⁺-release are reproduced in a modified 0-dimensional model, where 3 gates in series control the ryanodine receptor (RyR2) conductance. Two gates control RyR2 activation and inactivation and sense cytosolic Ca²⁺ while a third gate senses luminal junctional SR Ca²⁺. The model predicts EADs in bradycardia and low extracellular [K⁺] and cessation of SR Ca²⁺-release terminate salvos of EADs. Ca²⁺-waves, systolic cell-synchronous Ca²⁺-release, and multifocal diastolic Ca²⁺ release seen in subcellular Ca²⁺-mapping experiments are observed in the 2-dimensional version of the model. These results support the role of SR Ca²⁺-overload, abnormal SR Ca²⁺-release, and the subsequent activation of the electrogenic Na⁺/Ca²⁺-exchanger as the mechanism of TdP. The model offers new insights into the genesis of cardiac arrhythmia and new therapeutic strategies.

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.

Oxidative stress and ca(2+) release events in mouse cardiomyocytes

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

A human ventricular myocyte model with a refined representation of excitation-contraction coupling

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

Visualization of Ca2+ Filling Mechanisms upon Synaptic Inputs in the Endoplasmic Reticulum of Cerebellar Purkinje Cells

The endoplasmic reticulum (ER) plays crucial roles in intracellular Ca(2+) signaling, serving as both a source and sink of Ca(2+), and regulating a variety of physiological and pathophysiological events in neurons in the brain. However, spatiotemporal Ca(2+) dynamics within the ER in central neurons remain to be characterized. In this study, we visualized synaptic activity-dependent ER Ca(2+) dynamics in mouse cerebellar Purkinje cells (PCs) using an ER-targeted genetically encoded Ca(2+) indicator, G-CEPIA1er. We used brief parallel fiber stimulation to induce a local decrease in the ER luminal Ca(2+) concentration ([Ca(2+)]ER) in dendrites and spines. In this experimental system, the recovery of [Ca(2+)]ER takes several seconds, and recovery half-time depends on the extent of ER Ca(2+) depletion. By combining imaging analysis and numerical simulation, we show that the intraluminal diffusion of Ca(2+), rather than Ca(2+) reuptake, is the dominant mechanism for the replenishment of the local [Ca(2+)]ER depletion immediately following the stimulation. In spines, the ER filled almost simultaneously with parent dendrites, suggesting that the ER within the spine neck does not represent a significant barrier to Ca(2+) diffusion. Furthermore, we found that repetitive climbing fiber stimulation, which induces cytosolic Ca(2+) spikes in PCs, cumulatively increased [Ca(2+)]ER. These results indicate that the neuronal ER functions both as an intracellular tunnel to redistribute stored Ca(2+) within the neurons, and as a leaky integrator of Ca(2+) spike-inducing synaptic inputs.

Ca2+ Diffusion through Endoplasmic Reticulum Supports Elevated Intraterminal Ca2+ Levels Needed to Sustain Synaptic Release from Rods in Darkness

In addition to vesicle release at synaptic ribbons, rod photoreceptors are capable of substantial slow release at non-ribbon release sites triggered by Ca(2+)-induced Ca(2+) release (CICR) from intracellular stores. To maintain CICR as rods remain depolarized in darkness, we hypothesized that Ca(2+) released into the cytoplasm from terminal endoplasmic reticulum (ER) can be replenished continuously by ions diffusing within the ER from the soma. We measured [Ca(2+)] changes in cytoplasm and ER of rods from Ambystoma tigrinum retina using various dyes. ER [Ca(2+)] changes were measured by loading ER with fluo-5N and then washing dye from the cytoplasm with a dye-free patch pipette solution. Small dye molecules diffused within ER between soma and terminal showing a single continuous ER compartment. Depolarization of rods to -40 mV depleted Ca(2+) from terminal ER, followed by a decline in somatic ER [Ca(2+)]. Local activation of ryanodine receptors in terminals with a spatially confined puff of ryanodine caused a decline in terminal ER [Ca(2+)], followed by a secondary decrease in somatic ER. Localized photolytic uncaging of Ca(2+) from o-nitrophenyl-EGTA in somatic ER caused an abrupt Ca(2+) increase in somatic ER, followed by a slower Ca(2+) increase in terminal ER. These data suggest that, during maintained depolarization, a soma-to-terminal [Ca(2+)] gradient develops within the ER that promotes diffusion of Ca(2+) ions to resupply intraterminal ER Ca(2+) stores and thus sustain CICR-mediated synaptic release. The ability of Ca(2+) to move freely through the ER may also promote bidirectional communication of Ca(2+) changes between soma and terminal.

SIGNIFICANCE STATEMENT

Vertebrate rod and cone photoreceptors both release vesicles at synaptic ribbons, but rods also exhibit substantial slow release at non-ribbon sites triggered by Ca(2+)-induced Ca(2+) release (CICR). Blocking CICR inhibits >50% of release from rods in darkness. How do rods maintain sufficiently high [Ca(2+)] in terminal endoplasmic reticulum (ER) to support sustained CICR-driven synaptic transmission? We show that maintained depolarization creates a [Ca(2+)] gradient within the rod ER lumen that promotes soma-to-terminal diffusion of Ca(2+) to replenish intraterminal ER stores. This mechanism allows CICR-triggered synaptic release to be sustained indefinitely while rods remain depolarized in darkness. Free diffusion of Ca(2+) within the ER may also communicate synaptic Ca(2+) changes back to the soma to influence other critical cell processes.

Na⁺ ions as spatial intracellular messengers for co-ordinating Ca²⁺ signals during pH heterogeneity in cardiomyocytes

Aims

Contraction of the heart is regulated by electrically evoked Ca²⁺ transients (CaTs). H⁺ ions, the end products of metabolism, modulate CaTs through direct interactions with Ca²⁺-handling proteins and via Na⁺-mediated coupling between acid-extruding proteins (e.g. Na⁺/H⁺ exchange, NHE1) and Na⁺/Ca²⁺ exchange. Restricted H⁺ diffusivity in cytoplasm predisposes pH-sensitive Ca²⁺ signalling to becoming non-uniform, but the involvement of readily diffusible intracellular Na⁺ ions may provide a means for combatting this.

Methods and results

CaTs were imaged in fluo3-loaded rat ventricular myocytes paced at 2 Hz. Cytoplasmic [Na⁺] ([Na⁺]_i) was imaged using SBFI. Intracellular acidification by acetate exposure raised diastolic and systolic [Ca²⁺] (also observed with acid-loading by ammonium prepulse or CO₂ exposure). The systolic [Ca²⁺] response correlated with a rise in [Na⁺]_i and sarcoplasmic reticulum Ca²⁺ load, and was blocked by the NHE1 inhibitor cariporide (CO₂/HCO₃⁻-free media). Exposure of one half of a myocyte to acetate using dual microperfusion (CO₂/HCO₃⁻-free media) raised diastolic [Ca²⁺] locally in the acidified region. Systolic [Ca²⁺] and CaT amplitude increased more uniformly along the length of the cell, but only when NHE1 was functional. Cytoplasmic Na⁺ diffusivity (D_Na) was measured in quiescent cells, with strophanthidin present to inhibit the Na⁺/K⁺ pump. With regional acetate exposure to activate a local NHE-driven Na⁺-influx, D_Na was found to be sufficiently fast (680 µm²/s) for transmitting the pH–systolic Ca²⁺ interaction over long distances.

Conclusions

Na⁺ ions are rapidly diffusible messengers that expand the spatial scale of cytoplasmic pH–CaT interactions, helping to co-ordinate global Ca²⁺ signalling during conditions of intracellular pH non-uniformity.

Spatially defined InsP3-mediated signaling in embryonic stem cell-derived cardiomyocytes

The functional role of inositol 1,4,5-trisphosphate (InsP3) signaling in cardiomyocytes is not entirely understood but it was linked to an increased propensity for triggered activity. The aim of this study was to determine how InsP3 receptors can translate Ca(2+) release into a depolarization of the plasma membrane and consequently arrhythmic activity. We used embryonic stem cell-derived cardiomyocytes (ESdCs) as a model system since their spontaneous electrical activity depends on InsP3-mediated Ca(2+) release. [InsP3]i was monitored with the FRET-based InsP3-biosensor FIRE-1 (Fluorescent InsP3 Responsive Element) and heterogeneity in sub-cellular [InsP3]i was achieved by targeted expression of FIRE-1 in the nucleus (FIRE-1nuc) or expression of InsP3 5-phosphatase (m43) localized to the plasma membrane. Spontaneous activity of ESdCs was monitored simultaneously as cytosolic Ca(2+) transients (Fluo-4/AM) and action potentials (current clamp). During diastole, the diastolic depolarization was paralleled by an increase of [Ca(2+)]i and spontaneous activity was modulated by [InsP3]i. A 3.7% and 1.7% increase of FIRE-1 FRET ratio and 3.0 and 1.5 fold increase in beating frequency was recorded upon stimulation with endothelin-1 (ET-1, 100 nmol/L) or phenylephrine (PE, 10 µmol/L), respectively. Buffering of InsP3 by FIRE-1nuc had no effect on the basal frequency while attenuation of InsP3 signaling throughout the cell (FIRE-1), or at the plasma membrane (m43) resulted in a 53.7% and 54.0% decrease in beating frequency. In m43 expressing cells the response to ET-1 was completely suppressed. Ca(2+) released from InsP3Rs is more effective than Ca(2+) released from RyRs to enhance INCX. The results support the hypothesis that in ESdCs InsP3Rs form a functional signaling domain with NCX that translates Ca(2+) release efficiently into a depolarization of the membrane potential.

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.

Calcium movements inside the sarcoplasmic reticulum of cardiac myocytes

Sarcoplasmic reticulum (SR) Ca content ([Ca]SRT) is critical to both normal cardiac function and electrophysiology, and changes associated with pathology contribute to systolic and diastolic dysfunction and arrhythmias. The intra-SR free [Ca] ([Ca]SR) dictates the [Ca]SRT, the driving force for Ca release and regulates release channel gating. We discuss measurement of [Ca]SR and [Ca]SRT, how [Ca]SR regulates activation and termination of release, and how Ca diffuses within the SR and influences SR Ca release during excitation-contraction coupling, Ca sparks and Cac waves. The entire SR network is connected and its lumen is also continuous with the nuclear envelope. Rapid Ca diffusion within the SR could stabilize and balance local [Ca]SR within the myocyte, but restrictions to diffusion can create spatial inhomogeneities. Experimental measurements and mathematical models of [Ca]SR to date have greatly enriched our understanding of these [Ca]SR dynamics, but controversies exist and may stimulate new measurements and analysis.

Extraction of sub-microscopic Ca fluxes from blurred and noisy fluorescent indicator images with a detailed model fitting approach

The release of Ca from intracellular stores is key to cardiac muscle function; however, the molecular control of intracellular Ca release remains unclear. Depletion of the intracellular Ca store (sarcoplasmic reticulum, SR) may play an important role, but the ability to measure local SR Ca with fluorescent Ca indicators is limited by the microscope optical resolution and properties of the indicator. This leads to an uncertain degree of spatio-temporal blurring, which is not easily corrected (by deconvolution methods) due to the low signal-to-noise ratio of the recorded signals. In this study, a 3D computer model was constructed to calculate local Ca fluxes and consequent dye signals, which were then blurred by a measured microscope point spread function. Parameter fitting was employed to adjust a release basis function until the model output fitted recorded (2D) Ca spark data. This ‘forward method’ allowed us to obtain estimates of the time-course of Ca release flux and depletion within the sub-microscopic local SR associated with a number of Ca sparks. While variability in focal position relative to Ca spark sites causes more out-of-focus events to have smaller calculated fluxes (and less SR depletion), the average SR depletion was to 20±10% (s.d.) of the resting level. This focus problem implies that the actual SR depletion is likely to be larger and the five largest depletions analyzed were to 8±6% of the resting level. This profound depletion limits SR release flux during a Ca spark, which peaked at 8±3 pA and declined with a half time of 7±2 ms. By comparison, RyR open probability declined more slowly, suggesting release termination is dominated by neither SR Ca depletion nor intrinsic RyR gating, but results from an interaction of these processes.

Calcium levels inside myocytes regulate the heart's force of contraction. Calcium is released from the primary intracellular store called the sarcoplasmic reticulum. Calcium release was directly observed as ‘calcium sparks’ using fluorescent calcium indicators inside the cell. More recently, calcium levels inside the store have been measured as calcium ‘blinks’. These suggest that some depletion of store calcium occurs during cell excitation; however, the actual extent of depletion is made uncertain by the complex sarcoplasmic reticulum shape, dye saturation and optical properties of the microscope. While previous studies have assumed idealized microscope properties, we measured microscope blurring and applied it to a computer model of calcium movements inside the cell. In this model, calcium release was adjusted to match the simulated blurred calcium signals to experimental results. The calculations show that the depth of local sarcoplasmic reticulum calcium depletion is much greater than inferred from calcium blinks and that the time-course of calcium release is affected by this depletion. An estimate for the time-course of gating of the ion channels that regulate calcium release inside the cell was also calculated. We suggest that the time-course of SR Ca release arises from a complex interaction of Ca depletion and channel gating.

Alterations in T-tubule and dyad structure in heart disease: challenges and opportunities for computational analyses

Compelling recent experimental results make clear that sub-cellular structures are altered in ventricular myocytes during the development of heart failure, in both human samples and diverse experimental models. These alterations can include, but are not limited to, changes in the clusters of sarcoplasmic reticulum (SR) Ca(2+)-release channels, ryanodine receptors, and changes in the average distance between the cell membrane and ryanodine receptor clusters. In this review, we discuss the potential consequences of these structural alterations on the triggering of SR Ca(2+) release during excitation-contraction coupling. In particular, we describe how mathematical models of local SR Ca(2+) release can be used to predict functional changes resulting from diverse modifications that occur in disease states. We review recent studies that have used simulations to understand the consequences of sub-cellular structural changes, and we discuss modifications that will allow for future modelling studies to address unresolved questions. We conclude with a discussion of improvements in both experimental and mathematical modelling techniques that will be required to provide a stronger quantitative understanding of the functional consequences of changes in sub-cellular structure in heart disease.

Store-dependent deactivation: cooling the chain-reaction of myocardial calcium signaling

In heart cells, Ca(2+) released from the internal storage unit, the sarcoplasmic reticulum (SR) through ryanodine receptor (RyR2) channels is the predominant determinant of cardiac contractility. Evidence obtained in recent years suggests that SR Ca(2+) release is tightly regulated not only by cytosolic Ca(2+) but also by intra-store Ca(2+) concentration. Specifically, Ca(2+)-induced Ca(2+) release (CICR) that relies on auto-catalytic action of Ca(2+) at the cytosolic side of RyR2s is precisely balanced and counteracted by RyR2 deactivation dependent on a reciprocal decrease of Ca(2+) at the luminal side of RyR2s. Dysregulation of this inherently unstable Ca(2+) signaling is considered to be an underlying cause of triggered arrhythmias, and is associated with genetic and acquired forms of sudden cardiac death. In this article, we present an overview of recent advances in our understanding of the regulatory role luminal Ca(2+) plays in Ca(2+) handling, with a particular emphasis on the role of Ca(2+)release refractoriness in aberrant Ca(2+) release.

β-Adrenergic stimulation increases the intra-sarcoplasmic reticulum Ca2+ threshold for Ca2+ wave generation

β-Adrenergic signalling induces positive inotropic effects on the heart that associate with pro-arrhythmic spontaneous Ca(2+) waves. A threshold level of sarcoplasmic reticulum (SR) Ca(2+) ([Ca(2+)](SR)) is necessary to trigger Ca(2+) waves, and whether the increased incidence of Ca(2+) waves during β-adrenergic stimulation is due to an alteration in this threshold remains controversial. Using the low-affinity Ca(2+) indicator fluo-5N entrapped within the SR of rabbit ventricular myocytes, we addressed this controversy by directly monitoring [Ca(2+)](SR) and Ca(2+) waves during β-adrenergic stimulation. Electrical pacing in elevated extracellular Ca(2+) ([Ca(2+)](o) = 7 mM) was used to increase [Ca(2+)](SR) to the threshold where Ca(2+) waves were consistently observed. The β-adrenergic agonist isoproterenol (ISO; 1 μM) increased [Ca(2+)](SR) well above the control threshold and consistently triggered Ca(2+) waves. However, when [Ca(2+)](SR) was subsequently lowered in the presence of ISO (by lowering [Ca(2+)](o) to 1 mM and partially inhibiting sarcoplasmic/endoplasmic reticulum calcium ATPase with cyclopiazonic acid or thapsigargin), Ca(2+) waves ceased to occur at a [Ca(2+)](SR) that was higher than the control threshold. Furthermore, for a set [Ca(2+)](SR) level the refractoriness of wave occurrence (Ca(2+) wave latency) was prolonged during β-adrenergic stimulation, and was highly dependent on the extent that [Ca](SR) exceeded the wave threshold. These data show that acute β-adrenergic stimulation increases the [Ca(2+)](SR) threshold for Ca(2+) waves, and therefore the primary cause of Ca(2+) waves is the robust increase in [Ca(2+)](SR) above this higher threshold level. Elevation of the [Ca(2+)](SR) wave threshold and prolongation of wave latency represent potentially protective mechanisms against pro-arrhythmogenic Ca(2+) release during β-adrenergic stimulation.

Instabilities of the resting state in a mathematical model of calcium handling in cardiac myocytes

We analyze a recently published model of calcium handling in cardiac myocytes in order to find conditions for the presence of instabilities in the resting state of the model. Such instabilities can create calcium waves which in turn may be able to initiate cardiac arrhythmias. The model was developed by Swietach, Spitzer and Vaughan-Jones in order to study the effect, on calcium waves, of varying ryanodine receptor (RyR)-permeability, sarco/endoplasmic reticulum calcium ATPase (SERCA) and calcium diffusion. We study the model using the extracellular calcium concentration c(e) and the maximal velocity of the SERCA-pump v(SERCA) as control parameters. In the (c(e),v(SERCA))-domain we derive an explicit function v∗=v∗(c(e)), and we claim that any resting state based on parameters that lie above the curve (i.e. any pair (c(e),v(SERCA)) such that with v(SERCA) > v∗(c(e))) is unstable in the sense that small perturbations will grow and can eventually turn into a calcium wave. And conversely; any pair (c(e),v(SERCA)) below the curve is stable in the sense that small perturbations to the resting state will decay to rest. This claim is supported by analyzing the stability of the system in terms of computing the eigenmodes of the linearized model. Furthermore, the claim is supported by direct simulations based on the non-linear model. Since the curve separating stable from unstable states is given as an explicit function, we can show how stability depends on other parameters of the model.

Dynamic local changes in sarcoplasmic reticulum calcium: physiological and pathophysiological roles

Evidence obtained in recent years indicates that, in cardiac myocytes, release of Ca(2+) from the sarcoplasmic reticulum (SR) is regulated by changes in the concentration of Ca(2+) within the SR. In this review, we summarize recent advances in our understanding of this regulatory role, with a particular emphasis on dynamic and local changes in SR [Ca(2+)]. We focus on five important questions that are to some extent unresolved and controversial. These questions concern: (1) the importance of SR [Ca(2+)] depletion in the termination of Ca(2+) release; (2) the quantitative extent of depletion during local release events such as Ca(2+) sparks; (3) the influence of SR [Ca(2+)] refilling on release refractoriness and the propensity for pathological Ca(2+) release; (4) dynamic changes in SR [Ca(2+)] during propagating Ca(2+) waves; and (5) the speed of Ca(2+) diffusion within the SR. With each issue, we discuss data supporting alternative viewpoints, and we identify fundamental questions that are being actively investigated. We conclude with a discussion of experimental and computational advances that will help to resolve controversies. This article is part of a special issue entitled "Local Signaling in Myocytes."

Dynamic calcium movement inside cardiac sarcoplasmic reticulum during release

RATIONALE

Intra-sarcoplasmic reticulum (SR) free [Ca] ([Ca](SR)) provides the driving force for SR Ca release and is a key regulator of SR Ca release channel gating during normal SR Ca release or arrhythmogenic spontaneous Ca release events. However, little is known about [Ca](SR) spatiotemporal dynamics.

OBJECTIVE

To directly measure local [Ca](SR) with subsarcomeric spatiotemporal resolution during both normal global SR Ca release and spontaneous Ca sparks and to evaluate the quantitative implications of spatial [Ca](SR) gradients.

METHODS AND RESULTS

Intact and permeabilized rabbit ventricular myocytes were subjected to direct simultaneous measurement of cytosolic [Ca] and [Ca](SR) and FRAP (fluorescence recovery after photobleach). We found no detectable [Ca](SR) gradients between SR release sites (junctional SR) and Ca uptake sites (free SR) during normal global Ca release, clear spatiotemporal [Ca](SR) gradients during isolated Ca blinks, faster intra-SR diffusion in the longitudinal versus transverse direction, 3- to 4-fold slower diffusion of fluorophores in the SR than in cytosol, and that intra-SR Ca diffusion varies locally, dependent on local SR connectivity. A computational model clarified why spatiotemporal gradients are more detectable in isolated local releases versus global releases and provides a quantitative framework for understanding intra-SR Ca diffusion.

CONCLUSIONS

Intra-SR Ca diffusion is rapid, limiting spatial [Ca](SR) gradients during excitation-contraction coupling. Spatiotemporal [Ca](SR) gradients are apparent during Ca sparks, and these observations constrain models of dynamic Ca movement inside the SR. This has important implications for myocyte SR Ca handling, synchrony, and potentially arrhythmogenic spontaneous contraction.

Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca²+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca²+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca²+ dynamics: 1) the biphasic increment during the upstroke of the Ca²+ transient resulting from the delay between the peripheral and central SR Ca²+ release, and 2) the relative contribution of SL Ca²+ current and SR Ca²+ release to the Ca²+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca²+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca²+ release sites define the interface between Ca²+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca²+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca²+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca²+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.

Predicting local SR Ca(2+) dynamics during Ca(2+) 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(2+)) wave propagation, and 2), speed of Ca(2+) diffusion within the SR. Ca(2+) waves are initiated when a spontaneous local SR Ca(2+) release event triggers additional release from neighboring clusters of SR release channels (ryanodine receptors (RyRs)). A lack of consensus regarding the effective Ca(2+) diffusion constant in the SR (D(Ca,SR)) severely complicates our understanding of whether dynamic local changes in SR [Ca(2+)] can influence wave propagation. To address this problem, we have implemented a computational model of cytosolic and SR [Ca(2+)] during Ca(2+) waves. Simulations have investigated how dynamic local changes in SR [Ca(2+)] are influenced by 1), D(Ca,SR); 2), the distance between RyR clusters; 3), partial inhibition or stimulation of SR Ca(2+) pumps; 4), SR Ca(2+) pump dependence on cytosolic [Ca(2+)]; 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(2+)] in nearby regions. Specifically, our results show that local increases in SR [Ca(2+)] ahead of the wave can potentially facilitate Ca(2+) wave propagation, but only if SR diffusion is relatively slow. These simulations help to delineate what changes in [Ca(2+)] are possible during SR Ca(2+)release, and they broaden our understanding of the regulatory role played by dynamic changes in [Ca(2+)](SR).

Calsequestrin mediates changes in spontaneous calcium release profiles

Calsequestrin (CSQ) is the primary calcium buffer within the sarcoplasmic reticulum (SR) of cardiac cells. It has also been identified as a regulator of Ryanodine receptor (RyR) calcium release channels by serving as a SR luminal sensor. When calsequestrin is free and unbound to calcium, it can bind to RyR and desensitize the channel from cytoplasmic calcium activation. In this paper, we study the role of CSQ as a buffer and RyR luminal sensor using a mechanistic model of RyR-CSQ interaction. By using various asymptotic approximations and mean first exit time calculation, we derive a minimal model of a calcium release unit which includes CSQ dependence. Using this model, we then analyze the effect of changing CSQ expression on the calcium release profile and the rate of spontaneous calcium release. We show that because of its buffering capability, increasing CSQ increases the spark duration and size. However, because of luminal sensing effects, increasing CSQ depresses the basal spark rate and increases the critical SR level for calcium release termination. Finally, we show that with increased bulk cytoplasmic calcium concentration, the CRU model exhibits deterministic oscillations.

Models of cardiac excitation-contraction coupling in ventricular myocytes

Mathematical and computational modeling of cardiac excitation-contraction coupling has produced considerable insights into how the heart muscle contracts. With the increase in biophysical and physiological data available, the modeling has become more sophisticated with investigations spanning in scale from molecular components to whole cells. These modeling efforts have provided insight into cardiac excitation-contraction coupling that advanced and complemented experimental studies. One goal is to extend these detailed cellular models to model the whole heart. While this has been done with mechanical and electrophysiological models, the complexity and fast time course of calcium dynamics have made inclusion of detailed calcium dynamics in whole heart models impractical. Novel methods such as the probability density approach and moment closure technique which increase computational efficiency might make this tractable.

Excitation-contraction coupling in human heart failure examined by action potential clamp in rat cardiac myocytes

The effect of the loss of the notch in the human action potential (AP) during heart failure was examined by voltage clamping rat ventricular myocytes with human APs and recording intracellular Ca(2+) with fluorescent dyes. Loss of the notch resulted in about a 50% reduction in the initial phase of the Ca(2+) transient due to reduced ability of the L-type Ca(2+) channel to trigger release. The failing human AP increased non-uniformity of cytosolic Ca(2+), with some cellular regions failing to elicit Ca(2+)-induced Ca(2+) release from the sarcoplasmic reticulum. In addition, there was an increase in the occurrence of late Ca(2+) sparks. Monte-Carlo simulations of spark activation by L-type Ca(2+) current supported the idea that the decreased synchrony of Ca(2+) spark production associated with the loss of the notch could be explained by reduced Ca(2+) influx from open Ca(2+) channels. We conclude that the notch of the AP is critical for efficient and synchronous EC coupling and that the loss of the notch will reduce the SR Ca(2+) release in heart failure, without changes in (for example) SR Ca(2+)-ATPase uptake.

Cardiac cell modelling: observations from the heart of the cardiac physiome project

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field.

Modeling calcium waves in cardiac myocytes: importance of calcium diffusion

Under certain conditions, cardiac myocytes engage in a mode of calcium signaling in which calcium release from the sarcoplasmic reticulum (SR) to myoplasm occurs in self-propagating succession along the length of the cell. This event is called a calcium wave and is fundamentally a diffusion-reaction phenomenon. We present a simple, continuum mathematical model that simulates calcium waves. The framework features calcium diffusion within the SR and myoplasm, and dual modulation of ryanodine receptor (RyR) release channels by myoplasmic and SR calcium. The model is used to illustrate the effect of varying RyR permeability, sarco/endoplasmic reticulum Ca2+ ATPase (SERCA) activity and calcium ion mobility in myoplasm and SR on wave velocity. The model successfully reproduces calcium waves using experimentally-derived variables. It also supports the proposal for wave propagation driven by the diffusive spread of myoplasmic calcium, and highlights the importance of SR calcium load on wave propagation.

Ca2+: is there something new for the cardiovascular anesthesiologist?

PURPOSE OF REVIEW

Anesthesiologists are frequently called upon to treat abnormalities of heart rhythm or pumping ability. Intracellular Ca is crucial for normal excitation-contraction coupling in the heart and plays a major role in the sequence of events that starts with an electrical signal generated in the atria and ends with myocardial contraction.

RECENT FINDINGS

From controlled diffusion within the cell to a potential role as a biological clock, intracellular Ca is receiving a great deal of attention. For example, the pacemaking electrical signal is known to originate in the sinoatrial node myocyte, but exactly what role Ca plays is controversial despite the fact that the sinoatrial node was discovered over 100 years ago. Basic mechanisms involved in disease processes such as atrial fibrillation and new interventions for heart rate control are beginning to emerge. New discoveries in ventricular myocytes are also stimulating the development of promising therapeutic interventions to safely increase the pumping ability of the heart.

SUMMARY

As our understanding of cardiac physiology and pharmacology progresses at the subcellular and molecular levels, new therapies will continue to emerge and the practice of anesthesia will benefit greatly.

Termination of cardiac Ca2+ sparks: role of intra-SR [Ca2+], release flux, and intra-SR Ca2+ diffusion

Ca(2+) release from cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) is regulated by dyadic cleft [Ca(2+)] and intra-SR free [Ca(2+)] ([Ca(2+)](SR)). Robust SR Ca(2+) release termination is important for stable excitation-contraction coupling, and partial [Ca(2+)](SR) depletion may contribute to release termination. Here, we investigated the regulation of SR Ca(2+) release termination of spontaneous local SR Ca(2+) release events (Ca(2+) sparks) by [Ca(2+)](SR), release flux, and intra-SR Ca(2+) diffusion. We simultaneously measured Ca(2+) sparks and Ca(2+) blinks (localized elementary [Ca(2+)](SR) depletions) in permeabilized ventricular cardiomyocytes over a wide range of SR Ca(2+) loads and release fluxes. Sparks terminated via a [Ca(2+)](SR)-dependent mechanism at a fixed [Ca(2+)](SR) depletion threshold independent of the initial [Ca(2+)](SR) and release flux. Ca(2+) blink recovery depended mainly on intra-SR Ca(2+) diffusion rather than SR Ca(2+) uptake. Therefore, the large variation in Ca(2+) blink recovery rates at different release sites occurred because of differences in the degree of release site interconnection within the SR network. When SR release flux was greatly reduced, long-lasting release events occurred from well-connected junctions. These junctions could sustain release because local SR Ca(2+) release and [Ca(2+)](SR) refilling reached a balance, preventing [Ca(2+)](SR) from depleting to the termination threshold. Prolonged release events eventually terminated at a steady [Ca(2+)](SR), indicative of a slower, [Ca(2+)](SR)-independent termination mechanism. These results demonstrate that there is high variability in local SR connectivity but that SR Ca(2+) release terminates at a fixed [Ca(2+)](SR) termination threshold. Thus, reliable SR Ca(2+) release termination depends on tight RyR regulation by [Ca(2+)](SR).