Despite decades of research, we lack a mechanistic framework capable of predicting how movement-related signals are transformed into the diversity of muscle spindle afferent firing patterns observed experimentally, particularly in naturalistic behaviors. Here, a biophysical model demonstrates that well-known firing characteristics of mammalian muscle spindle la afferents-including movement history dependence, and nonlinear scaling with muscle stretch velocity-emerge from first principles of muscle contractile mechanics. Further, mechanical interactions of the muscle spindle with muscle-tendon dynamics reveal how motor commands to the muscle (alpha drive) versus muscle spindle (gamma drive) can cause highly variable and complex activity during active muscle contraction and muscle stretch that defy simple explanation. Depending on the neuromechanical conditions, the muscle spindle model output appears to ‘encode’ aspects of muscle force, yank, length, stiffness, velocity, and/or acceleration, providing an extendable, multiscale, biophysical framework for understanding and predicting proprioceptive sensory signals in health and disease.

Optimal control simulations have shown that both musculoskeletal dynamics and physiological noise are important determinants of movement. However, due to the limited efficiency of available computational tools, deterministic simulations of movement focus on accurately modelling the musculoskeletal system while neglecting physiological noise, and stochastic simulations account for noise while simplifying the dynamics. We took advantage of recent approaches where stochastic optimal control problems are approximated using deterministic optimal control problems, which can be solved efficiently using direct collocation. We were thus able to extend predictions of stochastic optimal control as a theory of motor coordination to include muscle coordination and movement patterns emerging from non-linear musculoskeletal dynamics. In stochastic optimal control simulations of human standing balance, we demonstrated that the inclusion of muscle dynamics can predict muscle co-contraction as minimal effort strategy that complements sensorimotor feedback control in the presence of sensory noise. In simulations of reaching, we demonstrated that nonlinear multi-segment musculoskeletal dynamics enables complex perturbed and unperturbed reach trajectories under a variety of task conditions to be predicted. In both behaviors, we demonstrated how interactions between task constraint, sensory noise, and the intrinsic properties of muscle influence optimal muscle coordination patterns, including muscle co-contraction, and the resulting movement trajectories. Our approach enables a true minimum effort solution to be identified as task constraints, such as movement accuracy, can be explicitly imposed, rather than being approximated using penalty terms in the cost function. Our approximate stochastic optimal control framework predicts complex features, not captured by previous simulation approaches, providing a generalizable and valuable tool to study how musculoskeletal dynamics and physiological noise may alter neural control of movement in both healthy and pathological movements.

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.

Recent evidence suggests that complex spatiotemporal patterns of muscle activity can be explained with a low-dimensional set of muscle synergies or M-modes. While it is clear that both spatial and temporal aspects of muscle coordination may be low dimensional, constraints on spatial versus temporal features of muscle coordination likely involve different neural control mechanisms. We hypothesized that the low-dimensional spatial and temporal features of muscle coordination are independent of each other. We further hypothesized that in reactive feedback tasks, spatially fixed muscle coordination patterns—or muscle synergies—are hierarchically recruited via time-varying neural commands based on delayed task-level feedback. We explicitly compared the ability of spatially fixed (SF) versus temporally fixed (TF) muscle synergies to reconstruct the entire time course of muscle activity during postural responses to anterior-posterior support-surface translations. While both SF and TF muscle synergies could account for EMG variability in a postural task, SF muscle synergies produced more consistent and physiologically interpretable results than TF muscle synergies during postural responses to perturbations. Moreover, a majority of SF muscle synergies were consistent in structure when extracted from epochs throughout postural responses. Temporal patterns of SF muscle synergy recruitment were well-reconstructed by delayed feedback of center of mass (CoM) kinematics and reproduced EMG activity of multiple muscles. Consistent with the idea that independent and hierarchical low-dimensional neural control structures define spatial and temporal patterns of muscle activity, our results suggest that CoM kinematics are a task variable used to recruit SF muscle synergies for feedback control of balance.

Postural stability depends on interactions between the musculoskeletal system and neural control mechanisms. We present a frontal plane model stabilized by delayed feedback to analyze the effects of altered stance width on postural responses to perturbations. We hypothesized that changing stance width alters the mechanical dynamics of the body and limits the range of delayed feedback gains that produce stable postural behaviors. Surprisingly, mechanical stability was found to decrease as stance width increased due to decreased effective inertia. Furthermore, due to sensorimotor delays and increased leverage of hip joint torque on center-of-mass motion, the magnitudes of the stabilizing delayed feedback gains decreased as stance width increased. Moreover, the ranges of the stable feedback gains were nonoverlapping across different stance widths such that using a single neural feedback control strategy at both narrow and wide stances could lead to instability. The set of stable feedback gains was further reduced by constraints on foot lift-off and perturbation magnitude. Simulations were fit to experimentally measured kinematics, and the identified feedback gains corroborated model predictions. In addition, analytical gain margin of the linearized system was found to predict step transitions without the need for simulation. In conclusion, this model offers a method to dissociate the complex interactions between postural configuration, delayed sensorimotor feedback, and nonlinear foot lift-off constraints. The model demonstrates that stability at wide stances can only be achieved if delayed neural feedback gains decrease. This model may be useful in explaining both expected and paradoxical changes in stance width in healthy and neurologically impaired individuals.

Little is known about the integration of neural mechanisms for balance and locomotion. Muscle synergies have been studied independently in standing balance and walking, but not compared. Here, we hypothesized that reactive balance and walking are mediated by a common set of lower-limb muscle synergies. In humans, we examined muscle activity during multidirectional support-surface perturbations during standing and walking, as well as unperturbed walking at two speeds. We show that most muscle synergies used in perturbations responses during standing were also used in perturbation responses during walking, suggesting common neural mechanisms for reactive balance across different contexts. We also show that most muscle synergies using in reactive balance were also used during unperturbed walking, suggesting that neural circuits mediating locomotion and reactive balance recruit a common set of muscle synergies to achieve task-level goals. Differences in muscle synergies across conditions reflected differences in the biomechanical demands of the tasks. For example, muscle synergies specific to walking perturbations may reflect biomechanical challenges associated with single limb stance, and muscle synergies used during sagittal balance recovery in standing but not walking were consistent with maintaining the different desired center of mass motions in standing vs. walking. Thus, muscle synergies specifying spatial organization of muscle activation patterns may define a repertoire of biomechanical subtasks available to different neural circuits governing walking and reactive balance and may be recruited based on task-level goals. Muscle synergy analysis may aid in dissociating deficits in spatial vs. temporal organization of muscle activity in motor deficits. Muscle synergy analysis may also provide a more generalizable assessment of motor function by identifying whether common modular mechanisms are impaired across the performance of multiple motor tasks.

Optimality principles have been proposed as a general framework for understanding motor control in animals and humans largely based on their ability to predict general features movement in idealized motor tasks. However, generalizing these concepts past proof-of-principle to understand the neuromechanical transformation from task-level control to detailed execution-level muscle activity and forces during behaviorally-relevant motor tasks has proved difficult. In an unrestrained balance task in cats, we demonstrate that achieving task-level constraints center of mass forces and moments while minimizing control effort predicts detailed patterns of muscle activity and ground reaction forces in an anatomically-realistic musculoskeletal model. Whereas optimization is typically used to resolve redundancy at a single level of the motor hierarchy, we simultaneously resolved redundancy across both muscles and limbs and directly compared predictions to experimental measures across multiple perturbation directions that elicit different intra- and interlimb coordination patterns. Further, although some candidate task-level variables and cost functions generated indistinguishable predictions in a single biomechanical context, we identified a common optimization framework that could predict up to 48 experimental conditions per animal (n = 3) across both perturbation directions and different biomechanical contexts created by altering animals' postural configuration. Predictions were further improved by imposing experimentally-derived muscle synergy constraints, suggesting additional task variables or costs that may be relevant to the neural control of balance. These results suggested that reduced-dimension neural control mechanisms such as muscle synergies can achieve similar kinetics to the optimal solution, but with increased control effort (≈2×) compared to individual muscle control. Our results are consistent with the idea that hierarchical, task-level neural control mechanisms previously associated with voluntary tasks may also be used in automatic brainstem-mediated pathways for balance.

The mechanisms underlying associations between cognitive set shifting impairments and balance dysfunction are unclear. Cognitive set shifting refers to the ability to flexibly adjust behavior to changes in task rules or contexts, which could be involved in flexibly adjusting balance recovery behavior to different contexts, such as the direction the body is falling. Prior studies found associations between cognitive set shifting impairments and severe balance dysfunction in populations experiencing frequent falls. The objective of this study was to test whether cognitive set shifting ability is expressed in successful balance recovery behavior in older adults with high clinical balance ability (N = 19, 71 ± 7 years, 6 female). We measured cognitive set shifting ability using the Trail Making Test and clinical balance ability using the miniBESTest. For most participants, cognitive set shifting performance (Trail Making Test B-A = 37 ± 20 s) was faster than normative averages (46 s for comparable age and education levels), and balance ability scores (miniBESTest = 25 ± 2/28) were above the threshold for fall risk (23 for people between 70 and 80 years). Reactive balance recovery in response to support-surface translations in anterior and posterior directions was assessed in terms of body motion, muscle activity, and brain activity. Across participants, lower cognitive set shifting ability was associated with smaller peak center of mass displacement during balance recovery, lower directional specificity of late phase balance-correcting muscle activity (i.e., greater antagonist muscle activity 200–300 ms after perturbation onset), and larger cortical N1 responses (100–200 ms). None of these measures were associated with clinical balance ability. Our results suggest that cognitive set shifting ability is expressed in balance recovery behavior even in the absence of profound clinical balance disability. Specifically, our results suggest that lower flexibility in cognitive task performance is associated with lower ability to incorporate the directional context into the cortically mediated later phase of the motor response. The resulting antagonist activity and stiffer balance behavior may help explain associations between cognitive set shifting impairments and frequent falls.