Dendritic and Synaptic Effects in Systems of Coupled Cortical Oscillators

We explore the influence of synaptic location and form on the behavior of networks of coupled cortical oscillators. First, we develop a model of two coupled somatic oscillators that includes passive dendritic cables. Using a phase model approach, we show that the synchronous solution can change from a stable solution to an unstable one as the cable lengthens and the synaptic position moves further from the soma. We confirm this prediction using a system of coupled compartmental models. We also demonstrate that when the synchronous solution becomes unstable, a bifurcation occurs and a pair of asynchronous stable solutions appear, causing a phase lag between the cells in the system. Then using a variety of coupling functions and different synaptic positions, we show that distal connections and broad synaptic time courses encourage phase lags that can be reduced, eliminated, or enhanced by the presence of active currents in the dendrite. This mechanism may appear in neural systems where proximal connections could be used to encourage synchrony, and distal connections and broad synaptic time courses could be used to produce phase lags that can be modulated by active currents.     

Rapid synchronization through fast threshold modulation

Synchronization properties of locally coupled neural oscillators were investigated analytically and by computer simulation. When coupled in a manner that mimics excitatory chemical synapses, oscillators having more than one time scale (relaxation oscillators) are shown to approach synchrony using mechanisms very different from that of oscillators with a more sinusoidal waveform. The relaxation oscillators make critical use of fast modulations of their thresholds, leading to a rate of synchronization relatively independent of coupling strength within some basin of attraction; this rate is faster for oscillators that have conductance-based features than for neural caricatures such as the FitzHugh-Nagumo equations that lack such features. Computer simulations of one-dimensional arrays show that oscillators in the relaxation regime synchronize much more rapidly than oscillators with the same equations whose parameters have been modulated to yield a more sinusoidal waveform. We present a heuristic explanation of this effect based on properties of the coupling mechanisms that can affect the way the synchronization scales with array length. These results suggest that the emergent synchronization behavior of oscillating neural networks can be dramatically influenced by the intrinsic properties of the network components. Possible implications for perceptual feature binding and attention are discussed.Supported in part by NASA (NGT-50497)Supported in part by NSF (DMS-8901913), and NIMH-47150 Present address and address for correspondence: Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, E25-618, Cambridge, MA 02139, USA     

Symmetry-breaking bifurcation: A possible mechanism for 2:1 frequency-locking in animal locomotion

The generation and control of animal locomotion is believed to involve central pattern generators — networks of neurons which are capable of producing oscillatory behavior. In the present work, the quadrupedal locomotor central pattern generator is modelled as four distinct but symmetrically coupled nonlinear oscillators. We show that the typical patterns for two such networks of oscillators include 2:1 frequency-locked oscillations. These patterns, which arise through symmetry-breaking Hopf bifurcation, correspond in part to observed patterns of 2:1 frequency-locking of limb movements during electrically elicited locomotion of decerebrate and spinal quadrupeds. We briefly describe how our theoretical predictions could be tested experimentally.     

Hexapodal gaits and coupled nonlinear oscillator models

The general, model-independent features of different networks of six symmetrically coupled nonlinear oscillators are investigated. These networks are considered as possible models for locomotor central pattern generators (CPGs) in insects. Numerical experiments with a specific oscillator network model are briefly described. It is shown that some generic phase-locked oscillation-patterns for various systems of six symmetrically coupled nonlinear oscillators correspond to the common forward-walking gaits adopted by insects. It is also demonstrated that transitions observed in insect gaits can be modelled as standard symmetry-breaking bifurcations occurring in such systems. The present analysis, which leads to a natural classification of hexapodal gaits by symmetry and to natural sequences of gait bifurcations, relates observed gaits to the overall organizational structure of the underlying CPG. The implications of the present results for the development of simplified control systems for hexapodal walking robots are discussed.     

Pulse coupled oscillators and the phase resetting curve

Limit cycle oscillators that are coupled in a pulsatile manner are referred to as pulse coupled oscillators. In these oscillators, the interactions take the form of brief pulses such that the effect of one input dies out before the next is received. A phase resetting curve (PRC) keeps track of how much an input advances or delays the next spike in an oscillatory neuron depending upon where in the cycle the input is applied. PRCs can be used to predict phase locking in networks of pulse coupled oscillators. In some studies of pulse coupled oscillators, a specific form is assumed for the interactions between oscillators, but a more general approach is to formulate the problem assuming a PRC that is generated using a perturbation that approximates the input received in the real biological network. In general, this approach requires that circuit architecture and a specific firing pattern be assumed. This allows the construction of discrete maps from one event to the next. The fixed points of these maps correspond to periodic firing modes and are easier to locate and analyze for stability compared to locating and analyzing periodic modes in the original network directly. Alternatively, maps based on the PRC have been constructed that do not presuppose a firing order. Specific circuits that have been analyzed under the assumption of pulsatile coupling include one to one lockings in a periodically forced oscillator or an oscillator forced at a fixed delay after a threshold event, two bidirectionally coupled oscillators with and without delays, a unidirectional N-ring of oscillators, and N all-to-all networks.     

Multi-Stability and Pattern-Selection in Oscillatory Networks with Fast Inhibition and Electrical Synapses

A model or hybrid network consisting of oscillatory cells interconnected by inhibitory and electrical synapses may express different stable activity patterns without any change of network topology or parameters, and switching between the patterns can be induced by specific transient signals. However, little is known of properties of such signals. In the present study, we employ numerical simulations of neural networks of different size composed of relaxation oscillators, to investigate switching between in-phase (IP) and anti-phase (AP) activity patterns. We show that the time windows of susceptibility to switching between the patterns are similar in 2-, 4- and 6-cell fully-connected networks. Moreover, in a network (N = 4, 6) expressing a given AP pattern, a stimulus with a given profile consisting of depolarizing and hyperpolarizing signals sent to different subpopulations of cells can evoke switching to another AP pattern. Interestingly, the resulting pattern encodes the profile of the switching stimulus. These results can be extended to different network architectures. Indeed, relaxation oscillators are not only models of cellular pacemakers, bursting or spiking, but are also analogous to firing-rate models of neural activity. We show that rules of switching similar to those found for relaxation oscillators apply to oscillating circuits of excitatory cells interconnected by electrical synapses and cross-inhibition. Our results suggest that incoming information, arriving in a proper time window, may be stored in an oscillatory network in the form of a specific spatio-temporal activity pattern which is expressed until new pertinent information arrives.     

A model for intersegmental coordination in the leech nerve cord

The neuronal circuits that generate swimming movements in the leech were simulated by a chain of coupled harmonic oscillators. Our model incorporates a gradient of rostrocaudally decreasing cycle periods along the oscillator chain, a finite conduction delay for coupling signals, and multiple coupling channels connecting each pair of oscillators. The interactions mediated by these channels are characterized by sinusoidal phase response curves. Investigations of this model were carried out with the aid of a digital computer and the results of a variety of manipulations were compared with data from analogous physiological experiments. The simulations reproduced many aspects of intersegmental coordination in the leech, including the findings that: 1) phase lags between adjacent ganglia are larger near the caudal than the rostral end of the leech nerve cord; 2) intersegmental phase lags increase as the number of ganglia in nervecord preparations is reduced; 3) severing one of the paired lateral connective nerves can reverse the phase lag across the lesion and 4) blocking synaptic transmission in midganglia of the ventral nerve cord reduces phase lags across the block.     

Neural mechanisms potentially contributing to the intersegmental phase lag in lamprey

Factors contributing to the production of a phase lag along chains of oscillatory networks consisting of Hodgkin-Huxley type neurons are analyzed by means of simulations. Simplified network configurations are explored consisting of the basic building blocks of the spinal central pattern generator (CPG) generating swimming in the lamprey. It consists of reciprocally coupled crossed inhibitory C interneurons and ipsilateral excitatory E interneurons that activate C neurons and other E neurons. Oscillatory activity in the model network can, in the simplest case, be produced by a pair of reciprocally coupled C interneurons oscillating through an escape mechanism. Different levels of tonic excitation drive the network over a wide burst frequency range. In this type of network, powerful frequency-regulating factors are the effective inhibition produced by the active side, in combination with the tendency of the inactive side to escape from the inhibition. These two mechanisms can be affected by several factors, e.g. spike frequency adaptation (calcium-dependent K(+) channels), N-methyl-D-aspartate membrane properties as well as presence of low-voltage activated calcium channels. A rostrocaudal phase lag can be produced either by extending the contralateral inhibitory projections or the ipsilateral excitatory projections relatively more in the caudal than the rostral direction, since both an increased inhibition and a phasic excitation slow down the receiving network. The phase lag becomes decreased if the length of the intersegmental projections is increased or if the projections are extended symmetrically in both the rostral and the caudal directions. The simulations indicate that the conditions in the ends of an oscillator chain may significantly affect sign, magnitude and constancy of the phase lag. Also, with short and relatively weak intersegmental connections, the network remains robust against perturbations as well as intrinsic frequency differences along the chain. The phase lag (percentage of cycle duration) increases, however, with burst frequency also when the coupling strength is comparatively weak. The results are discussed and compared with previous "phase pulling" models as well as relaxation oscillators.     

Mathematical models for the swimming pattern of a lamprey

A significant characteristic in a swimming pattern of a lamprey is the generation of a constant phase lag along its body in spite of the wide range of undulation frequencies. In this paper, we discuss a mathematical treatment for coupled oscillators with time-delayed interaction and propose a model for the central pattern generator (CPG) of a lamprey to account for the generation of a constant phase relation, with consideration of the signal conduction time. From this model, it is suggested that the desired phase relation can be produced by long ascending connections from the tail to the neck region of the CPG.     

Evaluating functional roles of phase resetting in generation of adaptive human bipedal walking with a physiologically based model of the spinal pattern generator

The central pattern generators (CPGs) in the spinal cord strongly contribute to locomotor behavior. To achieve adaptive locomotion, locomotor rhythm generated by the CPGs is suggested to be functionally modulated by phase resetting based on sensory afferent or perturbations. Although phase resetting has been investigated during fictive locomotion in cats, its functional roles in actual locomotion have not been clarified. Recently, simulation studies have been conducted to examine the roles of phase resetting during human bipedal walking, assuming that locomotion is generated based on prescribed kinematics and feedback control. However, such kinematically based modeling cannot be used to fully elucidate the mechanisms of adaptation. In this article we proposed a more physiologically based mathematical model of the neural system for locomotion and investigated the functional roles of phase resetting. We constructed a locomotor CPG model based on a two-layered hierarchical network model of the rhythm generator (RG) and pattern formation (PF) networks. The RG model produces rhythm information using phase oscillators and regulates it by phase resetting based on foot-contact information. The PF model creates feedforward command signals based on rhythm information, which consists of the combination of five rectangular pulses based on previous analyses of muscle synergy. Simulation results showed that our model establishes adaptive walking against perturbing forces and variations in the environment, with phase resetting playing important roles in increasing the robustness of responses, suggesting that this mechanism of regulation may contribute to the generation of adaptive human bipedal locomotion.     

Spatially Organized Dynamical States in Chemical Oscillator Networks: Synchronization,Dynamical Differentiation,and Chimera Patterns

Dynamical processes in many engineered and living systems take place on complex networks of discrete dynamical units. We present laboratory experiments with a networked chemical system of nickel electrodissolution in which synchronization patterns are recorded in systems with smooth periodic, relaxation periodic, and chaotic oscillators organized in networks composed of up to twenty dynamical units and 140 connections. The reaction system formed domains of synchronization patterns that are strongly affected by the architecture of the network. Spatially organized partial synchronization could be observed either due to densely connected network nodes or through the ‘chimera’ symmetry breaking mechanism. Relaxation periodic and chaotic oscillators formed structures by dynamical differentiation. We have identified effects of network structure on pattern selection (through permutation symmetry and coupling directness) and on formation of hierarchical and ‘fuzzy’ clusters. With chaotic oscillators we provide experimental evidence that critical coupling strengths at which transition to identical synchronization occurs can be interpreted by experiments with a pair of oscillators and analysis of the eigenvalues of the Laplacian connectivity matrix. The experiments thus provide an insight into the extent of the impact of the architecture of a network on self-organized synchronization patterns.     

Traveling electrical waves in cortex: insights from phase dynamics and speculation on a computational role

The theory of coupled phase oscillators provides a framework to understand the emergent properties of networks of neuronal oscillators. When the architecture of the network is dominated by short-range connections, the pattern of electrical output is predicted to correspond to traveling plane and rotating waves, in addition to synchronized output. We argue that this theory provides the foundation for understanding the traveling electrical waves that are observed across olfactory, visual, and visuomotor areas of cortex in a variety of species. The waves are typically present during periods outside of stimulation, while synchronous activity typically dominates in the presence of a strong stimulus. We suggest that the continuum of phase shifts during epochs with traveling waves provides a means to scan the incoming sensory stream for novel features. Experiments to test our theoretical approach are presented.     

Sensory and central mechanisms control intersegmental coordination

Recent experiments on the sensory and central mechanisms that coordinate animal locomotory movements have advanced our understanding of the relative importance of these two components and overturned some previously held notions. In different experimental preparations, sensory inputs and central pattern generators have now been shown to play different roles in setting intersegmental phase lags.     

Analysis of clustered firing patterns in synaptically coupled networks of oscillators

Oscillators in networks may display a variety of activity patterns. This paper presents a geometric singular perturbation analysis of clustering, or alternate firing of synchronized subgroups, among synaptically coupled oscillators. We consider oscillators in two types of networks: mutually coupled, with all-to-all inhibitory connections, and globally inhibitory, with one excitatory and one inhibitory population of oscillators, each of arbitrary size. Our analysis yields existence and stability conditions for clustered states, along with formulas for the periods of such firing patterns. By using two different approaches, we derive complementary conditions, the first set stated in terms of time lengths determined by intrinsic and synaptic properties of the oscillators and their coupling and the second set stated in terms of model parameters and phase space structures directly linked to parameters. These results suggest how biological components may interact to produce the spindle sleep rhythm in thalamocortical networks. Received: 9 September 1999 / Revised version: 7 July 2000 / Published online: 24 November 2000     

Synaptic mechanisms in invertebrate pattern generation

Although individual neurons can be intrinsically oscillatory and can be network pacemakers, motor patterns are often generated in a more distributed manner. Synaptic connections with other neurons are important because they either modify the rhythm of the pacemaker cell or are essential for pattern generation in the first place. Computational studies of half-center oscillators have made much progress in describing how neurons make transitions between active and inactive phases in these simple networks. In addition to characterizing phase transitions, recent studies have described the synaptic mechanisms that are important for the initiation and maintenance of activity in half-center oscillators.     

Analysis of the gait generation principle by a simulated quadruped model with a CPG incorporating vestibular modulation

This study aims to understand the principles of gait generation in a quadrupedal model. It is difficult to determine the essence of gait generation simply by observation of the movement of complicated animals composed of brains, nerves, muscles, etc. Therefore, we build a planar quadruped model with simplified nervous system and mechanisms, in order to observe its gaits under simulation. The model is equipped with a mathematical central pattern generator (CPG), consisting of four coupled neural oscillators, basically producing a trot pattern. The model also contains sensory feedback to the CPG, measuring the body tilt (vestibular modulation). This spontaneously gives rise to an unprogrammed lateral walk at low speeds, a transverse gallop while running, in addition to trotting at a medium speed. This is because the body oscillation exhibits a double peak per leg frequency at low speeds, no peak (little oscillation) at medium speeds, and a single peak while running. The body oscillation autonomously adjusts the phase differences between the neural oscillators via the feedback. We assume that the oscillations of the four legs produced by the CPG and the body oscillation varying according to the current speed are synchronized along with the varied phase differences to keep balance during locomotion through postural adaptation via the vestibular modulation, resulting in each gait. We succeeded in determining a single simple principle that accounts for gait transition from walking to trotting to galloping, even without brain control, complicated leg mechanisms, or a flexible trunk.     

Phase Resetting Curves Allow for Simple and Accurate Prediction of Robust N:1 Phase Locking for Strongly Coupled Neural Oscillators

Existence and stability criteria for harmonic locking modes were derived for two reciprocally pulse coupled oscillators based on their first and second order phase resetting curves. Our theoretical methods are general in the sense that no assumptions about the strength of coupling, type of synaptic coupling, and model are made. These methods were then tested using two reciprocally inhibitory Wang and Buzsáki model neurons. The existence of bands of 2:1, 3:1, 4:1, and 5:1 phase locking in the relative frequency parameter space was predicted correctly, as was the phase of the slow neuron's spike within the cycle of the fast neuron in which it occurred. For weak coupling the bands are very narrow, but strong coupling broadens the bands. The predictions of the pulse coupled method agreed with weak coupling methods in the weak coupling regime, but extended predictability into the strong coupling regime. We show that our prediction method generalizes to pairs of neural oscillators coupled through excitatory synapses, and to networks of multiple oscillatory neurons. The main limitation of the method is the central assumption that the effect of each input dies out before the next input is received.     

Mutually synchronized relaxation oscillators as prototypes of oscillating systems in biology

Summary This paper discusses the analogy between phenomena in populations of coupled biological oscillators and the behaviour of systems of synchronized mathematical oscillators. Frequency entrainment in a set of coupled relaxation oscillators is investigated with perturbation methods. This analysis leads to quantitative results for entrainment and explains phenomena such as travelling waves in systems of spatially distributed oscillators.     

Activity-dependent modulation of adaptation produces a constant burst proportion in a model of the lamprey spinal locomotor generator

The neuronal network underlying lamprey swimming has stimulated extensive modelling on different levels of abstraction. The lamprey swims with a burst frequency ranging from 0.3 to 8–10 Hz with a rostro-caudal lag between bursts in each segment along the spinal cord. The swimming motor pattern is characterized by a burst proportion that is independent of burst frequency and lasts around 30%–40% of the cycle duration. This also applies in preparations in which the reciprocal inhibition in the spinal cord between the left and right side is blocked. A network of coupled excitatory neurons producing hemisegmental oscillations may form the basis of the lamprey central pattern generator (CPG). Here we explored how such networks, in principle, could produce a large frequency range with a constant burst proportion. The computer simulations of the lamprey CPG use simplified, graded output units that could represent populations of neurons and that exhibit adaptation. We investigated the effect of an active modulation of the degree of adaptation of the CPG units to accomplish a constant burst proportion over the whole frequency range when, in addition, each hemisegment is assumed to be self-oscillatory. The degree of adaptation is increased with the degree of stimulation of the network. This will make the bursts terminate earlier at higher burst rates, allowing for a constant burst proportion. Without modulated adaptation the network operates in a limited range of swimming frequencies due to a progressive increase of burst duration with increasing background stimulation. By introducing a modulation of the adaptation, a broad burst frequency range can be produced. The reciprocal inhibition is thus not the primary burst terminating factor, as in many CPG models, and it is mainly responsible for producing alternation between the left and right sides. The results are compared with the Morris-Lecar oscillator model with parameters set to produce a type A and type B oscillator, in which the burst durations stay constant or increase, respectively, when the background stimulation is increased. Here as well, burst duration can be controlled by modulation of the slow variable in a similar way as above. When oscillatory hemisegmental networks are coupled together in a chain a phase lag is produced. The production of a phase lag in chains of such oscillators is compared with chains of Morris-Lecar relaxation oscillators. Models relating to the intact versus isolated spinal cord preparation are discussed, as well as the role of descending inhibition. Received: 1 April 1997 / Accepted in revised form: 20 March 1998     

The role of phase shifts of sensory inputs in walking revealed by means of phase reduction

Detailed neural network models of animal locomotion are important means to understand the underlying mechanisms that control the coordinated movement of individual limbs. Daun-Gruhn and Tóth, Journal of Computational Neuroscience 31(2), 43–60 (2011) constructed an inter-segmental network model of stick insect locomotion consisting of three interconnected central pattern generators (CPGs) that are associated with the protraction-retraction movements of the front, middle and hind leg. This model could reproduce the basic locomotion coordination patterns, such as tri- and tetrapod, and the transitions between them. However, the analysis of such a system is a formidable task because of its large number of variables and parameters. In this study, we employed phase reduction and averaging theory to this large network model in order to reduce it to a system of coupled phase oscillators. This enabled us to analyze the complex behavior of the system in a reduced parameter space. In this paper, we show that the reduced model reproduces the results of the original model. By analyzing the interaction of just two coupled phase oscillators, we found that the neighboring CPGs could operate within distinct regimes, depending on the phase shift between the sensory inputs from the extremities and the phases of the individual CPGs. We demonstrate that this dependence is essential to produce different coordination patterns and the transition between them. Additionally, applying averaging theory to the system of all three phase oscillators, we calculate the stable fixed points - they correspond to stable tripod or tetrapod coordination patterns and identify two ways of transition between them.