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

Optogenetic Modulation of Arrhythmia Triggers: Proof-of-Concept from Computational Modeling

Introduction

Early afterdepolarizations (EADs) are secondary voltage depolarizations associated with reduced repolarization reserve (RRR) that can trigger lethal arrhythmias. Relating EADs to triggered activity is difficult to study, so the ability to suppress or provoke EADs would be experimentally useful. Here, we use computational simulations to assess the feasibility of subthreshold optogenetic stimulation modulating the propensity for EADs (cell-scale) and EAD-associated ectopic beats (organ-scale).

Methods

We modified a ventricular ionic model by reducing rapid delayed rectifier potassium (0.25-0.1 × baseline) and increasing L-type calcium (1.0-3.5 × baseline) currents to create RRR conditions with varying severity. We ran simulations in models of single cardiomyocytes and left ventricles from post-myocardial infarction patient MRI scans. Optogenetic stimulation was simulated using either ChR2 (depolarizing) or GtACR1 (repolarizing) opsins.

Results

In cell-scale simulations without illumination, EADs were seen for 164 of 416 RRR conditions. Subthreshold stimulation of GtACR1 reduced EAD incidence by up to 84.8% (25/416 RRR conditions; 0.1 μW/mm2); in contrast, subthreshold ChR2 excitation increased EAD incidence by up to 136.6% (388/416 RRR conditions; 50 μW/mm2). At the organ scale, we assumed simultaneous, uniform illumination of the epicardial and endocardial surfaces. GtACR1-mediated suppression (10-50 μW/mm2) and ChR2-mediated unmasking (50-100 μW/mm2) of EAD-associated ectopic beats were feasible in three distinct ventricular models.

Conclusions

Our findings suggest that optogenetics could be used to silence or provoke both EADs and EAD-associated ectopic beats. Validation in animal models could lead to exciting new experimental regimes and potentially to novel anti-arrhythmia treatments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-023-00781-z.

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

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

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

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

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.

Predicting drug-induced arrhythmias by multiscale modeling

Drugs often have undesired side effects. In the heart, they can induce lethal arrhythmias such as torsades de pointes. The risk evaluation of a new compound is costly and can take a long time, which often hinders the development of new drugs. Here, we establish a high-resolution, multiscale computational model to quickly assess the cardiac toxicity of new and existing drugs. The input of the model is the drug-specific current block from single cell electrophysiology; the output is the spatio-temporal activation profile and the associated electrocardiogram. We demonstrate the potential of our model for a low-risk drug, ranolazine, and a high-risk drug, quinidine: For ranolazine, our model predicts a prolonged QT interval of 19.4% compared with baseline and a regular sinus rhythm at 60.15 beats per minute. For quinidine, our model predicts a prolonged QT interval of 78.4% and a spontaneous development of torsades de pointes both in the activation profile and in the electrocardiogram. Our model reveals the mechanisms by which electrophysiological abnormalities propagate across the spatio-temporal scales, from specific channel blockage, via altered single cell action potentials and prolonged QT intervals, to the spontaneous emergence of ventricular tachycardia in the form of torsades de pointes. Our model could have important implications for researchers, regulatory agencies, and pharmaceutical companies on rationalizing safe drug development and reducing the time-to-market of new drugs.

Early afterdepolarisation tendency as a simulated pro-arrhythmic risk indicator

Drug-induced Torsades de Pointes (TdP) arrhythmia is of major interest in predictive toxicology. Drugs which cause TdP block the hERG cardiac potassium channel. However, not all drugs that block hERG cause TdP. As such, further understanding of the mechanistic route to TdP is needed. Early afterdepolarisations (EADs) are a cell-level phenomenon in which the membrane of a cardiac cell depolarises a second time before repolarisation, and EADs are seen in hearts during TdP. Therefore, we propose a method of predicting TdP using induced EADs combined with multiple ion channel block in simulations using biophysically-based mathematical models of human ventricular cell electrophysiology. EADs were induced in cardiac action potential models using interventions based on diseases that are known to cause EADs, including: increasing the conduction of the L-type calcium channel, decreasing the conduction of the hERG channel, and shifting the inactivation curve of the fast sodium channel. The threshold of intervention that was required to cause an EAD was used to classify drugs into clinical risk categories. The metric that used L-type calcium induced EADs was the most accurate of the EAD metrics at classifying drugs into the correct risk categories, and increased in accuracy when combined with action potential duration measurements. The EAD metrics were all more accurate than hERG block alone, but not as predictive as simpler measures such as simulated action potential duration. This may be because different routes to EADs represent risk well for different patient subgroups, something that is difficult to assess at present.