Different gaits of locomotion in mammals are based on the appropriate coordination of neuronal activity in the spinal cord controlling movements of left and right limbs. This left-right coordination is provided in the spinal cord by the commissural interneurons (CINs) whose axons cross the midline and affect neural circuits on the contralateral side of the cord. Several types of CINs have been genetically identified, including the excitatory V3 CINs and the inhibitory (V0D) and excitatory (V0V) V0 CINs. Talpalar et al. [1] recently demonstrated that (a) ablation of both V0 CIN types leads to a left-right synchronized, “hopping” activity at all locomotor frequencies, whereas (b) selective ablation of the excitatory V0V CINs maintains alternation at low frequencies but switches to synchronized activity at high frequencies while (c) ablation of only the inhibitory V0D CINs leads to a lack of left–right alternation at low frequencies, but maintains alternation at high frequencies. The genetically identified, ipsilaterally projecting excitatory V2a interneurons are recruited with an increase in locomotor speed [2] and contribute to left-right alternation at high locomotor frequencies [3,4]. Our objective was to construct and analyze a computational model of the bilaterally interacting central pattern generators (CPGs) that could reproduce and explain these findings. In our model, the CPG on each side of the cord consisted of flexor and extensor half-centers. Each neural population, including the half-centers and CIN populations, consisted of 50-200 neurons modeled in the Hodgkin-Huxley style. The intrinsic bursting of CPG neurons was based on a persistent sodium current in these neurons. During model construction, we assumed that left-right coordination of activity depends on the balance between the three CIN pathways providing interactions between the left and right CPGs: the V3-mediated pathway that supports left-right synchronization and the V0D- and V0V- mediated pathways that provide left-right alternation. The activity of each (left and right) inhibitory V0D population was driven by the ipsilateral flexor half-center. The recruitment of V0D neurons was progressively reduced with an increase in locomotor speed, because of the reduction of burst amplitude. The left and right V0V pathways could be organized in two ways: (1) the V0V activity on each side was driven by the ipsilateral flexor half-center and its action on contralateral circuits was mediated by an inhibitory population, or (2) the V0V activity was driven by the ipsilateral extensor half-center and it excited the contralateral circuits. In any case, the V0V activation was mediated by the ipsilateral V2a neurons progressively recruited with increasing locomotor speed. The model demonstrates: (1) a left-right alternating pattern under control conditions; (2) a synchronized hopping pattern at any frequency after removing both the V0 populations; (3) a synchronized pattern at low frequencies with alternation at high frequencies after removing the V0D populations; (4) an alternating pattern at low frequencies with synchronized hopping at high frequencies after removing either V0V or V2a populations. The model closely reproduces and suggests an explanation for the experimental data of Talpalar et al. [1], Zhong et al. [2], and Crone et al. [3,4], proposes the organization of commissural interactions in the spinal cord defining the left-right alternation at different locomotor speeds, and generates predictions for future experimental investigations.
Read full abstract