Peak Tibiofemoral Contact Forces Estimated Using IMU-Based Approaches Are Not Significantly Different from Motion Capture-Based Estimations in Patients with Knee Osteoarthritis.

Altered tibiofemoral contact forces represent a risk factor for osteoarthritis onset and progression, making optimization of the knee force distribution a target of treatment strategies. Musculoskeletal model-based simulations are a state-of-the-art method to estimate joint contact forces, but they typically require laboratory-based input and skilled operators. To overcome these limitations, ambulatory methods, relying on inertial measurement units, have been proposed to estimated ground reaction forces and, consequently, knee contact forces out-of-the-lab. This study proposes the use of a full inertial-capture-based musculoskeletal modelling workflow with an underlying probabilistic principal component analysis model trained on 1787 gait cycles in patients with knee osteoarthritis. As validation, five patients with knee osteoarthritis were instrumented with 17 inertial measurement units and 76 opto-reflective markers. Participants performed multiple overground walking trials while motion and inertial capture methods were synchronously recorded. Moderate to strong correlations were found for the inertial capture-based knee contact forces compared to motion capture with root mean square error between 0.15 and 0.40 of body weight. The results show that our workflow can inform and potentially assist clinical practitioners to monitor knee joint loading in physical therapy sessions and eventually assess long-term therapeutic effects in a clinical context.

In Silico Biomarkers of Motor Function to Inform Musculoskeletal Rehabilitation and Orthopedic Treatment.

In this review, we elaborate on how musculoskeletal (MSK) modeling combined with dynamic movement simulation is gradually evolving from a research tool to a promising in silico tool to assist medical doctors and physical therapists in decision making by providing parameters relating to dynamic MSK function and loading. This review primarily focuses on our own and related work to illustrate the framework and the interpretation of MSK model-based parameters in patients with 3 different conditions, that is, degenerative joint disease, cerebral palsy, and adult spinal deformities. By selecting these 3 clinical applications, we also aim to demonstrate the differing levels of clinical readiness of the different simulation frameworks introducing in silico model-based biomarkers of motor function to inform MSK rehabilitation and treatment, with the application for adult spinal deformities being the most recent of the 3. Based on these applications, barriers to clinical integration and positioning of these in silico technologies within standard clinical practice are discussed in the light of specific challenges related to model assumptions, required level of complexity and personalization, and clinical implementation.

Inertial Sensor-to-Segment Calibration for Accurate 3D Joint Angle Calculation for Use in OpenSim.

Inertial capture (InCap) systems combined with musculoskeletal (MSK) models are an attractive option for monitoring 3D joint kinematics in an ecological context. However, the primary limiting factor is the sensor-to-segment calibration, which is crucial to estimate the body segment orientations. Walking, running, and stair ascent and descent trials were measured in eleven healthy subjects with the Xsens InCap system and the Vicon 3D motion capture (MoCap) system at a self-selected speed. A novel integrated method that combines previous sensor-to-segment calibration approaches was developed for use in a MSK model with three degree of freedom (DOF) hip and knee joints. The following were compared: RMSE, range of motion (ROM), peaks, and R2 between InCap kinematics estimated with different calibration methods and gold standard MoCap kinematics. The integrated method reduced the RSME for both the hip and the knee joints below 5°, and no statistically significant differences were found between MoCap and InCap kinematics. This was consistent across all the different analyzed movements. The developed method was integrated on an MSK model workflow, and it increased the sensor-to-segment calibration accuracy for an accurate estimate of 3D joint kinematics compared to MoCap, guaranteeing a clinical easy-to-use approach.

A multi-scale modelling framework combining musculoskeletal rigid-body simulations with adaptive finite element analyses, to evaluate the impact of femoral geometry on hip joint contact forces and femoral bone growth.

Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations.

Machine learning methods to support personalized neuromusculoskeletal modelling.

Many biomedical, orthopaedic, and industrial applications are emerging that will benefit from personalized neuromusculoskeletal models. Applications include refined diagnostics, prediction of treatment trajectories for neuromusculoskeletal diseases, in silico design, development, and testing of medical implants, and human-machine interfaces to support assistive technologies. This review proposes how physics-based simulation, combined with machine learning approaches from big data, can be used to develop high-fidelity personalized representations of the human neuromusculoskeletal system. The core neuromusculoskeletal model features requiring personalization are identified, and big data/machine learning approaches for implementation are presented together with recommendations for further research.