Anatomically-induced Fibrillation in a 3D model of the Human Atria

International audienceAtrial fibrillation (AF) requires both a trigger and a sub-strate that can maintain a complex reentrant activity. In patients and in experimental models this substrate is provided by both electrical and structural remodeling. Since these processes overlap in time it is impossible to assess their individual contributions to AF maintenance experimentally. Therefore we studied the effect of electrical re-modeling alone on AF initiation in a realistic numerical model of the human atria. We attempted to initiate AF by rapid pacing in 10 different locations, both with and without electrical remodeling. The protocols were repeated twice, with small variations in calcium conductivity, so that in total 30 simulations with and 30 simulations without remodeling were performed. In models with electrical remodeling, functional conduction block at structural in-homogeneities induced AF in 27 % of the simulations. In models without electrical remodeling, AF could not be induced. We conclude that in the complex anatomy of the atria electrical remodeling alone increases the probability of AF substantially. This finding supports a mechanism whereby electrical remodeling, which occurs relatively fast, accelerates the slower but irreversible structural remodeling process

Role of Fiber Orientation in Atrial Arrythmogenesis

Electrical wave-front propagation in the atria is determined largely by local fiber orientation. Recent study suggests that atrial fibrillation (AF) progresses with enhanced anisotropy. In this work, a 3D rabbit atrial anatomical model at 20 × 20 × 20 μm3 resolution with realistic fiber orientation was constructed based on the novel contrast-enhanced micro-CT imaging. The Fenton-Karma cellular activation model was adapted to reproduce rabbit atrial action potential period of 80 ms. Diffusivities were estimated for longitudinal and transverse directions of the fiber orientation respectively. Pacing was conducted in the 3D anisotropic atrial model with a reducing S2 interval to facilitate initiation of atrial arrhythmia. Multiple simulations were conducted with varying values of diffusion anisotropy and stimulus locations to evaluate the role of anisotropy in initiating AF. Under physiological anisotropy conditions, a rapid right atrial activation was followed by the left atrial activation. Excitation waves reached the atrio-ventricular border where they terminated. Upon reduction of conduction heterogeneity, re-entry was initiated by the rapid pacing and the activation of both atrial chambers was almost simultaneous. Myofiber orientation is an effective mechanism for regulating atrial activation. Modification of myoarchitecture is proarrhythmic

Directed networks as a novel way to describe and analyze cardiac excitation : directed graph mapping

Networks provide a powerful methodology with applications in a variety of biological, technological and social systems such as analysis of brain data, social networks, internet search engine algorithms, etc. To date, directed networks have not yet been applied to characterize the excitation of the human heart. In clinical practice, cardiac excitation is recorded by multiple discrete electrodes. During (normal) sinus rhythm or during cardiac arrhythmias, successive excitation connects neighboring electrodes, resulting in their own unique directed network. This in theory makes it a perfect fit for directed network analysis. In this study, we applied directed networks to the heart in order to describe and characterize cardiac arrhythmias. Proof-of-principle was established using in-silico and clinical data. We demonstrated that tools used in network theory analysis allow determination of the mechanism and location of certain cardiac arrhythmias. We show that the robustness of this approach can potentially exceed the existing state-of-the art methodology used in clinics. Furthermore, implementation of these techniques in daily practice can improve the accuracy and speed of cardiac arrhythmia analysis. It may also provide novel insights in arrhythmias that are still incompletely understood

Cardiac re-entry dynamics & self-termination in DT-MRI based model of Human Foetal Heart

The effect of heart geometry and anisotropy on cardiac re-entry dynamics and self-termination is studied here in anatomically realistic computer simulations of human foetal heart. 20 weeks of gestational age human foetal heart isotropic and anisotropic anatomy models from diffusion tensor MRI data sets are used in the computer simulations. The fibre orientation angles of the heart were obtained from the DT-MRI primary eigenvalues. In a spatially homogeneous electrophysiological mono domain model with the DT-MRI based heart geometries, we initiate simplified Fitz-Hugh-Nagumo kinetics cardiac re-entry at a prescribed location in a 2D slice, and in the full 3D anatomy model. In a slice of the heart, the MRI based fibre anisotropy changes the re-entry dynamics from pinned to anatomical re-entry. In the full 3D MRI based model, the foetal heart fibre anisotropy changes the re-entry dynamics from a persistent re-entry to the re-entry self-termination

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in IK1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on IK1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of IK1 and IKr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas IKr block alone may be a safer alternative for the E299V mutation. Combined inhibition of IKr and IKur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome:Insights from virtual human atria

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in IK1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on IK1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of IK1 and IKr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas IKr block alone may be a safer alternative for the E299V mutation. Combined inhibition of IKr and IKur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF

Cardiac re-entry dynamics and self-termination in DT-MRI based model of Human Foetal Heart

The effect of human fetal heart geometry and anisotropy on anatomy induced drift and self-termination of cardiac re-entry is studied here in MRI based 2D slice and 3D whole heart computer simulations. Isotropic and anisotropic models of 20 weeks of gestational age human fetal heart obtained from 100 μm voxel diffusion tensor MRI data sets were used in the computer simulations. The fiber orientation angles of the heart were obtained from the orientation of the DT-MRI primary eigenvectors. In a spatially homogeneous electrophysiological monodomain model with the DT-MRI based heart geometries, cardiac re-entry was initiated at a prescribed location in a 2D slice, and in the 3D whole heart anatomy models. Excitation was described by simplified FitzHugh-Nagumo kinetics. In a slice of the heart, with propagation velocity twice as fast along the fibers than across the fibers, DT-MRI based fiber anisotropy changes the re-entry dynamics from pinned to an anatomical re-entry. In the 3D whole heart models, the fiber anisotropy changes cardiac re-entry dynamics from a persistent re-entry to the re-entry self-termination. The self-termination time depends on the re-entry's initial position. In all the simulations with the DT-MRI based cardiac geometry, the anisotropy of the myocardial tissue shortens the time to re-entry self-termination several folds. The numerical simulations depend on the validity of the DT-MRI data set used. The ventricular wall showed the characteristic transmural rotation of the helix angle of the developed mammalian heart, while the fiber orientation in the atria was irregula

Influence of cardiac tissue anisotropy on re-entrant activation in computational models of ventricular fibrillation

The aim of this study was to establish the role played by anisotropic diffusion in (i) the number of filaments and epicardial phase singularities that sustain ventricular fibrillation in the heart, (ii) the lifetimes of filaments and phase singularities, and (iii) the creation and annihilation dynamics of filaments and phase singularities. A simplified monodomain model of cardiac tissue was used, with membrane excitation described by a simplified 3-variable model. The model was configured so that a single re-entrant wave was unstable, and fragmented into multiple re-entrant waves. Re-entry was then initiated in tissue slabs with varying anisotropy ratio. The main findings of this computational study are: (i) anisotropy ratio influenced the number of filaments Sustaining simulated ventricular fibrillation, with more filaments present in simulations with smaller values of transverse diffusion coefficient, (ii) each re-entrant filament was associated with around 0.9 phase singularities on the surface of the slab geometry, (iii) phase singularities were longer lived than filaments, and (iv) the creation and annihilation of filaments and phase singularities were linear functions of the number of filaments and phase singularities, and these relationships were independent of the anisotropy ratio. This study underscores the important role played by tissue anisotropy in cardiac ventricular fibrillation