Modeling gain and gradedness of Ca2+ release in the functional unit of the cardiac diadic space

A possible path to persistent re-entry waves at the outlet of the left pulmonary vein

Atrial fibrillation (AF) is the most common form of cardiac arrhythmia, often evolving from paroxysmal episodes to persistent stages over an extended timeframe. While various factors contribute to this progression, the precise biophysical mechanisms driving it remain unclear. Here we explore how rapid firing of cardiomyocytes at the outlet of the pulmonary vein of the left atria can create a substrate for a persistent re-entry wave. This is grounded in a recently formulated mathematical model of the regulation of calcium ion channel density by intracellular calcium concentration. According to the model, the number of calcium channels is controlled by the intracellular calcium concentration. In particular, if the concentration increases above a certain target level, the calcium current is weakened to restore the target level of calcium. During rapid pacing, the intracellular calcium concentration of the cardiomyocytes increases leading to a substantial reduction of the calcium current across the membrane of the myocytes, which again reduces the action potential duration. In a spatially resolved cell-based model of the outlet of the pulmonary vein of the left atria, we show that the reduced action potential duration can lead to re-entry. Initiated by rapid pacing, often stemming from paroxysmal AF episodes lasting several days, the reduction in calcium current is a critical factor. Our findings illustrate how such episodes can foster a conducive environment for persistent AF through electrical remodeling, characterized by diminished calcium currents. This underscores the importance of promptly addressing early AF episodes to prevent their progression to chronic stages.

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Do calcium channel blockers applied to cardiomyocytes cause increased channel expression resulting in reduced efficacy?

In the initial hours following the application of the calcium channel blocker (CCB) nifedipine to microtissues consisting of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), we observe notable variations in the drug's efficacy. Here, we investigate the possibility that these temporal changes in CCB effects are associated with adaptations in the expression of calcium ion channels in cardiomyocyte membranes. To explore this, we employ a recently developed mathematical model that delineates the regulation of calcium ion channel expression by intracellular calcium concentrations. According to the model, a decline in intracellular calcium levels below a certain target level triggers an upregulation of calcium ion channels. Such an upregulation, if instigated by a CCB, would then counteract the drug's inhibitory effect on calcium currents. We assess this hypothesis using time-dependent measurements of hiPSC-CMs dynamics and by refining an existing mathematical model of myocyte action potentials incorporating the dynamic nature of the number of calcium ion channels. The revised model forecasts that the CCB-induced reduction in intracellular calcium concentrations leads to a subsequent increase in calcium ion channel expression, thereby attenuating the drug's overall efficacy. The data and fit models suggest that dynamic changes in cardiac cells in the presence of CCBs may be explainable by induced changes in protein expression, and that this may lead to challenges in understanding calcium based drug effects on the heart unless timings of applications are carefully considered.

Local Control Model of a Human Ventricular Myocyte: An Exploration of Frequency-Dependent Changes and Calcium Sparks

Calcium (Ca2+) sparks are the elementary events of excitation-contraction coupling, yet they are not explicitly represented in human ventricular myocyte models. A stochastic ventricular cardiomyocyte human model that adapts to intracellular Ca2+ ([Ca2+]i) dynamics, spark regulation, and frequency-dependent changes in the form of locally controlled Ca2+ release was developed. The 20,000 CRUs in this model are composed of 9 individual LCCs and 49 RyRs that function as couplons. The simulated action potential duration at 1 Hz steady-state pacing is ~0.280 s similar to human ventricular cell recordings. Rate-dependence experiments reveal that APD shortening mechanisms are largely contributed by the L-type calcium channel inactivation, RyR open fraction, and [Ca2+]myo concentrations. The dynamic slow-rapid-slow pacing protocol shows that RyR open probability during high pacing frequency (2.5 Hz) switches to an adapted "nonconducting" form of Ca2+-dependent transition state. The predicted force was also observed to be increased in high pacing, but the SR Ca2+ fractional release was lower due to the smaller difference between diastolic and systolic [Ca2+]SR. Restitution analysis through the S1S2 protocol and increased LCC Ca2+-dependent activation rate show that the duration of LCC opening helps modulate its effects on the APD restitution at different diastolic intervals. Ultimately, a longer duration of calcium sparks was observed in relation to the SR Ca2+ loading at high pacing rates. Overall, this study demonstrates the spontaneous Ca2+ release events and ion channel responses throughout various stimuli.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

Critical Requirements for the Initiation of a Cardiac Arrhythmia in Rat Ventricle: How Many Myocytes?

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca²⁺ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca²⁺ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca²⁺ sparks can combine to form Ca²⁺ waves. These propagating elevations of (Ca²⁺)_i can activate an inward Na⁺-Ca²⁺ exchanger current (I_NCX) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.

Cardiac Alternans Occurs through the Synergy of Voltage- and Calcium-Dependent Mechanisms

Cardiac alternans is characterized by alternating weak and strong beats of the heart. This signaling at the cellular level may appear as alternating long and short action potentials (APs) that occur in synchrony with alternating large and small calcium transients, respectively. Previous studies have suggested that alternans manifests itself through either a voltage dependent mechanism based upon action potential restitution or as a calcium dependent mechanism based on refractoriness of calcium release. We use a novel model of cardiac excitation-contraction (EC) coupling in the rat ventricular myocyte that includes 20,000 calcium release units (CRU) each with 49 ryanodine receptors (RyR2s) and 7 L-type calcium channels that are all stochastically gated. The model suggests that at the cellular level in the case of alternans produced by rapid pacing, the mechanism requires a synergy of voltage- and calcium-dependent mechanisms. The rapid pacing reduces AP duration and magnitude reducing the number of L-type calcium channels activating individual CRUs during each AP and thus increases the population of CRUs that can be recruited stochastically. Elevated myoplasmic and sarcoplasmic reticulum (SR) calcium, [Ca2+]myo and [Ca2+]SR respectively, increases ryanodine receptor open probability (Po) according to our model used in this simulation and this increased the probability of activating additional CRUs. A CRU that opens in one beat is less likely to open the subsequent beat due to refractoriness caused by incomplete refilling of the junctional sarcoplasmic reticulum (jSR). Furthermore, the model includes estimates of changes in Na+ fluxes and [Na+]i and thus provides insight into how changes in electrical activity, [Na+]i and sodium-calcium exchanger activity can modulate alternans. The model thus tracks critical elements that can account for rate-dependent changes in [Na+]i and [Ca2+]myo and how they contribute to the generation of Ca2+ signaling alternans in the heart.

Modulations of Cardiac Functions and Pathogenesis by Reactive Oxygen Species and Natural Antioxidants

Homeostasis in the level of reactive oxygen species (ROS) in cardiac myocytes plays a critical role in regulating their physiological functions. Disturbance of balance between generation and removal of ROS is a major cause of cardiac myocyte remodeling, dysfunction, and failure. Cardiac myocytes possess several ROS-producing pathways, such as mitochondrial electron transport chain, NADPH oxidases, and nitric oxide synthases, and have endogenous antioxidation mechanisms. Cardiac Ca²⁺-signaling toolkit proteins, as well as mitochondrial functions, are largely modulated by ROS under physiological and pathological conditions, thereby producing alterations in contraction, membrane conductivity, cell metabolism and cell growth and death. Mechanical stresses under hypertension, post-myocardial infarction, heart failure, and valve diseases are the main causes for stress-induced cardiac remodeling and functional failure, which are associated with ROS-induced pathogenesis. Experimental evidence demonstrates that many cardioprotective natural antioxidants, enriched in foods or herbs, exert beneficial effects on cardiac functions (Ca²⁺ signal, contractility and rhythm), myocytes remodeling, inflammation and death in pathological hearts. The review may provide knowledge and insight into the modulation of cardiac pathogenesis by ROS and natural antioxidants.

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.

Improved Computational Identification of Drug Response Using Optical Measurements of Human Stem Cell Derived Cardiomyocytes in Microphysiological Systems

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) hold great potential for drug screening applications. However, their usefulness is limited by the relative immaturity of the cells' electrophysiological properties as compared to native cardiomyocytes in the adult human heart. In this work, we extend and improve on methodology to address this limitation, building on previously introduced computational procedures which predict drug effects for adult cells based on changes in optical measurements of action potentials and Ca2+ transients made in stem cell derived cardiac microtissues. This methodology quantifies ion channel changes through the inversion of data into a mathematical model, and maps this response to an adult phenotype through the assumption of functional invariance of fundamental intracellular and membrane channels during maturation. Here, we utilize an updated action potential model to represent both hiPSC-CMs and adult cardiomyocytes, apply an IC50-based model of dose-dependent drug effects, and introduce a continuation-based optimization algorithm for analysis of dose escalation measurements using five drugs with known effects. The improved methodology can identify drug induced changes more efficiently, and quantitate important metrics such as IC50 in line with published values. Consequently, the updated methodology is a step towards employing computational procedures to elucidate drug effects in adult cardiomyocytes for new drugs using stem cell-derived experimental tissues.

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

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

Myocyte-fibroblast communication in cardiac fibrosis and arrhythmias: Mechanisms and model systems

Development of cardiac fibrosis and arrhythmias is controlled by the activity of and communication between cardiomyocytes and fibroblasts in the heart. Myocyte-fibroblast interactions occur via both direct and indirect means including paracrine mediators, extracellular matrix interactions, electrical modulators, mechanical junctions, and membrane nanotubes. In the diseased heart, cardiomyocyte and fibroblast ratios and activity, and thus myocyte-fibroblast interactions, change and are thought to contribute to the course of disease including development of fibrosis and arrhythmogenic activity. Fibroblasts have a developing role in modulating cardiomyocyte electrical and hypertrophic activity, however gaps in knowledge regarding these interactions still exist. Research in this field has necessitated the development of unique approaches to isolate and control myocyte-fibroblast interactions. Numerous methods for 2D and 3D co-culture systems have been developed, while a growing part of this field is in the use of better tools for in vivo systems including cardiomyocyte and fibroblast specific Cre mouse lines for cell type specific genetic ablation. This review will focus on (i) mechanisms of myocyte-fibroblast communication and their effects on disease features such as cardiac fibrosis and arrhythmias as well as (ii) methods being used and currently developed in this field.

Ryanodine receptor sensitivity governs the stability and synchrony of local calcium release during cardiac excitation-contraction coupling

Calcium-induced calcium release is the principal mechanism that triggers the cell-wide [Ca(2+)]i transient that activates muscle contraction during cardiac excitation-contraction coupling (ECC). Here, we characterize this process in mouse cardiac myocytes with a novel mathematical action potential (AP) model that incorporates realistic stochastic gating of voltage-dependent L-type calcium (Ca(2+)) channels (LCCs) and sarcoplasmic reticulum (SR) Ca(2+) release channels (the ryanodine receptors, RyR2s). Depolarization of the sarcolemma during an AP stochastically activates the LCCs elevating subspace [Ca(2+)] within each of the cell's 20,000 independent calcium release units (CRUs) to trigger local RyR2 opening and initiate Ca(2+) sparks, the fundamental unit of triggered Ca(2+) release. Synchronization of Ca(2+) sparks during systole depends on the nearly uniform cellular activation of LCCs and the likelihood of local LCC openings triggering local Ca(2+) sparks (ECC fidelity). The detailed design and true SR Ca(2+) pump/leak balance displayed by our model permits investigation of ECC fidelity and Ca(2+) spark fidelity, the balance between visible (Ca(2+) spark) and invisible (Ca(2+) quark/sub-spark) SR Ca(2+) release events. Excess SR Ca(2+) leak is examined as a disease mechanism in the context of "catecholaminergic polymorphic ventricular tachycardia (CPVT)", a Ca(2+)-dependent arrhythmia. We find that that RyR2s (and therefore Ca(2+) sparks) are relatively insensitive to LCC openings across a wide range of membrane potentials; and that key differences exist between Ca(2+) sparks evoked during quiescence, diastole, and systole. The enhanced RyR2 [Ca(2+)]i sensitivity during CPVT leads to increased Ca(2+) spark fidelity resulting in asynchronous systolic Ca(2+) spark activity. It also produces increased diastolic SR Ca(2+) leak with some prolonged Ca(2+) sparks that at times become "metastable" and fail to efficiently terminate. There is a huge margin of safety for stable Ca(2+) handling within the cell and this novel mechanistic model provides insight into the molecular signaling characteristics that help maintain overall Ca(2+) stability even under the conditions of high SR Ca(2+) leak during CPVT. Finally, this model should provide tools for investigators to examine normal and pathological Ca(2+) signaling characteristics in the heart.

Stretch-induced increase in cardiac contractility is independent of myocyte Ca2+ while block of stretch channels by streptomycin improves contractility after ischemic stunning

Stretching the cardiac left ventricle (LV) enhances contractility but its effect on myoplasmic [Ca²⁺] is controversial. We measured LV pressure (LVP) and [Ca²⁺] as a function of intra-LV stretch in guinea pig intact hearts before and after 15 min global stunning ± perfusion with streptomycin (STM), a stretch-activated channel blocker. LV wall [Ca²⁺] was measured by indo-1 fluorescence and LVP by a saline-filled latex balloon inflated in 50 μL steps to stretch the LV. We implemented a mathematical model to interpret cross-bridge dynamics and myofilament Ca²⁺ responsiveness from the instantaneous relationship between [Ca²⁺] and LVP ± stretching. We found that: (1) stretch enhanced LVP but not [Ca²⁺] before and after stunning in either control (CON) and STM groups, (2) after stunning [Ca²⁺] increased in both groups although higher in STM versus CON (56% vs. 39%), (3) STM-enhanced LVP after stunning compared to CON (98% vs. 76% of prestunning values), and (4) stretch-induced effects on LVP were independent of [Ca²⁺] before or after stunning in both groups. Mathematical modeling suggested: (1) cooperativity in cross-bridge kinetics and myofilament Ca²⁺ handling is reduced after stunning in the unstretched heart, (2) stunning results in depressed myofilament Ca²⁺ sensitivity in the presence of attached cross-bridges regardless of stretch, and (3) the initial mechanism responsible for increased contractility during stretch may be enhanced formation of cross-bridges. Thus stretch-induced enhancement of contractility is not due to increased [Ca²⁺], whereas enhanced contractility after stunning in STM versus CON hearts results from improved Ca²⁺ handling and/or enhanced actinomyosin cross-bridge cycling.

Cytosolic Ca²⁺ buffering determines the intra-SR Ca²⁺ concentration at which cardiac Ca²⁺ sparks terminate

Single ryanodine receptor (RyR) Ca(2+) flux amplitude (i(Ca-RyR)) decreases as intra-sarcoplasmic reticulum (SR) Ca(2+) levels fall during a cardiac Ca(2+) spark. Since i(Ca-RyR) drives the inter-RyR Ca(2+)-induced Ca(2+) release (CICR) that underlies the spark, decreasing i(Ca-RyR) may contribute to spark termination because RyRs that spontaneously close may stay closed. To test this possibility, we simultaneously measured local cytosolic and intra-SR ([Ca(2+)]cyto and [Ca(2+)]SR) during Ca(2+) sparks in permeabilized rabbit ventricular myocytes. Local cytosolic or intra-SR Ca(2+) dynamics were manipulated using Ca(2+) buffers. Buffer manipulations applied in cells had no effect on individual RyR channels reconstituted in planar lipid bilayers. Presence of a fast cytosolic Ca(2+) buffer (BAPTA) significantly suppressed Ca(2+) spark activity and sparks terminated earlier at a higher than usual [Ca(2+)]SR level (∼80% vs. ∼62%). When cytosolic Ca(2+) buffer power was reduced (i.e. cytosolic EGTA level decreased), sparks terminated later and at a lower than usual [Ca(2+)]SR level (∼45% vs. ∼62%). When intra-SR Ca(2+) buffer power was increased, sparks also terminated later and at a lower than usual [Ca(2+)]SR (∼48% vs. ∼62%). These results suggest that cytosolic local control of inter-RyR CICR by i(Ca-RyR) plays a substantial role during the spark termination process. Thus, alterations in local cytosolic Ca(2+) handling dynamics in the dyadic cleft (Ca(2+) buffering, extrusion, etc.) likely influence Ca(2+) spark termination.

Calcium homeostasis in a local/global whole cell model of permeabilized ventricular myocytes with a Langevin description of stochastic calcium release

Population density approaches to modeling local control of Ca2+-induced Ca2+ release in cardiac myocytes can be used to construct minimal whole cell models that accurately represent heterogeneous local Ca2+ signals. Unfortunately, the computational complexity of such “local/global” whole cell models scales with the number of Ca2+ release unit (CaRU) states, which is a rapidly increasing function of the number of ryanodine receptors (RyRs) per CaRU. Here we present an alternative approach based on a Langevin description of the collective gating of RyRs coupled by local Ca2+ concentration ([Ca2+]). The computational efficiency of this approach no longer depends on the number of RyRs per CaRU. When the RyR model is minimal, Langevin equations may be replaced by a single Fokker-Planck equation, yielding an extremely compact and efficient local/global whole cell model that reproduces and helps interpret recent experiments that investigate Ca2+ homeostasis in permeabilized ventricular myocytes. Our calculations show that elevated myoplasmic [Ca2+] promotes elevated network sarcoplasmic reticulum (SR) [Ca2+] via SR Ca2+-ATPase-mediated Ca2+ uptake. However, elevated myoplasmic [Ca2+] may also activate RyRs and promote stochastic SR Ca2+ release, which can in turn decrease SR [Ca2+]. Increasing myoplasmic [Ca2+] results in an exponential increase in spark-mediated release and a linear increase in nonspark-mediated release, consistent with recent experiments. The model exhibits two steady-state release fluxes for the same network SR [Ca2+] depending on whether myoplasmic [Ca2+] is low or high. In the later case, spontaneous release decreases SR [Ca2+] in a manner that maintains robust Ca2+ sparks.

Mathematical modeling of physiological systems: an essential tool for discovery

Mathematical models are invaluable tools for understanding the relationships between components of a complex system. In the biological context, mathematical models help us understand the complex web of interrelations between various components (DNA, proteins, enzymes, signaling molecules etc.) in a biological system, gain better understanding of the system as a whole, and in turn predict its behavior in an altered state (e.g. disease). Mathematical modeling has enhanced our understanding of multiple complex biological processes like enzyme kinetics, metabolic networks, signal transduction pathways, gene regulatory networks, and electrophysiology. With recent advances in high throughput data generation methods, computational techniques and mathematical modeling have become even more central to the study of biological systems. In this review, we provide a brief history and highlight some of the important applications of modeling in biological systems with an emphasis on the study of excitable cells. We conclude with a discussion about opportunities and challenges for mathematical modeling going forward. In a larger sense, the review is designed to help answer a simple but important question that theoreticians frequently face from interested but skeptical colleagues on the experimental side: "What is the value of a model?"

Parameter sensitivity analysis of stochastic models provides insights into cardiac calcium sparks

We present a parameter sensitivity analysis method that is appropriate for stochastic models, and we demonstrate how this analysis generates experimentally testable predictions about the factors that influence local Ca(2+) release in heart cells. The method involves randomly varying all parameters, running a single simulation with each set of parameters, running simulations with hundreds of model variants, then statistically relating the parameters to the simulation results using regression methods. We tested this method on a stochastic model, containing 18 parameters, of the cardiac Ca(2+) spark. Results show that multivariable linear regression can successfully relate parameters to continuous model outputs such as Ca(2+) spark amplitude and duration, and multivariable logistic regression can provide insight into how parameters affect Ca(2+) spark triggering (a probabilistic process that is all-or-none in a single simulation). Benchmark studies demonstrate that this method is less computationally intensive than standard methods by a factor of 16. Importantly, predictions were tested experimentally by measuring Ca(2+) sparks in mice with knockout of the sarcoplasmic reticulum protein triadin. These mice exhibit multiple changes in Ca(2+) release unit structures, and the regression model both accurately predicts changes in Ca(2+) spark amplitude (30% decrease in model, 29% decrease in experiments) and provides an intuitive and quantitative understanding of how much each alteration contributes to the result. This approach is therefore an effective, efficient, and predictive method for analyzing stochastic mathematical models to gain biological insight.

Mechanism of calsequestrin regulation of single cardiac ryanodine receptor in normal and pathological conditions

Release of Ca²⁺ from the sarcoplasmic reticulum (SR) drives contractile function of cardiac myocytes. Luminal Ca²⁺ regulation of SR Ca²⁺ release is fundamental not only in physiology but also in physiopathology because abnormal luminal Ca²⁺ regulation is known to lead to arrhythmias, catecholaminergic polymorphic ventricular tachycardia (CPVT), and/or sudden cardiac arrest, as inferred from animal model studies. Luminal Ca²⁺ regulates ryanodine receptor (RyR)2-mediated SR Ca²⁺ release through mechanisms localized inside the SR; one of these involves luminal Ca²⁺ interacting with calsequestrin (CASQ), triadin, and/or junctin to regulate RyR2 function.

CASQ2-RyR2 regulation was examined at the single RyR2 channel level. Single RyR2s were incorporated into planar lipid bilayers by the fusion of native SR vesicles isolated from either wild-type (WT), CASQ2 knockout (KO), or R33Q-CASQ2 knock-in (KI) mice. KO and KI mice have CPVT-like phenotypes. We show that CASQ2(WT) action on RyR2 function (either activation or inhibition) was strongly influenced by the presence of cytosolic MgATP. Function of the reconstituted CASQ2(WT)–RyR2 complex was unaffected by changes in luminal free [Ca²⁺] (from 0.1 to 1 mM). The inhibition exerted by CASQ2(WT) association with the RyR2 determined a reduction in cytosolic Ca²⁺ activation sensitivity. RyR2s from KO mice were significantly more sensitive to cytosolic Ca²⁺ activation and had significantly longer mean open times than RyR2s from WT mice. Sensitivity of RyR2s from KI mice was in between that of RyR2 channels from KO and WT mice. Enhanced cytosolic RyR2 Ca²⁺ sensitivity and longer RyR2 open times likely explain the CPVT-like phenotype of both KO and KI mice.

Pernicious attrition and inter-RyR2 CICR current control in cardiac muscle

In cardiac muscle cells, ryanodine receptor (RyR) mediated Ca(2+) release from the sarcoplasmic reticulum (SR) drives the contractile apparatus. Spontaneous bouts of inter-RyR Ca(2+) induced Ca(2+) release (CICR) generate an elemental unit of SR Ca(2+) release called a spark. Sparks are localized events that terminate soon after they begin. The local control of sparks is not clearly understood. In this article, we review the potential regulatory role that the changing single RyR Ca(2+) current may play. Moreover, we aggregate RyR data into a working scheme of inter-RyR CICR current control of sparks and a potential inter-RyR CICR termination mechanism that we call pernicious attrition.

Decoding myocardial Ca²⁺ signals across multiple spatial scales: a role for sensitivity analysis

Numerous studies have employed mathematical modeling to quantitatively understand release of Ca(2+) from the sarcoplasmic reticulum (SR) in the heart. Models have been used to investigate physiologically important phenomena such as triggering of SR Ca(2+) release by Ca(2+) entry across the cell membrane and spontaneous leak of Ca(2+) from the SR in quiescent heart cells. In this review we summarize studies that have modeled myocardial Ca(2+) at different spatial scales: the sub-cellular level, the cellular level, and the multicellular level. We discuss each category of models from the standpoint of parameter sensitivity analysis, a common simulation procedure that can generate quantitative, comprehensive predictions about how changes in conditions influence model output. We propose that this is a useful perspective for conceptualizing models, in part because a sensitivity analysis requires the investigator to define the relevant parameters and model outputs. This procedure therefore helps to illustrate the capabilities and limitations of each model. We further suggest that in future studies, sensitivity analyses will aid in simplifying complex models and in suggesting experiments to differentiate between competing models built with different assumptions. We conclude with a discussion of unresolved questions that are likely to be addressed over the next several years.

Toward an integrative computational model of the Guinea pig cardiac myocyte

The local control theory of excitation-contraction (EC) coupling asserts that regulation of calcium (Ca²⁺) release occurs at the nanodomain level, where openings of single L-type Ca²⁺ channels (LCCs) trigger openings of small clusters of ryanodine receptors (RyRs) co-localized within the dyad. A consequence of local control is that the whole-cell Ca²⁺ transient is a smooth continuous function of influx of Ca²⁺ through LCCs. While this so-called graded release property has been known for some time, its functional importance to the integrated behavior of the cardiac ventricular myocyte has not been fully appreciated. We previously formulated a biophysically based model, in which LCCs and RyRs interact via a coarse-grained representation of the dyadic space. The model captures key features of local control using a low-dimensional system of ordinary differential equations. Voltage-dependent gain and graded Ca²⁺ release are emergent properties of this model by virtue of the fact that model formulation is closely based on the sub-cellular basis of local control. In this current work, we have incorporated this graded release model into a prior model of guinea pig ventricular myocyte electrophysiology, metabolism, and isometric force production. The resulting integrative model predicts the experimentally observed causal relationship between action potential (AP) shape and timing of Ca²⁺ and force transients, a relationship that is not explained by models lacking the graded release property. Model results suggest that even relatively subtle changes in AP morphology that may result, for example, from remodeling of membrane transporter expression in disease or spatial variation in cell properties, may have major impact on the temporal waveform of Ca²⁺ transients, thus influencing tissue level electromechanical function.

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

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

How the Hodgkin-Huxley equations inspired the Cardiac Physiome Project

Early modelling of cardiac cells (1960-1980) was based on extensions of the Hodgkin-Huxley nerve axon equations with additional channels incorporated, but after 1980 it became clear that processes other than ion channel gating were also critical in generating electrical activity. This article reviews the development of models representing almost all cell types in the heart, many different species, and the software tools that have been created to facilitate the cardiac Physiome Project.

Models of excitation-contraction coupling in cardiac ventricular myocytes

Excitation-contraction coupling describes the processes relating to electrical excitation through force generation and contraction in the heart. It occurs at multiple levels from the whole heart, to single myocytes and down to the sarcomere. A central process that links electrical excitation to contraction is calcium mobilization. Computational models that are well grounded in experimental data have been an effective tool to understand the complex dynamics of the processes involved in excitation-contraction coupling. Presented here is a summary of some computational models that have added to the understanding of the cellular and subcellular mechanisms that control ventricular myocyte calcium dynamics. Models of cardiac ventricular myocytes that have given insight into termination of calcium release and interval-force relations are discussed in this manuscript. Computational modeling of calcium sparks, the elementary events in cardiac excitation-contraction coupling, has given insight into mechanism governing their dynamics and termination as well as their role in excitation-contraction coupling and is described herein.

A model of the physiological basis of a multivariate phenotype that is mediated by Ca(2+) signaling and controlled by ryanodine receptor composition

Calcium-signals occur in a wide variety of tissue types - from skeletal, smooth and cardiac muscle to pancreatic and brain tissues. Ca(2+) signals regulate diverse processes including muscle contraction, hormone secretion, neural communication and gene expression. Together these different tissues and processes form the basis of a multivariate trait. Calcium signals are characterized by Ca(2+) transients, which are sharp increases in Ca(2+) concentration over a short period of time. In this paper we derive and analyze a model of Ca(2+) transients for skeletal muscle, neurons and cardiac tissue based on underlying biophysical principles. Tissue differentiation in our model and in nature comes about by varying the ryanodine receptor (RyR) channel composition of tissues. In vertebrates, there are typically three types of RyR channels (labeled RyR1, RyR2 and RyR3 in mammals and α-RyR, cardiac-RyR and β-RyR in birds, amphibians and fish). Different compositions of these three RyR channels generate different Ca(2+) transient properties. There are four Ca(2+) transient properties that we measure: maximum amplitude, duration, half duration (D(50)) and integrated concentration. In agreement with experimental work, our results find that the addition of RyR3 amplifies Ca(2+) transients in skeletal muscle. An important consequence of shared molecular components between tissue types in a multivariate setting is that the shared components cause individual traits of a multivariate trait to be correlated in function. Here we show how correlations in Ca(2+) transient properties between tissues can be predicted using an underlying biophysical model.

Integrative systems models of cardiac excitation-contraction coupling

Excitation-contraction coupling in the cardiac myocyte is mediated by a number of highly integrated mechanisms of intracellular Ca²(+) transport. The complexity and integrative nature of heart cell electrophysiology and Ca²(+) cycling has led to an evolution of computational models that have played a crucial role in shaping our understanding of heart function. An important emerging theme in systems biology is that the detailed nature of local signaling events, such as those that occur in the cardiac dyad, have important consequences at higher biological scales. Multiscale modeling techniques have revealed many mechanistic links between microscale events, such as Ca²(+) binding to a channel protein, and macroscale phenomena, such as excitation-contraction coupling gain. Here, we review experimentally based multiscale computational models of excitation-contraction coupling and the insights that have been gained through their application.

Modeling CICR in rat ventricular myocytes: voltage clamp studies

Background

The past thirty-five years have seen an intense search for the molecular mechanisms underlying calcium-induced calcium-release (CICR) in cardiac myocytes, with voltage clamp (VC) studies being the leading tool employed. Several VC protocols including lowering of extracellular calcium to affect Ca²⁺ loading of the sarcoplasmic reticulum (SR), and administration of blockers caffeine and thapsigargin have been utilized to probe the phenomena surrounding SR Ca²⁺ release. Here, we develop a deterministic mathematical model of a rat ventricular myocyte under VC conditions, to better understand mechanisms underlying the response of an isolated cell to calcium perturbation. Motivation for the study was to pinpoint key control variables influencing CICR and examine the role of CICR in the context of a physiological control system regulating cytosolic Ca²⁺ concentration ([Ca²⁺]_myo).

Methods

The cell model consists of an electrical-equivalent model for the cell membrane and a fluid-compartment model describing the flux of ionic species between the extracellular and several intracellular compartments (cell cytosol, SR and the dyadic coupling unit (DCU), in which resides the mechanistic basis of CICR). The DCU is described as a controller-actuator mechanism, internally stabilized by negative feedback control of the unit's two diametrically-opposed Ca²⁺ channels (trigger-channel and release-channel). It releases Ca²⁺ flux into the cyto-plasm and is in turn enclosed within a negative feedback loop involving the SERCA pump, regulating[Ca²⁺]_myo.

Results

Our model reproduces measured VC data published by several laboratories, and generates graded Ca²⁺ release at high Ca²⁺ gain in a homeostatically-controlled environment where [Ca²⁺]_myo is precisely regulated. We elucidate the importance of the DCU elements in this process, particularly the role of the ryanodine receptor in controlling SR Ca²⁺ release, its activation by trigger Ca²⁺, and its refractory characteristics mediated by the luminal SR Ca²⁺ sensor. Proper functioning of the DCU, sodium-calcium exchangers and SERCA pump are important in achieving negative feedback control and hence Ca²⁺ homeostasis.

Conclusions

We examine the role of the above Ca²⁺ regulating mechanisms in handling various types of induced disturbances in Ca²⁺ levels by quantifying cellular Ca²⁺ balance. Our model provides biophysically-based explanations of phenomena associated with CICR generating useful and testable hypotheses.

Ca2+ alternans in a cardiac myocyte model that uses moment equations to represent heterogeneous junctional SR Ca2+

Multiscale whole-cell models that accurately represent local control of Ca2+-induced Ca2+ release in cardiac myocytes can reproduce high-gain Ca2+ release that is graded with changes in membrane potential. Using a recently introduced formalism that represents heterogeneous local Ca2+ using moment equations, we present a model of cardiac myocyte Ca2+ cycling that exhibits alternating sarcoplasmic reticulum (SR) Ca2+ release when periodically stimulated by depolarizing voltage pulses. The model predicts that the distribution of junctional SR [Ca2+] across a large population of Ca2+ release units is distinct on alternating cycles. Load-release and release-uptake functions computed from this model give insight into how Ca2+ fluxes and stimulation frequency combine to determine the presence or absence of Ca2+ alternans. Our results show that the conditions for the onset of Ca2+ alternans cannot be explained solely by the steepness of the load-release function, but that changes in the release-uptake process also play an important role. We analyze the effect of the junctional SR refilling time constant on Ca2+ alternans and conclude that physiologically realistic models of defective Ca2+ cycling must represent the dynamics of heterogeneous junctional SR [Ca2+] without assuming rapid equilibration of junctional and network SR [Ca2+].

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.

Physiologic gating properties of unitary cardiac L-type Ca2+ channels

The contraction of adult mammalian ventricular cardiomyocytes is triggered by the influx of Ca(2+) ions through sarcolemmal L-type Ca(2+) channels (LCCs). However, the gating properties of unitary LCCs under physiologic conditions have remained elusive. Towards this end, we investigated the voltage-dependence of the gating kinetics of unitary LCCs, with a physiologic concentration of Ca(2+) ions permeating the channel. Unitary LCC currents were recorded with 2mM external Ca(2+) ions (in the absence of LCC agonists), using cell-attached patches on K-depolarized adult rat ventricular myocytes. The voltage-dependence of the peak probability of channel opening (Po vs. Vm) displayed a maximum value of 0.3, a midpoint of -12 mV, and a slope factor of 8.5. The maximum value for Po of the unitary LCC was significantly higher than previously assumed, under physiologic conditions. We also found that the mean open dwell time of the unitary LCC increased twofold with depolarization, ranging from 0.53+/-0.02 ms at -30 mV to 1.08+/-0.03 ms at 0 mV. The increase in mean LCC open time with depolarization counterbalanced the decrease in the single LCC current amplitude; the latter due to the decrease in driving force for Ca(2+) ion entry. Thus, the average amount of Ca(2+) ions entering through an individual LCC opening ( approximately 300-400 ions) remained relatively constant over this range of potentials. These novel results establish the voltage-dependence of unitary LCC gating kinetics using a physiologic Ca(2+) ion concentration. Moreover, they provide insight into local Ca(2+)-induced Ca(2+) release and a more accurate basis for mathematical modeling of excitation-contraction coupling in cardiac myocytes.

Ryanodine receptor-mediated arrhythmias and sudden cardiac death

The cardiac ryanodine receptor-Ca²⁺ release channel (RyR2) is an essential sarcoplasmic reticulum (SR) transmembrane protein that plays a central role in excitation–contraction coupling (ECC) in cardiomyocytes. Aberrant spontaneous, diastolic Ca²⁺ leak from the SR due to dysfunctional RyR2 contributes to the formation of delayed after-depolarisations, which are thought to underlie the fatal arrhythmia that occurs in both heart failure (HF) and in catecholaminergic polymorphic ventricular tachycardia (CPVT). CPVT is an inherited disorder associated with mutations in either the RyR2 or a SR luminal protein, calsequestrin. RyR2 shows normal function at rest in CPVT but the RyR2 dysfunction is unmasked by physical exercise or emotional stress, suggesting abnormal RyR2 activation as an underlying mechanism. Several potential mechanisms have been advanced to explain the dysfunctional RyR2 observed in HF and CPVT, including enhanced RyR2 phosphorylation status, altered RyR2 regulation at luminal/cytoplasmic sites and perturbed RyR2 intra/inter-molecular interactions. This review considers RyR2 dysfunction in the context of the structural and functional modulation of the channel, and potential therapeutic strategies to stabilise RyR2 function in cardiac pathology.

Computational modelling identifies the impact of subtle anatomical variations between amphibian and mammalian skeletal muscle on spatiotemporal calcium dynamics

The physical sites of calcium entry and exit in the skeletal muscle cell are distinct and highly organised in space. It was investigated whether the highly structured spatial organisation of sites of Ca(2+) release, uptake and action in skeletal muscle cells substantially impacts the dynamics of cytosolic Ca(2+) handling and thereby the physiology of the cell. Hereto, the spatiotemporal dynamics of the free calcium distribution in a fast-twitch (FT) muscle sarcomere was studied using a reaction-diffusion computational model for two genotypes with known anatomical differences. A computational model of a murine FT muscle sarcomere is developed, de novo including a closed calcium mass balance to simulate spatiotemporal high stimulation frequency calcium dynamics at 35 degrees C. Literature data on high-frequency calcium dye measurements were used as a first step towards model validation. The murine and amphibian sarcomere models were phenotypically distinct to capture known differences in positions of troponin C, actin-myosin overlap and calcium release within the sarcomere between frog and mouse. The models predicted large calcium gradients throughout the myoplasm as well as differences in calcium concentrations near the mitochondria of frog and mouse. Furthermore, the predicted Ca(2+) concentration was high at positions where Ca(2+) has a regulatory function, close to the mitochondria and troponin C.

Alternans of cardiac calcium cycling in a cluster of ryanodine receptors: a simulation study

Mechanical alternans in cardiac muscle is associated with intracellular Ca(2+) alternans. Mechanisms underlying intracellular Ca(2+) alternans are unclear. In previous experimental studies, we produced alternans of systolic Ca(2+) under voltage clamp, either by partially inhibiting the Ca(2+) release mechanism, or by applying small depolarizing pulses. In each case, alternans relied on propagating waves of Ca(2+) release. The aim of this study is to investigate by computer modeling how alternans of systolic Ca(2+) is produced. A mathematical model of a cardiac cell with 75 coupled elements is developed, with each element contains L-type Ca(2+) current, a subspace into which Ca release takes place, a cytoplasmic space, sarcoplasmic reticulum (SR) release channels [ryanodine receptor (RyR)], and uptake sites (SERCA). Interelement coupling is via Ca(2+) diffusion between neighboring subspaces via cytoplasmic spaces and network SR spaces. Small depolarizing pulses were simulated by step changes of cell membrane potential (20 mV) with random block of L-type channels. Partial inhibition of the release mechanism is mimicked by applying a reduction of RyR open probability in response to full stimulation by L-type channels. In both cases, systolic alternans follow, consistent with our experimental observations, being generated by propagating waves of Ca(2+) release and sustained through alternation of SR Ca(2+) content. This study provides novel and fundamental insights to understand mechanisms that may underlie intracellular Ca(2+) alternans without the need for refractoriness of L-type Ca or RyR channels under rapid pacing.

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 Ca2+-induced Ca2+ 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) [Ca2+] conditioned on Ca2+ release unit (CaRU) state. When coupled to ordinary differential equations (ODEs) for the bulk myoplasmic and network SR [Ca2+], a realistic but minimal model of cardiac excitation-contraction coupling was produced that avoids the computationally demanding task of resolving spatial aspects of global Ca2+ signaling, while accurately representing heterogeneous local Ca2+ 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 [Ca2+] are much faster than those of junctional SR [Ca2+]. 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 [Ca2+] 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 [Ca2+] 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 [Ca2+], this moment-closure approach to simulating local control of excitation-contraction coupling produces high-gain Ca2+ 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 Ca2+-induced Ca2+ release during simulated two-pulse voltage-clamp protocols.

A calcium-induced calcium release mechanism mediated by calsequestrin

Calcium (Ca(2+))-induced Ca(2+) release (CICR) is widely accepted as the principal mechanism linking electrical excitation and mechanical contraction in cardiac cells. The CICR mechanism has been understood mainly based on binding of cytosolic Ca(2+) with ryanodine receptors (RyRs) and inducing Ca(2+) release from the sarcoplasmic reticulum (SR). However, recent experiments suggest that SR lumenal Ca(2+) may also participate in regulating RyR gating through calsequestrin (CSQ), the SR lumenal Ca(2+) buffer. We investigate how SR Ca(2+) release via RyR is regulated by Ca(2+) and calsequestrin (CSQ). First, a mathematical model of RyR kinetics is derived based on experimental evidence. We assume that the RyR has three binding sites, two cytosolic sites for Ca(2+) activation and inactivation, and one SR lumenal site for CSQ binding. The open probability (P(o)) of the RyR is found by simulation under controlled cytosolic and SR lumenal Ca(2+). Both peak and steady-state P(o) effectively increase as SR lumenal Ca(2+) increases. Second, we incorporate the RyR model into a CICR model that has both a diadic space and the junctional SR (jSR). At low jSR Ca(2+) loads, CSQs are more likely to bind with the RyR and act to inhibit jSR Ca(2+) release, while at high SR loads CSQs are more likely to detach from the RyR, thereby increasing jSR Ca(2+) release. Furthermore, this CICR model produces a nonlinear relationship between fractional jSR Ca(2+) release and jSR load. These findings agree with experimental observations in lipid bilayers and cardiac myocytes.

Macroscopic consequences of calcium signaling in microdomains: a first-passage-time approach

Calcium (Ca) plays an important role in regulating various cellular processes. In a variety of cell types, Ca signaling occurs within microdomains where channels deliver localized pulses of Ca activating a nearby collection of Ca-sensitive receptors. The small number of channels involved ensures that the signaling process is stochastic. The aggregate response of several thousand of these microdomains yields a whole-cell response which dictates the cell behavior. Here, we study the statistical properties of a population of these microdomains in response to a trigger signal. We use a first-passage-time approach to show analytically how Ca release in the whole cell depends on properties of Ca channels within microdomains. Using these results we explain for the first time the underlying mechanism for the graded relationship between Ca influx and Ca release in cardiac cells.

Excitation-contraction coupling and mitochondrial energetics

Cardiac excitation-contraction (EC) coupling consumes vast amounts of cellular energy, most of which is produced in mitochondria by oxidative phosphorylation. In order to adapt the constantly varying workload of the heart to energy supply, tight coupling mechanisms are essential to maintain cellular pools of ATP, phosphocreatine and NADH. To our current knowledge, the most important regulators of oxidative phosphorylation are ADP, Pi, and Ca2+. However, the kinetics of mitochondrial Ca2+-uptake during EC coupling are currently a matter of intense debate. Recent experimental findings suggest the existence of a mitochondrial Ca2+ microdomain in cardiac myocytes, justified by the close proximity of mitochondria to the sites of cellular Ca2+ release, i. e., the ryanodine receptors of the sarcoplasmic reticulum. Such a Ca2+ microdomain could explain seemingly controversial results on mitochondrial Ca2+ uptake kinetics in isolated mitochondria versus whole cardiac myocytes. Another important consideration is that rapid mitochondrial Ca2+ uptake facilitated by microdomains may shape cytosolic Ca2+ signals in cardiac myocytes and have an impact on energy supply and demand matching. Defects in EC coupling in chronic heart failure may adversely affect mitochondrial Ca2+ uptake and energetics, initiating a vicious cycle of contractile dysfunction and energy depletion. Future therapeutic approaches in the treatment of heart failure could be aimed at interrupting this vicious cycle.

Multiscale modeling of calcium signaling in the cardiac dyad

Calcium (Ca(2+))-induced Ca(2+)-release (CICR) takes place in spatially restricted microdomains known as dyads. The length scale over which CICR occurs is on the order of nanometers and relevant time scales range from micro- to milliseconds. Quantitative understanding of CICR therefore requires development of models that are applicable over a range of spatio-temporal scales. We will present several new approaches for multiscale modeling of CICR. First, we present a model of dyad Ca(2+) dynamics in which the Fokker-Planck equation (FPE) is solved for the probability P(x, t) that a Ca(2+) ion is located at dyad position x at time t. Using this model, we demonstrate that (a) Ca(2+) signaling in the dyad is mediated by approximately tens of Ca(2+) ions; (b) these signaling events are noisy due to the small number of ions involved; and (c) the geometry of the RyR (ryanodine receptors) protein may function to restrict the diffusion of and to "funnel" Ca(2+) ions to activation-binding sites on the RyR, thus increasing RyR open probability and excitation-contraction (EC) coupling gain. Simplification of this model to one in which the dyadic space is represented using a single compartment yields the stochastic local-control model of CICR developed previously. We have shown that this model captures fundamental properties of CICR, such as graded release and voltage-dependent gain, may be integrated within a model of the myocyte and may be simulated in reasonable times using a combination of efficient numerical methods and parallel computing, but remains too complex for general use in cell simulations. To address this problem, we show how separation of time scales may be used to formulate a model in which nearby L-type Ca(2+) channels (LCCs) and RyRs gate as a coupled system that may be described using low-dimensional systems of ordinary differential equations, thus reducing computational complexity while capturing fundamentally important properties of CICR. The simplified model may be solved many orders of magnitude faster than can either of the more detailed models, thus enabling incorporation into tissue-level simulations.

Modelling calcium microdomains using homogenisation

Microdomains of calcium (i.e., areas on the nanometer scale that have qualitatively different calcium concentrations from that in the bulk cytosol) are known to be important in many situations. In cardiac cells, for instance, a calcium microdomain between the L-type channels and the ryanodine receptors, the so-called diadic cleft, is where the majority of the control of calcium release occurs. In other cell types that exhibit calcium oscillations and waves, the importance of microdomains in the vicinity of clusters of inositol trisphosphate receptors, or between the endoplasmic reticulum (ER) and other internal organelles or the plasma membrane, is clear. Given the limits of computational power, it is not currently realistic to model an entire cellular cytoplasm by incorporating detailed structural information about the ER throughout the entire cytoplasm. Hence, most models use a homogenised approach, assuming that both cytoplasm and ER coexist at each point of the domain. Conversely, microdomain models can be constructed, in which detailed structural information can be incorporated, but, until now, methods have not been developed for linking such a microdomain model to a model at the level of the entire cell. Using the homogenisation approach we developed in an earlier paper [Goel, P., Friedman, A., Sneyd, J., 2006. Homogenization of the cell cytoplasm: the calcium bidomain equations. SIAM J. Multiscale Modeling Simulation, in press] we show how a multiscale model of a calcium microdomain can be constructed. In this model a detailed model of the microdomain (in which the ER and the cytoplasm are separate compartments) is coupled to a homogenised model of the entire cell in a rigorous way. Our method is illustrated by a simple model of the diadic cleft of a cardiac half-sarcomere.

Protein geometry and placement in the cardiac dyad influence macroscopic properties of calcium-induced calcium release

In cardiac ventricular myocytes, events crucial to excitation-contraction coupling take place in spatially restricted microdomains known as dyads. The movement and dynamics of calcium (Ca2+) ions in the dyad have often been described by assigning continuously valued Ca2+ concentrations to one or more dyadic compartments. However, even at its peak, the estimated number of free Ca2+ ions present in a single dyad is small (approximately 10-100 ions). This in turn suggests that modeling dyadic calcium dynamics using laws of mass action may be inappropriate. In this study, we develop a model of stochastic molecular signaling between L-type Ca2+ channels (LCCs) and ryanodine receptors (RyR2s) that describes: a), known features of dyad geometry, including the space-filling properties of key dyadic proteins; and b), movement of individual Ca2+ ions within the dyad, as driven by electrodiffusion. The model enables investigation of how local Ca2+ signaling is influenced by dyad structure, including the configuration of key proteins within the dyad, the location of Ca2+ binding sites, and membrane surface charges. Using this model, we demonstrate that LCC-RyR2 signaling is influenced by both the stochastic dynamics of Ca2+ ions in the dyad as well as the shape and relative positioning of dyad proteins. Results suggest the hypothesis that the relative placement and shape of the RyR2 proteins helps to "funnel" Ca2+ ions to RyR2 binding sites, thus increasing excitation-contraction coupling gain.

A probability density approach to modeling local control of calcium-induced calcium release in cardiac myocytes

We present a probability density approach to modeling localized Ca2+ influx via L-type Ca2+ channels and Ca2+-induced Ca2+ release mediated by clusters of ryanodine receptors during excitation-contraction coupling in cardiac myocytes. Coupled advection-reaction equations are derived relating the time-dependent probability density of subsarcolemmal subspace and junctional sarcoplasmic reticulum [Ca2+] conditioned on "Ca2+ release unit" state. When these equations are solved numerically using a high-resolution finite difference scheme and the resulting probability densities are coupled to ordinary differential equations for the bulk myoplasmic and sarcoplasmic reticulum [Ca2+], a realistic but minimal model of cardiac excitation-contraction coupling is produced. Modeling Ca2+ release unit activity using this probability density approach avoids the computationally demanding task of resolving spatial aspects of global Ca2+ signaling, while accurately representing heterogeneous local Ca2+ signals in a population of diadic subspaces and junctional sarcoplasmic reticulum depletion domains. The probability density approach is validated for a physiologically realistic number of Ca2+ release units and benchmarked for computational efficiency by comparison to traditional Monte Carlo simulations. In simulated voltage-clamp protocols, both the probability density and Monte Carlo approaches to modeling local control of excitation-contraction coupling produce high-gain Ca2+ release that is graded with changes in membrane potential, a phenomenon not exhibited by so-called "common pool" models. However, a probability density calculation can be significantly faster than the corresponding Monte Carlo simulation, especially when cellular parameters are such that diadic subspace [Ca2+] is in quasistatic equilibrium with junctional sarcoplasmic reticulum [Ca2+] and, consequently, univariate rather than multivariate probability densities may be employed.

The dynamics of luminal depletion and the stochastic gating of Ca2+-activated Ca2+ channels and release sites

Single channel models of intracellular calcium (Ca(2+)) channels such as the 1,4,5-trisphosphate receptor and ryanodine receptor often assume that Ca(2+)-dependent transitions are mediated by constant background cytosolic [Ca(2+)]. This assumption neglects the fact that Ca(2+) released by open channels may influence subsequent gating through the processes of Ca(2+)-activation or inactivation. Similarly, the influence of the dynamics of luminal depletion on the stochastic gating of intracellular Ca(2+) channels is often neglected, in spite of the fact that the sarco/endoplasmic reticulum [Ca(2+)] near the luminal face of intracellular Ca(2+) channels influences the driving force for Ca(2+), the rate of Ca(2+) release, and the magnitude and time course of the consequent increase in cytosolic domain [Ca(2+)]. Here we analyze how the steady-state open probability of several minimal Ca(2+)-regulated Ca(2+) channel models depends on the conductance of the channel and the time constants for the relaxation of elevated cytosolic [Ca(2+)] and depleted luminal [Ca(2+)] to the bulk [Ca(2+)] of both compartments. Our approach includes Monte Carlo simulation as well as numerical solution of a system of advection-reaction equations for the multivariate probability density of elevated cytosolic [Ca(2+)] and depleted luminal [Ca(2+)] conditioned on each state of the stochastically gating channel. Both methods are subsequently used to study the role of luminal depletion in the dynamics of Ca(2+) puff/spark termination in release sites composed of Ca(2+) channels that are activated, but not inactivated, by cytosolic Ca(2+). The probability density approach shows that such minimal Ca(2+) release site models may exhibit puff/spark-like dynamics in either of two distinct parameter regimes. In one case, puffs/spark termination is due to the process of stochastic attrition and facilitated by rapid Ca(2+) domain collapse [cf. DeRemigio, H., Smith, G., 2005. The dynamics of stochastic attrition viewed as an absorption time on a terminating Markov chain. Cell Calcium 38, 73-86]. In the second case, puff/spark termination is promoted by the local depletion of luminal Ca(2+).

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

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

Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions

In this paper, we formulate a model for human ventricular cells that is efficient enough for whole organ arrhythmia simulations yet detailed enough to capture the effects of cell level processes such as current blocks and channelopathies. The model is obtained from our detailed human ventricular cell model by using mathematical techniques to reduce the number of variables from 19 to nine. We carefully compare our full and reduced model at the single cell, cable and 2D tissue level and show that the reduced model has a very similar behaviour. Importantly, the new model correctly produces the effects of current blocks and channelopathies on AP and spiral wave behaviour, processes at the core of current day arrhythmia research. The new model is well over four times more efficient than the full model. We conclude that the new model can be used for efficient simulations of the effects of current changes on arrhythmias in the human heart.