Abstract
Estrogen therapy was used to treat advanced breast cancer in postmenopausal women for decades until the introduction of tamoxifen. Resistance to long-term estrogen deprivation (LTED) with tamoxifen and aromatase inhibitors used as a treatment of breast cancer inevitably occurs, but unexpectedly low-dose estrogen can cause regression of breast cancer and increase disease-free survival in some patients. This therapeutic effect is attributed to estrogen-induced apoptosis in LTED breast cancer. Here, we describe modulation of the estrogen receptor (ER) liganded with antiestrogens (endoxifen and 4-hydroxytamoxifen) and an estrogenic triphenylethylene (TPE), ethoxytriphenylethylene (EtOXTPE), on estrogen-induced apoptosis in LTED breast cancer cells. Our results show that the angular TPE estrogen (EtOXTPE) is able to induce the ER-mediated apoptosis only at a later time compared with planar estradiol in these cells. Using real-time polymerase chain reaction, chromatin immunoprecipitation, western blotting, molecular modeling, and X-ray crystallography techniques, we report novel conformations of the ER complex with an angular estrogen EtOXTPE and endoxifen. We propose that alteration of the conformation of the ER complexes, with changes in coactivator binding, governs estrogen-induced apoptosis through the protein kinase regulated by RNA-like endoplasmic reticulum kinase sensor system to trigger an unfolded protein response.
Introduction
High-dose synthetic estrogen therapy was the first successful chemical treatment of any cancer (Haddow et al., 1944). Estrogen was used clinically to treat metastatic breast cancer (MBC) for 30 years prior to the successful introduction of the antiestrogen tamoxifen (Jordan, 2003). Tamoxifen is a structural derivative of the synthetic estrogen triphenylethylene (TPE). During the time before tamoxifen, clinical rules were defined for the successful application of high-dose estrogen treatment. Therapy was only effective in 30% of MBC and only if applied more than 5 years past menopause (Haddow, 1970). Mechanisms were unknown.
The advantage of tamoxifen compared with high-dose estrogen to treat MBC was not an increase in response rate, but a lower incidence of side effects (Cole et al., 1971; Ingle et al., 1981). Targeting tamoxifen to patients with estrogen receptor (ER) positive breast cancer, and the safety of the medicine, permitted implementation of translational research (Jordan, 2008, 2014b) to establish the value of tamoxifen as a long-term adjuvant therapy (Early Breast Cancer Trialists’ Collaborative Group, 1998). The proposition in the 1970s (Jordan et al., 1979; Jordan and Allen, 1980) that longer adjuvant therapy would be superior to shorter therapy also demanded an investigation of acquired resistance to long-term tamoxifen treatment in breast cancer. Nothing was known.
Acquired resistance to tamoxifen is unique. Initially, laboratory models in vivo demonstrated that tamoxifen actually stimulated tumor growth within 1 to 2 years (Gottardis and Jordan, 1988; Gottardis et al., 1989b). Nevertheless, low-dose estrogen also stimulated tumor growth. Mechanisms have subsequently been deciphered using breast cancer cell models in vitro (Fan et al., 2014a,b,c). A new steroidal pure antiestrogen was subsequently developed following proof of efficacy in vivo, to prevent tumor growth in tumors with acquired resistance (Gottardis et al., 1989a). Fulvestrant is now approved for the first line and second line treatment of MBC (Howell et al., 2002; Osborne et al., 2002; Moscetti et al., 2017).
It, therefore, came as some surprise to find that low-dose estrogen therapy triggered breast tumor regression following the development of acquired resistance to tamoxifen treatment of 5 years (Yao et al., 2000). This experimental model mimics the 5 years of adjuvant tamoxifen therapy that was the standard of care at the time. The rules described by Haddow (1970) applied to adjuvant tamoxifen therapy: 5 years of estrogen deprivation in ER-positive breast cancer is necessary to sensitize selected cell populations to undergo estrogen-induced apoptosis (Jordan, 2015). This created a general principle of ER positive breast cancer cell biology. Most importantly, low -dose estrogen salvage therapy is effective in producing a 30% clinical benefit rate in patients failing long-term adjuvant therapy with aromatase inhibitors (Ellis et al., 2009). Indeed, the science of estrogen-induced apoptosis has also been linked to the antitumor effects of estrogen therapy alone in the Women’s Health Initiative (Abderrahman and Jordan, 2016) and responsible for the “carry over effect” that maintains patients recurrence free after adjuvant therapy is terminated at 5 years (Jordan, 2014a).
The clinical significance of estrogen-induced apoptosis requires an understanding of molecular mechanisms to decipher appropriate applications for clinical care. We address the modulation of estrogen-induced apoptosis using estrogens of different shapes and related synthetic nonsteroidal antiestrogens.
Haddow et al. (1944) documented antitumor activity with both planar synthetic estrogens (diethylestibestrol) and angular estrogens (TPEs). In recent years, these synthetic estrogens were classified based on the resulting ER complex to activate or block an estrogen-responsive gene. Class I estrogens are planar estrogens and class II estrogens are angular estrogens (Jordan et al., 2001), based on their estrogenic or antiestrogenic activity at the transforming growth factor α gene. One explanation for the TPE estrogens having early antiestrogenic actions on estradiol-induced apoptosis (Maximov et al., 2011) was that the class II estrogen complex initially induces an antiestrogenic ER:class II estrogen complex (Obiorah et al., 2014), which then evolves to an estrogenic ER:class I complex to trigger apoptosis (Obiorah et al., 2014; Obiorah and Jordan, 2014).
Here, we address the hypothesis using X-ray crystallography of a novel type II estrogen known as ethoxytriphenylethylene (EtOXTPE) (Maximov et al., 2010, 2011) and compare and contrast the ER conformation with the TPE antiestrogen endoxifen, which is the major biologically active secondary metabolite of tamoxifen. The biology of these ER complexes to modulate estrogen-induced apoptosis now opens up new opportunities to examine the elasticity of unfolded protein response (UPR) to trigger or block apoptosis through the protein kinase regulated by RNA-like endoplasmic reticulum kinase (PERK) sensor in long-term estrogen deprivation (LTED) in breast cancer cell lines.
Materials and Methods
Cell Culture and Reagents.
The test compound was synthesized and the details of the synthesis have been reported previously (Maximov et al., 2010). The 17β-estradiol (E2) and 4-hydroxytamoxifen (4OHT) were acquired from Sigma-Aldrich (St. Louis, MO). Endoxifen (Z-isomer) was purchased from Santa Crus Biotechnology (Santa Cruz, CA). The ER positive and LTED breast cancer cells MCF7:5C were derived from MCF7 cells as reported previously (Jiang et al., 1992; Lewis et al., 2005b). MCF7:5C cells were maintained in phenol-red free RPMI 1610 media containing 10% charcoal dextran–treated fetal bovine serum, 6 ng/ml bovine insulin, L-glutamine, penicillin, and streptomycin, and were incubated at 37°C with 5% CO2. The cells were treated with indicated compounds for a specified time and then harvested.
Growth Assays.
For the assays, 10,000 MCF7:5C cells were used for the 7-day growth assay and 5000 MCF7:5C cells were used for the 14-day assays. The cells were plated in each well of 24-well plates and then treated the next day with specific concentrations of indicated compounds. The assays were carried out as described previously (Maximov et al., 2011, 2014) using a fluorescent DNA quantitation kit purchased from Bio-Rad (Hercules, CA) with sonication of samples after harvesting in isotonic buffer. All growth assays were performed in triplicate; the results represent the average of all replicates and error bars represent the S.D. in each treatment. One-way analysis of variance (ANOVA) was used with a follow-up Tukey’s test to determine the statistical significance of the treatments.
Immunoblotting.
The MCF-7:5C cells were seeded on 5-cm Petri dishes at a density of 2 million cells per plate and incubated overnight. The cells were treated for specified times with the indicated compounds. Protein isolation and immunoblotting were performed as previously described (Maximov et al., 2014). The primary antibodies used were anti-ERα clone G-20 (Santa Cruz Biotechnology), anti-eIF2α (9722) and anti-peIF2α (D9G8) (Santa Cruz Biotechnology), and goat anti-β-actin antibodies (Santa Cruz Biotechnology) diluted in 5% dry milk in Tris-buffered saline/Tween 20 blocking buffer at ratios recommended by the supplier at 4°C. All secondary antibodies were horseradish peroxidase linked (Santa Cruz Biotechnology). The signals were visualized by enhanced chemiluminescence. All immunoblots were performed in three replicates; data presented represent one the biologic replicates. Analysis was validated by densitometry using Image J (National Institutes of Health, Bethesda, MD), and the densitometry data are presented in Supplemental Tables 1–3.
Chromatin Immunoprecipitation (ChIP) Assay.
ChIP was performed as described previously (Sengupta et al., 2010; Obiorah et al., 2014) with minor modifications. The DNA fragments were purified using Qiaquick PCR purification kit (Qiagen, Germantown, MD). Then, 2 μl of eluted DNA was used for real-time polymerase chain reaction (RT-PCR) analysis. The primer sequences used are as follows: TFF1 promoter: 5′TGGGCTTCATGAGCTCCTTC3′ (forward); 5′TTCATAGTGAGAGATGGCCGG3′ (reverse). The data are expressed as percent input of starting chromatin material after subtracting the percent input pull down of the negative control (normal rabbit IgG). The assay was performed in triplicate; error bars represent the S.D. in each treatment. One-way ANOVA was used with a follow-up Tukey’s test to determine the statistical significance of the treatments.
Annexin V Staining.
MCF-7:5C cells were seeded at 300,000 cells per 10-cm Petri dishes and treated the next day with test compounds. Cells were treated with test compounds for 6 days, and for 3 days with 1 nM E2. Cells were harvested by aspirating media and washing cells with warm phosphate-buffered saline twice, and subsequently treated with accutase solution (Life Technologies, Grand Island, NY) for 4 minutes at 37°C. Cells were then harvested by pipetting after addition of phosphate-buffered saline and then transferred to centrifuge tubes and precipitated. Cells were put on ice afterward and stained using FITC Annexin V Apoptosis Detection Kit I (BD Pharmingen, San Diego, CA) according to the manufacturer’s instructions. The assay was performed in triplicate; data represent the average of the biologic replicates, and error bars represent the S.D. in each treatment. One-way ANOVA was used with a follow-up Tukey’s test to determine the statistical significance of the treatments.
RT-PCR.
MCF-7:5C cells, depending on the duration of treatment, were seeded at a density of 100,000–300,000 per well into six-well plates. Cells were treated the next day with test compounds for specified time points. RNA isolation, cDNA synthesis, and RT-PCR were performed as previously described (Obiorah et al., 2014). The primer sequences that were used for human TFF1 cDNA amplification are: 5′-CATCGACGTCCCTCCAGAAGA-3′ sense, and 5′-CTCTGGGACTAATCACCGTGCTG-3′ antisense; human GREB1 gene: 5′-CAAAGAATAACCTGTTGGCCCTGC-3′sense, 5′-GACATGCCTGCGCTCTCATACTTA-3′ antisense; human BCL2L11 gene: 5′-TCGGACTGAGAAACGCAAG-3′ sense, 5′-CTCGGTCACACTCAGAACTTAC-3′ antisense; human TP63 gene: 5′-TTCGGACAGTACAAAGAACGG-3′ sense, 5′-GCATTTCATAAGTCTCACGGC-3′ antisense; human HMOX1 gene: 5′- TCAGGCAGAGGGTGATAGAAG-3′ sense, 5′-TTGGTGTCATGGGTCAGC-3′ antisense; human TNFα gene: 5′-ACTTTGGAGTGATCGGCC-3′ sense, 5′-GCTTGAGGGTTTGCTACAAC-3′ antisense; and the reference gene RPLP0: 5′-GTGTTCGACAATGGCAGCAT-3′ sense, 5′-GACACCCTCCAGGAAGCGA-3′ antisense. All primers were obtained from Integrated DNA Technologies Inc. (Coralville, IA). All treatments were performed in triplicate; data represent average of the replicates, and error bars represent the S.D. in each treatment. One-way ANOVA was used with a follow-up Tukey’s test to determine the statistical significance of the treatments.
X-Ray Crystallography.
ERα ligand-binding domain (LBD) was incubated with a mixture of both cis- and trans-EtOXTPE isomers prior to crystallization. Separation of geometric isomers was not undertaken. Inclusion of the glucocorticoid receptor interacting protein 1 peptide was necessary to obtain diffraction quality crystals as was the use of Y537S mutation (Nettles et al., 2008). This mutation favors the agonist state of the receptor by forming a hydrogen bond with D351. Sitting drop was used for this crystallization. Clear rectangular crystals were observed after 1 week in 0.15 M KBr, 30% polyethylene glycol (PEG) methyl ether 2000, Tris pH 8.3. Paratone-N was used as the cryoprotectant. The structure was solved to 2.10 Å using molecular replacement (1GWR was used as the starting model) and one dimer was observed in the asymmetric unit (ASU). All EtOXTPE molecules were well ordered in the hormone-binding pocket. Only the trans-isomer of EtOXTPE was present. Coordinates have been deposited in the Protein Data Bank (PDB) under accession code 5T1Z. For estradiol complex, ERα LBD Y537S was incubated with 1 mM ligand and 2 mM glucocorticoid receptor interacting protein 1 peptide overnight at 4°C. Hanging drop vapor diffusion with a Hampton VDX plate (Hampton Research, Aliso Viejo, CA) was used for crystallization. Clear rectangular crystals were observed after 2 days in 20% PEG 3350, 100 mM MgCl2, and Tris pH 8.0. Paratone-N was used as the cryoprotectant. The structure was solved to 1.65 Å with one dimer in the ASU by molecular replacement with 1GWR as the input model. All estradiol ligands were resolved in the hormone-binding pocket. Coordinates were deposited in the PDB under accession code 6CBZ.
An ERα LBD construct with mutations C381S, C417S, C530S, and L536S was used to obtain co-crystal structures with endoxifen or 4OHT. Protein was expressed and purified as previously described (Fanning et al., 2016). Protein was incubated with 2 mM ligand at 4°C overnight prior to crystallization. Hanging drop at room temperature was used for these crystallizations. Clear hexagonal crystals were observed for the endoxifen co-crystals in Tris pH 8.0, 2 mM MgCl2, and 25% PEG 8000 after 5 days. For the 4OHT co-crystal structure, clear hexagonal crystals were observed after 1 week in Tris pH 6.5, 2 mM MgCl2, and 30% PEG 8000. All 4OHT ligands are resolved in the hormone-binding pocket. The endoxifen structure was solved with molecular replacement (using 5ACC as a starting model) to 1.65 Å with one dimer in the ASU. The 4OHT structure was solved by molecular replacement (with 5ACC as the starting model) to 1.80 Å with two dimers in the ASU. All endoxifen ligands were resolved in the hormone-binding pocket. Both structures were deposited in the PBD with accession codes 5W9C for the 4OHT structure and 5W9D for the endoxifen structure.
While structures of ERα LBD with E2 or 4OHT have been published, they have not been published using the same construct that we were using for EtOXTPE or endoxifen (with the 537S mutation for EtOXTPE or 381S, 417S, and 530S for endoxifen). Obtaining these structures enabled a more appropriate comparison between the ligands within the hormone-binding pocket. Furthermore, the previously published ER:4OHT structure (PDB:3ERT) possesses a crystal contact just after helix 11 that perturbs the loop connecting helices 11 and 12 and appears to alter how helix 12 sits in the activation function 2 cleft. Our structure does not have this crystal contact. Therefore, to properly compare endoxifen and 4OHT it was especially important to obtain this new 4OHT structure. The omit maps are shown in Supplemental Fig. 1.
Results
Effect of Nonplanar Estrogen EtOXTPE Alone on MCF-7:5C Cells.
To assess the activity of the test compounds on the viability of the MCF-7:5C cells (Fig. 1) we employed the proliferation assay as described in Materials and Methods. Following the increases in DNA levels at low concentrations (Fig. 2A), E2 started to induce apoptosis in cells at a concentration of 10−11 M, reducing the amount of viable cells in the wells by 60% after 7 days of treatment (Fig. 2A) (P < 0.05 vs. vehicle control). At a concentration of 10−10 M, E2 inhibited cell growth further by more than 90% (Fig. 2A) after 1 week of treatment. All other test compounds were used within a 10−12 to 10−6 M concentration range. Compound EtOXTPE produced partial agonist activity, inhibiting cell growth on average by 30% at the highest concentration of 10−6 M (Fig. 2A) after 1 week of treatment. Antiestrogen endoxifen produced no apoptotic activity with any statistically significant differences at any concentrations tested (Fig. 2A).
Chemical structures of compounds used in this study. EtOXPTE test compound is a 1:1 mixture or geometric isomers.
Cell proliferation assays in MCF-7:5C antihormone-resistant breast cancer cells. (A) Effects of test compounds alone after 7 days of treatment. Results show that the only compound able to completely inhibit the growth of the cells is E2, and the other test compounds possess no or only minor inhibitory effects. (B) Antiestrogenic effects of test compounds in combination with 1 nM E2 after 7 days of treatment. Results show that all compounds exhibit antiestrogenic effects after 7 days of treatment with the nonsteroidal antiestrogens 4OHT and endoxifen completely inhibiting the effect of E2 and angular estrogen EtOXTPE being able to inhibit E2 according to its intrinsic activity alone. (C) Effects of test compounds alone on cells after 14 days of treatment. Results show that besides E2, EtOXTPE can inhibit the cell growth after 14 days of treatment. All nonsteroid antiestrogenic compounds (4OHT and endoxifen) did not produce any inhibitory growth effect. (D) Reversal of EtOXTPE-induced inhibitory effect on MCF-7:5C cells by 4OHT after pretreatment of the cells with the test compound for various durations. All treatments were performed in triplicate; data represent the average of the replicates.
Effect of Nonplanar Estrogen EtOXTPE in Combination with 1 nM E2 on MCF7:5C Cells.
To assess the antiestrogenic properties of the test compounds to block E2-induced apoptosis, MCF-7:5C cells were treated with the compounds in combination with 1 nM E2 for 7 days. The results of the DNA quantification show that endoxifen completely inhibits 1 nM E2-induced apoptosis in the 10−7 to 10−6 M range with no statistical significance between the DNA values at these concentrations and vehicle control (P > 0.05) (Fig. 2B). Compound EtOXTPE has only partial antiestrogenic properties inhibiting 1 nM E2 according to its intrinsic activity alone at 10−6 M, and does not completely block E2 action (P < 0.05 compared with vehicle control) (Fig. 2B).
Effect of Nonplanar Estrogen EtOXTPE on MCF-7:5C Cells after an Extended Treatment.
To assess the activity of the test compounds of the cellular viability alone after an extended treatment we performed the same treatments as described previously for 14 days. As a result of the longer treatment, the DNA quantification assay showed that nonplanar estrogen is able to induce apoptosis in MCF-7:5C cells. Compound EtOXTPE is able to reduce the amount of cells by more than 90% starting at a 10−7 M concentration (Fig. 2C) (P < 0.05) with an IC50 value of approximately 2 × 10−9 M (Fig. 2C). Antiestrogen endoxifen is completely inactive alone and is not able to produce any reduction of cells at any of the concentration points (Fig. 2C).
Reversal of Nonplanar Estrogen EtOXTPE Effects by an Antiestrogen in MCF-7:5C Cells.
To test the possibility of reversal of the proapoptotic actions of the nonplanar estrogen EtOXTPE we treated MCF-7:5C cells with the compound for various durations (Fig. 2D), after which antiestrogen 4OHT was added at the concentration of 10−6 M. The treatments with E2, EtOXTPE, and 4OHT alone were used as controls. The cells were harvested for DNA fluorescent quantification assay as described in the Materials and Methods section after a total of 14 days of treatment. The results show that the test compound EtOXTPE is able to induce apoptosis in the cells at 10−6 M concentration after 14 days of treatment alone (Fig. 2D) by more than 90% compared with the vehicle control (P < 0.05), which is consistent with the dose-response curves with the same compound for equal duration of treatment (Fig. 2D). Addition of 4OHT at an equimolar concentration at early time points almost completely reversed the apoptotic actions of EtOXTPE after 1 and 2 days of treatment (Fig. 2D). However, after 4 days of treatment with EtOXTPE, 4OHT was able to reverse the apoptotic actions of the estrogen only by 50% of vehicle control (Fig. 2D). After 5 days of treatment with EtOXTPE, 4OHT was able to consistently reverse the apoptotic action of the compound by 30% on average throughout the remaining days of treatment in the experiment (Fig. 2D). This result is consistent with previously observed results with another TPE called bisphenol TPE (Obiorah and Jordan, 2014). Compared with planar E2, both EtOXTPE and biphenol TPE (Obiorah and Jordan, 2014) induce apoptosis consistently in MCF-7:5C cells later than E2 (Obiorah and Jordan, 2014).
The Apoptotic Effect of Nonplanar Estrogen EtOXTPE on MCF-7:5C Cells Is Delayed.
Cells undergo visual morphologic changes during treatment with EtOXTPE at a delayed rate compared with E2 (Fig. 3). The results of annexin V staining demonstrate that EtOXTPE is able to induce positive annexin V staining after 6 days of treatment compared with vehicle control (P < 0.05) (Fig. 4A), whereas endoxifen did not produce statistically significant changes in annexin V staining compared with vehicle control (Fig. 4A), or any morphologic changes (Fig. 3). Besides the positive annexin V staining results with the test compound, nonplanar estrogen is also able to induce activation of the proapaptotic genes, such as TNFα and BCL2L11 in MCF-7:5C cells after 120 hours of treatment compared with the vehicle control (P < 0.05) (Fig. 4, B and C). All these data indicate that nonplanar estrogen reduces the number of viable cells in the growth assay due to apoptosis and is delayed compared with E2. To assess the role of PERK signaling in EtOXTPE-induced apoptosis, we performed annexin V staining of the treated MCF-7:5C cells with the compound alone and in combination with PERK inhibitor GSK2606414 (Fig. 4D). We show that the inhibition of PERK with a selective inhibitor completely abrogates the EtOXTPE-induced apoptosis when compared with the vehicle control or to the PERK inhibitor alone treatment (Fig. 4D) (P > 0.05 in both cases), which is consistent with previously published data on E2-induced apoptosis in the same cell line (Fan et al., 2013).
Bright-field microscopy photographs of MCF-7:5C cell morphology after being treated with test compounds for indicated durations.
(A) Annexin V staining of MCF-7:5C cells after 3 days of treatment with 1 nM E2 and 6 days of treatment with EtOXTPE, 4OHT, and endoxifen at 1 μM concentrations. (B) effects of 1 μM EtOXTPE on TNFα gene mRNA in MCF7:5C cells after 120 hours of treatment. (C) Effects of 1 μM EtOXTPE on BCL2L11 gene mRNA in MCF7:5C cells after 120 hours of treatment. (D) Inhibition of EtOXTPE-induced apoptosis measured by annexin V staining in MCF7:5C cells by PERK inhibitor. All treatments were performed in triplicate; data represent the average of the replicates.
Activation of Estrogen-Regulated Genes by Nonplanar Estrogen EtOXTPE.
To test the impact of the test compounds on the transcription of the some of the estrogen-regulated genes in MCF-7:5C cells, we treated cells with the test compounds at the 10−6 M concentration for 24 hours and used quantitative RT-PCR (qRT-PCR) as described in Materials and Methods. The results of the experiments show that 1 nM E2 is able to induce mRNA production of TFF1 and GREB1 genes compared with vehicle control (P < 0.05) (Fig. 5). Treatment with test compound EtOXTPE produced mRNA production in both TFF1 and GREB1 genes when compared with vehicle controls (P < 0.05), but only partially compared with E2 treatment (P < 0.05) at 24 hours (Fig. 5).
Effects of well-documented estrogen-responsive gene mRNA expression in MCF-7:5C cells after treatment with test compounds at different time points at 1 nM concentration for E2 and 1 µM for other test compounds. (A) TFF1 gene after 24 hours of treatment; (B) GREB1 gene after 24 hours of treatment. All treatments were performed in triplicate; data represent the average of the replicates.
Delayed Effect on the Regulation of ERα mRNA and Protein by Nonplanar Estrogen EtOXTPE in MCF-7:5C Cells.
To assess the regulation of ERα protein and mRNA levels we employed immunoblotting and qRT-PCR. The results of immunoblotting show that E2 can downregulate ERα protein level considerably as early as 12 hours of treatment with subsequent downregulation maintained (Fig. 6A). Interestingly, nonplanar estrogen EtOXTPE can downregulate ERα protein levels to levels comparable with E2 only after 48 hours of treatment (Fig. 6A). Antiestrogens endoxifen did downregulate the levels of ERα protein compared with vehicle control, but the results were equivalent to the 4OHT treatment (Fig. 6A). Interestingly, 1 μM 4OHT can reverse the effect of 100 nM EtOXTPE downregulation of ERα protein after 24 hours of treatment (Fig. 6B). As for the regulation of ERα mRNA levels with test compounds, we employed qRT-PCR. The results demonstrate that E2 downregulated the mRNA levels as soon as 6 hours of treatment, although not statistically significantly (Fig. 6C). However, after 12 hours of treatment E2 downregulated the mRNA levels by 50% compared with vehicle control (P < 0.05) and stayed consistently low for the rest of the time points (Fig. 6C). Interestingly, EtOXTPE downregulated the ERα mRNA expression only after 24 hours of treatment (P < 0.05 compared with vehicle control) (down by 50%) and even further downregulated the mRNA levels after 36 hours of treatment (P < 0.05 compared with vehicle control; down by 60%) (Fig. 6C). Fulvestrant (ICI 182,780) was used as a positive control for ERα protein degradation in the immunoblotting experiments (Fig. 6, A and B) and as a negative control in the qRT-PCR experiments (Fig. 6C). The results indicate that nonplanar estrogen downregulates the ERα protein and mRNA levels as well as E2 in MCF-7:5C cells; however, this effect was delayed over time.
(A) Effects of treatments with indicated compounds on ERα protein levels in MCF-7:5C cells after different durations of treatments at 1 nM concentration for E2 and 1 µM for other test compounds. (B) Inhibition of ERα protein degradation induced by EtOXTPE and E2 treatment after 36 hours of treatment in MCF-7:5C cells. E2 was used at 1 nM concentration and EtOXTPE was used at 100 nM concentration, which is the minimum concentration at which the compound is able to still downregulate the ERα protein levels in the cells after 36 hours of treatment. 4OHT was used at 1 μM concentration. (C) Time course of treatments with indicated compounds and their effect on ERα mRNA gene expression in MCF-7:5C cells at 10 nM concentration for E2 and 1 µM for other test compounds. (C) Reversal of ERα protein degradation by estrogens after a 48-hour treatment, including EtOXTPE at suboptimal concentration by 1 μM 4OHT. (D) Phosphorylation levels of eIF2α (p-eIF2α) after various treatment durations with indicated compounds or their combinations. E2 and EtOXTPE both induce phosphorylation of eIF2α; however, EtOXTPE has a delayed response. Antiestrogens, such as endoxifen and ICI 182,780, do not induce phosphorylation of eIF2α at any time points, and inhibit activation of eIF2α by E2 or EtOXTPE. At a 24-hour time point EtOXTPE acts as an antiestrogen and inhibits E2-induced phosphorylation of eIF2α. In earlier work it was determined that PERK signaling protein in UPR was crucial for facilitating estrogen-induced apoptosis in MCF-7:5C cells via phosphorylation of eIF2α (Fan et al., 2013). All immunoblots were performed in three replicates; data presented represent one the biologic replicates. Analysis was validated by densitometry using Image J (National Institutes of Health) and the densitometry data are presented in Supplemental Tables 1–3.
Induction of Unfolded Protein Response through the PERK Pathway by Nonplanar Estrogen EtOXTPE in MCF-7:5C Cells.
One of the hallmarks and triggers of estrogen-induced apoptosis in MCF-7:5C cells is the induction of the UPR at early time points of treatment with estrogens (Fan et al., 2013). The unfolded protein response initiates endoplasmic reticulum stress that activates PERK, which subsequently activates eIF2α by phosphorylation. The PERK pathway is essential for triggering estrogen-induced apoptosis. To assess the UPR through the phosphorylation of eIF2α we employed immunoblotting. Our results show that E2 is able to induce phosphorylation of eIF2α as soon as 24 hours when compared with vehicle control (Fig. 6D); however, nonplanar EtOXTPE increases the phosphorylation levels of eIF2α equal to E2 treatment only after 36 hours of treatment when compared with vehicle control (Fig. 6D). Additionally, antiestrogens endoxifen and ICI 182,780 were used as negative controls and did not induce activation of eIF2α and were able to reverse the E2- and EtOXTPE-induced phosphorylation of eIF2α (Fig. 6D). Interestingly, EtOXTPE was also able to reverse E2-induced activation of eIF2α but only at 24 hours of treatment, and at later time points both compounds induced the phosphorylation of eIF2α (Fig. 6D).
Recruitment of ERα and SRC3 on GREB1 Promoter by Nonplanar Estrogen EtOXTPE in MCF-7:5C Cells.
To assess and compare the effects of test compounds on the recruitment of ERα protein and SRC3 coactivator in MCF-7:5C cells we have employed ChIP assay. The ChIP assays were done after 45-minute and 36-hour time points. At the 45-minute time point, EtOXTPE recruits only 25% of ER to the GREB1 promoter compared with E2 (Fig. 7A) and even less for SRC3 after 45 minutes of treatment (Fig. 7B). After 36 hours of treatment, EtOXTPE recruits about one-half of ERα compared with E2 (Fig. 7C). Interestingly, 4OHT also recruited as much ER as EtOXTPE after 36 hours of treatment (Fig. 7C), but after 45 minutes of treatment 4OHT recruited only twice as much as vehicle control (Fig. 7A) but no SRC3 was associated with the ER (Fig. 7, B and D). EtOXTPE recruited very little of the SRC3 coactivator when compared with vehicle control at the 36-hour time point (P > 0.05) and not more than 4OHT (P > 0.05) (Fig. 7D). In summary, EtOXTPE behaves more like an antiestrogen in terms of coactivator recruitment to the estrogen-responsive gene GREB1; however, biologically EtOXTPE is a partial agonist as demonstrated by growth assays and RT-PCR. This can indicate a different mechanism of ER activation for transcription based on a different surface conformation of the ER bound with the EtOXTPE compound.
ChIP assay showing recruitment of ERα and SRC3 at TFF1 estrogen-responsive element promoter after treatment (45 minutes and 36 hours) with test compounds in MCF7:5C cells at 10 nM concentration for E2 and 1 µM for other test compounds. Recruitment of ERα (A) and SRC3 (B) after 45 minutes of treatment with indicated ligands. Recruitment of ERα (C) and SRC3 (D) after 36 hours of treatment with indicated ligands. Recruitment of ERα and SRC3 was calculated as percentage of the total input after subtracting the IgG recruitment. All treatments were performed in triplicate; data represent the average of the replicates.
Novel Conformation of ERα LBD:EtOXTPE Compared with Estrogen, Endoxifen, and 4OHT Resolved by X-Ray Crystallography.
To investigate the structural basis for EtOXTPE binding to the ERα LBD we obtained an X-ray crystal structure of the co-crystallized complex. Clear density was observed in the ligand binding pocket corresponding to the trans-EtOXTPE isomer postrefinement (PBD code 5T1Z) (Fig. 8A). As such, this model shows that the trans-isomer of EtOXTPE is preferred over the cis-isomer in the ERα LBD binding pocket. The ethoxy moiety is not well ordered in the map. This is likely due to the absence of any stabilizing interactions with the LBD, making it free to adopt multiple conformations in the binding pocket. EtOXTPE appears to mainly form nonpolar interactions within the binding pocket. However, the hydrogen bond network formed between the A ring of estradiol, carboxylate group of E353, guanidinium group of R394, and a water molecule is conserved for the analogous phenol of EtOXTPE (Fig. 8, A and B).
Overall structures and close-up views of the binding site of ERα LBD in complex with EtOXTPE (A), E2 (B), endoxifen (C), and 4OHT (D). On the left-hand side of each panel is shown the overall structure of ERα LBD in complex with the ligands represented as spheres with helix 12 depicted in yellow, activation function 2 site colored in blue, and coactivator proteins in the agonist conformations (A and B) shown in green. On the right-hand side is presented a close-up view of the binding site of ERα LBD displaying the ligands and amino acids involved in the formation of hydrogen bonds and hydrophobic contacts. Hydrogen bonds are depicted as blue dashed lines.
To investigate structural differences between the EtOXTPE structure and the planar estrogen E2, an X-ray crystal structure of the ERα LBD Y537S mutant in complex with E2 and glucocorticoid receptor interacting protein 1 was solved to 1.65 Å using molecular replacement (PDB: 5DTA) (Fig. 8B). One dimer was observed in the ASU. Chain B of each structure was chosen for all comparisons because it is not influenced by crystal packing, while chain A shows crystal contacts at helices 11 and 12. Overall, both structures adopt the agonist conformation of the receptor with helix 12 (shown in yellow in Fig. 8, A and B) closing over the opening of the binding pocket. However, clear structural differences are apparent between the E2 and EtOXTPE structures (Fig. 8, A and B; Supplemental Fig. 2A). Notably, H524 is forced outside of the binding pocket by the bulkier phenol group of EtOXTPE and faces the solvent rather than adopting a conformation similar to that for E2 (Supplemental Fig. 2A). Consequently, the formation of a hydrogen bond with the ligand is no longer possible. In the ERα:E2 complex, this hydrogen bond is part of an extended hydrogen bond network, starting at the H524 side chain (helix 11) and terminating at E339 (helix 3) and K531 (helix 11) via E419. Thus, E2 is arrested in the binding pocket and this network induces stability to the E2 agonist conformation of the receptor. This network is no longer formed in the EtOXTPE structure and the receptor’s stability is affected. The imidazole group of H524 forms hydrogen bonds with the backbone of G521 and K520, while the amino group of K520 is involved in formation of a bridge salt with the carboxylate group of E523 (Supplemental Fig. 2A). The side chains of L525 and L540 are shifted to accommodate and participate in Van Der Waals interactions with the ethoxybenzene group of EtOXTPE (Fig. 8A; Supplemental Fig. 2A). Together, these changes propagate to alter the vector of helix 11 such that it is now closer to helix 12 by 1.876 Å (Cα to Cα at C530) (Supplemental Fig. 2A). As a result, the carbonyl oxygen of Y526 breaks its hydrogen bond with the amide nitrogen of C530, placing the phi and psi angles of C530 outside the range for an alpha helix. Together, these structures show that EtOXTPE induces an alternative conformation of the ERα LBD compared with E2.
EtOXTPE and endoxifen are structurally similar except that EtOXTPE is missing the terminal methylamine sidechain and additionally it has a phenolic hydroxyl group. Clear differences are observed between the X-ray crystal structures. Specifically, the bulkier terminal methylamine group of endoxifen forms hydrogen bonds with D351 and V533 (Fig. 8C) and forces helix 12 into the activation function 2 cleft to block coregulator binding, similar to 4OHT structure (Fig. 8, C and D; Supplemental Fig. 2B). The sum of the smaller size and absence of a hydrogen bond with D351 in the EtOXTPE arm, combined with its binding orientation that favors hydrophobic contacts with L536 and L540 of helix 12, enable the compound to stabilize the agonist conformation of the receptor (Fig. 8A), but to a lesser extent that E2. This conformation is not possible for endoxifen, which similar to 4-OHT contains a large ethoxyamino substituent that protrudes out of the binding pocket (Fig. 8D) and precludes the binding to the agonist conformation of ERα LBD, due to steric clashes with helix 12, especially L536 and L540 (Supplemental Fig. 2C). All X-ray crystal structure data collection and refinement statistics are given in Supplemental Table 4.
Discussion
An understanding of the modulation of estrogen-induced apoptosis through the ER has important implications for physiology through bone remodeling in osteoporosis and trophoblast turnover during implantation of the fertilized egg, in the prepared uterine lining (Jordan et al., 2016). These physiologic programs are also exposed during LTED therapy used to treat breast cancer (Jordan, 2015). We and others have previously reported (Song et al., 2001; Lewis et al., 2005a,b) that antiestrogen-resistant breast cancer cells can trigger apoptosis with low-dose estrogen treatment. This laboratory model has clinical significance since low- and high-dose estrogen treatments have antitumor actions following LTED for breast cancer (Jordan, 2014a; Coelingh Bennink et al., 2017).
Here, we address the hypothesis that a class II (angular) TPE-derived estrogen (Jordan et al., 2001) has an initial antiestrogenic response at the LTED breast tumor cell ER (Maximov et al., 2011), but subsequently triggers apoptosis days later. We have observed that despite the slowed triggering of ER-mediated apoptosis by EtOXTPE, X-ray crystallographic studies demonstrate that helix 12 seals the ligand within the LBD in a similar way as observed with E2 (Fig. 8). Nevertheless, the complex with EtOXTPE, referred to as 5T1Z, is distinct from the E2 complex.
The delay in EtOXTPE-induced apoptosis is illustrated in Fig. 2, A and C and compared with endoxifen, an important antiestrogenic metabolite of tamoxifen (Jordan, 2017), recently reported to be effective for salvage therapy following the failure of LTED therapy in breast cancer (Goetz et al., 2017). Endoxifen produces no apoptosis over the 14-day time course, whereas EtOXTPE produced complete apoptosis within the same time period. The X-ray crystallography of endoxifen is reported for the first time and compared and contrasted with 4OHT and EtOXTPE (Fig. 8; Supplemental Fig. 1). Estradiol, the natural ER binding ligand, induces complete closure of the LBD (Fig. 8A) via helix 12, which facilitates binding of coactivators and the formation of the transcriptional complex (Fig. 7). It is important to note, however, that the initial biologic response in the LTED breast cancer cells is growth stimulation with either E2 or EtOXTPE (Fig. 2A). We have addressed this decision-making mechanism to either grow or die in LTED breast cancer cells in an earlier publication (Fan et al., 2015). The interesting observation is that the class II estrogen EtOXTPE has a longer decision-making mechanism (Fig. 2, A and C).
EtOXTPE causes closure of the LBD with helix 12 locking the ligand inside; however, the complex is different than the E2:ER complex, which may account for the rapidity of estradiol-induced apoptosis (Fig. 2C). Coactivator SRC3 was not rapidly recruited to the EtOXTPE complex (Fig. 7). Indeed, the fact that EtOXTPE recruited ER and SRC3 to the TFF1 target gene promoter more like the antiestrogen 4OHT than E2 provides an explanation for the delayed partial agonist actions of EtOXTPE at TFF1 and GREB1 (Fig. 7). Indeed, the related class II (angular) triphenylethylene bisphenol also has impaired recruitment of the ER to the TFF1 promoter (Sengupta et al., 2013), and 25% of SRC3 binding compared with E2 or the planar class I estrogen bisphenol A. This accounts for the partial agonist actions at estrogen target genes in normal cell of the rat anterior pituitary gland (Jordan and Lieberman, 1984; Jordan et al., 1984) and MCF-7 breast cancer cells (Sengupta et al., 2013; Maximov et al., 2014; Obiorah et al., 2014).
Estrogen-induced apoptosis is the ultimate result of a protective mechanism embedded in the human genome to preserve fidelity of reproducing normal cells. The UPR response that triggers PERK signaling through eIF2α is delayed by EtOXTPE (Fig. 6D). In contrast, apoptosis is not triggered by endoxifen (Fig. 2C) over 2-week exposure. Additionally, apoptosis with EtOXTPE is rescued by addition of 4OHT within 48 hours after the addition of estrogenic EtOXTPE (Fig. 2C). These data are consistent with previous results observed with the class II synthetic angular estrogen bisphenol (Obiorah and Jordan, 2014).
In summary, we provide the first studies to report the X-ray crystallography of the potent antiestrogenic metabolite endoxifen (Fig. 8C) and a class II angular estrogen EtOXTPE (Fig. 8A). Both have related structures based on the estrogen TPE (Fig. 1), but contrasting pharmacologic actions in the LTED cells MCF-7:5C; endoxifen does not induce apoptosis or trigger estrogen-responsive gene transcription, whereas EtOXTPE is a partial agonist but is unable to mobilize full agonist responses. Nevertheless, it is the relentless activation of the UPR sensor PERK that eventually initiates delayed apoptosis. The key to immediate efficient estrogen-induced apoptosis in LTED cells is the efficient recruitment of coregulators to the external surface of the TPE:ER complex. This mechanism has its origins in earlier work by the McDonnell group 20 years ago (Paige et al., 1999). A comparison of class II synthetic estrogens, apoptosis, and coregulator recruitment, as previously reported by Han et al. (2016), will aid further understanding of the molecular modulation of estrogen-induced apoptosis via the ER.
These data illustrate the plasticity of the UPR system, which is dependent on the shape of the ligand:ER complex to trigger apoptosis. However, it is the translation of these findings to clinical utility for patient care that holds the most promise. There is a concern that 5 years of adjuvant therapy (Pan et al., 2017) is not sufficient to control disease recurrence for high-risk patients, i.e., in relation to large primary tumors and/or multiple lymph nodes. There is a proposal (Abderrahman and Jordan, 2018) to deploy estrogen-induced apoptosis and other precision medicines as a preemptive salvage therapy to lower micrometastatic tumor burden in high-risk patients. However, medicinal chemists have discovered (Xiong et al., 2016) raloxifene analogs that can occupy the ER LBD and trigger apoptosis without the collateral estrogenic activity. Studies are ongoing to address this hypothesis and a phase I trial is planned.
Authorship Contributions
Participated in research design: Maximov, Jordan.
Conducted experiments: Maximov, Abderrahman, Fanning, Sengupta, Fan, Curpan, Rincon, Greenland, Rajan, Greene.
Performed data analysis: Maximov, Fanning, Curpan, Jordan.
Wrote or contributed to the writing of the manuscript: Maximov, Fanning, Curpan, Jordan.
Footnotes
- Received December 27, 2017.
- Accepted May 4, 2018.
This work was supported by the Department of Defense Breast Center of Excellence Program [Award W81XWH-06-1-0590]; the Susan G. Komen for the Cure Foundation [Award SAC100009]; the National Institutes of Health Lombardi Comprehensive Cancer Center [Grant P30 CA051008]; the National Institutes of Health MD Anderson’s Cancer Center [Grant CA016672]; the Institute of Chemistry Timisoara of the Romanian Academy [Project 1.1/2015]; the Romanian National Authority for Scientific Research and Innovation, CNCS–UEFISCDI [Project PN-II-RU-TE-2014-4-0422]; the Department of Defense BCRP [Breakthrough Award W81XWH-14-1-0360]; the Susan G Komen Postdoctoral Fellowship [Grant PDF14301382]; and the Virginia and D. K. Ludwig Fund for Cancer Research.
↵
This article has supplemental material available at molpharm.aspetjournals.org.
Abbreviations
- 4OHT
- 4-hydroxytamoxifen
- ANOVA
- analysis of variance
- ASU
- asymmetric unit
- ChIP
- chromatin immunoprecipitation
- E2
- 17β-estradiol
- ER
- estrogen receptor
- EtOXTPE
- ethoxytriphenylethylene
- LBD
- ligand-binding domain
- LTED
- long-term estrogen deprivation
- MBC
- metastatic breast cancer
- PDB
- Protein Data Bank
- PEG
- polyethylene glycol
- PERK
- protein kinase regulated by RNA-like endoplasmic reticulum kinase
- qRT-PCR
- quantitative real-time polymerase chain reaction
- RT-PCR
- real-time polymerase chain reaction
- TPE
- triphenylethylene
- UPR
- unfolded protein response
- Copyright © 2018 by The Author(s)
This is an open access article distributed under the CC BY-NC Attribution 4.0 International license.