Model of intracellular calcium cycling in ventricular myocytes

A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in silico model, which we call the 'Hernandez-Hernandez model', of electrical and Ca2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular Ca2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type Ca2+ channels are predicted to drive sex-specific differences in intracellular Ca2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used Ca2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of antihypertensive drugs.

Intracellular ion accumulation in the genesis of complex action potential dynamics under cardiac diseases

Intracellular ions, including sodium (Na^{+}), calcium (Ca^{2+}), and potassium (K^{+}), etc., accumulate slowly after a change of the state of the heart, such as a change of the heart rate. The goal of this study is to understand the roles of slow ion accumulation in the genesis of cardiac memory and complex action-potential duration (APD) dynamics that can lead to lethal cardiac arrhythmias. We carry out numerical simulations of a detailed action potential model of ventricular myocytes under normal and diseased conditions, which exhibit memory effects and complex APD dynamics. We develop a low-dimensional iterated map (IM) model to describe the dynamics of Na^{+}, Ca^{2+}, and APD and use it to uncover the underlying dynamical mechanisms. The development of the IM model is informed by simulation results under the normal condition. We then use the IM model to perform linear stability analyses and computer simulations to investigate the bifurcations and complex APD dynamics, which depend on the feedback loops between APD and intracellular Ca^{2+} and Na^{+} concentrations and the steepness of the APD response to the ion concentrations. When the feedback between APD and Ca^{2+} concentration is positive, a Hopf bifurcation leading to periodic oscillatory behavior occurs as the steepness of the APD response to the ion concentrations increases. The negative feedback loop between APD and Na^{+} concentration is required for the Hopf bifurcation. When the feedback between APD and Ca^{2+} concentration is negative, period-doubling bifurcations leading to high periodicity and chaos occurs. In this case, Na^{+} accumulation plays little role in the dynamics. Finally, we carry out simulations of the detailed action potential model under two diseased conditions, which exhibit steep APD responses to ion concentrations. Under both conditions, Hopf bifurcations leading to slow oscillations or period-doubling bifurcations leading to high periodicity and chaotic APD dynamics occur, depending on the strength of the ion pump-Na^{+}-Ca^{2+} exchanger. Using functions reconstructed from the simulation data, the IM model accurately captures the bifurcations and dynamics under the two diseased conditions. In conclusion, besides using computer simulations of a detailed high-dimensional action-potential model to investigate the effects of slow ion accumulation and short-term memory on bifurcations and genesis of complex APD dynamics in cardiac myocytes under diseased conditions, this study also provides a low-dimensional mathematical tool, i.e., the IM model, to allow stability analyses for uncovering the underlying mechanisms.

A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in-silico model, which we call the "Hernandez-Hernandez model", of electrical and C a 2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular C a 2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type C a 2+ channels are predicted to drive sex-specific differences in intracellular C a 2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used C a 2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of anti-hypertensive drugs.

The Calcium Transient Coupled to the L-Type Calcium Current Attenuates Action Potential Alternans

Background

Action potential (AP) alternans are linked to increased arrhythmogenesis. It is suggested that calcium (Ca ²⁺ ) transient (CaT) alternans cause AP alternans through bi-directional coupling feedback mechanisms because CaT alternans can precede AP alternans and develop in AP alternans free conditions. However, the CaT is an emergent response to intracellular Ca ²⁺ handling, and the mechanisms linking AP and CaT alternans are still a topic of investigation. This study investigated the development of AP alternans in the absence of CaT.

Methods

AP (patch clamp) and intracellular Ca ²⁺ (Fluo-4 epifluorescence) were recorded simultaneously from isolated rabbit ventricle myocytes perfused with the intracellular Ca ²⁺ buffer BAPTA (10-20 mM) to abolish CaT and/or the L-type Ca2+ channel activator Bay K 8644 (25 nM).

Results

After a rate change, alternans were critically damped and stable, overdamped and ceased over seconds, underdamped with longer scale harmonics, or unstably underdamped progressing to 2:1 capture. Alternans in control cells were predominantly critically damped, but after CaT ablation with 10 or 20 mM BAPTA, exhibited respectively increased overdamping or increased underdamping. Alternans were easier to induce in CaT free cells as evidenced by a higher alternans threshold (ALT-TH: at least 7 pairs of alternating beats) relative to control cells. Alternans in Bay K 8644 treated cells were often underdamped, but the ALT-TH was similar to control. In CaT ablated cells, Bay K 8644 prolonged AP duration (APD) leading predominantly to unstably underdamped alternans.

Conclusions

AP alternans occur more readily in the absence of CaT suggesting that the CaT dampens the development of AP alternans. The data further demonstrate that agonizing the L-type calcium current without the negative feedback of the CaT leads to unstable alternans. This negative feedback mechanism may be important for understanding treatments aimed at reducing CaT or its dynamic response to prevent 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.

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.

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.

Optical Ultrastructure of Large Mammalian Hearts Recovers Discordant Alternans by In Silico Data Assimilation

Understanding and predicting the mechanisms promoting the onset and sustainability of cardiac arrhythmias represent a primary concern in the scientific and medical communities still today. Despite the long-lasting effort in clinical and physico-mathematical research, a critical aspect to be fully characterized and unveiled is represented by spatiotemporal alternans patterns of cardiac excitation. The identification of discordant alternans and higher-order alternating rhythms by advanced data analyses as well as their prediction by reliable mathematical models represents a major avenue of research for a broad and multidisciplinary scientific community. Current limitations concern two primary aspects: 1) robust and general-purpose feature extraction techniques and 2) in silico data assimilation within reliable and predictive mathematical models. Here, we address both aspects. At first, we extend our previous works on Fourier transformation imaging (FFI), applying the technique to whole-ventricle fluorescence optical mapping. Overall, we identify complex spatial patterns of voltage alternans and characterize higher-order rhythms by a frequency-series analysis. Then, we integrate the optical ultrastructure obtained by FFI analysis within a fine-tuned electrophysiological mathematical model of the cardiac action potential. We build up a novel data assimilation procedure demonstrating its reliability in reproducing complex alternans patterns in two-dimensional computational domains. Finally, we prove that the FFI approach applied to both experimental and simulated signals recovers the same information, thus closing the loop between the experiment, data analysis, and numerical simulations.

The development of L-type Ca2+ current mediated alternans does not depend on the restitution slope in canine ventricular myocardium

Cardiac alternans have crucial importance in the onset of ventricular fibrillation. The early explanation for alternans development was the voltage-driven mechanism, where the action potential (AP) restitution steepness was considered as crucial determining factor. Recent results suggest that restitution slope is an inadequate predictor for alternans development, but several studies still claim the role of membrane potential as underlying mechanism of alternans. These controversial data indicate that the relationship of restitution and alternans development is not completely understood. APs were measured by conventional microelectrode technique from canine right ventricular papillary muscles. Ionic currents combined with fluorescent measurements were recorded by patch-clamp technique. APs combined with fluorescent measurements were monitored by sharp microelectrodes. Rapid pacing evoked restitution-independent AP duration (APD) alternans. When non-alternating AP voltage command was used, Ca²⁺_i-transient (CaT) alternans were not observed. When alternating rectangular voltage pulses were applied, CaT alternans were proportional to I_CaL amplitude alternans. Selective I_CaL inhibition did not influence the fast phase of APD restitution. In this study we found that I_CaL has minor contribution in shaping the fast phase of restitution curve suggesting that I_CaL—if it plays important role in the alternans mechanism—could be an additional factor that attenuates the reliability of APD restitution slope to predict alternans.

Targeted Drug Delivery for Cardiovascular Disease: Modeling of Modulated Extracellular Vesicle Release Rates

Invasive and medical therapy has led to major improvements in cardiovascular disease management, but important challenges remain open. The discovery of a nano-sized system of extracellular vesicles (EVs) is opening new possibilities for reprogramming malfunctioning cells and indicates that EVs can be employed in therapeutic biomedical applications as engineered drug vehicles. Molecular communication (MC) has applications for treating cells with directed drug delivery, employing special targeting transmembrane proteins. In this paper, we propose a novel drug delivery system for cardiovascular diseases using an EV-mediated MC platform and exemplify the potential use in hypertrophic cardiomyopathy. We utilize intracellular calcium signaling as a natural mediator of EVs released from synthetic cells and model the release rate. We propose to use the cells as a therapeutic release system with a control signal input which modulates the EVs release rate as the output signal. We also study the frequency domain of the proposed model and estimate the transfer function of the therapeutic release system model numerically where the root-mean-square error for two separate estimated output signals are 0.0353 and 0.0124. The proposed EV-mediated targeted drug delivery system can make breakthroughs in future healthcare, in cardiovascular and other diseases where targeting is required.

A space-fractional bidomain framework for cardiac electrophysiology: 1D alternans dynamics

Cardiac electrophysiology modeling deals with a complex network of excitable cells forming an intricate syncytium: the heart. The electrical activity of the heart shows recurrent spatial patterns of activation, known as cardiac alternans, featuring multiscale emerging behavior. On these grounds, we propose a novel mathematical formulation for cardiac electrophysiology modeling and simulation incorporating spatially non-local couplings within a physiological reaction-diffusion scenario. In particular, we formulate, a space-fractional electrophysiological framework, extending and generalizing similar works conducted for the monodomain model. We characterize one-dimensional excitation patterns by performing an extended numerical analysis encompassing a broad spectrum of space-fractional derivative powers and various intra- and extracellular conductivity combinations. Our numerical study demonstrates that (i) symmetric properties occur in the conductivity parameters' space following the proposed theoretical framework, (ii) the degree of non-local coupling affects the onset and evolution of discordant alternans dynamics, and (iii) the theoretical framework fully recovers classical formulations and is amenable for parametric tuning relying on experimental conduction velocity and action potential morphology.

Dynamical analysis of early afterdepolarization patterns in a biophysically detailed cardiac model

Arrhythmogenic early afterdepolarizations (EADs) are investigated in a biophysically detailed mathematical model of a rabbit ventricular myocyte, providing their location in the parameter phase space and describing their dynamical mechanisms. Simulations using the Sato model, defined by 27 state variables and 177 parameters, are conducted to generate electrical action potentials (APs) for different values of the pacing cycle length and other parameters related to sodium and calcium concentrations. A detailed study of the different AP patterns with or without EADs is carried out, showing the presence of a high variety of temporal AP configurations with chaotic and quasiperiodic behaviors. Regions of bistability are identified and, importantly, linked to transitions between different behaviors. Using sweeping techniques, one-, two-, and three-parameter phase spaces are provided, allowing ascertainment of the role of the selected parameters as well as location of the transition regions. A Devil's staircase, with symbolic sequence analysis, is proposed to describe transitions in the ratio between the number of voltage (EAD and AP) peaks and the number of APs. To conclude, the obtained results are linked to recent studies for low-dimensional models and a conjecture is made for the internal dynamical structure of the transition region from non-EAD to EAD behavior using fold and cusp bifurcations and maximal canards.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Stress-driven cardiac calcium mishandling via a kinase-to-kinase crosstalk

Calcium homeostasis in the cardiomyocyte is critical to the regulation of normal cardiac function. Abnormal calcium dynamics such as altered uptake by the sarcoplasmic reticulum (SR) Ca²⁺-ATPase and increased diastolic SR calcium leak are involved in the development of maladaptive cardiac remodeling under pathological conditions. Ca²⁺/calmodulin-dependent protein kinase II-δ (CaMKIIδ) is a well-recognized key molecule in calcium dysregulation in cardiomyocytes. Elevated cellular stress is known as a common feature during pathological remodeling, and c-jun N-terminal kinase (JNK) is an important stress kinase that is activated in response to intrinsic and extrinsic stress stimuli. Our lab recently identified specific actions of JNK isoform 2 (JNK2) in CaMKIIδ expression, activation, and CaMKIIδ-dependent SR Ca²⁺ mishandling in the stressed heart. This review focuses on the current understanding of cardiac SR calcium handling under physiological and pathological conditions as well as the newly identified contribution of the stress kinase JNK2 in CaMKIIδ-dependent SR Ca²⁺ abnormal mishandling. The new findings identifying dual roles of JNK2 in CaMKIIδ expression and activation are also discussed in this review.

Controllability of voltage- and calcium-driven cardiac alternans in a map model

Certain cardiac arrhythmias are preceded by electrical alternans, a state characterized by beat-to-beat alternation in cellular action potential duration. Cardiac alternans may arise from different mechanisms including instabilities in voltage or intracellular calcium cycling. Although a number of techniques have been proposed to suppress alternans, these methods have mainly been tested using models that do not support calcium-driven alternans. Therefore, it is important to understand how control methods may perform when alternans is driven by instabilities in calcium cycling. In this study, we applied controllability analysis to a discrete map of alternans dynamics in a cardiac cell. We compared two different controllability measures to determine to what extent different control strategies could suppress alternans and tested these predictions using three feedback controllers. We found a modal controllability measure, unlike the minimum singular value of the controllability matrix, consistently indicated the control strategies requiring the least control effort and yielding the smallest closed-loop eigenvalue. In addition, action potential duration was identified as the most effective variable through which control can be applied, regardless of alternans mechanism, although sarcoplasmic reticulum calcium load was also useful for the calcium-driven alternans cases.

Mechanisms of ranolazine pretreatment in preventing ventricular tachyarrhythmias in diabetic db/db mice with acute regional ischemia-reperfusion injury

Studies have demonstrated that diabetic (db/db) mice have increased susceptibility to myocardial ischemia-reperfusion (IR) injury and ventricular tachyarrhythmias (VA). We aimed to investigate the antiarrhythmic and molecular mechanisms of ranolazine in db/db mouse hearts with acute IR injury. Ranolazine was administered for 1 week before coronary artery ligation. Diabetic db/db and control db/+ mice were divided into ranolazine-given and -nongiven groups. IR model was created by 15-min left coronary artery ligation and 10-min reperfusion. In vivo electrophysiological studies showed that the severity of VA inducibility was higher in db/db mice than control (db/ +) mice. Ranolazine suppressed the VA inducibility and severity. Optical mapping studies in Langendorff-perfused hearts showed that ranolazine significantly shortened action potential duration, Ca_i transient duration, Ca_i decay time, ameliorated conduction inhomogeneity, and suppressed arrhythmogenic alternans induction. Western blotting studies showed that the expression of pThr¹⁷-phospholamban, calsequestrin 2 and voltage-gated sodium channel in the IR zone was significantly downregulated in db/db mice, which was ameliorated with ranolazine pretreatment and might play a role in the anti-arrhythmic actions of ranolazine in db/db mouse hearts with IR injury.

Effect of constant-DI pacing on single cell pacing dynamics

Cardiac alternans, beat-to-beat alternations in action potential duration, is a precursor to fatal arrhythmias such as ventricular fibrillation. Previous research has shown that voltage driven alternans can be suppressed by application of a constant diastolic interval (DI) pacing protocol. However, the effect of constant-DI pacing on cardiac cell dynamics and its interaction with the intracellular calcium cycle remains to be determined. Therefore, we aimed to examine the effects of constant-DI pacing on the dynamical behavior of a single-cell numerical model of cardiac action potential and the influence of voltage-calcium (V-Ca) coupling on it. Single cell dynamics were analyzed in the vicinity of the bifurcation point using a hybrid pacing protocol, a combination of constant-basic cycle length (BCL) and constant-DI pacing. We demonstrated that in a small region beneath the bifurcation point, constant-DI pacing caused the cardiac cell to remain alternans-free after switching to the constant-BCL pacing, thus introducing a region of bistability (RB). The size of the RB increased with stronger V-Ca coupling and was diminished with weaker V-Ca coupling. Overall, our findings demonstrate that the application of constant-DI pacing on cardiac cells with strong V-Ca coupling may induce permanent changes to cardiac cell dynamics increasing the utility of constant-DI pacing.

Periplocymarin Plays an Efficacious Cardiotonic Role via Promoting Calcium Influx

Periplocymarin, which belongs to cardiac glycosides, is an effective component extracted from Periplocae Cortex. However, its cardiovascular effects remain unidentified. In the present study, injection of periplocymarin (5 mg/kg) through external jugular vein immediately increased the mean arterial pressure (MAP) in anesthetized C57BL/6 mice. Ex vivo experiments using mouse mesenteric artery rings were conducted to validate the role of periplocymarin on blood vessels. However, periplocymarin failed to induce vasoconstriction directly, and had no effects on vasoconstriction induced by phenylephrine (Phe) and angiotensin II (Ang II). In addition, vasodilatation induced by acetylcholine (Ach) was insusceptible to periplocymarin. Echocardiography was used to evaluate the effects of periplocymarin on cardiac function. The results showed that the injection of periplocymarin significantly increase the ejection fraction (EF) in mice without changing the heart rate. In vitro studies using isolated neonatal rat ventricular myocytes (NRVMs) revealed that periplocymarin transiently increased the intracellular Ca²⁺ concentration observed by confocal microscope. But in Ca²⁺-free buffer, this phenomenon vanished. Besides, inhibition of sodium potassium-activated adenosine triphosphatase (Na⁺-K⁺-ATPase) by digoxin significantly suppressed the increase of MAP and EF in mice, and the influx of Ca²⁺ in cardiomyocytes, mediated by periplocymarin. Collectively, these findings demonstrated that periplocymarin increased the contractility of myocardium by promoting the Ca²⁺ influx of cardiomyocytes via targeting on Na⁺-K⁺-ATPase, which indirectly led to the instantaneous rise of blood pressure.

Impact of ISK Voltage and Ca2+/Mg2+-Dependent Rectification on Cardiac Repolarization

Cardiac small conductance Ca²⁺-activated K⁺ (SK) channels are activated solely by Ca²⁺, but the SK current (I_SK) is inwardly rectified. However, the impact of inward rectification in shaping action potentials (APs) in ventricular cardiomyocytes under β-adrenergic stimulation or in disease states remains undefined. Two processes underlie this inward rectification: an intrinsic rectification caused by an electrostatic energy barrier from positively charged amino acids at the inner pore and a voltage-dependent Ca²⁺/Mg²⁺ block. Thus, Ca²⁺ has a biphasic effect on I_SK, activating at low [Ca²⁺] yet inhibiting I_SK at high [Ca²⁺]. We examined the effect of I_SK rectification on APs in rat cardiomyocytes by simultaneously recording whole-cell apamin-sensitive currents and Ca²⁺ transients during an AP waveform and developed a computer model of SK channels with rectification features. The typical profile of I_SK during AP clamp included an initial peak (mean 1.6 pA/pF) followed by decay to the point that submembrane [Ca²⁺] reached ∼10 μM. During the rest of the AP stimulus, I_SK either plateaued or gradually increased as the cell repolarized and submembrane [Ca²⁺] decreased further. We used a six-state gating model combined with intrinsic and Ca²⁺/Mg²⁺-dependent rectification to simulate I_SK and investigated the relative contributions of each type of rectification to AP shape. This SK channel model replicates key features of I_SK recording during AP clamp showing that intrinsic rectification limits I_SK at high V_m during the early and plateau phase of APs. Furthermore, the initial rise of Ca²⁺ transients activates, but higher [Ca²⁺] blocks SK channels, yielding a transient outward-like I_SK trajectory. During the decay phase of Ca²⁺, the Ca²⁺-dependent block is released, causing I_SK to rise again and contribute to repolarization. Therefore, I_SK is an important repolarizing current, and the rectification characteristics of an SK channel determine its impact on early, plateau, and repolarization phases of APs.

Sensitivity of a data-assimilation system for reconstructing three-dimensional cardiac electrical dynamics

Modelling of cardiac electrical behaviour has led to important mechanistic insights, but important challenges, including uncertainty in model formulations and parameter values, make it difficult to obtain quantitatively accurate results. An alternative approach is combining models with observations from experiments to produce a data-informed reconstruction of system states over time. Here, we extend our earlier data-assimilation studies using an ensemble Kalman filter to reconstruct a three-dimensional time series of states with complex spatio-temporal dynamics using only surface observations of voltage. We consider the effects of several algorithmic and model parameters on the accuracy of reconstructions of known scroll-wave truth states using synthetic observations. In particular, we study the algorithm's sensitivity to parameters governing different parts of the process and its robustness to several model-error conditions. We find that the algorithm can achieve an acceptable level of error in many cases, with the weakest performance occurring for model-error cases and more extreme parameter regimes with more complex dynamics. Analysis of the poorest-performing cases indicates an initial decrease in error followed by an increase when the ensemble spread is reduced. Our results suggest avenues for further improvement through increasing ensemble spread by incorporating additive inflation or using a parameter or multi-model ensemble. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

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.

Feedback control of calcium driven alternans in cardiac myocytes

Cardiac alternans is a beat-to-beat alternation of the action potential duration (APD), which has been implicated as a possible cause of ventricular fibrillation. Previous studies have shown that alternans can originate via a period doubling bifurcation caused by the nonlinear dependence of the APD on the previous diastolic interval. In this case, it has been demonstrated that alternans can be eliminated by applying feedback control on the pacing cycle length. However, studies have shown that alternans can also originate due to unstable calcium (Ca) cycling in cardiac myocytes. In this study, we explore the effectiveness of APD feedback control to suppress alternans when the underlying instability is due to unstable Ca cycling. In particular, we explore the role of the bi-directional coupling between Ca and voltage and determine the effectiveness of feedback control under a wide range of conditions. We also analyze the applicability of feedback control on a coupled two cell system and show that APD control induces spatially out-of-phase alternans. We analyze the onset and the necessary conditions for the emergence of these out-of-phase patterns and assess the effectiveness of feedback control to suppress Ca driven alternans in a multi-cellular system.

Calcium Handling Defects and Cardiac Arrhythmia Syndromes

Calcium ions (Ca²⁺) play a major role in the cardiac excitation-contraction coupling. Intracellular Ca²⁺ concentration increases during systole and falls in diastole thereby determining cardiac contraction and relaxation. Normal cardiac function also requires perfect organization of the ion currents at the cellular level to drive action potentials and to maintain action potential propagation and electrical homogeneity at the tissue level. Any imbalance in Ca²⁺ homeostasis of a cardiac myocyte can lead to electrical disturbances. This review aims to discuss cardiac physiology and pathophysiology from the elementary membrane processes that can cause the electrical instability of the ventricular myocytes through intracellular Ca²⁺ handling maladies to inherited and acquired arrhythmias. Finally, the paper will discuss the current therapeutic approaches targeting cardiac arrhythmias.

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.

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.

The Ca2+ transient as a feedback sensor controlling cardiomyocyte ionic conductances in mouse populations

Conductances of ion channels and transporters controlling cardiac excitation may vary in a population of subjects with different cardiac gene expression patterns. However, the amount of variability and its origin are not quantitatively known. We propose a new conceptual approach to predict this variability that consists of finding combinations of conductances generating a normal intracellular Ca²⁺ transient without any constraint on the action potential. Furthermore, we validate experimentally its predictions using the Hybrid Mouse Diversity Panel, a model system of genetically diverse mouse strains that allows us to quantify inter-subject versus intra-subject variability. The method predicts that conductances of inward Ca²⁺ and outward K⁺ currents compensate each other to generate a normal Ca²⁺ transient in good quantitative agreement with current measurements in ventricular myocytes from hearts of different isogenic strains. Our results suggest that a feedback mechanism sensing the aggregate Ca²⁺ transient of the heart suffices to regulate ionic conductances.

A model for cooperative gating of L-type Ca2+ channels and its effects on cardiac alternans dynamics

In ventricular myocytes, membrane depolarization during the action potential (AP) causes synchronous activation of multiple L-type CaV1.2 channels (LTCCs), which trigger the release of calcium (Ca2+) from the sarcoplasmic reticulum (SR). This results in an increase in intracellular Ca2+ (Cai) that initiates contraction. During pulsus alternans, cardiac contraction is unstable, going from weak to strong in successive beats despite a constant heart rate. These cardiac alternans can be caused by the instability of membrane potential (Vm) due to steep AP duration (APD) restitution (Vm-driven alternans), instability of Cai cycling (Ca2+-driven alternans), or both, and may be modulated by functional coupling between clustered CaV1.2 (e.g. cooperative gating). Here, mathematical analysis and computational models were used to determine how changes in the strength of cooperative gating between LTCCs may impact membrane voltage and intracellular Ca2+ dynamics in the heart. We found that increasing the degree of coupling between LTCCs increases the amplitude of Ca2+ currents (ICaL) and prolongs AP duration (APD). Increased AP duration is known to promote cardiac alternans, a potentially arrhythmogenic substrate. In addition, our analysis shows that increasing the strength of cooperative activation of LTCCs makes the coupling of Ca2+ on the membrane voltage (Cai→Vm coupling) more positive and destabilizes the Vm-Cai dynamics for Vm-driven alternans and Cai-driven alternans, but not for quasiperiodic oscillation. These results suggest that cooperative gating of LTCCs may have a major impact on cardiac excitation-contraction coupling, not only by prolonging APD, but also by altering Cai→Vm coupling and potentially promoting cardiac arrhythmias.

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 model of cardiac ryanodine receptor gating predicts experimental Ca2+-dynamics and Ca2+-triggered arrhythmia in the long QT syndrome

Abnormal Ca²⁺ handling is well-established as the trigger of cardiac arrhythmia in catecholaminergic polymorphic ventricular tachycardia and digoxin toxicity, but its role remains controversial in Torsade de Pointes (TdP), the arrhythmia associated with the long QT syndrome (LQTS). Recent experimental results show that early afterdepolarizations (EADs) that initiate TdP are caused by spontaneous (non-voltage-triggered) Ca²⁺ release from Ca²⁺-overloaded sarcoplasmic reticulum (SR) rather than the activation of the L-type Ca²⁺-channel window current. In bradycardia and long QT type 2 (LQT2), a second, non-voltage triggered cytosolic Ca²⁺ elevation increases gradually in amplitude, occurs before overt voltage instability, and then precedes the rise of EADs. Here, we used a modified Shannon-Puglisi-Bers model of rabbit ventricular myocytes to reproduce experimental Ca²⁺ dynamics in bradycardia and LQT2. Abnormal systolic Ca²⁺-oscillations and EADs caused by SR Ca²⁺-release are reproduced in a modified 0-dimensional model, where 3 gates in series control the ryanodine receptor (RyR2) conductance. Two gates control RyR2 activation and inactivation and sense cytosolic Ca²⁺ while a third gate senses luminal junctional SR Ca²⁺. The model predicts EADs in bradycardia and low extracellular [K⁺] and cessation of SR Ca²⁺-release terminate salvos of EADs. Ca²⁺-waves, systolic cell-synchronous Ca²⁺-release, and multifocal diastolic Ca²⁺ release seen in subcellular Ca²⁺-mapping experiments are observed in the 2-dimensional version of the model. These results support the role of SR Ca²⁺-overload, abnormal SR Ca²⁺-release, and the subsequent activation of the electrogenic Na⁺/Ca²⁺-exchanger as the mechanism of TdP. The model offers new insights into the genesis of cardiac arrhythmia and new therapeutic strategies.

Distinguishing mechanisms for alternans in cardiac cells using constant-diastolic-interval pacing

Alternans, a proarrhythmic dynamical state in which cardiac action potentials alternate between long and short durations despite a constant pacing period, traditionally has been explained at the cellular level using nonlinear dynamics principles under the assumption that the action potential duration (APD) is determined solely by the time elapsed since the end of the previous action potential, called the diastolic interval (DI). In this scenario, APDs at a steady state should be the same provided that the preceding DIs are the same. Nevertheless, experiments attempting to eliminate alternans by dynamically adjusting the timing of pacing stimuli to keep the DI constant showed that alternans persisted, contradicting the traditional theory. It is now widely known that alternans also can arise from a different mechanism associated with intracellular calcium cycling. Our goal is to determine whether intracellular calcium dynamics can explain the experimental findings regarding the persistence of alternans despite a constant DI. For this, we use mathematical models capable of producing alternans through both voltage- and calcium-mediated mechanisms. We show that for voltage-driven alternans, action potentials elicited from a constant-DI protocol are always the same. However, in the case of calcium-driven alternans, the constant-DI protocol can result in alternans. Reducing the strength of the calcium instability progressively reduces and finally eliminates constant-DI alternans. Our findings suggest that screening for the presence of alternans using a constant-DI protocol has the potential for differentiating between voltage-driven and calcium-driven alternans.

Mechanism for Triggered Waves in Atrial Myocytes

Excitation-contraction coupling in atrial cells is mediated by calcium (Ca) signaling between L-type Ca channels and Ryanodine receptors that occurs mainly at the cell boundary. This unique architecture dictates essential aspects of Ca signaling under both normal and diseased conditions. In this study we apply laser scanning confocal microscopy, along with an experimentally based computational model, to understand the Ca cycling dynamics of an atrial cell subjected to rapid pacing. Our main finding is that when an atrial cell is paced under Ca overload conditions, Ca waves can then nucleate on the cell boundary and propagate to the cell interior. These propagating Ca waves are referred to as "triggered waves" because they are initiated by L-type Ca channel openings during the action potential. These excitations are distinct from spontaneous Ca waves originating from random fluctuations of Ryanodine receptor channels, and which occur after much longer waiting times. Furthermore, we argue that the onset of these triggered waves is a highly nonlinear function of the sarcoplasmic reticulum Ca load. This strong nonlinearity leads to aperiodic response of Ca at rapid pacing rates that is caused by the complex interplay between paced Ca release and triggered waves. We argue further that this feature of atrial cells leads to dynamic instabilities that may underlie atrial arrhythmias. These studies will serve as a starting point to explore the nonlinear dynamics of atrial cells and will yield insights into the trigger and maintenance of atrial fibrillation.

Dynamics of spatiotemporal line defects and chaos control in complex excitable systems

Spatiotemporal pattern formation governs dynamics and functions in various biological systems. In the heart, excitable waves can form complex oscillatory and chaotic patterns even at an abnormally higher frequency than normal heart beats, which increase the risk of fatal heart conditions by inhibiting normal blood circulation. Previous studies suggested that line defects (nodal lines) play a critical role in stabilizing those undesirable patterns. However, it remains unknown if the line defects are static or dynamically changing structures in heart tissue. Through in vitro experiments of heart tissue observation, we reveal the spatiotemporal dynamics of line defects in rotating spiral waves. We combined a novel signaling over-sampling technique with a multi-dimensional Fourier analysis, showing that line defects can translate, merge, collapse and form stable singularities with even and odd parity while maintaining a stable oscillation of the spiral wave in the tissue. These findings provide insights into a broad class of complex periodic systems, with particular impact to the control and understanding of heart diseases.

Alternans in atria: Mechanisms and clinical relevance

Atrial fibrillation is the most common sustained arrhythmia and its prevalence is rapidly rising with the aging of the population. Cardiac alternans, defined as cyclic beat-to-beat alternations in contraction force, action potential (AP) duration and intracellular Ca²⁺ release at constant stimulation rate, has been associated with the development of ventricular arrhythmias. Recent clinical data also provide strong evidence that alternans plays a central role in arrhythmogenesis in atria. The aim of this article is to review the mechanisms that are responsible for repolarization alternans and contribute to the transition from spatially concordant alternans to the more arrhythmogenic spatially discordant alternans in atria.

Dynamical effects of calcium-sensitive potassium currents on voltage and calcium alternans

KEY POINTS

A mathematical model of a small conductance Ca²⁺ -activated potassium (SK) channel was developed and incorporated into a physiologically detailed ventricular myocyte model. Ca²⁺ -sensitive K⁺ currents promote negative intracellular Ca²⁺ to membrane voltage (CA_i²⁺ → V_m ) coupling. Increase of Ca²⁺ -sensitive K⁺ currents can be responsible for electromechanically discordant alternans and quasiperiodic oscillations at the cellular level. At the tissue level, Turing-type instability can occur when Ca²⁺ -sensitive K⁺ currents are increased.

ABSTRACT

Cardiac alternans is a precursor to life-threatening arrhythmias. Alternans can be caused by instability of the membrane voltage (V_m ), instability of the intracellular Ca²⁺ ( Ca i2+) cycling, or both. V_m dynamics and Ca i2+ dynamics are coupled via Ca²⁺ -sensitive currents. In cardiac myocytes, there are several Ca²⁺ -sensitive potassium (K⁺ ) currents such as the slowly activating delayed rectifier current (I_Ks ) and the small conductance Ca²⁺ -activated potassium (SK) current (I_SK ). However, the role of these currents in the development of arrhythmias is not well understood. In this study, we investigated how these currents affect voltage and Ca²⁺ alternans using a physiologically detailed computational model of the ventricular myocyte and mathematical analysis. We define the coupling between V_m and Ca i2+ cycling dynamics ( Ca i2+→V_m coupling) as positive (negative) when a larger Ca²⁺ transient at a given beat prolongs (shortens) the action potential duration (APD) of that beat. While positive coupling predominates at baseline, increasing I_Ks and I_SK promote negative Ca i2+→V_m coupling at the cellular level. Specifically, when alternans is Ca²⁺ -driven, electromechanically (APD-Ca²⁺ ) concordant alternans becomes electromechanically discordant alternans as I_Ks or I_SK increase. These cellular level dynamics lead to different types of spatially discordant alternans in tissue. These findings help to shed light on the underlying mechanisms of cardiac alternans especially when the relative strength of these currents becomes larger under pathological conditions or drug administrations.

New Insights in Cardiac Calcium Handling and Excitation-Contraction Coupling

Excitation-contraction (EC) coupling denotes the conversion of electric stimulus in mechanic output in contractile cells. Several studies have demonstrated that calcium (Ca²⁺) plays a pivotal role in this process. Here we present a comprehensive and updated description of the main systems involved in cardiac Ca²⁺ handling that ensure a functional EC coupling and their pathological alterations, mainly related to heart failure.

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.

Rabbit-specific computational modelling of ventricular cell electrophysiology: Using populations of models to explore variability in the response to ischemia

Computational modelling, combined with experimental investigations, is a powerful method for investigating complex cardiac electrophysiological behaviour. The use of rabbit-specific models, due to the similarities of cardiac electrophysiology in this species with human, is especially prevalent. In this paper, we first briefly review rabbit-specific computational modelling of ventricular cell electrophysiology, multi-cellular simulations including cellular heterogeneity, and acute ischemia. This mini-review is followed by an original computational investigation of variability in the electrophysiological response of two experimentally-calibrated populations of rabbit-specific ventricular myocyte action potential models to acute ischemia. We performed a systematic exploration of the response of the model populations to varying degrees of ischemia and individual ischemic parameters, to investigate their individual and combined effects on action potential duration and refractoriness. This revealed complex interactions between model population variability and ischemic factors, which combined to enhance variability during ischemia. This represents an important step towards an improved understanding of the role that physiological variability may play in electrophysiological alterations during acute ischemia.

Targeting the late component of the cardiac L-type Ca2+ current to suppress early afterdepolarizations

Early afterdepolarizations (EADs) associated with prolongation of the cardiac action potential (AP) can create heterogeneity of repolarization and premature extrasystoles, triggering focal and reentrant arrhythmias. Because the L-type Ca(2+) current (ICa,L) plays a key role in both AP prolongation and EAD formation, L-type Ca(2+) channels (LTCCs) represent a promising therapeutic target to normalize AP duration (APD) and suppress EADs and their arrhythmogenic consequences. We used the dynamic-clamp technique to systematically explore how the biophysical properties of LTCCs could be modified to normalize APD and suppress EADs without impairing excitation-contraction coupling. Isolated rabbit ventricular myocytes were first exposed to H2O2 or moderate hypokalemia to induce EADs, after which their endogenous ICa,L was replaced by a virtual ICa,L with tunable parameters, in dynamic-clamp mode. We probed the sensitivity of EADs to changes in the (a) amplitude of the noninactivating pedestal current; (b) slope of voltage-dependent activation; (c) slope of voltage-dependent inactivation; (d) time constant of voltage-dependent activation; and (e) time constant of voltage-dependent inactivation. We found that reducing the amplitude of the noninactivating pedestal component of ICa,L effectively suppressed both H2O2- and hypokalemia-induced EADs and restored APD. These results, together with our previous work, demonstrate the potential of this hybrid experimental-computational approach to guide drug discovery or gene therapy strategies by identifying and targeting selective properties of LTCC.

Cardiac alternans and intracellular calcium cycling

Cardiac alternans refers to a condition in which there is a periodic beat-to-beat oscillation in electrical activity and the strength of cardiac muscle contraction at a constant heart rate. Clinically, cardiac alternans occurs in settings that are typical for cardiac arrhythmias and has been causally linked to these conditions. At the cellular level, alternans is defined as beat-to-beat alternations in contraction amplitude (mechanical alternans), action potential duration (APD; electrical or APD alternans) and Ca(2+) transient amplitude (Ca(2+) alternans). The cause of alternans is multifactorial; however, alternans always originate from disturbances of the bidirectional coupling between membrane voltage (Vm ) and intracellular calcium ([Ca(2+) ]i ). Bidirectional coupling refers to the fact that, in cardiac cells, Vm depolarization and the generation of action potentials cause the elevation of [Ca(2+) ]i that is required for contraction (a process referred to as excitation-contraction coupling); conversely, changes of [Ca(2+) ]i control Vm because important membrane currents are Ca(2+) dependent. Evidence is mounting that alternans is ultimately caused by disturbances of cellular Ca(2+) signalling. Herein we review how two key factors of cardiac cellular Ca(2+) cycling, namely the release of Ca(2+) from internal stores and the capability of clearing the cytosol from Ca(2+) after each beat, determine the conditions under which alternans occurs. The contributions from key Ca(2+) -handling proteins (i.e. surface membrane channels, ion pumps and transporters and internal Ca(2+) release channels) are discussed.

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

Electromechanical alternans is a beat-to-beat alternation in the strength of contraction of a cardiac cell, which can be caused by an instability of calcium cycling. Using a distributed model of subcellular calcium we show that alternans occurs via an order-disorder phase transition which exhibits critical slowing down and a diverging correlation length. We apply finite size scaling along with a mapping to a stochastic coupled map model, to show that this transition in two dimensions is characterized by critical exponents consistent with the Ising universality class. These findings highlight the important role of cooperativity in biological cells, and suggest novel approaches to investigate the onset of the alternans instability in the heart.

The mechanisms of calcium cycling and action potential dynamics in cardiac alternans

RATIONALE

Alternans is a risk factor for cardiac arrhythmia, including atrial fibrillation. At the cellular level alternans manifests as beat-to-beat alternations in contraction, action potential duration (APD), and magnitude of the Ca(2+) transient (CaT). Electromechanical and CaT alternans are highly correlated, however, it has remained controversial whether the primary cause of alternans is a disturbance of cellular Ca(2+) signaling or electrical membrane properties.

OBJECTIVE

To determine whether a primary failure of intracellular Ca(2+) regulation or disturbances in membrane potential and AP regulation are responsible for the occurrence of alternans in atrial myocytes.

METHODS AND RESULTS

Pacing-induced APD and CaT alternans were studied in single rabbit atrial and ventricular myocytes using combined [Ca(2+)]i and electrophysiological measurements. In current-clamp experiments, APD and CaT alternans strongly correlated in time and magnitude. CaT alternans was observed without alternation in L-type Ca(2+) current, however, elimination of intracellular Ca(2+) release abolished APD alternans, indicating that [Ca(2+)]i dynamics have a profound effect on the occurrence of CaT alternans. Trains of 2 distinctive voltage commands in form of APs recorded during large and small alternans CaTs were applied to voltage-clamped cells. CaT alternans was observed with and without alternation in the voltage command shape. During alternans AP-clamp large CaTs coincided with both long and short AP waveforms, indicating that CaT alternans develop irrespective of AP dynamics.

CONCLUSIONS

The primary mechanism underlying alternans in atrial cells, similarly to ventricular cells, resides in a disturbance of Ca(2+) signaling, whereas APD alternans are a secondary consequence, mediated by Ca(2+)-dependent AP modulation.

Disrupted calcium release as a mechanism for atrial alternans associated with human atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia, but our knowledge of the arrhythmogenic substrate is incomplete. Alternans, the beat-to-beat alternation in the shape of cardiac electrical signals, typically occurs at fast heart rates and leads to arrhythmia. However, atrial alternans have been observed at slower pacing rates in AF patients than in controls, suggesting that increased vulnerability to arrhythmia in AF patients may be due to the proarrythmic influence of alternans at these slower rates. As such, alternans may present a useful therapeutic target for the treatment and prevention of AF, but the mechanism underlying alternans occurrence in AF patients at heart rates near rest is unknown. The goal of this study was to determine how cellular changes that occur in human AF affect the appearance of alternans at heart rates near rest. To achieve this, we developed a computational model of human atrial tissue incorporating electrophysiological remodeling associated with chronic AF (cAF) and performed parameter sensitivity analysis of ionic model parameters to determine which cellular changes led to alternans. Of the 20 parameters tested, only decreasing the ryanodine receptor (RyR) inactivation rate constant (ki_Ca) produced action potential duration (APD) alternans seen clinically at slower pacing rates. Using single-cell clamps of voltage, fluxes, and state variables, we determined that alternans onset was Ca²⁺-driven rather than voltage-driven and occurred as a result of decreased RyR inactivation which led to increased steepness of the sarcoplasmic reticulum (SR) Ca²⁺ release slope. Iterated map analysis revealed that because SR Ca²⁺ uptake efficiency was much higher in control atrial cells than in cAF cells, drastic reductions in ki_Ca were required to produce alternans at comparable pacing rates in control atrial cells. These findings suggest that RyR kinetics may play a critical role in altered Ca²⁺ homeostasis which drives proarrhythmic APD alternans in patients with AF.

Atrial fibrillation is an irregular heart rhythm affecting millions of people worldwide. Effective treatment of this cardiac disorder relies upon our detailed knowledge and understanding of the mechanisms that lead to arrhythmia. Recent clinical observations have suggested that alternans, a phenomenon where the shape of the electrical signal in the heart alternates from beat to beat, may play an important role in this process, but the underlying mechanisms remain unknown. In this study, we use computational models to conduct a detailed examination of the causes and contributors to alternans associated with human atrial fibrillation. We find that in atria remodeled by atrial fibrillation, alternans appears near resting heart rates because several aspects of calcium cycling are disrupted in the atrial cells. In particular, the release and uptake of calcium from the cellular storage compartment, the sarcoplasmic reticulum, becomes imbalanced, leading to alternation in calcium signals from beat to beat. These findings provide important insights into the mechanisms of proarrhythmic alternans in human atrial fibrillation which may be used to develop novel therapeutic targets and treatment strategies in the future.

Coexisting chaotic and multi-periodic dynamics in a model of cardiac alternans

The spatiotemporal dynamics of cardiac tissue is an active area of research for biologists, physicists, and mathematicians. Of particular interest is the study of period-doubling bifurcations and chaos due to their link with cardiac arrhythmogenesis. In this paper, we study the spatiotemporal dynamics of a recently developed model for calcium-driven alternans in a one dimensional cable of tissue. In particular, we observe in the cable coexistence of regions with chaotic and multi-periodic dynamics over wide ranges of parameters. We study these dynamics using global and local Lyapunov exponents and spatial trajectory correlations. Interestingly, near nodes-or phase reversals-low-periodic dynamics prevail, while away from the nodes, the dynamics tend to be higher-periodic and eventually chaotic. Finally, we show that similar coexisting multi-periodic and chaotic dynamics can also be observed in a detailed ionic model.

Data analysis in cardiac arrhythmias

Cardiac arrhythmias are an increasingly present in developed countries and represent a major health and economic burden. The occurrence of cardiac arrhythmias is closely linked to the electrical function of the heart. Consequently, the analysis of the electrical signal generated by the heart tissue, either recorded invasively or noninvasively, provides valuable information for the study of cardiac arrhythmias. In this chapter, novel cardiac signal analysis techniques that allow the study and diagnosis of cardiac arrhythmias are described, with emphasis on cardiac mapping which allows for spatiotemporal analysis of cardiac signals.Cardiac mapping can serve as a diagnostic tool by recording cardiac signals either in close contact to the heart tissue or noninvasively from the body surface, and allows the identification of cardiac sites responsible of the development or maintenance of arrhythmias. Cardiac mapping can also be used for research in cardiac arrhythmias in order to understand their mechanisms. For this purpose, both synthetic signals generated by computer simulations and animal experimental models allow for more controlled physiological conditions and complete access to the organ.

Cellular mechanism of cardiac alternans: an unresolved chicken or egg problem

T-wave alternans, a specific form of cardiac alternans, has been associated with the increased susceptibility to cardiac arrhythmias and sudden cardiac death (SCD). Plenty of evidence has related cardiac alternans at the tissue level to the instability of voltage kinetics or Ca(2+) handling dynamics at the cellular level. However, to date, none of the existing experiments could identify the exact cellular mechanism of cardiac alternans due to the bi-directional coupling between voltage kinetics and Ca(2+) handling dynamics. Either of these systems could be the origin of alternans and the other follows as a secondary change, therefore making the cellular mechanism of alternans a difficult chicken or egg problem. In this context, theoretical analysis combined with experimental techniques provides a possibility to explore this problem. In this review, we will summarize the experimental and theoretical advances in understanding the cellular mechanism of alternans. We focus on the roles of action potential duration (APD) restitution and Ca(2+) handling dynamics in the genesis of alternans and show how the theoretical analysis combined with experimental techniques has provided us a new insight into the cellular mechanism of alternans. We also discuss the possible reasons of increased propensity for alternans in heart failure (HF) and the new possible therapeutic targets. Finally, according to the level of electrophysiological recording techniques and theoretical strategies, we list some critical experimental or theoretical challenges which may help to determine the origin of alternans and to find more effective therapeutic targets in the future.

T-tubule disruption promotes calcium alternans in failing ventricular myocytes: mechanistic insights from computational modeling

In heart failure (HF), T-tubule (TT) disruption contributes to dyssynchronous calcium (Ca) release and impaired contraction, but its role in arrhythmogenesis remains unclear. In this study, we investigate the effects of TT disruption and other HF remodeling factors on Ca alternans in ventricular myocytes using computer modeling. A ventricular myocyte model with detailed spatiotemporal Ca cycling modeled by a coupled Ca release unit (CRU) network was used, in which the L-type Ca channels and the ryanodine receptor (RyR) channels were simulated by random Markov transitions. TT disruption, which removes the L-type Ca channels from the associated CRUs, results in "orphaned" RyR clusters and thus provides increased opportunity for spark-induced Ca sparks to occur. This effect combined with other HF remodeling factors promoted alternans by two distinct mechanisms: 1) for normal sarco-endoplasmic reticulum Ca ATPase (SERCA) activity, alternans was caused by both CRU refractoriness and coupling. The increased opportunity for spark-induced sparks by TT disruption combined with the enhanced CRU coupling by Ca elevation in the presence or absence of increased RyR leakiness facilitated spark synchronization on alternate beats to promote Ca alternans; 2) for down-regulated SERCA, alternans was caused by the sarcoplasmic reticulum (SR) Ca load-dependent mechanism, independent of CRU refractoriness. TT disruption and increased RyR leakiness shifted and steepened the SR Ca release-load relationship, which combines with down-regulated SERCA to promote Ca alternans. In conclusion, the mechanisms of Ca alternans for normal and down-regulated SERCA are different, and TT disruption promotes Ca alternans by both mechanisms, which may contribute to alternans at different stages of HF.