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.

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.

Ligaments provide stability to the human knee joint and play an essential role in restraining motion during daily activities. Compression-tension nonlinearity is a well-known characteristic of ligaments. Moreover, simpler material representations without this feature might give reasonable results because ligaments are primarily in tension during loading. However, the biomechanical role of different constitutive representations and their fibril-reinforced poroelastic properties is unknown. A numerical knee model which considers geometric and material nonlinearities of meniscus and cartilages was applied. Five different constitutive models for the ligaments (spring, elastic, hyperelastic, porohyperelastic, and fibril-reinforced porohyperelastic (FRPHE)) were implemented. Knee joint forces for the models with elastic, hyperelastic and porohyperelastic properties showed similar behavior throughout the stance, while the model with FRPHE properties exhibited lower joint forces during the last 50% of the stance phase. The model with ligaments as springs produced the lowest joint forces at this same stance phase. The results also showed that the fibril network contributed substantially to the knee joint forces, while the nonfibrillar matrix and fluid had small effects. Our results indicate that simpler material models of ligaments with similar properties in compression and tension can be used when the loading is directed primarily along the ligament axis in tension.

Post-traumatic osteoarthritis (PTOA) is associated with cartilage degradation, ultimately leading to disability and decrease of quality of life. Two key mechanisms have been suggested to occur in PTOA: tissue inflammation and abnormal biomechanical loading. Both mechanisms have been suggested to result in loss of cartilage proteoglycans, the source of tissue fixed charge density (FCD). In order to predict the simultaneous effect of these degrading mechanisms on FCD content, a computational model has been developed. We simulated spatial and temporal changes of FCD content in injured cartilage using a novel finite element model that incorporates (1) diffusion of the pro-inflammatory cytokine interleukin-1 into tissue, and (2) the effect of excessive levels of shear strain near chondral defects during physiologically relevant loading. Cytokine-induced biochemical cartilage explant degradation occurs near the sides, top, and lesion, consistent with the literature. In turn, biomechanically-driven FCD loss is predicted near the lesion, in accordance with experimental findings: regions near lesions showed significantly more FCD depletion compared to regions away from lesions (p<0.01). Combined biochemical and biomechanical degradation is found near the free surfaces and especially near the lesion, and the corresponding bulk FCD loss agrees with experiments. We suggest that the presence of lesions plays a role in cytokine diffusion-driven degradation, and also predisposes cartilage for further biomechanical degradation. Models considering both these cartilage degradation pathways concomitantly are promising in silico tools for predicting disease progression, recognizing lesions at high risk, simulating treatments, and ultimately optimizing treatments to postpone the development of PTOA.