Predicting gait adaptations due to ankle plantarflexor muscle weakness and contracture using physics-based musculoskeletal simulations

The role of the plantarflexor muscles on the daily life functioning in typically developing children

BACKGROUND

The plantarflexor muscles are important during walking. However, walking is not the only task that is used during daily life, and not much is known about the role of the plantarflexor muscles during other dynamic tasks. Insight into the demand of the plantarflexor muscles during different dynamic tasks might help to establish training or functional assessment for patients who experience plantarflexor weakness.

RESEARCH QUESTION

This study aimed to examine the role of the plantarflexor functioning during dynamic tasks in children of 5-9 years old, by identifying differences in ankle kinetics and plantarflexor muscle forces during dynamic tasks.

METHODS

In this observational cross-sectional study, 28 typically developing children from 5 to 9 years old were included. Three-dimensional movement analysis was performed to extract the peak ankle moment, power and plantarflexor muscle forces during the stance phase of walking, heel walking, toe walking, running and one-leg-hopping. To examine differences in peak ankle moment, ankle power and muscle forces, Friedman tests were used in combination with Wilcoxon signed-rank tests to determine differences in peak ankle moment, power and muscle forces of the various dynamic tasks compared to walking.

RESULTS

Based on the differences in peak ankle moment, power and muscle forces compared to walking, heel walking was the least demanding dynamic task. Running and one-leg-hopping were the most demanding tasks.

SIGNIFICANCE

Dynamic tasks with a high demand on the plantarflexor muscles, such as running and one-leg-hopping could be used to quantify deficits in patient populations with plantarflexor weakness that are not visible during walking. Furthermore, these tasks could be used during functional evaluation or training in these patient populations.

Active Arm Swing During Running Improves Rotational Stability of the Upper Body and Metabolic Energy Efficiency

PURPOSE

The kinematic benefits of arm swing during running for upper body stability have been previously investigated, while its role in metabolic energy efficiency remains controversial. To address this, this study aimed to test the hypothesis that active arm swing during running reduces both torso angular motion around the longitudinal axis and metabolic energy consumption.

METHODS

We employed forward dynamics musculoskeletal running simulations with different arm conditions to investigate the hypothesis. Full-body musculoskeletal running models, incorporating 150 muscles, were developed using artificial neural network-based running controllers. Three arm conditions were simulated using the running models and controllers: active arm swing, passive arm swing, and fixed arms.

RESULTS

Our results revealed that the active arm swing model demonstrated the lowest total metabolic energy consumption per traveling distance. The costs of transport were 5.52, 5.73, and 5.82 J/kg-m for active, passive, and fixed arm models, respectively. Interestingly, while metabolic energy consumption in the upper limb muscles was higher during active arm swing, the total energy consumption was lower. Additionally, the longitudinal rotation of the torso was minimal in the active arm swing condition.

CONCLUSION

These findings support our hypothesis, demonstrating that active arm swing during running reduces the angular motion of the torso and the metabolic energy consumption. This study provides evidence that arm swing during running is performed actively as an energy-saving mechanism. These results contribute to understanding of running biomechanics and may have implications for performance optimization in sports and rehabilitation settings.

How knee muscles and ground reaction forces shape knee buckling and ankle push-off in neuromuscular simulations of human walking

Ankle push-off is important for efficient, human-like walking, and many prosthetic devices mimic push-off using motors or elastic elements. The knee is extended throughout the stance phase and begins to buckle just before push-off, with timing being crucial. However, the exact mechanisms behind this buckling are still unclear. We use a predictive neuromuscular simulation to investigate whether active muscles are required for knee buckling and to what extent ground reaction forces (GRFs) drive it. In a systematic parameter search, we tested how long the knee muscles vastus (VAS), gastrocnemius (GAS), and hamstrings could be deactivated while maintaining a stable gait with impulsive push-off. VAS deactivation up to 35% of the gait cycle resulted in a dynamic gait with increased ankle peak power. GAS deactivation up to 20% of the gait cycle was detrimental to gait efficiency and showed reduced ankle peak power. At the start of knee buckling, the GRF vector is positioned near the knee joint's neutral axis, assisting in knee flexion. However, this mechanism is likely not enough to drive knee flexion independently. Our findings contribute to the biomechanical understanding of ankle push-off, with applications in prosthetic and bipedal robotic design, and fundamental research on human gait mechanics.

Analyzing gait data measured by wearable cyborg hybrid assistive limb during assisted walking: gait pattern clustering

Introduction

The wearable cyborg Hybrid Assistive Limb (HAL) is a therapeutic exoskeletal device that provides voluntary gait assistance using kinematic/kinetic gait data and bioelectrical signals. By utilizing the gait data automatically measured by HAL, we are developing a system to analyze the wearer's gait during the intervention, unlike conventional evaluations that compare pre- and post-treatment gait test results. Despite the potential use of the gait data from the HAL's sensor information, there is still a lack of analysis using such gait data and knowledge of gait patterns during HAL use. This study aimed to cluster gait patterns into subgroups based on the gait data that the HAL automatically collected during treatment and to investigate their characteristics.

Methods

Gait data acquired by HAL, including ground reaction forces, joint angles, trunk angles, and HAL joint torques, were analyzed in individuals with progressive neuromuscular diseases. For each measured item, principal component analysis was applied to the gait time-series data to extract the features of the gait patterns, followed by hierarchical cluster analysis to generate subgroups based on the principal component scores. Bayesian regression analysis was conducted to identify the influence of the wearer's attributes on the clustered gait patterns.

Results

The gait patterns of 13,710 gait cycles from 457 treatments among 48 individuals were divided into 5-10 clusters for each measured item. The clusters revealed a variety of gait patterns when wearing the HAL and identified the characteristics of multiple sub-group types. Bayesian regression models explained the influence of the wearer's disease type and gait ability on the distribution of gait patterns to subgroups.

Discussion

These results revealed key differences in gait patterns related to the wearer's condition, demonstrating the importance of monitoring HAL-assisted walking to provide appropriate interventions. Furthermore, our approach highlights the usefulness of the gait data that HAL automatically measures during the intervention. We anticipate that the HAL, designed as a therapeutic device, will expand its role as a data measurement device for analysis and evaluation that provides gait data simultaneously with interventions, creating a novel cybernics treatment system that facilitates a multi-faceted understanding of the wearer's gait.

Novel neuromuscular controllers with simplified muscle model and enhanced reflex modulation: A comparative study in hip exoskeletons

Neuromuscular controllers (NMCs) offer a promising approach to adaptive and task-invariant control of exoskeletons for walking assistance, leveraging the bioinspired models based on the peripheral nervous system. This article expands on our previous development of a novel structure for NMCs with modifications to the virtual muscle model and reflex modulation strategy. The modifications consist firstly of simplifications to the Hill-type virtual muscle model, resulting in a more straightforward formulation and reduced number of parameters; and second, using a finer division of gait subphases in the reflex modulation state machine, allowing for a higher degree of control over the shape of the assistive profile. Based on the proposed general structure, we present two controller variants for hip exoskeletons, with four- and five-state reflex modulations (NMC-4 and NMC-5). We used an iterative data-driven approach with two tuning stages (i.e., muscle parameters and reflex gains) to determine the controller parameters. Biological joint torque profiles and optimal torque profiles for metabolic cost reduction were used as references for the final tuning outcome. Experimental testing under various walking conditions demonstrated the capability of both variants for adapting to the locomotion task with minimal parameter adjustments, mostly in terms of timing. Furthermore, NMC-5 exhibited better alignment with biological and optimised torque profiles in terms of timing characteristics and relative magnitudes, resulting in less negative mechanical work. These findings firstly validate the adequacy of the simplified muscle model for assistive controllers, and demonstrate the utility of a more nuanced reflex modulation in improving the assistance quality.

Kinematic patterns in performing trunk flexion tasks influenced by various mechanical optimization targets: A simulation study

BACKGROUND

Low back pain is the most prevalent and disabling condition worldwide, with a high recurrence rate in the general adult population.

METHODS

A set of open-sourced trunk musculoskeletal models was used to investigate trunk flexion kinematics under different motor control strategies, including minimizing shearing or compressive loads at the L4/L5 or L5/S1 level.

FINDINGS

A control strategy that minimizes the load on the lower lumbar intervertebral disc can result in two kinematic patterns-the "restricted lumbar spine" and the "overflexed lumbar spine"-in performing the trunk flexion task. The "restricted" pattern can reduce the overall load on the lower lumbar levels, whereas the "overflexed" pattern can reduce the shearing force only at the L4/L5 level and increase the compressive and shearing forces at the L5/S1 level and the compressive force at the L4/L5 level.

INTERPRETATION

This study investigated the relationships between specific trunk kinematics in patients with low back pain and lumbar intervertebral loading via musculoskeletal modelling and simulation. The results provide insight into individualized treatment for patients with low back pain.

The Neuromusculoskeletal Modeling Pipeline: MATLAB-based Model Personalization and Treatment Optimization Functionality for OpenSim

Neuromusculoskeletal injuries including osteoarthritis, stroke, spinal cord injury, and traumatic brain injury affect roughly 19% of the U.S. adult population. Standardized interventions have produced suboptimal functional outcomes due to the unique treatment needs of each patient. Strides have been made to utilize computational models to develop personalized treatments, but researchers and clinicians have yet to cross the "valley of death" between fundamental research and clinical usefulness. This article introduces the Neuromusculoskeletal Modeling (NMSM) Pipeline, two MATLAB-based toolsets that add Model Personalization and Treatment Optimization functionality to OpenSim. The two toolsets facilitate computational design of individualized treatments for neuromusculoskeletal impairments through the use of personalized neuromusculoskeletal models and predictive simulation. The Model Personalization toolset contains four tools for personalizing 1) joint structure models, 2) muscle-tendon models, 3) neural control models, and 4) foot-ground contact models. The Treatment Optimization toolset contains three tools for predicting and optimizing a patient's functional outcome for different treatment options using a patient's personalized neuromusculoskeletal model with direct collocation optimal control methods. Support for user-defined cost functions and model modification functions facilitate simulation of a vast number of possible treatments. An NMSM Pipeline use case is presented for an individual post-stroke with impaired walking function, where the goal was to predict how the subject's neural control could be changed to improve walking speed without increasing metabolic cost. First the Model Personalization toolset was used to develop a personalized neuromusculoskeletal model of the subject starting from a generic OpenSim full-body model and experimental walking data (video motion capture, ground reaction, and electromyography) collected from the subject at his self-selected speed. Next the Treatment Optimization toolset was used with the personalized model to predict how the subject could recruit existing muscle synergies more effectively to reduce muscle activation disparities between the paretic and non-paretic legs. The software predicted that the subject could increase his walking speed by 60% without increasing his metabolic cost per unit time by modifying existing muscle synergy recruitment. This hypothetical treatment demonstrates how NMSM Pipeline tools could allow researchers working collaboratively with clinicians to develop personalized neuromusculoskeletal models of individual patients and to perform predictive simulations for the purpose of designing personalized treatments that maximize a patient's post-treatment functional outcome.

Perturbation-based estimation of within-stride cycle metabolic cost

Metabolic cost greatly impacts trade-offs within a variety of human movements. Standard respiratory measurements only obtain the mean cost of a movement cycle, preventing understanding of the contributions of different phases in, for example, walking. We present a method that estimates the within-stride cost of walking by leveraging measurements under different force perturbations. The method reproduces time series with greater consistency (r = 0.55 and 0.80 in two datasets) than previous model-based estimations (r = 0.29). This perturbation-based method reveals how the cost of push-off (10%) is much smaller than would be expected from positive mechanical work (~ 70%). This work elucidates the costliest phases during walking, offering new targets for assistive devices and rehabilitation strategies.

Exploring the Role of Asymmetric Auditory and Tactile Stimulation on Modulating Gait Kinetics

This study explores the influence of Auditory Rhythmic Asymmetric Cueing (A-RAC), Tactile Rhythmic Asymmetric Cueing (T-RAC), and their combination (AT) on key kinetic gait parameters in gait rehabilitation: Vertical Ground Reaction Force Asymmetry (GRF), Push-off Force Asymmetry (POF), and Braking Force Asymmetry (BRK). Utilizing the Computer-Assisted Rehabilitation Environment (CAREN) with 18 participants, this research examines these interventions' effectiveness in generating asymmetric gait. While the results during adaptation indicate that BRK was significantly affected by both A-RAC (p = 0.001) and AT (p = 0.003), only A-RAC had a significant effect on GRF (p = 0.002) during adaptation. None of the interventions significantly altered POF, suggesting its resistance to sensory cue modification. These findings provide valuable insights for enhancing gait rehabilitation strategies, particularly in addressing vertical load distribution, controlled deceleration, and overall walking safety.

The effect of including a mobile arch, toe joint, and joint coupling on predictive neuromuscular simulations of human walking

Predictive neuromuscular simulations are a powerful tool for studying the biomechanics of human walking, and deriving design criteria for technical devices like prostheses or biorobots. Good agreement between simulation and human data is essential for transferability to the real world. The human foot is often modeled with a single rigid element, but knowledge of how the foot model affects gait prediction is limited. Standardized procedures for selecting appropriate foot models are lacking. We performed 2D predictive neuromuscular simulations with six different foot models of increasing complexity to answer two questions: What is the effect of a mobile arch, a toe joint, and the coupling of toe and arch motion through the plantar fascia on gait prediction? and How much of the foot's anatomy do we need to model to predict sagittal plane walking kinematics and kinetics in good agreement with human data? We found that the foot model had a significant impact on ankle kinematics during terminal stance, push-off, and toe and arch kinematics. When focusing only on hip and knee kinematics, rigid foot models are sufficient. We hope our findings will help guide the community in modeling the human foot according to specific research goals and improve neuromuscular simulation accuracy.

A planar neuromuscular controller to simulate compensation strategies in the sit-to-walk movement

Standing up from a chair is a key daily life activity that is sensitive to functional limitations as we age and associated with falls, frailty, and institutional living. Predictive neuromusculoskeletal models can potentially shed light on the interconnectivity and interdependency of age-related changes in neuromuscular capacity, reinforcement schemes, sensory integration, and adaptation strategies during stand-up. Most stand-up movements transfer directly into walking (sit-to-walk). The aim of this study was to develop and validate a neuromusculoskeletal model with reflex-based muscle control that enables simulation of the sit-to-walk movement under various conditions (seat height, foot placement). We developed a planar sit-to-walk musculoskeletal model (11 degrees-of-freedom, 20 muscles) and neuromuscular controller, consisting of a two-phase stand-up controller and a reflex-based gait controller. The stand-up controller contains generic neural pathways of delayed proprioceptive feedback from muscle length, force, velocity, and upper-body orientation (vestibular feedback) and includes both monosynaptic an antagonistic feedback pathways. The control parameters where optimized using a shooting-based optimization method, based on a high-level optimization criterium. Simulations were compared to recorded kinematics, ground reaction forces, and muscle activation. The simulated kinematics resemble the measured kinematics and muscle activations. The adaptation strategies that resulted from alterations in seat height, are comparable to those observed in adults. The simulation framework and model are publicly available and allow to study age-related compensation strategies, including reduced muscular capacity, reduced neural capacity, external perturbations, and altered movement objectives.

Rethinking Exoskeleton Simulation-Based Design: The Effect of Using Different Cost Functions

Designing an exoskeleton that can improve user capabilities is a challenging task, and most designs rely on experiments to achieve this goal. A different approach is to use simulation-based designs to determine optimal device parameters. Most of these simulations use full trajectory tracking limb kinematics during a natural gait as a reference. However, exoskeletons typically change the natural gait kinematics of the user. Other types of simulations assume that human gait is optimized for a cost function that combines several objectives, such as the cost of transport, injury prevention, and stabilization. In this study, we use a 2D OpenSim model consisting of 10 degrees of freedom and considering 18 muscles, together with the Moco optimization tool, to investigate the differences between these two approaches with respect to running with a passive knee exoskeleton. Utilizing this model, we test the effect of a full trajectory tracking objective with different weights (representing the importance of the objective in the optimization cost function) and show that when using weights that are typically used in the literature, there is no deviation from the experimental data. Next, we develop a multi-objective cost function with foot clearance term based on peak knee angle during swing, that achieves trajectories similar (RMSE=7.4 deg) to experimental running data. Finally, we investigate the effect of different parameters in the design of a clutch-based passive knee exoskeleton (1.5 kg at each leg) and find that a design that utilizes a 2.5 Nm/deg spring achieves an improvement of up to 8% in net metabolic energy. Our results show that tracking objectives in the cost function, even with a low weight, hinders the simulation's ability to change the gait trajectory. Thus, there is a need for other predictive simulation methods for exoskeletons.

A dynamic foot model for predictive simulations of human gait reveals causal relations between foot structure and whole-body mechanics

The unique structure of the human foot is seen as a crucial adaptation for bipedalism. The foot's arched shape enables stiffening the foot to withstand high loads when pushing off, without compromising foot flexibility. Experimental studies demonstrated that manipulating foot stiffness has considerable effects on gait. In clinical practice, altered foot structure is associated with pathological gait. Yet, experimentally manipulating individual foot properties (e.g. arch height or tendon and ligament stiffness) is hard and therefore our understanding of how foot structure influences gait mechanics is still limited. Predictive simulations are a powerful tool to explore causal relationships between musculoskeletal properties and whole-body gait. However, musculoskeletal models used in three-dimensional predictive simulations assume a rigid foot arch, limiting their use for studying how foot structure influences three-dimensional gait mechanics. Here, we developed a four-segment foot model with a longitudinal arch for use in predictive simulations. We identified three properties of the ankle-foot complex that are important to capture ankle and knee kinematics, soleus activation, and ankle power of healthy adults: (1) compliant Achilles tendon, (2) stiff heel pad, (3) the ability to stiffen the foot. The latter requires sufficient arch height and contributions of plantar fascia, and intrinsic and extrinsic foot muscles. A reduced ability to stiffen the foot results in walking patterns with reduced push-off power. Simulations based on our model also captured the effects of walking with anaesthetised intrinsic foot muscles or an insole limiting arch compression. The ability to reproduce these different experiments indicates that our foot model captures the main mechanical properties of the foot. The presented four-segment foot model is a potentially powerful tool to study the relationship between foot properties and gait mechanics and energetics in health and disease.

Minimization of metabolic cost of transport predicts changes in gait mechanics over a range of ankle-foot orthosis stiffnesses in individuals with bilateral plantar flexor weakness

Neuromuscular disorders often lead to ankle plantar flexor muscle weakness, which impairs ankle push-off power and forward propulsion during gait. To improve walking speed and reduce metabolic cost of transport (mCoT), patients with plantar flexor weakness are provided dorsal-leaf spring ankle-foot orthoses (AFOs). It is widely believed that mCoT during gait depends on the AFO stiffness and an optimal AFO stiffness that minimizes mCoT exists. The biomechanics behind why and how an optimal stiffness exists and benefits individuals with plantar flexor weakness are not well understood. We hypothesized that the AFO would reduce the required support moment and, hence, metabolic cost contributions of the ankle plantar flexor and knee extensor muscles during stance, and reduce hip flexor metabolic cost to initiate swing. To test these hypotheses, we generated neuromusculoskeletal simulations to represent gait of an individual with bilateral plantar flexor weakness wearing an AFO with varying stiffness. Predictions were based on the objective of minimizing mCoT, loading rates at impact and head accelerations at each stiffness level, and the motor patterns were determined via dynamic optimization. The predictive gait simulation results were compared to experimental data from subjects with bilateral plantar flexor weakness walking with varying AFO stiffness. Our simulations demonstrated that reductions in mCoT with increasing stiffness were attributed to reductions in quadriceps metabolic cost during midstance. Increases in mCoT above optimum stiffness were attributed to the increasing metabolic cost of both hip flexor and hamstrings muscles. The insights gained from our predictive gait simulations could inform clinicians on the prescription of personalized AFOs. With further model individualization, simulations based on mCoT minimization may sufficiently predict adaptations to an AFO in individuals with plantar flexor weakness.

Telecoaching: a potential new training model for Charcot-Marie-Tooth patients: a systematic review

Introduction

Charcot-Marie-Tooth disease (CMT) is an inherited neuropathy that affects the sensory and motor nerves. It can be considered the most common neuromuscular disease, with a prevalence of 1/2500.

Methods

Considering the absence of a specific medical treatment and the benefits shown by physical activity in this population, a systematic review was completed using several search engines (Scopus, PubMed, and Web of Science) to analyze the use, effectiveness, and safety of a training program performed in telecoaching (TC). TC is a new training mode that uses mobile devices and digital technology to ensure remote access to training.

Results

Of the 382 studies identified, only 7 met the inclusion criteria. The effects of a TC training program included improvements in strength, cardiovascular ability, and functional abilities, as well as gait and fatigue. However, the quality of the studies was moderate, the size of the participants in each study was small, and the outcome measured was partial.

Discussion

Although many studies have identified statistically significant changes following the administration of the TC training protocol, further studies are needed, with appropriate study power, better quality, and a higher sample size.

Pole balancing on the fingertip: model-motivated machine learning forecasting of falls

Introduction: There is increasing interest in developing mathematical and computational models to forecast adverse events in physiological systems. Examples include falls, the onset of fatal cardiac arrhythmias, and adverse surgical outcomes. However, the dynamics of physiological systems are known to be exceedingly complex and perhaps even chaotic. Since no model can be perfect, it becomes important to understand how forecasting can be improved, especially when training data is limited. An adverse event that can be readily studied in the laboratory is the occurrence of stick falls when humans attempt to balance a stick on their fingertips. Over the last 20 years, this task has been extensively investigated experimentally, and presently detailed mathematical models are available. Methods: Here we use a long short-term memory (LTSM) deep learning network to forecast stick falls. We train this model to forecast stick falls in three ways: 1) using only data generated by the mathematical model (synthetic data), 2) using only stick balancing recordings of stick falls measured using high-speed motion capture measurements (human data), and 3) using transfer learning which combines a model trained using synthetic data plus a small amount of human balancing data. Results: We observe that the LTSM model is much more successful in forecasting a fall using synthetic data than it is in forecasting falls for models trained with limited available human data. However, with transfer learning, i.e., the LTSM model pre-trained with synthetic data and re-trained with a small amount of real human balancing data, the ability to forecast impending falls in human data is vastly improved. Indeed, it becomes possible to correctly forecast 60%-70% of real human stick falls up to 2.35 s in advance. Conclusion: These observations support the use of model-generated data and transfer learning techniques to improve the ability of computational models to forecast adverse physiological events.

Muscle strength and foot pressure vary depending on the type of foot pain

This study compared muscle strength and foot pressure among patients with metatarsalgia, patients with plantar fasciitis, and healthy controls. A total of 31 patients with foot pain (14 metatarsalgia and 17 plantar fasciitis) and 29 healthy controls participated in the study. The strengths of the plantar flexor and hip muscles were measured using isokinetic and handheld dynamometers, respectively. Foot pressure parameters, including the pressure-time integral (PTI) and foot arch index (AI), were assessed using pedobarography. Compared with the healthy control group, plantar flexor strength was significantly reduced in the affected feet of the metatarsalgia and plantar fasciitis groups (F = 0.083, all p < 0.001); however, hip strength was significantly decreased only in the affected feet of the metatarsalgia group (F = 20.900, p < 0.001). Plantar flexor (p < 0.001) and hip (p = 0.004) strength were significantly lower in the metatarsalgia group than in the plantar fasciitis group. The PTI was lower in the forefeet of the affected feet in the metatarsalgia (p < 0.001) and plantar fasciitis (p = 0.004) groups. Foot AI (p < 0.001) was significantly reduced only in the metatarsalgia group. These results suggest the need to consider the evaluation of muscle strength and foot pressure in both feet for the diagnosis and treatment of foot pain.

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.

Scalable musculoskeletal model for dynamic simulations of lower body movement

A musculoskeletal (MSK) model is an important tool for analysing human motions, calculating joint torques during movement, enhancing sports activity, and developing exoskeletons and prostheses. To enable biomechanical investigation of human motion, this work presents an open-source lower body MSK model. The MSK model of the lower body consists of 7 body segments (pelvis, left/right thigh, left/right leg, and left/right foot). The model has 20 degrees of freedom (DoFs) and 28 muscle torque generators (MTGs), which are developed from experimental data. The model can be modified for different anthropometric measurements and subject body characteristics, including sex, age, body mass, height, physical activity, and skin temperature. The model is validated by simulating the torque within the range of motion (ROM) of isolated movements; all simulation findings exhibit a good level of agreement with the literature.

NSF DARE-transforming modeling in neurorehabilitation: a patient-in-the-loop framework

In 2023, the National Science Foundation (NSF) and the National Institute of Health (NIH) brought together engineers, scientists, and clinicians by sponsoring a conference on computational modelling in neurorehabiilitation. To facilitate multidisciplinary collaborations and improve patient care, in this perspective piece we identify where and how computational modelling can support neurorehabilitation. To address the where, we developed a patient-in-the-loop framework that uses multiple and/or continual measurements to update diagnostic and treatment model parameters, treatment type, and treatment prescription, with the goal of maximizing clinically-relevant functional outcomes. This patient-in-the-loop framework has several key features: (i) it includes diagnostic and treatment models, (ii) it is clinically-grounded with the International Classification of Functioning, Disability and Health (ICF) and patient involvement, (iii) it uses multiple or continual data measurements over time, and (iv) it is applicable to a range of neurological and neurodevelopmental conditions. To address the how, we identify state-of-the-art and highlight promising avenues of future research across the realms of sensorimotor adaptation, neuroplasticity, musculoskeletal, and sensory & pain computational modelling. We also discuss both the importance of and how to perform model validation, as well as challenges to overcome when implementing computational models within a clinical setting. The patient-in-the-loop approach offers a unifying framework to guide multidisciplinary collaboration between computational and clinical stakeholders in the field of neurorehabilitation.

Predictive simulations identify potential neuromuscular contributors to idiopathic toe walking

BACKGROUND

Most cases of toe walking in children are idiopathic. We used pathology-specific neuromusculoskeletal predictive simulations to identify potential underlying neural and muscular mechanisms contributing to idiopathic toe walking.

METHODS

A musculotendon contracture was added to the ankle plantarflexors of a generic musculoskeletal model to represent a pathology-specific contracture model, matching the reduced ankle dorsiflexion range-of-motion in a cohort of children with idiopathic toe walking. This model was employed in a forward dynamic simulation controlled by reflexes and supraspinal drive, governed by a multi-objective cost function to predict gait patterns with the contracture model. We validated the predicted gait using experimental gait data from children with idiopathic toe walking with ankle contracture, by calculating the root mean square errors averaged over all biomechanical variables.

FINDINGS

A predictive simulation with the pathology-specific model with contracture approached experimental ITW data (root mean square error = 1.37SD). Gastrocnemius activation was doubled from typical gait simulations, but lacked a peak in early stance as present in electromyography. This synthesised idiopathic toe walking was more costly for all cost function criteria than typical gait simulation. Also, it employed a different neural control strategy, with increased length- and velocity-based reflex gains to the plantarflexors in early stance and swing than typical gait simulations.

INTERPRETATION

The simulations provide insights into how a musculotendon contracture combined with altered neural control could contribute to idiopathic toe walking. Insights into these neuromuscular mechanisms could guide future computational and experimental studies to gain improved insight into the cause of idiopathic toe walking.

Identifying essential factors for energy-efficient walking control across a wide range of velocities in reflex-based musculoskeletal systems

Humans can generate and sustain a wide range of walking velocities while optimizing their energy efficiency. Understanding the intricate mechanisms governing human walking will contribute to the engineering applications such as energy-efficient biped robots and walking assistive devices. Reflex-based control mechanisms, which generate motor patterns in response to sensory feedback, have shown promise in generating human-like walking in musculoskeletal models. However, the precise regulation of velocity remains a major challenge. This limitation makes it difficult to identify the essential reflex circuits for energy-efficient walking. To explore the reflex control mechanism and gain a better understanding of its energy-efficient maintenance mechanism, we extend the reflex-based control system to enable controlled walking velocities based on target speeds. We developed a novel performance-weighted least squares (PWLS) method to design a parameter modulator that optimizes walking efficiency while maintaining target velocity for the reflex-based bipedal system. We have successfully generated walking gaits from 0.7 to 1.6 m/s in a two-dimensional musculoskeletal model based on an input target velocity in the simulation environment. Our detailed analysis of the parameter modulator in a reflex-based system revealed two key reflex circuits that have a significant impact on energy efficiency. Furthermore, this finding was confirmed to be not influenced by setting parameters, i.e., leg length, sensory time delay, and weight coefficients in the objective cost function. These findings provide a powerful tool for exploring the neural bases of locomotion control while shedding light on the intricate mechanisms underlying human walking and hold significant potential for practical engineering applications.

Investigating the roles of reflexes and central pattern generators in the control and modulation of human locomotion using a physiologically plausible neuromechanical model

Objective.Studying the neural components regulating movement in human locomotion is obstructed by the inability to perform invasive experimental recording in the human neural circuits. Neuromechanical simulations can provide insights by modeling the locomotor circuits. Past neuromechanical models proposed control of locomotion either driven by central pattern generators (CPGs) with simple sensory commands or by a purely reflex-based network regulated by state-machine mechanisms, which activate and deactivate reflexes depending on the detected gait cycle phases. However, the physiological interpretation of these state machines remains unclear. Here, we present a physiologically plausible model to investigate spinal control and modulation of human locomotion.Approach.We propose a bio-inspired controller composed of two coupled CPGs that produce the rhythm and pattern, and a reflex-based network simulating low-level reflex pathways and Renshaw cells. This reflex network is based on leaky-integration neurons, and the whole system does not rely on changing reflex gains according to the gait cycle state. The musculoskeletal model is composed of a skeletal structure and nine muscles per leg generating movement in sagittal plane.Main results.Optimizing the open parameters for effort minimization and stability, human kinematics and muscle activation naturally emerged. Furthermore, when CPGs were not activated, periodic motion could not be achieved through optimization, suggesting the necessity of this component to generate rhythmic behavior without a state machine mechanism regulating reflex activation. The controller could reproduce ranges of speeds from 0.3 to 1.9 m s-1. The results showed that the net influence of feedback on motoneurons (MNs) during perturbed locomotion is predominantly inhibitory and that the CPGs provide the timing of MNs' activation by exciting or inhibiting muscles in specific gait phases.Significance.The proposed bio-inspired controller could contribute to our understanding of locomotor circuits of the intact spinal cord and could be used to study neuromotor disorders.

A spinal cord neuroprosthesis for locomotor deficits due to Parkinson's disease

People with late-stage Parkinson's disease (PD) often suffer from debilitating locomotor deficits that are resistant to currently available therapies. To alleviate these deficits, we developed a neuroprosthesis operating in closed loop that targets the dorsal root entry zones innervating lumbosacral segments to reproduce the natural spatiotemporal activation of the lumbosacral spinal cord during walking. We first developed this neuroprosthesis in a non-human primate model that replicates locomotor deficits due to PD. This neuroprosthesis not only alleviated locomotor deficits but also restored skilled walking in this model. We then implanted the neuroprosthesis in a 62-year-old male with a 30-year history of PD who presented with severe gait impairments and frequent falls that were medically refractory to currently available therapies. We found that the neuroprosthesis interacted synergistically with deep brain stimulation of the subthalamic nucleus and dopaminergic replacement therapies to alleviate asymmetry and promote longer steps, improve balance and reduce freezing of gait. This neuroprosthesis opens new perspectives to reduce the severity of locomotor deficits in people with PD.

Simulating human walking: a model-based reinforcement learning approach with musculoskeletal modeling

Introduction

Recent advancements in reinforcement learning algorithms have accelerated the development of control models with high-dimensional inputs and outputs that can reproduce human movement. However, the produced motion tends to be less human-like if algorithms do not involve a biomechanical human model that accounts for skeletal and muscle-tendon properties and geometry. In this study, we have integrated a reinforcement learning algorithm and a musculoskeletal model including trunk, pelvis, and leg segments to develop control modes that drive the model to walk.

Methods

We simulated human walking first without imposing target walking speed, in which the model was allowed to settle on a stable walking speed itself, which was 1.45 m/s. A range of other speeds were imposed for the simulation based on the previous self-developed walking speed. All simulations were generated by solving the Markov decision process problem with covariance matrix adaptation evolution strategy, without any reference motion data.

Results

Simulated hip and knee kinematics agreed well with those in experimental observations, but ankle kinematics were less well-predicted.

Discussion

We finally demonstrated that our reinforcement learning framework also has the potential to model and predict pathological gait that can result from muscle weakness.

Modeling Human Suboptimal Control: A Review

This review paper provides an overview of the approaches to model neuromuscular control, focusing on methods to identify nonoptimal control strategies typical of populations with neuromuscular disorders or children. Where possible, the authors tightened the description of the methods to the mechanisms behind the underlying biomechanical and physiological rationale. They start by describing the first and most simplified approach, the reductionist approach, which splits the role of the nervous and musculoskeletal systems. Static optimization and dynamic optimization methods and electromyography-based approaches are summarized to highlight their limitations and understand (the need for) their developments over time. Then, the authors look at the more recent stochastic approach, introduced to explore the space of plausible neural solutions, thus implementing the uncontrolled manifold theory, according to which the central nervous system only controls specific motions and tasks to limit energy consumption while allowing for some degree of adaptability to perturbations. Finally, they explore the literature covering the explicit modeling of the coupling between the nervous system (acting as controller) and the musculoskeletal system (the actuator), which may be employed to overcome the split characterizing the reductionist approach.

Predicting Gait Patterns of Children With Spasticity by Simulating Hyperreflexia

Spasticity is a common impairment within pediatric neuromusculoskeletal disorders. How spasticity contributes to gait deviations is important for treatment selection. Our aim was to evaluate the pathophysiological mechanisms underlying gait deviations seen in children with spasticity, using predictive simulations. A cluster analysis was performed to extract distinct gait patterns from experimental gait data of 17 children with spasticity to be used as comparative validation data. A forward dynamic simulation framework was employed to predict gait with either velocity- or force-based hyperreflexia. This framework entailed a generic musculoskeletal model controlled by reflexes and supraspinal drive, governed by a multiobjective cost function. Hyperreflexia values were optimized to enable the simulated gait to best match experimental gait patterns. Three experimental gait patterns were extracted: (1) increased knee flexion, (2) increased ankle plantar flexion, and (3) increased knee flexion and ankle plantar flexion when compared with typical gait. Overall, velocity-based hyperreflexia outperformed force-based hyperreflexia. The first gait pattern could mostly be explained by rectus femoris and hamstrings velocity-based hyperreflexia, the second by gastrocnemius velocity-based hyperreflexia, and the third by gastrocnemius, soleus, and hamstrings velocity-based hyperreflexia. This study shows how velocity-based hyperreflexia from specific muscles contributes to different spastic gait patterns, which may help in providing targeted treatment.

Novel Design and Implementation of a Neuromuscular Controller on a Hip Exoskeleton for Partial Gait Assistance

Exoskeletons intended for partial assistance of walking should be able to follow the gait pattern of their users, via online adaptive control strategies rather than imposing predefined kinetic or kinematic profiles. NeuroMuscular Controllers (NMCs) are adaptive strategies inspired by the neuromuscular modeling methods that seek to mimic and replicate the behavior of the human nervous system and skeletal muscles during gait. This study presents a novel design of a NMC, applied for the first time to partial assistance hip exoskeletons. Rather than the two-phase (stance/swing) division used in previous designs for the modulation of reflexes, a 5-state finite state machines is designed for gait phase synchronisation. The common virtual muscle model is also modified by assuming a stiff tendon, allowing for a more analytical computation approach for the muscle state resolution. As a first validation, the performance of the controller was tested with 9 healthy subjects walking at different speeds and slopes on a treadmill. The generated torque profiles show similarity to biological torques and optimal assistance profiles in the literature. Power output profiles of the exoskeleton indicate good synchronization with the users' intended movements, reflected in predominantly positive work by the assistance. The results also highlight the adaptability of the controller to different users and walking conditions, without the need for extensive parameter tuning.

Muscle strength and spatiotemporal gait parameters in people with Parkinson´s disease. A pilot study

INTRODUCTION

Muscle weakness in persons with Parkinson disease (PD) has been frequently recognized as a nonspecific symptom. In other neurological conditions, lower limb weakness, specifically quadriceps weakness, is the factor that causes greater gait disability. Little research has evaluated the relationship between lower limb muscle strength, using objective tools, in PD persons and gait performance. The aim of this study was to analyze the correlation between lower limb muscle strength, using an isokinetic dynamometer, and the spatiotemporal gait parameters in PD, compared with age- and sex- matched healthy controls.

SUBJECTS AND METHODS

The study was conducted with 7 persons with PD -Hoehn and Yahr (HY) between II-III- and 7 healthy controls. Isokinetic knee and ankle tests at 60 and 120°/s and the 10-meter walking test at comfortable and fast walking speed, were performed on all recruited subjects.

RESULTS

Significant differences in lower limb strength-related measures and gait parameters were observed between persons with PD and controls. Gait parameters showed excellent correlations (rho = 0.7) for both lower limb: ankle plantar flexion work/body wearing at 180°/s with number of steps (indirect) and stride (direct) at both speeds, and between the ankle plantar flexion peak torque/ body wearing at 180°/s with number of steps (indirect) and stride (direct) at maximum speed; and between knee extension work/body wearing at 60°/s) with stride (direct) at self-selected speed.

CONCLUSIONS

Persons with PD (HY II-III stages) lower limb muscle strength correlates excellently with gait pattern, showing lower isokinetic strength than healthy subjects of the same age and sex. This protocol showed safety to be performed in a larger sample.

Simulating the effect of ankle plantarflexion and inversion-eversion exoskeleton torques on center of mass kinematics during walking

Walking balance is central to independent mobility, and falls due to loss of balance are a leading cause of death for people 65 years of age and older. Bipedal gait is typically unstable, but healthy humans use corrective torques to counteract perturbations and stabilize gait. Exoskeleton assistance could benefit people with neuromuscular deficits by providing stabilizing torques at lower-limb joints to replace lost muscle strength and sensorimotor control. However, it is unclear how applied exoskeleton torques translate to changes in walking kinematics. This study used musculoskeletal simulation to investigate how exoskeleton torques applied to the ankle and subtalar joints alter center of mass kinematics during walking. We first created muscle-driven walking simulations using OpenSim Moco by tracking experimental kinematics and ground reaction forces recorded from five healthy adults. We then used forward integration to simulate the effect of exoskeleton torques applied to the ankle and subtalar joints while keeping muscle excitations fixed based on our previous tracking simulation results. Exoskeleton torque lasted for 15% of the gait cycle and was applied between foot-flat and toe-off during the stance phase, and changes in center of mass kinematics were recorded when the torque application ended. We found that changes in center of mass kinematics were dependent on both the type and timing of exoskeleton torques. Plantarflexion torques produced upward and backward changes in velocity of the center of mass in mid-stance and upward and smaller forward velocity changes near toe-off. Eversion and inversion torques primarily produced lateral and medial changes in velocity in mid-stance, respectively. Intrinsic muscle properties reduced kinematic changes from exoskeleton torques. Our results provide mappings between ankle plantarflexion and inversion-eversion torques and changes in center of mass kinematics which can inform designers building exoskeletons aimed at stabilizing balance during walking. Our simulations and software are freely available and allow researchers to explore the effects of applied torques on balance and gait.

Integration of feedforward and feedback control in the neuromechanics of vertebrate locomotion: a review of experimental, simulation and robotic studies

Animal locomotion is the result of complex and multi-layered interactions between the nervous system, the musculo-skeletal system and the environment. Decoding the underlying mechanisms requires an integrative approach. Comparative experimental biology has allowed researchers to study the underlying components and some of their interactions across diverse animals. These studies have shown that locomotor neural circuits are distributed in the spinal cord, the midbrain and higher brain regions in vertebrates. The spinal cord plays a key role in locomotor control because it contains central pattern generators (CPGs) - systems of coupled neuronal oscillators that provide coordinated rhythmic control of muscle activation that can be viewed as feedforward controllers - and multiple reflex loops that provide feedback mechanisms. These circuits are activated and modulated by descending pathways from the brain. The relative contributions of CPGs, feedback loops and descending modulation, and how these vary between species and locomotor conditions, remain poorly understood. Robots and neuromechanical simulations can complement experimental approaches by testing specific hypotheses and performing what-if scenarios. This Review will give an overview of key knowledge gained from comparative vertebrate experiments, and insights obtained from neuromechanical simulations and robotic approaches. We suggest that the roles of CPGs, feedback loops and descending modulation vary among animals depending on body size, intrinsic mechanical stability, time required to reach locomotor maturity and speed effects. We also hypothesize that distal joints rely more on feedback control compared with proximal joints. Finally, we highlight important opportunities to address fundamental biological questions through continued collaboration between experimentalists and engineers.

Dysfunctional neuro-muscular mechanisms explain gradual gait changes in prodromal spastic paraplegia

BACKGROUND

In Hereditary Spastic Paraplegia (HSP) type 4 (SPG4) a length-dependent axonal degeneration in the cortico-spinal tract leads to progressing symptoms of hyperreflexia, muscle weakness, and spasticity of lower extremities. Even before the manifestation of spastic gait, in the prodromal phase, axonal degeneration leads to subtle gait changes. These gait changes - depicted by digital gait recording - are related to disease severity in prodromal and early-to-moderate manifest SPG4 participants.

METHODS

We hypothesize that dysfunctional neuro-muscular mechanisms such as hyperreflexia and muscle weakness explain these disease severity-related gait changes of prodromal and early-to-moderate manifest SPG4 participants. We test our hypothesis in computer simulation with a neuro-muscular model of human walking. We introduce neuro-muscular dysfunction by gradually increasing sensory-motor reflex sensitivity based on increased velocity feedback and gradually increasing muscle weakness by reducing maximum isometric force.

RESULTS

By increasing hyperreflexia of plantarflexor and dorsiflexor muscles, we found gradual muscular and kinematic changes in neuro-musculoskeletal simulations that are comparable to subtle gait changes found in prodromal SPG4 participants.

CONCLUSIONS

Predicting kinematic changes of prodromal and early-to-moderate manifest SPG4 participants by gradual alterations of sensory-motor reflex sensitivity allows us to link gait as a directly accessible performance marker to emerging neuro-muscular changes for early therapeutic interventions.

Ten steps to becoming a musculoskeletal simulation expert: A half-century of progress and outlook for the future

Over the past half-century, musculoskeletal simulations have deepened our knowledge of human and animal movement. This article outlines ten steps to becoming a musculoskeletal simulation expert so you can contribute to the next half-century of technical innovation and scientific discovery. We advocate looking to the past, present, and future to harness the power of simulations that seek to understand and improve mobility. Instead of presenting a comprehensive literature review, we articulate a set of ideas intended to help researchers use simulations effectively and responsibly by understanding the work on which today's musculoskeletal simulations are built, following established modeling and simulation principles, and branching out in new directions.

The neuromechanics of animal locomotion: From biology to robotics and back

Robotics and neuroscience are sister disciplines that both aim to understand how agile, efficient, and robust locomotion can be achieved in autonomous agents. Robotics has already benefitted from neuromechanical principles discovered by investigating animals. These include the use of high-level commands to control low-level central pattern generator-like controllers, which, in turn, are informed by sensory feedback. Reciprocally, neuroscience has benefited from tools and intuitions in robotics to reveal how embodiment, physical interactions with the environment, and sensory feedback help sculpt animal behavior. We illustrate and discuss exemplar studies of this dialog between robotics and neuroscience. We also reveal how the increasing biorealism of simulations and robots is driving these two disciplines together, forging an integrative science of autonomous behavioral control with many exciting future opportunities.

Minimum effort simulations of split-belt treadmill walking exploit asymmetry to reduce metabolic energy expenditure

Walking on a split-belt treadmill elicits an adaptation response that changes baseline step length asymmetry. The underlying causes of this adaptation, however, are difficult to determine. It has been proposed that effort minimization may drive this adaptation, based on the idea that adopting longer steps on the fast belt, or positive step length asymmetry (SLA), can cause the treadmill to exert net-positive mechanical work on a bipedal walker. However, humans walking on split-belt treadmills have not been observed to reproduce this behavior when allowed to freely adapt. To determine if an effort-minimization motor control strategy would result in experimentally observed adaptation patterns, we conducted simulations of walking on different combinations of belt speeds with a human musculoskeletal model that minimized muscle excitations and metabolic rate. The model adopted increasing amounts of positive SLA and decreased its net metabolic rate with increasing belt speed difference, reaching +42.4% SLA and -5.7% metabolic rate relative to tied-belt walking at our maximum belt speed ratio of 3:1. These gains were primarily enabled by an increase in braking work and a reduction in propulsion work on the fast belt. The results suggest that a purely effort minimization driven split-belt walking strategy would involve substantial positive SLA, and that the lack of this characteristic in human behavior points to additional factors influencing the motor control strategy, such as aversion to excessive joint loads, asymmetry, or instability.NEW & NOTEWORTHY Behavioral observations of split-belt treadmill adaptation have been inconclusive toward its underlying causes. To estimate gait patterns when driven exclusively by one of these possible underlying causes, we simulated split-belt treadmill walking with a musculoskeletal model that minimized its summed muscle excitations. Our model took significantly longer steps on the fast belt and reduced its metabolic rate below tied-belt walking, unlike experimental observations. This suggests that asymmetry is energetically optimal, but human adaptation involves additional factors.

Robust walking control of a lower limb rehabilitation exoskeleton coupled with a musculoskeletal model via deep reinforcement learning

BACKGROUND

Few studies have systematically investigated robust controllers for lower limb rehabilitation exoskeletons (LLREs) that can safely and effectively assist users with a variety of neuromuscular disorders to walk with full autonomy. One of the key challenges for developing such a robust controller is to handle different degrees of uncertain human-exoskeleton interaction forces from the patients. Consequently, conventional walking controllers either are patient-condition specific or involve tuning of many control parameters, which could behave unreliably and even fail to maintain balance.

METHODS

We present a novel, deep neural network, reinforcement learning-based robust controller for a LLRE based on a decoupled offline human-exoskeleton simulation training with three independent networks, which aims to provide reliable walking assistance against various and uncertain human-exoskeleton interaction forces. The exoskeleton controller is driven by a neural network control policy that acts on a stream of the LLRE's proprioceptive signals, including joint kinematic states, and subsequently predicts real-time position control targets for the actuated joints. To handle uncertain human interaction forces, the control policy is trained intentionally with an integrated human musculoskeletal model and realistic human-exoskeleton interaction forces. Two other neural networks are connected with the control policy network to predict the interaction forces and muscle coordination. To further increase the robustness of the control policy to different human conditions, we employ domain randomization during training that includes not only randomization of exoskeleton dynamics properties but, more importantly, randomization of human muscle strength to simulate the variability of the patient's disability. Through this decoupled deep reinforcement learning framework, the trained controller of LLREs is able to provide reliable walking assistance to patients with different degrees of neuromuscular disorders without any control parameter tuning.

RESULTS AND CONCLUSION

A universal, RL-based walking controller is trained and virtually tested on a LLRE system to verify its effectiveness and robustness in assisting users with different disabilities such as passive muscles (quadriplegic), muscle weakness, or hemiplegic conditions without any control parameter tuning. Analysis of the RMSE for joint tracking, CoP-based stability, and gait symmetry shows the effectiveness of the controller. An ablation study also demonstrates the strong robustness of the control policy under large exoskeleton dynamic property ranges and various human-exoskeleton interaction forces. The decoupled network structure allows us to isolate the LLRE control policy network for testing and sim-to-real transfer since it uses only proprioception information of the LLRE (joint sensory state) as the input. Furthermore, the controller is shown to be able to handle different patient conditions without the need for patient-specific control parameter tuning.

Predictive simulation of single-leg landing scenarios for ACL injury risk factors evaluation

The Anterior Cruciate Ligament (ACL) rupture is a very common knee injury during sport activities. Landing after jump is one of the most prominent human body movements that can lead to such an injury. The landing-related ACL injury risk factors have been in the spotlight of research interest. Over the years, researchers and clinicians acquire knowledge about human movement during daily-life activities by organizing complex in vivo studies that feature high complexity, costs and technical and most importantly physical challenges. In an attempt to overcome these limitations, this paper introduces a computational modeling and simulation pipeline that aims to predict and identify key parameters of interest that are related to ACL injury during single-leg landings. We examined the following conditions: a) landing height, b) hip internal and external rotation, c) lumbar forward and backward leaning, d) lumbar medial and lateral bending, e) muscle forces permutations and f) effort goal weight. Identified on related research studies, we evaluated the following risk factors: vertical Ground Reaction Force (vGRF), knee joint Anterior force (AF), Medial force (MF), Compressive force (CF), Abduction moment (AbdM), Internal rotation moment (IRM), quadricep and hamstring muscle forces and Quadriceps/Hamstrings force ratio (Q/H force ratio). Our study clearly demonstrated that ACL injury is a rather complicated mechanism with many associated risk factors which are evidently correlated. Nevertheless, the results were mostly in agreement with other research studies regarding the ACL risk factors. The presented pipeline showcased promising potential of predictive simulations to evaluate different aspects of complicated phenomena, such as the ACL injury.

Age-related changes in gait biomechanics and their impact on the metabolic cost of walking: Report from a National Institute on Aging workshop

Changes in old age that contribute to the complex issue of an increased metabolic cost of walking (mass-specific energy cost per unit distance traveled) in older adults appear to center at least in part on changes in gait biomechanics. However, age-related changes in energy metabolism, neuromuscular function and connective tissue properties also likely contribute to this problem, of which the consequences are poor mobility and increased risk of inactivity-related disease and disability. The U.S. National Institute on Aging convened a workshop in September 2021 with an interdisciplinary group of scientists to address the gaps in research related to the mechanisms and consequences of changes in mobility in old age. The goal of the workshop was to identify promising ways to move the field forward toward improving gait performance, decreasing energy cost, and enhancing mobility for older adults. This report summarizes the workshop and brings multidisciplinary insight into the known and potential causes and consequences of age-related changes in gait biomechanics. We highlight how gait mechanics and energy cost change with aging, the potential neuromuscular mechanisms and role of connective tissue in these changes, and cutting-edge interventions and technologies that may be used to measure and improve gait and mobility in older adults. Key gaps in the literature that warrant targeted research in the future are identified and discussed.

Tolerable degree of muscle sacrifice when harvesting a vastus lateralis or myocutaneous anterolateral thigh flap

The myocutaneous anterolateral thigh (ALT) and vastus lateralis (VL) flaps include a large muscle mass and a sufficient vascular pedicle, and they have been used for decades to reconstruct traumatic and acquired defects of the head and neck and extremities. In spite of these benefits, musculoskeletal dysfunction was reported in nearly 1 out of 20 patients at follow-up. It is unclear whether the recently proposed muscle-sparing flap-raising approach could preserve VL muscle function and whether patients at increased risk could benefit from such an approach. Therefore, we performed a predictive dynamic gait simulation based on a biological motion model with gradual weakening of the VL during a self-selected and fast walking speed to determine the compensable degree of VL muscle reduction. Muscle force, joint angle, and joint moment were measured. Our study showed that VL muscle reduction could be compensated up to a certain degree, which could explain the observed incidence of musculoskeletal dysfunction. In elderly or fragile patients, the VL muscle should not be reduced by 50% or more, which could be achieved by muscle-sparing flap-raising of the superficial partition only. In young or athletic patients, a VL muscle reduction of 10%, which corresponds to a muscle cuff, has no relevant effect. Yet, a reduction of more than 30% leads to relevant weakening of the quadriceps. Therefore, in this patient population with the need for a large portion of muscle, alternative flaps should be considered. This study can serve as the first basis for further investigations of human locomotion after flap-raising.

Understanding altered contractile properties in advanced age: insights from a systematic muscle modelling approach

Age-related alterations of skeletal muscle are numerous and present inconsistently, and the effect of their interaction on contractile performance can be nonintuitive. Hill-type muscle models predict muscle force according to well-characterised contractile phenomena. Coupled with simple, yet reasonably realistic activation dynamics, such models consist of parameters that are meaningfully linked to fundamental aspects of muscle excitation and contraction. We aimed to illustrate the utility of a muscle model for elucidating relevant mechanisms and predicting changes in output by simulating the individual and combined effects on isometric force of several known ageing-related adaptations. Simulating literature-informed reductions in free Ca2+ concentration and Ca2+ sensitivity generated predictions at odds qualitatively with the characteristic slowing of contraction speed. Conversely, incorporating slower Ca2+ removal or a fractional increase in type I fibre area emulated expected changes; the former was required to simulate slowing of the twitch measured experimentally. Slower Ca2+ removal more than compensated for force loss arising from a large reduction in Ca2+ sensitivity or moderate reduction in Ca2+ release, producing realistic age-related shifts in the force-frequency relationship. Consistent with empirical data, reductions in free Ca2+ concentration and Ca2+ sensitivity reduced maximum tetanic force only slightly, even when acting in concert, suggesting a modest contribution to lower specific force. Lower tendon stiffness and slower intrinsic shortening speed slowed and prolonged force development in a compliance-dependent manner without affecting force decay. This work demonstrates the advantages of muscle modelling for exploring sources of variation and identifying mechanisms underpinning the altered contractile properties of aged muscle.

High-level motor planning allows flexible walking at different gait patterns in a neuromechanical model

Humans can freely adopt gait parameters like walking speed, step length, or cadence on the fly when walking. Planned movement that can be updated online to account for changes in the environment rather than having to rely on habitual, reflexive control that is adapted over long timescales. Here we present a neuromechanical model that accounts for this flexibility by combining movement goals and motor plans on a kinematic task level with low-level spinal feedback loops. We show that the model can walk at a wide range of different gait patterns by choosing a small number of high-level control parameters representing a movement goal. A larger number of parameters governing the low-level reflex loops in the spinal cord, on the other hand, remain fixed. We also show that the model can generalize the learned behavior by recombining the high-level control parameters and walk with gait patterns that it had not encountered before. Furthermore, the model can transition between different gaits without the loss of balance by switching to a new set of control parameters in real time.

Patterns of asymmetry and energy cost generated from predictive simulations of hemiparetic gait

Hemiparesis, defined as unilateral muscle weakness, often occurs in people post-stroke or people with cerebral palsy, however it is difficult to understand how this hemiparesis affects movement patterns as it often presents alongside a variety of other neuromuscular impairments. Predictive musculoskeletal modeling presents an opportunity to investigate how impairments affect gait performance assuming a particular cost function. Here, we use predictive simulation to quantify the spatiotemporal asymmetries and changes to metabolic cost that emerge when muscle strength is unilaterally reduced and how reducing spatiotemporal symmetry affects metabolic cost. We modified a 2-D musculoskeletal model by uniformly reducing the peak isometric muscle force unilaterally. We then solved optimal control simulations of walking across a range of speeds by minimizing the sum of the cubed muscle excitations. Lastly, we ran additional optimizations to test if reducing spatiotemporal asymmetry would result in an increase in metabolic cost. Our results showed that the magnitude and direction of effort-optimal spatiotemporal asymmetries depends on both the gait speed and level of weakness. Also, the optimal speed was 1.25 m/s for the symmetrical and 20% weakness models but slower (1.00 m/s) for the 40% and 60% weakness models, suggesting that hemiparesis can account for a portion of the slower gait speed seen in people with hemiparesis. Modifying the cost function to minimize spatiotemporal asymmetry resulted in small increases (~4%) in metabolic cost. Overall, our results indicate that spatiotemporal asymmetry may be optimal for people with hemiparesis. Additionally, the effect of speed and the level of weakness on spatiotemporal asymmetry may help explain the well-known heterogenous distribution of spatiotemporal asymmetries observed in the clinic. Future work could extend our results by testing the effects of other neuromuscular impairments on optimal gait strategies, and therefore build a more comprehensive understanding of the gait patterns observed in clinical populations.

Policy Design for an Ankle-Foot Orthosis Using Simulated Physical Human-Robot Interaction via Deep Reinforcement Learning

This paper presents a novel approach for designing a robotic orthosis controller considering physical human-robot interaction (pHRI). Computer simulation for this human-robot system can be advantageous in terms of time and cost due to the laborious nature of designing a robot controller that effectively assists humans with the appropriate magnitude and phase. Therefore, we propose a two-stage policy training framework based on deep reinforcement learning (deep RL) to design a robot controller using human-robot dynamic simulation. In Stage 1, the optimal policy of generating human gaits is obtained from deep RL-based imitation learning on a healthy subject model using the musculoskeletal simulation in OpenSim-RL. In Stage 2, human models in which the right soleus muscle is weakened to a certain severity are created by modifying the human model obtained from Stage 1. A robotic orthosis is then attached to the right ankle of these models. The orthosis policy that assists walking with optimal torque is then trained on these models. Here, the elastic foundation model is used to predict the pHRI in the coupling part between the human and robotic orthosis. Comparative analysis of kinematic and kinetic simulation results with the experimental data shows that the derived human musculoskeletal model imitates a human walking. It also shows that the robotic orthosis policy obtained from two-stage policy training can assist the weakened soleus muscle. The proposed approach was validated by applying the learned policy to ankle orthosis, conducting a gait experiment, and comparing it with the simulation results.

Fall Prediction Based on Instrumented Measures of Gait and Turning in Daily Life in People with Multiple Sclerosis

This study investigates the potential of passive monitoring of gait and turning in daily life in people with multiple sclerosis (PwMS) to identify those at future risk of falls. Seven days of passive monitoring of gait and turning were carried out in a pilot study of 26 PwMS in home settings using wearable inertial sensors. The retrospective fall history was collected at the baseline. After gait and turning data collection in daily life, PwMS were followed biweekly for a year and were classified as fallers if they experienced >1 fall. The ability of short-term passive monitoring of gait and turning, as well as retrospective fall history to predict future falls were compared using receiver operator curves and regression analysis. The history of retrospective falls was not identified as a significant predictor of future falls in this cohort (AUC = 0.62, p = 0.32). Among quantitative monitoring measures of gait and turning, the pitch at toe-off was the best predictor of falls (AUC = 0.86, p < 0.01). Fallers had a smaller pitch of their feet at toe-off, reflecting less plantarflexion during the push-off phase of walking, which can impact forward propulsion and swing initiation and can result in poor foot clearance and an increased metabolic cost of walking. In conclusion, our cohort of PwMS showed that objective monitoring of gait and turning in daily life can identify those at future risk of falls, and the pitch at toe-off was the single most influential predictor of future falls. Therefore, interventions aimed at improving the strength of plantarflexion muscles, range of motion, and increased proprioceptive input may benefit PwMS at future fall risk.

Characterizing stroke-induced changes in the variability of lower limb kinematics using multifractal detrended fluctuation analysis

Movement variability reflects the adaptation of the neuromuscular control system to internal or external perturbations, but its relationship to stroke-induced injury is still unclear. In this study, the multifractal detrended fluctuation analysis was used to explore the stroke-induced changes in movement variability by analyzing the joint angles in a treadmill-walking task. Eight healthy subjects and ten patients after stroke participated in the experiment, performing a treadmill-walking task at a comfortable speed. The kinematics data of the lower limbs were collected by the motion-capture system, and two indicators, the degree of multifractality (α) and degree of correlation [h(2)], were used to investigate the mechanisms underlying neuromuscular control. The results showed that the knee and ankle joint angles were multifractal and persistent at various scales, and there was a significant difference in the degree of multifractality and the degree of correlation at the knee and ankle joint angles among the three groups, with the values being ranked in the following order: healthy subjects < non-paretic limb < paretic limb. These observations highlighted increased movement variability and multifractal strength in patients after stroke due to neuromotor defects. This study provided evidence that multifractal detrended analysis of the angles of the knee and ankle joints is useful to investigate the changes in movement variability and multifractal after stroke. Further research is needed to verify and promote the clinical applications.

Development and Validation of a Closed-Loop Functional Electrical Stimulation-Based Controller for Gait Rehabilitation Using a Finite State Machine Model

Functional electrical stimulation (FES) can be used to initiate lower limb muscle contractions and has been widely applied in gait rehabilitation. Establishing the correct timing of FES activation during each phase of the gait (walking) cycle remains challenging as most FES systems rely on open-loop control, whereby the controller receives no feedback about joint kinematics and instead relies on predetermined/timed muscle stimulation. The objective of this study was to develop and validate a closed-loop FES-based control solution for gait rehabilitation using a finite state machine (FSM) model. A two-phased study approach was taken: (1) Experimentally-Informed Study: A neuromuscular-derived FSM model was developed to drive closed-loop FES-based control for gait rehabilitation. The finite states were determined using electromyography and joint kinematics data of 12 non-disabled adults, collected during treadmill walking. The gait cycles were divided into four states, namely: swing-to-stance, push off, pre-swing, and toe up. (2) Simulation Study: A closed-loop FES-based control solution that employed the resulting FSM model, was validated through comparisons of neuro-musculo-skeletal computer simulations of impaired versus healthy gait. This closed-loop controller yielded steadier simulated impaired gait, in comparison to an open-loop alternative. The simulation results confirmed that accurate timing of FES activation during the gait cycle, as informed by kinematics data, is important to natural gait retraining. The closed-loop FES-based solution, introduced in this study, contributes to the repository of gait rehabilitation control options and offers the advantage of being simplistic to implement. Furthermore, this control solution is expected to integrate well with powered exoskeleton technologies.

Adaptive ankle exoskeleton gait training demonstrates acute neuromuscular and spatiotemporal benefits for individuals with cerebral palsy: A pilot study

BACKGROUND

Gait abnormalities from neuromuscular conditions like cerebral palsy (CP) limit mobility and negatively affect quality of life. Increasing walking speed and stride length are essential clinical goals in the treatment of gait disorders from CP.

RESEARCH QUESTION

How does over-ground gait training with an untethered ankle exoskeleton providing adaptive assistance affect mobility-related spatiotemporal outcomes and lower-extremity muscle activity in people with CP?

METHODS

A diverse cohort of individuals with CP (n = 6, age 9-31, Gross Motor Function Classification System Level I - III) completed four over-ground training sessions (98 ± 17 min of assisted walking) and received pre- and post-training assessments. On both assessments, participants walked over-ground with and without the exoskeleton while we recorded spatiotemporal outcomes and muscle activity. We used two-tailed paired t-tests to compare all parameters pre- and post-training, and between assisted and unassisted conditions.

RESULTS

Following training, walking speed increased 0.24 m/s (p = 0.006) and stride length increased 0.17 m (p = 0.013) during unassisted walking, while walking speed increased 0.28 m/s (p = 0.023) and stride length increased 0.15 m (p = 0.002) during exoskeleton-assisted walking. Exoskeleton training improved stride-to-stride repeatability of soleus and vastus lateralis muscle activation by up to 51 % (p ≤ 0.046), while the amount of integrated stance-phase muscle activity was similar across visits and conditions. Relative to baseline, post-training walking with the exoskeleton resulted in a soleus activity pattern that was 39 % more similar to the typical pattern from unimpaired individuals (p < 0.001).

SIGNIFICANCE

This study demonstrates acute spatiotemporal and neuromuscular benefits from over-ground training with adaptive ankle exoskeleton assistance, and provides rationale for completion of a longer randomized controlled training protocol.

An approximate stochastic optimal control framework to simulate nonlinear neuro-musculoskeletal models in the presence of noise

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.

A neuromuscular model of human locomotion combines spinal reflex circuits with voluntary movements

Existing models of human walking use low-level reflexes or neural oscillators to generate movement. While appropriate to generate the stable, rhythmic movement patterns of steady-state walking, these models lack the ability to change their movement patterns or spontaneously generate new movements in the specific, goal-directed way characteristic of voluntary movements. Here we present a neuromuscular model of human locomotion that bridges this gap and combines the ability to execute goal directed movements with the generation of stable, rhythmic movement patterns that are required for robust locomotion. The model represents goals for voluntary movements of the swing leg on the task level of swing leg joint kinematics. Smooth movements plans towards the goal configuration are generated on the task level and transformed into descending motor commands that execute the planned movements, using internal models. The movement goals and plans are updated in real time based on sensory feedback and task constraints. On the spinal level, the descending commands during the swing phase are integrated with a generic stretch reflex for each muscle. Stance leg control solely relies on dedicated spinal reflex pathways. Spinal reflexes stimulate Hill-type muscles that actuate a biomechanical model with eight internal joints and six free-body degrees of freedom. The model is able to generate voluntary, goal-directed reaching movements with the swing leg and combine multiple movements in a rhythmic sequence. During walking, the swing leg is moved in a goal-directed manner to a target that is updated in real-time based on sensory feedback to maintain upright balance, while the stance leg is stabilized by low-level reflexes and a behavioral organization switching between swing and stance control for each leg. With this combination of reflex-based stance leg and voluntary, goal-directed control of the swing leg, the model controller generates rhythmic, stable walking patterns in which the swing leg movement can be flexibly updated in real-time to step over or around obstacles.

Design of a control framework for lower limb exoskeleton rehabilitation robot based on predictive assessment

BACKGROUND

Patients suffering from lower limb dyskinesia, especially in early stages of rehabilitation, have weak residual muscle strength in affected limb and require passive training by the lower limb rehabilitation robot. Anatomy indicates that the biceps femoris short head muscle has a strong influence on knee motion at the swing phase of walking. We sought to explore how it would influence on gait cycle in optimization framework. However, the training trajectory of conventional rehabilitation robots performing passive training usually follows gait planning based on general human gait data, which cannot simultaneously ensure both effective rehabilitation of affected limbs with varying severity pathological gait and comfort of the wearer within a safe motion trajectory.

METHODS

To elucidate the effects of weakness and contracture, we systematically introduced isolated defects into the musculoskeletal model and generated walking simulations to predict the gait adaptation due to these defects. An impedance control model of the rehabilitation robot is developed. Knee joint parameters optimized by predictive forward dynamics simulation are adopted as the expected values for the robot controller to achieve customized adjustment of the robot motion trajectory.

FINDINGS

Severe muscle contracture leads to severe knee flexion; severe muscle weakness induces a significant posterior tilt of the upper trunk, which hinders walking speed.

INTERPRETATION

Our simulation results attempt to reveal pathological gait features, which may help to reproduce the simulation of pathological gait. Furthermore, the robot simulation results show that the robot system achieves a speedy tracking by setting a larger stiffness value. The model also allows the implementation of different levels of damping or elasticity effects.

TRIAL REGISTRATION

The method proposed in this paper is an initial basic study that did not reach clinical trials and therefore retains retrospectively registered.