An efficient finite element approach for modeling fibrotic clefts in the heart

Quantifying anatomically-based in-silico electrocardiogram variability for cardiac digital twins

Human cardiac Cardiac digital twins (CDTs) are digital replicas of patient hearts, designed to match clinical observations precisely. The electro-cardiogram (ECG), as the most common non-invasive electrophysiology (EP) measurement, has been recently successfully employed for calibrating CDT. However, ECG-based calibration methods often fail to account for the inherent uncertainties in clinical data acquisition and CDT anatomical generation workflows. As a result, discrepancies inevitably arise between the actual physical and simulated patient EP and ECG. In this study, we aim to qualitatively and quantitatively analyze the impact of these uncertainties on ECG morphology and diagnostic markers, and therefore to assess the reliability of ECG-based CDT calibration. We analyze residual beat-to-beat variability in ECG recordings obtained from three datasets, including healthy subjects and patients treated for ventricular tachycardia and atrial fibrillation. Using a biophysically detailed and anatomically accurate computational model of whole-heart EP combined with a detailed torso model calibrated to closely replicate measured ECG signals, we vary anatomical factors (heart location, orientation, size), heterogeneity in electrical conductivities in the heart and torso, and electrode placements across ECG leads to assess their qualitative impact on ECG morphology. Our study demonstrates that diagnostically relevant ECG features and overall morphology remain close to the ground through ECG independently of the investigated uncertainties. This resilience is consistent with the narrow distribution of ECG due to residual beat-to-beat variability observed in both healthy subjects and patients. Overall, our results suggest that observation uncertainties do not impede an accurate calibration of the CDT.

Fibroblast Density is a Risk Factor for Drug-induced Arrhythmias

A recent study by Kawatou et al . has shown that the local heterogeneity of ion channel conductance is a critical substrate for focal or reentrant arrhythmias. However, the role of fibroblasts with repolarization heterogeneity in the initiation and maintenance of arrhythmias remains unknown. In this study, we investigated how diffuse fibrosis contributes to the formation of focal and reentrant arrhythmias under drug-induced heterogeneity using physiologically detailed mathematical models of the human heart. To simulate drug-induced heterogeneity, we varied the maximum conductance of transmembrane potassium and calcium currents, leading to heterogeneity in action potential duration (APD). Then, we assessed the effects of different fibrosis densities (FD) on the occurrence of premature ventricular complexes (PVCs). Fibroblasts were randomly and evenly inserted into the tissue, and various FD levels ranging from 0 to 35% were examined. We found a biphasic relationship between FD and drug-induced PVCs. Within a certain range of FD, FD positively correlated with PVC susceptibility. However, excessively high fibrosis levels were associated with reduced susceptibility to PVCs. In addition, the self-sustainability of arrhythmias exhibited a positive correlation with FD. This study demonstrates the interplay between the diffuse fibrosis and the drug-induced heterogeneity of APD in the genesis of ventricular arrhythmias.

Author summary

Sudden cardiac death remains a leading cause of death worldwide. Understanding the mechanisms underlying arrhythmia and its precursors is critical for the development of effective therapies and drugs. Repolarization heterogeneity plays a crucial role in both the initiation and maintenance of arrhythmias. Fibroblasts constitute a vital component of cardiac structure, originating from the remodeling of ventricular wall cells or the transformation of injured myocardial cells. Fibroblasts are known to couple with and alter the electrical properties of myocardial cells. However, our understanding of the role of fibroblasts in the development of arrhythmia remains limited. In this study, we employed a physiologically detailed mathematical model of cardiac tissue to investigate the roles of drug-induced heterogeneity and diffuse fibrosis in the initiation and maintenance of arrhythmias. We used 2D and 3D computational models to simulate various levels of drug-induced heterogeneity conditions with normal to pathological levels of fibroblast density (FD). We found that within a certain range of FD, fibroblasts promote PVCs under drug-induced heterogeneity. However, if FD exceeds 30%, the occurrence of PVCs decreases (biphasic relationship). On the other hand, the self-sustainability of VF (ventricular fibrillation) consistently increases with FD. This study implies that fibroblasts in cardiac tissue may play different roles in the initiation and maintenance of arrhythmia.

Re-entry in models of cardiac ventricular tissue with scar represented as a Gaussian random field

Introduction: Fibrotic scar in the heart is known to act as a substrate for arrhythmias. Regions of fibrotic scar are associated with slowed or blocked conduction of the action potential, but the detailed mechanisms of arrhythmia formation are not well characterised and this can limit the effective diagnosis and treatment of scar in patients. The aim of this computational study was to evaluate different representations of fibrotic scar in models of 2D 10 × 10 cm ventricular tissue, where the region of scar was defined by sampling a Gaussian random field with an adjustable length scale of between 1.25 and 10.0 mm. Methods: Cellular electrophysiology was represented by the Ten Tusscher 2006 model for human ventricular cells. Fibrotic scar was represented as a spatially varying diffusion, with different models of the boundary between normal and fibrotic tissue. Dispersion of activation time and action potential duration (APD) dispersion was assessed in each sample by pacing at an S1 cycle length of 400 ms followed by a premature S2 beat with a coupling interval of 323 ms. Vulnerability to reentry was assessed with an aggressive pacing protocol. In all models, simulated fibrosis acted to delay activation, to increase the dispersion of APD, and to generate re-entry. Results: A higher incidence of re-entry was observed in models with simulated fibrotic scar at shorter length scale, but the type of model used to represent fibrotic scar had a much bigger influence on the incidence of reentry. Discussion: This study shows that in computational models of fibrotic scar the effects that lead to either block or propagation of the action potential are strongly influenced by the way that fibrotic scar is represented in the model, and so the results of computational studies involving fibrotic scar should be interpreted carefully.

Fibrosis modeling choice affects morphology of ventricular arrhythmia in non-ischemic cardiomyopathy

Introduction: Patients with non-ischemic cardiomyopathy (NICM) are at risk for ventricular arrhythmias, but diagnosis and treatment planning remain a serious clinical challenge. Although computational modeling has provided valuable insight into arrhythmic mechanisms, the optimal method for simulating reentry in NICM patients with structural disease is unknown. Methods: Here, we compare the effects of fibrotic representation on both reentry initiation and reentry morphology in patient-specific cardiac models. We investigate models with heterogeneous networks of non-conducting structures (cleft models) and models where fibrosis is represented as a dense core with a surrounding border zone (non-cleft models). Using segmented cardiac magnetic resonance with late gadolinium enhancement (LGE) of five NICM patients, we created 185 3D ventricular electrophysiological models with different fibrotic representations (clefts, reduced conductivity and ionic remodeling). Results: Reentry was induced by electrical pacing in 647 out of 3,145 simulations. Both cleft and non-cleft models can give rise to double-loop reentries meandering through fibrotic regions (Type 1-reentry). When accounting for fibrotic volume, the initiation sites of these reentries are associated with high local fibrotic density (mean LGE in cleft models: p< 0.001, core volume in non-cleft models: p = 0.018, negative binomial regression). In non-cleft models, Type 1-reentries required slow conduction in core tissue (non-cleftsc models) as opposed to total conduction block. Incorporating ionic remodeling in fibrotic regions can give rise to single- or double-loop rotors close to healthy-fibrotic interfaces (Type 2-reentry). Increasing the cleft density or core-to-border zone ratio in cleft and non-cleftc models, respectively, leads to increased inducibility and a change in reentry morphology from Type 2 to Type 1. Conclusions: By demonstrating how fibrotic representation affects reentry morphology and location, our findings can aid model selection for simulating arrhythmogenesis in NICM.

Homogenisation for the monodomain model in the presence of microscopic fibrotic structures

Computational models in cardiac electrophysiology are notorious for long runtimes, restricting the numbers of nodes and mesh elements in the numerical discretisations used for their solution. This makes it particularly challenging to incorporate structural heterogeneities on small spatial scales, preventing a full understanding of the critical arrhythmogenic effects of conditions such as cardiac fibrosis. In this work, we explore the technique of homogenisation by volume averaging for the inclusion of non-conductive micro-structures into larger-scale cardiac meshes with minor computational overhead. Importantly, our approach is not restricted to periodic patterns, enabling homogenised models to represent, for example, the intricate patterns of collagen deposition present in different types of fibrosis. We first highlight the importance of appropriate boundary condition choice for the closure problems that define the parameters of homogenised models. Then, we demonstrate the technique's ability to correctly upscale the effects of fibrotic patterns with a spatial resolution of 10 µm into much larger numerical mesh sizes of 100- 250 µm . The homogenised models using these coarser meshes correctly predict critical pro-arrhythmic effects of fibrosis, including slowed conduction, source/sink mismatch, and stabilisation of re-entrant activation patterns. As such, this approach to homogenisation represents a significant step towards whole organ simulations that unravel the effects of microscopic cardiac tissue heterogeneities.

A personalized real-time virtual model of whole heart electrophysiology

Computer models capable of representing the intrinsic personal electrophysiology (EP) of the heart in silico are termed virtual heart technologies. When anatomy and EP are tailored to individual patients within the model, such technologies are promising clinical and industrial tools. Regardless of their vast potential, few virtual technologies simulating the entire organ-scale EP of all four-chambers of the heart have been reported and widespread clinical use is limited due to high computational costs and difficulty in validation. We thus report on the development of a novel virtual technology representing the electrophysiology of all four-chambers of the heart aiming to overcome these limitations. In our previous work, a model of ventricular EP embedded in a torso was constructed from clinical magnetic resonance image (MRI) data and personalized according to the measured 12 lead electrocardiogram (ECG) of a single subject under normal sinus rhythm. This model is then expanded upon to include whole heart EP and a detailed representation of the His-Purkinje system (HPS). To test the capacities of the personalized virtual heart technology to replicate standard clinical morphological ECG features under such conditions, bundle branch blocks within both the right and the left ventricles under two different conduction velocity settings are modeled alongside sinus rhythm. To ensure clinical viability, model generation was completely automated and simulations were performed using an efficient real-time cardiac EP simulator. Close correspondence between the measured and simulated 12 lead ECG was observed under normal sinus conditions and all simulated bundle branch blocks manifested relevant clinical morphological features.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

An automated near-real time computational method for induction and treatment of scar-related ventricular tachycardias

Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.

A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms

Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Etiology-Specific Remodeling in Ventricular Tissue of Heart Failure Patients and Its Implications for Computational Modeling of Electrical Conduction

With an estimated 64.3 million cases worldwide, heart failure (HF) imposes an enormous burden on healthcare systems. Sudden death from arrhythmia is the major cause of mortality in HF patients. Computational modeling of the failing heart provides insights into mechanisms of arrhythmogenesis, risk stratification of patients, and clinical treatment. However, the lack of a clinically informed approach to model cardiac tissues in HF hinders progress in developing patient-specific strategies. Here, we provide a microscopy-based foundation for modeling conduction in HF tissues. We acquired 2D images of left ventricular tissues from HF patients (n = 16) and donors (n = 5). The composition and heterogeneity of fibrosis were quantified at a sub-micrometer resolution over an area of 1 mm². From the images, we constructed computational bidomain models of tissue electrophysiology. We computed local upstroke velocities of the membrane voltage and anisotropic conduction velocities (CV). The non-myocyte volume fraction was higher in HF than donors (39.68 ± 14.23 vs. 22.09 ± 2.72%, p < 0.01), and higher in ischemic (IC) than nonischemic (NIC) cardiomyopathy (47.2 ± 16.18 vs. 32.16 ± 6.55%, p < 0.05). The heterogeneity of fibrosis within each subject was highest for IC (27.1 ± 6.03%) and lowest for donors (7.47 ± 1.37%) with NIC (15.69 ± 5.76%) in between. K-means clustering of this heterogeneity discriminated IC and NIC with an accuracy of 81.25%. The heterogeneity in CV increased from donor to NIC to IC tissues. CV decreased with increasing fibrosis for longitudinal (R ² = 0.28, p < 0.05) and transverse conduction (R ² = 0.46, p < 0.01). The tilt angle of the CV vectors increased 2.1° for longitudinal and 0.91° for transverse conduction per 1% increase in fibrosis. Our study suggests that conduction fundamentally differs in the two etiologies due to the characteristics of fibrosis. Our study highlights the importance of the etiology-specific modeling of HF tissues and integration of medical history into electrophysiology models for personalized risk stratification and treatment planning.

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.

3D Electrophysiological Modeling of Interstitial Fibrosis Networks and Their Role in Ventricular Arrhythmias in Non-Ischemic Cardiomyopathy

OBJECTIVE

Interstitial fibrosis is a pathological expansion of the heart's inter-cellular collagen matrix. It is a potential complication of nonischemic cardiomyopathy (NICM), a class of diseases involving electrical and or mechanical dysfunction of cardiac tissue not caused by atherosclerosis. Patients with NICM and interstitial fibrosis often suffer from life threatening arrhythmias, which we aim to simulate in this study.

METHODS

Our methodology builds on an efficient discrete finite element (DFE) method which allows for the representation of fibrosis as infinitesimal splits in a mesh. We update the DFE method with a local connectivity analysis which creates a consistent topology in the fibrosis network. This is particularly important in nonischemic disease due to the potential presence of large and contiguous fibrotic regions and therefore potentially complex fibrosis networks.

RESULTS

In experiments with an image-based model, we demonstrate that our methodology is able to simulate reentrant electrical events associated with cardiac arrhythmias. These reentries depended crucially upon sufficient fibrosis density, which was marked by conduction slowing at high pacing rates. We also created a 2D test-case which demonstrated that fibrosis topologies can modulate transient conduction block, and thereby reentrant activations.

CONCLUSION

Ventricular arrhythmias due to interstitial fibrosis in NICM can be efficiently simulated using our methods in medical image based geometries. Furthermore, fibrosis topology modulates transient conduction block, and should be accounted for in electrophysiological simulations with interstitial fibrosis.

SIGNIFICANCE

Our study provides methodology which has the potential to predict arrhythmias and to optimize treatments non-invasively for nonischemic cardiomyopathies.

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.

Myocardial Fibrosis in a 3D Model: Effect of Texture on Wave Propagation

Non-linear electrical waves propagate through the heart and control cardiac contraction. Abnormal wave propagation causes various forms of the heart disease and can be lethal. One of the main causes of abnormality is a condition of cardiac fibrosis, which, from mathematical point of view, is the presence of multiple non-conducting obstacles for wave propagation. The fibrosis can have different texture which varies from diffuse (e.g., small randomly distributed obstacles), patchy (e.g., elongated interstitional stria), and focal (e.g., post-infarct scars) forms. Recently, Nezlobinsky et al. (2020) used 2D biophysical models to quantify the effects of elongation of obstacles (fibrosis texture) and showed that longitudinal and transversal propagation differently depends on the obstacle length resulting in anisotropy for wave propagation. In this paper, we extend these studies to 3D tissue models. We show that 3D consideration brings essential new effects; for the same obstacle length in 3D systems, anisotropy is about two times smaller compared to 2D, however, wave propagation is more stable with percolation threshold of about 60% (compared to 35% in 2D). The percolation threshold increases with the obstacle length for the longitudinal propagation, while it decreases for the transversal propagation. Further, in 3D, the dependency of velocity on the obstacle length for the transversal propagation disappears.

Anisotropic conduction in the myocardium due to fibrosis: the effect of texture on wave propagation

Cardiac fibrosis occurs in many forms of heart disease. It is well established that the spatial pattern of fibrosis, its texture, substantially affects the onset of arrhythmia. However, in most modelling studies fibrosis is represented by multiple randomly distributed short obstacles that mimic only one possible texture, diffuse fibrosis. An important characteristic feature of other fibrosis textures, such as interstitial and patchy textures, is that fibrotic inclusions have substantial length, which is suggested to have a pronounced effect on wave propagation. In this paper, we study the effect of the elongation of inexcitable inclusions (obstacles) on wave propagation in a 2D model of cardiac tissue described by the TP06 model for human ventricular cells. We study in detail how the elongation of obstacles affects various characteristics of the waves. We quantify the anisotropy induced by the textures, its dependency on the obstacle length and the effects of the texture on the shape of the propagating wave. Because such anisotropy is a result of zig-zag propagation we show, for the first time, quantification of the effects of geometry and source-sink relationship, on the zig-zag nature of the pathway of electrical conduction. We also study the effect of fibrosis in the case of pre-existing anisotropy and introduce a procedure for scaling of the fibrosis texture. We show that fibrosis can decrease or increase the preexisting anisotropy depending on its scaled texture.

Scar shape analysis and simulated electrical instabilities in a non-ischemic dilated cardiomyopathy patient cohort

This paper presents a morphological analysis of fibrotic scarring in non-ischemic dilated cardiomyopathy, and its relationship to electrical instabilities which underlie reentrant arrhythmias. Two dimensional electrophysiological simulation models were constructed from a set of 699 late gadolinium enhanced cardiac magnetic resonance images originating from 157 patients. Areas of late gadolinium enhancement (LGE) in each image were assigned one of 10 possible microstructures, which modelled the details of fibrotic scarring an order of magnitude below the MRI scan resolution. A simulated programmed electrical stimulation protocol tested each model for the possibility of generating either a transmural block or a transmural reentry. The outcomes of the simulations were compared against morphological LGE features extracted from the images. Models which blocked or reentered, grouped by microstructure, were significantly different from one another in myocardial-LGE interface length, number of components and entropy, but not in relative area and transmurality. With an unknown microstructure, transmurality alone was the best predictor of block, whereas a combination of interface length, transmurality and number of components was the best predictor of reentry in linear discriminant analysis.

Acetylcholine Delays Atrial Activation to Facilitate Atrial Fibrillation

Background

Acetylcholine (ACh) shortens action potential duration (APD) in human atria. APD shortening facilitates atrial fibrillation (AF) by reducing the wavelength for reentry. However, the influence of ACh on electrical conduction in human atria and its contribution to AF are unclear, particularly when combined with impaired conduction from interstitial fibrosis.

Objective

To investigate the effect of ACh on human atrial conduction and its role in AF with computational, experimental, and clinical approaches.

Methods

S1S2 pacing (S1 = 600 ms and S2 = variable cycle lengths) was applied to the following human AF computer models: a left atrial appendage (LAA) myocyte to quantify the effects of ACh on APD, maximum upstroke velocity (V _max ), and resting membrane potential (RMP); a monolayer of LAA myocytes to quantify the effects of ACh on conduction; and 3) an intact left atrium (LA) to determine the effects of ACh on arrhythmogenicity. Heterogeneous ACh and interstitial fibrosis were applied to the monolayer and LA models. To corroborate the simulations, APD and RMP from isolated human atrial myocytes were recorded before and after 0.1 μM ACh. At the tissue level, LAAs from AF patients were optically mapped ex vivo using Di-4-ANEPPS. The difference in total activation time (AT) was determined between AT initially recorded with S1 pacing, and AT recorded during subsequent S1 pacing without (n = 6) or with (n = 7) 100 μM ACh.

Results

In LAA myocyte simulations, S1 pacing with 0.1 μM ACh shortened APD by 41 ms, hyperpolarized RMP by 7 mV, and increased V _max by 27 mV/ms. In human atrial myocytes, 0.1 μM ACh shortened APD by 48 ms, hyperpolarized RMP by 3 mV, and increased V _max by 6 mV/ms. In LAA monolayer simulations, S1 pacing with ACh hyperpolarized RMP to delay total AT by 32 ms without and 35 ms with fibrosis. This led to unidirectional conduction block and sustained reentry in fibrotic LA with heterogeneous ACh during S2 pacing. In AF patient LAAs, S1 pacing with ACh increased total AT from 39.3 ± 26 ms to 71.4 ± 31.2 ms (p = 0.036) compared to no change without ACh (56.7 ± 29.3 ms to 50.0 ± 21.9 ms, p = 0.140).

Conclusion

In fibrotic atria with heterogeneous parasympathetic activation, ACh facilitates AF by shortening APD and slowing conduction to promote unidirectional conduction block and reentry.

Factors Promoting Conduction Slowing as Substrates for Block and Reentry in Infarcted Hearts

The development of effective and safe therapies for scar-related ventricular tachycardias requires a detailed understanding of the mechanisms underlying the conduction block that initiates electrical re-entries associated with these arrhythmias. Conduction block has been often associated with electrophysiological changes that prolong action potential duration (APD) within the border zone (BZ) of chronically infarcted hearts. However, experimental evidence suggests that remodeling processes promoting conduction slowing as opposed to APD prolongation mark the chronic phase. In this context, the substrate for the initial block at the mouth of an isthmus/diastolic channel leading to ventricular tachycardia is unclear. The goal of this study was to determine whether electrophysiological parameters associated with conduction slowing can cause block and re-entry in the BZ. In silico experiments were conducted on two-dimensional idealized infarct tissue as well as on a cohort of postinfarction porcine left ventricular models constructed from ex vivo magnetic resonance imaging scans. Functional conduction slowing in the BZ was modeled by reducing sodium current density, whereas structural conduction slowing was represented by decreasing tissue conductivity and including fibrosis. The arrhythmogenic potential of APD prolongation was also tested as a basis for comparison. Within all models, the combination of reduced sodium current with structural remodeling more often degenerated into re-entry and, if so, was more likely to be sustained for more cycles. Although re-entries were also detected in experiments with prolonged APD, they were often not sustained because of the subsequent block caused by long-lasting repolarization. Functional and structural conditions associated with slow conduction rather than APD prolongation form a potent substrate for arrhythmogenesis at the isthmus/BZ of chronically infarcted hearts. Reduced excitability led to block while slow conduction shortened the wavelength of propagation, facilitating the sustenance of re-entries. These findings provide important insights for models of patient-specific risk stratification and therapy planning.

Extended Bidomain Modeling of Defibrillation: Quantifying Virtual Electrode Strengths in Fibrotic Myocardium

Defibrillation is a well-established therapy for atrial and ventricular arrhythmia. Here, we shed light on defibrillation in the fibrotic heart. Using the extended bidomain model of electrical conduction in cardiac tissue, we assessed the influence of fibrosis on the strength of virtual electrodes caused by extracellular electrical current. We created one-dimensional models of rabbit ventricular tissue with a central patch of fibrosis. The fibrosis was incorporated by altering volume fractions for extracellular, myocyte and fibroblast domains. In our prior work, we calculated these volume fractions from microscopic images at the infarct border zone of rabbit hearts. An average and a large degree of fibrosis were modeled. We simulated defibrillation by application of an extracellular current for a short duration (5 ms). We explored the effects of myocyte-fibroblast coupling, intra-fibroblast conductivity and patch length on the strength of the virtual electrodes present at the borders of the normal and fibrotic tissue. We discriminated between effects on myocyte and fibroblast membranes at both borders of the patch. Similarly, we studied defibrillation in two-dimensional models of fibrotic tissue. Square and disk-like patches of fibrotic tissue were embedded in control tissue. We quantified the influence of the geometry and fibrosis composition on virtual electrode strength. We compared the results obtained with a square and disk shape of the fibrotic patch with results from the one-dimensional simulations. Both, one- and two-dimensional simulations indicate that extracellular current application causes virtual electrodes at boundaries of fibrotic patches. A higher degree of fibrosis and larger patch size were associated with an increased strength of the virtual electrodes. Also, patch geometry affected the strength of the virtual electrodes. Our simulations suggest that increased fibroblast-myocyte coupling and intra-fibroblast conductivity reduce virtual electrode strength. However, experimental data to constrain these modeling parameters are limited and thus pinpointing the magnitude of the reduction will require further understanding of electrical coupling of fibroblasts in native cardiac tissues. We propose that the findings from our computational studies are important for development of patient-specific protocols for internal defibrillators.

Fibrosis Microstructure Modulates Reentry in Non-ischemic Dilated Cardiomyopathy: Insights From Imaged Guided 2D Computational Modeling

Aims: Patients who present with non-ischemic dilated cardiomyopathy (NIDCM) and enhancement on late gadolinium magnetic resonance imaging (LGE-CMR), are at high risk of sudden cardiac death (SCD). Further risk stratification of these patients based on LGE-CMR may be improved through better understanding of fibrosis microstructure. Our aim is to examine variations in fibrosis microstructure based on LGE imaging, and quantify the effect on reentry inducibility and mechanism. Furthermore, we examine the relationship between transmural activation time differences and reentry. Methods and Results: 2D Computational models were created from a single short axis LGE-CMR image, with 401 variations in fibrosis type (interstitial, replacement) and density, as well as presence or absence of reduced conductivity (RC). Transmural activation times (TAT) were measured, as well as reentry incidence and mechanism. Reentries were inducible above specific density thresholds (0.8, 0.6 for interstitial, replacement fibrosis). RC reduced these thresholds (0.3, 0.4 for interstitial, replacement fibrosis) and increased reentry incidence (48 no RC vs. 133 with RC). Reentries were classified as rotor, micro-reentry, or macro-reentry and depended on fibrosis micro-structure. Differences in TAT at coupling intervals 210 and 500ms predicted reentry in the models (sensitivity 89%, specificity 93%). A sensitivity analysis of TAT and reentry incidence showed that these quantities were robust to small changes in the pacing location. Conclusion: Computational models of fibrosis micro-structure underlying areas of LGE in NIDCM provide insight into the mechanisms and inducibility of reentry, and their dependence upon the type and density of fibrosis. Transmural activation times, measured at the central extent of the scar, can potentially differentiate microstructures which support reentry.

Quantitative descriptions of extended volume cardiac tissue architecture from multiple large hearts

The arrangement of cardiac cells into strand and sheet-like structures within the heart wall, confers important electrical properties onto heart tissue. Unraveling cardiomyocyte architecture in both healthy and diseased hearts is fundamental to understanding the mechanisms generating normal rhythm and arrhythmia. We analyzed five extended volume serial image stacks of normal pig left ventricular tissue. Analysis included: (1) reconstruction of original tissue volume and shape with non-linear correction maps; (2) segmentation and higher-order descriptions, areas and orientations of laminar structures through the heart wall; (3)computation of fiber directions; (4) computation of tissue connectivity using a shell filter. These measures contributed to a deeper and more objective understanding of cardiac tissue structures and their spatial variation than previously possible.

Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics

Success rates for catheter ablation of persistent atrial fibrillation patients are currently low; however, there is a subset of patients for whom electrical isolation of the pulmonary veins alone is a successful treatment strategy. It is difficult to identify these patients because there are a multitude of factors affecting arrhythmia susceptibility and maintenance, and the individual contributions of these factors are difficult to determine clinically. We hypothesised that the combination of pulmonary vein (PV) electrophysiology and atrial body fibrosis determine driver location and effectiveness of pulmonary vein isolation (PVI). We used bilayer biatrial computer models based on patient geometries to investigate the effects of PV properties and atrial fibrosis on arrhythmia inducibility, maintenance mechanisms, and the outcome of PVI. Short PV action potential duration (APD) increased arrhythmia susceptibility, while longer PV APD was found to be protective. Arrhythmia inducibility increased with slower conduction velocity (CV) at the LA/PV junction, but not for cases with homogeneous CV changes or slower CV at the distal PV. Phase singularity (PS) density in the PV region for cases with PV fibrosis was increased. Arrhythmia dynamics depend on both PV properties and fibrosis distribution, varying from meandering rotors to PV reentry (in cases with baseline or long APD), to stable rotors at regions of high fibrosis density. Measurement of fibrosis and PV properties may indicate patient specific susceptibility to AF initiation and maintenance. PV PS density before PVI was higher for cases in which AF terminated or converted to a macroreentry; thus, high PV PS density may indicate likelihood of PVI success.

Atrial fibrillation is the most commonly encountered cardiac arrhythmia, affecting a significant portion of the population. Currently, ablation is the most effective treatment but success rates are less than optimal, being 70% one-year post-treatment. There is a large effort to find better ablation strategies to permanently cure the condition. Pulmonary vein isolation by ablation is more or less the standard of care, but many questions remain since pulmonary vein ectopy by itself does not explain all of the clinical successes or failures. We used computer simulations to investigate how electrophysiological properties of the pulmonary veins can affect rotor formation and maintenance in patients suffering from atrial fibrillation. We used complex, biophysical representations of cellular electrophysiology in highly detailed geometries constructed from patient scans. We heterogeneously varied electrophysiological and structural properties to see their effects on rotor initiation and maintenance. Our study suggests a metric for indicating the likelihood of success of pulmonary vein isolation. Thus either measuring this clinically, or running patient-specific simulations to estimate this metric may suggest whether ablation in addition to pulmonary vein isolation should be performed. Our study provides motivation for a retrospective clinical study or experimental study into this metric.

Modeling dynamics in diseased cardiac tissue: Impact of model choice

Cardiac arrhythmias have been traditionally simulated using continuous models that assume tissue homogeneity and use a relatively large spatial discretization. However, it is believed that the tissue fibrosis and collagen deposition, which occur on a micron-level, are critical factors in arrhythmogenesis in diseased tissues. Consequently, it remains unclear how well continuous models, which use averaged electrical properties, are able to accurately capture complex conduction behaviors such as re-entry in fibrotic tissues. The objective of this study was to compare re-entrant behavior in discrete microstructural models of fibrosis and in two types of equivalent continuous models, a homogenous continuous model and a hybrid continuous model with distinct heterogeneities. In the discrete model, increasing levels of tissue fibrosis lead to a substantial increase in the re-entrant cycle length which is inadequately reflected in the homogenous continuous models. These cycle length increases appear to be primarily due to increases in the tip path length and to altered restitution behavior, and suggest that it is critical to consider the discrete effects of fibrosis on conduction when studying arrhythmogenesis in fibrotic myocardium. Hybrid models are able to accurately capture some aspects of re-entry and, if carefully tuned, may provide a framework for simulating conduction in diseased tissues with both accuracy and efficiency.

Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms

AIMS

Catheter ablation is an effective technique for terminating atrial arrhythmia. However, given a high atrial fibrillation (AF) recurrence rate, optimal ablation strategies have yet to be defined. Computer modelling can be a powerful aid but modelling of fibrosis, a major factor associated with AF, is an open question. Several groups have proposed methodologies based on imaging data, but no comparison to determine which methodology best corroborates clinically observed reentrant behaviour has been performed. We examined several methodologies to determine the best method for capturing fibrillation dynamics.

METHODS AND RESULTS

Patient late gadolinium-enhanced magnetic resonance imaging data were transferred onto a bilayer atrial computer model and used to assign fibrosis distributions. Fibrosis was modelled as conduction disturbances (lower conductivity, edge splitting, or percolation), transforming growth factor-β1 ionic channel effects, myocyte-fibroblast coupling, and combinations of the preceding. Reentry was induced through pulmonary vein ectopy and the ensuing rotor dynamics characterized. Non-invasive electrocardiographic imaging data of the patients in AF was used for comparison. Electrograms were computed and the fractionation durations measured over the surface. Edge splitting produced more phase singularities from wavebreaks than the other representations. The number of phase singularities seen with percolation was closer to the clinical values. Addition of fibroblast coupling had an organizing effect on rotor dynamics. Simple tissue conductivity changes with ionic changes localized rotors over fibrosis which was not observed with clinical data.

CONCLUSION

The specific representation of fibrosis has a large effect on rotor dynamics and needs to be carefully considered for patient specific modelling.

Bidomain Predictions of Virtual Electrode-Induced Make and Break Excitations around Blood Vessels

Introduction and background

Virtual electrodes formed by field stimulation during defibrillation of cardiac tissue play an important role in eliciting activations. It has been suggested that the coronary vasculature is an important source of virtual electrodes, especially during low-energy defibrillation. This work aims to further the understanding of how virtual electrodes from the coronary vasculature influence defibrillation outcomes.

Methods

Using the bidomain model, we investigated how field stimulation elicited activations from virtual electrodes around idealized intramural blood vessels. Strength–interval curves, which quantify the stimulus strength required to elicit wavefront propagation from the vessels at different states of tissue refractoriness, were computed for each idealized geometry.

Results

Make excitations occurred at late diastolic intervals, originating from regions of depolarization around the vessel. Break excitations occurred at early diastolic intervals, whereby the vessels were able to excite surrounding refractory tissue due to the local restoration of excitability by virtual electrode-induced hyperpolarizations. Overall, strength–interval curves had similar morphologies and underlying excitation mechanisms compared with previous experimental and numerical unipolar stimulation studies of cardiac tissue. Including the presence of the vessel wall increased the field strength required for make excitations but decreased the field strength required for break excitations, and the field strength at which break excitations occurred was generally greater than 5 V/cm. Finally, in a more realistic ventricular slice geometry, the proximity of virtual electrodes around subepicardial vessels was seen to cause break excitations in the form of propagating unstable wavelets to the subepicardial layer.

Conclusion

Representing the blood vessel wall microstructure in computational bidomain models of defibrillation is recommended as it significantly alters the electrophysiological response of the vessel to field stimulation. Although vessels may facilitate excitation of relatively refractory tissue via break excitations, the field strength required for this is generally greater than those used in the literature on low-energy defibrillation. However, the high-intensity shocks used in standard defibrillation may elicit break excitation propagation from the coronary vasculature.

Reentry and Ectopic Pacemakers Emerge in a Three-Dimensional Model for a Slab of Cardiac Tissue with Diffuse Microfibrosis near the Percolation Threshold

Arrhythmias in cardiac tissue are generally associated with irregular electrical wave propagation in the heart. Cardiac tissue is formed by a discrete cell network, which is often heterogeneous. Recently, it was shown in simulations of two-dimensional (2D) discrete models of cardiac tissue that a wave crossing a fibrotic, heterogeneous region may produce reentry and transient or persistent ectopic activity provided the fraction of conducting connections is just above the percolation threshold. Here, we investigate the occurrence of these phenomena in three-dimensions by simulations of a discrete model representing a thin slab of cardiac tissue. This is motivated (i) by the necessity to study the relevance and properties of the percolation-related mechanism for the emergence of microreentries in three dimensions and (ii) by the fact that atrial tissue is quite thin in comparison with ventricular tissue. Here, we simplify the model by neglecting details of tissue anatomy, e. g. geometries of atria or ventricles and the anisotropy in the conductivity. Hence, our modeling study is confined to the investigation of the effect of the tissue thickness as well as to the comparison of the dynamics of electrical excitation in a 2D layer with the one in a 3D slab. Our results indicate a strong and non-trivial effect of the thickness even for thin tissue slabs on the probability of microreentries and ectopic beat generation. The strong correlation of the occurrence of microreentry with the percolation threshold reported earlier in 2D layers persists in 3D slabs. Finally, a qualitative agreement of 3D simulated electrograms in the fibrotic region with the experimentally observed complex fractional atrial electrograms (CFAE) as well as strong difference between simulated electrograms in 2D and 3D were found for the cases where reentry and ectopic activity were triggered by the micro-fibrotic region.

Novel Radiofrequency Ablation Strategies for Terminating Atrial Fibrillation in the Left Atrium: A Simulation Study

Pulmonary vein isolation (PVI) with radiofrequency ablation (RFA) is the cornerstone of atrial fibrillation (AF) therapy, but few strategies exist for when it fails. To guide RFA, phase singularity (PS) mapping locates reentrant electrical waves (rotors) that perpetuate AF. The goal of this study was to test existing and develop new RFA strategies for terminating rotors identified with PS mapping. It is unsafe to test experimental RFA strategies in patients, so they were evaluated in silico using a bilayer computer model of the human atria with persistent AF (pAF) electrical (ionic) and structural (fibrosis) remodeling. pAF was initiated by rapidly pacing the right (RSPV) and left (LSPV) superior pulmonary veins during sinus rhythm, and rotor dynamics quantified by PS analysis. Three RFA strategies were studied: (i) PVI, roof, and mitral lines; (ii) circles, perforated circles, lines, and crosses 0.5-1.5 cm in diameter/length administered near rotor locations/pathways identified by PS mapping; and (iii) 4-8 lines streamlining the sequence of electrical activation during sinus rhythm. As in pAF patients, 2 ± 1 rotors with cycle length 185 ± 4 ms and short PS duration 452 ± 401 ms perpetuated simulated pAF. Spatially, PS density had weak to moderate positive correlations with fibrosis density (RSPV: r = 0.38, p = 0.35, LSPV: r = 0.77, p = 0.02). RFA PVI, mitral, and roof lines failed to terminate pAF, but RFA perforated circles and lines 1.5 cm in diameter/length terminated meandering rotors from RSPV pacing when placed at locations with high PS density. Similarly, RFA circles, perforated circles, and crosses 1.5 cm in diameter/length terminated stationary rotors from LSPV pacing. The most effective strategy for terminating pAF was to streamline the sequence of activation during sinus rhythm with >4 RFA lines. These results demonstrate that co-localizing 1.5 cm RFA lesions with locations of high PS density is a promising strategy for terminating pAF rotors. For patients immune to PVI, roof, mitral, and PS guided RFA strategies, streamlining patient-specific activation sequences during sinus rhythm is a robust but challenging alternative.

Patient-derived models link re-entrant driver localization in atrial fibrillation to fibrosis spatial pattern

AIMS

The mechanisms underlying persistent atrial fibrillation (AF) in patients with atrial fibrosis are poorly understood. The goal of this study was to use patient-derived atrial models to test the hypothesis that AF re-entrant drivers (RDs) persist only in regions with specific fibrosis patterns.

METHODS AND RESULTS

Twenty patients with persistent AF (PsAF) underwent late gadolinium-enhanced MRI to detect the presence of atrial fibrosis. Segmented images were used to construct personalized 3D models of the fibrotic atria with biophysically realistic atrial electrophysiology. In each model, rapid pacing was applied to induce AF. AF dynamics were analysed and RDs were identified using phase mapping. Fibrosis patterns in RD regions were characterized by computing maps of fibrosis density (FD) and entropy (FE). AF was inducible in 13/20 models and perpetuated by few RDs (2.7 ± 1.5) that were spatially confined (trajectory of phase singularities: 7.6 ± 2.3 mm). Compared with the remaining atrial tissue, regions where RDs persisted had higher FE (IQR: 0.42-0.60 vs. 0.00-0.40, P < 0.05) and FD (IQR: 0.59-0.77 vs. 0.00-0.33, P < 0.05). Machine learning classified RD and non-RD regions based on FD and FE and identified a subset of fibrotic boundary zones present in 13.8 ± 4.9% of atrial tissue where 83.5 ± 2.4% of all RD phase singularities were located.

CONCLUSION

Patient-derived models demonstrate that AF in fibrotic substrates is perpetuated by RDs persisting in fibrosis boundary zones characterized by specific regional fibrosis metrics (high FE and FD). These results provide new insights into the mechanisms that sustain PsAF and could pave the way for personalized, MRI-based management of PsAF.

Optogenetics-enabled assessment of viral gene and cell therapy for restoration of cardiac excitability

Multiple cardiac pathologies are accompanied by loss of tissue excitability, which leads to a range of heart rhythm disorders (arrhythmias). In addition to electronic device therapy (i.e. implantable pacemakers and cardioverter/defibrillators), biological approaches have recently been explored to restore pacemaking ability and to correct conduction slowing in the heart by delivering excitatory ion channels or ion channel agonists. Using optogenetics as a tool to selectively interrogate only cells transduced to produce an exogenous excitatory ion current, we experimentally and computationally quantify the efficiency of such biological approaches in rescuing cardiac excitability as a function of the mode of application (viral gene delivery or cell delivery) and the geometry of the transduced region (focal or spatially-distributed). We demonstrate that for each configuration (delivery mode and spatial pattern), the optical energy needed to excite can be used to predict therapeutic efficiency of excitability restoration. Taken directly, these results can help guide optogenetic interventions for light-based control of cardiac excitation. More generally, our findings can help optimize gene therapy for restoration of cardiac excitability.

Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization

Impulse propagation in biological tissues is known to be modulated by structural heterogeneity. In cardiac muscle, improved understanding on how this heterogeneity influences electrical spread is key to advancing our interpretation of dispersion of repolarization. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a means of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, analysed against in vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies relevant characteristics of cardiac electrical propagation at tissue level. These include conduction effects on action potential (AP) morphology, the shortening of AP duration along the activation pathway and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media.

Effects of macroscopic heterogeneity on propagation in a computationally inexpensive 2D model of the heart

We have developed a computationally inexpensive, two-dimensional, bidomain model of the heart to demonstrate the effect of tissue heterogeneity on propagation of cardiac impulses generated by the sino-atrial node (SAN). The geometry consists of a thin sheet of cardiac tissue with designated areas that represent the SAN and atria. The SAN auto-generates continuous impulses that result in waves of normal propagation throughout the tissue. On the introduction of heterogeneous patches with low tissue conductivities, the rhythm of the waveform becomes irregular. The study suggests that simplified and computationally inexpensive models can be insightful tools to better understand the mechanisms that cause atrial fibrillation (AF) and hence more effective treatment methods.