Advanced network pharmacology study reveals multi-pathway and multi-gene regulatory molecular mechanism of Bacopa monnieri in liver cancer based on data mining, molecular modeling, and microarray data analysis

Liver cancer is a malignant tumor that grows on the surface or inside the liver. The leading cause is a viral infection with hepatitis B or C virus. Natural products and their structural analogues have historically made a major contribution to pharmacotherapy, especially for cancer. A list of studies evidences the therapeutic efficacy of Bacopa monnieri against liver cancer, but the precise molecular mechanism is yet to be discovered. This study combines data mining, network pharmacology, and molecular docking analysis to potentially revolutionize liver cancer treatment by identifying effective phytochemicals. Initially, the information on active constituents of B. monnieri and target genes of both liver cancer and B. monnieri were retrieved from literature as well as from publicly available databases. Based on the matching results between B. monnieri potential targets and liver cancer targets, the protein-protein interaction (PPI) network was constructed using the STRING database and imported into Cytoscape for screening of hub genes based on their degree of connectivity. Later, the interactions network between compounds and overlapping genes was constructed using Cytoscape software to analyze the network pharmacological prospective effects of B. monnieri on liver cancer. Gene Ontology (GO) and KEGG pathway analysis of hub genes revealed that these genes are involved in the cancer-related pathway. Lastly, the expression level of core targets was analyzed using microarray data (GSE39791, GSE76427, GSE22058, GSE87630, and GSE112790). Further, the GEPIA server and PyRx software were used for survival and molecular docking analysis, respectively. In summary, we proposed that quercetin, luteolin, apigenin, catechin, epicatechin, stigmasterol, beta-sitosterol, celastrol, and betulic acid inhibit tumor growth by affecting tumor protein 53 (TP53), interleukin 6 (IL6), RAC-alpha serine/threonine protein kinases 1 (AKT1), caspase-3 (CASP3), tumor necrosis factor (TNF), jun proto-oncogene (JUN), heat shot protein 90 AA1 (HSP90AA1), vascular endothelial growth factor A (VEGFA), epidermal growth factor receptor (EGFR), and SRC proto-oncogene (SRC). Through, microarray data analysis, the expression level of JUN and IL6 were found to be upregulated while the expression level of HSP90AA1 was found to be downregulated. Kaplan–Meier survival analysis indicated that HSP90AA1 and JUN are promising candidate genes that can serve as diagnostic and prognostic biomarkers for liver cancer. Moreover, the molecular docking and molecular dynamic simulation of 60ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between the compound and binding pockets of HSP90AA1 and JUN. Despite that, in vivo and in vitro studies are mandatory to unveil pharmacokinetics and biosafety profiles to completely track the candidature status of B. monnieri in liver cancer.