AJP - Heart Journal of Applied Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 289: H1002-H1012, 2005. First published April 22, 2005; doi:10.1152/ajpheart.00252.2005
0363-6135/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Movie
Right arrow All Versions of this Article:
289/3/H1002    most recent
00252.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (9)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zou, R.
Right arrow Articles by Nattel, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zou, R.
Right arrow Articles by Nattel, S.

Substrate size as a determinant of fibrillatory activity maintenance in a mathematical model of canine atrium

Renqiang Zou,1 James Kneller,1,2 L. Joshua Leon,3 and Stanley Nattel1,2

1Department of Medicine and Research Center, University of Montreal and Montreal Heart Institute, and 2Department of Pharmacology and Therapeutics, McGill University, Montreal, Quebec; and 3Department of Electrical Engineering, University of Calgary, Calgary, Alberta, Canada

Submitted 15 March 2005 ; accepted in final form 15 April 2005


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Tissue size has been considered an important determinant of atrial fibrillation (AF), but recent work has questioned the critical size hypothesis. Here, we use a previously developed mathematical model of the two-dimensional canine atrium with realistic action potential, ionic, and conduction properties to address substrate size effects on the maintenance of fibrillatory activity. Cholinergic AF was simulated at different acetylcholine (ACh) concentrations ([ACh]) and distributions, with substrate area varied 11.1-fold. Automated phase singularity detection was used to facilitate the analysis of arrhythmic activity. The duration of activity induced by a single extrastimulus increased with increasing substrate dimensions. Two general mechanisms underlying activity were observed and were differentially affected by substrate size. For large mean [ACh], single primary rotors anchored in low-[ACh] zones maintained activity and substrate dimensions were not critical. At lower mean [ACh], extensive spiral wave meander prevented the emergence of single stable rotors. Prolonged activity was favored when substrate size permitted a sufficiently large number of simultaneous longer-lasting rotors that extinction of all was unlikely. Thus either single dominant rotor or multiple reentrant spiral generator mechanisms could maintain fibrillatory activity in this model and were differentially dependent on substrate size. These results speak to recent debates about the role in AF of single driver rotors versus multiple reentrant circuit mechanisms by suggesting that either may maintain fibrillatory atrial activity depending on atrial size and electrophysiological properties.

heart arrhythmia mechanisms; phase singularity; spiral waves; bioelectricity computer modeling; cardiac action potential


ATRIAL FIBRILLATION (AF) is the most common sustained arrhythmia in humans, is an important contributor to human morbidity and possibly mortality, and presents a variety of therapeutic challenges (29). The mechanisms underlying AF have been a subject of long-standing debate and recent intense investigation (31). Some of the basic properties of fibrillation induced by faradic stimulation were described by MacWilliam (23) in the late 1800s, along with the discoordinating effects of vagal nerve activation on atrial fluttering caused by faradic currents. MacWilliam also showed that ventricular fibrillation could persist in isolated sections of ventricular tissue. Garrey (8) extended MacWilliam's approach and noted that the ability of isolated tissue to fibrillate depends on its mass, introducing the notion of a "critical mass" for fibrillation related to conditions required for multiple circuit reentry. This idea was later incorporated in Moe's "multiple wavelet hypothesis" (26). An alternative conceptual framework championed by Lewis (20) was the idea of a primary reentrant circuit (a "mother wave"), with fibrillatory activity due to an irregular response of the atrial substrate.

Experimental support has been presented for both of these paradigms. AF persistence across species varies with their cardiac size, with AF being much longer lasting in larger hearts (28). On the other hand, the ability of the mouse heart to sustain fibrillation has been presented as a challenge to the critical mass hypothesis (41), and an analysis of dominant frequencies and phase singularity (PS) properties during AF in isolated sheep hearts has led to the suggestion that one or a small number of reentrant generators rather than multiple wavelet reentry underlies AF (4).

Testing the "critical mass" notion in biological tissues is difficult because paradigms that allow for the examination of a range of tissue sizes (e.g., different species, models producing atrial dilation, developmental changes, etc.) are also associated with discrepancies in atrial electrophysiological properties. Mathematical models of atrial activity do not share this limitation: the size of the computational substrate can be varied while the underlying electrophysiological properties are kept constant.

We (18) have previously analyzed the properties of AF in a mathematical model of cholinergic AF with a two-dimensional substrate having realistic cellular and tissue electrical properties for the dog. In a 5 x 10-cm substrate (taken to represent normal canine atrial dimensions), AF is generally maintained by a single dominant spiral wave with wavebreak in a heterogeneous substrate (18). The role of substrate dimensions was not addressed, and we are unaware of other studies that assessed the role of substrate dimensions specifically. In the present investigation, we analyzed the relationship between atrial substrate size and AF persistence, along with associated mechanisms, in the model. We evaluated activity under 12 discrete conditions of acetylcholine (ACh) concentration ([ACh]) and distribution in each of five rectangular atrial substrates over an 11.1-fold range of surface areas. Induced activity was followed for as long as it was maintained, to a maximum of 20 s.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Model description and implementation. We used the Ramirez-Nattel-Courtemanche (RNC) canine atrial action potential (AP) model, as previously described (17, 34). The cell membrane was modeled as an electrical capacitor-resistor system with ionic current (Iion) given by

(1)
where Iion includes contributions of Na+ current (INa); inward rectifier K+ current (IK1); transient outward current (Ito); ultrarapid (IKur,d), rapid (IKr), and slow (IKs) delayed rectifier K+ currents; L-type Ca2+ current (ICa); Ca2+-activated Cl current (IClCa); ACh-dependent K+ current (IKACh); sarcolemmal Ca2+ pump current (IpCa); Na+/Ca2+ exchanger current (INaCa); Na+-K+-ATPase current (INaK); and background Na+ (Ib,Na), Ca2+ (Ib,Ca), and Cl currents (Ib,Cl). APs were calculated by numerical integration of the relation

(2)
where V is voltage, t is time, Istim is the stimulus current, and Cm is membrane capacitance.

A two-dimensional rectangular atrial substrate was modeled as previously described (16, 18). The substrate consists of discreet cables representing muscle fibers, with each composed of a chain of myocytes that form a longitudinal syncytium. Fixed transverse resistors (10 M{Omega}) represent lateral coupling. Fiber resistivity and interfiber resistances were chosen to match experimental propagation characteristics. The discretization factor (segment length per length constant) was kept <0.1 to prevent numerical distortions (15, 36). Voltage-dependent expressions involving exponentials were precomputed with 50-µV resolution and stored in a table to increase computational efficiency (42). Activity was analyzed in substrates of 3 x 6, 4 x 8, 5 x 10, 7.5 x 15, and 10 x 20 cm, including two smaller and two larger substrates than in our previous studies (16, 18). To change tissue size, we varied both the number of cables (180, 240, 300, 450, and 600 cables, respectively) and the number of cells per cable (360, 480, 600, 900, and 1,200, correspondingly). The atrial conduction velocity (CV) anisotropy ratio was 2.4:1, and longitudinal CV was 80 cm/s, compatible with experimental measurements of CV in the canine atrium (21).

Two-dimensional simulations. Reentry was initiated with a cross-shock protocol (19). A basic (S1) stimulus was applied along the uppermost cable, creating a transversely propagating wavefront. Before the recovery front reached halfway, stimulus S2 was applied as used experimentally (46), creating a PS and initiating reentry. Stimuli were 1 ms and 180 µA/cm2. A previously described IKACh formulation was used (18). [ACh] was varied across the tissue as previously between a minimum and maximum according to a sine wave function (18). Four [ACh] distributions (1.67, 2.0, 2.4, and 3.3 cm peak to trough) and three mean [ACh] densities (1.875, 3.75, and 7.5 nM) were studied, based on values previously found to permit AF perpetuation (18). Concentration fluctuations were between 0 and twice the mean [ACh], such that for a given [ACh] distribution, larger mean [ACh] implied larger differences in maximum and minimum [ACh] and AP duration (APD) values. Note that the lowest [ACh] was 0 nM for all conditions, so that varying mean [ACh] corresponded to varying the maximum concentration and the [ACh] gradient, with minimum values fixed. To assess APD and APD distribution resulting from different [ACh] conditions, the 5 x 10-cm substrate was paced at 6 Hz. APD to 75% of full repolarization (APD75) of the 10th cycle was then determined at all sites in the grid with custom-written software. Simulations (20 s) were run for each concentration and spatial distribution in each substrate size, providing 12 conditions/size.

For display and PS analysis purposes, the transmembrane potential (TMP) was subsampled at the center of each 6 x 6 cell square on the computational grid after each millisecond, which we showed caused no loss of resolution. The model was implemented in C language on a Sun Enterprise 10000 computer running on 16 parallel processors. Data was output every millisecond.

Data analysis. To analyze AF quantitatively, we used an automated algorithm for PS identification and tracking (51). The algorithm proved essential for analysis of complex activation in larger grids. PS candidates were located with Bray's convolution method (3) and an image-analysis approach. We identified PSs that satisfied both criteria, along with tracking considerations, followed by manual confirmation to exclude artifacts. To facilitate analysis of activation, we simultaneously displayed animated PS trajectories and TMP maps. Illustrative examples are provided in on-line movies as indicated (see http://ajpheart.physiology.org/cgi/content/full/00252.2005/DC1). Mechanisms underlying fibrillatory activity were classified as 1) single generator (SG), when one dominant rotor with consistent position (core tip remaining in a single low-[ACh] zone)-maintained activity; 2) a series of single generators (SSG), when activity was maintained by single dominant rotors but individual rotors eventually died out and were replaced by rotors arising in other positions; and 3) multiple generator zones (MGZ), when activity resulted from multiple simultaneous rotors of varying location and duration. To explore the frequency content of simulated AF, spectral analysis was performed with fast Fourier transformations (FFTs), with activity sampled at 1 kHz for 20,000 frames (20 s), providing a maximal resolution of 0.05 Hz. Algorithms were implemented with Matlab version 6.5 for Windows on a Pentium 4 personal computer.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Substrate dimensions and AF maintenance. We ran simulations for 60 size/[ACh]/[ACh] distribution conditions. For each condition, we first defined a S1S2 interval that induced at least five cycles of spontaneous activity and then analyzed the activity induced. This was possible for all but the smallest substrate (3 x 6 cm), in which only one [ACh]/[ACh] distribution condition allowed for the induction of repetitive activity. Substrate size was clearly a strong determinant of AF duration (Fig. 1A). The roles of [ACh] and [ACh] distribution are illustrated in Fig. 1, B and C, respectively. The results for each mean [ACh] (Fig. 1B) are the average for all [ACh] distributions at that concentration, whereas the results for each [ACh] distribution (Fig. 1C) are the average for all [ACh]s at that distribution. For any mean [ACh] or [ACh] distribution, AF duration increased progressively as substrate dimension increased. Whereas mean [ACh] significantly determined AF duration, the spatial distribution of [ACh] was not a systematic determinant of AF duration. To verify that arrhythmic behavior was independent of the initial S1S2 conditions, we varied the S1S2 interval in the average-sized substrate (5 x 10 cm) by 5-ms increments and observed up to 5 s of activity initiated by three different S1S2 intervals for each ACh concentration/distribution condition. The results (Table 1) indicate that AF duration was relatively consistent for each ACh concentration/distribution condition and showed no systematic dependence on initial conditions.



View larger version (32K):
[in this window]
[in a new window]
 
Fig. 1. A: relationship between atrial fibrillation (AF) duration in each simulation and substrate size, along with means ± SD. B: AF duration in each simulation as a function of substrate dimensions and mean [ACh]. C: AF duration in each simulation as a function of substrate dimensions and [ACh] distribution (distance between peak and trough [ACh]). (Results are not shown for the 3 x 6-cm substrate because AF was only inducible in 1 of 12 [ACh]/[ACh] distribution conditions for this substrate size.)

 

View this table:
[in this window]
[in a new window]
 
Table 1. Relation of AF duration to starting (S1S2) conditions and ACh concentration/distribution for average-sized (5 x 10 cm) computational substrate

 
Mechanisms of perpetuation and termination of fibrillatory activity. Figure 1 indicates that substrate size is indeed a determinant of the maintenance of fibrillatory activity. To evaluate potential underlying mechanisms, we studied the activation pattern during each simulation.

In our previous analysis of simulated AF in a 5 x 10-cm substrate (18), most cases of longer-lasting AF were maintained by single dominant rotors. In the present study, activity triggered by S2 was always initially maintained by a single clockwise rotor with its core tip located in a zone of low [ACh]. The fate of this rotor depended largely on [ACh]. For the highest mean [ACh] (corresponding to the largest [ACh] differences between peak and trough [ACh]), the initial rotor tended to remain strongly anchored in the same low-[ACh] zone. Figure 2 shows an example in a 4 x 8-cm grid. Figure 2A is an atrial TMP map at a single time point during sustained fibrillatory activity. A star indicates the position of the dominant rotor core tip, the dashed arrow shows the direction of rotation, and solid arrows show the directions of wavefront propagation. An animation of TMP and PS maps is shown in on-line movie 1, with an evolving TMP map at the top and corresponding PS trajectories at the bottom (see http://ajpheart.physiology.org/cgi/content/full/00252.2005/DC1). PS trajectories are color coded, with PSs lasting >1 s in green, PSs lasting >1 cycle but <1 s in blue, and PSs lasting <1 cycle in red. Instantaneous PS points are shown as bright color-coded dots, with the preceding 50 ms of each trajectory as a thick line and the complete previous PS trajectory history as thin lines. Running simulation time is shown in digital format. For the example illustrated, PS meander occurred but was never sufficient for the primary rotor core tip to escape its initial low-[ACh] zone location. Wavebreak occurred when the emanating impulse impinged on other low-[ACh] zones, generally producing two short-lived counterrotating daughter wavelets that mutually annihilated or drifted off the substrate surface. Mutual annihilation of counterrotating secondary waves resulting from wavebreak was strongly favored because of the sharp rotation of such waves, which restricted them to the same low-[ACh] zone.



View larger version (47K):
[in this window]
[in a new window]
 
Fig. 2. Atrial activity in a substrate of size 4 x 8 cm with mean [ACh] of 7.5 nM and [ACh] distribution of 1.67 cm. A: atrial transmembrane potential (TMP) map at 2,551 ms of the simulation (voltage scale at right). The large white star indicates the primary generator core tip phase singularity (PS) location, with the dashed arrow indicating the direction of rotation. Solid arrows indicate other wavefront propagation directions. B: PS trajectory map superimposed on the [ACh] distribution grid ([ACh] scale to right). C: PS lifespan histogram. D: action potentials (APs) from the core of the primary spiral (position shown by * in B). E: fast Fourier transform (FFT) of activity in D. F: APs from the periphery of the primary rotor (position shown by # in B). G: FFT of activity in F. (See also on-line movie 1.)

 
Figure 2B shows all PS trajectories superimposed on an [ACh] map. The location of the primary rotor tip in a single low-[ACh] zone is evident, as are daughter wavelets, which self-annihilated by colliding with each other. The PS lifespan histogram (Fig. 2C) shows a large number of short-lived PSs and one PS at the primary rotor wave tip that remained stable throughout the simulation. APs in the low-[ACh] core region (Fig. 2D) show great amplitude variability and a lower dominant frequency (Fig. 2E) than in the adjacent high-[ACh] zone (Fig. 2, F and G). Both display the irregular periodicity of fibrillatory activity.

Single rotor-based AF, as shown in Fig. 2, occurred typically at the highest mean [ACh] (7.5 nM). Of 17 simulations at 7.5 nM mean [ACh], AF was sustained by single rotors for the full simulation in 14 simulations (82%; Table 2). In 10 of these, a single stable rotor maintained AF throughout the simulation (SG-type activity, as in Fig. 2), whereas in 4 the initial rotor eventually terminated but was immediately replaced by another long-lasting rotor generated by a daughter wave. Fibrillatory activity was then maintained by SSG. Figure 3 and on-line movie 2 (see http://ajpheart.physiology.org/cgi/content/full/00252.2005/DC1) show an example in the 10 x 20-cm substrate. Figure 3A shows APs in the low-[ACh] core of the predominant rotor. Activation is clearly irregular. The PS trajectory map (Fig. 3B) shows the principal generator rotor in the upper center as well as two transient rotor zones lasting >1 s (green trajectories) in other low-[ACh] regions. Figure 3C shows TMP maps at four time points corresponding to the extinction and replacement of a primary rotor. The clockwise rotor is shown at 8,401 ms, with its core tip PS at the star. At 8,891 ms, two counterrotating daughter wavelets have formed with core tip PSs (small stars) in the same low-[ACh] zone. The primary rotor collided with one of these wavelets (dashed red line), causing mutual annihilation. However, the remaining daughter wave formed a stable rotor that took over primary generator function, as shown at 9,076 and 9,401 ms. The sequence of events can be visualized in corresponding running time in on-line movie 2. SSG activity tended to occur with larger [ACh] distributions, which allowed for the transient formation of multiple wavefronts with core tips within individual (relatively large) low-[ACh] zones.


View this table:
[in this window]
[in a new window]
 
Table 2. AF mechanism and duration in relation to substrate size, [ACh], and [ACh] distribution

 


View larger version (56K):
[in this window]
[in a new window]
 
Fig. 3. Atrial activity in a substrate of size 10 x 20 cm with mean [ACh] of 7.5 nM and [ACh] distribution of 3.3 cm. A: atrial APs from the core region of the primary generator rotor. B: PS trajectory map superimposed on the [ACh] distribution grid ([ACh] scale to right). C: TMP maps at times indicated during the simulation. The large white star indicates PS at the core tip of the initial primary rotor. The dashed red line at 8,891 ms indicates annihilation between initial primary rotor and a daughter wave. The white star with the internal red asterisk indicates a generator rotor arising from a daughter wave at 8,891 ms. Dashed arrows show the direction of primary rotor rotation. Other arrows indicate directions of other wavefronts. For a discussion, see text (see also on-line movie 2).

 
In three cases, daughter waves succeeded in creating independent secondary rotors. For one of these, single rotors eventually resumed domination. In another, MGZ activity developed and remained stable for the entire simulation in a 5 x 10-cm substrate (for an example of MGZ activity, see Fig. 5 and the associated discussion below). In one case, MGZ activity terminated rapidly in the 4 x 8-cm substrate.



View larger version (58K):
[in this window]
[in a new window]
 
Fig. 5. Atrial activity in a substrate of size 10 x 20 cm with mean [ACh] of 1.875 nM and [ACh] distribution of 2.0 cm. A: atrial TMP snapshot (voltage scale to right). Arrows indicate wavefront propagation directions. B: PS trajectory map superimposed on the [ACh] distribution grid ([ACh] scale to right). C: PS lifespan histogram. D: APs from a low-[ACh] zone (top, position shown by * in B) and from a high-[ACh] zone (bottom, position shown by # in B). E: FFT of activity from the high-[ACh] zone corresponding to APs in D, bottom. (See also on-line movie 4.)

 
A consistent feature of the highest mean [ACh] simulations was strong anchoring in low-[ACh] regions, which tended to prevent drift of the primary rotor away from its initial location, as well as sharp curvature of daughter waves that result from wavebreak against other low-[ACh] zones, which favored mutual annihilation and prevented the production of long-lasting secondary rotors. Primary rotors were unlikely to terminate by collision with boundaries or with secondary waves, because of strong anchoring. Consequently, inducible fibrillatory activity was long lasting at the highest mean [ACh] for the majority of simulations at all substrate sizes, and substrate size was not a major determinant of AF persistence under this condition.

In contrast to the well-anchored primary rotors observed at the highest mean [ACh], at lower mean [ACh]s (with smaller spatial [ACh] gradients), primary rotors tended to drift away from their initial locations. Primary rotor drift favored annihilation against a boundary or a daughter wave. In smaller sheets, primary rotor extinction due to such drift usually terminated AF. Figure 4 and on-line movie 3 (see http://ajheart.physiology.org/cgi/content/full/00252.2005/DC1) show an example of this type of activity in a 4 x 8-cm substrate with the lowest mean [ACh] (1.875 nM) and a 2-cm [ACh] distribution distance. APs from a low-[ACh] zone are shown in Fig. 4A, and the PS trajectory map is shown in Fig. 4B. TMP maps showing the progression of activity are provided in Fig. 4C. S2 initiated a single clockwise-rotating rotor toward the middle of the substrate (Fig. 4C, top left). The PS at the core tip of this spiral showed much more meander (green trajectory in Fig. 4B) than for the example shown in Fig. 2, moving from the initial low-[ACh] area toward a low-[ACh] area to the upper left of the substrate at ~880 ms, returning toward the middle of the substrate (Fig. 4C, top right and bottom left), and then moving to another low-[ACh] area to the lower right part of the substrate at ~2,300 ms (Fig. 4C, bottom right). The latter area is adjacent to a border of the grid, and AF terminated at ~2,400 ms when the primary rotor drifted to the lower right corner and extinguished. In this example, several intermediate-lasting rotors established themselves temporarily in the lower left and far left portions of the grid, as reflected by the blue trajectories in Fig. 4B, but were annihilated by collisions with other waves or with borders. Because no longer-lasting secondary rotors were available to assume generator function at the time of primary generator rotor extinction, AF terminated.



View larger version (46K):
[in this window]
[in a new window]
 
Fig. 4. A: atrial activity recorded from a low-[ACh] zone in a substrate of size 4 x 8 cm with mean [ACh] of 1.875 nM and [ACh] distribution of 2.0 cm. B: PS trajectory map superimposed on the [ACh] distribution grid ([ACh] scale to right). C: snapshots of TMP at times indicated. For a discussion, see text (see also on-line movie 3).

 
The most important role of substrate size in determining AF maintenance was noted for the lowest mean [ACh]s. Figure 5 and on-line movie 4 (see http://ajheart.physiology.org/cgi/content/full/00252.2005/DC1) illustrate a simulation with the same [ACh] and spatial distribution as in Fig. 4 but in the largest computational substrate (10 x 20 cm). Because pairs of counterrotating secondary wavelets were not restricted in location by strong functional anchoring, they could avoid annihilation by separating from each other, leading to longer-lasting daughter waves that took over generator function upon initial primary rotor extinction. Reentry began with a single rotor in a central low-[ACh] area (see on-line movie 4). About 1 s after the initiation of reentry, a pair of counterrotating daughter waves was generated in an adjacent low-[ACh] zone and failed to annihilate against each other. One of them drifted to an adjacent, third low-[ACh] zone. Similar events continued to occur over time, with intermediate-lasting daughter waves avoiding annihilation by detaching from the low-[ACh] areas in which they were formed and drifting to adjacent areas. The number of intermediate- and long-lasting (blue and green PS) rotors was much larger than in the smaller substrate shown in Fig. 4. In addition, more low-[ACh] zones were located at a distance from boundaries, making extinction of spiral wave generators at boundaries less likely. Consequently, AF was sustained by multiple coexisting spiral wave generators that formed, lasted for varying periods of time and then annihilated against invading wavefronts, only to be replaced by new generators.

TMP maps during this MGZ type of activity, as illustrated in Fig. 5A, resembled classical multiple circuit reentry. PS trajectory maps (Fig. 5B) showed numerous zones of long-lasting (green) and intermediate-lasting (blue) generator activity. The PS histogram in this case (Fig. 5C) showed 39 PSs lasting >1 s, reflecting a large number of generator rotors over the course of 20 s of activity. However, the longest-lasting PS was present for 1.79 s (in contrast to 20 s for examples of sustained SG activity), indicating the tendency of even long-lasting generators to be annihilated and the need for multiple generators to ensure continued activity. APs showed highly irregular activity (Fig. 5D), and the power spectrum was much more broad based (Fig. 5E) than for activity maintained by one or successive single generators (e.g., compare with Fig. 2, E and G). Sustained MGZ activity clearly did not depend on the persistence of a single primary generator but rather on the presence of multiple shorter-lasting generators able to maintain activity at all times.

In contrast to the highest mean [ACh] substrate, for which MGZ activity occurred in only three simulations, MGZ was observed for seven simulations at the lowest mean [ACh] and nine simulations for the intermediate mean [ACh] (Table 2). The duration of activity maintained exclusively by MGZ reentry showed a strong dependency on substrate size at the two lower mean [ACh] values, with such activity lasting under 1 s in the 4 x 8-cm substrate, <4 s in the 5 x 10-cm substrate, up to 10 s in the 7.5 x 15-cm substrate, and between 15 and 20 s in six of seven cases in the 10 x 20-cm substrate.

Properties associated with different mechanisms of AF maintenance. Two important contributors to perpetuation of activity were noted. The occurrence of stable SG activity was closely related to mean [ACh], with longer-lasting SG and SSG activity being very likely at high [ACh]s irrespective of substrate size. Low mean [ACh]s tended to be associated with unstable, meandering primary generators, with AF persistence relying on MGZ mechanisms linked to substrate size.

Given the critical importance of mean [ACh] in governing the mechanisms underlying persistent activity in our model, we assessed the effects of [ACh] and [ACh] distribution on APD values and heterogeneity. Figure 6 shows the relationships between mean [ACh] and spatial [ACh] distribution on one hand and mean APD75 and APD gradient (difference between largest and smallest APD75 at a frequency of 6 Hz). The primary determinant of mean APD75 was clearly mean [ACh] (Fig. 6A). Larger spatial [ACh] gradients were associated with slightly larger mean APD75 values, likely because of reduced electrotonic influences of short APD75 (high [ACh]) regions on the longer APD75 zones. The APD gradient was predominantly determined by [ACh] distribution (Fig. 6B), presumably because larger peak to trough [ACh] distances resulted in reduced APD smoothing by electrotonic interaction. Mean [ACh] was a smaller determinant, tending to produce larger APD gradients as mean concentration increased. Mean [ACh] was also a significant determinant of activation rate during AF. The dominant frequency in the FFT histogram in the low-[ACh] core zone and in the high-[ACh] periphery of spiral wave generators averaged 8.7 ± 0.1 and 10.4 ± 0.1 Hz, respectively, at 1.75 nM mean [ACh]; 9.3 ± 0.1 and 11.4 ± 0.1 Hz at 3.5 nM; and 9.6 ± 0.1 and 12.1 ± 0.1 Hz at 7.5 nM. The larger dominant frequencies at higher [ACh]s reflect the shorter APD values, resulting in sharper rotation of reentrant wavefronts and smaller rotor cores.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 6. A: mean AP duration at 75% repolarization (APD75) in each simulation as a function of [ACh] distribution and mean [ACh]. B: APD gradient (Grad = minimum – maximum APD75) in each simulation as a function of [ACh] distribution and mean [ACh].

 
The number of PSs reflects the number of wavebreaks that can potentially form daughter wave generators to sustain MGZ reentry. Figure 7 shows an analysis of the relationships among mean [ACh], [ACh] distribution, and substrate size on one hand and the number of PSs normalized to observation time (in ms). Larger mean [ACh] tended to decrease PS numbers (Fig. 7A), presumably by limiting drift, and smaller [ACh] distribution distances increased PS numbers modestly by increasing the number of low-[ACh] zones per substrate (Fig. 7B). However, substrate size was clearly the strongest determinant of the number of PSs.



View larger version (30K):
[in this window]
[in a new window]
 
Fig. 7. Number of PSs per unit time as a function of substrate dimensions and mean [ACh] (A), substrate dimensions and [ACh] gradient (B), and relationship between PSs/time in each simulation and substrate size along with means ± SD (C). (Results are not shown for the 3 x 6-cm substrate because AF was only inducible in 1 of 12 [ACh]/[ACh] distribution conditions for this substrate size.)

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
We used a mathematical model of canine atrial tissue to evaluate the role of substrate dimensions in determining the persistence of fibrillatory activity. Continuous activity could be maintained either under conditions that produced very stable primary generator rotors or conditions that resulted in substantial rotor meander and a high probability of eventual extinction of individual generator sources. The role of substrate dimension was most important for conditions associated with unstable meandering rotors, because the maintenance of continuous activity required the existence of a large enough number of longer-lasting rotors that extinction of all simultaneously was unlikely. For stable generator rotor conditions, fibrillatory activity, once successfully induced, tended to remain stable irrespective of substrate size.

Relationship to previous conceptual models of AF. There are essentially three conceptual schemas of AF, all rooted in notions developed in the early 20th century: 1) multiple circuit reentry, 2) single circuit reentry; and 3) sustained rapid ectopic activity (31). Garrey (8) and Mines (25) were the first proponents of multiple circuit reentry as the mechanism underlying AF. Moe (27) refined the multiple reentry circuit schema with his "multiple wavelet hypothesis" and investigated the conditions promoting fibrillatory activity in a computer model. Lewis (20) held that AF is due to a single primary rapidly rotating wave, with wavefront break-up against variably refractory tissue producing fibrillatory activity. For many years, the multiple circuit model was preeminent, with substantial support from clinical and experimental evidence (reviewed in Refs. 30 and 31). Lewis' mother wave concept has recently received support from a range of experimental observations that point to single dominant rotors underlying AF in a variety of experimental paradigms (4, 5, 14, 24, 45).

When we initially developed the model described in the present study, we noted that AF was generally maintained by a single primary generator rotor (18). In the present study, we extended this work to a range of atrial substrate dimensions and to a longer simulation time, applying a recently developed method for automated PS identification and tracking (51) that permits the analysis of the complex activity that occurs in larger atrial substrates. This allowed us to refine the ideas that emerged from our previous analysis, by indicating that either dominant rotors (similar to Lewis' mother wave) or multiple unstable reentrant generators (essentially multiple wavelet reentry) can underlie fibrillatory activity in the model, depending on the functional properties and dimensions of the atrial substrate. To our knowledge, the present analysis is the first to identify in a mathematical tissue model a spiral wave counterpart to Moe's multiple wavelet hypothesis of AF.

We found that the MGZ mechanism, corresponding to multiple circuit reentry, was less likely than single rotor mechanisms (SG and SSG) to maintain fibrillatory activity for all conditions but the largest substrate. This finding supports the evolving notion that single rotor reentry with fibrillatory conduction is an important, and possibly predominant mechanism, for the maintenance of AF (4, 14, 24, 40). However, in larger, diseased atria, mechanisms corresponding to MGZ activity may become more important and may account for greater difficulty in AF suppression and prevention.

Relationship to experimental evidence regarding AF mechanisms. Detailed experimental evidence for multiple wavelet reentry in cholinergic AF was first presented in 1985 (1). Subsequent studies provided epicardial mapping evidence supporting multiple circuit reentry in vagal (47) and tachycardia remodeling AF (9). The role of atrial wavelength as a determinant of AF provided further support for multiple wavelet reentry (35). However, other evidence has called into question the role of multiple circuit reentry. Rapid left atrial activity, often with an apparent dominant generator zone, underlies cholinergic AF in isolated sheep hearts (4, 14, 24, 40). In these studies, most wavelets were generated by passive wave break (4). In the present analysis, either single dominant rotors (Figs. 2 and 3 and on-line movies 1 and 2) or continuous unstable generator activity (Fig. 5 and on-line movie 4) sustained fibrillatory activity. Our observation of SG AF in very small substrates supports observations of fibrillation in mouse hearts (41) and the notion that such fibrillation may not require multiple circuit reentry.

Several of our observations parallel findings in experimental studies of cholinergic AF. ACh alters cardiac electrophysiology primarily by activating IKACh (18). IKACh is particularly important in atrial tissue and, when activated, can produce atrial hyperpolarization (removing voltage-dependent Na+ current inactivation) and cause striking AP shortening. These effects provide greater source current for impulse conduction and allow for more rapid spiral wave rotation, accelerating and stabilizing atrial reentry (18). The importance of ACh actions are compounded by the spatial heterogeneity of effects associated with parasympathetic stimulation of the heart (22), which produces zones of lower [ACh]/effect that serve as functional barriers, and therefore anchoring points, in close proximity to areas of higher [ACh]/action. These properties are analogous to the changes induced by altered amplitude of K+ current, IK1, which appears to be a significant determinant of the stability of reentry and fibrillation in both atrial (33) and ventricular (37) tissues. We noted in the present study that lower mean [ACh]s were associated with meandering unstable rotors, whereas at higher concentrations single rotors stabilized and dominated activity. Schuessler et al. (39) observed multiple reentry circuits at moderately high [ACh] but single, dominant rotors that controlled the rhythm at the highest [ACh]s in canine atrial preparations. Ikeda et al. (12) noted meandering unstable wavefronts in the canine right atrium exposed to ACh. They subsequently described rapid nonstationary wavefronts during cholinergic AF in the canine right atrium, but stationary wavefronts associated with the shortest APDs, <30 ms (48).

Role of tissue size in AF maintenance. Garrey (8) first pointed to tissue size as a determinant of the likelihood of fibrillation, a notion supported by the persistence of AF in animals with different cardiac sizes (28). Atrial size is a significant determinant of clinical AF (10). Surgical procedures that reduce left atrial size may prevent recurrent AF (7, 38), and some of the benefit from linear ablation in AF may be due to decreased left atrial size (2). Smaller left atria are associated with higher conversion rates from AF (6). Our study supports the role of atrial size as a determinant of AF maintenance and provides potential mechanistic insights.

Consideration of the model. Compared with Moe's original computer model of AF (27), the model we used is much more realistic. Moe's model used a 992-unit array, with each unit having 5 discrete excitability states governed by an algebraic expression. The relative refractory period was much longer than in intact tissue, and there were no electrotonic interactions between functional units. Our model is based on cells of physiological size with APs controlled by realistic mathematical representations of canine ionic currents, with APs responding to various conditions similarly to experimental observations (34). Coupling, conduction, and anisotropic properties are appropriate for the dog (18).

On the other hand, the cellular substrate geometry is limited to a two-dimensional rectangular array. Three-dimensional complexities of in vivo atria, such as holes (venae cavae, pulmonary veins, and mitral and tricuspid valves), structural complexity (pectinate muscles and crista terminalis), and tissues with specific properties (septum and pulmonary veins) are not included. A more geometrically realistic three-dimensional model of the atria has been described (43) and provides insights into the potential involvement of various geometric structures in AF. However, other recent modeling studies of AF have similarly used two-dimensional planar arrays (32, 50), which remain instructive for asking specific questions such as the ones we posed in the present study. It will be interesting to determine in future work whether the principles defined in the present analysis continue to hold in more geometrically realistic models. Despite the simplifications inherent in our model, it has been shown in previous studies to agree closely with a range of experimental observations (16, 18) and with the results of more complex three-dimensional models (44). Of interest, the frequency peaks in FFT spectra in the present study are quite comparable to previous observations during vagal AF in the anesthetized dog (22, 46, 47).

We have previously observed that high, homogeneous [ACh]s stabilize reentry and greatly reduce meander of the spiral wave core tip (18). In the present study, we note that this anchoring effect of high [ACh]s results in very stable dominant generator-based AF in a heterogeneous substrate. It is likely that similar stabilization could occur with other factors that anchor primary generator spiral waves, such as anatomic obstacles and pectinate muscles, as suggested by previous experimental observations (13, 49).

We used a sinusoidal distribution to govern local [ACh] and produce spatially heterogeneous vagal action. Spatial heterogeneity in vagal effects is well recognized and is an important contributor to AF promotion (22, 46). The distribution of extracellular [ACh] during vagal stimulation in vivo has not been defined, but a sinusoidal distribution is clearly an oversimplification. During vagal AF in vivo, there are substantial differences in local frequency over distances of ~1–2 cm (22, 46), of the order of [ACh] distributions used in the present study. Our model has discrete boundary conditions. Boundary conditions do exist in the intact cardiac atria, notably at the septum, great veins, and atrioventricular valves. Ikeda et al. (11) noted spontaneous AF termination in isolated canine right atrial preparations by movement to the tissue border. Similar to findings in the present study, Skanes et al. (40) reported rotor termination upon collision with other wavefronts during AF in isolated sheep hearts.

Our study was designed to address a very specific question: how substrate dimensions affect the persistence of fibrillatory activity in a two-dimensional model of canine atrium. Many of our findings are similar to those observed in experimental contexts, as outlined above. However, we are in no way claiming that our simplified model reflects the full biological complexity of AF, particularly as observed clinically with the complicating features of cardiac disease, geometric complexities, autonomic factors, etc. In addition, our model deals only with fibrillatory activity of the reentrant type and is not relevant to forms of AF that may be due predominantly to rapid ectopic activity as believed to arise from the pulmonary veins in some patient populations (31).

In conclusion, we found that atrial substrate size is a potentially important determinant of the maintenance of fibrillatory activity in our model. Unlike our previous study, which found a predominance of single rotor mechanisms in a substrate of fixed, 5 x 10-cm dimensions, the present analysis shows that either single dominant rotors or multiple unstable rotor-based reentry can underlie persistent fibrillatory activity in an atrial model. Although AF can be maintained in relatively small substrates by single dominant rotor mechanisms, larger substrates allow for persistent AF more easily because they can support a sufficient number of meandering rotors that the probability of all being extinguished simultaneously is very low, a mechanism that cannot occur in smaller substrates. These observations support the notion that larger substrates favor AF persistence but indicate that the importance of this phenomenon depends on the mechanisms underlying AF. Furthermore, these findings speak to recent debates regarding the relative importance of single driver rotors versus multiple reentrant circuit mechanisms by suggesting that either may maintain AF depending on atrial size and electrophysiological properties.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This study was supported by the Canadian Institutes of Health Research, the Mathematics for Information Technology and Complex Systems Network of Centers of Excellence, and the Quebec Heart and Stroke Foundation. R. Zou is a Heart and Stroke Foundation of Canada research fellow.


    ACKNOWLEDGMENTS
 
The authors thank France Thériault for excellent secretarial help with the manuscript and the Réseau de Quebec de Calcul de Haute Performance for the use of sophisticated parallel processing facilities.


    FOOTNOTES
 

Address for reprint requests and other correspondence: S. Nattel, Research Center, Montreal Heart Institute, 5000 Belanger St., Montreal, Quebec, Canada H1T 1C8 (E-mail: stanley.nattel{at}icm-mhi.org)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Allessie MA, Lammers WJEP, Bonke FIM, and Hollen J. Experimental evaluation of Moe's multiple wavelet hypothesis of atrial fibrillation. In: Cardiac Arrhythmias, edited by Zipes DP and Jalife J. Philadelphia, PA: Grune & Stratton, 1985, p. 265–276.
  2. Avitall B, Urbonas A, Urboniene D, Millard S, and Helms R. Time course of left atrial mechanical recovery after linear lesions: normal sinus rhythm versus a chronic atrial fibrillation dog model. J Cardiovasc Electrophysiol 11: 1397–1406, 2000.[CrossRef][Web of Science][Medline]
  3. Bray MA, Lin SF, Aliev RR, Roth BJ, and Wikswo JP Jr. Experimental and theoretical analysis of phase singularity dynamics in cardiac tissue. J Cardiovasc Electrophysiol 12: 716–722, 2001.[CrossRef][Web of Science][Medline]
  4. Chen J, Mandapati R, Berenfeld O, Skanes AC, Gray RA, and Jalife J. Dynamics of wavelets and their role in atrial fibrillation in the isolated sheep heart. Cardiovasc Res 48: 220–232, 2000.[Abstract/Free Full Text]
  5. Derakhchan K, Li D, Courtemanche M, Smith B, Brouillette J, Pagé PL, and Nattel S. Method for simultaneous epicardial and endocardial mapping of in vivo canine heart: application to atrial conduction properties and arrhythmia mechanisms. J Cardiovasc Electrophysiol 12: 548–555, 2001.[CrossRef][Web of Science][Medline]
  6. Eversole A, Hancock W, Johns T, Lopez LM, and Conti CR. Ibutilide: efficacy and safety in atrial fibrillation and atrial flutter in a general cardiology practice. Clin Cardiol 24: 521–525, 2001.[Web of Science][Medline]
  7. Garcia-Villarreal OA, Rodriguez H, Trevino A, Gouveia AB, and Arguero R. Left atrial reduction and mitral valve surgery: the "functional-anatomic unit" concept. Ann Thorac Surg 71: 1044–1045, 2001.[Abstract/Free Full Text]
  8. Garrey WE. The nature of fibrillatory contraction of the heart–its relation to tissue mass and form. Am J Physiol 33: 397–414, 1914.[Free Full Text]
  9. Gaspo R, Bosch RF, Talajic M, and Nattel S. Functional mechanisms underlying tachycardia-induced sustained atrial fibrillation in a chronic dog model. Circulation 96: 4027–4035, 1997.[Abstract/Free Full Text]
  10. Henry WL, Morganroth J, Pearlman AS, Clark CE, Redwood DR, Itscoitz SB, and Epstein SE. Relation between echocardiographically determined left atrial size and atrial fibrillation. Circulation 53: 273–279, 1976.[Abstract/Free Full Text]
  11. Ikeda T, Uchida T, Hough D, Lee JJ, Fishbein MC, Mandel WJ, Chen PS, and Karagueuzian HS. Mechanism of spontaneous termination of functional reentry in isolated canine right atrium. Evidence for the presence of an excitable but nonexcited core. Circulation 94: 1962–1973, 1996.[Abstract/Free Full Text]
  12. Ikeda T, Wu TJ, Uchida T, Hough D, Fishbein MC, Mandel WJ, Chen PS, and Karagueuzian HS. Meandering and unstable reentrant wave fronts induced by acetylcholine in isolated canine right atrium. Am J Physiol Heart Circ Physiol 273: H356–H370, 1997.[Abstract/Free Full Text]
  13. Ikeda T, Yashima M, Uchida T, Hough D, Fishbein MC, Mandel WJ, Chen PS, and Karagueuzian HS. Attachment of meandering reentrant wave fronts to anatomic obstacles in the atrium. Role of the obstacle size. Circ Res 81: 753–764, 1997.[Abstract/Free Full Text]
  14. Jalife J, Berenfeld O, Skanes A, and Mandapati R. Mechanisms of atrial fibrillation: mother rotors or multiple daughter wavelets, or both? J Cardiovasc Electrophysiol 9,Suppl8: S2–S12, 1998.
  15. Joyner RW. Effects of the discrete pattern of electrical coupling on propagation through an electrical syncytium. Circ Res 50: 192–200, 1982.[Abstract/Free Full Text]
  16. Kneller J, Kalifa J, Zou R, Zaitsev AV, Warren M, Berenfeld O, Vigmond EJ, Leon LJ, Nattel S, and Jalife J. Mechanisms of AF termination by pure sodium channel blockade in an ionically-realistic mathematical model. Circ Res 96: e35–e47, 2005.[Abstract/Free Full Text]
  17. Kneller J, Ramirez RJ, Chartier D, Courtemanche M, and Nattel S. Time-dependent transients in an ionically based mathematical model of the canine atrial action potential. Am J Physiol Heart Circ Physiol 282: H1437–H1451, 2002.[Abstract/Free Full Text]
  18. Kneller J, Zou R, Vigmond EJ, Wang Z, Leon LJ, and Nattel S.Atrial fibrillation in a computer model of a 2-dimensional sheet of canine atrial cells with realistic ionic properties. Circ Res 90: 73e–87e, 2002.[Abstract/Free Full Text]
  19. Leon LJ, Roberge FA, and Vinet A. Simulation of two-dimensional anisotropic cardiac reentry: effects of the wavelength on the reentry characteristics. Ann Biomed Eng 22: 592–609, 1994.[CrossRef][Web of Science][Medline]
  20. Lewis T. The Mechanism and Graphic Registration of the Heart Beat (3rd ed.). London: Shaw, 1925, p. 319–374.
  21. Li D, Fareh S, Leung TK, and Nattel S. Promotion of atrial fibrillation by heart failure in dogs: atrial remodeling of a different sort. Circulation 100: 87–95, 1999.[Abstract/Free Full Text]
  22. Liu L and Nattel S. Differing sympathetic and vagal effects on atrial fibrillation in dogs: role of refractoriness heterogeneity. Am J Physiol Heart Circ Physiol 273: H805–H816, 1997.[Abstract/Free Full Text]
  23. MacWilliam JA. Fibrillar contraction of the heart. J Physiol 8: 296–310, 1887.
  24. Mandapati R, Skanes A, Chen J, Berenfeld O, and Jalife J. Stable micro-reentrant sources as a mechanism of atrial fibrillation in the isolated sheep heart. Circulation 101: 194–199, 2000.[Abstract/Free Full Text]
  25. Mines GR. On circulating excitations in heart muscles and their possible relation to tachycardia and fibrillation. Proc Trans R Soc Can 8: 43–53, 1914.
  26. Moe GK and Abildskov JA. Atrial fibrillation as a self-sustaining arrhythmia independent of focal discharge. Am Heart J 58: 59–70, 1959.[CrossRef][Web of Science][Medline]
  27. Moe GK, Rheinboldt WC, and Abildskov JA. A computer model of atrial fibrillation. Am Heart J 67: 200–220, 1964.[CrossRef][Web of Science][Medline]
  28. Moore EN and Spear JF. Electrophysiological studies on atrial fibrillation. Heart Vessels Suppl 2: 32–39, 1987.[Medline]
  29. Nattel S, Khairy P, Roy D, Thibault B, Guerra P, Talajic M, and Dubuc M. New approaches to atrial fibrillation management–a critical review of a rapidly-evolving field. Drugs 62: 2377–2397, 2002.[CrossRef][Web of Science][Medline]
  30. Nattel S, Li D, and Yue L. Basic mechanisms of atrial fibrillation–very new insights into very old ideas. Annu Rev Physiol 62: 51–77, 2000.[CrossRef][Web of Science][Medline]
  31. Nattel S. New ideas about atrial fibrillation 50 years on. Nature 415: 219–226, 2002.[CrossRef][Medline]
  32. Nygren A, Leon LJ, and Giles WR. Simulations of the human atrial action potential. Philos Trans R Soc Lond 359: 1111–1125, 2001.
  33. Pandit SV, Berenfeld O, Anumonwo J, Zaritski R, Kneller J, Nattel S, and Jalife J. Ionic determinants of functional reentry in a 2-D model of human atrial cells during simulated chronic atrial fibrillation. Biophys J 88: 3806–3821, 2005.[CrossRef][Web of Science][Medline]
  34. Ramirez RJ, Nattel S, and Courtemanche M. Mathematical analysis of canine atrial action potentials: rate, regional factors, and electrical remodeling. Am J Physiol Heart Circ Physiol 279: H1767–H1785, 2000.[Abstract/Free Full Text]
  35. Rensma PL, Allessie MA, Lammers WJ, Bonke FI, and Schalij MJ. Length of excitation wave and susceptibility to reentrant atrial arrhythmias in normal conscious dogs. Circ Res 62: 395–410, 1988.[Abstract/Free Full Text]
  36. Roberge FA, Vinet A, and Victorri B. Reconstruction of propagated electrical activity with a two-dimensional model of anisotropic heart muscle. Circ Res 58: 461–475, 1986.[Abstract/Free Full Text]
  37. Samie FH, Berenfeld O, Anumonwo J, Mironov SF, Udassi S, Beaumont J, Taffet S, Pertsov AM, and Jalife J. Rectification of the background potassium current: a determinant of rotor dynamics in ventricular fibrillation. Circ Res 89: 1216–1223, 2001.[Abstract/Free Full Text]
  38. Scherer M, Dzemali O, Aybek T, Wimmer-Greinecker G, and Moritz A. Impact of left atrial size reduction on chronic atrial fibrillation in mitral valve surgery. J Heart Valve Dis 12: 469–474, 2003.[Web of Science][Medline]
  39. Schuessler RB, Grayson TM, Bromberg BI, Cox JL, and Boineau JP. Cholinergically mediated tachyarrhythmias induced by a single extrastimulus in the isolated canine right atrium. Circ Res 71: 1254–1267, 1992.[Abstract/Free Full Text]
  40. Skanes AC, Mandapati R, Berenfeld O, Davidenko JM, and Jalife J. Spatio-temporal periodicity during atrial fibrillation in the isolated sheep heart. Circulation 98: 1236–1248, 1998.[Abstract/Free Full Text]
  41. Vaidya D, Morley GE, Samie FH, and Jalife J. Reentry and fibrillation in the mouse heart. A challenge to the critical mass hypothesis. Circ Res 85: 174–181, 1999.[Abstract/Free Full Text]
  42. Victorri B, Vinet A, Roberge FA, and Drouhard JP. Numerical integration in the reconstruction of cardiac action potentials using Hodgkin-Huxley-type models. Comput Biomed Res 18: 10–23, 1985.[CrossRef][Web of Science][Medline]
  43. Vigmond EJ, Ruckdeschel R, and Trayanova N. Reentry in a morphologically realistic atrial model. J Cardiovasc Electrophysiol 12: 1046–1054, 2001.[CrossRef][Web of Science][Medline]
  44. Vigmond EJ, Tsoi V, Kuo S, Arevalo H, Kneller J, Nattel S, and Trayanova N. The effect of vagally induced APD dispersion on atrial arrhythmogenesis. Heart Rhythm 1: 334–344, 2004.
  45. Waldo AL. Mechanisms of atrial flutter and atrial fibrillation: distinct entities or two sides of a coin? Cardiovasc Res 54: 217–229, 2002.[Free Full Text]
  46. Wang J, Liu L, Feng J, and Nattel S. Regional and functional factors determining induction and maintenance of atrial fibrillation in dogs. Am J Physiol Heart Circ Physiol 271: H148–H158, 1996.[Abstract/Free Full Text]
  47. Wang Z, Pagé P, and Nattel S. Mechanism of flecainide's antiarrrhythmic action in experimental atrial fibrillation. Circ Res 71: 271–287, 1992.[Abstract/Free Full Text]
  48. Wu TJ, Kim YH, Yashima M, Athill CA, Ting CT, Karagueuzian HS, and Chen PS. Progressive action potential duration shortening and the conversion from atrial flutter to atrial fibrillation in the isolated canine right atrium. J Am Coll Cardiol 38: 1757–1765, 2001.[Abstract/Free Full Text]
  49. Wu TJ, Yashima M, Xie F, Athill CA, Kim YH, Fishbein MC, Qu Z, Garfinkel A, Weiss JN, Karagueuzian HS, and Chen PS. Role of pectinate muscle bundles in the generation and maintenance of intra-atrial reentry: potential implications for the mechanism of conversion between atrial fibrillation and atrial flutter. Circ Res 83: 448–462, 1998.[Abstract/Free Full Text]
  50. Xie F, Qu Z, Garfinkel A, and Weiss JN. Electrical refractory period restitution and spiral wave reentry in simulated cardiac tissue. Am J Physiol Heart Circ Physiol 283: H448–H460, 2002.[Abstract/Free Full Text]
  51. Zou R, Kneller J, Leon LJ, and Nattel S. Development of a computer algorithm for the detection of phase singularities and initial application to analyze simulations of atrial fibrillation. Chaos 12: 764–778, 2002.[CrossRef][Web of Science][Medline]



This article has been cited by other articles:


Home page
EuropaceHome page
K. Nishida, G. Michael, D. Dobrev, and S. Nattel
Animal models for atrial fibrillation: clinical insights and scientific opportunities
Europace, October 29, 2009; (2009) eup328v1.
[Abstract] [Full Text] [PDF]


Home page
Circ Arrhythm ElectrophysiolHome page
S. Nattel, B. Burstein, and D. Dobrev
Atrial Remodeling and Atrial Fibrillation: Mechanisms and Implications
Circ Arrhythm Electrophysiol, April 1, 2008; 1(1): 62 - 73.
[Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
G. Laurent, G. Moe, X. Hu, H. Leong-Poi, K. A. Connelly, P. P.-S. So, A. Ramadeen, L. Doumanovskaia, A. Konig, J. Trogadis, et al.
Experimental studies of atrial fibrillation: a comparison of two pacing models
Am J Physiol Heart Circ Physiol, March 1, 2008; 294(3): H1206 - H1215.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
N. Ono, H. Hayashi, A. Kawase, S.-F. Lin, H. Li, J. N. Weiss, P.-S. Chen, and H. S. Karagueuzian
Spontaneous atrial fibrillation initiated by triggered activity near the pulmonary veins in aged rats subjected to glycolytic inhibition
Am J Physiol Heart Circ Physiol, January 1, 2007; 292(1): H639 - H648.
[Abstract] [Full Text] [PDF]


Home page
Exp PhysiolHome page
M. P. Nash, C. P. Bradley, P. M. Sutton, R. H. Clayton, P. Kallis, M. P. Hayward, D. J. Paterson, and P. Taggart
Whole heart action potential duration restitution properties in cardiac patients: a combined clinical and modelling study
Exp Physiol, March 1, 2006; 91(2): 339 - 354.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
I. R. Efimov and V. V. Fedorov
Chessboard of atrial fibrillation: reentry or focus? Single or multiple source(s)? Neurogenic or myogenic?
Am J Physiol Heart Circ Physiol, September 1, 2005; 289(3): H977 - H979.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Movie
Right arrow All Versions of this Article:
289/3/H1002    most recent
00252.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (9)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zou, R.
Right arrow Articles by Nattel, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zou, R.
Right arrow Articles by Nattel, S.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2005 by the American Physiological Society.