- 1Department of Computer Science, University of Oxford, Oxford, United Kingdom
- 2Departamento de Ciência da Computação, Universidade Federal De São João Del-Rei, São João del Rei, Brazil
Atrial fibrillation (AF) inducibility, sustainability and response to pharmacological treatment of individual patients are expected to be determined by their ionic current properties, especially in structurally-healthy atria. Mechanisms underlying AF and optimal cardioversion are however still unclear. In this study, in-silico drug trials were conducted using a population of human structurally-healthy atria models to 1) identify key ionic current properties determining AF inducibility, maintenance and pharmacological cardioversion, and 2) compare the prognostic value for predicting individual AF cardioversion of ionic current properties and electrocardiogram (ECG) metrics. In the population of structurally-healthy atria, 477 AF episodes were induced in ionic current profiles with both steep action potential duration (APD) restitution (eliciting APD alternans), and high excitability (enabling propagation at fast rates that transformed alternans into discordant). High excitability also favored 211 sustained AF episodes, so its decrease, through prolonged refractoriness, explained pharmacological cardioversion. In-silico trials over 200 AF episodes, 100 ionic profiles and 10 antiarrhythmic compounds were consistent with previous clinical trials, and identified optimal treatments for individual electrophysiological properties of the atria. Algorithms trained on 211 simulated AF episodes exhibited >70% accuracy in predictions of cardioversion for individual treatments using either ionic current profiles or ECG metrics. In structurally-healthy atria, AF inducibility and sustainability are enabled by discordant alternans, under high excitability and steep restitution conditions. Successful pharmacological cardioversion is predicted with 70% accuracy from either ionic or ECG properties, and it is optimal for treatments maximizing refractoriness (thus reducing excitability) for the given ionic current profile of the atria.
1 Introduction
In-silico trials with human modeling and simulation represent an effective means for drug safety and efficacy assessment (Passini et al., 2017; Musuamba et al., 2021; Margara et al., 2022). Besides yielding accurate predictions of pro-arrhythmic toxicity, superior to that obtained through animal studies (Passini et al., 2017), they provide mechanistic explanations to potential adverse events. Accordingly, in-silico trials constitute a central paradigm for drug safety (Passini et al., 2021), with unexploited capabilities for drug efficacy evaluation at higher dimensional levels (i.e., whole-organ scale) (Margara et al., 2022). The latter could support the introduction of human-based multiscale modeling and simulation into precision medicine, as shown by Roney et al., 2022 for the prediction of atrial fibrillation (AF) recurrence after catheter ablation. Similarly, in-silico trials offer the possibility of identifying key modulators of successful pharmacological treatment, that could guide tailored AF therapies.
Recent studies of AF pathophysiology have highlighted the low incidence (i.e., 3% of cases) of AF in the absence of heart disease (Wyse et al., 2014). Newly identified forms of heart disease, however, coincide with the absence of atrial structural abnormalities. Channelopathies (Wyse et al., 2014), genetic variants (Ragab et al., 2020) metabolic disorders, imbalanced autonomic nervous system (Hanna et al., 2021), endocrine control (Aguilar et al., 2021) or the inflammatory response (Scott et al., 2021), might increase AF susceptibility by transiently altering the electrical substrate (Heijman, Linz and Schotten, 2021). These dynamic substrates, as opposed to static structural inhomogeneities, might underlie the high variability observed across AF paroxysms. While patients presenting the above-mentioned disorders cannot be categorized as healthy, they might have structurally-healthy atria, in which AF inducibility has proven to be common, especially in males (Kumar et al., 2012).
In-vivo studies conducted in human subjects (Kim et al., 2002; Narayan et al., 2008; Krummen et al., 2012) observed that dynamic substrates favoring AF initiation resulted from steep action potential duration (APD) restitution (APDR). Steep APDR favored APD alternans and dynamic heterogeneities in repolarization when the atria were paced at rapid rates. Simulation studies (Gong et al., 2007) showed that discordant APD alternans (i.e., APD alternating in a beat-to-beat short-long pattern, with opposing phase in neighboring regions) could derive from steep APDR and preceded AF initiation. While discordant alternans was only observed for elevated L-type Ca2+ current, the authors did not analyze further alterations. Thus, other ionic current properties, mechanisms or channel dysregulation, resulting from genetic polymorphisms (Ragab et al., 2020), imbalanced regulatory systems (Heijman, Linz and Schotten, 2021) or natural inter-subject variability in ionic densities (Muszkiewicz et al., 2018), may also engage with steep electrical restitution and facilitate AF induction in patients with structurally-healthy atria. The mechanisms are however unknown.
Early detection of AF in patients with structurally-healthy atria is crucial, given their higher eligibility for antiarrhythmic drug treatment (Andrade et al., 2020; Hindricks et al., 2021). Indeed, absence of structural heart disease and AF duration <24 h are independent predictors of successful pharmacological cardioversion to normal sinus rhythm (Boriani et al., 2004). Restoring sinus rhythm limits the atrial structural remodeling (Andrade et al., 2020) and reduces the risk of stroke and cardiovascular death compared to rate control strategies (Kirchhof et al., 2020). Nevertheless, the success rate of pharmacological treatment is currently suboptimal (Heijman, Linz and Schotten, 2021), and little is known about how inter-subject variability in channel densities (Muszkiewicz et al., 2018) determines the response of individual AF patients to antiarrhythmic drugs (Capucci et al., 2018).
This lack of knowledge might be attributed to the unfeasible, non-ethical and risky nature of studying the effect of inter-subject variability in human electrophysiology on drug efficacy. The latter would require the patient to undergo multiple treatments and invasive screening procedures. Modeling and simulation studies, on the other hand, can provide strong evidence under high spatio-temporal resolution and perfect control over the parameters of interest, avoiding all three barriers and overcoming experimental limitations (Heijman et al., 2021).
Therefore, the aim of this study is to exploit human-based multiscale modeling and simulation to 1) identify the ionic profiles that determine electrical restitution and excitability properties, and favor AF in structurally-healthy atria, 2) understand their relationship with pharmacological cardioversion, and 3) quantify the prognostic significance for successful pharmacological cardioversion of ionic current densities versus electrocardiogram (ECG) metrics. To address these goals, simulated AF episodes were induced in a population of human-based whole-atria models with varying ionic densities and steep electrical restitution properties, as observed in-vivo. Pharmacological cardioversion was subsequently attempted by simulating the effects of the antiarrhythmic drugs used clinically for restoring sinus rhythm. Simulation results were consistently evaluated against experimental and clinical data. Ultimately, predictive algorithms for AF cardioversion, in the form of clinical decision support systems, were built using the information obtained from ECG metrics and ionic current profiles.
2 Methods
2.1 Population of human atrial cardiomyocyte models
A population of human atrial cardiomyocyte models was generated as described in Muszkiewicz et al., 2018. The CRN model (Courtemanche, Ramirez and Nattel, 1998) was considered as baseline, since it reproduces the long plateau observed in paroxysmal AF patients with steep APDR (Narayan et al., 2008; Krummen et al., 2012) and it has been extensively used for developing and testing antiarrhythmic drug modeling (Loewe et al., 2014, 2015; Sutanto et al., 2019). Effects of drug modeling on the CRN model, moreover, were compared to the effects on the Grandi-Bers (Grandi et al., 2011) and Maleckar-Trayanova (Maleckar et al., 2009) models to assess model independence of key results.
Key conductances and permeabilities were sampled up to ±70% of their control ranges using Latin Hypercube sampling, including the ultrarapid, rapid and slow delayed-rectifier K+ current density (GKur, GKr and GKs), transient outward K+ current density (Gto), inward rectifier K+ current density (GK1), L-type Ca2+ current density (GCaL), fast Na+ current density (GNa), Na+/K+ pump (GNaK), Ca2+/Na+ exchanger (GNCX) and the sarcoplasmic reticulum Ca2+ release (Grel), leak (Gleak) and uptake (Gup) currents. Experimental calibration was subsequently applied using action potential characteristics obtained from human cells of patients in sinus rhythm (dataset and further explanation available in Sánchez et al., 2014; Table 1). No AF remodeling was applied to the calibrated population of atrial cardiomyocyte models.
TABLE 1. Prediction accuracy (mean and confidence interval). Most relevant features selected by the elastic-net logistic regression. Abbreviations: DF: Dominant frequency; SaE: Sample entropy. RHE: Relative harmonic energy.
The APDR curves of the experimentally-calibrated population were computed by both standard S1S2 and dynamic restitution protocols. The standard S1S2 protocol consisted of a train of stimuli applied at a fixed S1 cycle length (CL), followed by a single extrastimulus with coupling interval S2. For the dynamic restitution protocol, the train was applied until steady-state at constant S1, for decreasing CLs. Conduction velocity (CV) restitution (CVR) was measured for models yielding APDR slope greater than unity for both restitution protocols (Koller, Riccio and Gilmour, 1998; Watanabe et al., 2001), and an APDR slope greater than unity over a diastolic interval range (Qu, Weiss and Garfinkel, 1999) of at least 30 milliseconds (ms). A detailed explanation of the APDR and CVR calculation can be found in the Supplementary Material.
2.2 Populations of whole-atria models with homogeneous vs. heterogeneous ionic current properties
Ninety-seven atrial cardiomyocyte models presenting steep APDR and CVR were used to populate two human-based whole-atria models, one with homogeneous and one with heterogeneous ionic current properties (Supplementary Figure S1). The population of 97 homogeneous whole-atria models was constructed by assigning one atrial cardiomyocyte model to all regions of the atria (to isolate the contribution of ionic densities to AF). The population of 97 heterogeneous models was similarly constructed with each atrial cardiomyocyte model, but in this case, the ionic profile was modified in specific atrial regions, as previously detailed (Seemann et al., 2006; Tobón et al., 2013; Sánchez et al., 2017). Specifically, the single-cell properties of the cardiomyocyte model were assigned to the left atrial tissue, so that the restitution properties interfered with burst pacing at the pulmonary veins (see below), and scaled in the right atrium, crista terminalis, pectinate muscles, left atrial appendage and atrio-ventricular rings, as detailed in Supplementary Table S1.
Regional heterogeneities in CV and anisotropy ratio were included in both populations of whole-atria models as in (Tobón et al., 2013; Sánchez et al., 2017) (Detailed in Supplementary Table S1). Simulations of electrical propagation were conducted using the monodomain equation with the MonoAlg3D software (Sachetto Oliveira et al., 2018; available at https://github.com/rsachetto/MonoAlg3D_C).
2.3 Atrial fibrillation inducibility
The protocol for simulated AF induction mimicked the spontaneous ectopic beats reported in human patients with frequent episodes of AF (Haïssaguerre et al., 1999). As such, a burst of five periodic stimuli was applied at the left pulmonary veins, the most common foci origin site. To account for repetitive discharges with irregular CL, different CLs, ranging from 170 to 260 ms (Haïssaguerre et al., 1999; mean CL of focal discharges 175 ± 30 msec) with 10 ms increment, were tested for the burst. AF episodes longer than 3 s were considered sustained.
2.4 Atrial fibrillation pharmacological management: Drug modeling
Sustained AF episodes were subjected to 10 simulated pharmacological treatments: five drugs currently approved for restoring sinus rhythm (i.e., flecainide, propafenone, vernakalant, amiodarone, and ibutilide; Hindricks et al., 2021), dronedarone (an analogue of amiodarone), ranolazine (as a promising compound for sinus rhythm maintenance; Capucci et al., 2018) and digoxin (a rate control agent previously used for rhythm control; The Digitalis in Acute Atrial Fibrillation (DAAF) Trial Group, 1997). Drug action was simulated as simple pore-block models according to previously reported effective plasma drug concentration, and their 50% inhibitory concentration and hill coefficient profiles (Supplementary Table S2, Supplementary Figure S2). For computational feasibility, considering the large number and complexity of the simulations and the deterministic nature of the models, drug administration was modelled 3 s from AF induction (Sánchez et al., 2017), and the episode was recorded for another 4 s. AF was considered successfully cardioverted by the drug if the atria were free of arrhythmic activity before the completion of the 7 s (Matene and Jacquemet, 2012). Drug efficacy was defined as the number of AF episodes cardioverted by a given drug over the total number of AF episodes.
Moreover, to evaluate whether pharmacological AF prevention and cardioversion followed the same restitution mechanisms, the AF-induction protocol was repeated in sinus rhythm after drug application. For this, only the two drugs that steepened/flattened the APDR curve the most were considered.
2.5 Simulated electrocardiogram
Simulated 8-lead (leads I, II, V1-V6) ECGs were computed for sustained fibrillation as in Gima and Rudy, 2002. The ECG was recorded during the 3-s AF episode prior to drug application. The resulting signals were standardized and processed using pass-band filtering with cut-off frequency of [0.5–50] Hz (Alcaraz et al., 2011). Five biomarkers, namely, dominant frequency, organization index, Shannon’s spectral entropy, sample entropy and relative harmonic energy, were extracted from each lead as in (Alcaraz et al., 2011; Zeemering et al., 2018). These biomarkers provided a characterization of AF complexity in the frequency and temporal domain.
2.6 Clinical decision support systems
An algorithm for predicting the likelihood of pharmacological cardioversion (equivalent to a clinical decision support system) was built using ECG metrics, ionic current profiles and cellular electrophysiological properties as input. For this, a dataset was constructed for each antiarrhythmic drug, combining the information of the AF episodes and the outcome of drug application. Every AF episode was characterized by M = 55 features, including 15 derived from the electrophysiological properties of the atria (i.e., baseline CV, APD90, resting membrane potential (RMP) and 12 ionic current densities) and 40 from the ECG analysis (i.e., 5 biomarkers from each lead, 8-lead ECG). Each episode was labelled either as 1, if AF was successfully cardioverted, or 0, otherwise.
The whole dataset was partitioned in sub-datasets to determine the predictive power of each set of biomarkers. Clinical decision support systems were built with each sub-dataset by elastic-net logistic regression (Zeemering et al., 2018), using the Glmnet toolbox (available at http://hastie.su.domains/glmnet_matlab/). Parameter selection was done by stepwise logistic regression with alpha equal 0.5 and 100 values for lambda. The performance of the predictive models was validated through 5-fold cross-validation, repeated 20 times, randomly allocating 70–30% of data for training and testing, respectively. Further details of the algorithms and the different sub-datasets can be found in the Supplementary Material.
3 Results
3.1 Evaluation of human atrial cardiomyocyte action potential models shows that steep restitution is associated with low triangulation and high cellular excitability
A calibrated subpopulation of human atrial cardiomyocyte models, matching experimental action potential characteristics, was selected from the initial candidate population. Four hundred fifty-one cardiomyocyte models in the calibrated subpopulation showed steep APDR for both restitution protocols, and 97 additionally presented steep CVR. Figure 1 illustrates action potential morphologies (Figure 1A), representative APDR and CVR curves of steep restitution properties (Figures 1B–E) and ionic density distribution (Figure 1F) of the subpopulations of atrial cardiomyocyte models considered.
FIGURE 1. (A) Simulated action potential traces for human atrial cardiomyocyte models in the population before (grey) and after (green) experimental calibration. (B–E) Comparison of the electrical restitution between a representative cardiomyocyte model characterized by steep action potential duration restitution (APDR) alone (blue) and additional steep conduction velocity restitution (CVR, pink). (B) Dynamic APDR plotted against the cycle length, showing transient alternation, and (C) against the diastolic interval. (D) CVR. (E) Standard S1S2 APDR. (F) Ionic density distribution (color-code as per panel (A). Data analyzed using Wilcoxon test. *p < 0.05, **p < 0.01, ***p < 0.001.
Steep APDR curves occurred for longer plateau phases (triangulation (APD90-APD50) 72.6 ± 18.1 vs 126.9 ± 61.3 ms, steep vs flat APDR curves). Models additionally yielding steep CVR curves had a more negative resting membrane potential (RMP: 82.4 ± 1.7 vs. -74.5 ± 3.3 mV; steep APDR and CVR vs steep APDR alone) (Figure 1A, pink versus blue traces). The hyperpolarization of the RMP allowed propagation for shorter CLs (307 ± 15 vs 538 ± 152 ms; steep APDR and CVR vs steep APDR alone), increasing the magnitude of APD alternans (Figure 1B) and favoring a greater modulation of CV (Figure 1D).
Steep APDR alone was found in models with significantly high GCaL and low GKur compared to flat restitution (Figure 1F, blue vs green boxes). Steep APDR and CVR curves (Figure 1F, pink boxes) were associated with even greater GCaL, elevated GNa, and high density of repolarization currents, especially GNaK, GK1, GKr and Gto. The 97 atrial cardiomyocyte models exhibiting steep APDR and CVR were used to develop the populations of human whole-atria models (Supplementary Figure S1).
3.2 Atrial fibrillation is facilitated by discordant APD alternans enabled by propagation at fast rates due to high excitability, for high L-type Ca2+ and inward rectifier K+ currents
212 and 265 AF episodes were induced in 75 out of the 97 homogeneous and heterogeneous whole-atria models, respectively, following burst stimulation at 10 different CLs. Whole-atria models with 22 ionic profiles failed to induce AF for both homogeneous and heterogeneous ionic current properties. Figure 2 illustrates a representative AF episode, showing multiple rotors in the whole-atria model (Figure 2A, left), the transmembrane voltage from 1 cell (Figure 2A, right), and similarities between the simulated ECG (atrial activity, Figure 2B) and the clinical ECG (atrial activity) obtained from a human AF patient (Figure 2C). More AF episodes were induced with the heterogeneous ionic substrate, but the 75 ionic profiles of the models ensuring AF inducibility were the same in both populations.
FIGURE 2. Representative atrial fibrillation (AF) episode. (A) Simulated AF in the whole atria model (left), and time course of the transmembrane voltage (Vm) in 1 cell (right). (B) Simulated atrial electrocardiogram (normalized ECG; no units). (C) Clinical ECG obtained in AF patient (adapted from Lankveld et al., 2016 under permission).
Figure 3 illustrates the distribution of ionic densities (Figure 3A), effective refractory period (ERP) (Figure 3B), and APDR curves (Figure 3C) between the 75 atrial cardiomyocyte models favoring AF and the 22 non-inducible models. The dynamic restitution curves in Figure 3C are plotted against the CLs used for burst pacing (170–260 ms), for two consecutive beats to show magnitude of APD alternans (i.e., difference between long and short APD). Figures 3D,E compare the response of the whole-atria models to burst stimulation between AF inducible (Figure 3D, top and Figure 3E) and non-inducible models (Figure 3D, middle and bottom), to illustrate the dynamics underlying AF initiation. The analysis presented in Figure 3 considers the population with homogeneous ionic current properties.
FIGURE 3. Ionic properties of AF-inducible (black) and non-inducible models (grey): (A) Ionic density distributions. (B) Effective refractory period (ERP). (C) Dynamic action potential duration (APD) restitution curves plotted against the cycle length used for burst pacing (170–260 ms). Two consecutive beats are considered (i.e., long and short APD), showing magnitude of APD alternans. Number (N) of models able to propagate all stimuli. (D) Comparison of propagation dynamics between an AF-inducible model (top), and two non-inducible models: one showing a 1:2 conduction pattern (middle row), the other propagating all stimuli without AF induction (bottom row). Burst stimulus applied at a cycle length of 190 ms. The blue star represents the location of the burst. (E) Discordant alternans favoring conduction block and re-entry. The snapshots in (E) are an expansion of the red box in (D). Abbreviations. S-ICV: Superior and inferior cava vein; ri-rs-li-ls-PV: Right inferior, right superior, left inferior and left superior pulmonary veins. Data are expressed as medians and IQR and analyzed using Wilcoxon Test. *p < 0.05, **p < 0.01. ***p < 0.001.
Analysis of whole-atria simulations revealed AF inducibility for ionic profiles presenting high GCaL, which increased APD alternans, and GK1, which conferred high cellular excitability (Figures 3A–C). Indeed, all whole-atria models susceptible to AF exhibited propagation at CL = 260 ms, while this only occurred in 11 models (50%) free of arrhythmia (Figure 3C). Likewise, 56 (75%) AF inducible and 5 (23%) non-inducible whole-atria models enabled propagation at 170 ms (i.e., the shortest CL tested in the burst).
Figure 3D compares the response to burst stimulation between one whole-atria model susceptible to AF (top row) and two non-inducible models (one of them belonging to the 11 models that failed to propagate burst stimulation at CL = 260 ms (middle row), the other belonging to the 5 models that propagated at 170 ms without inducing AF (bottom row)). Whole-atria models unable to propagate all stimuli in the burst exhibited a 1:2 conduction pattern (i.e., propagation of one every two stimuli), which prevented engagement with the steep restitution properties, as illustrated in Figure 3D (middle): the stimulus at t = 190 ms propagates following propagation failure at t = 0 ms.
However, even when all burst stimuli propagated, some models were resistant to AF (Figure 3D, bottom). Analysis of their associated APDR curves revealed that AF-inducible models presented greater magnitude of APD alternans (67.3 ± 92.5 vs 8.9 ± 6.6 ms, AF-inducible vs non-inducible models) (Figure 3C). The latter developed into discordant alternans (i.e., APD alternating in a beat-to-beat short-long pattern, with opposing phase in neighboring regions) during the application of rapid burst stimulation, causing unidirectional block and re-entry at different sites of the left atrium (Figure 3E). An extensive explanation of discordant alternans formation after burst pacing is provided in the Supplementary Figure S3. Discordant alternans were not observed in whole-atria models free of AF.
Restoring the baseline density (i.e., ionic conductance) of the L-type Ca2+ current in AF-inducible models recapitulated the magnitude of APD alternans observed in non-inducible models (67.3 ± 92.5 vs 8.3 ± 6.0 ms, AF-inducible models with elevated vs control L-type Ca2+ current density). Similarly, alternans magnitude was reduced after decreasing the inward rectifier K+ current density to control conditions (67.3 ± 92.5 vs 18.2 ± 12.1; AF-inducible models with elevated vs control inward rectifier K+ current). The latter effect was due to a reduction in the number of models propagating at short CLs (explanation given in the Supplementary Figure S4).
3.3 Analysis of 477 AF episodes reveals atrial fibrillation sustainability supported by high excitability, with enhanced up-regulation of the inward rectifier K+ current and elevated Na+/K+ pump
AF sustainability was evaluated in the 75 whole-atria models ensuring AF inducibility considering homogeneous versus heterogeneous ionic current properties, to compare the role of the ionic profile and spatial gradients. The heterogeneous ionic substrate favored AF maintenance: 154 (58%) of the 265 AF episodes induced in heterogeneous models were sustained for over 3 s, versus 57 (27%) of 212 AF episodes induced in homogeneous models. In the former, the left atrial appendage and atrioventricular rings had a shorter APD90 than the bulk tissue, which favored rotor stabilization. Together with the APD heterogeneities within the right atrium (i.e., crista terminalis, pectinate muscles and bulk tissue) and between the right and left atrium, the heterogeneous ionic substrate presented inhomogeneous restitution properties that reinforced AF maintenance.
Overall, 211 (44%) AF episodes sustained in both populations of whole-atria models. In order to isolate the role of ionic profile rather than gradients in AF sustainability, the properties of the homogeneous models yielding sustained (>3s) and un-sustained fibrillation were analyzed. Thus, Figure 4 shows ionic density distribution (Figure 4A), ERP (Figure 4B), APDR curves (Figure 4C), and Na+ current dynamics (Figure 4D).
FIGURE 4. Comparison between atrial cardiomyocyte models favoring sustained vs un-sustained fibrillation in homogeneous whole-atria models. (A) Ionic density distribution. (B) Effective refractory period (ERP). (C) Dynamic action potential duration (APD) restitution (as in Figure 3C). Number (N) of models able to propagate all stimuli. (D) Normalized fast Na+ current (INa) peak and normalized product of INa inactivation gates. Data are expressed as medians and IQR and analyzed using Wilcoxon Test. *p < 0.05, **p < 0.01, ***p < 0.001.
Higher GNaK and GK1, and lower GNa were associated with models sustaining AF. Analysis of Na+ inactivation gates revealed that, in spite of a lower GNa, Na+ availability was higher in the sustained AF group (Figure 4D, right). This was accentuated as the pacing rate increased: the maximum Na+ current peak, which was greater in the un-sustained AF group for lower rates, became higher in models yielding sustained AF when stimulated at short CLs (Figure 4D, left). This was due to the higher GNaK and GK1, which hyperpolarized the RMP (-82.6 ± 0.7 vs. -81.5 ± 0.6 mV, sustained and un-sustained AF, respectively). Accordingly, these atrial cardiomyocyte models had higher cellular excitability (and shorter effective refractory period, Figure 4B).
At slow rates (CL = 700 ms), lower CVs were found in whole-atria models leading to sustained AF (60.3 ± 5.7 vs 62.5 ± 4.5 cm/s), as a result of lower GNa. At fast rates (CL = 260 ms), however, the opposite trend was observed (56.4 ± 9.8 vs 50.7 ± 5.7 cm/s; sustained vs unsustained AF). The magnitude of APD alternans was similar across groups.
3.4 In-silico drug trials in over 200 atrial fibrillation episodes consistently reproduce cardioversion efficacy observed in clinical trials and identify key properties determining treatment success
Pharmacological cardioversion with 10 compounds was attempted in the 211 sustained AF episodes induced in both whole-atria populations (homogeneous, 58 episodes; heterogeneous, 153 episodes). Figure 5 displays the efficacy (Figure 5A) of each treatment obtained in-silico and reported in human clinical trials. In-silico antiarrhythmic drugs are displayed according to their preferential ionic channel target. Cardioversion efficacy was higher for homogeneous versus heterogeneous whole-atria models as illustrated by dashed and solid bars, respectively in Figure 5A (described below). Figure 5B compares the percentage (%) of APD90 variation obtained after drug action with respect to control conditions in-silico and in-vitro. Corresponding results were obtained for APD90 change between simulated drug action and experimental data (Figure 5B), with further validation, illustrating rate-dependence alterations in APD, CV, and maximum upstroke velocity, provided in Supplementary Table S3.
FIGURE 5. (A) Comparison between the cardioversion efficacy (Eff.) obtained in-silico (color bars) in whole-atria models with homogeneous (Homo.) and heterogeneous (Hete.) electrophysiology (EP) and human clinical trials (black bars). Drugs are grouped according to the ionic channels they target. A threshold of 25% block at therapeutic plasma concentration has been considered to denote ion current block. The arrows point towards greater ionic density block. (B) Percentage (%) of APD90 variation after drug modeling or administration with respect to control conditions, in-silico (color bars) and in-vitro (black bars, error bars in-vitro not available, see Supplementary Table S3). The animal silhouettes represent the species analyzed in-vitro.
Moreover, a comparison of APD90 variation after drug action between the CRN, Grandi-Bers and Maleckar-Trayanova models can be found in the Supplementary Figure S5). The latter models showed a higher sensitivity to IKur inhibition than the CRN model, which presented a greater APD prolongation after IKr block. Thus, repeating the whole-atria simulations with the Grandi-Bers or Maleckar-Trayanova models is expected to produce even higher efficacy for IKur blockers than that already shown in Figure 5.
Higher cardioversion efficacy derived from class III (i.e., dronedarone and vernakalant) and class V (i.e., digoxin) antiarrhythmic drugs. Targeting IKr, IKur, INaK and IK1 yielded a greater prolongation of refractoriness than inward current block (Supplementary Figure S6). Antiarrhythmic drugs resulting in IKr inhibition (i.e, ranolazine, ibutilide, and dronedarone) showed an increase in cardioversion efficacy proportional to the extent of IKr block. However, higher success rates resulted from blocking multiple K+ channels, even for lower percentages of IKr inhibition. For instance, a 25% INaK and 40% IKr block (i.e., digoxin) was far more effective than 45% IKr (i.e., ibutilide, efficacy in-silico matching clinical trials) or 70% IKr block (i.e., dronedarone).
The high efficacy of digoxin, in agreement with some clinical trials, owed to INaK inhibition. Blocking INaK or IK1 alone raised the RMP and decreased INa availability. The subsequent prolongation of post-repolarization refractoriness was sufficient to terminate AF regardless of a subtle APD change (Supplementary Figure S6). Indeed, single IK1 block and digoxin were the pharmacological approaches that provided the highest cardioversion rates.
Combined IKr and IKur block (i.e., flecainide and vernakalant, whose cardioversion efficacies matched the efficacy reported in clinical trials) was also more effective than IKr inhibition alone, especially in the population with heterogeneous ionic current properties. In general, the heterogeneous substrate hampered pharmacological cardioversion, since it facilitated AF maintenance (Section 3.3) and antiarrhythmic drugs had an inhomogeneous effect on different atrial regions. In particular, the heterogeneous ionic substrate was characterized by a lower GKr in the right atrium compared to the left atrium (Supplementary Table S1), so that IKr inhibition lost cardioversion efficacy in this setting. Thus, the efficacy of flecainide and dronedarone, but not of vernakalant, dropped significantly in the heterogeneous ionic substrate.
Low cardioversion rates stemmed from targeting the inward currents, as a result of drug-induced APD reduction (Figures 5A,B). INa block was effective when complemented with IKur and IKr inhibition (i.e., flecainide and vernakalant). However, combined INa and ICaL block (i.e., amiodarone and propafenone) shortened the APD, yielding the formation of focal re-entrant sources (Supplementary Figures S7, S8). The cardioversion efficacy decreased proportionally to ICaL block (equivalently, resulting in larger APD abbreviation, Figures 5A,B), with amiodarone (i.e., 50% ICaL and 30% INa block) showing the lowest cardioversion rate. The efficacy of amiodarone, and especially of propafenone, differed from the efficacy reported in clinical trials (see Discussion).
Nevertheless, amiodarone-induced APD shortening significantly flattened the APDR curved (Supplementary Figure S6). Investigations on cardio-protection instead of cardioversion (i.e., inducing the AF-initiation protocol after in-silico drug application) showed that amiodarone prevented 100% of recurrent AF episodes, with only atrial flutter (i.e., regular and periodic activation) remaining in 2% of cases. In contrast, dronedarone prolonged the APD, steepening the APDR. Therefore, while effective in preventing AF induced at very short CLs (i.e., 170–200 ms), 49% of AF episodes were still inducible in dronedarone-pretreated substrates.
3.5 Optimal pharmacological treatments for individual ionic profiles are those that maximize prolonged refractoriness thus minimizing excitability
Forty-nine AF episodes, sustained by 40 whole-atria models, were only successfully cardioverted by a single drug. Eleven episodes sustained in 9 models were only cardioverted by dronedarone, 14 episodes in 13 models by flecainide, and 24 episodes in 18 models by vernakalant. Figure 6 illustrates the ionic density distribution of the whole-atrial models that only responded favorably to dronedarone, flecainide, or vernakalant, and how the ERP of these models, grouped according to successful treatment, was modified after applying each drug.
FIGURE 6. (A) Comparison of the ionic density distribution between atrial cardiomyocyte models showing favorable response only to dronedarone, flecainide or vernakalant. (B) Effect of the three antiarrhythmic drugs on the effective refractory period (ERP), stratifying the atrial cardiomyocyte models according to the pharmacological treatment that produced successful cardioversion. Colors on the bars as in (A). Data are expressed as medians and IQR and analyzed using Wilcoxon Test. *p < 0.1, **p < 0.05, ***p < 0.01. ****p < 0.001.
The extent of ERP prolongation caused by each treatment in a particular model explained the ability for AF cardioversion. Thus, models showing a favorable response to dronedarone (N = 9) also presented a greater ERP prolongation after dronedarone application. Similarly, flecainide prolonged the ERP more than dronedarone or vernakalant in whole-atria models only cardioverted by flecainide (N = 13), and the same applied for vernakalant (Figure 6B).
Analysis of the ionic current profiles revealed that atrial cardiomyocyte models showing greater ERP prolongation and thus, successful cardioversion after vernakalant application, had high GKur and low GKr. This agrees with vernakalant action, with the strongest IKur inhibitory profile within the tested antiarrhythmic drugs. By contrast, dronedarone cardioverted AF in models presenting the opposite situation: higher GKr and lower GKur. Flecainide resulted more effective for an intermediate scenario, since it blocked GKur and GKr, but to a lesser extent than vernakalant and dronedarone, respectively (Figure 6A). Whole-atrial models presenting high GK1 and GNaK, and therefore high excitability (Figure 4), responded better to vernakalant and flecainide than to dronedarone. While none of them caused INaK or IK1 inhibition, the former yielded greater inward current block (Supplementary Table S2, Supplementary Figure S2). Thus, vernakalant- and flecainide-induced INa and ICaL inhibition facilitated pharmacological cardioversion in those models presenting higher excitability.
3.6 Accuracy of prediction of AF cardioversion with specific treatments is 70% using both electrocardiogram metrics and ionic current profiles
Clinical decision support systems were built based on the information of the atrial model sustaining AF (i.e., ionic current properties) and the ECG metrics obtained during AF. Figure 7 and Table 1 show the prediction accuracy and receiver operating characteristics (ROC) curves of successful cardioversion for the drugs illustrated in Figure 6.
FIGURE 7. Receiver operating characteristics (ROC) curves for the prediction of successful cardioversion with dronedarone, flecainide and vernakalant using the ECG, ionic properties and both. The confidence interval around each ROC is obtained repeating 5-fold cross-validation 20 times. Abbreviations: AUC: Areas under the curve.
Both the information of the ionic densities and the ECG biomarkers yielded high-accuracy predictions of successful AF cardioversion, with the former being slightly superior to the latter. Higher accuracies were obtained for dronedarone, due to a higher specificity (Figure 7). Dronedarone cardioverted around 50% of episodes, so that the dataset was more balanced (i.e., similar proportion of episodes cardioverted and not cardioverted). By contrast, the lowest specificity was obtained for vernakalant (70% cardioversion efficacy), with random oversampling slightly improving the prediction.
Three ionic densities were consistently selected as relevant features: GK1, GKur, and GNaK. As illustrated in Figure 6, these ionic densities differed to a greater extent in atrial cardiomyocyte models that responded better to flecainide, vernakalant or dronedarone, and thus, were revealing of the most favorable treatment for each episode. Conversely, the dominant frequency, sample entropy and the relative harmonic energy (directly related to the dominant frequency) were selected as relevant ECG metrics, suggesting that AF dynamics also had a crucial role on cardioversion efficacy.
4 Discussion
Key features underlying AF inducibility, maintenance and pharmacological cardioversion are identified in a population of structurally-healthy atria with varying ionic current properties. In-silico trials for 10 antiarrhythmic treatments are consistent with previously-reported clinical trials, and identify mechanisms and predictors of successful cardioversion. Specific findings are: 1) AF inducibility in structurally-healthy atria might be explained by discordant APD alternans resulting from steep electrical restitution and high excitability (occurring for elevated L-type Ca2+ and inward rectifier K+ currents); 2) AF maintenance is also enabled by high excitability, for profiles with high inward rectifier K+ and Na+/K+ pump currents; 3) pharmacological cardioversion is optimal for treatments minimizing excitability, through maximal prolongation of refractoriness, which is dependent on individual ionic profiles; 4) successful pharmacological cardioversion is predicted in virtual human atria with >70% accuracy using both ionic current profiles or ECG metrics. Methodologically, our study demonstrates the power of in-silico drug trials conducted in human whole-atria models constructed based on hundreds of ionic current profiles, with perfect control over the parameters studied and high spatio-temporal resolution.
4.1 Precision medicine for atrial fibrillation management: In-silico drug trials
Currently, rhythm control therapy for AF follows a “one-size fits all” approach (Kany et al., 2021). However, the response of individual patients to rhythm control therapy is highly variable, and predicting its success is still a challenge (Kany et al., 2021). Previous studies have demonstrated that analyzing AF complexity through the ECG might help predicting flecainide cardioversion efficacy (Zeemering et al., 2018). Others assert that pharmacogenetics would bring the best tailored approach (Capucci et al., 2018). In our study, we have shown that pharmacological cardioversion can be predicted with over 70% accuracy using either ionic or ECG metrics in virtual human atria. This is important as it points towards the possibility of predicting the response of pharmacological treatment using both non-invasive markers and ionic current properties. Moreover, it highlights the integral part of multi-scale modeling and simulation for AF treatment personalization (Roney et al., 2022).
In our in-silico drug trial, GK1 was selected as one of the most informative properties for predicting AF cardioversion, together with GNaK and GKur. Inhibiting the Na+/K+ pump reduced the cellular excitability and hampered AF maintenance, consistent with its potential for rhythm control (Bueno-Orovio et al., 2014) and rate control strategy (i.e., digoxin; Hindricks et al., 2021). Similarly, the importance of the inward rectifier K+ current in regulating AF inducibility and maintenance has been proven by previous experimental and computational studies (Noujaim et al., 2007; Sánchez et al., 2017). Results obtained with the cellular model used, however, should be carefully interpreted, given the strong model-dependent effects. This is especially important when evaluating the antiarrhythmic efficacy of IK1 inhibition, since a higher sensitivity for IK1 block is reported for Courtemanche compared to other ionic models (Sánchez et al., 2014; Sutanto, 2022). Moreover, IK1 inhibition might have opposing effects on triggers and substrates. While effective in terminating fibrillatory activity, as shown in this study, blocking IK1 might promote triggered activity and focal excitations (Filgueiras-Rama et al., 2012). In this sense, safe restoration of sinus rhythm (i.e., absence of ventricular arrhythmias) resulted from IK1 inhibition in goats with persistent AF (Ji et al., 2017). Moreover, as in our simulations, higher cardioversion efficacies resulted from applying chloroquine (IK1 blocker, among others) than flecainide in sheep models (Filgueiras-Rama et al., 2012).
Flecainide cardioverted 56% of the simulated AF episodes, matching the efficacy observed in clinical trials with AF patients (Boriani et al., 1998; Reisinger et al., 2004). Compared to flecainide, a higher proportion of AF episodes (i.e., 65%) were terminated for vernakalant simulation. Similarly, a non-randomized retrospective study (Pohjantähti-Maaroos et al., 2019) observed higher number of patients cardioverted to sinus rhythm with vernakalant than flecainide (67 vs 46%, respectively).
Dronedarone also proved to be effective in terminating a great number of simulated AF episodes. Dronedarone has shown a positive impact on clinical endpoints for paroxysmal AF patients and compared to amiodarone, it delays and reduces AF occurrence (Singh et al., 2008). However, amiodarone has proven more effective than dronedarone in rhythm control (Hindricks et al., 2021).
In our study, amiodarone showed the lowest cardioversion efficacy. Nevertheless, cardioversion rates as high as 87% are observed after 24 h from its administration (Bash et al., 2012). Amiodarone inhibits multiple ionic channels, including ICaL and INa. Blocking both channels reduces cellular excitability (Loewe et al., 2014; Liberos et al., 2016) and amiodarone-induced INa inhibition has been suggested as the main mechanism underlying rotor termination (Wilhelms et al., 2013). In our study, however, the baseline atrial cardiomyocyte models yielding AF presented increased GNa and GCaL compared to the baseline model. Consequently, the inward current block induced by amiodarone was not sufficient to reduce excitability and terminate AF. This is in agreement with previous simulations studies, that showed that the antiarrhythmic effects of INa and ICaL block depended on their basal ionic densities (Liberos et al., 2016). Of note, amiodarone is primarily recommended in the context of severe heart failure (Andrade et al., 2020), in which ICaL is expected to be down-regulated (Denham et al., 2018). Thus, amiodarone-induced ICaL inhibition might have a lower impact on the APD and higher on excitability than in our study. It has been also suggested that amiodarone metabolites, such as desethylamiodarone, which accumulates after long-term therapy with amiodarone, might be the actual responsible for the increased atrial refractory period (Talajic, DeRoode and Nattel, 1987).
Furthermore, amiodarone-induced ICaL inhibition flattened the APDR curve, reduced the frequency and amplitude of APD alternans and prevented AF inducibility in the structurally-healthy atria. Therefore, while ineffective for AF termination, amiodarone hindered the formation of dynamic substrates derived from steep electrical restitution.
4.2 Ionic current dysregulation and dynamic electrophysiological substrates in structurally-healthy atria
Steep electrical restitution was the mechanism adopted in this study to trigger dynamic substrates that facilitate AF, as observed in human AF patients (Narayan et al., 2008; Krummen et al., 2012). Unlike previous simulation works that changed single ionic currents (Garfinkel et al., 2000), steep electrical restitution in our study resulted from specific multivariate ionic density distributions. These distributions were obtained by scaling uniformly all ionic densities of the cellular model, which mimicked inter-patient variability in ionic densities and preserved the high non-linearity between ion channel conductances and the action potential morphology (Elliott et al., 2021). Moreover, we considered a detailed atrial geometry with realistic anisotropy ratios and CV, as opposed to using monolayers (Virag et al., 2002) or isotropic conduction (Gong et al., 2007), and replicate accurate AF-initiation conditions.
As such, we observed that increased GCaL, GK1 and GNaK characterized those atrial cardiomyocyte models favoring AF inducibility and maintenance in structurally-healthy atria. Previous simulation studies observed that increasing GCaL steepens the APD restitution slope, potentially leading to wave break-up and fibrillation (Garfinkel et al., 2000). In our study, however, we found that steep restitution was only engaged for high cellular excitability. The latter was enabled by GK1 up-regulation, and could be aggravated by increased acetylcholine-activated inward rectifier K+ (IK,ACh) current, not considered in the present study.
Likewise, elevated GK1 and GNaK ensured AF maintenance. In agreement with our simulations, increased Na+/K+ pump was linked to higher susceptibility to postoperative AF (Tran et al., 2009). Increased inward rectifier K+ current has been extensively associated with the AF-induced electrophysiological remodeling (Dobrev and Ravens, 2003) and proven essential for rotor stabilization in the electrically remodeled atria (Pandit et al., 2005). In our population of structurally-healthy atria, however, high GK1 was additionally required for AF inducibility. Similarly, other computational studies have observed that in the absence of fibrosis, the inhibition of the inward rectifier K+ current considerably hindered AF initiation (Gharaviri et al., 2021). A genetic analysis performed in a population of kindreds (Xia et al., 2005) identified a gain-of-function mutation in KCNJ2, encoding the Kir2. x channels of the inward rectifier K+ current, underlying familial AF. Missense mutations in KCNJ2 were also reported in human patients with paroxysmal AF (Deo et al., 2013). Therefore, as shown in this study, increased inward rectifier K+ current constitutes an AF predisposing condition that additionally favors AF perpetuation.
4.3 Limitations and future perspectives
Drug action was simulated as simple pore-block models as the data required were available for the compounds investigated. However, it is possible that this might have increased the model sensitivity to repolarization prolonging effects, explaining the grater efficacy of class III antiarrhythmic drugs compared to class Ic. Neglecting the rate dependence of class Ic agents and thus, potentially underestimating INa inhibition, might explain the different results obtained for propafenone and flecainide. Similarly, increasing the resting membrane potential through IK1 inhibition, therefore hampering INa (in)activation, resulted more effective than direct INa block. This might reflect the importance of class Ic agents for binding with different affinities to the different states of the Na+ channel, which was not considered in this study.
Moreover, pharmacological treatment and the performance of the predictive model should be further evaluated considering, not only electrophysiological, but also anatomical and structural variability. Even the electrophysiological variability considered was small. The 97 models selected, which presented the ionic current properties required to induce AF in structurally-healthy atria, had similar morphologies (i.e., a rapid phase 3 repolarization and relatively hyperpolarized resting membrane potential). Therefore, a wider range of action potentials has to be taken into account in future studies.
Restricting the study to similar action potential morphologies could also explain why the accuracy of the clinical decision support system was not higher than 70%. A better stratification of the ionic current properties that respond to specific antiarrhythmic drugs would result from considering a wider ionic current distribution in the first place (i.e., wider range of action potential morphologies). The latter could improve the prediction accuracy and boost associated prediction metrics. In this sense, employing other machine learning algorithms might also add extra predictive power.
In the clinic, obtaining the ionic density of the patient presents a challenge, especially without undertaking an invasive procedure (i.e., biopsy during surgery). However, the finding that similar accuracies are obtained with the ECG and ionic metrics for predicting successful AF cardioversion, suggests that the ECG already contain information of the ionic current properties. Thus, further studies should investigate whether the atrial ionic profile of a patient could be characterized non-invasively through the ECG.
5 Conclusion
In structurally-healthy atria, AF inducibility and sustainability are enabled by high excitability and steep restitution, due to elevated L-type Ca2+ and inward rectifier K+ currents. Accordingly, decreasing excitability, through prolonged refractoriness, results in pharmacological cardioversion. However, maximal prolongation of refractoriness depends on the interaction between the ionic distribution of the atria and the ionic currents targeted by the antiarrhythmic drug, highlighting optimal pharmacological treatments for individual ionic current profiles. Therefore, successful cardioversion is predicted in virtual human atria with 70% accuracy from both ionic current and ECG properties.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
All authors equally contributed to the conception of the work, revising it critically for important intellectual content, and final approval of the version to be published, ensuring that questions related to the accuracy or integrity of any part of the work were appropriately investigated and resolved. BR and AB-O acted as supervisors of the research, AR, RS and JC mainly contributed in methodological aspects, and AD was responsible for conducting all the simulations, analysis of results, and drafting the work.
Funding
This work was supported by a British Heart Foundation (BHF) Intermediate Fellowship (FS/17/22/32,644 to AB-O), a Wellcome Trust Senior Fellowship in Basic Biomedical Sciences (214290/Z/18/Z to BR), and the CompBioMed Centre of Excellence in Computational Biomedicine (European Union’s Horizon 2020; grant agreements 675,451, 823,712). This project has also received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No.860974 (to AD). The authors acknowledge additional support from an Infrastructure for Impact Award from the National Centre for the Replacement, Refinement and Reduction of Animals in Research (NC/P001076/1), the Oxford BHF Centre of Research Excellence (RE/13/1/30,181), PRACE, Piz Daint at the Swiss National Supercomputing Centre, Switzerland (ICEI-PRACE grants icp005 and icp013) and Fapemig.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.966046/full#supplementary-material
References
Aguilar M., Rose A., Takawale S., Nattel S., Reilly S. (2021). New aspects of endocrine control of atrial fibrillation and possibilities for clinical translation. Cardiovasc. Res. 117 (7), 1645–1661. doi:10.1093/cvr/cvab080
Alcaraz R., Sandberg L., Sörnmo J. J., Rieta J. J. (2011). Classification of paroxysmal and persistent atrial fibrillation in ambulatory ECG recordings. IEEE Trans. Biomed. Eng. 58 (5), 1441–1449. doi:10.1109/TBME.2011.2112658
Andrade J. G., Aguilar M., Atzema C., Bell A., Cairns J. A., Cheung C. C., et al. (2020). The 2020 Canadian cardiovascular society/Canadian heart rhythm society comprehensive guidelines for the management of atrial fibrillation. Can. J. Cardiol. 36 (12), 1847–1948. doi:10.1016/j.cjca.2020.09.001
Bash L. D., Buono G. M., Davies A., Martin K., Fahrbach H., et al. (2012). Systematic review and meta-analysis of the efficacy of cardioversion by vernakalant and comparators in patients with atrial fibrillation. Cardiovasc. Drugs Ther. 26 (2), 167–179. doi:10.1007/s10557-012-6374-4
Boriani G., Bifft A., Capucci G., Botto T., BroffoniM. T., et al. (1998). Conversion of recent-onset Atrial Fibrillation to sinus rhythm: Effects of different drug protocols. Pacing Clin. Electro 21 (11 II), 2470–2474. doi:10.1111/j.1540-8159.1998.tb01203.x
Boriani G., Diemberger M., Biffi C., Martignani A., Branzi A. (2004). Pharmacological cardioversion of atrial fibrillation. Drugs 64 (24), 2741–2762. doi:10.2165/00003495-200464240-00003
Bueno-Orovio A., Sánchez E., Pueyo B., Rodriguez B. (2014). Na/K pump regulation of cardiac repolarization: Insights from a systems biology approach. Pflugers Arch. - Eur. J. Physiol. 466 (2), 183–193. doi:10.1007/s00424-013-1293-1
Capucci A., Cipolletta L., Guerra F., Giannini I. (2018). Emerging pharmacotherapies for the treatment of atrial fibrillation. Expert Opin. Emerg. Drugs 23 (1), 25–36. doi:10.1080/14728214.2018.1446941
Courtemanche M., Ramirez R. J., Nattel S. (1998). Ionic mechanisms underlying human atrial action potential properties: Insights from a mathematical model. Am. J. Physiology-Heart Circulatory Physiology 275 (1 44-1), H301–H321. doi:10.1152/ajpheart.1998.275.1.h301
Denham N. C., Pearman J. L., Caldwell G. W. P., Madders D. A., Eisner A. W., et al. (2018). Calcium in the pathophysiology of atrial fibrillation and heart failure. Front. Physiol. 9, 1380. doi:10.3389/fphys.2018.01380
Deo M., Ruan S. V., Pandit K., Shah O., Berenfeld A., et al. (2013). KCNJ2 mutation in short QT syndrome 3 results in atrial fibrillation and ventricular proarrhythmia. Proc. Natl. Acad. Sci. U.S.A. 110 (11), 4291–4296. doi:10.1073/pnas.1218154110
Dobrev D., Ravens U. (2003). Remodeling of cardiomyocyte ion channels in human atrial fibrillation. Basic Res. Cardiol. 98 (3), 137–148. doi:10.1007/s00395-003-0409-8
Elliott J., Belen L., Mainardi J. F., Rodriguez Matas J. F. (2021). A comparison of regional classification strategies implemented for the population based approach to modelling atrial fibrillation. Mathematics 9 (14), 1686. doi:10.3390/math9141686
Filgueiras-Rama D., Martins S., Mironov M., Yamazaki C. J., Calvo S. R., et al. (2012). Chloroquine terminates stretch-induced atrial fibrillation more effectively than flecainide in the sheep heart. Circ Arrhythmia Electrophysiol. 5 (3), 561–570. doi:10.1161/CIRCEP.111.966820
Garfinkel A., Kim Z., Voroshilovsky J. R., Qu M. H., et al. (2000). Preventing ventricular fibrillation by flattening cardiac restitution. Proc. Natl. Acad. Sci. U.S.A. 97 (11), 6061–6066. doi:10.1073/pnas.090492697
Gharaviri A., Pezzuto M., Potse G., Conte S., Zeemering V., et al. (2021). Synergistic antiarrhythmic effect of inward rectifier current inhibition and pulmonary vein isolation in a 3D computer model for atrial fibrillation. Europace. 23, I161–I168. doi:10.1093/europace/euaa413
Gima K., Rudy Y. (2002). Ionic current basis of electrocardiographic waveforms. Circulation Res. 90 (8), 889–896. doi:10.1161/01.RES.0000016960.61087.86
Gong Y., Xie K. M., Stein A., Garfinkel C. A., Culianu B. B., et al. (2007). Mechanism underlying initiation of paroxysmal atrial flutter/atrial fibrillation by ectopic foci. Circulation 115 (16), 2094–2102. doi:10.1161/CIRCULATIONAHA.106.656504
Grandi E., Pandit N., Voigt A. J., Workman D., Dobrev J., et al. (2011). Human atrial action potential and Ca 2+ model. Circ. Res. 109 (9), 1055–1066. doi:10.1161/CIRCRESAHA.111.253955
Haïssaguerre M., Jaïs P., Shah D. C., Takahashi A., Hocini M., Quiniou G., et al. (1999). Spontaneous initiation of atrial fibrillation by ectopic impulses originating in the pulmonary veins. Cardiol. Rev. 7 (2), 65. doi:10.1097/00045415-199903000-00006
Hanna P., Buch S., Stavrakis C., Meyer J. D., Tompkins J. L., et al. (2021). Neuroscientific therapies for atrial fibrillation. Cardiovasc. Res. 117 (7), 1732–1745. doi:10.1093/cvr/cvab172
Heijman J., Linz D., Schotten U. (2021). Dynamics of atrial fibrillation mechanisms and comorbidities. Annu. Rev. Physiol. 83, 83–106. doi:10.1146/annurev-physiol-031720-085307
Hindricks G., Potpara N., Dagres E., Arbelo J. J., Bax C., et al. (2021). 2020 ESC Guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European Association for Cardio-Thoracic Surgery (EACTS): The Task Force for the diagnosis and management of atrial fibrillation of the European Society of Cardiology (ESC) Developed with the special contribution of the European Heart Rhythm Association (EHRA) of the ESC. Eur. Heart J. 42 (5), 373–498. doi:10.1093/eurheartj/ehaa612
Ji Y., Varkevisser D., Opacic A., Bossu M., Kuiper J. D. M., et al. (2017). The inward rectifier current inhibitor PA-6 terminates atrial fibrillation and does not cause ventricular arrhythmias in goat and dog models. Br. J. Pharmacol. 174 (15), 2576–2590. doi:10.1111/bph.13869
Kany S., Reissmann A., Metzner P., Kirchhof D., Darbar R. B., Schnabel R. B. (2021). Genetics of atrial fibrillation-practical applications for clinical management: If not now, when and how? Cardiovasc. Res. 117 (7), 1718–1731. doi:10.1093/cvr/cvab153
Kim B. S., Kim Y.-H., Hwang G.-S., Pak H.-N., Lee S. C., Shim W. J., et al. (2002). Action potential duration restitution kinetics in human atrial fibrillation. J. Am. Coll. Cardiol. 39 (8), 1329–1336. doi:10.1016/S0735-1097(02)01760-6
Kirchhof P., Camm A., Goette A., Brandes L., Eckardt A., et al. (2020). Early rhythm-control therapy in patients with atrial fibrillation. N. Engl. J. Med. 383 (14), 1305–1316. doi:10.1056/nejmoa2019422
Koller M. L., Riccio M. L., Jr. R. F. (1998). Dynamic restitution of action potential duration during electrical alternans and ventricular fibrillation. Am. J. Physiology-Heart Circulatory Physiology 275 (5 44-5), H1635–H1642. doi:10.1152/ajpheart.1998.275.5.h1635
Krummen D. E., Bayer J., Ho G., Ho M. R., Smetak P., et al. (2012). Mechanisms of human atrial fibrillation initiation. Circ Arrhythmia Electrophysiol. 5 (6), 1149–1159. doi:10.1161/CIRCEP.111.969022
Kumar S., Kalman F., Sutherland S. J., Spence S., Finch P. B., Sparks P. B. (2012). Atrial fibrillation inducibility in the absence of structural heart disease or clinical atrial fibrillation. Circ Arrhythmia Electrophysiol. 5 (3), 531–536. doi:10.1161/CIRCEP.111.968859
Lankveld T., de Vos C. B., Limantoro I., Zeemering S., Dudink E., Crijns H. J., et al. (2016). Systematic analysis of ECG predictors of sinus rhythm maintenance after electrical cardioversion for persistent atrial fibrillation. Heart rhythm. 13 (5), 1020–1027. doi:10.1016/j.hrthm.2016.01.004
Liberos A., Bueno-Orovio A., Rodrigo M., Ravens U., Hernandez-Romero I., Fernandez-Aviles F., et al. (2016). Balance between sodium and calcium currents underlying chronic atrial fibrillation termination: An in silico intersubject variability study. Heart rhythm. 13 (12), 2358–2365. doi:10.1016/j.hrthm.2016.08.028
Loewe A., Lutz M., Wilhelms D., Sinnecker P., Barthel E. P., et al. (2014). In-silico assessment of the dynamic effects of amiodarone and dronedarone on human atrial patho-electrophysiology. Europace 16, iv30–iv38. doi:10.1093/europace/euu230
Loewe A., Xu E. P., Scholz O., Dössel G., Seemann G. (2015). Understanding the cellular mode of action of vernakalant using a computational model: Answers and new questions. Curr. Dir. Biomed. Eng. 1 (1), 418–422. doi:10.1515/cdbme-2015-0101
Maleckar M. M., Greenstein W. R., Giles N. A., Trayanova N. A. (2009). K+current changes account for the rate dependence of the action potential in the human atrial myocyte. Am. J. Physiology-Heart Circulatory Physiology 297 (4), H1398–H1410. doi:10.1152/ajpheart.00411.2009
Margara F., Psaras Y., Jenny Wang Z. (2022). Human iPSC-CMs and in-silico technologies define mechanisms and accelerate targeted pharmacogenetics in hypertrophic cardiomyopathy. bioRxiv 2022, 495324. doi:10.1101/2022.06.08.495324
Matene E., Jacquemet V. (2012). Fully automated initiation of simulated episodes of atrial arrhythmias. Europace 14 (5), v17–v24. doi:10.1093/europace/eus271
Musuamba F. T., Skottheim Rusten R., Lesage G., Russo R., Bursi L., et al. (2021). Scientific and regulatory evaluation of mechanistic in silico drug and disease models in drug development: Building model credibility. CPT Pharmacometrics Syst. Pharmacol. 10 (8), 804–825. doi:10.1002/psp4.12669
Muszkiewicz A., Liu A., Bueno-Orovio B. A. J., Lawson K., Burrage B., et al. (2018). From ionic to cellular variability in human atrial myocytes: An integrative computational and experimental study. Am. J. Physiology-Heart Circulatory Physiology 314 (5), H895–H916. doi:10.1152/ajpheart.00477.2017
Narayan S. M., Kazi D. E., Krummen W. J., Rappel W.-J. (2008). Repolarization and activation restitution near human pulmonary veins and atrial fibrillation initiation. J. Am. Coll. Cardiol. 52 (15), 1222–1230. doi:10.1016/j.jacc.2008.07.012
Noujaim S. F., Pandit O., Berenfeld K., Vikstrom M., Cerrone S., et al. (2007). Up-regulation of the inward rectifier K+current (Ik1) in the mouse heart accelerates and stabilizes rotors. J. Physiol. 578 (1), 315–326. doi:10.1113/jphysiol.2006.121475
Pandit S. V., Berenfeld J. M. B., Anumonwo R. M., Zaritski J., Kneller S., et al. (2005). Ionic determinants of functional reentry in a 2-D model of human atrial cells during simulated chronic atrial fibrillation. Biophysical J. 88 (6), 3806–3821. doi:10.1529/biophysj.105.060459
Passini E., Britton H. R., Lu J., Rohrbacher A. N., Hermans D. J., et al. (2017). Human in silico drug trials demonstrate higher accuracy than animal models in predicting clinical pro-arrhythmic cardiotoxicity. Front. Physiol. 8 (SEP), 668–715. doi:10.3389/fphys.2017.00668
Passini E., Zhou X., Trovato C., Britton O. J., Bueno-Orovio A., Rodriguez B. (2021). The virtual assay software for human in silico drug trials to augment drug cardiac testing. J. Comput. Sci. 52, 101202. doi:10.1016/j.jocs.2020.101202
Pohjantähti-Maaroos H., Hyppölä M., Lekkala E., Sinisalo A., Heikkola J., Hartikainen J. (2019). Intravenous vernakalant in comparison with intravenous flecainide in the cardioversion of recent-onset atrial fibrillation. Eur. Heart J. Acute Cardiovasc. Care 8 (2), 114–120. doi:10.1177/2048872617728558
Qu Z., Weiss J. N., Garfinkel A. (1999). Cardiac electrical restitution properties and stability of reentrant spiral waves: A simulation study. Am. J. Physiology-Heart Circulatory Physiology 276 (1 45-1), H269–H283. doi:10.1152/ajpheart.1999.276.1.h269
Ragab A. A. Y., Sitorus B. B. J. J. M., Brundel N. M. S., Groot N. M. S. d. (2020). The genetic puzzle of familial atrial fibrillation. Front. Cardiovasc. Med. 7, 14–18. doi:10.3389/fcvm.2020.00014
Reisinger J., Gatterer E., Lang W., Vanicek T., Eisserer G., Bachleitner T., et al. (2004). Flecainide versus ibutilide for immediate cardioversion of atrial fibrillation of recent onset. Eur. Heart J. 25 (15), 1318–1324. doi:10.1016/j.ehj.2004.04.030
Roney C. H., Sim J., Yu M., Beach A., Mehta J., et al. (2022). Predicting atrial fibrillation recurrence by combining population data and virtual cohorts of patient-specific left atrial models. Circ Arrhythmia Electrophysiol. 15 (2), e010253. doi:10.1161/CIRCEP.121.010253
Sachetto Oliveira R., Martins Rocha D., Burgarelli W., Meira C., Constantinides R., Weber dos Santos R. (2018). Performance evaluation of GPU parallelization, space-time adaptive algorithms, and their combination for simulating cardiac electrophysiology. Int. J. Numer. Meth Biomed. Engng 34 (2), e2913–17. doi:10.1002/cnm.2913
Sánchez C., Bueno-Orovio E., Pueyo B., Rodríguez B. (2017). Atrial fibrillation dynamics and ionic block effects in six Heterogeneous human 3D virtual atria with distinct repolarization dynamics. Front. Bioeng. Biotechnol. 5 (MAY), 29–13. doi:10.3389/fbioe.2017.00029
Sánchez C., Bueno-Orovio E., Wettwer S., Loose J., Simon U., et al. (2014). Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation. PLoS ONE 9 (8), e105897. doi:10.1371/journal.pone.0105897
Scott Jr L., Fender A., Saljic L., Li X., Chen X., et al. (2021). NLRP3 inflammasome is a key driver of obesity-induced atrial arrhythmias. Cardiovasc. Res. 117 (7), 1746–1759. doi:10.1093/cvr/cvab024
Seemann G., Höper F. B., Sachse O., Dössel A. V., Holden H., Zhang H. (2006). Heterogeneous three-dimensional anatomical and electrophysiological model of human atria. Phil. Trans. R. Soc. A 364 (1843), 1465–1481. doi:10.1098/rsta.2006.1781
Singh B. N., Connolly S. J., Crijns H. J., Roy D., Kowey P. R., Capucci A., et al. (2008). Efficacy of dronedarone for maintenance of sinus rhythm in atrial fibrillation or flutter. Results of the EURIDIS and ADONIS. Kardiologiia 48 (7), 58–59.
Sutanto H. (2022). Individual contributions of cardiac ion channels on atrial repolarization and reentrant waves: A multiscale in-silico study. Jcdd 9 (1), 28. doi:10.3390/jcdd9010028
Sutanto H., Laudy L., Clerx M., Dobrev D., Crijns H. J. G. M., Heijman J. (2019). Maastricht antiarrhythmic drug evaluator (MANTA): A computational tool for better understanding of antiarrhythmic drugs. Pharmacol. Res. 148, 104444. doi:10.1016/j.phrs.2019.104444
Talajic M., DeRoode M. R., Nattel S. (1987). Comparative electrophysiologic effects of intravenous amiodarone and desethylamiodarone in dogs: Evidence for clinically relevant activity of the metabolite. Circulation 75 (1), 265–271. doi:10.1161/01.CIR.75.1.265
The Digitalis in Acute Atrial Fibrillation (DAAF) Trial Group (1997). Intravenous digoxin in acute atrial fibrillation: Results of a randomized, placebo-controlled multicentre trial in 239 patients. European Heart Journal 18 (4), 649–654. doi:10.1093/oxfordjournals.eurheartj.a015311
Tobón C., Ruiz-Villa E., Heidenreich L., Romero F., Hornero J., Saiz J. (2013). A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship. PLoS ONE 8 (2), e50883. doi:10.1371/journal.pone.0050883
Tran C. T., Schmidt J. B., Christensen K., Kjeldsen K. (2009). Atrial Na,K-ATPase increase and potassium dysregulation accentuate the risk of postoperative atrial fibrillation. Cardiology 114 (1), 1–7. doi:10.1159/000209264
Virag N., Jacquemet C. S., Henriquez S., Zozor O., Blanc J. M., et al. (2002). Study of atrial arrhythmias in a computer model based on magnetic resonance images of human atria. Chaos 12 (3), 754–763. doi:10.1063/1.1483935
Watanabe M. A., Fenton S. J., Evans H. M., Hastings A., Karma A. (2001). Mechanisms for discordant alternans. J. Cardiovasc. Electrophysiol. 12 (2), 196–206. doi:10.1046/j.1540-8167.2001.00196.x
Wilhelms M., Hettmann M. M., Maleckar J. T., Koivumäki O., Dössel G., Seemann G. (2013). Benchmarking electrophysiological models of human atrial myocytes. Front. Physio. 3, 487–516. JAN(January). doi:10.3389/fphys.2012.00487
Wyse D. G., Van Gelder P. T., Ellinor A. S., Go J. M., Kalman S. M., et al. (2014). Lone atrial fibrillation. J. Am. Coll. Cardiol. 63 (17), 1715–1723. doi:10.1016/j.jacc.2014.01.023
Xia M., Jin S., Bendahhou Y., He M. M., Larroque Y., et al. (2005). A Kir2.1 gain-of-function mutation underlies familial atrial fibrillation. Biochem. Biophysical Res. Commun. 332 (4), 1012–1019. doi:10.1016/j.bbrc.2005.05.054
Keywords: atrial fibrillation, in-silico, drug trials, cardioversion, electrocardiogram
Citation: Dasí A, Roy A, Sachetto R, Camps J, Bueno-Orovio A and Rodriguez B (2022) In-silico drug trials for precision medicine in atrial fibrillation: From ionic mechanisms to electrocardiogram-based predictions in structurally-healthy human atria. Front. Physiol. 13:966046. doi: 10.3389/fphys.2022.966046
Received: 10 June 2022; Accepted: 23 August 2022;
Published: 15 September 2022.
Edited by:
Omer Berenfeld, University of Michigan, United StatesReviewed by:
Eric A. Sobie, Icahn School of Medicine at Mount Sinai, United StatesWayne Rodney Giles, University of Calgary, Canada
Copyright © 2022 Dasí, Roy, Sachetto, Camps, Bueno-Orovio and Rodriguez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Blanca Rodriguez, YmxhbmNhLnJvZHJpZ3VlekBjcy5veC5hYy51aw==; Albert Dasí, YWxiZXJ0LmRhc2lpbWFydGluZXpAY3Mub3guYWMudWs=