Filament behavior in a computational model of ventricular fibrillation in the canine heart

The aim of this paper was to quantify the behavior of filaments in a computational model of re-entrant ventricular fibrillation. We simulated cardiac activation in an anisotropic monodomain with excitation described by the Fenton-Karma model with Beeler-Reuter restitution, and geometry by the Auckland canine ventricle. We initiated re-entry in the left and right ventricular free walls, as well as the septum. The number of filaments increased during the first 1.5 s before reaching a plateau with a mean value of about 36 in each simulation. Most re-entrant filaments were between 10 and 20 mm long. The proportion of filaments touching the epicardial surface was 65%, but most of these were visible for much less than one period of re-entry. This paper shows that useful information about filament dynamics can be gleaned from models of fibrillation in complex geometries, and suggests that the interplay of filament creation and destruction may offer a target for antifibrillatory therap

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

Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity

It has become widely accepted that the most dangerous cardiac arrhythmias are due to re- entrant waves, i.e., electrical wave(s) that re-circulate repeatedly throughout the tissue at a higher frequency than the waves produced by the heart's natural pacemaker (sinoatrial node). However, the complicated structure of cardiac tissue, as well as the complex ionic currents in the cell, has made it extremely difficult to pinpoint the detailed mechanisms of these life-threatening reentrant arrhythmias. A simplified ionic model of the cardiac action potential (AP), which can be fitted to a wide variety of experimentally and numerically obtained mesoscopic characteristics of cardiac tissue such as AP shape and restitution of AP duration and conduction velocity, is used to explain many different mechanisms of spiral wave breakup which in principle can occur in cardiac tissue. Some, but not all, of these mechanisms have been observed before using other models; therefore, the purpose of this paper is to demonstrate them using just one framework model and to explain the different parameter regimes or physiological properties necessary for each mechanism (such as high or low excitability, corresponding to normal or ischemic tissue, spiral tip trajectory types, and tissue structures such as rotational anisotropy and periodic boundary conditions). Each mechanism is compared with data from other ionic models or experiments to illustrate that they are not model-specific phenomena. The fact that many different breakup mechanisms exist has important implications for antiarrhythmic drug design and for comparisons of fibrillation experiments using different species, electromechanical uncoupling drugs, and initiation protocols.Comment: 128 pages, 42 figures (29 color, 13 b&w

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that are triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media and their application to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact in cardiac arrhythmias.Peer ReviewedPreprin

The role of M cells and the long QT syndrome in cardiac arrhythmias: simulation studies of reentrant excitations using a detailed electrophysiological model

In this numerical study, we investigate the role of intrinsic heterogeneities of cardiac tissue due to M cells in the generation and maintenance of reentrant excitations using the detailed Luo-Rudy dynamic model. This model has been extended to include a description of the long QT 3 syndrome, and is studied in both one dimension, corresponding to a cable traversing the ventricular wall, and two dimensions, representing a transmural slice. We focus on two possible mechanisms for the generation of reentrant events. We first investigate if early-after-depolarizations occurring in M cells can initiate reentry. We find that, even for large values of the long QT strength, the electrotonic coupling between neighboring cells prevents early-after-depolarizations from creating a reentry. We then study whether M cell domains, with their slow repolarization, can function as wave blocks for premature stimuli. We find that the inclusion of an M cell domain can result in some cases in reentrant excitations and we determine the lifetime of the reentry as a function of the size and geometry of the domain and of the strength of the long QT syndrome

Effect of Heart Structure on Ventricular Fibrillation in the Rabbit: A Simulation Study

Ventricular fibrillation (VF) is a lethal condition that affects millions worldwide. The mechanism underlying VF is unstable reentrant electrical waves rotating around lines called filaments. These complex spatio-temporal patterns can be studied using both experimental and numerical methods. Computer simulations provide unique insights including high resolution dynamics throughout the heart and systematic control of quantities such as fiber orientation and cellular kinetics that are not feasible experimentally. Here we study filament dynamics using two bi-ventricular 3-D high-resolution rabbit heart geometries, one with detailed fine structure and another without fine structure. We studied filament dynamics using anisotropic and isotropic conductivities, and with four cellular action potential models with different recovery kinetics. Spiral wave dynamics observed in isotropic two-dimensional sheets were not predictive of the behavior in the whole heart. In 2-D the four cell models exhibited stable reentry, meandering spiral waves, and spiral-wave breakup. In the whole heart with fine structure, all simulation results exhibited complex dynamics reminiscent of fibrillation observed experimentally. In the whole heart without fine structure, anisotropy acted to destabilize filament dynamics although the number of filaments was reduced compared to the heart with structure. In addition, in isotropic hearts without structure the two cell models that exhibited meandering spiral waves in 2-D, stabilized into figure-of-eight surface patterns. We also studied the sensitivity of filament dynamics to computer system configuration and initial conditions. After large simulation times, different macroscopic results sometimes occurred across different system configurations, likely due to a lack of bitwise reproducibility. The study conclusions were insensitive to initial condition perturbations, however, the exact number of filaments over time and their trends were altered by these changes. In summary, we present the following new results. First, we provide a new cell model that resembles the surface patterns of VF in the rabbit heart both qualitatively and quantitatively. Second, filament dynamics in the whole heart cannot be predicted from spiral wave dynamics in 2-D and we identified anisotropy as one destabilizing factor. Third, the exact dynamics of filaments are sensitive to a variety of factors, so we suggest caution in their interpretation and their quantitative analyses

Anisotropic Cardiac Conduction.

Anisotropy is the property of directional dependence. In cardiac tissue, conduction velocity is anisotropic and its orientation is determined by myocyte direction. Cell shape and size, excitability, myocardial fibrosis, gap junction distribution and function are all considered to contribute to anisotropic conduction. In disease states, anisotropic conduction may be enhanced, and is implicated, in the genesis of pathological arrhythmias. The principal mechanism responsible for enhanced anisotropy in disease remains uncertain. Possible contributors include changes in cellular excitability, changes in gap junction distribution or function and cellular uncoupling through interstitial fibrosis. It has recently been demonstrated that myocyte orientation may be identified using diffusion tensor magnetic resonance imaging in explanted hearts, and multisite pacing protocols have been proposed to estimate myocyte orientation and anisotropic conduction in vivo. These tools have the potential to contribute to the understanding of the role of myocyte disarray and anisotropic conduction in arrhythmic states

Cardiac cell modelling: Observations from the heart of the cardiac physiome project

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field

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