1
|
Zhang Y, Toyoda F, Himeno Y, Noma A, Amano A. Cell-specific models of hiPSC-CMs developed by the gradient-based parameter optimization method fitting two different action potential waveforms. Sci Rep 2024; 14:13086. [PMID: 38849433 PMCID: PMC11161598 DOI: 10.1038/s41598-024-63413-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 05/28/2024] [Indexed: 06/09/2024] Open
Abstract
Parameter optimization (PO) methods to determine the ionic current composition of experimental cardiac action potential (AP) waveform have been developed using a computer model of cardiac membrane excitation. However, it was suggested that fitting a single AP record in the PO method was not always successful in providing a unique answer because of a shortage of information. We found that the PO method worked perfectly if the PO method was applied to a pair of a control AP and a model output AP in which a single ionic current out of six current species, such as IKr, ICaL, INa, IKs, IKur or IbNSC was partially blocked in silico. When the target was replaced by a pair of experimental control and IKr-blocked records of APs generated spontaneously in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), the simultaneous fitting of the two waveforms by the PO method was hampered to some extent by the irregular slow fluctuations in the Vm recording and/or sporadic alteration in AP configurations in the hiPSC-CMs. This technical problem was largely removed by selecting stable segments of the records for the PO method. Moreover, the PO method was made fail-proof by running iteratively in identifying the optimized parameter set to reconstruct both the control and the IKr-blocked AP waveforms. In the lead potential analysis, the quantitative ionic mechanisms deduced from the optimized parameter set were totally consistent with the qualitative view of ionic mechanisms of AP so far described in physiological literature.
Collapse
Affiliation(s)
- Yixin Zhang
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Futoshi Toyoda
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan
- Central Research Laboratory, Shiga University of Medical Science, Otsu, Japan
| | - Yukiko Himeno
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan.
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan.
| | - Akinori Noma
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| |
Collapse
|
2
|
Himeno Y, Zhang Y, Enomoto S, Nomura H, Yamamoto N, Kiyokawa S, Ujihara M, Muangkram Y, Noma A, Amano A. Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model. Int J Mol Sci 2023; 24:15378. [PMID: 37895058 PMCID: PMC10607672 DOI: 10.3390/ijms242015378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 10/13/2023] [Indexed: 10/29/2023] Open
Abstract
Although repolarization has been suggested to propagate in cardiac tissue both theoretically and experimentally, it has been challenging to estimate how and to what extent the propagation of repolarization contributes to relaxation because repolarization only occurs in the course of membrane excitation in normal hearts. We established a mathematical model of a 1D strand of 600 myocytes stabilized at an equilibrium potential near the plateau potential level by introducing a sustained component of the late sodium current (INaL). By applying a hyperpolarizing stimulus to a small part of the strand, we succeeded in inducing repolarization which propagated along the strand at a velocity of 1~2 cm/s. The ionic mechanisms responsible for repolarization at the myocyte level, i.e., the deactivation of both the INaL and the L-type calcium current (ICaL), and the activation of the rapid component of delayed rectifier potassium current (IKr) and the inward rectifier potassium channel (IK1), were found to be important for the propagation of repolarization in the myocyte strand. Using an analogy with progressive activation of the sodium current (INa) in the propagation of excitation, regenerative activation of the predominant magnitude of IK1 makes the myocytes at the wave front start repolarization in succession through the electrical coupling via gap junction channels.
Collapse
Affiliation(s)
- Yukiko Himeno
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan; (Y.Z.); (A.N.); (A.A.)
| | | | | | | | | | | | | | | | | | | |
Collapse
|
3
|
Agrawal A, Wang K, Polonchuk L, Cooper J, Hendrix M, Gavaghan DJ, Mirams GR, Clerx M. Models of the cardiac L-type calcium current: A quantitative review. WIREs Mech Dis 2023; 15:e1581. [PMID: 36028219 PMCID: PMC10078428 DOI: 10.1002/wsbm.1581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Revised: 06/16/2022] [Accepted: 07/19/2022] [Indexed: 01/31/2023]
Abstract
The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.
Collapse
Affiliation(s)
- Aditi Agrawal
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | - Ken Wang
- Pharma Research and Early Development, Innovation Center BaselF. Hoffmann‐La Roche Ltd.BaselSwitzerland
| | - Liudmila Polonchuk
- Pharma Research and Early Development, Innovation Center BaselF. Hoffmann‐La Roche Ltd.BaselSwitzerland
| | - Jonathan Cooper
- Centre for Advanced Research ComputingUniversity College LondonLondonUK
| | - Maurice Hendrix
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
- Digital Research Service, Information SciencesUniversity of NottinghamNottinghamUK
| | - David J. Gavaghan
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | - Gary R. Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| | - Michael Clerx
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| |
Collapse
|
4
|
Qu Z, Yan D, Song Z. Modeling Calcium Cycling in the Heart: Progress, Pitfalls, and Challenges. Biomolecules 2022; 12:1686. [PMID: 36421700 PMCID: PMC9687412 DOI: 10.3390/biom12111686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Revised: 11/08/2022] [Accepted: 11/11/2022] [Indexed: 11/16/2022] Open
Abstract
Intracellular calcium (Ca) cycling in the heart plays key roles in excitation-contraction coupling and arrhythmogenesis. In cardiac myocytes, the Ca release channels, i.e., the ryanodine receptors (RyRs), are clustered in the sarcoplasmic reticulum membrane, forming Ca release units (CRUs). The RyRs in a CRU act collectively to give rise to discrete Ca release events, called Ca sparks. A cell contains hundreds to thousands of CRUs, diffusively coupled via Ca to form a CRU network. A rich spectrum of spatiotemporal Ca dynamics is observed in cardiac myocytes, including Ca sparks, spark clusters, mini-waves, persistent whole-cell waves, and oscillations. Models of different temporal and spatial scales have been developed to investigate these dynamics. Due to the complexities of the CRU network and the spatiotemporal Ca dynamics, it is challenging to model the Ca cycling dynamics in the cardiac system, particularly at the tissue sales. In this article, we review the progress of modeling of Ca cycling in cardiac systems from single RyRs to the tissue scale, the pros and cons of the current models and different modeling approaches, and the challenges to be tackled in the future.
Collapse
Affiliation(s)
- Zhilin Qu
- Department of Medicine, David Geffen School of Medicine, University of California, A2-237 CHS, 650 Charles E. Young Drive South, Los Angeles, CA 90095, USA
- Department of Computational Medicine, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
| | - Dasen Yan
- Peng Cheng Laboratory, Shenzhen 518066, China
| | - Zhen Song
- Peng Cheng Laboratory, Shenzhen 518066, China
| |
Collapse
|
5
|
Kohjitani H, Koda S, Himeno Y, Makiyama T, Yamamoto Y, Yoshinaga D, Wuriyanghai Y, Kashiwa A, Toyoda F, Zhang Y, Amano A, Noma A, Kimura T. Gradient-based parameter optimization method to determine membrane ionic current composition in human induced pluripotent stem cell-derived cardiomyocytes. Sci Rep 2022; 12:19110. [PMID: 36351955 PMCID: PMC9646722 DOI: 10.1038/s41598-022-23398-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Accepted: 10/31/2022] [Indexed: 11/11/2022] Open
Abstract
Premature cardiac myocytes derived from human induced pluripotent stem cells (hiPSC-CMs) show heterogeneous action potentials (APs), probably due to different expression patterns of membrane ionic currents. We developed a method for determining expression patterns of functional channels in terms of whole-cell ionic conductance (Gx) using individual spontaneous AP configurations. It has been suggested that apparently identical AP configurations can be obtained using different sets of ionic currents in mathematical models of cardiac membrane excitation. If so, the inverse problem of Gx estimation might not be solved. We computationally tested the feasibility of the gradient-based optimization method. For a realistic examination, conventional 'cell-specific models' were prepared by superimposing the model output of AP on each experimental AP recorded by conventional manual adjustment of Gxs of the baseline model. Gxs of 4-6 major ionic currents of the 'cell-specific models' were randomized within a range of ± 5-15% and used as an initial parameter set for the gradient-based automatic Gxs recovery by decreasing the mean square error (MSE) between the target and model output. Plotting all data points of the MSE-Gx relationship during optimization revealed progressive convergence of the randomized population of Gxs to the original value of the cell-specific model with decreasing MSE. The absence of any other local minimum in the global search space was confirmed by mapping the MSE by randomizing Gxs over a range of 0.1-10 times the control. No additional local minimum MSE was obvious in the whole parameter space, in addition to the global minimum of MSE at the default model parameter.
Collapse
Affiliation(s)
- Hirohiko Kohjitani
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Shigeya Koda
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Yukiko Himeno
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeru Makiyama
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yuta Yamamoto
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Daisuke Yoshinaga
- grid.258799.80000 0004 0372 2033Department Pediatrics, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yimin Wuriyanghai
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Asami Kashiwa
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Futoshi Toyoda
- grid.410827.80000 0000 9747 6806Department of Physiology, Shiga University of Medical Science, Otsu, Japan
| | - Yixin Zhang
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akinori Noma
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeshi Kimura
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| |
Collapse
|
6
|
Kim N, Pronto JD, Nickerson DP, Taberner AJ, Hunter PJ. A novel modular modeling approach for understanding different electromechanics between left and right heart in rat. Front Physiol 2022; 13:965054. [PMID: 36176770 PMCID: PMC9513479 DOI: 10.3389/fphys.2022.965054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 08/22/2022] [Indexed: 12/01/2022] Open
Abstract
While ion channels and transporters involved in excitation-contraction coupling have been linked and constructed as comprehensive computational models, validation of whether each individual component of a model can be reused has not been previously attempted. Here we address this issue while using a novel modular modeling approach to investigate the underlying mechanism for the differences between left ventricle (LV) and right ventricle (RV). Our model was developed from modules constructed using the module assembly principles of the CellML model markup language. The components of three existing separate models of cardiac function were disassembled as to create smaller modules, validated individually, and then the component parts were combined into a new integrative model of a rat ventricular myocyte. The model was implemented in OpenCOR using the CellML standard in order to ensure reproducibility. Simulated action potential (AP), Ca2+ transient, and tension were in close agreement with our experimental measurements: LV AP showed a prolonged duration and a more prominent plateau compared with RV AP; Ca2+ transient showed prolonged duration and slow decay in LV compared to RV; the peak value and relaxation of tension were larger and slower, respectively, in LV compared to RV. Our novel approach of module-based mathematical modeling has established that the ionic mechanisms underlying the APs and Ca2+ handling play a role in the variation in force production between ventricles. This simulation process also provides a useful way to reuse and elaborate upon existing models in order to develop a new model.
Collapse
Affiliation(s)
- Nari Kim
- NLRL for Innovative Cardiovascular Engineering, Department of Physiology, College of Medicine, Inje University, Busan, South Korea
- Cardiovascular and Metabolic Disease Center, Inje University, Busan, South Korea
- *Correspondence: Nari Kim,
| | - Julius D. Pronto
- NLRL for Innovative Cardiovascular Engineering, Department of Physiology, College of Medicine, Inje University, Busan, South Korea
- Cardiovascular and Metabolic Disease Center, Inje University, Busan, South Korea
| | - David P. Nickerson
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Andrew J. Taberner
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Peter J. Hunter
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
7
|
Niarakis A, Waltemath D, Glazier J, Schreiber F, Keating SM, Nickerson D, Chaouiya C, Siegel A, Noël V, Hermjakob H, Helikar T, Soliman S, Calzone L. Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology. Brief Bioinform 2022; 23:bbac212. [PMID: 35671510 PMCID: PMC9294410 DOI: 10.1093/bib/bbac212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 04/20/2022] [Accepted: 05/06/2022] [Indexed: 11/14/2022] Open
Abstract
Computational models are often employed in systems biology to study the dynamic behaviours of complex systems. With the rise in the number of computational models, finding ways to improve the reusability of these models and their ability to reproduce virtual experiments becomes critical. Correct and effective model annotation in community-supported and standardised formats is necessary for this improvement. Here, we present recent efforts toward a common framework for annotated, accessible, reproducible and interoperable computational models in biology, and discuss key challenges of the field.
Collapse
Affiliation(s)
- Anna Niarakis
- Université Paris-Saclay, Laboratoire Européen de Recherche pour la Polyarthrite rhumatoïde - Genhotel, Univ Evry, Evry, France
- Lifeware Group, Inria, Saclay-île de France, 91120 Palaiseau, France
| | - Dagmar Waltemath
- Department of Medical Informatics, University Medicine Greifswald, Greifswald, Germany
| | - James Glazier
- Biocomplexity Institute and Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | - Falk Schreiber
- Department of Computer and Information Science, University of Konstanz, Konstanz, Germany
- Faculty of Information Technology, Monash University, Clayton, Australia
| | | | - David Nickerson
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
| | | | - Anne Siegel
- Univ Rennes, CNRS, Inria - IRISA lab. Rennes
| | - Vincent Noël
- Institut Curie, PSL Research University, Paris, France
- INSERM, U900, Paris, France
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, Paris, France
| | - Henning Hermjakob
- EMBL-European Bioinformatics Institute, Wellcome Genome Campus, Cambridge, UK
| | - Tomáš Helikar
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, 68588, USA
| | - Sylvain Soliman
- Lifeware Group, Inria, Saclay-île de France, 91120 Palaiseau, France
| | - Laurence Calzone
- Institut Curie, PSL Research University, Paris, France
- INSERM, U900, Paris, France
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, Paris, France
| |
Collapse
|
8
|
Louch WE, Perdreau-Dahl H, Edwards AG. Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going? Front Physiol 2022; 13:834211. [PMID: 35356084 PMCID: PMC8959215 DOI: 10.3389/fphys.2022.834211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 01/31/2022] [Indexed: 11/24/2022] Open
Abstract
Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation–contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.
Collapse
Affiliation(s)
- William E. Louch
- Institute for Experimental Medical Research, Oslo University Hospital and University of Oslo, Oslo, Norway
- K.G. Jebsen Centre for Cardiac Research, University of Oslo, Oslo, Norway
| | - Harmonie Perdreau-Dahl
- Institute for Experimental Medical Research, Oslo University Hospital and University of Oslo, Oslo, Norway
- K.G. Jebsen Centre for Cardiac Research, University of Oslo, Oslo, Norway
| | - Andrew G. Edwards
- Simula Research Laboratory, Lysaker, Norway
- *Correspondence: Andrew G. Edwards,
| |
Collapse
|
9
|
Peng Z, Resnick A, Young YN. Primary cilium: a paradigm for integrating mathematical modeling with experiments and numerical simulations in mechanobiology. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:1215-1237. [PMID: 33757184 PMCID: PMC8552149 DOI: 10.3934/mbe.2021066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Primary cilia are non-motile, solitary (one per cell) microtubule-based organelles that emerge from the mother centriole after cells have exited the mitotic cycle. Identified as a mechanosensing organelle that responds to both mechanical and chemical stimuli, the primary cilium provides a fertile ground for integrative investigations of mathematical modeling, numerical simulations, and experiments. Recent experimental findings revealed considerable complexity to the underlying mechanosensory mechanisms that transmit extracellular stimuli to intracellular signaling many of which include primary cilia. In this invited review, we provide a brief survey of experimental findings on primary cilia and how these results lead to various mathematical models of the mechanics of the primary cilium bent under an external forcing such as a fluid flow or a trap. Mathematical modeling of the primary cilium as a fluid-structure interaction problem highlights the importance of basal anchorage and the anisotropic moduli of the microtubules. As theoretical modeling and numerical simulations progress, along with improved state-of-the-art experiments on primary cilia, we hope that details of ciliary regulated mechano-chemical signaling dynamics in cellular physiology will be understood in the near future.
Collapse
Affiliation(s)
- Zhangli Peng
- Department of Bioengineering, University of Illinois at Chicago, 851 S. Morgan St., Chicago, IL 60607, USA
| | - Andrew Resnick
- Department of Physics, Center for Gene Regulation in Health and Disease, Cleveland State University, Cleveland, OH 44115, USA
| | - Y.-N. Young
- Department of Mathematical Sciences, New Jersey Institute of Technology, University Heights, Newark, NJ 07102, USA
| |
Collapse
|
10
|
Corrado C, Avezzù A, Lee AWC, Mendoca Costa C, Roney CH, Strocchi M, Bishop M, Niederer SA. Using cardiac ionic cell models to interpret clinical data. WIREs Mech Dis 2020; 13:e1508. [PMID: 33027553 DOI: 10.1002/wsbm.1508] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 08/27/2020] [Accepted: 09/04/2020] [Indexed: 01/24/2023]
Abstract
For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.
Collapse
|
11
|
Marchena M, Echebarria B, Shiferaw Y, Alvarez-Lacalle E. Buffering and total calcium levels determine the presence of oscillatory regimes in cardiac cells. PLoS Comput Biol 2020; 16:e1007728. [PMID: 32970668 PMCID: PMC7537911 DOI: 10.1371/journal.pcbi.1007728] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 10/06/2020] [Accepted: 07/07/2020] [Indexed: 12/16/2022] Open
Abstract
Calcium oscillations and waves induce depolarization in cardiac cells which are believed to cause life-threathening arrhythimas. In this work, we study the conditions for the appearance of calcium oscillations in both a detailed subcellular model of calcium dynamics and a minimal model that takes into account just the minimal ingredients of the calcium toolkit. To avoid the effects of homeostatic changes and the interaction with the action potential we consider the somewhat artificial condition of a cell without pacing and with no calcium exchange with the extracellular medium. Both the full subcellular model and the minimal model present the same scenarios depending on the calcium load: two stationary states, one with closed ryanodine receptors (RyR) and most calcium in the cell stored in the sarcoplasmic reticulum (SR), and another, with open RyRs and a depleted SR. In between, calcium oscillations may appear. The robustness of these oscillations is determined by the amount of calsequestrin (CSQ). The lack of this buffer in the SR enhances the appearance of oscillations. The minimal model allows us to relate the stability of the oscillating state to the nullcline structure of the system, and find that its range of existence is bounded by a homoclinic and a Hopf bifurcation, resulting in a sudden transition to the oscillatory regime as the cell calcium load is increased. Adding a small amount of noise to the RyR behavior increases the parameter region where oscillations appear and provides a gradual transition from the resting state to the oscillatory regime, as observed in the subcellular model and experimentally. In cardiac cells, calcium plays a very important role. An increase in calcium levels is the trigger used by the cell to initiate contraction. Besides, calcium modulates several transmembrane currents, affecting the cell transmembrane potential. Thus, dysregulations in calcium handling have been associated with the appearance of arrhythmias. Often, this dysregulation results in the appearance of periodic calcium waves or global oscillations, providing a pro-arrhythmic substrate. In this paper, we study the onset of calcium oscillations in cardiac cells using both a detailed subcellular model of calcium dynamics and a minimal model that takes into account the essential ingredients of the calcium toolkit. Both reproduce the main experimental results and link this behavior with the presence of different steady-state solutions and bifurcations that depend on the total amount of calcium in the cell and in the level of buffering present. We expect that this work will help to clarify the conditions under which calcium oscillations appear in cardiac myocytes and, therefore, will represent a step further in the understanding of the origin of cardiac arrhythmias.
Collapse
Affiliation(s)
- Miquel Marchena
- Departament de Física, Universitat Politècnica de Catalunya-BarcelonaTech, Barcelona, Spain
| | - Blas Echebarria
- Departament de Física, Universitat Politècnica de Catalunya-BarcelonaTech, Barcelona, Spain
| | - Yohannes Shiferaw
- Physics Department, California State University, Northridge, California 91330, USA
| | | |
Collapse
|
12
|
Ca 2+ Release via IP 3 Receptors Shapes the Cardiac Ca 2+ Transient for Hypertrophic Signaling. Biophys J 2020; 119:1178-1192. [PMID: 32871099 DOI: 10.1016/j.bpj.2020.08.001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Revised: 06/16/2020] [Accepted: 08/04/2020] [Indexed: 12/27/2022] Open
Abstract
Calcium (Ca2+) plays a central role in mediating both contractile function and hypertrophic signaling in ventricular cardiomyocytes. L-type Ca2+ channels trigger release of Ca2+ from ryanodine receptors for cellular contraction, whereas signaling downstream of G-protein-coupled receptors stimulates Ca2+ release via inositol 1,4,5-trisphosphate receptors (IP3Rs), engaging hypertrophic signaling pathways. Modulation of the amplitude, duration, and duty cycle of the cytosolic Ca2+ contraction signal and spatial localization have all been proposed to encode this hypertrophic signal. Given current knowledge of IP3Rs, we develop a model describing the effect of functional interaction (cross talk) between ryanodine receptor and IP3R channels on the Ca2+ transient and examine the sensitivity of the Ca2+ transient shape to properties of IP3R activation. A key result of our study is that IP3R activation increases Ca2+ transient duration for a broad range of IP3R properties, but the effect of IP3R activation on Ca2+ transient amplitude is dependent on IP3 concentration. Furthermore we demonstrate that IP3-mediated Ca2+ release in the cytosol increases the duty cycle of the Ca2+ transient, the fraction of the cycle for which [Ca2+] is elevated, across a broad range of parameter values and IP3 concentrations. When coupled to a model of downstream transcription factor (NFAT) activation, we demonstrate that there is a high correspondence between the Ca2+ transient duty cycle and the proportion of activated NFAT in the nucleus. These findings suggest increased cytosolic Ca2+ duty cycle as a plausible mechanism for IP3-dependent hypertrophic signaling via Ca2+-sensitive transcription factors such as NFAT in ventricular cardiomyocytes.
Collapse
|
13
|
Guidry ME, Nickerson DP, Crampin EJ, Nash MP, Loiselle DS, Tran K. Insights From Computational Modeling Into the Contribution of Mechano-Calcium Feedback on the Cardiac End-Systolic Force-Length Relationship. Front Physiol 2020; 11:587. [PMID: 32547426 PMCID: PMC7273927 DOI: 10.3389/fphys.2020.00587] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Accepted: 05/11/2020] [Indexed: 11/23/2022] Open
Abstract
In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca2+ transients produced under equivalent force conditions. We found that the duration of the simulated Ca2+ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca2+ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca2+ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca2+ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.
Collapse
Affiliation(s)
- Megan E Guidry
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - David P Nickerson
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Edmund J Crampin
- Systems Biology Laboratory, School of Mathematics and Statistics, The University of Melbourne, Melbourne, VIC, Australia.,ARC Centre of Excellence in Convergent Bio-Nano Science and Technology, School of Chemical and Biomedical Engineering, The University of Melbourne, Melbourne, VIC, Australia
| | - Martyn P Nash
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.,Department of Engineering Science, The University of Auckland, Auckland, New Zealand
| | - Denis S Loiselle
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.,Department of Physiology, The University of Auckland, Auckland, New Zealand
| | - Kenneth Tran
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
14
|
Khokhlova A, Konovalov P, Iribe G, Solovyova O, Katsnelson L. The Effects of Mechanical Preload on Transmural Differences in Mechano-Calcium-Electric Feedback in Single Cardiomyocytes: Experiments and Mathematical Models. Front Physiol 2020; 11:171. [PMID: 32256377 PMCID: PMC7091561 DOI: 10.3389/fphys.2020.00171] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Accepted: 02/13/2020] [Indexed: 11/13/2022] Open
Abstract
Transmural differences in ventricular myocardium are maintained by electromechanical coupling and mechano-calcium/mechano-electric feedback. In the present study, we experimentally investigated the influence of preload on the force characteristics of subendocardial (Endo) and subepicardial (Epi) single ventricular cardiomyocytes stretched by up to 20% from slack sarcomere length (SL) and analyzed the results with the help of mathematical modeling. Mathematical models of Endo and Epi cells, which accounted for regional heterogeneity in ionic currents, Ca2+ handling, and myofilament contractile mechanisms, showed that a greater slope of the active tension–length relationship observed experimentally in Endo cardiomyocytes could be explained by greater length-dependent Ca2+ activation in Endo cells compared with Epi ones. The models also predicted that greater length dependence of Ca2+ activation in Endo cells compared to Epi ones underlies, via mechano-calcium-electric feedback, the reduction in the transmural gradient in action potential duration (APD) at a higher preload. However, the models were unable to reproduce the experimental data on a decrease of the transmural gradient in the time to peak contraction between Endo and Epi cells at longer end-diastolic SL. We hypothesize that preload-dependent changes in viscosity should be involved alongside the Frank–Starling effects to regulate the transmural gradient in length-dependent changes in the time course of contraction of Endo and Epi cardiomyocytes. Our experimental data and their analysis based on mathematical modeling give reason to believe that mechano-calcium-electric feedback plays a critical role in the modulation of electrophysiological and contractile properties of myocytes across the ventricular wall.
Collapse
Affiliation(s)
- Anastasia Khokhlova
- Institute of Immunology and Physiology, Russian Academy of Sciences, Yekaterinburg, Russia.,Institute of Natural Sciences and Mathematics, Ural Federal University, Yekaterinburg, Russia
| | - Pavel Konovalov
- Institute of Immunology and Physiology, Russian Academy of Sciences, Yekaterinburg, Russia
| | - Gentaro Iribe
- Department of Physiology, Asahikawa Medical University, Hokkaido, Japan.,Department of Cardiovascular Physiology, Okayama University, Okayama, Japan
| | - Olga Solovyova
- Institute of Immunology and Physiology, Russian Academy of Sciences, Yekaterinburg, Russia.,Institute of Natural Sciences and Mathematics, Ural Federal University, Yekaterinburg, Russia
| | - Leonid Katsnelson
- Institute of Immunology and Physiology, Russian Academy of Sciences, Yekaterinburg, Russia.,Institute of Natural Sciences and Mathematics, Ural Federal University, Yekaterinburg, Russia
| |
Collapse
|
15
|
Tomek J, Bueno-Orovio A, Passini E, Zhou X, Minchole A, Britton O, Bartolucci C, Severi S, Shrier A, Virag L, Varro A, Rodriguez B. Development, calibration, and validation of a novel human ventricular myocyte model in health, disease, and drug block. eLife 2019; 8:48890. [PMID: 31868580 PMCID: PMC6970534 DOI: 10.7554/elife.48890] [Citation(s) in RCA: 100] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2019] [Accepted: 12/18/2019] [Indexed: 12/19/2022] Open
Abstract
Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.
Collapse
Affiliation(s)
- Jakub Tomek
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Alfonso Bueno-Orovio
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Elisa Passini
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Xin Zhou
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Ana Minchole
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Oliver Britton
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| | - Chiara Bartolucci
- Department of Electrical, Electronic, and Information Engineering "Guglielmo Marconi", University of Bologna, Bologna, Italy
| | - Stefano Severi
- Department of Electrical, Electronic, and Information Engineering "Guglielmo Marconi", University of Bologna, Bologna, Italy
| | - Alvin Shrier
- Department of Physiology, McGill University, Montreal, Canada
| | - Laszlo Virag
- Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Andras Varro
- Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Blanca Rodriguez
- Department of Computer Science, British Heart Foundation Centre of Research Excellence, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
16
|
Umehara S, Tan X, Okamoto Y, Ono K, Noma A, Amano A, Himeno Y. Mechanisms Underlying Spontaneous Action Potential Generation Induced by Catecholamine in Pulmonary Vein Cardiomyocytes: A Simulation Study. Int J Mol Sci 2019; 20:ijms20122913. [PMID: 31207916 PMCID: PMC6628582 DOI: 10.3390/ijms20122913] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2019] [Revised: 06/11/2019] [Accepted: 06/12/2019] [Indexed: 12/24/2022] Open
Abstract
Cardiomyocytes and myocardial sleeves dissociated from pulmonary veins (PVs) potentially generate ectopic automaticity in response to noradrenaline (NA), and thereby trigger atrial fibrillation. We developed a mathematical model of rat PV cardiomyocytes (PVC) based on experimental data that incorporates the microscopic framework of the local control theory of Ca2+ release from the sarcoplasmic reticulum (SR), which can generate rhythmic Ca2+ release (limit cycle revealed by the bifurcation analysis) when total Ca2+ within the cell increased. Ca2+ overload in SR increased resting Ca2+ efflux through the type II inositol 1,4,5-trisphosphate (IP3) receptors (InsP3R) as well as ryanodine receptors (RyRs), which finally triggered massive Ca2+ release through activation of RyRs via local Ca2+ accumulation in the vicinity of RyRs. The new PVC model exhibited a resting potential of −68 mV. Under NA effects, repetitive Ca2+ release from SR triggered spontaneous action potentials (APs) by evoking transient depolarizations (TDs) through Na+/Ca2+ exchanger (APTDs). Marked and variable latencies initiating APTDs could be explained by the time courses of the α1- and β1-adrenergic influence on the regulation of intracellular Ca2+ content and random occurrences of spontaneous TD activating the first APTD. Positive and negative feedback relations were clarified under APTD generation.
Collapse
Affiliation(s)
- Shohei Umehara
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan.
| | - Xiaoqiu Tan
- Institute of Cardiovascular Research, Southwest Medical University, Luzhou 640000, China.
| | - Yosuke Okamoto
- Department of Cell Physiology, Graduate School of Medicine, Akita University, Akita 010-8543, Japan.
| | - Kyoichi Ono
- Department of Cell Physiology, Graduate School of Medicine, Akita University, Akita 010-8543, Japan.
| | - Akinori Noma
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan.
| | - Akira Amano
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan.
| | - Yukiko Himeno
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan.
| |
Collapse
|
17
|
Muangkram Y, Noma A, Amano A. A new myofilament contraction model with ATP consumption for ventricular cell model. J Physiol Sci 2018; 68:541-554. [PMID: 28770433 PMCID: PMC10717283 DOI: 10.1007/s12576-017-0560-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2017] [Accepted: 07/14/2017] [Indexed: 01/14/2023]
Abstract
A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.
Collapse
Affiliation(s)
- Yuttamol Muangkram
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, 525-8577, Japan
| | - Akinori Noma
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, 525-8577, Japan
| | - Akira Amano
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, 525-8577, Japan.
| |
Collapse
|
18
|
Starý T, Biktashev VN. Fast-slow asymptotics for a Markov chain model of fast sodium current. CHAOS (WOODBURY, N.Y.) 2017; 27:093937. [PMID: 28964114 DOI: 10.1063/1.5003013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We explore the feasibility of using fast-slow asymptotics to eliminate the computational stiffness of discrete-state, continuous-time deterministic Markov chain models of ionic channels underlying cardiac excitability. We focus on a Markov chain model of fast sodium current, and investigate its asymptotic behaviour with respect to small parameters identified in different ways.
Collapse
Affiliation(s)
- Tomáš Starý
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Harrison Building, North Park Road, Exeter EX4 4QF, United Kingdom
| | - Vadim N Biktashev
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Harrison Building, North Park Road, Exeter EX4 4QF, United Kingdom
| |
Collapse
|
19
|
Concise Whole-Cell Modeling of BK Ca-CaV Activity Controlled by Local Coupling and Stoichiometry. Biophys J 2017; 112:2387-2396. [PMID: 28591611 DOI: 10.1016/j.bpj.2017.04.035] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2017] [Revised: 04/12/2017] [Accepted: 04/25/2017] [Indexed: 11/23/2022] Open
Abstract
Large-conductance Ca2+-dependent K+ (BKCa) channels are important regulators of electrical activity. These channels colocalize and form ion channel complexes with voltage-dependent Ca2+ (CaV) channels. Recent stochastic simulations of the BKCa-CaV complex with 1:1 stoichiometry have given important insight into the local control of BKCa channels by fluctuating nanodomains of Ca2+. However, such Monte Carlo simulations are computationally expensive, and are therefore not suitable for large-scale simulations of cellular electrical activity. In this work we extend the stochastic model to more realistic BKCa-CaV complexes with 1:n stoichiometry, and analyze the single-complex model with Markov chain theory. From the description of a single BKCa-CaV complex, using arguments based on timescale analysis, we derive a concise model of whole-cell BKCa currents, which can readily be analyzed and inserted into models of cellular electrical activity. We illustrate the usefulness of our results by inserting our BKCa description into previously published whole-cell models, and perform simulations of electrical activity in various cell types, which show that BKCa-CaV stoichiometry can affect whole-cell behavior substantially. Our work provides a simple formulation for the whole-cell BKCa current that respects local interactions in BKCa-CaV complexes, and indicates how local-global coupling of ion channels may affect cell behavior.
Collapse
|
20
|
Calcium dynamics in cardiac excitatory and non-excitatory cells and the role of gap junction. Math Biosci 2017; 289:51-68. [PMID: 28457965 DOI: 10.1016/j.mbs.2017.04.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2016] [Revised: 11/12/2016] [Accepted: 04/26/2017] [Indexed: 11/21/2022]
Abstract
Calcium ions aid in the generation of action potential in myocytes and are responsible for the excitation-contraction coupling of heart. The heart muscle has specialized patches of cells, called excitatory cells (EC) such as the Sino-atrial node cells capable of auto-generation of action potential and cells which receive signals from the excitatory cells, called non-excitatory cells (NEC) such as cells of the ventricular and auricular walls. In order to understand cardiac calcium homeostasis, it is, therefore, important to study the calcium dynamics taking into account both types of cardiac cells. Here we have developed a model to capture the calcium dynamics in excitatory and non-excitatory cells taking into consideration the gap junction mediated calcium ion transfer from excitatory cell to non-excitatory cell. Our study revealed that the gap junctional coupling between excitatory and non-excitatory cells plays important role in the calcium dynamics. It is observed that any reduction in the functioning of gap junction may result in abnormal calcium oscillations in NEC, even when the calcium dynamics is normal in EC cell. Sensitivity of gap junction is observed to be independent of the pacing rate and hence a careful monitoring is required to maintain normal cardiomyocyte condition. It also highlights that sarcoplasmic reticulum may not be always able to control the amount of cytoplasmic calcium under the condition of calcium overload.
Collapse
|
21
|
Maleckar MM, Edwards AG, Louch WE, Lines GT. Studying dyadic structure-function relationships: a review of current modeling approaches and new insights into Ca 2+ (mis)handling. CLINICAL MEDICINE INSIGHTS-CARDIOLOGY 2017; 11:1179546817698602. [PMID: 28469494 PMCID: PMC5392018 DOI: 10.1177/1179546817698602] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Accepted: 12/19/2016] [Indexed: 11/25/2022]
Abstract
Excitation–contraction coupling in cardiac myocytes requires calcium influx through L-type calcium channels in the sarcolemma, which gates calcium release through sarcoplasmic reticulum ryanodine receptors in a process known as calcium-induced calcium release, producing a myoplasmic calcium transient and enabling cardiomyocyte contraction. The spatio-temporal dynamics of calcium release, buffering, and reuptake into the sarcoplasmic reticulum play a central role in excitation–contraction coupling in both normal and diseased cardiac myocytes. However, further quantitative understanding of these cells’ calcium machinery and the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease requires accurate knowledge of cardiac ultrastructure, protein distribution and subcellular function. As current imaging techniques are limited in spatial resolution, limiting insight into changes in calcium handling, computational models of excitation–contraction coupling have been increasingly employed to probe these structure–function relationships. This review will focus on the development of structural models of cardiac calcium dynamics at the subcellular level, orienting the reader broadly towards the development of models of subcellular calcium handling in cardiomyocytes. Specific focus will be given to progress in recent years in terms of multi-scale modeling employing resolved spatial models of subcellular calcium machinery. A review of the state-of-the-art will be followed by a review of emergent insights into calcium-dependent etiologies in heart disease and, finally, we will offer a perspective on future directions for related computational modeling and simulation efforts.
Collapse
Affiliation(s)
- Mary M Maleckar
- Simula Research Laboratory, Center for Cardiological Innovation and Center for Biomedical Computing, Lysaker, Norway
| | - Andrew G Edwards
- Simula Research Laboratory, Center for Cardiological Innovation and Center for Biomedical Computing, Lysaker, Norway.,University of Oslo, Oslo, Norway
| | - William E Louch
- Institute for Experimental Medical Research (IEMR), Oslo University Hospital and the University of Oslo, Oslo, Norway
| | - Glenn T Lines
- Simula Research Laboratory, Center for Cardiological Innovation and Center for Biomedical Computing, Lysaker, Norway
| |
Collapse
|
22
|
Gattoni S, Røe ÅT, Frisk M, Louch WE, Niederer SA, Smith NP. The calcium-frequency response in the rat ventricular myocyte: an experimental and modelling study. J Physiol 2016; 594:4193-224. [PMID: 26916026 DOI: 10.1113/jp272011] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2015] [Accepted: 02/22/2016] [Indexed: 12/11/2022] Open
Abstract
KEY POINTS In the majority of species, including humans, increased heart rate increases cardiac contractility. This change is known as the force-frequency response (FFR). The majority of mammals have a positive force-frequency relationship (FFR). In rat the FFR is controversial. We derive a species- and temperature-specific data-driven model of the rat ventricular myocyte. As a measure of the FFR, we test the effects of changes in frequency and extracellular calcium on the calcium-frequency response (CFR) in our model and three altered models. The results show a biphasic peak calcium-frequency response, due to biphasic behaviour of the ryanodine receptor and the combined effect of the rapid calmodulin buffer and the frequency-dependent increase in diastolic calcium. Alterations to the model reveal that inclusion of Ca(2+) /calmodulin-dependent protein kinase II (CAMKII)-mediated L-type channel and transient outward K(+) current activity enhances the positive magnitude calcium-frequency response, and the absence of CAMKII-mediated increase in activity of the sarco/endoplasmic reticulum Ca(2+) -ATPase induces a negative magnitude calcium-frequency response. ABSTRACT An increase in heart rate affects the strength of cardiac contraction by altering the Ca(2+) transient as a response to physiological demands. This is described by the force-frequency response (FFR), a change in developed force with pacing frequency. The majority of mammals, including humans, have a positive FFR, and cardiac contraction strength increases with heart rate. However, the rat and mouse are exceptions, with the majority of studies reporting a negative FFR, while others report either a biphasic or a positive FFR. Understanding the differences in the FFR between humans and rats is fundamental to interpreting rat-based experimental findings in the context of human physiology. We have developed a novel model of rat ventricular electrophysiology and calcium dynamics, derived predominantly from experimental data recorded under physiological conditions. As a measure of FFR, we tested the effects of changes in stimulation frequency and extracellular calcium concentration on the simulated Ca(2+) transient characteristics and showed a biphasic peak calcium-frequency relationship, consistent with recent observations of a shift from negative to positive FFR when approaching the rat physiological frequency range. We tested the hypotheses that (1) inhibition of Ca(2+) /calmodulin-dependent protein kinase II (CAMKII)-mediated increase in sarco/endoplasmic reticulum Ca(2+) -ATPase (SERCA) activity, (2) CAMKII modulation of SERCA, L-type channel and transient outward K(+) current activity and (3) Na(+) /K(+) pump dynamics play a significant role in the rat FFR. The results reveal a major role for CAMKII modulation of SERCA in the peak Ca(2+) -frequency response, driven most significantly by the cytosolic calcium buffering system and changes in diastolic Ca(2+) .
Collapse
Affiliation(s)
- Sara Gattoni
- King's College London, Department of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, London, UK
| | - Åsmund Treu Røe
- Institute for Experimental Medical Research, Oslo University Hospital and University of Oslo, Oslo, Norway.,K. G. Jebsen Cardiac Research Center and Center for Heart Failure Research, University of Oslo, Oslo, Norway
| | - Michael Frisk
- Institute for Experimental Medical Research, Oslo University Hospital and University of Oslo, Oslo, Norway.,K. G. Jebsen Cardiac Research Center and Center for Heart Failure Research, University of Oslo, Oslo, Norway
| | - William E Louch
- Institute for Experimental Medical Research, Oslo University Hospital and University of Oslo, Oslo, Norway.,K. G. Jebsen Cardiac Research Center and Center for Heart Failure Research, University of Oslo, Oslo, Norway
| | - Steven A Niederer
- King's College London, Department of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, London, UK
| | - Nicolas P Smith
- King's College London, Department of Biomedical Engineering and Imaging Sciences, St. Thomas' Hospital, London, UK.,University of Auckland, Engineering School Block 1, Level 5, 20 Symonds St, Auckland, 101, New Zealand
| |
Collapse
|
23
|
Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention. Ann Biomed Eng 2016; 44:2642-60. [PMID: 27138523 DOI: 10.1007/s10439-016-1628-0] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2015] [Accepted: 04/22/2016] [Indexed: 12/19/2022]
Abstract
Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.
Collapse
|
24
|
Himeno Y, Asakura K, Cha CY, Memida H, Powell T, Amano A, Noma A. A human ventricular myocyte model with a refined representation of excitation-contraction coupling. Biophys J 2016. [PMID: 26200878 DOI: 10.1016/j.bpj.2015.06.017] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Cardiac Ca(2+)-induced Ca(2+) release (CICR) occurs by a regenerative activation of ryanodine receptors (RyRs) within each Ca(2+)-releasing unit, triggered by the activation of L-type Ca(2+) channels (LCCs). CICR is then terminated, most probably by depletion of Ca(2+) in the junctional sarcoplasmic reticulum (SR). Hinch et al. previously developed a tightly coupled LCC-RyR mathematical model, known as the Hinch model, that enables simulations to deal with a variety of functional states of whole-cell populations of a Ca(2+)-releasing unit using a personal computer. In this study, we developed a membrane excitation-contraction model of the human ventricular myocyte, which we call the human ventricular cell (HuVEC) model. This model is a hybrid of the most recent HuVEC models and the Hinch model. We modified the Hinch model to reproduce the regenerative activation and termination of CICR. In particular, we removed the inactivated RyR state and separated the single step of RyR activation by LCCs into triggering and regenerative steps. More importantly, we included the experimental measurement of a transient rise in Ca(2+) concentrations ([Ca(2+)], 10-15 μM) during CICR in the vicinity of Ca(2+)-releasing sites, and thereby calculated the effects of the local Ca(2+) gradient on CICR as well as membrane excitation. This HuVEC model successfully reconstructed both membrane excitation and key properties of CICR. The time course of CICR evoked by an action potential was accounted for by autonomous changes in an instantaneous equilibrium open probability of couplons. This autonomous time course was driven by a core feedback loop including the pivotal local [Ca(2+)], influenced by a time-dependent decay in the SR Ca(2+) content during CICR.
Collapse
Affiliation(s)
- Yukiko Himeno
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan
| | - Keiichi Asakura
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan; Nippon Shinyaku Co., Ltd., Kyoto, Japan
| | - Chae Young Cha
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan; Oxford Centre for Diabetes Endocrinology and Metabolism, University of Oxford, Oxford, UK
| | - Hiraku Memida
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan
| | - Trevor Powell
- Department of Pharmacology, University of Oxford, Oxford, UK
| | - Akira Amano
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan
| | - Akinori Noma
- Biosimulation Research Center, College of Life Sciences, Ritsumeikan University, Shiga, Japan.
| |
Collapse
|
25
|
Wescott AP, Jafri MS, Lederer WJ, Williams GSB. Ryanodine receptor sensitivity governs the stability and synchrony of local calcium release during cardiac excitation-contraction coupling. J Mol Cell Cardiol 2016; 92:82-92. [PMID: 26827896 PMCID: PMC4807626 DOI: 10.1016/j.yjmcc.2016.01.024] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Revised: 01/13/2016] [Accepted: 01/27/2016] [Indexed: 11/27/2022]
Abstract
Calcium-induced calcium release is the principal mechanism that triggers the cell-wide [Ca(2+)]i transient that activates muscle contraction during cardiac excitation-contraction coupling (ECC). Here, we characterize this process in mouse cardiac myocytes with a novel mathematical action potential (AP) model that incorporates realistic stochastic gating of voltage-dependent L-type calcium (Ca(2+)) channels (LCCs) and sarcoplasmic reticulum (SR) Ca(2+) release channels (the ryanodine receptors, RyR2s). Depolarization of the sarcolemma during an AP stochastically activates the LCCs elevating subspace [Ca(2+)] within each of the cell's 20,000 independent calcium release units (CRUs) to trigger local RyR2 opening and initiate Ca(2+) sparks, the fundamental unit of triggered Ca(2+) release. Synchronization of Ca(2+) sparks during systole depends on the nearly uniform cellular activation of LCCs and the likelihood of local LCC openings triggering local Ca(2+) sparks (ECC fidelity). The detailed design and true SR Ca(2+) pump/leak balance displayed by our model permits investigation of ECC fidelity and Ca(2+) spark fidelity, the balance between visible (Ca(2+) spark) and invisible (Ca(2+) quark/sub-spark) SR Ca(2+) release events. Excess SR Ca(2+) leak is examined as a disease mechanism in the context of "catecholaminergic polymorphic ventricular tachycardia (CPVT)", a Ca(2+)-dependent arrhythmia. We find that that RyR2s (and therefore Ca(2+) sparks) are relatively insensitive to LCC openings across a wide range of membrane potentials; and that key differences exist between Ca(2+) sparks evoked during quiescence, diastole, and systole. The enhanced RyR2 [Ca(2+)]i sensitivity during CPVT leads to increased Ca(2+) spark fidelity resulting in asynchronous systolic Ca(2+) spark activity. It also produces increased diastolic SR Ca(2+) leak with some prolonged Ca(2+) sparks that at times become "metastable" and fail to efficiently terminate. There is a huge margin of safety for stable Ca(2+) handling within the cell and this novel mechanistic model provides insight into the molecular signaling characteristics that help maintain overall Ca(2+) stability even under the conditions of high SR Ca(2+) leak during CPVT. Finally, this model should provide tools for investigators to examine normal and pathological Ca(2+) signaling characteristics in the heart.
Collapse
Affiliation(s)
- Andrew P Wescott
- Center for Biomedical Engineering and Technology & Department of Physiology, University of Maryland, Baltimore, Baltimore, MD, United States
| | - M Saleet Jafri
- Center for Biomedical Engineering and Technology & Department of Physiology, University of Maryland, Baltimore, Baltimore, MD, United States; Molecular Neuroscience Department, George Mason University, Fairfax, VA, United States
| | - W J Lederer
- Center for Biomedical Engineering and Technology & Department of Physiology, University of Maryland, Baltimore, Baltimore, MD, United States
| | - George S B Williams
- Center for Biomedical Engineering and Technology & Department of Physiology, University of Maryland, Baltimore, Baltimore, MD, United States.
| |
Collapse
|
26
|
Song Z, Karma A, Weiss JN, Qu Z. Long-Lasting Sparks: Multi-Metastability and Release Competition in the Calcium Release Unit Network. PLoS Comput Biol 2016; 12:e1004671. [PMID: 26730593 PMCID: PMC4701461 DOI: 10.1371/journal.pcbi.1004671] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2015] [Accepted: 11/23/2015] [Indexed: 11/20/2022] Open
Abstract
Calcium (Ca) sparks are elementary events of biological Ca signaling. A normal Ca spark has a brief duration in the range of 10 to 100 ms, but long-lasting sparks with durations of several hundred milliseconds to seconds are also widely observed. Experiments have shown that the transition from normal to long-lasting sparks can occur when ryanodine receptor (RyR) open probability is either increased or decreased. Here, we demonstrate theoretically and computationally that long-lasting sparks emerge as a collective dynamical behavior of the network of diffusively coupled Ca release units (CRUs). We show that normal sparks occur when the CRU network is monostable and excitable, while long-lasting sparks occur when the network dynamics possesses multiple metastable attractors, each attractor corresponding to a different spatial firing pattern of sparks. We further highlight the mechanisms and conditions that produce long-lasting sparks, demonstrating the existence of an optimal range of RyR open probability favoring long-lasting sparks. We find that when CRU firings are sparse and sarcoplasmic reticulum (SR) Ca load is high, increasing RyR open probability promotes long-lasting sparks by potentiating Ca-induced Ca release (CICR). In contrast, when CICR is already strong enough to produce frequent firings, decreasing RyR open probability counter-intuitively promotes long-lasting sparks by decreasing spark frequency. The decrease in spark frequency promotes intra-SR Ca diffusion from neighboring non-firing CRUs to the firing CRUs, which helps to maintain the local SR Ca concentration of the firing CRUs above a critical level to sustain firing. In this setting, decreasing RyR open probability further suppresses long-lasting sparks by weakening CICR. Since a long-lasting spark terminates via the Kramers’ escape process over a potential barrier, its duration exhibits an exponential distribution determined by the barrier height and noise strength, which is modulated differently by different ways of altering the Ca release flux strength. Calcium (Ca) sparks, resulting from Ca-induced Ca release, are elementary events of biological Ca signaling. Sparks are normally brief, but long-lasting sparks have been widely observed experimentally under various conditions. The underlying mechanisms of spark duration or termination and the corresponding determinants remain a topic of debate. In this study, we demonstrate theoretically and computationally that normal brief sparks are excitable transients, while long-lasting sparks are multiple metastable states emerging in the diffusively coupled Ca release unit network, as a result of cooperativity and release competition among the Ca release units. Termination of a long-lasting spark is a Kramers’ escape process over a potential barrier, and the spark duration is the first-passage time, exhibiting an exponential distribution.
Collapse
Affiliation(s)
- Zhen Song
- The UCLA Cardiovascular Research Laboratory, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Alain Karma
- Department of Physics, Northeastern University, Boston, Massachusetts, United States of America
| | - James N. Weiss
- The UCLA Cardiovascular Research Laboratory, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
- Department of Physiology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Zhilin Qu
- The UCLA Cardiovascular Research Laboratory, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
- Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, California, United States of America
- * E-mail:
| |
Collapse
|
27
|
Semantics-Based Composition of Integrated Cardiomyocyte Models Motivated by Real-World Use Cases. PLoS One 2015; 10:e0145621. [PMID: 26716837 PMCID: PMC4696653 DOI: 10.1371/journal.pone.0145621] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Accepted: 11/06/2015] [Indexed: 11/19/2022] Open
Abstract
Semantics-based model composition is an approach for generating complex biosimulation models from existing components that relies on capturing the biological meaning of model elements in a machine-readable fashion. This approach allows the user to work at the biological rather than computational level of abstraction and helps minimize the amount of manual effort required for model composition. To support this compositional approach, we have developed the SemGen software, and here report on SemGen's semantics-based merging capabilities using real-world modeling use cases. We successfully reproduced a large, manually-encoded, multi-model merge: the "Pandit-Hinch-Niederer" (PHN) cardiomyocyte excitation-contraction model, previously developed using CellML. We describe our approach for annotating the three component models used in the PHN composition and for merging them at the biological level of abstraction within SemGen. We demonstrate that we were able to reproduce the original PHN model results in a semi-automated, semantics-based fashion and also rapidly generate a second, novel cardiomyocyte model composed using an alternative, independently-developed tension generation component. We discuss the time-saving features of our compositional approach in the context of these merging exercises, the limitations we encountered, and potential solutions for enhancing the approach.
Collapse
|
28
|
Winslow RL, Walker MA, Greenstein JL. Modeling calcium regulation of contraction, energetics, signaling, and transcription in the cardiac myocyte. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2015; 8:37-67. [PMID: 26562359 DOI: 10.1002/wsbm.1322] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Revised: 09/29/2015] [Accepted: 09/30/2015] [Indexed: 12/11/2022]
Abstract
Calcium (Ca(2+)) plays many important regulatory roles in cardiac muscle cells. In the initial phase of the action potential, influx of Ca(2+) through sarcolemmal voltage-gated L-type Ca(2+) channels (LCCs) acts as a feed-forward signal that triggers a large release of Ca(2+) from the junctional sarcoplasmic reticulum (SR). This Ca(2+) drives heart muscle contraction and pumping of blood in a process known as excitation-contraction coupling (ECC). Triggered and released Ca(2+) also feed back to inactivate LCCs, attenuating the triggered Ca(2+) signal once release has been achieved. The process of ECC consumes large amounts of ATP. It is now clear that in a process known as excitation-energetics coupling, Ca(2+) signals exert beat-to-beat regulation of mitochondrial ATP production that closely couples energy production with demand. This occurs through transport of Ca(2+) into mitochondria, where it regulates enzymes of the tricarboxylic acid cycle. In excitation-signaling coupling, Ca(2+) activates a number of signaling pathways in a feed-forward manner. Through effects on their target proteins, these interconnected pathways regulate Ca(2+) signals in complex ways to control electrical excitability and contractility of heart muscle. In a process known as excitation-transcription coupling, Ca(2+) acting primarily through signal transduction pathways also regulates the process of gene transcription. Because of these diverse and complex roles, experimentally based mechanistic computational models are proving to be very useful for understanding Ca(2+) signaling in the cardiac myocyte.
Collapse
Affiliation(s)
- Raimond L Winslow
- Institute for Computational Medicine and Department of Biomedical Engineering, The Johns Hopkins University School of Medicine and Whiting School of Engineering, Baltimore, MD, USA
| | - Mark A Walker
- Institute for Computational Medicine and Department of Biomedical Engineering, The Johns Hopkins University School of Medicine and Whiting School of Engineering, Baltimore, MD, USA
| | - Joseph L Greenstein
- Institute for Computational Medicine and Department of Biomedical Engineering, The Johns Hopkins University School of Medicine and Whiting School of Engineering, Baltimore, MD, USA
| |
Collapse
|
29
|
Bugenhagen SM, Beard DA. Computational analysis of the regulation of Ca(2+) dynamics in rat ventricular myocytes. Phys Biol 2015; 12:056008. [PMID: 26358004 DOI: 10.1088/1478-3975/12/5/056008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Force-frequency relationships of isolated cardiac myocytes show complex behaviors that are thought to be specific to both the species and the conditions associated with the experimental preparation. Ca(2+) signaling plays an important role in shaping the force-frequency relationship, and understanding the properties of the force-frequency relationship in vivo requires an understanding of Ca(2+) dynamics under physiologically relevant conditions. Ca(2+) signaling is itself a complicated process that is best understood on a quantitative level via biophysically based computational simulation. Although a large number of models are available in the literature, the models are often a conglomeration of components parameterized to data of incompatible species and/or experimental conditions. In addition, few models account for modulation of Ca(2+) dynamics via β-adrenergic and calmodulin-dependent protein kinase II (CaMKII) signaling pathways even though they are hypothesized to play an important regulatory role in vivo. Both protein-kinase-A and CaMKII are known to phosphorylate a variety of targets known to be involved in Ca(2+) signaling, but the effects of these pathways on the frequency- and inotrope-dependence of Ca(2+) dynamics are not currently well understood. In order to better understand Ca(2+) dynamics under physiological conditions relevant to rat, a previous computational model is adapted and re-parameterized to a self-consistent dataset obtained under physiological temperature and pacing frequency and updated to include β-adrenergic and CaMKII regulatory pathways. The necessity of specific effector mechanisms of these pathways in capturing inotrope- and frequency-dependence of the data is tested by attempting to fit the data while including and/or excluding those effector components. We find that: (1) β-adrenergic-mediated phosphorylation of the L-type calcium channel (LCC) (and not of phospholamban (PLB)) is sufficient to explain the inotrope-dependence; and (2) that CaMKII-mediated regulation of neither the LCC nor of PLB is required to explain the frequency-dependence of the data.
Collapse
Affiliation(s)
- Scott M Bugenhagen
- Department of Physiology, Medical College of Wisconsin, Milwaukee, WI 53226, USA
| | | |
Collapse
|
30
|
Aronsen JM, Skogestad J, Lewalle A, Louch WE, Hougen K, Stokke MK, Swift F, Niederer S, Smith NP, Sejersted OM, Sjaastad I. Hypokalaemia induces Ca²⁺ overload and Ca²⁺ waves in ventricular myocytes by reducing Na⁺,K⁺-ATPase α₂ activity. J Physiol 2014; 593:1509-21. [PMID: 25772299 DOI: 10.1113/jphysiol.2014.279893] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2014] [Accepted: 10/02/2014] [Indexed: 11/08/2022] Open
Abstract
KEY POINTS Hypokalaemia is a risk factor for development of ventricular arrhythmias. In rat ventricular myocytes, low extracellular K(+) (corresponding to clinical moderate hypokalaemia) increased Ca(2+) wave probability, Ca(2+) transient amplitude, sarcoplasmic reticulum (SR) Ca(2+) load and induced SR Ca(2+) leak. Low extracellular K(+) reduced Na(+),K(+)-ATPase (NKA) activity and hyperpolarized the resting membrane potential in ventricular myocytes. Both experimental data and modelling indicate that reduced NKA activity and subsequent Na(+) accumulation sensed by the Na(+), Ca(2+) exchanger (NCX) lead to increased Ca(2+) transient amplitude despite concomitant hyperpolarization of the resting membrane potential. Low extracellular K(+) induced Ca(2+) overload by lowering NKA α2 activity. Triggered ventricular arrhythmias in patients with hypokalaemia may therefore be attributed to reduced NCX forward mode activity linked to an effect on the NKA α2 isoform. ABSTRACT Hypokalaemia is a risk factor for development of ventricular arrhythmias. The aim of this study was to determine the cellular mechanisms leading to triggering of arrhythmias in ventricular myocytes exposed to low Ko. Low Ko, corresponding to moderate hypokalaemia, increased Ca(2+) transient amplitude, sarcoplasmic reticulum (SR) Ca(2+) load, SR Ca(2+) leak and Ca(2+) wave probability in field stimulated rat ventricular myocytes. The mechanisms leading to Ca(2+) overload were examined. Low Ko reduced Na(+),K(+)-ATPase (NKA) currents, increased cytosolic Na(+) concentration and increased the Na(+) level sensed by the Na(+), Ca(2+) exchanger (NCX). Low Ko also hyperpolarized the resting membrane potential (RMP) without significant alterations in action potential duration. Experiments in voltage clamped and field stimulated ventricular myocytes, along with mathematical modelling, suggested that low Ko increases the Ca(2+) transient amplitude by reducing NKA activity despite hyperpolarization of the RMP. Selective inhibition of the NKA α2 isoform by low dose ouabain abolished the ability of low Ko to reduce NKA currents, to increase Na(+) levels sensed by NCX and to increase the Ca(2+) transient amplitude. We conclude that low Ko, within the range of moderate hypokalaemia, increases Ca(2+) levels in ventricular myocytes by reducing the pumping rate of the NKA α2 isoform with subsequent Na(+) accumulation sensed by the NCX. These data highlight reduced NKA α2 -mediated control of NCX activity as a possible mechanism underlying triggered ventricular arrhythmias in patients with hypokalaemia.
Collapse
Affiliation(s)
- J M Aronsen
- Institute for Experimental Medical Research, Oslo University Hospital Ullevål and University of Oslo, Oslo, Norway; Bjørknes College, Oslo, Norway
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
31
|
EAD and DAD mechanisms analyzed by developing a new human ventricular cell model. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2014; 116:11-24. [PMID: 25192800 DOI: 10.1016/j.pbiomolbio.2014.08.008] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2014] [Accepted: 08/09/2014] [Indexed: 12/31/2022]
Abstract
It has long been suggested that the Ca(2+)-mechanisms are largely involved in generating the early afterdepolarization (EAD) as well as the delayed afterdepolarization (DAD). This view was examined in a quantitative manner by applying the lead potential analysis to a new human ventricular cell model. In this ventricular cell model, the tight coupled LCC-RyR model (CaRU) based on local control theory (Hinch et al. 2004) and ion channel models mostly based on human electrophysiological data were included to reproduce realistic Ca(2+) dynamics as well as the membrane excitation. Simultaneously, the Ca(2+) accumulation near the Ca(2+) releasing site was incorporated as observed in real cardiac myocytes. The maximum rate of ventricular repolarization (-1.02 mV/ms) is due to IK1 (-0.55 mV/ms) and the rest is provided nearly equally by INCX (-0.20 mV/ms), INaL (-0.16 mV/ms) and INaT (-0.13 mV/ms). These INaL and INaT components are due to closure of the voltage gate, which remains partially open during the plateau potential. DADs could be evoked by applying high-frequency stimulations supplemented by a partial Na(+)/K(+) pump inhibition, or by a microinjection of Ca(2+). EADs was evoked by retarding the inactivation of INaL. The lead potential (VL) analysis revealed that IK1 and IKr played the primary role to reverse the AP repolarization to depolarizing limb of EAD. ICaL and INCX amplified EAD, while the remaining currents partially antagonized dVL/dt. The maximum rate of rise of EAD was attributable to the rapid activation of both ICaL (45.5%) and INCX (54.5%).
Collapse
|
32
|
Maybank PJ, Whiteley JP. Automatic simplification of systems of reaction–diffusion equations by a posteriori analysis. Math Biosci 2014; 248:146-57. [DOI: 10.1016/j.mbs.2013.12.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2013] [Revised: 12/20/2013] [Accepted: 12/31/2013] [Indexed: 01/02/2023]
|
33
|
Liu K, Yao G, Yu Z. Parallel acceleration for modeling of calcium dynamics in cardiac myocytes. Biomed Mater Eng 2013; 24:1417-24. [PMID: 24212039 DOI: 10.3233/bme-130946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Spatial-temporal calcium dynamics due to calcium release, buffering, and re-uptaking plays a central role in studying excitation-contraction (E-C) coupling in both healthy and defected cardiac myocytes. In our previous work, partial differential equations (PDEs) had been used to simulate calcium dynamics with realistic geometries extracted from electron microscopic imaging data. However, the computational costs of such simulations are very high on a single processor. To alleviate this problem, we have accelerated the numerical simulations of calcium dynamics by using graphics processing units (GPUs). Computational performance and simulation accuracy are compared with those based on a single CPU and another popular parallel computing technique, OpenMP.
Collapse
Affiliation(s)
- Ke Liu
- Department of Computer Science, University of Wisconsin-Milwaukee, WI, USA
| | | | | |
Collapse
|
34
|
Shen C, Chen H. Coarse‐graining Calcium Dynamics on Stochastic Reaction‐diffusion Lattice Model. CHINESE J CHEM PHYS 2013. [DOI: 10.1063/1674-0068/26/02/181-184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
35
|
Markhasin VS, Balakin AA, Katsnelson LB, Konovalov P, Lookin ON, Protsenko Y, Solovyova O. Slow force response and auto-regulation of contractility in heterogeneous myocardium. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2012; 110:305-18. [DOI: 10.1016/j.pbiomolbio.2012.08.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2012] [Accepted: 08/09/2012] [Indexed: 11/25/2022]
|
36
|
Lee YS, Liu OZ, Sobie EA. Decoding myocardial Ca²⁺ signals across multiple spatial scales: a role for sensitivity analysis. J Mol Cell Cardiol 2012; 58:92-9. [PMID: 23026728 DOI: 10.1016/j.yjmcc.2012.09.009] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/21/2012] [Accepted: 09/20/2012] [Indexed: 11/28/2022]
Abstract
Numerous studies have employed mathematical modeling to quantitatively understand release of Ca(2+) from the sarcoplasmic reticulum (SR) in the heart. Models have been used to investigate physiologically important phenomena such as triggering of SR Ca(2+) release by Ca(2+) entry across the cell membrane and spontaneous leak of Ca(2+) from the SR in quiescent heart cells. In this review we summarize studies that have modeled myocardial Ca(2+) at different spatial scales: the sub-cellular level, the cellular level, and the multicellular level. We discuss each category of models from the standpoint of parameter sensitivity analysis, a common simulation procedure that can generate quantitative, comprehensive predictions about how changes in conditions influence model output. We propose that this is a useful perspective for conceptualizing models, in part because a sensitivity analysis requires the investigator to define the relevant parameters and model outputs. This procedure therefore helps to illustrate the capabilities and limitations of each model. We further suggest that in future studies, sensitivity analyses will aid in simplifying complex models and in suggesting experiments to differentiate between competing models built with different assumptions. We conclude with a discussion of unresolved questions that are likely to be addressed over the next several years.
Collapse
Affiliation(s)
- Young-Seon Lee
- Pharmacology and Systems Therapeutics, Mount Sinai School of Medicine, New York, NY, USA
| | | | | |
Collapse
|
37
|
Gauthier LD, Greenstein JL, Winslow RL. Toward an integrative computational model of the Guinea pig cardiac myocyte. Front Physiol 2012; 3:244. [PMID: 22783206 PMCID: PMC3389778 DOI: 10.3389/fphys.2012.00244] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2012] [Accepted: 06/14/2012] [Indexed: 11/22/2022] Open
Abstract
The local control theory of excitation-contraction (EC) coupling asserts that regulation of calcium (Ca2+) release occurs at the nanodomain level, where openings of single L-type Ca2+ channels (LCCs) trigger openings of small clusters of ryanodine receptors (RyRs) co-localized within the dyad. A consequence of local control is that the whole-cell Ca2+ transient is a smooth continuous function of influx of Ca2+ through LCCs. While this so-called graded release property has been known for some time, its functional importance to the integrated behavior of the cardiac ventricular myocyte has not been fully appreciated. We previously formulated a biophysically based model, in which LCCs and RyRs interact via a coarse-grained representation of the dyadic space. The model captures key features of local control using a low-dimensional system of ordinary differential equations. Voltage-dependent gain and graded Ca2+ release are emergent properties of this model by virtue of the fact that model formulation is closely based on the sub-cellular basis of local control. In this current work, we have incorporated this graded release model into a prior model of guinea pig ventricular myocyte electrophysiology, metabolism, and isometric force production. The resulting integrative model predicts the experimentally observed causal relationship between action potential (AP) shape and timing of Ca2+ and force transients, a relationship that is not explained by models lacking the graded release property. Model results suggest that even relatively subtle changes in AP morphology that may result, for example, from remodeling of membrane transporter expression in disease or spatial variation in cell properties, may have major impact on the temporal waveform of Ca2+ transients, thus influencing tissue level electromechanical function.
Collapse
Affiliation(s)
- Laura Doyle Gauthier
- Department of Biomedical Engineering, Institute for Computational Medicine, The Johns Hopkins University School of Medicine and Whiting School of Engineering Baltimore, MD, USA
| | | | | |
Collapse
|
38
|
Yang JH, Saucerman JJ. Phospholemman is a negative feed-forward regulator of Ca2+ in β-adrenergic signaling, accelerating β-adrenergic inotropy. J Mol Cell Cardiol 2012; 52:1048-55. [PMID: 22289214 PMCID: PMC3327824 DOI: 10.1016/j.yjmcc.2011.12.015] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/28/2011] [Revised: 11/21/2011] [Accepted: 12/29/2011] [Indexed: 01/20/2023]
Abstract
Sympathetic stimulation enhances cardiac contractility by stimulating β-adrenergic signaling and protein kinase A (PKA). Recently, phospholemman (PLM) has emerged as an important PKA substrate capable of regulating cytosolic Ca(2+) transients. However, it remains unclear how PLM contributes to β-adrenergic inotropy. Here we developed a computational model to clarify PLM's role in the β-adrenergic signaling response. Simulating Na(+) and sarcoplasmic reticulum (SR) Ca(2+) clamps, we identify an effect of PLM phosphorylation on SR unloading as the key mechanism by which PLM confers cytosolic Ca(2+) adaptation to long-term β-adrenergic receptor (β-AR) stimulation. Moreover, we show that phospholamban (PLB) opposes and overtakes these actions on SR load, forming a negative feed-forward loop in the β-adrenergic signaling cascade. This network motif dominates the negative feedback conferred by β-AR desensitization and accelerates β-AR-induced inotropy. Model analysis therefore unmasks key actions of PLM phosphorylation during β-adrenergic signaling, indicating that PLM is a critical component of the fight-or-flight response.
Collapse
Affiliation(s)
- Jason H. Yang
- Department of Biomedical Engineering, University of Virginia; Robert M. Berne Cardiovascular Research Center, University of Virginia
| | - Jeffrey J. Saucerman
- Department of Biomedical Engineering, University of Virginia; Robert M. Berne Cardiovascular Research Center, University of Virginia
| |
Collapse
|
39
|
Hake J, Edwards AG, Yu Z, Kekenes-Huskey PM, Michailova AP, McCammon JA, Holst MJ, Hoshijima M, McCulloch AD. Modelling cardiac calcium sparks in a three-dimensional reconstruction of a calcium release unit. J Physiol 2012; 590:4403-22. [PMID: 22495592 DOI: 10.1113/jphysiol.2012.227926] [Citation(s) in RCA: 74] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Johan Hake
- Department of Bioengineering, University of California San Diego, CA, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|
40
|
Noble D, Garny A, Noble PJ. How the Hodgkin-Huxley equations inspired the Cardiac Physiome Project. J Physiol 2012; 590:2613-28. [PMID: 22473779 DOI: 10.1113/jphysiol.2011.224238] [Citation(s) in RCA: 79] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Early modelling of cardiac cells (1960-1980) was based on extensions of the Hodgkin-Huxley nerve axon equations with additional channels incorporated, but after 1980 it became clear that processes other than ion channel gating were also critical in generating electrical activity. This article reviews the development of models representing almost all cell types in the heart, many different species, and the software tools that have been created to facilitate the cardiac Physiome Project.
Collapse
Affiliation(s)
- Denis Noble
- Department of Physiology, Anatomy & Genetics, University of Oxford, Oxford OX1 3PT, UK.
| | | | | |
Collapse
|
41
|
Cheng Y, Kekenes-Huskey P, Hake J, Holst M, McCammon J, Michailova A. Multi-Scale Continuum Modeling of Biological Processes: From Molecular Electro-Diffusion to Sub-Cellular Signaling Transduction. ACTA ACUST UNITED AC 2012; 5. [PMID: 23505398 DOI: 10.1088/1749-4699/5/1/015002] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
This article provides a brief review of multi-scale modeling at the molecular to cellular scale, with new results for heart muscle cells. A finite element-based simulation package (SMOL) was used to investigate the signaling transduction at molecular and sub-cellular scales (http://mccammon.ucsd.edu/smol/, http://FETK.org) by numerical solution of time-dependent Smoluchowski equations and a reaction-diffusion system. At the molecular scale, SMOL has yielded experimentally-validated estimates of the diffusion-limited association rates for the binding of acetylcholine to mouse acetylcholinesterase using crystallographic structural data. The predicted rate constants exhibit increasingly delayed steady-state times with increasing ionic strength and demonstrate the role of an enzyme's electrostatic potential in influencing ligand binding. At the sub-cellular scale, an extension of SMOL solves a non-linear, reaction-diffusion system describing Ca2+ ligand buffering and diffusion in experimentally-derived rodent ventricular myocyte geometries. Results reveal the important role for mobile and stationary Ca2+ buffers, including Ca2+ indicator dye. We found that the alterations in Ca2+-binding and dissociation rates of troponin C (TnC) and total TnC concentration modulate subcellular Ca2+ signals. Model predicts that reduced off-rate in whole troponin complex (TnC, TnI, TnT) versus reconstructed thin filaments (Tn, Tm, actin) alters cytosolic Ca2+ dynamics under control conditions or in disease-linked TnC mutations. The ultimate goal of these studies is to develop scalable methods and theories for integration of molecular-scale information into simulations of cellular-scale systems.
Collapse
Affiliation(s)
- Y Cheng
- Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA
| | | | | | | | | | | |
Collapse
|
42
|
Noble D. A theory of biological relativity: no privileged level of causation. Interface Focus 2012; 2:55-64. [PMID: 23386960 PMCID: PMC3262309 DOI: 10.1098/rsfs.2011.0067] [Citation(s) in RCA: 172] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2011] [Accepted: 10/17/2011] [Indexed: 12/29/2022] Open
Abstract
Must higher level biological processes always be derivable from lower level data and mechanisms, as assumed by the idea that an organism is completely defined by its genome? Or are higher level properties necessarily also causes of lower level behaviour, involving actions and interactions both ways? This article uses modelling of the heart, and its experimental basis, to show that downward causation is necessary and that this form of causation can be represented as the influences of initial and boundary conditions on the solutions of the differential equations used to represent the lower level processes. These insights are then generalized. A priori, there is no privileged level of causation. The relations between this form of 'biological relativity' and forms of relativity in physics are discussed. Biological relativity can be seen as an extension of the relativity principle by avoiding the assumption that there is a privileged scale at which biological functions are determined.
Collapse
Affiliation(s)
- Denis Noble
- Department of Physiology, Anatomy and Genetics, University of Oxford, Parks Road, Oxford OX1 3PT, UK
| |
Collapse
|
43
|
Yao G, Yu Z. A localized meshless approach for modeling spatial-temporal calcium dynamics in ventricular myocytes. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2012; 28:187-204. [PMID: 22408720 PMCID: PMC3294421 DOI: 10.1002/cnm.1453] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Spatial–temporal calcium dynamics due to calcium release, buffering and re-uptaking plays a central role in studying excitation–contraction (E–C) coupling in both normal and diseased cardiac myocytes. In this paper, we employ a meshless method, namely, the local radial basis function collocation method (LRBFCM), to model such calcium behaviors by solving a nonlinear system of reaction–diffusion partial differential equations. In particular, a simplified structural unit containing a single transverse tubule (T-tubule) and its surrounding half sarcomeres is investigated using the meshless method. Numerical results are compared with those generated by finite element methods, showing the capability and efficiency of the LRBFCM in modeling calcium dynamics in ventricular myocytes. The single T-tubule model is also extended to the whole-cell scale with T-tubules excluded to demonstrate the scalability of the proposed meshless method in handling very large domains. The experiments have shown that the LRBFCM is suitable to multiscale modeling of calcium dynamics in ventricular myocytes with high accuracy and efficiency.
Collapse
Affiliation(s)
| | - Zeyun Yu
- Correspondence to: Department of Computer Science, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA.
| |
Collapse
|
44
|
Tuckwell HC. Quantitative aspects of L-type Ca2+ currents. Prog Neurobiol 2012; 96:1-31. [DOI: 10.1016/j.pneurobio.2011.09.010] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2010] [Revised: 09/16/2011] [Accepted: 09/23/2011] [Indexed: 12/24/2022]
|
45
|
Subramaniam S, Fahy E, Gupta S, Sud M, Byrnes RW, Cotter D, Dinasarapu AR, Maurya MR. Bioinformatics and systems biology of the lipidome. Chem Rev 2011; 111:6452-90. [PMID: 21939287 PMCID: PMC3383319 DOI: 10.1021/cr200295k] [Citation(s) in RCA: 123] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Shankar Subramaniam
- Department of Bioengineering, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
- San Diego Supercomputer Center, 9500 Gilman Drive, La Jolla, California, 92093, USA
- Departments of Chemistry and Biochemistry, and Department of Cellular and Molecular Medicine, University of California at San Diego, La Jolla, California 92093, USA
| | - Eoin Fahy
- Department of Bioengineering, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| | - Shakti Gupta
- Department of Bioengineering, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| | - Manish Sud
- San Diego Supercomputer Center, 9500 Gilman Drive, La Jolla, California, 92093, USA
| | - Robert W. Byrnes
- San Diego Supercomputer Center, 9500 Gilman Drive, La Jolla, California, 92093, USA
| | - Dawn Cotter
- San Diego Supercomputer Center, 9500 Gilman Drive, La Jolla, California, 92093, USA
| | - Ashok Reddy Dinasarapu
- Department of Bioengineering, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| | - Mano Ram Maurya
- Department of Bioengineering, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| |
Collapse
|
46
|
Schendel T, Thul R, Sneyd J, Falcke M. How does the ryanodine receptor in the ventricular myocyte wake up: by a single or by multiple open L-type Ca2+ channels? EUROPEAN BIOPHYSICS JOURNAL: EBJ 2011; 41:27-39. [DOI: 10.1007/s00249-011-0755-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/02/2011] [Revised: 08/23/2011] [Accepted: 09/14/2011] [Indexed: 02/07/2023]
|
47
|
Winslow RL, Greenstein JL. Cardiac myocytes and local signaling in nano-domains. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2011; 107:48-59. [PMID: 21718716 DOI: 10.1016/j.pbiomolbio.2011.06.005] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2011] [Accepted: 06/14/2011] [Indexed: 10/18/2022]
Abstract
It is well known that calcium-induced calcium-release in cardiac myocytes takes place in spatially restricted regions known as dyads, where discrete patches of junctional sarcoplasmic reticulum tightly associate with the t-tubule membrane. The dimensions of a dyad are so small that it contains only a few Ca²⁺ ions at any given time. Ca²⁺ signaling in the dyad is therefore noisy, and dominated by the Brownian motion of Ca²⁺ ions in a potential field. Remarkably, from this complexity emerges the integrated behavior of the myocyte in which, under normal conditions, precise control of Ca²⁺ release and muscle contraction is maintained over the life of the cell. This is but one example of how signal processing within the cardiac myocyte and other cells often occurs in small "nano-domains" where proteins and protein complexes interact at spatial dimensions on the order of ∼1-10 nm and at time-scales on the order of nanoseconds to perform the functions of the cell. In this article, we will review several examples of local signaling in nano-domains, how it contributes to the integrative behavior of the cardiac myocyte, and present computational methods for modeling signal processing within these domains across differing spatio-temporal scales.
Collapse
Affiliation(s)
- Raimond L Winslow
- The Institute for Computational Medicine & Department of Biomedical Engineering, The Johns Hopkins University, School of Medicine & Whiting School of Engineering, Baltimore, MD 21218, USA.
| | | |
Collapse
|
48
|
Katsnelson LB, Solovyova O, Balakin A, Lookin O, Konovalov P, Protsenko Y, Sulman T, Markhasin VS. Contribution of mechanical factors to arrhythmogenesis in calcium overloaded cardiomyocytes: model predictions and experiments. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2011; 107:81-9. [PMID: 21699912 DOI: 10.1016/j.pbiomolbio.2011.06.001] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2011] [Accepted: 06/06/2011] [Indexed: 10/18/2022]
Abstract
It is well-known that Ca²⁺ overload in cardiomyocytes may underlie arrhythmias. However, the possible contribution of mechanical factors to rhythm disturbances in Ca²⁺ overloaded myocytes has not been sufficiently investigated. We used a mathematical model of the electrical and mechanical activity of cardiomyocytes to reveal an essential role of the mechanisms of cardiac mechano-electric feedback in arrhythmogenesis in Ca²⁺ overloaded myocardium. In the model, the following mechanical factors increased Ca²⁺ overload in contracting cardiomyocytes and promoted rhythm disturbances: i) a decrease in the mechanical load for afterloaded contractions; and ii) a decrease in the initial length of sarcomeres for isometric twitches. In exact accordance with the model predictions, in experiments on papillary muscles from the right ventricle of guinea pigs with Ca²⁺ overloaded cardiomyocytes (using 0.5-1 μM of ouabain), we found that emergence of rhythm disturbances and extrasystoles depends on the mechanical conditions of muscle contraction.
Collapse
Affiliation(s)
- Leonid B Katsnelson
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences, 106 Pervomayskaya str, Ekaterinburg 620049, Russian Federation.
| | | | | | | | | | | | | | | |
Collapse
|
49
|
Yu Z, Yao G, Hoshijima M, Michailova A, Holst M. Multiscale modeling of calcium dynamics in ventricular myocytes with realistic transverse tubules. IEEE Trans Biomed Eng 2011; 58:2947-51. [PMID: 21632291 DOI: 10.1109/tbme.2011.2158316] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Zeyun Yu
- Department of Computer Science, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA.
| | | | | | | | | |
Collapse
|
50
|
Abstract
Excitation-contraction coupling in the cardiac myocyte is mediated by a number of highly integrated mechanisms of intracellular Ca²(+) transport. The complexity and integrative nature of heart cell electrophysiology and Ca²(+) cycling has led to an evolution of computational models that have played a crucial role in shaping our understanding of heart function. An important emerging theme in systems biology is that the detailed nature of local signaling events, such as those that occur in the cardiac dyad, have important consequences at higher biological scales. Multiscale modeling techniques have revealed many mechanistic links between microscale events, such as Ca²(+) binding to a channel protein, and macroscale phenomena, such as excitation-contraction coupling gain. Here, we review experimentally based multiscale computational models of excitation-contraction coupling and the insights that have been gained through their application.
Collapse
Affiliation(s)
- Joseph L Greenstein
- Center for Cardiovascular Bioinformatics and Modeling, Whitaker Biomedical Engineering Institute, The Johns Hopkins University, Baltimore, MD, USA
| | | |
Collapse
|