Towards personalized computational modelling of the fibrotic substrate for atrial arrhythmia

Influence of genetic mutations to atria vulnerability to atrial fibrillation: An in-silico 3D human atria study

BACKGROUND AND OBJECTIVE

Personalized 3D computer models of atria have been extensively implemented in the last yearsas a tool to facilitate the understanding of the mechanisms underlying different forms of arrhythmia, such as atrial fibrillation (AF). Meanwhile, genetic mutations acting on potassium channel dynamics were demonstrated to induce fibrillatory episodes in asymptomatic patients. This research study aims at assessing the effects and the atrial susceptibility to AF of three gain-of-function mutations - namely, KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M - associated with AF outbreaks, using highly detailed 3D atrial models with realistic wall thickness and heterogenous histological properties.

METHODS

The 3D atrial model was generated by reconstructing segmented anatomical structures from CT scans of an AF patient. Modified versions of the Courtemanche human atrial myocyte model were used to reproduce the electrophysiological activity of the WT and of the three mutant cells. Ectopic foci (EF) were simulated in sixteen locations across the atrial mesh using an S1-S2 protocol with two S2 basic cycle lengths (BCL) and eleven coupling intervals in order to induce arrhythmias.

RESULTS

The three genetic mutations at 3D level reduced the APD90. The KCNE3-V17M mutation provoked the highest shortening (55 % in RA and LA with respect to WT), followed by KCNH2 T895M (14 % in RA and 18 % LA with respect to WT)and KCNH2 T436M (7 % in RA and 9 % LA with respect to WT). The KCNE3-V17M mutation led to arrhythmia in 67 % of the cases simulated and in 94 % of ectopic foci considered, at S2 BCL equal to 100 ms. The KCNH2 T436M and KCNH2 T895M mutations increased the vulnerability to AF in a similar way, leading to arrhythmic episodes in 7 % of the simulated conditions, at S2 BCL set to 160 ms. Overall, 60 % of the arrhythmic events generated arise in the left atrium. Spiral waves, multiple rotors and disordered electrical pattern were elicited in the presence of the KCNE3-V17M mutation, exhibiting an instantaneous mean frequency of 7.6 Hz with a mean standard deviation of 1.12 Hz. The scroll waves induced in the presence of the KCNH2 T436M and KCNH2 T895M mutations showed steadiness and regularity with an instantaneous mean frequencies in the range of 4.9 - 5.1 Hz and a mean standard deviation within 0.19 - 0.53 Hz.

CONCLUSIONS

The pro-arrhythmogenicity of the KCNE3-V17M, KCNH2 T895M and KCNH2 T436M mutations was studied and proved on personalized 3D cardiac models. The three genetic mutations were demonstrated to increase the predisposition of atrial tissue to the formation of AF-susceptible substrate in different ways based on their effects on electrophysiological properties of the atria.

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.

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.

25 years of basic and translational science in EP Europace: novel insights into arrhythmia mechanisms and therapeutic strategies

In the last 25 years, EP Europace has published more than 300 basic and translational science articles covering different arrhythmia types (ranging from atrial fibrillation to ventricular tachyarrhythmias), different diseases predisposing to arrhythmia formation (such as genetic arrhythmia disorders and heart failure), and different interventional and pharmacological anti-arrhythmic treatment strategies (ranging from pacing and defibrillation to different ablation approaches and novel drug-therapies). These studies have been conducted in cellular models, small and large animal models, and in the last couple of years increasingly in silico using computational approaches. In sum, these articles have contributed substantially to our pathophysiological understanding of arrhythmia mechanisms and treatment options; many of which have made their way into clinical applications. This review discusses a representative selection of EP Europace manuscripts covering the topics of pacing and ablation, atrial fibrillation, heart failure and pro-arrhythmic ventricular remodelling, ion channel (dys)function and pharmacology, inherited arrhythmia syndromes, and arrhythmogenic cardiomyopathies, highlighting some of the advances of the past 25 years. Given the increasingly recognized complexity and multidisciplinary nature of arrhythmogenesis and continued technological developments, basic and translational electrophysiological research is key advancing the field. EP Europace aims to further increase its contribution to the discovery of arrhythmia mechanisms and the implementation of mechanism-based precision therapy approaches in arrhythmia management.

Atrial fibrillation: Insights from animal models, computational modeling, and clinical studies

Atrial fibrillation (AF) is the most common human arrhythmia, affecting millions of patients worldwide. A combination of risk factors and comorbidities results in complex atrial remodeling, which increases AF vulnerability and persistence. Insights from animal models, clinical studies, and computational modeling have advanced the understanding of the mechanisms and pathophysiology of AF. Areas of heterogeneous pathological remodeling, as well as altered electrophysiological properties, serve as a substrate for AF drivers and spontaneous activations. The complex and individualized presentation of this arrhythmia suggests that mechanisms-based personalized approaches will likely be needed to overcome current challenges in AF management. In this paper, we review the insights on the mechanisms of AF obtained from animal models and clinical studies and how computational models integrate this knowledge to advance AF clinical management. We also assess the challenges that need to be overcome to implement these mechanistic models in clinical practice.

Thermal ablation effects on rotors that characterize functional re-entry cardiac arrhythmia

Thermal ablation is a well-established successful treatment for cardiac arrhythmia, but it still presents limitations that require further studies and developments. In the rotor-driven functional re-entry arrhythmia, tissue heterogeneity results on the generation of spiral/scroll waves and wave break dynamics that may cause dangerous sustainable fibrillation. The selection of the target region to perform thermal ablation to mitigate this type of arrhythmia is challenging, since it considerably affects the local electrophysiology dynamics. This work deals with the numerical simulation of the thermal ablation of a cardiac muscle tissue and its effects on the dynamics of rotor-driven functional re-entry arrhythmia. A non-homogeneous two-dimensional rectangular region is used in the present numerical analysis, where radiofrequency ablation is performed. The electrophysiology problem for the propagation of the action potential in the cardiac tissue is simulated with the Fenton-Karma model. Thermal damage caused to the tissue by the radiofrequency heating is modeled by the Arrhenius equation. The effects of size and position of a heterogeneous region in the original muscle tissue were first analyzed, in order to verify the possible existence of the functional re-entry arrhythmia during the time period considered in the simulations. For each case that exhibited re-entry arrhythmia, six different ablation procedures were analyzed, depending on the position of the radiofrequency electrode and heating time. The obtained results revealed the effects of different model parameters on the existence and possible mitigation of the functional re-entry arrhythmia.

Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions

Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100  μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.

Individual Contributions of Cardiac Ion Channels on Atrial Repolarization and Reentrant Waves: A Multiscale In-Silico Study

The excitation, contraction, and relaxation of an atrial cardiomyocyte are maintained by the activation and inactivation of numerous cardiac ion channels. Their collaborative efforts cause time-dependent changes of membrane potential, generating an action potential (AP), which is a surrogate marker of atrial arrhythmias. Recently, computational models of atrial electrophysiology emerged as a modality to investigate arrhythmia mechanisms and to predict the outcome of antiarrhythmic therapies. However, the individual contribution of atrial ion channels on atrial action potential and reentrant arrhythmia is not yet fully understood. Thus, in this multiscale in-silico study, perturbations of individual atrial ionic currents (I_Na, I_to, I_CaL, I_Kur, I_Kr, I_Ks, I_K1, I_NCX and I_NaK) in two in-silico models of human atrial cardiomyocyte (i.e., Courtemanche-1998 and Grandi-2011) were performed at both cellular and tissue levels. The results show that the inhibition of I_CaL and I_NCX resulted in AP shortening, while the inhibition of I_Kur, I_Kr, I_Ks, I_K1 and I_NaK prolonged AP duration (APD). Particularly, in-silico perturbations (inhibition and upregulation) of I_Kr and I_Ks only minorly affected atrial repolarization in the Grandi model. In contrast, in the Courtemanche model, the inhibition of I_Kr and I_Ks significantly prolonged APD and vice versa. Additionally, a 50% reduction of I_to density abbreviated APD in the Courtemanche model, while the same perturbation prolonged APD in the Grandi model. Similarly, a strong model dependence was also observed at tissue scale, with an observable I_K1-mediated reentry stabilizing effect in the Courtemanche model but not in the Grandi atrial model. Moreover, the Grandi model was highly sensitive to a change on intracellular Ca²⁺ concentration, promoting a repolarization failure in I_CaL upregulation above 150% and facilitating reentrant spiral waves stabilization by I_CaL inhibition. Finally, by incorporating the previously published atrial fibrillation (AF)-associated ionic remodeling in the Courtemanche atrial model, in-silico modeling revealed the antiarrhythmic effect of I_Kr inhibition in both acute and chronic settings. Overall, our multiscale computational study highlights the strong model-dependent effects of ionic perturbations which could affect the model’s accuracy, interpretability, and prediction. This observation also suggests the need for a careful selection of in-silico models of atrial electrophysiology to achieve specific research aims.

Simulation of diffuse and stringy fibrosis in a bilayer interconnected cable model of the left atrium

AIMS

The aim of this study is to design a computer model of the left atrium for investigating fibre-orientation-dependent microstructure such as stringy fibrosis.

METHODS AND RESULTS

We developed an approach for automatic construction of bilayer interconnected cable models from left atrial geometry and epi- and endocardial fibre orientation. The model consisted of two layers (epi- and endocardium) of longitudinal and transverse cables intertwined-like fabric threads, with a spatial discretization of 100 µm. Model validation was performed by comparison with cubic volumetric models in normal conditions. Then, diffuse (n = 2904), stringy (n = 3600), and mixed fibrosis patterns (n = 6840) were randomly generated by uncoupling longitudinal and transverse connections in the interconnected cable model. Fibrosis density was varied from 0% to 40% and mean stringy obstacle length from 0.1 to 2 mm. Total activation time, apparent anisotropy ratio, and local activation time jitter were computed during normal rhythm in each pattern. Non-linear regression formulas were identified for expressing measured propagation parameters as a function of fibrosis density and obstacle length (stringy and mixed patterns). Longer obstacles (even below tissue space constant) were independently associated with prolonged activation times, increased anisotropy, and local fluctuations in activation times. This effect was increased by endo-epicardial dissociation and mitigated when fibrosis was limited to the epicardium.

CONCLUSION

Interconnected cable models enable the study of microstructure in organ-size models despite limitations in the description of transmural structures.

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

AIMS

Computationally guided persistent atrial fibrillation (PsAF) ablation has emerged as an alternative to conventional treatment planning. To make this approach scalable, computational cost and the time required to conduct simulations must be minimized while maintaining predictive accuracy. Here, we assess the sensitivity of the process to finite-element mesh resolution. We also compare methods for pacing site distribution used to evaluate inducibility arrhythmia sustained by re-entrant drivers (RDs).

METHODS AND RESULTS

Simulations were conducted in low- and high-resolution models (average edge lengths: 400/350 µm) reconstructed from PsAF patients' late gadolinium enhancement magnetic resonance imaging scans. Pacing was simulated from 80 sites to assess RD inducibility. When pacing from the same site led to different outcomes in low-/high-resolution models, we characterized divergence dynamics by analysing dissimilarity index over time. Pacing site selection schemes prioritizing even spatial distribution and proximity to fibrotic tissue were evaluated. There were no RD sites observed in low-resolution models but not high-resolution models, or vice versa. Dissimilarity index analysis suggested that differences in simulation outcome arising from differences in discretization were the result of isolated conduction block incidents in one model but not the other; this never led to RD sites unique to one mesh resolution. Pacing site selection based on fibrosis proximity led to the best observed trade-off between number of stimulation locations and predictive accuracy.

CONCLUSION

Simulations conducted in meshes with 400 µm average edge length and ∼40 pacing sites proximal to fibrosis are sufficient to reveal the most comprehensive possible list of RD sites, given feasibility constraints.

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.

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.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

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.

Specific Electrogram Characteristics Identify the Extra-Pulmonary Vein Arrhythmogenic Sources of Persistent Atrial Fibrillation - Characterization of the Arrhythmogenic Electrogram Patterns During Atrial Fibrillation and Sinus Rhythm

Identification of atrial sites that perpetuate atrial fibrillation (AF), and ablation thereof terminates AF, is challenging. We hypothesized that specific electrogram (EGM) characteristics identify AF-termination sites (AFTS). Twenty-one patients in whom low-voltage-guided ablation after pulmonary vein isolation terminated clinical persistent AF were included. Patients were included if short RF-delivery for <8sec at a given atrial site was associated with acute termination of clinical persistent AF. EGM-characteristics at 21 AFTS, 105 targeted sites without termination and 105 non-targeted control sites were analyzed. Alteration of EGM-characteristics by local fibrosis was evaluated in a three-dimensional high resolution (100 µm)-computational AF model. AFTS demonstrated lower EGM-voltage, higher EGM-cycle-length-coverage, shorter AF-cycle-length and higher pattern consistency than control sites (0.49 ± 0.39 mV vs. 0.83 ± 0.76 mV, p < 0.0001; 79 ± 16% vs. 59 ± 22%, p = 0.0022; 173 ± 49 ms vs. 198 ± 34 ms, p = 0.047; 80% vs. 30%, p < 0.01). Among targeted sites, AFTS had higher EGM-cycle-length coverage, shorter local AF-cycle-length and higher pattern consistency than targeted sites without AF-termination (79 ± 16% vs. 63 ± 23%, p = 0.02; 173 ± 49 ms vs. 210 ± 44 ms, p = 0.002; 80% vs. 40%, p = 0.01). Low voltage (0.52 ± 0.3 mV) fractionated EGMs (79 ± 24 ms) with delayed components in sinus rhythm (‘atrial late potentials’, respectively ‘ALP’) were observed at 71% of AFTS. EGMs recorded from fibrotic areas in computational models demonstrated comparable EGM-characteristics both in simulated AF and sinus rhythm. AFTS may therefore be identified by locally consistent, fractionated low-voltage EGMs with high cycle-length-coverage and rapid activity in AF, with low-voltage, fractionated EGMs with delayed components/ ‘atrial late potentials’ (ALP) persisting in sinus rhythm.

Association of right atrial structure with incident atrial fibrillation: a longitudinal cohort cardiovascular magnetic resonance study from the Multi-Ethnic Study of Atherosclerosis (MESA)

BACKGROUND

While studies of the left atrium (LA) have demonstrated associations between volumes and emptying fraction with atrial fibrillation (AF), the contribution of right atrial (RA) abnormalities to incident AF remains poorly understood.

OBJECTIVES

Assess the association between RA structure and function with incident AF using feature-tracking cardiovascular magnetic resonance (CMR).

METHODS

This is a prospective cohort study of all participants in the Multi-Ethnic Study of Atherosclerosis with baseline CMR, sinus rhythm, and free of clinical cardiovascular disease at study initiation. RA volume, strain, and emptying fraction in participants with incident AF (n = 368) were compared against AF-free (n = 2779). Cox proportional-hazards models assessed association between variables.

RESULTS

Participants were aged 60 ± 10 yrs., 55% female, and followed an average 11.2 years. Individuals developing AF had higher baseline RA maximum volume index (mean ± standard deviation [SD]: 24 ± 9 vs 22 ± 8 mL/m2, p = 0.002) and minimum volume index (13 ± 7 vs 12 ± 6 mL/m2, p < 0.001), and lower baseline RA emptying fraction (45 ± 15% vs 47 ± 15%, p = 0.02), peak global strain (34 ± 17% vs 36 ± 19%, p < 0.001), and peak free-wall strain (40 ± 23% vs 42 ± 26%, p = 0.049) compared with the AF-free population. After adjusting for traditional cardiovascular risk factors and LA volume and function, we found RA maximum volume index (hazards ratio [HR]: 1.13 per SD, p = 0.041) and minimum volume index (HR: 1.12 per SD, p = 0.037) were independently associated with incident AF.

CONCLUSIONS

In a large multiethnic population, higher RA volume indices were independently associated with incident AF after adjustment for conventional cardiovascular risk factors and LA parameters. It is unclear if this predictive value persists when additional adjustment is made for ventricular parameters.

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.

Computationally guided personalized targeted ablation of persistent atrial fibrillation

Atrial fibrillation (AF) — the most common arrhythmia — significantly increases the risk of stroke and heart failure. Although catheter ablation can restore normal heart rhythms, patients with persistent AF who develop atrial fibrosis often undergo multiple failed ablations and thus increased procedural risks. Here, we present personalized computational modelling for the reliable predetermination of ablation targets, which are then used to guide the ablation procedure in patients with persistent AF and atrial fibrosis. We first show that a computational model of the atria of patients identifies fibrotic tissue that if ablated will not sustain AF. We then integrated the target-ablation sites in a clinical-mapping system, and tested its feasibility in 10 patients with persistent AF. The computational prediction of ablation targets avoids lengthy electrical mapping and could improve the accuracy and efficacy of targeted AF ablation in patients whilst eliminating the need for repeat procedures.

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.

A comprehensive, multiscale framework for evaluation of arrhythmias arising from cell therapy in the whole post-myocardial infarcted heart

Direct remuscularization approaches to cell-based heart repair seek to restore ventricular contractility following myocardial infarction (MI) by introducing new cardiomyocytes (CMs) to replace lost or injured ones. However, despite promising improvements in cardiac function, high incidences of ventricular arrhythmias have been observed in animal models of MI injected with pluripotent stem cell-derived cardiomyocytes (PSC-CMs). The mechanisms of arrhythmogenesis remain unclear. Here, we present a comprehensive framework for computational modeling of direct remuscularization approaches to cell therapy. Our multiscale 3D whole-heart modeling framework integrates realistic representations of cell delivery and transdifferentiation therapy modalities as well as representation of spatial distributions of engrafted cells, enabling simulation of clinical therapy and the prediction of emergent electrophysiological behavior and arrhythmogenensis. We employ this framework to explore how varying parameters of cell delivery and transdifferentiation could result in three mechanisms of arrhythmogenesis: focal ectopy, heart block, and reentry.

The Association Between Diabetes Mellitus and Atrial Fibrillation: Clinical and Mechanistic Insights

A number of clinical studies have reported that diabetes mellitus (DM) is an independent risk factor for Atrial fibrillation (AF). After adjustment for other known risk factors including age, sex, and cardiovascular risk factors, DM remains a significant if modest risk factor for development of AF. The mechanisms underlying the increased susceptibility to AF in DM are incompletely understood, but are thought to involve electrical, structural, and autonomic remodeling in the atria. Electrical remodeling in DM may involve alterations in gap junction function that affect atrial conduction velocity due to changes in expression or localization of connexins. Electrical remodeling can also occur due to changes in atrial action potential morphology in association with changes in ionic currents, such as sodium or potassium currents, that can affect conduction velocity or susceptibility to triggered activity. Structural remodeling in DM results in atrial fibrosis, which can alter conduction patterns and susceptibility to re-entry in the atria. In addition, increases in atrial adipose tissue, especially in Type II DM, can lead to disruptions in atrial conduction velocity or conduction patterns that may affect arrhythmogenesis. Whether the insulin resistance in type II DM activates unique intracellular signaling pathways independent of obesity requires further investigation. In addition, the relationship between incident AF and glycemic control requires further study.

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.

Simulated P wave morphology in the presence of endo-epicardial activation delay

Aims

Evidences of asynchrony between epicardial and endocardial activation in the atrial wall have been reported. We used a computer model of the atria and torso to investigate the consequences of such activation delay on P wave morphology, while controlling for P wave duration.

Methods and results

We created 390 models of the atria based on the same geometry. These models differed by atrial wall thickness (from 2 to 3 mm), transmural coupling, and tissue conductivity in the endocardial and epicardial layers. Among them, 18 were in baseline, 186 had slower conduction in the epicardium layer and 186 in the endocardial layer. Conduction properties were adjusted in such a way that total activation time was the same in all models. P waves on a 16-lead system were simulated during sinus rhythm. Activation maps were similar in all cases. Endo-epicardial delay varied between -5.5 and 5.5 ms vs. 0 ± 0.5 ms in baseline. All P waves had the same duration but variability in their morphology was observed. With slower epicardial conduction, P wave amplitude was reduced by an average of 20% on leads V3-V5 and P wave area decreased by 50% on leads V1-V2 and by 40% on lead V3. Reversed, lower magnitude effects were observed with slower endocardial conduction.

Conclusion

An endo-epicardial delay of a few milliseconds is sufficient to significantly alter P wave morphology, even if the activation map remains the same.

Computational Modeling of Electrophysiology and Pharmacotherapy of Atrial Fibrillation: Recent Advances and Future Challenges

The pathophysiology of atrial fibrillation (AF) is broad, with components related to the unique and diverse cellular electrophysiology of atrial myocytes, structural complexity, and heterogeneity of atrial tissue, and pronounced disease-associated remodeling of both cells and tissue. A major challenge for rational design of AF therapy, particularly pharmacotherapy, is integrating these multiscale characteristics to identify approaches that are both efficacious and independent of ventricular contraindications. Computational modeling has long been touted as a basis for achieving such integration in a rapid, economical, and scalable manner. However, computational pipelines for AF-specific drug screening are in their infancy, and while the field is progressing quite rapidly, major challenges remain before computational approaches can fill the role of workhorse in rational design of AF pharmacotherapies. In this review, we briefly detail the unique aspects of AF pathophysiology that determine requirements for compounds targeting AF rhythm control, with emphasis on delimiting mechanisms that promote AF triggers from those providing substrate or supporting reentry. We then describe modeling approaches that have been used to assess the outcomes of drugs acting on established AF targets, as well as on novel promising targets including the ultra-rapidly activating delayed rectifier potassium current, the acetylcholine-activated potassium current and the small conductance calcium-activated potassium channel. Finally, we describe how heterogeneity and variability are being incorporated into AF-specific models, and how these approaches are yielding novel insights into the basic physiology of disease, as well as aiding identification of the important molecular players in the complex AF etiology.

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.

Atrial Fibrosis Hampers Non-invasive Localization of Atrial Ectopic Foci From Multi-Electrode Signals: A 3D Simulation Study

Introduction: Focal atrial tachycardia is commonly treated by radio frequency ablation with an acceptable long-term success. Although the location of ectopic foci tends to appear in specific hot-spots, they can be located virtually in any atrial region. Multi-electrode surface ECG systems allow acquiring dense body surface potential maps (BSPM) for non-invasive therapy planning of cardiac arrhythmia. However, the activation of the atria could be affected by fibrosis and therefore biomarkers based on BSPM need to take these effects into account. We aim to analyze the effect of fibrosis on a BSPM derived index, and its potential application to predict the location of ectopic foci in the atria.

Methodology: We have developed a 3D atrial model that includes 5 distributions of patchy fibrosis in the left atrium at 5 different stages. Each stage corresponds to a different amount of fibrosis that ranges from 2 to 40%. The 25 resulting 3D models were used for simulation of Focal Atrial Tachycardia (FAT), triggered from 19 different locations described in clinical studies. BSPM were obtained for all simulations, and the body surface potential integral maps (BSPiM) were calculated to describe atrial activations. A machine learning (ML) pipeline using a supervised learning model and support vector machine was developed to learn the BSPM patterns of each of the 475 activation sequences and relate them to the origin of the FAT source.

Results: Activation maps for stages with more than 15% of fibrosis were greatly affected, producing conduction blocks and delays in propagation. BSPiMs did not always cluster into non-overlapped groups since BSPiMs were highly altered by the conduction blocks. From stage 3 (15% fibrosis) the BSPiMs showed differences for ectopic beats placed around the area of the pulmonary veins. Classification results were mostly above 84% for all the configurations studied when a large enough number of electrodes were used to map the torso. However, the presence of fibrosis increases the area of the ectopic focus location and therefore decreases the utility for the electrophysiologist.

Conclusions: The results indicate that the proposed ML pipeline is a promising methodology for non-invasive ectopic foci localization from BSPM signal even when fibrosis is present.

Comparing Reentrant Drivers Predicted by Image-Based Computational Modeling and Mapped by Electrocardiographic Imaging in Persistent Atrial Fibrillation

Electrocardiographic mapping (ECGI) detects reentrant drivers (RDs) that perpetuate arrhythmia in persistent AF (PsAF). Patient-specific computational models derived from late gadolinium-enhanced magnetic resonance imaging (LGE-MRI) identify all latent sites in the fibrotic substrate that could potentially sustain RDs, not just those manifested during mapped AF. The objective of this study was to compare RDs from simulations and ECGI (RD_sim/RD_ECGI) and analyze implications for ablation. We considered 12 PsAF patients who underwent RD_ECGI ablation. For the same cohort, we simulated AF and identified RD_sim sites in patient-specific models with geometry and fibrosis distribution from pre-ablation LGE-MRI. RD_sim- and RD_ECGI-harboring regions were compared, and the extent of agreement between macroscopic locations of RDs identified by simulations and ECGI was assessed. Effects of ablating RD_ECGI/RD_sim were analyzed. RD_sim were predicted in 28 atrial regions (median [inter-quartile range (IQR)] = 3.0 [1.0; 3.0] per model). ECGI detected 42 RD_ECGI-harboring regions (4.0 [2.0; 5.0] per patient). The number of regions with RD_sim and RD_ECGI per individual was not significantly correlated (R = 0.46, P = ns). The overall rate of regional agreement was fair (modified Cohen's κ₀ statistic = 0.11), as expected, based on the different mechanistic underpinning of RD_sim- and RD_ECGI. nineteen regions were found to harbor both RD_sim and RD_ECGI, suggesting that a subset of clinically observed RDs was fibrosis-mediated. The most frequent source of differences (23/32 regions) between the two modalities was the presence of RD_ECGI perpetuated by mechanisms other than the fibrotic substrate. In 6/12 patients, there was at least one region where a latent RD was observed in simulations but was not manifested during clinical mapping. Ablation of fibrosis-mediated RD_ECGI (i.e., targets in regions that also harbored RD_sim) trended toward a higher rate of positive response compared to ablation of other RD_ECGI targets (57 vs. 41%, P = ns). Our analysis suggests that RDs in human PsAF are at least partially fibrosis-mediated. Substrate-based ablation combining simulations with ECGI 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.

Efficient parameterization of cardiac action potential models using a genetic algorithm

Finding appropriate values for parameters in mathematical models of cardiac cells is a challenging task. Here, we show that it is possible to obtain good parameterizations in as little as 30-40 s when as many as 27 parameters are fit simultaneously using a genetic algorithm and two flexible phenomenological models of cardiac action potentials. We demonstrate how our implementation works by considering cases of "model recovery" in which we attempt to find parameter values that match model-derived action potential data from several cycle lengths. We assess performance by evaluating the parameter values obtained, action potentials at fit and non-fit cycle lengths, and bifurcation plots for fidelity to the truth as well as consistency across different runs of the algorithm. We also fit the models to action potentials recorded experimentally using microelectrodes and analyze performance. We find that our implementation can efficiently obtain model parameterizations that are in good agreement with the dynamics exhibited by the underlying systems that are included in the fitting process. However, the parameter values obtained in good parameterizations can exhibit a significant amount of variability, raising issues of parameter identifiability and sensitivity. Along similar lines, we also find that the two models differ in terms of the ease of obtaining parameterizations that reproduce model dynamics accurately, most likely reflecting different levels of parameter identifiability for the two models.

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.

Fibrosis and Atrial Fibrillation: Computerized and Optical Mapping; A View into the Human Atria at Submillimeter Resolution

Recent studies strongly suggest that the majority of atrial fibrillation (AF) patients with diagnosed or subclinical cardiac diseases have established or even pre-existing fibrotic structural remodeling, which may lead to conduction abnormalities and reentrant activity that sustain AF. As conventional treatments fail to treat AF in far too many cases, an urgent need exists to identify specific structural arrhythmogenic fibrosis patterns, which may maintain AF, in order to identify effective ablation targets for AF treatment. However, the existing challenge is to define what exact structural remodeling within the complex 3D human atrial wall is arrhythmogenic, as well as linking arrhythmogenic fibrosis to an underlying mechanism of AF maintenance in the clinical setting. This review is focused on the role of 3D fibrosis architecture in the mechanisms of AF maintenance revealed by submillimeter, high-resolution ex-vivo imaging modalities directly of human atria, as well as from in-silico 3D computational techniques that can be able to overcome in-vivo clinical limitations. The systematic integration of functional and structural imaging ex-vivo may inform the necessary integration of electrode and structural mapping in-vivo. A holistic view of AF driver mechanisms may begin to identify the defining characteristics or "fingerprints" of reentrant AF drivers, such as 3D fibrotic architecture, in order to design optimal patient-specific ablation strategies.

Demystifying rotors and their place in clinical translation of atrial fibrillation mechanisms

Treatment of atrial fibrillation (AF), the most common arrhythmia in clinical practice, remains challenging. Improved understanding of underlying mechanisms is needed to improve therapy. Functional re-entry is central to AF maintenance. The first detailed, quantitative theory of functional re-entry, the 'leading circle' model, was developed 40 years ago. Subsequently, an alternative paradigm based on 'spiral waves' has evolved. Spiral-wave generators, or 'rotors', have been identified using advanced mapping methods in experimental and clinical AF. A central tool in the analysis of spiral-wave rotors is the phase transformation, allowing for easier visualization of rotors and tracking of 'phase singularity' points at the rotor tip. In contrast to leading circle theory, which is expressed in terms familiar to (and easily understood by) cardiologists, the ideas needed to understand rotors are much more theoretical and harder for clinicians to apply. In this Review, we summarize the basic notions of phase mapping and spiral-wave rotors, and the ways in which rotor sources might be involved in AF maintenance. We discuss competing observations about the role of spatially confined rotors, short-lived rotors clustered at the edge of fibrotic zones, endocardial-epicardial interactive breeder properties and transmural re-entry, as well as studies underway to resolve them. We conclude with consideration of the clinical relevance of the issues discussed and their potential implications for the management of patients with AF.