# Identification of potential biomarkers related to the Bitong mixture in osteoarthritis based on bioinformatics and network pharmacology, and exploration of the mechanism involved

**Authors:** Ziwei Shen, Suming Li, Teng Li, Lining Wang

PMC · DOI: 10.3389/fimmu.2026.1739355 · Frontiers in Immunology · 2026-02-13

## TL;DR

This study identifies MMP9, MMP2, and SPP1 as potential biomarkers for Bitong Mixture in treating osteoarthritis and explores their role in immune and matrix regulation.

## Contribution

The study introduces a novel integrative approach combining bioinformatics and experimental validation to identify OA biomarkers for Bitong Mixture.

## Key findings

- MMP9, MMP2, and SPP1 were identified as potential OA biomarkers related to Bitong Mixture.
- BM treatment reduced OA symptoms and expression of these biomarkers in rats.
- Luteolin from BM showed strong binding affinity to MMP9, suggesting a possible mechanism.

## Abstract

Bitong Mixture (BM) has shown efficacy in alleviating pain in knee osteoarthritis (OA) in clinical practice; however, the molecular mechanisms underlying its therapeutic effects remain to be fully elucidated. This study aimed to identified BM-related OA biomarkers and explore their functional implications.

An integrative strategy combining bioinformatics prediction and experimental validation was used. Biomarkers were screened from public OA transcriptomic data using differential expression analysis, network pharmacology, and machine learning. Their functions were explored via enrichment and immune infiltration analyses. Molecular docking predicted interactions between herbal compounds and targets. Single-cell analysis characterized biomarker expression in chondrocyte subsets. A rat OA model and reverse transcription quantitative polymerase chain reaction (RT-qPCR) were employed for in vivo validation.

Bioinformatic prediction identified three potential biomarkers: MMP9, MMP2, and SPP1. They demonstrated certain diagnostic performance for OA and were implicated in pathways related to extracellular matrix organization and immune regulation. Immune analysis revealed significant correlations, notably between MMP2 and activated dendritic cells (cor = 0.66) and between SPP1 and CD4+ central memory T cells (cor = -0.75). Molecular docking suggested strong binding affinity between luteolin (a BM component) and MMP9. Single-cell analysis indicated high expression of these potential biomarkers in hypertrophic chondrocytes, inflammatory chondrocytes, and fibrochondrocytes. In vivo validation confirmed that BM alleviated OA symptoms and histopathological damage in rats. RT-qPCR results showed that BM treatment alleviated the OA-induced upregulation of MMP9, MMP2, and SPP1 expression.

MMP9, MMP2, and SPP1 are potential therapeutic biomarkers for BM in OA. The efficacy of BM may be attributed to its regulation of extracellular matrix remodeling and immune responses, which provides a possible mechanistic explanation for its clinical use.

## Linked entities

- **Genes:** MMP9 (matrix metallopeptidase 9) [NCBI Gene 4318], MMP2 (matrix metallopeptidase 2) [NCBI Gene 4313], SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696]
- **Chemicals:** luteolin (PubChem CID 5280445)
- **Diseases:** osteoarthritis (MONDO:0005178)
- **Species:** Rattus norvegicus (taxon 10116)

## Full-text entities

- **Genes:** SERPINE1 (serpin family E member 1) [NCBI Gene 5054] {aka PAI, PAI-1, PAI1, PLANH1}, LEP (leptin) [NCBI Gene 3952] {aka LEPD, OB, OBS}, Cd4 (Cd4 molecule) [NCBI Gene 24932] {aka W3/25, p55}, TAC1 (tachykinin precursor 1) [NCBI Gene 6863] {aka Hs.2563, NK2, NKNA, NPK, TAC2}, DOCK2 (dedicator of cytokinesis 2) [NCBI Gene 1794] {aka IMD40}, Mmp2 (matrix metallopeptidase 2) [NCBI Gene 81686], RENBP (renin binding protein) [NCBI Gene 5973] {aka RBP, RNBP}, RASSF1 (Ras association domain family member 1) [NCBI Gene 11186] {aka 123F2, NORE2A, RASSF1A, RDA32, REH3P21}, Mmp9 (matrix metallopeptidase 9) [NCBI Gene 81687], Spp1 (secreted phosphoprotein 1) [NCBI Gene 25353] {aka OSP}, NUF2 (NUF2 component of NDC80 kinetochore complex) [NCBI Gene 83540] {aka CDCA1, CT106, NUF2R}, IL6 (interleukin 6) [NCBI Gene 3569] {aka BSF-2, BSF2, CDF, HGF, HSF, IFN-beta-2}, MMP9 (matrix metallopeptidase 9) [NCBI Gene 4318] {aka CLG4B, GELB, MANDP2, MMP-9}, FN1 (fibronectin 1) [NCBI Gene 2335] {aka CIG, ED-B, FINC, FN, FNZ, GFND}, S100A12 (S100 calcium binding protein A12) [NCBI Gene 6283] {aka CAAF1, CAGC, CGRP, ENRAGE, MRP-6, MRP6}, TIMP1 (TIMP metallopeptidase inhibitor 1) [NCBI Gene 7076] {aka CLGI, EPA, EPO, HCI, TIMP, TIMP-1}, ADIPOQ (adiponectin, C1Q and collagen domain containing) [NCBI Gene 9370] {aka ACDC, ACRP30, ADIPQTL1, ADPN, APM-1, APM1}, MMP2 (matrix metallopeptidase 2) [NCBI Gene 4313] {aka CLG4, CLG4A, MMP-2, MMP-II, MONA, TBE-1}, NAIPP2 (NAIP pseudogene 2) [NCBI Gene 728519] {aka BIRC, NAIP1B}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696] {aka BNSP, BSPI, ETA-1, OPN}, IL1B (interleukin 1 beta) [NCBI Gene 3553] {aka IL-1, IL1-BETA, IL1F2, IL1beta}, Gapdh (glyceraldehyde-3-phosphate dehydrogenase) [NCBI Gene 24383] {aka BARS-38, Gapd}, AGER (advanced glycosylation end-product specific receptor) [NCBI Gene 177] {aka RAGE, SCARJ1, sRAGE}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, ITGA5 (integrin subunit alpha 5) [NCBI Gene 3678] {aka CD49e, FNRA, VLA-5, VLA5A}, PTGS2 (prostaglandin-endoperoxide synthase 2) [NCBI Gene 5743] {aka COX-2, COX2, GRIPGHS, PGG/HS, PGHS-2, PHS-2}
- **Diseases:** hyperplasia (MESH:D006965), MIA (MESH:C562377), renal impairment (MESH:D007674), musculoskeletal conditions (MESH:D009140), inflammatory cytokines (MESH:D000080424), toxicity (MESH:D064420), joint damage (MESH:D007592), Arthralgia (MESH:D018771), osteolysis (MESH:D010014), knee osteoarthritis (MESH:D020370), peptic ulcers (MESH:D010437), synovitis (MESH:D013585), immune dysfunction (MESH:D007154), cardiovascular complications (MESH:D002318), subchondral bone sclerosis (MESH:D001845), osteophyte formation (MESH:D054850), joint deformity (MESH:D016916), diabetic complications (MESH:D048909), cartilage (MESH:D002357), SLE (MESH:D008180), metabolic dysfunction (MESH:D008659), OA (MESH:D010003), arthritic (MESH:D015535), obesity (MESH:D009765), tumor (MESH:D009369), calcification (MESH:D002114), pain (MESH:D010146), inflammation (MESH:D007249), degenerative joint disease (MESH:D019636)
- **Chemicals:** alcohol (MESH:D000438), luteolin (MESH:D047311), eosin (MESH:D004801), paraformaldehyde (MESH:C003043), leukotrienes (MESH:D015289), CB (MESH:C063451), hematoxylin (MESH:D006416), H&amp;E (MESH:D006371), Safranin (MESH:C009195), Pericarpium Citri Reticulatae (-), hyaluronic acid (MESH:D006820), insulin (MESH:D007328), ethanol (MESH:D000431), hydrochloric acid (MESH:D006851), TRIzol (MESH:C411644), isoflurane (MESH:D007530), MIA (MESH:D019807), water (MESH:D014867), quercetin (MESH:D011794), EDTA (MESH:D004492), nobiletin (MESH:C008661), xylene (MESH:D014992), oxygen (MESH:D010100), saline (MESH:D012965), histamine (MESH:D006632), paraffin (MESH:D010232)
- **Species:** Curcuma longa (turmeric, species) [taxon 136217], Homo sapiens (human, species) [taxon 9606], Rattus norvegicus (brown rat, species) [taxon 10116], Paeonia lactiflora (Chinese peony, species) [taxon 35924], Aconitum kusnezoffii (species) [taxon 239685], Eucommia ulmoides (species) [taxon 4392], Aconitum carmichaelii (species) [taxon 85363]
- **Cell lines:** HTC — Homo sapiens (Human), Transformed cell line (CVCL_B0BV)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12946070/full.md

## Figures

11 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12946070/full.md

## References

86 references — full list in the complete paper: https://tomesphere.com/paper/PMC12946070/full.md

---
Source: https://tomesphere.com/paper/PMC12946070