Determinants of early afterdepolarization properties in ventricular myocyte models

Mechanism of Arrhythmogenesis Driven by Early After Depolarizations in Cardiac Tissue

Early-after depolarizations (EADs) are changes in the action potential plateau that can lead to cardiac arrhythmia. At the cellular level, these oscillations are irregular and change from beat to beat due to the sensitivity of voltage repolarization to subcellular stochastic processes. However, the behavior of EADs in tissue, where cells are strongly coupled by gap junctions, is less understood. In this study, we develop a computational model of EADs caused by a reduction in the rate of calcium-induced inactivation of the L-type calcium channel. We find that, as inactivation decreases EADs occur with durations varying randomly from beat to beat. In cardiac tissue, however, gap junction coupling between cells dampens these fluctuations, and it is unclear what dictates the formation of EADs. In this study we show that EADs in cardiac tissue can be modeled by the deterministic limit of a stochastic single-cell model. Analysis of this deterministic model reveals that EADs emerge in tissue after an abrupt transition to alternans, where large populations of cells suddenly synchronize, causing EADs on every other beat. We analyze this transition and show that it is due to a discontinuous bifurcation that leads to a large change in the action potential duration in response to very small changes in pacing rate. We further demonstrate that this transition is highly arrhythmogenic, as the sudden onset of EADs in cardiac tissue promotes conduction block and reentry. Our results highlight the importance of EAD alternans in arrhythmogenesis and suggests that ectopic beats are not required.

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Curvature-mediated source and sink effects on the genesis of premature ventricular complexes in long QT syndrome

Premature ventricular complexes (PVCs) are spontaneous excitations occurring in the ventricles of the heart that are associated with ventricular arrhythmias and sudden cardiac death. Under long QT conditions, PVCs can be mediated by repolarization gradient (RG) and early afterdepolarizations (EADs), yet the effects of heterogeneities or geometry of the RG or EAD regions on PVC genesis remain incompletely understood. In this study, we use computer simulation to systematically investigate the effects of the curvature of the RG border region on PVC genesis under long QT conditions. We show that PVCs can be either promoted or suppressed by negative or positive RG border curvature depending on the source and sink conditions. When the origin of oscillation is in the source region and the source is too strong, a positive RG border curvature can promote PVCs by causing the source area to oscillate. When the origin of oscillation is in the sink region, a negative RG border curvature can promote PVCs by causing the sink area to oscillate. Furthermore, EAD-mediated PVCs are also promoted by negative border curvature. We also investigate the effects of wavefront curvature and show that PVCs are promoted by convex but suppressed by concave wavefronts; however, the effect of wavefront curvature is much smaller than that of RG border curvature. In conclusion, besides the increase of RG and occurrence of EADs caused by QT prolongation, the geometry of the RG border plays important roles in PVC genesis, which can greatly increase the risk of arrhythmias in cardiac diseases.NEW & NOTEWORTHY The effects of the curvature or geometry of the repolarization gradient region and wavefront curvature on the genesis of premature ventricular complexes are systematically investigated using computer modeling and simulation. Premature ventricular complexes can be promoted by either positive or negative curvature of the gradient region depending on the source and sink conditions. The underlying mechanisms of the curvature effects are revealed, which provides mechanistic insights into arrhythmogenesis in cardiac diseases.

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.

Obesity and atrial fibrillation: a narrative review from arrhythmogenic mechanisms to clinical significance

The prevalence of obesity and atrial fibrillation (AF), which are inextricably linked, is rapidly increasing worldwide. Obesity rates are higher among patients with AF than healthy individuals. Some epidemiological data indicated that obese patients were more likely to develop AF, but others reported no significant correlation. Obesity-related hypertension, diabetes, and obstructive sleep apnea are all associated with AF. Additionally, increased epicardial fat, systemic inflammation, and oxidative stress caused by obesity can induce atrial enlargement, inflammatory activation, local myocardial fibrosis, and electrical conduction abnormalities, all of which led to AF and promoted its persistence. Weight loss reduced the risk and reversed natural progression of AF, which may be due to its anti-fibrosis and inflammation effect. However, fluctuations in weight offset the benefits of weight loss. Therefore, the importance of steady weight loss urges clinicians to incorporate weight management interventions in the treatment of patients with AF. In this review, we discuss the epidemiology of obesity and AF, summarize the mechanisms by which obesity triggers AF, and explain how weight loss improves the prognosis of AF.

Bifurcations to transient and oscillatory excitations in inhomogeneous excitable media: Insights into arrhythmogenesis in long QT syndrome

Ventricular arrhythmias are the leading cause of sudden cardiac death. Understanding the mechanisms of arrhythmia initiation is important for developing effective therapeutics for prevention. Arrhythmias can be induced via premature external stimuli or occur spontaneously via dynamical instabilities. Computer simulations have shown that a large repolarization gradient due to regional prolongation of the action potential duration can result in instabilities leading to premature excitations and arrhythmias, but the bifurcation remains to be elucidated. In this study we carry out numerical simulations and linear stability analyses using a one-dimensional heterogeneous cable consisting of the FitzHugh-Nagumo model. We show that a Hopf bifurcation leads to local oscillations, which, once their amplitudes are large enough, lead to spontaneous propagating excitations. Depending on the degree of heterogeneities, these excitations can range from one to many and to be sustained oscillations, manifesting as premature ventricular contractions (PVCs) and sustained arrhythmias. The dynamics depends on the repolarization gradient and the length of the cable. Complex dynamics is also induced by the repolarization gradient. The mechanistic insights from the simple model may help in the understanding of the genesis of PVCs and arrhythmias in long QT syndrome.

Spatiotemporal patterns of early afterdepolarizations underlying abnormal T-wave morphologies in a tissue model of the Purkinje-ventricular system

Sudden cardiac death (SCD) is a leading cause of death worldwide, and the majority of SCDs are caused by acute ventricular arrhythmias (VAs). Early afterdepolarizations (EADs) are an important trigger of VA under pathological conditions, e.g., inherited or acquired long QT syndrome (LQTS). However, it remains unclear how EAD events at the cellular level are spatially organized at the tissue level to induce and maintain ventricular arrhythmias and whether the spatial-temporal patterns of EADs at the tissue level are associated with abnormal T-wave morphologies that are often observed in LQTS, such as broad-based, notched or bifid; late appearance; and pointed T-waves. Here, a tissue model of the Purkinje-ventricular system (PVS) was developed to quantitatively investigate the complex spatial-temporal dynamics of EADs during T-wave abnormalities. We found that (1) while major inhibition of ICaL can substantially reduce the excitability of the PVS leading to conduction failures, moderate ICaL inhibition can promote occurrences of AP alternans at short cycle lengths (CLs), and EAD events preferentially occur with a major reduction of IKr (>50%) at long CLs; (2) with a minor reduction of ICaL, spatially synchronized steady-state EAD events with inverted and biphasic T-waves can be "weakened" into beat-to-beat concurrences of spatially synchronized EADs and T-wave alternans, and as pacing CLs increase, beat-to-beat concurrences of localized EADs with late-appearing and pointed T-wave morphologies can be observed; (3) under certain conditions, localized EAD events in the midmyocardium may trigger slow uni-directional electric propagation with inverted (antegrade) or upright (retrograde) broad-based T-waves; (4) spatially discordant EADs were typically characterized by desynchronized spontaneous onsets of EAD events between two groups of PVS tissues with biphasic T-wave morphologies, and they can evolve into spatially discordant oscillating EAD patterns with sustained or self-terminated alternating EAD and electrocardiogram (ECG) patterns. Our results provide new insights into the spatiotemporal aspects of the onset and development of EADs and suggest possible mechanistic links between the complex spatial dynamics of EADs and T-wave morphologies.

Dynamical mechanism for generation of arrhythmogenic early afterdepolarizations in cardiac myocytes: Insights from in silico electrophysiological models

We analyze the dynamical mechanisms underlying the formation of arrhythmogenic early afterdepolarizations (EADs) in two mathematical models of cardiac cellular electrophysiology: the Sato et al. biophysically detailed model of a rabbit ventricular myocyte of dimension 27 and a reduced version of the Luo-Rudy mammalian myocyte model of dimension 3. Based on a comparison of the two models, with detailed bifurcation analysis using spike-counting techniques and continuation methods in the simple model and numerical explorations in the complex model, we locate the point where the first EAD originates in an unstable branch of periodic orbits. These results serve as a basis to propose a conjectured scheme involving a hysteresis mechanism with the creation of alternans and EADs in the unstable branch. This theoretical scheme fits well with electrophysiological experimental data on EAD generation and hysteresis phenomena. Our findings open the door to the development of novel methods for pro-arrhythmia risk prediction related to EAD generation without actual induction of EADs.

The Impact of Mental Stress on Cardiovascular Health—Part II

Endothelial dysfunction is one of the earliest manifestations of atherosclerosis, contributing to its development and progression. Mental stress induces endothelial dysfunction through increased activity of the sympathetic nervous system, release of corticotropin-releasing hormone from the hypothalamus, inhibition of nitric oxide (NO) synthesis by cortisol, and increased levels of pro-inflammatory cytokines. Mental-stress-induced increased output of the sympathetic nervous system and concomitant withdrawal of the parasympathetic inflammatory reflex results in systemic inflammation and activation of a neural–hematopoietic–arterial axis. This includes the brainstem and subcortical regions network, bone marrow activation, release of leukocytes into the circulation and their migration to the arterial wall and atherosclerotic plaques. Low-grade, sterile inflammation is involved in all steps of atherogenesis, from coronary plaque formation to destabilisation and rupture. Increased sympathetic tone may cause arterial smooth-muscle-cell proliferation, resulting in vascular hypertrophy, thus contributing to the development of hypertension. Emotional events also cause instability of cardiac repolarisation due to brain lateralised imbalance of cardiac autonomic nervous stimulation, which may lead to asymmetric repolarisation and arrhythmia. Acute emotional stress can also provoke severe catecholamine release, leading to direct myocyte injury due to calcium overload, known as myocytolysis, coronary microvascular vasoconstriction, and an increase in left ventricular afterload. These changes can trigger a heart failure syndrome mimicking acute myocardial infarction, characterised by transient left ventricular dysfunction and apical ballooning, known as stress (Takotsubo) cardiomyopathy. Women are more prone than men to develop mental-stress-induced myocardial ischemia (MSIMI), probably reflecting gender differences in brain activation patterns during mental stress. Although guidelines on CV prevention recognise psychosocial factors as risk modifiers to improve risk prediction and decision making, the evidence that their assessment and treatment will prevent CAD needs further evaluation.

Rate dependence of early afterdepolarizations in the His-Purkinje system: A simulation study

BACKGROUND AND OBJECTIVE

Early afterdepolarizations (EADs) are associated with a variety of arrhythmias and have the property of rate dependence. EADs can occur in Purkinje cells while the effect of rate dependence of EADs in the His-Purkinje system has not been fully investigated. In order to reveal the rate dependence of EADs in the His-Purkinje system and its effect on ventricular electrical activities, the simulation research was carried out in this manuscript.

METHODS

This manuscript first studied the relationship between the occurrence of EADs and stimulation cycle length on the DiFranNoble cell model. Then, the relationship between the rate dependence of EADs and the conduction block of the His-Purkinje system at slow heart rates was studied on the rabbit whole ventricular model including the His-Purkinje system, and its mechanism was analyzed from multiple angles.

RESULTS

① The rate dependence of EADs is related to the inconsistency of EADs occurrence in the His-Purkinje system. When the stimulation cycle length is long or short enough, EADs either occur or not occur stably in the His-Purkinje system, while in a certain stimulation cycle length window, the chaotic state of EADs will be observed. ② The key subcellular factors x-gate is an important mechanism involved to the rate dependence of EADs in the His-Purkinje system. ③ The discrete distribution of x-gate values and the "source-sink" mechanism lead to the inconsistency of EADs in the His-Purkinje system. The prolonged action potential duration caused by EADs can lead to conduction block at slow heart rates.

CONCLUSION

The rate dependence of EADs in Purkinje system can lead to disordered ventricular electrical activity.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

Quantitative roles of ion channel dynamics on ventricular action potential

Mathematical models for the action potential (AP) generation of the electrically excitable cells including the heart are involved different mechanisms including the voltage-dependent currents with nonlinear time- and voltage-gating properties. From the shape of the AP waveforms to the duration of the refractory periods or heart rhythms are greatly affected by the functions describing the features or the quantities of these ion channels. In this work, a mathematical measure to analyze the regional contributions of voltage-gated channels is defined by dividing the AP into phases, epochs, and intervals of interest. The contribution of each time-dependent current for the newly defined cardiomyocyte model is successfully calculated and it is found that the contribution of dominant ion channels changes substantially not only for each phase but also for different regions of the cardiac AP. Besides, the defined method can also be applied in all Hodgkin-Huxley types of electrically excitable cell models to be able to understand the underlying dynamics better.

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.

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.

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.

KEY FACTORS TO THE RATE DEPENDENCE OF EARLY AFTERDEPOLARIZATIONS: A SIMULATION STUDY

Early afterdepolarizations (EADs) in cardiac myocytes have been reported to be associated with a series of cardiac arrhythmias. The generation of EADs has a high correlation with the excitation period, leading to repolarization dispersion over the cardiac tissue. However, the mechanism of EAD rate dependence has not been thoroughly revealed. In this study, the simulation approach was used to investigate the mechanism underlying EAD rate dependence. The results indicated that the gating variable of the delayed rectifier potassium current ([Formula: see text]-gate) and the intracellular sodium ion concentration ([Na[Formula: see text]][Formula: see text] were key factors contributing to EAD rate dependence. Also, different mathematical models showed different types of EAD rate dependence, which needs to be considered in the future simulation research related to EADs.

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.

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.

Targeting Ca2 + Handling Proteins for the Treatment of Heart Failure and Arrhythmias

Diseases of the heart, such as heart failure and cardiac arrhythmias, are a growing socio-economic burden. Calcium (Ca²⁺) dysregulation is key hallmark of the failing myocardium and has long been touted as a potential therapeutic target in the treatment of a variety of cardiovascular diseases (CVD). In the heart, Ca²⁺ is essential for maintaining normal cardiac function through the generation of the cardiac action potential and its involvement in excitation contraction coupling. As such, the proteins which regulate Ca²⁺ cycling and signaling play a vital role in maintaining Ca²⁺ homeostasis. Changes to the expression levels and function of Ca²⁺-channels, pumps and associated intracellular handling proteins contribute to altered Ca²⁺ homeostasis in CVD. The remodeling of Ca²⁺-handling proteins therefore results in impaired Ca²⁺ cycling, Ca²⁺ leak from the sarcoplasmic reticulum and reduced Ca²⁺ clearance, all of which contributes to increased intracellular Ca²⁺. Currently, approved treatments for targeting Ca²⁺ handling dysfunction in CVD are focused on Ca²⁺ channel blockers. However, whilst Ca²⁺ channel blockers have been successful in the treatment of some arrhythmic disorders, they are not universally prescribed to heart failure patients owing to their ability to depress cardiac function. Despite the progress in CVD treatments, there remains a clear need for novel therapeutic approaches which are able to reverse pathophysiology associated with heart failure and arrhythmias. Given that heart failure and cardiac arrhythmias are closely associated with altered Ca²⁺ homeostasis, this review will address the molecular changes to proteins associated with both Ca²⁺-handling and -signaling; their potential as novel therapeutic targets will be discussed in the context of pre-clinical and, where available, clinical data.

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.