Transverse tubular network structures in the genesis of intracellular calcium alternans and triggered activity in cardiac cells

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Mechanistic insights into spontaneous transition from cellular alternans to ventricular fibrillation

T-wave alternans (TWA) has been used for predicting the risk of malignant cardiac arrhythmias and sudden cardiac death (SCD) in multiple clinical settings; however, possible mechanism(s) underlying the spontaneous transition from cellular alternans reflected by TWA to arrhythmias in impaired repolarization remains unclear. The healthy guinea pig ventricular myocytes under E-4031 blocking IKr (0.1 μM, N = 12; 0.3 μM, N = 10; 1 μM, N = 10) were evaluated using whole-cell patch-clamp. The electrophysiological properties of isolated perfused guinea pig hearts under E-4031 (0.1 μM, N = 5; 0.3 μM, N = 5; 1 μM, N = 5) were evaluated using dual- optical mapping. The amplitude/threshold/restitution curves of action potential duration (APD) alternans and potential mechanism(s) underlying the spontaneous transition of cellular alternans to ventricular fibrillation (VF) were examined. There were longer APD80 and increased amplitude and threshold of APD alternans in E-4031 group compared with baseline group, which was reflected by more pronounced arrhythmogenesis at the tissue level, and were associated with steep restitution curves of the APD and the conduction velocity (CV). Conduction of AP alternans augmented tissue's functional spatiotemporal heterogeneity of regional AP/Ca alternans, as well as the AP/Ca dispersion, leading to localized uni-directional conduction block that spontaneous facilitated the formation of reentrant excitation waves without the need for additional premature stimulus. Our results provide a possible mechanism for the spontaneous transition from cardiac electrical alternans in cellular action potentials and intercellular conduction without the involvement of premature excitations, and explain the increased susceptibility to ventricular arrhythmias in impaired repolarization. In this study, we implemented voltage-clamp and dual-optical mapping approaches to investigate the underlying mechanism(s) for the arrhythmogenesis of cardiac alternans in the guinea pig heart at cellular and tissue levels. Our results demonstrated a spontaneous development of reentry from cellular alternans, arising from a combined actions of restitution properties of action potential duration, conduction velocity of excitation wave and interplay between alternants of action potential and the intracellular Ca handling. We believe this study provides new insights into underlying the mechanism, by which cellular cardiac alternans spontaneously evolves into cardiac arrhythmias.

Cardiac Alternans: From Bedside to Bench and Back

Cardiac alternans arises from dynamical instabilities in the electrical and calcium cycling systems of the heart, and often precedes ventricular arrhythmias and sudden cardiac death. In this review, we integrate clinical observations with theory and experiment to paint a holistic portrait of cardiac alternans: the underlying mechanisms, arrhythmic manifestations and electrocardiographic signatures. We first summarize the cellular and tissue mechanisms of alternans that have been demonstrated both theoretically and experimentally, including 3 voltage-driven and 2 calcium-driven alternans mechanisms. Based on experimental and simulation results, we describe their relevance to mechanisms of arrhythmogenesis under different disease conditions, and their link to electrocardiographic characteristics of alternans observed in patients. Our major conclusion is that alternans is not only a predictor, but also a causal mechanism of potentially lethal ventricular and atrial arrhythmias across the full spectrum of arrhythmia mechanisms that culminate in functional reentry, although less important for anatomic reentry and focal arrhythmias.

Increased SERCA2a sub-cellular heterogeneity in right-ventricular heart failure inhibits excitation-contraction coupling and modulates arrhythmogenic dynamics

The intracellular calcium handling system of cardiomyocytes is responsible for controlling excitation-contraction coupling (ECC) and has been linked to pro-arrhythmogenic cellular phenomena in conditions such as heart failure (HF). SERCA2a, responsible for intracellular uptake, is a primary regulator of calcium homeostasis, and remodelling of its function has been proposed as a causal factor underlying cellular and tissue dysfunction in disease. Whereas adaptations to the global (i.e. whole-cell) expression of SERCA2a have been previously investigated in the context of multiple diseases, the role of its spatial profile in the sub-cellular volume has yet to be elucidated. We present an approach to characterize the sub-cellular heterogeneity of SERCA2a and apply this approach to quantify adaptations to the length-scale of heterogeneity (the distance over which expression is correlated) associated with right-ventricular (RV)-HF. These characterizations informed simulations to predict the functional implications of this heterogeneity, and its remodelling in disease, on ECC, the dynamics of calcium-transient alternans and the emergence of spontaneous triggered activity. Image analysis reveals that RV-HF is associated with an increase in length-scale and its inter-cellular variability; simulations predict that this increase in length-scale can reduce ECC and critically modulate the vulnerability to both alternans and triggered activity. This article is part of the theme issue 'The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease'.

Modeling Calcium Cycling in the Heart: Progress, Pitfalls, and Challenges

Intracellular calcium (Ca) cycling in the heart plays key roles in excitation-contraction coupling and arrhythmogenesis. In cardiac myocytes, the Ca release channels, i.e., the ryanodine receptors (RyRs), are clustered in the sarcoplasmic reticulum membrane, forming Ca release units (CRUs). The RyRs in a CRU act collectively to give rise to discrete Ca release events, called Ca sparks. A cell contains hundreds to thousands of CRUs, diffusively coupled via Ca to form a CRU network. A rich spectrum of spatiotemporal Ca dynamics is observed in cardiac myocytes, including Ca sparks, spark clusters, mini-waves, persistent whole-cell waves, and oscillations. Models of different temporal and spatial scales have been developed to investigate these dynamics. Due to the complexities of the CRU network and the spatiotemporal Ca dynamics, it is challenging to model the Ca cycling dynamics in the cardiac system, particularly at the tissue sales. In this article, we review the progress of modeling of Ca cycling in cardiac systems from single RyRs to the tissue scale, the pros and cons of the current models and different modeling approaches, and the challenges to be tackled in the future.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca²⁺ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca²⁺ -handling proteins such as the L-type Ca²⁺ channel and Na⁺ -Ca²⁺ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca²⁺ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca²⁺ -handling and Ca²⁺ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca²⁺ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca²⁺ release. We found that varying TATS density and thus the localization of key Ca²⁺ -handling proteins has profound effects on Ca²⁺ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca²⁺ concentration through decreased Ca²⁺ extrusion. This elevated Ca²⁺ increases RyR open probability causing spontaneous Ca²⁺ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca²⁺ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca²⁺ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca²⁺ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca²⁺ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca²⁺ handling governed by the TATS. Simulated TATS loss causes diastolic Ca²⁺ and voltage instabilities through reduced Na⁺ -Ca²⁺ exchanger-mediated Ca²⁺ removal, cleft Ca²⁺ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca²⁺ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca²⁺ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca²⁺ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

Synchronization of spatially discordant voltage and calcium alternans in cardiac tissue

The heart is an excitable medium which is excited by membrane potential depolarization and propagation. Membrane potential depolarization brings in calcium (Ca) through the Ca channels to trigger intracellular Ca release for contraction of the heart. Ca also affects voltage via Ca-dependent ionic currents, and thus, voltage and Ca are bidirectionally coupled. It has been shown that the voltage subsystem or the Ca subsystem can generate its own dynamical instabilities which are affected by their bidirectional couplings, leading to complex dynamics of action potential and Ca cycling. Moreover, the dynamics become spatiotemporal in tissue in which cells are diffusively coupled through voltage. A widely investigated spatiotemporal dynamics is spatially discordant alternans (SDA) in which action potential duration (APD) or Ca amplitude exhibits temporally period-2 and spatially out-of-phase patterns, i.e., APD-SDA and Ca-SDA patterns, respectively. However, the mechanisms of formation, stability, and synchronization of APD-SDA and Ca-SDA patterns remain incompletely understood. In this paper, we use cardiac tissue models described by an amplitude equation, coupled iterated maps, and reaction-diffusion equations with detailed physiology (the ionic model) to perform analytical and computational investigations. We show that, when the Ca subsystem is stable, the Ca-SDA pattern always follows the APD-SDA pattern, and thus, they are always synchronized. When the Ca subsystem is unstable, synchronization of APD-SDA and Ca-SDA patterns depends on the stabilities of both subsystems, their coupling strengths, and the spatial scales of the initial Ca-SDA patterns. Spontaneous (initial condition-independent) synchronization is promoted by enhancing APD instability and reducing Ca instability as well as stronger Ca-to-APD and APD-to-Ca coupling, a pattern formation caused by dynamical instabilities. When Ca is more unstable and APD is less unstable or APD-to-Ca coupling is weak, synchronization of APD-SDA and Ca-SDA patterns is promoted by larger initially synchronized Ca-SDA clusters, i.e., initial condition-dependent synchronization. The synchronized APD-SDA and Ca-SDA patterns can be locked in-phase, antiphase, or quasiperiodic depending on the coupling relationship between APD and Ca. These theoretical and simulation results provide mechanistic insights into the APD-SDA and Ca-SDA dynamics observed in experimental studies.

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation

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.

Mitochondrial Contributions in the Genesis of Delayed Afterdepolarizations in Ventricular Myocytes

Mitochondria fulfill the cell's energy demand and affect the intracellular calcium (Ca2+) dynamics via direct Ca2+ exchange, the redox effect of reactive oxygen species (ROS) on Ca2+ handling proteins, and other signaling pathways. Recent experimental evidence indicates that mitochondrial depolarization promotes arrhythmogenic delayed afterdepolarizations (DADs) in cardiac myocytes. However, the nonlinear interactions among the Ca2+ signaling pathways, ROS, and oxidized Ca2+/calmodulin-dependent protein kinase II (CaMKII) pathways make it difficult to reveal the mechanisms. Here, we use a recently developed spatiotemporal ventricular myocyte computer model, which consists of a 3-dimensional network of Ca2+ release units (CRUs) intertwined with mitochondria and integrates mitochondrial Ca2+ signaling and other complex signaling pathways, to study the mitochondrial regulation of DADs. With a systematic investigation of the synergistic or competing factors that affect the occurrence of Ca2+ waves and DADs during mitochondrial depolarization, we find that the direct redox effect of ROS on ryanodine receptors (RyRs) plays a critical role in promoting Ca2+ waves and DADs under the acute effect of mitochondrial depolarization. Furthermore, the upregulation of mitochondrial Ca2+ uniporter can promote DADs through Ca2+-dependent opening of mitochondrial permeability transition pores (mPTPs). Also, due to much slower dynamics than Ca2+ cycling and ROS, oxidized CaMKII activation and the cytosolic ATP do not appear to significantly impact the genesis of DADs during the acute phase of mitochondrial depolarization. However, under chronic conditions, ATP depletion suppresses and enhanced CaMKII activation promotes Ca2+ waves and DADs.

The Physiology and Pathophysiology of T-Tubules in the Heart

In cardiomyocytes, invaginations of the sarcolemmal membrane called t-tubules are critically important for triggering contraction by excitation-contraction (EC) coupling. These structures form functional junctions with the sarcoplasmic reticulum (SR), and thereby enable close contact between L-type Ca²⁺ channels (LTCCs) and Ryanodine Receptors (RyRs). This arrangement in turn ensures efficient triggering of Ca²⁺ release, and contraction. While new data indicate that t-tubules are capable of exhibiting compensatory remodeling, they are also widely reported to be structurally and functionally compromised during disease, resulting in disrupted Ca²⁺ homeostasis, impaired systolic and/or diastolic function, and arrhythmogenesis. This review summarizes these findings, while highlighting an emerging appreciation of the distinct roles of t-tubules in the pathophysiology of heart failure with reduced and preserved ejection fraction (HFrEF and HFpEF). In this context, we review current understanding of the processes underlying t-tubule growth, maintenance, and degradation, underscoring the involvement of a variety of regulatory proteins, including junctophilin-2 (JPH2), amphiphysin-2 (BIN1), caveolin-3 (Cav3), and newer candidate proteins. Upstream regulation of t-tubule structure/function by cardiac workload and specifically ventricular wall stress is also discussed, alongside perspectives for novel strategies which may therapeutically target these mechanisms.

CaMKII inhibition has dual effects on spontaneous Ca2+ release and Ca2+ alternans in ventricular cardiomyocytes from mice with a gain-of-function RyR2 mutation

In conditions with abnormally increased activity of the cardiac ryanodine receptor (RyR2), Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) can contribute to a further destabilization of RyR2 that results in triggered arrhythmias. Therefore, inhibition of CaMKII in such conditions has been suggested as a strategy to suppress RyR2 activity and arrhythmias. However, suppression of RyR2 activity can lead to the development of arrhythmogenic Ca²⁺ alternans. The aim of this study was to test whether the suppression of RyR2 activity caused by inhibition of CaMKII increases propensity for Ca²⁺ alternans. We studied spontaneous Ca²⁺ release events and Ca²⁺ alternans in isolated left ventricular cardiomyocytes from mice carrying the gain-of-function RyR2 mutation RyR2-R2474S and from wild-type mice. CaMKII inhibition by KN-93 effectively decreased the frequency of spontaneous Ca²⁺ release events in RyR2-R2474S cardiomyocytes exposed to the β-adrenoceptor agonist isoprenaline. However, KN-93-treated RyR2-R2474S cardiomyocytes also showed increased propensity for Ca²⁺ alternans and increased Ca²⁺ alternans ratio compared with both an inactive analog of KN-93 and with vehicle-treated controls. This increased propensity for Ca²⁺ alternans was explained by prolongation of Ca²⁺ release refractoriness. Importantly, the increased propensity for Ca²⁺ alternans in KN-93-treated RyR2-R2474S cardiomyocytes did not surpass that of wild type. In conclusion, inhibition of CaMKII efficiently reduces spontaneous Ca²⁺ release but promotes Ca²⁺ alternans in RyR2-R2474S cardiomyocytes with a gain-of-function RyR2 mutation. The dominant effect in RyR2-R2474S is to reduce spontaneous Ca²⁺ release, which supports this intervention as a therapeutic strategy in this specific condition. However, future studies on CaMKII inhibition in conditions with increased propensity for Ca²⁺ alternans should include investigation of both phenomena.NEW & NOTEWORTHY Genetically increased RyR2 activity promotes arrhythmogenic Ca²⁺ release. Inhibition of CaMKII suppresses RyR2 activity and arrhythmogenic Ca²⁺ release. Suppression of RyR2 activity prolongs refractoriness of Ca²⁺ release. Prolonged refractoriness of Ca²⁺ release leads to arrhythmogenic Ca²⁺ alternans. CaMKII inhibition promotes Ca²⁺ alternans by prolonging Ca²⁺ release refractoriness.

Electrophysiological Mechanisms Underlying T-Wave Alternans and Their Role in Arrhythmogenesis

T-wave alternans (TWA) reflects every-other-beat alterations in the morphology of the electrocardiogram ST segment or T wave in the setting of a constant heart rate, hence, in the absence of heart rate variability. It is believed to be associated with the dispersion of repolarization and has been used as a non-invasive marker for predicting the risk of malignant cardiac arrhythmias and sudden cardiac death as numerous studies have shown. This review aims to provide up-to-date review on both experimental and simulation studies in elucidating possible mechanisms underlying the genesis of TWA at the cellular level, as well as the genesis of spatially concordant/discordant alternans at the tissue level, and their transition to cardiac arrhythmia. Recent progress and future perspectives in antiarrhythmic therapies associated with TWA are also discussed.

Mechanistic link between CaM-RyR2 interactions and the genesis of cardiac arrhythmia

In this study, we develop a computational model of the interaction between ryanodine receptor type 2 (RyR2) and calmodulin (CaM) to explore the mechanistic link between CaM-RyR2 interactions and cardiac arrhythmia. Our starting point is a biophysically based computational model of CaM binding to a single RyR2 subunit, which reproduces single-channel RyR2 measurements in lipid bilayers. We then integrate this CaM-RyR2 model into a spatially distributed whole-cell model of Ca cycling, which is used to investigate the relationship between CaM and Ca cycling homeostasis. We show that a reduction in CaM concentration leads to a substantial increase in the rate of spontaneous Ca sparks, and this induces a marked reduction in sarcoplasmic reticulum Ca load during steady-state pacing. Also, we show that a reduction in CaM modifies the RyR2 open probability, which makes the cell more prone to Ca wave propagation. These results indicate that aberrant Ca cycling activity during pacing is determined by the interplay between sarcoplasmic reticulum load reduction and the threshold for Ca wave propagation. Based on these results, we show that when CaM is reduced, Ca waves can occur in a cell and induce action potential perturbations that are arrhythmogenic. Thus, this study outlines a novel, to our knowledge, mechanistic link between CaM-RyR2 binding kinetics and the induction of arrhythmias in the heart.

Mitochondrial depolarization promotes calcium alternans: Mechanistic insights from a ventricular myocyte model

Mitochondria are vital organelles inside the cell and contribute to intracellular calcium (Ca2+) dynamics directly and indirectly via calcium exchange, ATP generation, and production of reactive oxygen species (ROS). Arrhythmogenic Ca2+ alternans in cardiac myocytes has been observed in experiments under abnormal mitochondrial depolarization. However, complex signaling pathways and Ca2+ cycling between mitochondria and cytosol make it difficult in experiments to reveal the underlying mechanisms of Ca2+ alternans under abnormal mitochondrial depolarization. In this study, we use a newly developed spatiotemporal ventricular myocyte computer model that integrates mitochondrial Ca2+ cycling and complex signaling pathways to investigate the mechanisms of Ca2+ alternans during mitochondrial depolarization. We find that elevation of ROS in response to mitochondrial depolarization plays a critical role in promoting Ca2+ alternans. Further examination reveals that the redox effect of ROS on ryanodine receptors and sarco/endoplasmic reticulum Ca2+-ATPase synergistically promote alternans. Upregulation of mitochondrial Ca2+ uniporter promotes Ca2+ alternans via Ca2+-dependent mitochondrial permeability transition pore opening. Due to their relatively slow kinetics, oxidized Ca2+/calmodulin-dependent protein kinase II activation and ATP do not play significant roles acutely in the genesis of Ca2+ alternans after mitochondrial depolarization, but their roles can be significant in the long term, mainly through their effects on sarco/endoplasmic reticulum Ca2+-ATPase activity. In conclusion, mitochondrial depolarization promotes Ca2+ alternans acutely via the redox effect of ROS and chronically by ATP reduction. It suppresses Ca2+ alternans chronically through oxidized Ca2+/calmodulin-dependent protein kinase II activation.

A Novel Computational Model of the Rabbit Atrial Cardiomyocyte With Spatial Calcium Dynamics

Models of cardiac electrophysiology are widely used to supplement experimental results and to provide insight into mechanisms of cardiac function and pathology. The rabbit has been a particularly important animal model for studying mechanisms of atrial pathophysiology and atrial fibrillation, which has motivated the development of models for the rabbit atrial cardiomyocyte electrophysiology. Previously developed models include detailed representations of membrane currents and intracellular ionic concentrations, but these so-called “common-pool” models lack a spatially distributed description of the calcium handling system, which reflects the detailed ultrastructure likely found in cells in vivo. Because of the less well-developed T-tubular system in atrial compared to ventricular cardiomyocytes, spatial gradients in intracellular calcium concentrations may play a more significant role in atrial cardiomyocyte pathophysiology, rendering common-pool models less suitable for investigating underlying electrophysiological mechanisms. In this study, we developed a novel computational model of the rabbit atrial cardiomyocyte incorporating detailed compartmentalization of intracellular calcium dynamics, in addition to a description of membrane currents and intracellular processes. The spatial representation of calcium was based on dividing the intracellular space into eighteen different compartments in the transversal direction, each with separate systems for internal calcium storage and release, and tracking ionic fluxes between compartments in addition to the dynamics driven by membrane currents and calcium release. The model was parameterized employing a population-of-models approach using experimental data from different sources. The parameterization of this novel model resulted in a reduced population of models with inherent variability in calcium dynamics and electrophysiological properties, all of which fall within the range of observed experimental values. As such, the population of models may represent natural variability in cardiomyocyte electrophysiology or inherent uncertainty in the underlying experimental data. The ionic model population was also able to reproduce the U-shaped waveform observed in line-scans of triggered calcium waves in atrial cardiomyocytes, characteristic of the absence of T-tubules, resulting in a centripetal calcium wave due to subcellular calcium diffusion. This novel spatial model of the rabbit atrial cardiomyocyte can be used to integrate experimental findings, offering the potential to enhance our understanding of the pathophysiological role of calcium-handling abnormalities under diseased conditions, such as atrial fibrillation.

Delayed global feedback in the genesis and stability of spatiotemporal excitation patterns in paced biological excitable media

Biological excitable media, such as cardiac or neural cells and tissue, exhibit memory in which a change in the present excitation may affect the behaviors in the next excitation. For example, a change in calcium (Ca²⁺) concentration in a cell in the present excitation may affect the Ca²⁺ dynamics in the next excitation via bi-directional coupling between voltage and Ca²⁺, forming a delayed feedback loop. Since the Ca²⁺ dynamics inside the excitable cells are spatiotemporal while the membrane voltage is a global signal, the feedback loop is then a delayed global feedback (DGF) loop. In this study, we investigate the roles of DGF in the genesis and stability of spatiotemporal excitation patterns in periodically-paced excitable media using mathematical models with different levels of complexity: a model composed of coupled FitzHugh-Nagumo units, a 3-dimensional physiologically-detailed ventricular myocyte model, and a coupled map lattice model. We investigate the dynamics of excitation patterns that are temporal period-2 (P2) and spatially concordant or discordant, such as subcellular concordant or discordant Ca²⁺alternans in cardiac myocytes or spatially concordant or discordant Ca²⁺ and repolarization alternans in cardiac tissue. Our modeling approach allows both computer simulations and rigorous analytical treatments, which lead to the following results and conclusions. When DGF is absent, concordant and discordant P2 patterns occur depending on initial conditions with the discordant P2 patterns being spatially random. When the DGF is negative, only concordant P2 patterns exist. When the DGF is positive, both concordant and discordant P2 patterns can occur. The discordant P2 patterns are still spatially random, but they satisfy that the global signal exhibits a temporal period-1 behavior. The theoretical analyses of the coupled map lattice model reveal the underlying instabilities and bifurcations for the genesis, selection, and stability of spatiotemporal excitation patterns.

Understanding the mechanisms of pattern formation in biological systems is of great importance. Here we investigate the dynamical mechanisms by which delayed global feedback affects excitation pattern formation and stability in periodically-paced biological excitable media, such as cardiac or neural cells and tissue. We focus on the formation and stability of the temporal period-2 and spatially in-phase and out-of-phase excitation patterns. Using models of different levels of complexity, we show that when the delayed global feedback is negative, only the spatially in-phase patterns are stable. When the feedback is positive, both spatially in-phase and out-of-phase patterns are stable, and the out-of-phase patterns are spatially random but satisfy that the global signals are temporal period-1 solutions.

Dual regulation by subcellular calcium heterogeneity and heart rate variability on cardiac electromechanical dynamics

Heart rate constantly varies under physiological conditions, termed heart rate variability (HRV), and in clinical studies, low HRV is associated with a greater risk of cardiac arrhythmias. Prior work has shown that HRV influences the temporal patterns of electrical activity, specifically the formation of pro-arrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca) levels. We previously showed that HRV may be anti-arrhythmic by disrupting APD and Ca alternations in a homogeneous cardiac myocyte. Here, we expand on our previous work, incorporating variation in subcellular Ca handling (also known to influence alternans) into a nonlinear map model of a cardiac myocyte composed of diffusively coupled Ca release units (CRUs). Ca-related parameters and initial conditions of each CRU are varied to mimic subcellular Ca heterogeneity, and a stochastic pacing sequence reproduces HRV. We find that subcellular Ca heterogeneity promotes the formation of spatially discordant subcellular alternans patterns, which decreases whole cell Ca and APD alternation for low and moderate HRV, while high subcellular Ca heterogeneity and HRV both promote electromechanical desynchronization. Finally, we find that for low and moderate HRV, both the specific subcellular Ca-related parameters and the pacing sequences influence measures of electromechanical dynamics, while for high HRV, these measures depend predominantly on the pacing sequence. Our results suggest that pro-arrhythmic subcellular discordant alternans tend to form for low levels of HRV, while high HRV may be anti-arrhythmic due to mitigated influence from subcellular Ca heterogeneity and desynchronization of APD from Ca instabilities.

Influence of the tubular network on the characteristics of calcium transients in cardiac myocytes

Transverse and axial tubules (TATS) are an essential ingredient of the excitation-contraction machinery that allow the effective coupling of L-type Calcium Channels (LCC) and ryanodine receptors (RyR2). They form a regular network in ventricular cells, while their presence in atrial myocytes is variable regionally and among animal species We have studied the effect of variations in the TAT network using a bidomain computational model of an atrial myocyte with variable density of tubules. At each z-line the t-tubule length is obtained from an exponential distribution, with a given mean penetration length. This gives rise to a distribution of t-tubules in the cell that is characterized by the fractional area (F.A.) occupied by the t-tubules. To obtain consistent results, we average over different realizations of the same mean penetration length. To this, in some simulations we add the effect of a network of axial tubules. Then we study global properties of calcium signaling, as well as regional heterogeneities and local properties of sparks and RyR2 openings. In agreement with recent experiments in detubulated ventricular and atrial cells, we find that detubulation reduces the calcium transient and synchronization in release. However, it does not affect sarcoplasmic reticulum (SR) load, so the decrease in SR calcium release is due to regional differences in Ca²⁺ release, that is restricted to the cell periphery in detubulated cells. Despite the decrease in release, the release gain is larger in detubulated cells, due to recruitment of orphaned RyR2s, i.e, those that are not confronting a cluster of LCCs. This probably provides a safeguard mechanism, allowing physiological values to be maintained upon small changes in the t-tubule density. Finally, we do not find any relevant change in spark properties between tubulated and detubulated cells, suggesting that the differences found in experiments could be due to differential properties of the RyR2s in the membrane and in the t-tubules, not incorporated in the present model. This work will help understand the effect of detubulation, that has been shown to occur in disease conditions such as heart failure (HF) in ventricular cells, or atrial fibrillation (AF) in atrial cells.

Multi-scale approaches for the simulation of cardiac electrophysiology: I - Sub-cellular and stochastic calcium dynamics from cell to organ

Computational models of the heart at multiple spatial scales, from sub-cellular nanodomains to the whole-organ, are a powerful tool for the simulation of cardiac electrophysiology. Application of these models has provided remarkable insight into the normal and pathological functioning of the heart. In these two articles, we present methods for modelling cardiac electrophysiology at all of these spatial scales. In part one, presented here, we discuss methods and approaches for modelling sub-cellular calcium dynamics at the whole-cell and organ scales, valuable for modelling excitation-contraction coupling and mechanisms of arrhythmia triggers.

Remodeling Promotes Proarrhythmic Disruption of Calcium Homeostasis in Failing Atrial Myocytes

It is well known that heart failure (HF) typically coexists with atrial fibrillation (AF). However, until now, no clear mechanism has been established that relates HF to AF. In this study, we apply a multiscale computational framework to establish a mechanistic link between atrial myocyte structural remodeling in HF and AF. Using a spatially distributed model of calcium (Ca) signaling, we show that disruption of the spatial relationship between L-type Ca channels (LCCs) and ryanodine receptors results in markedly increased Ca content of the sarcoplasmic reticulum (SR). This increase in SR load is due to changes in the balance between Ca entry via LCCs and Ca extrusion due to the sodium-calcium exchanger after an altered spatial relationship between these signaling proteins. Next, we show that the increased SR load in atrial myocytes predisposes these cells to subcellular Ca waves that occur during the action potential (AP) and are triggered by LCC openings. These waves are common in atrial cells because of the absence of a well-developed t-tubule system in most of these cells. This distinct spatial architecture allows for the presence of a large pool of orphaned ryanodine receptors, which can fire and sustain Ca waves during the AP. Finally, we incorporate our atrial cell model in two-dimensional tissue simulations and demonstrate that triggered wave generation in cells leads to electrical waves in tissue that tend to fractionate to form wavelets of excitation. This fractionation is driven by the underlying stochasticity of subcellular Ca waves, which perturbs AP repolarization and consequently induces localized conduction block in tissue. We outline the mechanism for this effect and argue that it may explain the propensity for atrial arrhythmias in HF.

Arrhythmia mechanisms and spontaneous calcium release: Bi-directional coupling between re-entrant and focal excitation

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation.

The Subcellular Distribution of Ryanodine Receptors and L-Type Ca2+ Channels Modulates Ca2+-Transient Properties and Spontaneous Ca2+-Release Events in Atrial Cardiomyocytes

Spontaneous Ca²⁺-release events (SCaEs) from the sarcoplasmic reticulum play crucial roles in the initiation of cardiac arrhythmias by promoting triggered activity. However, the subcellular determinants of these SCaEs remain incompletely understood. Structural differences between atrial and ventricular cardiomyocytes, e.g., regarding the density of T-tubular membrane invaginations, may influence cardiomyocyte Ca²⁺-handling and the distribution of cardiac ryanodine receptors (RyR2) has recently been shown to undergo remodeling in atrial fibrillation. These data suggest that the subcellular distribution of Ca²⁺-handling proteins influences proarrhythmic Ca²⁺-handling abnormalities. Here, we employ computational modeling to provide an in-depth analysis of the impact of variations in subcellular RyR2 and L-type Ca²⁺-channel distributions on Ca²⁺-transient properties and SCaEs in a human atrial cardiomyocyte model. We incorporate experimentally observed RyR2 expression patterns and various configurations of axial tubules in a previously published model of the human atrial cardiomyocyte. We identify an increased SCaE incidence for larger heterogeneity in RyR2 expression, in which SCaEs preferentially arise from regions of high local RyR2 expression. Furthermore, we show that the propagation of Ca²⁺ waves is modulated by the distance between RyR2 bands, as well as the presence of experimentally observed RyR2 clusters between bands near the lateral membranes. We also show that incorporation of axial tubules in various amounts and locations reduces Ca²⁺-transient time to peak. Furthermore, selective hyperphosphorylation of RyR2 around axial tubules increases the number of spontaneous waves. Finally, we present a novel model of the human atrial cardiomyocyte with physiological RyR2 and L-type Ca²⁺-channel distributions that reproduces experimentally observed Ca²⁺-handling properties. Taken together, these results significantly enhance our understanding of the structure-function relationship in cardiomyocytes, identifying that RyR2 and L-type Ca²⁺-channel distributions have a major impact on systolic Ca²⁺ transients and SCaEs.