Abstract
Objective: To explore the effective targets of Celecoxib in the treatment of heterotopic ossification using network pharmacology methods. Methods: Potential molecules related to heterotopic ossification were obtained by retrieving the GEO and CTD databases and intersecting them. Potential binding targets of Celecoxib were acquired from the STITCH database. A protein-protein interaction network was constructed between potential binding targets of Celecoxib and potential related molecules of heterotopic ossification using the STRING database. Molecules in the protein-protein interaction network were further analyzed using GO and KEGG enrichment analysis in R software, followed by enrichment analysis of active molecules in the Celecoxib-heterotopic ossification target dataset. Hub genes were selected based on the "degree" value and enrichment within the protein-protein interaction network. The binding affinity of hub genes to Celecoxib was observed using molecular docking techniques. Finally, in vitro experiments were conducted to validate the effectiveness of hub genes and explore their regulatory role in the progression of heterotopic ossification. Additionally, the therapeutic effect of Celecoxib, which modulates the expression of the hub genes, was investigated in the treatment of heterotopic ossification. Results: 568 potential molecules related to heterotopic ossification and 76 potential binding targets of Celecoxib were identified. After intersection, 13 potential functional molecules in Celecoxib's treatment of heterotopic ossification were obtained. KEGG analysis suggested pathways such as Rheumatoid arthritis, NF-kappa B signaling pathway, Pathways in cancer, Antifolate resistance, MicroRNAs in cancer play a role in the treatment of heterotopic ossification by Celecoxib. Further enrichment analysis of the 13 potential functional molecules identified 5 hub genes: IL6, CCND1, PTGS2, IGFBP3, CDH1. Molecular docking results indicated that Celecoxib displayed excellent binding affinity with CCND1 among the 5 hub genes. Experimental validation found that CCND1 is highly expressed in the progression of heterotopic ossification, promoting heterotopic ossification in the early stages and inhibiting it in the later stages, with Celecoxib's treatment of heterotopic ossification depending on CCND1. Conclusion: In the process of treating heterotopic ossification with Celecoxib, immune and inflammatory signaling pathways play a significant role. The therapeutic effect of Celecoxib on heterotopic ossification depends on the hub gene CCND1, which plays different roles at different stages of the progression of heterotopic ossification, ultimately inhibiting the occurrence of heterotopic ossification.