Image-Based Computational Evaluation of the Effects of Atrial Wall Thickness and Fibrosis on Re-entrant Drivers for Atrial Fibrillation

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.

Systematic in-silico evaluation of fibrosis effects on re-entrant wave dynamics in atrial tissue

Despite the key role of fibrosis in atrial fibrillation (AF), the effects of different spatial distributions and textures of fibrosis on wave propagation mechanisms in AF are not fully understood. To clarify these aspects, we performed a systematic computational study to assess fibrosis effects on the characteristics and stability of re-entrant waves in electrically-remodelled atrial tissues. A stochastic algorithm, which generated fibrotic distributions with controlled overall amount, average size, and orientation of fibrosis elements, was implemented on a monolayer spheric atrial model. 245 simulations were run at changing fibrosis parameters. The emerging propagation patterns were quantified in terms of rate, regularity, and coupling by frequency-domain analysis of correspondent synthetic bipolar electrograms. At the increase of fibrosis amount, the rate of reentrant waves significantly decreased and higher levels of regularity and coupling were observed (p < 0.0001). Higher spatial variability and pattern stochasticity over repetitions was observed for larger amount of fibrosis, especially in the presence of patchy and compact fibrosis. Overall, propagation slowing and organization led to higher stability of re-entrant waves. These results strengthen the evidence that the amount and spatial distribution of fibrosis concur in dictating re-entry dynamics in remodeled tissue and represent key factors in AF maintenance.

Chamber-specific wall thickness features in human atrial fibrillation

Persistent atrial fibrillation (AF) is not effectively treated due to a lack of adequate tools for identifying patient-specific AF substrates. Recent studies revealed that in 30-50% of patients, persistent AF is maintained by localized drivers not only in the left atrium (LA) but also in the right atrium (RA). The chamber-specific atrial wall thickness (AWT) features underlying AF remain elusive, though the important role of AWT in AF is widely acknowledged. We aimed to provide direct evidence of the existence of distinguished RA and LA AWT features underlying AF drivers by analysing functionally and structurally mapped human hearts ex vivo. Coronary-perfused intact human atria (n = 7, 47 ± 14 y.o.; two female) were mapped using panoramic near-infrared optical mapping during pacing-induced AF. Then the hearts were imaged at approximately 170 µm3 resolution by 9.4 T gadolinium-enhanced MRI. The heart was segmented, and 3D AWT throughout atrial chambers was estimated and analysed. Optical mapping identified six localized RA re-entrant drivers in four hearts and four LA drivers in three hearts. All RA AF drivers were anchored to the pectinate muscle junctions with the crista terminalis or atrial walls. The four LA AF drivers were in the posterior LA. RA (n = 4) with AF drivers were thicker with greater AWT variation than RA (n = 3) without drivers (5.4 ± 2.6 mm versus 5.0 ± 2.4 mm, T-test p < 0.05; F-test p < 0.05). Furthermore, AWT in RA driver regions was thicker and varied more than in RA non-driver regions (5.1 ± 2.5 mm versus 4.4 ± 2.2 mm, T-test p < 0.05; F-test p < 0.05). On the other hand, LA (n = 3) with drivers was thinner than the LA (n = 4) without drivers. In particular, LA driver regions were thinner than the rest of LA regions (3.4 ± 1.0 mm versus 4.2 ± 1.0 mm, T-test p < 0.05). This study demonstrates chamber-specific AWT features of AF drivers. In RA, driver regions are thicker and have more variable AWT than non-driver regions. By contrast, LA drivers are thinner than non-drivers. Robust evaluation of patient-specific AWT features should be considered for chamber-specific targeted ablation.

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

Exploring interpretability in deep learning prediction of successful ablation therapy for atrial fibrillation

Background: Radiofrequency catheter ablation (RFCA) therapy is the first-line treatment for atrial fibrillation (AF), the most common type of cardiac arrhythmia globally. However, the procedure currently has low success rates in dealing with persistent AF, with a reoccurrence rate of ∼50% post-ablation. Therefore, deep learning (DL) has increasingly been applied to improve RFCA treatment for AF. However, for a clinician to trust the prediction of a DL model, its decision process needs to be interpretable and have biomedical relevance.

Aim: This study explores interpretability in DL prediction of successful RFCA therapy for AF and evaluates if pro-arrhythmogenic regions in the left atrium (LA) were used in its decision process.

Methods: AF and its termination by RFCA have been simulated in MRI-derived 2D LA tissue models with segmented fibrotic regions (n = 187). Three ablation strategies were applied for each LA model: pulmonary vein isolation (PVI), fibrosis-based ablation (FIBRO) and a rotor-based ablation (ROTOR). The DL model was trained to predict the success of each RFCA strategy for each LA model. Three feature attribution (FA) map methods were then used to investigate interpretability of the DL model: GradCAM, Occlusions and LIME.

Results: The developed DL model had an AUC (area under the receiver operating characteristic curve) of 0.78 ± 0.04 for predicting the success of the PVI strategy, 0.92 ± 0.02 for FIBRO and 0.77 ± 0.02 for ROTOR. GradCAM had the highest percentage of informative regions in the FA maps (62% for FIBRO and 71% for ROTOR) that coincided with the successful RFCA lesions known from the 2D LA simulations, but unseen by the DL model. Moreover, GradCAM had the smallest coincidence of informative regions of the FA maps with non-arrhythmogenic regions (25% for FIBRO and 27% for ROTOR).

Conclusion: The most informative regions of the FA maps coincided with pro-arrhythmogenic regions, suggesting that the DL model leveraged structural features of MRI images to identify such regions and make its prediction. In the future, this technique could provide a clinician with a trustworthy decision support tool.

Identifying locations susceptible to micro-anatomical reentry using a spatial network representation of atrial fibre maps

Micro-anatomical reentry has been identified as a potential driver of atrial fibrillation (AF). In this paper, we introduce a novel computational method which aims to identify which atrial regions are most susceptible to micro-reentry. The approach, which considers the structural basis for micro-reentry only, is based on the premise that the accumulation of electrically insulating interstitial fibrosis can be modelled by simulating percolation-like phenomena on spatial networks. Our results suggest that at high coupling, where micro-reentry is rare, the micro-reentrant substrate is highly clustered in areas where the atrial walls are thin and have convex wall morphology, likely facilitating localised treatment via ablation. However, as transverse connections between fibres are removed, mimicking the accumulation of interstitial fibrosis, the substrate becomes less spatially clustered, and the bias to forming in thin, convex regions of the atria is reduced, possibly restricting the efficacy of localised ablation. Comparing our algorithm on image-based models with and without atrial fibre structure, we find that strong longitudinal fibre coupling can suppress the micro-reentrant substrate, whereas regions with disordered fibre orientations have an enhanced risk of micro-reentry. With further development, these methods may be useful for modelling the temporal development of the fibrotic substrate on an individualised basis.

A Mathematical Model for Electrical Activity in Pig Atrial Tissue

State of the art mathematical models are currently used to bridge the gap between basic research conducted in the laboratory and preclinical research conducted on large animals, which ultimately paves the way for clinical translation. In this regard, there is a great need for models that can be used alongside experiments for in-depth investigation and validation. One such experimental model is the porcine atrium, which is commonly used to study the mechanisms of onset and control of atrial fibrillation in the context of its surgical management. However, a mathematical model of pig atria is lacking. In this paper, we present the first ionically detailed mathematical model of porcine atrial electrophysiology, at body temperature. The model includes 12 ionic currents, 4 of which were designed based on experimental patch-clamp data directly obtained from literature. The formulations for the other currents are adopted from the human atrial model, and modified for porcine specificity based on our measured restitution data for different action potential characteristics: resting membrane potential, action potential amplitude, maximum upstroke velocity and action potential duration and different levels of membrane voltage repolarization. The intracellular Ca²⁺ dynamics follows the Luo-Rudy formulation for guinea pig ventricular cardiomyocytes. The resulting model represents “normal” cells which are formulated as a system of ordinary differential equations. We extend our model to two dimensions to obtain plane wave propagation in tissue with a velocity of 0.58 m/s and a wavelength of 8 cm. The wavelength reduces to 5 cm when the tissue is paced at 200 ms. Using S1-S2 cross-field protocol, we demonstrate in an 11.26 cm square simulation domain, the ability to initiate single spiral waves (rotation period ≃ 180 ms) that remain stable for more than 40 s. The spiral tip exhibits hypermeander. In agreement with previous experimental results using pig atria, our model shows that early repolarization is primarily driven by a calcium-mediated chloride current, I_ClCa, which is completely inactivated at high pacing frequencies. This is a condition that occurs only in porcine atria. Furthermore, the model shows spatiotemporal chaos with reduced repolarization.

Restitution Slope Affects the Outcome of Dominant Frequency Ablation in Persistent Atrial Fibrillation: CUVIA-AF2 Post-Hoc Analysis Based on Computational Modeling Study

Introduction

Although the dominant frequency (DF) localizes the reentrant drivers and the maximal slope of the action potential duration (APD) restitution curve (Smax) reflects the tendency of the wave-break, their interaction has never been studied. We hypothesized that DF ablation has different effects on atrial fibrillation (AF) depending on Smax.

Methods

We studied the DF and Smax in 25 realistic human persistent AF model samples (68% male, 60 ± 10 years old). Virtual AF was induced by ramp pacing measuring Smax, followed by spatiotemporal DF evaluation for 34 s. We assessed the DF ablation effect depending on Smax in both computational modeling and a previous clinical trial, CUVIA-AF (170 patients with persistent AF, 70.6% male, 60 ± 11 years old).

Results

Mean DF had an inverse relationship with Smax regardless of AF acquisition timing (p &lt; 0.001). Virtual DF ablations increased the defragmentation rate compared to pulmonary vein isolation (PVI) alone (p = 0.015), especially at Smax &lt;1 (61.5 vs. 7.7%, p = 0.011). In post-DF ablation defragmentation episodes, DF was significantly higher (p = 0.002), and Smax was lower (p = 0.003) than in episodes without defragmentation. In the post-hoc analysis of CUVIA-AF2, we replicated the inverse relationship between Smax and DF (r = −0.47, p &lt; 0.001), and we observed better rhythm outcomes of clinical DF ablations in addition to a PVI than of empirical PVI at Smax &lt;1 [hazard ratio 0.45, 95% CI (0.22–0.89), p = 0.022; log-rank p = 0.021] but not at ≥ 1 (log-rank p = 0.177).

Conclusion

We found an inverse relationship between DF and Smax and the outcome of DF ablation after PVI was superior at the condition with Smax &lt;1 in both in-silico and clinical trials.

EP-PINNs: Cardiac Electrophysiology Characterisation Using Physics-Informed Neural Networks

Accurately inferring underlying electrophysiological (EP) tissue properties from action potential recordings is expected to be clinically useful in the diagnosis and treatment of arrhythmias such as atrial fibrillation. It is, however, notoriously difficult to perform. We present EP-PINNs (Physics Informed Neural Networks), a novel tool for accurate action potential simulation and EP parameter estimation from sparse amounts of EP data. We demonstrate, using 1D and 2D in silico data, how EP-PINNs are able to reconstruct the spatio-temporal evolution of action potentials, whilst predicting parameters related to action potential duration (APD), excitability and diffusion coefficients. EP-PINNs are additionally able to identify heterogeneities in EP properties, making them potentially useful for the detection of fibrosis and other localised pathology linked to arrhythmias. Finally, we show EP-PINNs effectiveness on biological in vitro preparations, by characterising the effect of anti-arrhythmic drugs on APD using optical mapping data. EP-PINNs are a promising clinical tool for the characterisation and potential treatment guidance of arrhythmias.

Predicting Atrial Fibrillation Recurrence by Combining Population Data and Virtual Cohorts of Patient-Specific Left Atrial Models

BACKGROUND

Current ablation therapy for atrial fibrillation is suboptimal, and long-term response is challenging to predict. Clinical trials identify bedside properties that provide only modest prediction of long-term response in populations, while patient-specific models in small cohorts primarily explain acute response to ablation. We aimed to predict long-term atrial fibrillation recurrence after ablation in large cohorts, by using machine learning to complement biophysical simulations by encoding more interindividual variability.

METHODS

Patient-specific models were constructed for 100 atrial fibrillation patients (43 paroxysmal, 41 persistent, and 16 long-standing persistent), undergoing first ablation. Patients were followed for 1 year using ambulatory ECG monitoring. Each patient-specific biophysical model combined differing fibrosis patterns, fiber orientation maps, electrical properties, and ablation patterns to capture uncertainty in atrial properties and to test the ability of the tissue to sustain fibrillation. These simulation stress tests of different model variants were postprocessed to calculate atrial fibrillation simulation metrics. Machine learning classifiers were trained to predict atrial fibrillation recurrence using features from the patient history, imaging, and atrial fibrillation simulation metrics.

RESULTS

We performed 1100 atrial fibrillation ablation simulations across 100 patient-specific models. Models based on simulation stress tests alone showed a maximum accuracy of 0.63 for predicting long-term fibrillation recurrence. Classifiers trained to history, imaging, and simulation stress tests (average 10-fold cross-validation area under the curve, 0.85±0.09; recall, 0.80±0.13; precision, 0.74±0.13) outperformed those trained to history and imaging (area under the curve, 0.66±0.17) or history alone (area under the curve, 0.61±0.14).

CONCLUSION

A novel computational pipeline accurately predicted long-term atrial fibrillation recurrence in individual patients by combining outcome data with patient-specific acute simulation response. This technique could help to personalize selection for atrial fibrillation ablation.

The impact of the atrial wall thickness in normal/mild late‐gadolinium enhancement areas on atrial fibrillation rotors in persistent atrial fibrillation patients

Background

Some of atrial fibrillation (AF) drivers are found in normal/mild late‐gadolinium enhancement (LGE) areas, as well as moderate ones. The atrial wall thickness (AWT) has been reported to be important as a possible AF substrate. However, the AWT and degree of LGEs as an AF substrate has not been fully validated in humans.

Objective

The purpose of this study was to evaluate the impact of the AWT in normal/mild LGE areas on AF drivers.

Methods

A total of 287 segments in 15 persistent AF patients were assessed. AF drivers were defined as non‐passively activated areas (NPAs), where rotational activation was frequently observed, and were detected by the novel real‐time phase mapping (ExTRa Mapping), mild LGE areas were defined as areas with a volume ratio of the enhancement voxel of 0% to <10%. The AWT was defined as the minimum distance from the manually determined endocardium to the epicardial border on the LGE‐MRI.

Results

NPAs were found in 20 (18.0%) of 131 normal/mild LGE areas where AWT was significantly thicker than that in the passively activated areas (PAs) (2.5 ± 0.3 vs. 2.2 ± 0.3 mm, p < .001). However, NPAs were found in 41 (26.3%) of 156 moderate LGE areas where AWT was thinner than that of PAs (2.1 ± 0.2 mm vs. 2.23 ± 0.3 mm, p = .02). An ROC curve analysis yielded an optimal cutoff value of 2.2 mm for predicting the presence of an NPA in normal/mild LGE areas.

Conclusion

The location of AF drivers in normal/mild LGE areas might be more accurately identified by evaluating AWT.

Atrial slow conduction develops and dynamically expands during premature stimulation in an animal model of persistent atrial fibrillation

Slow conduction areas and conduction block in the atria are considered pro-arrhythmic conditions. Studies examining the size and distribution of slow conduction regions in the context of persistent atrial fibrillation (AF) may help to develop improved therapeutic strategies for patients with AF. In this work, we studied the differences of size and number in slow conduction areas between control and persistent AF goats and the influence of propagation direction on the development of these pathological conduction areas. Epicardial atrial electrical activations from the left atrial roof were optically mapped with physiological pacing cycle lengths and for the shortest captured cycle lengths. The recordings were converted to local activation times and conduction velocity measures. Regions with slow conduction velocity (less than [Formula: see text]) were identified. The size of the connected regions and the number of non-connected regions were counted for propagation from different orthogonal directions. We found that regions of slow conduction significantly increases in our 15 persistent AF goat recordings in response to premature stimulation (24.4±4.3% increase to 36.6±4.4%, p < 0.001). This increase is driven by an increase of size from (3.70±0.89[mm2] to 6.36±0.91[mm2], p = 0.014) for already existing regions and not by generation of new slow conduction regions (11.6±1.8 vs. 13±1.9, p = 0.242). In 12 control goat recordings, no increase from baseline pacing to premature pacing was found. Similarly, size of the slow conduction areas and the count did not change significantly in control animals.

Time-Averaged Wavefront Analysis Demonstrates Preferential Pathways of Atrial Fibrillation, Predicting Pulmonary Vein Isolation Acute Response

Electrical activation during atrial fibrillation (AF) appears chaotic and disorganised, which impedes characterisation of the underlying substrate and treatment planning. While globally chaotic, there may be local preferential activation pathways that represent potential ablation targets. This study aimed to identify preferential activation pathways during AF and predict the acute ablation response when these are targeted by pulmonary vein isolation (PVI). In patients with persistent AF (n = 14), simultaneous biatrial contact mapping with basket catheters was performed pre-ablation and following each ablation strategy (PVI, roof, and mitral lines). Unipolar wavefront activation directions were averaged over 10 s to identify preferential activation pathways. Clinical cases were classified as responders or non-responders to PVI during the procedure. Clinical data were augmented with a virtual cohort of 100 models. In AF pre-ablation, pathways originated from the pulmonary vein (PV) antra in PVI responders (7/7) but not in PVI non-responders (6/6). We proposed a novel index that measured activation waves from the PV antra into the atrial body. This index was significantly higher in PVI responders than non-responders (clinical: 16.3 vs. 3.7%, p = 0.04; simulated: 21.1 vs. 14.1%, p = 0.02). Overall, this novel technique and proof of concept study demonstrated that preferential activation pathways exist during AF. Targeting patient-specific activation pathways that flowed from the PV antra to the left atrial body using PVI resulted in AF termination during the procedure. These PV activation flow pathways may correspond to the presence of drivers in the PV regions.

Reinforcement Learning to Improve Image-Guidance of Ablation Therapy for Atrial Fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia and currently affects more than 650,000 people in the United Kingdom alone. Catheter ablation (CA) is the only AF treatment with a long-term curative effect as it involves destroying arrhythmogenic tissue in the atria. However, its success rate is suboptimal, approximately 50% after a 2-year follow-up, and this high AF recurrence rate warrants significant improvements. Image-guidance of CA procedures have shown clinical promise, enabling the identification of key patient anatomical and pathological (such as fibrosis) features of atrial tissue, which require ablation. However, the latter approach still suffers from a lack of functional information and the need to interpret structures in the images by a clinician. Deep learning plays an increasingly important role in biomedicine, facilitating efficient diagnosis and treatment of clinical problems. This study applies deep reinforcement learning in combination with patient imaging (to provide structural information of the atria) and image-based modelling (to provide functional information) to design patient-specific CA strategies to guide clinicians and improve treatment success rates. To achieve this, patient-specific 2D left atrial (LA) models were derived from late-gadolinium enhancement (LGE) MRI scans of AF patients and were used to simulate patient-specific AF scenarios. Then a reinforcement Q-learning algorithm was created, where an ablating agent moved around the 2D LA, applying CA lesions to terminate AF and learning through feedback imposed by a reward policy. The agent achieved 84% success rate in terminating AF during training and 72% success rate in testing. Finally, AF recurrence rate was measured by attempting to re-initiate AF in the 2D atrial models after CA with 11% recurrence showing a great improvement on the existing therapies. Thus, reinforcement Q-learning algorithms can predict successful CA strategies from patient MRI data and help to improve the patient-specific guidance of CA therapy.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to ICaL block. A combined in vitro and in silico approach

The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (I_CaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after I_CaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when I_CaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, I_CaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of I_CaL block.

Understanding PITX2-Dependent Atrial Fibrillation Mechanisms through Computational Models

Atrial fibrillation (AF) is a common arrhythmia. Better prevention and treatment of AF are needed to reduce AF-associated morbidity and mortality. Several major mechanisms cause AF in patients, including genetic predispositions to AF development. Genome-wide association studies have identified a number of genetic variants in association with AF populations, with the strongest hits clustering on chromosome 4q25, close to the gene for the homeobox transcription PITX2. Because of the inherent complexity of the human heart, experimental and basic research is insufficient for understanding the functional impacts of PITX2 variants on AF. Linking PITX2 properties to ion channels, cells, tissues, atriums and the whole heart, computational models provide a supplementary tool for achieving a quantitative understanding of the functional role of PITX2 in remodelling atrial structure and function to predispose to AF. It is hoped that computational approaches incorporating all we know about PITX2-related structural and electrical remodelling would provide better understanding into its proarrhythmic effects leading to development of improved anti-AF therapies. In the present review, we discuss advances in atrial modelling and focus on the mechanistic links between PITX2 and AF. Challenges in applying models for improving patient health are described, as well as a summary of future perspectives.

Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals With a Hybrid in silico and in vivo Dataset

In patients with atrial fibrillation, intracardiac electrogram signal amplitude is known to decrease with increased structural tissue remodeling, referred to as fibrosis. In addition to the isolation of the pulmonary veins, fibrotic sites are considered a suitable target for catheter ablation. However, it remains an open challenge to find fibrotic areas and to differentiate their density and transmurality. This study aims to identify the volume fraction and transmurality of fibrosis in the atrial substrate. Simulated cardiac electrograms, combined with a generalized model of clinical noise, reproduce clinically measured signals. Our hybrid dataset approach combines in silico and clinical electrograms to train a decision tree classifier to characterize the fibrotic atrial substrate. This approach captures different in vivo dynamics of the electrical propagation reflected on healthy electrogram morphology and synergistically combines it with synthetic fibrotic electrograms from in silico experiments. The machine learning algorithm was tested on five patients and compared against clinical voltage maps as a proof of concept, distinguishing non-fibrotic from fibrotic tissue and characterizing the patient's fibrotic tissue in terms of density and transmurality. The proposed approach can be used to overcome a single voltage cut-off value to identify fibrotic tissue and guide ablation targeting fibrotic areas.

Computational models of atrial fibrillation: achievements, challenges, and perspectives for improving clinical care

Despite significant advances in its detection, understanding and management, atrial fibrillation (AF) remains a highly prevalent cardiac arrhythmia with a major impact on morbidity and mortality of millions of patients. AF results from complex, dynamic interactions between risk factors and comorbidities that induce diverse atrial remodelling processes. Atrial remodelling increases AF vulnerability and persistence, while promoting disease progression. The variability in presentation and wide range of mechanisms involved in initiation, maintenance and progression of AF, as well as its associated adverse outcomes, make the early identification of causal factors modifiable with therapeutic interventions challenging, likely contributing to suboptimal efficacy of current AF management. Computational modelling facilitates the multilevel integration of multiple datasets and offers new opportunities for mechanistic understanding, risk prediction and personalized therapy. Mathematical simulations of cardiac electrophysiology have been around for 60 years and are being increasingly used to improve our understanding of AF mechanisms and guide AF therapy. This narrative review focuses on the emerging and future applications of computational modelling in AF management. We summarize clinical challenges that may benefit from computational modelling, provide an overview of the different in silico approaches that are available together with their notable achievements, and discuss the major limitations that hinder the routine clinical application of these approaches. Finally, future perspectives are addressed. With the rapid progress in electronic technologies including computing, clinical applications of computational modelling are advancing rapidly. We expect that their application will progressively increase in prominence, especially if their added value can be demonstrated in clinical trials.

Toward Patient-Specific Prediction of Ablation Strategies for Atrial Fibrillation Using Deep Learning

Atrial fibrillation (AF) is a common cardiac arrhythmia that affects 1% of the population worldwide and is associated with high levels of morbidity and mortality. Catheter ablation (CA) has become one of the first line treatments for AF, but its success rates are suboptimal, especially in the case of persistent AF. Computational approaches have shown promise in predicting the CA strategy using simulations of atrial models, as well as applying deep learning to atrial images. We propose a novel approach that combines image-based computational modelling of the atria with deep learning classifiers trained on patient-specific atrial models, which can be used to assist in CA therapy selection. Therefore, we trained a deep convolutional neural network (CNN) using a combination of (i) 122 atrial tissue images obtained by unfolding patient LGE-MRI datasets, (ii) 157 additional synthetic images derived from the patient data to enhance the training dataset, and (iii) the outcomes of 558 CA simulations to terminate several AF scenarios in the corresponding image-based atrial models. Four CNN classifiers were trained on this patient-specific dataset balanced using several techniques to predict three common CA strategies from the patient atrial images: pulmonary vein isolation (PVI), rotor-based ablation (Rotor) and fibrosis-based ablation (Fibro). The training accuracy for these classifiers ranged from 96.22 to 97.69%, while the validation accuracy was from 78.68 to 86.50%. After training, the classifiers were applied to predict CA strategies for an unseen holdout test set of atrial images, and the results were compared to outcomes of the respective image-based simulations. The highest success rate was observed in the correct prediction of the Rotor and Fibro strategies (100%), whereas the PVI class was predicted in 33.33% of the cases. In conclusion, this study provides a proof-of-concept that deep neural networks can learn from patient-specific MRI datasets and image-derived models of AF, providing a novel technology to assist in tailoring CA therapy to a patient.

Using machine learning to identify local cellular properties that support re-entrant activation in patient-specific models of atrial fibrillation

AIMS

Atrial fibrillation (AF) is sustained by re-entrant activation patterns. Ablation strategies have been proposed that target regions of tissue that may support re-entrant activation patterns. We aimed to characterize the tissue properties associated with regions that tether re-entrant activation patterns in a validated virtual patient cohort.

METHODS AND RESULTS

Atrial fibrillation patient-specific models (seven paroxysmal and three persistent) were generated and validated against local activation time (LAT) measurements during an S1-S2 pacing protocol from the coronary sinus and high right atrium, respectively. Atrial models were stimulated with burst pacing from three locations in the proximity of each pulmonary vein to initiate re-entrant activation patterns. Five atria exhibited sustained activation patterns for at least 80 s. Models with short maximum action potential durations (APDs) were associated with sustained activation. Phase singularities were mapped across the atria sustained activation patterns. Regions with a low maximum conduction velocity (CV) were associated with tethering of phase singularities. A support vector machine (SVM) was trained on maximum local conduction velocity and action potential duration to identify regions that tether phase singularities. The SVM identified regions of tissue that could support tethering with 91% accuracy. This accuracy increased to 95% when the SVM was also trained on surface area.

CONCLUSION

In a virtual patient cohort, local tissue properties, that can be measured (CV) or estimated (APD; using effective refractory period as a surrogate) clinically, identified regions of tissue that tether phase singularities. Combing CV and APD with atrial surface area further improved the accuracy in identifying regions that tether phase singularities.

Comprehensive evaluation of electrophysiological and 3D structural features of human atrial myocardium with insights on atrial fibrillation maintenance mechanisms

Atrial fibrillation (AF) occurrence and maintenance is associated with progressive remodeling of electrophysiological (repolarization and conduction) and 3D structural (fibrosis, fiber orientations, and wall thickness) features of the human atria. Significant diversity in AF etiology leads to heterogeneous arrhythmogenic electrophysiological and structural substrates within the 3D structure of the human atria. Since current clinical methods have yet to fully resolve the patient-specific arrhythmogenic substrates, mechanism-based AF treatments remain underdeveloped. Here, we review current knowledge from in-vivo, ex-vivo, and in-vitro human heart studies, and discuss how these studies may provide new insights on the synergy of atrial electrophysiological and 3D structural features in AF maintenance. In-vitro studies on surgically acquired human atrial samples provide a great opportunity to study a wide spectrum of AF pathology, including functional changes in single-cell action potentials, ion channels, and gene/protein expression. However, limited size of the samples prevents evaluation of heterogeneous AF substrates and reentrant mechanisms. In contrast, coronary-perfused ex-vivo human hearts can be studied with state-of-the-art functional and structural technologies, such as high-resolution near-infrared optical mapping and contrast-enhanced MRI. These imaging modalities can resolve atrial arrhythmogenic substrates and their role in reentrant mechanisms maintaining AF and validate clinical approaches. Nonetheless, longitudinal studies are not feasible in explanted human hearts. As no approach is perfect, we suggest that combining the strengths of direct human atrial studies with high fidelity approaches available in the laboratory and in realistic patient-specific computer models would elucidate deeper knowledge of AF mechanisms. We propose that a comprehensive translational pipeline from ex-vivo human heart studies to longitudinal clinically relevant AF animal studies and finally to clinical trials is necessary to identify patient-specific arrhythmogenic substrates and develop novel AF treatments.

A survey of pathways for mechano-electric coupling in the atria

Mechano-electric coupling (MEC) in atrial tissue has received sparse investigation to date, despite the well-known association between chronic atrial dilation and atrial fibrillation (AF). Of note, no fewer than six different mechanisms pertaining to stretch-activated channels, cellular capacitance and geometric effects have been identified in the literature as potential players. In this mini review, we briefly survey each of these pathways to MEC. We then perform computational simulations using single cell and tissue models in presence of various stretch regimes and MEC pathways. This allows us to assess the relative significance of each pathway in determining action potential duration, conduction velocity and rotor stability. For chronic atrial stretch, we find that stretch-induced alterations in membrane capacitance decrease conduction velocity and increase action potential duration, in agreement with experimental findings. In the presence of time-dependent passive atrial stretch, stretch-activated channels play the largest role, leading to after-depolarizations and rotor hypermeandering. These findings suggest that physiological atrial stretches, such as passive stretch during the atrial reservoir phase, may play an important part in the mechanisms of atrial arrhythmogenesis.

•

Passive strains caused by ventricular contraction need to be considered when incorporating mechano-electro feedback in atrial electrophysiology models.

•

In chronic stretch, stretch-induced capacitance changes dominate.

•

Chronic stretch leads to an increase in action potential duration and a reduction in conduction velocity, consistent with experimental studies.

•

In the presence of passive stretch, stretch-activated channels can induce delayed after-depolarisations and lead to rotor hypermeandering.

•

Mechano-electro feedback is thus likely to have implications for the genesis and maintenance of atrial arrhythmias.

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

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

Identifying locations of re-entrant drivers from patient-specific distribution of fibrosis in the left atrium

Clinical evidence suggests a link between fibrosis in the left atrium (LA) and atrial fibrillation (AF), the most common sustained arrhythmia. Image-derived fibrosis is increasingly used for patient stratification and therapy guidance. However, locations of re-entrant drivers (RDs) sustaining AF are unknown and therapy success rates remain suboptimal. This study used image-derived LA models to explore the dynamics of RD stabilization in fibrotic regions and generate maps of RD locations. LA models with patient-specific geometry and fibrosis distribution were derived from late gadolinium enhanced magnetic resonance imaging of 6 AF patients. In each model, RDs were initiated at multiple locations, and their trajectories were tracked and overlaid on the LA fibrosis distributions to identify the most likely regions where the RDs stabilized. The simulations showed that the RD dynamics were strongly influenced by the amount and spatial distribution of fibrosis. In patients with fibrosis burden greater than 25%, RDs anchored to specific locations near large fibrotic patches. In patients with fibrosis burden below 25%, RDs either moved near small fibrotic patches or anchored to anatomical features. The patient-specific maps of RD locations showed that areas that harboured the RDs were much smaller than the entire fibrotic areas, indicating potential targets for ablation therapy. Ablating the predicted locations and connecting them to the existing pulmonary vein ablation lesions was the most effective in-silico ablation strategy.

Cardiac computational modelling

Cardiovascular diseases currently have a major social and economic impact, constituting one of the leading causes of mortality and morbidity. Personalized computational models of the heart are demonstrating their usefulness both to help understand the mechanisms underlying cardiac disease, and to optimize their treatment and predict the patient's response. Within this framework, the Spanish Research Network for Cardiac Computational Modelling (VHeart-SN) has been launched. The general objective of the VHeart-SN network is the development of an integrated, modular and multiscale multiphysical computational model of the heart. This general objective is addressed through the following specific objectives: a) to integrate the different numerical methods and models taking into account the specificity of patients; b) to assist in advancing knowledge of the mechanisms associated with cardiac and vascular diseases; and c) to support the application of different personalized therapies. This article presents the current state of cardiac computational modelling and different scientific works conducted by the members of the network to gain greater understanding of the characteristics and usefulness of these models.

Arrhythmogenic propensity of the fibrotic substrate after atrial fibrillation ablation: a longitudinal study using magnetic resonance imaging-based atrial models

AIMS

Inadequate modification of the atrial fibrotic substrate necessary to sustain re-entrant drivers (RDs) may explain atrial fibrillation (AF) recurrence following failed pulmonary vein isolation (PVI). Personalized computational models of the fibrotic atrial substrate derived from late gadolinium enhanced (LGE)-magnetic resonance imaging (MRI) can be used to non-invasively determine the presence of RDs. The objective of this study is to assess the changes of the arrhythmogenic propensity of the fibrotic substrate after PVI.

METHODS AND RESULTS

Pre- and post-ablation individualized left atrial models were constructed from 12 AF patients who underwent pre- and post-PVI LGE-MRI, in six of whom PVI failed. Pre-ablation AF sustained by RDs was induced in 10 models. RDs in the post-ablation models were classified as either preserved or emergent. Pre-ablation models derived from patients for whom the procedure failed exhibited a higher number of RDs and larger areas defined as promoting RD formation when compared with atrial models from patients who had successful ablation, 2.6 ± 0.9 vs. 1.8 ± 0.2 and 18.9 ± 1.6% vs. 13.8 ± 1.5%, respectively. In cases of successful ablation, PVI eliminated completely the RDs sustaining AF. Preserved RDs unaffected by ablation were documented only in post-ablation models of patients who experienced recurrent AF (2/5 models); all of these models had also one or more emergent RDs at locations distinct from those of pre-ablation RDs. Emergent RDs occurred in regions that had the same characteristics of the fibrosis spatial distribution (entropy and density) as regions that harboured RDs in pre-ablation models.

CONCLUSION

Recurrent AF after PVI in the fibrotic atria may be attributable to both preserved RDs that sustain AF pre- and post-ablation, and the emergence of new RDs following ablation. The same levels of fibrosis entropy and density underlie the pro-RD propensity in both pre- and post-ablation substrates.

An audit of uncertainty in multi-scale cardiac electrophysiology models

Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Epicardial Fibrosis Explains Increased Endo-Epicardial Dissociation and Epicardial Breakthroughs in Human Atrial Fibrillation

Background

Atrial fibrillation (AF) is accompanied by progressive epicardial fibrosis, dissociation of electrical activity between the epicardial layer and the endocardial bundle network, and transmural conduction (breakthroughs). However, causal relationships between these phenomena have not been demonstrated yet. Our goal was to test the hypothesis that epicardial fibrosis suffices to increase endo–epicardial dissociation (EED) and breakthroughs (BT) during AF.

Methods

We simulated the effect of fibrosis in the epicardial layer on EED and BT in a detailed, high-resolution, three-dimensional model of the human atria with realistic electrophysiology. The model results were compared with simultaneous endo–epicardial mapping in human atria. The model geometry, specifically built for this study, was based on MR images and histo-anatomical studies. Clinical data were obtained in four patients with longstanding persistent AF (persAF) and three patients without a history of AF.

Results

The AF cycle length (AFCL), conduction velocity (CV), and EED were comparable in the mapping studies and the simulations. EED increased from 24.1 ± 3.4 to 56.58 ± 6.2% (p < 0.05), and number of BTs per cycle from 0.89 ± 0.55 to 6.74 ± 2.11% (p < 0.05), in different degrees of fibrosis in the epicardial layer. In both mapping data and simulations, EED correlated with prevalence of BTs. Fibrosis also increased the number of fibrillation waves per cycle in the model.

Conclusion

A realistic 3D computer model of AF in which epicardial fibrosis was increased, in the absence of other pathological changes, showed increases in EED and epicardial BT comparable to those in longstanding persAF. Thus, epicardial fibrosis can explain both phenomena.

A robust computational framework for estimating 3D Bi-Atrial chamber wall thickness

Atrial fibrillation (AF) is the most prevalent form of cardiac arrhythmia. The atrial wall thickness (AWT) can potentially improve our understanding of the mechanism underlying atrial structure that drives AF and provides important clinical information. However, most existing studies for estimating AWT rely on ruler-based measurements performed on only a few selected locations in 2D or 3D using digital calipers. Only a few studies have developed automatic approaches to estimate the AWT in the left atrium, and there are currently no methods to robustly estimate the AWT of both atrial chambers. Therefore, we have developed a computational pipeline to automatically calculate the 3D AWT across bi-atrial chambers and extensively validated our pipeline on both ex vivo and in vivo human atria data. The atrial geometry was first obtained by segmenting the atrial wall from the MRIs using a novel machine learning approach. The epicardial and endocardial surfaces were then separated using a multi-planar convex hull approach to define boundary conditions, from which, a Laplace equation was solved numerically to automatically separate bi-atrial chambers. To robustly estimate the AWT in each atrial chamber, coupled partial differential equations by coupling the Laplace solution with two surface trajectory functions were formulated and solved. Our pipeline enabled the reconstruction and visualization of the 3D AWT for bi-atrial chambers with a relative error of 8% and outperformed existing algorithms by >7%. Our approach can potentially lead to improved clinical diagnosis, patient stratification, and clinical guidance during ablation treatment for patients with AF.

Heterogeneous Effects of Fibroblast-Myocyte Coupling in Different Regions of the Human Atria Under Conditions of Atrial Fibrillation

Background: Atrial fibrillation (AF), the most common cardiac arrhythmia, is characterized by alteration of the action potential (AP) propagation. Under persistent AF, myocytes undergo electrophysiological and structural remodeling, which involves fibroblast proliferation and differentiation, modifying the substrate for AP propagation. The aim of this study was to analyze the effects on the AP of fibroblast-myocyte coupling during AF and its propagation in different regions of the atria.

Methods: Isolated myocytes were coupled to different numbers of fibroblasts using the established AP models and tissue simulations were performed by randomly distributing fibroblasts. Fibroblast formulations were updated to match recent experimental data. Major ion current conductances of the myocyte model were modified to simulate AP heterogeneity in four different atrial regions (right atrium posterior wall, crista terminalis, left atrium posterior wall, and pulmonary vein) according to experimental and computational studies.

Results: The results of the coupled myocyte-fibroblast simulations suggest that a more depolarized membrane potential and higher fibroblast membrane capacitance have a greater impact on AP duration and myocyte maximum depolarization velocity. The number of coupled fibroblasts and the stimulation frequency are determining factors in altering myocyte AP. Strand simulations show that conduction velocity tends to homogenize in all regions, while the left atrium is more likely to be affected by fibroblast and AP propagation block is more likely to occur. The pulmonary vein is the most affected region, even at low fibroblast densities. In 2D sheets with randomly placed fibroblasts, wavebreaks are observed in the low density (10%) central fibrotic zone and when fibroblast density increases (40%) propagation in the fibrotic region is practically blocked. At densities of 10 and 20% the width of the vulnerable window increases with respect to control but is decreased at 40%.

Conclusion: Myocyte-fibroblast coupling characteristics heterogeneously affect AP propagation and features in the different atrial zones, and myocytes from the left atria are more sensitive to fibroblast coupling.

Understanding AF Mechanisms Through Computational Modelling and Simulations

AF is a progressive disease of the atria, involving complex mechanisms related to its initiation, maintenance and progression. Computational modelling provides a framework for integration of experimental and clinical findings, and has emerged as an essential part of mechanistic research in AF. The authors summarise recent advancements in development of multi-scale AF models and focus on the mechanistic links between alternations in atrial structure and electrophysiology with AF. Key AF mechanisms that have been explored using atrial modelling are pulmonary vein ectopy; atrial fibrosis and fibrosis distribution; atrial wall thickness heterogeneity; atrial adipose tissue infiltration; development of repolarisation alternans; cardiac ion channel mutations; and atrial stretch with mechano-electrical feedback. They review modelling approaches that capture variability at the cohort level and provide cohort-specific mechanistic insights. The authors conclude with a summary of future perspectives, as envisioned for the contributions of atrial modelling in the mechanistic understanding of AF.

The role of personalized atrial modeling in understanding atrial fibrillation mechanisms and improving treatment

Atrial fibrillation is the most common arrhythmia in humans and is associated with high morbidity, mortality and health-related expenses. Computational approaches have been increasingly utilized in atrial electrophysiology. In this review we summarize the recent advancements in atrial fibrillation modeling at the organ scale. Multi-scale atrial models now incorporate high level detail of atrial anatomy, tissue ultrastructure and fibrosis distribution. We provide the state-of-the art methodologies in developing personalized atrial fibrillation models with realistic geometry and tissue properties. We then focus on the use of multi-scale atrial models to gain mechanistic insights in AF. Simulations using atrial models have provided important insight in the mechanisms underlying AF, showing the importance of the atrial fibrotic substrate and altered atrial electrophysiology in initiation and maintenance of AF. Last, we summarize the translational evidence that supports incorporation of computational modeling in clinical practice for development of personalized treatment strategies for patients with AF. In early-stages clinical studies, AF models successfully identify patients where pulmonary vein isolation alone is not adequate for treatment of AF and suggest novel targets for ablation. We conclude with a summary of the future developments envisioned for the field of atrial computational electrophysiology.