Daurisoline inhibits ESCC cell proliferation in vitro
By drug screening, 40 μM daurisoline was found to have obvious cytotoxicity to KYSE450 cells (Fig. 1a). Daurisoline is a bis-benzylisoquinoline alkaloid (Fig. 1b). The IC50 of KYSE150 and KYSE450 cells were 20.733 and 13.560 μM at 24 h, and 11.0149 and 9.936 μM at 48 h, respectively. Under daurisoline treatment for 24 h and 48 h, IC50 of immortalized esophageal epithelial cells (SHEE cells) was 41.102 μM and 11.406 μM, respectively. With the increase in daurisoline concentration, the survival rate of cells decreased gradually (Fig. 1c). In the cell proliferation assay, we used the following drug concentrations: 0, 2.5, 5, 10, and 20 μM. Daurisoline inhibited the proliferation of KYSE150 and KYSE450 cells in a dose-dependent manner. The proliferation inhibition rates of KYSE150, KYSE450, and SHEE cells were 44.51%, 68.83%, and 25.50%, respectively, after treatment with 20 μM daurisoline for 96 h (Fig. 1d). These results indicated that daurisoline was less toxic to SHEE cells compared with ESCC cells. To verify this, we further performed the anchorage-independent cell growth assay using the same drug concentration. Daurisoline inhibited the clonal size and number of clones formed for KYSE150 and KYSE450 cells. At 20 μM, the clonal formation rates of KYSE150 and KYSE450 cells were 30.25% and 42.36%, respectively, compared with that of the control cells (Fig. 1e). These results indicated that daurisoline inhibited ESCC cell proliferation in vitro.
Daurisoline significantly downregulates the phosphorylation of ERK1/2 in ESCC cells
We used quantitative phosphoproteomics to investigate the molecular mechanism underlying daurisoline inhibiting ESCC cell proliferation. We treated KYSE150 cells with 20 μM daurisoline for 24 h. (Fig. 2a). The first-order mass error of most spectrograms was within 10 ppm, which conformed to the high precision characteristic of orbital well MS (Fig. S1a). Most peptides were 7-20 amino acids in length. The distribution of the peptide lengths identified by MS conformed to the quality control requirements (Fig. S1b). The number of secondary spectra obtained was 383,273. After screening, the number of available effective spectrograms was 94,520. We identified 12,965 phosphorylation sites on 3130 proteins, of which 7415 sites on 2500 proteins contained quantitative data. To ensure a high credibility of the results, p > 0.75 was used to filter the identification data. Consequently, we identified 8549 sites on 2861 proteins, of which 6601 sites on 2408 proteins included quantitative data (Fig. S1c). The protein quantification group was normalized to remove the influence of protein expression on the modified signal. The screening criteria for different sites were as follows: 1.5 times the change threshold and two-sample two-tailed t-test (p < 0.05). Based on these data and standards, we found that the phosphorylation of 176 sites was enhanced and that of 340 sites was reduced in the daurisoline treated cells (Fig. S1d).
A total of 340 down regulated phosphorylation sites were enriched in the Kyoto Encyclopedia of genes and genomes (KEGG) pathway. KEGG pathway was selected according to Fisher's exact test (p < 0.05). Ultimately, nine KEGG pathway targets were screened: viral carcinogenesis, spliceosome, RNA transport, adipocytokine signaling pathway, herpes simplex virus infection, synaptic vesicle cycle, Salmonella infection, chronic myeloid leukemia and Th17 cell differentiation (Fig. 2b). The heat map shows the levels of the phosphorylation sites in the viral carcinogenesis pathway, and most of the phosphorylation sites are downregulated (Fig. 2c). In six KEGG pathways, multiple proteins appeared repeatedly in the different pathways. Therefore, Wayne diagrams are made based on the six KEGG pathways to find the key proteins in these pathways (Fig. 2d). We found that NFκB1 was enriched in six pathways, ERK2 in four pathways, and STAT3 in three pathways. The original mass spectra of p-ERK2 T185/Y187 are shown in Supplement Figure (Fig. S1e). In the volcanic map, according to the -log10 ( P-value), phosphorylation of the four phosphorylation sites is ranked from top to bottom as p-ERK2 Y187 (3.0177), p-ERK2 T185 (2.4976), p-STAT3 S727 (1.7433), and p-NFκB1 S907 (1.5686) (Fig. 2e). In addition, the protein expression ratios of the above four phosphorylation sites are 0.504, 0.458, 0.566, and 0.503, respectively (Fig. 2f). In conclusion, we believe that p-ERK2 T185/Y187 are important phosphorylation sites in the above nine KEGG pathways. In addition, although ERK1 was not enriched in the KEGG pathway, the phosphorylation of ERK1 was also downregulated in the phosphoproteomics data. The phosphorylation level of ERK1/2 in KYSE150 and KYSE450 cells was significantly downregulated after daurisoline (20 μM) treatment for 24 h (Fig. 2g). These results indicated that daurisoline inhibited ERK1/2 phosphorylation in ESCC cells.
Daurisoline directly binds with MEK1 /2
We used IGPS1.0 software to predict upstream kinases for specific phosphorylation sites in the above nine KEGG pathways. The prediction results showed that the upstream kinase proteins of p-ERK2 T185 and p-ERK2 Y187 were MEK1-7 (MAP2K1-7), and ALK and LTK, respectively (Table. S1). In order to verify the accuracy of the upstream kinase prediction, we performed the CNBr-activated Sepharose 4B pull down assay (Fig. S2). The protein binding to the drug was detected by Mass spectrum. The results of mass spectrometry showed that the protein binding to daurisoline were MEK1, MEK2, MEK3, and MEK7 with Unused ProtScores of 38.69, 17.28, 15.76, and 4.83, respectively (Table S2). Previous studies on MAPK signaling pathway showed that the sole substrate of MEK1/2 protein kinase is ERK1/2, the substrate of MEK3 is p38MAPK, and the substrate of MEK7 is JNK [26]. Therefore, based on the Unused ProtScores, daurisoline was likely to affect the biological function of MEK1/2 and downregulate ERK1/2 phosphorylation.
Autodock 4.0 software was used for molecular docking in order to explore the specific binding sites of daurisoline with MEK1/2 proteins. Daurisoline could bind to MEK1 (binding energy: -11.06 kcal/mol). The binding sites were Asn78 and Lys97 (Fig. 3a). Daurisoline could also bind to MEK2 (binding energy: -8.50 kcal/mol). The binding sites were Asp194 and Asp212 (Fig. 3b). Pull down assays displayed that daurisoline could directly bind to MEK1/2 (Fig. 3c, d). In addition, daurisoline also binds directly to MEK1 and MEK2 in KYSE150 and KYSE450 cells (Fig. 3e, f). To verify the accuracy of the molecular simulation results, we substituted the above-mentioned amino acids of MEK1 and MEK2 proteins with alanine (Ala, A). The results showed that after Asn78 substitution on MEK1 protein, the binding efficiency of the protein with daurisoline was significantly reduced (Fig. 3g). Furthermore, the binding efficiency of MEK2 protein with daurisoline was significantly reduced after Asp194 substitution (Fig. 3h).
Daurisoline inhibits the MEK1/2-ERK1/2 signaling pathway in ESCC
We performed an in vitro kinase assay. The phosphorylation efficiency of active MEK1 and MEK2 against inactive ERK2 was in turn reduced at 2.5, 5, 10, and 20 μM daurisoline (Fig. 4a, b). These results indicated that daurisoline could inhibit ERK2 phosphorylation by inhibiting MEK1/2 activity. Moreover, Coomassie blue staining showed the location and expression of MEK1/2 and ERK2 proteins.
We used gene set enrichment analysis (GSEA) to enrich the KEGG pathway associated with altered expression of all proteins identified in the phosphoproteomics analysis. We found that ERK1/2 was enriched and significantly downregulated in the Erbb signaling pathway (Fig. 4c) and prostate cancer pathway (Fig. 4d). In addition, MEK2 (MAP2K2) was also enriched in the aforementioned two pathways. Western blotting showed that daurisoline inhibited ERK1/2 phosphorylation in KYSE150 and KYSE450 cells. However, daurisoline had no significant effect on MEK1/2 phosphorylation (Fig. 4c, d). The above results indicate that daurisoline inhibits the MEK1/2-ERK1/2 signaling pathway in ESCC cells.
MEK1/2 were significantly overexpressed in ESCC
MEK1/2 have been reported to be highly expressed in ESCC [27]. The Cancer Genome Atlas (TCGA) database revealed that MEK1 (Fig. 5a) and MEK2 (Fig. 5b) were both significantly overexpressed in a variety of cancers, including ESCA. In addition, MEK1 (Fig. 5c) and MEK2 (Fig. 5d) are also overexpressed in EAC and ESCC. We used ESCC data from the TCGA database to conduct a multi-gene correlation analysis. We observed a significant correlation between MEK1 and ERK2 in ESCC (Fig. 5f, g).
Daurisoline inhibits ESCC tumor growth in vivo
To investigate the anticancer effect of daurisoline in vivo, we established the PDX model of ESCC (Fig. 6a). There was no significant difference in the body weight between mice treated with daurisoline and control, indicating that daurisoline had no toxic or side effects on mice (Fig. 6b). The final tumor volume was measured at the end of the treatment (Fig. 6c). The tumor growth rate of daurisoline treated group was significantly lower than that of control group (Fig. 6d, g). The tumor weight also decreased with the increase of the dose of daurisoline (Fig. 6e). We calculated the tumor weight inhibition rate of daurisoline treated group. The growth inhibition rates in the 20 and 40 mg/kg daurisoline treatment groups were 52.77% and 84.74%, respectively (Fig. 6f). Immunohistochemical staining showed that the Ki67 and ERK1/2 phosphorylation decreased after daurisoline treatment. The MEK1/2 phosphorylation did not change significantly (Fig. 6h). In conclusion, daurisoline inhibited ESCC tumor proliferation and ERK1/2 phosphorylation in tumor cells in vivo.