Development and Evaluation of a Subject-Specific Lower Limb Model With an Eleven-Degrees-of-Freedom Natural Knee Model Using Magnetic Resonance and Biplanar X-Ray Imaging During a Quasi-Static Lunge

Methodology to identify subject-specific dynamic laxity tests to stretch individual parts of knee ligaments

The mechanical properties of ligaments are important for multiple applications and are often estimated from laxity tests. However, the typical laxity tests are not optimized for this application and, a potential exists to develop better laxity tests in this respect. Therefore, the purpose of this study was to develop a methodology to identify optimal, dynamic laxity tests that isolate the stretch of the individual ligaments from each other. To this end, we applied an existing rigid body-based knee model and a dataset of ∼100.000 random samples of applied forces (0-150 N), moments (0-10 Nm) and knee flexion angles (0-90°) through Monte Carlo Simulations. For each modelled ligament bundle, we identified ten load cases; one producing the highest force and nine equally spaced between the maximal and zero force, where the maximal force in all other ligament bundles were minimized. We compared these novel laxity tests to standard internal/external and varus/valgus laxity tests using an isolation metric. We found that no laxity test could stretch the anterior part of the posterior cruciate and medial cruciate ligaments (PCL and MCL), whereas for all other ligaments, except the posterior PCL, the new laxity tests isolated the ligament stretch 28 % to 450 % better than standard tests. From our study, we conclude that it is possible to define better laxity tests than currently exist and these may be highly relevant for determination of mechanical properties of ligaments in vivo. Future studies should generalize our results and translate them to modern laxity measurements technologies.

Numerical study of osteophyte effects on preoperative knee functionality in patients undergoing total knee arthroplasty

Osteophytes are routinely removed during total knee arthroplasty, yet the preoperative planning currently relies on preoperative computed tomography (CT) scans of the patient's osteoarthritic knee, typically including osteophytic features. This complicates the surgeon's ability to anticipate the exact biomechanical effects of osteophytes and the consequences of their removal before the operation. The aim of this study was to investigate the effect of osteophytes on ligament strains and kinematics, and ascertain whether the osteophyte volume and location determine the extent of this effect. We segmented preoperative CT scans of 21 patients, featuring different osteophyte severity, using image-based active appearance models trained to identify the osteophytic and preosteophytic bone geometries and estimate the cartilage thickness in the segmented surfaces. The patients' morphologies were used to scale a template musculoskeletal knee model. Osteophytes induced clinically relevant changes to the knee's functional behavior, but these were variable and patient-specific. Generally, severe osteophytic knees significantly strained the oblique popliteal ligament (OPL) and posterior capsule (PC) relative to the preosteophytic state. Furthermore, there was a marked effect on the lateral collateral ligament and anterolateral ligament (ALL) strains compared to mild and moderate osteophytic knees, and concurrent alterations in the tibial lateral-medial translation and external-internal rotation. We found a strong correlation between the OPL, PC, and ALL strains and posterolateral condylar and tibial osteophytes, respectively. Our findings may have implications for the preoperative planning in total knee arthroplasty, toward reproducing the physiological knee biomechanics as close as feasibly possible.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

Prediction of ground reaction forces and moments during walking in children with cerebral palsy

Introduction

Gait analysis is increasingly used to support clinical decision-making regarding diagnosis and treatment planning for movement disorders. As a key part of gait analysis, inverse dynamics can be applied to estimate internal loading conditions during movement, which is essential for understanding pathological gait patterns. The inverse dynamics calculation uses external kinetic information, normally collected using force plates. However, collection of external ground reaction forces (GRFs) and moments (GRMs) can be challenging, especially in subjects with movement disorders. In recent years, a musculoskeletal modeling-based approach has been developed to predict external kinetics from kinematic data, but its performance has not yet been evaluated for altered locomotor patterns such as toe-walking. Therefore, the goal of this study was to investigate how well this prediction method performs for gait in children with cerebral palsy.

Methods

The method was applied to 25 subjects with various forms of hemiplegic spastic locomotor patterns. Predicted GRFs and GRMs, in addition to associated joint kinetics derived using inverse dynamics, were statistically compared against those based on force plate measurements.

Results

The results showed that the performance of the predictive method was similar for the affected and unaffected limbs, with Pearson correlation coefficients between predicted and measured GRFs of 0.71–0.96, similar to those previously reported for healthy adults, despite the motor pathology and the inclusion of toes-walkers within our cohort. However, errors were amplified when calculating the resulting joint moments to an extent that could influence clinical interpretation.

Conclusion

To conclude, the musculoskeletal modeling-based approach for estimating external kinetics is promising for pathological gait, offering the possibility of estimating GRFs and GRMs without the need for force plate data. However, further development is needed before implementation within clinical settings becomes possible.

Increased Femoral Anteversion Does Not Lead to Increased Joint Forces During Gait in a Cohort of Adolescent Patients

Orthopedic complications were previously reported for patients with increased femoral anteversion. A more comprehensive analysis of the influence of increased femoral anteversion on joint loading in these patients is required to better understand the pathology and its clinical management. Therefore, the aim was to investigate lower-limb kinematics, joint moments and forces during gait in adolescent patients with increased, isolated femoral anteversion compared to typically developing controls. Secondly, relationships between the joint loads experienced by the patients and different morphological and kinematic features were investigated. Patients with increased femoral anteversion (n = 42, 12.8 ± 1.9 years, femoral anteversion: 39.6 ± 6.9°) were compared to typically developing controls (n = 9, 12.0 ± 3.0 years, femoral anteversion: 18.7 ± 4.1°). Hip and knee joint kinematics and kinetics were calculated using subject-specific musculoskeletal models. Differences between patients and controls in the investigated outcome variables (joint kinematics, moments, and forces) were evaluated through statistical parametric mapping with Hotelling T2 and t-tests (α = 0.05). Canonical correlation analyses (CCAs) and regression analyses were used to evaluate within the patients’ cohort the effect of different morphological and kinematic predictors on the outcome variables. Predicted compressive proximo-distal loads in both hip and knee joints were significantly reduced in patients compared to controls. A gait pattern characterized by increased knee flexion during terminal stance (KneeFlex_tSt) was significantly correlated with hip and knee forces, as well as with the resultant force exerted by the quadriceps on the patella. On the other hand, hip internal rotation and in-toeing, did not affect the loads in the joints. Based on the finding of the CCAs and linear regression analyses, patients were further divided into two subgroups based KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented a significantly higher femoral anteversion than those with normal KneeFlex_tSt. Patients with excessive KneeFlex_tSt presented significantly larger quadriceps forces on the patella and a larger posteriorly-oriented shear force at the knee, compared to patients with normal KneeFlex_tSt, but both patients’ subgroups presented only limited differences in terms of joint loading compared to controls. This study showed that an altered femoral morphology does not necessarily lead to an increased risk of joint overloading, but instead patient-specific kinematics should be considered.

Subject-Specific Modeling of Femoral Torsion Influences the Prediction of Hip Loading During Gait in Asymptomatic Adults

Hip osteoarthritis may be caused by increased or abnormal intra-articular forces, which are known to be related to structural articular cartilage damage. Femoral torsional deformities have previously been correlated with hip pain and labral damage, and they may contribute to the onset of hip osteoarthritis by exacerbating the effects of existing pathoanatomies, such as cam and pincer morphologies. A comprehensive understanding of the influence of femoral morphotypes on hip joint loading requires subject-specific morphometric and biomechanical data on the movement characteristics of individuals exhibiting varying degrees of femoral torsion. The aim of this study was to evaluate hip kinematics and kinetics as well as muscle and joint loads during gait in a group of adult subjects presenting a heterogeneous range of femoral torsion by means of personalized musculoskeletal models. Thirty-seven healthy volunteers underwent a 3D gait analysis at a self-selected walking speed. Femoral torsion was evaluated with low-dosage biplanar radiography. The collected motion capture data were used as input for an inverse dynamics analysis. Personalized musculoskeletal models were created by including femoral geometries that matched each subject’s radiographically measured femoral torsion. Correlations between femoral torsion and hip kinematics and kinetics, hip contact forces (HCFs), and muscle forces were analyzed. Within the investigated cohort, higher femoral antetorsion led to significantly higher anteromedial HCFs during gait (medial during loaded stance phase and anterior during swing phase). Most of the loads during gait are transmitted through the anterior/superolateral quadrant of the acetabulum. Correlations with hip kinematics and muscle forces were also observed. Femoral antetorsion, through altered kinematic strategies and different muscle activations and forces, may therefore lead to altered joint mechanics and pose a risk for articular damage. The method proposed in this study, which accounts for both morphological and kinematic characteristics, might help in identifying in a clinical setting patients who, as a consequence of altered femoral torsional alignment, present more severe functional impairments and altered joint mechanics and are therefore at a higher risk for cartilage damage and early onset of hip osteoarthritis.

Subject-specific muscle properties from diffusion tensor imaging significantly improve the accuracy of musculoskeletal models

Musculoskeletal modelling is an important platform on which to study the biomechanics of morphological structures in vertebrates and is widely used in clinical, zoological and palaeontological fields. The popularity of this approach stems from the potential to non-invasively quantify biologically important but difficult-to-measure functional parameters. However, while it is known that model predictions are highly sensitive to input values, it is standard practice to build models by combining musculoskeletal data from different sources resulting in 'generic' models for a given species. At present, there are little quantitative data on how merging disparate anatomical data in models impacts the accuracy of these functional predictions. This issue is addressed herein by quantifying the accuracy of both subject-specific human limb models containing individualised muscle force-generating properties and models built using generic properties from both elderly and young individuals, relative to experimental muscle torques obtained from an isokinetic dynamometer. The results show that subject-specific models predict isokinetic muscle torques to a greater degree of accuracy than generic models at the ankle (root-mean-squared error - 7.9% vs. 49.3% in elderly anatomy-based models), knee (13.2% vs. 57.3%) and hip (21.9% vs. 32.8%). These results have important implications for the choice of musculoskeletal properties in future modelling studies, and the relatively high level of accuracy achieved in the subject-specific models suggests that such models can potentially address questions about inter-subject variations of muscle functions. However, despite relatively high levels of overall accuracy, models built using averaged generic muscle architecture data from young, healthy individuals may lack the resolution and accuracy required to study such differences between individuals, at least in certain circumstances. The results do not wholly discourage the continued use of averaged generic data in musculoskeletal modelling studies but do emphasise the need for to maximise the accuracy of input values if studying intra-species form-function relationships in the musculoskeletal system.