Subject-specific biomechanical analysis to estimate locations susceptible to osteoarthritis-Finite element modeling and MRI follow-up of ACL reconstructed patients.

The aims of this case-control study were to: (1) Identify cartilage locations and volumes at risk of osteoarthritis (OA) using subject-specific finite element (FE) models; (2) Quantify the relationships between the simulated biomechanical parameters and T2 and T1ρ relaxation times of magnetic resonance imaging (MRI). We created subject-specific FE models for seven patients with anterior cruciate ligament (ACL) reconstruction and six controls based on a previous proof-of-concept study. We identified locations and cartilage volumes susceptible to OA, based on maximum principal stresses and absolute maximum shear strains in cartilage exceeding thresholds of 7 MPa and 32%, respectively. The locations and volumes susceptible to OA were compared qualitatively and quantitatively against 2-year longitudinal changes in T2 and T1ρ relaxation times. The degeneration volumes predicted by the FE models, based on excessive maximum principal stresses, were significantly correlated (r = 0.711, p < 0.001) with the degeneration volumes determined from T2 relaxation times. There was also a significant correlation between the predicted stress values and changes in T2 relaxation time (r = 0.649, p < 0.001). Absolute maximum shear strains and changes in T1ρ relaxation time were not significantly correlated. Five out of seven patients with ACL reconstruction showed excessive maximum principal stresses in either one or both tibial cartilage compartments, in agreement with follow-up information from MRI. Expectedly, for controls, the FE models and follow-up information showed no degenerative signs. Our results suggest that the presented modelling methodology could be applied to prospectively identify ACL reconstructed patients at risk of biomechanically driven OA, particularly by the analysis of maximum principal stresses of cartilage.

Identification of locations susceptible to osteoarthritis in patients with anterior cruciate ligament reconstruction: Combining knee joint computational modelling with follow-up T1ρ and T2 imaging.

BACKGROUND: Finite element modelling can be used to evaluate altered loading conditions and failure locations in knee joint tissues. One limitation of this modelling approach has been experimental comparison. The aims of this proof-of-concept study were: 1) identify areas susceptible to osteoarthritis progression in anterior cruciate ligament reconstructed patients using finite element modelling; 2) compare the identified areas against changes in T2 and T1ρ values between 1-year and 3-year follow-up timepoints. METHODS: Two patient-specific finite element models of knee joints with anterior cruciate ligament reconstruction were created. The knee geometry was based on clinical magnetic resonance imaging and joint loading was obtained via motion capture. We evaluated biomechanical parameters linked with cartilage degeneration and compared the identified risk areas against T2 and T1ρ maps. FINDINGS: The risk areas identified by the finite element models matched the follow-up magnetic resonance imaging findings. For Patient 1, excessive values of maximum principal stresses and shear strains were observed in the posterior side of the lateral tibial and femoral cartilage. For Patient 2, high values of maximum principal stresses and shear strains of cartilage were observed in the posterior side of the medial joint compartment. For both patients, increased T2 and T1ρ values between the follow-up times were observed in the same areas. INTERPRETATION: Finite element models with patient-specific geometries and motions and relatively simple material models of tissues were able to identify areas susceptible to post-traumatic knee osteoarthritis. We suggest that the methodology presented here may be applied in large cohort studies.

Comparison between kinetic and kinetic-kinematic driven knee joint finite element models.

Use of knee joint finite element models for diagnostic purposes is challenging due to their complexity. Therefore, simpler models are needed for studies where a high number of patients need to be analyzed, without compromising the results of the model. In this study, more complex, kinetic (forces and moments) and simpler, kinetic-kinematic (forces and angles) driven finite element models were compared during the stance phase of gait. Patella and tendons were included in the most complex model, while they were absent in the simplest model. The greatest difference between the most complex and simplest models was observed in the internal-external rotation and axial joint reaction force, while all other rotations, translations and joint reaction forces were similar to one another. In terms of cartilage stresses and strains, the simpler models behaved similarly with the more complex models in the lateral joint compartment, while minor differences were observed in the medial compartment at the beginning of the stance phase. We suggest that it is feasible to use kinetic-kinematic driven knee joint models with a simpler geometry in studies with a large cohort size, particularly when analyzing cartilage responses and failures related to potential overloads.