Calcium Ion Fluctuations Alter Channel Gating in a Stochastic Luminal Calcium Release Site Model

Calcium Ion Fluctuations Alter Channel Gating in a Stochastic Luminal Calcium Release Site Model

Stochasticity and small system size effects in complex biochemical reaction networks can greatly alter transient and steady-state system properties. A common approach to modeling reaction networks, which accounts for system size, is the chemical master equation that governs the dynamics of the joint probability distribution for molecular copy number. However, calculation of the stationary distribution is often prohibitive, due to the large state-space associated with most biochemical reaction networks. Here, we analyze a network representing a luminal calcium release site model and investigate to what extent small system size effects and calcium fluctuations, driven by ion channel gating, influx and diffusion, alter steady-state ion channel properties including open probability. For a physiological ion channel gating model and number of channels, the state-space may be between approximately 10⁶-10⁸ elements, and a novel modified block power method is used to solve the associated dominant eigenvector problem required to calculate the stationary distribution. We demonstrate that both small local cytosolic domain volume and a small number of ion channels drive calcium fluctuations that result in deviation from the corresponding model that neglects small system size effects.

Role of Cytosolic Calcium Diffusion in Murine Cardiac Purkinje Cells

Cardiac Purkinje cells (PCs) are morphologically and electrophysiologically different from ventricular myocytes and, importantly, exhibit distinct calcium (Ca²⁺) homeostasis. Recent studies suggest that PCs are more susceptible to action potential (AP) abnormalities than ventricular myocytes; however, the exact mechanisms are poorly understood. In this study, we utilized a detailed biophysical mathematical model of a murine PC to systematically examine the role of cytosolic Ca²⁺ diffusion in shaping the AP in PCs. A biphasic spatiotemporal Ca²⁺ diffusion process, as recorded experimentally, was implemented in the model. In this study, we investigated the role of cytosolic Ca²⁺ dynamics on AP and ionic current properties by varying the effective Ca²⁺ diffusion rate. It was observed that AP morphology, specifically the plateau, was affected due to changes in the intracellular Ca²⁺ dynamics. Elevated Ca²⁺ concentration in the sarcolemmal region activated inward sodium–Ca²⁺ exchanger (NCX) current, resulting in a prolongation of the AP plateau at faster diffusion rates. Artificially clamping the NCX current to control values completely reversed the alterations in the AP plateau, thus confirming the role of NCX in modifying the AP morphology. Our results demonstrate that cytosolic Ca²⁺ diffusion waves play a significant role in shaping APs of PCs and could provide mechanistic insights in the increased arrhythmogeneity of PCs.

Impaired Sarcoplasmic Reticulum Calcium Uptake and Release Promote Electromechanically and Spatially Discordant Alternans: A Computational Study

Cardiac electrical dynamics are governed by cellular-level properties, such as action potential duration (APD) restitution and intracellular calcium (Ca) handling, and tissue-level properties, including conduction velocity restitution and cell-cell coupling. Irregular dynamics at the cellular level can lead to instabilities in cardiac tissue, including alternans, a beat-to-beat alternation in the action potential and/or the intracellular Ca transient. In this study, we incorporate a detailed single cell coupled map model of Ca cycling and bidirectional APD-Ca coupling into a spatially extended tissue model to investigate the influence of sarcoplasmic reticulum (SR) Ca uptake and release properties on alternans and conduction block. We find that an intermediate SR Ca uptake rate and larger SR Ca release resulted in the widest range of stimulus periods that promoted alternans. However, both reduced SR Ca uptake and release promote arrhythmogenic spatially and electromechanically discordant alternans, suggesting a complex interaction between SR Ca handling and alternans characteristics at the cellular and tissue level.