Confounder-Adjusted Transcriptomic Signature Reversal and Multi-Database Mechanistic Triangulation of Phytochemical Perturbagens in Human Osteoarthritic Synovial Fibroblasts: An Integrative Study of GSE29746, LINCS L1000, HERB 2.0, and the Comparative Toxicogenomics Database

Authors

Rizki Ayu

Department of Drug Safety and Effectiveness, CMHC Research Center, Palembang (Indonesia)

Ramakhrisnand Ramakhrisnand

Department of Drug Development, Avant Medical Center, Avarua (Cook Islands)

Aleisha Wulandari

Department of Nutrition, JTC Medical Center, Jakarta (Indonesia)

Sudarto

Department of Botanical Sciences, CMHC Research Center, Palembang (Indonesia)

Tiffany Gabriele

Department of Psychiatry, Gabon State Hospital, Libreville (Gabon)

Brenda Jaleel

Department of Neuroscience, San Fernando General Hospital, San Fernando (Trinidad and Tobago)

Andhika Kurnianta Kusuma

Department of Oral Health and Dentistry, Phlox Institute: Indonesian Medical Research Organization, Palembang (Indonesia)

Tomiola Owkwulu

Department of Obstetrics and Gynecology, Mogadishu State Hospital, Mogadishu (Somalia)

Cindy Susanti

Department of Internal Medicine, Provita Hospital, Jayapura (Indonesia)

Winata Putri

Department of Pediatrics, CMHC Research Center, Palembang (Indonesia)

Article Information

DOI: 10.51584/IJRIAS.2026.11100061

Subject Category: Health Science

Volume/Issue: 11/10 | Page No: 955-980

Publication Timeline

Submitted: 2026-10-18

Accepted: 2026-10-23

Published: 2026-11-05

Abstract

Background: Synovial fibroblast activation drives synovitis, matrix remodeling, and pain in osteoarthritis (OA), yet connectivity screens rarely audit demographic confounding, resampling stability, permutation controls, or curated human directional gene reversal. Objective: To computationally prioritize phytochemical perturbagens reversing the human OA synovial fibroblast transcriptomic signature by integrating GSE29746, LINCS L1000, HERB 2.0, and the Comparative Toxicogenomics Database (CTD). Methods: Passage-4 synovial fibroblast microarrays from GSE29746 (N = 22: 11 OA vs. 11 healthy controls; 19,261 genes) were analyzed via limma with and without age/sex adjustment, audited by jackknife/bootstrap resampling and external synovial validation (GSE55235/GSE55457), and queried in L1000CDS2 against 10,000-permutation controls, curated HERB 2.0/CTD links, and ADME/docking benchmarks. Results: Sex was balanced (6F/5M; p = 1.000), whereas OA donors were older (73.0 ± 7.1 vs. 59.0 ± 19.0 years; Δ = +14.0 years, Welch p = 0.040). Age/sex adjustment attenuated FDR < 0.05 hits by 93.5% (from 46 to 3 downregulated genes: MMP27 [log2FC = −1.733, FDR = 0.033], IFIT3 [−1.238, FDR = 0.033], and MX2 [−2.057, FDR = 0.048]; Spearman ρ = 0.914; jackknife concordance = 98.4%; external p ≤ 0.0041). Among 52 unique LINCS compounds, 5 (9.6%) matched curated HERB 2.0 phytochemicals. Celastrol ranked first (Tier B; score = 0.0350, p = 0.0038; HME1, 10 μM, 3 h) with curated CTD reversal on MDH2 and TNFSF10 (p = 0.0019) and GRP94/IKKβ/MDH2 docking (RMSD ≤ 0.91 Å). Reserpine (0.0435), emetine (0.0435), and homoharringtonine (0.0400) replicated across models (Tier C) without CTD reversal; isoeugenitol (0.0300) ranked Tier D. Reactome rRNA processing was enriched (p = 0.028). Zero candidates reached Tier A. Conclusions: Triangulation prioritizes celastrol for in vitro (IL-1β/TNF-α-stimulated OASF measuring IL-6, MMP-13, and COX-2) and in vivo testing before clinical interpretation.

Keywords

celastrol; Comparative Toxicogenomics Database; connectivity mapping; HERB 2.0; LINCS L1000; osteoarthritis; synovial fibroblasts

Downloads

References

1. Steinmetz JD, Culbreth GT, Haile LM, et al. Global, regional, and national burden of osteoarthritis, 1990–2020 and projections to 2050: a systematic analysis for the Global Burden of Disease Study 2021. Lancet Rheumatol. 2023;5(9):e508–e522. https://doi.org/10.1016/S2665-9913(23)00163-7 [Google Scholar] [Crossref]

2. Knights AJ, Farrell EC, Ellis OM, et al. Synovial fibroblasts assume distinct functional identities and secrete R-spondin 2 in osteoarthritis. Ann Rheum Dis. 2023;82(2):272–282. https://doi.org/10.1136/ard-2022-222773 [Google Scholar] [Crossref]

3. Nanus DE, Badoume A, Wijesinghe SN, et al. Synovial tissue from sites of joint pain in knee osteoarthritis patients exhibits a differential phenotype with distinct fibroblast subsets. EBioMedicine. 2021;72:103618. https://doi.org/10.1016/j.ebiom.2021.103618 [Google Scholar] [Crossref]

4. Philpott HT, Birmingham TB, Blackler G, et al. Association of synovial innate immune exhaustion with worse pain in knee osteoarthritis. Arthritis Rheumatol. 2025;77(6):664–676. https://doi.org/10.1002/art.43089 [Google Scholar] [Crossref]

5. Huang ZY, Luo ZY, Cai YR, et al. Single cell transcriptomics in human osteoarthritis synovium and in silico deconvoluted bulk RNA sequencing. Osteoarthritis Cartilage. 2022;30(3):475–480. https://doi.org/10.1016/j.joca.2021.12.007 [Google Scholar] [Crossref]

6. Liu Z, Sun Y, Pan J, et al. Single-cell profiling uncovers synovial fibroblast subpopulations associated with chondrocyte injury in osteoarthritis. Front Endocrinol (Lausanne). 2024;15:1479909. https://doi.org/10.3389/fendo.2024.1479909 [Google Scholar] [Crossref]

7. Liu B, Xian Y, Chen X, et al. Inflammatory fibroblast-like synoviocyte-derived exosomes aggravate osteoarthritis via enhancing macrophage glycolysis. Adv Sci (Weinh). 2024;11(14):e2307338. https://doi.org/10.1002/advs.202307338 [Google Scholar] [Crossref]

8. Del Rey MJ, Usategui A, Izquierdo E, et al. Transcriptome analysis reveals specific changes in osteoarthritis synovial fibroblasts. Ann Rheum Dis. 2012;71(2):275–280. https://doi.org/10.1136/annrheumdis-2011-200281 [Google Scholar] [Crossref]

9. Yang G, Wang K, Song H, et al. Celastrol ameliorates osteoarthritis via regulating TLR2/NF-κB signaling pathway. Front Pharmacol. 2022;13:963506. https://doi.org/10.3389/fphar.2022.963506 [Google Scholar] [Crossref]

10. Duan Q, Reid SP, Clark NR, et al. L1000CDS2: LINCS L1000 characteristic direction signatures search engine. npj Syst Biol Appl. 2016;2:16015. https://doi.org/10.1038/npjsba.2016.15 [Google Scholar] [Crossref]

11. Subramanian A, Narayan R, Corsello SM, et al. A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell. 2017;171(6):1437–1452.e17. https://doi.org/10.1016/j.cell.2017.10.049 [Google Scholar] [Crossref]

12. Evangelista JE, Clarke DJB, Xie Z, et al. SigCom LINCS: data and metadata search engine for a million gene expression signatures. Nucleic Acids Res. 2022;50(W1):W697–W709. https://doi.org/10.1093/nar/gkac328 [Google Scholar] [Crossref]

13. Marino GB, Evangelista JE, Clarke DJB, et al. L2S2: chemical perturbation and CRISPR KO LINCS L1000 signature search engine. Nucleic Acids Res. 2025;53(W1):W338–W350. https://doi.org/10.1093/nar/gkaf373 [Google Scholar] [Crossref]

14. Gao K, Liu L, Lei S, et al. HERB 2.0: an updated database integrating clinical and experimental evidence for traditional Chinese medicine. Nucleic Acids Res. 2025;53(D1):D1404–D1414. https://doi.org/10.1093/nar/gkae1037 [Google Scholar] [Crossref]

15. Davis AP, Wiegers TC, Sciaky D, et al. Comparative Toxicogenomics Database's 20th anniversary: update 2025. Nucleic Acids Res. 2025;53(D1):D1328–D1334. https://doi.org/10.1093/nar/gkae883 [Google Scholar] [Crossref]

16. Clough E, Barrett T, Wilhite SE, et al. NCBI GEO: archive for gene expression and epigenomics data sets: 23-year update. Nucleic Acids Res. 2024;52(D1):D138–D144. https://doi.org/10.1093/nar/gkad965 [Google Scholar] [Crossref]

17. Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007 [Google Scholar] [Crossref]

18. Chen Y, Chen L, Lun ATL, et al. edgeR v4: powerful differential analysis of sequencing data with expanded functionality and improved support for small counts and larger datasets. Nucleic Acids Res. 2025;53(2):gkaf018. https://doi.org/10.1093/nar/gkaf018 [Google Scholar] [Crossref]

19. Kim S, Chen J, Cheng T, et al. PubChem 2023 update. Nucleic Acids Res. 2023;51(D1):D1373–D1380. https://doi.org/10.1093/nar/gkac956 [Google Scholar] [Crossref]

20. Zdrazil B, Felix E, Hunter F, et al. The ChEMBL Database in 2023: a drug discovery platform spanning multiple bioactivity data types and time periods. Nucleic Acids Res. 2024;52(D1):D1180–D1192. https://doi.org/10.1093/nar/gkad1004 [Google Scholar] [Crossref]

21. Kolberg L, Raudvere U, Kuzmin I, et al. g:Profiler—interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update). Nucleic Acids Res. 2023;51(W1):W207–W212. https://doi.org/10.1093/nar/gkad347 [Google Scholar] [Crossref]

22. Milacic M, Beavers D, Conley P, et al. The Reactome Pathway Knowledgebase 2024. Nucleic Acids Res. 2024;52(D1):D672–D678. https://doi.org/10.1093/nar/gkad1025 [Google Scholar] [Crossref]

23. Aleksander SA, Balhoff J, Carbon S, et al. The Gene Ontology knowledgebase in 2023. Genetics. 2023;224(1):iyad031. https://doi.org/10.1093/genetics/iyad031 [Google Scholar] [Crossref]

24. Kanehisa M, Furumichi M, Sato Y, et al. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–D592. https://doi.org/10.1093/nar/gkac963 [Google Scholar] [Crossref]

25. Szklarczyk D, Kirsch R, Koutrouli M, et al. The STRING database in 2023: protein–protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023;51(D1):D638–D646. https://doi.org/10.1093/nar/gkac1000 [Google Scholar] [Crossref]

26. Huang C, Zeng B, Zhou B, et al. Single-cell transcriptomic analysis of chondrocytes in cartilage and pathogenesis of osteoarthritis. Genes Dis. 2024;12(2):101241. https://doi.org/10.1016/j.gendis.2024.101241 [Google Scholar] [Crossref]

27. Atasoy-Zeybek A, Hawse GP, Nagelli CV, et al. Transcriptomic changes during the replicative senescence of human articular chondrocytes. Int J Mol Sci. 2024;25(22):12130. https://doi.org/10.3390/ijms252212130 [Google Scholar] [Crossref]

28. Chabronova A, van den Akker GGH, Housmans BAC, et al. Ribosomal RNA-based epitranscriptomic regulation of chondrocyte translation and proteome in osteoarthritis. Osteoarthritis Cartilage. 2023;31(3):374–385. https://doi.org/10.1016/j.joca.2022.12.010 [Google Scholar] [Crossref]

29. Chen HY, Swaroop M, Papal S, et al. Reserpine maintains photoreceptor survival in retinal ciliopathy by resolving proteostasis imbalance and ciliogenesis defects. Elife. 2023;12:e83205. https://doi.org/10.7554/eLife.83205 [Google Scholar] [Crossref]

30. Sun Y, Dai J, Jiao R, et al. Homoharringtonine inhibits fibroblasts proliferation, extracellular matrix production and reduces surgery-induced knee arthrofibrosis via PI3K/AKT/mTOR pathway-mediated apoptosis. J Orthop Surg Res. 2021;16(1):9. https://doi.org/10.1186/s13018-020-02150-2 [Google Scholar] [Crossref]

31. Kumar R, Afsar M, Khandelwal N, et al. Emetine suppresses SARS-CoV-2 replication by inhibiting interaction of viral mRNA with eIF4E. Antiviral Res. 2021;189:105056. https://doi.org/10.1016/j.antiviral.2021.105056 [Google Scholar] [Crossref]

32. Valipour M, Irannejad H, Emami S. Application of emetine in SARS-CoV-2 treatment: regulation of p38 MAPK signaling pathway for preventing emetine-induced cardiac complications. Cell Cycle. 2022;21(22):2379–2386. https://doi.org/10.1080/15384101.2022.2100575 [Google Scholar] [Crossref]

33. Han J, Zhan LN, Huang Y, et al. Moderate mechanical stress suppresses chondrocyte ferroptosis in osteoarthritis by regulating NF-κB p65/GPX4 signaling pathway. Sci Rep. 2024;14(1):5078. https://doi.org/10.1038/s41598-024-55629-x [Google Scholar] [Crossref]

34. Zou Z, Pan S, Sun C, et al. AM1241 inhibits chondrocyte inflammation and ECM degradation through the Nrf2/HO-1 and NF-κB pathways and alleviates osteoarthritis in mice. Mol Med. 2025;31(1):9. https://doi.org/10.1186/s10020-024-01012-5 [Google Scholar] [Crossref]

35. Chang YT, Huang KC, Pranata R, et al. Evaluation of the protective effects of chondroitin sulfate oligosaccharide against osteoarthritis via inactivation of NLRP3 inflammasome by in vivo and in vitro studies. Int Immunopharmacol. 2024;142(Pt A):113148. https://doi.org/10.1016/j.intimp.2024.113148 [Google Scholar] [Crossref]

36. Wang X, Yu H, Zhang Y, et al. Curcumin alleviates osteoarthritis through the p38MAPK pathway: network pharmacological prediction and experimental confirmation. J Inflamm Res. 2024;17:5039–5056. https://doi.org/10.2147/JIR.S459867 [Google Scholar] [Crossref]

37. Guo P, Li H, Wang X, et al. PG545 prevents osteoarthritis development by regulating PI3K/AKT/mTOR signaling and activating chondrocyte autophagy. Pharmacology. 2023;108(6):576–588. https://doi.org/10.1159/000532078 [Google Scholar] [Crossref]

38. Chen W, Xiao J, Zhou Y, et al. Curcumenol regulates Histone H3K27me3 demethylases KDM6B affecting Succinic acid metabolism to alleviate cartilage degeneration in knee osteoarthritis. Phytomedicine. 2024;133:155922. https://doi.org/10.1016/j.phymed.2024.155922 [Google Scholar] [Crossref]

39. Gao X, Liu H, Guo J, et al. OGT regulates histone demethylation and alleviates osteoarthritis progression by O-GlcNAcylation of KDM6B. J Adv Res. 2026;87:455–471. https://doi.org/10.1016/j.jare.2026.01.003 [Google Scholar] [Crossref]

40. Yang C, Dong W, Wang Y, et al. DDIT3 aggravates TMJOA cartilage degradation via Nrf2/HO-1/NLRP3-mediated autophagy. Osteoarthritis Cartilage. 2024;32(8):921–937. https://doi.org/10.1016/j.joca.2024.04.017 [Google Scholar] [Crossref]

41. Zhao Q, Xu T, Du Z, et al. Cartilage targeting hydrogel nanoplatform degrades BRD4 to alleviate osteoarthritis via Nav1.7 axis. Nat Commun. 2026;17(1):4573. https://doi.org/10.1038/s41467-026-71246-w [Google Scholar] [Crossref]

42. Sun J, Wang XH, Song FH, et al. Inhibition of Brd4 alleviates osteoarthritis pain via suppression of neuroinflammation and activation of Nrf2-mediated antioxidant signalling. Br J Pharmacol. 2023;180(24):3194–3214. https://doi.org/10.1111/bph.16195 [Google Scholar] [Crossref]

43. Alenazi F, Moursi S, Mahmoud MR, et al. Withaferin A alleviates inflammation in animal models of arthritis by inhibiting the NF-κB pathway and cytokine release. Chem Biol Interact. 2024;398:111114. https://doi.org/10.1016/j.cbi.2024.111114 [Google Scholar] [Crossref]

44. Liu T, Wang J, Tong Y, et al. Integrating network pharmacology and animal experimental validation to investigate the action mechanism of oleanolic acid in obesity. J Transl Med. 2024;22(1):86. https://doi.org/10.1186/s12967-023-04840-x [Google Scholar] [Crossref]

45. Peng W, Yang Y, Lu H, et al. Network pharmacology combines machine learning, molecular simulation dynamics and experimental validation to explore the mechanism of acetylbinankadsurin A in the treatment of liver fibrosis. J Ethnopharmacol. 2024;323:117682. https://doi.org/10.1016/j.jep.2023.117682 [Google Scholar] [Crossref]

Metrics

Views & Downloads

Similar Articles