ORIGINAL RESEARCH article
Exploring the Regulation Mechanism of Xihuang Pill, Olibanum and β-Boswellic Acid on the Biomolecular Network of Triple-Negative Breast Cancer Based on Transcriptomics and Chemical Informatics Methodology
- 1Galactophore Department, The First Hospital of Hunan University of Chinese Medicine, Changsha, China
- 2Department of Cardiac Surgery, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
- 3Graduate College, Capital Medical University, Beijing, China
- 4Graduate College, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
- 5School of Medicine, Hunan University of Chinese Medicine, Changsha, China
- 6Department of Geratology, Xiyuan Hospital, China Academy of Chinese Medical Sciences, Beijing, China
- 7School of Clinical Medicine (Xiyuan Hospital), Beijing University of Chinese Medicine, Beijing, China
Background: Xihuang Pill (XHP) is mainly used to treat “Ru Yan (breast cancer)”. Evidence-based medical evidence and showed that XHP improves the efficacy of chemotherapy and reduced chemotherapy-induced toxicity in breast cancer patients. However, the mechanism of XHP against breast cancer is not clear.
Methods: The effect of XHP extract on cell half-inhibitory concentration (IC50) and cell viability of MD-MB-231 cells was detected by CCK-8 method. The cell inhibition rate of MDA-MB-453 cells were detected by MTT method. Apoptosis was detected by flow cytometry, cell transfer ability was detected by Transwell method, and cell proliferation ability was detected by colony formation assay. The expression of Notch1, β-catenin and c-myc mRNA in MDA-MB-453 cells were detected by real-time fluorescence quantitative PCR. Then, chemical informatics and transcriptomics methodology was utilized to predict the potential compounds and targets of XHP, and collect triple negative breast cancer (TNBC) genes and the data of Olibanum and β-boswellic acid intervention MD-MB-231 cells (from GSE102891). The cytoscape software was utilized to undergo network construction and network analysis. Finally, the data from the network analysis was imported into the DAVID database for enrichment analysis of signaling pathways and biological processes.
Results: The IC50 was 15.08 g/L (for MD-MB-231 cells). After interfering with MD-MB-231 cells with 15.08 g/L XHP extract for 72 h, compared with the control group, the cell viability, migration and proliferation was significantly decreased, while early apoptosis and late apoptosis were significantly increased (P < 0.01). After interfering with MDA-MB-453 cells with 6 g/L XHP extract for 72 h, compared with the control group, the cell inhibition and apoptosis rate increased, while the expression of Notch1, β-catenin and c-myc mRNA decreased. (P < 0.05). The chemical informatics and transcriptomics analysis showed that four networks were constructed and analyzed: (1) potential compounds-potential targets network of XHP; (2) XHP-TNBC PPI network; (3) DEGs PPI network of Olibanum-treated MD-MB 231 cells; (4) DEGs PPI network of β-boswellic acid -treated MD-MB 231 cells. Several anti-TNBC biological processes, signaling pathways, targets and so on were obtained.
Conclusion: XHP may exert anti-TNBC effects through regulating biological processes, signaling pathways, targets found in this study.
Breast cancer is the leading cause of death among women worldwide, and the incidence rate has increased significantly in recent years, which seriously threatens women’s health (Fan et al., 2014; Desantis et al., 2017). Breast cancer is currently divided into five subtypes by coding sequence microarray technology (Prat et al., 2015; Cejalvo et al., 2018): (1) Luminal-A type; (2) Luminal-B type; (3) human epidermal growth factor receptor 2 (HER2) overexpression type; (4) base-like type; (5) normal type. Basal-like breast cancer is non-specific invasive ductal carcinoma (Jiang et al., 2019), and its ER, PR and HER2 are negative, also known as triple negative breast cancer (TNBC). At present, TNBC accounts for 10 to 17% of all breast cancers (Warner et al., 2015; Kulkarni et al., 2019). The vast majority of TNBCs are highly invasive ductal carcinomas with nuclear polymorphism, high mitotic rate, and minimal tubule formation (Kulkarni et al., 2019). TNBC is generally classified into basal cell-like type 1, basal cell-like type 2, immunoregulatory, interstitial and mesenchymal stem cell types, and luminal androgen receptor type (Shi et al., 2018; Jiang et al., 2019). TNBC is a poorly differentiated tumor with strong invasive and metastatic ability, easy to invade blood vessels, and increased recurrence rate (Jiang et al., 2019).
At present, the management and treatment measures for triple-negative breast cancer are mainly: (1) local treatment: surgery is still the first choice for local treatment (Li et al., 2019a); (2) systemic therapy: combination chemotherapy of taxane and anthracycline is currently the common choice for TNBC neoadjuvant chemotherapy, but anthracyclines have irreversible toxicity to the heart (Mangone et al., 2019); (3) targeted therapy (Jhan and Andrechek, 2017). However, due to the ineffectiveness of traditional endocrine therapy and targeted therapy, patients with TNBC sometimes have tumor metastasis very early, which seriously affects their physical and mental health (McCann et al., 2019). At present, in order to find new chemotherapy or sensitizing drugs, plants and natural products are gradually becoming the source of new TNBC drug development (Szarc Vel Szic et al., 2017). The current study found that some traditional Chinese medicine compounds and natural medicines can inhibit the proliferation, metastasis, and drug resistance of TNBC cells in a variety of ways (Baraya et al., 2017; Szarc Vel Szic et al., 2017).
Xihuang Pill (XHP) is from the Wai Ke Quan Sheng Ji by Wang Weide in 1740, which is mainly used to treat “Ru Yan (breast cancer)” and so on. XHP is compose of Myrrha, Bovis Calculus, Olibanum and Moschus. Systematic reviews and meta-analysis showed that XHP combined with chemotherapy significantly enhanced tumor response in breast cancer patients, improved Karnofsky performance scores and reduced chemotherapy-induced toxicity (Guo et al., 2018; Mao et al., 2019). He et al. found that XHP-containing serum increased TP53 and Bax (P < 0.05), and decreased the ratio of Bcl-2/Bax in MDA-MB-435 cells (He et al., 2018). The mechanism of anti-TNBC of XHP has been reported in many studies, such as improving the immunosuppressive state of the tumor microenvironment and reversing immune escape, thereby inhibiting tumor growth. XHP reduces the number of Treg cells by inhibiting the expression of PI3K and AKT and up-regulating the expression of AP-1 in Treg cells, thereby promoting Treg cell apoptosis (Li et al., 2018a). Other study found that the mechanism of XHP inhibition of tumors may be related to the up-regulation of gene and protein expression of MEKK1, SEK1, JNK1 and AP-1 in Treg cells in the tumor microenvironment (Su et al., 2018). Zheng et al. found that XHP can block the cell cycle of the Hs578T cell line and promote its apoptosis (Zheng et al., 2016).
Recent studies showed that the main anti-breast cancer herbs in XHP are Myrrha and Olibanum, especially Olibanum (Cheng et al., 2016; Hao et al., 2018). Although the above studies have described some of the mechanisms of XHP against TNBC, the mechanism remains unclear. In our previous studies, we successfully used multiple bioinformatics techniques and transcriptomics to analyze the mechanisms by which traditional Chinese Medicine interferes with different types of breast cancer (Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019). Therefore, this study will use a multi-directional pharmacology strategy based on chemical informatics and transcriptomics to clarify the mechanisms by which XHP and Olibanum treat TNBC. The research process is shown in Figure 1.
Material and Methods
Experimental Material Preparation
Xihuang Pill (XHP) was purchased from Tianjin Tianshili (Liaoning) Pharmaceutical Co., Ltd. (Batch number: 20140726; Specification: 0.1g * 30 bottles/box; The composition ratio of Moschus, Bovis Calculus, Myrrha and Olibanum is 15: 15: 550: 550). Reference substance: acetyl-11-keto-β-boswellic acid (batch number: 111760-201502, mass fraction >98%) was purchased from China Food and Drug Research Institute.
Preparation of XHP solution: XHP was immersed in DMEM medium pre-cooled at 4°C for 24 h (mass concentration 0.1 g/mL) in a sterile sealed container; use ultrasonic vibration to help dissolve for 2 h and continue to soak for 48 h at 4°C. The supernatant was filtered through a 0.22 μm micropore filter to obtain an XHP leaching solution. The XHP solution is stored at 4°C (or −20°C); during the experiment, it was diluted to the desired concentration with DMEM medium.
Preparation of XHP solution required for High Performance Liquid Chromatography (HPLC): 1.00 g of XHP powder was accurately weighed and placed in a 50 ml Erlenmeyer flask. Pipet 20 ml of methanol accurately, sonicate in an ice bath for 20 min, extract twice, and place at room temperature. Centrifuge at 4,000 r/min for 5 min. The supernatant was placed in a pear-shaped bottle and concentrated under reduced pressure. Reconstitute with methanol and transfer to a 25 ml volumetric flask. Finally, make up to volume with methanol and shake well.
Preparation of acetyl-11-keto-β-boswellic acid reference substance: Take an appropriate amount of acetyl-11-keto-β-boswellic acid, accurately weigh it, place it in a measuring flask, add methanol to volume, and make it to a mass concentration of 1.22 mg/mL.
Human triple-negative breast cancer (TNBC) cell line MDA-MB-231 and MDA-MB-453 were provided by the Cell Center of Xiangya School of Medicine, Central South University. While experimenting, the MD-MB-231 cells and MDA-MB-453 cells were cultured in high glucose DMEM medium containing 10% FBS in an incubator at 37°C, 5% CO2, and the medium was changed every other day. The group without XHP was the control group, and the group with XHP extract was the experimental group (XHP group), and each group had three duplicate wells.
Reagent and Instrument
DMEM medium and Transwell kit were purchased from Corning Inc.; Fetal bovine serun (FBS) was purchased from Ausbian Inc., Australia; CCK-8 kit and Giemsa dye solution were purchased from Sigma Inc., USA; The Annexin V-FITC/PI double-stained cell apoptosis assay kit was purchased from eBioscience Inc., USA.
Ultra-clean workbench (Sujing Antai Company), CO2 incubator (Hitachi Company). LC-20A HPLC, including SPD-M20A detector, DGU-20A quaternary pump (Shimadzu company); KQ-300DE ultrasonic cleaner (Kunshan Ultrasonic Instrument Co., Ltd.)
Shimadzu LC-20A HPLC, Waters Symmetry C18 column (250 mm × 4.6 mm, 5 μm); injection volume 20 μl; column temperature 30°C; volume flow 0.8 ml/min; mobile phase is methanol-0.5% acetic acid aqueous solution.
The elution procedure is 0 to 10 min, 75% methanol; 10 to 45 min, 75 to 87% methanol; 45 to 75 min, 87% methanol; 75 to 85 min, 87 to 90% methanol; 85 to 115 min, 90% methanol; 115 to 120 min, 90 to 97% methanol; 120 to 150 min, 97% methanol. The detector is an evaporative light scattering detector, the drift tube temperature is 45°C; the air volume flow is 1.5 L/min, and the gain value is 2.
Experimental Methods for MD-MB-231 Cell
Detection of XHP Half-Inhibitory Concentration (IC50) and Cell Viability
Some 96-well plates were seeded at 3,000 cells (100 μl) per well. Nine (9) groups were set according to the concentration of added XHP, and the concentrations of XHP were 0, 1, 2, 5, 10, 20, 50, 75 and 100 g/L. After treatment for 72 h, 10 μl of CCK-8 reagent was added 2 to 4 h before the termination of culture, and the OD value was detected by a microplate reader at 450 nm. After logarithmic processing, a scatter plot was prepared to calculate the IC50 value of XHP. Using this IC50 value as the drug concentration of XHP intervention in MD-MB-231 cells in subsequent experiments (including cell viability detection).
After treatment with XHP extract for 72 h, the culture was continued for 5 d with the control group. 10 μl of CCK-8 reagent was added 2 to 4 h before the termination of the culture, and the cell viability was measured daily for 5 days via the above procedure.
MD-MB-231 Cell Apoptosis Detection by Flow Cytometry
Some 6-well plates were seeded at a minimum of 5 × 105 cells per well (2 ml) and plated for 24 h. Apoptosis detection was performed by flow cytometry according to the procedure in the instructions of Annexin V-FITC/PI double-stained cell apoptosis assay kit.
The cells were cultured in 24-well plates. About 100 μl of serum-free medium was added to each well of the Transwell’s inner chamber, and 600 μl of medium containing 30% FBS was added to each well of the Transwell’s outer chamber, and plated for 18 h at a cell number of 1 × 105 per well. Transfer cells, fix, stain with Giemsa stain, photograph with fluorescence microscopy, and count cells at a magnification of 200×.
Cell Clone Detection
Some 6-well plates were seeded at 800 cells (2 ml) per well, and continue to culture for 10 days, and change the solution once every 3 days. Cell clones were photographed under a fluorescent microscope before termination of the experiment. Fix the cells with 4% paraformaldehyde, crystal violet staining, photograph.
Experimental Methods for MDA-MB-453 Cell
Determination of Cell Inhibition Rate by MTT Method
The MDA-MB-453 cells in the logarithmic growth phase were used for experiments. The cell concentration was adjusted to 1 × 105 cells/mL and inoculated in 96-well culture plates; 90 μl of cell suspension was added to each well, and then 10 μl of different concentrations of XHP (0, 4, 6, 8, 10, 12, 14 g/L) were added. Six (6) duplicate wells were set for each group, and after incubating in a cell incubator at 37°C and 5% CO2 for 72 h, 20 μl MTT (5mg/mL) was added to each well. After continuing the culture for 4 h, the supernatant was discarded, and DMSO solution (DMSO) 150 μl/well was added to each well. After shaking for 10 min to fully dissolve the crystals, the OD value of each well was measured with a microplate reader at a wavelength of 490 nm.
MDA-MB-453 Cell Apoptosis Detection by Flow Cytometry
Apoptosis detection was performed by flow cytometry according to the procedure in the instructions of Annexin V-FITC/PI double-stained cell apoptosis assay kit.
The Expression of Notch1, β-catenin and c-myc mRNA Detection by Real-Time Fluorescence Quantitative PCR
The total RNA of each group of cells was extracted with Trizol according to the kit instructions. The OD260/OD280 ratio is calculated, and the ratio ≥1.8 means the purity and concentration of the RNA meet the experimental requirements. After detecting the integrity of RNA by agarose gel electrophoresis, the primers of Notch1, β-catenin, c-myc and internal reference GAPDH were added to amplify the corresponding target fragments. Finally, the real-time fluorescence quantitative PCR program was performed on the machine. The primers were designed by Primer 3.0 software and synthesized by Yuantai Bio-Technology Co., Ltd., see Table 1.
Statistical analysis was performed using SPSS 19.0 software. The measurement data were normally distributed, expressed as mean ± standard deviation, and t-test was used for comparison between groups. The difference was statistically significant at P < 0.05.
Chemical Informatics Methods
XHP Active Compounds Prediction
With the development of computer technology, Chinese medicine-related laboratories have built several large-scale Chinese medicine databases, which contain the components of commonly used Chinese medicines. Traditional Chinese Medicine Systems Pharmacology Database (TCMSP) (http://tcmspw.com/tcmsp.php) (Chen et al., 2014), Traditional Chinese Medicine Database@Taiwan (TCM@Taiwan) (http://tcm.cmu.edu.tw/zh-tw/) (Ru et al., 2014), Traditional Chinese Medicines Integrated Database (TCMID) (http://126.96.36.199:8000/tcmid/) (Lin et al., 2018) is a commonly used database. Oral bioavailability (OB), Caco-2 permeability and drug-likeness (DL) were utilized to identify the potential bioactive compounds of XHP (Walters and Murcko, 2002; Ano et al., 2004; Hu et al., 2009; Xu et al., 2012; Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019). The compounds with OB ≥30%, Caco-2 > −0.4 and DL ≥0.18 were regard as oral absorbable compounds with biologically active (Walters and Murcko, 2002; Ano et al., 2004; Hu et al., 2009; Xu et al., 2012; Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019). Finally, a lot of compounds were collected: (13E,17E,21E)-8-hydroxypolypodo-13,17,21-trien-3-one, (13E,17E,21E)-polypodo-13,17,21-triene-3,18-diol, (16S, 20R)-dihydroxydammar-24-en-3-one, (20R)-3β-acetoxy-16β-dihydroxydammar-24-ene, (20S)-3β,12β,16β,25-pentahydroxydammar-23-ene, (20S)-3β-acetoxy-12β,16β,25-tetrahydroxydammar-23-ene, (3R,20S)-3,20-dihydroxydammar- 24-ene, (8R)-3-oxo-8-hydroxy-polypoda -13E,17E,21-triene, 11α-hydroxypregna-4,17(20)-trans-diene-3,16-dione, 15α-hydroxymansumbinone, 16-hydroperoxymansumbin-13(17)-en-3β-ol, 1α-acetoxy-9,19-cyclolanost-24-en-3β-ol, 28-acetoxy-15α-hydroxymansumbinone, 2-methoxyfuranoguaia-9-ene-8-one, 35833-62-6, 3-methoxyfuranoguaia-9- en-8-one, 3β- hydroxydammar-24-ene, 3β-acetoxy-16β,20(R)-dihydroxydammar-24-ene, 4,17(20)-(cis)-pregnadiene-3,16-dione, 7β,15β- dihydroxypregn-4-ene-3,16-dione, beta-Sitosterol, Cabraleadiol monoacetate, Cabraleone, Chondrillasterol, Diayangambin, Epimansumbinol, Guggulsterol IV, Guggulsterol VI, Guggulsterone, Isofouquierone, Mansumbin-13(17)-en- 3,16-dione, Mansumbinoic acid, MOL001019, MOL001164, Myrrhanol C, Myrrhanone A, Myrrhanones B, Naringenin, Pelargonidin, Petunidin, Phellamurin, Quercetin, Stigmasterol, 3-oxo-tirucallic acid, Acetyl-alpha-boswellic acid, alpha-Boswellic acid, beta-Boswellic acid, Incensole, O-acetyl-α-boswellic acid, Phyllocladene, Tirucallol.
Since the application of biological models to predict XHP compounds has limitations (Metodiewa et al., 1997), in order to avoid missing active compounds during the pre-screening process, we searched a large number of references and selected oral absorbable compounds with pharmacological activity. Combined with relevant references (Huwez et al., 1998; Chen et al., 2013; Morikawa et al., 2017; Kunnumakkara et al., 2018; Hao, 2018), the following compounds are included: beta-Elemene, delta-Elemene, gamma-Elemene, Limonene, Curzerene, Germacrene B, Myrcenol, beta-selinene, Spathulenol, β-elemonic acid, 11-keto-β-boswellic acid, acetyl-11-keto-β-boswellic acid, trans-Anethole, ZINC01280365, SCHEMBL17727656, Oleanolic acid, Methyl deoxycholate, Deoxycholic Acid, Cholesterol, Cherianoine, Chenodeoxycholic acid, n-Nonane, Muscopyridine, Muscol, Estragole, Dihydroxanthyletin, Dihydroagarofuran, Cholesteryl ferulate, Androst-4-ene-3,17-dione, Allantoin.
Potential Targets of XHP
Four databases were utilized to predict the potential targets of XHP: Similarity ensemble approach (SEA) (http://sea.bkslab.org/) (Keiser et al., 2007), STITCH Database (http://stitch.embl.de/) (Kuhn et al., 2008; Kuhn et al., 2012), Swiss Target Prediction (http://www.swisstargetprediction.ch/) (Gfeller et al., 2014), and PubChem (https://pubchem.ncbi.nlm.nih.gov/) (Kim S. et al., 2019). The UniProtKB (http://www.uniprot.org/) was used for the correction of protein’s names and the collection of official symbols with the species limited to “Homo sapiens”. The details are described in Table S1 (see Supplementary Material).
Triple Negative Breast Cancer Biological Network
To construct the biological network of TNBC, the TNBC-related genes were collected from OMIM database (http://omim.org/) and Genecards (http://www.genecards.org) (Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019). Finally, one thousand and two hundred and twenty-one (1220) TNBC-related genes were obtained. These TNBC-related genes will be used for subsequent biological network construction and network analysis (see Table S2).
Human Transcriptomics Data
Network Construction and Analysis Methods
The data of protein-protein interaction (PPI) of XHP targets and TNBC genes were obtained from the String database (http://string-db.org/) with the species limited to “Homo sapiens” (Missiuro et al., 2009; Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019). Then, these data were input into cytoscape software ver 3.7.0 (https://cytoscape.org/) for network construction (Bader and Hogue, 2003).
The networks were analyzed by the plugin MCODE to obtain cluster. The definition and the methodology of acquisition of clusters were described in our previous work (Zeng and Yang, 2017; Yang et al., 2018; Yang et al., 2019), such as “Exploring the pharmacological mechanism of Yanghe Decoction on HER2-positive breast cancer by a network pharmacology approach” (Zeng and Yang, 2017) and “Investigating the regulation mechanism of baicalin on triple negative breast cancer’s biological network by a systematic biological strategy” (Yang et al., 2019).
Results and Discussion
XHP and acetyl-11-keto-β-boswellic acid reference substance were analyzed and compared according to the steps of HPLC Detection. The fingerprint was shown in Figure S1.
Experimental Results for MD-MB-231 Cell
Inhibition Effect of XHP on the Growth of MD-MB-231 Cell
Different concentrations of XHP were applied to MD-MB-231 cells for 72 h, and the OD value of each group was detected by CCK-8 kit. The OD value could indirectly reflect the number of viable cells. The results are shown in Table 2. The XHP extract concentration value was logarithmically processed. The OD value of each concentration group was compared with the OD value of the non-medicated solvent group, and the cell inhibition rate was calculated. Draw a scatter plot with Log (concentration value) and cell inhibition rate (Figure 2). The IC50 value of XHP extract intervention in MD-MB-231 cells for 72 h was 15.08 g/L.
Effect of XHP on the Viability of MD-MB-231 Cells
The CCK-8 reagent contains WST-8, which is reduced in the cell mitochondria to a highly water-soluble yellow formazan product, the number of which is proportional to the number of viable cells. The OD value measured at a wavelength of 450 nm can indirectly reflect the cell proliferation. The cells were intervened for 72 h with XHP extract at a final concentration of 15.08 g/L. Compared with the control group, the cell viability was significantly decreased (P < 0.01), as shown in Figure 3.
Effect of XHP on apoptosis of MD-MB-231 cells
The cells were intervened for 72 h with XHP extract at a final concentration of 15.08 g/L. Compared with the control group, the early apoptosis and late apoptosis of the cells after XHP intervention were significantly increased (P < 0.01) (Table 3 and Figure 4).
Figure 4 Effect of XHP on apoptosis of MD-MB-231 cells (A, C) are apoptotic maps of a random set of three duplicate wells in the control group; (B, D) are apoptotic maps of a random set of three duplicate wells in the XHP group; (E) is a histogram of the proportion of apoptotic cells in the two groups, compared with the control group; *P < 0.01).
Effect of XHP on Migration Ability of MD-MB-231 Cells
The cells were intervened for 72 h with XHP extract at a final concentration of 15.08 g/L. The effect of XHP on cell metastatic ability was determined by observing the migration of cells in the Transwell chamber to serum-containing media. Compared with the control group, the migration ability of the cells after XHP intervention were significantly decreased (P < 0.01) (Figure 5).
Figure 5 Effect of XHP on migration ability of MD-MB-231 cells (A, B) are the Giemsa staining maps of a random set of three duplicate wells in the control group and the XHP group, respectively; (C) is a histogram of the cell counts of the two groups. Each group of cells was randomly counted for nine fields of view; and take the average of the number of cells; compared with the control group, *P < 0.01.
Effect of XHP on Cell Cloning Ability
The cells were intervened for 72 h with XHP extract at a final concentration of 15.08 g/L. After 10 days of continuous culture, the clone formation rate was calculated to quantitatively analyze the proliferation ability of breast cancer cells. Compared with the control group, the cloning ability of the cells after XHP intervention were significantly decreased (P < 0.01) (Figure 6).
Figure 6 Effect of XHP on cell cloning ability. (A) is the 6-well plate culture plates after staining with crystal violet; (B, C) are the cell clones of a random set of three duplicate wells in the XHP group and the control group, respectively; (D) is the histogram of the number of cell clones in both groups; compared with the control group, *P < 0.01.
Experimental Results for MDA-MB-453 Cell
Inhibition Effect of XHP on the Growth of MDA-MB-453 Cell
The results of MTT also showed that XHP had a proliferation-inhibiting effect on MDA-MB-453 cells in a dose-dependent manner. When the concentration of XHP is lower than 6g/L, the inhibition rate of MDA-MB-453 cells is ≤20%, and there is no cytotoxicity. Therefore, 6 g/L was selected as the non-cytotoxic concentration of XHP for subsequent MDA-MB-453 experiments. The results are shown in Figure 7.
Figure 7 Inhibition rates of XHP on MDA-MB-453 Cell (* Compared with control group (0 g/L), P < 0.05).
Effect of XHP on apoptosis of MDA-MB-453 Cell
After 72 h intervention with 6 g/L XHP, the results of flow cytometry showed that the apoptosis rate in the test group was significantly higher than that in the control group (normal saline serum group) (P < 0.05) (Figure 8).
Figure 8 Effect of XHP on Apoptosis of MDA-MB-453 Cell. (A) apoptotic maps of a random set of three duplicate wells in the control group; (B) are apoptotic maps of a random set of three duplicate wells in the XHP group; (C) histogram of the proportion of apoptotic cells in the two groups, compared with the control group; *P < 0.05.
Effect of XHP on the Expression of Notch1, β-Catenin and c-myc mRNA
Compared with the control group (saline saline group), the expressions of Notch1 mRNA, β-catenin mRNA and c-myc mRNA in the XHP group were significantly reduced (P < 0.01) (Figure 9).
This experiment showed that XHP can promote apoptosis, inhibit cell proliferation, metastasis and vitality, and can down-regulate the expression of Notch1 mRNA, β-catenin mRNA and c-myc mRNA at the cellular level. In addition, other studies have also demonstrated the anti-TNBC effect of XHP. Zheng et al. (2016) found that XHP significantly inhibited the viability of Hs578T cell line in a dose- and time-dependent manner. The mechanism may be that XHP induced apoptosis through the inherent Bcl-2 dependent pathway and cell cycle arrest. Su et al. (2018) found that XHP may promote Treg cell apoptosis in the tumor microenvironment and further inhibit the tumor growth of breast cancer in 4T1 mice. The mechanism may be that XHP up-regulated the MEKK1, SEK1, JNK1 and AP-1 gene and protein expression in Treg cells in the tumor microenvironment. Li et al. (2018a) further research indicates that XHP promotes apoptosis of Treg cells by inhibiting PI3K/AKT/AP-1 signaling pathway, thereby reducing the number of Treg cells, weakening the immunosuppressive state of the tumor microenvironment, reversing the mainstream immune escape, and inhibiting the growth of tumors. However, the molecular mechanisms of XHP’s anti-TNBC effects, such as regulatory pathways and targets, need further study. Based on this, this study will further study the mechanism of XHP reversing the multidrug resistance of tumors, and provide a theoretical basis for the development of new drugs and clinical combination drugs.
XHP-TNBC PPI Network Analysis
XHP’s Potential Targets and TNBC Genes
After the potential target prediction, a total of 1,178 potential targets were obtained. Myrrha contains 550 potential targets; Bovis Calculus contains 680 potential targets; Olibanum contains 390 potential targets; Moschus contains 845 potential targets. Some of the targets contained in different herbs overlap (Figure S2). In the outer circle, red, blue, green, purple stand for potential targets of Bovis Calculus, Moschus, Myrrha, Olibanum, respectively. In the inner circle, the greater the number of purple links and the longer the dark orange arc, the more overlap between the input target lists. The blue link indicates the amount of functional overlap between the input target lists.
The relationship among potential compounds and potential targets was shown in Figure 10. This network consists of 81 potential compounds, 1,175 potential targets and 6,057 edges. The targets near the center are regulated by more compounds than ones in the peripheral. For example, the targets in the center are: AR (55 edges), CYP19A1 (52 edges), ESR1 (36 edges), NR3C1 (36 edges), ESR2 (35 edges), PTPN1 (35 edges), MAPT (34 edges), HSD11B2 (31 edges), NR1H4 (29 edges); the targets in the peripheral (AXL, CAMK2B, CXCR1, NEK2, NEK6) are regulated only by one compound. In addition, some important targets related to TNBC are regulated by more XHP compounds: EGFR (four edges), PARP1 (eight edges), PARP10 (three edges), PARP14 (one edge), PARP15 (one edge), PARP2 (two edges), VEGFA (two edges), PIK3CA (11 edges), AKT1 (eight edges), AKT2 (one edge). After searching, a total of 1,221 TNBC-related genes were obtained. These TNBC genes will be combined with XHP’s predicted targets to construct an XHP-TNBC PPI network so as to observe the association between XHP and TNBC.
Figure 10 Potential compounds-potential targets network of XHP (The red circle stands for compound of Olibanum; the blue circle stands for compound of Moschus; the yellow circle stands for compound of Bovis Calculus; the orange circle stands for compound of Myrrha. The pink hexagon stands for potential targets.)
It can be seen from the potential compounds-potential targets network of XHP that the herbs in XHP have multi-component and multi-target effects in anti-TNBC. Many compounds can act on one or more targets at the same time, and some targets can be regulated by multiple compounds. They may be the main active ingredient and target of XHP against TNBC. Currently, the study found that TNBC can also be subdivided into six subtypes: two basal-like-related (BL1 and BL2), two mesenchymal related subtypes [mesenchymal (M) and mesenchymal stem-like (MSL)], an immunoregulatory subtype (IM) and a tubular androgen receptor type (LAR) (Waks and Winer, 2019). LAR subtype is sensitive to androgen receptor (AR) inhibitors due to high expression of AR. Enzalutamide (an AR inhibitor) is an advanced prostate cancer drug approved by the US FDA (Traina et al., 2018). The current study on Enzalutamide is a phase II clinical study of patients with locally advanced or metastatic AR+TNBC; The study showed that approximately 55% of patients with TNBC were AR-positive, and enzalutamide had a good effect on the LAR subtype TNBC (Traina et al., 2018).
Recent studies have shown that tamoxifen may be effective against certain subtypes of TNBC, which is associated with ESR2 (Mukhopadhyay et al., 2019). Studies at the Roswell Park Comprehensive Cancer Center in the United States have shown that TP53 status is a determining factor in the duality of estrogen receptor-beta (ESR2) function (Mukhopadhyay et al., 2019). ESR2 and mutant TP53 can be combined to predict survival in patients with TNBC (Mukhopadhyay et al., 2019). Current research shows that miRNA genetic variation is associated with the expression of important receptors such as ER and HER2, and the single nucleotide polymorphism (SNP) site present in the estrogen receptor alpha gene (ESR1) may be involved in the development and progression of TNBC (Zhang, 2011).
Overall, AR, vascular endothelial growth factor (VEGF), poly(ADP-ribose) polymerase (PARP) and epidermal growth factor receptor (EGFR), PI3K, AkT, microRNAs and lncRNAs and so on are potential therapeutic targets for TNBC (Nagini, 2017; Chan et al., 2019 Nakhjavani et al., 2019).
XHP-TNBC PPI Network
The relationship among XHP’s potential targets and TNBC genes were shown in XHP-TNBC PPI network. This network is composed of 247 XHP-TNBC targets, 907 XHP targets, 830 TNBC targets, and 73,021 edges. The top 20 targets of high-degree are selected and divided into three categories: (1) XHP targets: ALB (545 edges), CASP3 (420 edges); (2) TNBC genes: MYC (562 edges), EGF (468 edges), FN1 (460 edges), PTEN (446 edges); (3) XHP-TNBC: TP53 (743 edges), AKT1 (722 edges), GAPDH (662 edges), INS (608 edges), EGFR (594 edges), IL6 (530 edges), VEGFA (525 edges), MAPK3 (513 edges), SRC (460 edges), TNF (456 edges), HRAS (454 edges), MAPK1 (452 edges), STAT3 (441 edges), JUN (427 edges) (Figure 11).
Figure 11 XHP-TNBC PPI Network (Purple circle stands for XHP-TNBC; Blue circle stands for TNBC genes; Pink circle stands for XHP targets.)
Biological Processes of XHP-TNBC PPI Network
The XHP-TNBC PPI network were analyzed by MCODE to obtain the clusters. The clusters of this network were shown in Figure S3 and Table 4. The genes and targets in clusters 1–10 were put into DAVID database to undergo GO enrichment analysis as an example.
After the GO enrichment analysis, a lot of biological processes of each cluster were return. Cluster 1 is associated with chemokines and their receptor-mediated signaling pathways, immune responses, ERK1/2 signaling pathways, angiogenesis, T cell chemotaxis. Cluster 2 is involved in apoptosis and cell proliferation, hypoxia induction, cell cycle, estrogen-mediated biological processes, canonical Wnt signaling pathway, ERK1/2 signaling pathway, angiogenesis, T cell-mediated immune response. Cluster 3 is involved in ERK1/2 signaling pathway, MAPK signaling pathway, PI3K signaling pathway, angiogenesis, and inflammatory response. Cluster 4 is involved in cell migration and adhesion to angiogenesis, extracellular matrix, hypoxia-induced, PI3K signaling pathway, MAPK signaling pathway, ERK1/2 signaling pathway, and endoplasmic reticulum stress. Cluster 5 is associated with tumor necrosis factor-mediated signaling, apoptosis, cell cycle, angiogenesis, T cell receptor signaling pathway, Wnt signaling pathway, and endoplasmic reticulum-mediated endogenous apoptotic signaling pathway. Cluster 6 is involved in transcriptional regulation, apoptosis, cell matrix adhesion, T cell receptor signaling pathway, angiogenesis, Wnt signaling pathway, cell cycle, immune cell chemotaxis and T cell costimulation. Cluster 7 is associated with apoptosis, extracellular matrix, negative regulation of Wnt signaling pathway, NF-κB signaling pathway, and angiogenesis. Cluster 8 is related to redox, steroid metabolism, androgen and estrogen-mediated biological processes, oxidative stress. Cluster 9 is associated with redox, chemokine generation and mediated biological processes, steroid hormone-mediated signaling, angiogenesis, T cell selection, cell migration. Cluster 10 is involved in the transmission of intracellular signals. The details of each cluster and biological process was described in Table S3.
Since Cluster 1 contains many classic biological processes, bubble chart is created using the main biological process data contained in Cluster 1 (Figure 12A).
Figure 12 Bubble chart (A) biological processes; (B) signaling pathway. X-axis stand for fold enrichment.
Signaling Pathways of XHP-TNBC PPI Network
The XHP targets combining with TNBC genes were put into DAVID database for pathway enrichment analysis. After this, thirty-four (34) anti-TNBC-related signaling pathways were returned (Figure 13).
Figure 13 Signaling pathway of XHP-TNBC PPI network (Red diamond stands for signaling pathway; Purple circle stands for XHP-TNBC; Blue circle stands for TNBC genes; Pink circle stands for XHP targets).
These signaling pathways are ranked according to the degree of enrichment (negative correlation with P value) and count from large to small. The top 10 is: Pathways in cancer (P= 5.47 * 10−49; Count = 209), Neuroactive ligand–receptor interaction (P = 4.54 * 10−38; Count = 153), PI3K-Akt signaling pathway (P = 5.01 * 10−34; Count = 170), Proteoglycans in cancer (P = 6.38 * 10−30; Count = 114), Focal adhesion (P = 1.22 * 10−25; Count = 110), HIF-1 signaling pathway (P = 1.10 * 10−24; Count = 67), FoxO signaling pathway (P = 8.15 * 10−23; Count = 80), TNF signaling pathway (P = 1.43 * 10−21; Count = 68), Neurotrophin signaling pathway (P = 1.60 * 10−21; Count = 73), Rap1 signaling pathway (P = 3.91 * 10−18; Count = 99) (Figure 12B). The details of each signaling pathway was described in Table S4.
Through the above network analysis, we found that Myrrha, Bovis Calculus, Olibanum and Moschus in XHP has anti-TNBC effect. Among them, the main compounds of anti-breast cancer in Olibanum are boswellic acids, and their targets are basically the same. Among them, CYP19A1, ESR1, NR3C1, ESR2, PTPN1, MAPT, PTPN2, SHBG and so on are the main common targets. The main anti-breast cancer chemical compounds of Olibanum are pentacyclic triterpenoids, such as beta-Boswellic acid, Acetyl-11-keto-β-boswellic acid, 11-keto-β-boswellic acid, Acetyl-α-boswellic acid, O-acetyl-α-boswellic acid and so on (see Figure S4).
Meanwhile, the main component of Myrrha are Quercetin, beta-Sitosterol, Guggulsterone, Limonene, Myrcenol. the main component of Bovis Calculus are Chenodeoxycholic acid, Deoxycholic Acid, Oleanolic acid, Cholesterol, Methyl deoxycholate. The main component of Moschus are Androst-4-ene-3,17-dione, Allantoin, Muscol, Estragole, Cholesteryl ferulate. The targets mainly regulated by them also include CYP19A1, ESR1, NR3C1, ESR2, PTPN1, MAPT, PTPN2, SHBG and so on.
At present, the research direction of TNBC therapeutic drugs is mainly three aspects: new targeted therapy, immunotherapy and new endocrine therapy (Lehmann et al., 2016; Garrido-Castro et al., 2019). The first aspect is the new targeted therapy, which is developed to reduce toxicity, reduce the risk of disease progression, and improve patient prognosis. A variety of targeted drugs have entered the clinical trial phase. The targeted treatment of TNBC is divided into four major categories (Lee et al., 2019; Nakhjavani et al., 2019). (1) Targeted therapy for DNA repair: The main research direction is poly ADP-ribose polymerase (PARP), a key enzyme involved in DNA repair, which recognizes breakpoints in DNA single strands and repairs (Lord and Ashworth, 2017; Caulfield et al., 2019). The BRCA gene (BRCA gene mutation detected in 11.2% of patients with TNBC) is sensitive to PARP inhibitors (Evans et al., 2017; Densham and Morris, 2019). Among them, olrapani combined with chemotherapy drugs achieved a high objective response rate in clinical trials. However, due to its mechanism of inhibiting DNA repair, patients are more likely to develop a second primary tumor (Parker, 2011). (2) Targeted therapy for tyrosine kinase inhibition (TKIs): TKIs are a class of drugs that inhibit the activity of tyrosine kinases (Smith et al., 2004), which blocks the downstream signaling pathway by inhibiting the phosphorylation of protein tyrosine residues, thereby inhibiting tumor growth and metastasis (Arora and Scholar, 2005). Both epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor (VEGFR) have receptors for tyrosine kinase activity, which can control the proliferation, invasion and metastasis of tumor cells by regulating various signaling pathways (Ayeni et al., 2015; Jayson et al., 2016). However, existing studies have shown that EGFR inhibitors alone are difficult to achieve anti-tumor effects (Ueda et al., 2017). Our transcriptomics and chemical informatics studies have shown that Olibanum can target EGFR and VEGFR. (3) PI3K/AKT/mTOR pathway inhibitor: the PI3K/AKT/mTOR pathway plays an important role in tumor cell proliferation, angiogenesis and metastasis, which is activated in breast cancer. Inhibitors of this pathway can inhibit the growth of breast cancer tumors and cause apoptosis of cancer cells, which are expected to become new drugs for TNBC targeted therapy (Costa et al., 2018; Kontzoglou and Garmpis, 2019). Our research shows that XHP can target the PI3K/AKT/mTOR pathway. (4) Silk-threonine protein kinase-related targeted drugs: PIMI kinase is a serine protein kinase and clinical studies have shown that PIMI kinase inhibitors can be used as new targeted therapies for patients with TNBC (Horiuchi et al., 2016). Studies have confirmed that the expression level of PIMI in TNBC tissues are higher than that in normal breast tissue and hormone receptor-positive breast cancer tissues (Zhao et al., 2017). In addition, PIMI plays an important role in the growth and proliferation of TNBC cells expressing MYC (Brasó-Maristany et al., 2016). Our chemical informatics studies show that XHP can regulate serine–threonine protein kinase-mediated biological processes. In the future, XHP and Olibanum may be used as sensitizers, and combined with monoclonal antibodies or inhibitors can significantly increase the objective response rate or survival rate of patients with TNBC.
The second aspect is immunotherapy. Tumor immunotherapy controls and kills tumor cells by stimulating or modulating the immune system to enhance the anti-tumor immunity of the tumor microenvironment (Lee Ventola, 2017). In the future, the most promising immunotherapeutic drug for the treatment of TNBC is the immunological checkpoint inhibitor (ICPI) (Li et al., 2018b). In addition, TNBC has a unique immune microenvironment. TNBC overexpresses VEGF, other molecules that promote tumor cell growth and migration, and has more tumor infiltrating lymphocytes (TIL) and tumor-associated macrophages (TAMs) (Saraiva et al., 2017; Yu and Di, 2017; Romero-Cordoba et al., 2019). The immune microenvironment of tumor cells is a place where the immune escape of tumor cells and the immune surveillance of the human immune system compete against each other. Therefore, the tumor immune microenvironment is closely related to the occurrence, development and prognosis of TNBC. The results of biological processes and signaling pathway enrichment analysis show that XHP can regulate the immune process in TNBC. Li et al. (2018a) and Su et al. (2018) have confirmed that XHP can inhibit tumor cell growth by improving the immunosuppressive state of the tumor microenvironment, promoting apoptosis of Treg cells in tumor microenvironment, and reversing immune escape.
The third aspect is endocrine therapy. Endocrine therapy refers to inhibition of tumor cell growth by regulating hormone levels and effects in the body. However, since the ER and PR in patients with TNBC are negative, traditional endocrine therapy is ineffective for them (Márquez-Garbán et al., 2019). Current research showed that androgen receptor (AR) inhibitors are potential therapeutic agents for TNBC, and AR signaling cascades can be inhibited by multiple pathways (Kono et al., 2017). Epidemiological surveys showed that 10%-35% of patients with TNBC express AR, while AR overexpression usually indicates poor prognosis (Barton et al., 2015). Therefore, blocking the expression of AR becomes a viable endocrine therapy for TNBC. For example, bicalutamide, the anti-AR inhibitor, achieved clinical benefit in 19% of patients in phase III clinical trials (Abe et al., 2014; Venema et al., 2019). In addition, TNBC’s endocrine therapy also includes targets such as gonadotropin-releasing hormone and growth hormone-releasing hormone, but the related studies are few.
Transcriptome Analysis of Olibanum-Treated MD-MB 231 Cells
Differentially Expressed Gene of Olibanum-Treated MD-MB 231 Cells
The transcriptome data comes from GSE102891. When performing the analysis, select “Boswellia Serrata Extract 128 ug/ml” as the experimental group, and select “control” as the control group to obtain gene expression data. Gene with a P value of <0.05 and Log2FC>1 or <−1 is considered to be a differentially expressed gene (DEG) (Figure S5). PCA plot showed that the results of Olibanum group, β-boswellic acid group and control group were significantly different (Figure S6A). The different clustering between Olibanum group and β-boswellic acid group were shown in Figure S6B. A total of 227 genes were identified as DEGs and were used to construct DEGs PPI network of Olibanum-treated MD-MB 231 cells (Figure 14A). The details of each DEG was shown in Table S5.
Figure 14 DEGs PPI networks. (A) DEGs PPI network of Olibanum-treated MD-MB 231 cells; (B) DEGs PPI network of β-boswellic acid -treated MD-MB 231 cells.
In Figure 14, the size of each node is related to its Degree; the bigger nodes have the larger value of Degree. The width of line is associated with its Edge Betweenness; the wider lines have the larger value of Edge Betweenness. In this network, the top 21 targets are: HSPA5 (33 edges), DDIT3 (28 edges), DNAJC3 (25 edges), ATF3 (24 edges), XBP1 (21 edges), HSPA1A (21 edges), PPP1R15A (19 edges), DNAJB9 (18 edges), TRIB3 (18 edges), HYOU1 (18 edges), PDIA3 (18 edges), ASNS (17 edges), RPS27A (16 edges), HMOX1 (16 edges), HERPUD1 (16 edges), DNAJB1 (16 edges), PDIA4 (16 edges), TXNRD1 (16 edges), EGR1 (15 edges), DUSP1 (15 edges), ERN1 (15 edges). These targets are considered to be the core targets of this network.
In the central node of the network, HSPA5, DDIT3, DNAJC3, ATF3, XBP1, PPP1R15A, DNAJB9, TRIB3, PDIA3, ASNS, HERPUD1, DNAJB1, PDIA4, TXNRD1, DUSP1, ERN1 are genes related to endoplasmic reticulum stress, which is essential for regulating the control of protein quality by endoplasmic reticulum and maintaining the balance of redox state (Xu and Park, 2018; Amen et al., 2019; Gao et al., 2019; Kim E. K. et al., 2019). Current research shows that Endoplasmic Reticulum Associated Unfolded Protein Response (UPR) and endoplasmic reticulum stress can affect the migration and invasion characteristics of breast cancer cells. The mechanisms include extracellular matrix (ECM) remodeling, cell adhesion modification, chemoattraction, epithelial-mesenchymal transition (EMT), regulation of signaling pathways associated with cell migration, and cytoskeletal remodeling, which in turn promotes breast cancer cell migration and invasion (Cook et al., 2016; Han and Wan, 2018; Sisinni et al., 2019). Hence, in the future, targeting UPR and breast cancer stress may be potential targeted therapeutic strategies for the treatment of breast cancer.
HSPA5, as a gene related to endoplasmic reticulum stress, plays an important role in the biological processes of breast cancer; inhibition of HSPA5 can inhibit TNBC cell migration and invasion (Chen et al., 2015). The endoplasmic reticulum stress-related gene DDIT3 also participates in the development of TNBC by regulating autophagy and apoptosis of TNBC cells (Singha et al., 2013). In vitro and in vivo showed that LYN-1604 exerts an anti-TNBC effect by targeting ULK1-regulated autophagic death; its induced autophagic death is closely related to these key genes such as ATF3, RAD21 and caspase3 (Ouyang et al., 2017; Zhang et al., 2017). The heat shock protein HYOU1 (also known as Orp150) plays an important role in hypoxia/ischemia and angiogenesis, which is overexpressed in certain invasive breast cancers, and its overexpression appears to be associated with poor prognostic indicators (Li et al., 2019b). In metastatic breast and ovarian cancer, HSPA1 has a different lysine methylation, and unmethylated HSPA1 shows potential as a prognostic marker in potentially highly serous carcinomas (Jakobsson et al., 2015).
Enrichment Analysis of DEGs PPI Network of Olibanum-Treated MD-MB 231 Cells
After the enrichment analysis, several biological processes and signaling pathways are obtained. The results of enrichment in the DAVID database indicate that these DEGs are primarily involved in biological processes associated with endoplasmic reticulum stress response, such as cellular response to topologically affected protein, MAPK signaling pathway, FoxO signaling pathway and so on. The main biological processes and signaling pathways were shown in Figures 15A, B. The details of each biological processes and signaling pathway were shown in Table S6.
Figure 15 Bubble chart. (A) biological processes of Olibanum-treated MD-MB 231 cells; (B) signaling pathways of Olibanum-treated MD-MB 231 cells; (C) biological processes of β-boswellic acid -treated MD-MB 231 cells; (D) signaling pathways of β-boswellic acid -treated MD-MB 231 cells.
Metascape is a portal that provides genetic annotation and analysis resources to help biologists understand one or more gene lists (http://metascape.org/gp/index.html#/main/step1) (Zhou et al., 2019). Metascape’s enrichment analysis of DEGs adds new biological processes and signaling pathways; compared to DAVID, the Reactom pathway is added, and the background annotation genes are more complete (Figure S7). Interestingly, the results of pathway enrichments indicate that Ferroptosis has the highest enrichment, and current studies show that induction of Ferroptosis has become an important strategy for the treatment of tumors (Galluzzi et al., 2018); especially in breast cancer, ferroptosis induces tumor cell death by activating the -ATF4 pathway of GCN2-eIF2α in TNBC cells (Chen et al., 2017). This suggests that the Olibanum extract may contain components that potentially induce Ferroptosis. The details of the biological processes, signaling pathways and Reactome pathways were shown in Table S7.
Transcriptome Analysis of β-Boswellic Acid-Treated MD-MB 231 Cells
Differentially Expressed Gene of β-Boswellic Acid-Treated MD-MB 231 Cells
The transcriptome data comes from GSE102891. When performing the analysis, select “3-O-Acetyl-β-boswellic acid 46 ug/ml” as the experimental group, and select “control” as the control group to obtain gene expression data. Gene with a P value of <0.05 and Log2FC >1 or <−1 is considered to be a differentially expressed gene (DEG) (Figure S8). A total of 950 genes were identified as DEGs and were used to construct DEGs PPI network of β-boswellic acid -treated MD-MB 231 cells (Figure 14). In Figure 14B, the size of each node is related to its Degree; the bigger nodes have the larger value of Degree. The width of line is associated with its Edge Betweenness; the wider lines have the larger value of Edge Betweenness. In this network, the top 20 targets are: FOS (51 edges), CXCL8 (49 edges), CD44 (46 edges), HSPA5 (39 edges), SIRT1 (38 edges), PTGS2 (37 edges), HDAC3 (36 edges), TIMP1 (33 edges), H2AFX (30 edges), TGFB1 (29 edges), UBE2C (28 edges), EGR1 (28 edges), ATF3 (27 edges), SERPINE1 (27 edges), HIST2H2BE (27 edges), CEBPB (27 edges), STUB1 (26 edges), TUFM (26 edges), CLPP (25 edges), HSPA1A (25 edges). Compared with DEGs PPI network of Olibanum-treated MD-MB 231 cell, their common targets are: HSPA5, ATF3, HSPA1A, EGR1. The details of each biological processes and signaling pathway were shown in Table S8.
FOS, HSPA5, SIRT1, PTGS2, HDAC3, TIMP1, TGFB1, ATF3, CEBPB, CLPP are involved in endoplasmic reticulum-related biological processes such as endoplasmic reticulum stress response and endoplasmic reticulum unfolded protein response (Wilkinson, 2019; Limia et al., 2019; Omidkhoda et al., 2019), which mediate the development of breast cancer. Recent studies have shown that the cell surface adhesion receptor CD44 is a key positive regulator of PD-L1 expression in TNBC and non-small cell lung cancer (NSCLC) (Kong et al., 2019). The Notch-mediated tumor-interstitial-inflammatory network promotes TNBC invasiveness and CXCL8 expression (Liubomirski et al., 2019). High expression of UBE2C is a potential factor for poor prognosis of TNBC. Furthermore, loss of BRCA1 function results in increased expression of UBE2C and chemoresistance to doxorubicin in breast cancer cells (Qin et al., 2017). The up-regulation of Serpine2 promotes breast cancer cell metastasis and reduces patient survival (Jin et al., 2017). Meanwhile, recent studies have shown that CHIP/STUB1 ubiquitin ligase is a negative chaperone molecule of HSP90/HSC70, which is reduced or lost in breast cancer. The absence of CHIP reshapes the cellular transcriptome and releases key cancer-promoting factors, such as the matrix degrading enzymes of the cathepsin family (Luan et al., 2018). In addition, downregulation of TUFM can induce epithelial-mesenchymal transition, and analogs of resveratrol HS-1793 can down-regulate its expression and increase anticancer activity against MCF-7 cells (Jeong et al., 2012).
Enrichment Analysis of DEGs PPI Network of β-Boswellic Acid-Treated MD-MB 231 Cells
After the enrichment analysis, several biological processes and signaling pathways are obtained. The results of enrichment in the DAVID database indicate that these DEGs are primarily involved in biological processes associated with endoplasmic reticulum stress response, such as PERK-mediated unfolded protein response, negative regulation of apoptotic process, and so on. The signaling pathways are protein processing in endoplasmic reticulum, TNF signaling pathway, Antigen processing and presentation, Estrogen signaling pathway, ECM-receptor interaction. The main biological processes and signaling pathways were shown in Figures 15C, D.
The biological processes in the Metascape database are similar to those in the DAVID database, in a different order (Figure S9). The details of the biological processes, signaling pathways and Reactome pathways were shown in Table S9.
In summary, through further analysis of transcriptomics data, it is found that the TNBC-related biological processes regulated by Olibanum extracts and β-boswellic acid mainly include endoplasmic reticulum stress response, oxidative stress, angiogenesis, inflammatory response, cell migration and adhesion, hypoxia induction, and autophagy. Of particular importance, the endoplasmic reticulum stress response is thought to be the primary upstream mechanism by which Olibanum extracts and β-boswellic acid play a role in TNBC cells. Hence, endoplasmic reticulum stress is a potential target in cancer therapy due to its important role in cancer development.
XHP may exert anti-TNBC effects through regulating biological processes, signaling pathways, targets found in this study. Meanwhile, the ability of Olibanum extracts and β-boswellic acid to induce endoplasmic reticulum stress and subsequently activate tumor cell death programs confirms that they are a promising class of anticancer agents. In addition, the Olibanum extracts may be used as an inducer of TNBC’s Ferroptosis in the future.
Data Availability Statement
All datasets for this study are included in the article/Supplementary Material.
LZ and KY dominated the concept and carried out a comprehensive design. LL, XX, AG, and TB were participants in the concept and design. KY, LZ, AG, TB, and LL are responsible for data analysis and interpretation in the chemical informatics section. TX and XX, and LL are responsible for data analysis and interpretation in experiments. KY, LZ, AG, and TB drafted the paper. LL and XX supervised the study. All authors participated in the analysis and interpretation of data and approved the final paper. KY, LZ, AG and TB should be considered joint first author. LL is the first corresponding author because she supervised the study.
This research is supported by Scientific Research Project of Hunan Provincial Department of Education (No. 19B434) and Double First-Class University Project of Hunan Province (Xiangjiaotong 469).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.00825/full#supplementary-material
Abe, E., Hirakawa, H., Suzuki, T., Nemoto, N., Miyashita, M., Tamaki, K., et al. (2014). Androgenic pathways in the progression of triple-negative breast carcinoma: a comparison between aggressive and non-aggressive subtypes. Breast Cancer Res. Treat. 145 (2), 281–293. doi: 10.1007/s10549-014-2942-6
Amen, O. M., Sarker, S. D., Ghildyal, R., Arya, A. (2019). Endoplasmic Reticulum Stress Activates Unfolded Protein Response Signaling and Mediates Inflammation, Obesity, and Cardiac Dysfunction: Therapeutic and Molecular Approach. Front. Pharmacol. 10, 977. doi: 10.3389/fphar.2019.00977
Ano, R., Kimura, Y., Shima, M., Matsuno, R., Ueno, T., Akamatsu, M. (2004). Relationships between structure and high-throughput screening permeability of peptide derivatives and related compounds with artificial membranes: Application to prediction of Caco-2 cell permeability. Bioorg. Med. Chem. 12, 257–264. doi: 10.1016/j.bmc.2003.10.002
Baraya, Y. S., Wong, K. K., Yaacob, N. S. (2017). The Immunomodulatory Potential of Selected Bioactive Plant-Based Compounds in Breast Cancer: A Review. Anticancer Agents Med. Chem. 17 (6), 770–783. doi: 10.2174/1871520616666160817111242
Barton, V. N., D’Amato, N. C., Gordon, M. A., Christenson, J. L., Elias, A., Richer, J. K. (2015). Androgen Receptor Biology in Triple Negative Breast Cancer: a Case for Classification as AR+ or Quadruple Negative Disease. Horm. Cancer. 6 (5-6), 206–213. doi: 10.1007/s12672-015-0232-3
Brasó-Maristany, F., Filosto, S., Catchpole, S., Marlow, R., Quist, J., Francesch-Domenech, E. (2016). PIM1 kinase regulates cell death, tumor growth and chemotherapy response revealing a novel target in triple-negative breast cancer. Nat. Med. 22 (11), 1303–1313. doi: 10.1038/nm.4198
Caulfield, S. E., Davis, C. C., Byers, K. F. (2019). Olaparib: A Novel Therapy for Metastatic Breast Cancer in Patients With a BRCA1/2 Mutation. J. Adv. Pract. Oncol. 10 (2), 167–174. doi: 10.6004/jadpro.2019.10.2.6
Cejalvo, J. M., Pascual, T., Fernández-Martínez, A., Brasó-Maristany, F., Gomis, R. R., Perou, C. M., et al. (2018). Clinical implications of the non-luminal intrinsic subtypes in hormone receptor-positive breast cancer. Cancer Treat Rev. 67, 63–70. doi: 10.1016/j.ctrv.2018.04.015
Chan, J. J., Tan, T. J. Y., Dent, R. A. (2019). Novel therapeutic avenues in triple-negative breast cancer: PI3K/AKT inhibition, androgen receptor blockade, and beyond. Ther. Adv. Med. Oncol. 11, 1758835919880429. doi: 10.1177/1758835919880429
Chen, Y., Zhou, C., Ge, Z., Liu, Y., Liu, Y., Feng, W., et al. (2013). Composition and potential anticancer activities of essential oils obtained from myrrh and frankincense. Oncol. Lett. 6 (4), 1140–1146. doi: 10.3892/ol.2013.1520
Chen, F. P., Chang, C. M., Hwang, S. J., Chen, Y. C., Chen, F. J. (2014). Chinese herbal prescriptions for osteoarthritis inTaiwan: analysis of National Health Insurance dataset. BMC Complement Altern. Med. 14, 91. doi: 10.1186/1472-6882-14-91
Chen, H. A., Chang, Y. W., Tseng, C. F., Chiu, C. F., Hong, C. C., Wang, W., et al. (2015). E1A-mediated inhibition of HSPA5 suppresses cell migration and invasion in triple-negative breast cancer. Ann. Surg. Oncol. 22 (3), 889–898. doi: 10.1245/s10434-014-4061-3
Chen, M. S., Wang, S. F., Hsu, C. Y., Yin, P. H., Yeh, T. S., Lee, H. C., et al. (2017). CHAC1 degradation of glutathione enhances cystine-starvation-induced necroptosis and ferroptosis in human triple negative breast cancer cells via the GCN2-eIF2α-ATF4 pathway. Oncotarget 8 (70), 114588–114602. doi: 10.18632/oncotarget.23055
Cheng, S., Qin, X., Jin, S., Qian, N. (2016). The comparative study on the antitumous effects of Xihuang pills , Frankincense , Myrrh , and combination of Frankincense and Myrrh on diverse human malignant tumor cells in vitro. Pharmacol. Clinics Chin. Mater. Med. 32 (3), 102–104. doi: 10.13412/j.cnki.zyyl.2016.03.027
Cook, K. L., Soto-Pantoja, D. R., Clarke, P. A., Cruz, M. I., Zwart, A., Wärri, A., et al. (2016). Endoplasmic Reticulum Stress Protein GRP78 Modulates Lipid Metabolism to Control Drug Sensitivity and Antitumor Immunity in Breast Cancer. Cancer Res. 76 (19), 5657–5670. doi: 10.1158/0008-5472.CAN-15-2616
Costa, R. L. B., Han, H. S., Gradishar, W. J. (2018). Targeting the PI3K/AKT/mTOR pathway in triple-negative breast cancer: a review. Breast Cancer Res. Treat. 169 (3), 397–406. doi: 10.1007/s10549-018-4697-y
Desantis, C. E., Ma, J., Goding Sauer, A., Newman, L. A., Jemal, A. (2017). Breast cancer statistics 2017, racial disparity in mortality by state. CA Cancer J. Clin. 67 (6), 439–448. doi: 10.3322/caac.21412
Evans, K. W., Yuca, E., Akcakanat, A. (2017). A Population of Heterogeneous Breast Cancer Patient-Derived Xenografts Demonstrate Broad Activity of PARP Inhibitor in BRCA1/2 Wild-Type Tumors. Clin. Cancer Res. 23 (21), 6468–6477. doi: 10.1158/1078-0432.CCR-17-0615
Galluzzi, L., Vitale, I., Aaronson, S. A., Abrams, J. M., Adam, D., Agostinis, P., et al. (2018). Molecular mechanisms of cell death: recommendations of the Nomenclature Committee on Cell Death 2018. Cell Death Differ. 25 (3), 486–541. doi: 10.1038/s41418-017-0012-4
Gao, Y., Kim, S., Lee, Y. I., Lee, J. (2019). Cellular Stress-Modulating Drugs Can Potentially Be Identified by in Silico Screening with Connectivity Map (CMap). Int. J. Mol. Sci. 20 (22), E5601. doi: 10.3390/ijms20225601
Garrido-Castro, A. C., Lin, N. U., Polyak, K. (2019). Insights into Molecular Classifications of Triple-Negative Breast Cancer: Improving Patient Selection for Treatment. Cancer Discovery 9 (2), 176–198. doi: 10.1158/2159-8290.CD-18-1177
Gfeller, D., Grosdidier, A., Wirth, M., Daina, A., Michielin, O., Zoete, V. (2014). SwissTargetPrediction: a web server for target prediction of bioactive small molecules. Nucleic Acids Res. 42, W32–W38. doi: 10.1093/nar/gku293
Guo, Q., Xu, X., He, S., Yuan, Y., Chen, S., Hua, B. (2018). Xi huang pills enhance the tumor treatment efficacy when combined with chemotherapy: A meta-analysis and systematic review. J. Cancer Res. Ther. 14 (Supplement), S1012–S1018. doi: 10.4103/0973-1482.192795
Hao, J., Jin, Z., Zhu, H., Liu, X., Mao, Y., Yang, X., et al. (2018). Antiestrogenic Activity of the Xi-Huang Formula for Breast Cancer by Targeting the Estrogen Receptor α. Cell Physiol. Biochem. 47 (6), 2199–2215. doi: 10.1159/000491533
He, L. J., Li, J. S., Chen, X., Zhang, H. T., Zhu, X. G., Zhang, X. D., et al. (2018). Effect of serum containing Xihuang pill on proliferation of human breast cancer cell line MDA-MB-435 and MCF-7 Cells. Zhongguo Zhong Yao Za Zhi. 43 (13), 2784–2788. doi: 10.19540/j.cnki.cjcmm.20180408.001. Chinese.
Horiuchi, D., Camarda, R., Zhou, A. Y., Yau, C., Momcilovic, O., Balakrishnan, S. (2016). Pim1 kinase inhibition as a targeted therapy against triple-negative breast tumors with elevated myc expression. Nat. Med. 22 (11), 1321–1329. doi: 10.1038/nm.4213
Jakobsson, M. E., Moen, A., Davidson, B., Falnes, PØ. (2015). Hsp70 (HSPA1) Lysine Methylation Status as a Potential Prognostic Factor in Metastatic High-Grade Serous Carcinoma. PloS One 10 (10), e0140168. doi: 10.1371/journal.pone.0140168. 2015.
Jayson, G. C., Kerbel, R., Ellis, L. M., Harris, A. L. (2016). Antiangiogenic therapy in oncology: current status and future directions. Lancet 388 (10043), 518–529. doi: 10.1016/S0140-6736(15)01088-0
Jeong, S. H., Song, I. S., Kim, H. K., Lee, S. R., Song, S., Suh, H., et al. (2012). An analogue of resveratrol HS-1793 exhibits anticancer activity against MCF-7 cells via inhibition of mitochondrial biogenesis gene expression. Mol. Cells 34 (4), 357–365. doi: 10.1007/s10059-012-0081-7
Jiang, Y. Z., Ma, D., Suo, C., Shi, J., Xue, M., Hu, X., et al. (2019). Genomic andTranscriptomic Landscape of Triple-Negative Breast Cancers: Subtypes and Treatment Strategies. Cancer Cell. 35 (3), 428–440.e5. doi: 10.1016/j.ccell.2019.02.001
Jin, T., Suk Kim, H., Ki Choi, S., Hye Hwang, E., Woo, J., Suk Ryu, H., et al. (2017). microRNA-200c/141 upregulates SerpinB2 to promote breast cancer cell metastasis and reduce patient survival. Oncotarget 8 (20), 32769–32782. doi: 10.18632/oncotarget.15680
Keiser, M. J., Roth, B. L., Armbruster, B. N., Ernsberger, P., Irwin, J. J., Shoichet, B. K. (2007). Relating protein pharmacology by ligand chemistry. Nat. Biotech. 25 (2), 197–206. doi: 10.1038/nbt1284
Kong, T., Ahn, R., Yang, K., Zhu, X., Fu, Z., Morin, G., et al. (2019). CD44 promotes PD-L1 expression and its tumor-intrinsic function in breast and lung cancers. Cancer Res. 80, 444–457. doi: 10.1158/0008-5472. canres.1108.2019.
Kono, M., Fujii, T., Lim, B., Karuturi, M. S., Tripathy, D., Ueno, N. T. (2017). Androgen Receptor Function and Androgen Receptor-Targeted Therapies in Breast Cancer: A Review. JAMA Oncol. 3 (9), 1266–1273. doi: 10.1001/jamaoncol.2016.4975
Kuhn, M., von Mering, C., Campillos, M., Jensen, L. J., Bork, P. (2008). STITCH: interaction networks of chemicals and proteins. Nucleic Acids Res. 36 (Database issue), D684–D688. doi: 10.1093/nar/gkm795
Kuhn, M., Szklarczyk, D., Franceschini, A., Campillos, M., von Mering, C., Jensen, L. J., et al. (2012). STITCH 3: zooming in on protein-chemical interactions Nucleic Acids Res. 40, Database issue, D876–D880. doi: 10.1093/nar/gkr1011
Kulkarni, A., Stroup, A. M., Paddock, L. E., Hill, S. M., Plascak, J. J., Llanos, A. A. M. (2019). Breast Cancer Incidence and Mortality by Molecular Subtype: Statewide Age and Racial/Ethnic Disparities in New Jersey. Cancer Health Disparities 3, e1–e17. doi: 10.9777/chd.2019.1012
Kunnumakkara, A. B., Banik, K., Bordoloi, D., Harsha, C., Sailo, B. L., Padmavathi, G., et al. (2018). Googling the Guggul (Commiphora and Boswellia) for Prevention of Chronic Diseases. Front. Pharmacol. 9, 686. doi: 10.3389/fphar.2018.00686
Lee, K. L., Kuo, Y. C., Ho, Y. S., Huang, Y. H. (2019). Triple-Negative Breast Cancer: Current Understanding and Future Therapeutic Breakthrough Targeting Cancer Stemness. Cancers (Basel) 11 (9), E1334. doi: 10.3390/cancers11091334
Lehmann, B. D., Jovanović, B., Chen, X., Estrada, M. V., Johnson, K. N., Shyr, Y., et al. (2016). Refinement of Triple-Negative Breast Cancer Molecular Subtypes: Implications for Neoadjuvant Chemotherapy Selection. PloS One 11 (6), e0157368. doi: 10.1371/journal.pone.0157368
Li, X. Y., Su, L., Jiang, Y. M., Gao, W. B., Xu, C. W., Zeng, C. Q., et al. (2018a). The Antitumor Effect of Xihuang Pill on Treg Cells Decreased in Tumor Microenvironment of 4T1 Breast Tumor-Bearing Mice by PI3K/AKT~AP-1 Signaling Pathway. Evid Based Complement Alternat Med. 2018, 6714829. doi: 10.1155/2018/6714829
Li, H., Chen, Y., Wang, X., Tang, L., Guan, X. (2019a). T1-2N0M0 Triple-Negative Breast Cancer Treated With Breast-Conserving Therapy Has Better Survival Compared to Mastectomy: A SEER Population-Based Retrospective Analysis. Clin. Breast Cancer S1526-8209 105 (19), 30316–30317. doi: 10.1016/j.clbc.2019.05.011
Li, X., Zhang, N. X., Ye, H. Y., Song, P. P., Chang, W., Chen, L., et al. (2019b). HYOU1 promotes cell growth and metastasis via activating PI3K/AKT signaling in epithelial ovarian cancer and predicts poor prognosis. Eur. Rev. Med. Pharmacol. Sci. 23 (10), 4126–4135. doi: 10.26355/eurrev_201901_17914
Limia, C. M., Sauzay, C., Urra, H., Hetz, C., Chevet, E., Avril, T. (2019). Emerging Roles of the Endoplasmic Reticulum Associated Unfolded Protein Response in Cancer Cell Migration and Invasion. Cancers (Basel) 11 (5), E631. doi: 10.3390/cancers11050631
Liubomirski, Y., Lerrer, S., Meshel, T., Morein, D., Rubinstein-Achiasaf, L., Sprinzak, D. (2019). Notch-Mediated Tumor-Stroma-Inflammation Networks Promote Invasive Properties and CXCL8Expression in Triple-Negative Breast Cancer. Front. Immunol. 10, 804. doi: 10.3389/fimmu.2019.00804
Luan, H., Mohapatra, B., Bielecki, T. A., Mushtaq, I., Mirza, S., Jennings, T. A., et al. (2018). Loss of the Nuclear Pool of Ubiquitin Ligase CHIP/STUB1 in Breast Cancer Unleashes the MZF1-Cathepsin Pro-oncogenic Program. Cancer Res. 78 (10), 2524–2535. doi: 10.1158/0008-5472.CAN-16-2140.
Márquez-Garbán, D. C., Deng, G., Comin-Anduix, B., Garcia, A. J., Xing, Y., Chen, H. W., et al. (2019). Antiestrogens in combination with immune checkpoint inhibitors in breast cancerimmunotherapy. J. Steroid Biochem. Mol. Biol. 193, 105415. doi: 10.1016/j.jsbmb.2019.105415
Mangone, L., Mancuso, P., Tagliabue, G., Filiberti, R. A., Carrozzi, G., Iacovacci, S., et al. (2019). Neoadjuvant therapy for breast cancer. Tumori 105, 300891619869505. doi: 10.1177/0300891619869505
Mao, D., Feng, L., Huang, S., Zhang, S., Peng, W., Zhang, S. (2019). Meta-Analysis of Xihuang Pill Efficacy When Combined with Chemotherapy for Treatment of Breast Cancer. Evid Based Complement Alternat Med. 2019, 3502460. doi: 10.1155/2019/3502460
Mazzio, E. A., Lewis, C. A., Soliman, K. F. A. (2017). Transcriptomic Profiling of MDA-MB-231 Cells Exposed to Boswellia Serrata and 3-O-Acetyl-B-Boswellic Acid; ER/UPR Mediated Programmed Cell Death. Cancer Genomics Proteomics. 14 (6), 409–425. doi: 10.21873/cgp.20051
Metodiewa, D., Kochman, A., Karolczak, S. (1997). Evidence for antiradical and antioxidant properties of four biologically active N,N-Diethylaminoethyl ethers of flavaone oximes: A comparison with natural polyphenolic flavonoid rutin action. IUBMB Life 41, 1067–1075. doi: 10.1080/15216549700202141
Missiuro, P. V., Liu, K., Zou, L., Ross, B. C., Zhao, G., Liu, J. S., et al. (2009). Information flow analysis of interactome networks. PloS Comput. Biol. 5, e1000350. doi: 10.1371/journal.pcbi.1000350
Mukhopadhyay, U. K., Oturkar, C. C., Adams, C., Wickramasekera, N., Bansal, S., Medisetty, R., et al. (2019). TP53 Status as a Determinant of Pro- vs Anti-Tumorigenic Effects of Estrogen Receptor-Beta in Breast Cancer. J. Natl. Cancer Inst. 111 (11), 1202–1215. doi: 10.1093/jnci/djz051
Nakhjavani, M., Hardingham, J. E., Palethorpe, H. M., Price, T. J., Townsend, A. R. (2019). Druggable Molecular Targets for the Treatment of Triple Negative Breast Cancer. J. Breast Cancer. 22 (3), 341–361. doi: 10.4048/jbc.2019.22.e39
Omidkhoda, N., Wallace Hayes, A., Reiter, R. J., Karimi, G. (2019). The role of MicroRNAs on endoplasmic reticulum stress in myocardial ischemia and cardiac hypertrophy. Pharmacol. Res. 150, 104516. doi: 10.1016/j.phrs.2019.104516
Ouyang, L., Zhang, L., Fu, L., Liu, B. (2017). A small-molecule activator induces ULK1-modulating autophagy-associated cell death in triple negative breast cancer. Autophagy 13 (4), 777–778. doi: 10.1080/15548627.2017.1283470
Prat, A., Pineda, E., Adamo, B., Galván, P., Fernández, A., Gaba, L., et al. (2015). Clinical implications of the intrinsic molecular subtypes of breast cancer. Breast 24 Suppl 2, S26–S35. doi: 10.1016/j.breast.2015.07.008
Qin, T., Huang, G., Chi, L., Sui, S., Song, C., Li, N., et al. (2017). Exceptionally high UBE2C expression is a unique phenomenon in basal-like type breast cancer and is regulated by BRCA1. BioMed. Pharmacother. 95, 649–655. doi: 10.1016/j.biopha.2017.08.095
Romero-Cordoba, S., Meneghini, E., Sant, M., Iorio, M. V., Sfondrini, L., Paolini, B., et al. (2019). Decoding Immune Heterogeneity of Triple Negative Breast Cancer and Its Association with Systemic Inflammation. Cancers (Basel) 11 (7), E911. doi: 10.3390/cancers11070911
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 13. doi: 10.1186/1758-2946-6-13
Saraiva, D. P., Guadalupe Cabral, M., Jacinto, A., Braga, S. (2017). How many diseases is triple negative breast cancer: the protagonism of the immune microenvironment. ESMO Open 2 (4), e000208. doi: 10.1136/esmoopen-2017-000208
Singha, P. K., Pandeswara, S., Venkatachalam, M. A., Saikumar, P. (2013). Manumycin A inhibits triple-negative breast cancer growth through LC3-mediated cytoplasmic vacuolation death. Cell Death Dis. 4, e457. doi: 10.1038/cddis.2012.192
Sisinni, L., Pietrafesa, M., Lepore, S., Maddalena, F., Condelli, V., Esposito, F., et al. (2019). Endoplasmic Reticulum Stress and Unfolded Protein Response in Breast Cancer: The Balance between Apoptosis and Autophagy and Its Role in Drug Resistance. Int. J. Mol. Sci. 20 (4), E857. doi: 10.3390/ijms20040857
Smith, J. K., Mamoon, N. M., Duhé, R. J. (2004). Emerging roles of targeted small molecule protein-tyrosine kinase inhibitors incancer therapy. Oncol. Res. 14 (4-5), 175–225. doi: 10.3727/000000003772462298
Su, L., Jiang, Y., Xu, Y., Li, X., Gao, W., Xu, C., et al. (2018). Xihuang pill promotes apoptosis of Treg cells in the tumor microenvironment in 4T1 mouse breast cancer by upregulating MEKK1/SEK1/JNK1/AP-1 pathway. BioMed. Pharmacother. 102, 1111–1119. doi: 10.1016/j.biopha.2018.03.063
Szarc Vel Szic, K., Declerck, K., Crans, R. A. J., Diddens, J., Scherf, D. B., Gerhäuser, C., et al. (2017). Epigenetic silencing of triple negative breast cancer hallmarks by Withaferin A. Oncotarget 8 (25), 40434–40453. doi: 10.18632/oncotarget.17107
Traina, T. A., Miller, K., Yardley, D. A., Eakle, J., Schwartzberg, L. S., O’Shaughnessy, J., et al. (2018). Enzalutamide for the Treatment of Androgen Receptor-Expressing Triple-Negative Breast Cancer. J. Clin. Oncol. 36 (9), 884–890. doi: 10.1200/JCO.2016.71.3495
Ueda, S., Saeki, T., Osaki, A., Yamane, T., Kuji, I. (2017). bevacizumab induces acute hypoxia and cancer progression in patients with refractory breast cancer: multimodal functional imaging and multiplex cytokine analysis. Clin. Cancer Res. 23, 5769–5778. doi: 10.1158/1078-0432.CCR-17-0874
Venema, C. M., Bense, R. D., Steenbruggen, T. G., Nienhuis, H. H., Qiu, S. Q., van Kruchten, M., et al. (2019). Consideration of breast cancer subtype in targeting the androgen receptor. Pharmacol. Ther. 200, 135–147. doi: 10.1016/j.pharmthera.2019.05.005
Warner, E. T., Tamimi, R. M., Hughes, M. E., Ottesen, R. A., Wong, Y. N., Edge, S. B., et al. (2015). Racial and Ethnic Differences in Breast Cancer Survival: Mediating Effect of Tumor Characteristics and Sociodemographic and Treatment Factors. J. Clin. Oncol. 33 (20), 2254–2261. doi: 10.1200/JCO.2014.57.1349
Xu, X., Zhang, W., Huang, C., Li, Y., Yu, H., Wang, Y., et al. (2012). A novel chemometric method for the prediction of human oral bioavailability. Int. J. Mol. Sci. 13, 6964–6982. doi: 10.3390/ijms13066964
Yang, K., Zeng, L., Ge, J. (2018). Exploring the Pharmacological Mechanism of Danzhi Xiaoyao Powder on ER-Positive Breast Cancer by a Network Pharmacology Approach. Evid Based Complement Alternat Med. 2018, 5059743. doi: 10.1155/2018/5059743
Yang, K., Zeng, L., Ge, A., Chen, Z., Bao, T., Long, Z., et al. (2019). Investigating the regulation mechanism of baicalin on triple negative breast cancer’s biological network by a systematic biological strategy. BioMed. Pharmacother. 118, 109253. doi: 10.1016/j.biopha.2019.109253
Zeng, L., Yang, K. (2017). Exploring the pharmacological mechanism of Yanghe Decoction on HER2-positive breast cancer by a network pharmacology approach. J. Ethnopharmacol. 199, 68–85. doi: 10.1016/j.jep.2017.01.045
Zhang, L., Fu, L., Zhang, S., Zhang, J., Zhao, Y., Zheng, Y., et al. (2017). Discovery of a small molecule targeting ULK1-modulated cell death of triple negative breast cancer in vitro and in vivo. Chem. Sci. 8 (4), 2687–2701. doi: 10.1039/c6sc05368h
Zheng, W., Han, S., Jiang, S., Pang, L., Li, X., Liu, X., et al. (2016). Multiple effects of Xihuang pill aqueous extract on the Hs578T triple-negative breast cancer cell line. BioMed. Rep. 5 (5), 559–566. doi: 10.3892/br.2016.769
Keywords: Xihuang Pill, triple negative breast cancer, transcriptomics, chemical informatics, MDA-MB-231 cell, high-throughput omics, bioinformatics, Chinese medicine
Citation: Yang K, Zeng L, Ge A, Bao T, Xu T, Xie X and Liu L (2020) Exploring the Regulation Mechanism of Xihuang Pill, Olibanum and β-Boswellic Acid on the Biomolecular Network of Triple-Negative Breast Cancer Based on Transcriptomics and Chemical Informatics Methodology. Front. Pharmacol. 11:825. doi: 10.3389/fphar.2020.00825
Received: 06 December 2019; Accepted: 19 May 2020;
Published: 11 June 2020.
Edited by:Jiarui Wu, Beijing University of Chinese Medicine, China
Reviewed by:Haci Ahmet Deveci, University of Gaziantep, Turkey
Jia-bo Wang, Fifth Medical Center of the PLA General Hospital, China
Copyright © 2020 Yang, Zeng, Ge, Bao, Xu, Xie and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.