[Deciphering cellular processes responding to lethality of 17β-estradiol by quantitative phosphoproteomics]

17β-Estradiol (E2), an important endocrine hormone in the mammalian body, participates in the regulation of the physiological functions of the reproductive system, mammary glands, bone, and cardiovascular system, among others. Paradoxically, despite the physiological actions of endogenous E2 (0.2-1.0 nmol/L), numerous clinical and experimental studies have demonstrated that high-dose E2 treatment can cause tumor regression and exert pro-apoptotic actions in multiple cell types; however, the underlying mechanism remains undescribed. In particular, little information of the cellular processes responding to the lethality of E2 is available. In the present study, we attempted to characterize the cellular processes responding to high-dose (μmol/L) E2 treatment using quantitative phosphoproteomics to obtain a better understanding of the regulatory mechanism of E2-induced cell death. First, the cell phenotype induced by high-dose E2 was determined by performing Cell Counting Kit-8 assay (CCK8), cell cytotoxicity analysis by trypan blue staining, and microscopic imaging on HeLa cells treated with 1-10 μmol/L E2 or dimethyl sulfoxide (DMSO) for 1-3 d. E2 inhibited cell proliferation and induced cell death in a dose- and time-dependent manner. Compared with the DMSO-treated HeLa cells, the cells treated with 5 μmol/L E2 for 2 d demonstrated >74% growth inhibition and approximately 50% cell death. Thus, these cells were used for quantitative phosphoproteomic analysis. Next, a solid-phase extraction (SPE)-based immobilized titanium ion affinity chromatography (Ti⁴⁺-IMAC) phosphopeptide-enrichment method coupled with data-independent acquisition (DIA)-based quantitative proteomics was employed for the in-depth screening of high-dose E2-regulated phosphorylation sites to investigate the intracellular processes responding to high-dose E2 treatment. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) identified over 10000 phosphorylation sites regulated by E2 and DMSO in HeLa cells. In comparison with the DMSO-treated cells, the cells treated with 5 μmol/L E2 showed 537 upregulated phosphorylation sites and 387 downregulated phosphorylation sites, with a threshold of p<0.01 and |log₂(fold change)|≥1. A total of 924 phosphorylation sites on 599 proteins were significantly regulated by high-dose E2, and these sites were subjected to enrichment analysis. In addition, 453 differently regulated phosphorylation sites on 325 proteins were identified only in the E2- or DMSO-treated cell samples. These phosphorylation sites may be phosphorylated or dephosphorylated in response to high-dose E2 stimulation and were subjected to parallel enrichment analyses. Taken together, 1218 phosphorylation sites on 741 proteins were significantly regulated by high-dose E2 treatment. The functional phosphoproteins in these two groups were then analyzed using Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) to determine the biological processes in which they participate and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Consistent with the cell-phenotype data, cell cycle-related proteins were highly enriched in the two groups of E2-regulated phosphoproteins (p<0.05), indicating that high-dose E2 treatment can regulate cell proliferation. In addition, E2-regulated phosphoproteins were highly enriched in the cellular processes of ribosome biogenesis, nucleocytoplasmic transport, and messenger ribonucleic acid (mRNA) processing/splicing (p<0.05), indicating that the activation of these processes may contribute to high-dose E2-induced cell death. These results further confirm that high-dose E2 treatment inhibits protein translation and induces cell death. Furthermore, the significant upregulation of multiple phosphorylation sites associated with epidermal growth factor receptor (EGFR) and mitogen-activated protein kinases (MAPKs) MAPK1, MAPK4, and MAPK14 by high-dose E2 indicates that the EGFR and MAPK signaling pathways are likely involved in the regulation of E2-induced cell death. These phosphorylation sites likely play vital roles in E2-induced cell death in HeLa cells. Overall, our phosphoproteomic data could be a valuable resource for uncovering the regulatory mechanisms of E2 in the micromolar range.