Toward robust phase-locking in Melibe swim central pattern generator models

Pairing cellular and synaptic dynamics into building blocks of rhythmic neural circuits. A tutorial

In this study we focus on two subnetworks common in the circuitry of swim central pattern generators (CPGs) in the sea slugs, Melibe leonina and Dendronotus iris and show that they are independently capable of stably producing emergent network bursting. This observation raises the question of whether the coordination of redundant bursting mechanisms plays a role in the generation of rhythm and its regulation in the given swim CPGs. To address this question, we investigate two pairwise rhythm-generating networks and examine the properties of their fundamental components: cellular and synaptic, which are crucial for proper network assembly and its stable function. We perform a slow-fast decomposition analysis of cellular dynamics and highlight its significant bifurcations occurring in isolated and coupled neurons. A novel model for slow synapses with high filtering efficiency and temporal delay is also introduced and examined. Our findings demonstrate the existence of two modes of oscillation in bicellular rhythm-generating networks with network hysteresis: i) a half-center oscillator and ii) an excitatory-inhibitory pair. These 2-cell networks offer potential as common building blocks combined in modular organization of larger neural circuits preserving robust network hysteresis.

Error Function Optimization to Compare Neural Activity and Train Blended Rhythmic Networks

We present a novel set of quantitative measures for "likeness" (error function) designed to alleviate the time-consuming and subjective nature of manually comparing biological recordings from electrophysiological experiments with the outcomes of their mathematical models. Our innovative "blended" system approach offers an objective, high-throughput, and computationally efficient method for comparing biological and mathematical models. This approach involves using voltage recordings of biological neurons to drive and train mathematical models, facilitating the derivation of the error function for further parameter optimization. Our calibration process incorporates measurements such as action potential (AP) frequency, voltage moving average, voltage envelopes, and the probability of post-synaptic channels. To assess the effectiveness of our method, we utilized the sea slug Melibe leonina swim central pattern generator (CPG) as our model circuit and conducted electrophysiological experiments with TTX to isolate CPG interneurons. During the comparison of biological recordings and mathematically simulated neurons, we performed a grid search of inhibitory and excitatory synapse conductance. Our findings indicate that a weighted sum of simple functions is essential for comprehensively capturing a neuron's rhythmic activity. Overall, our study suggests that our blended system approach holds promise for enabling objective and high-throughput comparisons between biological and mathematical models, offering significant potential for advancing research in neural circuitry and related fields.

Central pattern generator based on self-sustained oscillator coupled to a chain of oscillatory circuits

We have proposed and studied both numerically and experimentally a multistable system based on a self-sustained Van der Pol oscillator coupled to passive oscillatory circuits. The number of passive oscillators determines the number of multistable oscillatory regimes coexisting in the proposed system. It is shown that our system can be used in robotics applications as a simple model for a central pattern generator (CPG). In this case, the amplitude and phase relations between the active and passive oscillators control a gait, which can be adjusted by changing the system control parameters. Variation of the active oscillator's natural frequency leads to hard switching between the regimes characterized by different phase shifts between the oscillators. In contrast, the external forcing can change the frequency and amplitudes of oscillations, preserving the phase shifts. Therefore, the frequency of the external signal can serve as a control parameter of the model regime and realize a feedback in the proposed CPG depending on the environmental conditions. In particular, it allows one to switch the regime and change the velocity of the robot's gate and tune the gait to the environment. We have also shown that the studied oscillatory regimes in the proposed system are robust and not affected by external noise or fluctuations of the system parameters. Moreover, using the proposed scheme, we simulated the type of bipedal locomotion, including walking and running.

Generalized half-center oscillators with short-term synaptic plasticity

How can we develop simple yet realistic models of the small neural circuits known as central pattern generators (CPGs), which contribute to generate complex multiphase locomotion in living animals? In this paper we introduce a new model (with design criteria) of a generalized half-center oscillator, (pools of) neurons reciprocally coupled by fast/slow inhibitory and excitatory synapses, to produce either alternating bursting or other rhythmic patterns, characterized by different phase lags, depending on the sensory or other external input. We also show how to calibrate its parameters, based on both physiological and functional criteria and on bifurcation analysis. This model accounts for short-term neuromodulation in a biophysically plausible way and is a building block to develop more realistic and functionally accurate CPG models. Examples and counterexamples are used to point out the generality and effectiveness of our design approach.

Design Principles for Central Pattern Generators With Preset Rhythms

This article is concerned with the design of synthetic central pattern generators (CPGs). Biological CPGs are neural circuits that determine a variety of rhythmic activities, including locomotion, in animals. A synthetic CPG is a network of dynamical elements (here called cells) properly coupled by various synapses to emulate rhythms produced by a biological CPG. We focus on CPGs for locomotion of quadrupeds and present our design approach, based on the principles of nonlinear dynamics, bifurcation theory, and parameter optimization. This approach lets us design the synthetic CPG with a set of desired rhythms and switch between them as the parameter representing the control actions from the brain is varied. The developed four-cell CPG can produce four distinct gaits: walk, trot, gallop, and bound, similar to the mouse locomotion. The robustness and adaptability of the network design principles are verified using different cell and synapse models.

Dynamics and bifurcations in multistable 3-cell neural networks

We disclose the generality of the intrinsic mechanisms underlying multistability in reciprocally inhibitory 3-cell circuits composed of simplified, low-dimensional models of oscillatory neurons, as opposed to those of a detailed Hodgkin-Huxley type [Wojcik et al., PLoS One 9, e92918 (2014)]. The computational reduction to return maps for the phase-lags between neurons reveals a rich multiplicity of rhythmic patterns in such circuits. We perform a detailed bifurcation analysis to show how such rhythms can emerge, disappear, and gain or lose stability, as the parameters of the individual cells and the synapses are varied.

Control strategies of 3-cell Central Pattern Generator via global stimuli

The study of the synchronization patterns of small neuron networks that control several biological processes has become an interesting growing discipline. Some of these synchronization patterns of individual neurons are related to some undesirable neurological diseases, and they are believed to play a crucial role in the emergence of pathological rhythmic brain activity in different diseases, like Parkinson’s disease. We show how, with a suitable combination of short and weak global inhibitory and excitatory stimuli over the whole network, we can switch between different stable bursting patterns in small neuron networks (in our case a 3-neuron network). We develop a systematic study showing and explaining the effects of applying the pulses at different moments. Moreover, we compare the technique on a completely symmetric network and on a slightly perturbed one (a much more realistic situation). The present approach of using global stimuli may allow to avoid undesirable synchronization patterns with nonaggressive stimuli.

Two interconnected kernels of reciprocally inhibitory interneurons underlie alternating left-right swim motor pattern generation in the mollusk Melibe leonina

The central pattern generator (CPG) underlying the rhythmic swimming behavior of the nudibranch Melibe leonina (Mollusca, Gastropoda, Heterobranchia) has been described as a simple half-center oscillator consisting of two reciprocally inhibitory pairs of interneurons called swim interneuron 1 (Si1) and swim interneuron 2 (Si2). In this study, we identified two additional pairs of interneurons that are part of the swim CPG: swim interneuron 3 (Si3) and swim interneuron 4 (Si4). The somata of Si3 and Si4 were both located in the pedal ganglion, near that of Si2, and both had axons that projected through the pedal commissure to the contralateral pedal ganglion. These neurons fulfilled the criteria for inclusion as members of the swim CPG: 1) they fired at a fixed phase in relation to Si1 and Si2, 2) brief changes in their activity reset the motor pattern, 3) prolonged changes in their activity altered the periodicity of the motor pattern, 4) they had monosynaptic connections with each other and with Si1 and Si2, and 5) their synaptic actions helped explain the phasing of the motor pattern. The results of this study show that the motor pattern has more complex internal dynamics than a simple left/right alternation of firing; the CPG circuit appears to be composed of two kernels of reciprocally inhibitory neurons, one consisting of Si1, Si2, and the contralateral Si4 and the other consisting of Si3. These two kernels interact with each other to produce a stable rhythmic motor pattern.

Key bifurcations of bursting polyrhythms in 3-cell central pattern generators

We identify and describe the key qualitative rhythmic states in various 3-cell network motifs of a multifunctional central pattern generator (CPG). Such CPGs are neural microcircuits of cells whose synergetic interactions produce multiple states with distinct phase-locked patterns of bursting activity. To study biologically plausible CPG models, we develop a suite of computational tools that reduce the problem of stability and existence of rhythmic patterns in networks to the bifurcation analysis of fixed points and invariant curves of a Poincaré return maps for phase lags between cells. We explore different functional possibilities for motifs involving symmetry breaking and heterogeneity. This is achieved by varying coupling properties of the synapses between the cells and studying the qualitative changes in the structure of the corresponding return maps. Our findings provide a systematic basis for understanding plausible biophysical mechanisms for the regulation of rhythmic patterns generated by various CPGs in the context of motor control such as gait-switching in locomotion. Our analysis does not require knowledge of the equations modeling the system and provides a powerful qualitative approach to studying detailed models of rhythmic behavior. Thus, our approach is applicable to a wide range of biological phenomena beyond motor control.

Introduction to focus issue: rhythms and dynamic transitions in neurological disease: modeling, computation, and experiment

Rhythmic neuronal oscillations across a broad range of frequencies, as well as spatiotemporal phenomena, such as waves and bumps, have been observed in various areas of the brain and proposed as critical to brain function. While there is a long and distinguished history of studying rhythms in nerve cells and neuronal networks in healthy organisms, the association and analysis of rhythms to diseases are more recent developments. Indeed, it is now thought that certain aspects of diseases of the nervous system, such as epilepsy, schizophrenia, Parkinson's, and sleep disorders, are associated with transitions or disruptions of neurological rhythms. This focus issue brings together articles presenting modeling, computational, analytical, and experimental perspectives about rhythms and dynamic transitions between them that are associated to various diseases.