Drift of scroll waves in thin layers caused by thickness features: asymptotic theory and numerical simulations

Chamber-specific wall thickness features in human atrial fibrillation

Persistent atrial fibrillation (AF) is not effectively treated due to a lack of adequate tools for identifying patient-specific AF substrates. Recent studies revealed that in 30-50% of patients, persistent AF is maintained by localized drivers not only in the left atrium (LA) but also in the right atrium (RA). The chamber-specific atrial wall thickness (AWT) features underlying AF remain elusive, though the important role of AWT in AF is widely acknowledged. We aimed to provide direct evidence of the existence of distinguished RA and LA AWT features underlying AF drivers by analysing functionally and structurally mapped human hearts ex vivo. Coronary-perfused intact human atria (n = 7, 47 ± 14 y.o.; two female) were mapped using panoramic near-infrared optical mapping during pacing-induced AF. Then the hearts were imaged at approximately 170 µm3 resolution by 9.4 T gadolinium-enhanced MRI. The heart was segmented, and 3D AWT throughout atrial chambers was estimated and analysed. Optical mapping identified six localized RA re-entrant drivers in four hearts and four LA drivers in three hearts. All RA AF drivers were anchored to the pectinate muscle junctions with the crista terminalis or atrial walls. The four LA AF drivers were in the posterior LA. RA (n = 4) with AF drivers were thicker with greater AWT variation than RA (n = 3) without drivers (5.4 ± 2.6 mm versus 5.0 ± 2.4 mm, T-test p < 0.05; F-test p < 0.05). Furthermore, AWT in RA driver regions was thicker and varied more than in RA non-driver regions (5.1 ± 2.5 mm versus 4.4 ± 2.2 mm, T-test p < 0.05; F-test p < 0.05). On the other hand, LA (n = 3) with drivers was thinner than the LA (n = 4) without drivers. In particular, LA driver regions were thinner than the rest of LA regions (3.4 ± 1.0 mm versus 4.2 ± 1.0 mm, T-test p < 0.05). This study demonstrates chamber-specific AWT features of AF drivers. In RA, driver regions are thicker and have more variable AWT than non-driver regions. By contrast, LA drivers are thinner than non-drivers. Robust evaluation of patient-specific AWT features should be considered for chamber-specific targeted ablation.

Thermal effects on cardiac alternans onset and development: A spatiotemporal correlation analysis

Alternans of cardiac action potential duration represent critical precursors for the development of life-threatening arrhythmias and sudden cardiac death. The system's thermal state affects these electrical disorders requiring additional theoretical and experimental efforts to improve a patient-specific clinical understanding. In such a scenario, we generalize a recent work from Loppini et al. [Phys. Rev. E 100, 020201(R) (2019)PREHBM2470-004510.1103/PhysRevE.100.020201] by performing an extended spatiotemporal correlation study. We consider high-resolution optical mapping recordings of canine ventricular wedges' electrical activity at different temperatures and pacing frequencies. We aim to recommend the extracted characteristic length as a potential predictive index of cardiac alternans onset and evolution within a wide range of system states. In particular, we show that a reduction of temperature results in a drop of the characteristic length, confirming the impact of thermal instabilities on cardiac dynamics. Moreover, we theoretically investigate the use of such an index to identify and predict different alternans regimes. Finally, we propose a constitutive phenomenological law linking conduction velocity, characteristic length, and temperature in view of future numerical investigations.

Spatiotemporal correlation uncovers characteristic lengths in cardiac tissue

Complex spatiotemporal patterns of action potential duration have been shown to occur in many mammalian hearts due to period-doubling bifurcations that develop with increasing frequency of stimulation. Here, through high-resolution optical mapping experiments and mathematical modeling, we introduce a characteristic spatial length of cardiac activity in canine ventricular wedges via a spatiotemporal correlation analysis, at different stimulation frequencies and during fibrillation. We show that the characteristic length ranges from 40 to 20 cm during one-to-one responses and it decreases to a specific value of about 3 cm at the transition from period-doubling bifurcation to fibrillation. We further show that during fibrillation, the characteristic length is about 1 cm. Another significant outcome of our analysis is the finding of a constitutive phenomenological law obtained from a nonlinear fitting of experimental data which relates the conduction velocity restitution curve with the characteristic length of the system. The fractional exponent of 3/2 in our phenomenological law is in agreement with the domain size remapping required to reproduce experimental fibrillation dynamics within a realistic cardiac domain via accurate mathematical models.

Robust approach for rotor mapping in cardiac tissue

The motion of and interaction between phase singularities that lie at the centers of spiral waves capture many qualitative and, in some cases, quantitative features of complex dynamics in excitable systems. Being able to accurately reconstruct their position is thus quite important, even if the data are noisy and sparse, as in electrophysiology studies of cardiac arrhythmias, for instance. A recently proposed global topological approach [Marcotte and Grigoriev, Chaos 27, 093936 (2017)] promises to meaningfully improve the quality of the reconstruction compared with traditional, local approaches. Indeed, we found that this approach is capable of handling noise levels exceeding the range of the signal with minimal loss of accuracy. Moreover, it also works successfully with data sampled on sparse grids with spacing comparable to the mean separation between the phase singularities for complex patterns featuring multiple interacting spiral waves.

Response function framework for the dynamics of meandering or large-core spiral waves and modulated traveling waves

In many oscillatory or excitable systems, dynamical patterns emerge which are stationary or periodic in a moving frame of reference. Examples include traveling waves or spiral waves in chemical systems or cardiac tissue. We present a unified theoretical framework for the drift of such patterns under small external perturbations, in terms of overlap integrals between the perturbation and the adjoint critical eigenfunctions of the linearized operator (i.e., response functions). For spiral waves, the finite radius of the spiral tip trajectory and spiral wave meander are taken into account. Different coordinate systems can be chosen, depending on whether one wants to predict the motion of the spiral-wave tip, the time-averaged tip path, or the center of the meander flower. The framework is applied to analyze the drift of a meandering spiral wave in a constant external field in different regimes.

Muscle Thickness and Curvature Influence Atrial Conduction Velocities

Electroanatomical mapping is currently used to provide clinicians with information about the electrophysiological state of the heart and to guide interventions like ablation. These maps can be used to identify ectopic triggers of an arrhythmia such as atrial fibrillation (AF) or changes in the conduction velocity (CV) that have been associated with poor cell to cell coupling or fibrosis. Unfortunately, many factors are known to affect CV, including membrane excitability, pacing rate, wavefront curvature, and bath loading, making interpretation challenging. In this work, we show how endocardial conduction velocities are also affected by the geometrical factors of muscle thickness and wall curvature. Using an idealized three-dimensional strand, we show that transverse conductivities and boundary conditions can slow down or speed up signal propagation, depending on the curvature of the muscle tissue. In fact, a planar wavefront that is parallel to a straight line normal to the mid-surface does not remain normal to the mid-surface in a curved domain. We further demonstrate that the conclusions drawn from the idealized test case can be used to explain spatial changes in conduction velocities in a patient-specific reconstruction of the left atrial posterior wall. The simulations suggest that the widespread assumption of treating atrial muscle as a two-dimensional manifold for electrophysiological simulations will not accurately represent the endocardial conduction velocities in regions of the heart thicker than 0.5 mm with significant wall curvature.

Image-Based Computational Evaluation of the Effects of Atrial Wall Thickness and Fibrosis on Re-entrant Drivers for Atrial Fibrillation

Introduction: Catheter ablation (CA) is a common treatment for atrial fibrillation (AF), but the knowledge of optimal ablation sites, and hence clinical outcomes, are suboptimal. Increasing evidence suggest that ablation strategies based on patient-specific substrates information, such as distributions of fibrosis and atrial wall thickness (AWT), may be used to improve therapy. We hypothesized that competing influences of large AWT gradients and fibrotic patches on conductive properties of atrial tissue can determine locations of re-entrant drivers (RDs) sustaining AF. Methods: Two sets of models were used: (1) a simple model of 3D atrial tissue slab with a step change in AWT and a synthetic fibrosis patch, and (2) 3D models based on patient-specific right atrial (RA) and left atrial (LA) geometries. The latter were obtained from four healthy volunteers and two AF patients, respectively, using magnetic resonance imaging (MRI). A synthetic fibrotic patch was added in the RA and fibrosis distributions in the LA were obtained from gadolinium-enhanced MRI of the same patients. In all models, 3D geometry was combined with the Fenton-Karma atrial cell model to simulate RDs. Results: In the slab, RDs drifted toward, and then along the AWT step. However, with additional fibrosis, the RDs were localized in regions between the step and fibrosis. In the RA, RDs drifted toward and anchored to a large AWT gradient between the crista terminalis (CT) region and the surrounding atrial wall. Without such a gradient, RDs drifted toward the superior vena cava (SVC) or the tricuspid valve (TSV). With additional fibrosis, RDs initiated away from the CT anchored to the fibrotic patch, whereas RDs initiated close to the CT region remained localized between the two structures. In the LA, AWT was more uniform and RDs drifted toward the pulmonary veins (PVs). However, with additional fibrotic patches, RDs either anchored to them or multiplied. Conclusion: In the RA, RD locations are determined by both fibrosis and AWT gradients at the CT region. In the LA, they are determined by fibrosis due to the absence of large AWT gradients. These results elucidate mechanisms behind the stabilization of RDs sustaining AF and can help guide ablation therapy.

Filament Tension and Phase Locking of Meandering Scroll Waves

Meandering spiral waves are often observed in excitable media such as the Belousov-Zhabotinsky reaction and cardiac tissue. We derive a theory for drift dynamics of meandering rotors in general reaction-diffusion systems and apply it to two types of external disturbances: an external field and curvature-induced drift in three dimensions. We find two distinct regimes: with small filament curvature, meandering scroll waves exhibit filament tension, whose sign determines the stability and drift direction. In the regimes of strong external fields or meandering motion close to resonance, however, phase locking of the meander pattern is predicted and observed.

Nonlinear diffusion and thermo-electric coupling in a two-variable model of cardiac action potential

This work reports the results of the theoretical investigation of nonlinear dynamics and spiral wave breakup in a generalized two-variable model of cardiac action potential accounting for thermo-electric coupling and diffusion nonlinearities. As customary in excitable media, the common Q₁₀ and Moore factors are used to describe thermo-electric feedback in a 10° range. Motivated by the porous nature of the cardiac tissue, in this study we also propose a nonlinear Fickian flux formulated by Taylor expanding the voltage dependent diffusion coefficient up to quadratic terms. A fine tuning of the diffusive parameters is performed a priori to match the conduction velocity of the equivalent cable model. The resulting combined effects are then studied by numerically simulating different stimulation protocols on a one-dimensional cable. Model features are compared in terms of action potential morphology, restitution curves, frequency spectra, and spatio-temporal phase differences. Two-dimensional long-run simulations are finally performed to characterize spiral breakup during sustained fibrillation at different thermal states. Temperature and nonlinear diffusion effects are found to impact the repolarization phase of the action potential wave with non-monotone patterns and to increase the propensity of arrhythmogenesis.

Role of atrial wall thickness in wave-dynamics of atrial fibrillation

Background/Aims

Atrial anatomy and thickness may affect the electrical wave-dynamics of atrial fibrillation (AF). We explored the relationship between left atrial (LA) wall thickness (LAWT) or LA geometry and AF wave-dynamics.

Methods

We included 15 patients with persistent AF (age, 62.3 ± 11.9 years) who underwent AF catheter ablation. We measured the LAWT, LA endocardial curvature, and SD-curvature (surface bumpiness) from preprocedural computed tomography images. We compared those anatomical characteristics with electrophysiologic parameters such as dominant frequency (DF), Shannon entropy (ShEn), or complex fractionated atrial electrogram (CFAE)-cycle length (CL), calculated from intracardiac bipolar electrograms (300–500 points, 5 s), acquired during ablation procedures.

Results

1. LAWT (excluding fat) varied widely among patients, locations, and types of AF. LAWT was inversely correlated with LA volume (r = -0.565, p = 0.028) and positively correlated with SD-curvature (r = 0.272, p < 0.001). 2. LAWT was positively correlated with ShEn (r = 0.233, p < 0.001) and negatively correlated with CFAE-CL (r = -0.107, p = 0.038). 3. In the multivariate linear regression analyses for AF wave-dynamics parameters, DF (β = -0.29 [95% CI -0.44–-0.14], p < 0.001), ShEn (β = 0.19 [95% CI 0.12–0.25], p < 0.001), and CFAE-CL (β = 7.49 [95% CI 0.65–14.34], p = 0.032) were independently associated with LAWT.

Conclusion

Regional LAWT is associated with LA structural features, and has significant correlations with the wave-dynamics parameters associated with electrical wavebreaks or rotors in patients with persistent AF.

BeatBox-HPC simulation environment for biophysically and anatomically realistic cardiac electrophysiology

The BeatBox simulation environment combines flexible script language user interface with the robust computational tools, in order to setup cardiac electrophysiology in-silico experiments without re-coding at low-level, so that cell excitation, tissue/anatomy models, stimulation protocols may be included into a BeatBox script, and simulation run either sequentially or in parallel (MPI) without re-compilation. BeatBox is a free software written in C language to be run on a Unix-based platform. It provides the whole spectrum of multi scale tissue modelling from 0-dimensional individual cell simulation, 1-dimensional fibre, 2-dimensional sheet and 3-dimensional slab of tissue, up to anatomically realistic whole heart simulations, with run time measurements including cardiac re-entry tip/filament tracing, ECG, local/global samples of any variables, etc. BeatBox solvers, cell, and tissue/anatomy models repositories are extended via robust and flexible interfaces, thus providing an open framework for new developments in the field. In this paper we give an overview of the BeatBox current state, together with a description of the main computational methods and MPI parallelisation approaches.

Stability of patterns on thin curved surfaces

We consider reaction-diffusion equations on a thin curved surface and obtain a set of effective reaction-diffusion (R-D) equations to O(ε^{2}), where ε is the surface thickness. We observe that the R-D systems on these curved surfaces can have space-dependent reaction kinetics. Further, we use linear stability analysis to study the Schnakenberg model on spherical and cylindrical geometries. The dependence of the steady state on the thickness is determined for both cases, and we find that a change in the thickness can stabilize the unstable modes, and vice versa. The combined effect of thickness and curvature can play an important role in the rearrangement of spatial patterns on thin curved surfaces.

Adjoint eigenfunctions of temporally recurrent single-spiral solutions in a simple model of atrial fibrillation

This paper introduces a numerical method for computing the spectrum of adjoint (left) eigenfunctions of spiral wave solutions to reaction-diffusion systems in arbitrary geometries. The method is illustrated by computing over a hundred eigenfunctions associated with an unstable time-periodic single-spiral solution of the Karma model on a square domain. We show that all leading adjoint eigenfunctions are exponentially localized in the vicinity of the spiral tip, although the marginal modes (response functions) demonstrate the strongest localization. We also discuss the implications of the localization for the dynamics and control of unstable spiral waves. In particular, the interaction with no-flux boundaries leads to a drift of spiral waves which can be understood with the help of the response functions.

A novel technique to initiate and investigate scroll waves in thin layers of the photosensitive Belousov-Zhabotinsky reaction

While free scroll rings are non-stationary objects that either grow or contract with time, spatial confinement can have a large impact on their evolution reaching from significant lifetime extension (J.F. Totz, H. Engel, O. Steinbock, New J. Phys. 17, 093043 (2015)) up to formation of stable stationary and breathing pacemakers (A. Azhand, J.F. Totz, H. Engel, EPL 108, 10004 (2014)). Here, we explore the parameter range in which the interaction between an axis-symmetric scroll ring and a confining planar no-flux boundary can be studied experimentally in transparent gel layers supporting chemical wave propagation in the photosensitive variant of the Belousov-Zhabotinsky medium. Based on full three-dimensional simulations of the underlying modified complete Oregonator model for experimentally realistic parameters, we determine the conditions for successful initiation of scroll rings in a phase diagram spanned by the layer thickness and the applied light intensity. Furthermore, we discuss whether the illumination-induced excitability gradient due to Lambert-Beer's law as well as a possible inclination of the filament plane with respect to the no-flux boundary can destabilize the scroll ring.

Drift of Scroll Wave Filaments in an Anisotropic Model of the Left Ventricle of the Human Heart

Scroll waves are three-dimensional vortices which occur in excitable media. Their formation in the heart results in the onset of cardiac arrhythmias, and the dynamics of their filaments determine the arrhythmia type. Most studies of filament dynamics were performed in domains with simple geometries and generic description of the anisotropy of cardiac tissue. Recently, we developed an analytical model of fibre structure and anatomy of the left ventricle (LV) of the human heart. Here, we perform a systematic study of the dynamics of scroll wave filaments for the cases of positive and negative tension in this anatomical model. We study the various possible shapes of LV and different degree of anisotropy of cardiac tissue. We show that, for positive filament tension, the final position of scroll wave filament is mainly determined by the thickness of the myocardial wall but, however, anisotropy attracts the filament to the LV apex. For negative filament tension, the filament buckles, and for most cases, tends to the apex of the heart with no or slight dependency on the thickness of the LV. We discuss the mechanisms of the observed phenomena and their implications for cardiac arrhythmias.

Scroll wave drift along steps, troughs, and corners

Three-dimensional excitable systems can create nonlinear scroll waves that rotate around one-dimensional phase singularities. Recent theoretical work predicts that these filaments drift along step-like height variations. Here, we test this prediction using experiments with thin layers of the Belousov-Zhabotinsky reaction. We observe that over short distances scroll waves are attracted towards the step and then rapidly commence a steady drift along the step line. The translating filaments always reside on the shallow side of the step near the edge. Accordingly, filaments in the deep domain initially collide with and shorten at the step wall. The drift speeds obey the predicted proportional dependence on the logarithm of the height ratio and the direction depends on the vortex chirality. We also observe drift along the perimeter of rectangular plateaus and find that the filaments perform sharp turns at the corners. In addition, we investigate rectangular troughs for which vortices of equal chirality can drift in different directions. The latter two effects are reproduced in numerical simulations with the Barkley model. The simulations show that narrow troughs instigate scroll wave encounters that induce repulsive interaction and symmetry breaking. Similar phenomena could exist in the geometrically complicated ventricles of the human heart where reentrant vortex waves cause tachycardia and fibrillation.