A Computational Framework to Benchmark Basket Catheter Guided Ablation in Atrial Fibrillation

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant "gold standard". Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system's ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a "silver standard" comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

Computational Analysis of Mapping Catheter Geometry and Contact Quality Effects on Rotor Detection in Atrial Fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia and catheter mapping has been proved to be an effective approach for detecting AF drivers to be targeted by ablation. Among drivers, the so-called rotors have gained the most attention: their identification and spatial location could help to understand which patient-specific mechanisms are acting, and thus to guide the ablation execution. Since rotor detection by multi-electrode catheters may be influenced by several structural parameters including inter-electrode spacing, catheter coverage, and endocardium-catheter distance, in this study we proposed a tool for testing the ability of different catheter shapes to detect rotors in different conditions. An approach based on the solution of the monodomain equations coupled with a modified Courtemanche ionic atrial model, that considers an electrical remodeling, was applied to simulate spiral wave dynamics on a 2D model for 7.75 s. The developed framework allowed the acquisition of unipolar signals at 2 KHz. Two high-density multipolar catheters were simulated (Advisor™ HD Grid and PentaRay^®) and placed in a 2D region in which the simulated spiral wave persists longer. The configuration of the catheters was then modified by changing the number of electrodes, inter-electrodes distance, position, and atrial-wall distance for assessing how they would affect the rotor detection. In contact with the wall and at 1 mm distance from it, all the configurations detected the rotor correctly, irrespective of geometry, coverage, and inter-electrode distance. In the HDGrid-like geometry, the increase of the inter-electrode distance from 3 to 6 mm caused rotor detection failure at 2 mm distance from the LA wall. In the PentaRay-like configuration, regardless of inter-electrode distance, rotor detection failed at 3 mm endocardium-catheter distance. The asymmetry of this catheter resulted in rotation-dependent rotor detection. To conclude, the computational framework we developed is based on realistic catheter shapes designed with parameter configurations which resemble clinical settings. Results showed it is well suited to investigate how mapping catheter geometry and location affect AF driver detection, therefore it is a reliable tool to design and test new mapping catheters.

Applications of multimodality imaging for left atrial catheter ablation

Atrial arrhythmias, including atrial fibrillation and atrial flutter, may be treated through catheter ablation. The process of atrial arrhythmia catheter ablation, which includes patient selection, pre-procedural planning, intra-procedural guidance, and post-procedural assessment, is typically characterized by the use of several imaging modalities to sequentially inform key clinical decisions. Increasingly, advanced imaging modalities are processed via specialized image analysis techniques and combined with intra-procedural electrical measurements to inform treatment approaches. Here, we review the use of multimodality imaging for left atrial ablation procedures. The article first outlines how imaging modalities are routinely used in the peri-ablation period. We then describe how advanced imaging techniques may inform patient selection for ablation and ablation targets themselves. Ongoing research directions for improving catheter ablation outcomes by using imaging combined with advanced analyses for personalization of ablation targets are discussed, together with approaches for their integration in the standard clinical environment. Finally, we describe future research areas with the potential to improve catheter ablation outcomes.

Fibrotic Remodeling during Persistent Atrial Fibrillation: In Silico Investigation of the Role of Calcium for Human Atrial Myofibroblast Electrophysiology

During atrial fibrillation, cardiac tissue undergoes different remodeling processes at different scales from the molecular level to the tissue level. One central player that contributes to both electrical and structural remodeling is the myofibroblast. Based on recent experimental evidence on myofibroblasts' ability to contract, we extended a biophysical myofibroblast model with Ca²⁺ handling components and studied the effect on cellular and tissue electrophysiology. Using genetic algorithms, we fitted the myofibroblast model parameters to the existing in vitro data. In silico experiments showed that Ca²⁺ currents can explain the experimentally observed variability regarding the myofibroblast resting membrane potential. The presence of an L-type Ca²⁺ current can trigger automaticity in the myofibroblast with a cycle length of 799.9 ms. Myocyte action potentials were prolonged when coupled to myofibroblasts with Ca²⁺ handling machinery. Different spatial myofibroblast distribution patterns increased the vulnerable window to induce arrhythmia from 12 ms in non-fibrotic tissue to 22 ± 2.5 ms and altered the reentry dynamics. Our findings suggest that Ca²⁺ handling can considerably affect myofibroblast electrophysiology and alter the electrical propagation in atrial tissue composed of myocytes coupled with myofibroblasts. These findings can inform experimental validation experiments to further elucidate the role of myofibroblast Ca²⁺ handling in atrial arrhythmogenesis.

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.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.

Novel mapping techniques for rotor core detection using simulated intracardiac electrograms

BACKGROUND

Catheter ablation is associated with limited success rates in patients with persistent atrial fibrillation (AF). Currently, existing mapping systems fail to identify critical target sites for ablation. Recently, we proposed and validated several techniques (multiscale frequency [MSF], Shannon entropy [SE], kurtosis [Kt], and multiscale entropy [MSE]) to identify pivot point of rotors using ex-vivo optical mapping animal experiments. However, the performance of these techniques is unclear for the clinically recorded intracardiac electrograms (EGMs), due to the different nature of the signals.

OBJECTIVE

This study aims to evaluate the performance of MSF, MSE, SE, and Kt techniques to identify the pivot point of the rotor using unipolar and bipolar EGMs obtained from numerical simulations.

METHODS

Stationary and meandering rotors were simulated in a 2D human atria. The performances of new approaches were quantified by comparing the "true" core of the rotor with the core identified by the techniques. Also, the performances of all techniques were evaluated in the presence of noise, scar, and for the case of the multielectrode multispline and grid catheters.

RESULTS

Our results demonstrate that all the approaches are able to accurately identify the pivot point of both stationary and meandering rotors from both unipolar and bipolar EGMs. The presence of noise and scar tissue did not significantly affect the performance of the techniques. Finally, the core of the rotors was correctly identified for the case of multielectrode multispline and grid catheter simulations.

CONCLUSION

The core of rotors can be successfully identified from EGMs using novel techniques; thus, providing motivation for future clinical implementations.

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.

Universal atrial coordinates applied to visualisation, registration and construction of patient specific meshes

•

We introduce a coordinate system for the atria based on anatomical landmarks.

•

We construct the coordinates from solutions to Laplace’s equation.

•

We demonstrate the mapping of both scalar and vector data between different atria.

•

The coordinate system was used for registration and 2D visualisation of multimodal data.

•

Patient specific meshes with atrial structures and fibre direction were constructed using just five landmark points.

Integrating spatial information about atrial physiology and anatomy in a single patient from multimodal datasets, as well as generalizing these data across patients, requires a common coordinate system. In the atria, this is challenging due to the complexity and variability of the anatomy. We aimed to develop and validate a Universal Atrial Coordinate (UAC) system for the following applications: combination and assessment of multimodal data; comparison of spatial data across patients; 2D visualization; and construction of patient specific geometries to test mechanistic hypotheses. Left and right atrial LGE-MRI data were segmented and meshed. Two coordinates were calculated for each atrium by solving Laplace’s equation, with boundary conditions assigned using five landmark points. The coordinate system was used to map spatial information between atrial meshes, including scalar fields measured using different mapping modalities, and atrial anatomic structures and fibre directions from a reference geometry. Average error in point transfer from a source mesh to a destination mesh and back again was less than 0.1 mm for the left atrium and 0.02 mm for the right atrium. Patient specific meshes were constructed using the coordinate system and phase singularity density maps from arrhythmia simulations were visualised in 2D. In conclusion, we have developed a universal atrial coordinate system allowing automatic registration of imaging and electroanatomic mapping data, 2D visualisation, and patient specific model creation.

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.

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

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