Alternans and higher-order rhythms in an ionic model of a sheet of ischemic ventricular muscle

Mechanism for Amplitude Alternans in Electrocardiograms and the Initiation of Spatiotemporal Chaos

It is widely believed that one major life-threatening transition to chaotic fibrillation occurs via spiral-wave breakup that is preceded by spatiotemporal dispersion of refractoriness due to alternations in the duration of the cardiac action potential (AP). However, recent clinical and experimental evidence suggests that other characteristics of the AP may contribute to, and perhaps drive, this dangerous dynamical instability. To identify the relative roles of AP characteristics, we performed experiments in rabbit hearts under conditions to minimize AP duration dynamics which unmasked pronounced AP amplitude alternans just before the onset of fibrillation. We used a simplified ionic cell model to derive a return map and a stability condition that elucidates a novel underlying mechanism for AP alternans and spiral breakup. We found that inactivation of the sodium current is key to developing amplitude alternans and is directly connected to conduction block and initiation of arrhythmias. Simulations in 2D where AP amplitude alternation led to turbulence confirm our hypothesis.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Bistable dynamics underlying excitability of ion homeostasis in neuron models

When neurons fire action potentials, dissipation of free energy is usually not directly considered, because the change in free energy is often negligible compared to the immense reservoir stored in neural transmembrane ion gradients and the long-term energy requirements are met through chemical energy, i.e., metabolism. However, these gradients can temporarily nearly vanish in neurological diseases, such as migraine and stroke, and in traumatic brain injury from concussions to severe injuries. We study biophysical neuron models based on the Hodgkin-Huxley (HH) formalism extended to include time-dependent ion concentrations inside and outside the cell and metabolic energy-driven pumps. We reveal the basic mechanism of a state of free energy-starvation (FES) with bifurcation analyses showing that ion dynamics is for a large range of pump rates bistable without contact to an ion bath. This is interpreted as a threshold reduction of a new fundamental mechanism of ionic excitability that causes a long-lasting but transient FES as observed in pathological states. We can in particular conclude that a coupling of extracellular ion concentrations to a large glial-vascular bath can take a role as an inhibitory mechanism crucial in ion homeostasis, while the Na⁺/K⁺ pumps alone are insufficient to recover from FES. Our results provide the missing link between the HH formalism and activator-inhibitor models that have been successfully used for modeling migraine phenotypes, and therefore will allow us to validate the hypothesis that migraine symptoms are explained by disturbed function in ion channel subunits, Na⁺/K⁺ pumps, and other proteins that regulate ion homeostasis.

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.

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.

Phase resetting, phase locking, and bistability in the periodically driven saline oscillator: experiment and model

The saline oscillator consists of an inner vessel containing salt water partially immersed in an outer vessel of fresh water, with a small orifice in the center of the bottom of the inner vessel. There is a cyclic alternation between salt water flowing downwards out of the inner vessel into the outer vessel through the orifice and fresh water flowing upwards into the inner vessel from the outer vessel through that same orifice. We develop a very stable (i.e., stationary) version of this saline oscillator. We first investigate the response of the oscillator to periodic forcing with a train of stimuli (period=Tp) of large amplitude. Each stimulus is the quick injection of a fixed volume of fresh water into the outer vessel followed immediately by withdrawal of that very same volume. For Tp sufficiently close to the intrinsic period of the oscillator (T0) , there is 1:1 synchronization or phase locking between the stimulus train and the oscillator. As Tp is decreased below T0 , one finds the succession of phase-locking rhythms: 1:1, 2:2, 2:1, 2:2, and 1:1. As Tp is increased beyond T0 , one encounters successively 1:1, 1:2, 2:4, 2:3, 2:4, and 1:2 phase-locking rhythms. We next investigate the phase-resetting response, in which injection of a single stimulus transiently changes the period of the oscillation. By systematically changing the phase of the cycle at which the stimulus is delivered (the old phase), we construct the new-phase--old-phase curve (the phase transition curve), from which we then develop a one-dimensional finite-difference equation ("map") that predicts the response to periodic stimulation. These predicted phase-locking rhythms are close to the experimental findings. In addition, iteration of the map predicts the existence of bistability between two different 1:1 rhythms, which was then searched for and found experimentally. Bistability between 1:1 and 2:2 rhythms is also encountered. Finally, with one exception, numerical modeling with a phenomenologically derived Rayleigh oscillator reproduces all of the experimental behavior.

Cardiac electrical dynamics: maximizing dynamical heterogeneity

The relationships between key features of the cardiac electrical activity, such as electrical restitution, discordant alternans, wavebreak, and reentry, and the onset of ventricular tachyarrhythmias have been characterized extensively under the condition of constant rapid pacing. However, it is unlikely that this scenario applies directly to the clinical situation, where the induction of ventricular tachycardia (VT) typically is associated with the interruption of normal cardiac rhythm by several premature beats. To address this issue, we have developed a general theory to explain why specific patterns of premature stimuli increase dynamic heterogeneity of repolarization and precipitate conduction block. The theory predicts that conduction block is caused by (1) creation of a spatial gradient in diastolic interval (DI) by waves traveling at slightly different velocities (ie, conduction velocity dispersion) and (2) amplification of the spatial gradient in DI over subsequent action potentials, secondary to a strong dependence of action potential duration on the preceding DI (ie, a steep action potential duration restitution function). Tests of this theory have been conducted in computer models of homogeneous tissue, where increased spatial dispersion of repolarization during premature stimulation can be attributed solely to the development of dynamical heterogeneity, and in a canine model exhibiting spontaneously occurring VT and sudden death. Our results thus far indicate that the probability of inducing ventricular fibrillation (VF) in the animal model is highest for those sequences predicted to cause conduction block in the computer model. An understanding of the mechanisms underlying these observations will help to identify key electrical phenomena in the onset of VT and fibrillation. Drug and electrical therapies can then be improved by targeting these specific phenomena.

Region of slowed conduction acts as core for spiral wave reentry in cardiac cell monolayers

Pathophysiological heterogeneity in cardiac tissue is related to the occurrence of arrhythmias. Of importance are regions of slowed conduction, which have been implicated in the formation of conduction block and reentry. Experimentally, it has been a challenge to produce local heterogeneity in a manner that is both reversible and well controlled. Consequently, we developed a dual-zone superfusion chamber that can dynamically create a small (5 mm) central island of heterogeneity in cultured cardiac cell monolayers. Three different conditions were studied to explore the effect of regionally slowed conduction on wave propagation and reentry: depolarization by elevated extracellular potassium, sodium channel inhibition with lidocaine, and cell-cell decoupling with palmitoleic acid. Using optical mapping of transmembrane voltage, we found that the central region of slowed conduction always served as the core region around which a spiral wave formed and then revolved following a period of rapid pacing. Because of the localized slowing in the core region, we observed experimentally for the first time an S shape of the spiral wave front near its tip. These results indicate that a small region of slowed conduction can play a crucial role in the formation, anchoring, and modulation of reentrant spiral waves.

Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence

Ventricular fibrillation is a lethal arrhythmia characterized by multiple wavelets usually starting from a single or figure-of-eight re-entrant circuit. Understanding the factors regulating vulnerability to the re-entry is essential for developing effective therapeutic strategies to prevent ventricular fibrillation. In this study, we investigated how pre-existing tissue heterogeneities and electrical restitution properties affect the initiation of re-entry by premature extrastimuli in two-dimensional cardiac tissue models. We studied two pacing protocols for inducing re-entry following the "sinus" rhythm (S1) beat: (1) a single premature (S2) extrastimulus in heterogeneous tissue; (2) two premature extrastimuli (S2 and S3) in homogeneous tissue. In the first case, the vulnerable window of re-entry is determined by the spatial dimension and extent of the heterogeneity, and is also affected by electrical restitution properties and the location of the premature stimulus. The vulnerable window first increases as the action potential duration (APD) difference between the inside and outside of the heterogeneous region increases, but then decreases as this difference increases further. Steeper APD restitution reduces the vulnerable window of re-entry. In the second case, electrical restitution plays an essential role. When APD restitution is flat, no re-entry can be induced. When APD restitution is steep, re-entry can be induced by an S3 over a range of S1S2 intervals, which is also affected by conduction velocity restitution. When APD restitution is even steeper, the vulnerable window is reduced due to collision of the spiral tips.

Rhythms of high-grade block in an ionic model of a strand of regionally ischemic ventricular muscle

Electrical alternans, a beat-to-beat alternation in the electrocardiogram or electrogram, is frequently seen during the first few minutes of acute myocardial ischemia, and is often immediately followed by malignant cardiac arrhythmias such as ventricular tachycardia and ventricular fibrillation. As ischemia progresses, higher-order periodic rhythms (e.g., period-4) can replace the period-2 alternans rhythm. This is also seen in modelling work on a two-dimensional (2-D) sheet of regionally ischemic ventricular muscle. In addition, in the experimental work, ventricular arrhythmias are overwhelmingly seen only after the higher-order rhythms arise. We investigate an ionic model of a strand of ischemic ventricular muscle, constructed as a 3-cm-long 1-D cable with a centrally located 1-cm-long segment exposed to an elevated extracellular potassium concentration ([K(+)](o)). As [K(+)](o) is raised in this "ischemic segment" to represent one major effect of ongoing ischemia, the sequence of rhythms {1:1-->2:2 (alternans)-->2:1} is seen. With further increase in [K(+)](o), one sees higher-order periodic 2N:M rhythms {2:1-->4:2-->4:1-->6:2-->6:1-->8:2-->8:1}. In a 2N:M cycle, only M of the 2N action potentials generated at the proximal end of the cable successfully traverse the ischemic segment, with the remaining ones being blocked within the ischemic segment. Finally, there is a transition to complete block {8:1-->2:0-->1:0} (in an n:0 rhythm, all action potentials die out within the ischemic segment). Changing the length of the ischemic segment results in different rhythms and transitions being seen: e.g., when the ischemic segment is 2 cm long, the period-6 rhythms are not seen; when it is 0.5 cm long, there is a 3:1 rhythm interposed between the 2:1 and 1:0 rhythms. We discuss the relevance of our results to the experimental observations on the higher-order rhythms that presage reentrant ischemic ventricular arrhythmias.

Vulnerability to reentry in a regionally ischemic tissue: a simulation study

Sudden cardiac death is mainly provoked by arrhythmogenic processes. During myocardial ischemia many malignant arrhythmias, such as reentry, take place and can degenerate into ventricular fibrillation. It is thus of great interest to unravel the intricate mechanisms underlying the initiation and maintenance of a reentry. In this computational study, we analyze the probability of reentry during different stages of the acute phase of ischemia. We also aimed at the understanding of the role of its main components: hypoxia, hyperkalemia, and acidosis analyzing the intricate ionic mechanisms responsible for reentry generation. We simulated the electrical activity of a ventricular tissue affected by regional ischemia based on a modified version of the Luo-Rudy model (LRd00). The ischemic conditions were varied to simulate different stages of this pathology. After premature stimulation, we evaluated the vulnerability to reentry. We obtained an unimodal behavior for the vulnerable window as ischemia progressed, peaking at the eighth minute after the onset of ischemia where the vulnerable window yielded 58 ms. Under more severe conditions the vulnerable window decreased and became zero for minute 8.75. The present work provides insight into the mechanisms of reentry generation during ischemia, highlighting the role of acidosis and hypoxia when hyperkalemia is present.

Spatial heterogeneity of the restitution portrait in rabbit epicardium

Spatial heterogeneity of repolarization can provide a substrate for reentry to occur in myocardium. This heterogeneity may result from spatial differences in action potential duration (APD) restitution. The restitution portrait (RP) measures many aspects of rate-dependent restitution: the dynamic restitution curve (RC), S1-S2 RC, and short-term memory response. We used the RP to characterize epicardial patterns of spatial heterogeneity of restitution that were repeatable across animals. New Zealand White rabbit ventricles were paced from the epicardial apex, midventricle, or base, and optical action potentials were recorded from the same three regions. A perturbed downsweep pacing protocol was applied that measured the RP over a range of cycle lengths from 1,000 to 140 ms. The time constant of short-term memory measured close to the stimulus was dependent on location. In the midventricle the mean time constant was 19.1 +/- 1.1 s, but it was 39% longer at the apex (P < 0.01) and 23% longer at the base (P = 0.03). The S1-S2 RC slope was dependent on pacing site (P = 0.015), with steeper slope when pacing from the apex than from the base. There were no significant repeatable spatial patterns in steady-state APD at all cycle lengths or in dynamic RC slope. These results indicate that transient patterns of epicardial heterogeneity of APD may occur after a change in pacing rate. Thus it may affect cardiac electrical stability at the onset of a tachycardia or during a series of ectopic beats. Differences in restitution with respect to pacing site suggest that vulnerability may be affected by the location of reentry or ectopic foci.

Action potential duration dispersion and alternans in simulated heterogeneous cardiac tissue with a structural barrier

Structural barriers to wave propagation in cardiac tissue are associated with a decreased threshold for repolarization alternans both experimentally and clinically. Using computer simulations, we investigated the effects of a structural barrier on the onset of spatially concordant and discordant alternans. We used two-dimensional tissue geometry with heterogeneity in selected potassium conductances to mimic known apex-base gradients. Although we found that the actual onset of alternans was similar with and without the structural barrier, the increase in alternans magnitude with faster pacing was steeper with the barrier--giving the appearance of an earlier alternans onset in its presence. This is consistent with both experimental structural barrier findings and the clinical observation of T-wave alternans occurring at slower pacing rates in patients with structural heart disease. In ionically homogeneous tissue, discordant alternans induced by the presence of the structural barrier arose at intermediate pacing rates due to a source-sink mismatch behind the barrier. In heterogeneous tissue, discordant alternans occurred during fast pacing due to a barrier-induced decoupling of tissue with different restitution properties. Our results demonstrate a causal relationship between the presence of a structural barrier and increased alternans magnitude and action potential duration dispersion, which may contribute to why patients with structural heart disease are at higher risk for ventricular tachyarrhythmias.

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.

Alternating conduction in the ischaemic border zone as precursor of reentrant arrhythmias: a simulation study

AIMS

Here, we investigate the mechanisms underlying the onset of conduction-related arrhythmias in a three-dimensional (3D) computational model of acute regional ischaemia.

METHODS

Ischaemia was introduced by realistic gradients of potassium, pH, oxygen and electrical coupling in a 3D slab of ventricular tissue using the LRd model. We focused on a specific stage (10-15 min after occlusion) at which an intramural non-conductive ischaemic core (IC) surrounded by a border zone (BZ) has formed.

RESULTS

At pacing frequencies greater than 4.5 Hz, we observed narrow areas (0.5 mm wide) of 2:1 conduction blocks at the periphery of the IC. As the pacing frequency increased, the area of block widened to 9 mm and gave rise to reentry at the periphery of the BZ. Alternating conduction blocks produced discordant action potential duration (APD) alternans throughout the slab and T-wave alternans in pseudo-ECG. Slowing the recovery of the calcium current broadened the range of pacing frequencies at which blocks were observed. Hyperkalaemia alone was sufficient to induce the alternating blocks.

CONCLUSION

Computer modelling predicts that ischaemia-related arrhythmias are triggered by calcium-mediated alternating conduction blocks in the ischaemic border zone. Alternating conduction blocks lead to intramural reentry and APD alternans.

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.

Suppression of alternans and conduction blocks despite steep APD restitution: electrotonic, memory, and conduction velocity restitution effects

We examine the utility of the action potential (AP) duration (APD) restitution curve slope in predicting the onset of electrical alternans when electrotonic and memory effects are considered. We develop and use two ionic cell models without memory that have the same restitution curve with slope >1 but different AP shapes and, therefore, different electrotonic effects. We also study a third cell model that incorporates short-term memory of previous cycle lengths, so that it has a family of S1-S2 restitution curves as well as a dynamic restitution curve with slope >1. Our results indicate that both electrotonic and memory effects can suppress alternans, even when the APD restitution curve is steep. In the absence of memory, electrotonic currents related to the shape of the AP, as well as conduction velocity restitution, can affect how alternans develops in tissue and, in some cases, can prevent its induction entirely, even when isolated cells exhibit alternans. When short-term memory is included, alternans may not occur in isolated cells, despite a steep APD restitution curve, and may or may not occur in tissue, depending on conduction velocity restitution. We show for the first time that electrotonic and memory effects can prevent conduction blocks and stabilize reentrant waves in two and three dimensions. Thus we find that the slope of the APD restitution curve alone does not always well predict the onset of alternans and that incorporating electrotonic and memory effects may provide a more useful alternans criterion. A Data Supplement containing movies and JAVA applets is available online at http://ajpheart.physiology.org/cgi/content/full/00747.2003/DC1 .

Effect of acute global ischemia on the upper limit of vulnerability: a simulation study

The goal of this modeling research is to provide mechanistic insight into the effect of altered membrane kinetics associated with 5-12 min of acute global ischemia on the upper limit of cardiac vulnerability (ULV) to electric shocks. We simulate electrical activity in a finite-element bidomain model of a 4-mm-thick slice through the canine ventricles that incorporates realistic geometry and fiber architecture. Global acute ischemia is represented by changes in membrane dynamics due to hyperkalemia, acidosis, and hypoxia. Two stages of acute ischemia are simulated corresponding to 5-7 min (stage 1) and 10-12 min (stage 2) after the onset of ischemia. Monophasic shocks are delivered in normoxia and ischemia over a range of coupling intervals, and their outcomes are examined to determine the highest shock strength that resulted in induction of reentrant arrhythmia. Our results demonstrate that acute ischemia stage 1 results in ULV reduction to 0.8A from its normoxic value of 1.4A. In contrast, no arrhythmia is induced regardless of shock strength in acute ischemia stage 2. An investigation of mechanisms underlying this behavior revealed that decreased postshock refractoriness resulting mainly from 1) ischemic electrophysiological substrate and 2) decrease in the extent of areas positively-polarized by the shock is responsible for the change in ULV during stage 1. In contrast, conduction failure is the main cause for the lack of vulnerability in acute ischemia stage 2. The insight provided by this study furthers our understanding of mechanisms by which acute ischemia-induced changes at the ionic level modulate cardiac vulnerability to electric shocks.

Conduction block in one-dimensional heart fibers

We present a nonlinear dynamical systems analysis of the transition to conduction block in one-dimensional cardiac fibers. We study a simple model of wave propagation in heart tissue that depends only on the recovery of action potential duration and conduction velocity. If the recovery function has slope >or=1 and the velocity recovery function is nonconstant, rapid activation causes dynamical heterogeneity and finally conduction block away from the activation site. This dynamical mechanism may play a role in the initiation and breakup of spiral waves in excitable media.

Triggered alternans in an ionic model of ischemic cardiac ventricular muscle

It has been known for several decades that electrical alternans occurs during myocardial ischemia in both clinical and experimental work. There are a few reports showing that this alternans can be triggered into existence by a premature ventricular contraction. Detriggering of alternans by a premature ventricular contraction, as well as pause-induced triggering and detriggering, have also been reported. We conduct a search for triggered alternans in an ionic model of ischemic ventricular muscle in which alternans has been described recently: a one-dimensional cable of length 3 cm, containing a central ischemic zone 1 cm long, with 1 cm segments of normal (i.e., nonischemic) tissue at each end. We use a modified form of the Luo-Rudy [Circ. Res. 68, 1501-1526 (1991)] ionic model to represent the ventricular tissue, modeling the effect of ischemia by raising the external potassium ion concentration ([K(+)](o)) in the central ischemic zone. As [K(+)](o) is increased at a fixed pacing cycle length of 400 ms, there is first a transition from 1:1 rhythm to alternans or 2:2 rhythm, and then a transition from 2:2 rhythm to 2:1 block. There is a range of [K(+)](o) over which there is coexistence of 1:1 and 2:2 rhythms, so that dropping a stimulus from the periodic drive train during 1:1 rhythm can result in the conversion of 1:1 to 2:2 rhythm. Within the bistable range, the reverse transition from 2:2 to 1:1 rhythm can be produced by injection of a well-timed extrastimulus. Using a stimulation protocol involving delivery of pre- and post-mature stimuli, we derive a one-dimensional map that captures the salient features of the results of the cable simulations, i.e., the {1:1-->2:2-->2:1} transitions with {1:1<-->2:2} bistability. This map uses a new index of the global activity in the cable, the normalized voltage integral. Finally, we put forth a simple piecewise linear map that replicates the {1:1<-->2:2} bistability observed in the cable simulations and in the normalized voltage integral map. (c) 2002 American Institute of Physics.

Resetting and annihilation of reentrant activity in a model of a one-dimensional loop of ventricular tissue

Resetting and annihilation of reentrant activity by a single stimulus pulse (S1) or a pair (S1-S2) of coupled pulses are studied in a model of one-dimensional loop of cardiac tissue using a Beeler-Reuter-type ionic model. Different modes of reentry termination are described. The classical mode of termination by unidirectional block, in which a stimulus produces only a retrograde front that collides with the activation front of the reentry, can be obtained for both S1 and S1-S2 applied over a small vulnerable window. We demonstrate that another scenario of termination-that we term collision block-can also be induced by the S1-S2 protocol. This scenario is obtained over a much wider range of S1-S2 coupling intervals than the one leading to a unidirectional block. In the collision block, S1 produces a retrograde front, colliding with the activation front of the pre-existing reentry, and an antegrade front propagating in the same direction as the initial reentry. Then, S2 also produces an antegrade and a retrograde front. However, the propagation of these fronts in the spatial profile of repolarization left by S1 leads to a termination of the reentrant activity. More complex behaviors also occur in which the antegrade fronts produced by S1 and S2 both persist for several turns, displaying a growing alternation in action potential duration ("alternans amplification") that may lead to the termination of the reentrant activity. The hypothesis that both collision block and alternans amplification depend on the interaction between the action potential duration restitution curve and the recovery curve of conduction velocity is supported by the fact that the dynamical behaviors were reproduced using an integro-delay equation based on these two properties. We thus describe two new mechanisms (collision block and alternans amplification) whereby electrical stimulation can terminate reentrant activity. (c) 2002 American Institute of Physics.

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 reentrant waves, i.e., electrical wave(s) that recirculate 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, have made it extremely difficult to pinpoint the detailed dynamics 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. Movies showing all the breakup mechanisms are available at http://arrhythmia.hofstra.edu/breakup and at ftp://ftp.aip.org/epaps/chaos/E-CHAOEH-12-039203/ INDEX.html. 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. (c) 2002 American Institute of Physics.

A computationally efficient electrophysiological model of human ventricular cells

Recent experimental and theoretical results have stressed the importance of modeling studies of reentrant arrhythmias in cardiac tissue and at the whole heart level. We introduce a six-variable model obtained by a reformulation of the Priebe-Beuckelmann model of a single human ventricular cell. The reformulated model is 4.9 times faster for numerical computations and it is more stable than the original model. It retains the action potential shape at various frequencies, restitution of action potential duration, and restitution of conduction velocity. We were able to reproduce the main properties of epicardial, endocardial, and M cells by modifying selected ionic currents. We performed a simulation study of spiral wave behavior in a two-dimensional sheet of human ventricular tissue and showed that spiral waves have a frequency of 3.3 Hz and a linear core of approximately 50-mm diameter that rotates with an average frequency of 0.62 rad/s. Simulation results agreed with experimental data. In conclusion, the proposed model is suitable for efficient and accurate studies of reentrant phenomena in human ventricular tissue.