A multiscale computational model of spatially resolved calcium cycling in cardiac myocytes: from detailed cleft dynamics to the whole cell concentration profiles

Mathematical modeling of excitation-contraction coupling (ECC) in ventricular cardiac myocytes is a multiscale problem, and it is therefore difficult to develop spatially detailed simulation tools. ECC involves gradients on the length scale of 100 nm in dyadic spaces and concentration profiles along the 100 μm of the whole cell, as well as the sub-millisecond time scale of local concentration changes and the change of lumenal Ca2+ content within tens of seconds. Our concept for a multiscale mathematical model of Ca2+ -induced Ca2+ release (CICR) and whole cardiomyocyte electrophysiology incorporates stochastic simulation of individual LC- and RyR-channels, spatially detailed concentration dynamics in dyadic clefts, rabbit membrane potential dynamics, and a system of partial differential equations for myoplasmic and lumenal free Ca2+ and Ca2+-binding molecules in the bulk of the cell. We developed a novel computational approach to resolve the concentration gradients from dyadic space to cell level by using a quasistatic approximation within the dyad and finite element methods for integrating the partial differential equations. We show whole cell Ca2+-concentration profiles using three previously published RyR-channel Markov schemes.

[1]  W. Wier,et al.  Voltage dependence of intracellular [Ca2+]i transients in guinea pig ventricular myocytes. , 1987, Circulation research.

[2]  Takumi Washio,et al.  A three-dimensional simulation model of cardiomyocyte integrating excitation-contraction coupling and metabolism. , 2011, Biophysical journal.

[3]  Michael D. Stern,et al.  Local Control Models of Cardiac Excitation–Contraction Coupling , 1999, The Journal of general physiology.

[4]  Donald M. Bers,et al.  Na+-Ca2+ Exchange Current and Submembrane [Ca2+] During the Cardiac Action Potential , 2002, Circulation research.

[5]  Frank B Sachse,et al.  A modified local control model for Ca2+ transients in cardiomyocytes: junctional flux is accompanied by release from adjacent non-junctional RyRs. , 2014, Journal of molecular and cellular cardiology.

[6]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[7]  Jens Lang,et al.  Adaptive Multilevel Solution of Nonlinear Parabolic PDE Systems - Theory, Algorithm, and Applications , 2001, Lecture Notes in Computational Science and Engineering.

[8]  A. Levchenko,et al.  A 3D Monte Carlo analysis of the role of dyadic space geometry in spark generation. , 2006, Biophysical journal.

[9]  Donald M Bers,et al.  Voltage Dependence of Cardiac Excitation–Contraction Coupling: Unitary Ca2+ Current Amplitude and Open Channel Probability , 2007, Circulation research.

[10]  Johan Hake,et al.  Stochastic Binding of Ca2+ Ions in the Dyadic Cleft; Continuous versus Random Walk Description of Diffusion , 2008, Biophysical journal.

[11]  Eric A Sobie,et al.  Models of cardiac excitation-contraction coupling in ventricular myocytes. , 2010, Mathematical biosciences.

[12]  Joseph L Greenstein,et al.  The role of stochastic and modal gating of cardiac L-type Ca2+ channels on early after-depolarizations. , 2005, Biophysical journal.

[13]  Zhilin Qu,et al.  Computational Modeling and Numerical Methods for Spatiotemporal Calcium Cycling in Ventricular Myocytes , 2012, Front. Physio..

[14]  P. Dan,et al.  Distribution of proteins implicated in excitation-contraction coupling in rat ventricular myocytes. , 2000, Biophysical journal.

[15]  Zeyun Yu,et al.  Modelling cardiac calcium sparks in a three‐dimensional reconstruction of a calcium release unit , 2012, The Journal of physiology.

[16]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part II: implementation and tests in DUNE , 2008, Computing.

[17]  R. Gilmour,et al.  Memory models for the electrical properties of local cardiac systems. , 1997, Journal of theoretical biology.

[18]  J. Restrepo,et al.  A rabbit ventricular action potential model replicating cardiac dynamics at rapid heart rates. , 2007, Biophysical journal.

[19]  Martin Falcke,et al.  How does the ryanodine receptor in the ventricular myocyte wake up: by a single or by multiple open L-type Ca2+ channels? , 2011, European Biophysics Journal.

[20]  Joseph L Greenstein,et al.  Protein geometry and placement in the cardiac dyad influence macroscopic properties of calcium-induced calcium release. , 2007, Biophysical journal.

[21]  Brian O'Rourke,et al.  Elevated Cytosolic Na+ Decreases Mitochondrial Ca2+ Uptake During Excitation–Contraction Coupling and Impairs Energetic Adaptation in Cardiac Myocytes , 2006, Circulation research.

[22]  A. Klimek,et al.  The Molecular Architecture of Calcium Microdomains in Rat Cardiomyocytes , 2002, Annals of the New York Academy of Sciences.

[23]  Eric A. Sobie,et al.  Dynamics of calcium sparks and calcium leak in the heart. , 2011, Biophysical journal.

[24]  Isuru D. Jayasinghe,et al.  Organization of ryanodine receptors, transverse tubules, and sodium-calcium exchanger in rat myocytes. , 2009, Biophysical journal.

[25]  A. Fabiato,et al.  Contractions induced by a calcium‐triggered release of calcium from the sarcoplasmic reticulum of single skinned cardiac cells. , 1975, The Journal of physiology.

[26]  David R L Scriven,et al.  Ca²⁺ channel and Na⁺/Ca²⁺ exchange localization in cardiac myocytes. , 2013, Journal of molecular and cellular cardiology.

[27]  Godfrey L. Smith,et al.  Calcium movement in cardiac mitochondria. , 2014, Biophysical journal.

[28]  M. Cannell,et al.  Control of sarcoplasmic reticulum Ca2+ release by stochastic RyR gating within a 3D model of the cardiac dyad and importance of induction decay for CICR termination. , 2013, Biophysical journal.

[29]  C. Nagaiah.  Whole-cell simulation of hybride stochastic and deterministic calcium dynamics in 3D geometry. , 2012 .

[30]  B. Herman,et al.  Measurement of intracellular calcium. , 1999, Physiological reviews.

[31]  Juan G Restrepo,et al.  Spatiotemporal intracellular calcium dynamics during cardiac alternans. , 2009, Chaos.

[32]  M. Falcke,et al.  Release currents of IP(3) receptor channel clusters and concentration profiles. , 2004, Biophysical journal.

[33]  Martin Falcke,et al.  Efficient and detailed model of the local Ca2+ release unit in the ventricular cardiac myocyte. , 2010, Genome informatics. International Conference on Genome Informatics.

[34]  M. Stern,et al.  Theory of excitation-contraction coupling in cardiac muscle. , 1992, Biophysical journal.

[35]  W. Lederer,et al.  Effect of membrane potential changes on the calcium transient in single rat cardiac muscle cells. , 1987, Science.

[36]  M. Morad,et al.  Two-dimensional confocal images of organization, density, and gating of focal Ca2+ release sites in rat cardiac myocytes. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Joseph L Greenstein,et al.  Superresolution Modeling of Calcium Release in the Heart , 2014, Biophysical journal.

[38]  James P. Keener,et al.  Mathematical physiology , 1998 .

[39]  Gerald Warnecke,et al.  Adaptive space and time numerical simulation of reaction-diffusion models for intracellular calcium dynamics , 2012, Appl. Math. Comput..

[40]  M Falcke,et al.  Quasi-steady approximation for ion channel currents. , 2007, Biophysical journal.

[41]  Chamakuri Nagaiah,et al.  Whole-cell simulations of hybrid stochastic and deterministic calcium dynamics in 3 D geometry , 2012 .

[42]  C. Soeller,et al.  Numerical simulation of local calcium movements during L-type calcium channel gating in the cardiac diad. , 1997, Biophysical journal.