Approximating complex musculoskeletal biomechanics using multidimensional autogenerating polynomials.

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 evidence for nonlinear feedforward modulation of fusimotor drive to antagonistic co-contracting muscles.

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.

Correction: Biomechanical Constraints Underlying Motor Primitives Derived from the Musculoskeletal Anatomy of the Human Arm.

[This corrects the article DOI: 10.1371/journal.pone.0164050.].

Biomechanical Constraints Underlying Motor Primitives Derived from the Musculoskeletal Anatomy of the Human Arm.

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.