Interplay of ionic and structural heterogeneity on functional action potential duration gradients: Implications for arrhythmogenesis

Atrial Fibrillation Mechanisms and Implications for Catheter Ablation

AF is a heterogeneous rhythm disorder that is related to a wide spectrum of etiologies and has broad clinical presentations. Mechanisms underlying AF are complex and remain incompletely understood despite extensive research. They associate interactions between triggers, substrate and modulators including ionic and anatomic remodeling, genetic predisposition and neuro-humoral contributors. The pulmonary veins play a key role in the pathogenesis of AF and their isolation is associated to high rates of AF freedom in patients with paroxysmal AF. However, ablation of persistent AF remains less effective, mainly limited by the difficulty to identify the sources sustaining AF. Many theories were advanced to explain the perpetuation of this form of AF, ranging from a single localized focal and reentrant source to diffuse bi-atrial multiple wavelets. Translating these mechanisms to the clinical practice remains challenging and limited by the spatio-temporal resolution of the mapping techniques. AF is driven by focal or reentrant activities that are initially clustered in a relatively limited atrial surface then disseminate everywhere in both atria. Evidence for structural remodeling, mainly represented by atrial fibrosis suggests that reentrant activities using anatomical substrate are the key mechanism sustaining AF. These reentries can be endocardial, epicardial, and intramural which makes them less accessible for mapping and for ablation. Subsequently, early interventions before irreversible remodeling are of major importance. Circumferential pulmonary vein isolation remains the cornerstone of the treatment of AF, regardless of the AF form and of the AF duration. No ablation strategy consistently demonstrated superiority to pulmonary vein isolation in preventing long term recurrences of atrial arrhythmias. Further research that allows accurate identification of the mechanisms underlying AF and efficient ablation should improve the results of PsAF ablation.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model

Accumulation of intracellular Na⁺ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na⁺ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na⁺ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na⁺]_i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na⁺]_i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na⁺]_i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na⁺]_i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na⁺]_i may play complex roles in cellular and tissue-level cardiac dynamics.

Simultaneous conduction mapping and intracellular membrane potential recording in isolated atria

We describe a novel approach for simultaneously determining regional differences in action potential (AP) morphology and tissue electrophysiological properties in isolated atria. The epicardial surface of rat atrial preparations was placed in contact with a multi-electrode array (9 × 10 silver chloride electrodes, 0.1 mm diameter and 0.1 mm pitch). A glass microelectrode (100 MΩ) was simultaneously inserted into the endocardial surface to record intracellular AP from either of 2 regions (A, B) during pacing from 2 opposite corners of the tissue. AP duration at 80% of repolarisation and its restitution curve was significantly different only in region A (p < 0.01) when AP was initiated at different stimulation sites. Alternans in AP duration and AP amplitude, and in conduction velocity were observed during 2 separate arrhythmic episodes. This approach of combining microelectrode array and intracellular membrane potential recording may provide new insights into arrhythmogenic mechanisms in animal models of cardiovascular disease.

Attraction of rotors to the pulmonary veins in paroxysmal atrial fibrillation: a modeling study

Maintenance of paroxysmal atrial fibrillation (AF) by fast rotors in the left atrium (LA) or at the pulmonary veins (PVs) is not fully understood. To gain insight into this dynamic and complex process, we studied the role of the heterogeneous distribution of transmembrane currents in the PVs and LA junction (PV-LAJ) in the localization of rotors in the PVs. We also investigated whether simple pacing protocols could be used to predict rotor drift in the PV-LAJ. Experimentally observed heterogeneities in IK1, IKs, IKr, Ito, and ICaL in the PV-LAJ were incorporated into two- and pseudo three-dimensional models of Courtemanche-Ramirez-Nattel-Kneller human atrial kinetics to simulate various conditions and investigate rotor drifting mechanisms. Spatial gradients in the currents resulted in shorter action potential duration, minimum diastolic potential that was less negative, and slower upstroke and conduction velocity for rotors in the PV region than in the LA. Rotors under such conditions drifted toward the PV and stabilized at the shortest action potential duration and less-excitable region, consistent with drift direction under intercellular coupling heterogeneities and regardless of the geometrical constraint in the PVs. Simulations with various IK1 gradient conditions and current-voltage relationships substantiated its major role in the rotor drift. In our 1:1 pacing protocol, we found that among various action potential properties, only the minimum diastolic potential gradient was a rate-independent predictor of rotor drift direction. Consistent with experimental and clinical AF studies, simulations in an electrophysiologically heterogeneous model of the PV-LAJ showed rotor attraction toward the PV. Our simulations suggest that IK1 heterogeneity is dominant compared to other currents in determining the drift direction through its impact on the excitability gradient. These results provide a believed novel framework for understanding the complex dynamics of rotors in AF.

Interactions between cardiac fibrosis spatial pattern and ionic remodeling on electrical wave propagation

Cardiac fibrosis is an important form of pathological tissue remodeling. Fibrosis can electrically-uncouple neighboring excitable cardiomyocytes thus acting as an obstacle to electrical propagation. In this study, we investigated the effects of fibrosis spatial pattern on electrical propagation in control, decreased maximum sodium conductance, and increased intracellular resistivity conditions. Simulations were performed with a monodomain approach and a realistic canine ionic model. We found that the propagation failure is highly dependent on the spatial pattern of fibrosis for all conditions studied with maximum sensitivity for patterns with combination of small and large clusters. However, the effect is particularly sensitive to reduced sodium current condition where conduction block occurred at lower fibrosis density.

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

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

A computer model of engineered cardiac monolayers

Engineered monolayers created using microabrasion and micropatterning methods have provided a simplified in vitro system to study the effects of anisotropy and fiber direction on electrical propagation. Interpreting the behavior in these culture systems has often been performed using classical computer models with continuous properties. However, such models do not account for the effects of random cell shapes, cell orientations, and cleft spaces inherent in these monolayers on the resulting wavefront conduction. This work presents a novel methodology for modeling a monolayer of cardiac tissue in which the factors governing cell shape, cell-to-cell coupling, and degree of cleft space are not constant but rather are treated as spatially random with assigned distributions. This modeling approach makes it possible to simulate wavefront propagation in a manner analogous to performing experiments on engineered monolayer tissues. Simulated results are compared to previously published measured data from monolayers used to investigate the role of cellular architecture on conduction velocities and anisotropy ratios. We also present an estimate for obtaining the electrical properties from these networks and demonstrate how variations in the discrete cellular architecture affect the macroscopic conductivities. The simulations support the common assumption that under normal ranges of coupling strength, tissues with relatively uniform distributions of cell shapes and connectivity can be represented using continuous models with conductivities derived from random discrete cellular architecture using either global or local estimates. The results also reveal that in the presence of abrupt changes in cell orientation, local estimates of tissue properties predict smoother changes in conductivity that may not adequately predict the discrete nature of propagation at the transition sites.

[Multiscale modeling of cardiac electrical activity]

Models of cardiac electrical activity cover a wide range of spatial scales, from the genesis of the ionic currents in individual cardiomyocytes to the generation of electrocardiograms on the torso. The level of detail that is appropriate and practicable depends on the problem investigated and the scope of the computations that are required. We briefly present three examples of modelling: the dynamics of the entrainment of a single cell, the impact of fibrosis on electrical propagation in a piece of tissue and the generation of ECG in human. In each case, the methods, results and limitations are discussed.

Pacing-induced spatiotemporal dynamics can be exploited to improve reentry termination efficacy

Some potentially fatal cardiac arrhythmias may be terminated by a series of premature stimuli. Monomorphic ventricular tachycardia, which may be modeled as an excitation wave traveling around in a ring, is one such arrhythmia. We investigated the mechanisms and requirements for termination of such reentry using an ionic cardiac ring model. Termination requires conduction block, which in turn is facilitated by spatial dispersion in repolarization and recovery time. When applying short series of two or three stimuli, we found that for conduction block to robustly occur, the magnitude of the spatial gradient in recovery time must exceed a critical value of 20 ms/cm. Importantly, the required spatial gradient can be induced in this homogeneous system by the dynamics of the stimulus-induced waves-we show analytically the necessary conditions. Finally, we introduce a type of pacing protocol, the "aggressive ramp," which increases the termination efficacy by exploiting such pacing-induced heterogeneities. This technique, which is straightforward to implement, may therefore have important clinical implications.

Mechanisms of transition from normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of tissue heterogeneity and anisotropy

Experimental evidence suggests that regional differences in action potential (AP) morphology can provide a substrate for initiation and maintenance of reentrant arrhythmias in the right atrium (RA), but the relationships between the complex electrophysiological and anatomical organization of the RA and the genesis of reentry are unclear. In this study, a biophysically detailed three-dimensional computer model of the right atrial tissue was constructed to study the role of tissue heterogeneity and anisotropy in arrhythmogenesis. The model of Lindblad et al. for a rabbit atrial cell was modified to incorporate experimental data on regional differences in several ionic currents (primarily, I(Na), I(CaL), I(K1), I(to), and I(sus)) between the crista terminalis and pectinate muscle cells. The modified model was validated by its ability to reproduce the AP properties measured experimentally. The anatomical model of the rabbit RA (including tissue geometry and fiber orientation) was based on a recent histological reconstruction. Simulations with the resultant electrophysiologically and anatomically detailed three-dimensional model show that complex organization of the RA tissue causes breakdown of regular AP conduction patterns at high pacing rates (>11.75 Hz): as the AP in the crista terminalis cells is longer, and electrotonic coupling transverse to fibers of the crista terminalis is weak, high-frequency pacing at the border between the crista terminalis and pectinate muscles results in a unidirectional conduction block toward the crista terminalis and generation of reentry. Contributions of the tissue heterogeneity and anisotropy to reentry initiation mechanisms are quantified by measuring action potential duration (APD) gradients at the border between the crista terminalis and pectinate muscles: the APD gradients are high in areas where both heterogeneity and anisotropy are high, such that intrinsic APD differences are not diminished by electrotonic interactions. Thus, our detailed computer model reconstructs complex electrical activity in the RA, and provides new insights into the mechanisms of transition from focal atrial tachycardia into reentry.

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.

Model of reentrant ventricular tachycardia based on infarct border zone geometry predicts reentrant circuit features as determined by activation mapping

BACKGROUND

Infarct border zone (IBZ) geometry likely affects inducibility and characteristics of postinfarction reentrant ventricular tachycardia, but the connection has not been established.

OBJECTIVE

The purpose of this study was to determine characteristics of postinfarction ventricular tachycardia in the IBZ.

METHODS

A geometric model describing the relationship between IBZ geometry and wavefront propagation in reentrant circuits was developed. Based on the formulation, slow conduction and block were expected to coincide with areas where IBZ thickness (T) is minimal and the local spatial gradient in thickness (DeltaT) is maximal, so that the degree of wavefront curvature rho proportional, variant DeltaT/T is maximal. Regions of fastest conduction velocity were predicted to coincide with areas of minimum DeltaT. In seven arrhythmogenic postinfarction canine heart experiments, tachycardia was induced by programmed stimulation, and activation maps were constructed from multichannel recordings. IBZ thickness was measured in excised hearts from histologic analysis or magnetic resonance imaging. Reentrant circuit properties were predicted from IBZ geometry and compared with ventricular activation maps after tachycardia induction.

RESULTS

Mean IBZ thickness was 231 +/- 140 microm at the reentry isthmus and 1440 +/- 770 microm in the outer pathway (P <0.001). Mean curvature rho was 1.63 +/- 0.45 mm(-1) at functional block line locations, 0.71 +/- 0.18 mm(-1) at isthmus entrance-exit points, and 0.33 +/- 0.13 mm(-1) in the outer reentrant circuit pathway. The mean conduction velocity about the circuit during reentrant tachycardia was 0.32 +/- 0.04 mm/ms at entrance-exit points, 0.42 +/- 0.13 mm/ms for the entire outer pathway, and 0.64 +/- 0.16 mm/ms at outer pathway regions with minimum DeltaT. Model sensitivity and specificity to detect isthmus location was 75.0% and 97.2%.

CONCLUSIONS

Reentrant circuit features as determined by activation mapping can be predicted on the basis of IBZ geometrical relationships.

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.

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.

Vulnerable window for conduction block in a one-dimensional cable of cardiac cells, 2: multiple extrasystoles

Unidirectional conduction block of premature extrasystoles can lead to initiation of cardiac reentry, causing lethal arrhythmias including ventricular fibrillation. Multiple extrasystoles are often more effective at inducing unidirectional conduction block and reentry than a single extrasystole. Since the substrate for conduction block is spatial dispersion of refractoriness, in this study we investigate how the first extrasystole modulates this dispersion to influence the "vulnerable window" for conduction block by subsequent extrasystoles, particularly in relation to action potential duration restitution and conduction velocity restitution properties. Using a kinematic model to represent wavefront-waveback interactions and simulations with the Luo-Rudy model in a one-dimensional cable of cardiac cells, we show that in homogeneous tissue, a premature extrasystole can create a large dispersion of refractoriness leading to conduction block of a subsequent extrasystole. In heterogeneous tissue, however, a premature extrasystole can either reduce or enhance the dispersion of refractoriness depending on its propagation direction with respect to the previous beat. With multiple extrasystoles at random coupling intervals, vulnerability to conduction block is proportional to their number. In general, steep action potential duration restitution and broad conduction velocity restitution promote dispersion of refractoriness in response to multiple extrasystoles, and thus enhance vulnerability to conduction block. These restitution properties also promote spatially discordant alternans, a setting which is particularly prone to conduction block. The equivalent dispersion of refractoriness created dynamically in homogeneous tissue by spatially discordant alternans is more likely to cause conduction block than a comparable degree of preexisting dispersion in heterogeneous tissue.

Vulnerable window for conduction block in a one-dimensional cable of cardiac cells, 1: single extrasystoles

Spatial dispersion of refractoriness, which is amplified by genetic diseases, drugs, and electrical and structural remodeling during heart disease, is recognized as a major factor increasing the risk of lethal arrhythmias and sudden cardiac death. Dispersion forms the substrate for unidirectional conduction block, which is required for the initiation of reentry by extrasystoles or rapid pacing. In this study, we examine theoretically and numerically how preexisting gradients in refractoriness control the vulnerable window for unidirectional conduction block by a single premature extrasystole. Using a kinematic model to represent wavefront-waveback interactions, we first analytically derived the relationship (under simplified conditions) between the vulnerable window and various electrophysiological parameters such as action potential duration gradients, refractoriness barriers, conduction velocity restitution, etc. We then compared these findings to numerical simulations using the kinematic model or the Luo-Rudy action potential model in a one-dimensional cable of cardiac cells. The results from all three methods agreed well. We show that a critical gradient in action potential duration for conduction block can be analytically derived, and once this critical gradient is exceeded, the vulnerable window increases proportionately with the refractory barrier and is modulated by conduction velocity restitution and gap junctional conductance. Moreover, the critical gradient for conduction block is higher for an extrasystole traveling in the opposite direction from the sinus beat than for one traveling in the same direction (e.g., an epicardial extrasystole versus an endocardial extrasystole).

Reentry wave formation in excitable media with stochastically generated inhomogeneities

Clinical research shows that the frequency of arrhythmia events depends on the number and area of the border zones of infarct scars. We investigate the possibility that arrhythmia is initiated by reentry waves generated by the inhomogeneity of conduction velocity at the border zone. The interaction of a plane wave with a spatially extended inhomogeneity is simulated in the FitzHugh- Nagumo model. The inhomogeneity is introduced into the model by modifying the spatial dependence of the diffusion coefficient in a stochastic manner. This results in a rich variety of spatial distributions of conductivity. A plane wave propagating through such a system may break up on the regions with low conductivity and produce numerous spiral waves. The frequency of reentry wave formation is studied as a function of the parameters of the inhomogeneity generation algorithm. Three main scenarios of reentry wave formation were found: unidirectional block, main wave-wavelet collision, and wave break up during collision, on a region in which a conduction velocity gradient occurs. These scenarios are likely candidates for the mechanisms of arrhythmia initiation in a damaged tissue, e.g., the border zone of an infarct scar.

Electrotonic influences on action potential duration dispersion in small hearts: a simulation study

Intrinsic spatial variations in repolarization currents in the heart can produce spatial gradients in action potential duration (APD) that serve as possible sites for conduction block and the initiation of reentrant activity. In well-coupled myocardium, however, electrotonic influences at the stimulus site and wavefront collision sites act to modulate any intrinsic heterogeneity in APD. These effects alter APD gradients over an extent larger than that suggested by the length constant associated with propagation and, thus, are hypothesized to play a greater role in smaller hearts used as experimental models of human disease. This study uses computer simulation to investigate how heart size, tissue properties, and the spatial assignment of cell types affect functional APD dispersion. Simulations were carried out using the murine ventricular myocyte model of Pandit et al. or the Luo-Rudy mammalian model in three-dimensional models of mouse and rabbit ventricular geometries. Results show that the spatial extent of the APD dispersion is related to the dynamic changes in transmembrane resistance during recovery. Also, because of the small dimensions of the mouse heart, electrotonic effects on APD primarily determine the functional dispersion of refractoriness, even in the presence of large intrinsic cellular heterogeneity and reduced coupling. APD dispersion, however, is found to increase significantly when the heart size increases to the size of a rabbit heart, unmasking intrinsic cell types.

Dispersion of cardiac action potential duration and the initiation of re-entry: a computational study

BACKGROUND

The initiation of re-entrant cardiac arrhythmias is associated with increased dispersion of repolarisation, but the details are difficult to investigate either experimentally or clinically. We used a computational model of cardiac tissue to study systematically the association between action potential duration (APD) dispersion and susceptibility to re-entry.

METHODS

We simulated a 60 x 60 mm2 D sheet of cardiac ventricular tissue using the Luo-Rudy phase 1 model, with maximal conductance of the K+ channel gKmax set to 0.004 mS mm(-2). Within the central 40 x 40 mm region we introduced square regions with prolonged APD by reducing gKmax to between 0.001 and 0.003 mS mm(-2). We varied (i) the spatial scale of these regions, (ii) the magnitude of gKmax in these regions, and (iii) cell-to-cell coupling.

RESULTS

Changing spatial scale from 5 to 20 mm increased APD dispersion from 49 to 102 ms, and the susceptible window from 31 to 86 ms. Decreasing gKmax in regions with prolonged APD from 0.003 to 0.001 mS mm-2 increased APD dispersion from 22 to 70 ms, and the susceptible window from <1 to 56 ms. Decreasing cell-to-cell coupling by changing the diffusion coefficient from 0.2 to 0.05 mm2 ms(-1) increased APD dispersion from 57 to 88 ms, and increased the susceptible window from 41 to 74 ms.

CONCLUSION

We found a close association between increased APD dispersion and susceptibility to re-entrant arrhythmias, when APD dispersion is increased by larger spatial scale of heterogeneity, greater electrophysiological heterogeneity, and weaker cell-to-cell coupling.

Introduction: Mapping and control of complex cardiac arrhythmias

This paper serves as an introduction to the Focus Issue on mapping and control of complex cardiac arrhythmias. We first introduce basic concepts of cardiac electrophysiology and describe the main clinical methods being used to treat arrhythmia. We then provide a brief summary of the main themes contained in the articles in this Focus Issue. In recent years there have been important advances in the ability to map the spread of excitation in intact hearts and in laboratory settings. This work has been combined with simulations that use increasingly realistic geometry and physiology. Waves of excitation and contraction in the heart do not always propagate with constant velocity but are often subject to instabilities that may lead to fluctuations in velocity and cycle time. Such instabilities are often treated best in the context of simple one- or two-dimensional geometries. An understanding of the mechanisms of propagation and wave stability is leading to the implementation of different stimulation protocols in an effort to modify or eliminate abnormal rhythms. (c) 2002 American Institute of Physics.