Cartilage provides low-friction properties and plays an essential role in diarthrodial joints. A hydrated ground substance composed mainly of proteoglycans (PGs) and a fibrillar collagen network are the main constituents of cartilage. Unfortunately, traumatic joint loading can destroy this complex structure and produce lesions in tissue, leading later to changes in tissue composition and, ultimately, to post-traumatic osteoarthritis (PTOA). Consequently, the fixed charge density (FCD) of PGs may decrease near the lesion. However, the underlying mechanisms leading to these tissue changes are unknown. Here, knee cartilage disks from bovine calves were injuriously compressed, followed by a physiologically relevant dynamic compression for twelve days. FCD content at different follow-up time points was assessed using digital densitometry. A novel cartilage degeneration model was developed by implementing deviatoric and maximum shear strain, as well as fluid velocity controlled algorithms to simulate the FCD loss as a function of time. Predicted loss of FCD was quite uniform around the cartilage lesions when the degeneration algorithm was driven by the fluid velocity, while the deviatoric and shear strain driven mechanisms exhibited slightly discontinuous FCD loss around cracks. Our degeneration algorithm predictions fitted well with the FCD content measured from the experiments. The developed model could subsequently be applied for prediction of FCD depletion around different cartilage lesions and for suggesting optimal rehabilitation protocols.

We present a novel algorithm combined with computational modeling to simulate the development of knee osteoarthritis. The degeneration algorithm was based on excessive and cumulatively accumulated stresses within knee joint cartilage during physiological gait loading. In the algorithm, the collagen network stiffness of cartilage was reduced iteratively if excessive maximum principal stresses were observed. The developed algorithm was tested and validated against experimental baseline and 4-year follow-up Kellgren-Lawrence grades, indicating different levels of cartilage degeneration at the tibiofemoral contact region. Test groups consisted of normal weight and obese subjects with the same gender and similar age and height without osteoarthritic changes. The algorithm accurately simulated cartilage degeneration as compared to the Kellgren-Lawrence findings in the subject group with excess weight, while the healthy subject group’s joint remained intact. Furthermore, the developed algorithm followed the experimentally found trend of cartilage degeneration in the obese group (R2 = 0.95, p < 0.05; experiments vs. model), in which the rapid degeneration immediately after initiation of osteoarthritis (0–2 years, p < 0.001) was followed by a slow or negligible degeneration (2–4 years, p > 0.05). The proposed algorithm revealed a great potential to objectively simulate the progression of knee osteoarthritis.

The purpose of this proof-of-concept study was to develop three-dimensional patient-specific mechanobiological knee joint models to simulate alterations in the fixed charged density (FCD) around cartilage lesions during the stance phase of the walking gait. Two patients with anterior cruciate ligament (ACL) reconstructed knees were imaged at 1 and 3 years after surgery. The magnetic resonance imaging (MRI) data were used for segmenting the knee geometries, including the cartilage lesions. Based on these geometries, finite element (FE) models were developed. The gait of the patients was obtained using a motion capture system. Musculoskeletal modeling was utilized to calculate knee joint contact and lower extremity muscle forces for the FE models. Finally, a cartilage adaptation algorithm was implemented in both FE models. In the algorithm, it was assumed that excessive maximum shear and deviatoric strains (calculated as the combination of principal strains), and fluid velocity, are responsible for the FCD loss. Changes in the longitudinal T 1ρ and T 2 relaxation times were postulated to be related to changes in the cartilage composition and were compared with the numerical predictions. In patient 1 model, both the excessive fluid velocity and strain caused the FCD loss primarily near the cartilage lesion. T 1ρ and T 2 relaxation times increased during the follow-up in the same location. In contrast, in patient 2 model, only the excessive fluid velocity led to a slight FCD loss near the lesion, where MRI parameters did not show evidence of alterations. Significance: This novel proof-of-concept study suggests mechanisms through which a local FCD loss might occur near cartilage lesions. In order to obtain statistical evidence for these findings, the method should be investigated with a larger cohort of subjects.

No abstract

Abnormal mechanical loading is essential in the onset and progression of knee osteoarthritis. Combined musculoskeletal (MS) and finite element (FE) modeling is a typical method to estimate load distribution and tissue responses in the knee joint. However, earlier combined models mostly utilize static-optimization based MS models and muscle force driven FE models typically use elastic materials for soft tissues or analyze specific time points of gait. Therefore, here we develop an electromyography-assisted muscle force driven FE model with fibril-reinforced poro(visco)elastic cartilages and menisci to analyze knee joint loading during the stance phase of gait. Moreover, since ligament pre-strains are one of the important uncertainties in joint modeling, we conducted a sensitivity analysis on the pre-strains of anterior and posterior cruciate ligaments (ACL and PCL) as well as medial and lateral collateral ligaments (MCL and LCL). The model produced kinematics and kinetics consistent with previous experimental data. Joint contact forces and contact areas were highly sensitive to ACL and PCL pre-strains, while those changed less cartilage stresses, fibril strains, and fluid pressures. The presented workflow could be used in a wide range of applications related to the aetiology of cartilage degeneration, optimization of rehabilitation exercises, and simulation of knee surgeries.

The intricate properties of articular cartilage and the complexity of the loading environment are some of the key challenges in developing models for biomechanical analysis of the knee joint. Fibril-reinforced poroelastic (FRPE) material models have been reported to accurately capture characteristic responses of cartilage during dynamic and static loadings. However, high computational and time costs associated with such advanced models limit applicability of FRPE models when multiple subjects need to be analyzed. If choosing simpler material models, it is important to show that they can still produce truthful predictions. Therefore, the aim of this study was to compare depth-dependent maximum principal stresses and strains within articular cartilage in the 3D knee joint between FRPE material models and simpler isotropic elastic (IE), isotropic poroelastic (IPE) and transversely isotropic poroelastic (TIPE) material models during simulated gait cycle. When cartilage-cartilage contact pressures were matched between the models (15% allowed difference), maximum principal stresses in the IE, IPE and TIPE models were substantially lower than those in the FRPE model (by more than 50%, TIPE model being closest to the FRPE model), and stresses occurred only in compression in the IE model. Additional simulations were performed to find material parameters for the TIPE model (due to its anisotropic nature) that would yield maximum principal stresses similar to the FRPE model. The modified homogeneous TIPE model was in a better agreement with the homogeneous FRPE model, and the average and maximum differences in maximum principal stresses throughout the depth of cartilage were 7% and 9%, respectively, in the lateral compartment and 9% and 11% in the medial compartment. This study revealed that it is possible to match simultaneously maximum principal stresses and strains of cartilage between non-fibril-reinforced and fibril-reinforced knee joint models during gait. Depending on the research question (such as analysis of fibril strain necessitates the use of fibril-reinforced material models) or clinical demand (fast simulations with simpler material models), the choice of the material model should be done carefully.

Associate Editor Eric M. Darling oversaw the review of this article.

Post-traumatic osteoarthritis (PTOA) is a common disease, where the mechanical integrity of articular cartilage is compromised. PTOA can be a result of chondral defects formed due to injurious loading. One of the first changes around defects is proteoglycan depletion. Since there are no methods to restore injured cartilage fully back to its healthy state, preventing the onset and progression of the disease is advisable. However, this is problematic if the disease progression cannot be predicted. Thus, we developed an algorithm to predict proteoglycan loss of injured cartilage by decreasing the fixed charge density (FCD) concentration. We tested several mechanisms based on the local strains or stresses in the tissue for the FCD loss. By choosing the degeneration threshold suggested for inducing chondrocyte apoptosis and cartilage matrix damage, the algorithm driven by the maximum shear strain showed the most substantial FCD losses around the lesion. This is consistent with experimental findings in the literature. We also observed that by using coordinate system-independent strain measures and selecting the degeneration threshold in an ad hoc manner, all the resulting FCD distributions would appear qualitatively similar, i.e., the greatest FCD losses are found at the tissue adjacent to the lesion. The proposed strain-based FCD degeneration algorithm shows a great potential for predicting the progression of PTOA via biomechanical stimuli. This could allow identification of high-risk defects with an increased risk of PTOA progression.