Characterizing the arrhythmogenic substrate in personalized models of atrial fibrillation: sensitivity to mesh resolution and pacing protocol in AF models

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

In Silico TRials guide optimal stratification of ATrIal FIbrillation patients to Catheter Ablation and pharmacological medicaTION: the i-STRATIFICATION study

AIMS

Patients with persistent atrial fibrillation (AF) experience 50% recurrence despite pulmonary vein isolation (PVI), and no consensus is established for secondary treatments. The aim of our i-STRATIFICATION study is to provide evidence for stratifying patients with AF recurrence after PVI to optimal pharmacological and ablation therapies, through in silico trials.

METHODS AND RESULTS

A cohort of 800 virtual patients, with variability in atrial anatomy, electrophysiology, and tissue structure (low-voltage areas, LVAs), was developed and validated against clinical data from ionic currents to electrocardiogram. Virtual patients presenting AF post-PVI underwent 12 secondary treatments. Sustained AF developed in 522 virtual patients after PVI. Second ablation procedures involving left atrial ablation alone showed 55% efficacy, only succeeding in the small right atria (<60 mL). When additional cavo-tricuspid isthmus ablation was considered, Marshall-PLAN sufficed (66% efficacy) for the small left atria (<90 mL). For the bigger left atria, a more aggressive ablation approach was required, such as anterior mitral line (75% efficacy) or posterior wall isolation plus mitral isthmus ablation (77% efficacy). Virtual patients with LVAs greatly benefited from LVA ablation in the left and right atria (100% efficacy). Conversely, in the absence of LVAs, synergistic ablation and pharmacotherapy could terminate AF. In the absence of ablation, the patient's ionic current substrate modulated the response to antiarrhythmic drugs, being the inward currents critical for optimal stratification to amiodarone or vernakalant.

CONCLUSION

In silico trials identify optimal strategies for AF treatment based on virtual patient characteristics, evidencing the power of human modelling and simulation as a clinical assisting tool.

Anti- and pro-fibrillatory effects of pulmonary vein isolation gaps in human atrial fibrillation digital twins

Although pulmonary vein isolation (PVI) gaps and extrapulmonary vein triggers contribute to recurrence after atrial fibrillation (AF) ablation, their precise mechanisms remain unproven. Our study assessed the impact of PVI gaps on rhythm outcomes using a human AF digital twin. We included 50 patients (76.0% with persistent AF) who underwent catheter ablation with a realistic AF digital twin by integrating computed tomography and electroanatomical mapping. We evaluated the final rhythm status, including AF and atrial tachycardia (AT), across 600 AF episodes, considering factors including PVI level, PVI gap number, and pacing locations. Our findings revealed that antral PVI had a significantly lower ratio of AF at the final rhythm (28% vs. 56%, p = 0.002) than ostial PVI. Increasing PVI gap numbers correlated with an increased ratio of AF at the final rhythm (p < 0.001). Extra-PV induction yielded a higher ratio of AF at the final rhythm than internal PV induction (77.5% vs. 59.0%, p < 0.001). In conclusion, our human AF digital twin model helped assess AF maintenance mechanisms. Clinical trial registration: https://www.clinicaltrials.gov ; Unique identifier: NCT02138695.

Computational Modeling and Simulation of the Fibrotic Human Atria

Patient-specific modeling of atrial electrical activity enables the execution of simulations that can provide mechanistic insights and provide novel solutions to vexing clinical problems. The geometry and fibrotic remodeling of the heart can be reconstructed from clinical-grade medical scans and used to inform personalized models with detail incorporated at the cell- and tissue-scale to represent changes in image-identified diseased regions. Here, we provide a rubric for the reconstruction of realistic atrial models from pre-segmented 3D renderings of the left atrium with fibrotic tissue regions delineated, which are the output from clinical-grade systems for quantifying fibrosis. We then provide a roadmap for using those models to carry out patient-specific characterization of the fibrotic substrate in terms of its potential to harbor reentrant drivers via cardiac electrophysiology simulations.

Fat infiltration in the infarcted heart as a paradigm for ventricular arrhythmias

Infiltrating adipose tissue (inFAT) has been recently found to co-localize with scar in infarcted hearts and may contribute to ventricular arrhythmias (VAs), a life-threatening heart rhythm disorder. However, the contribution of inFAT to VA has not been well-established. We investigated the role of inFAT versus scar in VA through a combined prospective clinical and mechanistic computational study. Using personalized computational heart models and comparing the results from simulations of VA dynamics with measured electrophysiological abnormalities during the clinical procedure, we demonstrate that inFAT, rather than scar, is a primary driver of arrhythmogenic propensity and is frequently present in critical regions of the VA circuit. We determined that, within the VA circuitry, inFAT, as opposed to scar, is primarily responsible for conduction slowing in critical sites, mechanistically promoting VA. Our findings implicate inFAT as a dominant player in infarct-related VA, challenging existing paradigms and opening the door for unexplored anti-arrhythmic strategies.

Fast Characterization of Inducible Regions of Atrial Fibrillation Models With Multi-Fidelity Gaussian Process Classification

Computational models of atrial fibrillation have successfully been used to predict optimal ablation sites. A critical step to assess the effect of an ablation pattern is to pace the model from different, potentially random, locations to determine whether arrhythmias can be induced in the atria. In this work, we propose to use multi-fidelity Gaussian process classification on Riemannian manifolds to efficiently determine the regions in the atria where arrhythmias are inducible. We build a probabilistic classifier that operates directly on the atrial surface. We take advantage of lower resolution models to explore the atrial surface and combine seamlessly with high-resolution models to identify regions of inducibility. We test our methodology in 9 different cases, with different levels of fibrosis and ablation treatments, totalling 1,800 high resolution and 900 low resolution simulations of atrial fibrillation. When trained with 40 samples, our multi-fidelity classifier that combines low and high resolution models, shows a balanced accuracy that is, on average, 5.7% higher than a nearest neighbor classifier. We hope that this new technique will allow faster and more precise clinical applications of computational models for atrial fibrillation. All data and code accompanying this manuscript will be made publicly available at: https://github.com/fsahli/AtrialMFclass.

Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation

Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

Computational modeling identifies embolic stroke of undetermined source patients with potential arrhythmic substrate

Cardiac magnetic resonance imaging (MRI) has revealed fibrosis in embolic stroke of undetermined source (ESUS) patients comparable to levels seen in atrial fibrillation (AFib). We used computational modeling to understand the absence of arrhythmia in ESUS despite the presence of putatively pro-arrhythmic fibrosis. MRI-based atrial models were reconstructed for 45 ESUS and 45 AFib patients. The fibrotic substrate's arrhythmogenic capacity in each patient was assessed computationally. Reentrant drivers were induced in 24/45 (53%) ESUS and 22/45 (49%) AFib models. Inducible models had more fibrosis (16.7 ± 5.45%) than non-inducible models (11.07 ± 3.61%; p<0.0001); however, inducible subsets of ESUS and AFib models had similar fibrosis levels (p=0.90), meaning that the intrinsic pro-arrhythmic substrate properties of fibrosis in ESUS and AFib are indistinguishable. This suggests that some ESUS patients have latent pre-clinical fibrotic substrate that could be a future source of arrhythmogenicity. Thus, our work prompts the hypothesis that ESUS patients with fibrotic atria are spared from AFib due to an absence of arrhythmia triggers.