Nonuniform Muscle Fiber Orientation Causes Spiral Wave Drift in a Finite Element Model of Cardiac Action Potential Propagation

Abstract

Reentrant tachyarrhythmias are thought to involve spiral waves of excitation and recovery that may be nonstationary. The effect of muscle fiber curvature on spiral and planar wave propagation was studied using a computational model. Two-dimensional anisotropic cardiac propagation was modeled using a finite element method to solve a modification of the FitzHugh-Nagumo equations. Spiral waves that propagated stably about a fixed core in tissue with a uniform fiber orientation were found to drift at an oblique angle to the fibers when the fibers curved. The drift velocity was linearly related to the fiber angle gradient and was 10% of the longitudinal propagation velocity with a gradient of 4 degree/cm. Planar wavefronts were also affected by fiber curvature. The maximal upstroke rate, propagation velocity, and the action potential amplitude all increased when the fibers curved toward the wavefront and decreased when they curved away. For example, when the fibers curved toward the wavefront with a moderate gradient of 15 degree/cm, maximal upstroke rate increased 74%, transverse propagation velocity increased 65%, and action potential amplitude increased 9%. This phenomenon caused the spiral wave drift: As a spiral wave traverses a cycle, the angle between the wavefront at the wavetip and the fibers changes periodically, thus altering the propagation parameters. These periodic changes affect the instantaneous radius of curvature of the wavetip path, which results in drift. Spiral and planar waves are affected by muscle fiber curvature. The resulting dynamics may be important in determining the lifetime and stability of reentrant arrhythmias.