In-vivo hip joint contact forces (HJCF) can be estimated using computational neuromusculoskeletal (NMS) modelling. However, different neural solutions can result in different HJCF estimations. NMS model predictions are also influenced by the selection of neuromuscular parameters, which are either based on cadaveric data or calibrated to the individual. To date, the best combination of neural solution and parameter calibration to obtain plausible estimations of HJCF have not been identified. The aim of this study was to determine the effect of three electromyography (EMG)-informed neural solution modes (EMG-driven, EMG-hybrid, and EMG-assisted) and static optimisation, each using three different parameter calibrations (uncalibrated, minimise joint moments error, and minimise joint moments error and peak HJCF), on the estimation of HJCF in a healthy population (n = 23) during walking. When compared to existing in-vivo data, the EMG-assisted mode and static optimisation produced the most physiologically plausible HJCF when using a NMS model calibrated to minimise joint moments error and peak HJCF. EMG-assisted mode produced first and second peaks of 3.55 times body weight (BW) and 3.97 BW during walking; static optimisation produced 3.75 BW and 4.19 BW, respectively. However, compared to static optimisation, EMG-assisted mode generated muscle excitations closer to recorded EMG signals (average across hip muscles R = 0.60 ± 0.37 versus R = 0.12 ± 0.14). Findings suggest that the EMG-assisted mode combined with minimise joint moments error and peak HJCF calibration is preferable for the estimation of HJCF and generation of realistic load distribution across muscles.

Abnormal hip joint contact forces (HJCF) are considered a primary mechanical contributor to the progression of hip osteoarthritis (OA). Compared to healthy controls, people with hip OA often present with altered muscle activation patterns and greater muscle cocontraction, both of which can influence HJCF. Neuromusculoskeletal (NMS) modelling is non-invasive approach to estimating HJCF, whereby different neural control solutions can be used to estimate muscle forces. Static optimisation, available within the popular NMS modelling software OpenSim, is a commonly used neural control solution, but may not account for an individual's unique muscle activation patterns and/or co-contraction that are often evident in pathological population. Alternatively, electromyography (EMG)-assisted neural control solutions, available within CEINMS software, have been shown to account for individual activation patterns in healthy people. Nonetheless, their application in people with hip OA, with conceivably greater levels of co-contraction, is yet to be explored. The aim of this study was to compare HJCF estimations using static optimisation (in OpenSim) and EMG-assisted (in CEINMS) neural control solutions during walking in people with hip OA. EMG-assisted neural control solution was more consistent with both EMG and joint moment data than static optimisation, and also predicted significantly higher HJCF peaks (p<0.001). The EMG-assisted neural control solution also accounted for more muscle co-contraction than static optimisation (p=0.03), which probably contributed to these higher HJCF peaks. Findings suggest that the EMG-assisted neural control solution may estimate more physiologically plausible HJCF than static optimisation in a population with high levels of cocontraction, such as hip OA.

Gait analysis together with musculoskeletal modeling is widely used for research. In the absence of medical images, surface marker locations are used to scale a generic model to the individual's anthropometry. Studies evaluating the accuracy and reliability of different scaling approaches in a pediatric and/or clinical population have not yet been conducted and, therefore, formed the aim of this study. Magnetic resonance images (MRI) and motion capture data were collected from 12 participants with cerebral palsy and 6 typically developed participants. Accuracy was assessed by comparing the scaled model's segment measures to the corresponding MRI measures, whereas reliability was assessed by comparing the model's segments scaled with the experimental marker locations from the first and second motion capture session. The inclusion of joint centers into the scaling process significantly increased the accuracy of thigh and shank segment length estimates compared to scaling with markers alone. Pelvis scaling approaches which included the pelvis depth measure led to the highest errors compared to the MRI measures. Reliability was similar between scaling approaches with mean ICC of 0.97. The pelvis should be scaled using pelvic width and height and the thigh and shank segment should be scaled using the proximal and distal joint centers.

Physics-based simulations of walking have the theoretical potential to support clinical decision-making by predicting the functional outcome of treatments in terms of walking performance. Yet before using such simulations in clinical practice, their ability to identify the main treatment targets in specific patients needs to be demonstrated. In this study, we generated predictive simulations of walking with a medical imaging based neuro-musculoskeletal model of a child with cerebral palsy presenting crouch gait. We explored the influence of altered muscle-tendon properties, reduced neuromuscular control complexity, and spasticity on gait dysfunction in terms of joint kinematics, kinetics, muscle activity, and metabolic cost of transport. We modeled altered muscle-tendon properties by personalizing Hill-type muscle-tendon parameters based on data collected during functional movements, simpler neuromuscular control by reducing the number of independent muscle synergies, and spasticity through delayed muscle activity feedback from muscle force and force rate. Our simulations revealed that, in the presence of aberrant musculoskeletal geometries, altered muscle-tendon properties rather than reduced neuromuscular control complexity and spasticity were the primary cause of the crouch gait pattern observed for this child, which is in agreement with the clinical examination. These results suggest that muscle-tendon properties should be the primary target of interventions aiming to restore an upright gait pattern for this child. This suggestion is in line with the gait analysis following muscle-tendon property and bone deformity corrections. Future work should extend this single case analysis to more patients in order to validate the ability of our physics-based simulations to capture the gait patterns of individual patients pre-and post-treatment. Such validation would open the door for identifying targeted treatment strategies with the aim of designing optimized interventions for neuro-musculoskeletal disorders.

Gait deficits in cerebral palsy (CP) are often treated with a single-event multi-level surgery (SEMLS). Selecting the treatment options (combination of bony and soft tissue corrections) for a specific patient is a complex endeavor and very often treatment outcome is not satisfying. A deterioration in 22.8% of the parameters describing gait performance has been reported and there is need for additional surgery in 11% of the patients. Computational simulations based on musculoskeletal models that allow clinicians to test the effects of different treatment options before surgery have the potential to drastically improve treatment outcome. However, to date, no such simulation and modeling method is available. Two important challenges are the development of methods to include patient-specific neuromechanical impairments into the models and to simulate the effect of different surgical procedures on post-operative gait performance. Therefore, we developed the SimCP framework that allows the evaluation of the effect of different simulated surgeries on gait performance of a specific patient and includes a graphical user interface (GUI) that enables performing virtual surgery on the models. We demonstrated the potential of our framework for two case studies. Models reflecting the patient-specific musculoskeletal geometry and muscle properties are generated based solely on data collected before the treatment. The patient's motor control is described based on muscle synergies derived from pre-operative EMG. The GUI is then used to modify the musculoskeletal properties according to the surgical plan. Since SEMLS does not affect motor control, the same motor control model is used to define gait performance pre- and post-operative. We use the capability gap (CG), i.e., the difference between the joint moments needed to perform healthy walking and the joint moments the personalized model can generate, to quantify gait performance. In both cases, the CG was smaller post- then pre-operative and this was in accordance with the measured change in gait kinematics after treatment.