1
|
Veron G, Maltsev VA, Stern MD, Maltsev AV. Elementary intracellular Ca signals approximated as a transition of release channel system from a metastable state. JOURNAL OF APPLIED PHYSICS 2023; 134:124701. [PMID: 37744735 PMCID: PMC10517864 DOI: 10.1063/5.0151255] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 09/04/2023] [Indexed: 09/26/2023]
Abstract
Cardiac muscle contraction is initiated by an elementary Ca signal (called Ca spark) which is achieved by collective action of Ca release channels in a cluster. The mechanism of this synchronization remains uncertain. We approached Ca spark activation as an emergent phenomenon of an interactive system of release channels. We constructed a weakly lumped Markov chain that applies an Ising model formalism to such release channel clusters and probable open channel configurations and demonstrated that spark activation is described as a system transition from a metastable to an absorbing state, analogous to the pressure required to overcome surface tension in bubble formation. This yielded quantitative estimates of the spark generation probability as a function of various system parameters. We performed numerical simulations to find spark probabilities as a function of sarcoplasmic reticulum Ca concentration, obtaining similar values for spark activation threshold as our analytic model, as well as those reported in experimental studies. Our parametric sensitivity analyses also showed that the spark activation threshold decreased as Ca sensitivity of RyR activation and RyR cluster size increased.
Collapse
Affiliation(s)
- Guillermo Veron
- Cellular Biophysics Section, Laboratory of Cardiovascular Science, National Institute on Aging, NIH, Baltimore, Maryland 21224, USA
| | - Victor A. Maltsev
- Cellular Biophysics Section, Laboratory of Cardiovascular Science, National Institute on Aging, NIH, Baltimore, Maryland 21224, USA
| | - Michael D. Stern
- Cellular Biophysics Section, Laboratory of Cardiovascular Science, National Institute on Aging, NIH, Baltimore, Maryland 21224, USA
| | - Anna V. Maltsev
- School of Mathematical Sciences, Queen Mary University of London, London E14NS, United Kingdom
| |
Collapse
|
2
|
Maltsev AV, Stern MD, Maltsev VA. Disorder in Ca2+ release unit locations confers robustness but cuts flexibility of heart pacemaking. J Gen Physiol 2022; 154:213390. [PMID: 35943725 PMCID: PMC9366202 DOI: 10.1085/jgp.202113061] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Revised: 05/04/2022] [Accepted: 06/21/2022] [Indexed: 12/12/2022] Open
Abstract
Excitation-contraction coupling kinetics is dictated by the action potential rate of sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the pace. LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While CRU distribution exhibits substantial heterogeneity, its functional importance remains unknown. Using numerical modeling, here we show that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad range of realistic ICaL densities. Allowing each CRU to deviate randomly from its lattice position allows sparks to propagate, as observed experimentally. As disorder increases, the CRU distribution exhibits larger empty spaces and simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing action potential rate and reviving pacemaker function of dormant/nonfiring cells. However, cells with fully disordered CRU positions could not reach low firing rates and their β-adrenergic-receptor stimulation effect was substantially decreased. Inclusion of Cav1.3, a low-voltage activation L-type Ca channel isoform into ICaL, strongly increases recruitment of CRUs to fire during diastolic depolarization, increasing robustness of pacemaking and complementing effects of CRU distribution. Thus, order/disorder in CRU locations along with Cav1.3 expression regulates pacemaker function via synchronization of CRU firing. Excessive CRU disorder and/or overexpression of Cav1.3 boosts pacemaker function in the basal state, but limits the rate range, which may contribute to heart rate range decline with age and disease.
Collapse
Affiliation(s)
- Anna V. Maltsev
- School of Mathematics, Queen Mary University of London, London, UK,Correspondence to Anna V. Maltsev:
| | - Michael D. Stern
- Laboratory of Cardiovascular Science, Intramural Research Program, National Institute on Aging, National Institutes of Health, Baltimore, MD
| | - Victor A. Maltsev
- Laboratory of Cardiovascular Science, Intramural Research Program, National Institute on Aging, National Institutes of Health, Baltimore, MD,Victor A. Maltsev:
| |
Collapse
|
3
|
Yang D, Morrell CH, Lyashkov AE, Tagirova Sirenko S, Zahanich I, Yaniv Y, Vinogradova TM, Ziman BD, Maltsev VA, Lakatta EG. Ca 2+ and Membrane Potential Transitions During Action Potentials Are Self-Similar to Each Other and to Variability of AP Firing Intervals Across the Broad Physiologic Range of AP Intervals During Autonomic Receptor Stimulation. Front Physiol 2021; 12:612770. [PMID: 34566668 PMCID: PMC8456031 DOI: 10.3389/fphys.2021.612770] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 06/02/2021] [Indexed: 12/02/2022] Open
Abstract
Ca2+ and V m transitions occurring throughout action potential (AP) cycles in sinoatrial nodal (SAN) cells are cues that (1) not only regulate activation states of molecules operating within criticality (Ca2+ domain) and limit-cycle (V m domain) mechanisms of a coupled-clock system that underlies SAN cell automaticity, (2) but are also regulated by the activation states of the clock molecules they regulate. In other terms, these cues are both causes and effects of clock molecular activation (recursion). Recently, we demonstrated that Ca2+ and V m transitions during AP cycles in single SAN cells isolated from mice, guinea pigs, rabbits, and humans are self-similar (obey a power law) and are also self-similar to trans-species AP firing intervals (APFIs) of these cells in vitro, to heart rate in vivo, and to body mass. Neurotransmitter stimulation of β-adrenergic receptor or cholinergic receptor-initiated signaling in SAN cells modulates their AP firing rate and rhythm by impacting on the degree to which SAN clocks couple to each other, creating the broad physiologic range of SAN cell mean APFIs and firing interval variabilities. Here we show that Ca2+ and V m domain kinetic transitions (time to AP ignition in diastole and 90% AP recovery) occurring within given AP, the mean APFIs, and APFI variabilities within the time series of APs in 230 individual SAN cells are self-similar (obey power laws). In other terms, these long-range correlations inform on self-similar distributions of order among SAN cells across the entire broad physiologic range of SAN APFIs, regardless of whether autonomic receptors of these cells are stimulated or not and regardless of the type (adrenergic or cholinergic) of autonomic receptor stimulation. These long-range correlations among distributions of Ca2+ and V m kinetic functions that regulate SAN cell clock coupling during each AP cycle in different individual, isolated SAN cells not in contact with each other. Our numerical model simulations further extended our perspectives to the molecular scale and demonstrated that many ion currents also behave self-similar across autonomic states. Thus, to ensure rapid flexibility of AP firing rates in response to different types and degrees of autonomic input, nature "did not reinvent molecular wheels within the coupled-clock system of pacemaker cells," but differentially engaged or scaled the kinetics of gears that regulate the rate and rhythm at which the "wheels spin" in a given autonomic input context.
Collapse
Affiliation(s)
- Dongmei Yang
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Christopher H. Morrell
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
- Department of Mathematics and Statistics, Loyola University Maryland, Baltimore, MD, United States
| | - Alexey E. Lyashkov
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Syevda Tagirova Sirenko
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Ihor Zahanich
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Yael Yaniv
- Biomedical Engineering Faculty, Technion–Israel Institute of Technology, Haifa, Israel
| | - Tatiana M. Vinogradova
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Bruce D. Ziman
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Victor A. Maltsev
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| | - Edward G. Lakatta
- Laboratory of Cardiovascular Science, National Institute on Aging, National Institutes of Health, Baltimore, MD, United States
| |
Collapse
|
4
|
Baier MJ, Noack J, Seitz MT, Maier LS, Neef S. Phosphorylation of RyR2 Ser-2814 by CaMKII mediates β1-adrenergic stress induced Ca 2+ -leak from the sarcoplasmic reticulum. FEBS Open Bio 2021; 11:2756-2762. [PMID: 34403217 PMCID: PMC8487045 DOI: 10.1002/2211-5463.13274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 08/03/2021] [Accepted: 08/16/2021] [Indexed: 11/26/2022] Open
Abstract
Adrenergic stimulation, while being the central mechanism of cardiac positive inotropy, is a universally acknowledged inductor of undesirable sarcoplasmic reticulum (SR) Ca2+ leak. However, the exact mechanisms for this remained unspecified so far. This study shows that Ca2+/calmodulin‐dependent protein kinase II (CaMKII)‐specific phosphorylation of ryanodine receptor type 2 at Ser‐2814 is the pivotal mechanism by which SR Ca2+ leak develops downstream of β1‐adrenergic stress by increase of the leak/load relationship. Cardiomyocytes with a Ser‐2814 phosphoresistant mutation (S2814A) were protected from isoproterenol‐induced SR Ca2+ leak and consequently displayed improved postrest potentiation of systolic Ca2+ release under adrenergic stress compared to littermate wild‐type cells.
Collapse
Affiliation(s)
- Maria J Baier
- Department of Internal Medicine II, University Medical Center Regensburg, Germany
| | - Jannis Noack
- Department of Internal Medicine II, University Medical Center Regensburg, Germany
| | - Mark Tilmann Seitz
- Department of Trauma Surgery, Orthopaedics and Plastic Surgery, University Medical Center Göttingen, Germany
| | - Lars S Maier
- Department of Internal Medicine II, University Medical Center Regensburg, Germany
| | - Stefan Neef
- Department of Internal Medicine II, University Medical Center Regensburg, Germany
| |
Collapse
|
5
|
Iaparov BI, Zahradnik I, Moskvin AS, Zahradníková A. In silico simulations reveal that RYR distribution affects the dynamics of calcium release in cardiac myocytes. J Gen Physiol 2021; 153:211900. [PMID: 33735373 PMCID: PMC7980188 DOI: 10.1085/jgp.202012685] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2020] [Accepted: 02/08/2021] [Indexed: 11/20/2022] Open
Abstract
The dyads of cardiac myocytes contain ryanodine receptors (RYRs) that generate calcium sparks upon activation. To test how geometric factors of RYR distribution contribute to the formation of calcium sparks, which cannot be addressed experimentally, we performed in silico simulations on a large set of models of calcium release sites (CRSs). Our models covered the observed range of RYR number, density, and spatial arrangement. The calcium release function of CRSs was modeled by RYR openings, with an open probability dependent on concentrations of free Ca2+ and Mg2+ ions, in a rapidly buffered system, with a constant open RYR calcium current. We found that simulations of spontaneous sparks by repeatedly opening one of the RYRs in a CRS produced three different types of calcium release events (CREs) in any of the models. Transformation of simulated CREs into fluorescence signals yielded calcium sparks with characteristics close to the observed ones. CRE occurrence varied broadly with the spatial distribution of RYRs in the CRS but did not consistently correlate with RYR number, surface density, or calcium current. However, it correlated with RYR coupling strength, defined as the weighted product of RYR vicinity and calcium current, so that CRE characteristics of all models followed the same state-response function. This finding revealed the synergy between structure and function of CRSs in shaping dyad function. Lastly, rearrangements of RYRs simulating hypothetical experiments on splitting and compaction of a dyad revealed an increased propensity to generate spontaneous sparks and an overall increase in calcium release in smaller and more compact dyads, thus underlying the importance and physiological role of RYR arrangement in cardiac myocytes.
Collapse
Affiliation(s)
- Bogdan I Iaparov
- Department of Cellular Cardiology, Institute of Experimental Endocrinology, Biomedical Research Center, Slovak Academy of Sciences, Bratislava, Slovakia.,Research Institute of Physics and Applied Mathematics, and Department of Theoretical and Mathematical Physics, Ural Federal University, Ekaterinburg, Russia
| | - Ivan Zahradnik
- Department of Cellular Cardiology, Institute of Experimental Endocrinology, Biomedical Research Center, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Alexander S Moskvin
- Research Institute of Physics and Applied Mathematics, and Department of Theoretical and Mathematical Physics, Ural Federal University, Ekaterinburg, Russia
| | - Alexandra Zahradníková
- Department of Cellular Cardiology, Institute of Experimental Endocrinology, Biomedical Research Center, Slovak Academy of Sciences, Bratislava, Slovakia
| |
Collapse
|
6
|
Fill M, Gillespie D. Simulating cardiac Ca 2+ release units: effects of RyR cluster size and Ca 2+ buffers on diastolic Ca 2+ leak. Pflugers Arch 2021; 473:435-446. [PMID: 33608799 DOI: 10.1007/s00424-021-02539-w] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Revised: 01/27/2021] [Accepted: 02/05/2021] [Indexed: 10/22/2022]
Abstract
Leak of Ca2+ out of the cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) during diastole is vital to regulate SR Ca2+ levels. This leak can become deleterious when large spontaneous RyR-mediated Ca2+ release events evoke proarrhythmic Ca2+ waves that can lead to delayed after-depolarizations. Here, we model diastolic SR Ca2+ leak at individual SR Ca2+ release sites using computer simulations of RyR arrays like those in the dyadic cleft. The results show that RyR arrays size has a significant effect on SR Ca2+ leak, with bigger arrays producing larger and more frequent Ca2+ release events. Moreover, big RyR arrays are more susceptible to small changes in the levels of dyadic Ca2+ buffers. Such changes in buffering shift Ca2+ leak from small Ca2+ release events (involving few open RyRs) to larger events (with many open RyRs). Moreover, by analyzing a large parameter space of possible buffering and SR Ca2+ loads, we find further evidence for the hypothesis that SR Ca2+ leak by RyR arrays can undergo a sudden phase transition.
Collapse
Affiliation(s)
- Michael Fill
- Department of Physiology and Biophysics, Rush University Medical Center, Chicago, IL, USA
| | - Dirk Gillespie
- Department of Physiology and Biophysics, Rush University Medical Center, Chicago, IL, USA.
| |
Collapse
|
7
|
Gillespie D. Recruiting RyRs to Open in a Ca 2+ Release Unit: Single-RyR Gating Properties Make RyR Group Dynamics. Biophys J 2019; 118:232-242. [PMID: 31839264 DOI: 10.1016/j.bpj.2019.11.021] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 11/01/2019] [Accepted: 11/19/2019] [Indexed: 01/01/2023] Open
Abstract
In cardiac myocytes, clusters of type-2 ryanodine receptors (RyR2s) release Ca2+ from the sarcoplasmic reticulum (SR) via a positive feedback mechanism in which fluxed Ca2+ activates nearby RyRs. Although the general principles of this are understood, less is known about how single-RyR gating properties define the RyR group dynamics in an array of many channels. Here, we examine this using simulations with three models of RyR gating that have identical open probabilities: the commonly used two-state Markov gating model, one that utilizes multiple exponentials to fit single-channel open time (OT) and closed time (CT) distributions, and an extension of this multiexponential model that also includes experimentally measured correlations between single-channel OTs and CTs. The simulations of RyR clusters that utilize the multiexponential gating model produce infrequent Ca2+ release events with relatively few open RyRs. Ca2+ release events become even smaller when OT/CT correlations are included. This occurs because the correlations produce a small but consistent bias against recruiting more RyRs to open during the middle of a Ca2+ release event, between the initiation and termination phases (which are unaltered compared to the uncorrelated simulations). In comparison, the two-state model produces frequent, large, and long Ca2+ release events because it had a recruitment bias in favor of opening more RyRs. This difference stems from the two-state model's single-RyR OT and CT distributions being qualitatively different from the experimental ones. Thus, the details of single-RyR gating can profoundly affect SR Ca2+ release even if open probability and mean OTs and CTs are identical. We also show that Ca2+ release events can terminate spontaneously without any reduction in SR [Ca2+], luminal regulation, Ca2+-dependent inactivation, or physical coupling between RyRs when Ca2+ flux is below a threshold value. This supports and extends the pernicious attrition/induction decay hypothesis that SR Ca2+ release events terminate below a threshold Ca2+ flux.
Collapse
Affiliation(s)
- Dirk Gillespie
- Department of Physiology and Biophysics, Rush University Medical Center, Chicago, Illinois.
| |
Collapse
|