Study on the mechanism of Danshen-Guizhi drug pair in the treatment of ovarian cancer based on network pharmacology and in vitro experiment

Our study aims to explore the active components and mechanisms of the Danshen-Guizhi drug pair in treating ovarian cancer by network pharmacology and in vitro experiment. The “component-target-pathway” diagram of the Danshen-Guizhi drug pair was established by network pharmacology, and the effective active components, important targets as well as potential mechanisms of the Danshen-Guizhi drug pair were analyzed. The predicted results were verified by molecular docking and in vitro experiments. The main active components of the Danshen-Guizhi drug pair in the treatment of ovarian cancer are salviolone, luteolin, β-sitosterol and tanshinone IIA. The main core target is PTGS2. The pathways involved mainly include the cancer pathway, PI3K-Akt signaling pathway, and IL-17 signaling pathway. The molecular docking results showed that salviolone and tanshinone IIA had good binding ability to the target. The expression of PTGS2 mRNA and PGE2 in ovarian cells were significantly inhibited by salviolone. The mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer may be regulating cell proliferation, apoptosis and tumor immunity. This provides a theoretical basis for the clinical development and application of the Danshen-Guizhi drug pair.


INTRODUCTION
Ovarian cancer is one of the most common gynecologic malignancies. Its incidence ranks third in female reproductive system malignant tumors, after cervical cancer and uterine cancer (Hunn & Rodriguez, 2012). The early symptoms of ovarian cancer are unobvious, and late symptoms are untypical, the recurrence rate and mortality of ovarian cancer is the first in all kinds of gynecological tumors (Gaona-Luviano, Medina-Gaona & Magana-Perez, 2020). The treatment of ovarian cancer is mainly the comprehensive application of surgery and chemotherapy. However, surgical treatment has relatively large physical damage to patients, and there will be recurrence or serious complications, and platinum drug resistance is also a major problem in the treatment of ovarian cancer (McMullen, Madariaga & Lheureux, 2020).
Many clinical studies have shown that both a single Chinese herbal medicine and Chinese herbal formula for cancer treatment has been proved effective (Mao et al., 2017). It is reported that traditional Chinese medicine for the treatment of ovarian cancer plays an anti-cancer role mainly by affecting the proliferation and apoptosis of cancer cells, generating of tumor vessel and inhibiting the drug resistance of platinum chemotherapeutic drugs . Salvia miltiorrhiza Bunge (Radix Salviae, Danshen) was firstly recorded in Shennong's Herbal Classics, which is a popular component of many traditional Chinese medicines. Previous studies have shown that the effective active ingredients of Danshen play a role in the development of ovarian cancer Yang et al., 2018). The water-soluble components of Danshen can promote blood circulation and remove blood stasis, and fat-soluble components have anti-inflammatory, antibacterial, antioxidant as well as anti-tumor effects (Su et al., 2015). Tanshinone, the fat-soluble component of Danshen, has antiproliferative effects on ovarian cancer cells by promoting apoptosis and down-regulating cisplatin resistance genes (Jiao & Wen, 2011). Neolitsea cassia (L.) Kosterm (Cinnamomi Ramulus, Guizhi) was first found in Newly Compiled Materia Medica of the Tang Dynasty. It can sweat and relieve external syndromes, warm meridians, improve blood circulation, and relieve pain . Guizhi mainly contains cinnamic acid and cinnamic aldehyde, and can inhibit the growth of cancer cells (Reddy et al., 2016). The main active components of Danshen and Guizhi can play a synergistic effect, and improve the efficacy of the two drugs, with effect of promoting blood circulation, dredging collaterals, tonifying qi (vital energy) and nourishing heart. Traditional Chinese medicine believes that the formation of ovarian cancer is due to blood stasis, the pathological mechanism is mainly the disharmony of qi and blood, qi stagnation and blood stasis, and then the formation of mass as proposed by Medical Required (Jiang & Lin, 2006). The Danshen-Guizhi drug pair can improve the blood circulation of patients, balance qi and blood, and remove blood stasis. However, its therapeutic mechanism has not been clearly clarified.
With the rapid development of bioinformatics, network pharmacology has become a new discipline, which mainly reveals the pharmacological mechanism of traditional Chinese medicine. The research and development mode has changed from the traditional "single target, high selective drug development, disease-single target-single drug" to "multi-direction pharmacology, drug action on biological networks, disease-multi-targetsingle drug", which has promoted the modernization of traditional Chinese medicine Zhou et al., 2020). Network pharmacology is mainly to obtain disease-related targets from diseases by text mining, obtain effective components of drugs in traditional Chinese medicine compound for target prediction, construct drug-targetdisease network, and then carry out network analysis to predict the disease mechanism . In this study, the network pharmacology method was used to explore the mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer with multi-component, multi-target and multi-pathway, to provide a theoretical basis for the further development and utilization of the Danshen-Guizhi drug pair in clinical practice.

Screening of effective chemical component and collecting targets of Danshen-Guizhi drug pair
The main effective components of the Danshen-Guizhi drug pair were collected from the TCMSP (https://tcmspw.com/tcmsp.php). All components of the drug pair were screened by TCMSP based on absorption, distribution, metabolism, and excretion (ADME), including oral bioavailability (OB), drug similarity (DL), and lipids-water partition coefficient (AlogP) (Lee et al., 2019). OB is the percentage of oral drugs absorbed into systemic circulation, DL represents the similarity between the active ingredients of Traditional Chinese medicine and drugs, AlogP is one of the binding forces between drugs and receptors, and the occurrence of drug pharmacokinetic process is related to the hydrophobicity of drugs . All components of the Danshen-Guizhi drug pair were screened under the following conditions: oral bioavailability (OB) ≥30%, drug-like property (DL) ≥0.18, and lipid-water partition coefficient (AlogP) ≤5 (Liu et al., 2016). The main active compounds of the Danshen-Guizhi drug pair and their respective protein targets were obtained, the collected targets were transformed into gene names in the Uniprot protein database (https://www.uniprot.org/uploadlists/), and invalid genes were deleted.

Acquisition of disease-related targets
Enter "Ovarian Cancer" in Genecards database (https://www.genecards.org/) and OMIM database (https://omim.org/) to search for potential therapeutic targets of ovarian cancer.

Construction of PPI network and screening of core targets
The ovarian cancer-related gene targets obtained in Genecards were crossed with the drug targets collected in TCMSP to obtain the potential targets for drug treatment of ovarian cancer, and the VENNY diagram was plotted. Then the intersection targets were imported into the STRING (https://www.string-db.org/) database to restrict the species to human ("Homo sapiens"), and the protein-protein interaction network function enrichment analysis was carried out to obtain the PPI network. The results were topologically analyzed by Cytoscape3.8.0 software plug-in for potential therapeutic targets, and PPI network diagram of potential therapeutic targets of the Danshen-Guizhi drug pair for ovarian cancer was constructed.
In order to explore the mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer, the obtained PPI of disease and drug were combined with the intersection network in Cytoscape3.8.0 software, and then the topology network was analyzed by Cytoscape3.8.0 software plug-in to screen the core targets.

Enrichment analysis of Danshen-Guizhi drug pair in the treatment of ovarian cancer
The Matescap database (http://metascape.org/gp/index.html#/main/step1) is a powerful online tool for gene enrichment analysis, which can be used for gene annotation and analysis resources. The potential targets of the Danshen-Guizhi drug pair in the treatment of ovarian cancer were input into the Matescape database, and the species was limited to "H. sapiens". P < 0.01 was set for enrichment analysis of GO molecular function, cell composition and KEGG pathway.

Danshen-Guizhi drug pair "component-target-pathway" network construction
The network diagrams of active components, targets and pathways of drugs were drawn by Cytoscape 3.8.0 software. The network topology parameters of active components and targets of drugs were analyzed by CytoScape 3.8.0 plug-in CytoNCA. The core targets and main active components of drugs were judged according to the network topology parameters of Degree, Betweenness and Closeness.

Molecular docking verification
The main components corresponding to the target protein were searched from the TCMSP database (https://tcmspw.com/tcmsp.php), downloaded and saved as Mol2 format. The "component-target-pathway" network of the Danshen-Guizhi drug pair was analyzed, and the core targets with degree greater than 20 were selected for docking. The corresponding PDB ID of the core targets was found in the PBD (http://www1.rcsb.org/) database and saved as pdb format. Then the target was dehydrated and hydrogenated to optimize ligand in Pymol (Ren et al., 2020). Finally, the molecular docking was performed by AutoDock software to calculate the binding ability of the main active components of the Danshen-Guizhi drug pair to the core targets.

Detection of PTGS2 mRNA expression by qRT-PCR
Molecular docking results showed that the docking results of salvianone with seven main targets were the best. In order to further explore the role of salviolone in ovarian cancer, we detected the effect of salviolone on PTGS2 mRNA expression in ovarian cells. SKOV3 cells and A2780 cells were seeded in a 6-well plate and added with 0 and 20 mm of salviolone for 24 h. Cells were collected, total RNA was extracted, and PTGS2 mRNA was quantitatively detected. Data were represented as the Mean ± SD of three independent experiments.

Detection of PGE2 expression by ELISA
SKOV3 cells were seeded in 6-well plates at 1 × 10 5 cells per well. Cells were incubated in serum-free medium. The cells were incubated with 0, 20 mm salviolone for 2 h. The cell supernatant was collected and PGE2 was measured by ELISA kit (ab133021; Abcam, Cambridge, UK). The absorbance (OD) of the sample was measured at 415 nm. The average value was calculated by three repeated experiments.

Statistical analysis
SPSS Statistics 20.0 and GraphPad Prism v8.0.2 were used for statistical analysis. The data are expressed as Mean ± SD. Student's t-test was used to analyze the difference in effect of luteolin, salviolone, β-sitosterol as well as tanshinone IIA on ovarian cancer cells, compared with the control group, and P < 0.05 was considered to be statistically significant.

Obtaining active components and potential targets of Danshen-Guizhi drug pair
Through searching the TCMSP database, 202 effective active ingredients of Danshen and 220 effective active ingredients of Guizhi were obtained. According to ADME, 54 effective active ingredients of Danshen and six effective active ingredients of Guizhi were obtained. A total of 916 targets of effective components of Danshen and 74 targets of Guizhi were obtained. After deleting duplicate targets, 60 effective active ingredients and 136 targets were obtained, as shown in Table 1. The collected targets were transformed into genes in the Uniprot database.

Acquisition of ovarian cancer-related targets
A total of 8,366 potential targets related to ovarian cancer were obtained by Genecards human gene database. According to the target targets with Score greater than the median, the potential targets for the treatment of ovarian cancer were determined. Finally, 1,119 therapeutic targets for ovarian cancer were obtained (Table S1).

Construction and analysis of PPI network
A total of 119 ovarian cancer-related gene targets obtained in Genecards were crossed with 136 drug targets, and 66 potential targets for drug treatment of ovarian cancer were obtained. VENNY map was drawn, as shown in Fig. 1. The PPI network was constructed by STRING software for 66 targets. The network nodes represent proteins, and the edges represent the combination of proteins and proteins. The number of nodes is 66,1109. Cytoscape 3.8.0 software was used to construct the PPI network diagram of potential targets for the treatment of ovarian cancer by the Danshen-Guizhi drug pair, and the depth of node color as well as the size of node were related to the degree value. Gene targets can be screened from the color of node size: PTGS1, TP53, PTGS2, AKT1, MAPK1, JUN, AR, EGFR, ESR1, etc. These can be considered to be the core targets in the target interaction network of the Danshen-Guizhi drug pair for the treatment of ovarian cancer (Fig. 2). In order to reveal the pharmacological mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer, the PPI of disease (Fig. 3A) and drug (Fig. 3B) were combined to obtain the intersection network (Fig. 3C). The degree value of all nodes was calculated according to the median of two times. All nodes with the degree value greater than two times the median were selected to extract the network, and the Hithubs1  network ( Fig. 3D) was obtained. The median and each attribute value of nodes in the Hithubs network were calculated. All nodes with the degree value greater than the median were selected to extract the network again, and the Hithubs2 network ( Fig. 3E) was obtained. Finally, 70 key targets (Fig. 3F) were obtained.

Enrichment analysis of potential targets for the treatment of ovarian cancer by Danshen-Guizhi drug pair
The potential targets of Danshen-Guizhi drug pair in the treatment of ovarian cancer were input into Matescape database for analysis, and the limited type was "H. sapiens". The personalized analysis was carried out, and the P < 0.01 was set. The enrichment analysis of potential targets of the Danshen-Guizhi drug pair in the treatment of ovarian cancer was obtained, as shown in Fig. 4. The biological processes of GO enrichment include: apoptotic signaling pathway, response to oxidative stress, response to reactive oxygen species, regulation of mitotic cell cycle, cellular response to lipid, response to organic cyclic compound and response to steroid hormone. The results of KEGG enrichment analysis pathway were shown in Table 2. Pathways in cancer, hepatitis B, endocrine resistance, proteoglycans in cancer, AGE-RAGE signaling pathway in diabetes complications, PI3K-Akt signaling pathway, apoptosis, platinum drug resistance, TNF signaling pathway, IL-17 signaling pathway, progesterone-mediated oocyte maturation, and the prolactin signaling pathway were included in Fig. 5. Construction of "component-target-pathway" network of Danshen-Guizhi drug pair The results of effective active components, potential targets and enrichment pathways of the Danshen-Guizhi drug pair were analyzed by Cytoscape3.8.0 built-in plug-in CytoNCA, the core components and core targets of Danshen-Guizhi drug pair in the treatment of ovarian cancer were obtained, as shown in Fig. 6. Cytoscape analysis results showed that the active ingredients DS7 luteolin, DS52 tanshinone IIA, DS48 salviolone, GZ2 beta-sitosterol, DS32 dihydrotanshinone lactone, DS28 cryptotanshinone and so on, among which the highest degree of luteolin was 58, the betweenness centrality was 0.204, and the closeness centrality was 0.476. It can be judged that luteolin was the core component of the Danshen-Guizhi drug pair in the treatment of ovarian cancer, followed by tanshinone (connectivity was 46, betweenness centrality was 0.048, and closeness centrality was 0.422), salviolone (connectivity was 39, betweenness centrality was 0.099, and closeness centrality was 0.414), and β-sitosterol (connectivity is 38, betweenness centrality was 0.083, closeness centrality was 0.431), which can be considered to be effective core components in the treatment of ovarian cancer.

Molecular docking analysis
In AutoDock docking, the score of the binding ability between active ingredients and core targets of the Danshen-Guizhi drug pair was greater than 4.25, indicating that there was a binding ability between active ingredients and targets, greater than 5.25, indicating that the molecule has a good relationship with the target, and greater than 7.0, indicating that the molecule has an excellent binding ability to the target (Hsin, Ghosh & Kitano, 2013).
The results of molecular docking analysis showed that the docking scores of salviolone and tanshinone IIA with the core targets were higher, and the docking scores of PTGS2, AR and PTGS1 with the active components were higher. In the total docking scores, 89% of the total docking scores were higher than 4.25, and 68% of the total docking scores were higher than 5.25, indicating that the predicted targets had good binding ability with the active components. The results of molecular docking were shown in Fig. 7. The higher the score was, the deeper the color was. The highest value was 8.16. Inhibition of core components on ovarian cancer cells The effects of luteolin, salviolone, β-sitosterol as well as tanshinone IIA on SKOV3 and A2780 cells were determined by MTS assay. Compared with the control group, the activity of SKOV3 cells (Fig. 8A) and A2780 cells (Fig. 8B) were significantly inhibited by 10 mm, 20 mm luteolin, salviolone, β-sitosterol as well as tanshinone IIA, as shown in Fig. 8.

Expression of PTGS2 mRNA and PGE2
The expression of PTGS2 mRNA was evaluated by qRT-PCR. The expression of PTGS2 in SKOV3 and A2780 cells were significantly inhibited by salviolone (Fig. 9A). The expression of PGE2 was determined by ELISA. The expression of PGE2 in SKOV3 cells was significantly down-regulated by salviolone (Fig. 9B).

DISCUSSION
This was the first study to explore the main effective components and potential mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer by network pharmacology. The "component-target-pathway" network of the Danshen-Guizhi drug pair was constructed to analyze the effective active components, important targets and potential mechanisms of the Danshen-Guizhi drug pair. Finally, the molecular docking and in vitro study were used to verify the network. Four important active components were screened from the network of "componenttarget-pathway" of the Danshen-Guizhi drug pair, which were luteolin, salviolone, β-sitosterol and tanshinone IIA. Luteolin is a polyphenol flavonoids. Luteolin has antiinflammatory, antiallergic, antioxidant and immuneregulatory pharmacological activities that can be used as antioxidants or prooxidants (Lin et al., 2008;Lopez-Lazaro, 2009; Wang  al., 2018). At the same time, luteolin has a strong role in cell apoptosis, cell cycle regulation, proliferation and differentiation of cancer cells (Tavsan & Kayali, 2019). Salviolone is an antioxidant component extracted from rosemary. Salviolone has many pharmacological effects such as anti-radiation, anti-oxidation, anti-cancer, and inhibiting the migration and invasion of cancer cells (Zhang, Li & Ji, 2018). Beta-sitosterol, a plant derived nutrient, studies have shown that β-sitosterol can interfere with a variety of cell signaling pathways, including apoptosis, cell proliferation and angiogenesis, metastasis and inflammation (Bin Sayeed & Ameen, 2015). Tanshinone IIA is an active component extracted from Danshen. Tanshinone IIA can induce apoptosis of ovarian cancer TOV-21G cells by up-regulating miR-205, which mainly reveals the potential therapeutic application of tanshinone IIA in ovarian cancer . In vitro experiments showed that the activity of SKOV3 and A2780 cells were significantly inhibited by luteolin, salviolone, β-sitosterol and tanshinone IIA. This indicated that the Danshen-Guizhi drug pair might treat ovarian cancer through these components. By analyzing the "component-target-pathway" network of the Danshen-Guizhi drug pair, the main selected targets were PTGS2, HSP90AA1, PTGS1, AKT1 and ESR1. PTGS2 enhanced the proliferation and migration of human ovarian cancer CAOV-3 cells mainly by activating phosphatidylinositol 3 kinase/protein kinase B (PI3-k/Akt) pathway (Gu et al., 2008). It was reported that HSP90AA1 could affect the survival of tumor cells, the invasion and migration of cancer cells, and was closely related to the poor prognosis of tumours. HSP90AA1 RNAi could inhibit the proliferation of ovarian cancer SKOV3 cells and increase its apoptosis (Calderwood et al., 2006;Chu et al., 2013). AKT1 (AKT) is a serine/threonine-specific protein kinase. The activation of PI3K/AKT target caused about 70% of ovarian cancer cells to be activated, which maked many signaling pathways related to cell growth, proliferation, survival, metabolism and angiogenesis overactive (Li, Zeng & Shen, 2014). ESR1 is an important component of estrogen receptor. It is generally believed that ESR1 regulated the expression of protein-coding genes closely related to cell proliferation, survival, differentiation, and participation in the process of cancer (Hu, Chen & Ding, 2018). Therefore, it can be speculated that the Danshen-Guizhi drug pair mainly plays biological functions in the treatment of ovarian cancer through the above targets.
Based on the enrichment results and the analysis of the Danshen-Guizhi drug pair "component-target-pathway" network, it can be considered that the main biological processes are apoptosis, cell proliferation and tumor immunity. The main active components involved were luteolin, salviolone, β-sitosterol and tanshinone IIA.
The pathways involved in apoptosis are mainly apoptosis pathway and PI3K-Akt signaling pathway. PI3K-Akt signal transduction is usually destroyed in human cancer, and AKT1 is an important part of this pathway, which affects multiple processes directly involved in tumorigenesis (Brown & Banerji, 2017). AKT1, CCND1, EGFR, HSP90AA1, MAPK1, RELA and TP53 are involved in the PI3K-Akt signal pathway. The main pathways involved in cell proliferation are PI3K-Akt signaling pathway and proteoglycan pathway in cancer. As an important macromolecule, proteoglycans can participate in the biological process of various cancers and plays an important role. It can affect the progression of tumors by affecting the proliferation, angiogenesis and metastasis of tumor cells (Goldoni et al., 2008). AKT1, CCND1, CASP3, CDKN1A, ESR1, MAPK1, TP53 and other proteins participate in the proteoglycan pathway in cancer. The main pathway involved in tumor immunity is IL-17 signaling pathway. IL-17A is a multi-effect cytokine, and its expression is related to the pathological differentiation level and lymph node metastasis of ovarian cancer patients . CASP3, HSP90AA1, MAPK1, PTGS2, MAPK1 and TP53 are involved in this pathway.
PTGS2 and PGE2 could promote cancer progression by affecting the differentiation and survival of tumor-associated immune cells, and shifting the tumor microenvironment into an immunosuppressive environment (Wang & Dubois, 2010). PTGS2 and PGE2 inhibition may be useful to augment cancer immunotherapy (Chen & Smyth, 2011). In our study, in vitro cell experiments showed that the expression of PTGS2 in SKOV3 cells as well as A2780 cells and PGE2 in SKOV3 cells were significantly inhibited by salviolone. This further indicates the mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer may be inhibiting the expression of PTGS2.

CONCLUSION
In summary, the effective components of the Danshen-Guizhi drug pair in the treatment of ovarian cancer may be salviolone, luteolin, β-sitosterol and tanshinone IIA. The most effective target is PTGS2. The potential mechanism of the Danshen-Guizhi drug pair in the treatment of ovarian cancer may be regulating the proliferation, apoptosis and tumor immunity. This can provide a theoretical basis for the future application of the Danshen-Guizhi drug pair in the treatment of ovarian cancer, and also provide a reference for the clinical application of the Danshen-Guizhi drug pair. Our results need to be validated by further researches.

ADDITIONAL INFORMATION AND DECLARATIONS Funding
This work was supported by foundation projects: National natural science foundation of China (No. 81703780); Science and Technology Project of Henan Province (Nos. 212102310348; 212102310357); and Key Research and Development Program of Hunan Province (Grant No. 2018SK2121). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.