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

Forward Euler Stability of the Bidomain Model of Cardiac Tissue

The bidomain model describes the anisotropic electrical properties of cardiac tissue. One common numerical technique for solving the bidomain equations is the explicit forward Euler method. In this communication we derive a relationship between the time and space steps that ensures the stability of this numerical method.

A model for multi-site pacing of fibrillation using nonlinear dynamics feedback

Traditionally, cardiac defibrillation requires a strong electric shock. Many unwanted side effects of this shock could be eliminated if defibrillation were performed using weak stimuli applied to several locations throughout the heart. Such multi-site pacing algorithms have been shown to defibrillate both experimentally (Pak et al., Am J Physiol 285:H2704-H2711, 2003) and theoretically (Puwal and Roth, J Biol Systems 14:101-112, 2006). Gauthier et al. (Chaos, 12:952-961, 2002) proposed a method to pace the heart using an algorithm based on nonlinear dynamics feedback applied through a single electrode. Our study applies a related but simpler algorithm, which essentially configures each electrode as a demand pacemaker, to simulate the multi-site pacing of fibrillating cardiac tissue. We use the numerical model developed by Fenton et al. (Chaos, 12:852-892, 2002) as the reaction term in a reaction-diffusion equation that we solve over a two-dimensional sheet of tissue. The defibrillation rate after pacing for 3 s is about 30%, which is significantly higher than the spontaneous defibrillation rate and is higher than observed in previous experimental and theoretical studies. Tuning the algorithm period can increase this rate to 45%.

Alternans and the influence of ionic channel modifications: Cardiac three-dimensional simulations and one-dimensional numerical bifurcation analysis

Cardiac propagation is investigated by simulations using a realistic three-dimensional (3D) geometry including muscle fiber orientation of the ventricles of a rabbit heart and the modified Beeler-Reuter ionic model. Electrical excitation is introduced by a periodic pacing of the lower septum. Depending on the pacing frequency, qualitatively different dynamics are observed, namely, normal heart beat, T-wave alternans, and 2:1 conduction block at small, intermediate, and large pacing frequencies, respectively. In a second step, we performed a numerical stability and bifurcation analysis of a pulse propagating in a one-dimensional (1D) ring of cardiac tissue. The precise onset of the alternans instability is obtained from computer-assisted linear stability analysis of the pulse and computation of the associated spectrum. The critical frequency at the onset of alternans and the profiles of the membrane potential agree well with the ones obtained in the 3D simulations. Next, we computed changes in the wave profiles and in the onset of alternans for the Beeler-Reuter model with modifications of the sodium, calcium, and potassium channels, respectively. For this purpose, we employ the method of numerical bifurcation and stability analysis. While blocking of calcium channels has a stabilizing effect, blocked sodium or potassium channels lead to the occurrence of alternans at lower pacing frequencies. The findings regarding channel blocking are verified within three-dimensional simulations. Altogether, we have found T-wave alternans and conduction block in 3D simulations of a realistic rabbit heart geometry. The onset of alternans has been analyzed by numerical bifurcation and stability analysis of 1D wave trains. By comparing the results of the two approaches, we find that alternans is not strongly influenced by ingredients such as 3D geometry and propagation anisotropy, but depends mostly on the frequency of pacing (frequency of subsequent action potentials). In addition, we have introduced numerical bifurcation and stability analysis as a tool into heart modeling and demonstrated its efficiency in scanning a large set of parameters in the case of models with reduced conductivity. Bifurcation analysis also provides an accurate test for analytical theories of alternans as is demonstrated for the case of the restitution hypothesis.

Negative filament tension in the Luo-Rudy model of cardiac tissue

Scroll waves are vortices that occur in three-dimensional excitable media. Scroll waves have been observed in a variety of systems including cardiac tissue, where they are associated with cardiac arrhythmias. The disorganization of scroll waves into chaotic behavior is thought to be the mechanism of ventricular fibrillation, which lethality is widely known. One of the possible mechanisms of scroll wave instability is negative filament tension, which was studied theoretically using low-dimensional models of excitable medium. In this article we perform a numerical study of negative filament tension using the Luo-Rudy phase 1 model, which is widely used in cardiac electrophysiology. We show that this instability exists in this model, study its manifestation and discuss its relation to cardiac arrhythmogenesis.

Turbulence control with local pacing and its implication in cardiac defibrillation

In this review article, we describe turbulence control in excitable systems by using a local periodic pacing method. The controllability conditions of turbulence suppression and the mechanisms underlying these conditions are analyzed. The local pacing method is applied to control Winfree turbulence (WT) and defect turbulence (DT) induced by spiral-wave breakup. It is shown that WT can always be suppressed by local pacing if the pacing amplitude and frequency are properly chosen. On the other hand, the pacing method can achieve suppression of DT induced by instabilities associated with the motions of spiral tips while failing to suppress DT induced by the instabilities of wave propagation far from tips. In the latter case, an auxiliary method of applying gradient field is suggested to improve the control effects. The implication of this local pacing method to realistic cardiac defibrillation is addressed.

Arrhythmogenesis in the heart: Multiscale modeling of the effects of defibrillation shocks and the role of electrophysiological heterogeneity

The mechanisms of initiation of ventricular arrhythmias as well as those behind the complex spatiotemporal wave dynamics and its filament organization during ventricular fibrillation (VF) are the topic of intense research and debate. Mechanistic inquiry into the various mechanisms that lead to arrhythmia initiation and VF maintenance is hampered by the inability of current experimental techniques to resolve, with sufficient accuracy, electrical behavior confined to the depth of the ventricles. The objective of this article is to demonstrate that realistic 3D simulations of electrical activity in the heart are capable of bringing a new level of understanding of the mechanisms that underlie arrhythmia initiation and subsequent organization. The article does this by presenting the results of two multiscale simulation studies of ventricular electrical behavior. The first study aims to uncover the mechanisms responsible for rendering the ventricles vulnerable to electric shocks during a specific interval of time, the vulnerable window. The second study focuses on elucidating the role of electrophysiological heterogeneity, and specifically, differences in action potential duration in various ventricular structures, in VF organization. Both studies share common multiscale modeling approaches and analysis, including characterization of scroll-wave filament dynamics.

An ionically based mapping model with memory for cardiac restitution

Many features of the sequence of action potentials produced by repeated stimulation of a patch of cardiac muscle can be modeled by a 1D mapping, but not the full behavior included in the restitution portrait. Specifically, recent experiments have found that (i) the dynamic and S1-S2 restitution curves are different (rate dependence) and (ii) the approach to steady state, which requires many action potentials (accommodation), occurs along a curve distinct from either restitution curve. Neither behavior can be produced by a 1D mapping. To address these shortcomings, ad hoc 2D mappings, where the second variable is a "memory" variable, have been proposed; these models exhibit qualitative features of the relevant behavior, but a quantitative fit is not possible. In this paper we introduce a new 2D mapping and determine a set of parameters for it that gives a quantitatively accurate description of the full restitution portrait measured from a bullfrog ventricle. The mapping can be derived as an asymptotic limit of an idealized ionic model in which a generalized concentration acts as a memory variable. This ionic basis clarifies how the present model differs from previous models. The ionic basis also provides the foundation for more extensive cardiac modeling: e.g., constructing a PDE model that may be used to study the effect of memory on propagation. The fitting procedure for the mapping is straightforward and can easily be applied to obtain a mathematical model for data from other experiments, including experiments on different species.

Theory of action potential wave block at-a-distance in the heart

Propagation failure of an action potential wave at a finite distance from its source (so-called type-II block) may cause spiral wave formation or wave breakup in the heart, phenomena that are believed to underlie lethal and nonlethal heart rhythm disorders. In this study, we develop a sufficient condition for this type of block in a homogeneous, spatially one-dimensional system. Using a topological argument, we find that type-II block of a wave will always occur when launched within a finite range of times if the velocity of the trailing edge of the preceding wave, as measured at the stimulus site, is smaller than the velocity of a wave launched with the minimum diastolic interval (DI) for which propagation is possible. This "blocking condition" is robust, remaining valid even when memory and waveback electrotonic effects are included. The condition suggests that type-II block is greatly facilitated when waves are initiated at irregular intervals in time such that (1) the velocities of consecutive waves are as different as possible and (2) the DIs preceding each wave fall on the steeply sloped portion of the action potential duration restitution curve as often as possible. The set of timing intervals between stimuli that are predicted by the blocking condition to produce block are found to be consistent with these guidelines, and also to agree well with a coupled-maps computer simulation model, for the case of waves launched by four rapidly and irregularly timed stimuli.

Spiral Wave Attachment to Millimeter-Sized Obstacles

Background— Functional reentry in the heart takes the form of spiral waves. Drifting spiral waves can become pinned to anatomic obstacles and thus attain stability and persistence. Lidocaine is an antiarrhythmic agent commonly used to treat ventricular tachycardia clinically. We examined the ability of small obstacles to anchor spiral waves and the effect of lidocaine on their attachment.

Methods and Results— Spiral waves were electrically induced in confluent monolayers of cultured, neonatal rat cardiomyocytes. Small, circular anatomic obstacles (0.6 to 2.6 mm in diameter) were situated in the center of the monolayers to provide an anchoring site. Eighty reentry episodes consisting of at least 4 revolutions were studied. In 36 episodes, the spiral wave attached to the obstacle and became stationary and sustained, with a shorter reentry cycle length and higher rate. Spiral waves could attach to obstacles as small as 0.6 mm, with a likelihood for attachment that increased with obstacle size. After attachment, both conduction velocity of the wave-front tip and wavelength near the obstacle adapted from their pre-reentry values and increased linearly with obstacle size. In contrast, reentry cycle length did not correlate significantly with obstacle size. Addition of lidocaine 90 μmol/L depressed conduction velocity, increased reentry cycle length, and caused attached spiral waves to become quasi- attached to the obstacle or terminate.

Conclusions— Anchored spiral waves exhibit properties of both unattached spiral waves and anatomic reentry. Their behavior may be representative of functional reentry dynamics in cardiac tissue, particularly in the setting of monomorphic tachyarrhythmias.

Negative-Tension Instability of Scroll Waves and Winfree Turbulence in the Oregonator Model

Excitable media support self-organized scroll waves which can be unstable and give rise to three-dimensional wave chaos. Winfree turbulence of scroll waves results from the negative-tension instability of scroll waves; it plays an important role in the cardiac tissue where it may lead to ventricular fibrillation. By numerical simulations of the Oregonator model, we show that this instability and, thus, the Winfree turbulence may also be observed in the Belousov-Zhabotinsky reaction. The region in the parameter space, where the instability takes place, is determined, and a relationship between the negative-tension instability and the meandering behavior of spiral waves is found. The application of global periodic forcing to control such turbulence in the Oregonator model is discussed.

A tale of two dogs: analyzing two models of canine ventricular electrophysiology

The extensive development of detailed mathematical models of cardiac myocyte electrophysiology in recent years has led to a proliferation of models, including many that model the same animal species and specific region of the heart and thus would be expected to have similar properties. In this paper we review and compare two recently developed mathematical models of the electrophysiology of canine ventricular myocytes. To clarify their similarities and differences, we also present studies using them in a range of preparations from single cells to two-dimensional tissue. The models are compared with each other and with new and previously published experimental results in terms of a number of their properties, including action potential morphologies; transmembrane currents during normal heart rates and during alternans; alternans onsets, magnitudes, and cessations; and reentry dynamics of spiral waves. Action potential applets and spiral wave movies for the two canine ventricular models are available online as supplemental material. We find a number of differences between the models, including their rate dependence, alternans dynamics, and reentry stability, and a number of differences compared with experiments. Differences between models of the same species and region of the heart are not unique to these canine models. Similar differences can be found in the behavior of two models of human ventricular myocytes and of human atrial myocytes. We provide several possible explanations for the differences observed in models of the same species and region of the heart and discuss the implications for the applicability of models in addressing questions of mechanism in cardiac electrophysiology.

Magnetic resonance imaging of chemical waves in porous media

Magnetic resonance imaging (MRI) provides a powerful tool for the investigation of chemical structures in optically opaque porous media, in which chemical concentration gradients can be visualized, and diffusion and flow properties are simultaneously determined. In this paper we give an overview of the MRI technique and review theory and experiments on the formation of chemical waves in a tubular packed bed reactor upon the addition of a nonlinear chemical reaction. MR images are presented of reaction-diffusion waves propagating in the three-dimensional (3D) network of channels in the reactor, and the 3D structure of stationary concentration patterns formed via the flow-distributed oscillation mechanism is demonstrated to reflect the local hydrodynamics in the packed bed. Possible future directions regarding the influence of heterogeneities on transport and reaction are discussed.

Spiral wave dynamics in excitable media with spherical geometries

We describe the spatial and temporal organization of spiral and scroll waves in spherical shells of different sizes and solid spheres. We present simulation results for the evolution of the dynamics and clustering of spiral waves as a function of the excitability of the medium. The excitability, topology, and size of the domain places restrictions on how single and multiarmed spiral waves are organized in space. The results in spherical geometries are compared with those in planar two-dimensional media. These studies are relevant to the dynamics of spiral waves in a variety of media including the heart, and chemical reactions on spherical surfaces.

Mechanisms of destabilization and early termination of spiral wave reentry in the ventricle by a class III antiarrhythmic agent, nifekalant

Nifekalant (NF) is a novel class III antiarrhythmic agent that is effective in preventing life-threatening ventricular tachycardia/fibrillation (VT/VF). We investigated mechanisms of destabilization and early termination of spiral-type reentrant VT by NF in a two-dimensional subepicardial myocardial layer of Langendorff-perfused rabbit hearts (n = 21) using a high-resolution optical action potential mapping system. During basic stimulation, NF (0.1 microM) caused uniform prolongation of action potential duration (APD) without affecting conduction velocity and an increase of APD restitution slope. VTs induced by direct current stimulation in the presence of NF were of shorter duration (VTs > 30 s: 2/54 NF vs. 19/93 control). During VTs in control (with visible rotors), the wave front chased its own tail with a certain distance (repolarized zone), and they seldom met each other. The average number of phase singularity (PS) points was 1.31 +/- 0.14 per 665 ms (n = 7). In the presence of NF, the wave front frequently encountered its own tail, causing a transient breakup of the spiral wave or sudden movement of the rotation center (spatial jump of PS). The average number of PS was increased to 1.63 +/- 0.22 per 665 ms (n = 7, P < 0.05) after NF. The mode of spontaneous termination of rotors in control was in most cases (9/10, 90.0%) the result of mutual annihilation of counterrotating wave fronts. With NF, rotors frequently terminated by wave front collision with the atrioventricular groove (12/19, 63.2%) or by trapping the spiral tip in a refractory zone (7/19, 36.8%). Destabilization and early termination of spiral wave reentry induced by NF are the result of a limited proportion of excitable tissue after modulation of repolarization.

Inwardly rotating spiral wave breakup in oscillatory reaction-diffusion media

The breakup of inwardly rotating spiral waves has been investigated in an oscillatory reaction-diffusion system near a Hopf bifurcation point. The breakup first occurred at the region far away from the core area, then gradually involved the whole medium by increasing the diffusion coefficient ratio between the two components of the oscillator system. With the approximation of the Complex-Ginzburg-Landau equation (CGLE), the criteria for the occurrence of the inwardly rotating spiral wave are examined theoretically. The analysis of the stability in the corresponding CGLE revealed that the breakup of the inward spiral wave was related to the Eckhaus instability.

Action potential duration restitution portraits of mammalian ventricular myocytes: role of calcium current

Construction of the action potential duration (APD) restitution portrait allows visualization of multiple aspects of the dynamics of periodically paced myocytes at various basic cycle lengths (BCLs). For the first time, we obtained the restitution portrait of isolated rabbit and guinea pig cardiac ventricular myocytes and analyzed the time constant, tau, of APD accommodation and the slopes of different types of restitution curves, Sdyn and S12, measured at varying BCLs. Our results indicate that both tau and the individual slopes are species and pacing dependent. In contrast, the mutual relationship between slopes Sdyn and S12 does not depend on pacing history, being a generic feature of the species. In addition, the maximum slope S12, measured in the restitution portrait at the lowest BCL, predicts the onset of alternans. Further, we investigated the role of the L-type calcium current, ICa-L, in the restitution portrait. We found that ICa-L dramatically affects APD accommodation, as well as the individual slopes Sdyn and S12 measured in the restitution portrait. However, peak calcium current plays a role only at small values of BCL. In conclusion, the results demonstrate that the restitution portrait is a powerful technique to investigate restitution properties of periodically paced cardiac myocytes and the onset of alternans, in particular. Moreover, the data also show that ICa-L plays a crucial role in multiple aspects of cardiac dynamics measured through the restitution portrait.

Dynamics of bound states of same-chirality spiral waves

We describe the dynamics of bound states of same-chirality spirals in a generic numerical model of an excitable medium. For each bound state, we analyze its tip trajectory patterns and determine its characteristic frequencies. We report two previously unidentified bound states: for spiral pairs, a state that exhibits alternating cycles of small and large distances between collisions (A2); for triplets, the first example of a meandering bound state (M3). In parameter space, A2 lies in between the previously described oscillating pairs (O2) and master-slave pairs (MS). We present numerical evidence that the transition O2-->A2 occurs via a supercritical period-doubling bifurcation, while the transition A2-->MS occurs via a symmetry breaking secondary Hopf bifurcation. A classification of all regimes according to dynamical systems theory exposes the wealth of phenomena exhibited by multiarmed spiral waves.

Mechanistic enquiry into the effect of increased pacing rate on the upper limit of vulnerability

The goal of this study is to investigate the mechanisms responsible for the increase in the upper limit of vulnerability (ULV; highest shock strength that induces arrhythmia) following the increase in pacing rate. To accomplish this goal, the study employs a three-dimensional bidomain finite element model of a slice through the canine ventricles. The preparation was paced eight times at a basic cycle length (BCL) of either 80 or 150ms followed by delivery of shocks of various strengths and timings. Our results demonstrate that the shock strength, which induced an arrhythmia 50% of the time, increased 20% for the faster pacing compared to the slower pacing. Analysis of the mechanisms underlying the increased vulnerability revealed that delayed post-shock activations originating in the tissue depths appear as breakthrough activations on the surfaces of the preparation following an isoelectric window (IW). However, the IW duration was consistently shorter in the faster-paced preparation. Consequently, breakthrough activations appeared on the surfaces of this preparation earlier, when the tissue was less recovered, resulting in higher probability of unidirectional block and reentry. This explains why shocks of the same strength were more likely to result in arrhythmia induction when delivered to a preparation that was rapidly paced.

Action potential morphology heterogeneity in the atrium and its effect on atrial reentry: a two-dimensional and quasi-three-dimensional study

Atrial fibrillation (AF) is believed to be perpetuated by recirculating spiral waves. Atrial structures are often characterized with action potentials of varying morphologies; however, the role of the structure-dependent atrial electrophysiological heterogeneity in spiral wave behaviour is not well understood. The purpose of this study is to determine the effect of action potential morphology heterogeneity associated with the major atrial structures in spiral wave maintenance. The present study also focuses on how this effect is further modulated by the presence of the inherent periodicity in atrial structure. The goals of the study are achieved through the simulation of electrical behaviour in a two-dimensional atrial tissue model that incorporates the representation of action potentials in various structurally distinct regions in the right atrium. Periodic boundary conditions are then imposed to form a cylinder (quasi three-dimensional), thus allowing exploration of the additional effect of structure periodicity on spiral wave behaviour. Transmembrane potential maps and phase singularity traces are analysed to determine effects on spiral wave behaviour. Results demonstrate that the prolonged refractoriness of the crista terminalis (CT) affects the pattern of spiral wave reentry, while the variation in action potential morphology of the other structures does not. The CT anchors the spiral waves, preventing them from drifting away. Spiral wave dynamics is altered when the ends of the sheet are spliced together to form a cylinder. The main effect of the continuous surface is the generation of secondary spiral waves which influences the primary rotors. The interaction of the primary and secondary spiral waves decreased as cylinder diameter increased.

Periodic forcing of scroll rings and control of Winfree turbulence in excitable media

By simulations of the Barkley model, action of uniform periodic nonresonant forcing on scroll rings and wave turbulence in three-dimensional excitable media is investigated. Sufficiently strong rapid forcing converts expanding scroll rings into the collapsing ones and suppresses the Winfree turbulence caused by the negative tension of wave filaments. Slow strong forcing has an opposite effect, leading to expansion of scroll rings and induction of the turbulence. These effects are explained in the framework of the phenomenological kinematic theory of scroll waves.

Proarrhythmic Potential of Mesenchymal Stem Cell Transplantation Revealed in an In Vitro Coculture Model

Background— Mesenchymal stem cells (MSCs) are bone marrow stromal cells that are in phase 1 clinical studies of cellular cardiomyoplasty. However, the electrophysiological effects of MSC transplantation have not been studied. Although improvement of ventricular function would represent a positive outcome of MSC transplantation, focal application of stem cells has the potential downside of creating inhomogeneities that may predispose the heart to reentrant arrhythmias. In the present study we use an MSC and neonatal rat ventricular myocyte (NRVM) coculture system to investigate potential proarrhythmic consequences of MSC transplantation into the heart.

Methods and Results— Human MSCs were cocultured with NRVMs in ratios of 1:99, 1:9, and 1:4 and optically mapped. We found that conduction velocity was decreased in cocultures compared with controls, but action potential duration (APD 80 ) was not affected. Reentrant arrhythmias were induced in 86% of cocultures containing 10% and 20% MSCs (n=36) but not in controls (n=7) or cocultures containing only 1% MSCs (n=4). Immunostaining, Western blot, and dye transfer revealed the presence of functional gap junctions involving MSCs.

Conclusions— Our results suggest that mixtures of MSCs and NRVMs can produce an arrhythmogenic substrate. The mechanism of reentry is probably increased tissue heterogeneity resulting from electric coupling of inexcitable MSCs with myocytes.

Whole heart action potential duration restitution properties in cardiac patients: a combined clinical and modelling study

Steep action potential duration (APD) restitution has been shown to facilitate wavebreak and ventricular fibrillation. The global APD restitution properties in cardiac patients are unknown. We report a combined clinical electrophysiology and computer modelling study to: (1) determine global APD restitution properties in cardiac patients; and (2) examine the interaction of the observed APD restitution with known arrhythmia mechanisms. In 14 patients aged 52-85 years undergoing routine cardiac surgery, 256 electrode epicardial mapping was performed. Activation-recovery intervals (ARI; a surrogate for APD) were recorded over the entire ventricular surface. Mono-exponential restitution curves were constructed for each electrode site using a standard S1-S2 pacing protocol. The median maximum restitution slope was 0.91, with 27% of all electrode sites with slopes<0.5, 29% between 0.5 and 1.0, and 20% between 1.0 and 1.5. Eleven per cent of restitution curves maintained slope>1 over a range of diastolic intervals of at least 30 ms; and 0.3% for at least 50 ms. Activation-recovery interval restitution was spatially heterogeneous, showing regional organization with multiple discrete areas of steep and shallow slope. We used a simplified computer model of 2-D cardiac tissue to investigate how heterogeneous APD restitution can influence vulnerability to, and stability of re-entry. Our model showed that heterogeneity of restitution can act as a potent arrhythmogenic substrate, as well as influencing the stability of re-entrant arrhythmias. Global epicardial mapping in humans showed that APD restitution slopes were organized into regions of shallow and steep slopes. This heterogeneous organization of restitution may provide a substrate for arrhythmia.

A normal form for excitable media

We present a normal form for traveling waves in one-dimensional excitable media in the form of a differential delay equation. The normal form is built around the well-known saddle-node bifurcation generically present in excitable media. Finite wavelength effects are captured by a delay. The normal form describes the behavior of single pulses in a periodic domain and also the richer behavior of wave trains. The normal form exhibits a symmetry preserving Hopf bifurcation which may coalesce with the saddle node in a Bogdanov-Takens point, and a symmetry-breaking spatially inhomogeneous pitchfork bifurcation. We verify the existence of these bifurcations in numerical simulations. The parameters of the normal form are determined and its predictions are tested against numerical simulations of partial differential equation models of excitable media with good agreement.

Role of repolarization restitution in the development of coarse and fine atrial fibrillation in the isolated canine right atria

INTRODUCTION

Although the role of action potential duration restitution (APD-R) in the initiation and maintenance of ventricular fibrillation (VF) has been the subject of numerous investigations, its role in the generation of atrial fibrillation (AF) is less well studied. The cellular and ionic basis for coarse versus fine AF is not well delineated.

METHODS AND RESULTS

We measured APD-R during acetylcholine-mediated AF as well as during pacing (standard and dynamic protocols) in crista teriminalis, pectinate muscle, superior vena cava, and appendage of isolated canine arterially perfused right atria (n = 15). Transmembrane action potential (TAP), pseudo-ECG, and isometric tension development were simultaneously recorded. Acetylcholine flattened APD-R measured by both standard and dynamic protocols, but promoted induction of AF. AF was initially coarse, converting to fine within 3-15 minutes of AF. Coarse, but not fine AF was associated with dramatic fluctuations in tension development, reflecting wide variations in intracellular calcium activity ([Ca(2+)](i)). During coarse AF, APD-R data displayed a cloud-like distribution pattern, with a wide range of maximum APD-R slope (from 1.21 to 0.35). A maximum APD-R slope >1 was observed only in crista terminalis (3/10). The APD-R relationship was relatively linear and flat during fine AF. Reduction of [Ca(2+)](i) was associated with fine AF whereas augmentation of [Ca(2+)](i) with coarse AF.

CONCLUSIONS

Our data indicate that while APD-R may have a limited role in the maintenance of coarse AF, it is unlikely to contribute to the maintenance of fine AF and that [Ca(2+)](i) dynamics determine the degree to which AF is coarse or fine.

Multistability of reentrant rhythms in an ionic model of a two-dimensional annulus of cardiac tissue

The dynamics of reentry in a model of a two-dimensional annulus of homogeneous cardiac tissue, with a Beeler-Reuter type formulation of the membrane ionic currents, is examined. The bifurcation structure of the sustained propagated solutions is described as a function of Rin and Rout, the inner and outer radii of the annulus. The transition from periodic to quasiperiodic reentry occurs at a critical Rin, which first diminishes and then saturates as Rout is increased. The reduction of the critical Rin is a consequence of the increase of the wave-front curvature. There is a range of Rin below the critical radius in which two distinct quasiperiodic solutions coexist. Each of these solutions disappears at its own specific value of Rin, and their annihilation is preceded by a new type of bifurcation leading to a regime of propagation with transient successive detachments of the wave front from the inner border of the annulus.

Regional differences in APD restitution can initiate wavebreak and re-entry in cardiac tissue: a computational study

Background

Regional differences in action potential duration (APD) restitution in the heart favour arrhythmias, but the mechanism is not well understood.

Methods

We simulated a 150 × 150 mm 2D sheet of cardiac ventricular tissue using a simplified computational model. We investigated wavebreak and re-entry initiated by an S1S2S3 stimulus protocol in tissue sheets with two regions, each with different APD restitution. The two regions had a different APD at short diastolic interval (DI), but similar APD at long DI. Simulations were performed twice; once with both regions having steep (slope > 1), and once with both regions having flat (slope < 1) APD restitution.

Results

Wavebreak and re-entry were readily initiated using the S1S2S3 protocol in tissue sheets with two regions having different APD restitution properties. Initiation occurred irrespective of whether the APD restitution slopes were steep or flat. With steep APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms with S1S2 of 250 ms, to 75 ms (S1S2 180 ms). With flat APD restitution, the range of S2S3 intervals resulting in wavebreak increased from 1 ms (S1S2 250 ms), to 21 ms (S1S2 340 ms) and then 11 ms (S1S2 400 ms).

Conclusion

Regional differences in APD restitution are an arrhythmogenic substrate that can be concealed at normal heart rates. A premature stimulus produces regional differences in repolarisation, and a further premature stimulus can then result in wavebreak and initiate re-entry. This mechanism for initiating re-entry is independent of the steepness of the APD restitution curve.

Reproducing Cardiac Restitution Properties Using the Fenton–Karma Membrane Model

Modern models of cardiac membranes are often highly complex and computationally expensive, particularly when used in long simulations of spatially extended models of cardiac tissue. Therefore, there is a need for simpler membrane models that preserve the features of the complex models deemed important. This communication describes an empirical procedure that was used to choose the parameters of the three-variable Fenton-Karma (FK3V) model to reproduce the restitution properties of the Courtemanche-Ramirez-Nattel model of atrial tissue. The resulting parameter values for the FK3V model and the sensitivity table for all its parameters are provided. Thus, this study gives insight into the behavior of the FK3V model and the effect of its parameters on restitution properties.

Suppress Winfree turbulence by local forcing excitable systems

The occurrence of Winfree turbulence is currently regarded as one of the principal mechanisms underlying cardiac fibrillation. We develop a local stimulation method that suppresses Winfree turbulence in three-dimensional excitable media. We find that Winfree turbulence can be effectively suppressed by locally injecting periodic signals to only a very small subset (around some surface region) of total space sites. Our method for the first time demonstrates the effectiveness of local low-amplitude periodic excitations in suppressing turbulence in 3D excitable media and has fundamental improvements in efficiency, convenience, and turbulence suppression speed compared with previous strategies. Therefore, it has great potential for developing into a practical low-amplitude defibrillation approach.

Dynamics of conduction blocks in a model of paced cardiac tissue

We study numerically the dynamics of conduction blocks using a detailed electrophysiological model. We find that this dynamics depends critically on the size of the paced region. Small pacing regions lead to stationary conduction blocks while larger pacing regions can lead to conduction blocks that travel periodically towards the pacing region. We show that this size-dependence dynamics can lead to a novel arrhythmogenic mechanism. Furthermore, we show that the essential phenomena can be captured in a much simpler coupled-map model.

Ionic determinants of functional reentry in a 2-D model of human atrial cells during simulated chronic atrial fibrillation

Recent studies suggest that atrial fibrillation (AF) is maintained by fibrillatory conduction emanating from a small number of high-frequency reentrant sources (rotors). Our goal was to study the ionic correlates of a rotor during simulated chronic AF conditions. We utilized a two-dimensional (2-D), homogeneous, isotropic sheet (5 x 5 cm(2)) of human atrial cells to create a chronic AF substrate, which was able to sustain a stable rotor (dominant frequency approximately 5.7 Hz, rosette-like tip meander approximately 2.6 cm). Doubling the magnitude of the inward rectifier K(+) current (I(K1)) increased rotor frequency ( approximately 8.4 Hz), and reduced tip meander (approximately 1.7 cm). This rotor stabilization was due to a shortening of the action potential duration and an enhanced cardiac excitability. The latter was caused by a hyperpolarization of the diastolic membrane potential, which increased the availability of the Na(+) current (I(Na)). The rotor was terminated by reducing the maximum conductance (by 90%) of the atrial-specific ultrarapid delayed rectifier K(+) current (I(Kur)), or the transient outward K(+) current (I(to)), but not the fast or slow delayed rectifier K(+) currents (I(Kr)/I(Ks)). Importantly, blockade of I(Kur)/I(to) prolonged the atrial action potential at the plateau, but not at the terminal phase of repolarization, which led to random tip meander and wavebreak, resulting in rotor termination. Altering the rectification profile of I(K1) also slowed down or abolished reentrant activity. In combination, these simulation results provide novel insights into the ionic bases of a sustained rotor in a 2-D chronic AF substrate.

New approaches for identifying antiarrhythmic drug targets

Sudden cardiac death, secondary to ventricular fibrillation (VF), remains the leading cause of death in many developed countries. Substantial experimental and theoretical support exists for the idea that VF is caused by spiral wave re-entry. The initiation and subsequent break-up of spiral waves have been linked to electrical alternans, a phenomenon typically associated with a steeply sloped restitution relationship. Interventions that reduce the slope of the restitution relationship have been shown to prevent the induction of VF and to terminate existing VF in experimental models. These results suggest that electrical restitution may be a promising new target for antiarrhythmic therapies.

Modeling wave propagation in realistic heart geometries using the phase-field method

We present a novel algorithm for modeling electrical wave propagation in anatomical models of the heart. The algorithm uses a phase-field approach that represents the boundaries between the heart muscle and the surrounding medium as a spatially diffuse interface of finite thickness. The chief advantage of this method is to automatically handle the boundary conditions of the voltage in complex geometries without the need to track the location of these boundaries explicitly. The algorithm is shown to converge accurately in nontrivial test geometries with no-flux (zero normal current) boundary conditions as the width of the diffuse interface becomes small compared to the width of the cardiac action potential wavefront. Moreover, the method is illustrated for anatomically realistic models of isolated rabbit and canine ventricles as well as human atria.

The sensitivity of the heart to static magnetic fields

Static magnetic fields induce flow potentials in arterial flows in and around the heart, that have been detected as distortions in the ECG. The resultant currents flowing through the myocardium could alter the rate or rhythm of the heart. No such changes have been seen in animal experiments, or with humans, in static fields up to 8 T. The possible effects of such currents induced by fields larger than 8 T on cardiac pacemaker rate, and arrhythmogenesis are reviewed, using virtual cardiac tissues-computational models of cardiac electrophysiology. Arrhythmogenesis can be by the initiation of ectopic beats, or by re-entry, whose probability of occurrence is increased by any increase in the electrical heterogeneity, in particular, the action potential duration heterogeneity of the ventricle. Focal ectopic activity would be readily detectable, but since re-entrant arrhythmias are very rare events, even a large increase in their probability of occurrence still leaves them unlikely to be observed. Both of these two arrhythmogenic mechanisms would show a steep sigmoidal, or threshold dependence on induced current intensity, with the threshold for increasing the vulnerability to re-entry less than the threshold for initiating activity. Failure to observe them at fields less than 8 T provides only a lower bound for any threshold for arrhythmogenesis.

Global organization of dynamics in oscillatory heterogeneous excitable media

An oscillatory heterogeneous excitable medium undergoes a transition from periodic target patterns to a bursting rhythm driven by the spontaneous initiation and termination of spiral waves as coupling or density is reduced. We illustrate these phenomena in monolayers of chick embryonic heart cells using calcium-sensitive fluorescent dyes. These results are modeled in a heterogeneous cellular automaton in which the neighborhood of interaction and cell density is modified. Parameters that give rise to bursting rhythms are organized in distinct zones in parameter space, leading to a global organization that should be applicable to the dynamics in a large class of excitable media.

A membrane model of electrically remodelled atrial myocardium derived from in vivo measurements*

Aims

Contemporary ionic-based membrane models are computationally expensive and are not intended to match the properties of a given experimental preparation. The aim of this work was to use measured restitution properties of electrically remodelled atrial tissue to develop a simplified membrane model based on the Fenton–Karma (FK) equations amenable to large-scale simulation of chronic atrial fibrillation (CAF).

Methods

Two membrane models, the FK-CAF and FK-CNTRL parameter sets, were developed to match action potential duration (APD) and conduction velocity (CV) restitution properties of rapid-pacing-induced electrically remodelled sheep atria and healthy atria, respectively. The models were tested by inducing reentry in a two-dimensional anisotropic monodomain and comparing the resulting cycle lengths (CL) with measured CLs.

Results

Parameters for the FK models were obtained that reproduced APD and CV restitution properties measured in the CAF and healthy sheep atria. Using the FK-CAF parameters, reentry was sustained in a 2.5 by 2.5 cm sheet with a CL = 91.0 ± 3.0 ms. Reentry (CL = 113.2 ± 5.2 ms) could only be sustained in the FK-CNTRL model after the tissue was first activated at a fast rate (136.5 ms).

Conclusions

The FK-CAF model is shown to approximate the restitution properties of remodelled sheep atria and can be used to simulate reentry with short CLs similar to those measured during AF episodes.

Of circles and spirals: Bridging the gap between the leading circle and spiral wave concepts of cardiac reentry

The “leading circle model” was the first detailed attempt at understanding the mechanisms of functional reentry, and remains a widely-used notion in cardiac electrophysiology. The “spiral wave” concept was developed more recently as a result of modern theoretical analysis and is the basis for consideration of reentry mechanisms in present biophysical theory. The goal of this paper is to present these models in a way that is comprehensible to both the biophysical and electrophysiology communities, with the idea of helping clinical and experimental electrophysiologists to understand better the spiral wave concept and of helping biophysicists to understand why the leading circle concept is so attractive and widely used by electrophysiologists. To this end, the main properties of the leading circle and spiral wave models of reentry are presented. Their basic assumptions and determinants are discussed and the predictions of the two concepts with respect to pharmacological responses of arrhythmias are reviewed. A major difference between them lies in the predicted responses to Na+-channel blockade, for which the spiral wave paradigm appears more closely to correspond to the results of clinical and experimental observations. The basis of this difference is explored in the context of the fundamental properties of the models.

Suppression of scroll wave turbulence by noise

Rotating scroll waves are dynamical spatiotemporal structures characteristic of three-dimensional active media. It is well known that, under low excitability conditions, scroll waves develop an intrinsically unstable dynamical regime that leads to a highly disorganized pattern of wave propagation. Such a "turbulent" state bears some resemblance to fibrillation states in cardiac tissue. We show here that this unstable regime can be controlled by using a spatially distributed random forcing superimposed on a control parameter of the system. Our results are obtained from numerical simulations but an explicit analytical argument that rationalizes our observations is also presented.

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 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, the wave-back velocity exhibits pronounced rate dependence and the wave-front velocity does not. Moreover, the discrepancy between wave-back velocity RCs is most significant for a small DI. We provide an analytical explanation of these results, using a system of coupled maps to relate the wave-front and wave-back velocities. Our calculations show that rate-dependent wave-back velocity can be present even if neither APD nor wave-front velocity exhibits rate dependence.

Preventing alternans-induced spiral wave breakup in cardiac tissue: an ion-channel-based approach

The detailed processes involved in spiral wave breakup, believed to be one major mechanism by which tachycardia evolves into fibrillation, are still poorly understood. This has rendered difficult the proper design of an efficient and practical control stimulus protocol to eliminate such events. In order to gain new insights into the underlying electrophysiological and dynamical mechanisms of breakup, we applied linear perturbation theory to a steadily rotating spiral wave in two spatial dimensions. The tissue was composed of cells modeled using the Fenton-Karma equations whose parameters were chosen to emphasize alternans as a primary mechanism for breakup. Along with one meandering mode, not just one but several unstable alternans modes were found with differing growth rates, frequencies, and spatial structures. As the conductance of the fast inward current was increased, the instability of the modes increased, consistent with increased meandering and propensity for spiral breakup in simulations. We also explored a promising new approach, based on the theory, for the design of an energy efficient electrical stimulus protocol to control spiral wave breakup. The novelty lies in addressing the problem directly at the ion channel level and taking advantage of the inherent two dimensional nature of the rotating wave. With the help of the eigenmode method, we were able to calculate the exact timing and amplitude of the stimulus, and locate it optimally to maximize efficiency. The analysis led to a special-case example that demonstrated that a single, properly timed stimulus can have a global effect, suppressing all growing alternans modes over the entire tissue, thus inhibiting spiral wave breakup.

Molecular mechanisms and global dynamics of fibrillation: an integrative approach to the underlying basis of vortex-like reentry

Art Winfree's scientific legacy has been particularly important to our laboratory whose major goal is to understand the mechanisms of ventricular fibrillation (VF). Here, we take an integrative approach to review recent studies on the manner in which nonlinear electrical waves organize to result in VF. We describe the contribution of specific potassium channel proteins and of the myocardial fiber structure to such organization. The discussion centers on data derived from a model of stable VF in the Langendorff-perfused guinea pig heart that demonstrates distinct patterns of organization in the left (LV) and right (RV) ventricles. Analysis of optical mapping data reveals that VF excitation frequencies are distributed throughout the ventricles in clearly demarcated domains. The highest frequency domains are found on the anterior wall of the LV at a location where sustained reentrant activity is present. The optical data suggest that a high frequency rotor that remains stationary in the LV is the mechanism that sustains VF in this model. Computer simulations predict that the inward rectifying potassium current (IK1) is an essential determinant of rotor stability and frequency, and patch-clamp results strongly suggest that the outward component of IK1 of cells in the LV is significantly larger than in the RV. Additional computer simulations and analytical procedures predict that the filaments of the reentrant activity (scroll waves) adopt a non-random configuration depending on fiber organization within the ventricular wall. Using the minimal principle we have concluded that filaments align with the trajectory of least resistance (i.e. the geodesic) between their endpoints. Overall, the data discussed have opened new and potentially exciting avenues of research on the possible role played by inward rectifier channels in the mechanism of VF, as well as the organization of its reentrant sources in three-dimensional cardiac muscle. Such an integrative approach may lead us toward an understanding of the molecular and structural basis of VF and hopefully to new preventative approaches.

Dynamics of virtual electrode-induced scroll-wave reentry in a 3D bidomain model

Functional reentry in the heart can be caused by a wave front of excitation rotating around its edge. Previous simulations on the basis of monodomain cable equations predicted the existence of self-sustained, vortex-like wave fronts (scroll waves) rotating around a filament in three dimensions. In our simulations, we used the more accurate bidomain model with modified Beeler-Reuter ionic kinetics to study the dynamics of scroll-wave filaments in a 16 x 8 x 1.5-mm slab of ventricular tissue with straight fibers. Wave fronts were identified as the areas with inward current. Their edges represented the filaments. Both transmural and intramural reentries with I- and U-shaped filaments, respectively, were obtained by the S1-S2 point stimulation protocol through the virtual electrode-induced phase singularity mechanism. The filaments meandered along elongated trajectories and tended to attach to the tissue boundaries exposed to air (no current flow) rather than to the bath (zero extracellular potential). They completely detached from electroporated (zero transmembrane potential) boundaries. In our simulations, the presence of the bath led to generation of only U-shaped filaments, which survived for the 1.5-mm-thick slab but not for the slabs of 0.5- or 3-mm thicknesses. Thus boundary conditions may be another determinant of the type and dynamics of reentry.

Expanding scroll rings and negative tension turbulence in a model of excitable media

Scroll waves in excitable media, described by the Barkley model, are studied. In the parameter region of weak excitability, negative tension of wave filaments is found. It leads to expansion of scroll rings and instability of wave filaments. A circular filament tends to stretch, bend, loop, and produce an expanding tangle that fills up the volume. The filament does not undergo fragmentation before it touches the boundaries. Statistical properties of such Winfree turbulence of scroll waves are numerically investigated.

Multiarm spirals in a two-dimensional cardiac substrate

A variety of chemical and biological nonlinear excitable media, including heart tissue, can support stable, self-organized waves of activity in a form of rotating single-arm spirals. In particular, heart tissue can support stationary and meandering spirals of electrical excitation, which have been shown to underlie different forms of cardiac arrhythmias. In contrast to single-arm spirals, stable multiarm spirals (multiple spiral waves that rotate in the same direction around a common organizing center) have not been demonstrated and studied yet in living excitable tissues. Here, we show that persistent multiarm spirals of electrical activity can be induced in monolayer cultures of neonatal rat heart cells by a short, rapid train of electrical point stimuli applied during single-arm-spiral activity. Stable formation is accomplished only in monolayers that show a relatively broad and steep dependence of impulse wavelength and propagation velocity on rate of excitation. The resulting multiarm spirals emit waves of electrical activity at rates faster than for single-arm spirals and exhibit two distinct behaviors, namely "arm-switching" and "tip-switching." The phenomenon of rate acceleration due to an increase in the number of spiral arms possibly may underlie the acceleration of functional reentrant tachycardias paced by a clinician or an antitachycardia device.

Non-Tikhonov asymptotic properties of cardiac excitability

Models of electric excitability of cardiac cells can be studied by singular perturbation techniques. To do this one should take into account parameters appearing in equations in nonstandard ways. The physical reason for this is near-perfect switch behavior of ionic current gates. This leads to a definition of excitability different from the currently accepted one. The asymptotic structure revealed by our analysis can be used to devise simplified caricature models, to obtain approximate analytical solutions, and to facilitate numerical simulations.

Spiral wave drift in an electric field and scroll wave instabilities

Here, I present the numerical computation of speed and direction of the drift of a spiral wave in an excitable medium in the presence of an electric field. The drift speed presents a strong variation close to the parameter value where the drift-speed component along the field direction from parallel becomes antiparallel. Using a simple phenomenological model and results from a numerical linear stability analysis of scroll waves, I show that this behavior can be attributed to a resonance of the meander modes with the translation modes of the spiral wave. Extending this phenomenological model to scroll waves also clarifies the link between the drift and long wavelength instabilities of scroll waves.

Determining the effects of memory and action potential duration alternans on cardiac restitution using a constant-memory restitution protocol

Restitution, the dependence of action potential duration (APD) on diastolic interval, may be causally linked to the vulnerability of cardiac tissue to certain types of arrhythmias. While a number of pacing protocols are commonly used to quantify the restitution relation, one of these, the dynamic protocol, may result in the occurrence of APD alternans. However, the effects of APD alternans, and the concomitant alternation in cardiac memory, on the restitution curve are currently not well understood. Alternans preceding a given action potential may cause that action potential to have a different duration from one preceded by action potentials of identical duration. This interaction of alternans and memory can result in a dynamic restitution curve that is not unique. To address this, we have developed a constant-memory restitution protocol that enables the experimenter or modeller to obtain unique, constant-memory restitution curves at all diastolic intervals. Using this protocol, we obtained unique restitution curves for two ionic models of the cardiac action potential in the absence of alternans at all diastolic intervals. A comparison of the unique constant-memory and non-unique dynamic restitution curves for the two models shows that the presence of alternans can significantly alter the shape of the restitution curve compared to when alternans is absent.

Dynamical effects of diffusive cell coupling on cardiac excitation and propagation: a simulation study

Cell coupling is considered to be important for cardiac action potential propagation and arrhythmogenesis. We carried out computer simulations to investigate the effects of stimulation strength and cell-to-cell coupling on action potential duration (APD) restitution, APD alternans, and stability of reentry in models of isolated cell, one-dimensional cable, and two-dimensional tissue. Phase I formulation of the Luo and Rudy action potential model was used. We found that stronger stimulation resulted in a shallower APD restitution curve and onset of APD alternans at a faster pacing rate. Reducing diffusive coupling between cells prolonged APD. Weaker diffusive currents along the direction of propagation steepened APD restitution and caused APD alternans to occur at a slower pacing rate in tissue. Diffusive current due to curvature changed APD but had little effect on APD restitution slope and onset of instability. Heterogeneous cell coupling caused APD inhomogeneities in space. Reduction in coupling strength either uniformly or randomly had little effect on the rotation period and stability of a reentry, but random cell decoupling slowed the rotation period and, thus, stabilized the reentry, preventing it from breaking up into multiple waves. Therefore, in addition to its effects on action potential conduction velocity, diffusive cell coupling also affects APD in a rate-dependent manner, causes electrophysiological heterogeneities, and thus modulates the dynamics of cardiac excitation. These effects are brought about by the modulation of ionic current activation and inactivation.