A model of cardiac ryanodine receptor gating predicts experimental Ca2+-dynamics and Ca2+-triggered arrhythmia in the long QT syndrome

Dissecting the roles of calcium cycling and its coupling with voltage in the genesis of early afterdepolarizations in cardiac myocyte models

Early afterdepolarizations (EADs) are abnormal depolarizations during the plateau phase of the action potential, which are known to be associated with lethal arrhythmias in the heart. There are two major hypotheses for EAD genesis based on experimental observations, i.e., the voltage (Vm)-driven and intracellular calcium (Ca)-driven mechanisms. In ventricular myocytes, Ca and Vm are bidirectionally coupled, which can affect each other's dynamics and result in new dynamics, however, the roles of Ca cycling and its coupling with Vm in the genesis of EADs have not been well understood. In this study, we use an action potential model that is capable of independent Vm and Ca oscillations to investigate the roles of Vm and Ca coupling in EAD genesis. Four different mechanisms of EADs are identified, which are either driven by Vm oscillations or Ca oscillations alone, or oscillations caused by their interactions. We also use 5 other ventricular action potential models to assess these EAD mechanisms and show that EADs in these models are mainly Vm-driven. These mechanistic insights from our simulations provide a theoretical base for understanding experimentally observed EADs and EAD-related arrhythmogenesis.

Activation of small conductance Ca2+ -activated K+ channels suppresses Ca2+ transient and action potential alternans in ventricular myocytes

At the cellular level, cardiac alternans is observed as beat-to-beat alternations in contraction strength, action potential (AP) morphology and Ca2+ transient (CaT) amplitude, and is a risk factor for cardiac arrhythmia. The (patho)physiological roles of small conductance Ca2+ -activated K+ (SK) channels in ventricles are poorly understood. We tested the hypothesis that in single rabbit ventricular myocytes pharmacological modulation of SK channels plays a causative role for the development of pacing-induced CaT and AP duration (APD) alternans. SK channel blockers (apamin, UCL1684) had only a minor effect on AP repolarization. However, SK channel activation by NS309 resulted in significant APD shortening, demonstrating that functional SK channels are well expressed in ventricular myocytes. The effects of NS309 were prevented or reversed by apamin and UCL1684, indicating that NS309 acted on SK channels. SK channel activation abolished or reduced the degree of pacing-induced CaT and APD alternans. Inhibition of KV 7.1 (with HMR1556) and KV 11.1 (with E4031) channels was used to mimic conditions of long QT syndromes type-1 and type-2, respectively. Both HMR1556 and E4031 enhanced CaT alternans that was prevented by SK channel activation. In AP voltage-clamped cells the SK channel activator had no effect on CaT alternans, confirming that suppression of CaT alternans was caused by APD shortening. APD shortening contributed to protection from alternans by lowering sarcoplasmic reticulum Ca2+ content and curtailing Ca2+ release. The data suggest that SK activation could be a potential intervention to avert development of alternans with important ramifications for arrhythmia prevention and therapy for patients with long QT syndrome. KEY POINTS: At the cellular level, cardiac alternans is observed as beat-to-beat alternations in contraction strength, action potential (AP) morphology and intracellular Ca2+ release amplitude, and is a risk factor for cardiac arrhythmia. The (patho)physiological roles of small conductance Ca2+ -activated K+ (SK) channels in ventricles are poorly understood. We investigated whether pharmacological modulation of SK channels affects the development of cardiac alternans in normal ventricular cells and in cells with drug-induced long QT syndrome (LQTS). While SK channel blockers have only a minor effect on AP morphology, their activation leads to AP shortening and abolishes or reduces the degree of pacing-induced Ca2+ and AP alternans. AP shortening contributed to protection against alternans by lowering sarcoplasmic reticulum Ca2+ content and curtailing Ca2+ release. The data suggest SK activation as a potential intervention to avert the development of alternans with important ramifications for arrhythmia prevention for patients with LQTS.

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.

Effects of Sacubitril/Valsartan on the Expression of CaMKII/Cav1.2 in Atrial Fibrillation Stimulation Rabbit Model

Background and Objective. Atrial fibrillation (AF) is linked to high morbidity and death rates throughout the world due to limited therapeutic options and thus presents a major challenge to the developed and developing countries. In this study, we aim to investigate the influence of sacubitril/valsartan (sac/val) treatment on the calmodulin-dependent protein kinase II (CaMKII)/Cav1.2 expression in AF models. Methods. Overall, 18 rabbits were randomly divided into control, pacing (600 beats/min), and pacing+sac/val groups. The rabbits in the pacing+sac/val cohort received oral sac/val (10 mg/kg twice daily) across the 21-day investigation period. After three weeks, the atrial effective refractory period (AERP) and AF induction rate were compared. HL-1 cultures were exposed to fast pacing (24 h) with and without LBQ657 (active sacubitril form)/valsartan. Western blots were used for detecting Cav1.2 and CaMKII expression within atrial muscles of the rabbits and HL-1 cultures of AF model. Results. In comparison to the sham cohort, the AF induction rate was markedly increased together with AERP reduction within pacing cohort. Such changes were markedly rescued through sac/val treatment in pacing+sac/val cohort. The proteomic expression profiles of CaMKII and Cav1.2 showed that the CaMKII expression was markedly upregulated, while Cav1.2 expression was downregulated in the pacing cohort. Importantly, these effects were absent in pacing+sac/val cohort. Conclusion. Results of this study show that sac/val treatment regulates the expression of CaMKII/Cav1.2 and could alter this pathway in atrial rapid electrical stimulation models. Therefore, this investigation could contribute to a novel strategy in AF therapeutics in clinical settings.

Mechanisms of Premature Ventricular Complexes Caused by QT Prolongation

QT prolongation, due to lengthening of the action potential duration in the ventricles, is a major risk factor of lethal ventricular arrhythmias. A widely known consequence of QT prolongation is the genesis of early afterdepolarizations (EADs), which are associated with arrhythmias through the generation of premature ventricular complexes (PVCs). However, the vast majority of the EADs observed experimentally in isolated ventricular myocytes are phase-2 EADs, and whether phase-2 EADs are mechanistically linked to PVCs in cardiac tissue remains an unanswered question. In this study, we investigate the genesis of PVCs using computer simulations with eight different ventricular action potential models of various species. Based on our results, we classify PVCs as arising from two distinct mechanisms: repolarization gradient (RG)-induced PVCs and phase-2 EAD-induced PVCs. The RG-induced PVCs are promoted by increasing RG and L-type calcium current and are insensitive to gap junction coupling. EADs are not required for this PVC mechanism. In a paced beat, a single or multiple PVCs can occur depending on the properties of the RG. In contrast, phase-2 EAD-induced PVCs occur only when the RG is small and are suppressed by increasing RG and more sensitive to gap junction coupling. Unlike with RG-induced PVCs, in each paced beat, only a single EAD-induced PVC can occur no matter how many EADs in an action potential. In the wide parameter ranges we explore, RG-induced PVCs can be observed in all models, but the EAD-induced PVCs can only be observed in five of the eight models. The links between these two distinct PVC mechanisms and arrhythmogenesis in animal experiments and clinical settings are discussed.

Canard analysis reveals why a large Ca2+ window current promotes early afterdepolarizations in cardiac myocytes

The pumping of blood through the heart is due to a wave of muscle contractions that are in turn due to a wave of electrical activity initiated at the sinoatrial node. At the cellular level, this wave of electrical activity corresponds to the sequential excitation of electrically coupled cardiac cells. Under some conditions, the normally-long action potentials of cardiac cells are extended even further by small oscillations called early afterdepolarizations (EADs) that can occur either during the plateau phase or repolarizing phase of the action potential. Hence, cellular EADs have been implicated as a driver of potentially lethal cardiac arrhythmias. One of the major determinants of cellular EAD production and repolarization failure is the size of the overlap region between Ca2+ channel activation and inactivation, called the window region. In this article, we interpret the role of the window region in terms of the fast-slow structure of a low-dimensional model for ventricular action potential generation. We demonstrate that the effects of manipulation of the size of the window region can be understood from the point of view of canard theory. We use canard theory to explain why enlarging the size of the window region elicits EADs and why shrinking the window region can eliminate them. We also use the canard mechanism to explain why some manipulations in the size of the window region have a stronger influence on cellular electrical behavior than others. This dynamical viewpoint gives predictive power that is beyond that of the biophysical explanation alone while also uncovering a common mechanism for phenomena observed in experiments on both atrial and ventricular cardiac cells.

Conditions for the genesis of early afterdepolarization in a model of a ventricular myocyte

Early afterdepolarization (EAD) is a major arrhythmogenic factor in the long QT syndrome (LQTS), whose conditions for genesis have puzzled people for several decades. Here, we employ the phase I Luo-Rudy ventricular myocyte model to investigate EAD using methods from nonlinear dynamics and provide valuable insights into EAD genesis from a physical perspective. Two major results are obtained: (i) Sufficient parametric conditions for EAD are analytically determined and then used to analyze in detail the effects of the physiological parameters. (ii) The normal form of the Hopf bifurcation that leads to EAD is derived and then used to determine whether the Hopf bifurcation is subcritical or supercritical for EAD genesis and the corresponding amplitude and period of the EAD oscillation. Our work here paves the way for further studies of more complicated multi-scale dynamics of EAD and may lead to effective treatments for LQTS arrhythmias.

Multiple Dynamical Mechanisms of Phase-2 Early Afterdepolarizations in a Human Ventricular Myocyte Model: Involvement of Spontaneous SR Ca 2+ Release

Early afterdepolarization (EAD) is known to cause lethal ventricular arrhythmias in long QT syndrome (LQTS). In this study, dynamical mechanisms of EAD formation in human ventricular myocytes (HVMs) were investigated using the mathematical model developed by ten Tusscher and Panfilov (Am J Physiol Heart Circ Physiol 291, 2006). We explored how the rapid (IKr) and slow (IKs) components of delayed-rectifier K+ channel currents, L-type Ca2+ channel current (ICa L), Na+/Ca2+ exchanger current (INCX), and intracellular Ca2+ handling via the sarcoplasmic reticulum (SR) contribute to initiation, termination and modulation of phase-2 EADs during pacing in relation to bifurcation phenomena in non-paced model cells. Parameter-dependent dynamical behaviors of the non-paced model cell were determined by calculating stabilities of equilibrium points (EPs) and limit cycles, and bifurcation points to construct bifurcation diagrams. Action potentials (APs) and EADs during pacing were reproduced by numerical simulations for constructing phase diagrams of the paced model cell dynamics. Results are summarized as follows: (1) A modified version of the ten Tusscher-Panfilov model with accelerated ICaL inactivation could reproduce bradycardia-related EADs in LQTS type 2 and β-adrenergic stimulation-induced EADs in LQTS type 1. (2) Two types of EADs with different initiation mechanisms, ICaL reactivation-dependent and spontaneous SR Ca2+ release-mediated EADs, were detected. (3) Termination of EADs (AP repolarization) during pacing depended on the slow activation of IKs. (4) Spontaneous SR Ca2+ releases occurred at higher Ca2+ uptake rates, attributable to the instability of steady-state intracellular Ca2+ concentrations. Dynamical mechanisms of EAD formation and termination in the paced model cell are closely related to stability changes (bifurcations) in dynamical behaviors of the non-paced model cell, but they are model-dependent. Nevertheless, the modified ten Tusscher-Panfilov model would be useful for systematically investigating possible dynamical mechanisms of EAD-related arrhythmias in LQTS.

R-From-T as a Common Mechanism of Arrhythmia Initiation in Long QT Syndromes

BACKGROUND

Long QT syndromes (LQTS) arise from many genetic and nongenetic causes with certain characteristic ECG features preceding polymorphic ventricular tachyarrhythmias (PVTs). However, how the many molecular causes result in these characteristic ECG patterns and how these patterns are mechanistically linked to the spontaneous initiation of PVT remain poorly understood.

METHODS

Anatomic human ventricle and simplified tissue models were used to investigate the mechanisms of spontaneous initiation of PVT in LQTS.

RESULTS

Spontaneous initiation of PVT was elicited by gradually ramping up ICa,L to simulate the initial phase of a sympathetic surge or by changing the heart rate, reproducing the different genotype-dependent clinical ECG features. In LQTS type 2 (LQT2) and LQTS type 3 (LQT3), T-wave alternans was observed followed by premature ventricular complexes (PVCs). Compensatory pauses occurred resulting in short-long-short sequences. As ICa,L increased further, PVT episodes occurred, always preceded by a short-long-short sequence. However, in LQTS type 1 (LQT1), once a PVC occurred, it always immediately led to an episode of PVT. Arrhythmias in LQT2 and LQT3 were bradycardia dependent, whereas those in LQT1 were not. In all 3 genotypes, PVCs always originated spontaneously from the steep repolarization gradient region and manifested on ECG as R-on-T. We call this mechanism R-from-T, to distinguish it from the classic explanation of R-on-T arrhythmogenesis in which an exogenous PVC coincidentally encounters a repolarizing region. In R-from-T, the PVC and the T wave are causally related, where steep repolarization gradients combined with enhanced ICa,L lead to PVCs emerging from the T wave. Since enhanced ICa,L was required for R-from-T to occur, suppressing window ICa,L effectively prevented arrhythmias in all 3 genotypes.

CONCLUSIONS

Despite the complex molecular causes, these results suggest that R-from-T is likely a common mechanism for PVT initiation in LQTS. Targeting ICa,L properties, such as suppressing window ICa,L or preventing excessive ICa,L increase, could be an effective unified therapy for arrhythmia prevention in LQTS.

Determinants of early afterdepolarization properties in ventricular myocyte models

Early afterdepolarizations (EADs) are spontaneous depolarizations during the repolarization phase of an action potential in cardiac myocytes. It is widely known that EADs are promoted by increasing inward currents and/or decreasing outward currents, a condition called reduced repolarization reserve. Recent studies based on bifurcation theories show that EADs are caused by a dual Hopf-homoclinic bifurcation, bringing in further mechanistic insights into the genesis and dynamics of EADs. In this study, we investigated the EAD properties, such as the EAD amplitude, the inter-EAD interval, and the latency of the first EAD, and their major determinants. We first made predictions based on the bifurcation theory and then validated them in physiologically more detailed action potential models. These properties were investigated by varying one parameter at a time or using parameter sets randomly drawn from assigned intervals. The theoretical and simulation results were compared with experimental data from the literature. Our major findings are that the EAD amplitude and takeoff potential exhibit a negative linear correlation; the inter-EAD interval is insensitive to the maximum ionic current conductance but mainly determined by the kinetics of I_Ca,L and the dual Hopf-homoclinic bifurcation; and both inter-EAD interval and latency vary largely from model to model. Most of the model results generally agree with experimental observations in isolated ventricular myocytes. However, a major discrepancy between modeling results and experimental observations is that the inter-EAD intervals observed in experiments are mainly between 200 and 500 ms, irrespective of species, while those of the mathematical models exhibit a much wider range with some models exhibiting inter-EAD intervals less than 100 ms. Our simulations show that the cause of this discrepancy is likely due to the difference in I_Ca,L recovery properties in different mathematical models, which needs to be addressed in future action potential model development.

Early afterdepolarizations (EADs) are abnormal depolarizations during the plateau phase of action potential in cardiac myocytes, arising from a dual Hopf-homoclinic bifurcation. The same bifurcations are also responsible for certain types of bursting behaviors in other cell types, such as beta cells and neuronal cells. EADs are known to play important role in the genesis of lethal arrhythmias and have been widely studied in both experiments and computer models. However, a detailed comparison between the properties of EADs observed in experiments and those from mathematical models have not been carried out. In this study, we performed theoretical analyses and computer simulations of different ventricular action potential models as well as different species to investigate the properties of EADs and compared these properties to those observed in experiments. While the EAD properties in the action potential models capture many of the EAD properties seen in experiments, the inter-EAD intervals in the computer models differ a lot from model to model, and some of them show very large discrepancy with those observed in experiments. This discrepancy needs to be addressed in future cardiac action potential model development.