Generalized half-center oscillators with short-term synaptic plasticity

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.

Dynamics of antiphase bursting modulated by the inhibitory synaptic and hyperpolarization-activated cation currents

Antiphase bursting related to the rhythmic motor behavior exhibits complex dynamics modulated by the inhibitory synaptic current (Isyn), especially in the presence of the hyperpolarization-activated cation current (Ih). In the present paper, the dynamics of antiphase bursting modulated by the Ih and Isyn is studied in three aspects with a theoretical model. Firstly, the Isyn and the slow Ih with strong strength are the identified to be the necessary conditions for the antiphase bursting. The dependence of the antiphase bursting on the two currents is different for low (escape mode) and high (release mode) threshold voltages (Vth) of the inhibitory synapse. Secondly, more detailed co-regulations of the two currents to induce opposite changes of the bursting period are obtained. For the escape mode, increase of the Ih induces elevated membrane potential of the silence inhibited by a strong Isyn and shortened silence duration to go beyond Vth, resulting in reduced bursting period. For the release mode, increase of the Ih induces elevated tough value of the former part of the burst modulated by a nearly zero Isyn and lengthen burst duration to fall below Vth, resulting in prolonged bursting period. Finally, the fast-slow dynamics of the antiphase bursting are acquired. Using one-and two-parameter bifurcations of the fast subsystem of a single neuron, the burst of the antiphase bursting is related to the stable limit cycle, and the silence modulated by a strong Isyn to the stable equilibrium to a certain extent. The Ih mainly modulates the dynamics within the burst and quiescent state. Furthermore, with the fast subsystem of the coupled neurons, the silence is associated with the unstable equilibrium point. The results present theoretical explanations to the changes in the bursting period and fast-slow dynamics of the antiphase bursting modulated by the Isyn and Ih, which is helpful for understanding the antiphase bursting and modulating rhythmic motor patterns.

Optimization strategies to obtain smooth gait transitions through biologically plausible central pattern generators

Central pattern generators are small networks that contribute to generating animal locomotion. The models used to study gait generation and gait transition mechanisms often require biologically accurate neuron and synapse models, with high dimensionality and complex dynamics. Tuning the parameters of these models to elicit network dynamics compatible with gait features is not a trivial task, due to the impossibility of inferring a priori the effects of each parameter on the nonlinear system's emergent dynamics. In this paper we explore the use of global optimization strategies for parameter optimization in multigait central pattern generator (CPG) models with complex cell dynamics and minimal topology. We first consider an existing quadruped CPG model as a test bed for the objective function formulation, then proceed to optimize the parameters of a newly proposed multigait, interlimb hexapod CPG model. We successfully obtain hexapod gaits and prompt gait transitions by varying only control currents, while all CPG parameters, once optimized, are kept fixed. This mechanism of gait transitions is compatible with short-term synaptic plasticity.

Towards more biologically plausible central-pattern-generator models

Central pattern generators (CPGs) are relatively small neural networks that play a fundamental role in the control of animal locomotion. In this paper we define a method for the systematic design of CPG models able to exhibit biologically plausible gait transitions by implementing short-term synaptic plasticity mechanisms. As a case study, we focus on a simple CPG for quadruped locomotion. By applying the proposed method, three of four standard quadruped gaits were correctly reproduced by the obtained CPG model, not only in terms of the alternating sequence of the limbs but also in terms of frequency, duty cycle, and phase lags.

Computational Modeling of Spinal Locomotor Circuitry in the Age of Molecular Genetics

Neuronal circuits in the spinal cord are essential for the control of locomotion. They integrate supraspinal commands and afferent feedback signals to produce coordinated rhythmic muscle activations necessary for stable locomotion. For several decades, computational modeling has complemented experimental studies by providing a mechanistic rationale for experimental observations and by deriving experimentally testable predictions. This symbiotic relationship between experimental and computational approaches has resulted in numerous fundamental insights. With recent advances in molecular and genetic methods, it has become possible to manipulate specific constituent elements of the spinal circuitry and relate them to locomotor behavior. This has led to computational modeling studies investigating mechanisms at the level of genetically defined neuronal populations and their interactions. We review literature on the spinal locomotor circuitry from a computational perspective. By reviewing examples leading up to and in the age of molecular genetics, we demonstrate the importance of computational modeling and its interactions with experiments. Moving forward, neuromechanical models with neuronal circuitry modeled at the level of genetically defined neuronal populations will be required to further unravel the mechanisms by which neuronal interactions lead to locomotor behavior.