Vulnerability for reentry in a three dimensional model of human atria: a simulation study

Cellular heterogeneity and repolarisation across the atria: an in silico study

Mechanisms of atrial fibrillation and the susceptibility to reentries can be impacted by the repolarization across the atria. Studies into atrial fibrillation ignore cell-to-cell heterogeneity due to electrotonic coupling. Recent studies show that cellular variability may have a larger impact on electrophysiological behaviour than assumed. This paper aims to determine the impact of cellular heterogeneity on the repolarization phase across the AF remodelled atria. Using a population of models approach, 10 anatomically identical atrial models were created to include cellular heterogeneity. Atrial models were compared with an equivalent homogenous model. Activation, APD90, and repolarization maps were used to compare models. The impact of electrotonic coupling in the tissue was determined through a comparison of RMP, APD20, APD50, APD90, and triangulation between regional atrial tissue and the single cell populations. After calibration, cellular heterogeneity does not impact atrial depolarization. Repolarization patterns were significantly impacted by cellular heterogeneity, with the APD90 across the LA increasing due to heterogeneity and the reverse occurring in the RA. Electrotonic coupling caused a reduction in variability across all biomarkers but did not fully remove variability. Electrotonic coupling resulted in an increase in APD20 and APD50, and reduced triangulation compared to isolated cell populations. Heterogeneity also caused a reduction in triangulation compared with regionally homogeneous atria.

Signal Reconstruction of Pulmonary Vein Recordings Using a Phenomenological Mathematical Model: Application to Pulmonary Vein Isolation Therapy

Atrial fibrillation (AF), the most prevalent cardiac arrhythmia, is commonly initiated by ectopic beats originating from a small myocardial sleeve extending over the pulmonary veins. Pulmonary vein isolation therapy attempts to isolate the pulmonary veins from the left atrium by ablating tissue, commonly by using radiofrequency ablation. During this procedure, the cardiologist records electrical activity using a lasso catheter, and the activation pattern recorded is used as a guide toward which regions to ablate. However, poor contact between electrode and tissue can lead to important regions of electrical activity not being recorded in clinic. We reproduce these signals through the use of a phenomenological model of the cardiac action potential on a cylinder, which we fit to post-AF atrial cells, and model the bipolar electrodes of the lasso catheter by an approximation of the surface potential. The resulting activation pattern is validated by direct comparison with those of clinical recordings. A potential application of the model is to reconstruct the missing electrical activity, minimizing the impact of the information loss on the clinical procedure, and we present results to demonstrate this.

In-silico modeling of atrial repolarization in normal and atrial fibrillation remodeled state

Atrial fibrillation (AF) is the most common cardiac arrhythmia, and the total number of AF patients is constantly increasing. The mechanisms leading to and sustaining AF are not completely understood yet. Heterogeneities in atrial electrophysiology seem to play an important role in this context. Although some heterogeneities have been used in in-silico human atrial modeling studies, they have not been thoroughly investigated. In this study, the original electrophysiological (EP) models of Courtemanche et al., Nygren et al. and Maleckar et al. were adjusted to reproduce action potentials in 13 atrial regions. The parameter sets were validated against experimental action potential duration data and ECG data from patients with AV block. The use of the heterogeneous EP model led to a more synchronized repolarization sequence in a variety of 3D atrial anatomical models. Combination of the heterogeneous EP model with a model of persistent AF-remodeled electrophysiology led to a drastic change in cell electrophysiology. Simulated Ta-waves were significantly shorter under the remodeling. The heterogeneities in cell electrophysiology explain the previously observed Ta-wave effects. The results mark an important step toward the reliable simulation of the atrial repolarization sequence, give a deeper understanding of the mechanism of atrial repolarization and enable further clinical investigations.

Computational modeling of the human atrial anatomy and electrophysiology

This review article gives a comprehensive survey of the progress made in computational modeling of the human atria during the last 10 years. Modeling the anatomy has emerged from simple "peanut"-like structures to very detailed models including atrial wall and fiber direction. Electrophysiological models started with just two cellular models in 1998. Today, five models exist considering e.g. details of intracellular compartments and atrial heterogeneity. On the pathological side, modeling atrial remodeling and fibrotic tissue are the other important aspects. The bridge to data that are measured in the catheter laboratory and on the body surface (ECG) is under construction. Every measurement can be used either for model personalization or for validation. Potential clinical applications are briefly outlined and future research perspectives are suggested.

Eikonal-based initiation of fibrillatory activity in thin-walled cardiac propagation models

Reentrant arrhythmias can be simulated in electrophysiological models of electrical impulse propagation governed by a reaction-diffusion system. To facilitate the initiation of a large number of independent episodes of simulated arrhythmias with controllable level of complexity, a new approach is proposed for thin-walled geometries in which depolarization wave dynamics is essentially two-dimensional. Points representing phase singularities are first randomly distributed over the epicardial surface and are assigned a topological charge (direction of rotation). A qualitatively-correct phase map is then reconstructed on the whole surface by interpolation. The eikonal-diffusion equation is used to iteratively regularize the phase map based on a priori information on wavefront propagation. An initial condition for the reaction-diffusion model is created from the resulting phase map with multiple functional/anatomical reentries. Results in an atrial model demonstrate the ability to generate statistical realizations of the same dynamics and to vary the level of complexity measured by the number of phase singularities. A library of 100 simulations with an average number of phase singularities ranging from 1 to 10 is created. An extension to volumetric patient-specific atrial models including fiber orientation and a fast conducting system is presented to illustrate possible applications.