Indeterminacy of Spatiotemporal Cardiac Alternans

Cardiac alternans, a beat-to-beat alternation in action potential duration (at the cellular level) or in ECG morphology (at the whole heart level), is a marker of ventricular fibrillation, a fatal heart rhythm that kills hundreds of thousands of people in the US each year. Investigating cardiac alternans may lead to a better understanding of the mechanisms of cardiac arrhythmias and eventually better algorithms for the prediction and prevention of such dreadful diseases. In paced cardiac tissue, alternans develops under increasingly shorter pacing period. Existing experimental and theoretical studies adopt the assumption that alternans in homogeneous cardiac tissue is exclusively determined by the pacing period. In contrast, we find that, when calcium-driven alternans develops in cardiac fibers, it may take different spatiotemporal patterns depending on the pacing history. Because there coexist multiple alternans solutions for a given pacing period, the alternans pattern on a fiber becomes unpredictable. Using numerical simulation and theoretical analysis, we show that the coexistence of multiple alternans patterns is induced by the interaction between electrotonic coupling and an instability in calcium cycling.Comment: 20 pages, 10 figures, to be published in Phys. Rev.

Complex patterns of spontaneous initiations and terminations of reentrant circulation in a loop of cardiac tissue

A two-component model is developed that consists of a discrete loop of cardiac cells that circulates action potentials together with a cardiac pacing mechanism. Physiological properties of cells such as restitutions of refractoriness and of conduction velocity are given via experimentally measured functions. The dynamics of circulating pulses and their interactions with the pacer are regulated by two threshold relations. Patterns of spontaneous initiations and terminations of reentry (SITR) generated by this system are studied through numerical simulations and analytical observations. These patterns can be regular or irregular; causes of irregularities are identified as the threshold bistability of reentrant circulation (T-bistability) and in some cases, also phase-resetting interactions with the pacer.Comment: 27 pages, 10 figures, 61 references; A version of this paper (same results) is to appear in the Journal of Theoretical Biology; arXiv V2 adds helpful commments to facilitate reading and corrects minor errors in presentatio

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

Rate-dependent propagation of cardiac action potentials in a one-dimensional fiber

Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting the onset of abnormal cardiac rhythms. However, it is known that different pacing protocols lead to different APD restitution curves (RCs). This phenomenon, known as APD rate-dependence, is a consequence of memory in the tissue. In addition to APD restitution, conduction velocity restitution also plays an important role in the spatiotemporal dynamics of cardiac tissue. We present new results concerning rate-dependent restitution in the velocity of propagating action potentials in a one-dimensional fiber. Our numerical simulations show that, independent of the amount of memory in the tissue, waveback velocity exhibits pronounced rate-dependence and the wavefront velocity does not. Moreover, the discrepancy between waveback velocity RCs is most significant for small DI. We provide an analytical explanation of these results, using a system of coupled maps to relate the wavefront and waveback velocities. Our calculations show that waveback velocity rate-dependence is due to APD restitution, not memory.Comment: 17 pages, 7 figure

Global alternans instability and its effect on non-linear wave propagation : dynamical Wenckebach block and self terminating spiral waves

The main mechanism of formation of reentrant cardiac arrhythmias is via formation of waveblocks at heterogeneities of cardiac tissue. We report that heterogeneity and the area of waveblock can extend itself in space and can result formation of new additional sources, or termination of existing sources of arrhythmias. This effect is based on a new form of instability, which we coin as global alternans instability (GAI). GAI is closely related to the so-called (discordant) alternans instability, however its onset is determined by the global properties of the APD-restitution curve and not by its slope. The APD-restitution curve relates the duration of the cardiac pulse (APD) to the time interval between the pulses, and can easily be measured in an experimental or even clinical setting. We formulate the conditions for the onset of GAI, study its manifestation in various 1D and 2D situations and discuss its importance for the onset of cardiac arrhythmias

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

ROLE OF CONDUCTION IN THE GENESIS OF ALTERNANS OF ACTION POTENTIAL DURATION IN A SIMULATED ONE DIMENSIONAL FIBER

Ventricular fibrillation is one of the leading causes for Sudden Cardiac Death and is characterized by multiple activation wavefronts. Multiple activation wavefronts originate from a reentrant circuit which requires the presence of a unidirectional block in the path of a propagating excitation wave. It has been proposed that at the cellular level beat to beat alternation in the action potential duration at rapid pacing rates can result in a conduction block. Various mechanisms have been postulated to show the mechanisms of alternans. We use simulated activation in a one dimensional tissue fiber to show the existence of a new mechanism via which alternans can result. We used a new pacing protocol to eliminate alternans at the pacing site, and thus eliminating restitution of action potential duration at this site to reveal existence of alternans down the fiber. Effects on alternans of manipulations of specific ionic currents such as the sodium current (INa), calcium current (ICaL), potassium current (Ikr) and of the diffusion co-efficient (Dx) which simulates reduced expression of connexin 43 were determined. Decrease in sodium conductance, i.e. in excitability by half caused the alternans to occur at the pacing site itself even though APD restitution was eliminated. An increase or decrease in calcium current (ICaL) eliminated alternans throughout the fiber. The use of a novel pacing approach in investigation of alternans, as in this study, furthers our understanding of the mechanism of alternans and may prove helpful in the development of better anti-arrhythmic drugs in the future