Recent experimental and computational studies suggest that linearly correlated sets of parameters (intrinsic and synaptic properties of neurons) allow central pattern-generating networks to produce and maintain their rhythmic activity regardless of changing internal and external conditions. To determine the role of correlated conductances in the robust maintenance of functional bursting activity, we used our existing database of half-center oscillator (HCO) model instances of the leech heartbeat CPG. From the database, we identified functional activity groups of burster (isolated neuron) and half-center oscillator model instances and realistic subgroups of each that showed burst characteristics (principally period and spike frequency) similar to the animal. To find linear correlations among the conductance parameters maintaining functional leech bursting activity, we applied Principal Component Analysis (PCA) to each of these four groups. PCA identified a set of three maximal conductances (leak current, [Formula: see text]Leak; a persistent K current, [Formula: see text]K2; and of a persistent Na+ current, [Formula: see text]P) that correlate linearly for the two groups of burster instances but not for the HCO groups. Visualizations of HCO instances in a reduced space suggested that there might be non-linear relationships between these parameters for these instances. Experimental studies have shown that period is a key attribute influenced by modulatory inputs and temperature variations in heart interneurons. Thus, we explored the sensitivity of period to changes in maximal conductances of [Formula: see text]Leak, [Formula: see text]K2, and [Formula: see text]P, and we found that for our realistic bursters the effect of these parameters on period could not be assessed because when varied individually bursting activity was not maintained.

Flexibility in neuronal circuits has its roots in the dynamical richness of their neurons. Depending on their membrane properties single neurons can produce a plethora of activity regimes including silence, spiking and bursting. What is less appreciated is that these regimes can coexist with each other so that a transient stimulus can cause persistent change in the activity of a given neuron. Such multistability of the neuronal dynamics has been shown in a variety of neurons under different modulatory conditions. It can play either a functional role or present a substrate for dynamical diseases. We considered a database of an isolated leech heart interneuron model that can display silent, tonic spiking and bursting regimes. We analyzed only the cases of endogenous bursters producing functional half-center oscillators (HCOs). Using a one parameter (the leak conductance (g(leak)) bifurcation analysis, we extended the database to include silent regimes (stationary states) and systematically classified cases for the coexistence of silent and bursting regimes. We showed that different cases could exhibit two stable depolarized stationary states and two hyperpolarized stationary states in addition to various spiking and bursting regimes. We analyzed all cases of endogenous bursters and found that 18% of the cases were multistable, exhibiting coexistences of stationary states and bursting. Moreover, 91% of the cases exhibited multistability in some range of g(leak). We also explored HCOs built of multistable neuron cases with coexisting stationary states and a bursting regime. In 96% of cases analyzed, the HCOs resumed normal alternating bursting after one of the neurons was reset to a stationary state, proving themselves robust against this perturbation.

Robust activity of some networks, such as central pattern generators, suggests the existence of physiological mechanisms that maintain the most important characteristics, for example, the period and spike frequency of the pattern. Whatever these mechanisms are, they change the appropriate model parameters to or along the isomanifolds on which the characteristics of the pattern are constant, while their sensitivities to parameters may be different. Setting synaptic connections to zero at the points of isomanifolds allows for dissecting the maintenance mechanisms into components involving synaptic transmission and components involving intrinsic currents. The physiological meaning of the intrinsic current changes might be revealed by analysis of their impact on endogenous neuronal dynamics. Here, we sought answers to two questions: (1) Do parameter variations in insensitive directions (along isomanifolds) change endogenous dynamics of the network neurons? (2) Do sensitive and insensitive directions for network pattern characteristics depend on endogenous dynamics of the network neurons? We considered a leech heartbeat half-center oscillator model network and analyzed isomanifolds on which the burst period or spike frequency of the model, or both, are constant. Based on our analysis, we hypothesize that the dependence on endogenous dynamics of the isolated neurons is the stronger the more characteristics of the pattern have to be maintained. We also found that in general, the network was more flexible when it consisted of endogenously tonically spiking rather than bursting or silent neurons. Finally, we discuss the physiological implications of our findings.

We developed an analog very large-scale integrated system of two mutually inhibitory silicon neurons that display several different stable oscillations. For example, oscillations can be synchronous with weak inhibitory coupling and alternating with relatively strong inhibitory coupling. All oscillations observed experimentally were predicted by bifurcation analysis of a corresponding mathematical model. The synchronous oscillations do not require special synaptic properties and are apparently robust enough to survive the variability and constraints inherent in this physical system. In biological experiments with oscillatory neuronal networks, blockade of inhibitory synaptic coupling can sometimes lead to synchronous oscillations. An example of this phenomenon is the transition from alternating to synchronous bursting in the swimming central pattern generator of lamprey when synaptic inhibition is blocked by strychnine. Our results suggest a simple explanation for the observed oscillatory transitions in the lamprey central pattern generator network: that inhibitory connectivity alone is sufficient to produce the observed transition.