Application of neural oscillators to study the effects of walking speed on rhythmic activations at the ankle
© Chong et al; licensee BioMed Central Ltd. 2013
Received: 15 May 2012
Accepted: 8 February 2013
Published: 13 February 2013
Spinal pattern generators (SPG) are neural networks in the spinal cord that do not require a central input from the brain to generate a motor output. We wanted to determine whether SPG can adapt to the changing motor demands from walking at different speeds, and performing silly walks.
An SPG model consisting of an oscillator made up of two neurons was utilised in this study; one neuron activates the soleus and the other activates the tibialis anterior. The outputs of the SPG model therefore represent the electromyographic measurements from each muscle. Seven healthy subjects were requested to perform silly walks, normal walking at self-selected speed (4.8 ± 0.5 km/h), 3.5 km/h, 4.0 km/h and 4.5 km/h on a treadmill. Loading and hip angles were used as inputs into the model.
No significant differences in the model parameters were found between normal walking at self-selected speed and other walking speeds. Only the adaptation time constant for the ankle flexor during silly walks was significantly different from the other normal walking trials.
We showed that SPG in the spinal cord can interpret and respond accordingly to velocity-dependent afferent information. Changes in walking speed do not require a different motor control mechanism provided there is no disruption to the alternating muscular activations generated at the ankle.
KeywordsLocomotion Walking speed Silly walks EMG Spinal pattern generator Neural network Sensory afferents Motor control
Afferents interact continuously with different parts of the nervous system so as to enable a smooth and efficient gait. As the human nervous system should coordinate efficiently, responding and adapting to the immediate environment, it is important that the plethora of signals coming from the central, sensory and peripheral systems be selected and modulated, so that the motor output fulfils the demands of the locomotor task.
Different control networks at different levels of the nervous system contribute to human motor control. The lowest level of neural control, which is responsible for generating the basic patterns of locomotion, is believed to come from spinal pattern generators (SPG) located in the spinal cord [1–3]. Brown  showed that decerebrate cats can produce locomotor-like electromyographic (EMG) patterns while walking on a treadmill. Similar results were also obtained from other vertebrates and invertebrates [1, 5]. A study by Maegele et al. in 2002  showed that clinically incompletely and completely spinal cord-injured patients can activate lower limb muscles after treadmill therapy. While these studies successfully showed the ability of the SPG to produce a motor output with no interference from the brain, they also demonstrated that the interactions between SPG and sensory inputs are important in generating a dynamic movement . Taga  had shown that a real-time dynamic interaction between the neural and mechanical system, together with sensory information from the environment, could influence the motor output of the lower limbs.
Walking at a slower or faster pace creates different motor demands on the neural system. A number of gait components such as stance and swing phase intervals and muscle activations change with increasing speed. However, in healthy humans, it is not known whether these changes result from sensory cues to the neural network in the spinal cord, since higher commands from the brain can intervene. The aim of our study was to determine whether neural networks in the spinal cord can adapt to changing sensory afferents, and directly influence muscular activity to meet the motor demands of walking at different speeds. We therefore studied the response of the SPG model in situations where gait components in a gait cycle will be different from normal walking: change in walking speed and in performing “silly walks”. Our study used an SPG model that is triggered only by sensory afferents with no interference from a cortical signal.
Seven healthy male subjects (28.0 ± 4.4 years, 1.8 ± 0.1 m, 76.4 ± 9.5 kg) volunteered to participate in this study. They were thoroughly informed of the procedures and gave their consent. Each subject was requested to walk at his normal self-selected speed (4.8 ± 0.5 km/h), at 3.5 km/h, 4.0 km/h and 4.5 km/h on a treadmill (Kinetics s3, Kettler, Germany). In addition, they were asked to perform movements unlike their normal walking, i.e. “silly walks”, at a speed of their own choice (3.8 ± 0.4 km/h). Data from six consecutive strides were collected during steady-state walking. Three trials were recorded for each subject for each walking speed and silly walks, i.e. a total of 210 trials (for both right and left limbs). Trials were ignored if there were missing data in any one stride. Therefore, a total of only 176 trials was analysed in this study.
Vertical force data were collected at 200Hz and calculated from in-shoe pressure sensors (Gesellschaft für Biomechanik Münster, Germany) as a summation of the pressure acting on the entire area of the insole. Hip angles were acquired from an Oqus 3D motion analysis system (Qualisys, Sweden) at 100 Hz. This system used six infra-red cameras, which tracked a total of fifteen retro-reflective markers attached to the following body landmarks: lateral and medial knee, and four tracking markers on the thigh of each leg, sacrum, and left and right anterior superior iliac spine. Segment definitions and kinematic data were processed using Visual3D (C-Motion Inc, Maryland, USA). Muscle activation from the m. soleus (SOL) and m. tibialis anterior (TA) were captured using bipolar surface electrodes (5–700 Hz, Biovision, Wehrheim, Germany) at 2000 Hz. The SOL and TA muscles were chosen because they are the principal monoarticular plantarflexor and dorsiflexor muscles, respectively. Electrodes were placed according to recommendations by Hermens et al., 1999 . The electromyographic (EMG) signals were centred, rectified and filtered using a fifth-order low-pass Butterworth filter with a cut-off frequency of 40Hz.
The parameters a, b, m, n, r, T, v, w from the above equations determined the pattern and frequency of the output. A nonlinear least-squares fitting algorithm was used to determine a set of parameters that would fit the output to experimental data i.e. the output produced by the neuron representing the SOL would be fitted to measured EMG data of the SOL. This was done simultaneously for the TA. Initial values for each neuron were taken from the first value of the measured EMG data so as to solve the differential equations numerically. The fitting algorithm terminated once the relative deviation between two iterations fell below 0.001. A correlation coefficient R between the model output and experimental EMG data was calculated in each trial.
The following gait components were analysed; Maximum normalised force and maximum range of hip flexion-extension angles were calculated for each stride. Stance and swing phases determined from force profiles of each stride were also calculated. For these gait components, analysis of variance (ANOVA) and Tukey’s post-hoc test were performed to determine the significant differences between all the different walking types. In analysing the rectified EMG signals for different speeds, we adopted the method by Murray et al., 1984 . Cumulative numerical integration (IEMG) for each EMG signal in each stride was calculated for all speeds. The maximum of the mean IEMG values was designated as 1.00, regardless of speed. The other mean values were normalised with respect to this maximum value . To determine significant differences in the model parameters (p < 0.05), multivariate analysis of variance (MANOVA) along with analysis of variance (ANOVA) and Tukey’s post-hoc test were performed.
Mean and standard deviation (std) of gait components at different speeds and during silly walks
Self-selected 4.8 ± 0.5 km/h
Silly walks 3.8 ± 0.4 km/h
p < 0.05
p < 0.05
Hip flexion-extension range (rad)
p < 0.05
p < 0.05
p < 0.05
p < 0.05
The results of this study showed that neural networks in the spinal cord can activate muscles at the ankle to generate stepping motion during steady-state walking. In humans, it is difficult to determine whether the elevated EMG patterns during walking result from supraspinal control, activations from sensory inputs, or an interaction between supraspinal and spinal control. However, the outputs generated by our model, consisting only of spinal neurons, suggest that muscle activations can be generated by sensory inputs from loading and hip angles at the spinal level .
While there were significant differences in both the inputs (loading and hip angles) and outputs (IEMG) between the different walking speeds and silly walks (Table 1), significant differences were only found in T 2 between silly walks and the other walking tests (Figure 6). For normal walking at different speeds, this might imply that an insignificant change in a parameter is sufficient to cause a significant change in the output. Since the control of these parameters, which determine the neuronal properties of the SPG, could come from interneurons, presynaptic inhibition , or through descending pathways from supraspinal structures, the insignificant changes might imply that no regulation by the brain or inter-spinal circuitry is required to modulate the activation patterns during walking. Comparing our study to split-belt treadmill locomotion, Morton and Bastian (2006) found that subjects with cerebellar damage were able to perform rapid reactive adjustments to stride length and stance time when their legs were operating at different speeds . Our study corroborates their findings that higher control is not needed to alter the motor output of moving limbs, but the corrections could instead be performed predominantly by spinal structures using available sensory information. It has also been suggested that the same neural circuitry is responsible for gait transitions between walking and running . Thus, while an input from the cerebral cortex is required to initiate a movement, higher command centres need not be recruited to regulate motor output during locomotion regardless of speed.
Grasso et al.  suggested that the nervous system attempts to meet motor demands by controlling posture or limb joint motion rather than regulating muscle activations. We agree with their arguments, since we successfully used loading and hip angles as inputs to the SPG model to generate muscle activations. In addition, provided the gait patterns do not result in changes to equilibrium, the same neural network will be utilised . Since the data were captured during steady-state walking, it might also be important that the alternating activations of the flexor and extensor are not disrupted. Perhaps, changes to gait components are secondary, and could result from changes in stride length rather than a different motor control mechanism.
It has been shown that cats with lesions in the motor cortex encountered no problems walking on a flat horizontal surface until they were required to cross obstacles or climb a ladder . In addition, Armstrong and Drew  found that pulses measured in the cat’s cortical neurons were unrelated to speed, though muscle activity increased significantly. Therefore, as in humans, no conscious effort is necessary during level walking regardless of walking pace, until an obstacle or a sudden change in the external environment is encountered, when corrective responses are required.
Additional file 1: One subject performing a silly walk.(AVI 17 MB)
A limitation in this study was the restricted array of silly walks the subjects could perform while walking on a treadmill at a constant speed (an example is shown in ‘Additional file 1’). The movements performed by the subjects still involved an on-going, uninterrupted rhythmic pattern of activation between the antagonistic muscles at the ankle. Since we now know the same neural network is responsible for normal walking at different speeds, future studies can give the subjects a freer choice of the types of silly walks they would like to perform (like those seen in Monty Python’s sketch, The Ministry of Silly Walks). In such studies, significant differences in more model parameters might be found.
We proposed that SPG in the spinal cord can interpret and respond accordingly to velocity-dependent afferent information. Changes in walking speed do not require a different motor control mechanism provided equilibrium is not affected and there is no disruption of the continuous rhythmic patterns produced at the ankle.
Spinal pattern generators
M. tibialis anterior
Determine time course of adaptation in neuron
Adaptation in neuron
Weight of hip angles input
Weight of change of hip angles
Smoothed hip angles
Determine time course of adaptation in neuron
Weight of change of force input
Weight of force input
Inner state of neuron
- Max F:
Cumulative numerical integration of EMG signals
Cumulative numerical integration of soleus EMG signals
Cumulative numerical integration of tibialis anterior EMG signals
Multivariate analysis of variance
Analysis of variance
Persistent inward currents
This study was supported by the German Federal Ministry of Education and Research (BMBF), project no. 01EC1003A. Technical assistance was provided by Hanno Focke and David Schinowski. We acknowledge support by the German Research Foundation (DFG) and Open Access Publication Fund of University of Muenster.
- Grillner S, Wallen P: Central pattern generators for locomotion, with special reference to verterbrates. Ann Rev Neurosci. 1985, 8: 233-61. 10.1146/annurev.ne.08.030185.001313.View ArticlePubMed
- Ijspeert A: Central pattern generators for locomotion control in animals and robots: a review. Neural Netw. 2008, 21: 642-53. 10.1016/j.neunet.2008.03.014.View ArticlePubMed
- Rossignol S, Frigon A: Recovery of locomotion after spinal cord injury: some facts and mechanisms. Annu Rev Neurosci. 2011, 34: 413-40. 10.1146/annurev-neuro-061010-113746.View ArticlePubMed
- Brown TG: On the nature of the fundamental activity of the nervous centres; together with an analysis of the conditioning of rhythmic activity in progression, and a theory of the evolution of function in the nervous system. J Physiol-London. 1914, 48: 18-46.PubMed CentralView ArticlePubMed
- Marder E, Bucher D: Central pattern generators and the control of rhythmic movements. Curr Biol. 2001, 11: R986-96. 10.1016/S0960-9822(01)00581-4.View ArticlePubMed
- Maegele M, Müller S, Wernig A, Edgerton VR, Harkema SJ: Recruitment of spinal motor pools during voluntary movements versus stepping after human spinal cord injury. J Neurotrauma. 2002, 19 (10): 1217-29. 10.1089/08977150260338010.View ArticlePubMed
- Baessler U: On the definition of central pattern generator and its sensory control. Biol Cybern. 1986, 54: 65-9. 10.1007/BF00337116.View Article
- Taga G: A model of the neuro-musculo-skeletal system for human locomotion. i. emergence of basic gait. Biol Cybern. 1995, 73: 97-111. 10.1007/BF00204048.View ArticlePubMed
- Hermens HJ, Freriks B, Merletti R, Stegeman DF, Blok J, Rau G, Disselhorst-Klug C, Haegg G: European Recommendations For Surface Electromyography, Results Of The SENIAM Project. 1999, Enschede: Roessingh Research and Development
- Matsuoka K: Sustained oscillations generated by mutually inhibiting neurons with adaptation. Biol Cybern. 1985, 52: 367-76. 10.1007/BF00449593.View ArticlePubMed
- Matsuoka K: Mechanisms of frequency and pattern control in the neural rhythm generators. Biol Cybern. 1987, 56: 345-53. 10.1007/BF00319514.View ArticlePubMed
- Murray MP, Mollinger LA, Gardner GM, Sepic SB: Kinematic and EMG patterns during slow, free, and fast walking. J Orthop Res. 1984, 2: 272-80. 10.1002/jor.1100020309.View ArticlePubMed
- Chong SY, Wagner H, Wulf A: Neural oscillators triggered by loading and hip orientation can generate activation patterns at the ankle during walking in humans. Med Biol Eng Comput. 2012, 50 (9): 917-23. 10.1007/s11517-012-0944-2.View ArticlePubMed
- Morton SM, Bastian AJ: Cerebellar contributions to locomotor adaptations during splitbelt treadmill walking. J Neuroscience. 2006, 26: 9107-16. 10.1523/JNEUROSCI.2622-06.2006.View ArticlePubMed
- Labini FS, Ivanenko YP, Cappellini G, Gravano S, Lacquaniti F: Smooth changes in the EMG patterns during gait transitions under body weight unloading. J Neurophysiol. 2011, 106: 1525-36. 10.1152/jn.00160.2011.View ArticlePubMed
- Grasso R, Bianchi L, Lacquaniti F: Motor patterns for human gait: backward versus forward locomotion. J Neurophysiol. 1998, 80: 1868-85.PubMed
- Lamb T, Yang JF: Could different directions of infant stepping be controlled by the same locomotor centre pattern generator?. J Neurophysiol. 2000, 83: 2814-24.PubMed
- Armstrong DM: The Supraspinal control of mammalian locomotion. J Physiol. 1988, 405: 1-37.PubMed CentralView ArticlePubMed
- Armstrong DM, Drew TL: Discharges of pyramidal tract and other motor cortical neurones during locomotion in the cat. J Physiol. 1984, 346: 471-95.PubMed CentralView ArticlePubMed
- Li Y, Gorassini MA, Bennett DJ: Role of persistent sodium and calcium currents in motoneuron firing and spasticity in chronic spinal rats. J Neurophysiol. 2004, 91: 767-83.View ArticlePubMed
- Cotel F, Antri M, Barthe JY, Orsal D: Identified ankle extensor and flexor motoneurons display different firing profiles in the neonatal rat. J Neurosci. 2009, 29: 2748-53. 10.1523/JNEUROSCI.3462-08.2009.View ArticlePubMed
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.