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.)

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.

Objective: Magnesium regulates a large number of cellular processes. Small changes in intracellular free Mg(2+) ([Mg(2+)](i)) may have important effects on cardiac excitability and contractility. We investigated the effects of [Mg(2+)](i) on cardiac excitation-contraction coupling., Methods: We used our ionic-metabolic model that incorporates equations for Ca(2+) and Mg(2+) buffering and transport by ATP and ADP and equations for MgATP regulation of ion transporters (Na(+)-K(+) pump, sarcolemmal and sarcoplasmic Ca(2+) pumps)., Results: Model results indicate that variations in cytosolic Mg(2+) level might sensitively affect diastolic and systolic Ca(2+), sarcoplasmic Ca(2+) content, Ca(2+) influx through L-type channels, efficiency of the Na(+)/Ca(2+) exchanger and action potential shape. The analysis suggests that the most important reason for the observed effects is a modified normal function of sarcoplasmic Ca(2+)-ATPase pump by altered diastolic MgATP levels., Conclusion: The model is able to reproduce qualitatively a sequence of events that correspond well with experimental observations during cardiac excitation-contraction coupling in mammalian ventricular myocytes.

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

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca2+ channel clustering and allosteric activation of Na+/Ca2+ exchanger by L-type Ca2+ current affects intracellular Ca2+ dynamics. Our model includes a realistic 3D geometry of a single t tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially-distributed membrane ion-transporters (L-type Ca2+ channel, Na+/Ca2+ exchanger, sarcolemmal Ca2+ pump, sarcolemmal Ca2+ leak), and stationary and mobile Ca2+buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca2+ . We obtained parameters from voltage-clamp protocols of L-type Ca2+ current and line-scan recordings of Ca2+ concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca2+ transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca2+ concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca2+ crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca2+ flux distribution. The model additionally predicts that local Ca2+ trigger fluxes are at least 3-fold to 8-fold higher than the whole-cell Ca2+ trigger flux. We found also that the activation of allosteric Ca2+ - binding sites on the Na+/Ca2+ xchanger could provide a mechanism for regulating global and local Ca2+ trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na+/Ca2+ exchanger fluxes to intracellular Ca2+ dynamics. [ABSTRACT FROM AUTHOR]