Intra and inter-rater variability in the construction of patient-specific musculoskeletal model

BACKGROUND

Variations observed in biomechanical studies might be attributed to errors made by operators during the construction of musculoskeletal models, rather than being solely attributed to patient-specific geometry.

RESEARCH QUESTION

What is the impact of operator errors on the construction of musculoskeletal models, and how does it affect the estimation of muscle moment arms and hip joint reaction forces?

METHODS

Thirteen independent operators participated in defining the muscle model, while a single operator performed 13 repetitions to define the muscle model based on 3D bone geometry. For each model, the muscle moment arms relative to the hip joint center of rotation was evaluated. Additionally, the hip joint reaction force during one-legged stance was assessed using static inverse optimization.

RESULTS

The results indicated high levels of consistency, as evidenced by the intra- rater and inter-rater agreement measured by the Intraclass Correlation Coefficient (ICC), which yielded values of 0.95 and 0.99, respectively. However, the estimated muscle moment arms exhibited an error of up to 16 mm compared to the reference musculoskeletal model. It was found that muscles attached to prominent anatomical landmarks were specified with greater accuracy than those attached over larger areas. Furthermore, the variability in estimated moment arms contributed to variations of up to 12% in the hip joint reaction forces.

SIGNIFICANCE

Both moment arm and muscle force demonstrated significantly lower variability when assessed by a single operator, suggesting the preference for employing a single operator in the creation of musculoskeletal models for clinical biomechanical studies.

Movement control during one-leg standing is important for the bone mineral density maintenance or improvement

INTRODUCTION

Both hip fracture and bone mineral density (BMD) decline on the non-fractured side are more likely to occur within 1 year. However, there are no longitudinal study reports on the factors associated with BMD maintenance or improvement within the first year after hip fracture. This study aimed to investigate the factors influencing the neck BMD maintenance or improvement in patients with hip fractures from within 2 weeks-6 months after surgery.

MATERIALS AND METHODS

Patients were hip fracture after surgery and were divided into two groups: Among neck BMD changes (6 months minus 2 weeks after surgery) were calculated. Based on among neck BMD change, patients were classified into the BMD maintenance or improvement (change ≥ 0) and the BMD decrease groups (change < 0). Propensity score matching was performed to adjust for confounding factors. To predict the factors affecting neck BMD, hierarchical logistic regression analysis was performed. The dependent variable was the BMD maintenance or improvement group and the BMD decrease group. The independent variables were basic and medical information, and physical functions.

RESULTS

The hierarchical logistic regression analysis results showed that movement control during one-leg standing affected femoral neck BMD independently from age, sex. The odds ratio for movement control during one-leg standing was 8.22. The discrimination rate of the model was 69.7%.

CONCLUSION

This study suggested that adequate movement control during one-leg standing is important to maintain or improve neck BMD.

Is Anterior Rotation of the Acetabulum Necessary to Normalize Joint Contact Pressure in Periacetabular Osteotomy? A Finite-element Analysis Study

BACKGROUND

Inappropriate sagittal plane correction can result in an increased risk of osteoarthritis progression after periacetabular osteotomy (PAO). Individual and postural variations in sagittal pelvic tilt, along with acetabular deformity, affect joint contact mechanics in dysplastic hips and may impact the direction and degree of acetabular correction. Finite-element analyses that account for physiologic pelvic tilt may provide valuable insight into the effect of PAO on the contact mechanics of dysplastic hips, which may lead to improved acetabular correction during PAO.

QUESTIONS/PURPOSES

We performed virtual PAO using finite-element models with reference to the standing pelvic position to clarify (1) whether lateral rotation of the acetabulum normalizes the joint contact pressure, (2) risk factors for abnormal contact pressure after lateral rotation of the acetabulum, and (3) whether additional anterior rotation of the acetabulum further reduces contact pressure.

METHODS

Between 2016 and 2020, 85 patients (92 hips) underwent PAO to treat hip dysplasia. Eighty-two patients with hip dysplasia (lateral center-edge angle < 20°) were included. Patients with advanced osteoarthritis, femoral head deformity, prior hip or spine surgery, or poor-quality images were excluded. Thirty-eight patients (38 hips) were eligible to participate in this study. All patients were women, with a mean age of 39 ± 10 years. Thirty-three women volunteers without a history of hip disease were reviewed as control participants. Individuals with a lateral center-edge angle < 25° or poor-quality images were excluded. Sixteen individuals (16 hips) with a mean age of 36 ± 7 years were eligible as controls. Using CT images, we developed patient-specific three-dimensional surface hip models with the standing pelvic position as a reference. The loading scenario was based on single-leg stance. Four patterns of virtual PAO were performed in the models. First, the acetabular fragment was rotated laterally in the coronal plane so that the lateral center-edge angle was 30°; then, anterior rotation in the sagittal plane was added by 0°, 5°, 10°, and 15°. We developed finite-element models for each acetabular position and performed a nonlinear contact analysis to calculate the joint contact pressure of the acetabular cartilage. The normal range of the maximum joint contact pressure was calculated to be < 4.1 MPa using a receiver operating characteristic curve. A paired t-test or Wilcoxon signed rank test with Bonferroni correction was used to compare joint contact pressures among acetabular positions. We evaluated the association of joint contact pressure with the patient-specific sagittal pelvic tilt and acetabular version and coverage using Pearson or Spearman correlation coefficients. An exploratory univariate logistic regression analysis was performed to identify which of the preoperative factors (CT measurement parameters and sagittal pelvic tilt) were associated with abnormal contact pressure after lateral rotation of the acetabulum. Variables with p values < 0.05 (anterior center-edge angle and sagittal pelvic tilt) were included in a multivariable model to identify the independent influence of each factor.

RESULTS

Lateral rotation of the acetabulum decreased the median maximum contact pressure compared with that before virtual PAO (3.7 MPa [range 2.2-6.7] versus 7.2 MPa [range 4.1-14 MPa], difference of medians 3.5 MPa; p < 0.001). The resulting maximum contact pressures were within the normal range (< 4.1 MPa) in 63% of the hips (24 of 38 hips). The maximum contact pressure after lateral acetabular rotation was negatively correlated with the standing pelvic tilt (anterior pelvic plane angle) (ρ = -0.52; p < 0.001) and anterior center-edge angle (ρ = -0.47; p = 0.003). After controlling for confounding variables such as the lateral center-edge angle and sagittal pelvic tilt, we found that a decreased preoperative anterior center-edge angle (per 1°; odds ratio 1.14 [95% CI 1.01-1.28]; p = 0.01) was independently associated with elevated contact pressure (≥ 4.1 MPa) after lateral rotation; a preoperative anterior center-edge angle < 32° in the standing pelvic position was associated with elevated contact pressure (sensitivity 57%, specificity 96%, area under the curve 0.77). Additional anterior rotation further decreased the joint contact pressure; the maximum contact pressures were within the normal range in 74% (28 of 38 hips), 76% (29 of 38 hips), and 84% (32 of 38 hips) of the hips when the acetabulum was rotated anteriorly by 5°, 10°, and 15°, respectively.

CONCLUSION

Via virtual PAO, normal joint contact pressure was achieved in 63% of patients by normalizing the lateral acetabular coverage. However, lateral acetabular rotation was insufficient to normalize the joint contact pressure in patients with more posteriorly tilted pelvises and anterior acetabular deficiency. In patients with a preoperative anterior center-edge angle < 32° in the standing pelvic position, additional anterior rotation is expected to be a useful guide to normalize the joint contact pressure.

CLINICAL RELEVANCE

This virtual PAO study suggests that biomechanics-based planning for PAO should incorporate not only the morphology of the hip but also the physiologic pelvic tilt in the weightbearing position in order to customize acetabular reorientation for each patient.

Hip stress distribution - Predictor of dislocation in hip arthroplasties. A retrospective study of 149 arthroplasties

Dislocation after hip arthroplasty is still a major concern. Recent study of the volumetric wear of the cup has suggested that stresses studied in a one-legged stance model could predispose arthroplasty dislocation. The aim of this work was to study whether biomechanical parameters of contact stress distribution in total hip arthroplasty during a neutral hip position can predict a higher possibility of the arthroplasty dislocating. Biomechanical parameters were determined using 3-dimensional mathematical models of the one-legged stance within the HIPSTRESS method. Geometrical parameters were measured from standard anteroposterior X-ray images of the pelvis and proximal femora. Fifty-five patients subjected to total hip arthroplasty that later suffered dislocation of the head and, for comparison, ninety-four total hip arthroplasties that were functional at least 10 years after the implantation, were included in the study. Arthroplasties that suffered dislocation had on average a 6% higher resultant hip force than the control group (p = 0.004), 11% higher peak stress on the load-bearing area (p = 0.001) and a 50% more laterally positioned stress pole (p = 0.026), all parameters being less favorable in the group of unstable arthroplasties. There was no statistically significant difference in the gradient index or in the functional angle of the weight bearing. Our study showed that arthroplasties that show a tendency to push the head out of the cup in the representative body position—the one-legged stance—are prone to dislocation. An unfavorable resultant hip force, peak stress on the load bearing and laterally positioned stress pole are predictors of arthroplasty dislocation.

Effect of cup inclination on predicted contact stress-induced volumetric wear in total hip replacement

In order to increase the lifetime of the total hip endoprosthesis, it is necessary to understand mechanisms leading to its failure. In this work, we address volumetric wear of the artificial cup, in particular the effect of its inclination with respect to the vertical. Volumetric wear was calculated by using mathematical models for resultant hip force, contact stress and penetration of the prosthesis head into the cup. Relevance of the dependence of volumetric wear on inclination of the cup (its abduction angle ϑA) was assessed by the results of 95 hips with implanted endoprosthesis. Geometrical parameters obtained from standard antero-posterior radiographs were taken as input data. Volumetric wear decreases with increasing cup abduction angle ϑA. The correlation within the population of 95 hips was statistically significant (P = 0.006). Large cup abduction angle minimises predicted volumetric wear but may increase the risk for dislocation of the artificial head from the cup in the one-legged stance. Cup abduction angle and direction of the resultant hip force may compensate each other to achieve optimal position of the cup with respect to wear and dislocation in the one-legged stance for a particular patient.

Polyethylene wear is related to patient-specific contact stress in THA

BACKGROUND

General numerical models of polyethylene wear and THA simulators suggest contact stresses influence wear. These models do not account for some patient-specific factors. Whether the relationship between patient-specific contact stress and wear apply in vivo is unclear.

QUESTIONS/PURPOSES

We therefore determined whether (1) contact stress distribution at the prosthesis-cup interface and (2) hip geometry and cup inclination are related to wear in vivo.

METHODS

We retrospectively reviewed the radiographs of 80 patients who had aseptic loosening of their THAs as determined by radiographic criteria. We determined linear penetration and volumetric wear using postoperative and last followup radiographs. Contact stress distribution was determined by the HIPSTRESS method. The biomechanical model was scaled to fit the patient's musculoskeletal geometry of the pelvis, trochanteric position, and cup inclination using the standard postoperative radiograph.

RESULTS

Linear penetration and volumetric wear correlated with peak contact stress. Polyethylene wear was greater in THAs with a medial position of the greater trochanter and smaller inclination of the acetabular cup.

CONCLUSIONS

Our observations suggest wear is specific to contact stresses in vivo.

CLINICAL RELEVANCE

Long-term wear in a THA can be estimated using contact stress analysis based on analysis of the postoperative AP radiograph.