# Mechanisms of cellular senescence combined with molecular docking strategies: A biomarker study of potential therapeutic targets for allergic rhinitis

**Authors:** Qingyong Chen, Han Zhang, Yue Wu, Jixin Liu, Dongqing Wang, Dezhong Sun, Zhipeng Chen, Liqiang Lin, Huaiqing Lv, Qiang Shao

PMC · DOI: 10.1371/journal.pone.0338309 · PLOS One · 2026-01-09

## TL;DR

This study identifies key genes linked to cellular senescence in allergic rhinitis, offering potential biomarkers for new treatments.

## Contribution

The novel integration of bioinformatics and molecular docking reveals new therapeutic targets for allergic rhinitis.

## Key findings

- 20 differential senescence genes were identified by intersecting AR-related genes with cellular senescence genes.
- 19 hub genes, including CCL2 and STAT1, were found to be central in AR-related senescence networks.
- Molecular docking validated drug interactions with hub genes like STAT1 and TLR2.

## Abstract

Bioinformatics and molecular docking methods were used to screen potential biomarkers of cellular senescence in allergic rhinitis (Allergic rhinitis AR), which provided a theoretical basis for revealing the mechanism of AR and exploring new therapeutic approaches.

Four AR-related gene chips (GSE19187, GSE43523, GSE44037, and GSE51392) were downloaded from the gene expression database (GEO) for data pooling. Screening differential genes (DEGs) were taken to intersect with cellular senescence-related genes (SRGs) to obtain differential senescence genes (DESRGs). The differential senescence genes were subjected to Gene Ontology Database (GO) functional analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and GSEA enrichment analysis. Protein-protein interaction (PPI) networks were constructed through the STRING database, MCODE plugin weights were analyzed to identify important gene cluster modules, and Hub genes were screened using the CytoHubba plugin topological network algorithm. Hub gene protein interactions network (GeneMANIA) was constructed by the GeneMANIA database. Predict Hub gene construct mRNA-miNA-lncRNA interactions by miRanda, miRDB, miRWalk, TargetScan, and spongeScan databases; construct Hub gene transcription factor regulatory networks by TRRUST database; analyze Hub gene-drug interactions by DGIdb database and select commonly used drugs in the clinic for molecular docking validation.

A total of 264 differential genes were screened in the training set with corrected P.adj < 0.05 and |log2FC| ≥ 1.2 as the filtering condition, and a total of 866 cellular senescence genes, and 20 differential senescence genes (DESRGs) were obtained by taking the intersection of the two. A total of 19 Hub genes were obtained after PPI analysis, which were CCL2, STAT1, TLR2, IGFBP3, TLR3, KLF4, IL1RN, IRF1, SERPINB2, DPP4, MME, NQO1, SAMHD1, XAF1, PHGDH, EIF4EBP1, CTH, HSPA2, AHR The gene-protein interaction network identified 19 Hub genes associated with 21 functional proteins. 5 of the Hub gene loci were associated with 29 miRNAs and 53 lncRNAs. The transcription factor regulatory network obtained 15 transcription factors capable of regulating Hub genes. The analysis of drug–gene interactions identified 489 drugs that target hub genes. For example, in the case of budesonide, the interacting genes STAT1, TLR2, TLR3, and AHR were selected for molecular docking. Similarly, for mometasone, the interacting genes TLR2 and CTH were chosen for molecular docking.

Mining AR-related Hub senescence genes by bioinformatics analysis, constructing PPI network, ceRNA network, transcription factor regulatory network, gene-drug interaction network and molecular docking validation, we screened 19 CCL2, STAT1, TLR2, IGFBP3, TLR3, KLF4, IL1RN, IRF1, SERPINB2, DPP4, MME, NQO1, SAMHD1, XAF1, PHGDH, EIF4EBP1, CTH, HSPA2, and AHR are expected to be Hub genes for potential diagnostic and therapeutic biomarkers, which will provide targets and new insights for further in-depth explorations of AR cellular senescence-related mechanisms of action and therapy.

## Linked entities

- **Genes:** CCL2 (C-C motif chemokine ligand 2) [NCBI Gene 6347], STAT1 (signal transducer and activator of transcription 1) [NCBI Gene 6772], TLR2 (toll like receptor 2) [NCBI Gene 7097], IGFBP3 (insulin like growth factor binding protein 3) [NCBI Gene 3486], TLR3 (toll like receptor 3) [NCBI Gene 7098], KLF4 (KLF transcription factor 4) [NCBI Gene 9314], IL1RN (interleukin 1 receptor antagonist) [NCBI Gene 3557], IRF1 (interferon regulatory factor 1) [NCBI Gene 3659], SERPINB2 (serpin family B member 2) [NCBI Gene 5055], DPP4 (dipeptidyl peptidase 4) [NCBI Gene 1803], MME (membrane metalloendopeptidase) [NCBI Gene 4311], NQO1 (NAD(P)H quinone dehydrogenase 1) [NCBI Gene 1728], SAMHD1 (SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 1) [NCBI Gene 25939], XAF1 (XIAP associated factor 1) [NCBI Gene 54739], PHGDH (phosphoglycerate dehydrogenase) [NCBI Gene 26227], EIF4EBP1 (eukaryotic translation initiation factor 4E binding protein 1) [NCBI Gene 1978], CTH (cystathionine gamma-lyase) [NCBI Gene 1491], HSPA2 (heat shock protein family A (Hsp70) member 2) [NCBI Gene 3306], AHR (aryl hydrocarbon receptor) [NCBI Gene 196]
- **Diseases:** allergic rhinitis (MONDO:0011786)

## Full-text entities

- **Genes:** ELAVL2 (ELAV like RNA binding protein 2) [NCBI Gene 1993] {aka HEL-N1, HELN1, HUB}, IRF1 (interferon regulatory factor 1) [NCBI Gene 3659] {aka IMD117, IRF-1, MAR}, TLR3 (toll like receptor 3) [NCBI Gene 7098] {aka CD283, IIAE2, IMD83}, SAMHD1 (SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 1) [NCBI Gene 25939] {aka CHBL2, DCIP, HDDC1, MOP-5, SBBI88, hSAMHD1}, AHR (aryl hydrocarbon receptor) [NCBI Gene 196] {aka FVH3, RP85, bHLHe76}, MME (membrane metalloendopeptidase) [NCBI Gene 4311] {aka CALLA, CD10, CMT2T, NEP, SCA43, SFE}, DPP4 (dipeptidyl peptidase 4) [NCBI Gene 1803] {aka ADABP, ADCP2, CD26, DPPIV, TP103}, IL1RN (interleukin 1 receptor antagonist) [NCBI Gene 3557] {aka CRMO2, DIRA, ICIL-1RA, IL-1RN, IL-1ra, IL-1ra3}, SERPINB2 (serpin family B member 2) [NCBI Gene 5055] {aka HsT1201, PAI, PAI-2, PAI2, PLANH2}, IGFBP3 (insulin like growth factor binding protein 3) [NCBI Gene 3486] {aka BP-53, IBP-3, IBP3, IGFBP-3}, VSIG2 (V-set and immunoglobulin domain containing 2) [NCBI Gene 23584] {aka 2210413P10Rik, CTH, CTXL}, EIF4EBP1 (eukaryotic translation initiation factor 4E binding protein 1) [NCBI Gene 1978] {aka 4E-BP1, 4EBP1, BP-1, PHAS-I}, CCL2 (C-C motif chemokine ligand 2) [NCBI Gene 6347] {aka GDCF-2, HC11, HSMCR30, MCAF, MCP-1, MCP1}, PHGDH (phosphoglycerate dehydrogenase) [NCBI Gene 26227] {aka 3-PGDH, 3PGDH, HEL-S-113, NLS, NLS1, PDG}, NQO1 (NAD(P)H quinone dehydrogenase 1) [NCBI Gene 1728] {aka DHQU, DIA4, DTD, NMOR1, NMORI, QR1}, KLF4 (KLF transcription factor 4) [NCBI Gene 9314] {aka EZF, GKLF}, XAF1 (XIAP associated factor 1) [NCBI Gene 54739] {aka BIRC4BP, HSXIAPAF1, XIAPAF1}, STAT1 (signal transducer and activator of transcription 1) [NCBI Gene 6772] {aka CANDF7, IMD31A, IMD31B, IMD31C, ISGF-3, STAT91}, TLR2 (toll like receptor 2) [NCBI Gene 7097] {aka CD282, TIL4}, HSPA2 (heat shock protein family A (Hsp70) member 2) [NCBI Gene 3306] {aka HSP70-2, HSP70-3}
- **Diseases:** Allergic rhinitis AR (MESH:D065631), AR (MESH:D013734)
- **Chemicals:** mometasone (MESH:D000068656), budesonide (MESH:D019819)

## Full text

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

## Figures

13 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12788676/full.md

## References

44 references — full list in the complete paper: https://tomesphere.com/paper/PMC12788676/full.md

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