Interaction between spiral and paced waves in cardiac tissue

For prevention of lethal arrhythmias, patients at risk receive implantable cardioverter-defibrillators, which use high-frequency antitachycardia pacing (ATP) to convert tachycardias to a normal rhythm. One of the suggested ATP mechanisms involves paced-induced drift of rotating waves followed by their collision with the boundary of excitable tissue. This study provides direct experimental evidence of this mechanism. In monolayers of neonatal rat cardiomyocytes in which rotating waves of activity were initiated by premature stimuli, we used the Ca(2+)-sensitive indicator fluo 4 to observe propagating wave patterns. The interaction of the spiral tip with a paced wave was then monitored at a high spatial resolution. In the course of the experiments, we observed spiral wave pinning to local heterogeneities within the myocyte layer. High-frequency pacing led, in a majority of cases, to successful termination of spiral activity. Our data show that 1) stable spiral waves in cardiac monolayers tend to be pinned to local heterogeneities or areas of altered conduction, 2) overdrive pacing can shift a rotating wave from its original site, and 3) the wave break, formed as a result of interaction between the spiral tip and a paced wave front, moves by a paced-induced drift mechanism to an area where it may become unstable or collide with a boundary. The data were complemented by numerical simulations, which was used to further analyze experimentally observed behavior.     

Dynamics of Scroll Wave in a Three-Dimensional System with Changing Gradient

The dynamics of a scroll wave in an excitable medium with gradient excitability is studied in detail. Three parameter regimes can be distinguished by the degree of gradient. For a small gradient, the system reaches a simple rotating synchronization. In this regime, the rigid rotating velocity of spiral waves is maximal in the layers with the highest filament twist. As the excitability gradient increases, the scroll wave evolutes into a meandering synchronous state. This transition is accompanied by a variation in twisting rate. Filament twisting may prevent the breakup of spiral waves in the bottom layers with a low excitability with which a spiral breaks in a 2D medium. When the gradient is large enough, the twisted filament breaks up, which results in a semi-turbulent state where the lower part is turbulent while the upper part contains a scroll wave with a low twisting filament.     

Interaction of spiral and flat periodic autowaves in an active medium

Interaction between the rotating wave and a periodic external source in the model of Fitz Hugh--Nagumo type was computed. When the periods of the external source are longer than the rotation period of the spiral wave (T greater than Ts) the external source does not affect the spiral wave. At T less than Ts autowave synchronization effects are observed. The oscillation period predetermined by the external source is set in all the points of the medium except the neighbourhood of the spiral wave. The dislocation (wavebreak) persists in the medium drifting slowly at the angle to the wave vector of the flat waves. After the external source is eliminated, the spiral wave with the original period restores from this dislocation. When the dislocation reaches the interface, it disappears. In this case after the switching off of the external source the resting state is established. A theory of the drift is proposed which connects the drift velocity with the nucleus size and the rotation period of spiral wave.     

Spiral waves in two-dimensional models of ventricular muscle: formation of a stationary core.

Previous experimental studies have clearly demonstrated the existence of drifting and stationary electrical spiral waves in cardiac muscle and their involvement in cardiac arrhythmias. Here we present results of a study of reentrant excitation in computer simulations based on a membrane model of the ventricular cell. We have explored in detail the parameter space of the model, using tools derived from previous numerical studies in excitation-dynamics models. We have found appropriate parametric conditions for sustained stable spiral wave dynamics (1 s of activity or approximately 10 rotations) in simulations of an anisotropic (ratio in velocity 4:1) cardiac sheet of 2 cm x 2 cm. Initially, we used a model that reproduced well the characteristics of planar electrical waves exhibited by thin sheets of sheep ventricular epicardial muscle during rapid pacing at a cycle length of 300 ms. Under these conditions, the refractory period was 147 ms; the action potential duration (APD) was 120 ms; the propagation velocity along fibers was 33 cm/s; and the wavelength along fibers was 4.85 cm. Using cross-field stimulation in this model, we obtained a stable self-sustaining spiral wave rotating around an unexcited core of 1.75 mm x 7 mm at a period of 115 ms, which reproduced well the experimental results. Thus the data demonstrate that stable spiral wave activity can occur in small cardiac sheets whose wavelength during planar wave excitation in the longitudinal direction is larger than the size of the sheet. Analysis of the mechanism of this observation demonstrates that, during rotating activity, the core exerts a strong electrotonic influence that effectively abbreviates APD (and thus wavelength) in its immediate surroundings and is responsible for the stabilization and perpetuation of the activity. We conclude that appropriate adjustments in the kinetics of the activation front (i.e., threshold for activation and upstroke velocity of the initiating beat) of currently available models of the cardiac cell allow accurate reproduction of experimentally observed self-sustaining spiral wave activity. As such, the results set the stage for an understanding of functional reentry in terms of ionic mechanisms.     

Selection of Multiarmed Spiral Waves in a Regular Network of Neurons

Formation and selection of multiarmed spiral wave due to spontaneous symmetry breaking are investigated in a regular network of Hodgkin-Huxley neuron by changing the excitability and imposing spatial forcing currents on the neurons in the network. The arm number of the multiarmed spiral wave is dependent on the distribution of spatial forcing currents and excitability diversity in the network, and the selection criterion for supporting multiarmed spiral waves is discussed. A broken spiral segment is measured by a short polygonal line connected by three adjacent points (controlled nodes), and a double-spiral wave can be developed from the spiral segment. Multiarmed spiral wave is formed when a group of double-spiral waves rotate in the same direction in the network. In the numerical studies, a group of controlled nodes are selected and spatial forcing currents are imposed on these nodes, and our results show that l-arm stable spiral wave (l = 2, 3, 4,...8) can be induced to occupy the network completely. It is also confirmed that low excitability is critical to induce multiarmed spiral waves while high excitability is important to propagate the multiarmed spiral wave outside so that distinct multiarmed spiral wave can occupy the network completely. Our results confirm that symmetry breaking of target wave in the media accounts for emergence of multiarmed spiral wave, which can be developed from a group of spiral waves with single arm under appropriate condition, thus the potential formation mechanism of multiarmed spiral wave in the media is explained.     

Region [corrected] of slowed conduction acts as core for spiral wave reentry in cardiac cell monolayers

Pathophysiological heterogeneity in cardiac tissue is related to the occurrence of arrhythmias. Of importance are regions of slowed conduction, which have been implicated in the formation of conduction block and reentry. Experimentally, it has been a challenge to produce local heterogeneity in a manner that is both reversible and well controlled. Consequently, we developed a dual-zone superfusion chamber that can dynamically create a small (5 mm) central island of heterogeneity in cultured cardiac cell monolayers. Three different conditions were studied to explore the effect of regionally slowed conduction on wave propagation and reentry: depolarization by elevated extracellular potassium, sodium channel inhibition with lidocaine, and cell-cell decoupling with palmitoleic acid. Using optical mapping of transmembrane voltage, we found that the central region of slowed conduction always served as the core region around which a spiral wave formed and then revolved following a period of rapid pacing. Because of the localized slowing in the core region, we observed experimentally for the first time an S shape of the spiral wave front near its tip. These results indicate that a small region of slowed conduction can play a crucial role in the formation, anchoring, and modulation of reentrant spiral waves.     

A nonlinear single-mode structure in a plasma with anisotropic temperature

An exact solution is derived to the equations of vortex electron anisotropic hydrodynamics for a plasma that is unstable against the Weibel instability driven by the electron temperature anisotropy. This solution describes saturation of the Weibel instability in the single-mode regime with an arbitrary wavelength and corresponds to a standing helical wave of magnetic perturbations in which the amplitude of the generated magnetic field varies periodically over time. The longitudinal and transverse (with respect to the rotating anisotropy axis) plasma temperatures are subject to the same periodic variations. In this case, the maximum magnetic field energy can be on the order of the plasma thermal energy.

     

Ring-shaped assembly of microtubules shows preferential counterclockwise motion

In this paper, we reveal that microtubules (MTs), reconstructed from tubulin in vitro in the presence of guanosine-5'-triphosphate (GTP), have a ring or spiral shape on a motor protein-fixed surface, and these MTs show biased motion in the counterclockwise direction. By cross-linking these MTs during the sliding motion, we obtained large ring-shaped MT assemblies, 1 approximately 12.6 microm in diameter. The ratio of the rings rotating in the counterclockwise direction to those rotating in the clockwise direction was approximately 3/1. Under optimized conditions, the ratio was as high as 14/1. Thus, we successfully obtained aggregated MTs with a large hierarchic structure that shows a preferential motion, through a dynamic process in vitro.     

Inward Rotating Spiral Waves in Glycolysis

We report on the first observation of inward rotating spiral waves (antispirals) in a biochemical reaction-diffusion system. Experiments are performed with extracts from yeast cells in an open spatial reactor. By increasing the protein concentration of the extract we observe a transition from outward to inward propagating waves of glycolytic activity. Numerical simulations with an allosteric model for the phosphofructokinase can reproduce these inward propagating waves over a wide range of parameters if the octameric structure of yeast phosphofructokinase is taken into account.     

A discrete cell model with adaptive signalling for aggregation of Dictyostelium discoideum.

Dictyostelium discoideum (Dd) is a widely studied model system from which fundamental insights into cell movement, chemotaxis, aggregation and pattern formation can be gained. In this system aggregation results from the chemotactic response by dispersed amoebae to a travelling wave of the chemoattractant cAMP. We have developed a model in which the cells are treated as discrete points in a continuum field of the chemoattractant, and transduction of the extracellular cAMP signal into the intracellular signal is based on the G protein model developed by Tang & Othmer. The model reproduces a number of experimental observations and gives further insight into the aggregation process. We investigate different rules for cell movement the factors that influence stream formation the effect on aggregation of noise in the choice of the direction of movement and when spiral waves of chemoattractant and cell density are likely to occur. Our results give new insight into the origin of spiral waves and suggest that streaming is due to a finite amplitude instability.     

Spatially structured oscillations in a two-dimensional excitatory neuronal network with synaptic depression

We study the spatiotemporal dynamics of a two-dimensional excitatory neuronal network with synaptic depression. Coupling between populations of neurons is taken to be nonlocal, while depression is taken to be local and presynaptic. We show that the network supports a wide range of spatially structured oscillations, which are suggestive of phenomena seen in cortical slice experiments and in vivo. The particular form of the oscillations depends on initial conditions and the level of background noise. Given an initial, spatially localized stimulus, activity evolves to a spatially localized oscillating core that periodically emits target waves. Low levels of noise can spontaneously generate several pockets of oscillatory activity that interact via their target patterns. Periodic activity in space can also organize into spiral waves, provided that there is some source of rotational symmetry breaking due to external stimuli or noise. In the high gain limit, no oscillatory behavior exists, but a transient stimulus can lead to a single, outward propagating target wave.     

Spiral Wave Propagation in Communities with Spatially Correlated Heterogeneity

Many multicellular communities propagate signals in a directed manner via excitable waves. Cell-to-cell heterogeneity is a ubiquitous feature of multicellular communities, but the effects of heterogeneity on wave propagation are still unclear. Here, we use a minimal FitzHugh-Nagumo-type model to investigate excitable wave propagation in a two-dimensional heterogeneous community. The model shows three dynamic regimes in which waves either propagate directionally, die out, or spiral indefinitely, and we characterize how these regimes depend on the heterogeneity parameters. We find that in some parameter regimes, spatial correlations in the heterogeneity enhance directional propagation and suppress spiraling. However, in other regimes, spatial correlations promote spiraling, a surprising feature that we explain by demonstrating that these spirals form by a second, distinct mechanism. Finally, we characterize the dynamics using techniques from percolation theory. Despite the fact that percolation theory does not completely describe the dynamics quantitatively because it neglects the details of the excitable propagation, we find that it accounts for the transitions between the dynamic regimes and the general dependency of the spiral period on the heterogeneity and thus provides important insights. Our results reveal that the spatial structure of cell-to-cell heterogeneity can have important consequences for signal propagation in cellular communities.     

Learning to use an invisible visual signal for perception

How does the brain construct a percept from sensory signals? One approach to this fundamental question is to investigate perceptual learning as induced by exposure to statistical regularities in sensory signals [1-7]. Recent studies showed that exposure to novel correlations between sensory signals can cause a signal to have new perceptual effects [2, 3]. In those studies, however, the signals were clearly visible. The automaticity of the learning was therefore difficult to determine. Here we investigate whether learning of this sort, which causes new effects on appearance, can be low level and automatic by employing a visual signal whose perceptual consequences were made invisible-a vertical disparity gradient masked by other depth cues. This approach excluded high-level influences such as attention or consciousness. Our stimulus for probing perceptual appearance was a rotating cylinder. During exposure, we introduced a new contingency between the invisible signal and the rotation direction of the cylinder. When subsequently presenting an ambiguously rotating version of the cylinder, we found that the invisible signal influenced the perceived rotation direction. This demonstrates that perception can rapidly undergo "structure learning" by automatically picking up novel contingencies between sensory signals, thus automatically recruiting signals for novel uses during the construction of a percept.     

Attraction and repulsion of spiral waves by inhomogeneity of conduction anisotropy—a model of spiral wave interaction with electrical remodeling of heart tissue

Various forms of heart disease are associated with remodeling of the heart muscle, which results in a perturbation of cell-to-cell electrical coupling. These perturbations may alter the trajectory of spiral wave drift in the heart muscle. We investigate the effect of spatially extended inhomogeneity of transverse cell coupling on the spiral wave trajectory using a simple active media model. The spiral wave was either attracted or repelled from the center of inhomogeneity as a function of cell excitability and gradient of the cell coupling. High levels of excitability resulted in an attraction of the wave to the center of inhomogeneity, whereas low levels resulted in an escape and termination of the spiral wave. The spiral wave drift velocity was related to the gradient of the coupling and the initial position of the wave. In a diseased heart, a region of altered transverse coupling corresponds with local gap junction remodeling that may be responsible for stabilization-destabilization of spiral waves and hence reflect potentially important targets in the treatment of heart arrhythmias.     

The mathematics of gnomonic seashells

Parametric Cartesian vector-valued functions are constructed for the purpose of systematically describing various features of spiral shell geometry. The underlying geometrical hypothesis is that molluscan shell shapes can usually, to at least a good first approximation, be developed by rotating a generating curve about a fixed axis whilst simultaneously diminishing it by an “equiangular spiral” scale factor. A first-order symmetry equation is derived; then variational calculus is used to construct energy functionals which establish that Hooke's law is inherent in the formalism and that naturally occurring shell geometries are analogous to those of elastic spiral “clock springs”. The biological requirement that shelly structures must exist in a three-dimensional space is shown to be a sufficiently powerful mathematical constraint to ensure the existence of geometrical artifacts which can, perhaps, be likened to the conservation laws, pseudoforces, and fields of classical physics.     

Electromechanical model of excitable tissue to study reentrant cardiac arrhythmias

We introduce the concept of a contracting excitable medium that is capable of conducting non-linear waves of excitation that in turn initiate contraction. Furthermore, these kinematic deformations have a feedback effect on the excitation properties of the medium. Electrical characteristics resemble basic models of cardiac excitation that have been used to successfully study mechanisms of reentrant cardiac arrhythmias in electrophysiology. We present a computational framework that employs electromechanical and mechanoelectric feedback to couple a three-variable FitzHugh–Nagumo-type excitation-tension model to the non-linear stress equilibrium equations, which govern large deformation hyperelasticity. Numerically, the coupled electromechanical model combines a finite difference method approach to integrate the excitation equations, with a Galerkin finite element method to solve the equations governing tissue mechanics. We present example computations demonstrating various effects of contraction on stationary rotating spiral waves and spiral wave break. We show that tissue mechanics significantly contributes to the dynamics of electrical propagation, and that a coupled electromechanical approach should be pursued in future electrophysiological modelling studies.     

Generation and escape of local waves from the boundary of uncoupled cardiac tissue

We aim to understand the formation of abnormal waves of activity from myocardial regions with diminished cell-to-cell coupling. En route to this goal, we studied the behavior of a heterogeneous myocyte network in which a sharp coupling gradient was placed under conditions of increasing network automaticity. Experiments were conducted in monolayers of neonatal rat cardiomyocytes using heptanol and isoproterenol as means of altering cell-to-cell coupling and automaticity, respectively. Experimental findings were explained and expanded using a modified Beeler-Reuter numerical model. The data suggest that the combination of a heterogeneous substrate, a gradient of coupling, and an increase in oscillatory activity of individual cells creates a rich set of behaviors associated with self-generated spiral waves and ectopic sources. Spiral waves feature a flattened shape and a pin-unpin drift type of tip motion. These intercellular waves are action-potential based and can be visualized with either voltage or calcium transient measurements. A source/load mismatch on the interface between the boundary and well-coupled layers can lock wavefronts emanating from both ectopic sources and rotating waves within the inner layers of the coupling gradient. A numerical approach allowed us to explore how 1), the spatial distribution of cells, 2), the amplitude and dispersion of cell automaticity, and 3), the speed at which the coupling gradient moves in space affect wave behavior, including its escape into well-coupled tissue.     

Resonant Drift of Spiral Waves in the Complex Ginzburg-Landau Equation

Weak periodic external perturbations of an autowave medium can cause large-distance directed motion of the spiral wave. This happens when the period of the perturbation coincides with, or is close to the rotation period of a spiral wave, or its multiple. Such motion is called resonant or parametric drift. It may be used for low-voltage defibrillation of heart tissue. Theory of the resonant drift exists, but so far was used only qualitatively. In this paper, we show good quantitative agreement of the theory with direct numerical simulations. This is done for Complex Ginzburg-Landau Equation, one of the simplest autowave models.