![]()
|
|
||||||||
J. Biol. Chem., Vol. 281, Issue 26, 17909-17919, June 30, 2006
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

¶1




3
From the
Ben May Institute for Cancer Research, the
Department of Biochemistry and Molecular Biology, and the ¶Howard Hughes Medical Institute and Department of Chemistry, University of Chicago, Chicago, Illinois 60637
Received for publication, December 22, 2005 , and in revised form, April 28, 2006.
| ABSTRACT |
|---|
|
|
|---|
(ER
) and
(ER
) have distinct functions and differential expression in certain tissues. These differences have stimulated the search for subtype-selective ligands. Therapeutically, such ligands offer the potential to target specific tissues or pathways regulated by one receptor subtype without affecting the other. As reagents, they can be utilized to probe the physiological functions of the ER subtypes to provide information complementary to that obtained from knock-out animals. A fluorescence resonance energy transfer-based assay was used to screen a 10,000-compound chemical library for ER agonists. From the screen, we identified a family of ER
-selective agonists whose members contain bulky oxabicyclic scaffolds in place of the planar scaffolds common to most ER ligands. These agonists are 1050-fold selective for ER
in competitive binding assays and up to 60-fold selective in transactivation assays. The weak uterotrophic activity of these ligands in immature rats and their ability to stimulate expression of an ER
regulated gene in human U2OS osteosarcoma cells provides more physiological evidence of their ER
-selective nature. To provide insight into the molecular mechanisms of their activity and selectivity, we determined the crystal structures of the ER
ligand-binding domain (LBD) and a peptide from the glucocorticoid receptor-interacting protein 1 (GRIP1) coactivator complexed with the ligands OBCP-3M, OBCP-2M, and OBCP-1M. These structures illustrate how the bicyclic scaffolds of these ligands are accommodated in the flexible ligand-binding pocket of ER. A comparison of these structures with existing ER structures suggests that the ER
selectivity of OBCP ligands can be attributed to a combination of their interactions with Met-336 in ER
and Met-421 in ER
. These bicyclic ligands show promise as lead compounds that can target ER
. In addition, our understanding of the molecular determinants of their subtype selectivity provides a useful starting point for developing other ER modulators belonging to this relatively new structural class. | INTRODUCTION |
|---|
|
|
|---|
(ER
)4 and
(ER
) are ligand-inducible transcription factors that are involved in regulating cell growth, proliferation, and differentiation in various normal and cancerous tissues (13). Although the two subtypes of ER bind the endogenous estrogen, 17
-estradiol (E2), with similar affinity (4), they differ in size, share modest sequence identity (47%), and are encoded by different genes (5, 6). Studies of knock-out mice have also shown that the two subtypes have distinct functions and are differentially expressed in certain tissues (1). These differences have stimulated the search for subtype-specific ligands that can elicit tissue- or cell-specific ER activity. In particular, the dominance of ER
expression in the breast and uterus (7) suggests that ER
-selective ligands may offer some of the benefits of hormone replacement therapy such as a decrease in the risk of colorectal cancer (8) without increasing the risk of breast or uterine cancer. The recent discovery of ER
-selective ligands that display pathway-specific anti-inflammatory activity without classic estrogenic effects is an example of the possible therapeutic potential of subtype-selective ligands (9). Furthermore, subtype-selective ligands show promise as reagents to probe the physiological functions of ER
and ER
(1012), providing complementary information to studies in knock-out mice (1).
Progress towards the development of subtype-selective ligands was significantly advanced with the reports of crystal structures of ER
(13) and ER
(14). The ligand-binding pockets of the subtypes are similar but not identical. The ER
ligand-binding pocket is smaller (390 Å3 versus 490 Å3 for ER
) and differs in two residues from ER
; Leu-384 and Met-421 in ER
are replaced by Met-336 and Ile-373, respectively, in ER
(14). Notably, these two substitutions give rise to the selectivity of ligands for ER
or ER
. Although several ER
-selective ligands have been described (15), only a small number of crystal structures for these ligands complexed with ER
(14, 1620) and ER
(17, 20, 21) have been reported. These structures have provided valuable insights into possible molecular mechanisms of ER
selectivity. However, when one considers the structural diversity of the known ER
-selective ligands and the plasticity of the ER ligand-binding pocket, it is clear that more structures of ER subtype-selective ligands are needed to develop general mechanisms for subtype selectivity.
We describe here the identification and characterization of a family of ligands sharing an oxabicyclic scaffold that are ER agonists and variably selective for ER
. These ligands possess bulky bicyclic scaffolds that differ from the planar scaffolds common to most estrogenic compounds. Cell and animal studies demonstrate the promise of these ligands as agents that can target ER
for therapeutic or investigational purposes. In addition, we obtained crystal structures with three of these ligands, OBCP-3M, OBCP-2M, and OBCP-1M, in complex with the ER
ligand-binding domain (LBD) and a glucocorticoid receptor-interacting protein 1 (GRIP1) peptide. These structures provide evidence to suggest that the ER
selectivity of the OBCP ligands can be attributed to their interactions with only two residues: Met-336 in ER
and Met-421 in ER
.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
PlasmidsMany of the plasmids used in these studies were generous gifts from colleagues: pGST-ER
-LBD (ER
residues 282595) from Dr. Peter Kushner (University of California, San Francisco, CA), pCR3.1.SRC1 from Dr. Ming-Jer Tsai and Dr. Sophia Y. Tsai (Baylor College of Medicine, Houston, TX), p3xERE-Luc from Dr. Donald McDonnell (Duke University, Durham, NC), pCMV5-hER
(full-length ER
) and p6xHis-ER
LBD (residues 256505) from Dr. Benita Katzenellenbogen (University of Illinois, Champaign, IL), HEGO/pSG5 (full-length ER
) from Dr. Pierre Chambon (Institute of Genetics and Molecular and Cellular Biology, Strasbourg, France), pAR and p3xARE-Luc from Dr. Shutsung Liao (University of Chicago, Chicago, IL), and pMCSG7 from Dr. Frank Collart (Argonne National Laboratory, Argonne, IL). pCMV-
gal was obtained from Invitrogen.
The pFLAG-SRC1-His was constructed by cloning the steroid receptor coactivator-1 (SRC1) fragment (residues 568780) from pCR3.1.SRC1 into pFLAG-CMV2 (Sigma) and subsequently subcloning into pET28b (Novagen). pFLAG-ER
and pFLAG-ER
used in transient transfections were constructed by subcloning full-length hER
(from HEGO/pSG5) or hER
(from pCMV5-hER
) into pCMV-FLAG-2. p6xHis-ER
LBD used in the competitive binding assays was constructed by subcloning ER
LBD (residues 304554) into pET15b (Novagen). The ER
LBD (residues 298554; Y537S) used for crystallization was expressed from pMCSG7 (25) as a fusion protein with a 23-residue N-terminal tag. The Y537S mutation was introduced to increase the solubility and yield of the ER
LBD for crystallization. From comparisons of the crystal structures of this mutant ER with wild type ER complexed to the same ligand, the mutation does not appear to cause significant changes to the overall conformation of the protein.5
Protein Expression and Purification for in Vitro AssaysGlutathione S-transferase (GST)-ER
LBD was expressed in Escherichia coli Rosetta cells (Novagen) by induction with 0.2 mM isopropyl-
-D-thiogalactopyranoside for 3 h at 25°C.The protocol for the purification of the GST-ER
LBD was similar to that which has been described elsewhere (26). FLAG-SRC1, ER
LBD, and ER
LBD were all expressed and purified as generally described. Proteins were expressed in Rosetta cells by induction with 0.6 mM isopropyl-
-D-thiogalactopyranoside for 3 h at 25°C. Cell pellets were flash frozen in liquid nitrogen and stored at 80 °C. After thawing, they were resuspended in 5 volumes of lysis buffer (20 mM Tris buffer, pH 7.6, 500 mM NaCl, 10% (v/v) glycerol, 0.05%
-octyl glucoside, 10 mM imidazole, 10 mM
-mercaptoethanol (
-ME), 3 M urea (for ER LBD), protein inhibitor mixture (Calbiochem), 0.1 mg/ml lysozyme) and sonicated on ice three times for 30 s each. The lysates were centrifuged at 30,000 x g for 30 min to collect the supernatant, and the above procedure was repeated to extract additional soluble proteins from the cell debris. The two supernatant fractions were pooled and loaded on Ni2+-nitrilotriacetic acid resin (Qiagen) pre-equilibrated with lysis buffer. The resin was washed with 5 column volumes of wash buffer (50 mM Tris buffer, pH 7.6, 500 mM NaCl, 0.05%
-octyl glucoside, 20 mM imidazole, 10 mM
-ME, protein inhibitor mixture) and eluted with 10 column volumes of elution buffer (50 mM Tris, pH 7.6, 300 mM NaCl, 0.05%
-octyl glucoside, 250 mM imidazole, 4 mM tris(2-carboxyethyl)phosphine hydrochloride, protein inhibitor mixture). Collected fractions were analyzed by SDS-PAGE before selecting the purest protein fractions to concentrate. Additionally, ER
LBD and ER
LBD were dialyzed in binding buffer (50 mM Tris pH 8.0, 100 mM NaCl, 0.01%
-octyl glucoside, and 0.1% (v/v) glycerol) and FLAG-SRC1 in HTRF buffer (50 mM HEPES, pH 7.4, 50 mM KCl, 0.1% (v/v) Tween 20, and 1 mM EDTA).
Homogenous Time-resolved Fluorescence (HTRF) AssaysThe anti-FLAG-XL665 and anti-GST-Eu(K) antibodies were from Cisbio International (Gif/Yvette Cedex, France) and the Wallac 1420 MicroPlate Reader from PerkinElmer Life Sciences. FLAG-SRC1 and GST-ER
LBD were diluted in HTRF buffer that included 1 mM dithiothreitol and 1 mg/ml bovine serum albumin and preincubated with anti-FLAG-XL665 and anti-GST-Eu(K), respectively, for 1 h at 4°C.Thetwomixtures were then combined with test ligands (0.5%, v/v) and 2 M KF solution and incubated at 4 °C for 12 h with agitation. The plates were read on a Wallac 1420 MicroPlate Reader. The extent of the specific fluorescence resonance energy transfer was determined by taking a ratio of the emission intensities at 665 and 615 nm and scaled by 104. EC50 values in this and other assays were determined from the data using SigmaPlot (Systat Software).
Competitive Binding AssaysHis-tagged ER
LBD/ER
LBD was diluted with binding buffer to 810 nM and added to the wells of nickel-coated FlashPlates (PerkinElmer Life Sciences). The proteins were incubated for
2 h at 25 °C to allow the ER to bind. Plates were then washed five times with phosphate-buffered saline, before the addition of 12 nM [3H]E2 and test ligand diluted in binding buffer. The ligands were incubated for 4 h at 25 °Cor overnight at 4 °C, and the plates were analyzed with a MicroBeta Scintillation Counter (PerkinElmer Life Sciences). IC50 and Ki values were determined from the data using SigmaPlot (Systat Software).
Transient Transfection and Transactivation AssaysIn the transient transfection of COS-7 and SKBR3 cells, DNA was delivered to the cells using PolyFect (Qiagen). In the case of SKBR3, 150 ng of p3xERE-Luc, 30 ng of pFLAG-ER
/pFLAG-ER
, and 30 ng of normalization plasmid (pCMV-
gal) were used. For the COS-7 ER assays, 300 ng of p3xERE-Luc, 40 ng of pFLAG-ER
/pFLAG-ER
, and 40 ng of pCMV-
gal were used. In the COS-7 androgen receptor (AR) assays, 150 ng of p3xARE-Luc, 30 ng of pAR, and 30 ng of pCMV-
gal were used. Ligands that had been diluted 1:200 in phenol red-free Dulbecco's modified Eagle's medium supplemented with 10% charcoal-stripped fetal bovine serum were added to the cells 1824 h before cell harvest and lysis. Luciferase activities from lysed cells were analyzed on a MicroBeta (PerkinElmer Life Sciences) luminescence counter, and values were normalized with those from
-galactosidase activity (Promega).
|
Quantitative Real Time Reverse Transcription-PCR AssaysU2OS-ER
and U2OS-ER
cells (generous gifts from Drs. Thomas Spelsberg and David Monroe (Mayo Clinic, Rochester, MN)) were seeded in 6-well plates at 46 x 104 cells/well and pretreated for 24 h with 100 ng/ml doxycycline (Sigma). Cells were subsequently treated for an additional 24 h with ligands or vehicle (Me2SO) control in the presence of 100 ng/ml doxycycline to maintain ER expression. Each treatment group was performed in triplicate. Total RNA was prepared with TRIzol reagent (Invitrogen). 2 µg of RNA was treated with DNase I (Invitrogen) before being reverse transcribed using the Superscript III First-Strand Synthesis System (Invitrogen). A combination of oligo(dT) and random hexamers was used to prime the cDNA synthesis reaction. The resultant cDNA products were diluted to 40 µl, and 5 µl of cDNA was used in real time PCRs utilizing the QuantiTect SYBR Green PCR kit (Qiagen). The primer sequences used for real time PCR for angiotensinogen and
-actin have been previously reported by others (28). The reactions were carried out using the ABI 7300 Real-Time PCR System (Applied Biosystems) for 48 cycles (94 °C for 15 s, 5458 °C for 30 s, 72 °C for 40 s) after an initial 15-min incubation at 95 °C. RNA levels were determined for angiotensinogen and
-actin RNA by comparison with standard curves generated from reference RNA (Stratagene). Angiotensinogen expression was then normalized to the endogenous reference gene
-actin, and the relative expression was then determined by normalizing to the Me2SO-treated control.
Expression, Purification, and Crystallization of ER
LBDhER
LBD was expressed in BL21(DE3) magic strain E. coli (29). Following isopropyl-
-D-thiogalactopyranoside induction, the cell pellets in 0.1 M HEPES, pH 7.5, 0.5 M NaCl, 10 mM
-ME, 5% (v/v) glycerol, and 10 mM imidazole were incubated with lysozyme (1 mg/ml) and protein inhibitor mixture (Sigma) and then sonicated. The ER protein was purified on aNi2+-nitrilotriacetic acid column (Qiagen), and the N-terminal polyhistidine tag was cleaved using tobacco etch virus protease (30), which carries a noncleavable polyhistidine tag. Following its separation from the protease and the tag on the Ni2+-nitrilotriacetic acid column, the cleaved protein was diluted 1:10 in 25 mM Tris, pH 8.0, 0.1 M NaCl, 10 mM
-ME, and 5% (v/v) glycerol and purified further on a mono-Q anion exchange column (GE Healthcare) using an NaCl gradient of 25 mM to 0.6 M. The protein eluted at 0.3 M NaCl and was dialyzed in 4 mM HEPES, pH 7.5, 0.125 M NaCl, 10 mM
-ME, and 1.25% (v/v) glycerol. The protein was concentrated to
10 mg/ml, incubated overnight with a 1 mM concentration of both the ligand and the GRIP1 peptide (residues 686698) prior to screening for crystallization conditions in sitting drops (2:1 protein/reservoir) with the commercially available screens Index (Hampton Research) and Wizard (Emerald Biosystems) at both 16 and 25 °C. Optimization of initial screen conditions yielded diffraction quality crystals of ER in
1 week at 16 °C in 0.2 M sodium malonate, 20% (w/v) polyethylene glycol 3350, pH 7.0, for OBCP-3M and in 0.2 M ammonium sulfate, 0.1 M Tris-HCl, 25% (w/v) polyethylene glycol 3350, pH 8.5, for OBCP-2M and OBCP-1M. Crystals were soaked in the mother liquor and 10% (v/v) glycerol prior to freezing in liquid nitrogen.
Data Collection, Structure Determination, and RefinementSingle wavelength (0.979 Å) native data sets were collected on a MarCCD detector at the 19-BM beamline (OBCP-3M) and on a Quantum 135 detector at the 14-BM beamline (OBCP-2M and OBCP-1M) of the Advanced Photon Source (Argonne National Laboratory). Data were collected with an oscillation of 1.5, 1.0, and 0.5° and an exposure of 10, 12, and 3 s per image, to a resolution of 2.25, 2.10, and 1.80 Å with average redundancies of 3.6, 5.1, and 3.6 for OBCP-3M, OBCP-2M, and OBCP-1M structures, respectively. Data were integrated and merged using HKL2000 (31) and then used in molecular replacement using MOLREP (32) with the ER structure from Protein Data Bank entry 1L2I as the search model for the OBCP-3M structure and 1ZKY as the search model for both the OBCP-2M and OBCP-1M structures. XtalView (33), CNS (version 1.1) (34), and REFMAC (35) programs were used for model building and refinement. Missing residues and alternate conformations for some side chains were manually modeled and refined in subsequent cycles. The missing residues in the final models are at the termini and surface loops. Details of diffraction data collection and processing of the ER crystal complexes are shown in Table 5. The ribbon diagram was prepared with Swiss-PdbViewer (36) and rendered in POV-RAY (Persistence of Vision Ray Tracer; available on the World Wide Web at www.povray.org). Figures showing electron density maps were prepared with BOBSCRIPT (37) and rendered in Raster3D (38). Structural alignments were performed with Swiss-PdbViewer and presented using MOLSCRIPT (39) and Raster3D.
|
| RESULTS |
|---|
|
|
|---|
LBD and SRC1 NRD to Identify ER AgonistsSRC1 nuclear receptor interaction domain (SRC1 NRD) is a protein fragment encompassing all three nuclear receptor box motifs (40) known to interact with the LBDs of nuclear receptors (41). To find ER agonists, compounds that promote the interaction between the ER
LBD and SRC1 NRD were identified using HTRF (42), a technology based on fluorescence resonance energy transfer. Of the 10,000 compounds screened in mixtures of eight compounds each (to give working concentrations of
20 µM/compound), 40 mixtures increased HTRF signals to greater than 40% of the signal obtained with 10 nM E2 alone. We repeated the assay with each of the 320 compounds present in the original mixtures to select the 10 most active compounds. These compounds were tested over a range of concentrations to confirm dose-dependent promotion of ER
LBD/SRC1 NRD interactions. Compound 5474 (Fig. 1) was identified as the lead compound for structure activity relationship (SAR) studies.
SAR Studies on the 5474 Compound FamilyBased on the core scaffold of 5474, we purchased or synthesized a variety of chemical analogs (Fig. 1) to evaluate the contributions of different functional groups to their activity, as assessed by the HTRF assay. Concentrations ranging from 0.1 nM to 100 µM were used to obtain EC50 values for each compound (Table 1). EC50 values reflect the potency of each compound in promoting the ER
LBD interaction with SRC1 NRD. A comparison of compounds 5474 and 5474E showed that removal of the hydroxymethyl group attached to the oxabicyclic scaffold resulted in a loss of activity. However, a greater impact occurred when the phenolic OH group was moved from the ortho position (5474) to the para position (OBCP-3M), which led to a nearly 300-fold increase in potency. Modifications of the phenolic ring of OBCP-3M, including removal of the OH group (5474A), methylation of the OH group (5474C), or the addition of a methoxy group to the ring (5474H), all resulted in a loss of activity. In contrast, the addition or subtraction of methyl groups from the oxabicyclic scaffold had a smaller impact on the EC50 values. Successive removal of methyl groups from the oxabicyclic scaffold of OBCP-3M resulted in progressively less activity.
|
SelectivityThe three most potent ligands from the SAR studies, OBCP-3M, OBCP-2M, and OBCP-1M (OBCP compounds) were evaluated in a competitive binding assay using [3H]E2 to determine their affinities for the hER
and hER
LBDs. Each compound was able to displace [3H]E2 from recombinant ER
and ER
LBD, suggesting that these compounds bind in the ligand-binding pocket of the receptors. The Ki values of the compounds for ER
ranged from 560 to 1800 nM, and those for ER
ranged from 12 to 190 nM (Table 2). All three compounds showed selectivity for ER
, with preferences that ranged from
10- to 50-fold (Table 2). These affinities and selectivities for ER
are on the same order of magnitude as those for genistein, which binds ER
and ER
with affinities of 530 and 7 nM, respectively, and is 76-fold ER
-selective (Table 2).
|
and ER
in Mammalian Cell Lines and Retain ER
SelectivityThe OBCP compounds were also tested in transient transfection assays for their ability to activate ER transcription from a luciferase reporter under the control of a promoter containing three estrogen response elements (3xERE). In COS-7 (monkey kidney) cells, all three compounds tested were agonists on full-length ER
and ER
(Fig. 2, A and B). The selectivity of the OBCP compounds for ER
activation, determined from their transcriptional potencies on each receptor, ranged from
30- to 60-fold (Table 3). Efficacy on ER
transcription was comparable with E2, whereas for ER
, transcription efficacy was 6070% of the levels achieved by E2. To test for tissue-dependent effects, transient transfections using the same plasmids were carried out in SKBR3 (human breast cancer) cells with OBCP-3M and genistein. In this context, the ER
selectivity of OBCP-3M (16-fold) was less than in the COS-7 model (29-fold) described above, and the efficacy profile of OBCP-3M also changed. In SKBR3 cells, OBCP-3M activated ER
and ER
to 70% of the levels achieved by E2 (Fig. 2, C and D). Finally, in comparison with genistein, the ER
selectivity of OBCP-3M was also lower (16- versus 48-fold) (Table 4).
|
|
|
OBCP-3M and OBCP-1M Have Very Weak Uterotrophic Activity in Immature Sprague-Dawley RatsTo test OBCP compounds in vivo, we used a rat uterotrophic assay (27), which measures estrogenic activity based on the increase in uterine weight in immature rats (43). In contrast to E2, which at a dose of 1 nmol/animal/day increased uterine weight nearly 3-fold relative to vehicle control, both OBCP-3M and OBCP-1M were very weakly uterotrophic (
1.2- to 1.5-fold weight increase). These compounds stimulated significant uterine growth only at doses of 50 nmol/animal/day or greater (Fig. 3).
OBCP-3M Stimulates Expression of a Previously Identified Endogenous ER
Target GeneTo test the effect of OBCP compounds on an endogenous ER
gene target, we utilized two stably transfected U2OS cell lines containing doxycycline-inducible ER
(U2OS-ER
) or ER
(U2OS-ER
). These cell lines (28) and other similar U2OS cell lines (44, 45) have been previously used to identify genes regulated by specific ER subtypes. In these studies, angiotensinogen was shown to be a gene specifically up-regulated by E2 stimulation in the U2OS-ER
cells but not the U2OS-ER
cells (28). We monitored its expression to determine whether OBCP-3M had an effect on this ER
gene target. As controls, a non-subtype-selective agonist (E2), an ER
-selective agonist (PPT) (46), and an ER
-selective agonist (DPN) (47) were also tested on both cell lines. The expression of angiotensinogen in U2OS-ER
cells was unaffected by any of the treatments relative to the effects of the vehicle treatment (Fig. 4A), consistent with its identification as an ER
gene target. However, treatment of cells with E2, DPN, and OBCP-3M stimulated gene expression of angiotensinogen in the U2OS-ER
cells to a higher level than that found in vehicle treated cells (Fig. 4B). Unexpectedly, treatment with PPT caused a decrease in the expression of angiotensinogen in the U2OS-ER
cells relative to vehicle-treated cells.
Structure Determination of ER
LBD Complexes with GRIP1 Peptide and the OBCP-3M, OBCP-2M, and OBCP-1M LigandsThe crystal structures of the OBCP-3M-ER
LBD-GRIP1, OBCP-2M-ER
LBD-GRIP1, and OBCP-1M-ER
LBD-GRIP1 complexes were determined to identify conformational changes that could explain differences between the ER
binding selectivity of OBCP-3M (47-fold), OBCP-2M (29-fold), and OBCP-1M (9.5-fold). Details of x-ray diffraction data collection, processing, and refinement of all three complexes are presented in Table 5. These high resolution ER
LBD structures have overall conformations similar to those found in structures of ER
complexed with the agonists E2 and diethylstilbestrol (DES) (13, 48). Thus, helix 12 packs against helices 3, 5, 6, and 11 to form a hydrophobic groove that accommodates the binding of the
-helical GRIP1 coactivator peptide (Fig. 5). The OBCP-3M, OBCP-2M, and OBCP-1M complexes are very similar to one another and can be superimposed via their backbone atoms with root mean square deviations of <0.50 Å.
|
structures (Protein Data Bank code 1ERE
[PDB]
) are preserved in the OBCP-3M, OBCP-2M, and OBCP-1M structures. For example, the phenolic OH group in all three OBCP ligands participates in a hydrogen bond network that includes Glu-353, Arg-394, and an ordered water molecule. On the opposite side of the pocket, the OBCP hydroxymethyl group forms a hydrogen bond with His-524 (Fig. 6). The bulk of the oxabicyclic scaffolds of the OBCP ligands lies almost entirely above and below the areas that would be occupied by the C- and D-rings of E2 (Fig. 7B). The oxabicyclic scaffold is involved in extensive hydrophobic interactions with residues that line the ligand-binding pocket.
Different Diastereomers of OBCP-3M and OBCP-2M Bind to Each Monomer in the ER
LBD DimerThe OBCP-3M, OBCP-2M, and OBCP-1M compounds used in our studies are all mixtures of diastereomers. Interestingly, in the OBCP-3M and OBCP-2M complexes, two different diastereomers of the respective ligands are bound to each monomer of the ER
LBD dimer. In the OBCP-3M structure, the two diastereomers differ at the C-2 chiral center of the oxabicyclic scaffold, where the phenol ring is attached. Monomer A (MonA) of the ER
LBD dimer contains the C-2R diastereomer, and monomer B (MonB) contains the C-2S diastereomer. When MonA and MonB are superimposed, the two diastereomers can be seen to retain the same hydrogen bond contacts, despite different orientations of the OBCP-3M oxabicyclic scaffolds (Fig. 7A). This can occur because the phenolic OH and the hydroxymethyl OH of the two diastereomers are similarly positioned within each monomer. The major difference between the two monomers is that in MonA, the 6-methyl of OBCP-3M displaces the side chain of Met-343, which subsequently displaces that of Met-342. However, this change only locally perturbs the ER
LBD structure and does not affect the overall structure of the dimer. In contrast to findings in the OBCP-3M structure, the OBCP-2M oxabicyclic scaffold is bound to each monomer in a similar orientation. However, the electron density of the ligand in MonB appears to be produced by a combination of two OBCP-2M diastereomers that differ in the attachment of a methyl group at the C-9 chiral center. Thus, whereas MonA contains the C-9R diastereomer, Mon B is occupied by a mixture of the C-9S (Fig. 6C) and C-9R diastereomers. Interestingly, in the OBCP-1M complex, the same stereoisomer is bound to each monomer.
|
| DISCUSSION |
|---|
|
|
|---|
-selective agonists that have a nonplanar oxabicyclic scaffold. These compounds bind in the ligand-binding pocket of ER and induce an agonist conformation of the receptor (Fig. 5) that permits the recruitment of coactivators such as SRC1 or GRIP1 via their nuclear receptor box motifs. The most selective and potent of these compounds, OBCP-3M, exhibited a
50-fold selectivity for ER
in binding assays.
Cell-based experiments using full-length ERs show that OBCP-3M and its analogs retain their ER
-selective activity and are able to activate ER transcription through a promoter containing the consensus 3xERE motif. This activity is specific to ER, since no activity is seen with AR. In addition, OBCP-3M exhibits some degree of tissue-dependent activity. In COS-7 cells, OBCP-3M acts as a full agonist on ER
and as a partial agonist on ER
(Fig. 2, A and B) with a
30-fold ER
selectivity (Table 3). However, in SKBR3 cells, OBCP-3M is a partial agonist on both ER
and ER
(Fig. 2, C and D) and shows only 16-fold ER
selectivity (Table 4). The differences in the selectivity and efficacy of OBCP-3M in these two cell contexts may be due to the different tissue and species origins of the cell lines. Thus, differing compositions of coactivators in each cell line (49) and possible differences in coactivator recruitment by ER (50) are likely factors affecting the observed activity of OBCP-3M. Future studies to identify coactivators that preferentially interact with ER when OBCP-3M or its analogs are bound may help optimize these compounds for cell- or tissue-selective activity.
|
|
activity (10). The low level of uterotrophic activity exhibited at higher ligand doses is probably due to their residual ER
activity, an effect also seen with other ER
-selective ligands (11). However, our results do not preclude the possibility that the weaker activity compared with E2 can be due to differences in the absorption, distribution, metabolism, or excretion of the compounds.
Experiments in U2OS-ER
and U2OS-ER
cells provide further support for our in vivo findings and suggest that OBCP-3M can stimulate ER
activity. Treatment of U2OS-ER
cells with OBCP-3M showed that the compound can activate ER
to up-regulate the expression of angiotensinogen, a previously identified ER
target gene (28). Additional studies on the expression of other endogenous ER
target genes would better define the ER
-selective properties of the OBCP compounds.
SAR studies examining the different analogs within the 5474 compound family indicate that two functional groups on the OBCP ligands are important for binding to the ER
LBD: the phenol ring and the hydroxymethyl group at the C-2 and C-5 positions, respectively, of the oxabicyclic scaffold. Whereas the former is critical for strong interactions of the OBCP compounds with the ER
LBD, the hydroxymethyl group enhances the binding of these ligands. Based on the mechanisms by which E2 binds the ER
LBD (13), we predicted that the para phenolic OH group of the OBCP ligands probably mimics the binding mode of the E2 A-ring 3-OH, whereas the OH of the hydroxymethyl group binds in a similar mode as the 17
-OH of E2. Crystal complexes of ER
LBD-GRIP1 with OBCP-3M, OBCP-2M, and OBCP-1M confirmed these predictions (Figs. 5 and 6B). These structures also show that the oxabicyclic scaffolds of the OBCP ligands occupy some of the vacant space that would surround the C- and D-ring positions of E2 (13) in the ligand binding pocket (Fig. 7B).
|
LBD dimer. Another nuclear receptor, the pregnane X receptor, was shown previously to bind the compound SR12813 in a number of different orientations (51), but to our knowledge, this is the first report of a nuclear receptor binding two different stereoisomers in a dimer. Surprisingly, the OBCP-1M complex, which also crystallized in the presence of a mix of diastereomers and packed in the same space group with similar cell dimensions as the OBCP-3M and OBCP-2M complexes, bound the same stereoisomer in each monomer. In the OBCP-2M structure, the two diastereomers found in MonB bind in the same general orientation and differ only in the position of the methyl attached to the C-9 chiral center. In contrast, in the OBCP-3M structure, each diastereomer binds its monomer in a very different orientation (Fig. 6, A and B). In this case, the unique binding of OBCP-3M is probably due to the symmetry conferred on its six-membered ring (C1-C9-C5-C6-C7-C8) by the three methyl groups. This symmetry, which is lacking in OBCP-2M and OBCP-1M, allows the substitution of methyl groups in making ligand-protein contacts when different diastereomers bind the receptor. Comparative structural analysis of the two monomers (Fig. 7A) shows that the receptor contacts made by the 8-methyl and 9-methyl groups of OBCP-3M(MonA) are equivalent to those made by the 8-methyl and 6-methyl groups, respectively, of OBCP-3M(MonB). The third methyl group in each diastereomer, the 6-methyl of OBCP-3M(MonA) and the 9-methyl of OBCP-3M(MonB), make unique contacts with the receptor that have no equivalent in the other diastereomer. However, these interactions cause only local changes to the ligand-binding structure and do not affect the overall conformation of the two monomers. Selection for one diastereomer of OBCP-3M over the others in each monomer is an interesting crystallographic finding, because such microheterogeneity in the protein samples often interferes with crystallization or results in poor crystal order (52).
Comparative structural analyses of the three OBCP-ER
LBD structures with other ER structures suggest that the ER
selectivity of the OBCP ligands can be attributed to the two residues that differ in the ligand-binding pockets of the ER subtypes. Of these, the difference at Leu-384(ER
)/Met-336(ER
) probably contributes the most to the subtype selectivity of the OBCP ligands. Among the ER
structures, the Leu-384 side chain positions are essentially unchanged whether OBCP ligands or non-subtype-selective ligands, such as E2 or DES, are bound (Fig. 5, B and C). Leu-384 does not interact with E2 and interacts weakly with DES. Similarly, the OBCP ligands interact with Leu-384 through weak van der Waal contacts made by the 8-methyl groups of OBCP-3M(MonB), OBCP-2M, and OBCP-1M. However, superimposing ER
(Protein Data Bank code 1X7J) onto the OBCP-ER
structures shows that the change to Met-336 in ER
at the equivalent position to Leu-384 in ER
places the larger methionine close enough (
34 Å) to make van der Waal contacts with the OBCP ligands (Fig. 8). The flexibility observed for the Met-336 side chain (14) suggests that it can easily move to accommodate the oxabicyclic scaffold of the OBCP ligands. Thus, Met-336 could make extensive contacts with the hydrophobic surfaces of these ligands. These interactions would favor the binding of OBCP ligands to ER
. The additional contacts made by the extra methyl groups of OBCP-3M and OBCP-2M probably account for their increased affinity and improved ER
selectivity compared with OBCP-1M. In comparison, the bulk of the DES and E2 structures are too far away to make productive contacts with Met-336; E2 has only a 13
-methyl group and DES an ethyl group with which Met-336 could interact.
|
) also affects the ER
selectivity of the OBCP ligands. Superimposed structures of DES- or E2-ER
LBD with the OBCP structures show that the OBCP ligands would sterically clash with Met-421 if the residue is positioned as seen in structures with non-subtype-selective ligands (DES/E2) (Fig. 7, B and C). These steric clashes are with the ring atoms on the oxabicyclic scaffold (
2.53.5 Å) and the methyl of the hydroxymethyl group (<2 Å) of these compounds. To minimize such steric collisions with the OBCP ligands, the Met-421 side chains in the OBCP complexes are rotated away from the positions that they adopt in the DES/E2 structures (Fig. 7). However, in the OBCP-3M structure, the OH of the hydroxymethyl group is also positioned in close proximity to Met-421 (3.8 Å). This close approach probably results in electrostatic repulsion between the partial negative charges carried by the oxygen of the OH group and the sulfur of the methionine group. The hydroxymethyl group of OBCP-3M is limited in its ability to reduce this repulsion by rotating away because of steric constraints imposed by the two surrounding 6- and 9-methyl groups. In contrast, with only one or no methyl groups nearby, the hydroxymethyl of OBCP-2M and OBCP-1M can freely rotate to position itself away from the sulfur of Met-421 (5.0 Å) and toward His-524. Here, the oxygen of the hydroxymethyl experiences weaker electrostatic repulsions with the methionine and can hydrogen-bond with His-524. The replacement of Met-421 in ER
by Ile-373 in ER
would eliminate the repulsive forces experienced by the OBCP-3M OH because isoleucine is nonpolar and shorter and occupies less volume than methionine. This substitution relieves OBCP-3M of its electrostatic repulsion in ER
(but minimally affects OBCP-1M and OBCP-2M) and would contribute to the stronger selectivity of OBCP-3M for ER
.
ER
-selective agonists can use three major mechanisms to take advantage of the two residues that differ between the ligand-binding pockets of ER
and ER
. The first mechanism is to promote unfavorable interactions with ER
. This is the primary mechanism used by the benzoxazole/benzofuran compounds (17) (Fig. S1), in which acetonitrile/vinyl groups are positioned near Met-421, and by 8
-vinyl estradiol (11) (Fig. S1), in which a vinyl group clashes with Leu-384. The second mechanism is to promote favorable interactions with ER
. These are the proposed mechanisms used by DPN (53) (Fig. 1) and genistein (21), in which a polarizable nitrile group and an aromatic (B) ring, respectively, interact with Met-336. The third mechanism is to combine the aforementioned two mechanisms. We have observed this combination with OBCP-3M, and it has also been described as the secondary mechanism for the benzoxazole/benzofuran ligands (17). Our structural analyses suggest that the ER
selectivity of the OBCP ligands results mainly from favorable interactions with Met-336 in ER
, because this interaction is common to these ligands. However, the greater ER
selectivity of OBCP-3M compared with its analogs is probably achieved by the additional unfavorable interaction made with Met-421 in ER
. Although the ER
selectivity of OBCP-3M is modest, we believe that optimization of each of its interactions would improve its selectivity.
The identification of OBCP-3M and its analogs adds to the repertoire of bicyclic ligands that have been reported (5457), including a very recent report on similar oxabicyclic compounds (22). Promising results from our biological studies and those of others (22) suggest that the oxabicyclic compounds have potential as lead compounds in the development of subtype-selective ligands for investigative or therapeutic purposes. Our crystallographic studies show for the first time how ligands with bulky bicyclic scaffolds can be accommodated in the flexible ligand-binding pocket of ER and suggest how these ligands interact with both Met-421(ER
) and Met-336(ER
) to achieve ER
selectivity. These findings should stimulate interest in ligands that can utilize regions of the ER ligand-binding pocket that remain unoccupied when planar ligands are bound (13). Our studies should facilitate the structure-based design of more potent and subtype-selective ER modulators belonging to this relatively new structural class.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Fig. S1.
The atomic coordinates and structure factors (code 1ZKY, 2FAI, and 2B1V) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/). ![]()
1 Supported by the University of Chicago unendowed Medical Scientist Training Program. ![]()
2 To whom correspondence may be addressed: Dept. of Chemistry and Howard Hughes Medical Institute, University of Chicago, 929 E. 57th St., Chicago, IL 60637. Tel.: 773-702-1651; Fax: 773-702-1677; E-mail: mmrksich{at}uchicago.edu. 3 To whom correspondence may be addressed: Ben May Institute for Cancer Research, University of Chicago, W330, 929 E. 57th St., Chicago, IL 60637. Tel.: 773-702-6964; Fax: 773-834-9029; E-mail: ggreene{at}uchicago.edu.
4 The abbreviations used are: ER, estrogen receptor; hER, human ER; E2, 17
-estradiol; SAR, structure-activity relationship; LBD, ligand-binding domain; GRIP1, glucorticoid receptor-interacting protein 1; DPN, diarylpropionitrile; PPT, propylpyrazole triol; SRC1, steroid receptor coactivator 1; GST, glutathione S-transferase;
-ME,
-mercaptoethanol; OBCP-xM, oxabicyclic phenol with x methyls on the core scaffold; HTRF, homogenous time-resolved fluorescence; AR, androgen receptor; NRD, nuclear receptor interaction domain; ERE, estrogen response element; ARE, androgen response element; DES, diethlystilbestrol; MonA, monomer A; MonB, monomer B. ![]()
5 S. S. Rajan, K. W. Nettles, R. W. Hsieh, and Geoffrey L. Greene, unpublished results. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|