Dependence of phase-2 reentry and repolarization dispersion on epicardial and transmural ionic heterogeneity: a simulation study

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.

Inherited Arrhythmias in the Pediatric Population: An Updated Overview

Pediatric cardiomyopathies (CMs) and electrical diseases constitute a heterogeneous spectrum of disorders distinguished by structural and electrical abnormalities in the heart muscle, attributed to a genetic variant. They rank among the main causes of morbidity and mortality in the pediatric population, with an annual incidence of 1.1-1.5 per 100,000 in children under the age of 18. The most common conditions are dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). Despite great enthusiasm for research in this field, studies in this population are still limited, and the management and treatment often follow adult recommendations, which have significantly more data on treatment benefits. Although adult and pediatric cardiac diseases share similar morphological and clinical manifestations, their outcomes significantly differ. This review summarizes the latest evidence on genetics, clinical characteristics, management, and updated outcomes of primary pediatric CMs and electrical diseases, including DCM, HCM, arrhythmogenic right ventricular cardiomyopathy (ARVC), Brugada syndrome (BrS), catecholaminergic polymorphic ventricular tachycardia (CPVT), long QT syndrome (LQTS), and short QT syndrome (SQTS).

Phase 2 Re-Entry Without Ito: Role of Sodium Channel Kinetics in Brugada Syndrome Arrhythmias

BACKGROUND

In Brugada syndrome (BrS), phase 2 re-excitation/re-entry (P2R) induced by the transient outward potassium current (Ito) is a proposed arrhythmia mechanism; yet, the most common genetic defects are loss-of-function sodium channel mutations.

OBJECTIVES

The authors used computer simulations to investigate how sodium channel dysfunction affects P2R-mediated arrhythmogenesis in the presence and absence of Ito.

METHODS

Computer simulations were carried out in 1-dimensional cables and 2-dimensional tissue using guinea pig and human ventricular action potential models.

RESULTS

In the presence of Ito sufficient to generate robust P2R, reducing sodium current (INa) peak amplitude alone only slightly potentiated P2R. When INa inactivation kinetics were also altered to simulate reported effects of BrS mutations and sodium channel blockers, however, P2R occurred even in the absence of Ito. These effects could be potentiated by delaying L-type calcium channel activation or increasing ATP-sensitive potassium current, consistent with experimental and clinical findings. INa-mediated P2R also accounted for sex-related, day and night-related, and fever-related differences in arrhythmia risk in BrS patients.

CONCLUSIONS

Altered INa kinetics synergize powerfully with reduced INa amplitude to promote P2R-induced arrhythmias in BrS in the absence of Ito, establishing a robust mechanistic link between altered INa kinetics and the P2R-mediated arrhythmia mechanism.

The importance of mechanical conditions in the testing of excitation abnormalities in a population of electro-mechanical models of human ventricular cardiomyocytes

Background: Populations of in silico electrophysiological models of human cardiomyocytes represent natural variability in cell activity and are thoroughly calibrated and validated using experimental data from the human heart. The models have been shown to predict the effects of drugs and their pro-arrhythmic risks. However, excitation and contraction are known to be tightly coupled in the myocardium, with mechanical loads and stretching affecting both mechanics and excitation through mechanisms of mechano-calcium-electrical feedback. However, these couplings are not currently a focus of populations of cell models. Aim: We investigated the role of cardiomyocyte mechanical activity under different mechanical conditions in the generation, calibration, and validation of a population of electro-mechanical models of human cardiomyocytes. Methods: To generate a population, we assumed 11 input parameters of ionic currents and calcium dynamics in our recently developed TP + M model as varying within a wide range. A History matching algorithm was used to generate a non-implausible parameter space by calibrating the action potential and calcium transient biomarkers against experimental data and rejecting models with excitation abnormalities. The population was further calibrated using experimental data on human myocardial force characteristics and mechanical tests involving variations in preload and afterload. Models that passed the mechanical tests were validated with additional experimental data, including the effects of drugs with high or low pro-arrhythmic risk. Results: More than 10% of the models calibrated on electrophysiological data failed mechanical tests and were rejected from the population due to excitation abnormalities at reduced preload or afterload for cell contraction. The final population of accepted models yielded action potential, calcium transient, and force/shortening outputs consistent with experimental data. In agreement with experimental and clinical data, the models demonstrated a high frequency of excitation abnormalities in simulations of Dofetilide action on the ionic currents, in contrast to Verapamil. However, Verapamil showed a high frequency of failed contractions at high concentrations. Conclusion: Our results highlight the importance of considering mechanoelectric coupling in silico cardiomyocyte models. Mechanical tests allow a more thorough assessment of the effects of interventions on cardiac function, including drug testing.

Brugada Syndrome: From Molecular Mechanisms and Genetics to Risk Stratification

Brugada syndrome (BrS) is a rare hereditary arrhythmia disorder, with a distinctive ECG pattern, correlated with an increased risk of ventricular arrhythmias and sudden cardiac death (SCD) in young adults. BrS is a complex entity in terms of mechanisms, genetics, diagnosis, arrhythmia risk stratification, and management. The main electrophysiological mechanism of BrS requires further research, with prevailing theories centered on aberrant repolarization, depolarization, and current-load match. Computational modelling, pre-clinical, and clinical research show that BrS molecular anomalies result in excitation wavelength (k) modifications, which eventually increase the risk of arrhythmia. Although a mutation in the SCN5A (Sodium Voltage-Gated Channel Alpha Subunit 5) gene was first reported almost two decades ago, BrS is still currently regarded as a Mendelian condition inherited in an autosomal dominant manner with incomplete penetrance, despite the recent developments in the field of genetics and the latest hypothesis of additional inheritance pathways proposing a more complex mode of inheritance. In spite of the extensive use of the next-generation sequencing (NGS) technique with high coverage, genetics remains unexplained in a number of clinically confirmed cases. Except for the SCN5A which encodes the cardiac sodium channel NaV1.5, susceptibility genes remain mostly unidentified. The predominance of cardiac transcription factor loci suggests that transcriptional regulation is essential to the Brugada syndrome's pathogenesis. It appears that BrS is a multifactorial disease, which is influenced by several loci, each of which is affected by the environment. The primary challenge in individuals with a BrS type 1 ECG is to identify those who are at risk for sudden death, researchers propose the use of a multiparametric clinical and instrumental strategy for risk stratification. The aim of this review is to summarize the latest findings addressing the genetic architecture of BrS and to provide novel perspectives into its molecular underpinnings and novel models of risk stratification.

R-on-T and the initiation of reentry revisited: Integrating old and new concepts

Initiation of reentry requires 2 factors: (1) a triggering event, most commonly focal excitations such as premature ventricular complexes (PVCs); and (2) a vulnerable substrate with regional dispersion of refractoriness and/or excitability, such as occurs during the T wave of the electrocardiogram when some areas of the ventricle have repolarized and recovered excitability but others have not. When the R wave of a PVC coincides in time with the T wave of the previous beat, this timing can lead to unidirectional block and initiation of reentry, known as the R-on-T phenomenon. Classically, the PVC triggering reentry has been viewed as arising focally from 1 region and propagating into another region whose recovery is delayed, resulting in unidirectional conduction block and reentry initiation. However, more recent evidence indicates that PVCs also can arise from the T wave itself. In the latter case, the PVC initiating reentry is not a separate event from the T wave but rather is causally generated from the repolarization gradient that manifests as the T wave. We call the former an "R-to-T" mechanism and the latter an "R-from-T" mechanism, which are initiation mechanisms distinct from each other. Both are important components of the R-on-T phenomenon and need to be taken into account when designing antiarrhythmic strategies. Strategies targeting suppression of triggers alone or vulnerable substrate alone may be appropriate in some instances but not in others. Preventing R-from-T arrhythmias requires suppressing the underlying dynamic tissue instabilities responsible for producing both triggers and substrate vulnerability simultaneously. The same principles are likely to apply to supraventricular arrhythmias.

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.

Absence of Rgs5 Influences the Spatial and Temporal Fluctuation of Cardiac Repolarization in Mice

Aims

This study investigated the contribution of the regulator of G-protein signaling 5 (Rgs5) knockout to the alteration of the action potential duration (APD) restitution and repolarizing dispersion in ventricle.

Methods and Results

The effects of Rgs5^–/– were investigated by QT variance (QTv) and heart rate variability analysis of Rgs5^–/– mice. Monophasic action potential analysis was investigated in isolated Rgs5^–/– heart. Rgs5^–/– did not promote ventricular remodeling. The 24-h QTv and QT variability index (QTVI) of the Rgs5^–/– mice were higher than those of wild-type (WT) mice (P < 0.01). In WT mice, a positive correlation was found between QTv and the standard deviation of all NN intervals (r = 0.62; P < 0.01), but not in Rgs5^–/– mice (R = 0.01; P > 0.05). The absence of Rgs5 resulted in a significant prolongation of effective refractory period and APD in isolated ventricle. In addition, compared with WT mice, the knockout of Rgs5 significantly deepened the slope of the APD recovery curve at all 10 sites of the heart (P < 0.01) and increased the spatial dispersions of S_max (COV-S_max) (WT: 0.28 ± 0.03, Rgs5^–/–: 0.53 ± 0.08, P < 0.01). Compared with WT heart, Rgs5^–/– increased the induced S1–S2 interval at all sites of heart and widened the window of vulnerability of ventricular tachyarrhythmia (P < 0.05).

Conclusion

Our findings indicate that Rgs5^–/– is an important regulator of ventricular tachyarrhythmia in mice by prolonging ventricular repolarization and increasing spatial dispersion in ventricle.

The transient outward potassium current plays a key role in spiral wave breakup in ventricular tissue

Spiral wave reentry as a mechanism of lethal ventricular arrhythmias has been widely demonstrated in animal experiments and recordings from human hearts. It has been shown that in structurally normal hearts spiral waves are unstable, breaking up into multiple wavelets via dynamical instabilities. However, many of the second-generation action potential models give rise only to stable spiral waves, raising issues regarding the underlying mechanisms of spiral wave breakup. In this study, we carried out computer simulations of two-dimensional homogeneous tissues using five ventricular action potential models. We show that the transient outward potassium current (Ito), although it is not required, plays a key role in promoting spiral wave breakup in all five models. As the maximum conductance of Ito increases, it first promotes spiral wave breakup and then stabilizes the spiral waves. In the absence of Ito, speeding up the L-type calcium kinetics can prevent spiral wave breakup, however, with the same speedup kinetics, spiral wave breakup can be promoted by increasing Ito. Increasing Ito promotes single-cell dynamical instabilities, including action potential duration alternans and chaos, and increasing Ito further suppresses these action potential dynamics. These cellular properties agree with the observation that increasing Ito first promotes spiral wave breakup and then stabilizes spiral waves in tissue. Implications of our observations to spiral wave dynamics in the real hearts and action potential model improvements are discussed.NEW & NOTEWORTHY Spiral wave breakup manifesting as multiple wavelets is a mechanism of ventricular fibrillation. It has been known that spiral wave breakup in cardiac tissue can be caused by a steeply sloped action potential duration restitution curve, a property mainly determined by the recovery of L-type calcium current. Here, we show that the transient outward potassium current (Ito) is another current that plays a key role in spiral wave breakup, that is, spiral waves can be stable for low and high maximum Ito conductance but breakup occurs for intermediate maximum Ito conductance. Since Ito is present in normal hearts of many species and required for Brugada syndrome, it may play an important role in the spiral wave stability and arrhythmogenesis under both normal condition and Brugada syndrome.

Atrial Fibrillation in Inherited Channelopathies

Atrial fibrillation (AF), the common sustained arrhythmia in clinical practice, has major public health implications due to its associated morbidity and increased mortality. The AF epidemic is due to the burgeoning elderly population and the identification of novel risk factors, for example, genetics. Since the diagnosis of AF has a major impact on the clinical assessment and management of patients with inherited arrhythmia syndromes, improved understanding of the cause and pathogenesis of AF has provided important insights into the underlying pathophysiological mechanisms of this common arrhythmia and identified potential mechanism-based therapies.

Small-conductance Ca2+-activated K+ channels promote J-wave syndrome and phase 2 reentry

BACKGROUND

Small-conductance Ca²⁺-activated potassium (SK) channels play complex roles in cardiac arrhythmogenesis. SK channels colocalize with L-type Ca²⁺ channels, yet how this colocalization affects cardiac arrhythmogenesis is unknown.

OBJECTIVE

The purpose of this study was to investigate the role of colocalization of SK channels with L-type Ca²⁺ channels in promoting J-wave syndrome and ventricular arrhythmias.

METHODS

We carried out computer simulations of single-cell and tissue models. SK channels in the model were assigned to preferentially sense Ca²⁺ in the bulk cytosol, subsarcolemmal space, or junctional cleft.

RESULTS

When SK channels sense Ca²⁺ in the bulk cytosol, the SK current (I_SK) rises and decays slowly during an action potential, the action potential duration (APD) decreases as the maximum conductance increases, no complex APD dynamics and phase 2 reentry can be induced by I_SK. When SK channels sense Ca²⁺ in the subsarcolemmal space or junctional cleft, I_SK can rise and decay rapidly during an action potential in a spike-like pattern because of spiky Ca²⁺ transients in these compartments, which can cause spike-and-dome action potential morphology, APD alternans, J-wave elevation, and phase 2 reentry. Our results can account for the experimental finding that activation of I_SK induced J-wave syndrome and phase 2 reentry in rabbit hearts.

CONCLUSION

Colocalization of SK channels with L-type Ca²⁺ channels so that they preferentially sense Ca²⁺ in the subsarcolemmal or junctional space may result in a spiky I_SK, which can functionally play a similar role of the transient outward K⁺ current in promoting J-wave syndrome and ventricular arrhythmias.

Brugada syndrome: A comprehensive review of pathophysiological mechanisms and risk stratification strategies

Brugada syndrome (BrS) is an inherited ion channel channelopathy predisposing to ventricular arrhythmias and sudden cardiac death. Originally believed to be predominantly associated with mutations in SCN5A encoding for the cardiac sodium channel, mutations of 18 genes other than SCN5A have been implicated in the pathogenesis of BrS to date. Diagnosis is based on the presence of a spontaneous or drug-induced coved-type ST segment elevation. The predominant electrophysiological mechanism underlying BrS remains disputed, commonly revolving around the three main hypotheses based on abnormal repolarization, depolarization or current-load match. Evidence from computational modelling, pre-clinical and clinical studies illustrates that molecular abnormalities found in BrS lead to alterations in excitation wavelength (λ), which ultimately elevates arrhythmic risk. A major challenge for clinicians in managing this condition is the difficulty in predicting the subset of patients who will suffer from life-threatening ventricular arrhythmic events. Several repolarization risk markers have been used thus far, but these neglect the contributions of conduction abnormalities in the form of slowing and dispersion. Indices incorporating both repolarization and conduction based on the concept of λ have recently been proposed. These may have better predictive values than the existing markers. Current treatment options include pharmacological therapy to reduce the occurrence of arrhythmic events or to abort these episodes, and interventions such as implantable cardioverter-defibrillator insertion or radiofrequency ablation of abnormal arrhythmic substrate.

Icariin, a Novel Blocker of Sodium and Calcium Channels, Eliminates Early and Delayed Afterdepolarizations, As Well As Triggered Activity, in Rabbit Cardiomyocytes

Icariin, a flavonoid monomer from Herba Epimedii, has confirmed pharmacological and biological effects. However, its effects on arrhythmias and cardiac electrophysiology remain unclear. Here we investigate the effects of icariin on ion currents and action potentials (APs) in the rabbit myocardium. Furthermore, the effects of icariin on aconitine-induced arrhythmias were assessed in whole rabbits. Ion currents and APs were recorded in voltage-clamp and current-clamp mode in rabbit left ventricular myocytes (LVMs) and left atrial myocytes (LAMs), respectively. Icariin significantly shortened action potential durations (APDs) at 50 and 90% repolarization (APD₅₀ and APD₉₀) and reduced AP amplitude (APA) and the maximum upstroke velocity (V_max) of APs in LAMs and LVMs; however, icariin had no effect on resting membrane potential (RMP) in these cells. Icariin decreased the rate-dependence of the APD and completely abolished anemonia toxin II (ATX-II)-induced early afterdepolarizations (EADs). Moreover, icariin significantly suppressed delayed afterdepolarizations (DADs) and triggered activities (TAs) elicited by isoproterenol (ISO, 1 μM) and high extracellular calcium concentrations ([Ca²⁺]_o, 3.6 mM) in LVMs. Icariin also decreased I_NaT in a concentration-dependent manner in LAMs and LVMs, with IC₅₀ values of 12.28 ± 0.29 μM (n = 8 cells/4 rabbits) and 11.83 ± 0.92 μM (n = 10 cells/6 rabbits; p > 0.05 vs. LAMs), respectively, and reversed ATX-II-induced I_NaL in a concentration-dependent manner in LVMs. Furthermore, icariin attenuated I_CaL in a dose-dependent manner in LVMs. The corresponding IC₅₀ value was 4.78 ± 0.89 μM (n = 8 cells/4 rabbits), indicating that the aforementioned current in LVMs was 2.8-fold more sensitive to icariin than I_CaL in LAMs (13.43 ± 2.73 μM; n = 9 cells/5 rabbits). Icariin induced leftward shifts in the steady-state inactivation curves of I_NaT and I_CaL in LAMs and LVMs but did not have a significant effect on their activation processes. Moreover, icariin had no effects on I_K1 and I_Kr in LVMs or I_to and I_Kur in LAMs. These results revealed for the first time that icariin is a multichannel blocker that affects I_NaT, I_NaL and I_CaL in the myocardium and that the drug had significant inhibitory effects on aconitine-induced arrhythmias in whole rabbits. Therefore, icariin has potential as a class I and IV antiarrhythmic drug.

Memory-Induced Chaos in Cardiac Excitation

Excitable systems display memory, but how memory affects the excitation dynamics of such systems remains to be elucidated. Here we use computer simulation of cardiac action potential models to demonstrate that memory can cause dynamical instabilities that result in complex excitation dynamics and chaos. We develop an iterated map model that correctly describes these dynamics and show that memory converts a monotonic first return map of action potential duration into a nonmonotonic one, resulting in a period-doubling bifurcation route to chaos.

Electrophysiological Mechanisms of Brugada Syndrome: Insights from Pre-clinical and Clinical Studies

Brugada syndrome (BrS), is a primary electrical disorder predisposing affected individuals to sudden cardiac death via the development of ventricular tachycardia and fibrillation (VT/VF). Originally, BrS was linked to mutations in the SCN5A, which encodes for the cardiac Na+ channel. To date, variants in 19 genes have been implicated in this condition, with 11, 5, 3, and 1 genes affecting the Na+, K+, Ca2+, and funny currents, respectively. Diagnosis of BrS is based on ECG criteria of coved- or saddle-shaped ST segment elevation and/or T-wave inversion with or without drug challenge. Three hypotheses based on abnormal depolarization, abnormal repolarization, and current-load-mismatch have been put forward to explain the electrophysiological mechanisms responsible for BrS. Evidence from computational modeling, pre-clinical, and clinical studies illustrates that molecular abnormalities found in BrS lead to alterations in excitation wavelength (λ), which ultimately elevates arrhythmic risk. A major challenge for clinicians in managing this condition is the difficulty in predicting the subset of patients who will suffer from life-threatening VT/VF. Several repolarization risk markers have been used thus far, but these neglect the contributions of conduction abnormalities in the form of slowing and dispersion. Indices incorporating both repolarization and conduction and based on the concept of λ have recently been proposed. These may have better predictive values than the existing markers.

Perspective: a dynamics-based classification of ventricular arrhythmias

Despite key advances in the clinical management of life-threatening ventricular arrhythmias, culminating with the development of implantable cardioverter-defibrillators and catheter ablation techniques, pharmacologic/biologic therapeutics have lagged behind. The fundamental issue is that biological targets are molecular factors. Diseases, however, represent emergent properties at the scale of the organism that result from dynamic interactions between multiple constantly changing molecular factors. For a pharmacologic/biologic therapy to be effective, it must target the dynamic processes that underlie the disease. Here we propose a classification of ventricular arrhythmias that is based on our current understanding of the dynamics occurring at the subcellular, cellular, tissue and organism scales, which cause arrhythmias by simultaneously generating arrhythmia triggers and exacerbating tissue vulnerability. The goal is to create a framework that systematically links these key dynamic factors together with fixed factors (structural and electrophysiological heterogeneity) synergistically promoting electrical dispersion and increased arrhythmia risk to molecular factors that can serve as biological targets. We classify ventricular arrhythmias into three primary dynamic categories related generally to unstable Ca cycling, reduced repolarization, and excess repolarization, respectively. The clinical syndromes, arrhythmia mechanisms, dynamic factors and what is known about their molecular counterparts are discussed. Based on this framework, we propose a computational-experimental strategy for exploring the links between molecular factors, fixed factors and dynamic factors that underlie life-threatening ventricular arrhythmias. The ultimate objective is to facilitate drug development by creating an in silico platform to evaluate and predict comprehensively how molecular interventions affect not only a single targeted arrhythmia, but all primary arrhythmia dynamics categories as well as normal cardiac excitation-contraction coupling.