Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and sustenance of wave sources in a model of the PV tissue.

A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced electrical remodelling are introduced based on published experimental data, as changes of ion channel currents (I_CaL, I_K1, I_to, and I_Kur), the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5–12 ms) and anisotropy (conduction velocities of 0.75–1.25 and 0.21–0.31 m/s along and transverse to atrial fibres) can results in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wavelength (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulnerability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (I_K1 and I_Kur blockers) can also increase the heterogeneity (AP gradients of 26–27 ms) and the likelihood of wave breaks.

Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrangement enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological interventions.