Scroll-wave dynamics in human cardiac tissue: lessons from a mathematical model with inhomogeneities and fiber architecture

Cardiac arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), are among the leading causes of death in the industrialized world. These are associated with the formation of spiral and scroll waves of electrical activation in cardiac tissue; single spiral and scroll waves are believed to be associated with VT whereas their turbulent analogs are associated with VF. Thus, the study of these waves is an important biophysical problem. We present a systematic study of the combined effects of muscle-fiber rotation and inhomogeneities on scroll-wave dynamics in the TNNP (ten Tusscher Noble Noble Panfilov) model for human cardiac tissue. In particular, we use the three-dimensional TNNP model with fiber rotation and consider both conduction and ionic inhomogeneities. We find that, in addition to displaying a sensitive dependence on the positions, sizes, and types of inhomogeneities, scroll-wave dynamics also depends delicately upon the degree of fiber rotation. We find that the tendency of scroll waves to anchor to cylindrical conduction inhomogeneities increases with the radius of the inhomogeneity. Furthermore, the filament of the scroll wave can exhibit drift or meandering, transmural bending, twisting, and break-up. If the scroll-wave filament exhibits weak meandering, then there is a fine balance between the anchoring of this wave at the inhomogeneity and a disruption of wave-pinning by fiber rotation. If this filament displays strong meandering, then again the anchoring is suppressed by fiber rotation; also, the scroll wave can be eliminated from most of the layers only to be regenerated by a seed wave. Ionic inhomogeneities can also lead to an anchoring of the scroll wave; scroll waves can now enter the region inside an ionic inhomogeneity and can display a coexistence of spatiotemporal chaos and quasi-periodic behavior in different parts of the simulation domain. We discuss the experimental implications of our study.     

Scroll wave dynamics in a three-dimensional cardiac tissue model: roles of restitution, thickness, and fiber rotation

Scroll wave (vortex) breakup is hypothesized to underlie ventricular fibrillation, the leading cause of sudden cardiac death. We simulated scroll wave behaviors in a three-dimensional cardiac tissue model, using phase I of the Luo-Rudy (LR1) action potential model. The effects of action potential duration (APD) restitution, tissue thickness, filament twist, and fiber rotation were studied. We found that APD restitution is the major determinant of scroll wave behavior and that instabilities arising from APD restitution are the main determinants of scroll wave breakup in this cardiac model. We did not see a "thickness-induced instability" in the LR1 model, but a minimum thickness is required for scroll breakup in the presence of fiber rotation. The major effect of fiber rotation is to maintain twist in a scroll wave, promoting filament bending and thus scroll breakup. In addition, fiber rotation induces curvature in the scroll wave, which weakens conduction and further facilitates wave break.     

Dynamics of intramural scroll waves in three-dimensional continuous myocardium with rotational anisotropy.

It has been suggested that reentrant activity in three-dimensional cardiac muscle may be organized as a scroll wave rotating around a singularity line called the filament. Experimental studies indicate that filaments are often concealed inside the ventricular wall and consequently, scroll waves do not manifest reentrant activity on the surface. Here we analyse how such concealed scroll waves are affected by a twisted anisotropy resulting from rotation of layers of muscle fibers inside the ventricular wall. We used a computer model of a ventricular slab (15x15x15 mm(3)) with a fiber twist of 120 degrees from endocardium to epicardium. The action potential was simulated using FitzHugh-Nagumo equations. Scroll waves with rectilinear filaments were initiated at various depths of the slab and at different angles with respect to fiber orientation. The analysis shows that independent of initial conditions, after a certain transitional period, the filament aligns with the local fiber orientation. The alignment of the filament is determined by the directional variations in cell coupling due to fiber rotation and by boundary conditions. Our findings provide a mechanistic explanation for the prevalence of intramural reentry over transmural reentry during polymorphic ventricular tachycardia and fibrillation.     

Nonequilibrium Arrhythmic States and Transitions in a Mathematical Model for Diffuse Fibrosis in Human Cardiac Tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to fibroblasts via a gap-junctional conductance , reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.     

Effect of the form and anisotropy of the left ventricle on the drift of spiral waves

Three-dimensional spiral waves of electrical excitation in the myocardium are sources of dangerous cardiac arrhythmias. In this work, the dynamics of spiral waves of electrical excitation were studied in a symmetric anatomical model of the human heart left ventricle and a realistic ionic cell model of the human ventricular myocardium. Three factors that affect the drift waves in the heart were compared for the first time: the geometry of the heart wall, myocardial anisotropy, and wave chirality. Cardiac anisotropy was identified as a main factor in determining the drift of spiral waves. In the isotropic case, the dynamics were determined by the wall thickness, but did not depend on the wave chirality. In the anisotropic case, chirality was found to play a crucial role.

     

Theoretical analysis of the magnetocardiographic pattern for reentry wave propagation in a three-dimensional human heart model

We present a computational study of reentry wave propagation using electrophysiological models of human cardiac cells and the associated magnetic field map of a human heart. We examined the details of magnetic field variation and related physiological parameters for reentry waves in two-dimensional (2-D) human atrial tissue and a three-dimensional (3-D) human ventricle model. A 3-D mesh system representing the human ventricle was reconstructed from the surface geometry of a human heart. We used existing human cardiac cell models to simulate action potential (AP) propagation in atrial tissue and 3-D ventricular geometry, and a finite element method and the Galerkin approximation to discretize the 3-D domain spatially. The reentry wave was generated using an S1-S2 protocol. The calculations of the magnetic field pattern assumed a horizontally layered conductor for reentry wave propagation in the 3-D ventricle. We also compared the AP and magnetocardiograph (MCG) magnitudes during reentry wave propagation to those during normal wave propagation. The temporal changes in the reentry wave motion and magnetic field map patterns were also analyzed using two well-known MCG parameters: the current dipole direction and strength. The current vector in a reentry wave forms a rotating spiral. We delineated the magnetic field using the changes in the vector angle during a reentry wave, demonstrating that the MCG pattern can be helpful for theoretical analysis of reentry waves.     

A model for human ventricular tissue

The experimental and clinical possibilities for studying cardiac arrhythmias in human ventricular myocardium are very limited. Therefore, the use of alternative methods such as computer simulations is of great importance. In this article we introduce a mathematical model of the action potential of human ventricular cells that, while including a high level of electrophysiological detail, is computationally cost-effective enough to be applied in large-scale spatial simulations for the study of reentrant arrhythmias. The model is based on recent experimental data on most of the major ionic currents: the fast sodium, L-type calcium, transient outward, rapid and slow delayed rectifier, and inward rectifier currents. The model includes a basic calcium dynamics, allowing for the realistic modeling of calcium transients, calcium current inactivation, and the contraction staircase. We are able to reproduce human epicardial, endocardial, and M cell action potentials and show that differences can be explained by differences in the transient outward and slow delayed rectifier currents. Our model reproduces the experimentally observed data on action potential duration restitution, which is an important characteristic for reentrant arrhythmias. The conduction velocity restitution of our model is broader than in other models and agrees better with available data. Finally, we model the dynamics of spiral wave rotation in a two-dimensional sheet of human ventricular tissue and show that the spiral wave follows a complex meandering pattern and has a period of 265 ms. We conclude that the proposed model reproduces a variety of electrophysiological behaviors and provides a basis for studies of reentrant arrhythmias in human ventricular tissue.     

Effects of Mechano-Electric Feedback on Scroll Wave Stability in Human Ventricular Fibrillation

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, −60 mV and −10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of −60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of −10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of −10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.     

Multistability property in cardiac ionic models of mammalian and human ventricular cells

The underlying mechanisms of irregular cardiac rhythms are still poorly understood. Many experimental and modeling studies are aimed at identifying factors which cause cardiac arrhythmias. However, a lack of understanding of heart rhythm dynamical properties makes it difficult to uncover precise mechanisms of electrical instabilities, and hence to predict the onset of heart rhythm disorders. We review and compare the existing methods of studying cardiac dynamics, including restitution protocol (S1-S2), dynamic restitution protocol and multistability test protocol (S1-CI-S2). We focus on cardiac cell dynamics to elucidate regularities of heart rhythm. We demonstrate the advantages of our newly proposed systematic approach of analysis of cardiac cell dynamics using mammalian Luo Rudy 1991 and human ventricular Ten Tusscher 2006 single cell models under healthy and diseased conditions such as altered K⁺ or Ca²⁺ related currents. We investigate the role of ionic properties and the shape of an action potential on the nonlinear dynamics of electrical processes in periodically stimulated cardiac cells. We show the existence of multistability property for human ventricular cells. Moreover, the multistability is proposed to be an intrinsic property of cardiac cells, and is also suggested to be one of the mechanisms which could underlie the sudden triggering of life-threatening ventricular arrhythmias in the human heart.     

Effect of Twisted Fiber Anisotropy in Cardiac Tissue on Ablation with Pulsed Electric Fields

BackgroundAblation of cardiac tissue with pulsed electric fields is a promising alternative to current thermal ablation methods, and it critically depends on the electric field distribution in the heart.MethodsWe developed a model that incorporates the twisted anisotropy of cardiac tissue and computed the electric field distribution in the tissue. We also performed experiments in rabbit ventricles to validate our model. We find that the model agrees well with the experimentally determined ablation volume if we assume that all tissue that is exposed to a field greater than 3 kV/cm is ablated. In our numerical analysis, we considered how tissue thickness, degree of anisotropy, and electrode configuration affect the geometry of the ablated volume. We considered two electrode configurations: two parallel needles inserted into the myocardium (“penetrating needles” configuration) and one circular electrode each on epi- and endocardium, opposing each other (“epi-endo” configuration).ResultsFor thick tissues (10 mm) and moderate anisotropy ratio (a = 2), we find that the geometry of the ablated volume is almost unaffected by twisted anisotropy, i.e. it is approximately translationally symmetric from epi- to endocardium, for both electrode configurations. Higher anisotropy ratio (a = 10) leads to substantial variation in ablation width across the wall; these variations were more pronounced for the penetrating needle configuration than for the epi-endo configuration.For thinner tissues (4 mm, typical for human atria) and higher anisotropy ratio (a = 10), the epi-endo configuration yielded approximately translationally symmetric ablation volumes, while the penetrating electrodes configuration was much more sensitive to fiber twist.ConclusionsThese results suggest that the epi-endo configuration will be reliable for ablation of atrial fibrillation, independently of fiber orientation, while the penetrating electrode configuration may experience problems when the fiber orientation is not consistent across the atrial wall.     

On the Presence of Affine Fibril and Fiber Kinematics in the Mitral Valve Anterior Leaflet

In this study, we evaluated the hypothesis that the constituent fibers follow an affine deformation kinematic model for planar collagenous tissues. Results from two experimental datasets were utilized, taken at two scales (nanometer and micrometer), using mitral valve anterior leaflet (MVAL) tissues as the representative tissue. We simulated MVAL collagen fiber network as an ensemble of undulated fibers under a generalized two-dimensional deformation state, by representing the collagen fibrils based on a planar sinusoidally shaped geometric model. The proposed approach accounted for collagen fibril amplitude, crimp period, and rotation with applied macroscopic tissue-level deformation. When compared to the small angle x-ray scattering measurements, the model fit the data well, with an r² = 0.976. This important finding suggests that, at the homogenized tissue-level scale of ∼1 mm, the collagen fiber network in the MVAL deforms according to an affine kinematics model. Moreover, with respect to understanding its function, affine kinematics suggests that the constituent fibers are largely noninteracting and deform in accordance with the bulk tissue. It also suggests that the collagen fibrils are tightly bounded and deform as a single fiber-level unit. This greatly simplifies the modeling efforts at the tissue and organ levels, because affine kinematics allows a straightforward connection between the macroscopic and local fiber strains. It also suggests that the collagen and elastin fiber networks act independently of each other, with the collagen and elastin forming long fiber networks that allow for free rotations. Such freedom of rotation can greatly facilitate the observed high degree of mechanical anisotropy in the MVAL and other heart valves, which is essential to heart valve function. These apparently novel findings support modeling efforts directed toward improving our fundamental understanding of tissue biomechanics in healthy and diseased conditions.     

A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.     

Strain-based estimation of time-dependent transmural myocardial architecture in the ovine heart

Left ventricular myofibers are connected by an extensive extracellular collagen matrix to form myolaminar sheets. Histological cardiac tissue studies have previously observed a pleated transmural distribution of sheets in the ovine heart, alternating sign of the sheet angle from epicardium to endocardium. The present study investigated temporal variations in myocardial fiber and sheet architecture during the cardiac cycle. End-diastolic histological measurements made at subepicardium, midwall, and subendocardium at an anterior-basal and a lateral-equatorial region of the ovine heart, combined with transmural myocardial Lagrangian strains, showed that the sheet angle but not the fiber angle varied temporally throughout the cardiac cycle. The magnitude of the sheet angle decreased during systole at all transmural depths at the anterior-basal site and at midwall and subendocardium depths at the lateral-equatorial site, making the sheets more parallel to the radial axis. These results support a previously suggested accordion-like wall-thickening mechanism of the myocardial sheets.     

Circadian Rhythms of Early Afterdepolarizations and Ventricular Arrhythmias in a Cardiomyocyte Model

Sudden cardiac arrest is a malfunction of the heart’s electrical system, typically caused by ventricular arrhythmias, that can lead to sudden cardiac death (SCD) within minutes. Epidemiological studies have shown that SCD and ventricular arrhythmias are more likely to occur in the morning than in the evening, and laboratory studies indicate that these daily rhythms in adverse cardiovascular events are at least partially under the control of the endogenous circadian timekeeping system. However, the biophysical mechanisms linking molecular circadian clocks to cardiac arrhythmogenesis are not fully understood. Recent experiments have shown that L-type calcium channels exhibit circadian rhythms in both expression and function in guinea pig ventricular cardiomyocytes. We developed an electrophysiological model of these cells to simulate the effect of circadian variation in L-type calcium conductance. In our simulations, we found that there is a circadian pattern in the occurrence of early afterdepolarizations (EADs), which are abnormal depolarizations during the repolarization phase of a cardiac action potential that can trigger fatal ventricular arrhythmias. Specifically, the model produces EADs in the morning, but not at other times of day. We show that the model exhibits a codimension-2 Takens-Bogdanov bifurcation that serves as an organizing center for different types of EAD dynamics. We also simulated a two-dimensional spatial version of this model across a circadian cycle. We found that there is a circadian pattern in the breakup of spiral waves, which represents ventricular fibrillation in cardiac tissue. Specifically, the model produces spiral wave breakup in the morning, but not in the evening. Our computational study is the first, to our knowledge, to propose a link between circadian rhythms and EAD formation and suggests that the efficacy of drugs targeting EAD-mediated arrhythmias may depend on the time of day that they are administered.     

Electromechanical and structural alterations in the aging rabbit heart and aorta

Aging increases the risk for arrhythmias and sudden cardiac death (SCD). We aimed at elucidating aging-related electrical, functional, and structural changes in the heart and vasculature that account for this heightened arrhythmogenic risk. Young (5-9 mo) and old (3.5-6 yr) female New Zealand White (NZW) rabbits were subjected to in vivo hemodynamic, electrophysiological, and echocardiographic studies as well as ex vivo optical mapping, high-field magnetic resonance imaging (MRI), and histochemical experiments. Aging increased aortic stiffness (baseline pulse wave velocity: young, 3.54 ± 0.36 vs. old, 4.35 ± 0.28 m/s, P < 0.002) and diastolic (end diastolic pressure-volume relations: 3.28 ± 0.5 vs. 4.95 ± 1.5 mmHg/ml, P < 0.05) and systolic (end systolic pressure-volume relations: 20.56 ± 4.2 vs. 33.14 ± 8.4 mmHg/ml, P < 0.01) myocardial elastances in old rabbits. Electrophysiological and optical mapping studies revealed age-related slowing of ventricular and His-Purkinje conduction (His-to-ventricle interval: 23 ± 2.5 vs. 31.9 ± 2.9 ms, P < 0.0001), altered conduction anisotropy, and a greater inducibility of ventricular fibrillation (VF, 3/12 vs. 7/9, P < 0.05) in old rabbits. Histochemical studies confirmed an aging-related increased fibrosis in the ventricles. MRI showed a deterioration of the free-running Purkinje fiber network in ventricular and septal walls in old hearts as well as aging-related alterations of the myofibrillar orientation and myocardial sheet structure that may account for this slowed conduction velocity. Aging leads to parallel stiffening of the aorta and the heart, including an increase in systolic stiffness and contractility and diastolic stiffness. Increasingly, anisotropic conduction velocity due to fibrosis and altered myofibrillar orientation and myocardial sheet structure may contribute to the pathogenesis of VF in old hearts. The aging rabbit model represents a useful tool for elucidating age-related changes that predispose the aging heart to arrhythmias and SCD.     

A kinematic approach for efficient and robust simulation of the cardiac beating motion

Computer simulation techniques for cardiac beating motions potentially have many applications and a broad audience. However, most existing methods require enormous computational costs and often show unstable behavior for extreme parameter sets, which interrupts smooth simulation study and make it difficult to apply them to interactive applications. To address this issue, we present an efficient and robust framework for simulating the cardiac beating motion. The global cardiac motion is generated by the accumulation of local myocardial fiber contractions. We compute such local-to-global deformations using a kinematic approach; we divide a heart mesh model into overlapping local regions, contract them independently according to fiber orientation, and compute a global shape that satisfies contracted shapes of all local regions as much as possible. A comparison between our method and a physics-based method showed that our method can generate motion very close to that of a physics-based simulation. Our kinematic method has high controllability; the simulated ventricle-wall-contraction speed can be easily adjusted to that of a real heart by controlling local contraction timing. We demonstrate that our method achieves a highly realistic beating motion of a whole heart in real time on a consumer-level computer. Our method provides an important step to bridge a gap between cardiac simulations and interactive applications.     

Patient-specific Modeling of the Heart: Estimation of Ventricular Fiber Orientations

Patient-specific simulations of heart (dys)function aimed at personalizing cardiac therapy are hampered by the absence of in vivo imaging technology for clinically acquiring myocardial fiber orientations. The objective of this project was to develop a methodology to estimate cardiac fiber orientations from in vivo images of patient heart geometries. An accurate representation of ventricular geometry and fiber orientations was reconstructed, respectively, from high-resolution ex vivo structural magnetic resonance (MR) and diffusion tensor (DT) MR images of a normal human heart, referred to as the atlas. Ventricular geometry of a patient heart was extracted, via semiautomatic segmentation, from an in vivo computed tomography (CT) image. Using image transformation algorithms, the atlas ventricular geometry was deformed to match that of the patient. Finally, the deformation field was applied to the atlas fiber orientations to obtain an estimate of patient fiber orientations. The accuracy of the fiber estimates was assessed using six normal and three failing canine hearts. The mean absolute difference between inclination angles of acquired and estimated fiber orientations was 15.4 °. Computational simulations of ventricular activation maps and pseudo-ECGs in sinus rhythm and ventricular tachycardia indicated that there are no significant differences between estimated and acquired fiber orientations at a clinically observable level.The new insights obtained from the project will pave the way for the development of patient-specific models of the heart that can aid physicians in personalized diagnosis and decisions regarding electrophysiological interventions.