Re-entrant Excitation Waves Research Articles

PO-01-180 MECHANISMS OF PULMONARY ARTERIAL HYPERTENSION-INDUCED ATRIAL FIBRILLATION IN THE RIGHT HEART: INSIGHTS FROM MULTI-SCALE MODELS OF THE HUMAN ATRIA

Pulmonary arterial hypertension (PAH) is one of the common conditions primarily affecting the right heart and promoting atrial fibrillation (AF), but the mechanisms are poorly understood. This study aimed to use multi-scale atrial models to investigate mechanisms by which right atrial (RA) remodeling facilitates and perpetuates atrial arrhythmias. We have developed single-cell-to-tissue-level computer models to mimic control and disease conditions by incorporating the most recent experimental data on PAH-induced electrical and structural remodeling (Hiram et al., JACC, 2019). The effects of PAH-induced remodeling in both RA and left atria (LA) on action potential duration (APD), calcium transient (CaT), conduction velocity (CV), and dynamics of re-entrant excitation waves were assessed using the multi-scale computer models. Our computer simulation results suggested that at the single cell level, although PAH-induced remodeling prolonged action potential (ΔAPD: 49.6 ms in the RA vs. 41.6 ms in the LA) and increased calcium transient (ΔCaT: 7.5e-2 μM in the RA vs. 0.9e-3 μM in the LA), heterogeneous remodeling increased susceptibility to afterdepolarizations only in the RA (Fig. A). At the tissue level, dramatically reduced CV (ΔCV: -0.5 m/s in the RA vs. -0.05 m/s in the LA) abbreviated wavelength inthe RA but not inthe LA. In addition, afterdepolarizations in the RA increased repolarization dispersion and promoted unidirectional conduction block (Fig B). Finally, the upregulated fibrosis increased the chance of breakdown of excitation waves and sustained reentries in in the RA (Fig C). These results demonstrated that PAH-induced AF is mainly in RA not LA, due to a higher degree of cardiac remodeling. RA remodeling increases arrhythmias risk due to elevated susceptibility to afterdepolarizations, reduced wavelength, increased repolarization heterogeneity and upregulation of fibrosis, which together facilitate initiation and maintenance of re-entrant circuits in the right heart.

Mechanistic insights into spontaneous transition from cellular alternans to ventricular fibrillation.

T-wave alternans (TWA) has been used for predicting the risk of malignant cardiac arrhythmias and sudden cardiac death (SCD) in multiple clinical settings; however, possible mechanism(s) underlying the spontaneous transition from cellular alternans reflected by TWA to arrhythmias in impaired repolarization remains unclear. The healthy guinea pig ventricular myocytes under E-4031 blocking IKr (0.1μM, N=12; 0.3μM, N=10; 1μM, N=10) were evaluated using whole-cell patch-clamp. The electrophysiological properties of isolated perfused guinea pig hearts under E-4031 (0.1μM, N=5; 0.3μM, N=5; 1μM, N=5) were evaluated using dual- optical mapping. The amplitude/threshold/restitution curves of action potential duration (APD) alternans and potential mechanism(s) underlying the spontaneous transition of cellular alternans to ventricular fibrillation (VF) were examined. There were longer APD80 and increased amplitude and threshold of APD alternans in E-4031 group compared with baseline group, which was reflected by more pronounced arrhythmogenesis at the tissue level, and were associated with steep restitution curves of the APD and the conduction velocity (CV). Conduction of AP alternans augmented tissue's functional spatiotemporal heterogeneity of regional AP/Ca alternans, as well as the AP/Ca dispersion, leading to localized uni-directional conduction block that spontaneous facilitated the formation of reentrant excitation waves without the need for additional premature stimulus. Our results provide a possible mechanism for the spontaneous transition from cardiac electrical alternans in cellular action potentials and intercellular conduction without the involvement of premature excitations, and explain the increased susceptibility to ventricular arrhythmias in impaired repolarization. In this study, we implemented voltage-clamp and dual-optical mapping approaches to investigate the underlying mechanism(s) for the arrhythmogenesis of cardiac alternans in the guinea pig heart at cellular and tissue levels. Our results demonstrated a spontaneous development of reentry from cellular alternans, arising from a combined actions of restitution properties of action potential duration, conduction velocity of excitation wave and interplay between alternants of action potential and the intracellular Ca handling. We believe this study provides new insights into underlying the mechanism, by which cellular cardiac alternans spontaneously evolves into cardiac arrhythmias.

A Multi-Scale Computational Model for the Rat Ventricle: Construction, Parallelization, and Applications

BackgroundCardiovascular diseases are the top killer of human beings. The ventricular arrhythmia, as a type of malignant cardiac arrhythmias, typically leads to death if not treated within minutes. The multi-scale virtual heart provides an idealized tool for exploring the underlying mechanisms, by means of incorporating abundant experimental data at the level of ion channels and analyzing the subsequent pathological changes at organ levels. However, there are few studies on building a virtual heart model for rats—a species most widely used in experiments. ObjectiveTo build a multi-scale computational model for rats, with detailed methodology for the model construction, computational optimization, and its applications. MethodsFirst, approaches for building multi-scale models ranging from cellular to 3-D organ levels are introduced, with detailed descriptions of handling the ventricular myocardium heterogeneity, geometry processing, and boundary conditions, etc. Next, for dealing with the expensive computational costs of 3-D models, optimization approaches including an optimized representation and a GPU-based parallelization method are introduced. Finally, methods for reproducing of some key phenomenon (e.g., electrocardiograph, spiral/scroll waves) are demonstrated. ResultsThree types of heterogeneity, including the transmural heterogeneity, the interventricular heterogeneity, and the base-apex heterogeneity are incorporated into the model. The normal and reentrant excitation waves, as well as the corresponding pseudo-ECGs are reproduced by the constructed ventricle model. In addition, the temporal and spatial vulnerability to reentry arrhythmias are quantified based on the evaluation experiments of vulnerable window and the critical length. ConclusionsThe constructed multi-scale rat ventricle model is able to reproduce both the physiological and the pathological phenomenon in different scales. Evaluation experiments suggest that the apex is the most susceptible area to arrhythmias. The model can be a promising tool for the investigation of arrhythmogenesis and the screening of anti-arrhythmic drugs.

Proarrhythmia in the p.Met207Val PITX2c-Linked Familial Atrial Fibrillation-Insights From Modeling.

Functional analysis has shown that the p.Met207Val mutation was linked to atrial fibrillation and caused an increase in transactivation activity of PITX2c, which caused changes in mRNA synthesis related to ionic channels and intercellular electrical coupling. We assumed that these changes were quantitatively translated to the functional level. This study aimed to investigate the potential impact of the PITX2c p.Met207Val mutation on atrial electrical activity through multiscale computational models. The well-known Courtemanche-Ramirez-Nattel (CRN) model of human atrial cell action potentials (APs) was modified to incorporate experimental data on the expected p.Met207Val mutation-induced changes in ionic channel currents (INaL, IKs, and IKr) and intercellular electrical coupling. The cell models for wild-type (WT), heterozygous (Mutant/Wild type, MT/WT), and homozygous (Mutant, MT) PITX2c cases were incorporated into homogeneous multicellular 1D and 2D tissue models. Effects of this mutation-induced remodeling were quantified as changes in AP profile, AP duration (APD) restitution, conduction velocity (CV) restitution and wavelength (WL). Temporal and spatial vulnerabilities of atrial tissue to the genesis of reentry were computed. Dynamic behaviors of re-entrant excitation waves (Life span, tip trajectory and dominant frequency) in a homogeneous 2D tissue model were characterized. Our results suggest that the PITX2c p.Met207Val mutation abbreviated atrial APD and flattened APD restitution curves. It reduced atrial CV and WL that facilitated the conduction of high rate atrial excitation waves. It increased the tissue's temporal vulnerability by increasing the vulnerable window for initiating reentry and increased the tissue spatial vulnerability by reducing the substrate size necessary to sustain reentry. In the 2D models, the mutation also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained reentry. In conclusion, electrical and structural remodeling arising from the PITX2c p.Met207Val mutation may increase atrial susceptibility to arrhythmia due to shortened APD, reduced CV and increased tissue vulnerability, which, in combination, facilitate initiation and maintenance of re-entrant excitation waves.

Double-wave reentry in excitable media.

A monolayer of chick embryo cardiac cells grown in an annular geometry supports two simultaneous reentrant excitation waves that circulate as a doublet. We propose a mechanism that can lead to such behavior. The velocity restitution gives the instantaneous velocity of a wave as a function of the time since the passage of the previous wave at a given point in space. Nonmonotonic restitution relationships will lead to situations in which various spacings between circulating waves are possible. In cardiology, the situation in which two waves travel in an anatomically defined circuit is referred to as double-wave reentry. Since double-wave reentry may arise as a consequence of pacing during cardiac arrhythmias, understanding the dynamic features of double-wave reentry may be helpful in understanding the physiological properties of cardiac tissue and in the design of therapy.

Synergistic Anti-arrhythmic Effects in Human Atria with Combined Use of Sodium Blockers and Acacetin.

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Developing effective and safe anti-AF drugs remains an unmet challenge. Simultaneous block of both atrial-specific ultra-rapid delayed rectifier potassium (K+) current (IKur) and the Na+ current (INa) has been hypothesized to be anti-AF, without inducing significant QT prolongation and ventricular side effects. However, the antiarrhythmic advantage of simultaneously blocking these two channels vs. individual block in the setting of AF-induced electrical remodeling remains to be documented. Furthermore, many IKur blockers such as acacetin and AVE0118, partially inhibit other K+ currents in the atria. Whether this multi-K+-block produces greater anti-AF effects compared with selective IKur-block has not been fully understood. The aim of this study was to use computer models to (i) assess the impact of multi-K+-block as exhibited by many IKur blokers, and (ii) evaluate the antiarrhythmic effect of blocking IKur and INa, either alone or in combination, on atrial and ventricular electrical excitation and recovery in the setting of AF-induced electrical-remodeling. Contemporary mathematical models of human atrial and ventricular cells were modified to incorporate dose-dependent actions of acacetin (a multichannel blocker primarily inhibiting IKur while less potently blocking Ito, IKr, and IKs). Rate- and atrial-selective inhibition of INa was also incorporated into the models. These single myocyte models were then incorporated into multicellular two-dimensional (2D) and three-dimensional (3D) anatomical models of the human atria. As expected, application of IKur blocker produced pronounced action potential duration (APD) prolongation in atrial myocytes. Furthermore, combined multiple K+-channel block that mimicked the effects of acacetin exhibited synergistic APD prolongations. Synergistically anti-AF effects following inhibition of INa and combined IKur/K+-channels were also observed. The attainable maximal AF-selectivity of INa inhibition was greatly augmented by blocking IKur or multiple K+-currents in the atrial myocytes. This enhanced anti-arrhythmic effects of combined block of Na+- and K+-channels were also seen in 2D and 3D simulations; specially, there was an enhanced efficacy in terminating re-entrant excitation waves, exerting improved antiarrhythmic effects in the human atria as compared to a single-channel block. However, in the human ventricular myocytes and tissue, cellular repolarization and computed QT intervals were modestly affected in the presence of actions of acacetin and INa blockers (either alone or in combination). In conclusion, this study demonstrates synergistic antiarrhythmic benefits of combined block of IKur and INa, as well as those of INa and combined multi K+-current block of acacetin, without significant alterations of ventricular repolarization and QT intervals. This approach may be a valuable strategy for the treatment of AF.

Modelling the effects of chloroquine on KCNJ2-linked short QT syndrome.

A gain-of-function KCNJ2 D172N mutation in KCNJ2-encoded Kir2.1 channels underlies one form of short QT syndrome (SQT3), which is associated with increased susceptibility to arrhythmias and sudden death. Anti-malarial drug chloroquine was reported as an effective inhibitor of Kir2.1 channels. Using biophysically-detailed human ventricle computer models, this study assessed the effects of chloroquine on SQT3. The ten Tusscher et al. model of human ventricular cell action potential was modified to recapitulate functional changes in the inward rectifier K+ current (IK1) due to heterozygous and homozygous forms of the D172N mutation. Mutant formulations were incorporated into multi-scale models. The blocking effects of chloroquine on ionic currents were modelled using IC50 and Hill coefficient values from literatures. Effects of chloroquine on action potential duration (APD), effective refractory period (ERP) and pseudo-ECGs were quantified. It was shown that chloroquine caused a dose-dependent reduction in IK1, prolonged APD, and decreased the maximum voltage heterogeneity. Chloroquine prolonged QT interval and declined the T-wave amplitude. Although chloroquine reduced tissue’s temporal vulnerability, it increased the minimum substrate size necessary for sustaining re-entry. The actions of chloroquine decreased arrhythmia risk, due to the reduced tissue vulnerability, prolonged ERP and wavelength of re-entrant excitation waves, which in combination prevented and terminated re-entry in the tissue models. In conclusion, the results of this study provide new evidence that the anti-arrhythmic effects of chloroquine on SQT3 and, by extension, to the possibility that chloroquine may be a potential therapeutic agent for SQT3 treatment.

Computational Analysis of the Mode of Action of Disopyramide and Quinidine on hERG-Linked Short QT Syndrome in Human Ventricles.

The short QT syndrome (SQTS) is a rare cardiac disorder associated with arrhythmias and sudden death. Gain-of-function mutations to potassium channels mediating the rapid delayed rectifier current, IKr, underlie SQTS variant 1 (SQT1), in which treatment with Na+ and K+ channel blocking class Ia anti-arrhythmic agents has demonstrated some efficacy. This study used computational modeling to gain mechanistic insights into the actions of two such drugs, disopyramide and quinidine, in the setting of SQT1. The O'Hara-Rudy (ORd) human ventricle model was modified to incorporate a Markov chain formulation of IKr describing wild type (WT) and SQT1 mutant conditions. Effects of multi-channel block by disopyramide and quinidine, including binding kinetics and altered potency of IKr/hERG channel block in SQT1 and state-dependent block of sodium channels, were simulated on action potential and multicellular tissue models. A one-dimensional (1D) transmural ventricular strand model was used to assess prolongation of the QT interval, effective refractory period (ERP), and re-entry wavelength (WL) by both drugs. Dynamics of re-entrant excitation waves were investigated using a 3D human left ventricular wedge model. In the setting of SQT1, disopyramide, and quinidine both produced a dose-dependent prolongation in (i) the QT interval, which was primarily due to IKr block, and (ii) the ERP, which was mediated by a synergistic combination of IKr and INa block. Over the same range of concentrations quinidine was more effective in restoring the QT interval, due to more potent block of IKr. Both drugs demonstrated an anti-arrhythmic increase in the WL of re-entrant circuits. In the 3D wedge, disopyramide and quinidine at clinically-relevant concentrations decreased the dominant frequency of re-entrant excitations and exhibited anti-fibrillatory effects; preventing formation of multiple, chaotic wavelets which developed in SQT1, and could terminate arrhythmias. This computational modeling study provides novel insights into the clinical efficacy of disopyramide and quinidine in the setting of SQT1; it also dissects ionic mechanisms underlying QT and ERP prolongation. Our findings show that both drugs demonstrate efficacy in reversing the SQT1 phenotype, and indicate that disopyramide warrants further investigation as an alternative to quinidine in the treatment of SQT1.

In silico assessment of the effects of quinidine, disopyramide and E-4031 on short QT syndrome variant 1 in the human ventricles.

AimsShort QT syndrome (SQTS) is an inherited disorder associated with abnormally abbreviated QT intervals and an increased incidence of atrial and ventricular arrhythmias. SQT1 variant (linked to the rapid delayed rectifier potassium channel current, IKr) of SQTS, results from an inactivation-attenuated, gain-of-function mutation (N588K) in the KCNH2-encoded potassium channels. Pro-arrhythmogenic effects of SQT1 have been well characterized, but less is known about the possible pharmacological antiarrhythmic treatment of SQT1. Therefore, this study aimed to assess the potential effects of E-4031, disopyramide and quinidine on SQT1 using a mathematical model of human ventricular electrophysiology.MethodsThe ten Tusscher et al. biophysically detailed model of the human ventricular action potential (AP) was modified to incorporate IKr Markov chain (MC) formulations based on experimental data of the kinetics of the N588K mutation of the KCNH2-encoded subunit of the IKr channels. The modified ventricular cell model was then integrated into one-dimensional (1D) strand, 2D regular and realistic tissues with transmural heterogeneities. The channel-blocking effect of the drugs on ion currents in healthy and SQT1 cells was modeled using half-maximal inhibitory concentration (IC50) and Hill coefficient (nH) values from literatures. Effects of drugs on cell AP duration (APD), effective refractory period (ERP) and pseudo-ECG traces were calculated. Effects of drugs on the ventricular temporal and spatial vulnerability to re-entrant excitation waves were measured. Re-entry was simulated in both 2D regular and realistic ventricular tissue.ResultsAt the single cell level, the drugs E-4031 and disopyramide had hardly noticeable effects on the ventricular cell APD at 90% repolarization (APD90), whereas quinidine caused a significant prolongation of APD90. Quinidine prolonged and decreased the maximal transmural AP heterogeneity (δV); this led to the decreased transmural heterogeneity of APD across the 1D strand. Quinidine caused QT prolongation and a decrease in the T-wave amplitude, and increased ERP and decreased temporal susceptibility of the tissue to the initiation of re-entry and increased the minimum substrate size necessary to prevent re-entry in the 2D regular model, and further terminated re-entrant waves in the 2D realistic model. Quinidine exhibited significantly better therapeutic effects on SQT1 than E-4031 and disopyramide.ConclusionsThe simulated pharmacological actions of quinidine exhibited antiarrhythmic effects on SQT1. This study substantiates a causal link between quinidine and QT interval prolongation in SQT1 and suggests that quinidine may be a potential pharmacological agent for treating SQT1 patients.

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, IKur. 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 IKur 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 IKur 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 IKur contributes to atrial arrhythmias. In addition, as IKur is an atrial-specific channel and a number of IKur-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 IKur.

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria.

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in IK1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on IK1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of IK1 and IKr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas IKr block alone may be a safer alternative for the E299V mutation. Combined inhibition of IKr and IKur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF.

Pro-arrhythmogenic effects of CACNA1C G1911R mutation in human ventricular tachycardia: insights from cardiac multi-scale models.

Mutations in the CACNA1C gene are associated with ventricular tachycardia (VT). Although the CACNA1C mutations were well identified in patients with cardiac arrhythmias, mechanisms by which cardiac arrhythmias are generated in such genetic mutation conditions remain unclear. In this study, we identified a novel mechanism of VT resulted from enhanced repolarization dispersion which is a key factor for arrhythmias in the CACNA1C G1911R mutation using multi-scale computational models of the human ventricle. The increased calcium influx in the mutation prolonged action potential duration (APD), produced steepened action potential duration restitution (APDR) curves as well as augmented membrane potential differences among different cell types during repolarization, increasing transmural dispersion of repolarization (DOR) and the spatial and temporal heterogeneity of cardiac electrical activities. Consequentially, the vulnerability to unidirectional conduction block in response to a premature stimulus increased at tissue level in the G1911R mutation. The increased functional repolarization dispersion anchored reentrant excitation waves in tissue and organ models, facilitating the initiation and maintenance of VT due to less meandering rotor tip. Thus, the increased repolarization dispersion caused by the G1911R mutation is a primary factor that may primarily contribute to the genesis of cardiac arrhythmias in Timothy Syndrome.

Effects of Persistent Atrial Fibrillation-Induced Electrical Remodeling on Atrial Electro-Mechanics - Insights from a 3D Model of the Human Atria.

AimsAtrial stunning, a loss of atrial mechanical contraction, can occur following a successful cardioversion. It is hypothesized that persistent atrial fibrillation-induced electrical remodeling (AFER) on atrial electrophysiology may be responsible for such impaired atrial mechanics. This simulation study aimed to investigate the effects of AFER on atrial electro-mechanics.Methods and ResultsA 3D electromechanical model of the human atria was developed to investigate the effects of AFER on atrial electro-mechanics. Simulations were carried out in 3 conditions for 4 states: (i) the control condition, representing the normal tissue (state 1) and the tissue 2–3 months after cardioversion (state 2) when the atrial tissue recovers its electrophysiological properties after completion of reverse electrophysiological remodelling; (ii) AFER-SR condition for AF-remodeled tissue with normal sinus rhythm (SR) (state 3); and (iii) AFER-AF condition for AF-remodeled tissue with re-entrant excitation waves (state 4). Our results indicate that at the cellular level, AFER (states 3 & 4) abbreviated action potentials and reduced the Ca2+ content in the sarcoplasmic reticulum, resulting in a reduced amplitude of the intracellular Ca2+ transient leading to decreased cell active force and cell shortening as compared to the control condition (states 1 & 2). Consequently at the whole organ level, atrial contraction in AFER-SR condition (state 3) was dramatically reduced. In the AFER-AF condition (state 4) atrial contraction was almost abolished.ConclusionsThis study provides novel insights into understanding atrial electro-mechanics illustrating that AFER impairs atrial contraction due to reduced intracellular Ca2+ transients.

Computational analysis of the electromechanical consequences of short QT syndrome.

Computational analysis of the electromechanical consequences of short QT syndrome.

Effects of human atrial ionic remodelling by β-blocker therapy on mechanisms of atrial fibrillation: a computer simulation.

Atrial anti-arrhythmic effects of β-adrenoceptor antagonists (β-blockers) may involve both a suppression of pro-arrhythmic effects of catecholamines, and an adaptational electrophysiological response to chronic β-blocker use; so-called 'pharmacological remodelling'. In human atrium, such remodelling decreases the transient outward (Ito) and inward rectifier (IK1) K(+) currents, and increases the cellular action potential duration (APD) and effective refractory period (ERP). However, the consequences of these changes on mechanisms of genesis and maintenance of atrial fibrillation (AF) are unknown. Using mathematical modelling, we tested the hypothesis that the long-term adaptational decrease in human atrial Ito and IK1 caused by chronic β-blocker therapy, i.e. independent of acute electrophysiological effects of β-blockers, in an otherwise un-remodelled atrium, could suppress AF. Contemporarily, biophysically detailed human atrial cell and tissue models were used to investigate effects of the β-blocker-based pharmacological remodelling. Chronic β-blockade remodelling prolonged atrial cell APD and ERP. The incidence of small amplitude APD alternans in the CRN model was reduced. At the 1D tissue level, β-blocker remodelling decreased the maximum pacing rate at which APs could be conducted. At the three-dimensional organ level, β-blocker remodelling reduced the life span of re-entry scroll waves. This study improves our understanding of the electrophysiological mechanisms of AF suppression by chronic β-blocker therapy. Atrial fibrillation suppression may involve a reduced propensity for maintenance of re-entrant excitation waves, as a consequence of increased APD and ERP.

Pro‐arrhythmogenic effects of atrial fibrillation‐induced electrical remodelling: insights from the three‐dimensional virtual human atria

Chronic atrial fibrillation (AF) is associated with structural and electrical remodelling in the atria, which are associated with a high recurrence of AF. Through biophysically detailed computer modelling, this study investigated mechanisms by which AF-induced electrical remodelling promotes and perpetuates AF. A family of Courtemanche–Ramirez–Nattel variant models of human atrial cell action potentials (APs), taking into account of intrinsic atrial electrophysiological properties, was modified to incorporate various experimental data sets on AF-induced changes of major ionic channel currents (ICaL, IKur, Ito, IK1, IKs, INaCa) and on intracellular Ca2+ handling. The single cell models for control and AF-remodelled conditions were incorporated into multicellular three-dimensional (3D) atrial tissue models. Effects of the AF-induced electrical remodelling were quantified as the changes of AP profile, AP duration (APD) and its dispersion across the atria, and the vulnerability of atrial tissue to the initiation of re-entry. The dynamic behaviour of re-entrant excitation waves in the 3D models was characterised. In our simulations, AF-induced electrical remodelling abbreviated atrial APD non-uniformly across the atria; this resulted in relatively short APDs co-existing with marked regional differences in the APD at junctions of the crista terminalis/pectinate muscle, pulmonary veins/left atrium. As a result, the measured tissue vulnerability to re-entry initiation at these tissue junctions was increased. The AF-induced electrical remodelling also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained re-entry. Under the AF-remodelled condition, re-entrant scroll waves in the 3D model degenerated into persistent and erratic wavelets, leading to fibrillation. In conclusion, realistic 3D atrial tissue models indicate that AF-induced electrical remodelling produces regionally heterogeneous and shortened APD; these respectively facilitate initiation and maintenance of re-entrant excitation waves.

Pro‐arrhythmogenic effects of the S140G KCNQ1 mutation in human atrial fibrillation – insights from modelling

Functional analysis has shown that the missense gain-in-function KCNQ1 S140G mutation associated with familial atrial fibrillation produces an increase of the slow delayed rectifier potassium current (I(Ks)). Through computer modelling, this study investigated mechanisms by which the KCNQ1 S140G mutation promotes and perpetuates atrial fibrillation. In simulations, Courtemanche et al.'s model of human atrial cell action potentials (APs) was modified to incorporate experimental data on changes of I(Ks) induced by the KCNQ1 S140G mutation. The cell models for wild type (WT) and mutant type (MT) I(Ks) were incorporated into homogeneous multicellular 2D and 3D tissue models. Effects of the mutation were quantified on AP profile, AP duration (APD) restitution, effective refractory period (ERP) restitution, and conduction velocity (CV) restitution.Temporal and spatial vulnerabilities of atrial tissue to genesis of re-entry were computed. Dynamic behaviours of re-entrant excitation waves (lifespan (LS), tip meandering patterns and dominant frequency) in 2D and 3D models were characterised. It was shown that the KCNQ1 S140G mutation abbreviated atrial APD and ERP and flattened APD and ERP restitution curves. It reduced atrial CV at low excitation rates, but increased it at high excitation rates that facilitated the conduction of high rate atrial excitation waves. Although it increased slightly tissue temporal vulnerability for initiating re-entry, it reduced markedly the minimal substrate size necessary for sustaining re-entry (increasing the tissue spatial vulnerability). In the 2D and 3D models, the mutation also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained re-entry. In the 3D model, scroll waves under the mutation condition MT conditions also degenerated into persistent and erratic wavelets, leading to fibrillation. In conclusion, increased I(Ks) due to the KCNQ1 S140G mutation increases atrial susceptibility to arrhythmia due to increased tissue vulnerability, shortened ERP and altered atrial conduction velocity, which, in combination, facilitate initiation and maintenance of re-entrant excitation waves.

Proarrhythmia in KCNJ2-linked short QT syndrome: insights from modelling

One form of the short QT syndrome (SQT3) has been linked to the D172N gain-in-function mutation to Kir2.1, which preferentially increases outward current through channels responsible for inward rectifier K(+) current (I(K1)). This study investigated mechanisms by which the Kir2.1 D172N mutation facilitates and perpetuates ventricular arrhythmias. The ten Tusscher et al. model for human ventricular action potentials (APs) was modified to incorporate changes to I(K1) based on experimentally observed changes to Kir2.1 function: both heterozygous (WT-D172N) and homozygous (D172N) mutant scenarios were studied. Cell models were incorporated into heterogeneous one-dimensional (1D), 2D tissue, and 3D models to compute the restitution curves of AP duration (APD-R), effective refractory period (ERP-R), and conduction velocity (CV). Temporal and spatial vulnerability of ventricular tissue to re-entry was measured and dynamic behaviour of re-entrant excitation waves (lifespan and dominant frequency) in 2D and 3D models of the human ventricle was characterized. D172N 'mutant' I(K1) led to abbreviated APD and ERP, as well as steeper APD-R and ERP-R curves. It reduced tissue excitability at low excitation rates but increased it at high rates. It increased tissue temporal vulnerability for initiating re-entry, but reduced the minimal substrate size necessary to sustain re-entry. SQT3 'mutant' I(K1) also stabilized and accelerated re-entrant excitation waves, leading to sustained rapid re-entry. Increased I(K1) due to the Kir2.1 D172N mutation increases arrhythmia risk due to increased tissue vulnerability, shortened ERP, and altered excitability, which in combination facilitate initiation and maintenance of re-entrant circuits.

Increased Vulnerability of Human Ventricle to Re-entrant Excitation in hERG-linked Variant 1 Short QT Syndrome

The short QT syndrome (SQTS) is a genetically heterogeneous condition characterized by abbreviated QT intervals and an increased susceptibility to arrhythmia and sudden death. This simulation study identifies arrhythmogenic mechanisms in the rapid-delayed rectifier K+ current (IKr)-linked SQT1 variant of the SQTS. Markov chain (MC) models were found to be superior to Hodgkin-Huxley (HH) models in reproducing experimental data regarding effects of the N588K mutation on KCNH2-encoded hERG. These ionic channel models were then incorporated into human ventricular action potential (AP) models and into 1D and 2D idealised and realistic transmural ventricular tissue simulations and into a 3D anatomical model. In single cell models, the N588K mutation abbreviated ventricular cell AP duration at 90% repolarization (APD90) and decreased the maximal transmural voltage heterogeneity (δV) during APs. This resulted in decreased transmural heterogeneity of APD90 and of the effective refractory period (ERP): effects that are anticipated to be anti-arrhythmic rather than pro-arrhythmic. However, with consideration of transmural heterogeneity of IKr density in the intact tissue model based on the ten Tusscher-Noble-Noble-Panfilov ventricular model, not only did the N588K mutation lead to QT-shortening and increases in T-wave amplitude, but δV was found to be augmented in some local regions of ventricle tissue, resulting in increased tissue vulnerability for uni-directional conduction block and predisposing to formation of re-entrant excitation waves. In 2D and 3D tissue models, the N588K mutation facilitated and maintained re-entrant excitation waves due to the reduced substrate size necessary for sustaining re-entry. Thus, in SQT1 the N588K-hERG mutation facilitates initiation and maintenance of ventricular re-entry, increasing the lifespan of re-entrant spiral waves and the stability of scroll waves in 3D tissue.

Modelling changes in transmural propagation and susceptibility to arrhythmia induced by volatile anaesthetics in ventricular tissue

Volatile anaesthetics such as halothane, isoflurane and sevoflurane inhibit membrane currents contributing to the ventricular action potential. Transmural variation in the extent of current blockade induces differential effects on action potential duration (APD) in the endocardium and epicardium which may be pro-arrhythmic. Biophysical modelling techniques were used to simulate the functional impact of anaesthetic-induced blockade of membrane currents on APD and effective refractory period (ERP) in rat endocardial and epicardial cell models. Additionally, the transmural conduction of excitation waves in 1-dimensional cell arrays, the tissue's vulnerability to arrhythmogenesis and dynamic behaviour of re-entrant excitation in 2-dimensional cell arrays were studied. Simulated anaesthetic exposure reduced APD and ERP in both epicardial and endocardial cell models. The reduction in APD was greater in endocardial than epicardial cells, reducing transmural APD dispersion consistent with experimental data. However, the transmural ERP dispersion was augmented. All three anaesthetics increased the width of the tissue's vulnerable window during which a premature stimulus could induce unidirectional conduction block but only halothane reduced the critical size of ventricular substrates necessary to initiate and sustain re-entrant excitation. All three anaesthetics accelerated the rate of re-entrant excitation waves, but only halothane prolonged the lifespan of re-entry. These data illustrate in silico, that modest changes in ion channel conductance abbreviate rat ventricular APD and ERP, reduce transmural APD dispersion, but augment transmural ERP dispersion. These changes collectively enhance the propensity for arrhythmia generation and provide a substrate for re-entry circuits with a longer half life than in control conditions.