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

ArborSim: Articulated, branching, OpenSim routing for constructing models of multi-jointed appendages with complex muscle-tendon architecture

Computational models of musculoskeletal systems are essential tools for understanding how muscles, tendons, bones, and actuation signals generate motion. In particular, the OpenSim family of models has facilitated a wide range of studies on diverse human motions, clinical studies of gait, and even non-human locomotion. However, biological structures with many joints, such as fingers, necks, tails, and spines, have been a longstanding challenge to the OpenSim modeling community, especially because these structures comprise numerous bones and are frequently actuated by extrinsic muscles that span multiple joints-often more than three-and act through a complex network of branching tendons. Existing model building software, typically optimized for limb structures, makes it difficult to build OpenSim models that accurately reflect these intricacies. Here, we introduce ArborSim, customized software that efficiently creates musculoskeletal models of highly jointed structures and can build branched muscle-tendon architectures. We used ArborSim to construct toy models of articulated structures to determine which morphological features make a structure most sensitive to branching. By comparing the joint kinematics of models constructed with branched and parallel muscle-tendon units, we found that among various parameters-the number of tendon branches, the number of joints between branches, and the ratio of muscle fiber length to muscle tendon unit length-the number of tendon branches and the number of joints between branches are most sensitive to branching modeling method. Notably, the differences between these models showed no predictable pattern with increased complexity. As the proportion of muscle increased, the kinematic differences between branched and parallel models units also increased. Our findings suggest that stress and strain interactions between distal tendon branches and proximal tendon and muscle greatly affect the overall kinematics of a musculoskeletal system. By incorporating complex muscle-tendon branching into OpenSim models using ArborSim, we can gain deeper insight into the interactions between the axial and appendicular skeleton, model the evolution and function of diverse animal tails, and understand the mechanics of more complex motions and tasks.

Rethinking running biomechanics: a critical review of ground reaction forces, tibial bone loading, and the role of wearable sensors

This study presents a comprehensive review of the correlation between tibial acceleration (TA), ground reaction forces (GRF), and tibial bone loading, emphasizing the critical role of wearable sensor technology in accurately measuring these biomechanical forces in the context of running. This systematic review and meta-analysis searched various electronic databases (PubMed, SPORTDiscus, Scopus, IEEE Xplore, and ScienceDirect) to identify relevant studies. It critically evaluates existing research on GRF and tibial acceleration (TA) as indicators of running-related injuries, revealing mixed findings. Intriguingly, recent empirical data indicate only a marginal link between GRF, TA, and tibial bone stress, thus challenging the conventional understanding in this field. The study also highlights the limitations of current biomechanical models and methodologies, proposing a paradigm shift towards more holistic and integrated approaches. The study underscores wearable sensors' potential, enhanced by machine learning, in transforming the monitoring, prevention, and rehabilitation of running-related injuries.

Integrating an LSTM framework for predicting ankle joint biomechanics during gait using inertial sensors

The ankle joint plays a crucial role in gait, facilitating the articulation of the lower limb, maintaining foot-ground contact, balancing the body, and transmitting the center of gravity. This study aimed to implement long short-term memory (LSTM) networks for predicting ankle joint angles, torques, and contact forces using inertial measurement unit (IMU) sensors. Twenty-five healthy participants were recruited. Two IMU sensors were attached to the foot dorsum and the vertical axis of the distal anteromedial tibia in the right lower limb to record acceleration and angular velocity during running. We proposed a LSTM-MLP (multilayer perceptron) model for training time-series data from IMU sensors and predicting ankle joint biomechanics. The model underwent validation and testing using a custom nested k-fold cross-validation process. The average values of the coefficient of determination (R2), mean absolute error (MAE), and mean squared error (MSE) for ankle dorsiflexion joint and moment, subtalar inversion joint and moment, and ankle joint contact forces were 0.89 ± 0.04, 0.75 ± 1.04, and 2.96 ± 4.96 for walking, and 0.87 ± 0.07, 0.88 ± 1.26, and 4.1 ± 7.17 for running, respectively. This study demonstrates that IMU sensors, combined with LSTM neural networks, are invaluable tools for evaluating ankle joint biomechanics in lower limb pathological diagnosis and rehabilitation, offering a cost-effective and versatile alternative to traditional experimental settings.

A framework based on subject-specific musculoskeletal models and Monte Carlo simulations to personalize muscle coordination retraining

Excessive loads at lower limb joints can lead to pain and degenerative diseases. Altering joint loads with muscle coordination retraining might help to treat or prevent clinical symptoms in a non-invasive way. Knowing how much muscle coordination retraining can reduce joint loads and which muscles have the biggest impact on joint loads is crucial for personalized gait retraining. We introduced a simulation framework to quantify the potential of muscle coordination retraining to reduce joint loads for an individuum. Furthermore, the proposed framework enables to pinpoint muscles, which alterations have the highest likelihood to reduce joint loads. Simulations were performed based on three-dimensional motion capture data of five healthy adolescents (femoral torsion 10°-29°, tibial torsion 19°-38°) and five patients with idiopathic torsional deformities at the femur and/or tibia (femoral torsion 18°-52°, tibial torsion 3°-50°). For each participant, a musculoskeletal model was modified to match the femoral and tibial geometry obtained from magnetic resonance images. Each participant's model and the corresponding motion capture data were used as input for a Monte Carlo analysis to investigate how different muscle coordination strategies influence joint loads. OpenSim was used to run 10,000 simulations for each participant. Root-mean-square of muscle forces and peak joint contact forces were compared between simulations. Depending on the participant, altering muscle coordination led to a maximum reduction in hip, knee, patellofemoral and ankle joint loads between 5 and 18%, 4% and 45%, 16% and 36%, and 2% and 6%, respectively. In some but not all participants reducing joint loads at one joint increased joint loads at other joints. The required alteration in muscle forces to achieve a reduction in joint loads showed a large variability between participants. The potential of muscle coordination retraining to reduce joint loads depends on the person's musculoskeletal geometry and gait pattern and therefore showed a large variability between participants, which highlights the usefulness and importance of the proposed framework to personalize gait retraining.

NSF DARE-transforming modeling in neurorehabilitation: perspectives and opportunities from US funding agencies

In recognition of the importance and timeliness of computational models for accelerating progress in neurorehabilitation, the U.S. National Science Foundation (NSF) and the National Institutes of Health (NIH) sponsored a conference in March 2023 at the University of Southern California that drew global participation from engineers, scientists, clinicians, and trainees. This commentary highlights promising applications of computational models to understand neurorehabilitation ("Using computational models to understand complex mechanisms in neurorehabilitation" section), improve rehabilitation care in the context of digital twin frameworks ("Using computational models to improve delivery and implementation of rehabilitation care" section), and empower future interdisciplinary workforces to deliver higher-quality clinical care using computational models ("Using computational models in neurorehabilitation requires an interdisciplinary workforce" section). The authors describe near-term gaps and opportunities, all of which encourage interdisciplinary team science. Four major opportunities were identified including (1) deciphering the relationship between engineering figures of merit-a term commonly used by engineers to objectively quantify the performance of a device, system, method, or material relative to existing state of the art-and clinical outcome measures, (2) validating computational models from engineering and patient perspectives, (3) creating and curating datasets that are made publicly accessible, and (4) developing new transdisciplinary frameworks, theories, and models that incorporate the complexities of the nervous and musculoskeletal systems. This commentary summarizes U.S. funding opportunities by two Federal agencies that support computational research in neurorehabilitation. The NSF has funding programs that support high-risk/high-reward research proposals on computational methods in neurorehabilitation informed by theory- and data-driven approaches. The NIH supports the development of new interventions and therapies for a wide range of nervous system injuries and impairments informed by the field of computational modeling. The conference materials can be found at https://dare2023.usc.edu/ .

On the prediction of tibiofemoral contact forces for healthy individuals and osteoarthritis patients during gait: a comparative study of regression methods

Knee osteoarthritis (OA) is a public health problem affecting millions of people worldwide. The intensity of the tibiofemoral contact forces is related to cartilage degeneration, and so is the importance of quantifying joint loads during daily activities. Although simulation with musculoskeletal models has been used to calculate joint loads, it demands high-cost equipment and a very time-consuming process. This study aimed to evaluate consolidated machine learning algorithms to predict tibiofemoral forces during gait analysis of healthy individuals and knee OA patients. Also, we evaluated three different datasets to train each model, considering different combinations of primary kinematic and kinetic data, and post-processing data. We evaluated 14 patients with severe unilateral knee OA and 14 healthy individuals during 3-5 gait trials. Data were split into 70% and 30% of the samples as training and test data. Test data was independently evaluated considering a mixture of pathological and healthy individuals, and only OA and Control patients. The main results showed that accurate predictions of the tibiofemoral contact forces were achieved using machine learning methods and that the predictions were sensitive to changes in the input data as training. The present study provided insights into the most promising regressions methods to predict knee contact forces representing an important starting point for the broader application of biomechanical analysis in clinical environments.