Emergence of spiral wave activity in a mechanically heterogeneous reaction-diffusion-mechanics system

A survey of pathways for mechano-electric coupling in the atria

Mechano-electric coupling (MEC) in atrial tissue has received sparse investigation to date, despite the well-known association between chronic atrial dilation and atrial fibrillation (AF). Of note, no fewer than six different mechanisms pertaining to stretch-activated channels, cellular capacitance and geometric effects have been identified in the literature as potential players. In this mini review, we briefly survey each of these pathways to MEC. We then perform computational simulations using single cell and tissue models in presence of various stretch regimes and MEC pathways. This allows us to assess the relative significance of each pathway in determining action potential duration, conduction velocity and rotor stability. For chronic atrial stretch, we find that stretch-induced alterations in membrane capacitance decrease conduction velocity and increase action potential duration, in agreement with experimental findings. In the presence of time-dependent passive atrial stretch, stretch-activated channels play the largest role, leading to after-depolarizations and rotor hypermeandering. These findings suggest that physiological atrial stretches, such as passive stretch during the atrial reservoir phase, may play an important part in the mechanisms of atrial arrhythmogenesis.

•

Passive strains caused by ventricular contraction need to be considered when incorporating mechano-electro feedback in atrial electrophysiology models.

•

In chronic stretch, stretch-induced capacitance changes dominate.

•

Chronic stretch leads to an increase in action potential duration and a reduction in conduction velocity, consistent with experimental studies.

•

In the presence of passive stretch, stretch-activated channels can induce delayed after-depolarisations and lead to rotor hypermeandering.

•

Mechano-electro feedback is thus likely to have implications for the genesis and maintenance of atrial arrhythmias.

Synchronization of viscoelastically coupled excitable oscillators

Viscoelastically coupled excitable oscillators are used to model individually beating spatially separated cardiomyocytes surrounded by an extracellular matrix (ECM). We investigate how mechanical coupling via the ECM can synchronize two such oscillators with excitation contraction coupling and electromechanical feedback and how this synchronization depends on the rheological properties of the ECM. Extending our study to a linear chain of coupled oscillators we find a transition to synchronization as the ECM becomes stiffer. In the case of purely elastic coupling we observe antiphase chimera states.

Uniaxial strain of cultured mouse and rat cardiomyocyte strands slows conduction more when its axis is parallel to impulse propagation than when it is perpendicular

AIM

Cardiac tissue deformation can modify tissue resistance, membrane capacitance and ion currents and hence cause arrhythmogenic slow conduction. Our aim was to investigate whether uniaxial strain causes different changes in conduction velocity (θ) when the principal strain axis is parallel vs perpendicular to impulse propagation.

METHODS

Cardiomyocyte strands were cultured on stretchable custom microelectrode arrays, and θ was determined during steady-state pacing. Uniaxial strain (5%) with principal axis parallel (orthodromic) or perpendicular (paradromic) to propagation was applied for 1 minute and controlled by imaging a grid of markers. The results were analysed in terms of cable theory.

RESULTS

Both types of strain induced immediate changes of θ upon application and release. In material coordinates, orthodromic strain decreased θ significantly more (P < .001) than paradromic strain (2.2 ± 0.5% vs 1.0 ± 0.2% in n = 8 mouse cardiomyocyte cultures, 2.3 ± 0.4% vs 0.9 ± 0.5% in n = 4 rat cardiomyocyte cultures, respectively). The larger effect of orthodromic strain can be explained by the increase in axial myoplasmic resistance, which is not altered by paradromic strain. Thus, changes in tissue resistance substantially contributed to the changes of θ during strain, in addition to other influences (eg stretch-activated channels). Besides these immediate effects, the application of strain also consistently initiated a slow progressive decrease in θ and a slow recovery of θ upon release.

CONCLUSION

Changes in cardiac conduction velocity caused by acute stretch do not only depend on the magnitude of strain but also on its orientation relative to impulse propagation. This dependence is due to different effects on tissue resistance.

Fast propagation regions cause self-sustained reentry in excitable media

Self-sustained waves of electrophysiological activity can cause arrhythmia in the heart. These reentrant excitations have been associated with spiral waves circulating around either an anatomically defined weakly conducting region or a functionally determined core. Recently, an ablation procedure has been clinically introduced that stops atrial fibrillation of the heart by destroying the electrical activity at the spiral core. This is puzzling because the tissue at the anatomically defined spiral core would already be weakly conducting, and a further decrease should not improve the situation. In the case of a functionally determined core, an ablation procedure should even further stabilize the rotating wave. The efficacy of the procedure thus needs explanation. Here, we show theoretically that fundamentally in any excitable medium a region with a propagation velocity faster than its surrounding can act as a nucleation center for reentry and can anchor an induced spiral wave. Our findings demonstrate a mechanistic underpinning for the recently developed ablation procedure. Our theoretical results are based on a very general and widely used two-component model of an excitable medium. Moreover, the important control parameters used to realize conditions for the discovered phenomena are applicable to quite different multicomponent models.

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that is triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media with applications to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact for cardiac arrhythmias.

Transient and periodic spatiotemporal structures in a reaction-diffusion-mechanics system

We study transient spatiotemporal structures induced by a weak space-time localized stimulus in an excitable contractile fiber within a two-component globally coupled reaction-diffusion model. The model which we develop allows us to analyze various regimes of excitation spreading and determine origin of the induced structures for various contraction types (defined by the fiber fixation) and global coupling strengths. One of the most notable effects we observed is the after-excitation effect. It leads to emergence of multiple excitation pulses excited by a single external stimulus and can result in long-lasting transient activity and appearance of new oscillatory attractor regimes, including the ones with multiple phase clusters.

Chiral selection and frequency response of spiral waves in reaction-diffusion systems under a chiral electric field

Chirality is one of the most fundamental properties of many physical, chemical, and biological systems. However, the mechanisms underlying the onset and control of chiral symmetry are largely understudied. We investigate possibility of chirality control in a chemical excitable system (the Belousov-Zhabotinsky reaction) by application of a chiral (rotating) electric field using the Oregonator model. We find that unlike previous findings, we can achieve the chirality control not only in the field rotation direction, but also opposite to it, depending on the field rotation frequency. To unravel the mechanism, we further develop a comprehensive theory of frequency synchronization based on the response function approach. We find that this problem can be described by the Adler equation and show phase-locking phenomena, known as the Arnold tongue. Our theoretical predictions are in good quantitative agreement with the numerical simulations and provide a solid basis for chirality control in excitable media.

Spatial symmetry breaking determines spiral wave chirality

Chirality represents a fundamental property of spiral waves. Introducing obstacles into cardiac monolayers leads to the initiation of clockwise-rotating, counterclockwise-rotating, and pairs of spiral waves. Simulations show that the precise location of the obstacle and the pacing frequency determine spiral wave chirality. Instabilities predicted by curves relating the action potential duration and the pacing frequency at different spatial locations predict sites of wave break initiation and, hence, spiral wave chirality.

Influences of periodic mechanical deformation on pinned spiral waves

In a generic model of excitable media, we study the behavior of spiral waves interacting with obstacles and their dynamics under the influences of simple periodic mechanical deformation (PMD). Depending on the characteristics of the obstacles, i.e., size and excitability, the rotation of a pinned spiral wave shows different scenarios, e.g., embedding into or anchoring on an obstacle. Three different drift phenomena induced by PMD are observed: scattering on small partial-excitable obstacles, meander-induced unpinning on big partial-excitable obstacles, and drifting around small unexcitable obstacles. Their underlying mechanisms are discussed. The dependence of the threshold amplitude of PMD on the characteristics of the obstacles to successfully remove pinned spiral waves on big partial-excitable obstacles is studied.

Electric-field-controlled unpinning of scroll waves

Three-dimensional excitation vortices exist in systems such as chemical reactions and the human heart. Their one-dimensional rotation backbone can pin to unexcitable heterogeneities, which greatly affect the structure, dynamics, and lifetime of the vortex. In experiments with the Belousov-Zhabotinsky reaction, we demonstrate vortex unpinning from a pair of inert and impermeable spheres using externally applied electric fields. Unpinning occurs abruptly but is preceded by a slow reorientation and deformation of the initially circular vortex loop. Our experimental findings are reproduced by numerical simulations of an excitable reaction-diffusion-advection model.

A discrete electromechanical model for human cardiac tissue: effects of stretch-activated currents and stretch conditions on restitution properties and spiral wave dynamics

We introduce an electromechanical model for human cardiac tissue which couples a biophysical model of cardiac excitation (Tusscher, Noble, Noble, Panfilov, 2006) and tension development (adjusted Niederer, Hunter, Smith, 2006 model) with a discrete elastic mass-lattice model. The equations for the excitation processes are solved with a finite difference approach, and the equations of the mass-lattice model are solved using Verlet integration. This allows the coupled problem to be solved with high numerical resolution. Passive mechanical properties of the mass-lattice model are described by a generalized Hooke's law for finite deformations (Seth material). Active mechanical contraction is initiated by changes of the intracellular calcium concentration, which is a variable of the electrical model. Mechanical deformation feeds back on the electrophysiology via stretch-activated ion channels whose conductivity is controlled by the local stretch of the medium. We apply the model to study how stretch-activated currents affect the action potential shape, restitution properties, and dynamics of spiral waves, under constant stretch, and dynamic stretch caused by active mechanical contraction. We find that stretch conditions substantially affect these properties via stretch-activated currents. In constantly stretched medium, we observe a substantial decrease in conduction velocity, and an increase of action potential duration; whereas, with dynamic stretch, action potential duration is increased only slightly, and the conduction velocity restitution curve becomes biphasic. Moreover, in constantly stretched medium, we find an increase of the core size and period of a spiral wave, but no change in rotation dynamics; in contrast, in the dynamically stretching medium, we observe spiral drift. Our results may be important to understand how altered stretch conditions affect the heart's functioning.

Pacemaker interactions induce reentrant wave dynamics in engineered cardiac culture

Pacemaker interactions can lead to complex wave dynamics seen in certain types of cardiac arrhythmias. We use experimental and mathematical models of pacemakers in heterogeneous excitable media to investigate how pacemaker interactions can be a mechanism for wave break and reentrant wave dynamics. Embryonic chick ventricular cells are cultured in vitro so as to create a dominant central pacemaker site that entrains other pacemakers in the medium. Exposure of those cultures to a potassium channel blocker, E-4031, leads to emergence of peripheral pacemakers that compete with each other and with the central pacemaker. Waves emitted by faster pacemakers break up over the slower pacemaker to form reentrant waves. Similar dynamics are observed in a modified FitzHugh-Nagumo model of heterogeneous excitable media with two distinct sites of pacemaking. These findings elucidate a mechanism of pacemaker-induced reentry in excitable media.