Computational Modeling and Numerical Methods for Spatiotemporal Calcium Cycling in Ventricular Myocytes

Intracellular calcium (Ca) cycling dynamics in cardiac myocytes is regulated by a complex network of spatially distributed organelles, such as sarcoplasmic reticulum (SR), mitochondria, and myofibrils. In this study, we present a mathematical model of intracellular Ca cycling and numerical and computational methods for computer simulations. The model consists of a coupled Ca release unit (CRU) network, which includes a SR domain and a myoplasm domain. Each CRU contains 10 L-type Ca channels and 100 ryanodine receptor channels, with individual channels simulated stochastically using a variant of Gillespie’s method, modified here to handle time-dependent transition rates. Both the SR domain and the myoplasm domain in each CRU are modeled by 5 × 5 × 5 voxels to maintain proper Ca diffusion. Advanced numerical algorithms implemented on graphical processing units were used for fast computational simulations. For a myocyte containing 100 × 20 × 10 CRUs, a 1-s heart time simulation takes about 10 min of machine time on a single NVIDIA Tesla C2050. Examples of simulated Ca cycling dynamics, such as Ca sparks, Ca waves, and Ca alternans, are shown.

[1]  Wei Chen,et al.  Role of coupled gating between cardiac ryanodine receptors in the genesis of triggered arrhythmias. , 2009, American journal of physiology. Heart and circulatory physiology.

[2]  G. Smith,et al.  Ryanodine receptor allosteric coupling and the dynamics of calcium sparks. , 2008, Biophysical journal.

[3]  Eric A Sobie,et al.  Moment closure for local control models of calcium-induced calcium release in cardiac myocytes. , 2008, Biophysical journal.

[4]  R. Winslow,et al.  An integrative model of the cardiac ventricular myocyte incorporating local control of Ca2+ release. , 2002, Biophysical journal.

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

[6]  H.J. Chizeck,et al.  The Physiome Projects and Multiscale Modeling [Life sciences] , 2008, IEEE Signal Processing Magazine.

[7]  Y. Shiferaw,et al.  Variability in Timing of Spontaneous Calcium Release in the Intact Rat Heart Is Determined by the Time Course of Sarcoplasmic Reticulum Calcium Load , 2010, Circulation research.

[8]  Alan Garfinkel,et al.  Spark-Induced Sparks As a Mechanism of Intracellular Calcium Alternans in Cardiac Myocytes , 2010, Circulation research.

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

[10]  Isuru D. Jayasinghe,et al.  Optical single-channel resolution imaging of the ryanodine receptor distribution in rat cardiac myocytes , 2009, Proceedings of the National Academy of Sciences.

[11]  M. Berridge,et al.  The versatility and universality of calcium signalling , 2000, Nature Reviews Molecular Cell Biology.

[12]  Peter Hunter,et al.  A strategy for integrative computational physiology. , 2005, Physiology.

[13]  Y. Tu,et al.  Excitation-contraction coupling gain and cooperativity of the cardiac ryanodine receptor: a modeling approach. , 2005, Biophysical journal.

[14]  Christian Soeller,et al.  Analysis of ryanodine receptor clusters in rat and human cardiac myocytes , 2007, Proceedings of the National Academy of Sciences.

[15]  R. Winslow,et al.  Cardiac Ca2+ dynamics: the roles of ryanodine receptor adaptation and sarcoplasmic reticulum load. , 1998, Biophysical journal.

[16]  D. Bers.  Cardiac excitation–contraction coupling , 2002, Nature.

[17]  A. Garfinkel,et al.  Alternans and Arrhythmias: From Cell to Heart , 2011, Circulation research.

[18]  Eric A Sobie,et al.  Termination of cardiac Ca(2+) sparks: an investigative mathematical model of calcium-induced calcium release. , 2002, Biophysical journal.

[19]  W. Lederer,et al.  Calcium sparks: elementary events underlying excitation-contraction coupling in heart muscle. , 1993, Science.

[20]  J. Ulrich.  [Physiology of the heart]. , 1950, Zeitschrift fur Kreislaufforschung.

[21]  F. Protasi,et al.  Shape, size, and distribution of Ca(2+) release units and couplons in skeletal and cardiac muscles. , 1999, Biophysical journal.

[22]  Sheng Wei,et al.  T-Tubule Remodeling During Transition From Hypertrophy to Heart Failure , 2010, Circulation research.

[23]  Ona Z Liu,et al.  Does the Goldilocks Principle apply to calcium release restitution in heart cells? , 2012, Journal of molecular and cellular cardiology.

[24]  S. Rush,et al.  A Practical Algorithm for Solving Dynamic Membrane Equations , 1978, IEEE Transactions on Biomedical Engineering.

[25]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[26]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[27]  Donald M. Bers,et al.  Allosteric Regulation of Na/Ca Exchange Current by Cytosolic Ca in Intact Cardiac Myocytes , 2001, The Journal of general physiology.

[28]  James D. Murray.  Mathematical Biology: I. An Introduction , 2007 .

[29]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[30]  Heping Cheng,et al.  Calcium sparks. , 2008, Physiological reviews.

[31]  Alan Garfinkel,et al.  Multi-scale modeling in biology: how to bridge the gaps between scales? , 2011, Progress in biophysics and molecular biology.

[32]  D. Bers.  Calcium cycling and signaling in cardiac myocytes. , 2008, Annual review of physiology.

[33]  X. Wehrens,et al.  Phosphorylation of RyR2 and shortening of RyR2 cluster spacing in spontaneously hypertensive rat with heart failure. , 2007, American journal of physiology. Heart and circulatory physiology.

[34]  A. Garfinkel,et al.  An advanced algorithm for solving partial differential equation in cardiac conduction , 1999, IEEE Transactions on Biomedical Engineering.

[35]  Gregory D Smith,et al.  Ca2+ alternans in a cardiac myocyte model that uses moment equations to represent heterogeneous junctional SR Ca2+. , 2010, Biophysical journal.

[36]  C. Soeller,et al.  Examination of the transverse tubular system in living cardiac rat myocytes by 2-photon microscopy and digital image-processing techniques. , 1999, Circulation research.

[37]  A Garfinkel,et al.  Model of intracellular calcium cycling in ventricular myocytes. , 2003, Biophysical journal.

[38]  Eric A Sobie,et al.  Predicting local SR Ca(2+) dynamics during Ca(2+) wave propagation in ventricular myocytes. , 2010, Biophysical journal.

[39]  J. R. Clay,et al.  Relationship between membrane excitability and single channel open-close kinetics. , 1983, Biophysical journal.

[40]  J. Keizer,et al.  Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations. , 1994, Biophysical journal.

[41]  Richard Haberman,et al.  Applied Partial Differential Equations with Fourier Series and Boundary Value Problems , 2012 .

[42]  Daisuke Sato,et al.  Acceleration of cardiac tissue simulation with graphic processing units , 2009, Medical & Biological Engineering & Computing.

[43]  Edward G Lakatta,et al.  A coupled SYSTEM of intracellular Ca2+ clocks and surface membrane voltage clocks controls the timekeeping mechanism of the heart's pacemaker. , 2010, Circulation research.

[44]  Donald M Bers,et al.  A mathematical treatment of integrated Ca dynamics within the ventricular myocyte. , 2004, Biophysical journal.

[45]  J. Shadid,et al.  Interplay of ryanodine receptor distribution and calcium dynamics. , 2006, Biophysical journal.

[46]  B. O’Rourke,et al.  Mitochondrial Ca2+ uptake: tortoise or hare? , 2009, Journal of molecular and cellular cardiology.

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