Identification of amino acids in the hormone binding domain of the human estrogen receptor important in estrogen binding

Nuclear receptor superfamily structural diversity in pacific oyster: In silico identification of estradiol binding candidates

The increasing presence of anthropogenic contaminants in aquatic environments poses challenges for species inhabiting contaminated sites. Due to their structural binding characteristics to ligands that inhibit or activate gene transcription, these xenobiotic compounds frequently target the nuclear receptor superfamily. The present work aims to understand the potential interaction between the hormone 17-β-estradiol, an environmental contaminant, and the nuclear receptors of Crassostrea gigas, the Pacific oyster. This filter-feeding, sessile oyster species is subject to environmental changes and exposure to contaminants. In the Pacific oyster, the estrogen-binding nuclear receptor is not able to bind this hormone as it does in vertebrates. However, another receptor may exhibit responsiveness to estrogen-like molecules and derivatives. We employed high-performance in silico methodologies, including three-dimensional modeling, molecular docking and atomistic molecular dynamics to identify likely binding candidates with the target moecule. Our approach revealed that among the C. gigas nuclear receptor superfamily, candidates with the most favorable interaction with the molecule of interest belonged to the NR1D, NR1H, NR1P, NR2E, NHR42, and NR0B groups. Interestingly, NR1H and NR0B were associated with planktonic/larval life cycle stages, while NR1P, NR2E, and NR0B were associated with sessile/adult life stages. The application of this computational methodological strategy demonstrated high performance in the virtual screening of candidates for binding with the target xenobiotic molecule and can be employed in other studies in the field of ecotoxicology in non-model organisms.

Species Differences in Response to Binding Interactions of Bisphenol A and its Analogs with the Modeled Estrogen Receptor 1 and In Vitro Reporter Gene Assay in Human and Zebrafish

Adverse impacts associated with the interactions of numerous endocrine-disruptor chemicals (EDCs) with estrogen receptor 1 play a pivotal role in reproductive dysfunction. The predictive studies on these interactions thus are crucial in the risk assessment of EDCs but rely heavily on the accuracy of specific protein structure in three dimensions. As the three-dimensional (3D) structure of zebrafish estrogen receptor 1 (zEsr1) is not available, the 3D structure of zEsr1 ligand-binding domain (zEsr1-LBD) was generated using MODELLER and its quality was assessed by the PROCHECK, ERRAT, ProSA, and Verify-3D tools. After the generated model was verified as reliable, bisphenol A and its analogs were docked on the zEsr1-LBD and human estrogen receptor 1 ligand-binding domain (hESR1-LBD) using the Discovery Studio and Autodock Vina programs. The molecular dynamics followed by molecular docking were simulated using the Nanoscale Molecular Dynamics program and compared to those of the in vitro reporter gene assays. Some chemicals were bound with an orientation similar to that of 17β-estradiol in both models and in silico binding energies showed moderate or high correlations with in vitro results (0.33 ≤ r²  ≤ 0.71). Notably, hydrogen bond occupancy during molecular dynamics simulations exhibited a high correlation with in vitro results (r²  ≥ 0.81) in both complexes. These results show that the combined in silico and in vitro approaches is a valuable tool for identifying EDCs in different species, facilitating the assessment of EDC-induced reproductive toxicity. Environ Toxicol Chem 2022;41:2431-2443. © 2022 SETAC.

Experimental, molecular docking and molecular dynamic studies of natural products targeting overexpressed receptors in breast cancer

Natural compounds are proper tools for inhibiting cancer cell proliferation. Hence, the search for these ligands of overexpressed receptors in breast cancer has been a competitive challenge recently and opens new avenues for drug discovery. In this research, we have investigated molecular interactions between natural products and overexpressed receptors in breast cancer using molecular docking and dynamic simulation approaches followed by extraction of the best ligand from Citrus limetta and developing for nanoscale encapsulation composed of soy lecithin using a sonicator machine. The encapsulation process was confirmed by DLS and TEM analyses. Anticancer activity was also examined using MTT method. Among the investigated natural compounds, hesperidin was found to bind to specific targets with stronger binding energy. The molecular dynamics results indicated that the hesperidin-MCL-1 complex is very stable at 310.15 K for 200 ns. The RP-HPLC analysis revealed that the purity of extracted hesperidin was 98.8% with a yield of 1.72%. The results of DLS and TEM showed a strong interaction between hesperidin and lecithin with an entrapped efficiency of 92.02 ± 1.08%. Finally, the cytotoxicity effect of hesperidin was increased against the MDA-MB-231 cell line with an IC50 value of 62.93 μg/mL after encapsulation, whereas no significant effect against the MCF10A cell line. We showed for the first time that hesperidin is a flexible and strong ligand for the MCL-1 receptor. Also, it has the in vitro ability to kill the MDA-MB-231 cell lines without having a significant effect on the MCF10A cell lines. Therefore, hesperidin could be used as a food ingredient to generate functional foods.

Monohaloacetonitriles induce cytotoxicity and exhibit different mode of action in endocrine disruption

Haloacetonitriles are emerging disinfection by-products that can be detected in various aquatic environments. They are cytotoxic, genotoxic, mutagenic, and tumorigenic in vitro and in vivo, but their endocrine-disrupting potency remains unknown. In this study, we examined the agonistic and antagonistic estrogenic and androgenic activities of haloacetonitriles, as well as their cytotoxicity, using a yeast-based reporter assay. We also investigated the interactions of haloacetonitriles with human estrogen receptor alpha (hERα) through molecular docking. We observed that iodoacetonitrile (median lethal dose: 1.96 × 10^-5 M) and bromoacetonitrile (median lethal dose: 1.97 × 10^-5 M) had similar cytotoxicities, which are higher than that of chloroacetonitrile (median lethal dose: 7.16 × 10^-5 M). We observed bromoacetonitrile and chloroacetonitrile elicited estrogenic activity with 10% effective concentrations of 3.30 × 10^-9 M and 2.36 × 10^-9 M, respectively. This finding indicates that bromoacetonitrile and chloroacetonitrile may mimic estrogen signaling through interaction with hERα. Consistent with that result, we identified bromoacetonitrile and chloroacetonitrile interacted with residues in the original estrogen recognition sites of hERα. Our results show that bromoacetonitrile and chloroacetonitrile affect the endocrine-disrupting potency mediated via estrogen receptors by using in vitro assay and molecular docking.

LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor

Protein dynamics is central to all biological processes, including signal transduction, cellular regulation and biological catalysis. Among them, in-depth exploration of ligand-driven protein dynamics contributes to an optimal understanding of protein function, which is particularly relevant to drug discovery. Hence, a wide range of computational tools have been designed to investigate the important dynamic information in proteins. However, performing and analyzing protein dynamics is still challenging due to the complicated operation steps, giving rise to great difficulty, especially for nonexperts. Moreover, there is a lack of web protocol to provide online facility to investigate and visualize ligand-driven protein dynamics. To this end, in this study, we integrated several bioinformatic tools to develop a protocol, named Ligand and Receptor Molecular Dynamics (LARMD, http://chemyang.ccnu.edu.cn/ccb/server/LARMD/ and http://agroda.gzu.edu.cn:9999/ccb/server/LARMD/), for profiling ligand-driven protein dynamics. To be specific, estrogen receptor (ER) was used as a case to reveal ERβ-selective mechanism, which plays a vital role in the treatment of inflammatory diseases and many types of cancers in clinical practice. Two different residues (Ile373/Met421 and Met336/Leu384) in the pocket of ERβ/ERα were the significant determinants for selectivity, especially Met336 of ERβ. The helix H8, helix H11 and H7-H8 loop influenced the migration of selective agonist (WAY-244). These computational results were consistent with the experimental results. Therefore, LARMD provides a user-friendly online protocol to study the dynamic property of protein and to design new ligand or site-directed mutagenesis.

Freshwater snail Parafossarulus striatulus estrogen receptor: Characteristics and expression profiles under lab and field exposure

The modes of action by which putative endocrine disrupting chemicals (EDCs) elicit toxicity in mollusks remains unclear due to our limited understanding of the molluscan endocrine system. We identified and partially characterised the estrogen receptor (ER) of the mollusk Parafossarulus striatulus. The full-length cDNA of the ER of P. striatulus (psER) was isolated and found to have an ORF of 1386 bp which corresponded to 461 amino acids. Phylogenetic analysis revealed that psER is an orthologue of ER of other mollusks. Moreover, the DNA-binding domain, ligand-binding domain, P-box, D-box, and AF2 domain were also identified in psER. Exposure of females and males to 17β-estradiol (E2, 100 ng/L) for 24 h and 72 h did not alter psER transcription, but exposure to 17α-methyltestosterone (MT, 100 μg/L) for 72 h significantly decreased ER transcription in females only (p < 0.05). psER transcription was surveyed in males and females seeded in different regions in Taihu Lake, China. psER transcription were elevated among females and males maintained at site ML. This elevation was statistically significant (p < 0.05) among male snails as compared to snails held at the more pristine site of SZ. This was different to the results from lab, implying that some unknown chemicals or other environmental factors in field could affect psER transcription level in snails. Furthermore, females and males held at site ML also exhibited a significant elevation in vitellogenin transcription as compared to snails held at site SZ, suggesting that vitellogenin production may be directly regulated by psER or co-regulated with psER in this species.

Molecular characterisation of oestrogen receptor ERα and the effects of bisphenol A on its expression during sexual development in the Chinese giant salamander (Andrias davidianus)

The aim of this study was to characterise the molecular structure of the oestrogen receptor ERα and to evaluate the effect of bisphenol A (BPA) on ERα expression during sexual development of the Chinese giant salamander (Andrias davidianus). The ERα cDNA of A. davidianus includes an open reading frame of 1755bp (encoding 584 amino acids), a 219-bp 5' untranslated region (UTR) and a 611-bp 3'UTR. A polyadenylation signal was not found in the 3'UTR. Amino acid sequence analysis showed high homology between ERα of A. davidianus and that of other amphibians, such as Andrias japonicas (99.66% identity) and Rana rugose (81.06% identity). In 3-year-old A. davidianus, highest ERα expression was observed in the liver and gonads. During different developmental stages in A. davidianus (from 1 to 3 years of age), ERα expression in the testes increased gradually. ERα was localised in the epithelial cells of seminiferous lobules and in interstitial cells. ERα-positive cells were more abundant in the interstitial tissue during testicular development. ERα was located in the nucleus of oocytes during ovary development. We found that the sex of 6-month-old A. davidianus larvae could not be distinguished anatomically. The sex ratio did not change after larvae were treated with 10μM BPA for 1 month. However, BPA treatment reduced bodyweight and ERα expression in the gonads in male larvae.

Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists

We present the energy profiles of estrogen receptor–agonist ligand interactions in atomic detail using a quantum biochemical approach.

Estrogenicity of halogenated bisphenol A: in vitro and in silico investigations

The binding interactions of bisphenol A (BPA) and its halogenated derivatives (halogenated BPAs) to human estrogen receptor α ligand binding domain (hERα-LBD) was investigated using a combined in vitro and in silico approach. First, the recombinant hERα-LBD was prepared as a soluble protein in Escherichia coli BL21(DE3)pLysS. A native fluorescent phytoestrogen, coumestrol, was employed as tracer for the fluorescence polarization assay. The results of the in vitro binding assay showed that bisphenol compounds could bind to hERα-LBD as the affinity ligands. All the tested halogenated BPAs exhibited weaker receptor binding than BPA, which might be explained by the steric effect of substituents. Molecular docking studies elucidated that the halogenated BPAs adopted different conformations in the flexible hydrophobic ligand binding pocket (LBP), which is mainly dependent on their distinct halogenation patterns. The compounds with halogen substituents on the phenolic rings and on the bridging alkyl moiety acted as agonists and antagonists for hERα, respectively. Interestingly, all the compounds in the agonist conformation of hERα formed a hydrogen bond with His524, while the compounds in the antagonist conformation formed a hydrogen bond with Thr347. These docking results suggested a pivotal role of His524/Thr347 in maintaining the hERα structure in the biologically active agonist/antagonist conformation. Comparison of the calculated binding energies vs. experimental binding affinities yielded a good correlation, which might be applicable for the structure-based design of novel bisphenol compounds with reduced toxicities and for environmental risk assessment. In addition, based on hERα-LBD as a recognition element, the proposed fluorescence polarization assay may offer an alternative to chromatographic techniques for the multi-residue determination of bisphenol compounds.

Potential mechanisms underlying estrogen-induced expression of the molluscan estrogen receptor (ER) gene

In vertebrates, estrogens and estrogen mimicking chemicals modulate gene expression mainly through a genomic pathway mediated by the estrogen receptors (ERs). Although the existence of an ER orthologue in the mollusc genome has been known for some time, its role in estrogen signalling has yet to be deciphered. This is largely due to its constitutive (ligand-independent) activation and a limited mechanistic understanding of its regulation. To fill this knowledge gap, we cloned and characterised an ER cDNA (sgER) and the 5'-flanking region of the gene from the Sydney rock oyster Saccostrea glomerata. The sgER cDNA is predicted to encode a 477-amino acid protein that contains a DNA-binding domain (DBD) and a ligand-binding domain (LBD) typically conserved among both vertebrate and invertebrate ERs. A comparison of the sgER LBD sequence with those of other ligand-dependent ERs revealed that the sgER LBD is variable at several conserved residues known to be critical for ligand binding and receptor activation. Ligand binding assays using fluorescent-labelled E2 and purified sgER protein confirmed that sgER is devoid of estrogen binding. In silico analysis of the sgER 5'-flanking sequence indicated the presence of three putative estrogen responsive element (ERE) half-sites and several putative sites for ER-interacting transcription factors, suggesting that the sgER promoter may be autoregulated by its own gene product. sgER mRNA is ubiquitously expressed in adult oyster tissues, with the highest expression found in the ovary. Ovarian expression of sgER mRNA was significantly upregulated following in vitro and in vivo exposure to 17β-estradiol (E2). Notably, the activation of sgER expression by E2 in vitro was abolished by the specific ER antagonist ICI 182, 780. To determine whether sgER expression is epigenetically regulated, the in vivo DNA methylation status of the putative proximal promoter in ovarian tissues was assessed using bisulfite genomic sequencing. The results showed that the promoter is predominantly hypomethylated (with 0-3.3% methylcytosines) regardless of sgER mRNA levels. Overall, our investigations suggest that the estrogen responsiveness of sgER is regulated by a novel ligand-dependent receptor, presumably via a non-genomic pathway(s) of estrogen signalling.

Characterization of agonist binding to His524 in the estrogen receptor α ligand binding domain

The bioactivities of the natural steroidal estrogen 17β-estradiol (E2), the synthetic estrogen diethylstilbestrol (DES), and the phytoestrogen genistein (GEN) are intimately associated with their binding to the estrogen receptor α ligand binding domain (ERα LBD) and accordingly allostery. Molecular modeling techniques have been performed on agonists in complex with the LBD, focusing on the pivotal role of His524 modeled as the ε-tautomer and the protonated form (depending on pH). It is found that E2 binds to the active LBD with the aid of Leu525, showing existing stable patterns of an H-binding network with Glu419 via His524 in all models. The main difference seen in the effect is that the full agonists E2 and DES have higher binding energies to the protonated His524 than the partial agonists GEN and Way-169916 (W), which is in line with noted experimental transcriptional activities. In conclusion, the study demonstrates that the phytoestrogen GEN interacts differently with the LBD than what E2 and DES do, which explains the observed signaling differences.

Freshwater mudsnail (Potamopyrgus antipodarum) estrogen receptor: identification and expression analysis under exposure to (xeno-)hormones

Molluscs are raising attention as ecotoxicological test organisms due to their high diversity and ecological importance. The ovoviviparous prosobranch gastropod Potamopyrgus antipodarum (freshwater mudsnail) responds very sensitively to xenobiotics and has therefore been proposed as OECD standard test organism. Endocrine disrupting chemicals influence the reproduction of P. antipodarum, which can be assessed by embryo numbers in the brood pouch. However, the knowledge about the endocrine system of P. antipodarum is rather limited. The aim of this study was to identify an estrogen receptor in the endocrine system of P. antipodarum and to investigate if this receptor is differentially expressed under exposure to (xeno-)hormones (17α-ethinylestradiol, bisphenol A and 17α-methyltestosterone). The DNA-binding domain of the identified ER-like transcript has an amino acid identity of 92 percent compared to the ER of the gastropod Nucella lapillus (84 percent to human ERα) and 83 percent in the ligand binding domain (38 percent to human ERα). Furthermore, the P. antipodarum ER is transcriptionally regulated as shown by quantitative real-time PCRs of (xeno-)hormone exposed snails. 17α-ethinylestradiol and bisphenol A exposure resulted in a transitory ER-mRNA increase while17α-methyltestosterone caused a transitory reduction of ER-mRNA. In addition the solvent dimethyl sulfoxide had also a modulating effect on the receptor.

RUNX3 acts as a tumor suppressor in breast cancer by targeting estrogen receptor α

Transcription factor RUNX3 is inactivated in a number of malignancies, including breast cancer, and is suggested to function as a tumor suppressor. How RUNX3 functions as a tumor suppressor in breast cancer remains undefined. Here, we show that about 20% of female Runx3(+/-) mice spontaneously developed ductal carcinoma at an average age of 14.5 months. Additionally, RUNX3 inhibits the estrogen-dependent proliferation and transformation potential of ERα-positive MCF-7 breast cancer cells in liquid culture and in soft agar and suppresses the tumorigenicity of MCF-7 cells in severe combined immunodeficiency mice. Furthermore, RUNX3 inhibits ERα-dependent transactivation by reducing the stability of ERα. Consistent with its ability to regulate the levels of ERα, expression of RUNX3 inversely correlates with the expression of ERα in breast cancer cell lines, human breast cancer tissues and Runx3(+/-) mouse mammary tumors. By destabilizing ERα, RUNX3 acts as a novel tumor suppressor in breast cancer.

Molecular Cloning and Characterization of the Estrogen Receptor from the Striped Bitterling (Acheilognathus yamatsutae)

OBJECTIVES

In order to identify the possibility of striped bitterling (SB) (Acheilognathus yamatsutae) being used as a test species for estrogenic endocrine disrupting chemicals (EEDCs), we carried out the cloning and sequence characterization of the estrogen receptor (ER).

METHODS

The ER from a striped bitterling was obtained by reverse transcriptase-polymerase chain reaction (RT-PCR), 5'- and 3'-rapid amplification of cDNA ends (5'-RACE and 3'-RACE) and T-vector cloning. The expression of ER mRNA was also analyzed in six tissues (brain, liver, kidney, gill, gonad, and intestines) by real-time PCR.

RESULTS

We obtained an ER from the striped bitterling. The SB ER cDNA was 2189 base pairs (bp) in length and contained a 1707 bp open reading frame that encoded 568 amino acid residues. The SB ER amino acid sequence clustered in a monophyletic group with the ERα of other fish, and was more closely related to zebrafish ERα (88% identity) than to the ERα of other fish. The SB ER cDNA was divided into A/B, C, D, E and F domains. The SB ER has conserved important sequences for ER functions, such as the DNA binding domain (D domain), which are consistent with those of other teleosts.

CONCLUSIONS

The ER of the striped bitterling could provide basic information in toxicological studies of EEDCs in the striped bitterling.

Molecular characterization of estrogen receptor genes in Gobiocypris rarus and their expression upon endocrine disrupting chemicals exposure in juveniles

Estrogens play an important role in many physiological processes of vertebrates, mediated by estrogen receptors (ERs). The full length of the cDNAs for ERα, ERβ1, and ERβ2 were isolated and characterized from Gobiocypris rarus. G. rarus ERs shared the highest amino acid identities with counterparts of three cyprinidae species (Pimephales promelas ERα: 91.1%, Rutilus rutilus ERβ1: 92.9%, Tanichthy albonubes ERβ2: 93.5%). The phylogenic tree of vertebrate ERs indicates G. rarus ER isoforms are more related to counterparts of cyprinidae species. The expression of ERα mRNA was high in gonad and liver. The ERβ1 transcript was the highest in the liver of female fish and was evenly high in the liver, testis and intestine in male. The ERβ2 transcript was high in liver, gonad, and intestine. G. rarus juvenile at 34 days post fertilization were exposed for 3 days to endocrine disrupting chemicals including 17α-ethynylestradiol (EE2), 4-nonylphenol (NP) and bisphenol A (BPA). ER mRNA expression following the xenoestrogens' exposure was analyzed by quantitative real-time PCR. EE2 exposure at 0.01, 0.1 and 1 nM significantly up-regulated ERα transcript. ERβ1 mRNA expression was suppressed by EE2 at all concentrations. However ERβ2 transcript had opposite response to EE2 at low and high concentrations (up-regulation at 0.1 nM, down-regulation at 1 nM). Except a weak increase of ERα at 10 nM EE2, varying decrease of three ER transcripts was resulted in by NP at 10, 100 and 1000 nM. ERα transcript was significantly up-regulated by BPA at 10 nM. A non-significant weak increase in ERβ1 mRNA expression was caused by 1 nM BPA. However 1 nM and 10 nM BPA exposures resulted in significant and non-significant decrease of ERβ2 transcript, respectively. The BPA exposures at other concentrations almost had no effect on the ER transcripts. Vitellogenin (Vtg) mRNA expression profiling following exposure to three xenoestrogens indicated that Vtg transcript is a sensitive biomarker of the juvenile G. rarus at 34 dpf to the EDCs, especially to EE2. These results combined suggest that the ER genes are not modulated in the same manner by EE2, NP, and BPA and that ERs may not contribute equally to the transcriptional regulation of genes involved in fish development and reproduction.

Nucleotide sequence, tissue expression patterns and phylogenetic analysis of estrogen receptor one mRNA in the Murray rainbowfish (Melanotaenia fluviatilis) (Atheriniformes, Actinopterygii)

Estrogens are steroidal hormones that control many physiological processes in both female and male vertebrates. Like other vertebrates, fish have two distinct estrogen receptors (Esr) subtypes, Esr1 and Esr2a that have been isolated in a number of species, as well as a third subtype, Esr2b. The mRNA encoding the Esr1 was isolated from the female liver of an Australian freshwater fish, the Murray rainbowfish, Melanotaenia fluviatilis. The rainbowfish esr1 cDNA was 2569 bp in length and with an open reading frame to encode a protein of 611 amino acids. Phylogenetic analysis and multiple amino acid sequence alignment indicated close relationship and high similarity with killifish (Fundulus heteroclitus) and gilthead sea bream (Sparus aurata). Expression of rainbowfish esr1 mRNA was abundant in the liver, gonads and intestine of adult female and male rainbowfish. This is the first isolation of the full-length nucleotide sequence of an estrogen receptor from rainbowfish. This sequence provides a valuable molecular tool that can be used in future studies investigating estrogen mechanisms, actions and tissue-specific expression in juvenile and adult rainbowfish.

Computational insights into the mechanism of ligand unbinding and selectivity of estrogen receptors

Estrogen receptors (ER) belong to the nuclear receptor superfamily, and two subtypes, ERalpha and ERbeta, have been identified to date. The differentiated functions and receptor expressions of ERalpha and ERbeta made it attracted to discover subtype-specified ligands with high selectivity. However, these two subtypes are highly homologous and only two residues differ in the ligand binding pocket. Therefore, the mechanism of ligand selectivity has become an important issue in searching selective ligands of ER subtypes. In this study, steered molecular dynamics simulations were carried out to investigate the unbinding pathways of two selective ERbeta ligands from the binding pocket of both ERalpha and ERbeta, which demonstrated that the pathway between the H11 helix and the H7 approximately H8 loop was the most probable for ligand escaping. Then potentials of mean force for ligands unbinding along this pathway were calculated in order to gain insights into the molecular basis for energetics of ligand unbinding and find clues of ligand selectivity. The results indicated that His524/475 in ERalpha/ERbeta acted as a "gatekeeper" during the ligand unbinding. Especially, the H7 approximately H8 loop of ERbeta acted as a polar "transmitter" that controlled the ligand unbinding from the binding site and contributed to the ligand selectivity. Finally, the mechanism of ligand selectivity of ER subtypes was discussed from a kinetic perspective and suggestions for improving the ligand selectivity of ERbeta were also presented. These findings could be helpful for rational design of highly selective ERbeta ligands.

An amphioxus orthologue of the estrogen receptor that does not bind estradiol: insights into estrogen receptor evolution

Background

The origin of nuclear receptors (NRs) and the question whether the ancestral NR was a liganded or an unliganded transcription factor has been recently debated. To obtain insight into the evolution of the ligand binding ability of estrogen receptors (ER), we comparatively characterized the ER from the protochordate amphioxus (Branchiostoma floridae), and the ER from lamprey (Petromyzon marinus), a basal vertebrate.

Results

Extensive phylogenetic studies as well as signature analysis allowed us to confirm that the amphioxus ER (amphiER) and the lamprey ER (lampER) belong to the ER group. LampER behaves as a "classical" vertebrate ER, as it binds to specific DNA Estrogen Responsive Elements (EREs), and is activated by estradiol (E₂), the classical ER natural ligand. In contrast, we found that although amphiER binds EREs, it is unable to bind E₂ and to activate transcription in response to E₂. Among the 7 natural and synthetic ER ligands tested as well as a large repertoire of 14 cholesterol derivatives, only Bisphenol A (an endocrine disruptor with estrogenic activity) bound to amphiER, suggesting that a ligand binding pocket exists within the receptor. Parsimony analysis considering all available ER sequences suggest that the ancestral ER was not able to bind E₂ and that this ability evolved specifically in the vertebrate lineage. This result does not support a previous analysis based on ancestral sequence reconstruction that proposed the ancestral steroid receptor to bind estradiol. We show that biased taxonomic sampling can alter the calculation of ancestral sequence and that the previous result might stem from a high proportion of vertebrate ERs in the dataset used to compute the ancestral sequence.

Conclusion

Taken together, our results highlight the importance of comparative experimental approaches vs ancestral reconstructions for the evolutionary study of endocrine systems: comparative analysis of extant ERs suggests that the ancestral ER did not bind estradiol and that it gained the ability to be regulated by estradiol specifically in the vertebrate lineage, before lamprey split.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17beta-estradiol (E(2)) and the selective ER modulator raloxifene (RAL)] from the human ERalpha ligand-binding domain in monomeric and dimeric forms. E(2) dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II'). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I'). Remarkably, ER dimerization strongly suppresses Paths II and II' for E(2) dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

A human estrogen receptor (ER)alpha mutation with differential responsiveness to nonsteroidal ligands: novel approaches for studying mechanism of ER action

Estrogens, acting through the estrogen receptors (ERs), play crucial roles in regulating the function of reproductive and other systems under physiological and pathological conditions. ER activity in regulating target genes is modulated by the binding of both steroidal and synthetic nonsteroidal ligands, with ligand binding inducing ERs to adopt various conformations that control their interactions with transcriptional coregulators. Previously, we developed an intramolecular folding sensor with a mutant form of ERalpha (ER(G521T)) that proved to be essentially unresponsive to the endogenous ligand 17beta-estradiol, yet responded very well to certain synthetic ligands. In this study, we have characterized this G521T-ER mutation in terms of the potency and efficacy of receptor response toward several steroidal and nonsteroidal ligands in two different ways: directly, by ligand effects on mutant ER conformation (by the split-luciferase complementation system), and indirectly, by ligand effects on mutant ER transactivation. Full-length G521T-ER shows no affinity for estradiol and does not activate an estrogen-responsive reporter gene. The synthetic pyrazole agonist ligand propyl-pyrazole-triol is approximately 100-fold more potent than estradiol in inducing intramolecular folding and reporter gene transactivation with the mutant ER, whereas both ligands have high potency on wild-type ER. This estradiol-unresponsive mutant ER can also specifically highlight the agonistic property of the selective ER modulator, 4-hydroxytamoxifen, by reporter gene transactivation, even in the presence of estradiol, and it can exert a dominant-negative effect on estrogen-stimulated wild-type ER. This system provides a model for ER-mutants that show differential ligand responsiveness to gene activation to gain insight into the phenomenon of hormone resistance observed in endocrine therapies of ER-positive breast cancers.

NFkappaB selectivity of estrogen receptor ligands revealed by comparative crystallographic analyses

Our understanding of how steroid hormones regulate physiological functions has been significantly advanced by structural biology approaches. However, progress has been hampered by misfolding of the ligand binding domains in heterologous expression systems and by conformational flexibility that interferes with crystallization. Here, we show that protein folding problems that are common to steroid hormone receptors are circumvented by mutations that stabilize well-characterized conformations of the receptor. We use this approach to present the structure of an apo steroid receptor that reveals a ligand-accessible channel allowing soaking of preformed crystals. Furthermore, crystallization of different pharmacological classes of compounds allowed us to define the structural basis of NFkappaB-selective signaling through the estrogen receptor, thus revealing a unique conformation of the receptor that allows selective suppression of inflammatory gene expression. The ability to crystallize many receptor-ligand complexes with distinct pharmacophores allows one to define structural features of signaling specificity that would not be apparent in a single structure.

Stimulation of transactivation of the largemouth bass estrogen receptors alpha, beta-a, and beta-b by methoxychlor and its mono- and bis-demethylated metabolites in HepG2 cells

The purpose of this study was to determine the mechanisms by which the pesticide, methoxychlor (MXC), acts as an environmental endocrine disruptor through interaction with the three largemouth bass (Micropterus salmoides) estrogen receptors (ERs) alpha, betaa, and betab. MXC is a less-environmentally persistent analog of DDT that behaves as a weak estrogen. Using transient transfection assays in HepG2 cells, we have previously shown that each receptor is responsive to the endogenous ligand 17beta-estradiol (E(2)) in a dose-dependent manner. The parent compound, MXC, showed dose-dependent stimulation of transcriptional activation through all three ERs. In addition to the parent molecule, each of the metabolites was also estrogenic with all three ERs. The order of potency for ERalpha and ERbetab was HPTE>OH-MXC>MXC, while the opposite order was seen for ERbetaa. HepG2 cells did not substantially metabolize MXC to the active metabolites, thus the activity of MXC was not due to metabolism. When examining the effects of increasing concentrations of MXC at a fixed concentration of E(2), all three ERs show increased activity compared to that with E(2) alone, showing that the effects of MXC and E(2) are additive. However, when this experiment was repeated with increasing concentrations of HPTE at a fixed concentration of E(2), the activity of ERalpha was decreased, that of ERbetab was increased, while that of ERbetaa was unaffected compared to E(2) alone. These experiments suggest that HPTE functions as an E(2) antagonist with ERalpha, an E(2) agonist with ERbetab and does not perturb E(2) stimulation of ERbetaa. While it is clear the ERbeta subtypes are the products of different genes (due to a gene duplication in teleosts) the differences in their responses to MXC and its metabolites indicate that their functions diverge, both in their in vivo molecular response to E(2), as well as in their interaction with endocrine disrupting compounds found in the wild.

A novel ligand-independent function of the estrogen receptor is essential for osteocyte and osteoblast mechanotransduction

Bone senses and adapts to meet mechanical needs by means of an extensive mechanotransduction network comprising osteocytes (former osteoblasts entrapped in mineral) and their cytoplasmic projections through which osteocytes communicate with osteoblasts and osteoclasts on the bone surface. Mechanical stimulation promotes osteocyte (and osteoblast) survival by activating the extracellular signal-regulated kinases, ERKs. Estrogens have similar effects and, intriguingly, the adaptive response of bone to mechanical forces is defective in mice lacking estrogen receptor (ER) alpha or ERbeta. We report that ERKs are not activated by stretching in osteocytic and osteoblastic cells in which both ERalpha and ERbeta have been knocked out or knocked down and this is reversed partially by transfection of either one of the two human ERs and fully by transfection of both receptors. ERK activation in response to stretching is also recovered by transfecting the ligand-binding domain (E) of either receptor or an ERalpha mutant that does not bind estrogens. Furthermore, mechano-responsiveness is restored by transfecting the Ealpha targeted to the plasma membrane, but not to the nucleus, whereas ERalpha mutants with impaired plasma membrane localization or binding to caveolin-1 fail to confer ERK activation in response to stretching. Lastly, the ER antagonist ICI 182,780 abrogates ERK activation and the anti-apoptotic effect of mechanical stimulation. We conclude that in addition to their role as ligand-dependent mediators of the effects of estrogens, the ERs participate in the transduction of mechanical forces into pro-survival signaling in bone cells, albeit in a ligand-independent manner.

Mass spectrometry identification of covalent attachment sites of two related estrogenic ligands on human estrogen receptor alpha

A purified preparation of human estrogen receptor alpha (hERalpha) ligand-binding domain (LBD) involving mainly the Ser(309)Ala(569) (approximately 30%) and Ser(309)Ala(571) (approximately 63%) ER portions was used to identify the covalent attachment sites of two closely related estrogenic ER affinity labels 17alpha-bromoacetamidopropylestradiol (17BAPE(2)) and 17alpha-bromoacetamidomethylestradiol (17BAME(2)). To identify and quantify the electrophile covalent attachment sites, [(14)C]17BAPE(2)- and [(14)C]17BAME(2)-alkylated hLBD preparations were trypsinized and submitted to HPLC. In each case, two radioactive fractions were obtained. Mass spectrometry analyses of the two fractions showed signals, which closely matched the molecular masses of alkylated Cys(530)Lys(531) and Cys(417)Arg(434) hLBD tryptic peptides. The covalent attachment of the two electrophiles on hLBD was assigned to the S atoms of Cys(530) and Cys(417). However, the balance between Cys(530) and Cys(417) labeling markedly differed according to the affinity label used, with the Cys(530)/Cys(417) ratio being 2.1 for 17BAPE(2), and 20 for 17BAME(2). We attempted to interpret the covalent attachment of electrophiles by molecular modeling using the crystallographic structure of LBD bound to E(2). In agreement with the different levels of Cys(417) alkylation, the LBD model with unchanged helices could not easily account for Cys(417) labeling by 17BAME(2), whereas favorable results were obtained through 17BAPE(2) docking. Moreover, labeling at Cys(530) by the two electrophiles could not be interpreted using the LBD model. This indicates that some states of solute LBD bound to the estrogenic E(2) 17alpha-derivatives differ from the structure of crystallized LBD bound to E(2).

Molecular Characterization of Estrogen Receptors 1, 2a, and 2b and Their Tissue and Ontogenic Expression Profiles in Fathead Minnow (Pimephales promelas)1

There are two estrogen receptor (ER) subtypes in fish, Esr1 and Esr2 (formerly ERalpha and ERbeta), and in some species the Esr2 subtype has two forms, Esr2b (formerly ERbeta1) and Esr2a (formerly ERbeta2 or ERgamma). There is little information, however, on the different characteristics and functional significance of the two receptor subtypes in fish, and this is especially relevant for understanding the disruption of ER signaling by chemicals with estrogenic activity. In this study, the full-length cDNAs for esr1 (3167 base pairs [bp]) and esr2b (2318 bp), and a partial-length (267 bp) cDNA for esr2a, were cloned and characterized in fathead minnow (fhm; Pimephales promelas), and their patterns of expression established during development and in adults. Real-time polymerase chain reaction revealed some clear distinctions in the ontogenic and tissue expression of fhm esr1, esr2b, and esr2a, suggesting different functions for each ER subtype. Fhm ERs were expressed in brain, pituitary, liver, gonad, intestine, and gill of male and female fish, esr2b and esr2a were also expressed in muscle. Fhm esr1 and esr2b were expressed predominantly in the liver, whereas fhm esr2a was expressed predominantly in intestine and was lowest expressed in liver. Responses of the different hepatic ERs in male fathead minnow exposed to 100 ng estradiol/L differed, with a significant induction (5-fold) of fhm esr1 but no effect on esr2b or esr2a expression, suggesting different mechanisms of regulation for the different ERs. The detailed characterization of ERs in fathead minnow provides the foundation for understanding the molecular basis of estrogenic disruption in fish.

Regulation of protein activity with small-molecule-controlled inteins

Inteins are the protein analogs of self-splicing RNA introns, as they post-translationally excise themselves from a variety of protein hosts. Intein insertion abolishes, in general, the activity of its host protein, which is subsequently restored upon intein excision. These protein elements therefore have the potential to be used as general molecular "switches" for the control of arbitrary target proteins. Based on rational design, an intein-based protein switch has been constructed whose splicing activity is conditionally triggered in vivo by the presence of thyroid hormone or synthetic analogs. This modified intein was used in Escherichia coli to demonstrate that a number of different proteins can be inactivated by intein insertion and then reactivated by the addition of thyroid hormone via ligand-induced splicing. This conditional activation was also found to occur in a dose-dependent manner. Rational protein engineering was then combined with genetic selection to evolve an additional intein whose activity is controlled by the presence of synthetic estrogen ligands. The ability to regulate protein function post-translationally through the use of ligand-controlled intein splicing will most likely find applications in metabolic engineering, drug discovery and delivery, biosensing, molecular computation, as well as many additional areas of biotechnology.

A novel pesticide-induced conformational state of the oestrogen receptor ligand-binding domain, detected by conformation-specific peptide binding

The diverse effects of different natural and synthetic oestrogen receptor ligands depend on induction of different receptor conformations, allowing differential interactions with other transcription factors. Different conformations of the oestrogen receptor ligand binding domains can be monitored by conformation-specific binding to peptides selected from phage-displayed peptide libraries. We now report that a group of chlorinated pesticides, including 2,4-dichlorodiphenyl-dichloroethylene, induces a peptide recognition pattern different from those induced by any one of the classical oestrogen receptor ligands. The pesticide-complexed oestrogen receptors recognized peptides reacting with the receptors complexed both with the natural oestrogen 17beta-oestradiol and with the synthetic partial antagonist 4-hydroxy-tamoxifen, respectively, indicating that the pesticide-induced conformation shares features with both the 17beta-oestradiol- and the 4-hydroxy-tamoxifen-induced conformations. The substitution H524A in the ligand binding domain conferred the pesticide-specific peptide recognition pattern and transactivation activity to the oestradiol- and the 4-hydroxy-tamoxifen-complexed receptors, indicating that one important determinant of the pesticide-induced conformation is a lack of stabilisation of any one particular receptor conformation by ligand interaction with H524, which is known to interact with both oestradiol and 4-hydroxy-tamoxifen. Thus, peptide binding analyses of oestrogen receptor conformations induced by environmental endocrine disruptors can give novel information about molecular mechanisms of oestrogen action in general.

Estrogen receptor mutations in human disease

As early as the 1800s, the actions of estrogen have been implicated in the development and progression of breast cancer. The estrogen receptor (ER) was identified in the late 1950s and purified a few years later. However, it was not until the 1980s that the first ER was molecularly cloned, and in the mid 1990s, a second ER was cloned. These two related receptors are now called ERalpha and ERbeta, respectively. Since their discovery, much research has focused on identifying alterations within the coding sequence of these receptors in clinical samples. As a result, a large number of naturally occurring splice variants of both ERalpha and ERbeta have been identified in normal epithelium and diseased or cancerous tissues. In contrast, only a few point mutations have been identified in human patient samples from a variety of disease states, including breast cancer, endometrial cancer, and psychiatric diseases. To elucidate the mechanism of action for these variant isoforms or mutant receptors, experimental mutagenesis has been used to analyze the function of distinct amino acid residues in the ERs. This review will focus on ERalpha and ERbeta alterations in breast cancer.

Requirements for reliable determination of binding affinity constants by saturation analysis approach

Accurate calculation of the equilibrium association constant (K) and binding site concentration (N) related to a receptor (R)/ligand (L) interaction, via R saturation analysis, requires exact determination of the specifically bound L concentration (B(S)) and the unbound L concentration (U) at equilibrium. However, most binding determinations involve a procedure for separation of bound and unbound L. In such situations, it was previously shown that correct calculation of B(S) and U from binding data requires prior determination of alpha, i.e. the procedure parameter representing the proportion of equilibrium B(S) recovered after running the separation process, and of kn, i.e. the equilibrium nonspecific binding coefficient. For the simplest model of R/L interaction, the consequences of alpha neglect and/or kn neglect on determination of K and N, via R saturation analysis, are investigated. When alpha but not kn has been determined, B(S) can be accurately calculated, whereas U is overestimated by factor (kn + 1). Consequently the type (linear or hyperbolic) of theoretic curves obtained by usual representations (such as the Scatchard, the Lineweaver-Burk or the Michaelis-Menten plot) of the R/L binding is unchanged; these curves afford correct N and underestimation of K by factor (kn + 1). When alpha (alpha < 1) has not been determined B(S) and U are underestimated and overestimated, respectively. Then erroneous representations of the R/L binding result (e.g. instead of regular straight line segments, Scatchard plot and Lineweaver-Burk plot involve convex-upward and convex-downward hyperbola portions, respectively, suggestive of positive cooperativity of L binding), which leads to incorrect N and K. Errors in N and K would depend on (i) the binding (K, N and kn) and method (alpha) parameters and (ii) the expressions used to calculate approximate B(S) and U values. Simulations involving variable alpha, KN and kn values indicate that: (1) the magnitude of error in N determination (mainly involving moderate underestimation) directly depends on the alpha value; (2) the magnitude of K underestimation mainly depends on the KN value; it is moderate (usually < two-fold) with KN values < 1, but could become very high (e.g. > 100-fold), when KN > 10(2). In this case, the K underestimation is modulated by the alpha and kn values. Practical situations which afford high KN and thus might result in very marked underestimation of K are discussed. A single R dilution method is proposed to assess the validity of K determinations using the R saturation analysis approach.

Cloning, characterisation, and expression of three oestrogen receptors (ERalpha, ERbeta1 and ERbeta2) in the European sea bass, Dicentrarchus labrax

Three oestrogen receptor [ER] subtypes have been described in teleost fish, namely ERalpha, and two ERbeta subtypes, called ERbeta1 and ERbeta2 (or ERbeta and ERgamma in Atlantic croaker). Their expression during embryonic development and gonadal growth has evoked interest in their potential role in sexual differentiation and gonadal development in fish. We cloned three oestrogen receptors from adult liver (sb-ERalpha cDNA) and ovary (partial sb-ERbeta1 and sb-ERbeta2 cDNAs) of the European sea bass, and according to their phylogenetic relatedness to other ERs in teleosts, named them sea bass [sb-] ERalpha, ERbeta1 and ERbeta2. Deduced amino acid numbers for sb-ERalpha, sb-ERbeta1 and sb-ERbeta2 were 639, 517 and 608, respectively, representing in the case of sb-ERbeta1 and sb-ERbeta2 about 90% of the open reading frame. Highest amino acid identities were found for sb-ERalpha with eelpout ERalpha (88.7%), for sb-ERbeta1 with Atlantic croaker ERgamma (85.8%), and for sb-ERbeta2 with Atlantic croaker ERbeta (90.1%). Southern analysis confirmed that all three sea bass oestrogen receptors (sb-ERs) are the products of three distinct genes. In adult sea bass, ERalpha was predominantly expressed in liver and pituitary, while sb-ERbeta1 and sb-ERbeta2 were more ubiquitously expressed, with highest expression levels in pituitary. In a mixed-sex population of juvenile sea bass, sb-ERalpha expression was significantly elevated in gonads at 200 days posthatch (dph), while for sb-ERbeta1 and sb-ERbeta2 highest expression levels were observed in gonads at 250 dph. For sb-ERbeta2, expression was also significantly higher in the brain at 250 dph. The cloning of these three ER subtypes in the European sea bass together with the results obtained on expression levels in adult and juvenile animals has given us the foundation to investigate their possible role in sexual differentiation and development in this species in future studies.

Estrogenic Activities of Nitrophenols in Diesel Exhaust Particles1

We recently isolated 3-methyl-4-nitrophenol (4-nitro-m-cresol; PNMC) and 4-nitro-3-phenylphenol (PNMPP) from diesel exhaust particles (DEP) and identified them as vasodilators. Because these compounds are alkylphenolic derivatives that might mimic hormones, we evaluated their estrogenic activity by using recombinant yeast screens, myometrial contractility assays, and in vivo uterotrophic assays. Recombinant yeast screen assays showed that both PNMC and PNMPP possess estrogenic activity. Furthermore, ovariectomized 25-day-old immature female rats injected with PNMC and PNMPP subcutaneously for 2 days showed significant increases in uterine weight among those receiving 100 mg/kg PNMC and 0.1 and 1.0 mg/kg PNMPP. To clarify further the estrogenic activity of PNMC and PNMPP, rat uterine horns were monitored in organ bath chambers for myometrial contractility in response to oxytocin (OT). Significant differences occurred in the initial and maximum contractilities to OT at 0.25 and 25 mIU/ml in uterine horns obtained from animals treated with 100 mg/kg PNMC and in the maximum contractilities to OT at 0.025, 0.25, and 25 mIU/ml in those from rats treated with 0.1 mg/kg PNMPP. These results clearly demonstrated that PNMC and PNMPP in DEP have estrogenic activity both in vitro and in vivo and might therefore be considered as endocrine-disrupting chemicals.

Differential expression of largemouth bass (Micropterus salmoides) estrogen receptor isotypes alpha, beta, and gamma by estradiol

The expression levels of three estrogen receptor (ER) isotypes alpha, beta, and gamma were quantified in female largemouth bass (Micropterus salmoides) (LMB) liver, ovary, brain, and pituitary tissues. ER alpha and beta expression predominated in the liver, while ERs beta and gamma predominated in the other tissues. Temporally in females, ER alpha was highly up-regulated, ER gamma was slightly up-regulated, and ER beta levels remained unchanged in the liver when plasma 17-beta estradiol (E2) and vitellogenin (Vtg) levels were elevated in the spring. In ovarian tissue from these same fish, all three ERs were maximally expressed in the fall, during early oocyte development and prior to peak plasma E2 levels. When males were injected with E2, ER alpha was highly inducible, ER gamma was moderately up-regulated, and ER beta levels were not affected. None of the ER isotypes were induced by E2 in gonadal tissues. These results combined suggest that the ERs themselves are not regulated in the same manner by E2, and furthermore, do not contribute equally to the transcriptional regulation of genes involved in fish reproduction such as Vtg.

Selective estrogen receptor modulators: discrimination of agonistic versus antagonistic activities by gene expression profiling in breast cancer cells

Selective estrogen receptor modulators (SERMs) such as tamoxifen are effective in the treatment of many estrogen receptor-positive breast cancers and have also proven to be effective in the prevention of breast cancer in women at high risk for the disease. The comparative abilities of tamoxifen versus raloxifene in breast cancer prevention are currently being compared in the Study of Tamoxifen and Raloxifene trial. To better understand the actions of these compounds in breast cancer, we have examined their effects on the expression of approximately 12,000 genes, using Affymetrix GeneChip microarrays, with quantitative PCR verification in many cases, categorizing their actions as agonist, antagonist, or partial agonist/antagonist. Analysis of gene stimulation and inhibition by the SERMs trans-hydroxytamoxifen (TOT) and raloxifene (Ral) or ICI 182,780 (ICI) and by estradiol (E2) in estrogen receptor-containing MCF-7 human breast cancer cells revealed that (a) TOT was the most E2-like of the three compounds, (b) all three compounds either partially or fully antagonized the action of E2 on most genes, with the order of antagonist activity being ICI > Ral > TOT, (c) TOT and Ral, but not ICI, displayed partial agonist/partial antagonist activity on a number of E2-regulated genes, (d) several stimulatory cell cycle-related genes were down-regulated exclusively by ICI, (e) the estrogen-like activity of Ral nearly always overlapped with that of TOT, indicating that Ral has little unique agonist activity different from that of TOT, and (f) some genes were specifically up-regulated by TOT but not Ral, ICI, or E2. Hence, gene expression profiling can discern fundamental differences among SERMs and provides insight into the distinct biologies of TOT, Ral, and ICI in breast cancer.

Therapeutic targeting in the estrogen receptor hormonal pathway

Estrogens work along with genetic changes to promote the development and growth of breast cancers. Because estrogenic hormones act via the estrogen receptors (ERs), ER-alpha and ER-beta, and the ER is present in more than half of breast tumors, this receptor has been the most widely targeted protein in breast cancer therapy. The presence of the ER in breast tumors predicts improved disease-free survival and response to selective ER modulators (SERMs), such as tamoxifen, or other forms of endocrine therapy. Suppression of ER activity by SERMs has proven to be a great benefit in the treatment of breast cancers and also in the prevention of breast cancer in women at high risk for the disease. The Study of Tamoxifen and Raloxifene trial comparing tamoxifen versus raloxifene effectiveness in breast cancer prevention is currently under way. To understand the balance of beneficial and undesirable effects of SERMs and to optimize their effectiveness, current investigations seek to characterize the genes activated or suppressed by these agents. Elucidation of the gene networks and cell signaling pathways under estrogen and SERM regulation and a clearer definition of the respective roles of ER-alpha and ER-beta and their coregulators in the actions of selective ER ligands, should enable the identification of new gene targets for therapeutic intervention and the development of novel drugs for the optimal treatment and prevention of breast cancer.

Molecular determinants of ER alpha and ER beta involved in selectivity of 16 alpha-iodo-7 beta estradiol

The two known estrogen receptors, ER alpha and ER beta, are hormone inducible transcription factors that have distinct roles in regulating cell proliferation and differentiation. The natural ligand, 17 beta-estradiol (E2), binds with high affinity to both ER alpha and ER beta. However, a close analogue, 16 alpha-iodo-17 beta-estradiol (16 alpha IE2) showed about 10-fold selectivity for ER alpha over ER beta. From X-ray studies, it has been shown that the ligand-binding domains (LBD) of the two receptors are strikingly similar, and that only two changes fall within the binding cavity (ER alpha Leu384 to ER beta Met336, and ER alpha Met421 to ER beta Ile373). To understand the molecular basis for the ER alpha selectivity of 16 alpha IE2, mutants and chimeras of ER alpha and ER beta were generated, and ligand-binding and transactivation functions were studied. The ER alpha Leu384 Met mutant behaved like ER alpha WT in the presence of 16 alpha IE2; whereas the profile of the ER alpha Met421 Ile mutant was similar to that of ER beta WT. The ER beta mutant Ile373 Met behaved like ER alpha with 16 alpha IE2. The results clearly demonstrate the role of ER alpha Met421 in the ER alpha selectivity of 16 alpha IE2.

The food contaminants bisphenol A and 4-nonylphenol act as agonists for estrogen receptor alpha in MCF7 breast cancer cells

Xenoestrogens are chemically distinct industrial products potentially able to disrupt the endocrine system by mimicking the action of endogenous steroid hormones. Among such compounds, the ubiquitous environmental contaminants bisphenol A (BPA) and 4-nonylphenol (NPH) may promote adverse effects in humans triggering estrogenic signals in target tissues. Following a research program on human exposure to endocrine disruptors, we found contamination of fresh food by BPA and NPH. More important, these contaminants were found to display estrogen-like activity using as a model system the estrogen-dependent MCF7 breast cancer cells (MCF7wt); its variant named MCF7SH, which is hormone-independent but still ERalpha-positive, and the steroid receptor-negative human cervical carcinoma HeLa cells. In transfection experiments BPA and NPH activated in a direct manner the endogenous ERalpha in MCF7wt and MCF7SH cells, as the antiestrogen hydroxytamoxifen was able to reverse both responses. Moreover, only the hormone-binding domains of ERalpha and ERbeta expressed by chimeric proteins in HeLa cells were sufficient to elicit the transcriptional activity upon BPA and NPH treatments. Transfecting the same cell line with ERalpha mutants, both contaminants triggered an estrogen-like response. These transactivation properties were interestingly supported in MCF7wt cells by the autoregulation of ERalpha which was assessed by RT-PCR for the mRNA evaluation and by immunoblotting and immunocytochemistry for the determination of protein levels. The ability of BPA and NPH to modulate gene expression was further confirmed by the upregulation of an estrogen target gene like pS2. As a biological counterpart, concentrations of xenoestrogens eliciting transcriptional activity were able to stimulate the proliferation of MCF7wt and MCFSH cells. Only NPH at a dose likely too high to be of any physiological relevance induced a severe cytotoxicity in an ERalpha-independent manner as ascertained in HeLa cells. The estrogenic effects of such industrial agents together with an increasing widespread human exposure should be taken into account for the potential influence also on hormone-dependent breast cancer disease.

Cloning and characterization of estrogen receptor alpha in mummichog, Fundulus heteroclitus

Developmental exposure to 17 beta-estradiol (E(2)) induced the death of embryos and fry, malformations, sex reversal, and incomplete ossification of vertebrae and cranial bones in the cyprinodont fish, the mummichog (Fundulus heteroclitus). To clarify the mechanism by which exogenous estrogens caused these developmental effects, we determined the sequence of an estrogen receptor (ER) coding region, encoded by 620 amino acid residues. This region shared 80% identity to that of ER alpha of medaka (Oryzias latipes). Northern blot analysis showed that two ER alpha mRNAs with 5.5 and 4 kb were expressed in the liver. These mRNAs were strongly induced by E(2) stimulation. The 4 kb mRNA was expressed 8 h after treatment, whereas the 5.5 kb mRNA was not induced until 12 h after E(2) stimulation. Vitellogenin (VTG) was expressed 8 h after E(2) stimulation in the male liver. Receptor binding assays using the protein of F. heteroclitus ER alpha (fh ER alpha) ligand binding domain showed that alkylphenols bind to fh ER alpha with a higher affinity (50 times or more) as compared with the human ER alpha. The present results demonstrate that the fh ER alpha has a sequence very similar to that of medaka, and the mRNA for this receptor was induced by E(2)-stimulation, followed subsequently by VTG expression. Furthermore, alkylphenols bind to fh ER alpha more efficiently than to human ER alpha.

Molecular basis for the subtype discrimination of the estrogen receptor-beta-selective ligand, diarylpropionitrile

Although the two subtypes of the human estrogen receptor (ER), ERalpha and ERbeta, share only 56% amino acid sequence identity in their ligand binding domain (LBD), the residues that surround the ligand are nearly identical; nevertheless, subtype-selective ligands are known. To understand the molecular basis by which diarylpropionitrile (DPN), an ERbeta-selective ligand, is able to discriminate between the two ERs, we examined its activity on ER mutants and chimeric constructs generated by DNA shuffling. The N-terminal region of the ERbeta LBD (through helix 6) appears to be fully responsible for the ERbeta selectivity of DPN. In fact, a single ERalpha point mutation (L384M) was largely sufficient to switch the DPN response of this ER to that of the ERbeta type, but residues in helix 3 are also important in achieving the full ERbeta selectivity of DPN. Using molecular modeling, we found an energetically favorable fit for the S-DPN enantiomer in ERbeta, in which the proximal phenol mimics the A ring of estradiol, and the nitrile engages in stabilizing interactions with residues in the ligand-binding pocket of ERbeta. Our findings highlight that a limited number of critical interactions of DPN with the ERbeta ligand-binding pocket underlie its ER subtype-selective character.

Transcriptional activation of estrogen receptor ERalpha and ERbeta by polycyclic musks is cell type dependent

In the past decade the list of chemicals in the environment that are able to mimic the natural hormone estrogen, thereby disrupting endocrine function, has grown rapidly. These chemicals are able to bind to estrogen receptors (ERs) and influence estrogen signaling pathways, although several of them have structures that differ substantially from the endogenous hormone 17beta-estradiol. Because of their polycyclic nature, the polycyclic musks AHTN (6-acetyl-1,1,2,4,4,7-hexamethyltetraline) and HHCB (1,2,4,6,7,8-hexahydro-4,6,6,7,8,8-hexamethylcyclopenta-gamma-2-benzopyran) were assessed for interaction with estrogen receptors. These compounds are ubiquitously present in surface waters and have been identified in human milk fat and blood. Using sensitive in vitro reporter gene assays, we found that AHTN and HHCB act as selective estrogen receptor modulators (SERMs), inducing both estrogenic and antiestrogenic activity dependent on the cell line and the ER subtype targeted. Weak estrogenic effects were observed only at relatively high concentrations (10 microM). Antiestrogenic effects were observed in various cell lines starting at concentrations of 0.1 microM. In comparison with the well-known SERM, 4-hydroxytamoxifen, AHTN and HHCB have a much lower potency in suppressing estradiol-induced transactivation. Results with various mutant ER receptor types suggest that binding sites of the musk compounds differ from 17beta-estradiol and 4-hydroxytamoxifen. The cell type dependency of ER transactivation shows that caution should be exerted when interpreting effects of estrogenic compounds using in vitro systems.

Functionally orthogonal ligand-receptor pairs for the selective regulation of gene expression generated by manipulation of charged residues at the ligand-receptor interface of ER alpha and ER beta

The reengineering of protein-small molecule interfaces represents a powerful tool of chemical biology. For many applications it is necessary to engineer receptors so that they do not interact with their endogenous ligands but are highly responsive to designed ligand analogues, which in turn do not interact with endogenous proteins. The chemical design strategy used to reengineer protein-small molecule interfaces is particularly challenging for interfaces involving relatively plastic receptor binding sites and therefore presents a unique challenge in molecular design. In this study we explore the scope and limitations of a new strategy for manipulating polar/charged residues across the ligand receptor interface of estradiol (E2) and the estrogen receptor (ER). Carboxylate-functionalized E2 analogues can activate ER alpha(Glu353-->Ala) and ER beta(Glu305-->Ala) with very large selectivites, demonstrating that this design strategy is extendable to other members of the steroid hormone receptor family. Neutral E2 analogues were found to complement ER alpha(E353A) with similar potencies but with generally lower selectivities. This suggests that the high selectivity observed with ligand-receptor pairs generated by exchanging charged residues across ligand-receptor interfaces is only due in part to their complementary shapes and that appropriate introduction of charged functionality on the ligand can provide substantial enhancement of selectivity by decreasing the engineered ligands affinity for the endogenous receptor. Attempts to modify the cationic residues by complementing Arg394-->Ala or Arg394-->Glu were not successful.

Alteration of large-scale chromatin structure by estrogen receptor

The estrogen receptor (ER), a member of the nuclear hormone receptor superfamily important in human physiology and disease, recruits coactivators which modify local chromatin structure. Here we describe effects of ER on large-scale chromatin structure as visualized in live cells. We targeted ER to gene-amplified chromosome arms containing large numbers of lac operator sites either directly, through a lac repressor-ER fusion protein (lac rep-ER), or indirectly, by fusing lac repressor with the ER interaction domain of the coactivator steroid receptor coactivator 1. Significant decondensation of large-scale chromatin structure, comparable to that produced by the approximately 150-fold-stronger viral protein 16 (VP16) transcriptional activator, was produced by ER in the absence of estradiol using both approaches. Addition of estradiol induced a partial reversal of this unfolding by green fluorescent protein-lac rep-ER but not by wild-type ER recruited by a lac repressor-SRC570-780 fusion protein. The chromatin decondensation activity did not require transcriptional activation by ER nor did it require ligand-induced coactivator interactions, and unfolding did not correlate with histone hyperacetylation. Ligand-induced coactivator interactions with helix 12 of ER were necessary for the partial refolding of chromatin in response to estradiol using the lac rep-ER tethering system. This work demonstrates that when tethered or recruited to DNA, ER possesses a novel large-scale chromatin unfolding activity.

Ligand-specific regulation of proteasome-mediated proteolysis of estrogen receptor-alpha

Proteasome-mediated proteolysis modulates the cellular concentration of estrogen receptor-alpha (ERalpha) and is induced by treatment of cells with 17beta-estradiol. Herein, we show that multiple receptor agonists, including 17alpha-estradiol and estriol as well as the antagonist ICI-182780, stimulate proteasome-dependent proteolysis of ERalpha in a process that requires ligand binding to the receptor. Proteolysis of receptor depends on ligand concentration, and there exists a direct correlation between ligand-binding affinity and the half-maximal dose of ligand required to stimulate receptor degradation. Furthermore, introduction of a point mutation into the receptor ligand-binding pocket yields a stable receptor resistant to proteolysis. Interestingly, although all ligands stimulate receptor degradation, the extent to which overall ER levels are affected varies with each ligand and is not related to ligand-binding affinity or activation of transcription. These results demonstrate ligand-specific regulation of ERalpha proteolysis, and they introduce the concept that cellular receptor concentration is governed not only at the level of induction of proteolysis but also by the efficiency with which the receptor is degraded.

Structure-function relationships in estrogen receptors and the characterization of novel selective estrogen receptor modulators with unique pharmacological profiles

This article summarizes recent research on the development of estrogen receptor alpha (ER alpha) and estrogen receptor beta (ER beta) subtype-selective ligands based on our understanding of structure-activity relationships in these two estrogen receptors and differences in their ligand binding domains and activation function domains. The use of these ligands should enable greater understanding of the unique biologies mediated by ER alpha versus ER beta and may, as well, provide selective estrogen receptor modulators having unique biological and pharmacological profiles optimal for prevention and treatment of breast cancer, for menopausal hormone replacement, for prevention of osteoporosis, and for potential cardiovascular benefit.