1
|
Eisner D, Neher E, Taschenberger H, Smith G. Physiology of intracellular calcium buffering. Physiol Rev 2023; 103:2767-2845. [PMID: 37326298 DOI: 10.1152/physrev.00042.2022] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 05/08/2023] [Accepted: 06/11/2023] [Indexed: 06/17/2023] Open
Abstract
Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.
Collapse
Affiliation(s)
- David Eisner
- Division of Cardiovascular Sciences, University of Manchester, Manchester, United Kingdom
| | - Erwin Neher
- Membrane Biophysics Laboratory, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
- Cluster of Excellence "Multiscale Bioimaging: from Molecular Machines to Networks of Excitable Cells" (MBExC), University of Göttingen, Göttingen, Germany
| | - Holger Taschenberger
- Department of Molecular Neurobiology, Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
| | - Godfrey Smith
- School of Cardiovascular and Metabolic Health, College of Medical, Veterinary, and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| |
Collapse
|
2
|
Mazhar F, Bartolucci C, Regazzoni F, Paci M, Dedè L, Quarteroni A, Corsi C, Severi S. A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics. J Physiol 2023. [PMID: 37641426 DOI: 10.1113/jp283974] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 08/02/2023] [Indexed: 08/31/2023] Open
Abstract
Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+ -handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.
Collapse
Affiliation(s)
- Fazeelat Mazhar
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena, Italy
| | - Chiara Bartolucci
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena, Italy
| | | | - Michelangelo Paci
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena, Italy
- BioMediTech, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
| | - Luca Dedè
- MOX - Dipartimento di Matematica, Politecnico di Milano, Milan, Italy
| | - Alfio Quarteroni
- MOX - Dipartimento di Matematica, Politecnico di Milano, Milan, Italy
- Mathematics Institute, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Cristiana Corsi
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena, Italy
| | - Stefano Severi
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena, Italy
| |
Collapse
|
3
|
Hoang-Trong TM, Ullah A, Lederer WJ, Jafri MS. A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation. MEMBRANES 2021; 11:989. [PMID: 34940490 PMCID: PMC8706945 DOI: 10.3390/membranes11120989] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 12/12/2021] [Accepted: 12/15/2021] [Indexed: 11/16/2022]
Abstract
Calcium (Ca2+) plays a central role in the excitation and contraction of cardiac myocytes. Experiments have indicated that calcium release is stochastic and regulated locally suggesting the possibility of spatially heterogeneous calcium levels in the cells. This spatial heterogeneity might be important in mediating different signaling pathways. During more than 50 years of computational cell biology, the computational models have been advanced to incorporate more ionic currents, going from deterministic models to stochastic models. While periodic increases in cytoplasmic Ca2+ concentration drive cardiac contraction, aberrant Ca2+ release can underly cardiac arrhythmia. However, the study of the spatial role of calcium ions has been limited due to the computational expense of using a three-dimensional stochastic computational model. In this paper, we introduce a three-dimensional stochastic computational model for rat ventricular myocytes at the whole-cell level that incorporate detailed calcium dynamics, with (1) non-uniform release site placement, (2) non-uniform membrane ionic currents and membrane buffers, (3) stochastic calcium-leak dynamics and (4) non-junctional or rogue ryanodine receptors. The model simulates spark-induced spark activation and spark-induced Ca2+ wave initiation and propagation that occur under conditions of calcium overload at the closed-cell condition, but not when Ca2+ levels are normal. This is considered important since the presence of Ca2+ waves contribute to the activation of arrhythmogenic currents.
Collapse
Affiliation(s)
- Tuan Minh Hoang-Trong
- School of Systems Biology, Krasnow Institute for Advanced Study, George Mason University, Fairfax, VA 22030, USA; (T.M.H.-T.); (A.U.)
| | - Aman Ullah
- School of Systems Biology, Krasnow Institute for Advanced Study, George Mason University, Fairfax, VA 22030, USA; (T.M.H.-T.); (A.U.)
| | - William Jonathan Lederer
- Center for Biomedical Engineering and Technology, University of Maryland School of Medicine, Baltimore, MD 21201, USA;
| | - Mohsin Saleet Jafri
- School of Systems Biology, Krasnow Institute for Advanced Study, George Mason University, Fairfax, VA 22030, USA; (T.M.H.-T.); (A.U.)
- Center for Biomedical Engineering and Technology, University of Maryland School of Medicine, Baltimore, MD 21201, USA;
| |
Collapse
|
4
|
Bazmi M, Escobar AL. Excitation-Contraction Coupling in the Goldfish ( Carassius auratus) Intact Heart. Front Physiol 2020; 11:1103. [PMID: 33041845 PMCID: PMC7518121 DOI: 10.3389/fphys.2020.01103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Accepted: 08/10/2020] [Indexed: 12/17/2022] Open
Abstract
Cardiac physiology of fish models is an emerging field given the ease of genome editing and the development of transgenic models. Several studies have described the cardiac properties of zebrafish (Denio rerio). The goldfish (Carassius auratus) belongs to the same family as the zebrafish and has emerged as an alternative model with which to study cardiac function. Here, we propose to acutely study electrophysiological and systolic Ca2+ signaling in intact goldfish hearts. We assessed the Ca2+ dynamics and the electrophysiological cardiac function of goldfish, zebrafish, and mice models, using pulsed local field fluorescence microscopy, intracellular microelectrodes, and flash photolysis in perfused hearts. We observed goldfish ventricular action potentials (APs) and Ca2+ transients to be significantly longer when compared to the zebrafish. The action potential half duration at 50% (APD50) of goldfish was 370.38 ± 8.8 ms long, and in the zebrafish they were observed to be only 83.9 ± 9.4 ms. Additionally, the half duration of the Ca2+ transients was also longer for goldfish (402.1 ± 4.4 ms) compared to the zebrafish (99.1 ± 2.7 ms). Also, blocking of the L-type Ca2+ channels with nifedipine revealed this current has a major role in defining the amplitude and the duration of goldfish Ca2+ transients. Interestingly, nifedipine flash photolysis experiments in the intact heart identified whether or not the decrease in the amplitude of Ca2+ transients was due to shorter APs. Moreover, an increase in temperature and heart rate had a strong shortening effect on the AP and Ca2+ transients of goldfish hearts. Furthermore, ryanodine (Ry) and thapsigargin (Tg) significantly reduced the amplitude of the Ca2+ transients, induced a prolongation in the APs, and altogether exhibited the degree to which the Ca2+ release from the sarcoplasmic reticulum contributed to the Ca2+ transients. We conclude that the electrophysiological properties and Ca2+ signaling in intact goldfish hearts strongly resembles the endocardial layer of larger mammals.
Collapse
Affiliation(s)
- Maedeh Bazmi
- Quantitative Systems Biology Program, School of Natural Sciences, University of California, Merced, Merced, CA, United States
| | - Ariel L Escobar
- Department of Bioengineering, School of Engineering, University of California, Merced, Merced, CA, United States
| |
Collapse
|
5
|
Abstract
Changes of intracellular Ca2+ concentration regulate many aspects of cardiac myocyte function. About 99% of the cytoplasmic calcium in cardiac myocytes is bound to buffers, and their properties will therefore have a major influence on Ca2+ signaling. This article considers the fundamental properties and identities of the buffers and how to measure them. It reviews the effects of buffering on the systolic Ca2+ transient and how this may change physiologically, and in heart failure and both atrial and ventricular arrhythmias, as well. It is concluded that the consequences of this strong buffering may be more significant than currently appreciated, and a fuller understanding is needed for proper understanding of cardiac calcium cycling and contractility.
Collapse
Affiliation(s)
- Godfrey L Smith
- Institute of Cardiovascular and Medical Sciences, College of Medical, Veterinary, and Life Sciences, University of Glasgow, UK (G.L.S.)
| | - David A Eisner
- Unit of Cardiac Physiology, Division of Cardiovascular Sciences, University of Manchester, UK (D.A.E.)
| |
Collapse
|
6
|
Denizot A, Arizono M, Nägerl UV, Soula H, Berry H. Simulation of calcium signaling in fine astrocytic processes: Effect of spatial properties on spontaneous activity. PLoS Comput Biol 2019; 15:e1006795. [PMID: 31425510 PMCID: PMC6726244 DOI: 10.1371/journal.pcbi.1006795] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Revised: 09/04/2019] [Accepted: 07/08/2019] [Indexed: 12/20/2022] Open
Abstract
Astrocytes, a glial cell type of the central nervous system, have emerged as detectors and regulators of neuronal information processing. Astrocyte excitability resides in transient variations of free cytosolic calcium concentration over a range of temporal and spatial scales, from sub-microdomains to waves propagating throughout the cell. Despite extensive experimental approaches, it is not clear how these signals are transmitted to and integrated within an astrocyte. The localization of the main molecular actors and the geometry of the system, including the spatial organization of calcium channels IP3R, are deemed essential. However, as most calcium signals occur in astrocytic ramifications that are too fine to be resolved by conventional light microscopy, most of those spatial data are unknown and computational modeling remains the only methodology to study this issue. Here, we propose an IP3R-mediated calcium signaling model for dynamics in such small sub-cellular volumes. To account for the expected stochasticity and low copy numbers, our model is both spatially explicit and particle-based. Extensive simulations show that spontaneous calcium signals arise in the model via the interplay between excitability and stochasticity. The model reproduces the main forms of calcium signals and indicates that their frequency crucially depends on the spatial organization of the IP3R channels. Importantly, we show that two processes expressing exactly the same calcium channels can display different types of calcium signals depending on the spatial organization of the channels. Our model with realistic process volume and calcium concentrations successfully reproduces spontaneous calcium signals that we measured in calcium micro-domains with confocal microscopy and predicts that local variations of calcium indicators might contribute to the diversity of calcium signals observed in astrocytes. To our knowledge, this model is the first model suited to investigate calcium dynamics in fine astrocytic processes and to propose plausible mechanisms responsible for their variability.
Collapse
Affiliation(s)
- Audrey Denizot
- INRIA, F-69603, Villeurbanne, France
- Univ Lyon, LIRIS, UMR5205 CNRS, F-69621, Villeurbanne, France
| | - Misa Arizono
- Interdisciplinary Institute for Neuroscience, Université de Bordeaux, Bordeaux, France
- Interdisciplinary Institute for Neuroscience, CNRS UMR 5297, Bordeaux, France
| | - U. Valentin Nägerl
- Interdisciplinary Institute for Neuroscience, Université de Bordeaux, Bordeaux, France
- Interdisciplinary Institute for Neuroscience, CNRS UMR 5297, Bordeaux, France
| | - Hédi Soula
- INRIA, F-69603, Villeurbanne, France
- Univ P&M Curie, CRC, INSERM UMRS 1138, F-75006, Paris, France
| | - Hugues Berry
- INRIA, F-69603, Villeurbanne, France
- Univ Lyon, LIRIS, UMR5205 CNRS, F-69621, Villeurbanne, France
| |
Collapse
|
7
|
Singh N, Adlakha N. A mathematical model for interdependent calcium and inositol 1,4,5-trisphosphate in cardiac myocyte. ACTA ACUST UNITED AC 2019. [DOI: 10.1007/s13721-019-0198-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
8
|
Finite element model to study two dimensional unsteady state calcium distribution in cardiac myocytes. ALEXANDRIA JOURNAL OF MEDICINE 2019. [DOI: 10.1016/j.ajme.2015.09.007] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
9
|
Singh N, Adlakha N. Three dimensional coupled reaction–diffusion modeling of calcium and inositol 1,4,5-trisphosphate dynamics in cardiomyocytes. RSC Adv 2019; 9:42459-42469. [PMID: 35542883 PMCID: PMC9076935 DOI: 10.1039/c9ra06929a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2019] [Accepted: 11/28/2019] [Indexed: 11/30/2022] Open
Abstract
Nanoparticles have shown great promise in improving cancer treatment efficacy by changing the intracellular calcium level through activation of intracellular mechanisms. One of the mechanisms of the killing of the cancerous cell by a nanoparticle is through elevation of the intracellular calcium level. Evidence accumulated over the past decade indicates a pivotal role for the IP3 receptor mediated Ca2+ release in the regulation of the cytosolic and the nuclear Ca2+ signals. There have been various studies done suggesting the role of IP3 receptors (IP3R) and IP3 production and degradation in cardiomyocytes. In the present work, we have proposed a three-dimensional unsteady-state mathematical model to describe the mechanism of cardiomyocytes which focuses on evaluation of various parameters that affect these coupled dynamics and elevate the cytosolic calcium concentration which can be helpful to search for novel therapies to cure these malignancies by targeting the complex calcium signaling process in cardiomyocytes. Our study suggests that there are other factors involved in this signaling which can increase the calcium level, which can help in finding treatment for cancer. The cytosolic calcium level may be controlled by IP3 signaling, leak, source influx of calcium (σ) and maximum production of IP3 (VP). We believe that the proposed model suggests new insight into finding treatment for cancer in cardiomyocytes through elevation of the cytosolic Ca2+ concentration by various parameters like leak, σ, VP and especially by other complex cell signaling dynamics, namely IP3 dynamics. We propose a three-dimensional unsteady-state mathematical model to describe the mechanism of cardiomyocytes.![]()
Collapse
Affiliation(s)
- Nisha Singh
- Applied Mathematics and Humanities Department
- SVNIT
- Surat
- India
| | - Neeru Adlakha
- Applied Mathematics and Humanities Department
- SVNIT
- Surat
- India
| |
Collapse
|
10
|
Yang F, Zhang L, Lu W, Zhang Y, Zuo W, Wang K, Zhang H. A composite visualization method for electrophysiology-morphous merging of human heart. Biomed Eng Online 2017; 16:70. [PMID: 28595607 PMCID: PMC5465514 DOI: 10.1186/s12938-017-0368-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2016] [Accepted: 06/01/2017] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Electrophysiological behavior is of great importance for analyzing the cardiac functional mechanism under cardiac physiological and pathological condition. Due to the complexity of cardiac structure and biophysiological function, visualization of a cardiac electrophysiological model compositively is still a challenge. The lack of either modality of the whole organ structure or cardiac electrophysiological behaviors makes analysis of the intricate mechanisms of cardiac dynamic function a difficult task. This study aims at exploring 3D conduction of stimulus and electrical excitation reactivity on the level of organ with the authentic fine cardiac anatomy structure. METHODS In this paper, a cardiac electrical excitation propagation model is established based on the human cardiac cross-sectional data to explore detailed cardiac electrical activities. A novel biophysical merging visualization method is then presented for biophysical integration of cardiac anatomy and electrophysiological properties in the form of the merging optical model, which provides the corresponding position, spatial relationship and the whole process in 3D space with the context of anatomical structure for representing the biophysical detailed electrophysiological activity. RESULTS The visualization result present the action potential propagation of the left ventricle within the excitation cycle with the authentic fine cardiac organ anatomy. In the visualized images, all vital organs are identified and distinguished without ambiguity. The three dimensional spatial position, relation and the process of cardiac excitation conduction and re-entry propagation in the anatomical structure during the phase of depolarization and repolarization is also shown in the result images, which exhibits the performance of a more detailed biophysical understanding of the electrophysiological kinetics of human heart in vivo. CONCLUSIONS Results suggest that the proposed merging optical model can merge cardiac electrophysiological activity with the anatomy structure. By specifying the respective opacity for the cardiac anatomy structure and the electrophysiological model in the merging attenuation function, the visualized images can provide an in-depth insight into the biophysical detailed cardiac functioning phenomena and the corresponding electrophysiological behavior mechanism, which is helpful for further speculating cardiac physiological and pathological responses and is fundamental to the cardiac research and clinical diagnoses.
Collapse
Affiliation(s)
- Fei Yang
- School of Mechanical, Electrical & Information Engineering, Shandong University, Weihai, China
| | - Lei Zhang
- School of Art and Design, Harbin University, Harbin, China
| | - Weigang Lu
- Department of Educational Technology, Ocean University of China, Qingdao, China.
| | - Yue Zhang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Wangmeng Zuo
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Kuanquan Wang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Henggui Zhang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China.,School of Physics and Astronomy, University of Manchester, Manchester, M139PL, UK
| |
Collapse
|
11
|
Greiser M, Kerfant BG, Williams GS, Voigt N, Harks E, Dibb KM, Giese A, Meszaros J, Verheule S, Ravens U, Allessie MA, Gammie JS, van der Velden J, Lederer WJ, Dobrev D, Schotten U. Tachycardia-induced silencing of subcellular Ca2+ signaling in atrial myocytes. J Clin Invest 2014; 124:4759-72. [PMID: 25329692 PMCID: PMC4347234 DOI: 10.1172/jci70102] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2013] [Accepted: 08/28/2014] [Indexed: 01/06/2023] Open
Abstract
Atrial fibrillation (AF) is characterized by sustained high atrial activation rates and arrhythmogenic cellular Ca2+ signaling instability; however, it is not clear how a high atrial rate and Ca2+ instability may be related. Here, we characterized subcellular Ca2+ signaling after 5 days of high atrial rates in a rabbit model. While some changes were similar to those in persistent AF, we identified a distinct pattern of stabilized subcellular Ca2+ signaling. Ca2+ sparks, arrhythmogenic Ca2+ waves, sarcoplasmic reticulum (SR) Ca2+ leak, and SR Ca2+ content were largely unaltered. Based on computational analysis, these findings were consistent with a higher Ca2+ leak due to PKA-dependent phosphorylation of SR Ca2+ channels (RyR2s), fewer RyR2s, and smaller RyR2 clusters in the SR. We determined that less Ca2+ release per [Ca2+]i transient, increased Ca2+ buffering strength, shortened action potentials, and reduced L-type Ca2+ current contribute to a stunning reduction of intracellular Na+ concentration following rapid atrial pacing. In both patients with AF and in our rabbit model, this silencing led to failed propagation of the [Ca2+]i signal to the myocyte center. We conclude that sustained high atrial rates alone silence Ca2+ signaling and do not produce Ca2+ signaling instability, consistent with an adaptive molecular and cellular response to atrial tachycardia.
Collapse
Affiliation(s)
- Maura Greiser
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Benoît-Gilles Kerfant
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - George S.B. Williams
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Niels Voigt
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Erik Harks
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Katharine M. Dibb
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Anne Giese
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Janos Meszaros
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Sander Verheule
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Ursula Ravens
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Maurits A. Allessie
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - James S. Gammie
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Jolanda van der Velden
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - W. Jonathan Lederer
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Dobromir Dobrev
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| | - Ulrich Schotten
- Department of Physiology, Maastricht University, Maastricht, the Netherlands. Center for Biomedical Engineering and Technology, Laboratory of Molecular Cardiology, and Department of Physiology, University of Maryland School of Medicine, Baltimore, Maryland, USA. Institute of Pharmacology, Faculty of Medicine, University Duisburg-Essen, Essen, Germany. Unit of Cardiac Physiology, Manchester Academic Health Sciences Centre, Manchester, United Kingdom. Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany. Division of Cardiac Surgery, University of Maryland Medical Center, Baltimore, Maryland, USA. Laboratory for Physiology, VU University Medical Center, Amsterdam, the Netherlands
| |
Collapse
|
12
|
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
|
13
|
Sensitivity of rabbit ventricular action potential and Ca²⁺ dynamics to small variations in membrane currents and ion diffusion coefficients. BIOMED RESEARCH INTERNATIONAL 2013; 2013:565431. [PMID: 24222910 PMCID: PMC3814049 DOI: 10.1155/2013/565431] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/28/2013] [Accepted: 08/19/2013] [Indexed: 12/19/2022]
Abstract
Little is known about how small variations in ionic currents and Ca²⁺ and Na⁺ diffusion coefficients impact action potential and Ca²⁺ dynamics in rabbit ventricular myocytes. We applied sensitivity analysis to quantify the sensitivity of Shannon et al. model (Biophys. J., 2004) to 5%-10% changes in currents conductance, channels distribution, and ion diffusion in rabbit ventricular cells. We found that action potential duration and Ca²⁺ peaks are highly sensitive to 10% increase in L-type Ca²⁺ current; moderately influenced by 10% increase in Na⁺-Ca²⁺ exchanger, Na⁺-K⁺ pump, rapid delayed and slow transient outward K⁺ currents, and Cl⁻ background current; insensitive to 10% increases in all other ionic currents and sarcoplasmic reticulum Ca²⁺ fluxes. Cell electrical activity is strongly affected by 5% shift of L-type Ca²⁺ channels and Na⁺-Ca²⁺ exchanger in between junctional and submembrane spaces while Ca²⁺-activated Cl⁻-channel redistribution has the modest effect. Small changes in submembrane and cytosolic diffusion coefficients for Ca²⁺, but not in Na⁺ transfer, may alter notably myocyte contraction. Our studies highlight the need for more precise measurements and further extending and testing of the Shannon et al. model. Our results demonstrate usefulness of sensitivity analysis to identify specific knowledge gaps and controversies related to ventricular cell electrophysiology and Ca²⁺ signaling.
Collapse
|
14
|
Yu Z, Wang J, Gao Z, Xu M, Hoshijima M. New software developments for quality mesh generation and optimization from biomedical imaging data. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2013; 113:226-40. [PMID: 24252469 PMCID: PMC3836056 DOI: 10.1016/j.cmpb.2013.08.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2013] [Revised: 08/16/2013] [Accepted: 08/16/2013] [Indexed: 06/02/2023]
Abstract
In this paper we present a new software toolkit for generating and optimizing surface and volumetric meshes from three-dimensional (3D) biomedical imaging data, targeted at image-based finite element analysis of some biomedical activities in a single material domain. Our toolkit includes a series of geometric processing algorithms including surface re-meshing and quality-guaranteed tetrahedral mesh generation and optimization. All methods described have been encapsulated into a user-friendly graphical interface for easy manipulation and informative visualization of biomedical images and mesh models. Numerous examples are presented to demonstrate the effectiveness and efficiency of the described methods and toolkit.
Collapse
Affiliation(s)
- Zeyun Yu
- Department of Computer Science, University of Wisconsin at Milwaukee, USA.
| | | | | | | | | |
Collapse
|
15
|
Haq KT, Daniels RE, Miller LS, Miura M, ter Keurs HEDJ, Bungay SD, Stuyvers BD. Evoked centripetal Ca(2+) mobilization in cardiac Purkinje cells: insight from a model of three Ca(2+) release regions. J Physiol 2013; 591:4301-19. [PMID: 23897231 DOI: 10.1113/jphysiol.2013.253583] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Abstract
Despite strong suspicion that abnormal Ca(2+) handling in Purkinje cells (P-cells) is implicated in life-threatening forms of ventricular tachycardias, the mechanism underlying the Ca(2+) cycling of these cells under normal conditions is still unclear. There is mounting evidence that P-cells have a unique Ca(2+) handling system. Notably complex spontaneous Ca(2+) activity was previously recorded in canine P-cells and was explained by a mechanistic hypothesis involving a triple layered system of Ca(2+) release channels. Here we examined the validity of this hypothesis for the electrically evoked Ca(2+) transient which was shown, in the dog and rabbit, to occur progressively from the periphery to the interior of the cell. To do so, the hypothesis was incorporated in a model of intracellular Ca(2+) dynamics which was then used to reproduce numerically the Ca(2+) activity of P-cells under stimulated conditions. The modelling was thus performed through a 2D computational array that encompassed three distinct Ca(2+) release nodes arranged, respectively, into three consecutive adjacent regions. A system of partial differential equations (PDEs) expressed numerically the principal cellular functions that modulate the local cytosolic Ca(2+) concentration (Cai). The apparent node-to-node progression of elevated Cai was obtained by combining Ca(2+) diffusion and 'Ca(2+)-induced Ca(2+) release'. To provide the modelling with a reliable experimental reference, we first re-examined the Ca(2+) mobilization in swine stimulated P-cells by 2D confocal microscopy. As reported earlier for the dog and rabbit, a centripetal Ca(2+) transient was readily visible in 22 stimulated P-cells from six adult Yucatan swine hearts (pacing rate: 0.1 Hz; pulse duration: 25 ms, pulse amplitude: 10% above threshold; 1 mm Ca(2+); 35°C; pH 7.3). An accurate replication of the observed centripetal Ca(2+) propagation was generated by the model for four representative cell examples and confirmed by statistical comparisons of simulations against cell data. Selective inactivation of Ca(2+) release regions of the computational array showed that an intermediate layer of Ca(2+) release nodes with an ~30-40% lower Ca(2+) activation threshold was required to reproduce the phenomenon. Our computational analysis was therefore fully consistent with the activation of a triple layered system of Ca(2+) release channels as a mechanism of centripetal Ca(2+) signalling in P-cells. Moreover, the model clearly indicated that the intermediate Ca(2+) release layer with increased sensitivity for Ca(2+) plays an important role in the specific intracellular Ca(2+) mobilization of Purkinje fibres and could therefore be a relevant determinant of cardiac conduction.
Collapse
Affiliation(s)
- Kazi T Haq
- B. D. Stuyvers: Memorial University, Faculty of Medicine, Division of BioMedical Sciences, 300 Prince Phillip Bd, St John's, NL, A1B 3V6, Canada.
| | | | | | | | | | | | | |
Collapse
|
16
|
Recruitment of the intracellular Ca2+ by ultrashort electric stimuli: the impact of pulse duration. Cell Calcium 2013; 54:145-50. [PMID: 23777980 DOI: 10.1016/j.ceca.2013.05.008] [Citation(s) in RCA: 83] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2013] [Revised: 05/20/2013] [Accepted: 05/24/2013] [Indexed: 01/25/2023]
Abstract
Nanosecond-duration electric stimuli are distinguished by the ability to permeabilize intracellular membranes and recruit Ca2+ from intracellular stores. We quantified this effect in non-excitable cells (CHO) using ratiometric Ca2+ imaging with Fura-2. In a Ca(2+)-free medium, 10-, 60-, and 300-ns stimuli evoked Ca2+ transients by mobilization of Ca2+ from the endoplasmic reticulum. With 2 mM external Ca2+, the transients included both extra- and intracellular components. The recruitment of intracellular Ca2+ increased as the stimulus duration decreased. At the threshold of 200-300 nM, the transients were amplified by calcium-induced calcium release. We conclude that nanosecond stimuli mimic Ca2+ signaling while bypassing the usual receptor- and channels-mediated cascades. The recruitment of the intracellular Ca2+ can be controlled by the duration of the stimulus.
Collapse
|
17
|
Do calcium buffers always slow down the propagation of calcium waves? J Math Biol 2012; 67:1587-632. [PMID: 23076831 DOI: 10.1007/s00285-012-0605-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2012] [Revised: 09/27/2012] [Indexed: 10/27/2022]
Abstract
Calcium buffers are large proteins that act as binding sites for free cytosolic calcium. Since a large fraction of cytosolic calcium is bound to calcium buffers, calcium waves are widely observed under the condition that free cytosolic calcium is heavily buffered. In addition, all physiological buffered excitable systems contain multiple buffers with different affinities. It is thus important to understand the properties of waves in excitable systems with the inclusion of buffers. There is an ongoing controversy about whether or not the addition of calcium buffers into the system always slows down the propagation of calcium waves. To solve this controversy, we incorporate the buffering effect into the generic excitable system, the FitzHugh-Nagumo model, to get the buffered FitzHugh-Nagumo model, and then to study the effect of the added buffer with large diffusivity on traveling waves of such a model in one spatial dimension. We can find a critical dissociation constant (K = K(a)) characterized by system excitability parameter a such that calcium buffers can be classified into two types: weak buffers (K ∈ (K(a), ∞)) and strong buffers (K ∈ (0, K(a))). We analytically show that the addition of weak buffers or strong buffers but with its total concentration b(0)(1) below some critical total concentration b(0,c)(1) into the system can generate a traveling wave of the resulting system which propagates faster than that of the origin system, provided that the diffusivity D1 of the added buffers is sufficiently large. Further, the magnitude of the wave speed of traveling waves of the resulting system is proportional to √D1 as D1 --> ∞. In contrast, the addition of strong buffers with the total concentration b(0)(1) > b(0,c)(1) into the system may not be able to support the formation of a biologically acceptable wave provided that the diffusivity D1 of the added buffers is sufficiently large.
Collapse
|
18
|
Kekenes-Huskey PM, Cheng Y, Hake JE, Sachse FB, Bridge JH, Holst MJ, McCammon JA, McCulloch AD, Michailova AP. Modeling effects of L-type ca(2+) current and na(+)-ca(2+) exchanger on ca(2+) trigger flux in rabbit myocytes with realistic T-tubule geometries. Front Physiol 2012; 3:351. [PMID: 23060801 PMCID: PMC3463892 DOI: 10.3389/fphys.2012.00351] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2012] [Accepted: 08/16/2012] [Indexed: 12/26/2022] Open
Abstract
The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca2+ channel (LCC) clustering, and allosteric activation of Na+/Ca2+ exchanger by L-type Ca2+ current affects intracellular Ca2+ dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na+/Ca2+ exchanger, sarcolemmal Ca2+ pump, and sarcolemmal Ca2+ leak), and stationary and mobile Ca2+ buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca2+. We obtained parameters from voltage-clamp protocols of L-type Ca2+ current and line-scan recordings of Ca2+ concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca2+ transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca2+ concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca2+ crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca2+ flux distribution. The model additionally predicts that local Ca2+ trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca2+ trigger flux. We found also that the activation of allosteric Ca2+-binding sites on the Na+/Ca2+ exchanger could provide a mechanism for regulating global and local Ca2+ trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na+/Ca2+ exchanger fluxes to intracellular Ca2+ dynamics.
Collapse
|
19
|
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
|
20
|
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
|
21
|
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
|
22
|
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
|
23
|
Koivumäki JT, Korhonen T, Tavi P. Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study. PLoS Comput Biol 2011; 7:e1001067. [PMID: 21298076 PMCID: PMC3029229 DOI: 10.1371/journal.pcbi.1001067] [Citation(s) in RCA: 87] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2010] [Accepted: 12/21/2010] [Indexed: 02/02/2023] Open
Abstract
Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca²+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca²+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca²+ dynamics: 1) the biphasic increment during the upstroke of the Ca²+ transient resulting from the delay between the peripheral and central SR Ca²+ release, and 2) the relative contribution of SL Ca²+ current and SR Ca²+ release to the Ca²+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca²+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca²+ release sites define the interface between Ca²+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca²+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca²+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca²+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.
Collapse
Affiliation(s)
- Jussi T. Koivumäki
- Department of Biotechnology and Molecular Medicine, A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Topi Korhonen
- Department of Biotechnology and Molecular Medicine, A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Pasi Tavi
- Department of Biotechnology and Molecular Medicine, A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
24
|
Cheng Y, Yu Z, Hoshijima M, Holst MJ, McCulloch AD, McCammon JA, Michailova AP. Numerical analysis of Ca2+ signaling in rat ventricular myocytes with realistic transverse-axial tubular geometry and inhibited sarcoplasmic reticulum. PLoS Comput Biol 2010; 6:e1000972. [PMID: 21060856 PMCID: PMC2965743 DOI: 10.1371/journal.pcbi.1000972] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2009] [Accepted: 09/23/2010] [Indexed: 12/21/2022] Open
Abstract
The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca2+). To investigate how the t-tubule microanatomy and the distribution of membrane Ca2+ flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca2+ signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca2+ signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca2+ flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca2+ buffers, including the Ca2+ indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca2+ transients was found when the Ca2+ flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca2+ signals are predicted. Even at modest elevation of Ca2+, reached during Ca2+ influx, large and steep Ca2+ gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca2+ flux density along the cell membrane support initiation and propagation of Ca2+ waves in rat myocytes. In cardiac muscle cells, calcium (Ca2+) is best known for its role in contraction activation. A remarkable amount of quantitative data on cardiac cell structure, ion-transporting protein distributions and intracellular Ca2+ dynamics has been accumulated. Various alterations in the protein distributions or cell ultra-structure are now recognized to be the primary mechanisms of cardiac dysfunction in a diverse range of common pathologies including cardiac arrhythmias and hypertrophy. Using a 3-D computational model, incorporating more realistic transverse-axial t-tubule geometry and considering geometric irregularities and inhomogeneities in the distribution of ion-transporting proteins, we analyze several important spatial and temporal features of Ca2+ signaling in rat ventricular myocytes. This study demonstrates that the computational models could serve as powerful tools for prediction and analyses of how the Ca2+ dynamics and cardiac excitation-contraction coupling are regulated under normal conditions or certain pathologies. The use of computational and mathematical approaches will help also to better understand aspects of cell functions that are not currently amenable to experimental investigation.
Collapse
Affiliation(s)
- Yuhui Cheng
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Zeyun Yu
- Department of Computer Science, University of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Masahiko Hoshijima
- Department of Medicine, University of California San Diego, La Jolla, California, United States of America
| | - Michael J. Holst
- Department of Mathematics, University of California San Diego, La Jolla, California, United States of America
| | - Andrew D. McCulloch
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - J. Andrew McCammon
- Department of Chemistry and Biochemistry, Department of Pharmacology, University of California San Diego, La Jolla, California, United States of America
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, California, United States of America
| | - Anushka P. Michailova
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
25
|
Niggli E. Ca2+ on the move: what is better - fire and forget, or fire-diffuse-fire? J Physiol 2010; 588:1391-2. [PMID: 20436043 DOI: 10.1113/jphysiol.2010.190306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Affiliation(s)
- Ernst Niggli
- Department of Physiology, University of Bern, Buehlplatz 5, 3012 Bern, Switzerland.
| |
Collapse
|
26
|
Fink M, Niederer SA, Cherry EM, Fenton FH, Koivumäki JT, Seemann G, Thul R, Zhang H, Sachse FB, Beard D, Crampin EJ, Smith NP. Cardiac cell modelling: observations from the heart of the cardiac physiome project. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2010; 104:2-21. [PMID: 20303361 DOI: 10.1016/j.pbiomolbio.2010.03.002] [Citation(s) in RCA: 108] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2009] [Revised: 12/06/2009] [Accepted: 03/04/2010] [Indexed: 10/19/2022]
Abstract
In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field.
Collapse
Affiliation(s)
- Martin Fink
- Department of Physiology, Anatomy and Genetics, University of Oxford, OX1 3JP, United Kingdom
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
27
|
Lu S, Michailova A, Saucerman J, Cheng Y, Yu Z, Kaiser T, Li W, Bank R, Holst M, McCammon J, Hayashi T, Hoshijima M, Arzberger P, McCulloch A. Multiscale modeling in rodent ventricular myocytes. ACTA ACUST UNITED AC 2009; 28:46-57. [PMID: 19349251 DOI: 10.1109/memb.2009.931787] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
There is a growing body of experimental evidence suggesting that the Ca(2+) signaling in ventricular myocytes is characterized by a high gradient near the cell membrane and a more uniform Ca(2+) distribution in the cell interior [1]--[7]. An important reason for this phenomenon might be that in these cells the t-tubular system forms a network of extracellular space, extending deep into the cell interior. This allows the electrical signal, that propagates rapidly along the cell membrane, to reach the vicinity of the sarcoplasmic reticulum (SR), where intracellular Ca(2+) required for myofilament activation is stored [1], [8]--[11]. Early studies of cardiac muscle showed that the t-tubules are found at intervals of about 2 lm along the longitudinal cell axis in close proximity to the Z-disks of the sarcomeres [12]. Subsequent studies have demonstrated that the t-tubular system has also longitudinal extensions [9]--[11], [13].
Collapse
Affiliation(s)
- Shaoying Lu
- Department of Bioengineering, University of California San Diego, La Jolla, 92093-0412, USA.
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Model of excitation-contraction coupling of rat neonatal ventricular myocytes. Biophys J 2009; 96:1189-209. [PMID: 19186154 DOI: 10.1016/j.bpj.2008.10.026] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2008] [Accepted: 10/27/2008] [Indexed: 11/22/2022] Open
Abstract
The neonatal rat ventricular myocyte culture is one of the most popular experimental cardiac cell models. To our knowledge, the excitation-contraction coupling (ECC) of these cells, i.e., the process linking the electrical activity to the cytosolic Ca2+ transient and contraction, has not been previously analyzed, nor has it been presented as a complete system in detail. Neonatal cardiomyocytes are in the postnatal developmental stage, and therefore, the features of their ECC differ vastly from those of adult ventricular myocytes. We present the first complete analysis of ECC in these cells by characterizing experimentally the action potential and calcium signaling and developing the first mathematical model of ECC in neonatal cardiomyocytes that we know of. We show that in comparison to adult cardiomyocytes, neonatal cardiomyocytes have long action potentials, heterogeneous cytosolic Ca2+ signals, weaker sarcoplasmic reticulum Ca2+ handling, and stronger sarcolemmal Ca2+ handling, with a significant contribution by the Na+/Ca2+ exchanger. The developed model reproduces faithfully the ECC of rat neonatal cardiomyocytes with a novel description of spatial cytosolic [Ca2+] signals. Simulations also demonstrate how an increase in the cell size (hypertrophy) affects the ECC in neonatal cardiomyocytes. This model of ECC in developing cardiomyocytes provides a platform for developing future models of cardiomyocytes at different developmental stages.
Collapse
|
29
|
Yu Z, Holst MJ, Hayashi T, Bajaj CL, Ellisman MH, McCammon JA, Hoshijima M. Three-dimensional geometric modeling of membrane-bound organelles in ventricular myocytes: bridging the gap between microscopic imaging and mathematical simulation. J Struct Biol 2008; 164:304-13. [PMID: 18835449 DOI: 10.1016/j.jsb.2008.09.004] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2008] [Revised: 09/09/2008] [Accepted: 09/10/2008] [Indexed: 10/21/2022]
Abstract
A general framework of image-based geometric processing is presented to bridge the gap between three-dimensional (3D) imaging that provides structural details of a biological system and mathematical simulation where high-quality surface or volumetric meshes are required. A 3D density map is processed in the order of image pre-processing (contrast enhancement and anisotropic filtering), feature extraction (boundary segmentation and skeletonization), and high-quality and realistic surface (triangular) and volumetric (tetrahedral) mesh generation. While the tool-chain described is applicable to general types of 3D imaging data, the performance is demonstrated specifically on membrane-bound organelles in ventricular myocytes that are imaged and reconstructed with electron microscopic (EM) tomography and two-photon microscopy (T-PM). Of particular interest in this study are two types of membrane-bound Ca(2+)-handling organelles, namely, transverse tubules (T-tubules) and junctional sarcoplasmic reticulum (jSR), both of which play an important role in regulating the excitation-contraction (E-C) coupling through dynamic Ca(2+) mobilization in cardiomyocytes.
Collapse
Affiliation(s)
- Zeyun Yu
- Department of Mathematics, University of California, San Diego, La Jolla, CA 92093, USA.
| | | | | | | | | | | | | |
Collapse
|
30
|
Zima AV, Picht E, Bers DM, Blatter LA. Termination of cardiac Ca2+ sparks: role of intra-SR [Ca2+], release flux, and intra-SR Ca2+ diffusion. Circ Res 2008; 103:e105-15. [PMID: 18787194 DOI: 10.1161/circresaha.107.183236] [Citation(s) in RCA: 127] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Ca(2+) release from cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) is regulated by dyadic cleft [Ca(2+)] and intra-SR free [Ca(2+)] ([Ca(2+)](SR)). Robust SR Ca(2+) release termination is important for stable excitation-contraction coupling, and partial [Ca(2+)](SR) depletion may contribute to release termination. Here, we investigated the regulation of SR Ca(2+) release termination of spontaneous local SR Ca(2+) release events (Ca(2+) sparks) by [Ca(2+)](SR), release flux, and intra-SR Ca(2+) diffusion. We simultaneously measured Ca(2+) sparks and Ca(2+) blinks (localized elementary [Ca(2+)](SR) depletions) in permeabilized ventricular cardiomyocytes over a wide range of SR Ca(2+) loads and release fluxes. Sparks terminated via a [Ca(2+)](SR)-dependent mechanism at a fixed [Ca(2+)](SR) depletion threshold independent of the initial [Ca(2+)](SR) and release flux. Ca(2+) blink recovery depended mainly on intra-SR Ca(2+) diffusion rather than SR Ca(2+) uptake. Therefore, the large variation in Ca(2+) blink recovery rates at different release sites occurred because of differences in the degree of release site interconnection within the SR network. When SR release flux was greatly reduced, long-lasting release events occurred from well-connected junctions. These junctions could sustain release because local SR Ca(2+) release and [Ca(2+)](SR) refilling reached a balance, preventing [Ca(2+)](SR) from depleting to the termination threshold. Prolonged release events eventually terminated at a steady [Ca(2+)](SR), indicative of a slower, [Ca(2+)](SR)-independent termination mechanism. These results demonstrate that there is high variability in local SR connectivity but that SR Ca(2+) release terminates at a fixed [Ca(2+)](SR) termination threshold. Thus, reliable SR Ca(2+) release termination depends on tight RyR regulation by [Ca(2+)](SR).
Collapse
Affiliation(s)
- Aleksey V Zima
- Department of Molecular Biophysics and Physiology, Rush University Medical Center, Chicago, IL 60612, USA
| | | | | | | |
Collapse
|
31
|
Ter Keurs HEDJ, Shinozaki T, Zhang YM, Wakayama Y, Sugai Y, Kagaya Y, Miura M, Boyden PA, Stuyvers BDM, Landesberg A. Sarcomere mechanics in uniform and nonuniform cardiac muscle: a link between pump function and arrhythmias. Ann N Y Acad Sci 2008; 1123:79-95. [PMID: 18375580 DOI: 10.1196/annals.1420.010] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Starling's law and the end-systolic pressure-volume relationship (ESPVR) reflect the effect of sarcomere length (SL) on the development of stress (sigma) and shortening by myocytes in the uniform ventricle. We show here that tetanic contractions of rat cardiac trabeculae exhibit a sigma-SL relationship at saturating [Ca2+] that depends on sarcomere geometry in a manner similar to that of skeletal sarcomeres and the existence of opposing forces in cardiac muscle shortened below slack length. The sigma-SL -[Ca2+](free) relationships (sigma-SL-Ca relationships) at submaximal [Ca2+] in intact and skinned trabeculae were similar, although the sensitivity for Ca2+ of intact muscle was higher. We analyzed the mechanisms underlying the sigma-SL-Ca relationship by using a kinetic model assuming that the rates of Tn-C Ca2+ binding and/or cross-bridge (XB) cycling are determined by either the SL, [Ca2+], or sigma. We analyzed the correlation between the model results and steady-state sigma measurements at varied SL at [Ca2+] from skinned rat cardiac trabeculae to test the hypotheses that the dominant feedback mechanism is SL-, sigma-, or [Ca2+]-dependent, and that the feedback mechanism regulates Tn-C Ca2+ affinity, XB kinetics, or the unitary XB force. The analysis strongly suggests that the feedback of the number of strong XBs to cardiac Tn-C Ca2+ affinity is the dominant mechanism regulating XB recruitment. Using this concept in a model of twitch-sigma accurately reproduced the sigma-SL-Ca relationship and the time courses of twitch sigma and the intracellular [Ca2+]i. The foregoing concept has equally important repercussions for the nonuniformly contracting heart, in which arrhythmogenic Ca2+ waves arise from weakened areas in the cardiac muscle. These Ca2+ waves can reversibly be induced with nonuniform excitation-contraction coupling (ECC) by the cycle of stretch and release in the border zone between the damaged and intact regions. Stimulus trains induced propagating Ca2+ waves and reversibly induced arrhythmias. We hypothesize that rapid force loss by the sarcomeres in the border zone during relaxation causes Ca2+ release from Tn-C and initiates Ca2+ waves propagated by the sarcoplasmic reticulum (SR). Modeling of the response of the cardiac twitch to rapid force changes using the feedback concept uniquely predicts the occurrence of [Ca2+]i transients as a result of accelerated Ca2+ dissociation from Tn-C. These results are consistent with the hypothesis that a force feedback to Ca2+ binding by Tn-C is responsible for Starling's law and the ESPVR in the uniform myocardium and leads to a surge of Ca2+ released by the myofilaments during relaxation in the nonuniform myocardium, which initiates arrhythmogenic propagating Ca2+ release by the SR.
Collapse
Affiliation(s)
- Henk E D J Ter Keurs
- Department of Physiology, School of Medicine, University of Calgary, 3330 Hospital Dr., N.W., Calgary, Alberta T2N 4N1, Canada.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
32
|
ter Keurs HE, Shinozaki T, Zhang YM, Zhang ML, Wakayama Y, Sugai Y, Kagaya Y, Miura M, Boyden PA, Stuyvers BD, Landesberg A. Sarcomere mechanics in uniform and non-uniform cardiac muscle: A link between pump function and arrhythmias. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2008; 97:312-31. [DOI: 10.1016/j.pbiomolbio.2008.02.013] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
33
|
Dan P, Lin E, Huang J, Biln P, Tibbits GF. Three-dimensional distribution of cardiac Na+-Ca2+ exchanger and ryanodine receptor during development. Biophys J 2007; 93:2504-18. [PMID: 17557789 PMCID: PMC1965441 DOI: 10.1529/biophysj.107.104943] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2007] [Accepted: 05/29/2007] [Indexed: 11/18/2022] Open
Abstract
Mechanisms of cardiac excitation-contraction coupling in neonates are still not clearly defined. Previous work in neonates shows reverse-mode Na(+)-Ca(2+) exchange to be the primary route of Ca(2+) entry during systole and the neonatal sarcoplasmic reticulum to have similar capability as that of adult in storing and releasing Ca(2+). We investigated Na(+)-Ca(2+) exchanger (NCX) and ryanodine receptor (RyR) distribution in developing ventricular myocytes using immunofluorescence, confocal microscopy, and digital image analysis. In neonates, both NCX and RyR clusters on the surface of the cell displayed a short longitudinal periodicity of approximately 0.7 microm. However, by adulthood, both proteins were also found in the interior. In the adult, clusters of NCX on the surface of the cell retained the approximately 0.7-microm periodicity whereas clusters of RyR adopted a longer longitudinal periodicity of approximately 2.0 microm. This suggests that neonatal myocytes also have a peri-M-line RyR distribution that is absent in adult myocytes. NCX and RyR colocalized voxel density was maximal in neonates and declined significantly with ontogeny. We conclude in newborns, Ca(2+) influx via NCX could potentially activate the dense network of peripheral Ca(2+) stores via peripheral couplings, evoking Ca(2+)-induced Ca(2+) release.
Collapse
Affiliation(s)
- Pauline Dan
- Cardiac Membrane Research Laboratory, Simon Fraser University, Burnaby, British Columbia, Canada
| | | | | | | | | |
Collapse
|
34
|
Mironova LA, Mironov SL. Approximate analytical time-dependent solutions to describe large-amplitude local calcium transients in the presence of buffers. Biophys J 2007; 94:349-58. [PMID: 17872951 PMCID: PMC2157246 DOI: 10.1529/biophysj.107.113340] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Local Ca(2+) signaling controls many neuronal functions, which is often achieved through spatial localization of Ca(2+) signals. These nanodomains are formed due to combined effects of Ca(2+) diffusion and binding to the cytoplasmic buffers. In this article we derived simple analytical expressions to describe Ca(2+) diffusion in the presence of mobile and immobile buffers. A nonlinear character of the reaction-diffusion problem was circumvented by introducing a logarithmic approximation of the concentration term. The obtained formulas reproduce free Ca(2+) levels up to 50 microM and their changes in the millisecond range. Derived equations can be useful to predict spatiotemporal profiles of large-amplitude [Ca(2+)] transients, which participate in various physiological processes.
Collapse
Affiliation(s)
- Lidia A Mironova
- DFG-Center of Molecular Physiology of the Brain, Department of Neuro- and Sensory Physiology, Georg-August-University, Göttingen, Germany
| | | |
Collapse
|
35
|
Michailova A, Lorentz W, McCulloch A. Modeling transmural heterogeneity of K(ATP) current in rabbit ventricular myocytes. Am J Physiol Cell Physiol 2007; 293:C542-57. [PMID: 17329404 DOI: 10.1152/ajpcell.00148.2006] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
To investigate the mechanisms regulating excitation-metabolic coupling in rabbit epicardial, midmyocardial, and endocardial ventricular myocytes we extended the LabHEART model (Puglisi JL and Bers DM. Am J Physiol Cell Physiol 281: C2049-C2060, 2001). We incorporated equations for Ca(2+) and Mg(2+) buffering by ATP and ADP, equations for nucleotide regulation of ATP-sensitive K(+) channel and L-type Ca(2+) channel, Na(+)-K(+)-ATPase, and sarcolemmal and sarcoplasmic Ca(2+)-ATPases, and equations describing the basic pathways (creatine and adenylate kinase reactions) known to communicate the flux changes generated by intracellular ATPases. Under normal conditions and during 20 min of ischemia, the three regions were characterized by different I(Na), I(to), I(Kr), I(Ks), and I(Kp) channel properties. The results indicate that the ATP-sensitive K(+) channel is activated by the smallest reduction in ATP in epicardial cells and largest in endocardial cells when cytosolic ADP, AMP, PCr, Cr, P(i), total Mg(2+), Na(+), K(+), Ca(2+), and pH diastolic levels are normal. The model predicts that only K(ATP) ionophore (Kir6.2 subunit) and not the regulatory subunit (SUR2A) might differ from endocardium to epicardium. The analysis suggests that during ischemia, the inhomogeneous accumulation of the metabolites in the tissue sublayers may alter in a very irregular manner the K(ATP) channel opening through metabolic interactions with the endogenous PI cascade (PIP(2), PIP) that in turn may cause differential action potential shortening among the ventricular myocyte subtypes. The model predictions are in qualitative agreement with experimental data measured under normal and ischemic conditions in rabbit ventricular myocytes.
Collapse
Affiliation(s)
- Anushka Michailova
- Dept of Bioengineering, PFBH 241, University of California San Diego, La Jolla, CA 92093-0412, USA.
| | | | | |
Collapse
|
36
|
Santamaria F, Wils S, De Schutter E, Augustine GJ. Anomalous diffusion in Purkinje cell dendrites caused by spines. Neuron 2007; 52:635-48. [PMID: 17114048 PMCID: PMC1994115 DOI: 10.1016/j.neuron.2006.10.025] [Citation(s) in RCA: 137] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2005] [Revised: 02/15/2006] [Accepted: 10/18/2006] [Indexed: 11/28/2022]
Abstract
We combined local photolysis of caged compounds with fluorescence imaging to visualize molecular diffusion within dendrites of cerebellar Purkinje cells. Diffusion of a volume marker, fluorescein dextran, within spiny dendrites was remarkably slow in comparison to its diffusion in smooth dendrites. Computer simulations indicate that this retardation is due to a transient trapping of molecules within dendritic spines, yielding anomalous diffusion. We considered the influence of spine trapping on the diffusion of calcium ions (Ca(2+)) and inositol-1,4,5-triphospate (IP(3)), two synaptic second messengers. Diffusion of IP(3) was strongly influenced by the presence of dendritic spines, while Ca(2+) was removed so rapidly that it could not diffuse far enough to be trapped. We conclude that an important function of dendritic spines may be to trap chemical signals and thereby create slowed anomalous diffusion within dendrites.
Collapse
Affiliation(s)
- Fidel Santamaria
- Department of Neurobiology, Duke University Medical Center, PO Box 3209, Durham, North Carolina 27710, USA
| | | | | | | |
Collapse
|
37
|
Thomas RC, Postma M. Dynamic and static calcium gradients inside large snail (Helix aspersa) neurones detected with calcium-sensitive microelectrodes. Cell Calcium 2006; 41:365-78. [PMID: 16962659 PMCID: PMC1885971 DOI: 10.1016/j.ceca.2006.07.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2006] [Revised: 07/20/2006] [Accepted: 07/26/2006] [Indexed: 11/15/2022]
Abstract
We have used quartz Ca2+-sensitive microelectrodes (CASMs) in large voltage-clamped snail neurones to investigate the inward spread of Ca2+ after a brief depolarisation. Both steady state and [Ca2+]i transients changed with depth of penetration. When the CASM tip was within 20 μm of the far side of the cell the [Ca2+]i transient time to peak was 4.4 ± 0.5 s, rising to 14.7 ± 0.7 s at a distance of 80 μm. We estimate that the Ca2+ transients travelled centripetally at an average speed of 6 μm2 s−1 and decreased in size by half over a distance of about 45 μm. Cyclopiazonic acid had little effect on the size and time to peak of Ca2+ transients but slowed their recovery significantly. This suggests that the endoplasmic reticulum curtails rather than reinforces the transients. Injecting the calcium buffer BAPTA made the Ca2+ transients more uniform in size and increased their times to peak and rates of recovery near the membrane. We have developed a computational model for the transients, which includes diffusion, uptake and Ca2+ extrusion. Good fits were obtained with a rather large apparent diffusion coefficient of about 90 ± 20 μm2 s−1.This may assist fast recovery by extrusion.
Collapse
Affiliation(s)
- Roger C Thomas
- Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge CB2 3EG, UK.
| | | |
Collapse
|
38
|
Soler-Llavina GJ, Sabatini BL. Synapse-specific plasticity and compartmentalized signaling in cerebellar stellate cells. Nat Neurosci 2006; 9:798-806. [PMID: 16680164 DOI: 10.1038/nn1698] [Citation(s) in RCA: 127] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2006] [Accepted: 04/17/2006] [Indexed: 11/09/2022]
Abstract
Here we demonstrate that cerebellar stellate cells diffusionally isolate synaptically evoked signals in dendrites and are capable of input-specific synaptic plasticity. Sustained activity of parallel fibers induces a form of long-term depression that requires opening of calcium (Ca(2+))-permeable AMPA-type glutamate receptors (CP-AMPARs) and signaling through class 1 metabotropic glutamate receptors (mGluR1) and CB1 receptors. This depression is induced by postsynaptic increases in Ca(2+) concentration ([Ca(2+)]) and is limited to activated synapses. To understand how synapse-specific plasticity is induced by diffusible second messengers in aspiny dendrites, we examined diffusion of Ca(2+) and small molecules within stellate cell dendrites. Activation of a single parallel fiber opened CP-AMPARs, generating long-lived Ca(2+) transients that were confined to submicron dendritic stretches. The diffusion of Ca(2+) was severely retarded due to interactions with parvalbumin and a general restriction of small molecule mobility. Thus stellate cell dendrites spatially restrict signaling cascades that lead from CP-AMPAR activation to endocannabinoid production and trigger the selective regulation of active synapses.
Collapse
Affiliation(s)
- Gilberto J Soler-Llavina
- Department of Neurobiology, Harvard Medical School, 220 Longwood Avenue, Boston, Massachusetts 02115, USA
| | | |
Collapse
|
39
|
Ter Keurs HEDJ, Wakayama Y, Miura M, Stuyvers BD, Boyden PA, Landesberg A. Spatial Nonuniformity of Contraction Causes Arrhythmogenic Ca2+Waves in Rat Cardiac Muscle. Ann N Y Acad Sci 2006; 1047:345-65. [PMID: 16093510 DOI: 10.1196/annals.1341.031] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Landesberg and Sideman's four state model of the cardiac cross-bridge (XB) hypothesizes a feedback of force development to Ca(2+) binding by troponin C (TnC). We have further modeled this behavior and observed that the force (F)-Ca(2+) relationship as well as the F-sarcomere length (SL) relationship and the time course of F and Ca(2+) transients in cardiac muscle can be reproduced faithfully by a single effect of F on deformation of the TnC-Ca complex and, thereby, on the dissociation rate of Ca(2+). Furthermore, this feedback predicts that rapid decline of F in the activated sarcomere causes release of Ca(2+) from TnC-Ca(2+), which is sufficient to initiate arrhythmogenic Ca(2+) release from the sarcoplasmic reticulum (SR). This work investigated the initiation of Ca(2+) waves underlying triggered propagated contractions (TPCs) in rat cardiac trabeculae under conditions that simulate functional nonuniformity caused by mechanical or ischemic local damage of the myocardium. A mechanical discontinuity along the trabeculae was created by exposing the preparation to a small constant flow jet of solution that reduces excitation-contraction coupling in myocytes within that segment. Force was measured, and SL as well as [Ca(2+)](i) were measured regionally. When the jet contained caffeine, 2,3-butanedione monoxime or low-[Ca(2+)], muscle-twitch F decreased and the sarcomeres in the exposed segment were stretched by shortening the normal regions outside the jet. During relaxation, the sarcomeres in the exposed segment shortened rapidly. Short trains of stimulation at 2.5 Hz reproducibly caused Ca(2+) waves to rise from the borders exposed to the jet. Ca(2+) waves started during F relaxation of the last stimulated twitch and propagated into segments both inside and outside of the jet. Arrhythmias, in the form of nondriven rhythmic activity, were triggered when the amplitude of the Ca(2+) wave increased by raising [Ca(2+)](o). The arrhythmias disappeared when the muscle uniformity was restored by turning the jet off. These results show that nonuniform contraction can cause Ca(2+) waves underlying TPCs, and suggest that Ca(2+) dissociated from myofilaments plays an important role in the initiation of arrhythmogenic Ca(2+) waves.
Collapse
Affiliation(s)
- Henk E D J Ter Keurs
- Department of Medicine, Health Sciences Centre, University of Calgary, Calgary, Alberta T2N 4N1, Canada.
| | | | | | | | | | | |
Collapse
|
40
|
Shannon TR, Bers DM, Blatter LA, Niggli E. Confocal imaging of CICR events from isolated and immobilized SR vesicles. Cell Calcium 2005; 38:497-505. [PMID: 16122794 DOI: 10.1016/j.ceca.2005.07.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2004] [Revised: 05/12/2005] [Accepted: 07/01/2005] [Indexed: 10/25/2022]
Abstract
The Ca2+ concentration inside the sarcoplasmic reticulum ([Ca2+]SR) is a difficult parameter to measure in ventricular cardiac myocytes. Interference from Ca2+-sensitive dye loading into cellular compartments other than the SR interferes with free Ca2+ measurement. In addition, the composition of the cytosol surrounding the SR in intact cells cannot be easily controlled. We have developed a method to measure localized [Ca2+]SR in immobilized membrane vesicles during rapid solution switches. Ca2+ uptake and release in rat SR membrane vesicles was monitored using confocal microscopy. Vesicles were immobilized on a coverslip using an agarose matrix. Perfusion with a Ca2+-containing solution supplemented with ATP initiated SR Ca2+ uptake, causing a rise in intravesicular fluorescence in vesicles containing the low-affinity Ca2+ indicator fluo-5N. Perfusion with caffeine caused SR Ca2+ release and a decrease in intravesicular flourescence. Although caffeine-dependent release was readily visible with extravesicular Ca2+-green, Ca2+ which leaked from the SR was detected only indirectly as eventless release. We conclude that SR Ca2+ uptake and release can be selectively measured in functional SR vesicles using a confocal microscope. Caffeine-dependent release is directly measurable though SR Ca2+ leak can only be inferred as subresolution events, presumably because channels in separate vesicles were not close enough to result in concerted Ca2+-induced Ca2+ release.
Collapse
Affiliation(s)
- T R Shannon
- Department of Physiology, Loyola University, 2160 South First Avenue, Maywood, IL 60153, USA
| | | | | | | |
Collapse
|
41
|
McCulloch AD. Functionally and structurally integrated computational modeling of ventricular physiology. ACTA ACUST UNITED AC 2005; 54:531-9. [PMID: 15760485 DOI: 10.2170/jjphysiol.54.531] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Computational biology is integrative in several ways. Functionally, computational models are valuable for integrating the many interacting processes within biochemical networks and the many interacting physiological subsystems within the cell. Structurally detailed models provide a way of integrating across scales of biological organization from molecule to organism. Data integration across diverse laboratory and clinical measurements is another unique strength of computational biology. We describe examples of all three categories of integration by using recent advances in modeling cardiac excitation-contraction coupling and whole-heart electromechanics in health and disease.
Collapse
Affiliation(s)
- Andrew D McCulloch
- Department of Bioengineering, Whitaker Institute for Biomedical Engineering, University of California San Diego, La Jolla, CA 92093-0412, USA.
| |
Collapse
|
42
|
Mackenzie L, Roderick HL, Berridge MJ, Conway SJ, Bootman MD. The spatial pattern of atrial cardiomyocyte calcium signalling modulates contraction. J Cell Sci 2004; 117:6327-37. [PMID: 15561771 DOI: 10.1242/jcs.01559] [Citation(s) in RCA: 116] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
We examined the regulation of calcium signalling in atrial cardiomyocytes during excitation-contraction coupling, and how changes in the distribution of calcium impacts on contractility. Under control conditions, calcium transients originated in subsarcolemmal locations and showed local regeneration through activation of calcium-induced calcium release from ryanodine receptors. Despite functional ryanodine receptors being expressed at regular (∼2 μm) intervals throughout atrial myocytes, the subsarcolemmal calcium signal did not spread in a fully regenerative manner through the interior of a cell. Rather, there was a diminishing centripetal propagation of calcium. The lack of regeneration was due to mitochondria and SERCA pumps preventing the inward movement of calcium. Inhibiting these calcium buffering mechanisms allowed the globalisation of action potential-evoked responses. In addition, physiological positive inotropic agents, such as endothelin-1 and β-adrenergic agonists, as well as enhanced calcium current, calcium store loading and inositol 1,4,5-trisphosphate infusion also led to regenerative global responses. The consequence of globalising calcium signals was a significant increase in cellular contraction. These data indicate how calcium signals and their consequences are determined by the interplay of multiple subcellular calcium management systems.
Collapse
Affiliation(s)
- Lauren Mackenzie
- Laboratory of Molecular Signalling, The Babraham Institute, Babraham Hall, Babraham, Cambridge, CB2 4AT, UK
| | | | | | | | | |
Collapse
|
43
|
Michailova AP, Belik ME, McCulloch AD. Effects of Magnesium on Cardiac Excitation-Contraction Coupling. J Am Coll Nutr 2004; 23:514S-517S. [PMID: 15466954 DOI: 10.1080/07315724.2004.10719392] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
OBJECTIVE Magnesium regulates a large number of cellular processes. Small changes in intracellular free Mg(2+) ([Mg(2+)](i)) may have important effects on cardiac excitability and contractility. We investigated the effects of [Mg(2+)](i) on cardiac excitation-contraction coupling. METHODS We used our ionic-metabolic model that incorporates equations for Ca(2+) and Mg(2+) buffering and transport by ATP and ADP and equations for MgATP regulation of ion transporters (Na(+)-K(+) pump, sarcolemmal and sarcoplasmic Ca(2+) pumps). RESULTS Model results indicate that variations in cytosolic Mg(2+) level might sensitively affect diastolic and systolic Ca(2+), sarcoplasmic Ca(2+) content, Ca(2+) influx through L-type channels, efficiency of the Na(+)/Ca(2+) exchanger and action potential shape. The analysis suggests that the most important reason for the observed effects is a modified normal function of sarcoplasmic Ca(2+)-ATPase pump by altered diastolic MgATP levels. CONCLUSION The model is able to reproduce qualitatively a sequence of events that correspond well with experimental observations during cardiac excitation-contraction coupling in mammalian ventricular myocytes.
Collapse
Affiliation(s)
- Anushka P Michailova
- Department of Bioengineering, PFBH 241, University of California-San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0412, USA.
| | | | | |
Collapse
|
44
|
Puglisi JL, Wang F, Bers DM. Modeling the isolated cardiac myocyte. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2004; 85:163-78. [PMID: 15142742 DOI: 10.1016/j.pbiomolbio.2003.12.003] [Citation(s) in RCA: 38] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Computer modeling of cardiac myocytes has flourished in recent years. Models have evolved from mathematical descriptions of ionic channels alone to more sophisticated formulations that include calcium transport mechanisms, ATP production and metabolic pathways. The increased complexity is fueled by the new data available in the field. The continuous production of experimental data has led to the evolution of increasingly refined descriptions of the phenomena by modelers. Integrating the numerous systems involved in cardiac myocyte homeostasis makes the use of computer models necessary due to the unreliability of intuitive approaches. However the complexity of the model should not imply a cumbersome operation of the program. As with any tool, computer models have to be easy to operate or their strength will be diminished and potential users will not benefit fully from them. The contribution of the computer modeler to their respective biological fields will be more successful and enduring if modelers devote sufficient time to implement their equations into a model with user-friendly characteristics.
Collapse
Affiliation(s)
- Jose L Puglisi
- Department of Physiology, Loyola University Chicago, 2160 S. First Avenue, Maywood, IL 60153, USA.
| | | | | |
Collapse
|
45
|
Strier DE, Dawson SP. Role of complexing agents in the appearance of Turing patterns. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 69:066207. [PMID: 15244708 DOI: 10.1103/physreve.69.066207] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2003] [Indexed: 05/24/2023]
Abstract
In this paper we study a four-species reaction-diffusion system where Turing patterns are stabilized by the presence of fast reversible reactions between the morphogens and two different mobile complexing agents (CAs) that are not necessarily in excess. We provide a quantitative explanation of how the interaction with the CA changes the size of the Turing space making it possible to observe patterns even in a region where the free diffusion coefficients of the relevant species are equal, as is usually the case in real systems. Our analytical treatment gives a series of mathematical relations that can be helpful for those designing experiments where Turing patterns are expected to appear. We also show how the mobility of CAs affect the characteristic size of the pattern. Finally, we provide an example of biological interest in order to illustrate the main procedures and results.
Collapse
Affiliation(s)
- Damián E Strier
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, UBA, Ciudad Universitaria, Pabellón I, (1428) Buenos Aires, Argentina.
| | | |
Collapse
|
46
|
McHugh JM, Kenyon JL. An Excel-based model of Ca2+diffusion and fura 2 measurements in a spherical cell. Am J Physiol Cell Physiol 2004; 286:C342-8. [PMID: 14512292 DOI: 10.1152/ajpcell.00270.2003] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
We wrote a program that runs as a Microsoft Excel spreadsheet to calculate the diffusion of Ca2+in a spherical cell in the presence of a fixed Ca2+buffer and two diffusible Ca2+buffers, one of which is considered to be a fluorescent Ca2+indicator. We modeled Ca2+diffusion during and after Ca2+influx across the plasma membrane with parameters chosen to approximate amphibian sympathetic neurons, mammalian adrenal chromaffin cells, and rat dorsal root ganglion neurons. In each of these cell types, the model predicts that spatially averaged intracellular Ca2+activity ([Ca2+]avg) rises to a high peak and starts to decline promptly on the termination of Ca2+influx. We compared [Ca2+]avgwith predictions of ratiometric Ca2+measurements analyzed in two ways. Method 1 sums the fluorescence at each of the two excitation or emission wavelengths over the N compartments of the model, calculates the ratio of the summed signals, and converts this ratio to Ca2+([Ca2+]avg,M1). Method 2 sums the measured number of moles of Ca2+in each of the N compartments and divides by the volume of the cell ([Ca2+]avg,M2). [Ca2+]avg,M1peaks well after the termination of Ca2+influx at a value substantially less than [Ca2+]avgbecause the summed signals do not reflect the averaged free Ca2+if the signals come from compartments containing gradients in free Ca2+spanning nonlinear regions of the relationship between free Ca2+and the fluorescence signals. In contrast, [Ca2+]avg,M2follows [Ca2+]avgclosely.
Collapse
Affiliation(s)
- J M McHugh
- Department of Physiology & Cell Biology, University of Nevada School of Medicine, Reno, NV 89557, USA
| | | |
Collapse
|
47
|
Abstract
Dendritic spines receive excitatory synapses and serve as calcium compartments, which appear to be necessary for input-specific synaptic plasticity. Dendrites of GABAergic interneurons have few or no spines and thus do not possess a clear morphological basis for synapse-specific compartmentalization. We demonstrate using two-photon calcium imaging that activation of single synapses on aspiny dendrites of neocortical fast spiking (FS) interneurons creates highly localized calcium microdomains, often restricted to less than 1 microm of dendritic space. We confirm using ultrastructural reconstruction of imaged dendrites the absence of any morphological basis for this compartmentalization and show that it is dependent on the fast kinetics of calcium-permeable (CP) AMPA receptors and fast local extrusion via the Na+/Ca2+ exchanger. Because aspiny dendrites throughout the CNS express CP-AMPA receptors, we propose that CP-AMPA receptors mediate a spine-free mechanism of input-specific calcium compartmentalization.
Collapse
Affiliation(s)
- Jesse H Goldberg
- Department of Biological Sciences, Columbia University, New York, NY 10027, USA.
| | | | | | | |
Collapse
|