Two-dimensional variation of bursting properties in a silicon-neuron half-center oscillator

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.

Neuromorphic Instantiation of Spiking Half-Centered Oscillator Models for Central Pattern Generation

In both invertebrate and vertebrate animals, small networks called central pattern generators (CPGs) form the building blocks of the neuronal circuits involved in locomotion. Most CPGs contain a simple half-center oscillator (HCO) motif which consists of two neurons, or populations of neurons, connected by reciprocal inhibition. CPGs and HCOs are well characterized neuronal networks and have been extensively modeled at different levels of abstraction. In the past two decades, hardware implementation of spiking CPG and HCO models in neuromorphic hardware has opened up new applications in mobile robotics, computational neuroscience, and neuroprosthetics. Despite their relative simplicity, the parameter space of GPG and HCO models can become exhaustive when considering various neuron models and network topologies. Motivated by computational work in neuroscience that used a brute-force approach to generate a large database of millions of simulations of the heartbeat HCO of the leech, we have started to build a database of spiking chains of multiple HCOs for different neuron model types and network topologies. Here we present preliminary results using the Izhikevich and Morris-Lecar neuron models for single and pairs of HCOs with different inter-HCO coupling schemes.

Qualitative-Modeling-Based Silicon Neurons and Their Networks

The ionic conductance models of neuronal cells can finely reproduce a wide variety of complex neuronal activities. However, the complexity of these models has prompted the development of qualitative neuron models. They are described by differential equations with a reduced number of variables and their low-dimensional polynomials, which retain the core mathematical structures. Such simple models form the foundation of a bottom-up approach in computational and theoretical neuroscience. We proposed a qualitative-modeling-based approach for designing silicon neuron circuits, in which the mathematical structures in the polynomial-based qualitative models are reproduced by differential equations with silicon-native expressions. This approach can realize low-power-consuming circuits that can be configured to realize various classes of neuronal cells. In this article, our qualitative-modeling-based silicon neuron circuits for analog and digital implementations are quickly reviewed. One of our CMOS analog silicon neuron circuits can realize a variety of neuronal activities with a power consumption less than 72 nW. The square-wave bursting mode of this circuit is explained. Another circuit can realize Class I and II neuronal activities with about 3 nW. Our digital silicon neuron circuit can also realize these classes. An auto-associative memory realized on an all-to-all connected network of these silicon neurons is also reviewed, in which the neuron class plays important roles in its performance.

Robust synchronization of coupled neural oscillators using the derivative-free nonlinear Kalman Filter

A synchronizing control scheme for coupled neural oscillators of the FitzHugh-Nagumo type is proposed. Using differential flatness theory the dynamical model of two coupled neural oscillators is transformed into an equivalent model in the linear canonical (Brunovsky) form. A similar linearized description is succeeded using differential geometry methods and the computation of Lie derivatives. For such a model it becomes possible to design a state feedback controller that assures the synchronization of the membrane's voltage variations for the two neurons. To compensate for disturbances that affect the neurons' model as well as for parametric uncertainties and variations a disturbance observer is designed based on Kalman Filtering. This consists of implementation of the standard Kalman Filter recursion on the linearized equivalent model of the coupled neurons and computation of state and disturbance estimates using the diffeomorphism (relations about state variables transformation) provided by differential flatness theory. After estimating the disturbance terms in the neurons' model their compensation becomes possible. The performance of the synchronization control loop is tested through simulation experiments.

A modeling approach on why simple central pattern generators are built of irregular neurons

The crustacean pyloric Central Pattern Generator (CPG) is a nervous circuit that endogenously provides periodic motor patterns. Even after about 40 years of intensive studies, the rhythm genesis is still not rigorously understood in this CPG, mainly because it is made of neurons with irregular intrinsic activity. Using mathematical models we addressed the question of using a network of irregularly behaving elements to generate periodic oscillations, and we show some advantages of using non-periodic neurons with intrinsic behavior in the transition from bursting to tonic spiking (as found in biological pyloric CPGs) as building components. We studied two- and three-neuron model CPGs built either with Hindmarsh-Rose or with conductance-based Hodgkin-Huxley-like model neurons. By changing a model’s parameter we could span the neuron’s intrinsic dynamical behavior from slow periodic bursting to fast tonic spiking, passing through a transition where irregular bursting was observed. Two-neuron CPG, half center oscillator (HCO), was obtained for each intrinsic behavior of the neurons by coupling them with mutual symmetric synaptic inhibition. Most of these HCOs presented regular antiphasic bursting activity and the changes of the bursting frequencies was studied as a function of the inhibitory synaptic strength. Among all HCOs, those made of intrinsic irregular neurons presented a wider burst frequency range while keeping a reliable regular oscillatory (bursting) behavior. HCOs of periodic neurons tended to be either hard to change their behavior with synaptic strength variations (slow periodic burster neurons) or unable to perform a physiologically meaningful rhythm (fast tonic spiking neurons). Moreover, 3-neuron CPGs with connectivity and output similar to those of the pyloric CPG presented the same results.

Improving noise resistance of intrinsic rhythms in a square-wave burster model

The square-wave burster (Wang and Rinzel, 2003) is a class of autonomous bursting cells that share a bifurcation structure. It is known that this class of cells is involved in the generation of various life-supporting rhythms. In our research to realize an electronic circuit that mimics the rhythm generating mechanism in the square-wave burster, our circuit experimentally exhibited severe fluctuations in its rhythmic activity. We have found a noise-sensitive region in the phase portrait of the ideal model and have proposed modifications of the model that can reduce this fluctuation. A possible modification to ionic-conductance neuron models (Kohno and Aihara, 2011) was inspired by them. This modification, however, cannot be applied to a group of square-wave bursters, including the Butera-Rinzel-Smith model (Butera et al., 1999; Del Negro et al., 2001), which is a model of the pre-Bötzinger complex bursting neuron that plays a crucial role in the generation of respiration rhythms, because this modification premises that the slow dynamics originates from an activation gate variable of a hyperpolarizing ionic current. However, in some square-wave bursters, they are controlled by an inactivation gate variable of a depolarizing ionic current. In this study, we proposed a similar modification with a completely different mechanism that can be applied to this group of square-wave bursters. In the presence of noises, the modified Butera-Rinzel-Smith model can generate rhythmic activity that is more stable and similar to biological observations than the original model. The mechanisms underlying this modification are explained with noisy bifurcation diagrams.

Emergent central pattern generator behavior in gap-junction-coupled Hodgkin-Huxley style neuron model

Most models of central pattern generators (CPGs) involve two distinct nuclei mutually inhibiting one another via synapses. Here, we present a single-nucleus model of biologically realistic Hodgkin-Huxley neurons with random gap junction coupling. Despite no explicit division of neurons into two groups, we observe a spontaneous division of neurons into two distinct firing groups. In addition, we also demonstrate this phenomenon in a simplified version of the model, highlighting the importance of afterhyperpolarization currents (I(AHP)) to CPGs utilizing gap junction coupling. The properties of these CPGs also appear sensitive to gap junction conductance, probability of gap junction coupling between cells, topology of gap junction coupling, and, to a lesser extent, input current into our simulated nucleus.

Biophysical Neural Spiking, Bursting, and Excitability Dynamics in Reconfigurable Analog VLSI

We study a range of neural dynamics under variations in biophysical parameters underlying extended Morris-Lecar and Hodgkin-Huxley models in three gating variables. The extended models are implemented in NeuroDyn, a four neuron, twelve synapse continuous-time analog VLSI programmable neural emulation platform with generalized channel kinetics and biophysical membrane dynamics. The dynamics exhibit a wide range of time scales extending beyond 100 ms neglected in typical silicon models of tonic spiking neurons. Circuit simulations and measurements show transition from tonic spiking to tonic bursting dynamics through variation of a single conductance parameter governing calcium recovery. We similarly demonstrate transition from graded to all-or-none neural excitability in the onset of spiking dynamics through the variation of channel kinetic parameters governing the speed of potassium activation. Other combinations of variations in conductance and channel kinetic parameters give rise to phasic spiking and spike frequency adaptation dynamics. The NeuroDyn chip consumes 1.29 mW and occupies 3 mm × 3 mm in 0.5 μm CMOS, supporting emerging developments in neuromorphic silicon-neuron interfaces.

A library of analog operators based on the hodgkin-huxley formalism for the design of tunable, real-time, silicon neurons

In this paper, we present a library of analog operators used for the analog real-time computation of the Hodgkin-Huxley formalism. These operators make it possible to design a silicon (Si) neuron that is dynamically tunable, and that reproduces different kinds of neurons. We used an original method in neuromorphic engineering to characterize this Si neuron. In electrophysiology, this method is well known as the "voltage-clamp" technique. We also compare the features of an application-specific integrated circuit built with this library with results obtained from software simulations. We then present the complex behavior of neural membrane voltages and the potential applications of this Si neuron.

Sensory Feedback in a Half-Center Oscillator Model

We hypothesize that one role of sensorimotor feedback for rhythmic movements in biological organisms is to synchronize the frequency of movements to the mechanical resonance of the body. Our hypothesis is based on recent studies that have shown the advantage of moving at mechanical resonance and how such synchronization may be possible in biology. We test our hypothesis by developing a physical system that consists of a silicon-neuron central pattern generator (CPG), which controls the motion of a beam, and position sensors that provide feedback information to the CPG. The silicon neurons that we use are integrated circuits that generate neural signals based on the Hodgkin-Huxley dynamics. We use this physical system to develop a model of the interaction between the sensory feedback and the complex dynamics of the neurons to create the closed-loop system behavior. This model is then used to describe the conditions under which our hypothesis is valid and the general effects of sensorimotor feedback on the rhythmic movements of this system.