Elizabeth M. Cherry

Effects of Pacing Site and Stimulation History on Alternans Dynamics and the Development of Complex Spatiotemporal Patterns in Cardiac Tissue

Alternans of action potential duration has been associated with T wave alternans and the development of arrhythmias because it produces large gradients of repolarization. However, little is known about alternans dynamics in large mammalian hearts. Using optical mapping to record electrical activations simultaneously from the epicardium and endocardium of 9 canine right ventricles, we demonstrate novel arrhythmogenic complex spatiotemporal dynamics. (i) Alternans predominantly develops first on the endocardium. (ii) The postulated simple progression from normal rhythm to concordant to discordant alternans is not always observed; concordant alternans can develop from discordant alternans as the pacing period is decreased. (iii) In contrast to smaller tissue preparations, multiple stationary nodal lines may exist and need not be perpendicular to the pacing site or to each other. (iv) Alternans has fully three-dimensional dynamics and the epicardium and endocardium can show significantly different dynamics: multiple nodal surfaces can be transmural or intramural and can form concave/convex surfaces resulting in islands of discordant alternans. (v) The complex spatiotemporal patterns observed during alternans are very sensitive to both the site of stimulation and the stimulation history. Alternans in canine ventricles not only exhibit larger amplitudes and persist for longer cycle length regimes compared to those found in smaller mammalian hearts, but also show novel dynamics not previously described that enhance dispersion and show high sensitivity to initial conditions. This indicates some underlying predisposition to chaos and can help to guide the design of new drugs and devices controlling and preventing arrhythmic events.

Low-energy control of electrical turbulence in the heart

Controlling the complex spatio-temporal dynamics underlying life-threatening cardiac arrhythmias such as fibrillation is extremely difficult, because of the nonlinear interaction of excitation waves in a heterogeneous anatomical substrate. In the absence of a better strategy, strong, globally resetting electrical shocks remain the only reliable treatment for cardiac fibrillation. Here we establish the relationship between the response of the tissue to an electric field and the spatial distribution of heterogeneities in the scale-free coronary vascular structure. We show that in response to a pulsed electric field, E, these heterogeneities serve as nucleation sites for the generation of intramural electrical waves with a source density ρ(E) and a characteristic time, τ, for tissue depolarization that obeys the power law τu2009∝u2009Eα. These intramural wave sources permit targeting of electrical turbulence near the cores of the vortices of electrical activity that drive complex fibrillatory dynamics. We show in vitro that simultaneous and direct access to multiple vortex cores results in rapid synchronization of cardiac tissue and therefore, efficient termination of fibrillation. Using this control strategy, we demonstrate low-energy termination of fibrillation in vivo. Our results give new insights into the mechanisms and dynamics underlying the control of spatio-temporal chaos in heterogeneous excitable media and provide new research perspectives towards alternative, life-saving low-energy defibrillation techniques.

Verification of cardiac tissue electrophysiology simulators using an N-version benchmark

Ongoing developments in cardiac modelling have resulted, in particular, in the development of advanced and increasingly complex computational frameworks for simulating cardiac tissue electrophysiology. The goal of these simulations is often to represent the detailed physiology and pathologies of the heart using codes that exploit the computational potential of high-performance computing architectures. These developments have rapidly progressed the simulation capacity of cardiac virtual physiological human style models; however, they have also made it increasingly challenging to verify that a given code provides a faithful representation of the purported governing equations and corresponding solution techniques. This study provides the first cardiac tissue electrophysiology simulation benchmark to allow these codes to be verified. The benchmark was successfully evaluated on 11 simulation platforms to generate a consensus gold-standard converged solution. The benchmark definition in combination with the gold-standard solution can now be used to verify new simulation codes and numerical methods in the future.

Properties of two human atrial cell models in tissue: Restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.

Toward real-time simulation of cardiac dynamics

We show that through careful and model-specific optimizations of their GPU implementations, simulations of realistic, detailed cardiac-cell models now can be performed in 2D and 3D in times that are close to real time using a desktop computer. Previously, large-scale simulations of detailed mathematical models of cardiac cells were possible only using supercomputers. In our study, we consider five different models of cardiac electrophysiology that span a broad range of computational complexity: the two-variable Karma model, the four-variable Bueno-Orovio-Cherry-Fenton model, the eight-variable Beeler-Reuter model, the 19-variable Ten Tusscher-Panfilov model, and the 67-variable Iyer-Mazhari-Winslow model. For each of these models, we treat both their single- and double-precision versions and demonstrate linear or even sub-linear growth in simulation times with an increase in the size of the grid used to model cardiac tissue. We also show that our GPU implementations of these models can increase simulation speeds to near real-time for simulations of complex spatial patterns indicative of cardiac arrhythmic disorders, including spiral waves and spiral wave breakup. The achievement of real-time applications without the need for supercomputers may, in the near term, facilitate the adoption of modeling-based clinical diagnostics and treatment planning, including patient-specific electrophysiological studies.

Effects of boundaries and geometry on the spatial distribution of action potential duration in cardiac tissue

Increased dispersion of action potential duration across cardiac tissue has long been considered an important substrate for the development of most electrical arrhythmias. Although this dispersion has been studied previously by characterizing the static intrinsic gradients in cellular electrophysiology and dynamical gradients generated by fast pacing, few studies have concentrated on dispersions generated solely by structural effects. Here we show how boundaries and geometry can produce spatially dependent changes in action potential duration (APD) in homogeneous and isotropic tissue, where all the cells have the same APD in the absence of diffusion. Electrotonic currents due to coupling within the tissue and at the tissue boundaries can generate dispersion, and the profile of this dispersion can change dramatically depending on tissue size and shape, action potential morphology, tissue dimensionality, and stimulus frequency and location. The dispersion generated by pure geometrical effects can be on the order of tens of milliseconds, enough under certain conditions to produce conduction blocks and initiate reentrant waves.

Mechanisms of ventricular arrhythmias: a dynamical systems-based perspective

Defining the cellular electrophysiological mechanisms for ventricular tachyarrhythmias is difficult, given the wide array of potential mechanisms, ranging from abnormal automaticity to various types of reentry and kk activity. The degree of difficulty is increased further by the fact that any particular mechanism may be influenced by the evolving ionic and anatomic environments associated with many forms of heart disease. Consequently, static measures of a single electrophysiological characteristic are unlikely to be useful in establishing mechanisms. Rather, the dynamics of the electrophysiological triggers and substrates that predispose to arrhythmia development need to be considered. Moreover, the dynamics need to be considered in the context of a system, one that displays certain predictable behaviors, but also one that may contain seemingly stochastic elements. It also is essential to recognize that even the predictable behaviors of this complex nonlinear system are subject to small changes in the state of the system at any given time. Here we briefly review some of the short-, medium-, and long-term alterations of the electrophysiological substrate that accompany myocardial disease and their potential impact on the initiation and maintenance of ventricular arrhythmias. We also provide examples of cases in which small changes in the electrophysiological substrate can result in rather large differences in arrhythmia outcome. These results suggest that an interrogation of cardiac electrical dynamics is required to provide a meaningful assessment of the immediate risk for arrhythmia development and for evaluating the effects of putative antiarrhythmic interventions.

A Quantitative Comparison of the Behavior of Human Ventricular Cardiac Electrophysiology Models in Tissue

Numerical integration of mathematical models of heart cell electrophysiology provides an important computational tool for studying cardiac arrhythmias, but the abundance of available models complicates selecting an appropriate model. We study the behavior of two recently published models of human ventricular action potentials, the Grandi-Pasqualini-Bers (GPB) and the OHara-Virág-Varró-Rudy (OVVR) models, and compare the results with four previously published models and with available experimental and clinical data. We find the shapes and durations of action potentials and calcium transients differ between the GPB and OVVR models, as do the magnitudes and rate-dependent properties of transmembrane currents and the calcium transient. Differences also occur in the steady-state and S1–S2 action potential duration and conduction velocity restitution curves, including a maximum conduction velocity for the OVVR model roughly half that of the GPB model and well below clinical values. Between single cells and tissue, both models exhibit differences in properties, including maximum upstroke velocity, action potential amplitude, and minimum diastolic interval. Compared to experimental data, action potential durations for the GPB and OVVR models agree fairly well (although OVVR epicardial action potentials are shorter), but maximum slopes of steady-state restitution curves are smaller. Although studies show alternans in normal hearts, it occurs only in the OVVR model, and only for a narrow range of cycle lengths. We find initiated spiral waves do not progress to sustained breakup for either model. The dominant spiral wave period of the GPB model falls within clinically relevant values for ventricular tachycardia (VT), but for the OVVR model, the dominant period is longer than periods associated with VT. Our results should facilitate choosing a model to match properties of interest in human cardiac tissue and to replicate arrhythmia behavior more closely. Furthermore, by indicating areas where existing models disagree, our findings suggest avenues for further experimental work.

Termination of equine atrial fibrillation by quinidine: an optical mapping study.

OBJECTIVEnTo perform the first optical mapping studies of equine atrium to assess the spatiotemporal dynamics of atrial fibrillation (AF) and of its termination by quinidine.nnnANIMALSnIntact, perfused atrial preparations obtained from four horses with normal cardiovascular examinations.nnnMATERIALS AND METHODSnAF was induced by a rapid pacing protocol with or without acetylcholine perfusion, and optical mapping was used to determine spatial dominant frequency distributions, electrical activity maps, and single-pixel optical signals. Following induction of AF, quinidine gluconate was perfused into the preparation and these parameters were monitored during quinidine-induced termination of AF.nnnRESULTSnEquine AF develops in the context of spatial gradients in action potential duration (APD) and diastolic interval (DI) that produce alternans, conduction block, and Wenckebach conduction in different regions at fast pacing rates. Quinidine terminates AF and prevents subsequent reinduction by reducing the maximal frequency and increasing frequency homogeneity.nnnCONCLUSIONSnHeterogeneity of APD and DI promote alternans and conduction block at fast pacing rates in the equine atrium, predisposing to the development of AF. Quinidine terminates AF by reducing maximum frequency and increasing frequency homogeneity. Our results are consistent with the hypothesis that quinidine increases effective refractory period, thereby decreasing frequency.

Teaching Cardiac Electrophysiology Modeling to Undergraduate Students: Laboratory Exercises and GPU Programming for the Study of Arrhythmias and Spiral Wave Dynamics

As part of a 3-wk intersession workshop funded by a National Science Foundation Expeditions in Computing award, 15 undergraduate students from the City University of New York(1) collaborated on a study aimed at characterizing the voltage dynamics and arrhythmogenic behavior of cardiac cells for a broad range of physiologically relevant conditions using an in silico model. The primary goal of the workshop was to cultivate student interest in computational modeling and analysis of complex systems by introducing them through lectures and laboratory activities to current research in cardiac modeling and by engaging them in a hands-on research experience. The success of the workshop lay in the exposure of the students to active researchers and experts in their fields, the use of hands-on activities to communicate important concepts, active engagement of the students in research, and explanations of the significance of results as the students generated them. The workshop content addressed how spiral waves of electrical activity are initiated in the heart and how different parameter values affect the dynamics of these reentrant waves. Spiral waves are clinically associated with tachycardia, when the waves remain stable, and with fibrillation, when the waves exhibit breakup. All in silico experiments were conducted by simulating a mathematical model of cardiac cells on graphics processing units instead of the standard central processing units of desktop computers. This approach decreased the run time for each simulation to almost real time, thereby allowing the students to quickly analyze and characterize the simulated arrhythmias. Results from these simulations, as well as some of the background and methodology taught during the workshop, is presented in this article along with the programming code and the explanations of simulation results in an effort to allow other teachers and students to perform their own demonstrations, simulations, and studies.