Computational modeling of bursting pacemaker neurons in the pre-Bötzinger complex

Central control of interlimb coordination and speed-dependent gait expression in quadrupeds

KEY POINTS

Quadrupeds express different gaits depending on speed of locomotion. Central pattern generators (one per limb) within the spinal cord generate locomotor oscillations and control limb movements. Neural interactions between these generators define interlimb coordination and gait. We present a computational model of spinal circuits representing four rhythm generators with left-right excitatory and inhibitory commissural and fore-hind inhibitory interactions within the cord. Increasing brainstem drive to all rhythm generators and excitatory commissural interneurons induces an increasing frequency of locomotor oscillations accompanied by speed-dependent gait changes from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for multiple experimental data, including speed-dependent gait transitions in intact mice and changes in gait expression in mutants lacking certain types of commissural interneurons. The model suggests the possible circuit organization in the spinal cord and proposes predictions that can be tested experimentally.

ABSTRACT

As speed of locomotion is increasing, most quadrupeds, including mice, demonstrate sequential gait transitions from walk to trot and to gallop and bound. The neural mechanisms underlying these transitions are poorly understood. We propose that the speed-dependent expression of different gaits results from speed-dependent changes in the interactions between spinal circuits controlling different limbs and interlimb coordination. As a result, the expression of each gait depends on (1) left-right interactions within the spinal cord mediated by different commissural interneurons (CINs), (2) fore-hind interactions on each side of the spinal cord and (3) brainstem drives to rhythm-generating circuits and CIN pathways. We developed a computational model of spinal circuits consisting of four rhythm generators (RGs) with bilateral left-right interactions mediated by V0 CINs (V0_D and V0_V sub-types) providing left-right alternation, and conditional V3 CINs promoting left-right synchronization. Fore and hind RGs mutually inhibited each other. We demonstrate that linearly increasing excitatory drives to the RGs and V3 CINs can produce a progressive increase in the locomotor speed accompanied by sequential changes of gaits from walk to trot and to gallop and bound. The model closely reproduces and suggests explanations for the speed-dependent gait expression observed in vivo in intact mice and in mutants lacking V0_V or all V0 CINs. Specifically, trot is not expressed after removal of V0_V CINs, and only bound is expressed after removal of all V0 CINs. The model provides important insights into the organization of spinal circuits and neural control of locomotion.

Cycle-by-cycle assembly of respiratory network activity is dynamic and stochastic

Rhythmically active networks are typically composed of neurons that can be classified as silent, tonic spiking, or rhythmic bursting based on their intrinsic activity patterns. Within these networks, neurons are thought to discharge in distinct phase relationships with their overall network output, and it has been hypothesized that bursting pacemaker neurons may lead and potentially trigger cycle onsets. We used multielectrode recording from 72 experiments to test these ideas in rhythmically active slices containing the pre-Bötzinger complex, a region critical for breathing. Following synaptic blockade, respiratory neurons exhibited a gradient of intrinsic spiking to rhythmic bursting activities and thus defied an easy classification into bursting pacemaker and nonbursting categories. Features of their firing activity within the functional network were analyzed for correlation with subsequent rhythmic bursting in synaptic isolation. Higher firing rates through all phases of fictive respiration statistically predicted bursting pacemaker behavior. However, a cycle-by-cycle analysis indicated that respiratory neurons were stochastically activated with each burst. Intrinsically bursting pacemakers led some population bursts and followed others. This variability was not reproduced in traditional fully interconnected computational models, while sparsely connected network models reproduced these results both qualitatively and quantitatively. We hypothesize that pacemaker neurons do not act as clock-like drivers of the respiratory rhythm but rather play a flexible and dynamic role in the initiation and stabilization of each burst. Thus, at the behavioral level, each breath can be thought of as de novo assembly of a stochastic collaboration of network topology and intrinsic properties.

Toxin detection based on action potential shape analysis using a realistic mathematical model of differentiated NG108-15 cells

The NG108-15 neuroblastoma/glioma hybrid cell line has been frequently used for toxin detection, pharmaceutical screening and as a whole-cell biosensor. However, detailed analysis of its action potentials during toxin or drug administration has not been accomplished previously using patch clamp electrophysiology. In order to explore the possibility of identifying toxins based on their effect on the shape of intracellularly or extracellularly detected action potentials, we created a computer model of the action potential generation of this cell type. To generate the experimental data to validate the model, voltage dependent sodium, potassium and high-threshold calcium currents, as well as action potentials, were recorded from NG108-15 cells with conventional whole-cell patch-clamp methods. Based on the classic Hodgkin-Huxley formalism and the linear thermodynamic description of the rate constants, ion-channel parameters were estimated using an automatic fitting method. Utilizing the established parameters, action potentials were generated in the model and were optimized to represent the actual recorded action potentials to establish baseline conditions. To demonstrate the applicability of the method for toxin detection and discrimination, the effect of tetrodotoxin (a sodium channel blocker) and tefluthrin (a pyrethroid that is a sodium channel opener) were studied. The two toxins affected the shape of the action potentials differently and their respective effects were identified based on the changes in the fitted parameters. Our results represent one of the first steps to establish a complex model of NG108-15 cells for quantitative toxin detection based on action potential shape analysis of the experimental results.