A Multiscale Computational Model of Spatially Resolved Calcium Cycling in Cardiac Myocytes: from Detailed Cleft Dynamics to the Whole Cell Concentration Profiles

Overview

Journal Front Physiol

Date 2015 Oct 7

PMID 26441674

Citations 7

Authors

Janine Vierheller

Wilhelm Neubert

Martin Falcke

Stephen H Gilbert

Nagaiah Chamakuri

Affiliations

Soon will be listed here.

Abstract

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 Ca(2+) content within tens of seconds. Our concept for a multiscale mathematical model of Ca(2+) -induced Ca(2+) 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 Ca(2+) and Ca(2+)-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 Ca(2+)-concentration profiles using three previously published RyR-channel Markov schemes.

Citing Articles

T cell Ca microdomains through the lens of computational modeling.

Gil Montoya D, Ornelas-Guevara R, Diercks B, Guse A, Dupont G Front Immunol. 2023; 14:1235737.

PMID: 37860008 PMC: 10582754. DOI: 10.3389/fimmu.2023.1235737.

Ca Release via IP Receptors Shapes the Cardiac Ca Transient for Hypertrophic Signaling.

Hunt H, Tilunaite A, Bass G, Soeller C, Roderick H, Rajagopal V Biophys J. 2020; 119(6):1178-1192.

PMID: 32871099 PMC: 7499065. DOI: 10.1016/j.bpj.2020.08.001.

Multiscale Modeling of Dyadic Structure-Function Relation in Ventricular Cardiac Myocytes.

Cosi F, Giese W, Neubert W, Luther S, Chamakuri N, Parlitz U Biophys J. 2019; 117(12):2409-2419.

PMID: 31635789 PMC: 6990380. DOI: 10.1016/j.bpj.2019.09.023.

NCX-Mediated Subcellular Ca Dynamics Underlying Early Afterdepolarizations in LQT2 Cardiomyocytes.

Zhong M, Rees C, Terentyev D, Choi B, Koren G, Karma A Biophys J. 2018; 115(6):1019-1032.

PMID: 30173888 PMC: 6140381. DOI: 10.1016/j.bpj.2018.08.004.

Studying dyadic structure-function relationships: a review of current modeling approaches and new insights into Ca (mis)handling.

Maleckar M, Edwards A, Louch W, Lines G Clin Med Insights Cardiol. 2017; 11:1179546817698602.

PMID: 28469494 PMC: 5392018. DOI: 10.1177/1179546817698602.

References

Williams G, Chikando A, Tuan H, Sobie E, Lederer W, Jafri M . Dynamics of calcium sparks and calcium leak in the heart. Biophys J. 2011; 101(6):1287-96. PMC: 3177068. DOI: 10.1016/j.bpj.2011.07.021. View

Stern M, Song L, Cheng H, Sham J, Yang H, Boheler K . Local control models of cardiac excitation-contraction coupling. A possible role for allosteric interactions between ryanodine receptors. J Gen Physiol. 1999; 113(3):469-89. PMC: 2222895. DOI: 10.1085/jgp.113.3.469. View

Williams G, Smith G, Sobie E, Jafri M . Models of cardiac excitation-contraction coupling in ventricular myocytes. Math Biosci. 2010; 226(1):1-15. PMC: 5499386. DOI: 10.1016/j.mbs.2010.03.005. View

Altamirano J, Bers D . Voltage dependence of cardiac excitation-contraction coupling: unitary Ca2+ current amplitude and open channel probability. Circ Res. 2007; 101(6):590-7. DOI: 10.1161/CIRCRESAHA.107.152322. View

Cannell M, Kong C, Imtiaz M, Laver D . 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. Biophys J. 2013; 104(10):2149-59. PMC: 3660628. DOI: 10.1016/j.bpj.2013.03.058. View

Bentele K, Falcke M . Quasi-steady approximation for ion channel currents. Biophys J. 2007; 93(8):2597-608. PMC: 1989719. DOI: 10.1529/biophysj.107.104299. View

Cleemann L, Wang W, Morad M . Two-dimensional confocal images of organization, density, and gating of focal Ca2+ release sites in rat cardiac myocytes. Proc Natl Acad Sci U S A. 1998; 95(18):10984-9. PMC: 28007. DOI: 10.1073/pnas.95.18.10984. View

Schendel T, Falcke M . Efficient and detailed model of the local Ca2+ release unit in the ventricular cardiac myocyte. Genome Inform. 2010; 22:142-55. View

Koh X, Srinivasan B, Ching H, Levchenko A . A 3D Monte Carlo analysis of the role of dyadic space geometry in spark generation. Biophys J. 2006; 90(6):1999-2014. PMC: 1386779. DOI: 10.1529/biophysj.105.065466. View

10.

Soeller C, Cannell M . Numerical simulation of local calcium movements during L-type calcium channel gating in the cardiac diad. Biophys J. 1997; 73(1):97-111. PMC: 1180912. DOI: 10.1016/S0006-3495(97)78051-2. View

11.

Hake J, Lines G . Stochastic binding of Ca2+ ions in the dyadic cleft; continuous versus random walk description of diffusion. Biophys J. 2008; 94(11):4184-201. PMC: 2480677. DOI: 10.1529/biophysj.106.103523. View

12.

Tanskanen A, Greenstein J, ORourke B, Winslow R . The role of stochastic and modal gating of cardiac L-type Ca2+ channels on early after-depolarizations. Biophys J. 2004; 88(1):85-95. PMC: 1201453. DOI: 10.1529/biophysj.104.051508. View

13.

Mahajan A, Shiferaw Y, Sato D, Baher A, Olcese R, Xie L . A rabbit ventricular action potential model replicating cardiac dynamics at rapid heart rates. Biophys J. 2007; 94(2):392-410. PMC: 2157228. DOI: 10.1529/biophysj.106.98160. View

14.

Weber C, Piacentino 3rd V, Ginsburg K, Houser S, Bers D . Na(+)-Ca(2+) exchange current and submembrane [Ca(2+)] during the cardiac action potential. Circ Res. 2002; 90(2):182-9. DOI: 10.1161/hh0202.103940. View

15.

Scriven D, Moore E . Ca²⁺ channel and Na⁺/Ca²⁺ exchange localization in cardiac myocytes. J Mol Cell Cardiol. 2012; 58:22-31. DOI: 10.1016/j.yjmcc.2012.11.022. View

16.

Restrepo J, Karma A . Spatiotemporal intracellular calcium dynamics during cardiac alternans. Chaos. 2009; 19(3):037115. PMC: 2771706. DOI: 10.1063/1.3207835. View

17.

Maack C, Cortassa S, Aon M, Ganesan A, Liu T, ORourke B . Elevated cytosolic Na+ decreases mitochondrial Ca2+ uptake during excitation-contraction coupling and impairs energetic adaptation in cardiac myocytes. Circ Res. 2006; 99(2):172-82. PMC: 2711867. DOI: 10.1161/01.RES.0000232546.92777.05. View

18.

Scriven D, Klimek A, Lee K, Moore E . The molecular architecture of calcium microdomains in rat cardiomyocytes. Ann N Y Acad Sci. 2002; 976:488-99. DOI: 10.1111/j.1749-6632.2002.tb04783.x. View

19.

Takahashi A, Camacho P, Lechleiter J, Herman B . Measurement of intracellular calcium. Physiol Rev. 1999; 79(4):1089-125. DOI: 10.1152/physrev.1999.79.4.1089. View

20.

Walker M, Williams G, Kohl T, Lehnart S, Jafri M, Greenstein J . Superresolution modeling of calcium release in the heart. Biophys J. 2014; 107(12):3018-3029. PMC: 4269784. DOI: 10.1016/j.bpj.2014.11.003. View