Anomalous drift of spiral waves in heterogeneous excitable media

Drift of sparse and dense spiral waves under joint external forces

Many methods have been employed to investigate the drift behaviors of spiral waves in an effort to understand and control their dynamics. Drift behaviors of sparse and dense spirals induced by external forces have been investigated, yet they remain incompletely understood. Here we employ joint external forces to study and control the drift dynamics. First, sparse and dense spiral waves are synchronized by the suitable external current. Then, under another weak current or heterogeneity, the synchronized spirals undergo a directional drift, and the dependence of their drift velocity on the strength and frequency of the joint external force is studied.

Arrhythmia Mechanism and Dynamics in a Humanized Mouse Model of Inherited Cardiomyopathy Caused by Phospholamban R14del Mutation

BACKGROUND

Arginine (Arg) 14 deletion (R14del) in the calcium regulatory protein phospholamban (hPLN^R14del) has been identified as a disease-causing mutation in patients with an inherited cardiomyopathy. Mechanisms underlying the early arrhythmogenic phenotype that predisposes carriers of this mutation to sudden death with no apparent structural remodeling remain unclear.

METHODS

To address this, we performed high spatiotemporal resolution optical mapping of intact hearts from adult knock-in mice harboring the human PLN^WT (wildtype [WT], n=12) or the heterozygous human PLN^R14del mutation (R14del, n=12) before and after ex vivo challenge with isoproterenol and rapid pacing.

RESULTS

Adverse electrophysiological remodeling was evident in the absence of significant structural or hemodynamic changes. R14del hearts exhibited increased arrhythmia susceptibility compared with wildtype. Underlying this susceptibility was preferential right ventricular action potential prolongation that was unresponsive to β-adrenergic stimulation. A steep repolarization gradient at the left ventricular/right ventricular interface provided the substrate for interventricular activation delays and ultimately local conduction block during rapid pacing. This was followed by the initiation of macroreentrant circuits supporting the onset of ventricular tachycardia. Once sustained, these circuits evolved into high-frequency rotors, which in their majority were pinned to the right ventricle. These rotors exhibited unique spatiotemporal dynamics that promoted their increased stability in R14del compared with wildtype hearts.

CONCLUSIONS

Our findings highlight the crucial role of primary electric remodeling caused by the hPLN^R14del mutation. These inherently arrhythmogenic features form the substrate for adrenergic-mediated VT at early stages of PLN^R14del induced cardiomyopathy.

Anisotropic shortening in the wavelength of electrical waves promotes onset of electrical turbulence in cardiac tissue: An in silico study

Several pathological conditions introduce spatial variations in the electrical properties of cardiac tissue. These variations occur as localized or distributed gradients in ion-channel functionality over extended tissue media. Electrical waves, propagating through such affected tissue, demonstrate distortions, depending on the nature of the ionic gradient in the diseased substrate. If the degree of distortion is large, reentrant activity may develop, in the form of rotating spiral (2d) and scroll (3d) waves of electrical activity. These reentrant waves are associated with the occurrence of lethal cardiac rhythm disorders, known as arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), which are believed to be common precursors of sudden cardiac arrest. By using state-of-the-art mathematical models for generic, and ionically-realistic (human) cardiac tissue, we study the detrimental effects of these ionic gradients on electrical wave propagation. We propose a possible mechanism for the development of instabilities in reentrant wave patterns, in the presence of ionic gradients in cardiac tissue, which may explain how one type of arrhythmia (VT) can degenerate into another (VF). Our proposed mechanism entails anisotropic reduction in the wavelength of the excitation waves because of anisotropic variation in its electrical properties, in particular the action potential duration (APD). We find that the variation in the APD, which we induce by varying ion-channel conductances, imposes a spatial variation in the spiral- or scroll-wave frequency ω. Such gradients in ω induce anisotropic shortening of wavelength of the spiral or scroll arms and eventually leads to instabilitites.

Emergent dynamics of spatio-temporal chaos in a heterogeneous excitable medium

Self-organized activation patterns in excitable media such as spiral waves and spatio-temporal chaos underlie dangerous cardiac arrhythmias. While the interaction of single spiral waves with different types of heterogeneity has been studied extensively, the effect of heterogeneity on fully developed spatio-temporal chaos remains poorly understood. We investigate how the complexity and stability properties of spatio-temporal chaos in the Bär-Eiswirth model of excitable media depend on the heterogeneity of the underlying medium. We employ different measures characterizing the chaoticity of the system and find that the spatial arrangement of multiple discrete lower excitability regions has a strong impact on the complexity of the dynamics. Varying the number, shape, and spatial arrangement of the heterogeneities, we observe strong emergent effects ranging from increases in chaoticity to the complete cessation of chaos, contrasting the expectation from the homogeneous behavior. The implications of our findings for the development and treatment of arrhythmias in the heterogeneous cardiac muscle are discussed.

Effect of myocyte-fibroblast coupling on the onset of pathological dynamics in a model of ventricular tissue

Managing lethal cardiac arrhythmias is one of the biggest challenges in modern cardiology, and hence it is very important to understand the factors underlying such arrhythmias. While early afterdepolarizations (EAD) of cardiac cells is known to be one such arrhythmogenic factor, the mechanisms underlying the emergence of tissue level arrhythmias from cellular level EADs is not fully understood. Another known arrhythmogenic condition is fibrosis of cardiac tissue that occurs both due to aging and in many types of heart diseases. In this paper we describe the results of a systematic in-silico study, using the TNNP model of human cardiac cells and MacCannell model for (myo)fibroblasts, on the possible effects of diffuse fibrosis on arrhythmias occurring via EADs. We find that depending on the resting potential of fibroblasts (VFR), M-F coupling can either increase or decrease the region of parameters showing EADs. Fibrosis increases the probability of occurrence of arrhythmias after a single focal stimulation and this effect increases with the strength of the M-F coupling. While in our simulations, arrhythmias occur due to fibrosis induced ectopic activity, we do not observe any specific fibrotic pattern that promotes the occurrence of these ectopic sources.

How Does the Xenopus laevis Embryonic Cell Cycle Avoid Spatial Chaos?

Theoretical studies have shown that a deterministic biochemical oscillator can become chaotic when operating over a sufficiently large volume and have suggested that the Xenopus laevis cell cycle oscillator operates close to such a chaotic regime. To experimentally test this hypothesis, we decreased the speed of the post-fertilization calcium wave, which had been predicted to generate chaos. However, cell divisions were found to develop normally, and eggs developed into normal tadpoles. Motivated by these experiments, we carried out modeling studies to understand the prerequisites for the predicted spatial chaos. We showed that this type of spatial chaos requires oscillatory reaction dynamics with short pulse duration and postulated that the mitotic exit in Xenopus laevis is likely slow enough to avoid chaos. In systems with shorter pulses, chaos may be an important hazard, as in cardiac arrhythmias, or a useful feature, as in the pigmentation of certain mollusk shells.

Electrical wave propagation in an anisotropic model of the left ventricle based on analytical description of cardiac architecture

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation.

Attractive and repulsive contributions of localized excitability inhomogeneities and elimination of spiral waves in excitable media

The attracting and repelling of spiral waves in a two-dimensional excitable medium in the presence of localized excitability inhomogeneities are studied. The choice of two effects depends on the comparison of excitabilities inside and outside the localized obstacle. We inspect the changes in attracting and repelling behaviors with respect to the size of the obstacle and the initial distance between the center of the spiral core and the obstacle. To understand the occurrence of these phenomena, we investigated the small v-value areas near the tip and the function of the wave velocity as the excitability parameter ε. Considering the attributes of the attractive obstacle, an eliminating scheme of spiral waves is proposed in which the attractive obstacle is rapidly moved at several fixed times. This method can avoid the high-amplitude and high-frequency stimulus in most conventional methods.

Evolution of spiral and scroll waves of excitation in a mathematical model of ischaemic border zone

Abnormal electrical activity from the boundaries of ischemic cardiac tissue is recognized as one of the major causes in generation of ischemia-reperfusion arrhythmias. Here we present theoretical analysis of the waves of electrical activity that can rise on the boundary of cardiac cell network upon its recovery from ischaemia-like conditions. The main factors included in our analysis are macroscopic gradients of the cell-to-cell coupling and cell excitability and microscopic heterogeneity of individual cells. The interplay between these factors allows one to explain how spirals form, drift together with the moving boundary, get transiently pinned to local inhomogeneities, and finally penetrate into the bulk of the well-coupled tissue where they reach macroscopic scale. The asymptotic theory of the drift of spiral and scroll waves based on response functions provides explanation of the drifts involved in this mechanism, with the exception of effects due to the discreteness of cardiac tissue. In particular, this asymptotic theory allows an extrapolation of 2D events into 3D, which has shown that cells within the border zone can give rise to 3D analogues of spirals, the scroll waves. When and if such scroll waves escape into a better coupled tissue, they are likely to collapse due to the positive filament tension. However, our simulations have shown that such collapse of newly generated scrolls is not inevitable and that under certain conditions filament tension becomes negative, leading to scroll filaments to expand and multiply leading to a fibrillation-like state within small areas of cardiac tissue.