Early Afterdepolarizations with Growing Amplitudes via Delayed Subcritical Hopf Bifurcations and Unstable Manifolds of Saddle Foci in Cardiac Action Potential Dynamics

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.

Estimating the probability of early afterdepolarizations and predicting arrhythmic risk associated with long QT syndrome type 1 mutations

Early afterdepolarizations (EADs) are action potential (AP) repolarization abnormalities that can trigger lethal arrhythmias. Simulations using biophysically detailed cardiac myocyte models can reveal how model parameters influence the probability of these cellular arrhythmias; however, such analyses can pose a huge computational burden. We have previously developed a highly simplified approach in which logistic regression models (LRMs) map parameters of complex cell models to the probability of ectopic beats. Here, we extend this approach to predict the probability of EADs (P(EAD)) as a mechanistic metric of arrhythmic risk. We use the LRM to investigate how changes in parameters of the slow-activating delayed rectifier current (IKs) affect P(EAD) for 17 different long QT syndrome type 1 (LQTS1) mutations. In this LQTS1 clinical arrhythmic risk prediction task, we compared P(EAD) for these 17 mutations with two other recently published model-based arrhythmia risk metrics (AP morphology metric across populations of myocyte models and transmural repolarization prolongation based on a one-dimensional [1D] tissue-level model). These model-based risk metrics yield similar prediction performance; however, each fails to stratify clinical risk for a significant number of the 17 studied LQTS1 mutations. Nevertheless, an interpretable ensemble model using multivariate linear regression built by combining all of these model-based risk metrics successfully predicts the clinical risk of 17 mutations. These results illustrate the potential of computational approaches in arrhythmia risk prediction.

Fast-slow analysis as a technique for understanding the neuronal response to current ramps

The standard protocol for studying the spiking properties of single neurons is the application of current steps while monitoring the voltage response. Although this is informative, the jump in applied current is artificial. A more physiological input is where the applied current is ramped up, reflecting chemosensory input. Unsurprisingly, neurons can respond differently to the two protocols, since ion channel activation and inactivation are affected differently. Understanding the effects of current ramps, and changes in their slopes, is facilitated by mathematical models. However, techniques for analyzing current ramps are under-developed. In this article, we demonstrate how current ramps can be analyzed in single neuron models. The primary issue is the presence of gating variables that activate on slow time scales and are therefore far from equilibrium throughout the ramp. The use of an appropriate fast-slow analysis technique allows one to fully understand the neural response to ramps of different slopes. This study is motivated by data from olfactory bulb dopamine neurons, where both fast ramp (tens of milliseconds) and slow ramp (tens of seconds) protocols are used to understand the spiking profiles of the cells. The slow ramps generate experimental bifurcation diagrams with the applied current as a bifurcation parameter, thereby establishing asymptotic spiking activity patterns. The faster ramps elicit purely transient behavior that is of relevance to most physiological inputs, which are short in duration. The two protocols together provide a broader understanding of the neuron's spiking profile and the role that slowly activating ion channels can play.

Life and death saddles in the heart

Saddle points are responsible for threshold phenomena of many biological systems. In the heart, saddle points determine the normal excitability and conduction, but are also responsible for certain abnormal action potential behaviors associated with lethal arrhythmias. We investigate the dynamical mechanisms for the genesis of lethal extra heartbeats in heterogeneous cardiac tissue under two diseased conditions. For both conditions, the lethal events occur when the system is close to the saddle point, implying the pivotal role of the saddle point in cardiac arrhythmogenesis.

Mechanisms of phase-3 early afterdepolarizations and triggered activities in ventricular myocyte models

Early afterdepolarizations (EADs) are abnormal depolarizations during the repolarizing phase of the action potential, which are associated with cardiac arrhythmogenesis. EADs are classified into phase-2 and phase-3 EADs. Phase-2 EADs occur during phase 2 of the action potential, with takeoff potentials typically above -40 mV. Phase-3 EADs occur during phase 3 of the action potential, with takeoff potential between -70 and -50 mV. Since the amplitude of phase-3 EADs can be as large as that of a regular action potential, they are also called triggered activities (TAs). This also makes phase-3 EADs and TAs much more arrhythmogenic than phase-2 EADs since they can propagate easily in tissue. Although phase-2 EADs have been widely observed, phase-3 EADs and TAs have been rarely demonstrated in isolated ventricular myocytes. Here we carry out computer simulations of three widely used ventricular action potential models to investigate the mechanisms of phase-3 EADs and TAs. We show that when the T-type Ca2+ current (ICa,T ) is absent (e.g., in normal ventricular myocytes), besides the requirement of increasing inward currents and reducing outward currents as for phase-2 EADs, the occurrence of phase-3 EADs and TAs requires a substantially large increase of the L-type Ca2+ current and the slow component of the delayed rectifier K+ current. The presence of ICa,T (e.g., in neonatal and failing ventricular myocytes) can greatly reduce the thresholds of these two currents for phase-3 EADs and TAs. This implies that ICa,T may play an important role in arrhythmogenesis in cardiac diseases.

Deterministic and Stochastic Cellular Mechanisms Contributing to Carbon Monoxide Induced Ventricular Arrhythmias

Chronic exposure to low levels of Carbon Monoxide is associated with an increased risk of cardiac arrhythmia. Microelectrode recordings from rat and guinea pig single isolated ventricular myocytes exposed to CO releasing molecule CORM-2 and excited at 0.2/s show repolarisation changes that develop over hundreds of seconds: action potential prolongation by delayed repolarisation, EADs, multiple EADs and oscillations around the plateau, leading to irreversible repolarisation failure. The measured direct effects of CO on currents in these cells, and ion channels expressed in mammalian systems showed an increase in prolonged late Na+, and a decrease in the maximal T- and L-type Ca++. peak and late Na+, ultra-rapid delayed, delayed rectifier, and the inward rectifier K+ currents. Incorporation of these CO induced changes in maximal currents in ventricular cell models; (Gattoni et al., J. Physiol., 2016, 594, 4193-4224) (rat) and (Luo and Rudy, Circ. Res., 1994, 74, 1071-1096) (guinea-pig) and human endo-, mid-myo- and epi-cardial (O'Hara et al., PLoS Comput. Biol., 2011, 7, e1002061) models, by changes in maximal ionic conductance reproduces these repolarisation abnormalities. Simulations of cell populations with Gaussian distributions of maximal conductance parameters predict a CO induced increase in APD and its variability. Incorporation of these predicted CO induced conductance changes in human ventricular cell electrophysiology into ventricular tissue and wall models give changes in indices for the probability of the initiation of re-entrant arrhythmia.

Circadian Rhythms of Early Afterdepolarizations and Ventricular Arrhythmias in a Cardiomyocyte Model

Sudden cardiac arrest is a malfunction of the heart's electrical system, typically caused by ventricular arrhythmias, that can lead to sudden cardiac death (SCD) within minutes. Epidemiological studies have shown that SCD and ventricular arrhythmias are more likely to occur in the morning than in the evening, and laboratory studies indicate that these daily rhythms in adverse cardiovascular events are at least partially under the control of the endogenous circadian timekeeping system. However, the biophysical mechanisms linking molecular circadian clocks to cardiac arrhythmogenesis are not fully understood. Recent experiments have shown that L-type calcium channels exhibit circadian rhythms in both expression and function in guinea pig ventricular cardiomyocytes. We developed an electrophysiological model of these cells to simulate the effect of circadian variation in L-type calcium conductance. In our simulations, we found that there is a circadian pattern in the occurrence of early afterdepolarizations (EADs), which are abnormal depolarizations during the repolarization phase of a cardiac action potential that can trigger fatal ventricular arrhythmias. Specifically, the model produces EADs in the morning, but not at other times of day. We show that the model exhibits a codimension-2 Takens-Bogdanov bifurcation that serves as an organizing center for different types of EAD dynamics. We also simulated a two-dimensional spatial version of this model across a circadian cycle. We found that there is a circadian pattern in the breakup of spiral waves, which represents ventricular fibrillation in cardiac tissue. Specifically, the model produces spiral wave breakup in the morning, but not in the evening. Our computational study is the first, to our knowledge, to propose a link between circadian rhythms and EAD formation and suggests that the efficacy of drugs targeting EAD-mediated arrhythmias may depend on the time of day that they are administered.

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.

Bifurcations and Slow-Fast Analysis in a Cardiac Cell Model for Investigation of Early Afterdepolarizations

In this study, we teased out the dynamical mechanisms underlying the generation of arrhythmogenic early afterdepolarizations (EADs) in a three-variable model of a mammalian ventricular cell. Based on recently published studies, we consider a 1-fast, 2-slow variable decomposition of the system describing the cellular action potential. We use sweeping techniques, such as the spike-counting method, and bifurcation and continuation methods to identify parametric regions with EADs. We show the existence of isolas of periodic orbits organizing the different EAD patterns and we provide a preliminary classification of our fast–slow decomposition according to the involved dynamical phenomena. This investigation represents a basis for further studies into the organization of EAD patterns in the parameter space and the involved bifurcations.

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.

Why pacing frequency affects the production of early afterdepolarizations in cardiomyocytes: An explanation revealed by slow-fast analysis of a minimal model

Early afterdepolarizations (EADs) are pathological voltage oscillations in cardiomyocytes that have been observed in response to a number of pharmacological agents and disease conditions. Phase-2 EADs consist of small voltage fluctuations during the plateau of an action potential, typically under conditions in which the action potential is elongated. Although a single-cell behavior, EADs can lead to tissue-level arrhythmias. Much is currently known about the biophysical mechanisms (i.e., the roles of ion channels and intracellular Ca^{2+} stores) for the various forms of EADs, due partially to the development and analysis of mathematical models. This includes the application of slow-fast analysis, which takes advantage of timescale separation inherent in the system to simplify its analysis. We take this further, using a minimal three-dimensional model to demonstrate that phase-2 EADs are canards formed in the neighborhood of a folded node singularity. This allows us to predict the number of EADs that can be produced for a given parameter set, and provides guidance on parameter changes that facilitate or inhibit EAD production. With this approach, we demonstrate why periodic stimulation, as occurs in intact heart, preferentially facilitates EAD production when applied at low frequencies. We also explain the origin of complex alternan dynamics that can occur with intermediate-frequency stimulation, in which varying numbers of EADs are produced with each pulse. These revelations fall out naturally from an understanding of folded node singularities, but are difficult to glean from knowledge of the biophysical mechanism for EADs alone. Therefore, understanding the canard mechanism is a useful complement to understanding of the biophysical mechanism that has been developed over years of experimental and computational investigations.

Early afterdepolarizations in cardiac action potentials as mixed mode oscillations due to a folded node singularity

Early afterdepolarizations (EADs) are pathological voltage oscillations during the repolarization phase of cardiac action potentials. They are considered as potential precursors to cardiac arrhythmias and have recently gained much attention in the context of preclinical drug safety testing under the Comprehensive in vitro Proarrhythmia Assay (CiPA) paradigm. From the viewpoint of multiple time scales theory, the onset of EADs has previously been studied by means of mathematical action potential models with one slow ion channel gating variable. In this article, we for the first time associate EADs with mixed mode oscillations in dynamical systems with two slow gating variables and present a folded node singularity of the slow flow as a novel mechanism for EADs genesis. We derive regions of the pharmacology parameter space in which EADs occur using both the folded node analysis and a full system bifurcation analysis, and we suggest the normal distance to the boundary of the EADs region as a mechanism-based risk metric to computationally estimate a drug’s proarrhythmic liability.

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.

Identification of Ion Currents Components Generating Field Potential Recorded in MEA From hiPSC-CM

OBJECTIVE

Multi electrodes arrays (MEAs) combined with cardiomyocytes derived from human-induced pluripotent stem cells (hiPSC-CMs) can enable high- or medium-throughput drug screening in safety pharmacology. This technology has recently attracted a lot of attention, in particular from an international initiative named CiPA. But it is currently limited by the difficulty to analyze the measured signals. We propose a strategy to analyze the signals acquired by the MEA and to automatically deduce the channels affected by the drug.

METHODS

Our method is based on the bidomain equations, a model for the MEA electrodes, and an inverse problem strategy.

RESULTS

in silico MEA signals are obtained for two commercial devices and an example of early after depolarization is presented. Then, by processing real signals obtained for four different compounds, our algorithm was able to provide dose-response curves for potassium, sodium, and calcium channels. For ivabradine and moxifloxacin, the IC50 and dose-response curves are in very good agreement with known values.

SIGNIFICANCE

The proposed strategy offers a possible answer to a major question raised by the community of safety pharmacology. By allowing a more automated analysis of the signals, our approach could contribute to promote the technology based on MEA and hiPSC-CMs and, therefore, improve reliability and efficiency of drug screening.

Period doubling cascades of limit cycles in cardiac action potential models as precursors to chaotic early Afterdepolarizations

Background

Early afterdepolarizations (EADs) are pathological voltage oscillations during the repolarization phase of cardiac action potentials (APs). EADs are caused by drugs, oxidative stress or ion channel disease, and they are considered as potential precursors to cardiac arrhythmias in recent attempts to redefine the cardiac drug safety paradigm. The irregular behaviour of EADs observed in experiments has been previously attributed to chaotic EAD dynamics under periodic pacing, made possible by a homoclinic bifurcation in the fast subsystem of the deterministic AP system of differential equations.

Results

In this article we demonstrate that a homoclinic bifurcation in the fast subsystem of the action potential model is neither a necessary nor a sufficient condition for the genesis of chaotic EADs. We rather argue that a cascade of period doubling (PD) bifurcations of limit cycles in the full AP system paves the way to chaotic EAD dynamics across a variety of models including a) periodically paced and spontaneously active cardiomyocytes, b) periodically paced and non-active cardiomyocytes as well as c) unpaced and spontaneously active cardiomyocytes. Furthermore, our bifurcation analysis reveals that chaotic EAD dynamics may coexist in a stable manner with fully regular AP dynamics, where only the initial conditions decide which type of dynamics is displayed.

Conclusions

EADs are a potential source of cardiac arrhythmias and hence are of relevance both from the viewpoint of drug cardiotoxicity testing and the treatment of cardiomyopathies. The model-independent association of chaotic EADs with period doubling cascades of limit cycles introduced in this article opens novel opportunities to study chaotic EADs by means of bifurcation control theory and inverse bifurcation analysis. Furthermore, our results may shed new light on the synchronization and propagation of chaotic EADs in homogeneous and heterogeneous multicellular and cardiac tissue preparations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0422-4) contains supplementary material, which is available to authorized users.

Phase diagrams and dynamics of a computationally efficient map-based neuron model

We introduce a new map-based neuron model derived from the dynamical perceptron family that has the best compromise between computational efficiency, analytical tractability, reduced parameter space and many dynamical behaviors. We calculate bifurcation and phase diagrams analytically and computationally that underpins a rich repertoire of autonomous and excitable dynamical behaviors. We report the existence of a new regime of cardiac spikes corresponding to nonchaotic aperiodic behavior. We compare the features of our model to standard neuron models currently available in the literature.