A simplified local control model of calcium-induced calcium release in cardiac ventricular myocytes

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.

Numerical analysis of Ca2+ signaling in rat ventricular myocytes with realistic transverse-axial tubular geometry and inhibited sarcoplasmic reticulum

The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca²⁺). To investigate how the t-tubule microanatomy and the distribution of membrane Ca²⁺ flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca²⁺ signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca²⁺ signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca²⁺ flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca²⁺ buffers, including the Ca²⁺ indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca²⁺ transients was found when the Ca²⁺ flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca²⁺ signals are predicted. Even at modest elevation of Ca²⁺, reached during Ca²⁺ influx, large and steep Ca²⁺ gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca²⁺ flux density along the cell membrane support initiation and propagation of Ca²⁺ waves in rat myocytes.

In cardiac muscle cells, calcium (Ca²⁺) is best known for its role in contraction activation. A remarkable amount of quantitative data on cardiac cell structure, ion-transporting protein distributions and intracellular Ca²⁺ dynamics has been accumulated. Various alterations in the protein distributions or cell ultra-structure are now recognized to be the primary mechanisms of cardiac dysfunction in a diverse range of common pathologies including cardiac arrhythmias and hypertrophy. Using a 3-D computational model, incorporating more realistic transverse-axial t-tubule geometry and considering geometric irregularities and inhomogeneities in the distribution of ion-transporting proteins, we analyze several important spatial and temporal features of Ca²⁺ signaling in rat ventricular myocytes. This study demonstrates that the computational models could serve as powerful tools for prediction and analyses of how the Ca²⁺ dynamics and cardiac excitation-contraction coupling are regulated under normal conditions or certain pathologies. The use of computational and mathematical approaches will help also to better understand aspects of cell functions that are not currently amenable to experimental investigation.

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

Integrative modeling of the cardiac ventricular myocyte

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field. We then narrow our focus to models of the cardiac ventricular myocyte and describe these models in the context of their subcellular functional systems including dynamic models of voltage-gated ion channels, mitochondrial energy production, ATP-dependent and electrogenic membrane transporters, intracellular Ca dynamics, mechanical contraction, and regulatory signal transduction pathways. We describe key advances and limitations of the models as well as point to new directions for future modeling research. WIREs Syst Biol Med 2011 3 392-413 DOI: 10.1002/wsbm.122

Spontaneous Ca2+ sparks and Ca2+ homeostasis in a minimal model of permeabilized ventricular myocytes

Many issues remain unresolved concerning how local, subcellular Ca(2+) signals interact with bulk cellular concentrations to maintain homeostasis in health and disease. To aid in the interpretation of data obtained in quiescent ventricular myocytes, we present here a minimal whole cell model that accounts for both localized (subcellular) and global (cellular) aspects of Ca(2+) signaling. Using a minimal formulation of the distribution of local [Ca(2+)] associated with a large number of Ca(2+)-release sites, the model simulates both random spontaneous Ca(2+) sparks and the changes in myoplasmic and sarcoplasmic reticulum (SR) [Ca(2+)] that result from the balance between stochastic release and reuptake into the SR. Ca(2+)-release sites are composed of clusters of two-state ryanodine receptors (RyRs) that exhibit activation by local cytosolic [Ca(2+)] but no inactivation or regulation by luminal Ca(2+). Decreasing RyR open probability in the model causes a decrease in aggregate release flux and an increase in SR [Ca(2+)], regardless of whether RyR inhibition is mediated by a decrease in RyR open dwell time or an increase in RyR closed dwell time. The same balance of stochastic release and reuptake can be achieved, however, by either high-frequency/short-duration or low-frequency/long-duration Ca(2+) sparks. The results are well correlated with recent experimental observations using pharmacological RyR inhibitors and clarify those aspects of the release-reuptake balance that are inherent to the coupling between local and global Ca(2+) signals and those aspects that depend on molecular-level details. The model of Ca(2+) sparks and homeostasis presented here can be a useful tool for understanding changes in cardiac Ca(2+ )release resulting from drugs, mutations, or acquired diseases.

High-throughput cardiac science on the Grid

Cardiac electrophysiology is a mature discipline, with the first model of a cardiac cell action potential having been developed in 1962. Current models range from single ion channels, through very complex models of individual cardiac cells, to geometrically and anatomically detailed models of the electrical activity in whole ventricles. A critical issue for model developers is how to choose parameters that allow the model to faithfully reproduce observed physiological effects without over-fitting. In this paper, we discuss the use of a parametric modelling toolkit, called Nimrod, that makes it possible both to explore model behaviour as parameters are changed and also to tune parameters by optimizing model output. Importantly, Nimrod leverages computers on the Grid, accelerating experiments by using available high-performance platforms. We illustrate the use of Nimrod with two case studies, one at the cardiac tissue level and one at the cellular level.

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.

A mathematical model of the murine ventricular myocyte: a data-driven biophysically based approach applied to mice overexpressing the canine NCX isoform

Mathematical modeling of Ca(2+) dynamics in the heart has the potential to provide an integrated understanding of Ca(2+)-handling mechanisms. However, many previous published models used heterogeneous experimental data sources from a variety of animals and temperatures to characterize model parameters and motivate model equations. This methodology limits the direct comparison of these models with any particular experimental data set. To directly address this issue, in this study, we present a biophysically based model of Ca(2+) dynamics directly fitted to experimental data collected in left ventricular myocytes isolated from the C57BL/6 mouse, the most commonly used genetic background for genetically modified mice in studies of heart diseases. This Ca(2+) dynamics model was then integrated into an existing mouse cardiac electrophysiology model, which was reparameterized using experimental data recorded at consistent and physiological temperatures. The model was validated against the experimentally observed frequency response of Ca(2+) dynamics, action potential shape, dependence of action potential duration on cycle length, and electrical restitution. Using this framework, the implications of cardiac Na(+)/Ca(2+) exchanger (NCX) overexpression in transgenic mice were investigated. These simulations showed that heterozygous overexpression of the canine cardiac NCX increases intracellular Ca(2+) concentration transient magnitude and sarcoplasmic reticulum Ca(2+) loading, in agreement with experimental observations, whereas acute overexpression of the murine cardiac NCX results in a significant loss of Ca(2+) from the cell and, hence, depressed sarcoplasmic reticulum Ca(2+) load and intracellular Ca(2+) concentration transient magnitude. From this analysis, we conclude that these differences are primarily due to the presence of allosteric regulation in the canine cardiac NCX, which has not been observed experimentally in the wild-type mouse heart.

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.

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.

Multiscale modeling in rodent ventricular myocytes

There is a growing body of experimental evidence suggesting that the Ca(2+) signaling in ventricular myocytes is characterized by a high gradient near the cell membrane and a more uniform Ca(2+) distribution in the cell interior [1]--[7]. An important reason for this phenomenon might be that in these cells the t-tubular system forms a network of extracellular space, extending deep into the cell interior. This allows the electrical signal, that propagates rapidly along the cell membrane, to reach the vicinity of the sarcoplasmic reticulum (SR), where intracellular Ca(2+) required for myofilament activation is stored [1], [8]--[11]. Early studies of cardiac muscle showed that the t-tubules are found at intervals of about 2 lm along the longitudinal cell axis in close proximity to the Z-disks of the sarcomeres [12]. Subsequent studies have demonstrated that the t-tubular system has also longitudinal extensions [9]--[11], [13].

K+ current changes account for the rate dependence of the action potential in the human atrial myocyte

Ongoing investigation of the electrophysiology and pathophysiology of the human atria requires an accurate representation of the membrane dynamics of the human atrial myocyte. However, existing models of the human atrial myocyte action potential do not accurately reproduce experimental observations with respect to the kinetics of key repolarizing currents or rate dependence of the action potential and fail to properly enforce charge conservation, an essential characteristic in any model of the cardiac membrane. In addition, recent advances in experimental methods have resulted in new data regarding the kinetics of repolarizing currents in the human atria. The goal of this study was to develop a new model of the human atrial action potential, based on the Nygren et al. model of the human atrial myocyte and newly available experimental data, that ensures an accurate representation of repolarization processes and reproduction of action potential rate dependence and enforces charge conservation. Specifically, the transient outward K(+) current (I(t)) and ultrarapid rectifier K(+) current (I(Kur)) were newly formulated. The inwardly recitifying K(+) current (I(K1)) was also reanalyzed and implemented appropriately. Simulations of the human atrial myocyte action potential with this new model demonstrated that early repolarization is dependent on the relative conductances of I(t) and I(Kur), whereas densities of both I(Kur) and I(K1) underlie later repolarization. In addition, this model reproduces experimental measurements of rate dependence of I(t), I(Kur), and action potential duration. This new model constitutes an improved representation of excitability and repolarization reserve in the human atrial myocyte and, therefore, provides a useful computational tool for future studies involving the human atrium in both health and disease.

Stochastic simulation of calcium microdomains in the vicinity of an L-type calcium channel

We study numerically the local dynamics of the intracellular calcium concentration in the vicinity of a voltage- and calcium-dependent plasma membrane L-type calcium channel. To account for the low number of Ca(2+) ions and buffer molecules present in sub-femtoliter volumes, we use an exact stochastic simulation algorithm including diffusion. We present a novel, unified simulation method that implements reaction-diffusion events of Ca(2+) ions and buffer molecules, stochastic ion channel gating and channel conductance as a multivariate Markov process. For fixed-voltage dynamics, e.g. under voltage-clamp conditions, it is shown that voltage-sensitive channel-gating steps can be incorporated exactly. We compare multi- and single-voxel geometries and show that the single-voxel approach leads to almost identical first- and second-order moments, at much lower computation time. Numerical examples illustrate the variability in local Ca(2+) fluctuations as induced by bursts of channel openings in response to membrane depolarisations. Finally, by introducing calmodulin as a link, it is shown how this variability is passed on to downstream signalling pathways. The method may prove useful to study calcium microdomains and calcium-regulated processes triggered by membrane depolarisations as evoked by, e.g., viral channel-forming proteins during virus-host cell interactions.

Markov models for ion channels: versatility versus identifiability and speed

Markov models (MMs) represent a generalization of Hodgkin-Huxley models. They provide a versatile structure for modelling single channel data, gating currents, state-dependent drug interaction data, exchanger and pump dynamics, etc. This paper uses examples from cardiac electrophysiology to discuss aspects related to parameter estimation. (i) Parameter unidentifiability (found in 9 out of 13 of the considered models) results in an inability to determine the correct layout of a model, contradicting the idea that model structure and parameters provide insights into underlying molecular processes. (ii) The information content of experimental voltage step clamp data is discussed, and a short but sufficient protocol for parameter estimation is presented. (iii) MMs have been associated with high computational cost (owing to their large number of state variables), presenting an obstacle for multicellular whole organ simulations as well as parameter estimation. It is shown that the stiffness of models increases computation time more than the number of states. (iv) Algorithms and software programs are provided for steady-state analysis, analytical solutions for voltage steps and numerical derivation of parameter identifiability. The results provide a new standard for ion channel modelling to further the automation of model development, the validation process and the predictive power of these models.

Coupling contraction, excitation, ventricular and coronary blood flow across scale and physics in the heart

In this study, we review the development and application of multi-physics and multi-scale coupling in the construction of whole-heart physiological models. Through an examination of recent computational modelling developments, we analyse the significance of coupling mechanisms for the increased understanding of cardiac function in the areas of excitation-contraction, coronary blood flow and ventricular fluid mechanical coupling. Within these physiological domains, we demonstrate and discuss the importance of model parametrization, imaging-based model anatomy and computational implementation.

The role of the Frank-Starling law in the transduction of cellular work to whole organ pump function: a computational modeling analysis

We have developed a multi-scale biophysical electromechanics model of the rat left ventricle at room temperature. This model has been applied to investigate the relative roles of cellular scale length dependent regulators of tension generation on the transduction of work from the cell to whole organ pump function. Specifically, the role of the length dependent Ca(2+) sensitivity of tension (Ca(50)), filament overlap tension dependence, velocity dependence of tension, and tension dependent binding of Ca(2+) to Troponin C on metrics of efficient transduction of work and stress and strain homogeneity were predicted by performing simulations in the absence of each of these feedback mechanisms. The length dependent Ca(50) and the filament overlap, which make up the Frank-Starling Law, were found to be the two dominant regulators of the efficient transduction of work. Analyzing the fiber velocity field in the absence of the Frank-Starling mechanisms showed that the decreased efficiency in the transduction of work in the absence of filament overlap effects was caused by increased post systolic shortening, whereas the decreased efficiency in the absence of length dependent Ca(50) was caused by an inversion in the regional distribution of strain.

Modelling and measuring electromechanical coupling in the rat heart

Tension-dependent binding of Ca(2+) to troponin C in the cardiac myocyte has been shown to play an important role in the regulation of Ca(2+) and the activation of tension development. The significance of this regulatory mechanism is quantified experimentally by the quantity of Ca(2+) released following a rapid change in the muscle length. Using a computational, coupled, electromechanics cell model, we have confirmed that the tension dependence of Ca(2+) binding to troponin C, rather than cross-bridge kinetics or the rate of Ca(2+) uptake by the sarcoplasmic reticulum, determines the quantity of Ca(2+) released following a length step. This cell model has been successfully applied in a continuum model of the papillary muscle to analyse experimental data, suggesting the tension-dependent binding of Ca(2+) to troponin C as the likely pathway through which the effects of localized impaired tension generation alter the Ca(2+) transient. These experimental results are qualitatively reproduced using a three-dimensional coupled electromechanics model. Furthermore, the model predicts that changes in the Ca(2+) transient in the viable myocardium surrounding the impaired region are amplified in the absence of tension-dependent binding of Ca(2+) to troponin C.

Early afterdepolarisations and ventricular arrhythmias in cardiac tissue: a computational study

Afterdepolarisations are associated with arrhythmias in the heart, but are difficult to study experimentally. In this study we used a simplified computational model of 1D and 2D cardiac ventricular tissue, where we could control the size of the region generating afterdepolarisations, as well as the properties of the afterdepolarisation waveform. Provided the size of the afterdepolarisation region was greater than around 1 mm, propagating extrasystoles were produced in both 1D and 2D. The number of extrasystoles produced depended on the amplitude, period, and duration of the oscillatory EAD waveform. In 2D, re-entry was also initiated for specific combinations of EAD amplitude, period, and duration, with the afterdepolarisation region acting as a common pathway. The main finding from this modelling study is therefore that afterdepolarisations can act as potent sources of propagating extrasystoles, as well as a source of re-entrant activation.

Mechanisms of transmurally varying myocyte electromechanics in an integrated computational model

The mechanical properties of myocardium vary across the transmural aspect of the left ventricular wall. Some of these functional heterogeneities may be related to differences in excitation-contraction coupling characteristics that have been observed in cells isolated from the epicardial, mid-myocardial and endocardial regions of the left ventricle of many species, including canine. Integrative models of coupled myocyte electromechanics are reviewed and used here to investigate sources of heterogeneous electromechanical behaviour in these cells. The simulations (i) illustrate a previously unrecognized role of the transient outward potassium current in mechanical function and (ii) suggest that there may also exist additional heterogeneities affecting crossbridge cycling rates in cells from different transmural regions.

From mitochondrial ion channels to arrhythmias in the heart: computational techniques to bridge the spatio-temporal scales

Computer simulations of electrical behaviour in the whole ventricles have become commonplace during the last few years. The goals of this article are (i) to review the techniques that are currently employed to model cardiac electrical activity in the heart, discussing the strengths and weaknesses of the various approaches, and (ii) to implement a novel modelling approach, based on physiological reasoning, that lifts some of the restrictions imposed by current state-of-the-art ionic models. To illustrate the latter approach, the present study uses a recently developed ionic model of the ventricular myocyte that incorporates an excitation-contraction coupling and mitochondrial energetics model. A paradigm to bridge the vastly disparate spatial and temporal scales, from subcellular processes to the entire organ, and from sub-microseconds to minutes, is presented. Achieving sufficient computational efficiency is the key to success in the quest to develop multiscale realistic models that are expected to lead to better understanding of the mechanisms of arrhythmia induction following failure at the organelle level, and ultimately to the development of novel therapeutic applications.

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

The sarcoplasmic reticulum (SR) in ventricular myocytes contains releasable Ca(2+) for activating cellular contraction. Recent measurements of intra-SR (luminal) Ca(2+) suggest a high diffusive Ca(2+)-mobility constant (D(CaSR)). This could help spatially to unify SR Ca(2+)-content ([Ca(2+)](SRT)) and standardize Ca(2+)-release throughout the cell. But measurements of localized depletions of luminal Ca(2+) (Ca(2+)-blinks), associated with local Ca(2+)-release (Ca(2+)-sparks), suggest D(CaSR) may actually be low. Here we describe a novel method for measuring D(CaSR). Using a cytoplasmic Ca(2+)-fluorophore, we estimate regional [Ca(2+)](SRT) from localized, caffeine-induced SR Ca(2+)-release. Caffeine microperfusion of one end of a guinea pig or rat myocyte diffusively empties the whole SR at a rate indicating D(CaSR) is 8-9 microm(2)/s, up to tenfold lower than previous estimates. Ignoring background SR Ca(2+)-leakage in our measurement protocol produces an artifactually high D(CaSR) (>40 microm(2)/s), which may also explain the previous high values. Diffusion-reaction modeling suggests that a low D(CaSR) would be sufficient to support local SR Ca(2+)-signaling within sarcomeres during excitation-contraction coupling. Low D(CaSR) also implies that [Ca(2+)](SRT) may readily become spatially nonuniform, particularly under pathological conditions of spatially nonuniform Ca(2+)-release. Local control of luminal Ca(2+), imposed by low D(CaSR), may complement the well-established local control of SR Ca(2+)-release by Ca(2+)-channel/ryanodine receptor couplons.

Using Physiome standards to couple cellular functions for rat cardiac excitation-contraction

Scientific endeavour is reliant upon the extension and reuse of previous knowledge. The formalization of this process for computational modelling is facilitated by the use of accepted standards with which to describe and simulate models, ensuring consistency between the models and thus reducing the development and propagation of errors. CellML 1.1, an XML-based programming language, has been designed as a modelling standard which, by virtue of its import and grouping functions, facilitates model combination and reuse. Using CellML 1.1, we demonstrate the process of formalized model reuse by combining three separate models of rat cardiomyocyte function (an electrophysiology model, a model of cellular calcium dynamics and a mechanics model) which together make up the Pandit-Hinch-Niederer et al. cell model. Not only is this integrative model of rat electromechanics a useful tool for cardiac modelling but it is also an ideal framework with which to demonstrate both the power of model reuse and the challenges associated with this process. We highlight and classify a number of these issues associated with combining models and provide some suggested solutions.

Stochastic binding of Ca2+ ions in the dyadic cleft; continuous versus random walk description of diffusion

Ca(2+) signaling in the dyadic cleft in ventricular myocytes is fundamentally discrete and stochastic. We study the stochastic binding of single Ca(2+) ions to receptors in the cleft using two different models of diffusion: a stochastic and discrete Random Walk (RW) model, and a deterministic continuous model. We investigate whether the latter model, together with a stochastic receptor model, can reproduce binding events registered in fully stochastic RW simulations. By evaluating the continuous model goodness-of-fit for a large range of parameters, we present evidence that it can. Further, we show that the large fluctuations in binding rate observed at the level of single time-steps are integrated and smoothed at the larger timescale of binding events, which explains the continuous model goodness-of-fit. With these results we demonstrate that the stochasticity and discreteness of the Ca(2+) signaling in the dyadic cleft, determined by single binding events, can be described using a deterministic model of Ca(2+) diffusion together with a stochastic model of the binding events, for a specific range of physiological relevant parameters. Time-consuming RW simulations can thus be avoided. We also present a new analytical model of bimolecular binding probabilities, which we use in the RW simulations and the statistical analysis.

Physiological modulation of voltage-dependent inactivation in the cardiac muscle L-type calcium channel: A modelling study

The inactivation of the L-type Ca2+ current is composed of voltage-dependent and calcium-dependent mechanisms. The relative contribution of these processes is still under dispute and the idea that the voltage-dependent inactivation could be subject to further modulation by other physiological processes had been ignored. This study sought to model physiological modulation of inactivation of the current in cardiac ventricular myocytes, based upon the recent detailed experimental data that separated total and voltage-dependent inactivation (VDI) by replacing extracellular Ca2+ with Mg2+ and monitoring L-type Ca2+ channel behaviour by outward K+ current flowing through the channel in the absence of inward current flow. Calcium-dependent inactivation (CDI) was based upon Ca2+ influx and formulated from data that was recorded during beta-adrenergic stimulation of the myocytes. Ca2+ influx and its competition with non-selective monovalent cation permeation were also incorporated into channel permeation in the model. The constructed model could closely reproduce the experimental Ba2+ and Ca2+ current results under basal condition where no beta-stimulation was added after a slight reduction of the development of fast voltage-dependent inactivation with depolarization. The model also predicted that under beta-adrenergic stimulation voltage-dependent inactivation is lost and calcium-dependent inactivation largely compensates it. The developed model thus will be useful to estimate the respective roles of VDI and CDI of L-type Ca2+ channels in various physiological and pathological conditions of the heart which would otherwise be difficult to show experimentally.

A kinetic model for calcium dynamics in RAW 264.7 cells: 1. Mechanisms, parameters, and subpopulational variability

Calcium (Ca(2+)) is an important second messenger and has been the subject of numerous experimental measurements and mechanistic studies in intracellular signaling. Calcium profile can also serve as a useful cellular phenotype. Kinetic models of calcium dynamics provide quantitative insights into the calcium signaling networks. We report here the development of a complex kinetic model for calcium dynamics in RAW 264.7 cells stimulated by the C5a ligand. The model is developed using the vast number of measurements of in vivo calcium dynamics carried out in the Alliance for Cellular Signaling (AfCS) Laboratories. Ligand binding, phospholipase C-beta (PLC-beta) activation, inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) dynamics, and calcium exchange with mitochondria and extracellular matrix have all been incorporated into the model. The experimental data include data from both native and knockdown cell lines. Subpopulational variability in measurements is addressed by allowing nonkinetic parameters to vary across datasets. The model predicts temporal response of Ca(2+) concentration for various doses of C5a under different initial conditions. The optimized parameters for IP(3)R dynamics are in agreement with the legacy data. Further, the half-maximal effect concentration of C5a and the predicted dose response are comparable to those seen in AfCS measurements. Sensitivity analysis shows that the model is robust to parametric perturbations.

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.

A mathematical model of the slow force response to stretch in rat ventricular myocytes

We developed a model of the rat ventricular myocyte at room temperature to predict the relative effects of different mechanisms on the cause of the slow increase in force in response to a step change in muscle length. We performed simulations in the presence of stretch-dependent increases in flux through the Na(+)-H(+) exchanger (NHE) and Cl(-)-HCO(3)(-) exchanger (AE), stretch-activated channels (SAC), and the stretch-dependent nitric oxide (NO) induced increased open probability of the ryanodine receptors to estimate the capacity of each mechanism to produce the slow force response (SFR). Inclusion of stretch-dependent NHE & AE, SACs, and stretch-dependent NO effects caused an increase in tension following 15 min of stretch of 0.87%, 32%, and 0%, respectively. Comparing [Ca(2+)](i) dynamics before and after stretch in the presence of combinations of the three stretch-dependent elements, which produced significant SFR values (>20%), showed that the inclusion of stretch-dependent NO effects produced [Ca(2+)](i) transients, which were not consistent with experimental results. Further simulations showed that in the presence of SACs and the absence of stretch-dependent NHE & AE inhibition of NHE attenuated the SFR, such that reduced SFR in the presence of NHE blockers does not indicate a stretch dependence of NHE. Rather, a functioning NHE is responsible for a portion of the SFR. Based on our simulations we estimate that in rat cardiac myocytes at room temperature SACs play a significant role in producing the SFR, potentially in the presence of stretch-dependent NHE & AE and that NO effects, if any, must involve more mechanisms than just increasing the open probability of ryanodine receptors.

Functional Significance of Na+/Ca2+ Exchangers Co-localization with Ryanodine Receptors

Co-localization of Na+/Ca2+ exchangers (NCX) with ryanodine receptors (RyRs) is debated. We incorporate local NCX current in a biophysically detailed model of L-type Ca2+ channels (LCCs) and RyRs and study the effect of NCX on the regulation of Ca2+-induced Ca2+ release and the shape of the action potential. In canine ventricular cells, under pathological conditions, e.g., impaired LCCs, local NCXs become an enhancer of sarcoplasmic reticulum release. Under such conditions incorporation of local NCXs is critical to accurately capture mechanisms of excitation-contraction coupling.

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.

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

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

Contributions of sustained INa and IKv43 to transmural heterogeneity of early repolarization and arrhythmogenesis in canine left ventricular myocytes

The roles of sustained components of I(Na) and I(Kv43) in shaping the action potentials (AP) of myocytes isolated from the canine left ventricle (LV) have not been studied in detail. Here we investigate the hypothesis that these two currents can contribute substantially to heterogeneity of early repolarization and arrhythmic risk. Quantitative data from voltage-clamp and expression profiling experiments were used to complete meaningful modifications to an existing "local control" model of canine midmyocardial myocyte excitation-contraction coupling for epicardial and endocardial cells. We include 1) heterogeneous I(Kv43), I(Ks), and I(SERCA) density; 2) modulation of I(Kv43) by Kv channel interacting protein type 2 (KChIP2) channel subunits; 3) a possible Ca(2+)-dependent open-state inactivation of I(Kv43); and 4) a sustained component of the inward Na(+) current, I(NaL). The resulting simulations illustrate ways in which KChIP2- and Ca(2+)-dependent control of I(Kv43) can result in a sustained outward current that can neutralize I(NaL) in a rate- and myocyte subtype-dependent manner. Both these currents appear to play significant roles in modulating AP duration and rate dependence in midmyocardial myocytes. Furthermore, an increased ratio of I(Kv43) to I(NaL) is capable of protecting epicardial myocytes from the early afterdepolarizations resulting from the SCN5A-I1768V mutation-induced increase in I(NaL). Experimentally observed transmural differences in Ca(2+) handling, including greater sarcoplasmic reticulum Ca(2+) content and faster Ca(2+) transient decay rates on the epicardium, were recapitulated in our simulations. By design, these models allow upward integration into organ models or may be used as a basis for further investigations into cellular heterogeneities.

Inhibition of cAMP-dependent protein kinase under conditions occurring in the cardiac dyad during a Ca2+ transient

The space between the t-tubule invagination and the sarcoplasmic reticulum (SR) membrane, the dyad, in ventricular myocytes has been predicted to experience very high [Ca2+] for short periods of time during a Ca2+ transient. The dyadic space accommodates many protein kinases responsible for the regulation of Ca2+ handling proteins of the cell. We show in vitro that cAMP-dependent protein kinase (PKA) is inhibited by high [Ca2+] through a shift in the ratio of CaATP/MgATP toward CaATP. We further generate a three-dimensional mathematical model of Ca2+ and ATP diffusion within dyad. We use this model to predict the extent to which PKA would be inhibited by an increased CaATP/MgATP ratio during a Ca2+ transient in the dyad in vivo. Our results suggest that under normal physiological conditions a myocyte paced at 1 Hz would experience up to 55% inhibition of PKA within the cardiac dyad, with inhibition averaging 5% throughout the transient, an effect which becomes more pronounced as the myocyte contractile frequency increases (at 7 Hz, PKA inhibition averages 28% across the dyad throughout the duration of a Ca2+ transient).

A mathematical analysis of the action potential plateau duration of a human ventricular myocyte

The plateau phase of a human ventricular myocyte is analysed. The plateau duration is a function of the time required for a myocyte's transmembrane voltage to decrease by a certain voltage, DeltaV. The timing of the plateau is shown to be controlled by two slowly changing gate variables, the inactivation gate that controls the inward/depolarizing L-type calcium current and the inactivation gate that controls the outward/repolarizing slow rectifier potassium current. The amount of current controlled by these variables is a function of the net conductivity of the corresponding sodium and potassium channels. An equation is derived that relates action potential duration to these net conductivities and the time dependence of the slowly moving variables. This equation is used to estimate plateau duration for a given value of DeltaV. The initial conditions of the slowly moving inactivation variables are shown to affect plateau duration. These initial conditions depend on the amount of time that has elapsed between a previous repolarization and a current depolarization (diastolic interval). The analysis thus helps to quantify the characteristics of action potential duration restitution.

Mechanisms of excitation-contraction coupling in an integrative model of the cardiac ventricular myocyte

It is now well established that characteristic properties of excitation-contraction (EC) coupling in cardiac myocytes, such as high gain and graded Ca(2+) release, arise from the interactions that occur between L-type Ca(2+) channels (LCCs) and nearby ryanodine-sensitive Ca(2+) release channels (RyRs) in localized microdomains. Descriptions of Ca(2+)-induced Ca(2+) release (CICR) that account for these local mechanisms are lacking from many previous models of the cardiac action potential, and those that do include local control of CICR are able to reconstruct properties of EC coupling, but require computationally demanding stochastic simulations of approximately 10(5) individual ion channels. In this study, we generalize a recently developed analytical approach for deriving simplified mechanistic models of CICR to formulate an integrative model of the canine cardiac myocyte which is computationally efficient. The resulting model faithfully reproduces experimentally measured properties of EC coupling and whole cell phenomena. The model is used to study the role of local redundancy in L-type Ca(2+) channel gating and the role of dyad configuration on EC coupling. Simulations suggest that the characteristic steep rise in EC coupling gain observed at hyperpolarized potentials is a result of increased functional coupling between LCCs and RyRs. We also demonstrate mechanisms by which alterations in the early repolarization phase of the action potential, resulting from reduction of the transient outward potassium current, alters properties of EC coupling.