The 3D structures of the chemical ingredients in RXLD were imported into the PharmMapper database to construct the pharmacophore models. The gene targets of GBM were obtained from databases. The pharmacophore-targets network and the protein-protein interactions (PPI) were constructed using the String database and were visualized by using Cytoscape. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses were conducted using Bioconductor software. Cytoscape visualized the relationship of pathways and candidate genes to screen for key target genes. Software packages PyMOL, AutoDock, and Vina acquired the molecular docking results.
We identified 216 potential targets associated with GBM. The core components in RXLD were neochamaejasmin A, wikstrol A, isochamaejasmin, chamaejasmine, and subtoxin A. The undertaken GO enrichment analysis revealed that oxidative stress, cell proliferation, cell cycle, cell invasion, and cell migration were involved in the biological processes. The KEGG enrichment analysis revealed that the crucial pathway was MAPK pathway, while HRAS, PRKCB, MAPK9, CCND1, and TP53 were distributed in core locations. A total of seven RXLD pharmacophores demonstrated strong spontaneous docking activities with MAPK9.
We speculate that the RAS/MAPK pathway might be an upstream pathway through which the RXLD exerts its anti-GBM effects and might be able to regulate further the Wnt/β-catenin, the oxidative stress, and the ferroptosis pathways.