Calcium alternans is due to an order-disorder phase transition in cardiac cells

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'.

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.

An ensemble of parameters from a robust Markov-based model reproduces L-type calcium currents from different human cardiac myocytes

The development of modeling structures at the channel level that can integrate subcellular and cell models and properly reproduce different experimental data is of utmost importance in cardiac electrophysiology. In contrast to gate-based models, Markov Chain models are well suited to promote the integration of the subcellular level of the cardiomyocyte to the whole cell. In this paper, we develop Markov Chain models for the L-type Calcium current that can reproduce the electrophysiology of two established human models for the ventricular and Purkinje cells. In addition, instead of presenting a single set of parameters, we present a collection of set of parameters employing Differential Evolution algorithms that can properly reproduce very different protocol data. We show the importance of using an ensemble of a set of parameter values to obtain proper results when considering a second protocol that suppresses calcium inactivation and mimics a pathological condition. We discuss how model discrepancy, data availability, and parameter identifiability can influence the choice of the size of the collection. In summary, we have modified two cardiac models by proposing new Markov Chain models for the L-type Calcium. We keep the original whole-cell dynamics by reproducing the same characteristic action potential and calcium dynamics, whereas the Markov chain-based description of the L-type Calcium channels allows novel small spatial scale simulations of subcellular processes. Finally, the use of collections of parameters was crucial for addressing model discrepancy, identifiability issues, and avoiding fitting parameters overly precisely, i.e., overfitting.

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.

Voltage-mediated mechanism for calcium wave synchronization and arrhythmogenesis in atrial tissue

A wide range of atrial arrythmias are caused by molecular defects in proteins that regulate calcium (Ca) cycling. In many cases, these defects promote the propagation of subcellular Ca waves in the cell, which can perturb the voltage time course and induce dangerous perturbations of the action potential (AP). However, subcellular Ca waves occur randomly in cells and, therefore, electrical coupling between cells substantially decreases their effect on the AP. In this study, we present evidence that Ca waves in atrial tissue can synchronize in-phase owing to an order-disorder phase transition. In particular, we show that, below a critical pacing rate, Ca waves are desynchronized and therefore do not induce substantial AP fluctuations in tissue. However, above this critical pacing rate, Ca waves gradually synchronize over millions of cells, which leads to a dramatic amplification of AP fluctuations. We exploit an underlying Ising symmetry of paced cardiac tissue to show that this transition exhibits universal properties common to a wide range of physical systems in nature. Finally, we show that in the heart, phase synchronization induces spatially out-of-phase AP duration alternans which drives wave break and reentry. These results suggest that cardiac tissue exhibits a phase transition that is required for subcellular Ca cycling defects to induce a life-threatening arrhythmia.

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.

Stability of spatially discordant repolarization alternans in cardiac tissue

Cardiac alternans, a period-2 behavior of excitation and contraction of the heart, is a precursor of ventricular arrhythmias and sudden cardiac death. One form of alternans is repolarization or action potential duration alternans. In cardiac tissue, repolarization alternans can be spatially in-phase, called spatially concordant alternans, or spatially out-of-phase, called spatially discordant alternans (SDA). In SDA, the border between two out-of-phase regions is called a node in a one-dimensional cable or a nodal line in a two-dimensional tissue. In this study, we investigate the stability and dynamics of the nodes and nodal lines of repolarization alternans driven by voltage instabilities. We use amplitude equation and coupled map lattice models to derive theoretical results, which are compared with simulation results from the ionic model. Both conduction velocity restitution induced SDA and non-conduction velocity restitution induced SDA are investigated. We show that the stability and dynamics of the SDA nodes or nodal lines are determined by the balance of the tensions generated by conduction velocity restitution, convection due to action potential propagation, curvature of the nodal lines, and repolarization and coupling heterogeneities. Our study provides mechanistic insights into the different SDA behaviors observed in experiments.

Buffering and total calcium levels determine the presence of oscillatory regimes in cardiac cells

Calcium oscillations and waves induce depolarization in cardiac cells which are believed to cause life-threathening arrhythimas. In this work, we study the conditions for the appearance of calcium oscillations in both a detailed subcellular model of calcium dynamics and a minimal model that takes into account just the minimal ingredients of the calcium toolkit. To avoid the effects of homeostatic changes and the interaction with the action potential we consider the somewhat artificial condition of a cell without pacing and with no calcium exchange with the extracellular medium. Both the full subcellular model and the minimal model present the same scenarios depending on the calcium load: two stationary states, one with closed ryanodine receptors (RyR) and most calcium in the cell stored in the sarcoplasmic reticulum (SR), and another, with open RyRs and a depleted SR. In between, calcium oscillations may appear. The robustness of these oscillations is determined by the amount of calsequestrin (CSQ). The lack of this buffer in the SR enhances the appearance of oscillations. The minimal model allows us to relate the stability of the oscillating state to the nullcline structure of the system, and find that its range of existence is bounded by a homoclinic and a Hopf bifurcation, resulting in a sudden transition to the oscillatory regime as the cell calcium load is increased. Adding a small amount of noise to the RyR behavior increases the parameter region where oscillations appear and provides a gradual transition from the resting state to the oscillatory regime, as observed in the subcellular model and experimentally.

In cardiac cells, calcium plays a very important role. An increase in calcium levels is the trigger used by the cell to initiate contraction. Besides, calcium modulates several transmembrane currents, affecting the cell transmembrane potential. Thus, dysregulations in calcium handling have been associated with the appearance of arrhythmias. Often, this dysregulation results in the appearance of periodic calcium waves or global oscillations, providing a pro-arrhythmic substrate. In this paper, we study the onset of calcium oscillations in cardiac cells using both a detailed subcellular model of calcium dynamics and a minimal model that takes into account the essential ingredients of the calcium toolkit. Both reproduce the main experimental results and link this behavior with the presence of different steady-state solutions and bifurcations that depend on the total amount of calcium in the cell and in the level of buffering present. We expect that this work will help to clarify the conditions under which calcium oscillations appear in cardiac myocytes and, therefore, will represent a step further in the understanding of the origin of cardiac arrhythmias.

Two-variable nullcline analysis of ionic general equilibrium predicts calcium homeostasis in ventricular myocytes

Ventricular contraction is roughly proportional to the amount of calcium released from the Sarcoplasmic Reticulum (SR) during systole. While it is rather straightforward to measure calcium levels and contractibility under different physiological conditions, the complexity of calcium handling during systole and diastole has made the prediction of its release at steady state impossible. Here we approach the problem analyzing the evolution of intracellular and extracellular calcium fluxes during a single beat which is away from homeostatic balance. Using an in-silico subcellular model of rabbit ventricular myocyte, we show that the high dimensional nonlinear problem of finding the steady state can be reduced to a two-variable general equilibrium condition where pre-systolic calcium level in the cytosol and in the SR must fulfill simultaneously two different equalities. This renders calcium homeostasis as a problem that can be studied in terms of its equilibrium structure, leading to precise predictions of steady state from single-beat measurements. We show how changes in ion channels modify the general equilibrium, as shocks would do in general equilibrium macroeconomic models. This allows us to predict when an enhanced entrance of calcium in the cell reduces its contractibility and explain why SERCA gene therapy, a change in calcium handling to treat heart failure, might fail to improve contraction even when it successfully increases SERCA expression.

Cardiomyocytes, upon voltage excitation, release calcium, which leads to cell contraction. However, under some pathological conditions, calcium handling is impaired. Recently, SERCA gene therapy, whose aim is to improve Ca²⁺ sequestration by the Sarcoplasmic Reticulum (SR), has failed to improve the prognosis of patients with Heart Failure. This, together with recent counterintuitive results in calcium handling, has highlighted the need for a framework to understand calcium homeostasis across species and pathologies. We show here that the proper framework is a general equilibrium approach of two independent variables. The development of this framework allows us to find a possible mechanism for the failure of SERCA gene therapy even when it manages to increase Ca SERCA expression.

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.

Exploring Impaired SERCA Pump-Caused Alternation Occurrence in Ischemia

Impaired sarcoplasmic reticulum (SR) calcium transport ATPase (SERCA) gives rise to Ca²⁺ alternans and changes of the Ca2+release amount. These changes in Ca²⁺ release amount can reveal the mechanism underlying how the interaction between Ca²⁺ release and Ca²⁺ uptake induces Ca²⁺ alternans. This study of alternans by calculating the values of Ca²⁺ release properties with impaired SERCA has not been explored before. Here, we induced Ca²⁺ alternans by using an impaired SERCA pump under ischemic conditions. The results showed that the recruitment and refractoriness of the Ca²⁺ release increased as Ca²⁺ alternans occurred. This indicates triggering Ca waves. As the propagation of Ca waves is linked to the occurrence of Ca²⁺ alternans, the “threshold” for Ca waves reflects the key factor in Ca²⁺ alternans development, and it is still controversial nowadays. We proposed the ratio between the diastolic network SR (NSR) Ca content (Ca_nsr) and the cytoplasmic Ca content (Ca_i) (Ca_nsr/Ca_i) as the “threshold” of Ca waves and Ca²⁺ alternans. Diastolic Ca_nsr, Ca_i, and their ratio were recorded at the onset of Ca²⁺ alternans. Compared with certain Ca_nsr and Ca_i, the “threshold” of the ratio can better explain the comprehensive effects of the Ca²⁺ release and the Ca²⁺ uptake on Ca²⁺ alternans onset. In addition, these ratios are related with the function of SERCA pumps, which vary with different ischemic conditions. Thus, values of these ratios could be used to differentiate Ca²⁺ alternans from different ischemic cases. This agrees with some experimental results. Therefore, the certain value of diastolic Ca_nsr/Ca_i can be the better “threshold” for Ca waves and Ca²⁺ alternans.

Complex patterns of subcellular cardiac alternans

Cardiac alternans, in which the membrane potential and the intracellular calcium concentration exhibit alternating durations and peak amplitudes at consecutive beats, constitute a precursor to fatal cardiac arrhythmia such as sudden cardiac death. A crucial question therefore concerns the onset of cardiac alternans. Typically, alternans are only reported when they are fully developed. Here, we present a modelling approach to explore recently discovered microscopic alternans, which represent one of the earliest manifestations of cardiac alternans. In this case, the regular periodic dynamics of the local intracellular calcium concentration is already unstable, while the whole-cell behaviour suggests a healthy cell state. In particular, we use our model to investigate the impact of calcium diffusion in both the cytosol and the sarcoplasmic reticulum on the formation of microscopic calcium alternans. We find that for dominant cytosolic coupling, calcium alternans emerge via the traditional period doubling bifurcation. In contrast, dominant luminal coupling leads to a novel route to calcium alternans through a saddle-node bifurcation at the network level. Combining semi-analytical and computational approaches, we compute areas of stability in parameter space and find that as we cross from stable to unstable regions, the emergent patterns of the intracellular calcium concentration change abruptly in a fashion that is highly dependent upon position along the stability boundary. Our results demonstrate that microscopic calcium alternans may possess a much richer dynamical repertoire than previously thought and further strengthen the role of luminal calcium in shaping cardiac calcium dynamics.

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.

Stochastic coupled map model of subcellular calcium cycling in cardiac cells

In this study, we analyze a nonlinear map model of intracellular calcium (Ca) and voltage in cardiac cells. In this model, Ca release from the sarcoplasmic reticulum (SR) occurs at spatially distributed dyadic junctions that are diffusively coupled. At these junctions, release occurs with a probability that depends on key variables such as the SR load and the diastolic interval. Using this model, we explore how nonlinearity and stochasticity determine the spatial distribution of Ca release events within a cardiac cell. In particular, we identify a novel synchronization transition, which occurs at rapid pacing rates, in which the global Ca transient transitions from a period 2 response to a period 1 response. In the global period 2 response dyadic junctions fire in unison, on average, on alternate beats, while in the period 1 regime, Ca release at individual dyads is highly irregular. A close examination of the spatial distribution of Ca reveals that in the period 1 regime, the system coarsens into spatially out-of-phase regions with a length scale much smaller than the system size, but larger than the spacing between dyads. We have also explored in detail the coupling to membrane voltage. We study first the case of positive coupling, where a large Ca transient promotes a long action potential duration (APD). Here, the coupling to voltage synchronizes Ca release so that the system exhibits a robust period 2 response that is independent of initial conditions. On the other hand, in the case of negative coupling, where a large Ca transient tends to shorten the APD, we find a multitude of metastable states which consist of complex spatially discordant alternans patterns. Using an analogy to equilibrium statistical mechanics, we show that the spatial patterns observed can be explained by a mapping to the Potts model, with an additional term that accounts for a global coupling of spin states. Using this analogy, we argue that Ca cycling in cardiac cells exhibits complex spatiotemporal patterns that emerge via first or second order phase transitions. These results show that voltage and Ca can interact in order to induce complex subcellular responses, which can potentially lead to heart rhythm disorders.

Computational Model of Calcium Signaling in Cardiac Atrial Cells at the Submicron Scale

In cardiac cells, calcium is the mediator of excitation-contraction coupling. Dysfunctions in calcium handling have been identified as the origin of some cardiac arrhythmias. In the particular case of atrial myocytes, recent available experimental data has found links between these dysfunctions and structural changes in the calcium handling machinery (ryanodine cluster size and distribution, t-tubular network, etc). To address this issue, we have developed a computational model of an atrial myocyte that takes into account the detailed intracellular structure. The homogenized macroscopic behavior is described with a two-concentration field model, using effective diffusion coefficients of calcium in the sarcoplasmic reticulum (SR) and in the cytoplasm. The model reproduces the right calcium transients and dependence with pacing frequency. Under basal conditions, the calcium rise is mostly restricted to the periphery of the cell, with a large concentration ratio between the periphery and the interior. We have then studied the dependence of the speed of the calcium wave on cytosolic and SR diffusion coefficients, finding an almost linear relation with the former, in agreement with a diffusive and fire mechanism of propagation, and little dependence on the latter. Finally, we have studied the effect of a change in RyR cluster microstructure. We find that, under resting conditions, the spark frequency decreases slightly with RyR cluster spatial dispersion, but markedly increases when the RyRs are distributed in clusters of larger size, stressing the importance of RyR cluster organization to understand atrial arrhythmias, as recent experimental results suggest (Macquaide et al., 2015).

Locating Order-Disorder Phase Transition in a Cardiac System

To prevent sudden cardiac death, predicting where in the cardiac system an order-disorder phase transition into ventricular fibrillation begins is as important as when it begins. We present a computationally efficient, information-theoretic approach to predicting the locations of the wavebreaks. Such wavebreaks initiate fibrillation in a cardiac system where the order-disorder behavior is controlled by a single driving component, mimicking electrical misfiring from the pulmonary veins or from the Purkinje fibers. Communication analysis between the driving component and each component of the system reveals that channel capacity, mutual information and transfer entropy can locate the wavebreaks. This approach is applicable to interventional therapies to prevent sudden death, and to a wide range of systems to mitigate or prevent imminent phase transitions.

Minimal model for calcium alternans due to SR release refractoriness

In the heart, rapid pacing rates may induce alternations in the strength of cardiac contraction, termed pulsus alternans. Often, this is due to an instability in the dynamics of the intracellular calcium concentration, whose transients become larger and smaller at consecutive beats. This alternation has been linked experimentally and theoretically to two different mechanisms: an instability due to (1) a strong dependence of calcium release on sarcoplasmic reticulum (SR) load, together with a slow calcium reuptake into the SR or (2) to SR release refractoriness, due to a slow recovery of the ryanodine receptors (RyR2) from inactivation. The relationship between calcium alternans and refractoriness of the RyR2 has been more elusive than the corresponding SR Ca load mechanism. To study the former, we reduce a general calcium model, which mimics the deterministic evolution of a calcium release unit, to its most basic elements. We show that calcium alternans can be understood using a simple nonlinear equation for calcium concentration at the dyadic space, coupled to a relaxation equation for the number of recovered RyR2s. Depending on the number of RyR2s that are recovered at the beginning of a stimulation, the increase in calcium concentration may pass, or not, over an excitability threshold that limits the occurrence of a large calcium transient. When the recovery of the RyR2 is slow, this produces naturally a period doubling bifurcation, resulting in calcium alternans. We then study the effects of inactivation, calcium diffusion, and release conductance for the onset of alternans. We find that the development of alternans requires a well-defined value of diffusion while it is less sensitive to the values of inactivation or release conductance.

A computational model of spatio-temporal cardiac intracellular calcium handling with realistic structure and spatial flux distribution from sarcoplasmic reticulum and t-tubule reconstructions

Intracellular calcium cycling is a vital component of cardiac excitation-contraction coupling. The key structures responsible for controlling calcium dynamics are the cell membrane (comprising the surface sarcolemma and transverse-tubules), the intracellular calcium store (the sarcoplasmic reticulum), and the co-localisation of these two structures to form dyads within which calcium-induced-calcium-release occurs. The organisation of these structures tightly controls intracellular calcium dynamics. In this study, we present a computational model of intracellular calcium cycling in three-dimensions (3-D), which incorporates high resolution reconstructions of these key regulatory structures, attained through imaging of tissue taken from the sheep left ventricle using serial block face scanning electron microscopy. An approach was developed to model the sarcoplasmic reticulum structure at the whole-cell scale, by reducing its full 3-D structure to a 3-D network of one-dimensional strands. The model reproduces intracellular calcium dynamics during control pacing and reveals the high-resolution 3-D spatial structure of calcium gradients and intracellular fluxes in both the cytoplasm and sarcoplasmic reticulum. We also demonstrated the capability of the model to reproduce potentially pro-arrhythmic dynamics under perturbed conditions, pertaining to calcium-transient alternans and spontaneous release events. Comparison with idealised cell models emphasised the importance of structure in determining calcium gradients and controlling the spatial dynamics associated with calcium-transient alternans, wherein the probabilistic nature of dyad activation and recruitment was constrained. The model was further used to highlight the criticality in calcium spark propagation in relation to inter-dyad distances. The model presented provides a powerful tool for future investigation of structure-function relationships underlying physiological and pathophysiological intracellular calcium handling phenomena at the whole-cell. The approach allows for the first time direct integration of high-resolution images of 3-D intracellular structures with models of calcium cycling, presenting the possibility to directly assess the functional impact of structural remodelling at the cellular scale.

The organisation of the membrane and sub-cellular structures of cells in the heart closely controls the coupling between its electrical and mechanical function. Computational models of the cellular calcium handling system, which is responsible for this electro-mechanical coupling, have been developed in recent years to study underlying structure-function relationships. Previous models have been largely idealised in structure; we present a new model which incorporates experimental data describing the high-resolution organisation of the primary structures involved in calcium dynamics. Significantly, the structure of the intracellular calcium store is modelled for the first time. The model is shown to reproduce calcium dynamics in control cells in both normal and abnormal conditions, demonstrating its suitability for future investigation of structure-function relationships. Thus, the model presented provides a powerful tool for the direct integration of experimentally acquired structural data in healthy and diseased cells and assessment of the role of structure in regulating normal and abnormal calcium dynamics.

A unified theory of calcium alternans in ventricular myocytes

Intracellular calcium (Ca²⁺) alternans is a dynamical phenomenon in ventricular myocytes, which is linked to the genesis of lethal arrhythmias. Iterated map models of intracellular Ca²⁺ cycling dynamics in ventricular myocytes under periodic pacing have been developed to study the mechanisms of Ca²⁺ alternans. Two mechanisms of Ca²⁺ alternans have been demonstrated in these models: one relies mainly on fractional sarcoplasmic reticulum Ca²⁺ release and uptake, and the other on refractoriness and other properties of Ca²⁺ sparks. Each of the two mechanisms can partially explain the experimental observations, but both have their inconsistencies with the experimental results. Here we developed an iterated map model that is composed of two coupled iterated maps, which unifies the two mechanisms into a single cohesive mathematical framework. The unified theory can consistently explain the seemingly contradictory experimental observations and shows that the two mechanisms work synergistically to promote Ca²⁺ alternans. Predictions of the theory were examined in a physiologically-detailed spatial Ca²⁺ cycling model of ventricular myocytes.

Nonlinear onset of calcium wave propagation in cardiac cells

Spontaneous calcium (Ca) waves in cardiac myocytes are known to underlie a wide range of cardiac arrhythmias. However, it is not understood which physiological parameters determine the onset of waves. In this study, we explore the relationship between Ca signaling between ion channels and the nucleation of Ca waves. In particular, we apply a master equation approach to analyze the stochastic interaction between neighboring clusters of ryanodine receptor (RyR) channels. Using this analysis, we show that signaling between clusters can be described as a barrier hopping process with exponential sensitivity to system parameters. A consequence of this feature is that the probability that Ca release at a cluster induces release at a neighboring cluster exhibits a sigmoid dependence on the Ca content in the cell. This nonlinearity originates from the regulation of RyR opening due to more than one Ca ion binding site, in conjunction with Ca mediated cooperativity between RyR channels in clusters. We apply a spatially distributed stochastic model of Ca cycling to analyze the physiological consequences of this nonlinearity, and show that it explains the sharp onset of Ca wave nucleation in cardiac cells. Furthermore, we show that this sharp onset can serve as a mechanism for Ca alternans under physiologically relevant conditions. Thus our findings identify the nonlinear features of Ca signaling which potentially underlie the onset of Ca waves and Ca alternans in cardiac cells.

Long-Lasting Sparks: Multi-Metastability and Release Competition in the Calcium Release Unit Network

Calcium (Ca) sparks are elementary events of biological Ca signaling. A normal Ca spark has a brief duration in the range of 10 to 100 ms, but long-lasting sparks with durations of several hundred milliseconds to seconds are also widely observed. Experiments have shown that the transition from normal to long-lasting sparks can occur when ryanodine receptor (RyR) open probability is either increased or decreased. Here, we demonstrate theoretically and computationally that long-lasting sparks emerge as a collective dynamical behavior of the network of diffusively coupled Ca release units (CRUs). We show that normal sparks occur when the CRU network is monostable and excitable, while long-lasting sparks occur when the network dynamics possesses multiple metastable attractors, each attractor corresponding to a different spatial firing pattern of sparks. We further highlight the mechanisms and conditions that produce long-lasting sparks, demonstrating the existence of an optimal range of RyR open probability favoring long-lasting sparks. We find that when CRU firings are sparse and sarcoplasmic reticulum (SR) Ca load is high, increasing RyR open probability promotes long-lasting sparks by potentiating Ca-induced Ca release (CICR). In contrast, when CICR is already strong enough to produce frequent firings, decreasing RyR open probability counter-intuitively promotes long-lasting sparks by decreasing spark frequency. The decrease in spark frequency promotes intra-SR Ca diffusion from neighboring non-firing CRUs to the firing CRUs, which helps to maintain the local SR Ca concentration of the firing CRUs above a critical level to sustain firing. In this setting, decreasing RyR open probability further suppresses long-lasting sparks by weakening CICR. Since a long-lasting spark terminates via the Kramers’ escape process over a potential barrier, its duration exhibits an exponential distribution determined by the barrier height and noise strength, which is modulated differently by different ways of altering the Ca release flux strength.

Calcium (Ca) sparks, resulting from Ca-induced Ca release, are elementary events of biological Ca signaling. Sparks are normally brief, but long-lasting sparks have been widely observed experimentally under various conditions. The underlying mechanisms of spark duration or termination and the corresponding determinants remain a topic of debate. In this study, we demonstrate theoretically and computationally that normal brief sparks are excitable transients, while long-lasting sparks are multiple metastable states emerging in the diffusively coupled Ca release unit network, as a result of cooperativity and release competition among the Ca release units. Termination of a long-lasting spark is a Kramers’ escape process over a potential barrier, and the spark duration is the first-passage time, exhibiting an exponential distribution.

Predicting the onset of period-doubling bifurcations in noisy cardiac systems

Biological, physical, and social systems often display qualitative changes in dynamics. Developing early warning signals to predict the onset of these transitions is an important goal. The current work is motivated by transitions of cardiac rhythms, where the appearance of alternating features in the timing of cardiac events is often a precursor to the initiation of serious cardiac arrhythmias. We treat embryonic chick cardiac cells with a potassium channel blocker, which leads to the initiation of alternating rhythms. We associate this transition with a mathematical instability, called a period-doubling bifurcation, in a model of the cardiac cells. Period-doubling bifurcations have been linked to the onset of abnormal alternating cardiac rhythms, which have been implicated in cardiac arrhythmias such as T-wave alternans and various tachycardias. Theory predicts that in the neighborhood of the transition, the system's dynamics slow down, leading to noise amplification and the manifestation of oscillations in the autocorrelation function. Examining the aggregates' interbeat intervals, we observe the oscillations in the autocorrelation function and noise amplification preceding the bifurcation. We analyze plots--termed return maps--that relate the current interbeat interval with the following interbeat interval. Based on these plots, we develop a quantitative measure using the slope of the return map to assess how close the system is to the bifurcation. Furthermore, the slope of the return map and the lag-1 autocorrelation coefficient are equal. Our results suggest that the slope and the lag-1 autocorrelation coefficient represent quantitative measures to predict the onset of abnormal alternating cardiac rhythms.