Although the neural basis of balance control remains unknown, recent studies suggest that a feedback law on center-of-mass (CoM) kinematics determines the temporal patterning of muscle activity during human postural responses. We hypothesized that the same feedback law would also explain variations in muscle activity to support-surface translation as perturbation characteristics vary. Subject CoM motion was experimentally modulated using 34 different anterior–posterior support-surface translations of varying peak acceleration and velocity but the same total displacement. Electromyographic (EMG) recordings from several muscles of the lower limbs and trunk were compared to predicted EMG patterns from an inverted pendulum model under delayed feedback control. In both recorded and predicted EMG patterns, the initial burst of muscle activity scaled linearly with peak acceleration, whereas the tonic “plateau” region scaled with peak velocity. The relatively invariant duration of the initial burst was modeled by incorporating a transient, time-limited encoding of CoM acceleration inspired by muscle spindle primary afferent dynamic responses. The entire time course of recorded and predicted muscle activity compared favorably across all conditions, suggesting that the initial burst of muscle activity is not generated by feedforward neural mechanisms. Perturbation conditions were presented randomly and subjects maintained relatively constant feedback gains across all conditions. In contrast, an optimal feedback solution based on a trade-off between CoM stabilization and energy expenditure predicted that feedback gains should change with perturbation characteristics. These results suggest that an invariant feedback law was used to generate the entire time course of muscle activity across a variety of postural disturbances.

We developed wavelet-based functional ANOVA (wfANOVA) as a novel approach for comparing neurophysiological signals that are functions of time. Temporal resolution is often sacrificed by analyzing such data in large time bins, increasing statistical power by reducing the number of comparisons. We performed ANOVA in the wavelet domain because differences between curves tend to be represented by a few temporally localized wavelets, which we transformed back to the time domain for visualization. We compared wfANOVA and ANOVA performed in the time domain (tANOVA) on both experimental electromyographic (EMG) signals from responses to perturbation during standing balance across changes in peak perturbation acceleration (3 levels) and velocity (4 levels) and on simulated data with known contrasts. In experimental EMG data, wfANOVA revealed the continuous shape and magnitude of significant differences over time without a priori selection of time bins. However, tANOVA revealed only the largest differences at discontinuous time points, resulting in features with later onsets and shorter durations than those identified using wfANOVA (P < 0.02). Furthermore, wfANOVA required significantly fewer (∼¼×; P < 0.015) significant F tests than tANOVA, resulting in post hoc tests with increased power. In simulated EMG data, wfANOVA identified known contrast curves with a high level of precision (r2 = 0.94 ± 0.08) and performed better than tANOVA across noise levels (P < <0.01). Therefore, wfANOVA may be useful for revealing differences in the shape and magnitude of neurophysiological signals (e.g., EMG, firing rates) across multiple conditions with both high temporal resolution and high statistical power.

We hypothesized that motor outputs are hierarchically organized such that descending temporal commands based on desired task-level goals flexibly recruit muscle synergies that specify the spatial patterns of muscle coordination that allow the task to be achieved. According to this hypothesis, it should be possible to predict the patterns of muscle synergy recruitment based on task-level goals. We demonstrated that the temporal recruitment of muscle synergies during standing balance control was robustly predicted across multiple perturbation directions based on delayed sensorimotor feedback of center of mass (CoM) kinematics (displacement, velocity, and acceleration). The modulation of a muscle synergy's recruitment amplitude across perturbation directions was predicted by the projection of CoM kinematic variables along the preferred tuning direction(s), generating cosine tuning functions. Moreover, these findings were robust in biphasic perturbations that initially imposed a perturbation in the sagittal plane and then, before sagittal balance was recovered, perturbed the body in multiple directions. Therefore, biphasic perturbations caused the initial state of the CoM to differ from the desired state, and muscle synergy recruitment was predicted based on the error between the actual and desired upright state of the CoM. These results demonstrate that that temporal motor commands to muscle synergies reflect task-relevant error as opposed to sensory inflow. The proposed hierarchical framework may represent a common principle of motor control across motor tasks and levels of the nervous system, allowing motor intentions to be transformed into motor actions.

During postural responses to perturbations, horizontal plane forces generated by the cat hindlimb are stereotypically directed either towards or away from the animal’s center of mass, independent of perturbation direction. We used a static, three-dimensional musculoskeletal model of the hindlimb to investigate possible biomechanical determinants of this “force constraint strategy” (Macpherson, 1988). We hypothesized that directions in which the hindlimb can produce large forces are preferentially used in postural control. We computed feasible force sets (FFS) based on hindlimb configurations of three cats during postural equilibrium tasks (Jacobs and Macpherson, 1996) and compared them to horizontal plane postural force directions. The grand mean FFS was bimodal, with maxima near the posterior-anterior axis (−86±8° and 71±4°), and minima near the medial-lateral axis (177±8° and 8±8°). Postural force directions clustered near both maxima; there were no medial postural forces near the absolute minimum. However, the medians of the anterior and posterior postural force direction histograms in the right hindlimb were rotated counter-clockwise from the FFS maxima (p<0.05; Wilcoxon signed-rank test). Because the posterior-anterior alignment of the FFS is consistent with a hindlimb structure optimized for locomotion, we conclude that the biomechanical capabilities of the hindlimb strongly influence, but do not uniquely determine the force directions observed in the force constraint strategy. Forces used in postural control may reflect a balance between a neural preference for using forces in the directions of large feasible forces and other criteria, such as the stabilization of the center of mass, and muscular coordination strategies.

The pendulum test is a sensitive clinical assessment of spasticity where the lower leg is dropped from the horizontal position and features of limb motion are recorded. Three key kinematic features are associated with the degree of severity of spasticity in children with cerebral palsy: decreased initial limb excursion, reduced number of limb oscillations, and a non-vertical resting limb angle. While spasticity is attributed to increased velocity-dependent resistance to motion, prior models simulating increased sensorimotor feedback of muscle velocity fail to explain the key pendulum test kinematic outcomes in spastic individuals. Here we hypothesized that increased muscle tone, causing a transient increase in muscle force, i.e. short-range stiffness, could account for reduced first swing excursion and non-vertical resting limb angle. We further hypothesized that hyperreflexia modeled based on muscle fiber force, and not velocity, feedback would be necessary to reduce the number of oscillations because of its interaction with transiently increased muscle force due to short-range stiffness. We simulated the lower leg as a torque-driven single-link pendulum. Muscle tone was modeled as a constant baseline joint torque, short-range stiffness torque was dependent on the level of muscle tone, and delayed sensory feedback torque to simulate reflex activity was based on either muscle velocity or force. Muscle tone and transient short-range stiffness were necessary to simulate decreased initial swing excursion and non-vertical resting leg angle. Moreover, the reduction in the number of oscillations was best reproduced by simulating stretch reflex activity in terms of force, and not velocity, feedback. Varying only baseline muscle torque and reflex gain, we simulated a range of pendulum test kinematics observed across different levels of spasticity. Our model lends insight into physiological mechanisms of spasticity whose contributions can vary on an individual-specific basis, and potentially across different neurological disorders that manifest spasticity as a symptom.

Neuromusculoskeletal models solve the basic problem of determining how the body moves under the influence of external and internal forces. Existing biomechanical modeling programs often emphasize dynamics with the goal of finding a feed-forward neural program to replicate experimental data or of estimating force contributions or individual muscles. The computation of rigid-body dynamics, muscle forces, and activation of the muscles are often performed separately. We have developed an intrinsically forward computational platform (Neuromechanic, www.neuromechanic.com) that explicitly represents the interdependencies among rigid body dynamics, frictional contact, muscle mechanics, and neural control modules. This formulation has significant advantages for optimization and forward simulation, particularly with application to neural controllers with feedback or regulatory features. Explicit inclusion of all state dependencies allows calculation of system derivatives with respect to kinematic states and muscle and neural control states, thus affording a wealth of analytical tools, including linearization, stability analyses and calculation of initial conditions for forward simulations. In this review, we describe our algorithm for generating state equations and explain how they may be used in integration, linearization, and stability analysis tools to provide structural insights into the neural control of movement.

The musculoskeletal redundancy of the body provides multiple solutions for performing motor tasks. We have proposed that the nervous system solves this unconstrained problem through the recruitment of motor modules or functional muscle synergies that map motor intention to action. Consistent with this hypothesis, we showed that trial-by-trial variations in muscle activation for multidirectional balance control in humans were constrained by a small set of muscle synergies. However, apparent muscle synergy structures could arise from characteristic patterns of sensory input resulting from perturbations or from low-dimensional optimal motor solutions. Here we studied electromyographic (EMG) responses for balance control across a range of biomechanical contexts, which alter not only the sensory inflow generated by postural perturbations, but also the muscle activation patterns used to restore balance. Support-surface translations in 12 directions were delivered to subjects standing in six different postural configurations: one-leg, narrow, wide, very wide, crouched, and normal stance. Muscle synergies were extracted from each condition using nonnegative matrix factorization. In addition, muscle synergies from the normal stance condition were used to reconstruct muscle activation patterns across all stance conditions. A consistent set of muscle synergies were recruited by each subject across conditions. When balance demands were extremely different from the normal stance (e.g., one-legged or crouched stance), task-specific muscle synergies were recruited in addition to the preexisting ones, rather generating de novo muscle synergies. Taken together, our results suggest that muscle synergies represent consistent motor modules that map intention to action, regardless of the biomechanical context of the task.

Here we show that novel, energy-recycling stairs reduce the amount of work required for humans to both ascend and descend stairs. Our low-power, interactive, and modular steps can be placed on existing staircases, storing energy during stair descent and returning that energy to the user during stair ascent. Energy is recycled through event-triggered latching and unlatching of passive springs without the use of powered actuators. When ascending the energy-recycling stairs, naive users generated 17.4 ± 6.9% less positive work with their leading legs compared to conventional stairs, with the knee joint positive work reduced by 37.7 ± 10.5%. Users also generated 21.9 ± 17.8% less negative work with their trailing legs during stair descent, with ankle joint negative work reduced by 26.0 ± 15.9%. Our low-power energy-recycling stairs have the potential to assist people with mobility impairments during stair negotiation on existing staircases.

Simulating realistic musculoskeletal dynamics is critical to understanding neural control of muscle activity evoked in sensorimotor feedback responses that have inherent neural transmission delays. Thus, the initial mechanical response of muscles to perturbations in the absence of any change in muscle activity determines which corrective neural responses are required to stabilize body posture. Muscle short-range stiffness, a history-dependent property of muscle that causes a rapid and transient rise in muscle force upon stretch, likely affects musculoskeletal dynamics in the initial mechanical response to perturbations. Here we identified the contributions of short-range stiffness to joint torques and angles in the initial mechanical response to support surface translations using dynamic simulation. We developed a dynamic model of muscle short-range stiffness to augment a Hill-type muscle model. Our simulations show that short-range stiffness can provide stability against external perturbations during the neuromechanical response delay. Assuming constant muscle activation during the initial mechanical response, including muscle short-range stiffness was necessary to account for the rapid rise in experimental sagittal plane knee and hip joint torques that occurs simultaneously with very small changes in joint angles and reduced root mean square errors between simulated and experimental torques by 56% and 47%, respectively. Moreover, forward simulations lacking short-range stiffness produced unreasonably large joint angle changes during the initial response. Using muscle models accounting for short-range stiffness along with other aspects of history-dependent muscle dynamics may be important to advance our ability to simulate inherently unstable human movements based on principles of neural control and biomechanics.

The biomechanical principles underlying the organization of muscle activation patterns during standing balance are poorly understood. The goal of this study was to understand the influence of biomechanical inter-joint coupling on endpoint forces and accelerations induced by the activation of individual muscles during postural tasks. We calculated induced endpoint forces and accelerations of 31 muscles in a 7 degree-of-freedom, three-dimensional model of the cat hindlimb. To test the effects of inter-joint coupling, we systematically immobilized the joints (excluded kinematic degrees of freedom) and evaluated how the endpoint force and acceleration directions changed for each muscle in 7 different conditions. We hypothesized that altered inter-joint coupling due to joint immobilization of remote joints would substantially change the induced directions of endpoint force and acceleration of individual muscles. Our results show that for most muscles crossing the knee or the hip, joint immobilization altered the endpoint force or acceleration direction by more than 90° in the dorsal and sagittal planes. Induced endpoint forces were typically consistent with behaviorally observed forces only when the ankle was immobilized. We then activated a proximal muscle simultaneous with an ankle torque of varying magnitude, which demonstrated that the resulting endpoint force or acceleration direction is modulated by the magnitude of the ankle torque. We argue that this simple manipulation can lend insight into the functional effects of co-activating muscles. We conclude that inter-joint coupling may be an essential biomechanical principle underlying the coordination of proximal and distal muscles to produce functional endpoint actions during motor tasks.