Skip to content


You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Neuroscience

Open Access

Mathematical modeling and analysis of spinal circuits involved in locomotor pattern generation and frequency-dependent left-right coordination

  • Yaroslav I Molkov1Email author,
  • Bartholomew J Bacak2 and
  • Ilya A Rybak2
BMC Neuroscience201415(Suppl 1):P53

Published: 21 July 2014

Coordination between left and right neural activities in the spinal cord during locomotion is controlled by commissural interneurons (CINs). Several CIN types have been genetically identified, including the excitatory V3 and excitatory and inhibitory V0 types. Talpalar et al. [1] recently reported that genetic elimination of the V0 CINs caused switching from a normal left-right alternating pattern of motor activity to a left-right synchronized “hopping” pattern. Furthermore, ablation of only the inhibitory V0 neurons (V0D subtype) resulted in a lack of left-right alternation at low locomotor frequencies and maintaining this alternation at high frequencies, whereas selective ablation of the excitatory V0 neurons (V0V subtype) maintained the left–right alternation at low locomotor frequencies and switched the motor output to a left-right synchronized ("hopping") pattern at high frequencies.

To analyze and explain the above findings, we developed a simplified mathematical model of neural circuits consisted of four pacemaker neurons representing left (LF) and right (RF) flexor and left (LE) and right (RE) extensor half-centers interacting via commissural pathways representing V3, V0D, and V0V CINs (Fig. 1). The “locomotor” frequency in the model was controlled by a parameter defining the excitability of neurons (via the leak reversal potentials) and commissural pathways, whose changes represented the corresponding changes induced by changing the concentration of N-methyl-D-aspartate (NMDA) applied to control the locomotor frequency in the isolated rodent spinal cord preparations [1].
Figure 1

Model schematic.

The model demonstrated: (1) a typical left-right alternating pattern under control conditions; (2) switching to a synchronized hopping activity at any frequency after removing commissural connections representing both V0 (V0D and V0V) neurons; (3) a synchronized left-right pattern at low frequencies and left-right alternation of activity at high frequencies after removing the commissural connections representing V0D neurons; (4) an alternating left-right pattern at low frequencies with synchronized hopping at high frequencies after removing the commissural connections representing V0V neurons. We used the bifurcation and fast–slow decomposition methods to analyze the behavior of this network in the four above states/regimes and transitions between them. The model was able to reproduce and provide explanation to several important experimental phenomena and generated predictions that can be tested experimentally.



Supported by NIH grant R01NS7323.

Authors’ Affiliations

Department of Mathematical Sciences, Indiana University – Purdue University Indianapolis
Department of Neurobiology and Anatomy, Drexel University College of Medicine


  1. Talpalar AE, Bouvier J, Borgius L, Fortin G, Pierani A, Kiehn O: Dual-mode operation of neuronal networks involved in left-right alternation. Nature. 2013, 500 (7460): 85-93. 10.1038/nature12286.View ArticlePubMedGoogle Scholar


© Molkov et al; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.