Do Muscle Synergies Improve Optimization Prediction of Muscle Activations During Gait?

Positional Analysis of Assisting Muscles for Handling-Assisted Exoskeletons

In order to better design handling-assisted exoskeletons, it is necessary to analyze the biomechanics of human hand movements. In this study, Anybody Modeling System (AMS) simulation was used to analyze the movement state of muscles during human handling. Combined with surface electromyography (sEMG) experiments, specific analysis and verification were carried out to obtain the position of muscles that the human body needs to assist during handling. In this study, the simulation and experiment were carried out for the manual handling process. A treatment group and an experimental group were set up. This study found that the vastus medialis muscle, vastus lateralis muscle, latissimus dorsi muscle, trapezius muscle, deltoid muscle and triceps brachii muscle require more energy in the process of handling, and it is reasonable and effective to combine sEMG signals with the simulation of the musculoskeletal model to analyze the muscle condition of human movement.

Does crouch alter the effects of neuromuscular impairments on gait? A simulation study

Cerebral palsy (CP) is a neurologic injury that impacts control of movement. Individuals with CP also often develop secondary impairments like weakness and contracture. Both altered motor control and secondary impairments influence how an individual walks after neurologic injury. However, understanding the complex interactions between and relative effects of these impairments makes analyzing and improving walking capacity in CP challenging. We used a sagittal-plane musculoskeletal model and neuromuscular control framework to simulate crouch and nondisabled gait. We perturbed each simulation by varying the number of synergies controlling each leg (altered control), and imposed weakness and contracture. A Bayesian Additive Regression Trees (BART) model was also used to parse the relative effects of each impairment on the muscle activations required for each gait pattern. By using these simulations to evaluate gait-pattern specific effects of neuromuscular impairments, we identified some advantages of crouch gait. For example, crouch tolerated 13 % and 22 % more plantarflexor weakness than nondisabled gait without and with altered control, respectively. Furthermore, BART demonstrated that plantarflexor weakness had twice the effect on total muscle activity required during nondisabled gait than crouch gait. However, crouch gait was also disadvantageous in the presence of vasti weakness: crouch gait increased the effects of vasti weakness on gait without and with altered control. These simulations highlight gait-pattern specific effects and interactions between neuromuscular impairments. Utilizing computational techniques to understand these effects can elicit advantages of gait deviations, providing insight into why individuals may select their gait pattern and possible interventions to improve energetics.

Reinforcement-based processes actively regulate motor exploration along redundant solution manifolds

From a baby's babbling to a songbird practising a new tune, exploration is critical to motor learning. A hallmark of exploration is the emergence of random walk behaviour along solution manifolds, where successive motor actions are not independent but rather become serially dependent. Such exploratory random walk behaviour is ubiquitous across species' neural firing, gait patterns and reaching behaviour. The past work has suggested that exploratory random walk behaviour arises from an accumulation of movement variability and a lack of error-based corrections. Here, we test a fundamentally different idea-that reinforcement-based processes regulate random walk behaviour to promote continual motor exploration to maximize success. Across three human reaching experiments, we manipulated the size of both the visually displayed target and an unseen reward zone, as well as the probability of reinforcement feedback. Our empirical and modelling results parsimoniously support the notion that exploratory random walk behaviour emerges by utilizing knowledge of movement variability to update intended reach aim towards recently reinforced motor actions. This mechanism leads to active and continuous exploration of the solution manifold, currently thought by prominent theories to arise passively. The ability to continually explore muscle, joint and task redundant solution manifolds is beneficial while acting in uncertain environments, during motor development or when recovering from a neurological disorder to discover and learn new motor actions.

A procedure and model for the identification of uni- and biarticular structures passive contribution to inter-segmental dynamics

Inter-segmental moments come from muscles contractions, but also from passive moments, resulting from the resistance of the periarticular structures. To quantify the passive contribution of uni- and biarticular structures during gait, we propose an innovative procedure and model. 12 typically developed (TD) children and 17 with cerebral palsy (CP) participated in a passive testing protocol. The relaxed lower limb joints were manipulated through full ranges of motion while kinematics and applied forces were simultaneously measured. The relationships between uni-/biarticular passive moments/forces and joint angles/musculo-tendon lengths were modelled by a set of exponential functions. Then, subject specific gait joint angles/musculo-tendon lengths were input into the determined passive models to estimate joint moments and power attributable to passive structures. We found that passive mechanisms contribute substantially in both populations, mainly during push-off and swing phases for hip and knee and push-off for the ankle, with a distinction between uni- and biarticular structures. CP children showed comparable passive mechanisms but larger variability than the TD ones and higher contributions. The proposed procedure and model enable a comprehensive assessment of the passive mechanisms for a subject-specific treatment of the stiffness implying gait disorders by targeting when and how passive forces are impacting gait.

EMG-driven musculoskeletal model calibration with estimation of unmeasured muscle excitations via synergy extrapolation

Subject-specific electromyography (EMG)-driven musculoskeletal models that predict muscle forces have the potential to enhance our knowledge of internal biomechanics and neural control of normal and pathological movements. However, technical gaps in experimental EMG measurement, such as inaccessibility of deep muscles using surface electrodes or an insufficient number of EMG channels, can cause difficulties in collecting EMG data from muscles that contribute substantially to joint moments, thereby hindering the ability of EMG-driven models to predict muscle forces and joint moments reliably. This study presents a novel computational approach to address the problem of a small number of missing EMG signals during EMG-driven model calibration. The approach (henceforth called "synergy extrapolation" or SynX) linearly combines time-varying synergy excitations extracted from measured muscle excitations to estimate 1) unmeasured muscle excitations and 2) residual muscle excitations added to measured muscle excitations. Time-invariant synergy vector weights defining the contribution of each measured synergy excitation to all unmeasured and residual muscle excitations were calibrated simultaneously with EMG-driven model parameters through a multi-objective optimization. The cost function was formulated as a trade-off between minimizing joint moment tracking errors and minimizing unmeasured and residual muscle activation magnitudes. We developed and evaluated the approach by treating a measured fine wire EMG signal (iliopsoas) as though it were "unmeasured" for walking datasets collected from two individuals post-stroke-one high functioning and one low functioning. How well unmeasured muscle excitations and activations could be predicted with SynX was assessed quantitatively for different combinations of SynX methodological choices, including the number of synergies and categories of variability in unmeasured and residual synergy vector weights across trials. The two best methodological combinations were identified, one for analyzing experimental walking trials used for calibration and another for analyzing experimental walking trials not used for calibration or for predicting new walking motions computationally. Both methodological combinations consistently provided reliable and efficient estimates of unmeasured muscle excitations and activations, muscle forces, and joint moments across both subjects. This approach broadens the possibilities for EMG-driven calibration of muscle-tendon properties in personalized neuromusculoskeletal models and may eventually contribute to the design of personalized treatments for mobility impairments.

InverseMuscleNET: Alternative Machine Learning Solution to Static Optimization and Inverse Muscle Modeling

InverseMuscleNET, a machine learning model, is proposed as an alternative to static optimization for resolving the redundancy issue in inverse muscle models. A recurrent neural network (RNN) was optimally configured, trained, and tested to estimate the pattern of muscle activation signals. Five biomechanical variables (joint angle, joint velocity, joint acceleration, joint torque, and activation torque) were used as inputs to the RNN. A set of surface electromyography (EMG) signals, experimentally measured around the shoulder joint for flexion/extension, were used to train and validate the RNN model. The obtained machine learning model yields a normalized regression in the range of 88-91% between experimental data and estimated muscle activation. A sequential backward selection algorithm was used as a sensitivity analysis to discover the less dominant inputs. The order of most essential signals to least dominant ones was as follows: joint angle, activation torque, joint torque, joint velocity, and joint acceleration. The RNN model required 0.06 s of the previous biomechanical input signals and 0.01 s of the predicted feedback EMG signals, demonstrating the dynamic temporal relationships of the muscle activation profiles. The proposed approach permits a fast and direct estimation ability instead of iterative solutions for the inverse muscle model. It raises the possibility of integrating such a model in a real-time device for functional rehabilitation and sports evaluation devices with real-time estimation and tracking. This method provides clinicians with a means of estimating EMG activity without an invasive electrode setup.

A musculoskeletal model driven by muscle synergy-derived excitations for hand and wrist movements

OBJECTIVE

Musculoskeletal model (MM) driven by electromyography (EMG) signals has been identified as a promising approach to predicting human motions in the control of prostheses and robots. However, muscle excitations in MMs are generally derived from the EMG signals of the targeted sensor covering the muscle, inconsistent with the fact that signals of a sensor are from multiple muscles considering signal crosstalk in actual situation. To identify more accurate muscle excitations for MM in the presence of crosstalk, we proposed a novel excitation-extracting method inspired by muscle synergy for simultaneously estimating hand and wrist movements.

APPROACH

Muscle excitations were firstly extracted using a two-step muscle synergy-derived method. Specifically, we calculated subject-specific muscle weighting matrix and corresponding profiles according to contributions of different muscles for movements derived from synergistic motion relation. Then, the improved excitations were used to simultaneously estimate hand and wrist movements through musculoskeletal modeling. Moreover, the offline comparison among the proposed method, traditional MM and regression methods, and an online test of the proposed method were conducted.

MAIN RESULTS

The offline experiments demonstrated that the proposed approach outperformed the EMG envelope-driven MM and three regression models with higher R and lower NRMSE. Furthermore, the comparison of excitations of two MMs validated the effectiveness of the proposed approach in extracting muscle excitations in the presence of crosstalk. The online test further indicated the superior performance of the proposed method than the MM driven by EMG envelopes.

SIGNIFICANCE

The proposed excitation-extracting method identified more accurate neural commands for MMs, providing a promising approach in rehabilitation and robot control to model the transformation from surface EMG to joint kinematics.

A fair and EMG-validated comparison of recruitment criteria, musculotendon models and muscle coordination strategies, for the inverse-dynamics based optimization of muscle forces during gait

Experimental studies and EMG collections suggest that a specific strategy of muscle coordination is chosen by the central nervous system to perform a given motor task. A popular mathematical approach for solving the muscle recruitment problem is optimization. Optimization-based methods minimize or maximize some criterion (objective function or cost function) which reflects the mechanism used by the central nervous system to recruit muscles for the movement considered. The proper cost function is not known a priori, so the adequacy of the chosen function must be validated according to the obtained results. In addition of the many criteria proposed, several physiological representations of the musculotendon actuator dynamics (that prescribe constraints for the forces) along with different musculoskeletal models can be found in the literature, which hinders the selection of the best neuromusculotendon model for each application. Seeking to provide a fair base for comparison, this study measures the efficiency and accuracy of: (i) four different criteria within the static optimization approach (where the physiological character of the muscle, which affects the constraints of the forces, is not considered); (ii) three physiological representations of the musculotendon actuator dynamics: activation dynamics with elastic tendon, simplified activation dynamics with rigid tendon and rigid tendon without activation dynamics; (iii) a synergy-based method; all of them within the framework of inverse-dynamics based optimization. Motion/force/EMG gait analyses were performed on ten healthy subjects. A musculoskeletal model of the right leg actuated by 43 Hill-type muscles was scaled to each subject and used to calculate joint moments, musculotendon kinematics and moment arms. Muscle activations were then estimated using the different approaches, and these estimates were compared with EMG measurements. Although no significant differences were obtained with all the methods at statistical level, it must be pointed out that a higher complexity of the method does not guarantee better results, as the best correlations with experimental values were obtained with two simplified approaches: the static optimization and the physiological approach with simplified activation dynamics and rigid tendon, both using the sum of the squares of muscle forces as objective function.