Using personalized computer models to custom-tailor ablation procedures for atrial fibrillation patients: are we there yet?

A Review of Personalised Cardiac Computational Modelling Using Electroanatomical Mapping Data

Computational models of cardiac electrophysiology have gradually matured during the past few decades and are now being personalised to provide patient-specific therapy guidance for improving suboptimal treatment outcomes. The predictive features of these personalised electrophysiology models hold the promise of providing optimal treatment planning, which is currently limited in the clinic owing to reliance on a population-based or average patient approach. The generation of a personalised electrophysiology model entails a sequence of steps for which a range of activation mapping, calibration methods and therapy simulation pipelines have been suggested. However, the optimal methods that can potentially constitute a clinically relevant in silico treatment are still being investigated and face limitations, such as uncertainty of electroanatomical data recordings, generation and calibration of models within clinical timelines and requirements to validate or benchmark the recovered tissue parameters. This paper is aimed at reporting techniques on the personalisation of cardiac computational models, with a focus on calibrating cardiac tissue conductivity based on electroanatomical mapping data.

Towards Improved Human In Vitro Models for Cardiac Arrhythmia: Disease Mechanisms, Treatment, and Models of Atrial Fibrillation

Heart rhythm disorders, arrhythmias, place a huge economic burden on society and have a large impact on the quality of life of a vast number of people. Arrhythmias can have genetic causes but primarily arise from heart tissue remodeling during aging or heart disease. As current therapies do not address the causes of arrhythmias but only manage the symptoms, it is of paramount importance to generate innovative test models and platforms for gaining knowledge about the underlying disease mechanisms which are compatible with drug screening. In this review, we outline the most important features of atrial fibrillation (AFib), the most common cardiac arrhythmia. We will discuss the epidemiology, risk factors, underlying causes, and present therapies of AFib, as well as the shortcomings and opportunities of current models for cardiac arrhythmia, including animal models, in silico and in vitro models utilizing human pluripotent stem cell (hPSC)-derived cardiomyocytes.

A Review on Atrial Fibrillation (Computer Simulation and Clinical Perspectives)

Atrial fibrillation (AF), a heart condition, has been a well-researched topic for the past few decades. This multidisciplinary field of study deals with signal processing, finite element analysis, mathematical modeling, optimization, and clinical procedure. This article is focused on a comprehensive review of journal articles published in the field of AF. Topics from the age-old fundamental concepts to specialized modern techniques involved in today’s AF research are discussed. It was found that a lot of research articles have already been published in modeling and simulation of AF. In comparison to that, the diagnosis and post-operative procedures for AF patients have not yet been totally understood or explored by the researchers. The simulation and modeling of AF have been investigated by many researchers in this field. Cellular model, tissue model, and geometric model among others have been used to simulate AF. Due to a very complex nature, the causes of AF have not been fully perceived to date, but the simulated results are validated with real-life patient data. Many algorithms have been proposed to detect the source of AF in human atria. There are many ablation strategies for AF patients, but the search for more efficient ablation strategies is still going on. AF management for patients with different stages of AF has been discussed in the literature as well but is somehow limited mostly to the patients with persistent AF. The authors hope that this study helps to find existing research gaps in the analysis and the diagnosis of AF.

Synergistic antiarrhythmic effect of inward rectifier current inhibition and pulmonary vein isolation in a 3D computer model for atrial fibrillation

AIMS

Recent clinical studies showed that antiarrhythmic drug (AAD) treatment and pulmonary vein isolation (PVI) synergistically reduce atrial fibrillation (AF) recurrences after initially successful ablation. Among newly developed atrial-selective AADs, inhibitors of the G-protein-gated acetylcholine-activated inward rectifier current (IKACh) were shown to effectively suppress AF in an experimental model but have not yet been evaluated clinically. We tested in silico whether inhibition of inward rectifier current or its combination with PVI reduces AF inducibility.

METHODS AND RESULTS

We simulated the effect of inward rectifier current blockade (IK blockade), PVI, and their combination on AF inducibility in a detailed three-dimensional model of the human atria with different degrees of fibrosis. IK blockade was simulated with a 30% reduction of its conductivity. Atrial fibrillation was initiated using incremental pacing applied at 20 different locations, in both atria. IK blockade effectively prevented AF induction in simulations without fibrosis as did PVI in simulations without fibrosis and with moderate fibrosis. Both interventions lost their efficacy in severe fibrosis. The combination of IK blockade and PVI prevented AF in simulations without fibrosis, with moderate fibrosis, and even with severe fibrosis. The combined therapy strongly decreased the number of fibrillation waves, due to a synergistic reduction of wavefront generation rate while the wavefront lifespan remained unchanged.

CONCLUSION

Newly developed blockers of atrial-specific inward rectifier currents, such as IKAch, might prevent AF occurrences and when combined with PVI effectively supress AF recurrences in human.

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.

Arrhythmia dynamics in computational models of the atria following virtual ablation of re-entrant drivers

Aims

Efforts to improve ablation success rates in persistent atrial fibrillation (AF) patients by targeting re-entrant driver (RD) sites have been hindered by weak mechanistic understanding regarding emergent RDs localization following initial fibrotic substrate modification. This study aimed to systematically assess arrhythmia dynamics after virtual ablation of RD sites in computational models.

Methods and results

Simulations were conducted in 12 patient-specific atrial models reconstructed from pre-procedure late gadolinium-enhanced magnetic resonance imaging scans. In a previous study involving these same models, we comprehensively characterized pre-ablation RDs in simulations conducted with either 'average human AF'-based electrophysiology (i.e. EPavg) or ±10% action potential duration or conduction velocity (i.e. EPvar). Re-entrant drivers seen under the EPavg condition were virtually ablated and the AF initiation protocol was re-applied. Twenty-one emergent RDs were observed in 9/12 atrial models (1.75 ± 1.35 emergent RDs per model); these dynamically localized to boundary regions between fibrotic and non-fibrotic tissue. Most emergent RD locations (15/21, 71.4%) were within 0.1 cm of sites where RDs were seen pre-ablation in simulations under EPvar conditions. Importantly, this suggests that the level of uncertainty in our models' ability to predict patient-specific ablation targets can be substantially mitigated by running additional simulations that include virtual ablation of RDs. In 7/12 atrial models, at least one episode of macro-reentry around ablation lesion(s) was observed.

Conclusion

Arrhythmia episodes after virtual RD ablation are perpetuated by both emergent RDs and by macro-reentrant circuits formed around lesions. Custom-tailoring of ablation procedures based on models should take steps to mitigate these sources of AF recurrence.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Computational modeling: What does it tell us about atrial fibrillation therapy?

Atrial fibrillation (AF) is a complex cardiac arrhythmia with diverse etiology that negatively affects morbidity and mortality of millions of patients. Technological and experimental advances have provided a wealth of information on the pathogenesis of AF, highlighting a multitude of mechanisms involved in arrhythmia initiation and maintenance, and disease progression. However, it remains challenging to identify the predominant mechanisms for specific subgroups of AF patients, which, together with an incomplete understanding of the pleiotropic effects of antiarrhythmic therapies, likely contributes to the suboptimal efficacy of current antiarrhythmic approaches. Computer modeling of cardiac electrophysiology has advanced in parallel to experimental research and provides an integrative framework to attempt to overcome some of these challenges. Multi-scale cardiac modeling and simulation integrate structural and functional data from experimental and clinical work with knowledge of atrial electrophysiological mechanisms and dynamics, thereby improving our understanding of AF mechanisms and therapy. In this review, we describe recent advances in our quantitative understanding of AF through mathematical models. We discuss computational modeling of AF mechanisms and therapy using detailed, mechanistic cell/tissue-level models, including approaches to incorporate variability in patient populations. We also highlight efforts using whole-atria models to improve catheter ablation therapies. Finally, we describe recent efforts and suggest future extensions to model clinical concepts of AF using patient-level models.

Patient-Specific Identification of Atrial Flutter Vulnerability-A Computational Approach to Reveal Latent Reentry Pathways

Atypical atrial flutter (AFlut) is a reentrant arrhythmia which patients frequently develop after ablation for atrial fibrillation (AF). Indeed, substrate modifications during AF ablation can increase the likelihood to develop AFlut and it is clinically not feasible to reliably and sensitively test if a patient is vulnerable to AFlut. Here, we present a novel method based on personalized computational models to identify pathways along which AFlut can be sustained in an individual patient. We build a personalized model of atrial excitation propagation considering the anatomy as well as the spatial distribution of anisotropic conduction velocity and repolarization characteristics based on a combination of a priori knowledge on the population level and information derived from measurements performed in the individual patient. The fast marching scheme is employed to compute activation times for stimuli from all parts of the atria. Potential flutter pathways are then identified by tracing loops from wave front collision sites and constricting them using a geometric snake approach under consideration of the heterogeneous wavelength condition. In this way, all pathways along which AFlut can be sustained are identified. Flutter pathways can be instantiated by using an eikonal-diffusion phase extrapolation approach and a dynamic multifront fast marching simulation. In these dynamic simulations, the initial pattern eventually turns into the one driven by the dominant pathway, which is the only pathway that can be observed clinically. We assessed the sensitivity of the flutter pathway maps with respect to conduction velocity and its anisotropy. Moreover, we demonstrate the application of tailored models considering disease-specific repolarization properties (healthy, AF-remodeled, potassium channel mutations) as well as applicabiltiy on a clinical dataset. Finally, we tested how AFlut vulnerability of these substrates is modulated by exemplary antiarrhythmic drugs (amiodarone, dronedarone). Our novel method allows to assess the vulnerability of an individual patient to develop AFlut based on the personal anatomical, electrophysiological, and pharmacological characteristics. In contrast to clinical electrophysiological studies, our computational approach provides the means to identify all possible AFlut pathways and not just the currently dominant one. This allows to consider all relevant AFlut pathways when tailoring clinical ablation therapy in order to reduce the development and recurrence of AFlut.

The Fibrotic Substrate in Persistent Atrial Fibrillation Patients: Comparison Between Predictions From Computational Modeling and Measurements From Focal Impulse and Rotor Mapping

Focal impulse and rotor mapping (FIRM) involves intracardiac detection and catheter ablation of re-entrant drivers (RDs), some of which may contribute to arrhythmia perpetuation in persistent atrial fibrillation (PsAF). Patient-specific computational models derived from late gadolinium-enhanced magnetic resonance imaging (LGE-MRI) has the potential to non-invasively identify all areas of the fibrotic substrate where RDs could potentially be sustained, including locations where RDs may not manifest during mapped AF episodes. The objective of this study was to carry out multi-modal assessment of the arrhythmogenic propensity of the fibrotic substrate in PsAF patients by comparing locations of RD-harboring regions found in simulations and detected by FIRM (RD_sim and RD_FIRM) and analyze implications for ablation strategies predicated on targeting RDs. For 11 PsAF patients who underwent pre-procedure LGE-MRI and FIRM-guided ablation, we retrospectively simulated AF in individualized atrial models, with geometry and fibrosis distribution reconstructed from pre-ablation LGE-MRI scans, and identified RD_sim sites. Regions harboring RD_sim and RD_FIRM were compared. RD_sim were found in 38 atrial regions (median [inter-quartile range (IQR)] = 4 [3; 4] per model). RD_FIRM were identified and subsequently ablated in 24 atrial regions (2 [1; 3] per patient), which was significantly fewer than the number of RD_sim-harboring regions in corresponding models (p < 0.05). Computational modeling predicted RD_sim in 20 of 24 (83%) atrial regions identified as RD_FIRM-harboring during clinical mapping. In a large number of cases, we uncovered RD_sim-harboring regions in which RD_FIRM were never observed (18/22 regions that differed between the two modalities; 82%); we termed such cases “latent” RD_sim sites. During follow-up (230 [180; 326] days), AF recurrence occurred in 7/11 (64%) individuals. Interestingly, latent RD_sim sites were observed in all seven computational models corresponding to patients who experienced recurrent AF (2 [2; 2] per patient); in contrast, latent RD_sim sites were only discovered in two of four patients who were free from AF during follow-up (0.5 [0; 1.5] per patient; p < 0.05 vs. patients with AF recurrence). We conclude that substrate-based ablation based on computational modeling could improve outcomes.

Virtual medicine: Utilization of the advanced cardiac imaging patient avatar for procedural planning and facilitation

Advances in imaging technology have led to a paradigm shift from planning of cardiovascular procedures and surgeries requiring the actual patient in a "brick and mortar" hospital to utilization of the digitalized patient in the virtual hospital. Cardiovascular computed tomographic angiography (CCTA) and cardiovascular magnetic resonance (CMR) digitalized 3-D patient representation of individual patient anatomy and physiology serves as an avatar allowing for virtual delineation of the most optimal approaches to cardiovascular procedures and surgeries prior to actual hospitalization. Pre-hospitalization reconstruction and analysis of anatomy and pathophysiology previously only accessible during the actual procedure could potentially limit the intrinsic risks related to time in the operating room, cardiac procedural laboratory and overall hospital environment. Although applications are specific to areas of cardiovascular specialty focus, there are unifying themes related to the utilization of technologies. The virtual patient avatar computer can also be used for procedural planning, computational modeling of anatomy, simulation of predicted therapeutic result, printing of 3-D models, and augmentation of real time procedural performance. Examples of the above techniques are at various stages of development for application to the spectrum of cardiovascular disease processes, including percutaneous, surgical and hybrid minimally invasive interventions. A multidisciplinary approach within medicine and engineering is necessary for creation of robust algorithms for maximal utilization of the virtual patient avatar in the digital medical center. Utilization of the virtual advanced cardiac imaging patient avatar will play an important role in the virtual health care system. Although there has been a rapid proliferation of early data, advanced imaging applications require further assessment and validation of accuracy, reproducibility, standardization, safety, efficacy, quality, cost effectiveness, and overall value to medical care.

Sensitivity of reentrant driver localization to electrophysiological parameter variability in image-based computational models of persistent atrial fibrillation sustained by a fibrotic substrate

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia, causing morbidity and mortality in millions worldwide. The atria of patients with persistent AF (PsAF) are characterized by the presence of extensive and distributed atrial fibrosis, which facilitates the formation of persistent reentrant drivers (RDs, i.e., spiral waves), which promote fibrillatory activity. Targeted catheter ablation of RD-harboring tissues has shown promise as a clinical treatment for PsAF, but the outcomes remain sub-par. Personalized computational modeling has been proposed as a means of non-invasively predicting optimal ablation targets in individual PsAF patients, but it remains unclear how RD localization dynamics are influenced by inter-patient variability in the spatial distribution of atrial fibrosis, action potential duration (APD), and conduction velocity (CV). Here, we conduct simulations in computational models of fibrotic atria derived from the clinical imaging of PsAF patients to characterize the sensitivity of RD locations to these three factors. We show that RDs consistently anchor to boundaries between fibrotic and non-fibrotic tissues, as delineated by late gadolinium-enhanced magnetic resonance imaging, but those changes in APD/CV can enhance or attenuate the likelihood that an RD will anchor to a specific site. These findings show that the level of uncertainty present in patient-specific atrial models reconstructed without any invasive measurements (i.e., incorporating each individual's unique distribution of fibrotic tissue from medical imaging alongside an average representation of AF-remodeled electrophysiology) is sufficiently high that a personalized ablation strategy based on targeting simulation-predicted RD trajectories alone may not produce the desired result.