Dynamics of atrial arrhythmias modulated by time-dependent acetylcholine concentration: a simulation study

25 years of basic and translational science in EP Europace: novel insights into arrhythmia mechanisms and therapeutic strategies

In the last 25 years, EP Europace has published more than 300 basic and translational science articles covering different arrhythmia types (ranging from atrial fibrillation to ventricular tachyarrhythmias), different diseases predisposing to arrhythmia formation (such as genetic arrhythmia disorders and heart failure), and different interventional and pharmacological anti-arrhythmic treatment strategies (ranging from pacing and defibrillation to different ablation approaches and novel drug-therapies). These studies have been conducted in cellular models, small and large animal models, and in the last couple of years increasingly in silico using computational approaches. In sum, these articles have contributed substantially to our pathophysiological understanding of arrhythmia mechanisms and treatment options; many of which have made their way into clinical applications. This review discusses a representative selection of EP Europace manuscripts covering the topics of pacing and ablation, atrial fibrillation, heart failure and pro-arrhythmic ventricular remodelling, ion channel (dys)function and pharmacology, inherited arrhythmia syndromes, and arrhythmogenic cardiomyopathies, highlighting some of the advances of the past 25 years. Given the increasingly recognized complexity and multidisciplinary nature of arrhythmogenesis and continued technological developments, basic and translational electrophysiological research is key advancing the field. EP Europace aims to further increase its contribution to the discovery of arrhythmia mechanisms and the implementation of mechanism-based precision therapy approaches in arrhythmia management.

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: A computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and β-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD₉₀) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD₉₀ and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 μM. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD₉₀ prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility.

In silico Comparison of Left Atrial Ablation Techniques That Target the Anatomical, Structural, and Electrical Substrates of Atrial Fibrillation

Catheter ablation therapy for persistent atrial fibrillation (AF) typically includes pulmonary vein isolation (PVI) and may include additional ablation lesions that target patient-specific anatomical, electrical, or structural features. Clinical centers employ different ablation strategies, which use imaging data together with electroanatomic mapping data, depending on data availability. The aim of this study was to compare ablation techniques across a virtual cohort of AF patients. We constructed 20 paroxysmal and 30 persistent AF patient-specific left atrial (LA) bilayer models incorporating fibrotic remodeling from late-gadolinium enhancement (LGE) MRI scans. AF was simulated and post-processed using phase mapping to determine electrical driver locations over 15 s. Six different ablation approaches were tested: (i) PVI alone, modeled as wide-area encirclement of the pulmonary veins; PVI together with: (ii) roof and inferior lines to model posterior wall box isolation; (iii) isolating the largest fibrotic area (identified by LGE-MRI); (iv) isolating all fibrotic areas; (v) isolating the largest driver hotspot region [identified as high simulated phase singularity (PS) density]; and (vi) isolating all driver hotspot regions. Ablation efficacy was assessed to predict optimal ablation therapies for individual patients. We subsequently trained a random forest classifier to predict ablation response using (a) imaging metrics alone, (b) imaging and electrical metrics, or (c) imaging, electrical, and ablation lesion metrics. The optimal ablation approach resulting in termination, or if not possible atrial tachycardia (AT), varied among the virtual patient cohort: (i) 20% PVI alone, (ii) 6% box ablation, (iii) 2% largest fibrosis area, (iv) 4% all fibrosis areas, (v) 2% largest driver hotspot, and (vi) 46% all driver hotspots. Around 20% of cases remained in AF for all ablation strategies. The addition of patient-specific and ablation pattern specific lesion metrics to the trained random forest classifier improved predictive capability from an accuracy of 0.73 to 0.83. The trained classifier results demonstrate that the surface areas of pre-ablation driver regions and of fibrotic tissue not isolated by the proposed ablation strategy are both important for predicting ablation outcome. Overall, our study demonstrates the need to select the optimal ablation strategy for each patient. It suggests that both patient-specific fibrosis properties and driver locations are important for planning ablation approaches, and the distribution of lesions is important for predicting an acute response.

Heart failure-induced atrial remodelling promotes electrical and conduction alternans

Heart failure (HF) is associated with an increased propensity for atrial fibrillation (AF), causing higher mortality than AF or HF alone. It is hypothesized that HF-induced remodelling of atrial cellular and tissue properties promotes the genesis of atrial action potential (AP) alternans and conduction alternans that perpetuate AF. However, the mechanism underlying the increased susceptibility to atrial alternans in HF remains incompletely elucidated. In this study, we investigated the effects of how HF-induced atrial cellular electrophysiological (with prolonged AP duration) and tissue structural (reduced cell-to-cell coupling caused by atrial fibrosis) remodelling can have an effect on the generation of atrial AP alternans and their conduction at the cellular and one-dimensional (1D) tissue levels. Simulation results showed that HF-induced atrial electrical remodelling prolonged AP duration, which was accompanied by an increased sarcoplasmic reticulum (SR) Ca²⁺ content and Ca²⁺ transient amplitude. Further analysis demonstrated that HF-induced atrial electrical remodelling increased susceptibility to atrial alternans mainly due to the increased sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) Ca²⁺ reuptake, modulated by increased phospholamban (PLB) phosphorylation, and the decreased transient outward K⁺ current (I_to). The underlying mechanism has been suggested that the increased SR Ca²⁺ content and prolonged AP did not fully recover to their previous levels at the end of diastole, resulting in a smaller SR Ca²⁺ release and AP in the next beat. These produced Ca²⁺ transient alternans and AP alternans, and further caused AP alternans and Ca²⁺ transient alternans through Ca²⁺→AP coupling and AP→Ca²⁺ coupling, respectively. Simulation of a 1D tissue model showed that the combined action of HF-induced ion channel remodelling and a decrease in cell-to-cell coupling due to fibrosis increased the heart tissue’s susceptibility to the formation of spatially discordant alternans, resulting in an increased functional AP propagation dispersion, which is pro-arrhythmic. These findings provide insights into how HF promotes atrial arrhythmia in association with atrial alternans.

Atrial Fibrillation (AF) is the most common arrhythmia in adults, especially in the elderly, with the increased incidence of stroke being a major complication that increases morbidity and mortality. The occurrence of AF is often accompanied by heart failure (HF). AF and HF are also known to have the bidirectional relationship that AF worsens HF and HF promotes AF. HF can induce atrial remodelling, including electrical remodelling, atrial fibrosis, stretch and dilatation, and oxidative stress, in which many factors are associated with arrhythmogenic atrial alternans. HF-induced atrial remodelling varies during various stages and complications of HF, but possible mechanisms underlying their pro-susceptibility to alternans have not been completely elucidated. In this study, we investigated the effects of HF-induced atrial remodelling with prolonged action potential duration (APD) and decreased cell-to-cell coupling on susceptibility to atrial alternans. Simulation results showed that HF-induced an increase in sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) Ca²⁺ reuptake caused by increased phospholamban phosphorylation and a decrease in transient outward K⁺ current played significant roles in the genesis of Ca²⁺ transient alternans and action potential alternans at the single-cell level. The HF-induced decline of cell-to-cell coupling and APD prolongation promoted the genesis of spatially discordant alternans in atrial tissue. This provides insights into how HF facilitates atrial arrhythmia in relation to atrial alternans.

Acetylcholine Delays Atrial Activation to Facilitate Atrial Fibrillation

Background

Acetylcholine (ACh) shortens action potential duration (APD) in human atria. APD shortening facilitates atrial fibrillation (AF) by reducing the wavelength for reentry. However, the influence of ACh on electrical conduction in human atria and its contribution to AF are unclear, particularly when combined with impaired conduction from interstitial fibrosis.

Objective

To investigate the effect of ACh on human atrial conduction and its role in AF with computational, experimental, and clinical approaches.

Methods

S1S2 pacing (S1 = 600 ms and S2 = variable cycle lengths) was applied to the following human AF computer models: a left atrial appendage (LAA) myocyte to quantify the effects of ACh on APD, maximum upstroke velocity (V _max ), and resting membrane potential (RMP); a monolayer of LAA myocytes to quantify the effects of ACh on conduction; and 3) an intact left atrium (LA) to determine the effects of ACh on arrhythmogenicity. Heterogeneous ACh and interstitial fibrosis were applied to the monolayer and LA models. To corroborate the simulations, APD and RMP from isolated human atrial myocytes were recorded before and after 0.1 μM ACh. At the tissue level, LAAs from AF patients were optically mapped ex vivo using Di-4-ANEPPS. The difference in total activation time (AT) was determined between AT initially recorded with S1 pacing, and AT recorded during subsequent S1 pacing without (n = 6) or with (n = 7) 100 μM ACh.

Results

In LAA myocyte simulations, S1 pacing with 0.1 μM ACh shortened APD by 41 ms, hyperpolarized RMP by 7 mV, and increased V _max by 27 mV/ms. In human atrial myocytes, 0.1 μM ACh shortened APD by 48 ms, hyperpolarized RMP by 3 mV, and increased V _max by 6 mV/ms. In LAA monolayer simulations, S1 pacing with ACh hyperpolarized RMP to delay total AT by 32 ms without and 35 ms with fibrosis. This led to unidirectional conduction block and sustained reentry in fibrotic LA with heterogeneous ACh during S2 pacing. In AF patient LAAs, S1 pacing with ACh increased total AT from 39.3 ± 26 ms to 71.4 ± 31.2 ms (p = 0.036) compared to no change without ACh (56.7 ± 29.3 ms to 50.0 ± 21.9 ms, p = 0.140).

Conclusion

In fibrotic atria with heterogeneous parasympathetic activation, ACh facilitates AF by shortening APD and slowing conduction to promote unidirectional conduction block and reentry.

Fitting local repolarization parameters in cardiac reaction-diffusion models in the presence of electrotonic coupling

BACKGROUND

Repolarization gradients contribute to arrhythmogenicity. In reaction-diffusion models of cardiac tissue, heterogeneities in action potential duration (APD) can be created by locally modifying an intrinsic membrane kinetics parameter. Electrotonic coupling, however, acts as a confounding factor that modulates APD dispersion.

METHOD

We developed an algorithm based on a quasi-Newton method that iteratively adjusts the spatial distribution of a membrane parameter to reproduce a pre-defined target APD map in a coupled tissue. The method assumes that the relation between the adjustable parameter and APD is bijective in an isolated cell. Each iteration of the algorithm involved simulating the cardiac reaction-diffusion system with the updated parameter profile for one beat and extracting the APD map. The algorithm was extended to simultaneous estimation of two parameter profiles based on two APD maps at different repolarization thresholds.

RESULTS

The method was validated in 1D, 2D and 3D atrial tissues using synthetic target APD maps with controllable total variation and maximum APD gradient. The adjustable parameter was local acetylcholine concentration. The iterations converged provided that APD gradients were not too steep. Convergence was found to be faster 2-5 iterations) when the maximal gradient was less steep, when APD range was smaller and when tissue conductivity was reduced.

CONCLUSION

This algorithm provides a tool to automatically generate arrhythmogenic substrates with controllable repolarization gradients and possibly incorporate experimental APD maps into computer models.

Lessons from computer simulations of ablation of atrial fibrillation

This paper reviews the simulations of catheter ablation in computer models of the atria, from the first attempts to the most recent anatomical models. It describes how postulated substrates of atrial fibrillation can be incorporated into mathematical models, how modelling studies can be designed to test ablation strategies, what their current trade-offs and limitations are, and what clinically relevant lessons can be learnt from these simulations. Drawing a parallel between clinical and modelling studies, six ablation targets are considered: pulmonary vein isolation, linear ablation, ectopic foci, complex fractionated atrial electrogram, rotors and ganglionated plexi. The examples presented for each ablation target illustrate a major advantage of computer models, the ability to identify why a therapy is successful or not in a given atrial fibrillation substrate. The integration of pathophysiological data to create detailed models of arrhythmogenic substrates is expected to solidify the understanding of ablation mechanisms and to provide theoretical arguments supporting substrate-specific ablation strategies.

A simple analytical model of action potential duration profile in electrotonically-coupled cells

Electrotonic interactions between cardiac cells modulate the dispersion of action potential duration (APD). This paper provides a complete mathematical analysis of a simple model of exponential-shaped repolarization in a network of electrotonically-coupled cells with different intrinsic APDs. The forward problem consists in computing the APD map in the coupled system from the intrinsic APD map. A closed-form algebraic formula is derived for the forward problem. The inverse problem, inferring the intrinsic APDs from an APD map, is proved to have a unique solution (if any). Perturbation analysis leads to an efficient and accurate Newton-based solver for this specific inverse problem. Finally, an analytical expression is obtained for the convolution filter that solves the forward problem in one dimension. This mathematical framework forms a solid theoretical basis for future development and validation of repolarization parameter estimation techniques in detailed models of cardiac tissue.