The effect of fixed charge density and cartilage swelling on mechanics of knee joint cartilage during simulated gait

Baseline-to-loaded changes in regional tibial cartilage thickness, T1ρ and T2: Utilization of an MRI compatible loading device

The objective of the study was to evaluate tibial cartilage thickness (TCT), T1ρ and T2 values within both loaded and baseline configurations in a cadaveric knee model using a 3D bone based tibial coordinate system. Ten intact cadaveric knees were mounted into an magnetic resonance imaging (MRI) compatible loading device. Morphologic and quantitative MRI (qMRI) images were acquired with the knee in a baseline configuration and after application of 50% body weight. The morphologic images were evaluated for cartilage degeneration using a modified Noyes scoring system. A 3D bone-based tibial coordinate system was utilized to evaluate regional changes of tibial T1ρ, T2, and cartilage thickness values among regions covered and uncovered by the meniscus. Inter-regional differences in medial and lateral MRI outcomes were found between loaded and baseline configurations. Cartilage regions covered by the meniscus demonstrated disparate qMRI and TCT results as compared to cartilage regions not covered by the meniscus. The regions covered by meniscus experienced a ~3.5%, ~0.5%, and ~5.5% reduction of T1ρ (p < 0.05, medial and lateral compartments), T2 and TCT, respectively, in both compartments while regions not covered by the meniscus experienced larger reductions of ~10%, ~2%, and ~10.5% reduction of T1ρ (p < 0.05, medial and lateral compartments), T2 and TCT (p < 0.05, lateral compartment only), respectively, in both compartments. T1ρ and T2 decreases following application of 50% body weight load were substantially larger in the tibial regions with modified Noyes grade 3 (n = 2) compared to either healthy regions (n = 85, p < 0.0.003) or regions with modified Noyes grade 2 (n = 13, p < 0.004). Interregional differences in MRI outcomes reflect variations in structure and function, and largely followed a pattern in cartilage regions that were covered or not covered by the meniscus. Results of the current study suggest that ΔT1ρ and ΔT2 values may be sensitive to superficial fissuring, more than baseline or loaded T1ρ or T2 values, or TCT alone, however future studies with additional specimens, with greater variability in OA grade distribution, may further emphasize the current findings.

Towards a load bearing hydrogel: A proof of principle in the use of osmotic pressure for biomimetic cartilage constructs

Cartilage defects occur frequently and can lead to osteoarthritis. Hydrogels are a promising regenerative strategy for treating such defects, using their ability of mimicking the native extracellular matrix. However, commonly used hydrogels for tissue regeneration are too soft to resist load-bearing in the joint. To overcome this, an implant is being developed in which the mechanical loadbearing function originates from the osmotic pressure generated by the swelling potential of a charged hydrogel, which is restricted from swelling by a textile spacer fabric. This study aims to quantify the relationship between the swelling potential of the hydrogel and the compressive stiffness of the implant.

Creating a Functional Biomimetic Cartilage Implant Using Hydrogels Based on Methacrylated Chondroitin Sulfate and Hyaluronic Acid

The load-bearing function of articular cartilage tissue contrasts with the poor load-bearing capacity of most soft hydrogels used for its regeneration. The present study explores whether a hydrogel based on the methacrylated natural polymers chondroitin sulfate (CSMA) and hyaluronic acid (HAMA), injected into warp-knitted spacer fabrics, could be used to create a biomimetic construct with cartilage-like mechanical properties. The swelling ratio of the combined CSMA/HAMA hydrogels in the first 20 days was higher for hydrogels with a higher CSMA concentration, and these hydrogels also degraded quicker, whereas those with a 1.33 wt% of HAMA were stable for more than 120 days. When confined by a polyamide 6 (PA6) spacer fabric, the volumetric swelling of the combined CSMA/HAMA gels (10 wt%, 6.5 × CSMA:HAMA ratio) was reduced by ~53%. Both the apparent peak and the equilibrium modulus significantly increased in the PA6-restricted constructs compared to the free-swelling hydrogels after 28 days of swelling, and no significant differences in the moduli and time constant compared to native bovine cartilage were observed. Moreover, the cell viability in the CSMA/HAMA PA6 constructs was comparable to that in gelatin–methacrylamide (GelMA) PA6 constructs at one day after polymerization. These results suggest that using a HydroSpacer construct with an extracellular matrix (ECM)-like biopolymer-based hydrogel is a promising approach for mimicking the load-bearing properties of native cartilage.

An EMG-assisted muscle-force driven finite element analysis pipeline to investigate joint- and tissue-level mechanical responses in functional activities: towards a rapid assessment toolbox

Joint tissue mechanics (e.g., stress and strain) are believed to have a major involvement in the onset and progression of musculoskeletal disorders, e.g., knee osteoarthritis (KOA). Accordingly, considerable efforts have been made to develop musculoskeletal finite element (MS-FE) models to estimate highly detailed tissue mechanics that predict cartilage degeneration. However, creating such models is time-consuming and requires advanced expertise. This limits these complex, yet promising MS-FE models to research applications with few participants and makes the models impractical for clinical assessments. Also, these previously developed MS-FE models have not been used to assess activities other than gait. This study introduces and verifies a semi-automated rapid state-of-the-art MS-FE modeling and simulation toolbox incorporating an electromyography- (EMG) assisted MS model and a muscle-force driven FE model of the knee with fibril-reinforced poro(visco)elastic cartilages and menisci. To showcase the usability of the pipeline, we estimated joint- and tissue-level knee mechanics in 15 KOA individuals performing different daily activities. The pipeline was verified by comparing the estimated muscle activations and joint mechanics to existing experimental data. To determine the importance of EMG-assisted MS approach, results were compared to those from the same FE models but driven by static-optimization-based MS models. The EMG-assisted MS-FE pipeline bore a closer resemblance to experiments compared to the static-optimization-based MS-FE pipeline. Importantly, the developed pipeline showed great potential as a rapid MS-FE analysis toolbox to investigate multiscale knee mechanics during different activities of individuals with KOA. The template FE model of the study is freely available here.

Prediction of local fixed charge density loss in cartilage following ACL injury and reconstruction: A computational proof-of-concept study with MRI follow-up

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₂ 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₂ 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.

Computational evaluation of altered biomechanics related to articular cartilage lesions observed in vivo

Chondral lesions provide a potential risk factor for development of osteoarthritis. Despite the variety of in vitro studies on lesion degeneration, in vivo studies that evaluate relation between lesion characteristics and the risk for the possible progression of OA are lacking. Here, we aimed to characterize different lesions and quantify biomechanical responses experienced by surrounding cartilage tissue. We generated computational knee joint models with nine chondral injuries based on clinical in vivo arthrographic computed tomography images. Finite element models with fibril-reinforced poro(visco)elastic cartilage and menisci were constructed to simulate physiological loading. Systematically, the lesions experienced increased peak values of maximum principal strain, maximum shear strain, and minimum principal strain in the surrounding chondral tissue (p < 0.01) compared with intact tissue. Depth, volume, and area of the lesion correlated with the maximum shear strain (p < 0.05, Spearman rank correlation coefficient ρ = 0.733-0.917). Depth and volume of the lesion correlated also with the maximum principal strain (p < 0.05, ρ = 0.767, and ρ = 0.717, respectively). However, the lesion area had non-significant correlation with this strain parameter (p = 0.06, ρ = 0.65). Potentially, the introduced approach could be developed for clinical evaluation of biomechanical risks of a chondral lesion and planning an intervention. Statement of Clinical Relevance: In this study, we computationally characterized different in vivo chondral lesions and evaluated their risk of cartilage degeneration. This information is vital in decision-making for intervention in order to prevent post-traumatic osteoarthritis. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

Tibiofemoral Cartilage Contact Differences Between Level Walking and Downhill Running

Background:

Some studies have suggested that altered tibiofemoral cartilage contact behavior (arthrokinematics) may contribute to long-term cartilage degeneration, potentially leading to tibiofemoral osteoarthritis. However, few studies have assessed normal tibiofemoral arthrokinematics during dynamic activities.

Purpose:

To characterize tibiofemoral arthrokinematics during the impact phase of level walking and downhill running.

Study Design:

Descriptive laboratory study.

Methods:

Arthrokinematic data were collected on uninjured knees of 44 participants (mean age, 20.7 ± 6.6 years). Using a dynamic stereoradiographic imaging system with superimposed 3-dimensional bone models from computed tomography and magnetic resonance imaging of participant-specific tibiofemoral joints, arthrokinematics were assessed during the first 15% of the gait cycle during level walking and the first 10% of the gait cycle during downhill running.

Results:

During level walking and downhill running, the medial compartment had a greater cartilage contact area versus the lateral compartment. Both compartments had a significantly less cartilage contact area during running versus walking (medial compartment gait cycle affected: 8%-10%; lateral compartment gait cycle affected: 5%-10%). Further, medial and lateral compartment tibiofemoral contact paths were significantly more posterior and longer during downhill running.

Conclusion:

There was a decreased tibiofemoral cartilage contact area during downhill running compared with level walking, suggesting that underlying bone morphology may play a key role in determining the size of cartilage contact regions.

Clinical Relevance:

This study provides the first data characterizing tibiofemoral cartilage contact patterns during level walking and downhill running. These results provide evidence in support of performing biomechanical assessments during both level walking and downhill running to obtain a comprehensive picture of tibiofemoral cartilage behavior after clinical interventions.

The ionic contribution of proteoglycans to mechanical stiffness of the meniscus

Load transmission is an important function of the meniscus. In articular cartilage, proteoglycans help maintain hydration via negatively charged moieties. We aimed to investigate the influence of electrostatic effects on stiffness of meniscal tissue. Circular discs from bovine menisci of 8 mm diameter and 5 mm thickness were placed within a confined compression chamber. The apparatus was bathed in distilled water, 0.14 M PBS or 3 M PBS before being subjected to 5% ramp compressive strain and held for 300s. FEBio software was used to fit resultant relaxation curves to a non-linear poroviscoelastic model with strain dependent Holmes-Mow permeability. Analysis was conducted using one-way ANOVA with Tukey post-hoc analysis. 10 samples were tested in each solution. Significant differences (p < 0.05) were observed between the values for Young's modulus, zero strain dependent permeability and the viscoelastic coefficient for samples tested in 3 M PBS as compared to deionised water/0.14 M PBS. No significant differences were observed in the strain dependent/stiffening coefficients or the relaxation time. Approximately 79% of the stiffness of the meniscus appears attributable to ionic effects. Ionic effects play a significant role in the mechanical stiffness of the meniscus. It is important to include the influence of ionic effects when developing mathematical models of this tissue.

On incorporating osmotic prestretch/prestress in image-driven finite element simulations of cartilage

Medical imaging performed in vivo captures geometries under Donnan osmotic loading, even when the articulating joint is otherwise mechanically unloaded. Hence patient-specific finite element (FE) models constructed from such medical images of cartilage represent osmotically induced prestretched/prestressed states. When applying classical modeling approaches to patient-specific simulations of cartilage a theoretical inconsistency arises: the in-vivo imaged geometry (used to construct the model) is not an unloaded, stress-free reference configuration. Furthermore when fitting nonlinear constitutive models that include osmotic swelling (to obtain material parameters), if one assumes that experimental data-generated from osmotically loaded cartilage-begin from a stress-free reference configuration the fitted stress-stretch relationship (parameters) obtained will actually describe a different behavior. In this study we: (1) establish a practical computational method to include osmotically induced prestretch/prestress in image-driven simulations of cartilage; and (2) investigate the influence of considering the prestretched/prestressed state both when fitting fiber-reinforced, biphasic constitutive models of cartilage that include osmotic swelling and when simulating cartilage responses. Our results highlight the importance of determining the prestretched/prestressed state within cartilage induced by osmotic loading in the imaged configuration prior to solving boundary value problems of interest. With our new constitutive model and modeling methods, we aim to improve the fidelity of FE-based, patient-specific biomechanical simulations of joints and cartilage. Improved simulations can provide medical researchers with new information often unavailable in a clinical setting, information that may contribute to better insight into the pathophysiology of cartilage diseases.