Injury-related cell death and proteoglycan loss in articular cartilage: Numerical model combining necrosis, reactive oxygen species, and inflammatory cytokines

PLoS Comput Biol. 2023 Jan 26;19(1):e1010337. doi: 10.1371/journal.pcbi.1010337. eCollection 2023 Jan.

Abstract

Osteoarthritis (OA) is a common musculoskeletal disease that leads to deterioration of articular cartilage, joint pain, and decreased quality of life. When OA develops after a joint injury, it is designated as post-traumatic OA (PTOA). The etiology of PTOA remains poorly understood, but it is known that proteoglycan (PG) loss, cell dysfunction, and cell death in cartilage are among the first signs of the disease. These processes, influenced by biomechanical and inflammatory stimuli, disturb the normal cell-regulated balance between tissue synthesis and degeneration. Previous computational mechanobiological models have not explicitly incorporated the cell-mediated degradation mechanisms triggered by an injury that eventually can lead to tissue-level compositional changes. Here, we developed a 2-D mechanobiological finite element model to predict necrosis, apoptosis following excessive production of reactive oxygen species (ROS), and inflammatory cytokine (interleukin-1)-driven apoptosis in cartilage explant. The resulting PG loss over 30 days was simulated. Biomechanically triggered PG degeneration, associated with cell necrosis, excessive ROS production, and cell apoptosis, was predicted to be localized near a lesion, while interleukin-1 diffusion-driven PG degeneration was manifested more globally. Interestingly, the model also showed proteolytic activity and PG biosynthesis closer to the levels of healthy tissue when pro-inflammatory cytokines were rapidly inhibited or cleared from the culture medium, leading to partial recovery of PG content. The numerical predictions of cell death and PG loss were supported by previous experimental findings. Furthermore, the simulated ROS and inflammation mechanisms had longer-lasting effects (over 3 days) on the PG content than localized necrosis. The mechanobiological model presented here may serve as a numerical tool for assessing early cartilage degeneration mechanisms and the efficacy of interventions to mitigate PTOA progression.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Apoptosis
  • Cartilage, Articular* / metabolism
  • Cartilage, Articular* / pathology
  • Cytokines / metabolism
  • Humans
  • Interleukin-1 / metabolism
  • Interleukin-1 / pharmacology
  • Necrosis / metabolism
  • Necrosis / pathology
  • Osteoarthritis* / metabolism
  • Proteoglycans
  • Quality of Life
  • Reactive Oxygen Species / metabolism

Substances

  • Proteoglycans
  • Cytokines
  • Reactive Oxygen Species
  • Interleukin-1

Grants and funding

We acknowledge funding support from: The Doctoral Programme in Science, Technology and Computing (SCITECO) (JPK), Strategic funding of the University of Eastern Finland (PT), Academy of Finland (grant nos. 334773 – under the frame of ERA PerMed, 324529) (RKK), Novo Nordisk Foundation (grant no. NNF21OC0065373, the Center for Mathematical Modeling of Knee Osteoarthritis (MathKOA)) (RKK), Alfred Kordelin Foundation (grant no 190317) (PT), Maire Lisko Foundation (PT), Sigrid Juselius Foundation (RKK, DDA), Saastamoinen Foundation (PT), Instrumentarium Science Foundation (ASAE), the Swedish Research Council (2019-00953—under the frame of ERA PerMed). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.