Inferring and quantifying the role of an intrinsic current in a mechanism for a half-center bursting oscillation: A dominant scale and hybrid dynamical systems analysis

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.

A Single Central Pattern Generator for the Control of a Locomotor Rolling Wave in Mollusc Aplysia

Locomotion in mollusc Aplysia is implemented by a pedal rolling wave, a type of axial locomotion. Well-studied examples of axial locomotion (pedal waves in Drosophila larvae and body waves in leech, lamprey, and fish) are generated in a segmented nervous system via activation of multiple coupled central pattern generators (CPGs). Pedal waves in molluscs, however, are generated by a single pedal ganglion, and it is unknown whether there are single or multiple CPGs that generate rhythmic activity and phase shifts between different body parts. During locomotion in intact Aplysia, bursting activity in the parapedal commissural nerve (PPCN) was found to occur during tail contraction. A cluster of 20 to 30 P1 root neurons (P1Ns) on the ventral surface of the pedal ganglion, active during the pedal wave, were identified. Computational cluster analysis revealed that there are 2 phases to the motor program: phase I (centered around 168°) and phase II (centered around 357°). PPCN activity occurs during phase II. The majority of P1Ns are motoneurons. Coactive P1Ns tend to be electrically coupled. Two classes of pedal interneurons (PIs) were characterized. Class 1 (PI1 and PI2) is active during phase I. Their axons make a loop within the pedal ganglion and contribute to locomotor pattern generation. They are electrically coupled to P1Ns that fire during phase I. Class 2 (PI3) is active during phase II and innervates the contralateral pedal ganglion. PI3 may contribute to bilateral coordination. Overall, our findings support the idea that Aplysia pedal waves are generated by a single CPG.

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.

Estimating parameters and predicting membrane voltages with conductance-based neuron models

Recent results demonstrate techniques for fully quantitative, statistical inference of the dynamics of individual neurons under the Hodgkin-Huxley framework of voltage-gated conductances. Using a variational approximation, this approach has been successfully applied to simulated data from model neurons. Here, we use this method to analyze a population of real neurons recorded in a slice preparation of the zebra finch forebrain nucleus HVC. Our results demonstrate that using only 1,500 ms of voltage recorded while injecting a complex current waveform, we can estimate the values of 12 state variables and 72 parameters in a dynamical model, such that the model accurately predicts the responses of the neuron to novel injected currents. A less complex model produced consistently worse predictions, indicating that the additional currents contribute significantly to the dynamics of these neurons. Preliminary results indicate some differences in the channel complement of the models for different classes of HVC neurons, which accords with expectations from the biology. Whereas the model for each cell is incomplete (representing only the somatic compartment, and likely to be missing classes of channels that the real neurons possess), our approach opens the possibility to investigate in modeling the plausibility of additional classes of channels the cell might possess, thus improving the models over time. These results provide an important foundational basis for building biologically realistic network models, such as the one in HVC that contributes to the process of song production and developmental vocal learning in songbirds.

Determining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model

The Hodgkin-Huxley (HH) model is the basis for numerous neural models. There are two negative feedback processes in the HH model that regulate rhythmic spiking. The first is an outward current with an activation variable n that has an opposite influence to the excitatory inward current and therefore provides subtractive negative feedback. The other is the inactivation of an inward current with an inactivation variable h that reduces the amount of positive feedback and therefore provides divisive feedback. Rhythmic spiking can be obtained with either negative feedback process, so we ask what is gained by having two feedback processes. We also ask how the different negative feedback processes contribute to spiking. We show that having two negative feedback processes makes the HH model more robust to changes in applied currents and conductance densities than models that possess only one negative feedback variable. We also show that the contributions made by the subtractive and divisive feedback variables are not static, but depend on time scales and conductance values. In particular, they contribute differently to the dynamics in Type I versus Type II neurons.

Phase response properties of half-center oscillators

We examine the phase response properties of half-center oscillators (HCOs) that are modeled by a pair of Morris-Lecar-type neurons connected by strong fast inhibitory synapses. We find that the two basic mechanisms for half-center oscillations, "release" and "escape", give rise to strikingly different phase response curves (PRCs). Release-type HCOs are most sensitive to perturbations delivered to cells at times when they are about to transition from the active to the suppressed state, and PRCs are dominated by a large negative peak (phase delays) at corresponding phases. On the other hand, escape-type HCOs are most sensitive to perturbations delivered to cells at times when they are about to transition from the suppressed to the active state, and PRCs are dominated by a large positive peak (phase advances) at corresponding phases. By analyzing the phase space structure of Morris-Lecar-type HCO models with fast synaptic dynamics, we identify the dynamical mechanisms underlying the shapes of the PRCs. To demonstrate the significance of the different shapes of the PRCs for the release-type and escape-type HCOs, we link the shapes of the PRCs to the different frequency modulation properties of release-type and escape-type HCOs, and we show that the different shapes of the PRCs for the release-type and escape-type HCOs can lead to fundamentally different phase-locking dynamics.

Hybrid models and biological model reduction with PyDSTool

The PyDSTool software environment is designed to develop, simulate, and analyze dynamical systems models, particularly for biological applications. Unlike the engineering application focus and graphical specification environments of most general purpose simulation tools, PyDSTool provides a programmatic environment well suited to exploratory data- and hypothesis-driven biological modeling problems. In this work, we show how the environment facilitates the application of hybrid dynamical modeling to the reverse engineering of complex biophysical dynamics; in this case, of an excitable membrane. The example demonstrates how the software provides novel tools that support the inference and validation of mechanistic hypotheses and the inclusion of data constraints in both quantitative and qualitative ways. The biophysical application is broadly relevant to models in the biosciences. The open source and platform-independent PyDSTool package is freely available under the BSD license from http://sourceforge.net/projects/pydstool/. The hosting service provides links to documentation and online forums for user support.