Learning poly-synaptic paths with traveling waves

Traveling waves are commonly observed across the brain. While previous studies have suggested the role of traveling waves in learning, the mechanism remains unclear. We adopted a computational approach to investigate the effect of traveling waves on synaptic plasticity. Our results indicate that traveling waves facilitate the learning of poly-synaptic network paths when combined with a reward-dependent local synaptic plasticity rule. We also demonstrate that traveling waves expedite finding the shortest paths and learning nonlinear input/output mapping, such as exclusive or (XOR) function.     

Effect of Sharp Jumps at the Edges of Phase Response Curves on Synchronization of Electrically Coupled Neuronal Oscillators

We study synchronization phenomenon of coupled neuronal oscillators using the theory of weakly coupled oscillators. The role of sudden jumps in the phase response curve profiles found in some experimental recordings and models on the ability of coupled neurons to exhibit synchronous and antisynchronous behavior is investigated, when the coupling between the neurons is electrical. The level of jumps in the phase response curve at either end, spike width and frequency of voltage time course of the coupled neurons are parameterized using piecewise linear functional forms, and the conditions for stable synchrony and stable antisynchrony in terms of those parameters are computed analytically. The role of the peak position of the phase response curve on phase-locking is also investigated.     

Transitory behaviors in diffusively coupled nonlinear oscillators

We study collective behaviors of diffusively coupled oscillators which exhibit out-of-phase synchrony for the case of weakly interacting two oscillators. In large populations of such oscillators interacting via one-dimensionally nearest neighbor couplings, there appear various collective behaviors depending on the coupling strength, regardless of the number of oscillators. Among others, we focus on an intermittent behavior consisting of the all-synchronized state, a weakly chaotic state and some sorts of metachronal waves. Here, a metachronal wave means a wave with orderly phase shifts of oscillations. Such phase shifts are produced by the dephasing interaction which produces the out-of-phase synchronized states in two coupled oscillators. We also show that the abovementioned intermittent behavior can be interpreted as in-out intermittency where two saddles on an invariant subspace, the all-synchronized state and one of the metachronal waves play an important role.     

Three people can synchronize as coupled oscillators during sports activities

We experimentally investigated the synchronized patterns of three people during sports activities and found that the activity corresponded to spatiotemporal patterns in rings of coupled biological oscillators derived from symmetric Hopf bifurcation theory, which is based on group theory. This theory can provide catalogs of possible generic spatiotemporal patterns irrespective of their internal models. Instead, they are simply based on the geometrical symmetries of the systems. We predicted the synchronization patterns of rings of three coupled oscillators as trajectories on the phase plane. The interactions among three people during a 3 vs. 1 ball possession task were plotted on the phase plane. We then demonstrated that two patterns conformed to two of the three patterns predicted by the theory. One of these patterns was a rotation pattern (R) in which phase differences between adjacent oscillators were almost 2π/3. The other was a partial anti-phase pattern (PA) in which the two oscillators were anti-phase and the third oscillator frequency was dead. These results suggested that symmetric Hopf bifurcation theory could be used to understand synchronization phenomena among three people who communicate via perceptual information, not just physically connected systems such as slime molds, chemical reactions, and animal gaits. In addition, the skill level in human synchronization may play the role of the bifurcation parameter.     

The spatial origin of a perceptual transition in binocular rivalry

When the left and the right eye are simultaneously presented with incompatible images at overlapping retinal locations, an observer typically reports perceiving only one of the two images at a time. This phenomenon is called binocular rivalry. Perception during binocular rivalry is not stable; one of the images is perceptually dominant for a certain duration (typically in the order of a few seconds) after which perception switches towards the other image. This alternation between perceptual dominance and suppression will continue for as long the images are presented. A characteristic of binocular rivalry is that a perceptual transition from one image to the other generally occurs in a gradual manner: the image that was temporarily suppressed will regain perceptual dominance at isolated locations within the perceived image, after which its visibility spreads throughout the whole image. These gradual transitions from perceptual suppression to perceptual dominance have been labeled as traveling waves of perceptual dominance. In this study we investigate whether stimulus parameters affect the location at which a traveling wave starts. We varied the contrast, spatial frequency or motion speed in one of the rivaling images, while keeping the same parameter constant in the other image. We used a flash-suppression paradigm to force one of the rival images into perceptual suppression. Observers waited until the suppressed image became perceptually dominant again, and indicated the position at which this breakthrough from suppression occurred. Our results show that the starting point of a traveling wave during binocular rivalry is highly dependent on local stimulus parameters. More specifically, a traveling wave most likely started at the location where the contrast of the suppressed image was higher than that of the dominant one, the spatial frequency of the suppressed image was lower than that of the dominant one, and the motion speed of the suppressed image was higher than that of the dominant one. We suggest that a breakthrough from suppression to dominance occurs at the location where salience (the degree to which a stimulus element stands out relative to neighboring elements) of the suppressed image is higher than that of the dominant one. Our results further show that stimulus parameters affecting the temporal dynamics during continuous viewing of rival images described in other studies, also affect the spatial origin of traveling waves during binocular rivalry.     

Unfolding an electronic integrate-and-fire circuit

Many physical and biological phenomena involve accumulation and discharge processes that can occur on significantly different time scales. Models of these processes have contributed to understand excitability self-sustained oscillations and synchronization in arrays of oscillators. Integrate-and-fire (I+F) models are popular minimal fill-and-flush mathematical models. They are used in neuroscience to study spiking and phase locking in single neuron membranes, large scale neural networks, and in a variety of applications in physics and electrical engineering. We show here how the classical first-order I+F model fits into the theory of nonlinear oscillators of van der Pol type by demonstrating that a particular second-order oscillator having small parameters converges in a singular perturbation limit to the I+F model. In this sense, our study provides a novel unfolding of such models and it identifies a constructible electronic circuit that is closely related to I+F.     

Formation and Dynamics of Waves in a Cortical Model of Cholinergic Modulation

Acetylcholine (ACh) is a regulator of neural excitability and one of the neurochemical substrates of sleep. Amongst the cellular effects induced by cholinergic modulation are a reduction in spike-frequency adaptation (SFA) and a shift in the phase response curve (PRC). We demonstrate in a biophysical model how changes in neural excitability and network structure interact to create three distinct functional regimes: localized asynchronous, traveling asynchronous, and traveling synchronous. Our results qualitatively match those observed experimentally. Cortical activity during slow wave sleep (SWS) differs from that during REM sleep or waking states. During SWS there are traveling patterns of activity in the cortex; in other states stationary patterns occur. Our model is a network composed of Hodgkin-Huxley type neurons with a M-current regulated by ACh. Regulation of ACh level can account for dynamical changes between functional regimes. Reduction of the magnitude of this current recreates the reduction in SFA the shift from a type 2 to a type 1 PRC observed in the presence of ACh. When SFA is minimal (in waking or REM sleep state, high ACh) patterns of activity are localized and easily pinned by network inhomogeneities. When SFA is present (decreasing ACh), traveling waves of activity naturally arise. A further decrease in ACh leads to a high degree of synchrony within traveling waves. We also show that the level of ACh determines how sensitive network activity is to synaptic heterogeneity. These regimes may have a profound functional significance as stationary patterns may play a role in the proper encoding of external input as memory and traveling waves could lead to synaptic regularization, giving unique insights into the role and significance of ACh in determining patterns of cortical activity and functional differences arising from the patterns.     

Modeling of glycolytic wave propagation in an open spatial reactor with inhomogeneous substrate influx

The spatio-temporal dynamics of traveling waves in glycolysis as it occurs in yeast extract have been studied, both theoretically and experimentally. We describe this phenomenon with the distributed Selkov model that accounts for the reactions of phosphofructokinase, which is a key enzyme of the glycolytic reaction cascade. To describe the experimentally observed phase waves in an open spatial reactor we introduce a non-homogeneous flux of substrate in the model. The experimental observation that waves can change their direction of propagation during the experiment is considered in the model. The mechanism for such a change in wave direction is discussed.     

Gates and oscillators: a network model of the brain clock

The suprachiasmatic nuclei (SCN) control circadian oscillations of physiology and behavior. Measurements of electrical activity and of gene expression indicate that these heterogeneous structures are composed of both rhythmic and nonrhythmic cells. A fundamental question with regard to the organization of the circadian system is how the SCN achieve a coherent output while their constituent independent cellular oscillators express a wide range of periods. Previously, the consensus output of individual oscillators had been attributed to coupling among cells. The authors propose a model that incorporates nonrhythmic "gate" cells and rhythmic oscillator cells with a wide range of periods, that neither requires nor excludes a role for interoscillator coupling. The gate provides daily input to oscillator cells and is in turn regulated (directly or indirectly) by the oscillator cells. In the authors' model, individual oscillators with initial random phases are able to self-assemble so as to maintain cohesive rhythmic output. In this view, SCN circuits are important for self-sustained oscillation, and their network properties distinguish these nuclei from other tissues that rhythmically express clock genes. The model explains how individual SCN cells oscillate independently and yet work together to produce a coherent rhythm.     

Gap junctions modulate seizures in a mean-field model of general anesthesia for the cortex

During slow-wave sleep, general anesthesia, and generalized seizures, there is an absence of consciousness. These states are characterized by low-frequency large-amplitude traveling waves in scalp electroencephalogram. Therefore the oscillatory state might be an indication of failure to form coherent neuronal assemblies necessary for consciousness. A generalized seizure event is a pathological brain state that is the clearest manifestation of waves of synchronized neuronal activity. Since gap junctions provide a direct electrical connection between adjoining neurons, thus enhancing synchronous behavior, reducing gap-junction conductance should suppress seizures; however there is no clear experimental evidence for this. Here we report theoretical predictions for a physiologically-based cortical model that describes the general anesthetic phase transition from consciousness to coma, and includes both chemical synaptic and direct electrotonic synapses. The model dynamics exhibits both Hopf (temporal) and Turing (spatial) instabilities; the Hopf instability corresponds to the slow (≲8 Hz) oscillatory states similar to those seen in slow-wave sleep, general anesthesia, and seizures. We argue that a delicately balanced interplay between Hopf and Turing modes provides a canonical mechanism for the default non-cognitive rest state of the brain. We show that the Turing mode, set by gap-junction diffusion, is generally protective against entering oscillatory modes; and that weakening the Turing mode by reducing gap conduction can release an uncontrolled Hopf oscillation and hence an increased propensity for seizure and simultaneously an increased sensitivity to GABAergic anesthesia.     

Waves on Reissner's membrane: a mechanism for the propagation of otoacoustic emissions from the cochlea

Sound is detected and converted into electrical signals within the ear. The cochlea not only acts as a passive detector of sound, however, but can also produce tones itself. These otoacoustic emissions are a striking manifestation of the cochlea's mechanical active process. A controversy remains of how these mechanical signals propagate back to the middle ear, from which they are emitted as sound. Here, we combine theoretical and experimental studies to show that mechanical signals can be transmitted by waves on Reissner's membrane, an elastic structure within the cochlea. We develop a theory for wave propagation on Reissner's membrane and its role in otoacoustic emissions. Employing a scanning laser interferometer, we measure traveling waves on Reissner's membrane in the gerbil, guinea pig, and chinchilla. The results are in accord with the theory and thus support a role for Reissner's membrane in otoacoustic emissions.     

Gβ Regulates Coupling between Actin Oscillators for Cell Polarity and Directional Migration

For directional movement, eukaryotic cells depend on the proper organization of their actin cytoskeleton. This engine of motility is made up of highly dynamic nonequilibrium actin structures such as flashes, oscillations, and traveling waves. In Dictyostelium, oscillatory actin foci interact with signals such as Ras and phosphatidylinositol 3,4,5-trisphosphate (PIP₃) to form protrusions. However, how signaling cues tame actin dynamics to produce a pseudopod and guide cellular motility is a critical open question in eukaryotic chemotaxis. Here, we demonstrate that the strength of coupling between individual actin oscillators controls cell polarization and directional movement. We implement an inducible sequestration system to inactivate the heterotrimeric G protein subunit Gβ and find that this acute perturbation triggers persistent, high-amplitude cortical oscillations of F-actin. Actin oscillators that are normally weakly coupled to one another in wild-type cells become strongly synchronized following acute inactivation of Gβ. This global coupling impairs sensing of internal cues during spontaneous polarization and sensing of external cues during directional motility. A simple mathematical model of coupled actin oscillators reveals the importance of appropriate coupling strength for chemotaxis: moderate coupling can increase sensitivity to noisy inputs. Taken together, our data suggest that Gβ regulates the strength of coupling between actin oscillators for efficient polarity and directional migration. As these observations are only possible following acute inhibition of Gβ and are masked by slow compensation in genetic knockouts, our work also shows that acute loss-of-function approaches can complement and extend the reach of classical genetics in Dictyostelium and likely other systems as well.     

Effects of electrical and structural remodeling on atrial fibrillation maintenance: a simulation study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength--electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely.     

Traveling waves and pulses in a one-dimensional network of excitable integrate-and-fire neurons

 We study the existence and stability of traveling waves and pulses in a one-dimensional network of integrate-and-fire neurons with synaptic coupling. This provides a simple model of excitable neural tissue. We first derive a self-consistency condition for the existence of traveling waves, which generates a dispersion relation between velocity and wavelength. We use this to investigate how wave-propagation depends on various parameters that characterize neuronal interactions such as synaptic and axonal delays, and the passive membrane properties of dendritic cables. We also establish that excitable networks support the propagation of solitary pulses in the long-wavelength limit. We then derive a general condition for the (local) asymptotic stability of traveling waves in terms of the characteristic equation of the linearized firing time map, which takes the form of an integro-difference equation of infinite order. We use this to analyze the stability of solitary pulses in the long-wavelength limit. Solitary wave solutions are shown to come in pairs with the faster (slower) solution stable (unstable) in the case of zero axonal delays; for non-zero delays and fast synapses the stable wave can itself destabilize via a Hopf bifurcation. Received: 27 October 1998     

Oscillator decomposition of infant fNIRS data

The functional near-infrared spectroscopy (fNIRS) can detect hemodynamic responses in the brain and the data consist of bivariate time series of oxygenated hemoglobin (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) on each channel. In this study, we investigate oscillatory changes in infant fNIRS signals by using the oscillator decompisition method (OSC-DECOMP), which is a statistical method for extracting oscillators from time series data based on Gaussian linear state space models. OSC-DECOMP provides a natural decomposition of fNIRS data into oscillation components in a data-driven manner and does not require the arbitrary selection of band-pass filters. We analyzed 18-ch fNIRS data (3 minutes) acquired from 21 sleeping 3-month-old infants. Five to seven oscillators were extracted on most channels, and their frequency distribution had three peaks in the vicinity of 0.01-0.1 Hz, 1.6-2.4 Hz and 3.6-4.4 Hz. The first peak was considered to reflect hemodynamic changes in response to the brain activity, and the phase difference between oxy-Hb and deoxy-Hb for the associated oscillators was at approximately 230 degrees. The second peak was attributed to cardiac pulse waves and mirroring noise. Although these oscillators have close frequencies, OSC-DECOMP can separate them through estimating their different projection patterns on oxy-Hb and deoxy-Hb. The third peak was regarded as the harmonic of the second peak. By comparing the Akaike Information Criterion (AIC) of two state space models, we determined that the time series of oxy-Hb and deoxy-Hb on each channel originate from common oscillatory activity. We also utilized the result of OSC-DECOMP to investigate the frequency-specific functional connectivity. Whereas the brain oscillator exhibited functional connectivity, the pulse waves and mirroring noise oscillators showed spatially homogeneous and independent changes. OSC-DECOMP is a promising tool for data-driven extraction of oscillation components from biological time series data.     

The equation of motion for sperm flagella.

The equation of motion for sperm flagella, in which the elastic bending moment and the active contractile moment are balanced by the moment from the viscous resistance of the surrounding fluid, is solved for a wave solution that superimposes partial solutions. Substitution of the expression for the wave solution into the equation leads to an expression for the active contractile moment. This active moment can be decomposed into two parts. The first part describes an active moment that travels over the flagellum with the mechanical flagellar wave, the second part represents a moment in phase over the entire length of the flagellum, which decreases linearly towards the distal tip. The linear synchronous moment, to which an amount of traveling moment has been added as a perturbation, leads to wave solutions that closely resemble flagellar waves. Properties such as wavelength and wave amplitudes and also the shape of the waves in sea urchin sperm flagella at different frequencies are accurately described by the theory. The change in wave shape in sea urchin sperm flagella at raised viscosity is predicted well by the theory. The different wave properties caused in bull sperm flagella by different boundary conditions at the proximal junction are explained. When only a traveling active moment is present in a flagellum, the wave solutions describe waves of a small wave length in a long flagellum. Some properties of the wave motion of sperm flagella are derived from the theory and verified experimentally.     

Calcium-Actin Waves and Oscillations of Cellular Membranes

We propose a mechanism for the formation of membrane oscillations and traveling waves, which arise due to the coupling between the actin cytoskeleton and the calcium flux through the membrane. In our model, the fluid cell membrane has a mobile but constant population of proteins with a convex spontaneous curvature, which act as nucleators of actin polymerization and adhesion. Such a continuum model couples the forces of cell-substrate adhesion, actin polymerization, membrane curvature, and the flux of calcium through the membrane. Linear stability analysis shows that sufficiently strong coupling among the calcium, membrane, and protein dynamics may induce robust traveling waves on the membrane. This result was checked for a reduced feedback scheme and is compared to the results without the effects of calcium, where permanent phase separation without waves or oscillations is obtained. The model results are compared to the published observations of calcium waves in cell membranes, and a number of testable predictions are proposed.     

Phase computations and phase models for discrete molecular oscillators

Background

Biochemical oscillators perform crucial functions in cells, e.g., they set up circadian clocks. The dynamical behavior of oscillators is best described and analyzed in terms of the scalar quantity, phase. A rigorous and useful definition for phase is based on the so-called isochrons of oscillators. Phase computation techniques for continuous oscillators that are based on isochrons have been used for characterizing the behavior of various types of oscillators under the influence of perturbations such as noise.

Results

In this article, we extend the applicability of these phase computation methods to biochemical oscillators as discrete molecular systems, upon the information obtained from a continuous-state approximation of such oscillators. In particular, we describe techniques for computing the instantaneous phase of discrete, molecular oscillators for stochastic simulation algorithm generated sample paths. We comment on the accuracies and derive certain measures for assessing the feasibilities of the proposed phase computation methods. Phase computation experiments on the sample paths of well-known biological oscillators validate our analyses.

Conclusions

The impact of noise that arises from the discrete and random nature of the mechanisms that make up molecular oscillators can be characterized based on the phase computation techniques proposed in this article. The concept of isochrons is the natural choice upon which the phase notion of oscillators can be founded. The isochron-theoretic phase computation methods that we propose can be applied to discrete molecular oscillators of any dimension, provided that the oscillatory behavior observed in discrete-state does not vanish in a continuous-state approximation. Analysis of the full versatility of phase noise phenomena in molecular oscillators will be possible if a proper phase model theory is developed, without resorting to such approximations.     

The Spatial Origins of Cochlear Amplification Assessed by Stimulus-Frequency Otoacoustic Emissions

Cochlear amplification of basilar membrane traveling waves is thought to occur between a tone’s characteristic frequency (CF) place and within one octave basal of the CF. Evidence for this view comes only from the cochlear base. Stimulus-frequency otoacoustic emissions (SFOAEs) provide a noninvasive alternative to direct measurements of cochlear motion that can be measured across a wide range of CF regions. Coherent reflection theory indicates that SFOAEs arise mostly from the peak region of the traveling wave, but several studies using far-basal suppressor tones claimed that SFOAE components originate many octaves basal of CF. We measured SFOAEs while perfusing guinea pig cochleas from apex to base with salicylate or KCl solutions that reduced outer-hair-cell function and SFOAE amplification. Solution effects on inner hair cells reduced auditory nerve compound action potentials (CAPs) and provided reference times for when solutions reached the SFOAE-frequency CF region. As solution flowed from apex to base, SFOAE reductions generally occurred later than CAP reductions and showed that the effects of cochlear amplification usually peaked ∼1/2 octave basal of the CF region. For tones ≥2 kHz, cochlear amplification typically extended ∼1.5 octaves basal of CF, and the data are consistent with coherent reflection theory. SFOAE amplification did not extend to the basal end of the cochlea, even though reticular lamina motion is amplified in this region, which indicates that reticular lamina motion is not directly coupled to basilar membrane traveling waves. Previous reports of SFOAE-frequency residuals produced by suppressor frequencies far above the SFOAE frequency are most likely due to additional sources created by the suppressor. For some tones <2 kHz, SFOAE amplification extended two octaves apical of CF, which highlights that different vibratory motions produce SFOAEs and CAPs, and that the amplification region depends on the cochlear mode of motion considered. The concept that there is a single “cochlear amplification region” needs to be revised.     

Intrinsic regulation of spatiotemporal organization within the suprachiasmatic nucleus

The mammalian pacemaker in the suprachiasmatic nucleus (SCN) contains a population of neural oscillators capable of sustaining cell-autonomous rhythms in gene expression and electrical firing. A critical question for understanding pacemaker function is how SCN oscillators are organized into a coherent tissue capable of coordinating circadian rhythms in behavior and physiology. Here we undertake a comprehensive analysis of oscillatory function across the SCN of the adult PER2::LUC mouse by developing a novel approach involving multi-position bioluminescence imaging and unbiased computational analyses. We demonstrate that there is phase heterogeneity across all three dimensions of the SCN that is intrinsically regulated and extrinsically modulated by light in a region-specific manner. By investigating the mechanistic bases of SCN phase heterogeneity, we show for the first time that phase differences are not systematically related to regional differences in period, waveform, amplitude, or brightness. Furthermore, phase differences are not related to regional differences in the expression of arginine vasopressin and vasoactive intestinal polypeptide, two key neuropeptides characterizing functionally distinct subdivisions of the SCN. The consistency of SCN spatiotemporal organization across individuals and across planes of section suggests that the precise phasing of oscillators is a robust feature of the pacemaker important for its function.