Grand challenge competition to predict in vivo knee loads

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 Use of Wearable Inertial Sensors and Workplace-Based Exercises to Reduce Lateral Epicondylitis in the Workstation of a Textile Logistics Center

People whose jobs involve repetitive motions of the wrist and forearm can suffer from lateral epicondylitis, which is a significant burden on both the individual and the employer due to treatment costs, reduced productivity, and work absenteeism. This paper describes an ergonomic intervention to reduce lateral epicondylitis in the workstation of a textile logistics center. The intervention includes workplace-based exercise programs, evaluation of risk factors, and movement correction. An injury- and subject-specific score was calculated from the motion captured with wearable inertial sensors at the workplace to evaluate the risk factors of 93 workers. Then, a new working movement was adapted to the workplace, which limited the observed risk factors and took into account the subject-specific physical abilities. The movement was taught to the workers during personalized sessions. The risk factors of 27 workers were evaluated again after the intervention to validate the effectiveness of the movement correction. In addition, active warm-up and stretching programs were introduced as part of the workday to promote muscle endurance and improve resistance to repetitive stress. The present strategy offered good results at low cost, without any physical modification of the workplace and without any detriment to productivity.

Subject-specific tribo-contact conditions in total knee replacements: a simulation framework across scales

Fundamental knowledge about in vivo kinematics and contact conditions at the articulating interfaces of total knee replacements are essential for predicting and optimizing their behavior and durability. However, the prevailing motions and contact stresses in total knee replacements cannot be precisely determined using conventional in vivo measurement methods. In silico modeling, in turn, allows for a prediction of the loads, velocities, deformations, stress, and lubrication conditions across the scales during gait. Within the scope of this paper, we therefore combine musculoskeletal modeling with tribo-contact modeling. In the first step, we compute contact forces and sliding velocities by means of inverse dynamics approach and force-dependent kinematic solver based upon experimental gait data, revealing contact forces during healthy/physiological gait of young subjects. In a second step, the derived data are employed as input data for an elastohydrodynamic model based upon the finite element method full-system approach taking into account elastic deformation, the synovial fluid's hydrodynamics as well as mixed lubrication to predict and discuss the subject-specific pressure and lubrication conditions.

Overcoming Challenges and Innovations in Orthopedic Prosthesis Design: An Interdisciplinary Perspective

Recent advances in the orthopedic prostheses design have significantly improved the quality of life for individuals with orthopedic disabilities. However, there are still critical challenges that need to be addressed to further enhance the functionality of orthopedic prostheses improving biocompatibility to promote better integration with natural tissues, enhancing durability to withstand the demands of daily use, and improving sensory feedback for better control of movement are the most pressing issues. To address these challenges, promising emerging solutions such as smart prosthetics, 3D printing, regenerative medicine, and artificial intelligence have been developed. These innovative technologies hold the potential to significantly enhance the functionality of orthopedic prostheses. Realizing the full potential of these next-generation orthopedic prostheses requires addressing several critical factors. These include interdisciplinary collaboration between experts in orthopedics, materials science, biology, and engineering, increased investment in research and development, standardization of components to ensure quality and reliability, and improved access to prosthetics. A comprehensive review of these challenges and considerations for future orthopedic prosthesis design is s provided in this paper addressing the further advances to the field. By addressing these issues, we can continue to improve the lives of individuals with orthopedic disabilities and further enhance the field of orthopedic prosthetics.

Conformity design can change the effect of tibial component malrotation on knee biomechanics after total knee arthroplasty

BACKGROUND

Component alignment is essential to improve knee function and survival in total knee arthroplasty. However, it is still unclear whether the conformity design of tibiofemoral component can mitigate abnormal knee biomechanics caused by component malrotation. The purpose of this study was to investigate whether the sagittal/coronal conformity design of the tibial component could change the effect of the tibial component malrotation on knee biomechanics in total knee arthroplasty.

METHODS

A developed patient-specific musculoskeletal multi-body dynamics model of total knee arthroplasty was used to investigate the effects of the sagittal/coronal conformity of the tibial component on knee contact forces and kinematics caused by tibial component malrotation during the walking gait.

FINDINGS

Medial and lateral contact forces, internal-external rotation, and anterior-posterior translation were significantly affected by tibial component malrotation after total knee arthroplasty during the walking gait. The lower sagittal conformity of the tibial component can mitigate the abnormal internal-external rotation caused by tibial component malrotation in total knee arthroplasty, the higher coronal conformity of the tibial component can mitigate the abnormal medial-lateral translation caused by tibial component malrotation in total knee arthroplasty.

INTERPRETATION

This study highlights the importance of the tibiofemoral conformity designs on knee biomechanics caused by component malrotation in total knee arthroplasty. The optimization of the tibiofemoral conformity designs should be thoroughly considered in the design of new implants and in the planning of surgical procedures.

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.

The influence of induced gait asymmetry on joint reaction forces

Chronic injury- or disease-induced joint impairments result in asymmetric gait deviations that may precipitate changes in joint loading associated with pain and osteoarthritis. Understanding the impact of gait deviations on joint reaction forces (JRFs) is challenging because of concurrent neurological and/or anatomical changes and because measuring JRFs requires medically invasive instrumented implants. Instead, we investigated the impact of joint motion limitations and induced asymmetry on JRFs by simulating data recorded as 8 unimpaired participants walked with bracing to unilaterally and bilaterally restrict ankle, knee, and simultaneous ankle + knee motion. Personalized models, calculated kinematics, and ground reaction forces (GRFs) were input into a computed muscle control tool to determine lower limb JRFs and simulated muscle activations guided by electromyography-driven timing constraints. Unilateral knee restriction increased GRF peak and loading rate ipsilaterally but peak values decreased contralaterally when compared to walking without joint restriction. GRF peak and loading rate increased with bilateral restriction compared to the contralateral limb of unilaterally restricted conditions. Despite changes in GRFs, JRFs were relatively unchanged due to reduced muscle forces during loading response. Thus, while joint restriction results in increased limb loading, reductions in muscle forces counteract changes in limb loading such that JRFs were relatively unchanged.

Can static optimization detect changes in peak medial knee contact forces induced by gait modifications?

Medial knee contact force (MCF) is related to the pathomechanics of medial knee osteoarthritis. However, MCF cannot be directly measured in the native knee, making it difficult for therapeutic gait modifications to target this metric. Static optimization, a musculoskeletal simulation technique, can estimate MCF, but there has been little work validating its ability to detect changes in MCF induced by gait modifications. In this study, we quantified the error in MCF estimates from static optimization compared to measurements from instrumented knee replacements during normal walking and seven different gait modifications. We then identified minimum magnitudes of simulated MCF changes for which static optimization correctly identified the direction of change (i.e., whether MCF increased or decreased) at least 70% of the time. A full-body musculoskeletal model with a multi-compartment knee and static optimization was used to estimate MCF. Simulations were evaluated using experimental data from three subjects with instrumented knee replacements who walked with various gait modifications for a total of 115 steps. Static optimization underpredicted the first peak (mean absolute error = 0.16 bodyweights) and overpredicted the second peak (mean absolute error = 0.31 bodyweights) of MCF. Average root mean square error in MCF over stance phase was 0.32 bodyweights. Static optimization detected the direction of change with at least 70% accuracy for early-stance reductions, late-stance reductions, and early-stance increases in peak MCF of at least 0.10 bodyweights. These results suggest that a static optimization approach accurately detects the direction of change in early-stance medial knee loading, potentially making it a valuable tool for evaluating the biomechanical efficacy of gait modifications for knee osteoarthritis.

A novel method to accurately recreate in vivo loads and kinematics in computational models of the knee

Despite availability of in vivo knee loads and kinematics data, conventional load- and displacement-controlled configurations still can't accurately predict tibiofemoral loads from kinematics or vice versa. We propose a combined load- and displacement-control method for joint-level simulations of the knee to reliably reproduce in vivo contact mechanics. Prediction errors of the new approach were compared to those of conventional purely load- or displacement-controlled models using in vivo implant loads and kinematics for multiple subjects and activities (CAMS-Knee dataset). Our method reproduced both loads and kinematics more closely than conventional models and thus demonstrates clear advantages for investigating tibiofemoral contact or wear.

Machine Learning for Optical Motion Capture-Driven Musculoskeletal Modelling from Inertial Motion Capture Data

Marker-based Optical Motion Capture (OMC) systems and associated musculoskeletal (MSK) modelling predictions offer non-invasively obtainable insights into muscle and joint loading at an in vivo level, aiding clinical decision-making. However, an OMC system is lab-based, expensive, and requires a line of sight. Inertial Motion Capture (IMC) techniques are widely-used alternatives, which are portable, user-friendly, and relatively low-cost, although with lesser accuracy. Irrespective of the choice of motion capture technique, one typically uses an MSK model to obtain the kinematic and kinetic outputs, which is a computationally expensive tool increasingly well approximated by machine learning (ML) methods. Here, an ML approach is presented that maps experimentally recorded IMC input data to the human upper-extremity MSK model outputs computed from ('gold standard') OMC input data. Essentially, this proof-of-concept study aims to predict higher-quality MSK outputs from the much easier-to-obtain IMC data. We use OMC and IMC data simultaneously collected for the same subjects to train different ML architectures that predict OMC-driven MSK outputs from IMC measurements. In particular, we employed various neural network (NN) architectures, such as Feed-Forward Neural Networks (FFNNs) and Recurrent Neural Networks (RNNs) (vanilla, Long Short-Term Memory, and Gated Recurrent Unit) and a comprehensive search for the best-fit model in the hyperparameters space in both subject-exposed (SE) as well as subject-naive (SN) settings. We observed a comparable performance for both FFNN and RNN models, which have a high degree of agreement (ravg,SE,FFNN=0.90±0.19, ravg,SE,RNN=0.89±0.17, ravg,SN,FFNN=0.84±0.23, and ravg,SN,RNN=0.78±0.23) with the desired OMC-driven MSK estimates for held-out test data. The findings demonstrate that mapping IMC inputs to OMC-driven MSK outputs using ML models could be instrumental in transitioning MSK modelling from 'lab to field'.

Patella height influences patellofemoral contact and kinematics following cruciate-retaining total knee replacement

The role of patella height is discussed controversially in total knee arthroplasty (TKA). Therefore, this computational study aims to systematically analyze the biomechanical effect of different patella heights on patellofemoral (PF) forces and kinematics after cruciate-retaining (CR) TKA. We implemented a CR bicondylar TKA with a dome patellar button in a validated dynamic musculoskeletal multibody model of a male human knee joint. Retropatellar dynamics (contact force [N], shear force [N], patellar shift [mm], tilt [°], and rotation [°]) were evaluated during dual-limb squat motion (flexion from 0° to 90°) with simulated active muscle forces and the effects of different patella heights (Blackburne-Peel [BP] ratio of 0.39, 0.49, 0.65, 0.85, 1.01, and 1.1 were systematically examined). As active knee flexion increased, PF contact force also increased. Patella alta (BP = 1.1) resulted in higher PF contact forces compared to normal patella height (BP = 0.65) by up to 16%. Contrarily, patella baja was associated with decreased PF forces by 7%. Compared to patella baja (BP = 0.39), patella alta (BP = 1.1) considerably increased the contact force by up to 25%. Different patellar heights mainly affected PF shear forces during early knee flexion. Concerning PF kinematics, patella alta (BP = 1.1) yielded a greater lateral tilt of more than 4° and higher patellar rotation by up to 3° during deep knee flexion, compared to normal patella height (BP = 0.65). Our computational study indicates that patella alta is associated with the highest PF contact and shear force after the implantation of a CR bicondylar TKA. This should be considered in PF disorders following TKA.

Knee joint underloading does not evolve after a two-week reintroduction to running program after anterior cruciate ligament reconstruction

OBJECTIVES

Knee underloading patterns have been reported mid- and long-term after return to running post-ACLR, but changes in these patterns during the reintroduction to running are unknown. We evaluated knee biomechanics in individuals within 6 months of ACL-R at the start and completion of a reintroduction to running program.

DESIGN

Longitudinal laboratory study.

SETTING

Three-dimensional running biomechanics during instrumented treadmill running.

PARTICIPANTS

24 participants post-ACL-R with hamstring autograft and 24 healthy, matched controls.

MAIN OUTCOME MEASURES

Tibiofemoral joint (TFJ) and patellofemoral joint (PFJ) contact forces, peak knee extension moment and peak knee flexion angle.

RESULTS

Significant LIMB∗GROUP interactions (all p < 0.05) but no TIME effects were found. PFJ and TFJ contact forces, peak knee flexion angle and peak knee extensor moment were lower (all p < 0.001) on the injured-limb compared to both contralateral-limb and CONTROL. PFJ and TFJ contact forces and peak knee flexion, knee extension moment were greater (all p < 0.01) on the contralateral-limb of ACL-R compared to CONTROL. There was no change in knee biomechanics after two weeks of the reintroduction to running.

CONCLUSIONS

Clinicians should be aware that substantial and persistent knee underloading does not resolve upon reintroduction to running after ACL-R.

LEVEL OF EVIDENCE

Longitudinal observational study, level III.

Numerical instability of Hill-type muscle models

Hill-type muscle models are highly preferred as phenomenological models for musculoskeletal simulation studies despite their introduction almost a century ago. The use of simple Hill-type models in simulations, instead of more recent cross-bridge models, is well justified since computationally 'light-weight'-although less accurate-Hill-type models have great value for large-scale simulations. However, this article aims to invite discussion on numerical instability issues of Hill-type muscle models in simulation studies, which can lead to computational failures and, therefore, cannot be simply dismissed as an inevitable but acceptable consequence of simplification. We will first revisit the basic premises and assumptions on the force-length and force-velocity relationships that Hill-type models are based upon, and their often overlooked but major theoretical limitations. We will then use several simple conceptual simulation studies to discuss how these numerical instability issues can manifest as practical computational problems. Lastly, we will review how such numerical instability issues are dealt with, mostly in an ad hoc fashion, in two main areas of application: musculoskeletal biomechanics and computer animation.

Tibiofemoral compressive force during downhill walking in patients with primary total knee arthroplasty: A statistical parametric mapping approach

BACKGROUND

Downhill walking is a necessary part of daily life and an effective activity in post-operative rehabilitation following total knee arthroplasty. The purpose of this study was to determine differences in the behavior of total, medial, and lateral tibiofemoral compressive forces as well as knee extensor and flexor muscle forces between different limbs of patients with total knee arthroplasty (replaced, non-replaced) during downhill and level walking.

METHODS

Musculoskeletal modeling and simulation were implemented to determine muscle forces and tibiofemoral compressive forces in 25 patients with total knee arthroplasty. A 2 × 2 [Limb (replaced, non-replaced) × Slope (0°, 10°)] Statistical parametric mapping repeated measures analysis of variance was conducted on selected variables.

FINDINGS

Statistical parametric mapping did not identify any between-limb differences for compressive or muscle forces. Differences in joint compressive and muscle forces persisted throughout different intervals of stance-phase between level and downhill walking. Knee extensor muscle forces were distinctly greater during level walking for nearly all of stance phase. Knee flexor muscle force was greater during downhill walking for >60% of stance. Statistical parametric mapping did identify regions of significance between level and downhill walking that coincided temporally (near loading response and push off) with peak joint moment and joint compressive forces traditionally reported using discrete variable analyses.

INTERPRETATION

Downhill walking may be a safe and useful rehabilitation tool for post-knee arthroplasty rehabilitation that will not disproportionally load either the replaced or the non-replaced joint and where the quadriceps muscles can be strengthened during a gait-specific task.

Computational Modelling for Managing Pathways to Cartilage Failure

Over several decades the perception and therefore description of articular cartilage changed substantially. It has transitioned from being described as a relatively inert tissue with limited repair capacity, to a tissue undergoing continuous maintenance and even adaption, through a range of complex regulatory processes. Even from the narrower lens of biomechanics, the engagement with articular cartilage has changed from it being an interesting, slippery material found in the hostile mechanical environment between opposing long bones, to an intriguing example of mechanobiology in action. The progress revealing this complexity, where physics, chemistry, material science and biology are merging, has been described with increasingly sophisticated computational models. Here we describe how these computational models of cartilage as an integrated system can be combined with the approach of structural reliability analysis. That is, causal, deterministic models placed in the framework of the probabilistic approach of structural reliability analysis could be used to understand, predict, and mitigate the risk of cartilage failure or pathology. At the heart of this approach is seeing cartilage overuse and disease processes as a 'material failure', resulting in failure to perform its function, which is largely mechanical. One can then describe pathways to failure, for example, how homeostatic repair processes can be overwhelmed leading to a compromised tissue. To illustrate this 'pathways to failure' approach, we use the interplay between cartilage consolidation and lubrication to analyse the increase in expected wear rates associated with cartilage defects or meniscectomy.

The effect of subtalar joint axis location on muscle moment arms

Most dynamic musculoskeletal models define the subtalar joint (STJ) as a one degree of freedom (DOF) hinge with a tri-planar axis. The orientation of this axis of rotation is often determined as a combination of inclination and deviation angles measured from the ground and midline of the foot, respectively. In defining the location of the axis, often the origin is found at the distal aspect of the heel instead of at the articulation of the talus and calcaneus. Key musculoskeletal modeling definitions, such as muscle moment arms, are dependent on the distance and relative location of muscle insertion to the axis of rotation. Since the axis orientation and origin location affect calculations of muscle moment arm and joint dynamics, there is much need for accurate characterization of the STJ axis to understand the STJ's role in dynamic weight-bearing motion. The purpose of this study is to explore how the STJ origin location and axis orientation affect muscle moment arms surrounding the ankle. Datasets from the Grand Knee Challenge, posted on the open-source SimTK website, were modeled using OpenSim. Modifying the location of the STJ axis from the original location closer to the articulation between the talus and calcaneus resulted in significant differences in STJ muscle moment arms and peak STJ moments. The findings of this study conclude that the location of the STJ axis origin needs to be considered and accurately defined, especially if the inclination/deviation angles of the rotational axis will be modified to represent a more subject-specific definition.

Computer-based analysis of different component positions and insert thicknesses on tibio-femoral and patello-femoral joint dynamics after cruciate-retaining total knee replacement

BACKGROUND

Positioning of the implant components and tibial insert thickness constitute critical aspects of total knee replacement (TKR) that influence the postoperative knee joint dynamics. This study aimed to investigate the impact of implant component positioning (anterior-posterior and medio-lateral shift) and varying tibial insert thickness on the tibio-femoral (TF) and patello-femoral (PF) joint kinematics and contact forces after cruciate-retaining (CR)-TKR.

METHOD

A validated musculoskeletal multibody simulation (MMBS) model with a fixed-bearing CR-TKR during a squat motion up to 90° knee flexion was deployed to calculate PF and TF joint dynamics for varied implant component positions and tibial insert thicknesses. Evaluation was performed consecutively by comparing the respective knee joint parameters (e.g. contact force, quadriceps muscle force, joint kinematics) to a reference implant position.

RESULTS

The PF contact forces were mostly affected by the anterior-posterior as well as medio-lateral positioning of the femoral component (by 3 mm anterior up to 31 % and by 6 mm lateral up to 14 %). TF contact forces were considerably altered by tibial insert thickness (24 % in case of + 4 mm increase) and by the anterior-posterior position of the femoral component (by 3 mm posterior up to 16 %). Concerning PF kinematics, a medialised femoral component by 6 mm increased the lateral patellar tilt by more than 5°.

CONCLUSIONS

Our results indicate that regarding PF kinematics and contact forces the positioning of the femoral component was more critical than the tibial component. The positioning of the femoral component in anterior-posterior direction on and PF contact force was evident. Orthopaedic surgeons should strictly monitor the anterior-posterior as well as the medio-lateral position of the femoral component and the insert thickness.

On the estimation of hip joint loads through musculoskeletal modeling

Noninvasive estimation of joint loads is still an open challenge in biomechanics. Although musculoskeletal modeling represents a solid resource, multiple improvements are still necessary to obtain accurate predictions of joint loads and to translate such potential into practical utility. The present study, focused on the hip joint, is aimed at reviewing the state-of-the-art literature on the estimation of hip joint reaction forces through musculoskeletal modeling. Our literature inspection, based on well-defined selection criteria, returned seventeen works, which were compared in terms of methods and results. Deviations between predicted and in vivo measured hip joint loads, taken from the OrthoLoad database, were assessed through quantitative deviation indices. Despite the numerous modeling and computational improvements made over the last two decades, predicted hip joint loads still deviate from their experimental counterparts and typically overestimate them. Several critical aspects have emerged that affect muscle force estimation, hence joint loads. Among them, the physical fidelity of the musculoskeletal model, with its parameters and geometry, plays a crucial role. Also, predicted joint loads are markedly affected by the selected muscle recruitment strategy, which reflects the underlying motor control policy. Practical guidelines for researchers interested in noninvasive estimation of hip joint loads are also provided.

Estimation and Comparison of Knee Joint Contact Forces During Heel Contact and Heel Rise Deep Squatting

Background

Increased knee flexion is required for deep squatting in the daily life of the non-western population as well as in many sports activities. The purpose of this study was to estimate as well as to compare knee joint contact forces during heel contact (HC) and heel rise (HR) deep squatting in 10 healthy young Indian participants.

Materials and Methods

Kinematic data were captured using a 12-camera Motion Analysis system. Kinetic data were collected using two Kistler force plates. EMG of 6 lower limb muscles was monitored by Noraxon wireless EMG. OpenSim musculoskeletal model was customized to increase the maximum knee flexion capability of the existing model and knee joint contact forces were estimated.

Results

A significant difference in tibiofemoral (p < 0.001) as well as patellofemoral (p = 0.006) knee joint contact force was observed between HC and HR squatting. The resultant maximum tibiofemoral KJCF was 5.9 (± 0.54) times body weight (BW) and 5.3 (± 0.6) BW for the HC and HR, respectively. The resultant maximum patellofemoral KJCF was 7.8 (± 0.57) BW and 7.1 (± 0.73) BW for the HC and HR, respectively.

Conclusion

The findings can provide implications for physiotherapists to design rehabilitation exercise protocols, exercise professionals, and the development of high flexion knee implants.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s43465-022-00798-y.

Prediction of in vivo hip contact forces during common activities of daily living using a segment-based musculoskeletal model

Background: Quantifying in vivo hip muscle and contact forces during activities of daily living (ADL) provides valuable information for diagnosis and treatment of hip-related disorders. The objective of this study was to utilize Freebody, a segment-based musculoskeletal model, for the prediction of hip contact forces using a novel objective function during seven common ADLs and validate its performance against the publicly available HIP98 dataset. Methods: Marker data, ground reaction forces, and hip contact forces during slow, normal, and fast walking, stair ascent and descent, and standing up and sitting down were extracted for 3 subjects from the HIP98 dataset. A musculoskeletal anatomical dataset was scaled to match the dimensions of each subject, and muscle and hip contact forces were estimated by minimizing a novel objective function, which was the summation of the muscle stresses squared and body weight-normalised hip contact force. The accuracy of predictions were quantified using several metrics, and muscle forces were qualitatively compared to experimental EMGs in the literature. Results: FreeBody predicted the hip contact forces during the ADLs with encouraging accuracy: The root mean squared error of predictions were 44.0 ± 8.5, 47.4 ± 6.5, and 59.8 ± 7.1% BW during slow, normal, and fast walking, 44.2 ± 16.8% and 53.3 ± 12.2% BW for stair ascent and descent, and 31.8 ± 8.2% and 17.1 ± 5.0% BW for standing up and sitting down, respectively. The error in prediction of peak hip contact forces were 14-18%, 24-28%, 17-35% for slow, normal, and fast walking, 7-25% and 15-32% in stair ascent and descent, and around 10% for standing up and sitting down. The coefficient of determination was larger than 0.90 in all activities except in standing up (0.86 ± 0.08). Conclusion: This study has implemented a novel objective function in a segment-based musculoskeletal model, FreeBody, for the prediction of hip contact forces during a large range of ADLs. The model outputs compare favourably for all ADLs and are the best in standing up and sitting down, while muscle activation patterns are consistent with experimental EMGs from literature. This new objective function addresses one of the major limitations associated with musculoskeletal models in the literature, namely the high non-physiological predicted hip joint contact forces.

Contributions of individual muscle forces to hip, knee, and ankle contact forces during the stance phase of running: a model-based study

Knowledge of muscle forces' contributions to the joint contact forces can assist in the evaluation of muscle function, joint injury prevention, treatment of gait disorders, and arthroplasty planning. This study's objective was to evaluate the contributions of human lower limb muscles to the hip, knee, and ankle joint contact forces during the stance phase of running. A total of 25 muscles (or groups) were investigated based on the OpenSim framework along the anterior-posterior, superoinferior, and mediolateral components of each joint coordinate system. It was revealed that, during the running stance phase, the gluteus medius, gluteus maximus, and iliopsoas mainly contributed to the hip contact force. The soleus, vastus group, and rectus femoris primarily contributed to the knee contact force, while the peroneus, soleus, gluteus medius, and gastrocnemius mainly contributed to the ankle joint force; some muscles simultaneously offloaded the joints during the stance phase. The distributive pattern of the individual muscle functions contributing to the joint load may substantially differ during the running and walking stance phases. This study's findings may further provide suggestive information for the design of lower limb joint prosthesis, the study of the biomechanics of pathologic walking and running, and the progression of joint osteoarthritis.

Patellar ligament desmopathy in the horse – a review and comparison to human patellar tendinopathy (‘Jumper’s knee’)

Patellar ligament desmopathy in horses is regarded as an uncommon condition with unclear aetiology. Of the three patellar ligaments in the horse, the intermediate is the one most often diagnosed with desmopathy in horses presenting with chronic lameness. This structure corresponds to the patellar tendon in humans. As diagnostic imaging modalities continuously improve, changes in echogenicity of the patellar ligaments are identified ultrasonographically with increasing frequency. However, disruption of the normal fibre pattern may be present also in patellar ligaments in horses that show no signs of lameness. Similarly, there is a poor correlation between pain and diagnostic imaging findings in human patellar tendinopathy. Consequently, there appears to be a knowledge gap pertaining to normal ultrasonographic variation and diagnostic criteria for disease of the patellar ligaments in horses. Furthermore, local anaesthetic techniques to verify the diagnosis are poorly described, and due to the low number of treated cases, no specific treatment modality can be recommended on a scientific basis. The aim of this paper is to review the current knowledge regarding the pathogenesis, diagnosis and management of patellar ligament desmopathy in horses, compare this condition with patellar tendinopathy in humans, and identify areas for further research to increase the diagnostic accuracy in horses. We conclude that there is a profound need for better descriptions of ultrasonographic variation and pathological changes of the equine patellar ligaments. Identification of areas of maximal ligament strain and descriptions of early histopathological changes could render more information on the possible aetiology, preventive measurements and treatment options of desmopathy. Description of regional innervation could aid in development of methods for diagnostic anaesthesia to verify pain originating from the ligaments.

Multi-level personalization of neuromusculoskeletal models to estimate physiologically plausible knee joint contact forces in children

Neuromusculoskeletal models are a powerful tool to investigate the internal biomechanics of an individual. However, commonly used neuromusculoskeletal models are generated via linear scaling of generic templates derived from elderly adult anatomies and poorly represent a child, let alone children with a neuromuscular disorder whose musculoskeletal structures and muscle activation patterns are profoundly altered. Model personalization can capture abnormalities and appropriately describe the underlying (altered) biomechanics of an individual. In this work, we explored the effect of six different levels of neuromusculoskeletal model personalization on estimates of muscle forces and knee joint contact forces to tease out the importance of model personalization for normal and abnormal musculoskeletal structures and muscle activation patterns. For six children, with and without cerebral palsy, generic scaled models were developed and progressively personalized by (1) tuning and calibrating musculotendon units' parameters, (2) implementing an electromyogram-assisted approach to synthesize muscle activations, and (3) replacing generic anatomies with image-based bony geometries, and physiologically and physically plausible muscle kinematics. Biomechanical simulations of gait were performed in the OpenSim and CEINMS software on ten overground walking trials per participant. A mixed-ANOVA test, with Bonferroni corrections, was conducted to compare all models' estimates. The model with the highest level of personalization produced the most physiologically plausible estimates. Model personalization is crucial to produce physiologically plausible estimates of internal biomechanical quantities. In particular, personalization of musculoskeletal anatomy and muscle activation patterns had the largest effect overall. Increased research efforts are needed to ease the creation of personalized neuromusculoskeletal models.

Joint kinematics alone can distinguish hip or knee osteoarthritis patients from asymptomatic controls with high accuracy

Osteoarthritis (OA) is one of the leading musculoskeletal disabilities worldwide, and several interventions intend to change the gait pattern in OA patients to more healthy patterns. However, an accessible way to follow up the biomechanical changes in a clinical setting is still missing. Therefore, this study aims to evaluate whether we can use biomechanical data collected from a specific activity of daily living to help distinguish hip OA patients from controls and knee OA patients from controls using features that potentially could be measured in a clinical setting. To achieve this goal, we considered three different classes of statistical models with different levels of data complexity. Class 1 is kinematics based only (clinically applicable), class 2 includes joint kinetics (semi-applicable under the condition of access to a force plate or prediction models), and class 3 uses data from advanced musculoskeletal modeling (not clinically applicable). We used a machine learning pipeline to determine which classification model was best. We found 100% classification accuracy for KneeOA-vs-Asymptomatic and 93.9% for HipOA-vs-Asymptomatic using seven features derived from the lumbar spine and hip kinematics collected during ascending stairs. These results indicate that kinematical data alone can distinguish hip or knee OA patients from asymptomatic controls. However, to enable clinical use, we need to validate if the classifier also works with sensor-based kinematical data and whether the probabilistic outcome of the logistic regression model can be used in the follow-up of patients with OA.

Integrating wearables and modelling for monitoring rehabilitation following total knee joint replacement

BACKGROUND AND OBJECTIVE

Wearable inertial devices integrated with modelling and cloud computing have been widely adopted in the sports sector, however, their use in the health and medical field has yet to be fully realised. To date, there have been no reported studies concerning the use of wearables as a surrogate tool to monitor knee joint loading during recovery following a total knee joint replacement. The objective of this study is to firstly evaluate if peak tibial acceleration from wearables during gait is a good surrogate metric for computer modelling predicted functional knee loading; and secondly evaluate if traditional clinical patient related outcomes measures are consistent with wearable predictions.

METHODS

Following ethical approval, four healthy participants were used to establish the relationship between computer modelling predicted knee joint loading and wearable measured tibial acceleration. Following this, ten patients who had total knee joint replacements were then followed during their 6-week rehabilitation. Gait analysis, wearable acceleration, computer models of knee joint loading, and patient related outcomes measures including the Oxford knee score and range of motion were recorded.

RESULTS

A linear correlation (R² of 0.7-0.97) was observed between peak tibial acceleration (from wearables) and musculoskeletal model predicted knee joint loading during gait in healthy participants first. Whilst patient related outcome measures (Oxford knee score and patient range of motion) were observed to improve consistently during rehabilitation, this was not consistent with all patient's tibial acceleration. Only those patients that exhibited increasing peak tibial acceleration over 6-weeks rehabilitation were positively correlated with the Oxford knee score (R² of 0.51 to 0.97). Wearable predicted tibial acceleration revealed three patients with a consistent knee loading, five patients with improving knee loading, and two patients with declining knee loading during recovery. Hence, 20% of patients did not present with satisfactory joint loading following total knee joint replacement and this was not detected with current patient related outcome measures.

CONCLUSIONS

The use of inertial measurement units or wearables in this study provided additional insight into patients who were not exhibiting functional improvements in joint loading, and offers clinicians an 'off-site' early warning metric to identify potential complications during recovery and provide the opportunity for early intervention. This study has important implications for improving patient outcomes, equity, and for those who live in rural regions.

Study on the Role of MicroRNA-214 in the Rehabilitation of Cartilage in Mice with Exercise-Induced Traumatic Osteoarthritis

This study aimed to explore the possible relationship between the expression of Micro RNA-214 (miR-214) and the pathogenesis and recovery in mice with post-traumatic osteoarthritis (PTOA). In this study, 40 male C57BL/6 mice were randomly divided into five groups: model control (MC) group, model (M) group, rehabilitation control (RC) group, model + rehabilitation (M + R) group, and model + convalescent (M + C) group. Four weeks of high-intensity treadmill exercise (HITE) and 4 weeks of moderate-intensity treadmill exercise (MITE) were implemented for PTOA modeling and rehabilitation, respectively. In vitro, 10% elongation mechanical strain was used for IL-1β stimulated chondrocytes. We found that compared with the MC group, there was a significant increase in the aspect of inflammation and catabolism while a dramatic fall in miR-214 expression was observed in the M group. After the 4 weeks of MITE, the level of inflammation and metabolism, as well as miR-214 expression, was partially reversed in the M + R group compared with the M + C group. The expression of miR-214 decreased dramatically after chondrocyte stimulation by IL-1β and then increased significantly after 10% strain was applied to IL-1β-treated cells. These results suggest that a suitable mechanical load can increase the expression of miR-214, and that miR-214 may play a chondroprotective effect in the development of OA.

Can genetics guide exercise prescriptions in osteoarthritis?

Osteoarthritis (OA) is the most common form of arthritis and has a multifactorial etiology. Current management for OA focuses on minimizing pain and functional loss, typically involving pharmacological, physical, psychosocial, and mind-body interventions. However, there remain challenges in determining which patients will benefit most from which interventions. Although exercise-based interventions are recommended as first-line treatments and are known to be beneficial for managing both the disease and illness of OA, the optimal exercise “prescription” is unknown, due in part to our limited understanding of the precise mechanisms underlying its action. Here we present our perspective on the potential role of genetics in guiding exercise prescription for persons with OA. We describe key publications in the areas of exercise and OA, genetics and OA, and exercise and genetics, and point to a paucity of knowledge at the intersection of exercise, genetics, and OA. We suggest there is emerging evidence to support the use of genetics and epigenetics to explain the beneficial effects of exercise for OA. We identify missing links in the existing research relating to exercise, genetics, and OA, and highlight epigenetics as a promising mechanism through which environmental exposures such as exercise may impact OA outcomes. We anticipate future studies will improve our understanding of how genetic and epigenetic factors mediate exercise-based interventions to support implementation and ultimately improve OA patient care.

A subject-specific musculoskeletal model to predict the tibiofemoral contact forces during daily living activities

Accurate prediction of tibiofemoral contact force (TFCF) during daily living activities is significant for understanding the initiation, progression, and treatment of knee osteoarthritis (KOA). However, the diversity of target activities, prediction accuracy, and computational efficiency of the current musculoskeletal simulations need to be further improved. In this study, a subject-specific musculoskeletal model considered the tibiofemoral alignment, medial-lateral contact locations, secondary tibiofemoral and all patellofemoral motions, and knee ligaments was proposed to predict the TFCFs during the five daily activities (normal walking, sit-to-stand, stand-to-sit, stair ascent, and stair descent) in OpenSim software. The standing lower-limbs-full-length radiograph, local radiograph of knee joint, motion capture data, and force plate data of eighteen subjects were acquired as the input data of the musculoskeletal model. The results showed good agreements of TFCFs between the predictions based on our proposed musculoskeletal model and the in-vivo measurements based on instrumented knee implants during the five daily activities (RMSE: 0.16 ∼ 0.31 BW, R²: 0.88 ∼ 0.97, M: -0.11 ∼ -0.02, P: 0.03 ∼ 0.10, and C: 0.04 ∼ 0.14). Additionally, the order of the peak total and lateral TFCFs from low to high was normal walking, stair ascent and stand-to-sit, and stair descent and sit-to-stand (P < 0.05), and the peak medial TFCF was stand-to-sit, sit-to-stand, normal walking, stair ascent and stair descent (P < 0.05). The outcomes of this study are valuable for further understanding the knee biomechanics during daily living activities and providing theoretical guidance for the treatments of KOA.

Finding Emergent Gait Patterns May Reduce Progression of Knee Osteoarthritis in a Clinically Relevant Time Frame

A high contact force between the medial femoral condyle and the tibial plateau is the primary cause of medial compartment knee osteoarthritis (OA). A high medial contact force (MCF) during gait has been shown to be correlated to both the knee adduction moment (KAM) and knee flexion/extension moment (KFM). In this study, we used OpenSim Moco to find gait kinematics that reduced the peaks of the KAM, without increasing the peaks of the KFM, which could potentially reduce the MCF and, hence, the progression of knee OA. We used gait data from four knee OA participants. Our simulations decreased both peaks of the KAM without increasing either peak of the KFM. We found that increasing the step width was the primary mechanism, followed by simulations of all participants to reduce the frontal plane lever arm of the ground reaction force vector about the knee, in turn reducing the KAM. Importantly, each participant simulation followed different patterns of kinematic changes to achieve this reduction, which highlighted the need for participant-specific gait modifications. Moreover, we were able to simulate emerging gait patterns within 15 min, enhancing the relevance and potential for the application of developed methods in clinical settings.

Muscle coordination retraining inspired by musculoskeletal simulations reduces knee contact force

Humans typically coordinate their muscles to meet movement objectives like minimizing energy expenditure. In the presence of pathology, new objectives gain importance, like reducing loading in an osteoarthritic joint, but people often do not change their muscle coordination patterns to meet these new objectives. Here we use musculoskeletal simulations to identify simple changes in coordination that can be taught using electromyographic biofeedback, achieving the therapeutic goal of reducing joint loading. Our simulations predicted that changing the relative activation of two redundant ankle plantarflexor muscles—the gastrocnemius and soleus—could reduce knee contact force during walking, but it was unclear whether humans could re-coordinate redundant muscles during a complex task like walking. Our experiments showed that after a single session of walking with biofeedback of summary measures of plantarflexor muscle activation, healthy individuals reduced the ratio of gastrocnemius-to-soleus muscle activation by 25 ± 15% (p = 0.004, paired t test, n = 10). Participants who walked with this “gastrocnemius avoidance” gait pattern reduced late-stance knee contact force by 12 ± 12% (p = 0.029, paired t test, n = 8). Simulation-informed coordination retraining could be a promising treatment for knee osteoarthritis and a powerful tool for optimizing coordination for a variety of rehabilitation and performance applications.

Tibio-Femoral Contact Force Distribution of Knee Before and After Total Knee Arthroplasty: Combined Finite Element and Gait Analysis

Objective

To assess the tibio‐femoral contact forces before and after total knee arthroplasty (TKA) in patients with knee osteoarthritis (KOA) by three‐dimensional (3D) finite element analysis (FEA) models and gait analysis.

Methods

Two hospitalized patients with Kellgren–Lawrence grade IV varus KOA and two healthy subjects were enrolled in this study. Both patients underwent unilateral TKA. FEA models were established based on CT and MR images of the knees of the patients with KOA and healthy subjects. Gait analysis was performed using a three‐dimensional motion capture system with a force plate. Three direction forces at the ankle joints were calculated by inverse dynamic analysis, which provided the load for the FEA models. The total contact forces of the knee joints were also calculated by inverse dynamic analysis to enable comparisons with the results from the FEA models. The total knee contact forces, maximum von Mises stress, and stress distribution of the medial plateau were compared between the patients and healthy subjects. The distributions of the medial plateau force at 2 and 6 months postoperatively were compared with the distributions of the forces preoperatively and those in the healthy subjects.

Results

During static standing, the medial plateau bore the most of the total contact forces in the knees with varus KOA (90.78% for patient 1 and 93.53% for patient 2) compared with 64.75 ± 3.34% of the total force in the healthy knees. At the first and second peaks of the ground reaction force during the stance phase of a gait cycle, the medial plateau bore a much higher percentage of contact forces in patients with KOA (74.78% and 86.48%, respectively, for patient 1; 70.68% and 83.56%, respectively, for patient 2) than healthy subjects (61.06% ± 3.43% at the first peak and 72.09% ± 1.83% at the second peak). Two months after TKA, the percentages of contact forces on the medial tibial plateau were 79.65%–85.19% at the first and second peaks of ground reaction forces during the stance phase of a gait cycle, and the percentages decreased to 53.99% – 68.13% 6 months after TKA.

Conclusion

FEA showed that TKA effectively restored the distribution of tibio‐femoral contact forces during static standing and walking, especially 6 months after the surgery. The changes in the gait were consistent with the changes in the contact force distribution calculated by the FEA model.

Uncertainty in Muscle–Tendon Parameters can Greatly Influence the Accuracy of Knee Contact Force Estimates of Musculoskeletal Models

Understanding the sources of error is critical before models of the musculoskeletal system can be usefully translated. Using in vivo measured tibiofemoral forces, the impact of uncertainty in muscle–tendon parameters on the accuracy of knee contact force estimates of a generic musculoskeletal model was investigated following a probabilistic approach. Population variability was introduced to the routine musculoskeletal modeling framework by perturbing input parameters of the lower limb muscles around their baseline values. Using ground reaction force and skin marker trajectory data collected from six subjects performing body-weight squat, the knee contact force was calculated for the perturbed models. The combined impact of input uncertainties resulted in a considerable variation in the knee contact force estimates (up to 2.1 BW change in the predicted force), especially at larger knee flexion angles, hence explaining up to 70% of the simulation error. Although individual muscle groups exhibited different contributions to the overall error, variation in the maximum isometric force and pathway of the muscles showed the highest impacts on the model outcomes. Importantly, this study highlights parameters that should be personalized in order to achieve the best possible predictions when using generic musculoskeletal models for activities involving deep knee flexion.

Altering the strength of the muscles crossing the lower limb joints only affects knee joint reaction forces

BACKGROUND

Generic musculoskeletal models based on literature data are often used to estimate joint reaction forces (JRFs) that otherwise could only be measured invasively. Estimated JRFs are sensitive to changes in maximum isometric force (F_iso) of the muscles, but these are normally simply scaled using a multiplicative coefficient. The impact of varying F_iso, or strength, of muscles crossing each lower limb joint on estimated JRFs has not been systematically explored in musculoskeletal models of the lower limb.

RESEARCH QUESTION

How do alterations in the strength of joint-crossing muscles influence the lower limb JRF magnitudes computed through a generic musculoskeletal model?

METHODS

By modifying F_iso of muscles crossing hip, knee, ankle, or all joints at once up to ± 40% in 10% increments, thirty-two models were created to simulate the gait of a patient with an instrumented tibial prosthesis (5^th Grand Challenge dataset). A standard workflow (inverse kinematics, static optimization, joint reaction analysis) was utilized to calculate JRFs. Both alterations in JRF magnitudes due to joint crossing muscles' strength modifications and their accuracy against in vivo knee loading measurements were quantified.

RESULTS

The knee JRF was the most sensitive force to changes in the joint-crossing muscles' strength (variations ranging from -37.9 ± 0.5% to +37.9 ± 3.2%), while the hip and ankle JRFs were almost unaffected (maximum variation: +6.1%). Reducing the strength of knee and ankle-crossing muscles and intensifying the strength of hip-crossing muscles lowered the knee JRF. The knee JRF was best estimated (peak error: 0.42 ± 0.15 body weight, root mean squared error: 0.37 ± 0.06 body weight, coefficient of determination: 0.76 ± 0.10) by the model with -40% weakened knee-crossing muscles.

SIGNIFICANCE

Altering strengths mainly affects knee JRF estimated with generic musculoskeletal models, suggesting that personalization of strength of joint-crossing muscles is required for accurate knee JRF estimations. Rehabilitation regimes meant to strengthen muscles crossing a joint should be carefully designed to avoid undesired effects on the other joints.

Techniques to determine knee joint contact forces during squatting: A systematic review

This review article provides an overview of techniques used to determine human knee joint contact forces during squatting. The main two approaches are experimental and theoretical. Thigh calf contact has a significant effect on knee forces and should not be neglected. In this study, data were searched electronically and organized by techniques to find knee joint contact force during squatting theoretically and experimentally. There was a large variation in peak tibiofemoral (CV = 0.45) and patellofemoral (CV = 0.38) contact forces predicted theoretically. However, very little variation was observed between peak tibiofemoral contact forces (CV = 0.12) measured in vivo experimentally but measured knee joint force is available up to a limited knee flexion angle. There was a reduction in knee joint contact forces due to thigh calf contact. Literature of knee joint contact force prediction theoretically during squatting incorporating thigh calf contact force is very limited.

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.

Triceps surae strength balancing as a management option for early-stage knee osteoarthritis: A patient case

BACKGROUND

Knee osteoarthritis is a progressive disease that may require management for years before knee arthroplasty can be considered. Previously reported musculoskeletal models suggest that rebalancing the strength of the triceps surae muscles can reduce the joint loads.

METHODS

A single patient diagnosed with mild/moderate medial left knee osteoarthritis was treated with botulinum toxin injections in the gastrocnemius muscle of the calf, based on the hypothesis that this would rebalance the triceps surae load distribution and reduce tibiofemoral joint loads. Tests were performed before and 4 weeks after injection to record functional clinical scores and to obtain lower limb joint kinematic and kinetic data of walking, which were subsequently analyzed with a musculoskeletal simulation model.

FINDINGS

The patient experienced a clinically relevant improvement in self-reported pain levels in activities-of-daily-living, stair climbing, 6 minutes' treadmill test, range-of-motion, and in the functional knee questionnaire, KOOS. No improvement was seen when performing lunges. The musculoskeletal simulations showed the expected shift in loads between the muscles, reduced knee loads, and improvement of the load symmetry between the legs.

INTERPRETATION

The case corroborates the hypothesis, and this suggests further tests by randomized controlled trials. If confirmed, this simple and reversible medical intervention can improve the management of early-stage knee osteoarthritis.

[Osteotomy close to the knee joint in adults]

Osteotomy close to the knee joint in cases of axial deformities or in the context of corrective interventions is still an elementary component of the spectrum of orthopedic surgery. The precise localization and analysis of the pathology are crucial in order to restore approximate physiological loading conditions by means of a suitable osteotomy procedure. Thus, above all, the medial open-wedge tibial osteotomy has become established as a standard procedure for the treatment of knee osteoarthritis caused by extra-articular tibial varus deformity. Furthermore, corrective osteotomy of the sagittal axis of the tibia and of the distal femur also show good long-term results when the indications are correctly interpreted. The osteotomies described in this article are corrections that should be performed after careful consideration of the biological and mechanical aspects. Correct planning, attention to risk factors and the surgical technique are the main factors that determine the long-term success of the procedure.

Movement Quality Parameters during Gait Assessed by a Single Accelerometer in Subjects with Osteoarthritis and Following Total Joint Arthroplasty

This study’s aim is threefold: (I) Evaluate movement quality parameters of gait in people with hip or knee osteoarthritis (OA) compared to asymptomatic controls from a single trunk-worn 3D accelerometer. (II) Evaluate the sensitivity of these parameters to capture changes at 6-weeks, 3-, 6-, and 12-months following total knee arthroplasty (TKA). (III) Investigate whether observed changes in movement quality from 6-weeks and 12-months post-TKA relates to changes in patient-reported outcome measures (PROMs). We invited 20 asymptomatic controls, 20 people with hip OA, 18 people pre- and post-TKA to our movement lap. They wore a single trunk-worn accelerometer and walked at a self-selected speed. Movement quality parameters (symmetry, complexity, smoothness, and dynamic stability) were calculated from the 3D acceleration signal. Between groups and between timepoints comparisons were made, and changes in movement quality were correlated with PROMs. We found significant differences in symmetry and stability in both OA groups. Post-TKA, most parameters reflected an initial decrease in movement quality at 6-weeks post-TKA, which mostly normalised 6-months post-TKA. Finally, improved movement quality relates to improvements in PROMs. Thus, a single accelerometer can characterise movement quality in both OA groups and post-TKA. The correlation shows the potential to monitor movement quality in a clinical setting to inform objective, data-driven personalised rehabilitation.

Biomechanical Effect of Coronal Alignment and Ligament Laxity in Total Knee Arthroplasty: A Simulation Study

The purposes of this study were to develop a cruciate-retaining total knee arthroplasty musculoskeletal model, which enables the adjustment of ligament length and implant alignment; validate the model; and evaluate the effects of varus/valgus alignment adjustment and unbalanced medial/lateral ligament laxity during gait. A cruciate-retaining total knee arthroplasty musculoskeletal model was constructed and validated against the in vivo contact forces. This model was transformed to 2° varus/valgus alignment of femoral or tibial replacement models and 2° medial/lateral laxity models. The contact forces and ligament tensions of the adjusted models were calculated. The contact forces in the model showed good agreement with the in vivo contact forces. Valgus replacement alignment with balanced ligament models showed a lower contact force at the medial compartment than at the neutral alignment model, whereas the varus replacement alignment with balanced ligament models showed a greater contact force at the medial compartment and medial/posterior cruciate ligament tension. The medial laxity with neutral alignment model showed a similar contact force with decreased medial ligament tension compared to the balanced neutral alignment model, whereas the lateral laxity with the neutral alignment model showed a greater contact force and decreased lateral ligament tension. The cruciate-retaining total knee arthroplasty model was validated using in vivo contact forces (r = 0.939) Two degrees of valgus alignment adjustment with balanced ligament or neutral alignment with 2° of medial laxity can be safe without increasing contact force or ligament tension compared to neutral alignment with a balanced extension gap. However, 2° of varus alignment adjustment with balanced ligament or neutral alignment with 2° of lateral laxity may be unfavorable due to the overloading of the joints and knee ligaments.

Comparison of Optimization Strategies for Musculoskeletal Modeling of the Wrist for Therapy Planning in Case of Total Wrist Arthroplasty

The human wrist joint is an elegant mechanism. The wrist allows the positioning and orienting of the hand to the forearm. The computational modeling of the human hand, especially of the wrist joint, can reveal important information about biomechanical mechanisms and provide the basis for its dysfunction and pathologies. For instance, this could be used for therapy planning in total wrist arthroplasty (TWA). In this study, different optimization methods and sensitivity analyses of anatomical parameters for musculoskeletal modeling were presented. Optimization includes finding the best available value of an objective function, including a variety of different types of objective functions. In the simplest case, optimization consists of maximizing or minimizing a function by systematically choosing input values from within an allowed set and computing the value of the function. Optimization techniques are used in many facets, such as the model building of joints or joint systems such as the wrist. The purpose of this study is to show the variability and influence of the included information for modeling, investigating the biomechanical function and load situation of the joint in representative scenarios. These possibilities to take them into account by an optimization and seem essential for the application of computational modeling to joint pathologies.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

A Computationally Efficient Musculoskeletal Model of the Lower Limb for the Control of Rehabilitation Robots: Assumptions and Validation

We present and validate a computationally efficient lower limb musculoskeletal model for the control of a rehabilitation robot. It is a parametric model that allows the customization of joint kinematics, and it is able to operate in real time. Methods: Since the rehabilitation exercises corresponds to low-speed movements, a quasi-static model can be assumed, and then muscle force coefficients are position dependent. This enables their calculation in an offline stage. In addition, the concept of a single functional degree of freedom is used to minimize drastically the workspace of the stored coefficients. Finally, we have developed a force calculation process based on Lagrange multipliers that provides a closed-form solution; in this way, the problem of dynamic indeterminacy is solved without the need to use an iterative process. Results: The model has been validated by comparing muscle forces estimated by the model with the corresponding electromyography (EMG) values using squat exercise, in which the Spearman’s correlation coefficient is higher than 0.93. Its computational time is lower than 2.5 ms in a conventional computer using MATLAB. Conclusions: This procedure presents a good agreement with the experimental values of the forces, and it can be integrated into real time control systems.

Kinematics and kinetics comparison of ultra-congruent versus medial-pivot designs for total knee arthroplasty by multibody analysis

Nowadays, several configurations of total knee arthroplasty (TKA) implants are commercially available whose designs resulted from clinical and biomechanical considerations. Previous research activities led to the development of the so-called medial-pivot (MP) design. However, the actual benefits of the MP, with respect to other prosthesis designs, are still not well understood. The present work compares the impact of two insert geometries, namely the ultra-congruent (UC) and medial-pivot (MP), on the biomechanical behaviour of a bicondylar total knee endoprosthesis. For this purpose, a multibody model of a lower limb was created alternatively integrating the two implants having the insert geometry discretized. Joint dynamics and contact pressure distributions were evaluated by simulating a squat motion. Results showed a similar tibial internal rotation range of about 3.5°, but an early rotation occurs for the MP design. Furthermore, the discretization of the insert geometry allowed to efficiently derive the contact pressure distributions, directly within the multibody simulation framework, reporting peak pressure values of 33 MPa and 20 MPa for the UC and MP, respectively. Clinically, the presented findings confirm the possibility, through a MP design, to achieve a more natural joint kinematics, consequently improving the post-operative patient satisfaction and potentially reducing the occurrence of phenomena leading to the insert loosening.

Descending 13 real world steps: A dataset and analysis of stair descent

BACKGROUND

Stair descent analysis has been typically limited to laboratory staircases of 4 or 5 steps. To date there has been no report of gait parameters during unconstrained stair descent outside of the laboratory, and few motion capture datasets are publicly available.

RESEARCH QUESTION

We aim to collect a dataset and perform gait analysis for stair descent outside of the laboratory. We aim to measure basic kinematic and kinetic gait parameters and foot placement behavior.

METHODS

We present a public stair descent dataset from 101 unimpaired participants aged 18-35 on an unconstrained 13-step staircase collected using wearable sensors. The dataset consists of kinematics (full-body joint angle and position), kinetics (plantar normal forces, acceleration), and foot placement for 30,609 steps.

RESULTS

We report the lower limb joint angle ranges (30° and 8° for hip flexion and extension, 85° and -11° for knee flexion and extension, and 31° and 28° for ankle dorsi- and plantar-flexion). The self-selected speed was 0.79 ± 0.16 m/s, with cycle duration of 0.97 ± 0.18 s. Mean foot overhang as a percentage of foot length was 17.07 ± 6.66 %, and we calculate that foot size explains only 6% of heel placement variation, but 79% of toe placement variation. We also find a minor but significant asymmetry between left and right maximum hip flexion angle, though all other measured parameters were symmetrical.

SIGNIFICANCE

This is the first quantitative observation of gait data from a large number (n = 101) of participants descending an unconstrained staircase outside of a laboratory. This study enables analysis of gait characteristics including self-selected walking speed and foot placement to better understand typical stair gait behavior. The dataset is a public resource for understanding typical stair descent.

The effect of modelling parameters in the development and validation of knee joint models on ligament mechanics: A systematic review

BACKGROUND

The ligaments in the knee are prone to injury especially during dynamic activities. The resulting instability can have a profound impact on a patient's daily activities and functional capacity. Musculoskeletal knee modelling provides a non-invasive tool for investigating ligament force-strain behaviour in various dynamic scenarios, as well as potentially complementing existing pre-planning tools to optimise surgical reconstructions. However, despite the development and validation of many musculoskeletal knee models, the effect of modelling parameters on ligament mechanics has not yet been systematically reviewed.

OBJECTIVES

This systematic review aimed to investigate the results of the most recent studies using musculoskeletal modelling techniques to create models of the native knee joint, focusing on ligament mechanics and modelling parameters in various simulated movements.

DATA SOURCES

PubMed, ScienceDirect, Google Scholar, and IEEE Xplore.

ELIGIBILITY CRITERIA FOR SELECTING STUDIES

Databases were searched for articles containing any numerical ligament strain or force data on the intact, ACL-deficient, PCL-deficient, or lateral extra-articular reconstructed (LER) knee joints. The studies had to derive these results from musculoskeletal modelling methods. The dates of the publications were between 1 January 1995 and 30 November 2021.

METHOD

A customised data extraction form was created to extract each selected study's critical musculoskeletal model development parameters. Specific parameters of the musculoskeletal knee model development used in each eligible study were independently extracted, including the (1) musculoskeletal model definition (i.e., software used for modelling, knee type, source of geometry, the inclusion of cartilage and menisci, and articulating joints and joint boundary conditions (i.e., number of degrees of freedom (DoF), subjects, type of activity, collected data and type of simulation)), (2) specifically ligaments modelling techniques (i.e., ligament bundles, attachment points, pathway, wrapping surfaces and ligament material properties such as stiffness and reference length), (3) sensitivity analysis, (4) validation approaches, (5) predicted ligament mechanics (i.e., force, length or strain) and (6) clinical applications if available. The eligible papers were then discussed quantitatively and qualitatively with respect to the above parameters.

RESULTS AND DISCUSSION

From the 1004 articles retrieved by the initial electronic search, only 25 met all inclusion criteria. The results obtained by aggregating data reported in the eligible studies indicate that considerable variability in the predicted ligament mechanics is caused by differences in geometry, boundary conditions and ligament modelling parameters.

CONCLUSION

This systematic review revealed that there is currently a lack of consensus on knee ligament mechanics. Despite this lack of consensus, some papers highlight the potential of developing translational tools using musculoskeletal modelling. Greater consistency in model design, incorporation of sensitivity assessment of the model outcomes and more rigorous validation methods should lead to better agreement in predictions for ligament mechanics between studies. The resulting confidence in the musculoskeletal model outputs may lead to the development of clinical tools that could be used for patient-specific treatments.

Optimization Reduces Knee-Joint Forces During Walking and Squatting: Validating the Inverse Dynamics Approach for Full Body Movements on Instrumented Knee Prostheses

Because of the redundancy of our motor system, movements can be performed in many ways. While multiple motor control strategies can all lead to the desired behavior, they result in different joint and muscle forces. This creates opportunities to explore this redundancy, for example, for pain avoidance or reducing the risk of further injury. To assess the effect of different motor control optimization strategies, a direct measurement of muscle and joint forces is desirable, but problematic for medical and ethical reasons. Computational modeling might provide a solution by calculating approximations of these forces. In this study, we used a full-body computational musculoskeletal model to (a) predict forces measured in knee prostheses during walking and squatting and (b) study the effect of different motor control strategies (i.e., minimizing joint force vs. muscle activation) on the joint load and prediction error. We found that musculoskeletal models can accurately predict knee joint forces with a root mean squared error of <0.5 body weight (BW) in the superior direction and about 0.1 BW in the medial and anterior directions. Generally, minimization of joint forces produced the best predictions. Furthermore, minimizing muscle activation resulted in maximum knee forces of about 4 BW for walking and 2.5 BW for squatting. Minimizing joint forces resulted in maximum knee forces of 2.25 BW and 2.12 BW, that is, a reduction of 44% and 15%, respectively. Thus, changing the muscular coordination strategy can strongly affect knee joint forces. Patients with a knee prosthesis may adapt their neuromuscular activation to reduce joint forces during locomotion.

Moderate-intensity treadmill running relieves motion-induced post-traumatic osteoarthritis mice by up-regulating the expression of lncRNA H19

BACKGROUND

The purpose of this study was to explore whether moderate-intensity exercise can alleviate motion-induced post-traumatic osteoarthritis (PTOA) and the expression change of lncRNA H19 during this progression.

METHODS

Twenty-week-old male C57BL/6 mice were randomly divided into five groups: model control group (MC group, n = 6), treadmill model group (M group, n = 6), rehabilitation control group (RC group, n = 6), treadmill model + rehabilitation training group (M + R group, n = 6) and treadmill model + convalescent group (M + C group, n = 6). Paraffin sections were used to observe the pathological changes in the mouse knee joint in each group. A micro-CT was used to scan the knee joint to obtain the morphological indexes of the tibial plateau bone. Real-time PCR was used to detect the mRNA levels of inflammatory factors, synthetic and catabolic factors in cartilage.

RESULTS

After high-intensity exercise for 4 weeks, the inflammation and catabolism of the mouse knee cartilage were enhanced, and the anabolism was weakened. Further study showed that these results were partially reversed after 4-week moderate-intensity training. The results of hematoxylin-eosin staining confirmed this finding. Meanwhile, high-intensity exercise reduced the expression of lncRNA H19 in cartilage, while the expression of lncRNA H19 increased after 4 weeks of moderate-intensity exercise.

CONCLUSION

High-intensity treadmill running can cause injury to the knee cartilage in C57BL/6 mice which leads to PTOA and a decrease of lncRNA H19 expression in cartilage. Moderate-intensity exercise can relieve PTOA and partially reverse lncRNA H19 expression.