Wavelength and Fibrosis Affect Phase Singularity Locations During 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.

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.

A novel approach for quantitative electrogram analysis for driver identification: Implications for ablation in persistent atrial fibrillation

Objective

This study sought to study the feasibility, efficacy, and safety of using multiscale entropy (MSE) analysis to guide catheter ablation for persistent atrial fibrillation (PsAF) and predict ablation outcomes.

Methods

We prospectively enrolled 108 patients undergoing initial ablation for PsAF. MSE was calculated based on bipolar intracardiac electrograms (iEGMs) to measure the dynamical complexity of biological signals. The iEGMs data were exported after pulmonary vein isolation (PVI), then calculated in a customed platform, and finally re-annotated into the CARTO system. After PVI, regions of the highest mean MSE (mMSE) values were ablated in descending order until AF termination, or three areas had been ablated.

Results

Baseline characteristics were evenly distributed between the AF termination (n = 38, 35.19%) and the non-termination group. The RA-to-LA mean MSE (mMSE) gradient demonstrated a positive gradient in the non-termination group and a negative gradient in the termination group (0.105 ± 0.180 vs. −0.235 ± 0.256, P < 0.001). During a 12-month follow-up, 29 patients (26.9%) had arrhythmia recurrence after single ablation, and 18 of them had AF (62.1%). The termination group had lower rates of arrhythmia recurrence (15.79 vs. 32.86%, Log-Rank P = 0.053) and AF recurrence (10.53 vs. 20%, Log-Rank P = 0.173) after single ablation and a lower rate of arrhythmia recurrence (7.89 vs. 27.14%, Log-Rank P = 0.018) after repeated ablation. Correspondingly, subjects with negative RA-to-LA mMSE gradient had lower incidences of arrhythmia (16.67 vs. 35%, Log-Rank P = 0.028) and AF (16.67 vs. 35%, Log-Rank P = 0.032) recurrence after single ablation and arrhythmia recurrence after repeated ablation (12.5 vs. 26.67%, Log-Rank P = 0.062). Marginal peri-procedural safety outcomes were observed.

Conclusion

MSE analysis-guided driver ablation in addition to PVI for PsAF could be feasible, efficient, and safe. An RA < LA mMSE gradient before ablation could predict freedom from arrhythmia. The RA-LA MSE gradient could be useful for guiding ablation strategy selection.

Current progress of computational modeling for guiding clinical atrial fibrillation ablation

Atrial fibrillation (AF) is one of the most common arrhythmias, associated with high morbidity, mortality, and healthcare costs, and it places a significant burden on both individuals and society. Anti-arrhythmic drugs are the most commonly used strategy for treating AF. However, drug therapy faces challenges because of its limited efficacy and potential side effects. Catheter ablation is widely used as an alternative treatment for AF. Nevertheless, because the mechanism of AF is not fully understood, the recurrence rate after ablation remains high. In addition, the outcomes of ablation can vary significantly between medical institutions and patients, especially for persistent AF. Therefore, the issue of which ablation strategy is optimal is still far from settled. Computational modeling has the advantages of repeatable operation, low cost, freedom from risk, and complete control, and is a useful tool for not only predicting the results of different ablation strategies on the same model but also finding optimal personalized ablation targets for clinical reference and even guidance. This review summarizes three-dimensional computational modeling simulations of catheter ablation for AF, from the early-stage attempts such as Maze III or circumferential pulmonary vein isolation to the latest advances based on personalized substrate-guided ablation. Finally, we summarize current developments and challenges and provide our perspectives and suggestions for future directions.

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.

Atrial fibrosis and substrate based characterization in atrial fibrillation: Time to move forwards

Atrial fibrillation (AF) is the most commonly encountered cardiac arrhythmia in clinical practice. However, current therapeutic interventions for atrial fibrillation have limited clinical efficacy as a consequence of major knowledge gaps in the mechanisms sustaining atrial fibrillation. From a mechanistic perspective, there is increasing evidence that atrial fibrosis plays a central role in the maintenance and perpetuation of atrial fibrillation. Electrophysiologically, atrial fibrosis results in alterations in conduction velocity, cellular refractoriness, and produces conduction block promoting meandering, unstable wavelets and micro-reentrant circuits. Clinically, atrial fibrosis has also linked to poor clinical outcomes including AF-related thromboembolic complications and arrhythmia recurrences post catheter ablation. In this article, we review the pathophysiology behind the formation of fibrosis as AF progresses, the role of fibrosis in arrhythmogenesis, surrogate markers for detection of fibrosis using cardiac magnetic resonance imaging, echocardiography and electroanatomic mapping, along with their respective limitations. We then proceed to review the current evidence behind therapeutic interventions targeting atrial fibrosis, including drugs and substrate-based catheter ablation therapies followed by the potential future use of electro phenotyping for AF characterization to overcome the limitations of contemporary substrate-based methodologies.

Theoretical Models and Computational Analysis of Action Potential Dispersion for Cardiac Arrhythmia Risk Stratification

Reentrant cardiac arrhythmias such as atrial fibrillation (AF) and ventricular fibrillation (VF) are common cardiac arrhythmias that account for substantial morbidity and mortality throughout the world. However, the mechanisms and optimal ablation treatment strategies for such arrhythmias are still unclear. Using 2D optical mapping of a mouse model with AF and VF, we have identified regional heterogeneity of the action potential duration (APD) in the atria and ventricles of the heart as key drivers for the initiation and persistence of reentry. The purpose of this paper is to discuss theoretical patterns of dispersion, demonstrate patterns of dispersion seen in our mouse model and discuss the computational analysis of APD dispersion patterns. These analyses and discussions may lead to better understanding of dispersion patterns in patients with these arrhythmias, as well as help comprehend whether and how reducing dispersion can lead to arrhythmia risk stratification and treatment strategies for arrhythmias.

Translational applications of computational modelling for patients with cardiac arrhythmias

Cardiac arrhythmia is associated with high morbidity, and its underlying mechanisms are poorly understood. Computational modelling and simulation approaches have the potential to improve standard-of-care therapy for these disorders, offering deeper understanding of complex disease processes and sophisticated translational tools for planning clinical procedures. This review provides a clinician-friendly summary of recent advancements in computational cardiology. Organ-scale models automatically generated from clinical-grade imaging data are used to custom tailor our understanding of arrhythmia drivers, estimate future arrhythmogenic risk and personalise treatment plans. Recent mechanistic insights derived from atrial and ventricular arrhythmia simulations are highlighted, and the potential avenues to patient care (eg, by revealing new antiarrhythmic drug targets) are covered. Computational approaches geared towards improving outcomes in resynchronisation therapy have used simulations to elucidate optimal patient selection and lead location. Technology to personalise catheter ablation procedures are also covered, specifically preliminary outcomes form early-stage or pilot clinical studies. To conclude, future developments in computational cardiology are discussed, including improving the representation of patient-specific fibre orientations and fibrotic remodelling characterisation and how these might improve understanding of arrhythmia mechanisms and provide transformative tools for patient-specific therapy.

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.

Using cardiac ionic cell models to interpret clinical data

For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

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.

Preprocedure Application of Machine Learning and Mechanistic Simulations Predicts Likelihood of Paroxysmal Atrial Fibrillation Recurrence Following Pulmonary Vein Isolation

BACKGROUND

Pulmonary vein isolation (PVI) is an effective treatment strategy for patients with atrial fibrillation (AF), but many experience AF recurrence and require repeat ablation procedures. The goal of this study was to develop and evaluate a methodology that combines machine learning (ML) and personalized computational modeling to predict, before PVI, which patients are most likely to experience AF recurrence after PVI.

METHODS

This single-center retrospective proof-of-concept study included 32 patients with documented paroxysmal AF who underwent PVI and had preprocedural late gadolinium enhanced magnetic resonance imaging. For each patient, a personalized computational model of the left atrium simulated AF induction via rapid pacing. Features were derived from pre-PVI late gadolinium enhanced magnetic resonance images and from results of simulations of AF induction. The most predictive features were used as input to a quadratic discriminant analysis ML classifier, which was trained, optimized, and evaluated with 10-fold nested cross-validation to predict the probability of AF recurrence post-PVI.

RESULTS

In our cohort, the ML classifier predicted probability of AF recurrence with an average validation sensitivity and specificity of 82% and 89%, respectively, and a validation area under the curve of 0.82. Dissecting the relative contributions of simulations of AF induction and raw images to the predictive capability of the ML classifier, we found that when only features from simulations of AF induction were used to train the ML classifier, its performance remained similar (validation area under the curve, 0.81). However, when only features extracted from raw images were used for training, the validation area under the curve significantly decreased (0.47).

CONCLUSIONS

ML and personalized computational modeling can be used together to accurately predict, using only pre-PVI late gadolinium enhanced magnetic resonance imaging scans as input, whether a patient is likely to experience AF recurrence following PVI, even when the patient cohort is small.

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'.

Constructing a Human Atrial Fibre Atlas

Atrial anisotropy affects electrical propagation patterns, anchor locations of atrial reentrant drivers, and atrial mechanics. However, patient-specific atrial fibre fields and anisotropy measurements are not currently available, and consequently assigning fibre fields to atrial models is challenging. We aimed to construct an atrial fibre atlas from a high-resolution DTMRI dataset that optimally reproduces electrophysiology simulation predictions corresponding to patient-specific fibre fields, and to develop a methodology for automatically assigning fibres to patient-specific anatomies. We extended an atrial coordinate system to map the pulmonary veins, vena cava and appendages to standardised positions in the coordinate system corresponding to the average location across the anatomies. We then expressed each fibre field in this atrial coordinate system and calculated an average fibre field. To assess the effects of fibre field on patient-specific modelling predictions, we calculated paced activation time maps and electrical driver locations during AF. In total, 756 activation time maps were calculated (7 anatomies with 9 fibre maps and 2 pacing locations, for the endocardial, epicardial and bilayer surface models of the LA and RA). Patient-specific fibre fields had a relatively small effect on average paced activation maps (range of mean local activation time difference for LA fields: 2.67-3.60 ms, and for RA fields: 2.29-3.44 ms), but had a larger effect on maximum LAT differences (range for LA 12.7-16.6%; range for RA 11.9-15.0%). A total of 126 phase singularity density maps were calculated (7 anatomies with 9 fibre maps for the LA and RA bilayer models). The fibre field corresponding to anatomy 1 had the highest median PS density map correlation coefficient for LA bilayer simulations (0.44 compared to the other correlations, ranging from 0.14 to 0.39), while the average fibre field had the highest correlation for the RA bilayer simulations (0.61 compared to the other correlations, ranging from 0.37 to 0.56). For sinus rhythm simulations, average activation time is robust to fibre field direction; however, maximum differences can still be significant. Patient specific fibres are more important for arrhythmia simulations, particularly in the left atrium. We propose using the fibre field corresponding to DTMRI dataset 1 for LA simulations, and the average fibre field for RA simulations as these optimally predicted arrhythmia properties.

Time resolution for wavefront and phase singularity tracking using activation maps in cardiac propagation models

The dynamics of cardiac fibrillation can be described by the number, the trajectory, the stability, and the lifespan of phase singularities (PSs). Accurate PS tracking is straightforward in simple uniform tissues but becomes more challenging as fibrosis, structural heterogeneity, and strong anisotropy are combined. In this paper, we derive a mathematical formulation for PS tracking in two-dimensional reaction-diffusion models. The method simultaneously tracks wavefronts and PS based on activation maps at full spatiotemporal resolution. PS tracking is formulated as a linear assignment problem solved by the Hungarian algorithm. The cost matrix incorporates information about distances between PS, chirality, and wavefronts. A graph of PS trajectories is generated to represent the creations and annihilations of PS pairs. Structure-preserving graph transformations are applied to provide a simplified description at longer observation time scales. The approach is validated in 180 simulations of fibrillation in four different types of substrates featuring, respectively, wavebreaks, ionic heterogeneities, fibrosis, and breakthrough patterns. The time step of PS tracking is studied in the range from 0.1 to 10 ms. The results show the benefits of improving time resolution from 1 to 0.1 ms. The tracking error rate decreases by an order of magnitude because the occurrence of simultaneous events becomes less likely. As observed on PS survival curves, the graph-based analysis facilitates the identification of macroscopically stable rotors despite wavefront fragmentation by fibrosis.

In situ procedure for high-efficiency computational modeling of atrial fibrillation reflecting personal anatomy, fiber orientation, fibrosis, and electrophysiology

We previously reported the feasibility and efficacy of a simulation-guided clinical catheter ablation of atrial fibrillation (AF) in an in-silico AF model. We developed a highly efficient realistic AF model reflecting the patient endocardial voltage and local conduction and tested its clinical feasibility. We acquired > 500 endocardial bipolar electrograms during right atrial pacing at the beginning of the AF ablation procedures. Based on the clinical bipolar electrograms, we generated simulated voltage maps by applying fibrosis and local activation maps adjusted for the fiber orientation. The software's accuracy (CUVIA2.5) was retrospectively tested in 17 patients and feasibility prospectively in 10 during clinical AF ablation. Results: We found excellent correlations between the clinical and simulated voltage maps (R = 0.933, p < 0.001) and clinical and virtual local conduction (R = 0.958, p < 0.001). The proportion of virtual local fibrosis was 15.4, 22.2, and 36.9% in the paroxysmal AF, persistent AF, and post-pulmonary vein isolation (PVI) states, respectively. The reconstructed virtual bipolar electrogram exhibited a relatively good similarities of morphology to the local clinical bipolar electrogram (R = 0.60 ± 0.08, p < 0.001). Feasibility testing revealed an in situ procedural computing time from the clinical data acquisition to wave-dynamics analyses of 48.2 ± 4.9 min. All virtual analyses were successfully achieved during clinical PVI procedures. We developed a highly efficient, realistic, in situ procedural simulation model reflective of individual anatomy, fiber orientation, fibrosis, and electrophysiology that can be applied during AF ablation.

Challenges Associated with Interpreting Mechanisms of AF

Determining optimal treatment strategies for complex arrhythmogenesis in AF is confounded by the lack of consensus regarding the mechanisms causing AF. Studies report different mechanisms for AF, ranging from hierarchical drivers to anarchical multiple activation wavelets. Differences in the assessment of AF mechanisms are likely due to AF being recorded across diverse models using different investigational tools, spatial scales and clinical populations. The authors review different AF mechanisms, including anatomical and functional re-entry, hierarchical drivers and anarchical multiple wavelets. They then describe different cardiac mapping techniques and analysis tools, including activation mapping, phase mapping and fibrosis identification. They explain and review different data challenges, including differences between recording devices in spatial and temporal resolutions, spatial coverage and recording surface, and report clinical outcomes using different data modalities. They suggest future research directions for investigating the mechanisms underlying human AF.

Addressing challenges of quantitative methodologies and event interpretation in the study of atrial fibrillation

Atrial fibrillation (AF) is the commonest arrhythmia, yet the mechanisms of its onset and persistence are incompletely known. Although techniques for quantitative assessment have been investigated, there have been few attempts to integrate this information to advance disease treatment protocols. In this review, key quantitative methods for AF analysis are described, and suggestions are provided for the coordination of the available information, and to develop foci and directions for future research efforts. Quantitative biologists may have an interest in this topic in order to develop machine learning and tools for arrhythmia characterization, but they may perhaps have a minimal background in the clinical methodology and in the types of observed events and mechanistic hypotheses that have thus far been developed. We attempt to address these issues via exploration of the published literature. Although no new data is presented in this review, examples are shown of current lines of investigation, and in particular, how electrogram analysis and whole-chamber quantitative modeling of the left atrium may be useful to characterize fibrillatory patterns of activity, so as to propose avenues for more efficacious acquisition and interpretation of AF data.

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.

Heterogeneity of the action potential duration is required for sustained atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia and accounts for substantial morbidity and mortality. Recently, we created a mouse model with spontaneous and sustained AF caused by a mutation in the NaV1.5 channel (F1759A) that enhances persistent Na+ current, thereby enabling the investigation of molecular mechanisms that cause AF and the identification of novel treatment strategies. The mice have regional heterogeneity of action potential duration of the atria similar to observations in patients with AF. In these mice, we found that the initiation and persistence of the rotational reentrant AF arrhythmias, known as spiral waves or rotors, were dependent upon action potential duration heterogeneity. The centers of the rotors were localized to regions of greatest heterogeneity of the action potential duration. Pharmacologically attenuating the action potential duration heterogeneity reduced both spontaneous and pacing-induced AF. Computer-based simulations also demonstrated that the action potential duration heterogeneity is sufficient to generate rotors that manifest as AF. Taken together, these findings suggest that action potential duration heterogeneity in mice and humans is one mechanism by which AF is initiated and that reducing action potential duration heterogeneity can lessen the burden 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.