Prediction of the Que Target Genes
The chemical structure of Que was obtained from the PubChem database as shown in Figure 2A. Based on its structure, PharmMapper database, TCMSP database, and Swiss Target Prediction database were used to predict the potential target genes of Que. The targets obtained from these databases were combined and duplicates were removed to obtain 482 potential target genes.
Prediction of Precancerous lesions of breast cancer-Related Que Target Genes by System Pharmacology Approach
By GeneCards database (https://www.genecards.org/), with "Precancerous lesions," as the keyword search genes associated with breast cancer lesion before, get 1602 disease targets information. Intersecting Que 482 potential targets genes with 1602 disease targets resulted in 314 genes, as shown in Figure 2B, and these target genes were not only Precancerous lesions related genes but also the drug targets. The PPI network graph of target contains 157 nodes and 589 edges, in which the nodes represent targets and the edges represent the interaction between targets. Import the PPI data into Cytoscape 3.6.1 software for visualization. As shown in Figure 2C, the top 20 targets for the degree value of target are MAPK1, SRC, HRAS, AKT1, HSP90AA1, MAPK8, RHOA, MAPK14, ESR1, EGFR, IGF1, RXRA, JAK2, CASP3, AR, PTK2, IL-2, STAT1 and MMP9.
Functional Enrichment Analysis for the Que Target Genes
In order to explore the relationship between these 314 potential target genes, the GO and KEGG enrichment analysis were done by the STRING database. A variety of GO enrichment terms were enriched, including 1980 biological processes, 100 cellular components, and 223 molecular functions. We found that the biological processes such as cellular process and metabolic process, cellular components like binding and catalytic activity, and molecular function like cell and intracellular (Figure 2D), were enriched, which may be involved in the biological activity in the Que treatment process. In addition, 176 KEGG pathways were enriched (Figure 2E), among the major Que related pathways screened by KEGG, the pathways were mainly associated with cancer pathway and PI3K-Akt signal pathway, however, JAT/STAT ranked fourth among the pathways screened for immunity to Que (Figure 2F). It is suggested that Que may improve the progression of breast cancer from precancerous lesions to BC or improve the prognosis of breast cancer patients through JAK/STAT1 signaling pathway.
Cell proliferation experiments in MCF-10AT, MCF-7 and MDA-MB-231 cells
The effects of Que the proliferation of MCF-10AT, MCF-7 and MDA-MB-231 cells were determined by MTT assay and the results were given in Figure 3C. The half maximal inhibitory concentration (IC50) values of Que on MCF-10AT, MCF-7 and MDA-MB-231 cells were 52.39µM, 53.76µM and 64.23μM, respectively.
The apoptosis of MCF-10AT, MCF-7 and MDA-MB-231 cells was induced by Que at different concentrations at different time periods
To evaluate the apoptotic properties of Que, induced cell death, we performed the Annexin-V binding assay. We detected the apoptosis of MCF-10A, MCF-10AT, MCF-7, MDA-MB-231 cells treated with Que at 5μM, 20μM, 80μM, and 120μM at 24 h and 48 h, respectively (Figure 3A and 3B). First, Que was treated with mammary gland cells for 24h. We compared MCF-10AT with MCF-10A, MCF-7 and MDA-MB-231 cell lines at Que concentration of 5μM. The total apoptotic cell population % was determined as 2.85±0.21% and 3.27±0.14%, 2.84±0.21%, 2.85± 0.21% (P=0.006, P=0.109, P=0.004, respectively), there was significant statistical significance. MCF-10AT was compared with MCF-10A, MCF-7, MDA-MB-231 cell lines at 20μM concentration of Que. The total apoptotic cell population % was determined as 7.60±0.21% and 3.31±0.12%, 12.76±0.37%, 20.17± 2.90% (P=0.004, P=0.002, P=0.004, respectively), there was significant statistical significance. MCF-10AT was compared with MCF-10A, MCF-7, MDA-MB-231 at 80μM concentration of Que. The total apoptotic cell population % was determined as 35.63±0.60% and 1.64±0.05%, 25.34±5.06%, 29.80± 0.98% (P=0.002, P=0.004, P=0.004, respectively), there was significant statistical significance. MCF-10AT was compared with MCF-10A, MCF-7 and MDA-MB-231 cell lines at 120μM concentration of Que. The total apoptotic cell population % was determined as 40.25±0.71% and 3.40±0.07%, 29.38±2.50%, 29.40± 2.51% (P=0.002, P=0.002, P=0.002, respectively), there was significant statistical significance. Next, We found that MCF-10AT was compared with MCF-10A, MCF-7, MDA-MB-231 cell lines at 5μM concentration of Que, the total apoptotic cell population % was determined as 1.75±0.33% and 1.34±0.15%, 9.63±1.70%, 9.58± 0.32% (P=0.240, P=0.004, P=0.004, respectively). MCF-10AT was compared with MCF-10A, MCF-7, MDA-MB-231 cell lines at 20μM concentration of Que, the total apoptotic cell population % was determined as 8.44±1.54% and 0.84±0.10%, 11.12±1.38%, 14.23±1.62% (P=0.004, P=0.016, P=0.004, respectively), there was significant statistical significance. MCF-10AT was compared with MCF-10A, MCF-7 and MDA-MB-231 cell lines at 80μM concentration of Que, the total apoptotic cell population % was determined as 72.36±1.77% and 3.52±0.14%, 61.95±1.22%, 74.73±0.77% (P=0.004, P=0.004, P=0.009, respectively), there was significant statistical significance. MCF-10AT was compared with MCF-10A, MCF-7, MDA-MB-231 cell lines at 120μM concentration of Que, the total apoptotic cell population % was determined as 79.19±1.76% and 5.52±0.30%, 61.95±1.22%, 74.37±0.80% (P=0.004, P=0.004, P=0.002, respectively), there was significant statistical significance (Figure 3D). These results indicate that Que can induce apoptosis on MCF-10AT, MCF-7, MDA-MB-231 and other pre-breast cancer cells and BC cells at different times and at different concentrations, showing a time and concentration dependence. However, Que has a stronger apoptosis effect on MCF-10AT cells with pre-breast cancer lesions. However, MCF-10A cells had no apoptosis effect. Apoptosis trends at different concentrations at 24h and 48h were shown in the figure (Figure 3E).
Que can promote γδ T cell subsets and their amplification
In vitro induction of peripheral blood mononuclear cells from healthy volunteers in complete medium containing TCR γδ monoclonal antibody and cytokine IL-2. After adding 0μM, 2.5μM, 5μM, and 10μM Que on day 3, the cell morphology was observed under microscope and the proportion of γδ T cell subsets was determined by flow cytometry after 10-12 days. We found that the number of γδ T cells increased with the increase of the number of days in vitro. Flow cytometry showed that the γδ T cells could be effectively amplified, and the purity of γδ T cells could reach more than 60%, and the γδ T cells could be expanded to 90% at the concentration of 2.5μM to 5μM. Vδ2 T cell subsets were dominant at 2.5μM Que, Vδ1 T cell subsets were dominant at 5μM Que, and Vδ2 T cell killing subsets were dominant at 10μM Que (Figure 4A and B). These results indicate that γδ T cells have better killing and immunomodulatory effects.
Cytotoxicity of healthy human γδ T cells against MCF-10A, MCF-10AT, MCF-7, MDA-MB-231 cell lines
After amplification in vitro, the effector cell: target cell (E/T) was 10:1, the killing rates of γδ T cells against MCF-10A, MCF-10AT, MCF-7, MDA-MB-231 were 61.44±4.70, 55.52±3.10, 53.94±2.74, 53.28±1.73 (P=0.114, P=0.486, P=0.343, respectively) (Figure 4C). There was no significant difference between the groups. These results indicate that γδ T cells have a certain killing effect on both precancer and BC cells.
Que and γδ T cells have synergistic cytotoxic effects on mammary gland cell
In order to investigate the killing effect of Que and γδ T cells on mammary gland cell, we used 5μM Que to detect E/T (1:1, 5:1, 10:1) respectively to investigate the specific killing effect. We found that the cell killing rates of MCF-10A in E/T (1:1, 5:1, 10:1) were 24.12±4.34, 51.93±6.47, 64.94±3.61. The cytotoxicity rates of MCF-10AT in E/T (1:1, 5:1, 10:1) were 19.38±5.30, 33.45±5.49, 64.96±5.45,and MCF-10A > MCF-10AT (P=0.222, P=0.008, P=0.917, respectively). The cell killing rates of MCF-7 at E/T (1:1, 5:1, 10:1) were 13.23±2.68, 24.39±3.13, 55.59±5.98, compared with MCF-10AT, the killing rate of MCF-7 versus MCF-10AT was statistically significant at 5:1 (P=0.095, P=0.032, P=0.056, respectively). The cell killing rates of MDA-MB-231 at E/T: 1:1, 5:1, 10:1 were 12.77±3.64, 22.7±1.39, 59.04±5.67, respectively. Thus, compared with MDA-MB-231, MCF-10AT was also statistically significant at 5:1 (P=0.056, P=0.016, P=0.222, respectively). In addition, with the increase of effector cell proportion, Que concentration was still 5μM, and MCF-10A > MCF-10AT > MCF-7 > MDA-MB-231, 1:1 < 5:1 < 10:1 (Figure 4D). These results indicated that Que combined with γδ T cells had a certain specific killing effect on both precancerous breast cancer cells and BC cells, and the strongest killing effect on precancerous breast cancer cells and BC cells was found when the Que concentration was 5μM and E/T (10:1).
Effect of Que on IFNγ-R, Phospho-JAK2 (p-JAK2), Phospho-STAT1 (p-STAT1) and PD-L1 in MCF‐10AT and MCF-7 cell lines protein expression
Western blotting of IFNγ-R, p-JAK2, p-STAT1 and PD-L1protein were performed. The concentrations of Que were 0μM, 5μM, 20μM, 80μM and 120μM, and the cells were co-cultured with MCF-10AT and MCF-7 cells for 48 h, respectively. Our results showed that when MCF-10AT and MCF-7 cells were treated with Que concentration of 80μM and 120μM, respectively, a substantial increase in protein levels of IFNγ-R, p-JAK2 and p-STAT1 were significantly increased (P<0.0001), while the concomitant decrease protein levels of PD-L1 (P<0.0001, P=0.0005) (Figure 5A and B).