Neural control of movement can only be realized though the interaction between the mechanical properties of the limb and the environment. Thus, a fundamental question is whether anatomy has evolved to simplify neural control by shaping these interactions in a beneficial way. This inductive data-driven study analyzed the patterns of muscle actions across multiple joints using the musculoskeletal model of the human upper limb. This model was used to calculate muscle lengths across the full range of motion of the arm and examined the correlations between these values between all pairs of muscles. Musculoskeletal coupling was quantified using hierarchical clustering analysis. Muscle lengths between multiple pairs of muscles across multiple postures were highly correlated. These correlations broadly formed two proximal and distal groups, where proximal muscles of the arm were correlated with each other and distal muscles of the arm and hand were correlated with each other, but not between groups. Using hierarchical clustering, between 11 and 14 reliable muscle groups were identified. This shows that musculoskeletal anatomy does indeed shape the mechanical interactions by grouping muscles into functional clusters that generally match the functional repertoire of the human arm. Together, these results support the idea that the structure of the musculoskeletal system is tuned to solve movement complexity problem by reducing the dimensionality of available solutions.

the sensorimotor integration during unconstrained reaching movements in the presence of variable environmental forces remains poorly understood. the objective of this study was to quantify how much the primary afferent activity of muscle spindles can contribute to shaping muscle coactivation patterns during reaching movements with complex dynamics. to achieve this objective, we designed a virtual reality task that guided healthy human participants through a set of planar reaching movements with controlled kinematic and dynamic conditions that were accompanied by variable muscle co-contraction. Next, we approximated the Ia afferent activity using a phenomenological model of the muscle spindle and muscle lengths derived from a musculoskeletal model. the parameters of the spindle model were altered systematically to evaluate the effect of fusimotor drive on the shape of the temporal profile of afferent activity during movement. The experimental and simulated data were analyzed with hierarchical clustering. We found that the pattern of co-activation of agonistic and antagonistic muscles changed based on whether passive forces in each movement played assistive or resistive roles in limb dynamics. the reaching task with assistive limb dynamics was associated with the most muscle co-contraction. In contrast, the simulated Ia afferent profiles were not changing between tasks and they were largely reciprocal with homonymous muscle activity. Simulated physiological changes to the fusimotor drive were not sufficient to reproduce muscle co-contraction. These results largely rule out the static set and α-γ coactivation as the main types of fusimotor drive that transform the monosynaptic Ia afferent feedback into task-dependent co-contraction of antagonistic muscles. We speculate that another type of nonlinear transformation of Ia afferent signals that is independent of signals modulating the activity of α motoneurons is required for Ia afferent-based co-contraction. This transformation could either be applied through a complex nonlinear profile of fusimotor drive that is not yet experimentally observed or through presynaptic inhibition. Movement is the product of interactions between neural signals and the musculoskeletal dynamics that depends on limb anatomy 1-4. The motor control problem is then solved within a system with coupled neural and mechanical dynamical elements 5-7. Therefore, the relationship between neural signals driving muscle contraction and the resulting motion is nonlinear. Muscle contractions generate forces that sum into active moments defined by the agonistic or antagonistic relationships between the muscle's moment arms around a given axis of rotation of the joint. The components of these forces that sum to zero moment, such as forces produced by balanced

Computational models of the musculoskeletal system are scientific tools used to study human movement, quantify the effects of injury and disease, plan surgical interventions, or control realistic high-dimensional articulated prosthetic limbs. If the models are sufficiently accurate, they may embed complex relationships within the sensorimotor system. These potential benefits are limited by the challenge of implementing fast and accurate musculoskeletal computations. A typical hand muscle spans over 3 degrees of freedom (DOF), wrapping over complex geometrical constraints that change its moment arms and lead to complex posture-dependent variation in torque generation. Here, we report a method to accurately and efficiently calculate musculotendon length and moment arms across all physiological postures of the forearm muscles that actuate the hand and wrist. Then, we use this model to test the hypothesis that the functional similarities of muscle actions are embedded in muscle structure. The posture dependent muscle geometry, moment arms and lengths of modeled muscles were captured using autogenerating polynomials that expanded their optimal selection of terms using information measurements. The iterative process approximated 33 musculotendon actuators, each spanning up to 6 DOFs in an 18 DOF model of the human arm and hand, defined over the full physiological range of motion. Using these polynomials, the entire forearm anatomy could be computed in <10 μs, which is far better than what is required for real-time performance, and with low errors in moment arms (below 5%) and lengths (below 0.4%). Moreover, we demonstrate that the number of elements in these autogenerating polynomials does not increase exponentially with increasing muscle complexity; complexity increases linearly instead. Dimensionality reduction using the polynomial terms alone resulted in clusters comprised of muscles with similar functions, indicating the high accuracy of approximating models. We propose that this novel method of describing musculoskeletal biomechanics might further improve the applications of detailed and scalable models to describe human movement.

Computational models of the musculoskeletal system are scientific tools used to study human movement, quantify the effects of injury and disease, and plan surgical interventions. Additionally, these models could also be used to intuitively link biological control signals and realistic high-dimensional articulated prosthetic limbs. However, implementing fast and accurate musculoskeletal computations that can be used to control a prosthetic limb in real-time is a challenging problem. As muscles typically span multiple joints, the wrapping over complex geometrical constraints changes their moment arms and length as a function of joint angle and, thus, their ability to generate joint torques. As a result of these biomechanical complexities, calculating these muscle state variables in real-time is a difficult simulation problem. Here, we report a method to accurately and efficiently calculate these variables for the forearm muscles that actuate the hand and wrist across multiple postures. The posture dependent muscle geometry, moment arms and lengths of modeled muscles, were captured using autogenerating polynomials that expanded their optimal selection of terms using information measurements. The iterative process approximated 33 musculotendon actuators, each spanning up to 6 DOFs in an 18 DOF model of the human arm and hand, defined over the full physiological range of motion. Using these polynomials, the entire forearm anatomy could be computed in <10 µs, which is far better than what is required for real-time performance, and with low errors in moment arms (below 5%) and lengths (below 0.4%). Moreover, we demonstrate that the number of elements in these autogenerating polynomials does not increase exponentially with the increase in complexity of muscles, increasing linearly instead. The similar structure and function of muscles are represented with specific invariant polynomial terms. Dimensionality reduction using the polynomial terms alone resulted in clusters comprised of muscles with similar functions, suggesting that the polynomials themselves captured biologically relevant features of muscle structure and function. We propose that this novel method of describing musculoskeletal biomechanics might further improve the applications of detailed and scalable models for the description of human movement.