Background

The mechanisms involved in onset, maintenance, and termination of atrial fibrillation are not well understood. A biophysical model could be useful to determine how the events unfold.
Method

A two-dimensional cellular automaton consisting of 576 × 576 grid nodes was implemented to demonstrate the types of electrical activity that may occur in compromised atrial substrate. Electrical activation between nodes was made anisotropic (2:1), and the refractory period (RP) was adjusted from 74 to 192 ms in the spatial domain. Presence of collagen fibers were simulated as short lines of conduction block at many random grid sites, while ablation lesions were delineated as longer lines of block. An S1-S2 pulse from one grid corner was utilized to initiate simulated electrical activity. Simulations were done in which 1. no ablation lines, 2. random ablation lines, and 3. parallel ablation lines were added to the grid to determine how this affected the formation and annihilation of rotational activity after S1-S2 stimulation.
Results

As the premature (S2) wavefront traversed the grid, rotational activity formed near boundaries where wavefronts propagated from shorter to longer refractory regions, causing unidirectional block, and were anchored by fiber clusters. Multiple wavelets appeared when wavefronts originating from different driving rotational features collided, and/or by their encounter with RP discontinuities. With the addition of randomly orientated simulated ablation lesions, followed by reinduction of fibrillatory activity, mean activation interval (AI) prolonged from a baseline level of 144.2 ms–160.3 ms (p < 0.001 in most comparisons). During fibrillatory activity, when parallel ablation lines were added to short RP regions, AI prolonged to 150.4 ms (p < 0.001), and when added to long RP regions, AI prolonged to 185.3 ms (p < 0.001). In all cases, AI prolongation after simulated ablation resulted from reduced number and/or from the isolation of local drivers, so that distant drivers in short RP regions activated long RP regions N:1, while distant drivers in long RP regions activated short RP regions at a relatively slow rate.
Conclusions

An automaton model was found useful to generate and test hypotheses concerning fibrillatory activity, which can then be validated in the clinical electrophysiology laboratory.