Contribution of collagen degradation and proteoglycan depletion to cartilage degeneration in primary and secondary osteoarthritis: an in silico study.

OBJECTIVES: Current experimental approaches cannot elucidate the effect of maladaptive changes on the main cartilage constituents during the degeneration process in osteoarthritis (OA). In silico approaches, however, allow creating 'virtual knock-out' cases to elucidate these effects in a constituent-specific manner. We used such an approach to study the main mechanisms of cartilage degeneration in different mechanical loadings associated with the following OA etiologies: (1) physiological loading of degenerated cartilage, (2) injurious loading of healthy intact cartilage and (3) physiological loading of cartilage with a focal defect. METHODS: We used the recently developed Cartilage Adaptive REorientation Degeneration (CARED) framework to simulate cartilage degeneration associated with primary and secondary OA (OA cases (1)-(3)). CARED incorporates numerical description of tissue-level cartilage degeneration mechanisms in OA, namely, collagen degradation, collagen reorientation, fixed charged density loss and tissue hydration increase following mechanical loading. We created 'virtual knock-out' scenarios by deactivating these degenerative processes one at a time in each of the three OA cases. RESULTS: In the injurious loading of intact and physiological loading of degenerated cartilage, collagen degradation drives degenerative changes through fixed charge density loss and tissue hydration rise. In contrast, the two later mechanisms were more prominent in the focal defect cartilage model. CONCLUSION: The virtual knock-out models reveal that injurious loading to intact cartilage and physiological loading to degenerated cartilage induce initial degenerative changes in the collagen network, whereas, in the presence of a focal cartilage defect, mechanical loading initially causes proteoglycans (PG) depletion, before changes in the collagen fibril network occur.

Comparison of material models for anterior cruciate ligament in tension: from poroelastic to a novel fibril-reinforced nonlinear composite model.

Computational models of the knee joint are useful for evaluating stresses and strains within the joint tissues. However, the outcome of those models is sensitive to the material model and material properties chosen for ligaments, the collagen reinforced tissues connecting bone to bone. The purpose of this study was to investigate different compositionally motivated material models and further to develop a model that can accurately reproduce experimentally measured stress-relaxation data of bovine anterior cruciate ligament (ACL). Tensile testing samples were extracted from ACLs of bovine knee joints (N = 10) and subjected to a three-step stress-relaxation test at the toe region. Data from the experiments was averaged and one average finite element model was generated to replicate the experiment. Poroelastic and different fibril-reinforced poro(visco)elastic material models were applied, and their material parameters were optimized to reproduce the experimental force-time response. Material models with only fluid flow mediated relaxation were not able to capture the stress-relaxation behavior (R2 = 0.806, 0.803 and 0.938). The inclusion of the viscoelasticity of the fibrillar network improved the model prediction (R2 = 0.978 and 0.976), but the complex stress-relaxation behavior was best captured by a poroelastic model with a nonlinear two-relaxation-time strain-recruited viscoelastic fibrillar network (R2 = 0.997). The results suggest that in order to replicate the multi-step stress-relaxation behavior of ACL in tension, the fibrillar network formulation should include the complex nonlinear viscoelastic phenomena.

12 Degrees of Freedom Muscle Force Driven Fibril-Reinforced Poroviscoelastic Finite Element Model of the Knee Joint.

Accurate knowledge of the joint kinematics, kinetics, and soft tissue mechanical responses is essential in the evaluation of musculoskeletal (MS) disorders. Since in vivo measurement of these quantities requires invasive methods, musculoskeletal finite element (MSFE) models are widely used for simulations. There are, however, limitations in the current approaches. Sequentially linked MSFE models benefit from complex MS and FE models; however, MS model's outputs are independent of the FE model calculations. On the other hand, due to the computational burden, embedded (concurrent) MSFE models are limited to simple material models and cannot estimate detailed responses of the soft tissue. Thus, first we developed a MSFE model of the knee with a subject-specific MS model utilizing an embedded 12 degrees of freedom (DoFs) knee joint with elastic cartilages in which included both secondary kinematic and soft tissue deformations in the muscle force estimation (inverse dynamics). Then, a muscle-force-driven FE model with fibril-reinforced poroviscoelastic cartilages and fibril-reinforced poroelastic menisci was used in series to calculate detailed tissue mechanical responses (forward dynamics). Second, to demonstrate that our workflow improves the simulation results, outputs were compared to results from the same FE models which were driven by conventional MS models with a 1 DoF knee, with and without electromyography (EMG) assistance. The FE model driven by both the embedded and the EMG-assisted MS models estimated similar results and consistent with experiments from literature, compared to the results estimated by the FE model driven by the MS model with 1 DoF knee without EMG assistance.

EMG-Assisted Muscle Force Driven Finite Element Model of the Knee Joint with Fibril-Reinforced Poroelastic Cartilages and Menisci.

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.

Interrelationship of cartilage composition and chondrocyte mechanics after a partial meniscectomy in the rabbit knee joint - Experimental and numerical analysis.

Site-specific and depth-dependent properties of cartilage were implemented within a finite element (FE) model to determine if compositional or structural changes in the tissue could explain site-specific alterations of chondrocyte deformations due to cartilage loading in rabbit knee joints 3 days after a partial meniscectomy (PM). Depth-dependent proteoglycan (PG) content, collagen content and collagen orientation in the cartilage extracellular matrix (ECM), and PG content in the pericellular matrix (PCM) were assessed with microscopic and spectroscopic methods. Patellar, femoral groove and samples from both the lateral and medial compartments of the femoral condyle and tibial plateau were extracted from healthy controls and from the partial meniscectomy group. For both groups and each knee joint site, axisymmetric FE models with measured properties were generated. Experimental cartilage loading was applied in the simulations and chondrocyte volumes were compared to the experimental values. ECM and PCM PG loss occurred within the superficial cartilage layer in the PM group at all locations, except in the lateral tibial plateau. Collagen content and orientation were not significantly altered due to the PM. The FE simulations predicted similar chondrocyte volume changes and group differences as obtained experimentally. Loss of PCM fixed charge density (FCD) decreased cell volume loss, as observed in the medial femur and medial tibia, whereas loss of ECM FCD increased cell volume loss, as seen in the patella, femoral groove and lateral femur. The model outcome, cell volume change, was also sensitive to applied tissue geometry, collagen fibril orientation and loading conditions.

Tissue viscoelasticity is related to tissue composition but may not fully predict the apparent-level viscoelasticity in human trabecular bone - An experimental and finite element study.

Trabecular bone is viscoelastic under dynamic loading. However, it is unclear how tissue viscoelasticity controls viscoelasticity at the apparent-level. In this study, viscoelasticity of cylindrical human trabecular bone samples (n=11, male, age 18-78 years) from 11 proximal femurs were characterized using dynamic and stress-relaxation testing at the apparent-level and with creep nanoindentation at the tissue-level. In addition, bone tissue elasticity was determined using scanning acoustic microscope (SAM). Tissue composition and collagen crosslinks were assessed using Raman micro-spectroscopy and high performance liquid chromatography (HPLC), respectively. Values of material parameters were obtained from finite element (FE) models by optimizing tissue-level creep and apparent-level stress-relaxation to experimental nanoindentation and unconfined compression testing values, respectively, utilizing the second order Prony series to depict viscoelasticity. FE simulations showed that tissue-level equilibrium elastic modulus (Eeq) increased with increasing crystallinity (r=0.730, p=.011) while at the apparent-level it increased with increasing hydroxylysyl pyridinoline content (r=0.718, p=.019). In addition, the normalized shear modulus g1 (r=-0.780, p=.005) decreased with increasing collagen ratio (amide III/CH2) at the tissue-level, but increased (r=0.696, p=.025) with increasing collagen ratio at the apparent-level. No significant relations were found between the measured or simulated viscoelastic parameters at the tissue- and apparent-levels nor were the parameters related to tissue elasticity determined with SAM. However, only Eeq, g2 and relaxation time τ1 from simulated viscoelastic values were statistically different between tissue- and apparent-levels (p<.01). These findings indicate that bone tissue viscoelasticity is affected by tissue composition but may not fully predict the macroscale viscoelasticity in human trabecular bone.

Characterization of site-specific biomechanical properties of human meniscus-Importance of collagen and fluid on mechanical nonlinearities.

Meniscus adapts to joint loads by depth- and site-specific variations in its composition and structure. However, site-specific mechanical characteristics of intact meniscus under compression are poorly known. In particular, mechanical nonlinearities caused by different meniscal constituents (collagen and fluid) are not known. In the current study, in situ indentation testing was conducted to determine site-specific elastic, viscoelastic and poroelastic properties of intact human menisci. Lateral and medial menisci (n=26) were harvested from the left knee joint of 13 human cadavers. Indentation tests, using stress-relaxation and dynamic (sinusoidal) loading protocols, were conducted for menisci at different sites (anterior, middle, posterior, n=78). Sample- and site-specific axisymmetric finite element models with fibril-reinforced poroelastic properties were fitted to the corresponding stress-relaxation curves to determine the mechanical parameters. Elastic moduli, especially the instantaneous and dynamic moduli, showed site-specific variation only in the medial meniscus (p<0.05 between the sites). The instantaneous and dynamic elastic moduli of the anterior horn were significantly (p<0.05) greater in the medial than lateral meniscus. The phase angle showed no statistically significant variation between the sites (p>0.05). The values for the strain-dependent fibril network modulus (nonlinear behaviour of collagen) were significantly different (p<0.05) between all sites in the medial menisci. Additionally, there was a significant difference (p<0.01) in the strain-dependent fibril network modulus between the lateral and medial anterior horns. The initial permeability was significantly different (p<0.05) in the medial meniscus only between the middle and posterior sites. For the strain-dependent permeability coefficient, only anterior and middle sites showed a significant difference (p<0.05) in the medial meniscus. This parameter demonstrated a significant difference (p<0.05) between lateral and medial menisci at the anterior horns. Our results reveal that under in situ indentation loading, medial meniscus shows more site-dependent variation in the mechanical properties as compared to lateral meniscus. In particular, anterior horn of medial meniscus was the stiffest and showed the most nonlinear mechanical behaviour. The nonlinearity was related to both collagen fibrils and fluid.

The effect of collagen degradation on chondrocyte volume and morphology in bovine articular cartilage following a hypotonic challenge.

Collagen degradation is one of the early signs of osteoarthritis. It is not known how collagen degradation affects chondrocyte volume and morphology. Thus, the aim of this study was to investigate the effect of enzymatically induced collagen degradation on cell volume and shape changes in articular cartilage after a hypotonic challenge. Confocal laser scanning microscopy was used for imaging superficial zone chondrocytes in intact and degraded cartilage exposed to a hypotonic challenge. Fourier transform infrared microspectroscopy, polarized light microscopy, and mechanical testing were used to quantify differences in proteoglycan and collagen content, collagen orientation, and biomechanical properties, respectively, between the intact and degraded cartilage. Collagen content decreased and collagen orientation angle increased significantly (p < 0.05) in the superficial zone cartilage after collagenase treatment, and the instantaneous modulus of the samples was reduced significantly (p < 0.05). Normalized cell volume and height 20 min after the osmotic challenge (with respect to the original volume and height) were significantly (p < 0.001 and p < 0.01, respectively) larger in the intact compared to the degraded cartilage. These findings suggest that the mechanical environment of chondrocytes, specifically collagen content and orientation, affects cell volume and shape changes in the superficial zone articular cartilage when exposed to osmotic loading. This emphasizes the role of collagen in modulating cartilage mechanobiology in diseased tissue.