Search

Little is known about how small variations in ionic currents and Ca²⁺ and Na⁺ diffusion coefficients impact action potential and Ca²⁺ dynamics in rabbit ventricular myocytes. We applied sensitivity analysis to quantify the sensitivity of Shannon et al. model (Biophys. J., 2004) to 5%-10% changes in currents conductance, channels distribution, and ion diffusion in rabbit ventricular cells. We found that action potential duration and Ca²⁺ peaks are highly sensitive to 10% increase in L-type Ca²⁺ current; moderately influenced by 10% increase in Na⁺-Ca²⁺ exchanger, Na⁺-K⁺ pump, rapid delayed and slow transient outward K⁺ currents, and Cl⁻ background current; insensitive to 10% increases in all other ionic currents and sarcoplasmic reticulum Ca²⁺ fluxes. Cell electrical activity is strongly affected by 5% shift of L-type Ca²⁺ channels and Na⁺-Ca²⁺ exchanger in between junctional and submembrane spaces while Ca²⁺-activated Cl⁻-channel redistribution has the modest effect. Small changes in submembrane and cytosolic diffusion coefficients for Ca²⁺, but not in Na⁺ transfer, may alter notably myocyte contraction. Our studies highlight the need for more precise measurements and further extending and testing of the Shannon et al. model. Our results demonstrate usefulness of sensitivity analysis to identify specific knowledge gaps and controversies related to ventricular cell electrophysiology and Ca²⁺ signaling.

Intercellular calcium waves in cardiac myocytes are a well-recognized, if incompletely understood, phenomenon. In a variety of preparations, investigators have reported multi-cellular calcium waves or triggered propagated contractions, but the mechanisms of propagation and pathological importance of these events remain unclear. Here, we review existing experimental data and present a computational approach to investigate the mechanisms of multi-cellular calcium wave propagation. Over the past 50 years, the standard modeling paradigm for excitable cardiac tissue has seen increasingly detailed models of the dynamics of individual cells coupled in tissue solely by intercellular and interstitial current flow. Although very successful, this modeling regime has been unable to capture two important phenomena: 1) the slow intercellular calcium waves observed experimentally, and 2) how intercellular calcium events resulting in delayed after depolarizations at the cellular level could overcome a source-sink mismatch to initiate depolarization waves in tissue. In this paper, we introduce a mathematical model with subcellular spatial resolution, in which we allow both inter- and intracellular current flow and calcium diffusion. In simulations of coupled cells employing this model, we observe: a) slow inter-cellular calcium waves propagating at about 0.1 mm/s, b) faster Calcium-Depolarization-Calcium (CDC) waves, traveling at about 1 mm/s, and c) CDC-waves that can set off fast depolarization-waves (50 cm/s) in tissue with varying gap-junction conductivity., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

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.

Spatial-temporal Ca(2+) dynamics due to Ca(2+) release, buffering, and reuptaking plays a central role in studying excitation-contraction (E-C) coupling in both normal and diseased cardiac myocytes. In this paper, we employ two numerical methods, namely, the meshless method and the finite element method, to model such Ca(2+) behaviors by solving a nonlinear system of reaction-diffusion partial differential equations at two scales. In particular, a subcellular model containing several realistic transverse tubules (or t-tubules) is investigated and assumed to reside at different locations relative to the cell membrane. To this end, the Ca(2+) concentration calculated from the whole-cell modeling is adopted as part of the boundary constraint in the subcellular model. The preliminary simulations show that Ca(2+) concentration changes in ventricular myocytes are mainly influenced by calcium release from t-tubules.

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(2+)). To investigate how the t-tubule microanatomy and the distribution of membrane Ca(2+) flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca(2+) 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(2+) signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca(2+) flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca(2+) buffers, including the Ca(2+) 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(2+) transients was found when the Ca(2+) flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca(2+) signals are predicted. Even at modest elevation of Ca(2+), reached during Ca(2+) influx, large and steep Ca(2+) gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca(2+) flux density along the cell membrane support initiation and propagation of Ca(2+) waves in rat 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].

Changes in cytosolic free Mg(2+) and adenosine nucleotide phosphates affect cardiac excitability and contractility. To investigate how modulation by Mg(2+), ATP, and ADP of K(ATP) and L-type Ca(2+) channels influences excitation-contraction coupling, we incorporated equations for intracellular ATP and MgADP regulation of the K(ATP) current and MgATP regulation of the L-type Ca(2+) current in an ionic-metabolic model of the canine ventricular myocyte. The new model: 1), quantitatively reproduces a dose-response relationship for the effects of changes in ATP on K(ATP) current, 2), simulates effects of ADP in modulating ATP sensitivity of K(ATP) channel, 3), predicts activation of Ca(2+) current during rapid increase in MgATP, and 4), demonstrates that decreased ATP/ADP ratio with normal total Mg(2+) or increased free Mg(2+) with normal ATP and ADP activate K(ATP) current, shorten action potential, and alter ionic currents and intracellular Ca(2+) signals. The model predictions are in agreement with experimental data measured under normal and a variety of pathological conditions.

Ca(2+) signaling in cells is largely governed by Ca(2+) diffusion and Ca(2+) binding to mobile and stationary Ca(2+) buffers, including organelles. To examine Ca(2+) signaling in cardiac atrial myocytes, a mathematical model of Ca(2+) diffusion was developed which represents several subcellular compartments, including a subsarcolemmal space with restricted diffusion, a myofilament space, and the cytosol. The model was used to quantitatively simulate experimental Ca(2+) signals in terms of amplitude, time course, and spatial features. For experimental reference data, L-type Ca(2+) currents were recorded from atrial cells with the whole-cell voltage-clamp technique. Ca(2+) signals were simultaneously imaged with the fluorescent Ca(2+) indicator Fluo-3 and a laser-scanning confocal microscope. The simulations indicate that in atrial myocytes lacking T-tubules, Ca(2+) movement from the cell membrane to the center of the cells relies strongly on the presence of mobile Ca(2+) buffers, particularly when the sarcoplasmic reticulum is inhibited pharmacologically. Furthermore, during the influx of Ca(2+) large and steep concentration gradients are predicted between the cytosol and the submicroscopically narrow subsarcolemmal space. In addition, the computations revealed that, despite its low Ca(2+) affinity, ATP acts as a significant buffer and carrier for Ca(2+), even at the modest elevations of [Ca(2+)](i) reached during influx of Ca(2+).