Computational Design of FastFES Treatment to Improve Propulsive Force Symmetry During Post-stroke Gait: A Feasibility Study

Stroke is a leading cause of long-term disability worldwide and often impairs walking ability. To improve recovery of walking function post-stroke, researchers have investigated the use of treatments such as fast functional electrical stimulation (FastFES). During FastFES treatments, individuals post-stroke walk on a treadmill at their fastest comfortable speed while electrical stimulation is delivered to two muscles of the paretic ankle, ideally to improve paretic leg propulsion and toe clearance. However, muscle selection and stimulation timing are currently standardized based on clinical intuition and a one-size-fits-all approach, which may explain in part why some patients respond to FastFES training while others do not. This study explores how personalized neuromusculoskeletal models could potentially be used to enable individual-specific selection of target muscles and stimulation timing to address unique functional limitations of individual patients post-stroke. Treadmill gait data, including EMG, surface marker positions, and ground reactions, were collected from an individual post-stroke who was a non-responder to FastFES treatment. The patient's gait data were used to personalize key aspects of a full-body neuromusculoskeletal walking model, including lower-body joint functional axes, lower-body muscle force generating properties, deformable foot-ground contact properties, and paretic and non-paretic leg neural control properties. The personalized model was utilized within a direct collocation optimal control framework to reproduce the patient's unstimulated treadmill gait data (verification problem) and to generate three stimulated walking predictions that sought to minimize inter-limb propulsive force asymmetry (prediction problems). The three predictions used: (1) Standard muscle selection (gastrocnemius and tibialis anterior) with standard stimulation timing, (2) Standard muscle selection with optimized stimulation timing, and (3) Optimized muscle selection (soleus and semimembranosus) with optimized stimulation timing. Relative to unstimulated walking, the optimal control problems predicted a 41% reduction in propulsive force asymmetry for scenario (1), a 45% reduction for scenario (2), and a 64% reduction for scenario (3), suggesting that non-standard muscle selection may be superior for this patient. Despite these predicted improvements, kinematic symmetry was not noticeably improved for any of the walking predictions. These results suggest that personalized neuromusculoskeletal models may be able to predict personalized FastFES training prescriptions that could improve propulsive force symmetry, though inclusion of kinematic requirements would be necessary to improve kinematic symmetry as well.

The effect of plate design, bridging span, and fracture healing on the performance of high tibial osteotomy plates: An experimental and finite element study

Objectives

Opening wedge high tibial osteotomy (HTO) is an established surgical procedure for the treatment of early-stage knee arthritis. Other than infection, the majority of complications are related to mechanical factors – in particular, stimulation of healing at the osteotomy site. This study used finite element (FE) analysis to investigate the effect of plate design and bridging span on interfragmentary movement (IFM) and the influence of fracture healing on plate stress and potential failure.

Materials and Methods

A 10° opening wedge HTO was created in a composite tibia. Imaging and strain gauge data were used to create and validate FE models. Models of an intact tibia and a tibia implanted with a custom HTO plate using two different bridging spans were validated against experimental data. Physiological muscle forces and different stages of osteotomy gap healing simulating up to six weeks postoperatively were then incorporated. Predictions of plate stress and IFM for the custom plate were compared against predictions for an industry standard plate (TomoFix).

Results

For both plate types, long spans increased IFM but did not substantially alter peak plate stress. The custom plate increased axial and shear IFM values by up to 24% and 47%, respectively, compared with the TomoFix. In all cases, a callus stiffness of 528 MPa was required to reduce plate stress below the fatigue strength of titanium alloy.

Conclusion

We demonstrate that larger bridging spans in opening wedge HTO increase IFM without substantially increasing plate stress. The results indicate, however, that callus healing is required to prevent fatigue failure.

Cite this article: A. R. MacLeod, G. Serrancoli, B. J. Fregly, A. D. Toms, H. S. Gill. The effect of plate design, bridging span, and fracture healing on the performance of high tibial osteotomy plates: An experimental and finite element study. Bone Joint Res 2018;7:639–649. DOI: 10.1302/2046-3758.712.BJR-2018-0035.R1.

Can Measured Synergy Excitations Accurately Construct Unmeasured Muscle Excitations?

Accurate prediction of muscle and joint contact forces during human movement could improve treatment planning for disorders such as osteoarthritis, stroke, Parkinson's disease, and cerebral palsy. Recent studies suggest that muscle synergies, a low-dimensional representation of a large set of muscle electromyographic (EMG) signals (henceforth called "muscle excitations"), may reduce the redundancy of muscle excitation solutions predicted by optimization methods. This study explores the feasibility of using muscle synergy information extracted from eight muscle EMG signals (henceforth called "included" muscle excitations) to accurately construct muscle excitations from up to 16 additional EMG signals (henceforth called "excluded" muscle excitations). Using treadmill walking data collected at multiple speeds from two subjects (one healthy, one poststroke), we performed muscle synergy analysis on all possible subsets of eight included muscle excitations and evaluated how well the calculated time-varying synergy excitations could construct the remaining excluded muscle excitations (henceforth called "synergy extrapolation"). We found that some, but not all, eight-muscle subsets yielded synergy excitations that achieved >90% extrapolation variance accounted for (VAF). Using the top 10% of subsets, we developed muscle selection heuristics to identify included muscle combinations whose synergy excitations achieved high extrapolation accuracy. For 3, 4, and 5 synergies, these heuristics yielded extrapolation VAF values approximately 5% lower than corresponding reconstruction VAF values for each associated eight-muscle subset. These results suggest that synergy excitations obtained from experimentally measured muscle excitations can accurately construct unmeasured muscle excitations, which could help limit muscle excitations predicted by muscle force optimizations.

A computational framework for simultaneous estimation of muscle and joint contact forces and body motion using optimization and surrogate modeling

Concurrent estimation of muscle activations, joint contact forces, and joint kinematics by means of gradient-based optimization of musculoskeletal models is hindered by computationally expensive and non-smooth joint contact and muscle wrapping algorithms. We present a framework that simultaneously speeds up computation and removes sources of non-smoothness from muscle force optimizations using a combination of parallelization and surrogate modeling, with special emphasis on a novel method for modeling joint contact as a surrogate model of a static analysis. The approach allows one to efficiently introduce elastic joint contact models within static and dynamic optimizations of human motion. We demonstrate the approach by performing two optimizations, one static and one dynamic, using a pelvis-leg musculoskeletal model undergoing a gait cycle. We observed convergence on the order of seconds for a static optimization time frame and on the order of minutes for an entire dynamic optimization. The presented framework may facilitate model-based efforts to predict how planned surgical or rehabilitation interventions will affect post-treatment joint and muscle function.

Development of a Subject-Specific Foot-Ground Contact Model for Walking

Computational walking simulations could facilitate the development of improved treatments for clinical conditions affecting walking ability. Since an effective treatment is likely to change a patient's foot-ground contact pattern and timing, such simulations should ideally utilize deformable foot-ground contact models tailored to the patient's foot anatomy and footwear. However, no study has reported a deformable modeling approach that can reproduce all six ground reaction quantities (expressed as three reaction force components, two center of pressure (CoP) coordinates, and a free reaction moment) for an individual subject during walking. This study proposes such an approach for use in predictive optimizations of walking. To minimize complexity, we modeled each foot as two rigid segments-a hindfoot (HF) segment and a forefoot (FF) segment-connected by a pin joint representing the toes flexion-extension axis. Ground reaction forces (GRFs) and moments acting on each segment were generated by a grid of linear springs with nonlinear damping and Coulomb friction spread across the bottom of each segment. The stiffness and damping of each spring and common friction parameter values for all springs were calibrated for both feet simultaneously via a novel three-stage optimization process that used motion capture and ground reaction data collected from a single walking trial. The sequential three-stage process involved matching (1) the vertical force component, (2) all three force components, and finally (3) all six ground reaction quantities. The calibrated model was tested using four additional walking trials excluded from calibration. With only small changes in input kinematics, the calibrated model reproduced all six ground reaction quantities closely (root mean square (RMS) errors less than 13 N for all three forces, 25 mm for anterior-posterior (AP) CoP, 8 mm for medial-lateral (ML) CoP, and 2 N·m for the free moment) for both feet in all walking trials. The largest errors in AP CoP occurred at the beginning and end of stance phase when the vertical ground reaction force (vGRF) was small. Subject-specific deformable foot-ground contact models created using this approach should enable changes in foot-ground contact pattern to be predicted accurately by gait optimization studies, which may lead to improvements in personalized rehabilitation medicine.

Methodological Choices in Muscle Synergy Analysis Impact Differentiation of Physiological Characteristics Following Stroke

Muscle synergy analysis (MSA) is a mathematical technique that reduces the dimensionality of electromyographic (EMG) data. Used increasingly in biomechanics research, MSA requires methodological choices at each stage of the analysis. Differences in methodological steps affect the overall outcome, making it difficult to compare results across studies. We applied MSA to EMG data collected from individuals post-stroke identified as either responders (RES) or non-responders (nRES) on the basis of a critical post-treatment increase in walking speed. Importantly, no clinical or functional indicators identified differences between the cohort of RES and nRES at baseline. For this exploratory study, we selected the five highest RES and five lowest nRES available from a larger sample. Our goal was to assess how the methodological choices made before, during, and after MSA affect the ability to differentiate two groups with intrinsic physiologic differences based on MSA results. We investigated 30 variations in MSA methodology to determine which choices allowed differentiation of RES from nRES at baseline. Trial-to-trial variability in time-independent synergy vectors (SVs) and time-varying neural commands (NCs) were measured as a function of: (1) number of synergies computed; (2) EMG normalization method before MSA; (3) whether SVs were held constant across trials or allowed to vary during MSA; and (4) synergy analysis output normalization method after MSA. MSA methodology had a strong effect on our ability to differentiate RES from nRES at baseline. Across all 10 individuals and MSA variations, two synergies were needed to reach an average of 90% variance accounted for (VAF). Based on effect sizes, differences in SV and NC variability between groups were greatest using two synergies with SVs that varied from trial-to-trial. Differences in SV variability were clearest using unit magnitude per trial EMG normalization, while NC variability was less sensitive to EMG normalization method. No outcomes were greatly impacted by output normalization method. MSA variability for some, but not all, methods successfully differentiated intrinsic physiological differences inaccessible to traditional clinical or biomechanical assessments. Our results were sensitive to methodological choices, highlighting the need for disclosure of all aspects of MSA methodology in future studies.

Lower extremity EMG-driven modeling of walking with automated adjustment of musculoskeletal geometry

Neuromusculoskeletal disorders affecting walking ability are often difficult to manage, in part due to limited understanding of how a patient’s lower extremity muscle excitations contribute to the patient’s lower extremity joint moments. To assist in the study of these disorders, researchers have developed electromyography (EMG) driven neuromusculoskeletal models utilizing scaled generic musculoskeletal geometry. While these models can predict individual muscle contributions to lower extremity joint moments during walking, the accuracy of the predictions can be hindered by errors in the scaled geometry. This study presents a novel EMG-driven modeling method that automatically adjusts surrogate representations of the patient’s musculoskeletal geometry to improve prediction of lower extremity joint moments during walking. In addition to commonly adjusted neuromusculoskeletal model parameters, the proposed method adjusts model parameters defining muscle-tendon lengths, velocities, and moment arms. We evaluated our EMG-driven modeling method using data collected from a high-functioning hemiparetic subject walking on an instrumented treadmill at speeds ranging from 0.4 to 0.8 m/s. EMG-driven model parameter values were calibrated to match inverse dynamic moments for five degrees of freedom in each leg while keeping musculoskeletal geometry close to that of an initial scaled musculoskeletal model. We found that our EMG-driven modeling method incorporating automated adjustment of musculoskeletal geometry predicted net joint moments during walking more accurately than did the same method without geometric adjustments. Geometric adjustments improved moment prediction errors by 25% on average and up to 52%, with the largest improvements occurring at the hip. Predicted adjustments to musculoskeletal geometry were comparable to errors reported in the literature between scaled generic geometric models and measurements made from imaging data. Our results demonstrate that with appropriate experimental data, joint moment predictions for walking generated by an EMG-driven model can be improved significantly when automated adjustment of musculoskeletal geometry is included in the model calibration process.

Muscle Synergies Facilitate Computational Prediction of Subject-Specific Walking Motions

Researchers have explored a variety of neurorehabilitation approaches to restore normal walking function following a stroke. However, there is currently no objective means for prescribing and implementing treatments that are likely to maximize recovery of walking function for any particular patient. As a first step toward optimizing neurorehabilitation effectiveness, this study develops and evaluates a patient-specific synergy-controlled neuromusculoskeletal simulation framework that can predict walking motions for an individual post-stroke. The main question we addressed was whether driving a subject-specific neuromusculoskeletal model with muscle synergy controls (5 per leg) facilitates generation of accurate walking predictions compared to a model driven by muscle activation controls (35 per leg) or joint torque controls (5 per leg). To explore this question, we developed a subject-specific neuromusculoskeletal model of a single high-functioning hemiparetic subject using instrumented treadmill walking data collected at the subject's self-selected speed of 0.5 m/s. The model included subject-specific representations of lower-body kinematic structure, foot-ground contact behavior, electromyography-driven muscle force generation, and neural control limitations and remaining capabilities. Using direct collocation optimal control and the subject-specific model, we evaluated the ability of the three control approaches to predict the subject's walking kinematics and kinetics at two speeds (0.5 and 0.8 m/s) for which experimental data were available from the subject. We also evaluated whether synergy controls could predict a physically realistic gait period at one speed (1.1 m/s) for which no experimental data were available. All three control approaches predicted the subject's walking kinematics and kinetics (including ground reaction forces) well for the model calibration speed of 0.5 m/s. However, only activation and synergy controls could predict the subject's walking kinematics and kinetics well for the faster non-calibration speed of 0.8 m/s, with synergy controls predicting the new gait period the most accurately. When used to predict how the subject would walk at 1.1 m/s, synergy controls predicted a gait period close to that estimated from the linear relationship between gait speed and stride length. These findings suggest that our neuromusculoskeletal simulation framework may be able to bridge the gap between patient-specific muscle synergy information and resulting functional capabilities and limitations.

Computational Prediction of Muscle Moments During ARED Squat Exercise on the International Space Station

Prevention of muscle atrophy caused by reduced mechanical loading in microgravity conditions remains a challenge for long-duration spaceflight. To combat leg muscle atrophy, astronauts on the International Space Station (ISS) often perform squat exercise using the Advanced Resistive Exercise Device (ARED). While the ARED is effective at building muscle strength and volume on Earth, NASA researchers do not know how closely ARED squat exercise on the ISS replicates Earth-level squat muscle moments, or how small variations in exercise form affect muscle loading. This study used dynamic simulations of ARED squat exercise on the ISS to address these two questions. A multibody dynamic model of the complete astronaut-ARED system was constructed in OpenSim. With the ARED base locked to ground and gravity set to 9.81 m/s², we validated the model by reproducing muscle moments, ground reaction forces, and foot center of pressure (CoP) positions for ARED squat exercise on Earth. With the ARED base free to move relative to the ISS and gravity set to zero, we then used the validated model to simulate ARED squat exercise on the ISS for a reference squat motion and eight altered squat motions involving changes in anterior-posterior (AP) foot or CoP position on the ARED footplate. The reference squat motion closely reproduced Earth-level muscle moments for all joints except the ankle. For the altered squat motions, changing the foot position was more effective at altering muscle moments than was changing the CoP position. All CoP adjustments introduced an undesirable shear foot reaction force that could cause the feet to slip on the ARED footplate, while some foot and CoP adjustments introduced an undesirable sagittal plane foot reaction moment that would cause the astronaut to rotate off the ARED footplate without the use of some type of foot fixation. Our results provide potentially useful information for achieving desired increases or decreases in specific muscle moments during ARED squat exercise performed on the ISS.

Neuromusculoskeletal Model Calibration Significantly Affects Predicted Knee Contact Forces for Walking

Though walking impairments are prevalent in society, clinical treatments are often ineffective at restoring lost function. For this reason, researchers have begun to explore the use of patient-specific computational walking models to develop more effective treatments. However, the accuracy with which models can predict internal body forces in muscles and across joints depends on how well relevant model parameter values can be calibrated for the patient. This study investigated how knowledge of internal knee contact forces affects calibration of neuromusculoskeletal model parameter values and subsequent prediction of internal knee contact and leg muscle forces during walking. Model calibration was performed using a novel two-level optimization procedure applied to six normal walking trials from the Fourth Grand Challenge Competition to Predict In Vivo Knee Loads. The outer-level optimization adjusted time-invariant model parameter values to minimize passive muscle forces, reserve actuator moments, and model parameter value changes with (Approach A) and without (Approach B) tracking of experimental knee contact forces. Using the current guess for model parameter values but no knee contact force information, the inner-level optimization predicted time-varying muscle activations that were close to experimental muscle synergy patterns and consistent with the experimental inverse dynamic loads (both approaches). For all the six gait trials, Approach A predicted knee contact forces with high accuracy for both compartments (average correlation coefficient r = 0.99 and root mean square error (RMSE) = 52.6 N medial; average r = 0.95 and RMSE = 56.6 N lateral). In contrast, Approach B overpredicted contact force magnitude for both compartments (average RMSE = 323 N medial and 348 N lateral) and poorly matched contact force shape for the lateral compartment (average r = 0.90 medial and -0.10 lateral). Approach B had statistically higher lateral muscle forces and lateral optimal muscle fiber lengths but lower medial, central, and lateral normalized muscle fiber lengths compared to Approach A. These findings suggest that poorly calibrated model parameter values may be a major factor limiting the ability of neuromusculoskeletal models to predict knee contact and leg muscle forces accurately for walking.

Evaluation of Direct Collocation Optimal Control Problem Formulations for Solving the Muscle Redundancy Problem

Estimation of muscle forces during motion involves solving an indeterminate problem (more unknown muscle forces than joint moment constraints), frequently via optimization methods. When the dynamics of muscle activation and contraction are modeled for consistency with muscle physiology, the resulting optimization problem is dynamic and challenging to solve. This study sought to identify a robust and computationally efficient formulation for solving these dynamic optimization problems using direct collocation optimal control methods. Four problem formulations were investigated for walking based on both a two and three dimensional model. Formulations differed in the use of either an explicit or implicit representation of contraction dynamics with either muscle length or tendon force as a state variable. The implicit representations introduced additional controls defined as the time derivatives of the states, allowing the nonlinear equations describing contraction dynamics to be imposed as algebraic path constraints, simplifying their evaluation. Problem formulation affected computational speed and robustness to the initial guess. The formulation that used explicit contraction dynamics with muscle length as a state failed to converge in most cases. In contrast, the two formulations that used implicit contraction dynamics converged to an optimal solution in all cases for all initial guesses, with tendon force as a state generally being the fastest. Future work should focus on comparing the present approach to other approaches for computing muscle forces. The present approach lacks some of the major limitations of established methods such as static optimization and computed muscle control while remaining computationally efficient.

CEINMS: A toolbox to investigate the influence of different neural control solutions on the prediction of muscle excitation and joint moments during dynamic motor tasks

Personalized neuromusculoskeletal (NMS) models can represent the neurological, physiological, and anatomical characteristics of an individual and can be used to estimate the forces generated inside the human body. Currently, publicly available software to calculate muscle forces are restricted to static and dynamic optimisation methods, or limited to isometric tasks only. We have created and made freely available for the research community the Calibrated EMG-Informed NMS Modelling Toolbox (CEINMS), an OpenSim plug-in that enables investigators to predict different neural control solutions for the same musculoskeletal geometry and measured movements. CEINMS comprises EMG-driven and EMG-informed algorithms that have been previously published and tested. It operates on dynamic skeletal models possessing any number of degrees of freedom and musculotendon units and can be calibrated to the individual to predict measured joint moments and EMG patterns. In this paper we describe the components of CEINMS and its integration with OpenSim. We then analyse how EMG-driven, EMG-assisted, and static optimisation neural control solutions affect the estimated joint moments, muscle forces, and muscle excitations, including muscle co-contraction.

Residual Elimination Algorithm Enhancements to Improve Foot Motion Tracking During Forward Dynamic Simulations of Gait

Patient-specific gait optimizations capable of predicting post-treatment changes in joint motions and loads could improve treatment design for gait-related disorders. To maximize potential clinical utility, such optimizations should utilize full-body three-dimensional patient-specific musculoskeletal models, generate dynamically consistent gait motions that reproduce pretreatment marker measurements closely, and achieve accurate foot motion tracking to permit deformable foot-ground contact modeling. This study enhances an existing residual elimination algorithm (REA) Remy, C. D., and Thelen, D. G., 2009, “Optimal Estimation of Dynamically Consistent Kinematics and Kinetics for Forward Dynamic Simulation of Gait,” ASME J. Biomech. Eng., 131(3), p. 031005) to achieve all three requirements within a single gait optimization framework. We investigated four primary enhancements to the original REA: (1) manual modification of tracked marker weights, (2) automatic modification of tracked joint acceleration curves, (3) automatic modification of algorithm feedback gains, and (4) automatic calibration of model joint and inertial parameter values. We evaluated the enhanced REA using a full-body three-dimensional dynamic skeletal model and movement data collected from a subject who performed four distinct gait patterns: walking, marching, running, and bounding. When all four enhancements were implemented together, the enhanced REA achieved dynamic consistency with lower marker tracking errors for all segments, especially the feet (mean root-mean-square (RMS) errors of 3.1 versus 18.4 mm), compared to the original REA. When the enhancements were implemented separately and in combinations, the most important one was automatic modification of tracked joint acceleration curves, while the least important enhancement was automatic modification of algorithm feedback gains. The enhanced REA provides a framework for future gait optimization studies that seek to predict subject-specific post-treatment gait patterns involving large changes in foot-ground contact patterns made possible through deformable foot-ground contact models.

An Open-Source Toolbox for Surrogate Modeling of Joint Contact Mechanics

GOAL

Incorporation of elastic joint contact models into simulations of human movement could facilitate studying the interactions between muscles, ligaments, and bones. Unfortunately, elastic joint contact models are often too expensive computationally to be used within iterative simulation frameworks. This limitation can be overcome by using fast and accurate surrogate contact models that fit or interpolate input-output data sampled from existing elastic contact models. However, construction of surrogate contact models remains an arduous task. The aim of this paper is to introduce an open-source program called Surrogate Contact Modeling Toolbox (SCMT) that facilitates surrogate contact model creation, evaluation, and use.

METHODS

SCMT interacts with the third-party software FEBio to perform elastic contact analyses of finite-element models and uses MATLAB to train neural networks that fit the input-output contact data. SCMT features sample point generation for multiple domains, automated sampling, sample point filtering, and surrogate model training and testing.

RESULTS

An overview of the software is presented along with two example applications. The first example demonstrates creation of surrogate contact models of artificial tibiofemoral and patellofemoral joints and evaluates their computational speed and accuracy, while the second demonstrates the use of surrogate contact models in a forward dynamic simulation of an open-chain leg extension-flexion motion.

CONCLUSION

SCMT facilitates the creation of computationally fast and accurate surrogate contact models. Additionally, it serves as a bridge between FEBio and OpenSim musculoskeletal modeling software.

SIGNIFICANCE

Researchers may now create and deploy surrogate models of elastic joint contact with minimal effort.

Contribution of tibiofemoral joint contact to net loads at the knee in gait

Inverse dynamics analysis is commonly used to estimate the net loads at a joint during human motion. Most lower-limb models of movement represent the knee as a simple hinge joint when calculating muscle forces. This approach is limited because it neglects the contributions from tibiofemoral joint contact forces and may therefore lead to errors in estimated muscle forces. The aim of this study was to quantify the contributions of tibiofemoral joint contact loads to the net knee loads calculated from inverse dynamics for multiple subjects and multiple gait patterns. Tibiofemoral joint contact loads were measured in four subjects with instrumented implants as each subject walked at their preferred speed (normal gait) and performed prescribed gait modifications designed to treat medial knee osteoarthritis. Tibiofemoral contact loads contributed substantially to the net knee extension and knee adduction moments in normal gait with mean values of 16% and 54%, respectively. These findings suggest that knee-contact kinematics and loads should be included in lower-limb models of movement for more accurate determination of muscle forces. The results of this study may be used to guide the development of more realistic lower-limb models that account for the effects of tibiofemoral joint contact at the knee.

Patient-specific computer model of dynamic squatting after total knee arthroplasty

Knee forces are highly relevant to performance after total knee arthroplasty especially during high flexion activities such as squatting. We constructed subject-specific models of two patients implanted with instrumented knee prostheses that measured knee forces in vivo. In vivo peak forces ranged from 2.2 to 2.3 times bodyweight but peaked at different flexion angles based on the type of squatting activity. Our model predicted tibiofemoral contact force with reasonable accuracy in both subjects. This model can be a very useful tool to predict the effect of surgical techniques and component alignment on contact forces. In addition, this model could be used for implant design development, to enhance knee function, to predict forces generated during other activities, and for predicting clinical outcomes.

Muscle synergies may improve optimization prediction of knee contact forces during walking

The ability to predict patient-specific joint contact and muscle forces accurately could improve the treatment of walking-related disorders. Muscle synergy analysis, which decomposes a large number of muscle electromyographic (EMG) signals into a small number of synergy control signals, could reduce the dimensionality and thus redundancy of the muscle and contact force prediction process. This study investigated whether use of subject-specific synergy controls can improve optimization prediction of knee contact forces during walking. To generate the predictions, we performed mixed dynamic muscle force optimizations (i.e., inverse skeletal dynamics with forward muscle activation and contraction dynamics) using data collected from a subject implanted with a force-measuring knee replacement. Twelve optimization problems (three cases with four subcases each) that minimized the sum of squares of muscle excitations were formulated to investigate how synergy controls affect knee contact force predictions. The three cases were: (1) Calibrate+Match where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously matched, (2) Precalibrate+Predict where experimental knee contact forces were predicted using precalibrated muscle model parameters values from the first case, and (3) Calibrate+Predict where muscle model parameter values were calibrated and experimental knee contact forces were simultaneously predicted, all while matching inverse dynamic loads at the hip, knee, and ankle. The four subcases used either 44 independent controls or five synergy controls with and without EMG shape tracking. For the Calibrate+Match case, all four subcases closely reproduced the measured medial and lateral knee contact forces (R2 ≥ 0.94, root-mean-square (RMS) error < 66 N), indicating sufficient model fidelity for contact force prediction. For the Precalibrate+Predict and Calibrate+Predict cases, synergy controls yielded better contact force predictions (0.61 < R2 < 0.90, 83 N < RMS error < 161 N) than did independent controls (-0.15 < R2 < 0.79, 124 N < RMS error < 343 N) for corresponding subcases. For independent controls, contact force predictions improved when precalibrated model parameter values or EMG shape tracking was used. For synergy controls, contact force predictions were relatively insensitive to how model parameter values were calibrated, while EMG shape tracking made lateral (but not medial) contact force predictions worse. For the subject and optimization cost function analyzed in this study, use of subject-specific synergy controls improved the accuracy of knee contact force predictions, especially for lateral contact force when EMG shape tracking was omitted, and reduced prediction sensitivity to uncertainties in muscle model parameter values.

Neuromuscular complexity during gait is not responsive to medication in persons with Parkinson's disease

The purpose of this study was to investigate the effects of dopaminergic therapy on neuromuscular complexity during gait and on the relationship between neuromuscular complexity and gait speed in persons with Parkinson's disease (PD). Nine persons with PD walked at self-selected speed for 5 min after having withdrawn from dopaminergic medication for at least 12 h and while optimally-medicated. Electromyographic recordings were taken from eight leg muscles bilaterally. Non-negative matrix factorization was applied to reduce the dimensionality of the electromyographic signals into motor modules. We assessed neuromuscular complexity by investigating the number, structure, and timing of the modules. We also investigated the influence of dopaminergic medication on the relationships between neuromuscular complexity and gait speed. Though gait speed increased significantly after medication intake, medication did not affect neuromuscular complexity. Neuromuscular complexity was significantly associated with gait speed only while the participants were medicated. Thus, the supraspinal structures that govern neuromuscular complexity during gait do not appear to be solely dopaminergically-influenced in PD. The lack of dopaminergic influence on neuromuscular complexity may explain why persons with PD exhibit gait slowness even while medicated, and an intervention that restores neuromuscular complexity may result in gait speed improvement in PD.

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle-tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.

Are external knee load and EMG measures accurate indicators of internal knee contact forces during gait?

Mechanical loading is believed to be a critical factor in the development and treatment of knee osteoarthritis. However, the contact forces to which the knee articular surfaces are subjected during daily activities cannot be measured clinically. Thus, the ability to predict internal knee contact forces accurately using external measures (i.e., external knee loads and muscle electromyographic [EMG] signals) would be clinically valuable. We quantified how well external knee load and EMG measures predict internal knee contact forces during gait. A single subject with a force-measuring tibial prosthesis and post-operative valgus alignment performed four gait patterns (normal, medial thrust, walking pole, and trunk sway) to induce a wide range of external and internal knee joint loads. Linear regression analyses were performed to assess how much of the variability in internal contact forces was accounted for by variability in the external measures. Though the different gait patterns successfully induced significant changes in the external and internal quantities, changes in external measures were generally weak indicators of changes in total, medial, and lateral contact force. Our results suggest that when total contact force may be changing, caution should be exercised when inferring changes in knee contact forces based on observed changes in external knee load and EMG measures. Advances in musculoskeletal modeling methods may be needed for accurate estimation of in vivo knee contact forces.

Persons with Parkinson's disease exhibit decreased neuromuscular complexity during gait

OBJECTIVE

Individual muscle activation patterns may be controlled by motor modules constructed by the central nervous system to simplify motor control. This study compared modular control of gait between persons with Parkinson's disease (PD) and neurologically-healthy older adults (HOA) and investigated relationships between modular organization and gait parameters in persons with PD.

METHODS

Fifteen persons with idiopathic PD and fourteen HOA participated. Electromyographic recordings were made from eight leg muscles bilaterally while participants walked at their preferred walking speed for 10 min on an instrumented treadmill. Non-negative matrix factorization techniques decomposed the electromyographic signals, identifying the number and nature of modules accounting for 95% of variability in muscle activations during treadmill walking.

RESULTS

Generally, fewer modules were required to reconstruct muscle activation patterns during treadmill walking in PD compared to HOA (p < .05). Control of knee flexor and ankle plantar flexor musculature was simplified in PD. Activation timing was altered in PD while muscle weightings were unaffected. Simplified neuromuscular control was related to decreased walking speed in PD.

CONCLUSION

Neuromuscular control of gait is simplified in PD and may contribute to gait deficits in this population.

SIGNIFICANCE

Future studies of locomotor rehabilitation in PD should consider neuromuscular complexity to maximize intervention effectiveness.

Changes in in vivo knee contact forces through gait modification

Knee osteoarthritis (OA) commonly occurs in the medial compartment of the knee and has been linked to overloading of the medial articular cartilage. Gait modification represents a non-invasive treatment strategy for reducing medial compartment knee force. The purpose of this study was to evaluate the effectiveness of a variety of gait modifications that were expected to alter medial contact force. A single subject implanted with a force-measuring knee replacement walked using nine modified gait patterns, four of which involved different hiking pole configurations. Medial and lateral contact force at 25, 50, and 75% of stance phase, and the average value over all of stance phase (0-100%), were determined for each gait pattern. Changes in medial and lateral contact force values relative to the subject's normal gait pattern were determined by a Kruskal-Wallis test. Apart from early stance (25% of stance), medial contact force was most effectively reduced by walking with long hiking poles and wide pole placement, which significantly reduced medial and lateral contact force during stance phase by up to 34% (at 75% of stance) and 26% (at 50% of stance), respectively. Although this study is based on data from a single subject, the results provide important insight into changes in medial and lateral contact forces through gait modification. The results of this study suggest that an optimal configuration of bilateral hiking poles may significantly reduce both medial and lateral compartment knee forces in individuals with medial knee osteoarthritis.

Update on grand challenge competition to predict in vivo knee loads

Validation is critical if clinicians are to use musculoskeletal models to optimize treatment of individual patients with a variety of musculoskeletal disorders. This paper provides an update on the annual Grand Challenge Competition to Predict in Vivo Knee Loads, a unique opportunity for direct validation of knee contact forces and indirect validation of knee muscle forces predicted by musculoskeletal models. Three competitions (2010, 2011, and 2012) have been held at the annual American Society of Mechanical Engineers Summer Bioengineering Conference, and two more competitions are planned for the 2013 and 2014 conferences. Each year of the competition, a comprehensive data set collected from a single subject implanted with a force-measuring knee replacement is released. Competitors predict medial and lateral knee contact forces for two gait trials without knowledge of the experimental knee contact force measurements. Predictions are evaluated by calculating root-mean-square (RMS) errors and R(2) values relative to the experimentally measured medial and lateral contact forces. For the first three years of the competition, competitors used a variety of methods to predict knee contact and muscle forces, including static and dynamic optimization, EMG-driven models, and parametric numerical models. Overall, errors in predicted contact forces were comparable across years, with average RMS errors for the four competition winners ranging from 229 N to 312 N for medial contact force and from 238 N to 326 N for lateral contact force. Competitors generally predicted variations in medial contact force (highest R(2 )= 0.91) better than variations in lateral contact force (highest R(2 )= 0.70). Thus, significant room for improvement exists in the remaining two competitions. The entire musculoskeletal modeling community is encouraged to use the competition data and models for their own model validation efforts.

Implantable sensor technology: measuring bone and joint biomechanics of daily life in vivo

Stresses and strains are major factors influencing growth, remodeling and repair of musculoskeletal tissues. Therefore, knowledge of forces and deformation within bones and joints is critical to gain insight into the complex behavior of these tissues during development, aging, and response to injury and disease. Sensors have been used in vivo to measure strains in bone, intraarticular cartilage contact pressures, and forces in the spine, shoulder, hip, and knee. Implantable sensors have a high impact on several clinical applications, including fracture fixation, spine fixation, and joint arthroplasty. This review summarizes the developments in strain-measurement-based implantable sensor technology for musculoskeletal research.

Recent trends in assistive technology for mobility

Loss of physical mobility makes maximal participation in desired activities more difficult and in the worst case fully prevents participation. This paper surveys recent work in assistive technology to improve mobility for persons with a disability, drawing on examples observed during a tour of academic and industrial research sites in Europe. The underlying theme of this recent work is a more seamless integration of the capabilities of the user and the assistive technology. This improved integration spans diverse technologies, including powered wheelchairs, prosthetic limbs, functional electrical stimulation, and wearable exoskeletons. Improved integration is being accomplished in three ways: 1) improving the assistive technology mechanics; 2) improving the user-technology physical interface; and 3) sharing of control between the user and the technology. We provide an overview of these improvements in user-technology integration and discuss whether such improvements have the potential to be transformative for people with mobility impairments.