Three dimensional patient-specific collagen architecture modulates cartilage responses in the knee joint during gait

Revealing Detailed Cartilage Function Through Nanoparticle Diffusion Imaging: A Computed Tomography & Finite Element Study

The ability of articular cartilage to withstand significant mechanical stresses during activities, such as walking or running, relies on its distinctive structure. Integrating detailed tissue properties into subject-specific biomechanical models is challenging due to the complexity of analyzing these characteristics. This limitation compromises the accuracy of models in replicating cartilage function and impacts predictive capabilities. To address this, methods revealing cartilage function at the constituent-specific level are essential. In this study, we demonstrated that computational modeling derived individual constituent-specific biomechanical properties could be predicted by a novel nanoparticle contrast-enhanced computer tomography (CECT) method. We imaged articular cartilage samples collected from the equine stifle joint (n = 60) using contrast-enhanced micro-computed tomography (µCECT) to determine contrast agents' intake within the samples, and compared those to cartilage functional properties, derived from a fibril-reinforced poroelastic finite element model. Two distinct imaging techniques were investigated: conventional energy-integrating µCECT employing a cationic tantalum oxide nanoparticle (Ta2O5-cNP) contrast agent and novel photon-counting µCECT utilizing a dual-contrast agent, comprising Ta2O5-cNP and neutral iodixanol. The results demonstrate the capacity to evaluate fibrillar and non-fibrillar functionality of cartilage, along with permeability-affected fluid flow in cartilage. This finding indicates the feasibility of incorporating these specific functional properties into biomechanical computational models, holding potential for personalized approaches to cartilage diagnostics and treatment.

Triple contrast computed tomography reveals site-specific biomechanical differences in the human knee joint-A proof of concept study

Cartilage and synovial fluid are challenging to observe separately in native computed tomography (CT). We report the use of triple contrast agent (bismuth nanoparticles [BiNPs], CA4+, and gadoteridol) to image and segment cartilage in cadaveric knee joints with a clinical CT scanner. We hypothesize that BiNPs will remain in synovial fluid while the CA4+ and gadoteridol will diffuse into cartilage, allowing (1) segmentation of cartilage, and (2) evaluation of cartilage biomechanical properties based on contrast agent concentrations. To investigate these hypotheses, triple contrast agent was injected into both knee joints of a cadaver (N = 1), imaged with a clinical CT at multiple timepoints during the contrast agent diffusion. Knee joints were extracted, imaged with micro-CT (µCT), and biomechanical properties of the cartilage surface were determined by stress-relaxation mapping. Cartilage was segmented and contrast agent concentrations (CA4+ and gadoteridol) were compared with the biomechanical properties at multiple locations (n = 185). Spearman's correlation between cartilage thickness from clinical CT and reference µCT images verifies successful and reliable segmentation. CA4+ concentration is significantly higher in femoral than in tibial cartilage at 60 min and further timepoints, which corresponds to the higher Young's modulus observed in femoral cartilage. In this pilot study, we show that (1) large BiNPs do not diffuse into cartilage, facilitating straightforward segmentation of human knee joint cartilage in a clinical setting, and (2) CA4+ concentration in cartilage reflects the biomechanical differences between femoral and tibial cartilage. Thus, the triple contrast agent CT shows potential in cartilage morphology and condition estimation in clinical CT.

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 T₂ and T_1ρ 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 T₂ and T_1ρ 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 T₂ relaxation times. There was also a significant correlation between the predicted stress values and changes in T₂ relaxation time (r = 0.649, p < 0.001). Absolute maximum shear strains and changes in T_1ρ 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.

Towards Tailored Rehabilitation by Implementation of a Novel Musculoskeletal Finite Element Analysis Pipeline

Tissue-level mechanics (e.g., stress and strain) are important factors governing tissue remodeling and development of knee osteoarthritis (KOA), and hence, the success of physical rehabilitation. To date, no clinically feasible analysis toolbox has been introduced and used to inform clinical decision making with subject-specific in-depth joint mechanics of different activities. Herein, we utilized a rapid state-of-the-art electromyography-assisted musculoskeletal finite element analysis toolbox with fibril-reinforced poro(visco)elastic cartilages and menisci to investigate knee mechanics in different activities. Tissue mechanical responses, believed to govern collagen damage, cell death, and fixed charge density loss of proteoglycans, were characterized within 15 patients with KOA while various daily activities and rehabilitation exercises were performed. Results showed more inter-participant variation in joint mechanics during rehabilitation exercises compared to daily activities. Accordingly, the devised workflow may be used for designing subject-specific rehabilitation protocols. Further, results showed the potential to tailor rehabilitation exercises, or assess capacity for daily activity modifications, to optimally load knee tissue, especially when mechanically-induced cartilage degeneration and adaptation are of interest.

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.

Rapid CT-based Estimation of Articular Cartilage Biomechanics in the Knee Joint Without Cartilage Segmentation

Knee osteoarthritis (OA) is a painful joint disease, causing disabilities in daily activities. However, there is no known cure for OA, and the best treatment strategy might be prevention. Finite element (FE) modeling has demonstrated potential for evaluating personalized risks for the progression of OA. Current FE modeling approaches use primarily magnetic resonance imaging (MRI) to construct personalized knee joint models. However, MRI is expensive and has lower resolution than computed tomography (CT). In this study, we extend a previously presented atlas-based FE modeling framework for automatic model generation and simulation of knee joint tissue responses using contrast agent-free CT. In this method, based on certain anatomical dimensions measured from bone surfaces, an optimal template is selected and scaled to generate a personalized FE model. We compared the simulated tissue responses of the CT-based models with those of the MRI-based models. We show that the CT-based models are capable of producing similar tensile stresses, fibril strains, and fluid pressures of knee joint cartilage compared to those of the MRI-based models. This study provides a new methodology for the analysis of knee joint and cartilage mechanics based on measurement of bone dimensions from native CT scans.

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.

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.

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.

White Matter Anisotropy for Impact Simulation and Response Sampling in Traumatic Brain Injury

Advanced neuroimaging provides new opportunities to enhance head injury models, including the incorporation of white matter (WM) structural anisotropy. Information from high-resolution neuroimaging, however, usually has to be "down-sampled" to match a typically coarse brain mesh. To understand how this mesh-image resolution mismatch affects impact simulation and subsequent response sampling, we compared three competing anisotropy implementations (using either voxels, tractography, or a multiscale submodeling) and two response sampling strategies (element-wise or tractography-based, using brain mesh or neuroimaging for region segmentation, respectively). Using the combination of high resolution options as a baseline, we studied how the choice in each individual category affected the resulting injury metrics. By simulating a recorded loss of consciousness head impact, we found that injury metrics including peak strain and injury susceptibility in the deep WM regions based on fiber strain, but not on maximum principal strain, were sensitive to the anisotropy implementation, response sampling, and region segmentation. Overall, it was recommended to use tractography for anisotropy implementation and response sampling, and to employ neuroimaging for region segmentation, because they led to more accurate injury metrics. Further refining mesh locally via submodeling was unnecessary. Brain strain responses were also parametrically found to be closer to that from minimum fiber reinforcement, consistent with the fact that the majority of WM had a rather high degree of fiber dispersion. Finally, the upgraded Worcester Head Injury Model incorporating WM anisotropy was successfully re-validated against cadaveric impacts and an in vivo head rotation ("good" to "excellent" validation with an average Correlation Analysis score of 0.437 and 0.509, respectively). These investigations may facilitate further continual development of head injury models to better study traumatic brain injury.

The effects of geometric uncertainties on computational modelling of knee biomechanics

The geometry of the articular components of the knee is an important factor in predicting joint mechanics in computational models. There are a number of uncertainties in the definition of the geometry of cartilage and meniscus, and evaluating the effects of these uncertainties is fundamental to understanding the level of reliability of the models. In this study, the sensitivity of knee mechanics to geometric uncertainties was investigated by comparing polynomial-based and image-based knee models and varying the size of meniscus. The results suggested that the geometric uncertainties in cartilage and meniscus resulting from the resolution of MRI and the accuracy of segmentation caused considerable effects on the predicted knee mechanics. Moreover, even if the mathematical geometric descriptors can be very close to the imaged-based articular surfaces, the detailed contact pressure distribution produced by the mathematical geometric descriptors was not the same as that of the image-based model. However, the trends predicted by the models based on mathematical geometric descriptors were similar to those of the imaged-based models.

Quantitative Evaluation of the Mechanical Risks Caused by Focal Cartilage Defects in the Knee

Focal cartilage lesions can proceed to severe osteoarthritis or remain unaltered even for years. A method to identify high risk defects would be of utmost importance to guide clinical decision making and to identify the patients that are at the highest risk for the onset and progression of osteoarthritis. Based on cone beam computed tomography arthrography, we present a novel computational model for evaluating changes in local mechanical responses around cartilage defects. Our model, based on data obtained from a human knee in vivo, demonstrated that the most substantial alterations around the defect, as compared to the intact tissue, were observed in minimum principal (compressive) strains and shear strains. Both strain values experienced up to 3-fold increase, exceeding levels previously associated with chondrocyte apoptosis and failure of collagen crosslinks. Furthermore, defects at the central regions of medial tibial cartilage with direct cartilage-cartilage contact were the most vulnerable to loading. Also locations under the meniscus experienced substantially increased minimum principal strains. We suggest that during knee joint loading particularly minimum principal and shear strains are increased above tissue failure limits around cartilage defects which might lead to osteoarthritis. However, this increase in strains is highly location-specific on the joint surface.