Early afterdepolarizations in cardiac myocytes: beyond reduced repolarization reserve

COVID-19: Electrophysiological mechanisms underlying sudden cardiac death during exercise with facemasks

The mandatory use of facemasks is a public health measure implemented by various countries in response to the novel coronavirus disease 19 (COVID-19) pandemic. However, there have been case reports of sudden cardiac death (SCD) with the wearing of facemasks during exercise. In this paper, we hypothesize that exercise with facemasks may increase the risk of ventricular tachycardia/ventricular fibrillation (VT/VF) leading to SCD via the development of acute and/or intermittent hypoxia and hypercapnia. We discuss the potential underlying mechanisms including increases in adrenergic stimulation and oxidative stress leading to electrophysiological abnormalities that promote arrhythmias via non-reentrant and reentrant mechanisms. Given the interplay of multiple variables contributing to the increased arrhythmic risk, we advise avoidance of a facemask during high intensity exercise, or if wearing of a mask is mandatory, exercise intensity should remain low to avoid precipitation of lethal arrhythmias. However, we cannot exclude the possibility of an arrhythmic substrate even with low intensity exercise especially in those with established chronic cardiovascular disease in whom baseline electrophysiological abnormalities may be found.

Triggered Ca2+ Waves Induce Depolarization of Maximum Diastolic Potential and Action Potential Prolongation in Dog Atrial Myocytes

BACKGROUND

We have identified a novel form of abnormal Ca2+ wave activity in normal and failing dog atrial myocytes which occurs during the action potential (AP) and is absent during diastole. The goal of this study was to determine if triggered Ca2+ waves affect cellular electrophysiological properties.

METHODS

Simultaneous recordings of intracellular Ca2+ and APs allowed measurements of maximum diastolic potential and AP duration during triggered calcium waves (TCWs) in isolated dog atrial myocytes. Computer simulations then explored electrophysiological behavior arising from TCWs at the tissue scale.

RESULTS

At 3.3 to 5 Hz, TCWs occurred during the AP and often outlasted several AP cycles. Maximum diastolic potential was reduced, and AP duration was significantly prolonged during TCWs. All electrophysiological responses to TCWs were abolished by SEA0400 and ORM10103, indicating that Na-Ca exchange current caused depolarization. The time constant of recovery from inactivation of Ca2+ current was 40 to 70 ms in atrial myocytes (depending on holding potential) so this current could be responsible for AP activation during depolarization induced by TCWs. Modeling studies demonstrated that the characteristic properties of TCWs are potentially arrhythmogenic by promoting both conduction block and reentry arising from the depolarization induced by TCWs.

CONCLUSIONS

Triggered Ca2+ waves activate inward NCX and dramatically reduce atrial maximum diastolic potential and prolong AP duration, establishing the substrate for reentry which could contribute to the initiation and maintenance of atrial arrhythmias.

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.

Classifying Drugs by their Arrhythmogenic Risk Using Machine Learning

All medications have adverse effects. Among the most serious of these are cardiac arrhythmias. Current paradigms for drug safety evaluation are costly, lengthy, conservative, and impede efficient drug development. Here, we combine multiscale experiment and simulation, high-performance computing, and machine learning to create a risk estimator to stratify new and existing drugs according to their proarrhythmic potential. We capitalize on recent developments in machine learning and integrate information across 10 orders of magnitude in space and time to provide a holistic picture of the effects of drugs, either individually or in combination with other drugs. We show, both experimentally and computationally, that drug-induced arrhythmias are dominated by the interplay between two currents with opposing effects: the rapid delayed rectifier potassium current and the L-type calcium current. Using Gaussian process classification, we create a classifier that stratifies drugs into safe and arrhythmic domains for any combinations of these two currents. We demonstrate that our classifier correctly identifies the risk categories of 22 common drugs exclusively on the basis of their concentrations at 50% current block. Our new risk assessment tool explains under which conditions blocking the L-type calcium current can delay or even entirely suppress arrhythmogenic events. Using machine learning in drug safety evaluation can provide a more accurate and comprehensive mechanistic assessment of the proarrhythmic potential of new drugs. Our study paves the way toward establishing science-based criteria to accelerate drug development, design safer drugs, and reduce heart rhythm disorders.

Classifying Drugs by their Arrhythmogenic Risk Using Machine Learning

All medications have adverse effects. Among the most serious of these are cardiac arrhythmias. Current paradigms for drug safety evaluation are costly, lengthy, conservative, and impede efficient drug development. Here, we combine multiscale experiment and simulation, high-performance computing, and machine learning to create a risk estimator to stratify new and existing drugs according to their proarrhythmic potential. We capitalize on recent developments in machine learning and integrate information across 10 orders of magnitude in space and time to provide a holistic picture of the effects of drugs, either individually or in combination with other drugs. We show, both experimentally and computationally, that drug-induced arrhythmias are dominated by the interplay between two currents with opposing effects: the rapid delayed rectifier potassium current and the L-type calcium current. Using Gaussian process classification, we create a classifier that stratifies drugs into safe and arrhythmic domains for any combinations of these two currents. We demonstrate that our classifier correctly identifies the risk categories of 22 common drugs exclusively on the basis of their concentrations at 50% current block. Our new risk assessment tool explains under which conditions blocking the L-type calcium current can delay or even entirely suppress arrhythmogenic events. Using machine learning in drug safety evaluation can provide a more accurate and comprehensive mechanistic assessment of the proarrhythmic potential of new drugs. Our study paves the way toward establishing science-based criteria to accelerate drug development, design safer drugs, and reduce heart rhythm disorders.

Thioridazine Induces Cardiotoxicity via Reactive Oxygen Species-Mediated hERG Channel Deficiency and L-Type Calcium Channel Activation

Thioridazine (THIO) is a phenothiazine derivative that is mainly used for the treatment of psychotic disorders. However, cardiac arrhythmias especially QT interval prolongation associated with the application of this compound have received serious attention after its introduction into clinical practice, and the mechanisms underlying the cardiotoxicity induced by THIO have not been well defined. The present study was aimed at exploring the long-term effects of THIO on the hERG and L-type calcium channels, both of which are relevant to the development of QT prolongation. The hERG current (I _hERG) and the calcium current (I _Ca-L) were measured by patch clamp techniques. Protein levels were analyzed by Western blot, and channel-chaperone interactions were determined by coimmunoprecipitation. Reactive oxygen species (ROS) were determined by flow cytometry and laser scanning confocal microscopy. Our results demonstrated that THIO induced hERG channel deficiency but did not alter channel kinetics. THIO promoted ROS production and stimulated endoplasmic reticulum (ER) stress and the related proteins. The ROS scavenger N-acetyl cysteine (NAC) significantly attenuated hERG reduction induced by THIO and abolished the upregulation of ER stress marker proteins. Meanwhile, THIO increased the degradation of hERG channels via disrupting hERG-Hsp70 interactions. The disordered hERG proteins were degraded in proteasomes after ubiquitin modification. On the other hand, THIO increased I _Ca-L density and intracellular Ca²⁺ ([Ca²⁺]_i) in neonatal rat ventricular cardiomyocytes (NRVMs). The specific CaMKII inhibitor KN-93 attenuated the intracellular Ca²⁺ overload, indicating that ROS-mediated CaMKII activation promoted calcium channel activation induced by THIO. Optical mapping analysis demonstrated the slowing effects of THIO on cardiac repolarization in mouse hearts. THIO significantly prolonged APD₅₀ and APD₉₀ and increased the incidence of early afterdepolarizations (EADs). In human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), THIO also resulted in APD prolongation. In conclusion, dysfunction of hERG channel proteins and activation of L-type calcium channels via ROS production might be the ionic mechanisms for QT prolongation induced by THIO.

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.

Ginsenoside Rb1 exerts antiarrhythmic effects by inhibiting INa and ICaL in rabbit ventricular myocytes

Ginsenoside Rb1 exerts its pharmacological action by regulating sodium, potassium and calcium ion channels in the membranes of nerve cells. These ion channels are also present in cardiomyocytes, but no studies have been reported to date regarding the effects of Rb1 on cardiac sodium currents (I_Na), L-type calcium currents (I_CaL) and action potentials (APs). Additionally, the antiarrhythmic potential of Rb1 has not been assessed. In this study, we used a whole-cell patch clamp technique to assess the effect of Rb1 on these ion channels. The results showed that Rb1 inhibited I_Na and I_CaL, reduced the action potential amplitude (APA) and maximum upstroke velocity (V_max), and shortened the action potential duration (APD) in a concentration-dependent manner but had no effect on the inward rectifier potassium current (I_K1), delayed rectifier potassium current (I_K) or resting membrane potential (RMP). We also designed a pathological model at the cellular and organ level to verify the role of Rb1. The results showed that Rb1 abolished high calcium-induced delayed afterdepolarizations (DADs), depressed the increase in intracellular calcium ([Ca²⁺]_i), relieved calcium overload and protected cardiomyocytes. Rb1 can also reduce the occurrence of ventricular premature beats (VPBs) and ventricular tachycardia (VT) in ischemia-reperfusion (I-R) injury.

Allicin disrupts cardiac Cav1.2 channels via trafficking

CONTEXT

Allicin is a potential antiarrhythmic agent. The antiarrhythmic properties of allicin rely on its blockade of various cardiac ion channels. The l-type calcium (Cav1.2) channel provides a pivotal substrate for cardiac electrophysiologic activities. The mechanism of allicin on Cav1.2 remains unclear.

OBJECTIVE

This study evaluated the potential of allicin on the synthesis and trafficking of Cav1.2 channels.

MATERIALS AND METHODS

Primary cardiomyocytes (CMs) from neonatal Sprague-Dawley (SD) rats were exposed to allicin (0, 0.0001, 0.001, 0.01, 0.1, 1, 10, 100 μg/mL) for 24 and 48 h. The CellTiter-Glo assay was performed to measure CM viability. Western blot with grayscale analysis and confocal laser scanning microscopy were used to evaluate the effects of allicin on the synthesis and trafficking of Cav1.2 channel proteins in primary CMs.

RESULTS

There was no significant difference in apoptotic toxicity from the actual cell viability (p > 0.05) in any group (0, 0.0001, 0.001, 0.01, 0.1, 1, 10, 100 μg/mL allicin), except that viability in the 0.001 and 0.01 μg/mL groups at 24 h were significantly greater (137.37 and 135.96%) (p < 0.05). Western blot with grayscale analysis revealed no substantial inhibition by allicin of the synthesis of Cav1.2 proteins. Confocal laser scanning microscopy revealed trafficking dysfunction of Cav1.2 channels caused by allicin in primary CMs.

CONCLUSION

This study is the first to demonstrate that allicin inhibits cardiac Cav1.2 channels by disrupting trafficking, possibly mediating its antiarrhythmic benefits. Therefore, allicin might serve as a new antiarrhythmic agent in the future.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Bifurcations Caused by Feedback between Voltage and Intracellular Ion Concentrations in Ventricular Myocytes

We develop an iterated map model to describe the bifurcations and complex dynamics caused by the feedback between voltage and intracellular Ca^{2+} and Na^{+} concentrations in paced ventricular myocytes. Voltage and Ca^{2+} can form either a positive or a negative feedback loop, while voltage and Na^{+} form a negative feedback loop. Under certain diseased conditions, when the feedback between voltage and Ca^{2+} is positive, Hopf bifurcations occur, leading to periodic oscillatory behaviors. When this feedback is negative, period-doubling bifurcation routes to alternans and chaos occur.

Transmural Autonomic Regulation of Cardiac Contractility at the Intact Heart Level

The relationship between cardiac excitability and contractility depends on when Ca²⁺ influx occurs during the ventricular action potential (AP). In mammals, it is accepted that Ca²⁺ influx through the L-type Ca²⁺ channels occurs during AP phase 2. However, in murine models, experimental evidence shows Ca²⁺ influx takes place during phase 1. Interestingly, Ca²⁺ influx that activates contraction is highly regulated by the autonomic nervous system. Indeed, autonomic regulation exerts multiple effects on Ca²⁺ handling and cardiac electrophysiology. In this paper, we explore autonomic regulation in endocardial and epicardial layers of intact beating mice hearts to evaluate their role on cardiac excitability and contractility. We hypothesize that in mouse cardiac ventricles the influx of Ca²⁺ that triggers excitation–contraction coupling (ECC) does not occur during phase 2. Using pulsed local field fluorescence microscopy and loose patch photolysis, we show sympathetic stimulation by isoproterenol increased the amplitude of Ca²⁺ transients in both layers. This increase in contractility was driven by an increase in amplitude and duration of the L-type Ca²⁺ current during phase 1. Interestingly, the β-adrenergic increase of Ca²⁺ influx slowed the repolarization of phase 1, suggesting a competition between Ca²⁺ and K⁺ currents during this phase. In addition, cAMP activated L-type Ca²⁺ currents before SR Ca²⁺ release activated the Na⁺-Ca²⁺ exchanger currents, indicating Ca_v1.2 channels are the initial target of PKA phosphorylation. In contrast, parasympathetic stimulation by carbachol did not have a substantial effect on amplitude and kinetics of endocardial and epicardial Ca²⁺ transients. However, carbachol transiently decreased the duration of the AP late phase 2 repolarization. The carbachol-induced shortening of phase 2 did not have a considerable effect on ventricular pressure and systolic Ca²⁺ dynamics. Interestingly, blockade of muscarinic receptors by atropine prolonged the duration of phase 2 indicating that, in isolated hearts, there is an intrinsic release of acetylcholine. In addition, the acceleration of repolarization induced by carbachol was blocked by the acetylcholine-mediated K⁺ current inhibition. Our results reveal the transmural ramifications of autonomic regulation in intact mice hearts and support our hypothesis that Ca²⁺ influx that triggers ECC occurs in AP phase 1 and not in phase 2.

Predicting critical drug concentrations and torsadogenic risk using a multiscale exposure-response simulator

Torsades de pointes is a serious side effect of many drugs that can trigger sudden cardiac death, even in patients with structurally normal hearts. Torsadogenic risk has traditionally been correlated with the blockage of a specific potassium channel and a prolonged recovery period in the electrocardiogram. However, the precise mechanisms by which single channel block translates into heart rhythm disorders remain incompletely understood. Here we establish a multiscale exposure-response simulator that converts block-concentration characteristics from single cell recordings into three-dimensional excitation profiles and electrocardiograms to rapidly assess torsadogenic risk. For the drug dofetilide, we characterize the QT interval and heart rate at different drug concentrations and identify the critical concentration at the onset of torsades de pointes: For dofetilide concentrations of 2x, 3x, and 4x, as multiples of the free plasma concentration Cmax = 2.1 nM, the QT interval increased by +62.0%, +71.2%, and +82.3% compared to baseline, and the heart rate changed by -21.7%, -23.3%, and +88.3%. The last number indicates that, at the critical concentration of 4x, the heart spontaneously developed an episode of a torsades-like arrhythmia. Strikingly, this critical drug concentration is higher than the concentration estimated from early afterdepolarizations in single cells and lower than in one-dimensional cable models. Our results highlight the importance of whole heart modeling and explain, at least in part, why current regulatory paradigms often fail to accurately quantify the pro-arrhythmic potential of a drug. Our exposure-response simulator could provide a more mechanistic assessment of pro-arrhythmic risk and help establish science-based guidelines to reduce rhythm disorders, design safer drugs, and accelerate drug development.

Dynamic clamping human and rabbit atrial calcium current: narrowing ICaL window abolishes early afterdepolarizations

Key points

Early‐afterdepolarizations (EADs) are abnormal action potential oscillations and a known cause of cardiac arrhythmias. Ventricular EADs involve reactivation of a Ca²⁺ current (I_CaL) in its ‘window region’ voltage range. However, electrical mechanisms of atrial EADs, a potential cause of atrial fibrillation, are poorly understood.

Atrial cells were obtained from consenting patients undergoing heart surgery, as well as from rabbits. I_CaL was blocked with nifedipine and then a hybrid patch clamp/mathematical‐modelling technique, ‘dynamic clamping’, was used to record action potentials at the same time as injecting an artificial, modifiable, I_CaL (I_CaL,D‐C).

Progressively widening the I_CaL,D‐C window region produced EADs of various types, dependent on window width. EAD production was strongest upon moving the activation (vs. inactivation) side of the window.

EADs were then induced by a different method: increasing I_CaL,D‐C amplitude and/or K⁺ channel‐blockade (4‐aminopyridine). Narrowing of the I_CaL,D‐C window by ∼10 mV abolished these EADs.

Atrial I_CaL window narrowing is worthy of further testing as a potential anti‐atrial fibrillation drug mechanism.

Abstract

Atrial early‐afterdepolarizations (EADs) may contribute to atrial fibrillation (AF), perhaps involving reactivation of L‐type Ca²⁺ current (I_CaL) in its window region voltage range. The present study aimed (i) to validate the dynamic clamp technique for modifying the I_CaL contribution to atrial action potential (AP) waveform; (ii) to investigate the effects of widening the window I_CaL on EAD‐propensity; and (iii) to test whether EADs from increased I_CaL and AP duration are supressed by narrowing the window I_CaL. I_CaL and APs were recorded from rabbit and human atrial myocytes by whole‐cell‐patch clamp. During AP recording, I_CaL was inhibited (3 µm nifedipine) and replaced by a dynamic clamp model current, I_CaL,D‐C (tuned to native I_CaL characteristics), computed in real‐time (every 50 µs) based on myocyte membrane potential. I_CaL,D‐C‐injection restored the nifedipine‐suppressed AP plateau. Widening the window I_CaL,D‐C, symmetrically by stepwise simultaneous equal shifts of half‐voltages (V_0.5) of I_CaL,D‐C activation (negatively) and inactivation (positively), generated EADs (single, multiple or preceding repolarization failure) in a window width‐dependent manner, as well as AP alternans. A stronger EAD‐generating effect resulted from independently shifting activation V_0.5 (asymmetrical widening) than inactivation V_0.5; for example, a 15 mV activation shift produced EADs in nine of 17 (53%) human atrial myocytes vs. 0 of 18 from inactivation shift (P < 0.05). In 11 rabbit atrial myocytes in which EADs were generated either by increasing the conductance of normal window width I_CaL,D‐C or subsequent 4‐aminopyridine (2 mm), window I_CaL,D‐C narrowing (10 mV) abolished EADs of all types (P < 0.05). The present study validated the dynamic clamp for I_CaL, which is novel in atrial cardiomyocytes, and showed that EADs of various types are generated by widening (particularly asymmetrically) the window I_CaL, as well as abolished by narrowing it. Window I_CaL narrowing is a potential therapeutic mechanism worth pursuing in the search for improved anti‐AF drugs.

Early‐afterdepolarizations (EADs) are abnormal action potential oscillations and a known cause of cardiac arrhythmias. Ventricular EADs involve reactivation of a Ca²⁺ current (I_CaL) in its ‘window region’ voltage range. However, electrical mechanisms of atrial EADs, a potential cause of atrial fibrillation, are poorly understood.

Atrial cells were obtained from consenting patients undergoing heart surgery, as well as from rabbits. I_CaL was blocked with nifedipine and then a hybrid patch clamp/mathematical‐modelling technique, ‘dynamic clamping’, was used to record action potentials at the same time as injecting an artificial, modifiable, I_CaL (I_CaL,D‐C).

Progressively widening the I_CaL,D‐C window region produced EADs of various types, dependent on window width. EAD production was strongest upon moving the activation (vs. inactivation) side of the window.

EADs were then induced by a different method: increasing I_CaL,D‐C amplitude and/or K⁺ channel‐blockade (4‐aminopyridine). Narrowing of the I_CaL,D‐C window by ∼10 mV abolished these EADs.

Atrial I_CaL window narrowing is worthy of further testing as a potential anti‐atrial fibrillation drug mechanism.

Overexpression of the HCN2 channel increases the arrhythmogenicity induced by hypokalemia

Hypokalemia, an abnormally low level of potassium (K⁺), is a electrolyte imbalance that commonly occurs in heart failure patients. Hypokalemia is well known to induce lethal ventricular arrhythmia. However, the effects of hypokalemia in failing hearts that have undergone electrophysiological remodeling, i.e., the reactivation of fetal-type ion channels, remain unexplored. We have examined the effect of hypokalemia in the myocytes of transgenic mice overexpressing the hyperpolarization-activated, cyclic nucleotide-sensitive (HCN) channel in the heart (HCN2-Tg mice). Perfusion with a mild hypokalemic solution containing 3 mM K⁺ induced ectopic ventricular automaticity in 55.0% of HCN2-Tg mouse myocytes. In the remaining HCN2-Tg mouse myocytes, the resting membrane potential (RMP) was more depolarized than that of wild-type myocytes subjected to the same treatment and could also be hyperpolarized by an HCN channel blocker. We conclude that in hypokalemia in our mice model, the HCN2 channel was constitutively activated at the hyperpolarized RMP, thereby destabilizing the electrophysiological activity of ventricular myocytes.

Integrative and theoretical research on the architecture of a biological system and its disorder

An organism stems from assemblies of a variety of cells and proteins. This complex system serves as a unit, and it exhibits highly sophisticated functions in response to exogenous stimuli that change over time. The complete sequencing of the entire human genome has allowed researchers to address the enigmas of life and disease at the gene- or molecular-based level. The consequence of such studies is the rapid accumulation of a multitude of data at multiple levels, ranging from molecules to the whole body, that has necessitated the development of entirely new concepts, tools, and methodologies to analyze and integrate these data. This necessity has given birth to systems biology, an advanced theoretical and practical research framework that has totally changed the directions of not only basic life science but also medicine. During the symposium of the 95th Annual Meeting of The Physiological Society of Japan 2018, five researchers reported on their respective studies on systems biology. The topics included reactions of drugs, ion-transport architecture in an epithelial system, multi-omics in renal disease, cardiac electrophysiological systems, and a software platform for computer simulation. In this review article these authors have summarized recent achievements in the field and discuss next-generation studies on health and disease.

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

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

Impact of disease state on arrhythmic event detection by action potential modelling in cardiac safety pharmacology

INTRODUCTION

The use of in silico cardiac action potential simulations is one of the pillars of the CiPA initiative (Comprehensive in vitro Proarrhythmia Assay) currently under evaluation designed to detect more accurately proarrhythmic liabilities of new drug candidate. In order to take into account the variability of clinical situations, we propose to improve this method by studying the impact of various disease states on arrhythmic events induced by 30 torsadogenic or non-torsadogenic compounds.

METHOD

In silico modelling was done on the human myocytes using the Dutta revised O'Hara-Rudy algorithm. Results were analysed using a new metric based on the compound IC_50s against the seven cardiac ionic currents considered to be the most important by the CiPA initiative (I_Kr, I_Ks, I_Na, I_NaL, I_K1, I_to, I_CaL) and the minimal rate of action potential voltage decrease calculated at the early-afterdepolarization (EAD) take-off membrane voltage (V_min).

RESULTS

The specific threshold at which each torsadogenic compounds induced EAD, was exacerbated by the presence of cardiac risk factors ranked as follows: congestive heart failure > hypertrophic cardiomyopathy > cardiac pause > no risk factor. Non-torsadogenic compounds induced no EAD even in the presence of cardiac risk factors.

DISCUSSION

The present study highlighted the impact of pre-existing cardiovascular disease on arrhythmic event detection suggesting that disease state modelling may need to be incorporated in order to fully realize the goal of the CiPA paradigm in a more accurate predictability of proarrhythmic liabilities of new drug candidate.

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.

NCX-Mediated Subcellular Ca2+ Dynamics Underlying Early Afterdepolarizations in LQT2 Cardiomyocytes

Long QT syndrome type 2 (LQT2) is a congenital disease characterized by loss of function mutations in hERG potassium channels (I_Kr). LQT2 is associated with fatal ventricular arrhythmias promoted by triggered activity in the form of early afterdepolarizations (EADs). We previously demonstrated that intracellular Ca²⁺ handling is remodeled in LQT2 myocytes. Remodeling leads to aberrant late RyR-mediated Ca²⁺ releases that drive forward-mode Na⁺-Ca²⁺ exchanger (NCX) current and slow repolarization to promote reopening of L-type calcium channels and EADs. Forward-mode NCX was found to be enhanced despite the fact that these late releases do not significantly alter the whole-cell cytosolic calcium concentration during a vulnerable period of phase 2 of the action potential corresponding to the onset of EADs. Here, we use a multiscale ventricular myocyte model to explain this finding. We show that because the local NCX current is a saturating nonlinear function of the local submembrane calcium concentration, a larger number of smaller-amplitude discrete Ca²⁺ release events can produce a large increase in whole-cell forward-mode NCX current without increasing significantly the whole-cell cytosolic calcium concentration. Furthermore, we develop novel insights, to our knowledge, into how alterations of stochastic RyR activity at the single-channel level cause late aberrant Ca²⁺ release events. Experimental measurements in transgenic LTQ2 rabbits confirm the critical arrhythmogenic role of NCX and identify this current as a potential target for antiarrhythmic therapies in LQT2.

Transient Outward K+ Current (Ito) Underlies the Right Ventricular Initiation of Polymorphic Ventricular Tachycardia in a Transgenic Rabbit Model of Long-QT Syndrome Type 1

BACKGROUND

Sudden death in long-QT syndrome type 1 (LQT1), an inherited disease caused by loss-of-function mutations in KCNQ1, is triggered by early afterdepolarizations (EADs) that initiate polymorphic ventricular tachycardia (pVT). We investigated ionic mechanisms that underlie pVT in LQT1 using a transgenic rabbit model of LQT1.

METHODS

Optical mapping, cellular patch clamping, and computer modeling were used to elucidate the mechanisms of EADs in transgenic LQT1 rabbits.

RESULTS

The results showed that shorter action potential duration in the right ventricle (RV) was associated with focal activity during pVT initiation. RV cardiomyocytes demonstrated higher incidence of EADs under 50 nmol/L isoproterenol. Voltage-clamp studies revealed that the transient outward potassium current (Ito) magnitude was 28% greater in RV associated with KChiP2 but with no differences in terms of calcium-cycling kinetics and other sarcolemmal currents. Perfusing with the Ito blocker 4-aminopyridine changed the initial focal sites of pVT from the RV to the left ventricle, corroborating the role of Ito in pVT initiation. Computer modeling showed that EADs occur preferentially in the RV because of the larger conductance of the slow-inactivating component of Ito, which repolarizes the membrane potential sufficiently rapidly to allow reactivation of ICa,L before IKr has had sufficient time to activate.

CONCLUSIONS

Ito heterogeneity creates both triggers and an arrhythmogenic substrate in LQT1. In the absence of IKs, Ito interactions with ICa,L and IKr promote EADs in the RV while prolonging action potential duration in the left ventricle. This heterogeneity of action potential enhances dispersion of refractoriness and facilitates conduction blocks that initiate pVTs.

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.

Memory-induced nonlinear dynamics of excitation in cardiac diseases

Excitable cells, such as cardiac myocytes, exhibit short-term memory, i.e., the state of the cell depends on its history of excitation. Memory can originate from slow recovery of membrane ion channels or from accumulation of intracellular ion concentrations, such as calcium ion or sodium ion concentration accumulation. Here we examine the effects of memory on excitation dynamics in cardiac myocytes under two diseased conditions, early repolarization and reduced repolarization reserve, each with memory from two different sources: slow recovery of a potassium ion channel and slow accumulation of the intracellular calcium ion concentration. We first carry out computer simulations of action potential models described by differential equations to demonstrate complex excitation dynamics, such as chaos. We then develop iterated map models that incorporate memory, which accurately capture the complex excitation dynamics and bifurcations of the action potential models. Finally, we carry out theoretical analyses of the iterated map models to reveal the underlying mechanisms of memory-induced nonlinear dynamics. Our study demonstrates that the memory effect can be unmasked or greatly exacerbated under certain diseased conditions, which promotes complex excitation dynamics, such as chaos. The iterated map models reveal that memory converts a monotonic iterated map function into a nonmonotonic one to promote the bifurcations leading to high periodicity and chaos.

Effects of early afterdepolarizations on excitation patterns in an accurate model of the human ventricles

Early Afterdepolarizations, EADs, are defined as the reversal of the action potential before completion of the repolarization phase, which can result in ectopic beats. However, the series of mechanisms of EADs leading to these ectopic beats and related cardiac arrhythmias are not well understood. Therefore, we aimed to investigate the influence of this single cell behavior on the whole heart level. For this study we used a modified version of the Ten Tusscher-Panfilov model of human ventricular cells (TP06) which we implemented in a 3D ventricle model including realistic fiber orientations. To increase the likelihood of EAD formation at the single cell level, we reduced the repolarization reserve (RR) by reducing the rapid delayed rectifier Potassium current and raising the L-type Calcium current. Varying these parameters defined a 2D parametric space where different excitation patterns could be classified. Depending on the initial conditions, by either exciting the ventricles with a spiral formation or burst pacing protocol, we found multiple different spatio-temporal excitation patterns. The spiral formation protocol resulted in the categorization of a stable spiral (S), a meandering spiral (MS), a spiral break-up regime (SB), spiral fibrillation type B (B), spiral fibrillation type A (A) and an oscillatory excitation type (O). The last three patterns are a 3D generalization of previously found patterns in 2D. First, the spiral fibrillation type B showed waves determined by a chaotic bi-excitable regime, i.e. mediated by both Sodium and Calcium waves at the same time and in same tissue settings. In the parameter region governed by the B pattern, single cells were able to repolarize completely and different (spiral) waves chaotically burst into each other without finishing a 360 degree rotation. Second, spiral fibrillation type A patterns consisted of multiple small rotating spirals. Single cells failed to repolarize to the resting membrane potential hence prohibiting the Sodium channel gates to recover. Accordingly, we found that Calcium waves mediated these patterns. Third, a further reduction of the RR resulted in a more exotic parameter regime whereby the individual cells behaved independently as oscillators. The patterns arose due to a phase-shift of different oscillators as disconnection of the cells resulted in continuation of the patterns. For all patterns, we computed realistic 9 lead ECGs by including a torso model. The B and A type pattern exposed the behavior of Ventricular Tachycardia (VT). We conclude that EADs at the single cell level can result in different types of cardiac fibrillation at the tissue and 3D ventricle level.

Cardiac action potential repolarization revisited: early repolarization shows all-or-none behaviour

In healthy mammalian hearts the action potential (AP) waveform initiates and modulates each contraction, or heartbeat. As a result, AP height and duration are key physiological variables. In addition, rate-dependent changes in ventricular AP duration (APD), and variations in APD at a fixed heart rate are both reliable biomarkers of electrophysiological stability. Present guidelines for the likelihood that candidate drugs will increase arrhythmias rely on small changes in APD and Q-T intervals as criteria for safety pharmacology decisions. However, both of these measurements correspond to the final repolarization of the AP. Emerging clinical evidence draws attention to the early repolarization phase of the action potential (and the J-wave of the ECG) as an additional important biomarker for arrhythmogenesis. Here we provide a mechanistic background to this early repolarization syndrome by summarizing the evidence that both the initial depolarization and repolarization phases of the cardiac action potential can exhibit distinct time- and voltage-dependent thresholds, and also demonstrating that both can show regenerative all-or-none behaviour. An important consequence of this is that not all of the dynamics of action potential repolarization in human ventricle can be captured by data from single myocytes when these results are expressed as 'repolarization reserve'. For example, the complex pattern of cell-to-cell current flow that is responsible for AP conduction (propagation) within the mammalian myocardium can change APD and the Q-T interval of the electrocardiogram alter APD stability, and modulate responsiveness to pharmacological agents (such as Class III anti-arrhythmic drugs).

Hysteretic Dynamics of Multi-Stable Early Afterdepolarisations with Repolarisation Reserve Attenuation: A Potential Dynamical Mechanism for Cardiac Arrhythmias

Some cardiovascular and non-cardiovascular drugs frequently cause excessive prolongation of the cardiac action potential (AP) and lead to the development of early afterdepolarisations (EADs), which trigger lethal ventricular arrhythmias. Combining computer simulations in APs with numerical calculations based on dynamical system theory, we investigated stability changes of APs observed in a paced human ventricular myocyte model by decreasing and/or increasing the rapid (I_Kr) and slow (I_Ks) components of delayed rectifying K⁺ current. Upon reducing I_Kr, the APs without EADs (no-EAD response) showed gradual prolongation of AP duration (APD), and were annihilated without AP configuration changes due to the occurrence of saddle-node bifurcations. This annihilation caused a transition to an AP with EADs as a new stable steady state. Furthermore, reducing repolarisation currents (repolarisation reserve attenuation) evoked multi-stable states consisting of APs with different APDs, and caused multiple hysteretic dynamics. Depending on initial ion circumstances within ventricular myocytes, these multi-stable AP states might increase the local/global heterogeneity of AP repolarisations in the ventricle. Thus, the EAD-induced arrhythmias with repolarisation reserve attenuation might be attributed to the APD variability caused by multi-stability in cardiac AP dynamics.

Optimization of an In silico Cardiac Cell Model for Proarrhythmia Risk Assessment

Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive in vitro Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using in silico models and in vitro multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model). We also developed a novel metric that could separate drugs with different TdP liabilities at high concentrations based on total electronic charge carried by the major inward ionic currents during the action potential. In this study, we further optimized the IKr-dynamic ORd model by refining model parameters using published human cardiomyocyte experimental data under control and drug block conditions. Using this optimized model and manual patch clamp data, we developed an updated version of the metric that quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state action potential, which could classify the level of drug-induced TdP risk across a wide range of concentrations and pacing rates. We also established a framework to quantitatively evaluate a system's robustness against the induction of early afterdepolarizations (EADs), and demonstrated that the new metric is correlated with the cell's robustness to the pro-EAD perturbation of IKr conductance reduction. In summary, in this work we present an optimized model that is more consistent with experimental data, an improved metric that can classify drugs at concentrations both near and higher than clinical exposure, and a physiological framework to check the relationship between a metric and EAD. These findings provide a solid foundation for using in silico models for the regulatory assessment of TdP risk under the CiPA paradigm.

Quantitative Comparison of Effects of Dofetilide, Sotalol, Quinidine, and Verapamil between Human Ex vivo Trabeculae and In silico Ventricular Models Incorporating Inter-Individual Action Potential Variability

Background:In silico modeling could soon become a mainstream method of pro-arrhythmic risk assessment in drug development. However, a lack of human-specific data and appropriate modeling techniques has previously prevented quantitative comparison of drug effects between in silico models and recordings from human cardiac preparations. Here, we directly compare changes in repolarization biomarkers caused by dofetilide, dl-sotalol, quinidine, and verapamil, between in silico populations of human ventricular cell models and ex vivo human ventricular trabeculae. Methods and Results:Ex vivo recordings from human ventricular trabeculae in control conditions were used to develop populations of in silico human ventricular cell models that integrated intra- and inter-individual variability in action potential (AP) biomarker values. Models were based on the O'Hara-Rudy ventricular cardiomyocyte model, but integrated experimental AP variability through variation in underlying ionic conductances. Changes to AP duration, triangulation and early after-depolarization occurrence from application of the four drugs at multiple concentrations and pacing frequencies were compared between simulations and experiments. To assess the impact of variability in IC50 measurements, and the effects of including state-dependent drug binding dynamics, each drug simulation was repeated with two different IC50 datasets, and with both the original O'Hara-Rudy hERG model and a recently published state-dependent model of hERG and hERG block. For the selective hERG blockers dofetilide and sotalol, simulation predictions of AP prolongation and repolarization abnormality occurrence showed overall good agreement with experiments. However, for multichannel blockers quinidine and verapamil, simulations were not in agreement with experiments across all IC50 datasets and IKr block models tested. Quinidine simulations resulted in overprolonged APs and high incidence of repolarization abnormalities, which were not observed in experiments. Verapamil simulations showed substantial AP prolongation while experiments showed mild AP shortening. Conclusions: Results for dofetilide and sotalol show good agreement between experiments and simulations for selective compounds, however lack of agreement from simulations of quinidine and verapamil suggest further work is needed to understand the more complex electrophysiological effects of these multichannel blocking drugs.

Homeostatic enhancement of sensory transduction

Our sense of hearing boasts exquisite sensitivity, precise frequency discrimination, and a broad dynamic range. Experiments and modeling imply, however, that the auditory system achieves this performance for only a narrow range of parameter values. Small changes in these values could compromise hair cells' ability to detect stimuli. We propose that, rather than exerting tight control over parameters, the auditory system uses a homeostatic mechanism that increases the robustness of its operation to variation in parameter values. To slowly adjust the response to sinusoidal stimulation, the homeostatic mechanism feeds back a rectified version of the hair bundle's displacement to its adaptation process. When homeostasis is enforced, the range of parameter values for which the sensitivity, tuning sharpness, and dynamic range exceed specified thresholds can increase by more than an order of magnitude. Signatures in the hair cell's behavior provide a means to determine through experiment whether such a mechanism operates in the auditory system. Robustness of function through homeostasis may be ensured in any system through mechanisms similar to those that we describe here.

Role of abnormal repolarization in the mechanism of cardiac arrhythmia

In cardiac patients, life-threatening tachyarrhythmia is often precipitated by abnormal changes in ventricular repolarization and refractoriness. Repolarization abnormalities typically evolve as a consequence of impaired function of outward K⁺ currents in cardiac myocytes, which may be caused by genetic defects or result from various acquired pathophysiological conditions, including electrical remodelling in cardiac disease, ion channel modulation by clinically used pharmacological agents, and systemic electrolyte disorders seen in heart failure, such as hypokalaemia. Cardiac electrical instability attributed to abnormal repolarization relies on the complex interplay between a provocative arrhythmic trigger and vulnerable arrhythmic substrate, with a central role played by the excessive prolongation of ventricular action potential duration, impaired intracellular Ca²⁺ handling, and slowed impulse conduction. This review outlines the electrical activity of ventricular myocytes in normal conditions and cardiac disease, describes classical electrophysiological mechanisms of cardiac arrhythmia, and provides an update on repolarization-related surrogates currently used to assess arrhythmic propensity, including spatial dispersion of repolarization, activation-repolarization coupling, electrical restitution, TRIaD (triangulation, reverse use dependence, instability, and dispersion), and the electromechanical window. This is followed by a discussion of the mechanisms that account for the dependence of arrhythmic vulnerability on the location of the ventricular pacing site. Finally, the review clarifies the electrophysiological basis for cardiac arrhythmia produced by hypokalaemia, and gives insight into the clinical importance and pathophysiology of drug-induced arrhythmia, with particular focus on class Ia (quinidine, procainamide) and Ic (flecainide) Na⁺ channel blockers, and class III antiarrhythmic agents that block the delayed rectifier K⁺ channel (dofetilide).

Ca2+-activated Cl- current is antiarrhythmic by reducing both spatial and temporal heterogeneity of cardiac repolarization

The role of Ca²⁺-activated Cl^- current (I_Cl(Ca)) in cardiac arrhythmias is still controversial. It can generate delayed afterdepolarizations in Ca²⁺-overloaded cells while in other studies incidence of early afterdepolarization (EAD) was reduced by I_Cl(Ca). Therefore our goal was to examine the role of I_Cl(Ca) in spatial and temporal heterogeneity of cardiac repolarization and EAD formation. Experiments were performed on isolated canine cardiomyocytes originating from various regions of the left ventricle; subepicardial, midmyocardial and subendocardial cells, as well as apical and basal cells of the midmyocardium. I_Cl(Ca) was blocked by 0.5mmol/L 9-anthracene carboxylic acid (9-AC). Action potential (AP) changes were tested with sharp microelectrode recording. Whole-cell 9-AC-sensitive current was measured with either square pulse voltage-clamp or AP voltage-clamp (APVC). Protein expression of TMEM16A and Bestrophin-3, ion channel proteins mediating I_Cl(Ca), was detected by Western blot. 9-AC reduced phase-1 repolarization in every tested cell. 9-AC also increased AP duration in a reverse rate-dependent manner in all cell types except for subepicardial cells. Neither I_Cl(Ca) density recorded with square pulses nor the normalized expressions of TMEM16A and Bestrophin-3 proteins differed significantly among the examined groups of cells. The early outward component of I_Cl(Ca) was significantly larger in subepicardial than in subendocardial cells in APVC setting. Applying a typical subepicardial AP as a command pulse resulted in a significantly larger early outward component in both subepicardial and subendocardial cells, compared to experiments when a typical subendocardial AP was applied. Inhibiting I_Cl(Ca) by 9-AC generated EADs at low stimulation rates and their incidence increased upon beta-adrenergic stimulation. 9-AC increased the short-term variability of repolarization also. We suggest a protective role for I_Cl(Ca) against risk of arrhythmias by reducing spatial and temporal heterogeneity of cardiac repolarization and EAD formation.

Spatial Patterns of Excitation at Tissue and Whole Organ Level Due to Early Afterdepolarizations

Early after depolarizations (EAD) occur in many pathological conditions, such as congenital or acquired channelopathies, drug induced arrhythmias, and several other situations that are associated with increased arrhythmogenicity. In this paper we present an overview of the relevant computational studies on spatial EAD dynamics in 1D, 2D, and in 3D anatomical models and discuss the relation of EADs to cardiac arrhythmias. We also discuss unsolved problems and highlight new lines of research in this area.

In silico assessment of genetic variation in KCNA5 reveals multiple mechanisms of human atrial arrhythmogenesis

A recent experimental study investigating patients with lone atrial fibrillation identified six novel mutations in the KCNA5 gene. The mutants exhibited both gain- and loss-of-function of the atrial specific ultra-rapid delayed rectifier K⁺ current, I_Kur. The aim of this study is to elucidate and quantify the functional impact of these KCNA5 mutations on atrial electrical activity. A multi-scale model of the human atria was updated to incorporate detailed experimental data on I_Kur from both wild-type and mutants. The effects of the mutations on human atrial action potential and rate dependence were investigated at the cellular level. In tissue, we assessed the effects of the mutations on the vulnerability to unidirectional conduction patterns and dynamics of re-entrant excitation waves. Gain-of-function mutations shortened the action potential duration in single cells, and stabilised and accelerated re-entrant excitation in tissue. Loss-of-function mutations had heterogeneous effects on action potential duration and promoted early-after-depolarisations following beta-adrenergic stimulation. In the tissue model, loss-of-function mutations facilitated breakdown of excitation waves at more physiological excitation rates than the wild-type, and the generation of early-after-depolarisations promoted unidirectional patterns of excitation. Gain- and loss-of-function I_Kur mutations produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups of mutations. This study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

In a recent study, six mutations resulting in either gain-of-function or loss-of-function in the ultra-rapid delayed rectifier potassium current I_Kur, were identified to be associated with atrial fibrillation (AF). However, the causative link between the mutant I_Kur (either gain- or loss-of-function) and AF genesis, especially the difference and similarity between the two mutant groups, has not been elucidated. In our study, we used multiscale computational models to investigate the mechanism of arrhythmogenesis mediated by the two groups of mutations. The results suggest that the gain-of-function mutations shortened atrial action potential duration, stabilised and accelerated re-entrant excitation waves in tissue; the loss-of-function mutation promoted early-after-depolarisations following beta-adrenergic stimulation and thus wave breaks in tissue. We show these two groups of mutants carrying I_Kur produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups. Our study provides new insights into understanding the mechanisms by which mutant I_Kur contributes to atrial arrhythmias. In addition, as I_Kur is an atrial-specific channel and a number of I_Kur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking I_Kur.

Improving the In Silico Assessment of Proarrhythmia Risk by Combining hERG (Human Ether-à-go-go-Related Gene) Channel-Drug Binding Kinetics and Multichannel Pharmacology

BACKGROUND

The current proarrhythmia safety testing paradigm, although highly efficient in preventing new torsadogenic drugs from entering the market, has important limitations that can restrict the development and use of valuable new therapeutics. The CiPA (Comprehensive in vitro Proarrhythmia Assay) proposes to overcome these limitations by evaluating drug effects on multiple cardiac ion channels in vitro and using these data in a predictive in silico model of the adult human ventricular myocyte. A set of drugs with known clinical torsade de pointes risk was selected to develop and calibrate the in silico model.

METHODS AND RESULTS

Manual patch-clamp data assessing drug effects on expressed cardiac ion channels were integrated into the O'Hara-Rudy myocyte model modified to include dynamic drug-hERG channel (human Ether-à-go-go-Related Gene) interactions. Together with multichannel pharmacology data, this model predicts that compounds with high torsadogenic risk are more likely to be trapped within the hERG channel and show stronger reverse use dependency of action potential prolongation. Furthermore, drug-induced changes in the amount of electronic charge carried by the late sodium and L-type calcium currents was evaluated as a potential metric for assigning torsadogenic risk.

CONCLUSIONS

Modeling dynamic drug-hERG channel interactions and multi-ion channel pharmacology improves the prediction of torsadogenic risk. With further development, these methods have the potential to improve the regulatory assessment of drug safety models under the CiPA paradigm.

Age-dependent electrical and morphological remodeling of the Drosophila heart caused by hERG/seizure mutations

Understanding the cellular-molecular substrates of heart disease is key to the development of cardiac specific therapies and to the prevention of off-target effects by non-cardiac targeted drugs. One of the primary targets for therapeutic intervention has been the human ether a go-go (hERG) K⁺ channel that, together with the KCNQ channel, controls the rate and efficiency of repolarization in human myocardial cells. Neither of these channels plays a major role in adult mouse heart function; however, we show here that the hERG homolog seizure (sei), along with KCNQ, both contribute significantly to adult heart function as they do in humans. In Drosophila, mutations in or cardiac knockdown of sei channels cause arrhythmias that become progressively more severe with age. Intracellular recordings of semi-intact heart preparations revealed that these perturbations also cause electrical remodeling that is reminiscent of the early afterdepolarizations seen in human myocardial cells defective in these channels. In contrast to KCNQ, however, mutations in sei also cause extensive structural remodeling of the myofibrillar organization, which suggests that hERG channel function has a novel link to sarcomeric and myofibrillar integrity. We conclude that deficiency of ion channels with similar electrical functions in cardiomyocytes can lead to different types or extents of electrical and/or structural remodeling impacting cardiac output.

We have used the fruit fly cardiac model to show that seizure, the fly homolog of the human ether a go-go K⁺ channel hERG, is functional in the fly heart. This channel plays a major role in cardiac repolarization in humans but not in adult rodent hearts. Loss of channel function in the fly causes bradycardia, electrical arrhythmia and altered myofibrillar structure. Gene expression analysis indicates that Wnt signaling is affected and we show a genetic interaction between sei and pygopus, a Wnt pathway component, on heart function.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Enhanced Late Na and Ca Currents as Effective Antiarrhythmic Drug Targets

While recent advances clarified the molecular and cellular modes of action of antiarrhythmic drugs (AADs), their link to suppression of dynamical arrhythmia mechanisms remains only partially understood. The current classifications of AADs (Classes I, III, and IV) rely on blocking peak Na, K and L-type calcium currents (I_Ca,L), with Class II with dominant beta receptor blocking activity and Class V including drugs with diverse classes of actions. The discovery that the calcium and redox sensor, cardiac Ca/calmodulin-dependent protein kinase II (CaMKII) enhances both the late Na (I_Na-L) and the late I_Ca,L in patients at high risk of VT/VF provided a new and a rational AAD target. Pathological rise of either or both of I_Na-L and late I_Ca,L are demonstrated to promote cellular early afterdepolarizations (EADs) and EAD-mediated triggered activity that can initiate VT/VF in remodeled hearts. Selective inhibition of the I_Na-L without affecting their peak transients with the highly specific prototype drug, GS-967 suppresses these EAD-mediated VT/VFs. As in the case of I_Na-L, selective inhibition of the late I_Ca,L without affecting its peak with the prototype drug, roscovitine suppressed oxidative EAD-mediated VT/VF. These findings indicate that specific blockers of the late inward currents without affecting their peaks (gating modifiers), offer a new and effective AAD class action i.e., “Class VI.” The development of safe drugs with selective Class VI actions provides a rational and effective approach to treat VT/VF particularly in cardiac conditions associated with enhanced CaMKII activity such as heart failure.

Differential roles of two delayed rectifier potassium currents in regulation of ventricular action potential duration and arrhythmia susceptibility

KEY POINTS

Arrhythmias result from disruptions to cardiac electrical activity, although the factors that control cellular action potentials are incompletely understood. We combined mathematical modelling with experiments in heart cells from guinea pigs to determine how cellular electrical activity is regulated. A mismatch between modelling predictions and the experimental results allowed us to construct an improved, more predictive mathematical model. The balance between two particular potassium currents dictates how heart cells respond to perturbations and their susceptibility to arrhythmias.

ABSTRACT

Imbalances of ionic currents can destabilize the cardiac action potential and potentially trigger lethal cardiac arrhythmias. In the present study, we combined mathematical modelling with information-rich dynamic clamp experiments to determine the regulation of action potential morphology in guinea pig ventricular myocytes. Parameter sensitivity analysis was used to predict how changes in ionic currents alter action potential duration, and these were tested experimentally using dynamic clamp, a technique that allows for multiple perturbations to be tested in each cell. Surprisingly, we found that a leading mathematical model, developed with traditional approaches, systematically underestimated experimental responses to dynamic clamp perturbations. We then re-parameterized the model using a genetic algorithm, which allowed us to estimate ionic current levels in each of the cells studied. This unbiased model adjustment consistently predicted an increase in the rapid delayed rectifier K⁺ current and a drastic decrease in the slow delayed rectifier K⁺ current, and this prediction was validated experimentally. Subsequent simulations with the adjusted model generated the clinically relevant prediction that the slow delayed rectifier is better able to stabilize the action potential and suppress pro-arrhythmic events than the rapid delayed rectifier. In summary, iterative coupling of simulations and experiments enabled novel insight into how the balance between cardiac K⁺ currents influences ventricular arrhythmia susceptibility.

Dynamical mechanisms of phase-2 early afterdepolarizations in human ventricular myocytes: insights from bifurcation analyses of two mathematical models

Early afterdepolarization (EAD) is known as a cause of ventricular arrhythmias in long QT syndromes. We theoretically investigated how the rapid (I_Kr) and slow (I_Ks) components of delayed-rectifier K⁺ channel currents, L-type Ca²⁺ channel current (I_CaL), Na⁺/Ca²⁺ exchanger current (I_NCX), Na⁺-K⁺ pump current (I_NaK), intracellular Ca²⁺ (Ca_i) handling via sarcoplasmic reticulum (SR), and intracellular Na⁺ concentration (Na_i) contribute to initiation, termination, and modulation of phase-2 EADs, using two human ventricular myocyte models. Bifurcation structures of dynamical behaviors in model cells were explored by calculating equilibrium points, limit cycles (LCs), and bifurcation points as functions of parameters. EADs were reproduced by numerical simulations. The results are summarized as follows: 1) decreasing I_Ks and/or I_Kr or increasing I_CaL led to EAD generation, to which mid-myocardial cell models were especially susceptible; the parameter regions of EADs overlapped the regions of stable LCs. 2) Two types of EADs (termination mechanisms), I_Ks activation-dependent and I_CaL inactivation-dependent EADs, were detected; I_Ks was not necessarily required for EAD formation. 3) Inhibiting I_NCX suppressed EADs via facilitating Ca²⁺-dependent I_CaL inactivation. 4) Ca_i dynamics (SR Ca²⁺ handling) and Na_i strongly affected bifurcations and EAD generation in model cells via modulating I_CaL, I_NCX, and I_NaK Parameter regions of EADs, often overlapping those of stable LCs, shifted depending on Ca_i and Na_i in stationary and dynamic states. 5) Bradycardia-related induction of EADs was mainly due to decreases in Na_i at lower pacing rates. This study demonstrates that bifurcation analysis allows us to understand the dynamical mechanisms of EAD formation more profoundly.

NEW & NOTEWORTHY

We investigated mechanisms of phase-2 early afterdepolarization (EAD) by bifurcation analyses of human ventricular myocyte (HVM) models. EAD formation in paced HVMs basically depended on bifurcation phenomena in non-paced HVMs, but was strongly affected by intracellular ion concentrations in stationary and dynamic states. EAD generation did not necessarily require I_Ks.

Increase in Cardiac Ischemia-Reperfusion Injuries in Opa1+/- Mouse Model

Background

Recent data suggests the involvement of mitochondrial dynamics in cardiac ischemia/reperfusion (I/R) injuries. Whilst excessive mitochondrial fission has been described as detrimental, the role of fusion proteins in this context remains uncertain.

Objectives

To investigate whether Opa1 (protein involved in mitochondrial inner-membrane fusion) deficiency affects I/R injuries.

Methods and Results

We examined mice exhibiting Opa1^delTTAG mutations (Opa1^+/-), showing 70% Opa1 protein expression in the myocardium as compared to their wild-type (WT) littermates. Cardiac left-ventricular systolic function assessed by means of echocardiography was observed to be similar in 3-month-old WT and Opa1^+/- mice. After subjection to I/R, infarct size was significantly greater in Opa1^+/- than in WTs both in vivo (43.2±4.1% vs. 28.4±3.5%, respectively; p<0.01) and ex vivo (71.1±3.2% vs. 59.6±8.5%, respectively; p<0.05). No difference was observed in the expression of other main fission/fusion protein, oxidative phosphorylation, apoptotic markers, or mitochondrial permeability transition pore (mPTP) function. Analysis of calcium transients in isolated ventricular cardiomyocytes demonstrated a lower sarcoplasmic reticulum Ca²⁺ uptake, whereas cytosolic Ca²⁺ removal from the Na⁺/Ca²⁺ exchanger (NCX) was increased, whilst SERCA2a, phospholamban, and NCX protein expression levels were unaffected in Opa1^+/- compared to WT mice. Simultaneous whole-cell patch-clamp recordings of mitochondrial Ca²⁺ movements and ventricular action potential (AP) showed impairment of dynamic mitochondrial Ca²⁺ uptake and a marked increase in the AP late repolarization phase in conjunction with greater occurrence of arrhythmia in Opa1^+/- mice.

Conclusion

Opa1 deficiency was associated with increased sensitivity to I/R, imbalance in dynamic mitochondrial Ca²⁺ uptake, and subsequent increase in NCX activity.

Spontaneous initiation of premature ventricular complexes and arrhythmias in type 2 long QT syndrome

The occurrence of early afterdepolarizations (EADs) and increased dispersion of repolarization are two known factors for arrhythmogenesis in long QT syndrome. However, increased dispersion of repolarization tends to suppress EADs due to the source-sink effect, and thus how the two competing factors cause initiation of arrhythmias remains incompletely understood. Here we used optical mapping and computer simulation to investigate the mechanisms underlying spontaneous initiation of arrhythmias in type 2 long QT (LQT2) syndrome. In optical mapping experiments of transgenic LQT2 rabbit hearts under isoproterenol, premature ventricular complexes (PVCs) were observed to originate from the steep spatial repolarization gradient (RG) regions and propagated unidirectionally. The same PVC behaviors were demonstrated in computer simulations of tissue models of rabbits. Depending on the heterogeneities, these PVCs could lead to either repetitive focal excitations or reentry without requiring an additional vulnerable substrate. Systematic simulations showed that cellular phase 2 EADs were either suppressed or confined to the long action potential region due to the source-sink effect. Tissue-scale phase 3 EADs and PVCs occurred due to tissue-scale dynamical instabilities caused by RG and enhanced L-type calcium current (I_Ca,L), occurring under both large and small RG. Presence of cellular EADs was not required but potentiated PVCs when RG was small. We also investigated how other factors affect the dynamical instabilities causing PVCs. Our main conclusion is that tissue-scale dynamical instabilities caused by RG and enhanced I_Ca,L give rise to both the trigger and the vulnerable substrate simultaneously for spontaneous initiation of arrhythmias in LQT2 syndrome.

Quinine and the ABCs of Long QT: A Patient's Misfortune with Arthritis, (Alcoholic) Beverages, and Cramps

A 91-year-old woman presented to the emergency department by ambulance after her family found her minimally responsive. Telemetry monitoring demonstrated episodes of non-sustained polymorphic ventricular tachycardia (PMVT) associated with significantly prolonged repolarization. Her medical history revealed that she was taking quinine or a derivative in three different forms: hydroxychloroquine, quinine sulfate (for leg cramps), and her gin mixed with tonic water (containing quinine). The present case is illustrative of classic etiologies and findings of acquired long QT syndrome, and serves as an important reminder for providers to take a complete medication history, including use of duplicative and alternative medicines and type of alcohol consumption.

CLOCK-BMAL1 regulate the cardiac L-type calcium channel subunit CACNA1C through PI3K-Akt signaling pathway

The heterodimerized transcription factors CLOCK-BMAL1 regulate the cardiomyocyte circadian rhythms. The L-type calcium currents play important role in the cardiac electrogenesis and arrhythmogenesis. Whether and how the CLOCK-BMAL1 regulate the cardiac L-type calcium channels are yet to be determined. The functions of the L-type calcium channels were evaluated with patch clamping techniques. Recombinant adenoviruses of CLOCK and BMAL1 were used in the expression experiments. We reported that the expressions and functions of CACNA1C (the α-subunit of the L-type calcium channels) showed circadian rhythms, with the peak at zeitgeber time 3 (ZT3). The endocardial action potential durations 90 (APD90) were correspondingly longer at ZT3. The protein levels of the phosphorylated Akt at threonine 308 (pAkt T308) also showed circadian rhythms. Overexpressions of CLOCK-BMAL1 significantly reduced the levels of CACNA1C while increasing the levels of pAkt T308 and pik3r1. Furthermore, the inhibitory effects of CLOCK-BMAL1 on CACNA1C could be abolished by the Akt inhibitor MK2206 or the PDK1 inhibitor GSK2334470. Collectively, our findings suggested that the expressions of the cardiac CACNA1C were under the CLOCK-BMAL1 regulation, probably through the PI3K-Akt signal pathway.

Effects of allocryptopine on outward potassium current and slow delayed rectifier potassium current in rabbit myocardium

Objective

Allocryptopine (ALL) is an effective alkaloid of Corydalis decumbens (Thunb.) Pers. Papaveraceae and has proved to be anti-arrhythmic. The purpose of our study is to investigate the effects of ALL on transmural repolarizing ionic ingredients of outward potassium current (I_to) and slow delayed rectifier potassium current (I_Ks).

Methods

The monophasic action potential (MAP) technique was used to record the MAP duration of the epicardium (Epi), myocardium (M) and endocardium (Endo) of the rabbit heart and the whole cell patch clamp was used to record I_to and I_Ks in cardiomyocytes of Epi, M and Endo layers that were isolated from rabbit ventricles.

Results

The effects of ALL on MAP of Epi, M and Endo layers were disequilibrium. ALL could effectively reduce the transmural dispersion of repolarization (TDR) in rabbit transmural ventricular wall. ALL decreased the current densities of I_to and I_Ks in a voltage and concentration dependent way and narrowed the repolarizing differences among three layers. The analysis of gating kinetics showed ALL accelerated the channel activation of I_to in M layers and partly inhibit the channel openings of I_to in Epi, M and Endo cells. On the other hand, ALL mainly slowed channel deactivation of I_Ks channel in Epi and Endo layers without affecting its activation.

Conclusions

Our study gives partially explanation about the mechanisms of transmural inhibition of I_to and I_Ks channels by ALL in rabbit myocardium. These findings provide novel perspective regarding the anti-arrhythmogenesis application of ALL in clinical settings.

Ventricular Arrhythmias in Apparently Normal Hearts: Who Needs an Implantable Cardiac Defibrillator?

Idiopathic ventricular tachycardia is often considered a benign form of ventricular arrhythmia in patients without apparent structural heart disease. However, a subset of patients may develop malignant ventricular arrhythmias and present with syncope and sudden cardiac arrest. Survivors of cardiac arrest are candidates for implantable cardiac defibrillators (ICDs). The indications for ICDs in patients with less than a full-blown cardiac arrest presentation but with electrocardiographically high-risk ectopy features remain uncertain. This article addresses some of the uncertainties and pitfalls in ICD risk stratification in this patient group and explores potential mechanisms for malignant conversion of benign premature ventricular complexes to sustained arrhythmia.

Use of FDSS/μCell imaging platform for preclinical cardiac electrophysiology safety screening of compounds in human induced pluripotent stem cell-derived cardiomyocytes

FDSS/μCell is a high-speed acquisition imaging platform (Hamamatsu Ltd., Hamamatsu, Japan) that allows for simultaneous high-throughput reading under controlled conditions. We evaluated the Ca(2+) transients or optical membrane potential changes of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) (iCells) in the presence or absence of 44 pharmacological agents known to interfere with cardiac ion channels (e.g., hERG, IKs, NaV1.5, CaV1.2). We tested two Ca(2+)-sensitive fluorescence dyes (Codex ACTOne® and EarlyTox®) and a membrane potential dye (FLIPR® membrane potential dye). We were able to quantify and report drug-induced early-after depolarizations (EAD)-like waveforms, cardiomyocyte ectopic beats and changes in beating rate from a subgroup of pharmacological agents acting acutely (within a 1-hour period). Cardiovascular drugs, such as dofetilide and d,l-sotalol, exhibited EAD-like signals at 3nM and 10μM, respectively. CNS drugs, such as haloperidol and sertindole, exhibited EAD-like signals and ectopic beats at 30nM and 1μM, respectively. Other drugs, such as astemizole, solifenacin, and moxifloxacin, exhibited similar arrhythmias at 30nM, 3μM and 300μM, respectively. Our data suggest that the membrane potential and intracellular Ca(2+) signal are tightly coupled, supporting the idea that the EAD-like signals reported are the accurate representation of an EAD signal of the cardiac action potential. Finally, the EAD-like Ca(2+) signal was well correlated to clinically-relevant concentrations where Torsade de Pointes (TdPs) arrhythmias were noted in healthy volunteers treated orally with some of the compounds we tested, as reported in PharmaPendium®.

Beat-to-Beat Variability in Preload Unmasks Latent Risk of Torsade de Pointes in Anesthetized Chronic Atrioventricular Block Dogs

BACKGROUND

Beat-to-beat variability in ventricular repolarization (BVR) associates with increased arrhythmic risk. Proarrhythmic remodeling in the dog with chronic AV-block (CAVB) compromises repolarization reserve and associates with increased BVR, which further increases upon dofetilide infusion and correlates with Torsade de Pointes (TdP) arrhythmias. It was hypothesized that these pro-arrhythmia-associated increases in BVR are induced by beat-to-beat variability in preload.

METHODS AND RESULTS

Left ventricular monophasic action potential duration (LVMAPD) was recorded in acute (AAVB) and CAVB dogs, before and after dofetilide infusion. BVR was quantified as short-term variability of LVMAPD. The PQ-interval was controlled by pacing: either a constant or an alternating preload pattern was established, verified by PV-loop. The effect of the stretch-activated channel blocker, streptomycin, on BVR was evaluated in a second CAVB group. At alternating preload only, BVR was increased after proarrhythmic remodeling (0.45±0.14 ms AAVB vs. 2.2±1.1 ms CAVB, P<0.01). At CAVB, but not at AAVB, dofetilide induced significant proarrhythmia. Preload variability augmented the dofetilide-induced BVR increase at CAVB (+1.5±0.8 ms vs. +0.9±0.9 ms, P=0.058). In the second group, the increase in baseline BVR by alternating preload (0.3±0.03 ms to 1.0±0.8 ms, P<0.01) was abolished by streptomycin (0.5±0.2 ms, P<0.05).

CONCLUSIONS

In CAVB dogs, the inverse relation between BVR and repolarization reserve originates from an augmented sensitivity of ventricular repolarization to beat-to-beat preload changes. Stretch-activated channels appear to be involved in the mechanism of BVR. (Circ J 2016; 80: 1336-1345).

Atrial-Selective Potassium Channel Blockers

Atrial fibrillation (AF) is associated with increased morbidity and mortality. Atrial-selective potassium (K(+)) channel blockers may represent a novel therapeutic target. The best validated atrial-specific ion currents are the acetylcholine-activated inward-rectifier K(+) current IK,ACh and ultrarapidly activating delayed-rectifier K(+) current IKur. Two-pore domain and small-conductance Ca(2+)-activated K(+) channels and Kv1.1 channels may also contribute to the atrial repolarization. We review the molecular and electrophysiologic characteristics of atrial-selective K(+) channels and their potential pathophysiologic role in AF. We summarize currently available K(+) channel blockers focusing on the most important compounds.

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

Early afterdepolarizations (EADs) are pathological oscillations in cardiac action potentials during the repolarization phase and may be caused by drug side effects, ion channel disease or oxidative stress. The most widely observed EAD pattern is characterized by oscillations with growing amplitudes. So far, its occurence has been explained in terms of a supercritical Hopf bifurcation in the fast subsystem of the action potential dynamics from which stable limit cycles with growing amplitudes emerge. The novel contribution of this article is the introduction of two alternative explanations of EAD genesis with growing amplitudes that do not involve stable limit cycles in fast subsystems. In particular, we demonstrate that EAD patterns with growing amplitudes may alternatively arise due to a delayed subcritical Hopf bifurcation or an unstable manifold of a saddle focus fixed point in the full fast-slow system modelling the action potential. Our work extends the list of possible dynamical EAD mechanisms and may contribute to a classification of drug effects in preclinical cardiotoxicity testing.

Power-Law Dynamics of Membrane Conductances Increase Spiking Diversity in a Hodgkin-Huxley Model

We studied the effects of non-Markovian power-law voltage dependent conductances on the generation of action potentials and spiking patterns in a Hodgkin-Huxley model. To implement slow-adapting power-law dynamics of the gating variables of the potassium, n, and sodium, m and h, conductances we used fractional derivatives of order η≤1. The fractional derivatives were used to solve the kinetic equations of each gate. We systematically classified the properties of each gate as a function of η. We then tested if the full model could generate action potentials with the different power-law behaving gates. Finally, we studied the patterns of action potential that emerged in each case. Our results show the model produces a wide range of action potential shapes and spiking patterns in response to constant current stimulation as a function of η. In comparison with the classical model, the action potential shapes for power-law behaving potassium conductance (n gate) showed a longer peak and shallow hyperpolarization; for power-law activation of the sodium conductance (m gate), the action potentials had a sharp rise time; and for power-law inactivation of the sodium conductance (h gate) the spikes had wider peak that for low values of η replicated pituitary- and cardiac-type action potentials. With all physiological parameters fixed a wide range of spiking patterns emerged as a function of the value of the constant input current and η, such as square wave bursting, mixed mode oscillations, and pseudo-plateau potentials. Our analyses show that the intrinsic memory trace of the fractional derivative provides a negative feedback mechanism between the voltage trace and the activity of the power-law behaving gate variable. As a consequence, power-law behaving conductances result in an increase in the number of spiking patterns a neuron can generate and, we propose, expand the computational capacity of the neuron.

There is increasing evidence that the activity of individual membrane ion channels, conductances, and the firing rate of neurons are history dependent. In this work we studied how history dependent activation of membrane conductances affect the action potential activity of the Hodgkin-Huxley model, a widely used model of action potential generation. In order to implement history dependent activation, we made use of fractional order differential equations. This type of history dependent differential equations are increasingly being used in biomedical sciences to simulate complex phenomena. We use fractional order derivatives to model the kinetic dynamics of the gate variables for the potassium and sodium conductances of the Hodgkin-Huxley model. Our results show that power-law dynamics of the different gate variables result in a wide range of action potential shapes and spiking patterns, even in the case where the model was stimulated with constant current. As a consequence, power-law behaving conductances result in an increase in the number of spiking patterns a neuron can generate and, we propose, expand the computational capacity of the neuron.