Calcium-voltage coupling in the genesis of early and delayed afterdepolarizations in cardiac myocytes

Antiarrhythmic effects of metformin

Atrial fibrillation/flutter (AF) is a major public health problem and is associated with stroke, heart failure, dementia, and death. It is estimated that 20%-30% of Americans will develop AF at some point in their life. Current medications to prevent AF have limited efficacy and significant adverse effects. Newer and safer therapies to prevent AF are needed. Ventricular arrhythmias are less prevalent than AF but may have significant consequences including sudden cardiac death. Metformin is the most prescribed, first-line medication for treatment of diabetes mellitus (DM). It decreases hepatic glucose production but also reduces inflammation and oxidative stress. Experimental studies have shown that metformin improves metabolic, electrical, and histologic risk factors associated with AF and ventricular arrhythmias. Furthermore, in large clinical observational studies, metformin has been associated with a reduced risk of AF in people with DM. These data suggest that metformin may have antiarrhythmic properties and may be a candidate to be repurposed as a medication to prevent cardiac arrhythmias. In this article, we review the clinical observational and experimental evidence for the association between metformin and cardiac arrhythmias. We also discuss the potential antiarrhythmic mechanisms underlying this association. Repurposing a well-tolerated, safe, and inexpensive medication to prevent cardiac arrhythmias has significant positive public health implications.

Enhanced Ca2+-Driven Arrhythmias in Female Patients with Atrial Fibrillation: Insights from Computational Modeling

BACKGROUND AND AIMS

Substantial sex-based differences have been reported in atrial fibrillation (AF), with female patients experiencing worse symptoms, increased complications from drug side effects or ablation, and elevated risk of AF-related stroke and mortality. Recent studies revealed sex-specific alterations in AF-associated Ca2+ dysregulation, whereby female cardiomyocytes more frequently exhibit potentially proarrhythmic Ca2+-driven instabilities compared to male cardiomyocytes. In this study, we aim to gain a mechanistic understanding of the Ca2+-handling disturbances and Ca2+-driven arrhythmogenic events in males vs females and establish their responses to Ca2+-targeted interventions.

METHODS AND RESULTS

We incorporated known sex differences and AF-associated changes in the expression and phosphorylation of key Ca2+-handling proteins and in ultrastructural properties and dimensions of atrial cardiomyocytes into our recently developed 3D atrial cardiomyocyte model that couples electrophysiology with spatially detailed Ca2+-handling processes. Our simulations of quiescent cardiomyocytes show increased incidence of Ca2+ sparks in female vs male myocytes in AF, in agreement with previous experimental reports. Additionally, our female model exhibited elevated propensity to develop pacing-induced spontaneous Ca2+ releases (SCRs) and augmented beat-to-beat variability in action potential (AP)-elicited Ca2+ transients compared with the male model. Parameter sensitivity analysis uncovered precise arrhythmogenic contributions of each component that was implicated in sex and/or AF alterations. Specifically, increased ryanodine receptor phosphorylation in female AF cardiomyocytes emerged as the major SCR contributor, while reduced L-type Ca2+ current was protective against SCRs for male AF cardiomyocytes. Furthermore, simulations of tentative Ca2+-targeted interventions identified potential strategies to attenuate Ca2+-driven arrhythmogenic events in female atria (e.g., t-tubule restoration, and inhibition of ryanodine receptor and sarcoplasmic/endoplasmic reticulum Ca2+-ATPase), and revealed enhanced efficacy when applied in combination.

CONCLUSIONS

Our sex-specific computational models of human atrial cardiomyocytes uncover increased propensity to Ca2+-driven arrhythmogenic events in female compared to male atrial cardiomyocytes in AF, and point to combined Ca2+-targeted interventions as promising approaches to treat AF in female patients. Our study establishes that AF treatment may benefit from sex-dependent strategies informed by sex-specific mechanisms.

Ventricular Tachycardia Due to Triggered Activity: Role of Early and Delayed Afterdepolarizations

Most forms of sustained ventricular tachycardia (VT) are caused by re-entry, resulting from altered myocardial conduction and refractoriness secondary to underlying structural heart disease. In contrast, VT caused by triggered activity (TA) is unrelated to an abnormal structural substrate and is often caused by molecular defects affecting ion channel function or regulation of intracellular calcium cycling. This review summarizes the cellular and molecular bases underlying TA and exemplifies their clinical relevance with selective representative scenarios. The underlying basis of TA caused by delayed afterdepolarizations is related to sarcoplasmic reticulum calcium overload, calcium waves, and diastolic sarcoplasmic reticulum calcium leak. Clinical examples of TA caused by delayed afterdepolarizations include sustained right and left ventricular outflow tract tachycardia and catecholaminergic polymorphic VT. The other form of afterpotentials, early afterdepolarizations, are systolic events and inscribe early afterdepolarizations during phase 2 or phase 3 of the action potential. The fundamental defect is a decrease in repolarization reserve with associated increases in late plateau inward currents. Malignant ventricular arrhythmias in the long QT syndromes are initiated by early afterdepolarization-mediated TA. An understanding of the molecular and cellular bases of these arrhythmias has resulted in generally effective pharmacologic-based therapies, but these are nonspecific agents that have off-target effects. Therapeutic efficacy may need to be augmented with an implantable defibrillator. Next-generation therapies will include novel agents that rescue arrhythmogenic abnormalities in cellular signaling pathways and gene therapy approaches that transfer or edit pathogenic gene variants or silence mutant messenger ribonucleic acid.

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.

Ryanodine receptor inhibition with acute dantrolene treatment reduces arrhythmia susceptibility in human hearts

Ryanodine receptor 2 (RyR2) hyperactivity is observed in structural heart diseases that are a result of ischemia or heart failure. It causes abnormal calcium handling and calcium leaks that cause metabolic, electrical, and mechanical dysfunction, which can trigger arrhythmias. Here, we tested the antiarrhythmic potential of dantrolene (RyR inhibitor) in human hearts. Human hearts not used in transplantation were obtained, and right ventricular outflow tract (RVOT) wedges and left ventricular (LV) slices were prepared. Pseudo-ECGs were recorded to determine premature ventricular contraction (PVC) incidences. Optical mapping was performed to determine arrhythmogenic substrates. After baseline optical recordings, tissues were treated with 1) isoproterenol (250 nM), 2) caffeine (200 mM), and 3) dantrolene (2 or 10 mM). Optical recordings were obtained after each treatment. Isoproterenol and caffeine treatment increased PVC incidence, whereas dantrolene reduced the PVC burden. Isoproterenol shortened action potential duration (APD) in the RV, RVOT, and LV regions and shortened calcium transient duration (CaTD) in the LV. Caffeine further shortened APD in the RV, did not modulate APD in the RVOT, and prolonged APD in the LV. In addition, in the LV, CaTD prolongation was also observed. More importantly, adding dantrolene did not alter APD in the RV or RVOT regions but produced a trend toward APD prolongation and significant CaTD prolongation in the LV, restoring these parameters to baseline values. In conclusions, dantrolene treatment suppresses triggers and reverses arrhythmogenic substrates in the human heart and could be a novel antiarrhythmic therapy in patients with structural heart disease.NEW & NOTEWORTHY Ryanodine receptor 2 hyperactivity is observed in structural heart diseases caused by ischemia or heart failure. It causes abnormal calcium leaks, which can trigger arrhythmias. To prevent arrhythmias, we applied dantrolene in human hearts ex vivo. Isoproterenol and caffeine treatment increased PVC incidence, whereas dantrolene reduced the PVC burden. Dantrolene treatment suppresses triggers and reverses arrhythmogenic substrates and could be a novel antiarrhythmic therapy in patients with structural heart disease.

Naringin exerts antiarrhythmic effects by inhibiting channel currents in mouse cardiomyocytes

INTRODUCTION

Naringin, a flavonoid extracted from citrus plants, has a variety of biological effects. Studies have shown that increasing the consumption of flavonoid-rich foods can reduce the incidence of cardiac arrhythmia. Naringin has been reported to have beneficial cardiovascular effects and thus can be used to prevent cardiovascular diseases, but the electrophysiological mechanism through which it prevents arrhythmias has not been elucidated. This study was conducted to investigate the effect of naringin on the transmembrane ion channel currents in mouse ventricular myocytes and the antiarrhythmic effect of this compound on Langendorff-perfused mouse hearts.

METHODS

Action potentials (APs) and ionic currents were recorded in isolated ventricular myocytes using the whole-cell patch-clamp technique. Anemone toxin II (ATX II) and CaCl2 were used to induce early afterdepolarizations (EADs) and delayed afterdepolarizations (DADs), respectively. Electrocardiogram (ECG) recordings were conducted in Langendorff-perfused mouse hearts with a BL-420F biological signal acquisition and analysis system.

RESULTS

At the cellular level, naringin shortened the action potential duration (APD) of ventricular myocytes and decreased the maximum depolarization velocity (Vmax) of APs.Naringin inhibited the L-type calcium current (ICa.L) and ATX II enhanced the late sodium current (INa.L) in a concentration-dependent manner with IC50 values of 508.5 μmol/L (n = 9) and 311.6 μmol/L (n = 10), respectively. In addition, naringin also inhibited the peak sodium current (INa·P) and delayed the rectifier potassium current (IK) and the transient outward potassium current (Ito). Moreover, naringin reduced ATX II-induced APD prolongation and EADs and had a significant inhibitory effect on CaCl2-induced DADs as well. At the organ level, naringin reduced the incidence of ventricular tachycardia (VT) and ventricular fibrillation (VF) induced by ATX II and shortened the duration of both in isolated hearts.

CONCLUSION

Naringin can inhibit the occurrence of EADs and DADs at the cellular level; furthermore, it can inhibit INa.L, ICa.L, INa·P, IK, and Ito in ventricular myocytes. Naringin also inhibits arrhythmias induced by ATX II in hearts. By investigating naringin with this electrophysiological method for the first time, we determined that this flavonoid may be a multichannel blocker with antiarrhythmic effects.

Electrical and Structural Insights into Right Ventricular Outflow Tract Arrhythmogenesis

The right ventricular outflow tract (RVOT) is the major origin of ventricular arrhythmias, including premature ventricular contractions, idiopathic ventricular arrhythmias, Brugada syndrome, torsade de pointes, long QT syndrome, and arrhythmogenic right ventricular cardiomyopathy. The RVOT has distinct developmental origins and cellular characteristics and a complex myocardial architecture with high shear wall stress, which may lead to its high vulnerability to arrhythmogenesis. RVOT myocytes are vulnerable to intracellular sodium and calcium overload due to calcium handling protein modulation, enhanced CaMKII activity, ryanodine receptor phosphorylation, and a higher cAMP level activated by predisposing factors or pathological conditions. A reduction in Cx43 and Scn5a expression may lead to electrical uncoupling in RVOT. The purpose of this review is to update the current understanding of the cellular and molecular mechanisms of RVOT arrhythmogenesis.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca2+ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca2+ -handling proteins such as the L-type Ca2+ channel and Na+ -Ca2+ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca2+ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca2+ -handling and Ca2+ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca2+ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca2+ release. We found that varying TATS density and thus the localization of key Ca2+ -handling proteins has profound effects on Ca2+ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca2+ concentration through decreased Ca2+ extrusion. This elevated Ca2+ increases RyR open probability causing spontaneous Ca2+ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca2+ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca2+ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca2+ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca2+ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca2+ handling governed by the TATS. Simulated TATS loss causes diastolic Ca2+ and voltage instabilities through reduced Na+ -Ca2+ exchanger-mediated Ca2+ removal, cleft Ca2+ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca2+ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca2+ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca2+ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

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.

Uremic toxins in chronic kidney disease highlight a fundamental gap in understanding their detrimental effects on cardiac electrophysiology and arrhythmogenesis

Chronic kidney disease (CKD) and cardiovascular disease (CVD) have an estimated 700-800 and 523 million cases worldwide, respectively, with CVD being the leading cause of death in CKD patients. The pathophysiological interplay between the heart and kidneys is defined as the cardiorenal syndrome (CRS), in which worsening of kidney function is represented by increased plasma concentrations of uremic toxins (UTs), culminating in dialysis patients. As there is a high incidence of CVD in CKD patients, accompanied by arrhythmias and sudden cardiac death, knowledge on electrophysiological remodeling would be instrumental for understanding the CRS. While the interplay between both organs is clearly of importance in CRS, the involvement of UTs in pro-arrhythmic remodeling is only poorly investigated, especially regarding the mechanistic background. Currently, the clinical approach against potential arrhythmic events is mainly restricted to symptom treatment, stressing the need for fundamental research on UT in relation to electrophysiology. This review addresses the existing knowledge of UTs and cardiac electrophysiology, and the experimental research gap between fundamental research and clinical research of the CRS. Clinically, mainly absorbents like ibuprofen and AST-120 are studied, which show limited safe and efficient usability. Experimental research shows disturbances in cardiac electrical activation and conduction after inducing CKD or exposure to UTs, but are scarcely present or focus solely on already well-investigated UTs. Based on UTs data derived from CKD patient cohort studies, a clinically relevant overview of physiological and pathological UTs concentrations is created. Using this, future experimental research is stimulated to involve electrophysiologically translatable animals, such as rabbits, or in vitro engineered heart tissues.

Experimental hypothyroidism induces cardiac arrhythmias and ranolazine reverts and prevents the phenotype

AIMS

Hypothyroidism is associated with an increased risk of cardiovascular disease and enhanced susceptibility to arrhythmias. In our investigation, we evaluated the potential involvement of late sodium current (I_Na,late) in cardiac arrhythmias in an experimental murine model of hypothyroidism.

MAIN METHODS

Male Swiss mice were treated with methimazole (0.1 % w/vol, during 21 days) to induce experimental hypothyroidism before ECG, action potential (AP) and intracellular Ca²⁺ dynamics were evaluated. Susceptibility to arrhythmia was measured in vitro and in vivo.

KEY FINDINGS

The results revealed that hypothyroid animals presented ECG alterations (e.g. increased QTc) with the presence of spontaneous sustained ventricular tachycardia. These changes were associated with depolarized resting membrane potential in isolated cardiomyocytes and increased AP duration and dispersion at 90 % of the repolarization. Aberrant AP waveforms were related to increased Ca²⁺ sparks and out-of-pace Ca²⁺ waves. These changes were observed in a scenario of enhanced I_Na,late. Interestingly, ranolazine, a clinically used blocker of I_Na,late, restored the ECG alterations, reduced Ca²⁺ sparks and aberrant waves, decreased the in vitro events and the severity of arrhythmias observed in isolated cardiomyocytes from hypothyroid animals. Using the in vivo dobutamine + caffeine protocol, animals with hypothyroidism developed catecholaminergic bidirectional ventricular tachycardia, but pre-treatment with ranolazine prevented this.

SIGNIFICANCE

We concluded that animals with hypothyroidism have increased susceptibility to developing arrhythmias and ranolazine, a clinically used blocker of I_Na,late, is able to correct the arrhythmic phenotype.

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.

Cardiac Macrophages and Their Effects on Arrhythmogenesis

Cardiac electrophysiology is a complex system established by a plethora of inward and outward ion currents in cardiomyocytes generating and conducting electrical signals in the heart. However, not only cardiomyocytes but also other cell types can modulate the heart rhythm. Recently, cardiac macrophages were demonstrated as important players in both electrophysiology and arrhythmogenesis. Cardiac macrophages are a heterogeneous group of immune cells including resident macrophages derived from embryonic and fetal precursors and recruited macrophages derived from circulating monocytes from the bone marrow. Recent studies suggest antiarrhythmic as well as proarrhythmic effects of cardiac macrophages. The proposed mechanisms of how cardiac macrophages affect electrophysiology vary and include both direct and indirect interactions with other cardiac cells. In this review, we provide an overview of the different subsets of macrophages in the heart and their possible interactions with cardiomyocytes under both physiologic conditions and heart disease. Furthermore, we elucidate similarities and differences between human, murine and porcine cardiac macrophages, thus providing detailed information for researchers investigating cardiac macrophages in important animal species for electrophysiologic research. Finally, we discuss the pros and cons of mice and pigs to investigate the role of cardiac macrophages in arrhythmogenesis from a translational perspective.

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.

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

Mettl3 deficiency leads to the upregulation of Cav1.2 and increases arrhythmia susceptibility in mice

Methyltransferase-like 3 (Mettl3) is a component of methyltransferase complex that mediates mA modification of RNAs, and participates in multiple biological processes. However, the role of Mettl3 in cardiac electrophysiology remains unknown. This study aims to explore the ventricular arrhythmia susceptibility of Mettl3 mice and the underlying mechanisms. Mice were anesthetized with 2% avertin (0.1 mL/ body weight) for echocardiography and programmed electrical pacing. Whole-cell patch clamp technique was used to examine the electrophysiological property of cardiomyocytes. The expression of Cav1.2 was determined by qRT-PCR and western blot analysis. The mA medication of mRNA was examined by MeRIP-Seq and MeRIP-qPCR. No differences are found in the morphology and function of the hearts between Mettl3 mice and wild-type (WT) controls. The QT and QTc intervals of Mettl3 mice are significantly longer. High-frequency electrical stimulation showed that heterozygous knockout of Mettl3 increases ventricular arrhythmia susceptibility. The whole-cell patch-clamp recordings showed that the APD is prolonged in Mettl3 ventricular myocytes and more EADs were observed. The density of is substantially increased in ventricular myocytes of Mettl3 mice. The pore-forming subunit of L-type calcium channel Cav1.2 is upregulated in Mettl3 mice, while the mRNA of its coding gene does not change. MeRIP-Seq and MeRIP-qPCR showed that the mA methylation of mRNA is decreased in cultured Mettl3-knockdown cardiomyocytes and Mettl3 hearts. Collectively, deficiency of Mettl3 increases ventricular arrhythmia susceptibility due to the upregulation of Cav1.2 by reducing mA modification onmRNA in mice. This study highlights the role of mA modification in the regulation of cardiac electrophysiology.

Thermal modulation of epicardial Ca2+ dynamics uncovers molecular mechanisms of Ca2+ alternans

Ca²⁺ alternans (Ca-Alts) are alternating beat-to-beat changes in the amplitude of Ca²⁺ transients that frequently occur during tachycardia, ischemia, or hypothermia that can lead to sudden cardiac death. Ca-Alts appear to result from a variation in the amount of Ca²⁺ released from the sarcoplasmic reticulum (SR) between two consecutive heartbeats. This variable Ca²⁺ release has been attributed to the alternation of the action potential duration, delay in the recovery from inactivation of RYR Ca²⁺ release channel (RYR2), or an incomplete Ca²⁺ refilling of the SR. In all three cases, the RYR2 mobilizes less Ca²⁺ from the SR in an alternating manner, thereby generating an alternating profile of the Ca²⁺ transients. We used a new experimental approach, fluorescence local field optical mapping (FLOM), to record at the epicardial layer of an intact heart with subcellular resolution. In conjunction with a local cold finger, a series of images were recorded within an area where the local cooling induced a temperature gradient. Ca-Alts were larger in colder regions and occurred without changes in action potential duration. Analysis of the change in the enthalpy and Q₁₀ of several kinetic processes defining intracellular Ca²⁺ dynamics indicated that the effects of temperature change on the relaxation of intracellular Ca²⁺ transients involved both passive and active mechanisms. The steep temperature dependency of Ca-Alts during tachycardia suggests Ca-Alts are generated by insufficient SERCA-mediated Ca²⁺ uptake into the SR. We found that Ca-Alts are heavily dependent on intra-SR Ca²⁺ and can be promoted through partial pharmacologic inhibition of SERCA2a. Finally, the FLOM experimental approach has the potential to help us understand how arrhythmogenesis correlates with the spatial distribution of metabolically impaired myocytes along the myocardium.

Chronic Ethanol Exposure Induces Deleterious Changes in Cardiomyocytes Derived from Human Induced Pluripotent Stem Cells

Chronic alcohol consumption in adults can induce cardiomyopathy, arrhythmias, and heart failure. In newborns, prenatal alcohol exposure can increase the risk of congenital heart diseases. Understanding biological mechanisms involved in the long-term alcohol exposure-induced cardiotoxicity is pivotal to the discovery of therapeutic strategies. In this study, cardiomyocytes derived from human pluripotent stem cells (hiPSC-CMs) were treated with clinically relevant doses of ethanol for various durations up to 5 weeks. The treated cells were characterized for their cellular properties and functions, and global proteomic profiling was conducted to understand the molecular changes associated with long-term ethanol exposure. Increased cell death, oxidative stress, deranged Ca²⁺ handling, abnormal action potential, altered contractility, and suppressed structure development were observed in ethanol-treated cells. Many dysregulated proteins identified by global proteomic profiling were involved in apoptosis, heart contraction, and extracellular collagen matrix. In addition, several signaling pathways including the Wnt and TGFβ signaling pathways were affected due to long-term ethanol treatment. Therefore, chronic ethanol treatment of hiPSC-CMs induces cardiotoxicity, impairs cardiac functions, and alters protein expression and signaling pathways. This study demonstrates the utility of hiPSC-CMs as a novel model for chronic alcohol exposure study and provides the molecular mechanisms associated with long-term alcohol exposure in human cardiomyocytes.

Functional Calsequestrin-1 Is Expressed in the Heart and Its Deficiency Is Causally Related to Malignant Hyperthermia-Like Arrhythmia

Background: Calsequestrins (Casqs), comprising the Casq1 and Casq2 isoforms, buffer Ca²⁺ and regulate its release in the sarcoplasmic reticulum (SR) of skeletal and cardiac muscle, respectively. Human inherited diseases associated with mutations in CASQ1 or CASQ2 include malignant hyperthermia/environmental heat stroke (MH/EHS) and catecholaminergic polymorphic ventricular tachycardia. However, patients with an MH/EHS event often suffer from arrhythmia for which the underlying mechanism remains unknown. Methods: Working hearts from conventional (Casq1-KO) and cardiac-specific (Casq1-CKO) Casq1 knockout mice were monitored in vivo and ex vivo by electrocardiogram and electrical mapping, respectively. MH was induced by 2% isoflurane and treated intraperitoneally with dantrolene. Time-lapse imaging was used to monitor intracellular Ca²⁺ activity in isolated mouse cardiomyocytes or neonatal rat ventricular myocytes (NRVMs) with knockdown, over-expression or truncation of the Casq1 gene. Conformational change in both Casqs was determined by crosslinking Western blot analysis. Results: Like MH/EHS patients, Casq1-KO and Casq1-CKO mice had faster basal heart rate, and ventricular tachycardia upon exposure to 2% isoflurane, which could be relieved by dantrolene. Basal sinus tachycardia and ventricular ectopic electrical triggering also occurred in Casq1-KO hearts ex vivo. Accordingly, the ventricular cardiomyocytes from Casq1-CKO mice displayed dantrolene-sensitive increased Ca²⁺ waves and diastole premature Ca²⁺ transients/oscillations upon isoflurane. NRVMs with Casq1-knockdown had enhanced spontaneous Ca2+ sparks/transients upon isoflurane, while cells over-expressing Casq1 exhibited decreased Ca2+ sparks/transients that were absent in cells with truncation of 9 amino acids at the C-terminus of Casq1. Structural evaluation showed that most of the Casq1 protein was present as a polymer and physically interacted with RyR2 in the ventricular SR. The Casq1 isoform was also expressed in human myocardium. Mechanistically, exposure to 2% isoflurane or heating at 41ºC induced Casq1 oligomerization in mouse ventricular and skeletal muscle tissues, leading to a reduced Casq1/RyR2 interaction and increased RyR2 activity in the ventricle. Conclusions: Casq1 is expressed in the heart, where it regulates SR Ca²⁺ release and heart rate. Casq1 deficiency independently causes MH/EHS-like ventricular arrhythmia by trigger-induced Casq1 oligomerization and a relief of its inhibitory effect on RyR2-mediated Ca²⁺ release, thus revealing a new inherited arrhythmia and a novel mechanism for MH/EHS arrhythmogenesis.

Efficacy of RyR2 inhibitor EL20 in induced pluripotent stem cell-derived cardiomyocytes from a patient with catecholaminergic polymorphic ventricular tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited cardiac arrhythmia syndrome that often leads to sudden cardiac death. The most common form of CPVT is caused by autosomal-dominant variants in the cardiac ryanodine receptor type-2 (RYR2) gene. Mutations in RYR2 promote calcium (Ca2+ ) leak from the sarcoplasmic reticulum (SR), triggering lethal arrhythmias. Recently, it was demonstrated that tetracaine derivative EL20 specifically inhibits mutant RyR2, normalizes Ca2+ handling and suppresses arrhythmias in a CPVT mouse model. The objective of this study was to determine whether EL20 normalizes SR Ca2+ handling and arrhythmic events in induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) from a CPVT patient. Blood samples from a child carrying RyR2 variant RyR2 variant Arg-176-Glu (R176Q) and a mutation-negative relative were reprogrammed into iPSCs using a Sendai virus system. iPSC-CMs were derived using the StemdiffTM kit. Confocal Ca2+ imaging was used to quantify RyR2 activity in the absence and presence of EL20. iPSC-CMs harbouring the R176Q variant demonstrated spontaneous SR Ca2+ release events, whereas administration of EL20 diminished these abnormal events at low nanomolar concentrations (IC50  = 82 nM). Importantly, treatment with EL20 did not have any adverse effects on systolic Ca2+ handling in control iPSC-CMs. Our results show for the first time that tetracaine derivative EL20 normalized SR Ca2+ handling and suppresses arrhythmogenic activity in iPSC-CMs derived from a CPVT patient. Hence, this study confirms that this RyR2-inhibitor represents a promising therapeutic candidate for treatment of CPVT.

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.

Evolution of mathematical models of cardiomyocyte electrophysiology

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

Activation of TRPC (Transient Receptor Potential Canonical) Channel Currents in Iron Overloaded Cardiac Myocytes

BACKGROUND

Arrhythmias and heart failure are common cardiac complications leading to substantial morbidity and mortality in patients with hemochromatosis, yet mechanistic insights remain incomplete. We investigated the effects of iron (Fe) on electrophysiological properties and intracellular Ca²⁺ (Ca²⁺_i) handling in mouse left ventricular cardiomyocytes.

METHODS

Cardiomyocytes were isolated from the left ventricle of mouse hearts and were superfused with Fe³⁺/8-hydroxyquinoline complex (5-100 μM). Membrane potential and ionic currents including TRPC (transient receptor potential canonical) were recorded using the patch-clamp technique. Ca²⁺_i was evaluated by using Fluo-4. Cell contraction was measured with a video-based edge detection system. The role of TRPCs in the genesis of arrhythmias was also investigated by using a mathematical model of a mouse ventricular myocyte with the incorporation of the TRPC component.

RESULTS

We observed prolongation of the action potential duration and induction of early and delayed afterdepolarizations in myocytes superfused with 15 µmol/L Fe³⁺/8-hydroxyquinoline complex. Iron treatment decreased the peak amplitude of the L-type Ca²⁺ current and total K⁺ current, altered Ca²⁺_i dynamics, and decreased cell contractility. During the final phase of Fe treatment, sustained Ca²⁺_i waves and repolarization failure occurred and ventricular cells became unexcitable. Gadolinium abolished Ca²⁺_i waves and restored the resting membrane potential to the normal range. The involvement of TRPC activation was confirmed by TRPC channel current recordings in the absence or presence of functional TRPC channel antibodies. Computer modeling captured the same action potential and Ca²⁺_i dynamics and provided additional mechanistic insights.

CONCLUSIONS

We conclude that iron overload induces cardiac dysfunction that is associated with TRPC channel activation and alterations in membrane potential and Ca²⁺_i dynamics.

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 (I_to), although it is not required, plays a key role in promoting spiral wave breakup in all five models. As the maximum conductance of I_to increases, it first promotes spiral wave breakup and then stabilizes the spiral waves. In the absence of I_to, 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 I_to. Increasing I_to promotes single-cell dynamical instabilities, including action potential duration alternans and chaos, and increasing I_to further suppresses these action potential dynamics. These cellular properties agree with the observation that increasing I_to 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 (I_to) is another current that plays a key role in spiral wave breakup, that is, spiral waves can be stable for low and high maximum I_to conductance but breakup occurs for intermediate maximum I_to conductance. Since I_to 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.

Mitochondrial depolarization promotes calcium alternans: Mechanistic insights from a ventricular myocyte model

Mitochondria are vital organelles inside the cell and contribute to intracellular calcium (Ca2+) dynamics directly and indirectly via calcium exchange, ATP generation, and production of reactive oxygen species (ROS). Arrhythmogenic Ca2+ alternans in cardiac myocytes has been observed in experiments under abnormal mitochondrial depolarization. However, complex signaling pathways and Ca2+ cycling between mitochondria and cytosol make it difficult in experiments to reveal the underlying mechanisms of Ca2+ alternans under abnormal mitochondrial depolarization. In this study, we use a newly developed spatiotemporal ventricular myocyte computer model that integrates mitochondrial Ca2+ cycling and complex signaling pathways to investigate the mechanisms of Ca2+ alternans during mitochondrial depolarization. We find that elevation of ROS in response to mitochondrial depolarization plays a critical role in promoting Ca2+ alternans. Further examination reveals that the redox effect of ROS on ryanodine receptors and sarco/endoplasmic reticulum Ca2+-ATPase synergistically promote alternans. Upregulation of mitochondrial Ca2+ uniporter promotes Ca2+ alternans via Ca2+-dependent mitochondrial permeability transition pore opening. Due to their relatively slow kinetics, oxidized Ca2+/calmodulin-dependent protein kinase II activation and ATP do not play significant roles acutely in the genesis of Ca2+ alternans after mitochondrial depolarization, but their roles can be significant in the long term, mainly through their effects on sarco/endoplasmic reticulum Ca2+-ATPase activity. In conclusion, mitochondrial depolarization promotes Ca2+ alternans acutely via the redox effect of ROS and chronically by ATP reduction. It suppresses Ca2+ alternans chronically through oxidized Ca2+/calmodulin-dependent protein kinase II activation.

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.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Delayed global feedback in the genesis and stability of spatiotemporal excitation patterns in paced biological excitable media

Biological excitable media, such as cardiac or neural cells and tissue, exhibit memory in which a change in the present excitation may affect the behaviors in the next excitation. For example, a change in calcium (Ca²⁺) concentration in a cell in the present excitation may affect the Ca²⁺ dynamics in the next excitation via bi-directional coupling between voltage and Ca²⁺, forming a delayed feedback loop. Since the Ca²⁺ dynamics inside the excitable cells are spatiotemporal while the membrane voltage is a global signal, the feedback loop is then a delayed global feedback (DGF) loop. In this study, we investigate the roles of DGF in the genesis and stability of spatiotemporal excitation patterns in periodically-paced excitable media using mathematical models with different levels of complexity: a model composed of coupled FitzHugh-Nagumo units, a 3-dimensional physiologically-detailed ventricular myocyte model, and a coupled map lattice model. We investigate the dynamics of excitation patterns that are temporal period-2 (P2) and spatially concordant or discordant, such as subcellular concordant or discordant Ca²⁺alternans in cardiac myocytes or spatially concordant or discordant Ca²⁺ and repolarization alternans in cardiac tissue. Our modeling approach allows both computer simulations and rigorous analytical treatments, which lead to the following results and conclusions. When DGF is absent, concordant and discordant P2 patterns occur depending on initial conditions with the discordant P2 patterns being spatially random. When the DGF is negative, only concordant P2 patterns exist. When the DGF is positive, both concordant and discordant P2 patterns can occur. The discordant P2 patterns are still spatially random, but they satisfy that the global signal exhibits a temporal period-1 behavior. The theoretical analyses of the coupled map lattice model reveal the underlying instabilities and bifurcations for the genesis, selection, and stability of spatiotemporal excitation patterns.

Understanding the mechanisms of pattern formation in biological systems is of great importance. Here we investigate the dynamical mechanisms by which delayed global feedback affects excitation pattern formation and stability in periodically-paced biological excitable media, such as cardiac or neural cells and tissue. We focus on the formation and stability of the temporal period-2 and spatially in-phase and out-of-phase excitation patterns. Using models of different levels of complexity, we show that when the delayed global feedback is negative, only the spatially in-phase patterns are stable. When the feedback is positive, both spatially in-phase and out-of-phase patterns are stable, and the out-of-phase patterns are spatially random but satisfy that the global signals are temporal period-1 solutions.

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.

Multi-scale approaches for the simulation of cardiac electrophysiology: I - Sub-cellular and stochastic calcium dynamics from cell to organ

Computational models of the heart at multiple spatial scales, from sub-cellular nanodomains to the whole-organ, are a powerful tool for the simulation of cardiac electrophysiology. Application of these models has provided remarkable insight into the normal and pathological functioning of the heart. In these two articles, we present methods for modelling cardiac electrophysiology at all of these spatial scales. In part one, presented here, we discuss methods and approaches for modelling sub-cellular calcium dynamics at the whole-cell and organ scales, valuable for modelling excitation-contraction coupling and mechanisms of arrhythmia triggers.

Molecular mechanisms of arrhythmogenic cardiomyopathy

Arrhythmogenic cardiomyopathy is a genetic disorder characterized by the risk of life-threatening arrhythmias, myocardial dysfunction and fibrofatty replacement of myocardial tissue. Mutations in genes that encode components of desmosomes, the adhesive junctions that connect cardiomyocytes, are the predominant cause of arrhythmogenic cardiomyopathy and can be identified in about half of patients with the condition. However, the molecular mechanisms leading to myocardial destruction, remodelling and arrhythmic predisposition remain poorly understood. Through the development of animal, induced pluripotent stem cell and other models of disease, advances in our understanding of the pathogenic mechanisms of arrhythmogenic cardiomyopathy over the past decade have brought several signalling pathways into focus. These pathways include canonical and non-canonical WNT signalling, the Hippo-Yes-associated protein (YAP) pathway and transforming growth factor-β signalling. These studies have begun to identify potential therapeutic targets whose modulation has shown promise in preclinical models. In this Review, we summarize and discuss the reported molecular mechanisms underlying the pathogenesis of arrhythmogenic cardiomyopathy.

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.

A Spatiotemporal Ventricular Myocyte Model Incorporating Mitochondrial Calcium Cycling

Intracellular calcium (Ca2+) cycling dynamics in cardiac myocytes are spatiotemporally generated by stochastic events arising from a spatially distributed network of coupled Ca2+ release units that interact with an intertwined mitochondrial network. In this study, we developed a spatiotemporal ventricular myocyte model that integrates mitochondria-related Ca2+ cycling components into our previously developed ventricular myocyte model consisting of a three-dimensional Ca2+ release unit network. Mathematical formulations of mitochondrial membrane potential, mitochondrial Ca2+ cycling, mitochondrial permeability transition pore stochastic opening and closing, intracellular reactive oxygen species signaling, and oxidized Ca2+/calmodulin-dependent protein kinase II signaling were incorporated into the model. We then used the model to simulate the effects of mitochondrial depolarization on mitochondrial Ca2+ cycling, Ca2+ spark frequency, and Ca2+ amplitude, which agree well with experimental data. We also simulated the effects of the strength of mitochondrial Ca2+ uniporters and their spatial localization on intracellular Ca2+ cycling properties, which substantially affected diastolic and systolic Ca2+ levels in the mitochondria but exhibited only a small effect on sarcoplasmic reticulum and cytosolic Ca2+ levels under normal conditions. We show that mitochondrial depolarization can cause Ca2+ waves and Ca2+ alternans, which agrees with previous experimental observations. We propose that this new, to our knowledge, spatiotemporal ventricular myocyte model, incorporating properties of mitochondrial Ca2+ cycling and reactive-oxygen-species-dependent signaling, will be useful for investigating the effects of mitochondria on intracellular Ca2+ cycling and action potential dynamics in ventricular myocytes.

Frequency and Phase Characteristics of Candle Flame Oscillation

The combustion of candles exhibits a variety of dynamical behaviors. Binding several candles together will result in flickering of candle flames, which is generally described as a nonlinear oscillator. The impact on the frequency of the flame by several factors, such as the arrangement, the number and the asymmetry of the oscillators, is discussed. Experimental results show that the frequency gradually decreases as the number of candles increases in the case of an isolated oscillator, while alternation between the in-phase and the anti-phase synchronization appears in a coupled system of two oscillators. Moreover, envelopes in the amplitude of the oscillatory luminance are displayed when candles are coupled asymmetrically. Since the coupling between oscillators is dominated by thermal radiation, a “overlapped peaks model” is proposed to phenomenologically explain the relationship between temperature distribution, coupling strength and the collective behavior in coupled system of candle oscillators in both symmetric and asymmetric cases.

Action Potential Prolongation, β-Adrenergic Stimulation, and Angiotensin II as Co-factors in Sarcoplasmic Reticulum Instability

Introduction: Increases in action potential duration (APD), genetic or acquired, and arrhythmias are often associated; nonetheless, the relationship between the two phenomena is inconstant, suggesting coexisting factors. β-adrenergic activation increases sarcoplasmic reticulum (SR) Ca²⁺-content; angiotensin II (ATII) may increase cytosolic Ca²⁺ and ROS production, all actions stimulating RyRs opening. Here we test how APD interacts with β-adrenergic and AT-receptor stimulation in facilitating spontaneous Ca²⁺ release events (SCR).

Methods: Under “action potential (AP) clamp”, guinea-pig cardiomyocytes (CMs) were driven with long (200 ms), normal (150 ms), and short (100 ms) AP waveforms at a CL of 500 ms; in a subset of CMs, all the 3 waveforms could be tested within the same cell. SCR were detected as inward current transients (I_TI) following repolarization; I_TI incidence and repetition within the same cycle were measured under increasing isoprenaline concentration ([ISO]) alone, or plus 100 nM ATII (30 min incubation+superfusion).

Results: I_TI incidence and repetition increased with [ISO]; at longer APs the [ISO]-response curve was shifted upward and I_TI coupling interval was reduced. ATII increased I_TI incidence more at low [ISO] and under normal (as compared to long) APs. Efficacy of AP shortening in suppressing I_TI decreased in ATII-treated myocytes and at higher [ISO].

Conclusions: AP prolongation sensitized the SR to the destabilizing actions of ISO and ATII. Summation of ISO, ATII and AP duration effects had a “saturating” effect on SCR incidence, thus suggesting convergence on a common factor (RyRs stability) “reset” by the occurrence of spontaneous Ca²⁺ release events.

Computational Model of Calcium Signaling in Cardiac Atrial Cells at the Submicron Scale

In cardiac cells, calcium is the mediator of excitation-contraction coupling. Dysfunctions in calcium handling have been identified as the origin of some cardiac arrhythmias. In the particular case of atrial myocytes, recent available experimental data has found links between these dysfunctions and structural changes in the calcium handling machinery (ryanodine cluster size and distribution, t-tubular network, etc). To address this issue, we have developed a computational model of an atrial myocyte that takes into account the detailed intracellular structure. The homogenized macroscopic behavior is described with a two-concentration field model, using effective diffusion coefficients of calcium in the sarcoplasmic reticulum (SR) and in the cytoplasm. The model reproduces the right calcium transients and dependence with pacing frequency. Under basal conditions, the calcium rise is mostly restricted to the periphery of the cell, with a large concentration ratio between the periphery and the interior. We have then studied the dependence of the speed of the calcium wave on cytosolic and SR diffusion coefficients, finding an almost linear relation with the former, in agreement with a diffusive and fire mechanism of propagation, and little dependence on the latter. Finally, we have studied the effect of a change in RyR cluster microstructure. We find that, under resting conditions, the spark frequency decreases slightly with RyR cluster spatial dispersion, but markedly increases when the RyRs are distributed in clusters of larger size, stressing the importance of RyR cluster organization to understand atrial arrhythmias, as recent experimental results suggest (Macquaide et al., 2015).

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.

Calcium-Mediated Oscillation in Membrane Potentials and Atrial-Triggered Activity in Atrial Cells of Casq2R33Q/R33Q Mutation Mice

Aim: We investigated the underlying mechanisms in atrial fibrillation (AF) associated with R33Q mutation and Ca²⁺-triggered activity. Methods and Results: We examined AF susceptibility with intraesophageal burst pacing in the sarcoplasmic reticulum (SR) Ca²⁺ leak model calsequestrin 2 R33Q (Casq2^R33Q/R33Q) mice. Atrial trigger appeared in R33Q mice but not WT mice (17.24%, 5/29 vs. 0.00%, 0/32, P < 0.05). AF was induced by 25 Hz pacing in R33Q mice (48.27%, 14/29 vs. 6.25%, 2/32, P < 0.01). The mice were given 1.5 mg/kg isoproterenol (Iso), and the incidences of AF increased (65.51%, 19/29 vs. 9.21%, 3/32, P < 0.01). Electrophysiology experiments and the recording of intracellular Ca²⁺ indicated significant increases in the Ca²⁺ sparks (5.24 ± 0.75 100 μM^-1.s^-1 vs. 0.29 ± 0.04 100 μM^-1.s^-1, n = 20, P < 0.05), intracellular free Ca²⁺ (0.238 ± 0.009 μM vs. 0.172 ± 0.006 μM, n = 20, P < 0.05), Ca²⁺ wave (11.74% vs. 2.24%, n = 20, P < 0.05), transient inward current (ITi) (-0.56 ± 0.02 pA/pF vs. -0.42 ± 0.01 pA/pF, n = 10, P < 0.05), and oscillation in membrane potentials (10.71%, 3/28 vs. 4.16%, 1/24, P < 0.05) in the R33Q group, but there was no significant difference in the L-type calcium current. These effects were enhanced by Iso, and the inhibition of calmodulin-dependent protein kinase II (CaMKII) by 1 μM KN93 reversed the effects of Iso on Ca²⁺ sparks (5.01 ± 0.66 100 μm^-1.s^-1 vs. 11.33 ± 1.63 100 μm^-1.s^-1, P < 0.05), intracellular Ca²⁺ (0.245 ± 0.005 μM vs. 0.324 ± 0.008 μM, P < 0.05), Ca²⁺ wave (12.35% vs. 17.83%, P < 0.05), ITi (-0.61 ± 0.02 pA/pF vs. -0.78 ± 0.03 pA/pF, n = 10, P < 0.05), and oscillation in membrane potential (17.85% 5/28 vs. 32.17% 9/28, P < 0.05). The reduction of ryanodine receptor 2 (RyR2) stable subunits (Casq2, triadin, and junctin) rather than RYR2 and the increase in CaMKII, phosphor-CaMKII, phosphor-RyR2 (Ser 2814), SERCA, and NCX1.1 was reflected in the R33Q group. Conclusion: This study demonstrates that the increase in spontaneous calcium elevations corresponding to ITi that may trigger the oscillation in membrane potentials in the R33Q group, thereby increasing the risk of AF. The occurrence of spontaneous calcium elevations in R33Q atrial myocytes is due to the dysfunction of RyR2 stable subunits, CaMKII hyperactivity, and CaMKII-mediated RyR phosphorylation. An effective therapeutic strategy to intervene in Ca²⁺-induced AF associated with the R33Q mutation may be through CaMKII inhibition.

Mapping and Ablation of Idiopathic Ventricular Fibrillation

Idiopathic ventricular fibrillation (IVF) is the main cause of unexplained sudden cardiac death, particularly in young patients under the age of 35. IVF is a diagnosis of exclusion in patients who have survived a VF episode without any identifiable structural or metabolic causes despite extensive diagnostic testing. Genetic testing allows identification of a likely causative mutation in up to 27% of unexplained sudden deaths in children and young adults. In the majority of cases, VF is triggered by PVCs that originate from the Purkinje network. Ablation of VF triggers in this setting is associated with high rates of acute success and long-term freedom from VF recurrence. Recent studies demonstrate that a significant subset of IVF defined by negative comprehensive investigations, demonstrate in fact subclinical structural alterations. These localized myocardial alterations are identified by high density electrogram mapping, are of small size and are mainly located in the epicardium. As reentrant VF drivers are often colocated with regions of abnormal electrograms, this localized substrate can be shown to be mechanistically linked with VF. Such areas may represent an important target for ablation.

Mitochondrial Ca2+ Influx Contributes to Arrhythmic Risk in Nonischemic Cardiomyopathy

Background

Heart failure (HF) is associated with increased arrhythmia risk and triggered activity. Abnormal Ca²⁺ handling is thought to underlie triggered activity, and mitochondria participate in Ca²⁺ homeostasis.

Methods and Results

A model of nonischemic HF was induced in C57BL/6 mice by hypertension. Computer simulations were performed using a mouse ventricular myocyte model of HF. Isoproterenol‐induced premature ventricular contractions and ventricular fibrillation were more prevalent in nonischemic HF mice than sham controls. Isolated myopathic myocytes showed decreased cytoplasmic Ca²⁺ transients, increased mitochondrial Ca²⁺ transients, and increased action potential duration at 90% repolarization. The alteration of action potential duration at 90% repolarization was consistent with in vivo corrected QT prolongation and could be explained by augmented L‐type Ca²⁺ currents, increased Na⁺‐Ca²⁺ exchange currents, and decreased total K⁺ currents. Of myopathic ventricular myocytes, 66% showed early afterdepolarizations (EADs) compared with 17% of sham myocytes (P<0.05). Intracellular application of 1 μmol/L Ru360, a mitochondrial Ca²⁺ uniporter–specific antagonist, could reduce mitochondrial Ca²⁺ transients, decrease action potential duration at 90% repolarization, and ameliorate EADs. Furthermore, genetic knockdown of mitochondrial Ca²⁺ uniporters inhibited mitochondrial Ca²⁺ uptake, reduced Na⁺‐Ca²⁺ exchange currents, decreased action potential duration at 90% repolarization, suppressed EADs, and reduced ventricular fibrillation in nonischemic HF mice. Computer simulations showed that EADs promoted by HF remodeling could be abolished by blocking either the mitochondrial Ca²⁺ uniporter or the L‐type Ca²⁺ current, consistent with the experimental observations.

Conclusions

Mitochondrial Ca²⁺ handling plays an important role in EADs seen with nonischemic cardiomyopathy and may represent a therapeutic target to reduce arrhythmic risk in this condition.

Divergent effects of adrenaline in human induced pluripotent stem cell-derived cardiomyocytes obtained from hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a common inherited cardiac disease that affects the heart muscle with diverse clinical outcomes. HCM can cause sudden cardiac death (SCD) during or immediately after mild to rigorous physical activity in young patients. However, the mechanism causing SCD as a result of exercise remains unknown, but exercise-induced ventricular arrhythmias are thought to be responsible for this fatal consequence. To understand the disease mechanism behind HCM in a better way, we generated patient-specific induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) from HCM patients carrying either the MYBPC3-Gln1061X or TPM1-Asp175Asn mutation. We extensively investigated the effects of low to high concentrations of adrenaline on action potential characteristics, and the occurrence of arrhythmias in the presence of various concentrations of adrenaline and in wash-out condition. We classified and quantified different types of arrhythmias observed in hiPSC-CMs, and found that the occurrence of arrhythmias was dependent on concentrations of adrenaline and positions of mutations in genes causing HCM. In addition, we observed ventricular tachycardia types of arrhythmias in hiPSC-CMs carrying the TPM1-Asp175Asn mutation. We additionally examined the antiarrhythmic potency of bisoprolol in HCM-specific hiPSC-CMs. However, bisoprolol could not reduce the occurrence of arrhythmias during administration or during the wash-out condition of adrenaline in HCM-specific hiPSC-CMs. Our study demonstrates hiPSC-CMs as a promising tool for studying HCM. The experimental design used in this study could be suitable and beneficial for studying other components and drugs related to cardiac disease in general.

Summary: Different concentrations of adrenaline have divergent effects during and immediately after administration in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) obtained from hypertrophic cardiomyopathy (HCM) patients. Bisoprolol could not reduce the arrhythmias in HCM-specific hiPSC-CMs.

Transverse tubular network structures in the genesis of intracellular calcium alternans and triggered activity in cardiac cells

RATIONALE

The major role of a transverse-tubular (TT) network in a cardiac cell is to facilitate effective excitation-contraction coupling and signaling. The TT network structures are heterogeneous within a single cell, and vary between different types of cells and species. They are also remodeled in cardiac diseases. However, how different TT network structures predispose cardiac cells to arrhythmogenesis remains to be revealed.

OBJECTIVE

To systematically investigate the roles of TT network structure and the underlying mechanisms in the genesis of intracellular calcium (Ca2+) alternans and triggered activity (TA).

METHODS AND RESULTS

Based on recent experimental observations, different TT network structures, including uniformly and non-uniformly random TT distributions, were modeled in a cardiac cell model consisting of a three-dimensional network of Ca2+ release units (CRUs). Our simulations showed that both Ca2+ alternans and Ca2+ wave-mediated TA were promoted when the fraction of orphaned CRUs was in an intermediate range, but suppressed in cells exhibiting either well-organized TT networks or low TT densities. Ca2+ alternans and TA could be promoted by low TT densities when the cells were small or the CRU coupling was strong. Both alternans and TA occurred more easily in uniformly random TT networks than in non-uniformly random TT networks. Subcellular spatially discordant Ca2+ alternans was promoted by non-uniformly random TT networks but suppressed by increasing CRU coupling strength. These mechanistic insights provide a holistic understanding of the effects of TT network structure on the susceptibility to arrhythmogenesis.

CONCLUSIONS

The TT network plays important roles in promoting Ca2+ alternans and TA, and different TT network structures may predispose cardiac cells differently to arrhythmogenesis.

Atrial Fibrillation Initiated by Early Afterdepolarization-Mediated Triggered Activity during Acute Oxidative Stress: Efficacy of Late Sodium Current Blockade

Background

The mechanism of Atrial Fibrillation (AF) that emerges spontaneously during acute oxidative stress is poorly defined and its drug therapy remains suboptimal. We hypothesized that oxidative activation of Ca-calmodulin dependent protein kinase (CaMKII) promotes Early Afterdepolarization-(EAD)-mediated triggered AF in aged fibrotic atria that is sensitive to late Na current (I_Na-L) blockade.

Method and Results

High-resolution voltage optical mapping of the Left and Right Atrial (LA & RA) epicardial surfaces along with microelectrode recordings were performed in isolated-perfused male Fisher 344 rat hearts in Langendorff setting. Aged atria (23-24 months) manifested 10-fold increase in atrial tissue fibrosis compared to young/adult (2-4 months) atria (P<0001. Spontaneous AF arose in 39 out of 41 of the aged atria but in 0 out of 12 young/adult hearts (P<001) during arterial perfusion of with 0.1 mm of hydrogen peroxide (H₂O₂). Optical Action Potential (AP) activation maps showed that the AF was initiated by a focal mechanism in the LA suggestive of EAD-mediated triggered activity. Cellular AP recordings with glass microelectrodes from the LA epicardial sites showing focal activity confirmed optical AP recordings that the spontaneous AF was initiated by late phase 3 EAD-mediated triggered activity. Inhibition of CaMKII activity with KN-93 (1 μM) (N=6) or its downstream target, the enhanced I_Na-L with GS-967 (1 μM), a specific blocker of I_Na-L (N=6), potently suppressed the AF and prevented its initiation when perfused 15 min prior to H₂O₂ (n=6).

Conclusions

Increased atrial tissue fibrosis combined with acute oxidative activation of CaMK II Initiate AF by EAD-mediated triggered activity. Specific block of the I_Na-L with GS-967 effectively suppresses the AF. Drug therapy of oxidative AF in humans with traditional antiarrhythmic drugs remains suboptimal; suppressing I_Na-L offers a potential new strategy for effective suppression of oxidative human AF that remains suboptimal.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model

Accumulation of intracellular Na⁺ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na⁺ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na⁺ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na⁺]_i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na⁺]_i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na⁺]_i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na⁺]_i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na⁺]_i may play complex roles in cellular and tissue-level cardiac dynamics.

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

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