1
|
Rapid X-Ray-Based 3-D Finite Element Modeling of Medial Knee Joint Cartilage Biomechanics During Walking. Ann Biomed Eng 2022; 50:666-679. [PMID: 35262835 PMCID: PMC9079039 DOI: 10.1007/s10439-022-02941-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Accepted: 02/23/2022] [Indexed: 11/30/2022]
Abstract
Finite element (FE) modeling is becoming an increasingly popular method for analyzing knee joint mechanics and biomechanical mechanisms leading to osteoarthritis (OA). The most common and widely available imaging method for knee OA diagnostics is planar X-ray imaging, while more sophisticated imaging methods, e.g., magnetic resonance imaging (MRI) and computed tomography (CT), are seldom used. Hence, the capability to produce accurate biomechanical knee joint models directly from X-ray imaging would bring FE modeling closer to clinical use. Here, we extend our atlas-based framework by generating FE knee models from X-ray images (N = 28). Based on measured anatomical landmarks from X-ray and MRI, knee joint templates were selected from the atlas library. The cartilage stresses and strains of the X-ray-based model were then compared with the MRI-based model during the stance phase of the gait. The biomechanical responses were statistically not different between MRI- vs. X-ray-based models when the template obtained from X-ray imaging was the same as the MRI template. However, if this was not the case, the peak values of biomechanical responses were statistically different between X-ray and MRI models. The developed X-ray-based framework may pave the way for a clinically feasible approach for knee joint FE modeling.
Collapse
|
2
|
Mohammadi A, Myller KAH, Tanska P, Hirvasniemi J, Saarakkala S, Töyräs J, Korhonen RK, Mononen ME. Rapid CT-based Estimation of Articular Cartilage Biomechanics in the Knee Joint Without Cartilage Segmentation. Ann Biomed Eng 2020; 48:2965-2975. [PMID: 33179182 PMCID: PMC7723937 DOI: 10.1007/s10439-020-02666-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Accepted: 10/17/2020] [Indexed: 12/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Ali Mohammadi
- Department of Applied Physics, University of Eastern Finland, POB 1627, 70211, Kuopio, Finland.
| | - Katariina A H Myller
- Diagnostic Imaging Center, Kuopio University Hospital, Kuopio, Finland.,Department of Medical Physics, Turku University Central Hospital, 20500, Turku, Finland
| | - Petri Tanska
- Department of Applied Physics, University of Eastern Finland, POB 1627, 70211, Kuopio, Finland
| | - Jukka Hirvasniemi
- Department of Radiology & Nuclear Medicine, Erasmus University Medical Center, Rotterdam, The Netherlands
| | - Simo Saarakkala
- Research Unit of Medical Imaging, Physics and Technology, Faculty of Medicine, University of Oulu, Oulu, Finland.,Department of Diagnostic Radiology, Oulu University Hospital, Oulu, Finland
| | - Juha Töyräs
- Department of Applied Physics, University of Eastern Finland, POB 1627, 70211, Kuopio, Finland.,Diagnostic Imaging Center, Kuopio University Hospital, Kuopio, Finland.,School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia
| | - Rami K Korhonen
- Department of Applied Physics, University of Eastern Finland, POB 1627, 70211, Kuopio, Finland
| | - Mika E Mononen
- Department of Applied Physics, University of Eastern Finland, POB 1627, 70211, Kuopio, Finland
| |
Collapse
|
3
|
Eskelinen ASA, Tanska P, Florea C, Orozco GA, Julkunen P, Grodzinsky AJ, Korhonen RK. Mechanobiological model for simulation of injured cartilage degradation via pro-inflammatory cytokines and mechanical stimulus. PLoS Comput Biol 2020; 16:e1007998. [PMID: 32584809 PMCID: PMC7343184 DOI: 10.1371/journal.pcbi.1007998] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2020] [Revised: 07/08/2020] [Accepted: 05/28/2020] [Indexed: 01/12/2023] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Petri Tanska
- Department of Applied Physics, University of Eastern Finland, Finland
| | - Cristina Florea
- Department of Applied Physics, University of Eastern Finland, Finland
- Departments of Biological Engineering, Electrical Engineering and Computer Science and Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, United States of America
| | - Gustavo A. Orozco
- Department of Applied Physics, University of Eastern Finland, Finland
| | - Petro Julkunen
- Department of Applied Physics, University of Eastern Finland, Finland
- Department of Clinical Neurophysiology, Kuopio University Hospital, Finland
| | - Alan J. Grodzinsky
- Departments of Biological Engineering, Electrical Engineering and Computer Science and Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, United States of America
| | - Rami K. Korhonen
- Department of Applied Physics, University of Eastern Finland, Finland
| |
Collapse
|
4
|
Myller KAH, Korhonen RK, Töyräs J, Tanska P, Väänänen SP, Jurvelin JS, Saarakkala S, Mononen ME. Clinical Contrast-Enhanced Computed Tomography With Semi-Automatic Segmentation Provides Feasible Input for Computational Models of the Knee Joint. J Biomech Eng 2020; 142:051001. [PMID: 31647541 DOI: 10.1115/1.4045279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Indexed: 11/08/2022]
Abstract
Computational models can provide information on joint function and risk of tissue failure related to progression of osteoarthritis (OA). Currently, the joint geometries utilized in modeling are primarily obtained via manual segmentation, which is time-consuming and hence impractical for direct clinical application. The aim of this study was to evaluate the applicability of a previously developed semi-automatic method for segmenting tibial and femoral cartilage to serve as input geometry for finite element (FE) models. Knee joints from seven volunteers were first imaged using a clinical computed tomography (CT) with contrast enhancement and then segmented with semi-automatic and manual methods. In both segmentations, knee joint models with fibril-reinforced poroviscoelastic (FRPVE) properties were generated and the mechanical responses of articular cartilage were computed during physiologically relevant loading. The mean differences in the absolute values of maximum principal stress, maximum principal strain, and fibril strain between the models generated from semi-automatic and manual segmentations were <1 MPa, <0.72% and <0.40%, respectively. Furthermore, contact areas, contact forces, average pore pressures, and average maximum principal strains were not statistically different between the models (p >0.05). This semi-automatic method speeded up the segmentation process by over 90% and there were only negligible differences in the results provided by the models utilizing either manual or semi-automatic segmentations. Thus, the presented CT imaging-based segmentation method represents a novel tool for application in FE modeling in the clinic when a physician needs to evaluate knee joint function.
Collapse
Affiliation(s)
- Katariina A H Myller
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland; Diagnostic Imaging Center, Kuopio University Hospital, P.O. Box 100, Kuopio FI-70029, Finland
| | - Rami K Korhonen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland
| | - Juha Töyräs
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland; Diagnostic Imaging Center, Kuopio University Hospital, P.O. Box 100, Kuopio FI-70029, Finland; School of Information Technology and Electrical Engineering, The University of Queensland, St Lucia Qld, Brisbane 4072, Australia
| | - Petri Tanska
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland
| | - Sami P Väänänen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland; Diagnostic Imaging Center, Kuopio University Hospital, P.O. Box 100, Kuopio FI-70029, Finland; Central Finland Central Hospital, Department of Physics, Keskussairaalantie 19, Jyväskylä FI-40620, Finland
| | - Jukka S Jurvelin
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland
| | - Simo Saarakkala
- Department of Diagnostic Radiology, Oulu University Hospital, Kajaanintie 50, Oulu FI-90220, Finland; Research Unit of Medical Imaging, Physics and Technology, University of Oulu, P.O. Box 5000, Oulu FI-90014, Finland
| | - Mika E Mononen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, Kuopio FI-70211, Finland
| |
Collapse
|
5
|
A biphasic visco-hyperelastic damage model for articular cartilage: application to micromechanical modelling of the osteoarthritis-induced degradation behaviour. Biomech Model Mechanobiol 2019; 19:1055-1077. [DOI: 10.1007/s10237-019-01270-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 11/22/2019] [Indexed: 01/10/2023]
|
6
|
Maximum shear strain-based algorithm can predict proteoglycan loss in damaged articular cartilage. Biomech Model Mechanobiol 2019; 18:753-778. [PMID: 30631999 DOI: 10.1007/s10237-018-01113-1] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Accepted: 12/24/2018] [Indexed: 01/25/2023]
Abstract
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.
Collapse
|
7
|
Orozco GA, Tanska P, Florea C, Grodzinsky AJ, Korhonen RK. A novel mechanobiological model can predict how physiologically relevant dynamic loading causes proteoglycan loss in mechanically injured articular cartilage. Sci Rep 2018; 8:15599. [PMID: 30348953 PMCID: PMC6197240 DOI: 10.1038/s41598-018-33759-3] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Accepted: 10/02/2018] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Gustavo A Orozco
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland.
| | - Petri Tanska
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
| | - Cristina Florea
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
- Departments of Biological Engineering, Electrical Engineering and Computer Science and Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Alan J Grodzinsky
- Departments of Biological Engineering, Electrical Engineering and Computer Science and Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Rami K Korhonen
- Department of Applied Physics, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
8
|
Kaleem B, Maier F, Drissi H, Pierce DM. Low-energy impact of human cartilage: predictors for microcracking the network of collagen. Osteoarthritis Cartilage 2017; 25:544-553. [PMID: 27903450 DOI: 10.1016/j.joca.2016.11.009] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/31/2016] [Revised: 11/08/2016] [Accepted: 11/21/2016] [Indexed: 02/02/2023]
Abstract
OBJECTIVE We aimed to determine the minimum mechanical impact to cause microstructural damage in the network of collagen (microcracking) within human cartilage and hypothesized that energies below 0.1 J or 1 mJ/mm3 would suffice. DESIGN We completed 108 low-energy impact tests (0.05, 0.07, or 0.09 J; 0.75 or 1.0 m/s2) using healthy cartilage specimens from six male donors (30.2 ± 8.8 yrs old). Before and after impact we acquired, imaging the second harmonic generation (SHG), ten images from each specimen (50 μm depth, 5 μm step size), resulting in 2160 images. We quantified both the presence and morphology of microcracks. We then correlated test parameters (predictors) impact energy/energy dissipation density, nominal stress/stress rate, and strain/strain rate to microcracking and tested for significance. Where predictors significantly correlated with microstructural outcomes we fitted binary logistic regression plots with 95% confidence intervals (CIs). RESULTS No specimens presented visible damage following impact. We found that impact energy/energy dissipation density and nominal stress/stress rate were significant (P < 0.05) predictors of microcracking while both strain and strain rate were not. In our test configuration, an impact energy density of 2.93 mJ/mm3, an energy dissipation density of 1.68 mJ/mm3, a nominal stress of 4.18 MPa, and a nominal stress rate of 689 MPa/s all corresponded to a 50% probability of microcracking in the network of collagen. CONCLUSIONS An impact energy density of 1.0 mJ/mm3 corresponded to a ∼20% probability of microcracking. Such changes may initiate a degenerative cascade leading to post-traumatic osteoarthritis.
Collapse
Affiliation(s)
- B Kaleem
- University of Connecticut, Department of Biomedical Engineering, Storrs, CT, USA
| | - F Maier
- Department of Mechanical Engineering, University of Connecticut, Storrs, CT, USA
| | - H Drissi
- University of Connecticut Health Center, Orthopedic Surgery, Farmington, CT, USA
| | - D M Pierce
- University of Connecticut, Department of Biomedical Engineering, Storrs, CT, USA; Department of Mechanical Engineering, University of Connecticut, Storrs, CT, USA.
| |
Collapse
|
9
|
Klennert BJ, Ellis BJ, Maak TG, Kapron AL, Weiss JA. The mechanics of focal chondral defects in the hip. J Biomech 2016; 52:31-37. [PMID: 28041611 DOI: 10.1016/j.jbiomech.2016.11.056] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Revised: 11/20/2016] [Accepted: 11/22/2016] [Indexed: 11/24/2022]
Abstract
There is a mean incidence of osteoarthritis (OA) of the hip in 8% of the overall population. In the presence of focal chondral defects, defined as localized damage to the articular cartilage, there is an increased risk of symptomatic progression toward OA. This relationship between chondral defects and subsequent development of OA has led to substantial efforts to develop effective procedures for surgical cartilage repair. This study examined the effects of chondral defects and labral delamination on cartilage mechanics in the dysplastic hip during the gait cycle using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with simulated chondral defects at the chondrolabral junction on the posterior acetabulum during five distinct points in the gait cycle. Focal chondral defects increased maximum shear stress on the osteochondral surface of the acetabular cartilage, when compared to the intact case. This effect was amplified with labral delamination. Additionally, chondral defects increased the first principal Lagrange strain on the articular surface of the acetabular cartilage and labrum. Labral delamination relieved some of this tensile strain. As defect size was increased, contact stress increased in the medial zone of the acetabulum, while it decreased anteriorly. The results suggest that in the presence of chondral defects and labral delamination the cartilage experiences elevated tensile strains and shear and contact stress, which could lead to further damage of the cartilage, and subsequent arthritic progression. The framework presented here will serve as the procedure for future finite element studies on cartilage mechanics in hips with varying disease states with simulated chondral defects and labral tears.
Collapse
Affiliation(s)
- Brenden J Klennert
- Department of Bioengineering, University of Utah, 36 S. Wasatch Drive, Rm. 3100, Salt Lake City, UT 84112, United States; Scientific Computing and Imaging Institute, University of Utah, 72 S. Central Campus Drive, Salt Lake City, UT 84112, United States
| | - Benjamin J Ellis
- Department of Bioengineering, University of Utah, 36 S. Wasatch Drive, Rm. 3100, Salt Lake City, UT 84112, United States; Scientific Computing and Imaging Institute, University of Utah, 72 S. Central Campus Drive, Salt Lake City, UT 84112, United States
| | - Travis G Maak
- Department of Orthopaedics, University of Utah, 590 Wakara Way, Salt Lake City, UT 84108, United States
| | - Ashley L Kapron
- Department of Orthopaedics, University of Utah, 590 Wakara Way, Salt Lake City, UT 84108, United States
| | - Jeffrey A Weiss
- Department of Bioengineering, University of Utah, 36 S. Wasatch Drive, Rm. 3100, Salt Lake City, UT 84112, United States; Scientific Computing and Imaging Institute, University of Utah, 72 S. Central Campus Drive, Salt Lake City, UT 84112, United States; Department of Orthopaedics, University of Utah, 590 Wakara Way, Salt Lake City, UT 84108, United States.
| |
Collapse
|
10
|
Kaplan JT, Neu CP, Drissi H, Emery NC, Pierce DM. Cyclic loading of human articular cartilage: The transition from compaction to fatigue. J Mech Behav Biomed Mater 2016; 65:734-742. [PMID: 27756049 DOI: 10.1016/j.jmbbm.2016.09.040] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2016] [Revised: 09/26/2016] [Accepted: 09/30/2016] [Indexed: 11/19/2022]
Abstract
Osteoarthritis and articular cartilage injuries are common conditions in human joints and a frequent cause of pain and disability. Unfortunately, cartilage is avascular and has limited capabilities for self-repair. Despite the societal impact, there is little information on the dynamic process of cartilage degeneration. We performed a series of cyclic unconfined compression tests motivated by in vivo loading conditions and designed to generate mechanical fatigue. We examined the functional (both stress-stretch and creep) responses of the tissue after recovery from a specified number of loading cycles, as well as histology and second harmonic generation microscopy images. The effect of compaction was complimented by the effect of fatigue in our unconfined compression tests. A three-way, repeated-measures mixed model ANOVA showed significant differences between loading with a physiologically relevant low magnitude, and two more severe loading magnitudes, in terms of the resulting specimen stiffness, time to equilibrium and thickness. There was a statistically significant effect of loading frequency on a specimen's time to equilibrium and significant interaction of force and frequency on specimen thickness and time to equilibrium. Increasing the number of loading cycles significantly impacted a specimen's effective stiffness and resulting thickness. We attribute permanent loss of mechanical function under cyclic loading to rearrangement and disruption of the collagen network and resulting proteoglycan (PG) aggregation, as seen in histological and second harmonic generation images, as a result of induced mechanical fatigue.
Collapse
Affiliation(s)
- Jonathan T Kaplan
- Departments of Biomedical and Mechanical Engineering, University of Connecticut, Storrs, CT, United States; U.S. Army Natick Soldier Research, Development and Engineering Center, Natick, MA, United States
| | - Corey P Neu
- Department of Mechanical Engineering, University of Colorado, Boulder, CO, United States
| | - Hicham Drissi
- Department of Orthopaedic Surgery, University of Connecticut Health Center, Farmington, CT, United States
| | - Nancy C Emery
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| | - David M Pierce
- Departments of Biomedical and Mechanical Engineering, University of Connecticut, Storrs, CT, United States.
| |
Collapse
|
11
|
Stender ME, Carpenter RD, Regueiro RA, Ferguson VL. An evolutionary model of osteoarthritis including articular cartilage damage, and bone remodeling in a computational study. J Biomech 2016; 49:3502-3508. [PMID: 27717548 DOI: 10.1016/j.jbiomech.2016.09.024] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2015] [Revised: 09/09/2016] [Accepted: 09/16/2016] [Indexed: 11/26/2022]
Abstract
With osteoarthritis, a complex set of progressive chemical, biological, and mechanical changes occur in both cartilage and bone. The aim of this study is to develop a high-fidelity computational model of the complete bone-cartilage unit to study the evolution of osterarthritis-induced articular cartilage (AC) damage and remodeling of subchondral cortical bone (SCB) and subchondral trabecular bone (STB). A finite element model of spherical indentation was developed with a depth-dependent anisotropic model of degenerating articular cartilage, a calcified cartilage (CC) zone, and SCB and STB remodeling regions. Calcified tissue (CC, SCB, and STB) and AC material regions were integrated to form an evolutionary bone-cartilage unit model. Results indicate that with indentation loading, articular cartilage damage occurs at the articular surface. Furthermore, bone remodeling was predicted to occur with a net stiffening of the subchondral bone plate. Changes in indentation force were minimal (<2%) between initial and final peak indentation loading. However, additional degradation and wear of AC and/or alterations in loading may have more pronounced effects on the mechanical response of the bone-cartilage unit. Bone remodeling and articular cartilage damage predictions are consistent with experimental observations that cartilage damage begins at the articular surface and subchondral bone experiences a thickening (i.e., stiffening) response with osteoarthritis. Our results provide insight into the early-term initiation behavior of osteoarthritis; the potential consequences of evolutions in AC, SCB, and STB with disease progression; and may guide future experimental and computational studies to elucidate mechanisms of osteoarthritis progression.
Collapse
Affiliation(s)
- Michael E Stender
- Department of Mechanical Engineering, University of Colorado, Boulder, CO 80309, United States
| | - R Dana Carpenter
- Department of Mechanical Engineering, University of Colorado, Boulder, CO 80309, United States
| | - Richard A Regueiro
- Department of Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, CO 80309, United States
| | - Virginia L Ferguson
- Department of Mechanical Engineering, University of Colorado, Boulder, CO 80309, United States.
| |
Collapse
|