A Combined Geometric Morphometric and Discrete Element Modeling Approach for Hip Cartilage Contact Mechanics

Incorporating patient-specific hip orientation from weightbearing computed tomography affects discrete element analysis-computed regional joint contact mechanics in individuals treated with periacetabular osteotomy for hip dysplasia

Computational models of the hip often omit patient-specific functional orientation when placing imaging-derived bony geometry into anatomic landmark-based coordinate systems for application of joint loading schemes. The purpose of this study was to determine if this omission meaningfully alters computed contact mechanics. Discrete element analysis models were created from non-weightbearing (NWB) clinical CT scans of 10 hip dysplasia patients (11 hips) and oriented in the International Society of Biomechanics (ISB) coordinate system (NWB-ISB). Three additional models were generated for each hip by adding patient-specific stance information obtained via weightbearing CT (WBCT) to each ISB-oriented model: (1) patient-specific sagittal tilt added (WBCT-sagittal), (2) coronal and axial rotation from optical motion capture added to (1; WBCT-combo), and (3) WBCT-derived axial, sagittal, and coronal rotation added to (1; WBCT-original). Identical gait cycle loading was applied to all models for a given hip, and computed contact stress and contact area were compared between model initialization techniques. Addition of sagittal tilt did not significantly change whole-joint peak (p = 0.922) or mean (p = 0.871) contact stress or contact area (p = 0.638). Inclusion of motion-captured coronal and axial rotation (WBCT-combo) decreased peak contact stress (p = 0.014) and slightly increased average contact area (p = 0.071) from WBCT-sagittal models. Including all WBCT-derived rotations (WBCT-original) further reduced computed peak contact stress (p = 0.001) and significantly increased contact area (p = 0.001). Variably significant differences (p = 0.001-1.0) in patient-specific acetabular subregion mechanics indicate the importance of functional orientation incorporation for modeling applications in which local contact mechanics are of interest.

Neonatal Hip Loading in Developmental Dysplasia: Finite Element Simulation of Proximal Femur Growth and Treatment

Background

Abnormal prenatal hip joint loading can lead to compromised hip joint function. Early intervention is crucial for favorable outcomes.

Purpose

This study investigates the impact of treatment timing (initiation and duration) on cartilage growth and ossification in the proximal femur of infants with developmental dysplasia of the hip, a condition affecting newborns.

Methods

We used a mechanobiological model to simulate proximal femur growth during treatment durations of 3 months, 6 months, and a late-start treatment.

Results

The findings indicate that the timing of treatment initiation is crucial, while a longer treatment duration does not contribute to improved morphological development of the hip joint.

Conclusions

Mechanobiological models of growth can be used to develop treatments and therapies that correct loading conditions. Growing bone is particularly sensitive to loading conditions, and altered loading during growth can affect bone shape and functionality.

A Correspondence-Based Network Approach for Groupwise Analysis of Patient-Specific Spatiotemporal Data

Methods for statistically analyzing patient-specific data that vary both spatially and over time are currently either limited to summary statistics or require elaborate surface registration. We propose a new method, called correspondence-based network analysis, which leverages particle-based shape modeling to establish correspondence across a population and preserve patient-specific measurements and predictions through statistical analysis. Herein, we evaluated this method using three published datasets of the hip describing cortical bone thickness of the proximal femur, cartilage contact stress, and dynamic joint space between control and patient cohorts to evaluate activity- and group-based differences, as applicable, using traditional statistical parametric mapping (SPM) and our proposed spatially considerate correspondence-based network analysis approach. The network approach was insensitive to correspondence density, while the traditional application of SPM showed decreasing area of the region of significance with increasing correspondence density. In comparison to SPM, the network approach identified broader and more connected regions of significance for all three datasets. The correspondence-based network analysis approach identified differences between groups and activities without loss of subject and spatial specificity which could improve clinical interpretation of results.

Validation of a personalized ligament-constraining discrete element framework for computing ankle joint contact mechanics

BACKGROUND AND OBJECTIVE

Computer simulations of joint contact mechanics have great merit to improve our current understanding of articular ankle pathology. Owed to its computational simplicity, discrete element analysis (DEA) is an encouraging alternative to finite element analysis (FEA). However, previous DEA models lack subject-specific anatomy and may oversimplify the biomechanics of the ankle. The objective of this study was to develop and validate a personalized DEA framework that permits movement of the fibula and incorporates personalized cartilage thickness as well as ligamentous constraints.

METHODS

A linear and non-linear DEA framework, representing cartilage as compressive springs, was established, verified, and validated. Three-dimensional (3D) bony ankle models were constructed from cadaveric lower limb CT scans imaged during application of weight (85 kg) and/or torque (10 Nm). These 3D models were used to generate cartilage thickness and ligament insertion sites based on a previously validated statistical shape model. Ligaments were modelled as non-linear tension-only springs. Validation of contact stress prediction was performed using a simple, axially constrained tibiotalar DEA model against an equivalent FEA model. Validation of ligamentous constraints compared the final position of the ankle mortise to that of the cadaver after application of torque and sequential ligament sectioning. Finally, a combined ligamentous-constraining DEA model was validated for predicted contact stress against an equivalent ligament-constraining FEA model.

RESULTS

The linear and non-linear DEA model reproduced a mean articular contact stress within 0.36 MPa and 0.39 MPa of the FEA calculated stress, respectively. With respect to the ligamentous validation, the DEA ligament-balancing algorithm could reproduce the position of the distal fibula within the ankle mortise to within 0.97 mm of the experimental observed distal fibula. When combining the ligament-constraining and contact stress algorithm, DEA was able to reproduce a mean articular contact stress to within 0.50 MPa of the FEA calculated contact stress.

CONCLUSION

The DEA framework presented herein offers a computationally efficient alternative to FEA for the prediction of contact stress in the ankle joint, manifesting its potential to enhance the mechanical understanding of articular ankle pathologies on both a patient-specific and population-wide level. The novelty of this model lies in its personalized nature, inclusion of the distal tibiofibular joint and the use of non-linear ligament balancing to maintain the physiological ankle joint articulation.

Personalized statistical modeling of soft tissue structures in the knee

Background and Objective: As in vivo measurements of knee joint contact forces remain challenging, computational musculoskeletal modeling has been popularized as an encouraging solution for non-invasive estimation of joint mechanical loading. Computational musculoskeletal modeling typically relies on laborious manual segmentation as it requires reliable osseous and soft tissue geometry. To improve on feasibility and accuracy of patient-specific geometry predictions, a generic computational approach that can easily be scaled, morphed and fitted to patient-specific knee joint anatomy is presented.

Methods: A personalized prediction algorithm was established to derive soft tissue geometry of the knee, originating solely from skeletal anatomy. Based on a MRI dataset (n = 53), manual identification of soft-tissue anatomy and landmarks served as input for our model by use of geometric morphometrics. Topographic distance maps were generated for cartilage thickness predictions. Meniscal modeling relied on wrapping a triangular geometry with varying height and width from the anterior to the posterior root. Elastic mesh wrapping was applied for ligamentous and patellar tendon path modeling. Leave-one-out validation experiments were conducted for accuracy assessment.

Results: The Root Mean Square Error (RMSE) for the cartilage layers of the medial tibial plateau, the lateral tibial plateau, the femur and the patella equaled respectively 0.32 mm (range 0.14–0.48), 0.35 mm (range 0.16–0.53), 0.39 mm (range 0.15–0.80) and 0.75 mm (range 0.16–1.11). Similarly, the RMSE equaled respectively 1.16 mm (range 0.99–1.59), 0.91 mm (0.75–1.33), 2.93 mm (range 1.85–4.66) and 2.04 mm (1.88–3.29), calculated over the course of the anterior cruciate ligament, posterior cruciate ligament, the medial and the lateral meniscus.

Conclusion: A methodological workflow is presented for patient-specific, morphological knee joint modeling that avoids laborious segmentation. By allowing to accurately predict personalized geometry this method has the potential for generating large (virtual) sample sizes applicable for biomechanical research and improving personalized, computer-assisted medicine.

Joint contact stress improves in dysplastic hips after periacetabular osteotomy but remains higher than in normal hips

AIM

The purpose of this study was to use computational modeling to determine if surgical correction of hip dysplasia restores hip contact mechanics to those of asymptomatic, radiographically normal hips.

METHODS

Discrete element analysis (DEA) was used to compute joint contact stresses during the stance phase of normal walking gait for 10 individuals with radiographically normal, asymptomatic hips and 10 age- and weight-matched patients with acetabular dysplasia who underwent periacetabular osteotomy (PAO).

RESULTS

Mean and peak contact stresses were higher (p < 0.001 and p = 0.036, respectively) in the dysplastic hips than in the matched normal hips. PAO normalised standard radiographic measurements and medialised the location of computed contact stress within the joint. Mean contact stress computed in dysplastic hips throughout the stance phase of gait (median 5.5 MPa, [IQR 3.9-6.1 MPa]) did not significantly decrease after PAO (3.7 MPa, [IQR 3.2-4.8]; p = 0.109) and remained significantly (p < 0.001) elevated compared to radiographically normal hips (2.4 MPa, [IQR 2.2-2.8 MPa]). Peak contact stress demonstrated a similar trend. Joint contact area during the stance phase of gait in the dysplastic hips increased significantly (p = 0.036) after PAO from 395 mm2 (IQR 378-496 mm2) to 595 mm2 (IQR 474-660 mm2), but remained significantly smaller (p = 0.001) than that for radiographically normal hips (median 1120 mm2, IQR 853-1444 mm2).

CONCLUSIONS

While contact mechanics in dysplastic hips more closely resembled those of normal hips after PAO, the elevated contact stresses and smaller contact areas remaining after PAO indicate ongoing mechanical abnormalities should be expected even after radiographically successful surgical correction.

Wear patterns in knee OA correlate with native limb geometry

Background: To date, the amount of cartilage loss is graded by means of discrete scoring systems on artificially divided regions of interest (ROI). However, optimal statistical comparison between and within populations requires anatomically standardized cartilage thickness assessment. Providing anatomical standardization relying on non-rigid registration, we aim to compare morphotypes of a healthy control cohort and virtual reconstructed twins of end-stage knee OA subjects to assess the shape-related knee OA risk and to evaluate possible correlations between phenotype and location of cartilage loss.

Methods: Out of an anonymized dataset provided by the Medacta company (Medacta International SA, Castel S. Pietro, CH), 798 end-stage knee OA cases were extracted. Cartilage wear patterns were observed by computing joint space width. The three-dimensional joint space width data was translated into a two-dimensional pixel image, which served as the input for a principal polynomial autoencoder developed for non-linear encoding of wear patterns. Virtual healthy twin reconstruction enabled the investigation of the morphology-related risk for OA requiring joint arthroplasty.

Results: The polynomial autoencoder revealed 4 dominant, orthogonal components, accounting for 94% of variance in the latent feature space. This could be interpreted as medial (54.8%), bicompartmental (25.2%) and lateral (9.1%) wear. Medial wear was subdivided into anteromedial (11.3%) and posteromedial (10.4%) wear. Pre-diseased limb geometry had a positive predictive value of 0.80 in the prediction of OA incidence (r 0.58, p &lt; 0.001).

Conclusion: An innovative methodological workflow is presented to correlate cartilage wear patterns with knee joint phenotype and to assess the distinct knee OA risk based on pre-diseased lower limb morphology. Confirming previous research, both alignment and joint geometry are of importance in knee OA disease onset and progression.

Quantification of femoral position on radiographic acetabular coverage in children

An accurate assessment of the radiographic acetabular coverage is essential for clinical diagnosis or surgical decision-making in hip disorders. This study aimed to evaluate the effect of femoral position on acetabular coverage and to predict the actual acetabular coverage from nonstandard radiographs. A total of 21 children (34 hips) with normative acetabular coverage were screened in this retrospective study. The Mimics-based local-rotation fluoroscopy simulation method was used to tilt, incline, and rotate the femur in 4° increments within the range of femoral motion. The acetabular coverage, namely acetabular-head index (AHI) and center-edge angle (CEA), increased with femoral abduction but decreased with other motions. Compared to the femoral neutral position, no significant differences were identified in AHI with the rotation (range: 0°-16°) and in CEA with the tilt (range: -20°-4°), inclination (range: 0°-4°), or rotation (range: -8°-40°). The linear regression analysis showed that the CEA increased by about 0.20° for each 1° increase in femoral inclination and decreased by about 0.01°, 0.07°, 0.06°, or 0.07° for each 1° increase in internal rotation, external rotation, flexion, or extension, respectively. And a more significant change in AHI was observed. All femoral malpositions, especially the inclination, affected radiographic acetabular coverage in children. Therefore, each pelvic radiograph should assess potential femoral malpositioning before diagnosing hip disorders. This study will assist surgeons in predicting the acetabular coverage on nonstandard radiographs.

Personalised statistical modelling of soft tissue structures in the ankle

BACKGROUND AND OBJECTIVE

Revealing the complexity behind subject-specific ankle joint mechanics requires simultaneous analysis of three-dimensional bony and soft-tissue structures. 3D musculoskeletal models have become pivotal in orthopedic treatment planning and biomechanical research. Since manual segmentation of these models is time-consuming and subject to manual errors, (semi-) automatic methods could improve the accuracy and enlarge the sample size of personalised 'in silico' biomechanical experiments and computer-assisted treatment planning. Therefore, our aim was to automatically predict ligament paths, cartilage topography and thickness in the ankle joint based on statistical shape modelling.

METHODS

A personalised cartilage and ligamentous prediction algorithm was established using geometric morphometrics, based on an 'in-house' generated lower limb skeletal model (N = 542), tibiotalar cartilage (N = 60) and ankle ligament segmentations (N = 10). For cartilage, a population-averaged thickness map was determined by use of partial least-squares regression. Ligaments were wrapped around bony contours based on iterative shortest path calculation. Accuracy of ligament path and cartilage thickness prediction was quantified using leave-one-out experiments. The novel personalised thickness prediction was compared with a constant cartilage thickness of 1.50 mm by use of a paired sample T-test.

RESULTS

Mean distance error of cartilage and ligament prediction was 0.12 mm (SD 0.04 mm) and 0.54 mm (SD 0.05 mm), respectively. No significant differences were found between the personalised thickness cartilage and segmented cartilage of the tibia (p = 0.73, CI [-1.60 .10^-17, 1.13 .10^-17]) and talus (p = 0.95, CI[ -1.35 .10^-17, 1.28 .10^-17]). For the constant thickness cartilage, a statistically significant difference was found in 89% and 92% of the tibial (p < 0.001, CI [0.51, 0.58]) and talar (p < 0.001, CI [0.33, 0.40]) cartilage area.

CONCLUSIONS

In this study, we described a personalised prediction algorithm of cartilage and ligaments in the ankle joint. We were able to predict cartilage and main ankle ligaments with submillimeter accuracy. The proposed method has a high potential for generating large (virtual) sample sizes in biomechanical research and mitigates technological advances in computer-assisted orthopaedic surgery.

Accurate Arthroscopic Cam Resection Normalizes Contact Stresses in Patients With Femoroacetabular Impingement

BACKGROUND

Femoroacetabular impingement (FAI) is increasingly recognized as a cause of hip pain in young adults. The condition leads to chondrolabral separation and chondral delamination and eventually predisposes to osteoarthritis of the hip. FAI that inflicts cartilage damage has been observed in hips with abnormal morphological characteristics and is related to a long-term evolution toward osteoarthritis. Arthroscopic surgery, which allows for correction of morphological characteristics and restores impingement-free motions, is the current standard of treatment.

HYPOTHESIS

Arthroscopic cam resection can restore the normal mechanical environment of the hip joint in cam-type FAI.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Patient-specific discrete element models from 10 patients with cam-type FAI (all male; age, 18-40 years) were defined based on preoperative computed tomography scans and postoperative magnetic resonance imaging (MRI) scans. Complete cam resection postoperatively on MRI was confirmed with alpha angles <55°. The preoperative and postoperative peak contact stress findings during impingement testing were compared against a matched control group.

RESULTS

Peak contact stress was significantly elevated in patients with cam-type FAI during impingement testing, with increasing amounts of internal hip rotation (26.6 ± 11.64 MPa in cam patients preoperatively, 12.1 ± 4.62 MPa in those same patients postoperatively, and 11.4 ± 1.72 MPa in the virtual control group during impingement testing at 20° of internal hip rotation; P < .01). This effect was normalized after arthroscopic cam resection and loading patterns matched those of the control group.

CONCLUSION

Accurate arthroscopic cam resection restored the normal peak joint contact stresses in the hip joint. This highlights the importance of early and complete cam resections in the face of a positive diagnosis of cam-type FAI.

CLINICAL RELEVANCE

Treatment of cam-type FAI effectively normalizes hip joint contact mechanics.