![]()
|
|
||||||||
J. Biol. Chem., Vol. 282, Issue 32, 23231-23239, August 10, 2007
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
in Complex with a Synthetic Inverse Agonist Reveals Its Novel Molecular Mechanism*From the Novartis Institutes for BioMedical Research, CH-4002 Basel, Switzerland
Received for publication, April 20, 2007 , and in revised form, June 6, 2007.
| ABSTRACT |
|---|
|
|
|---|
(ERR
, NR3B1) are of potential interest for several disease indications (e.g. breast cancer, metabolic diseases, or osteoporosis). ERR
is constitutively active, because its ligand binding pocket (LBP) is practically filled with side chains (in particular with Phe328, which is replaced by Ala in ERR
and ERR
). We present here the crystal structure of the ligand binding domain of ERR
(containing the mutation C325S) in complex with the inverse agonist cyclohexylmethyl-(1-p-tolyl-1H-indol-3-ylmethyl)-amine (compound 1a), to a resolution of 2.3Å. The structure reveals the dramatic multiple conformational changes in the LBP, which create the necessary space for the ligand. As a consequence of the new side chain conformation of Phe328 (on helix H3), Phe510(H12) has to move away, and thus the activation helix H12 is displaced from its agonist position. This is a novel mechanism of H12 inactivation, different from ERR
, estrogen receptor (ER)
, and ER
. H12 binds (with a surprising binding mode) in the coactivator groove of its ligand binding domain, at a similar place as a coactivator peptide. This is in contrast to ERR
but resembles the situation for ER
(raloxifene or 4-hydroxytamoxifen complexes). Our results explain the novel molecular mechanism of an inverse agonist for ERR
and provide the basis for rational drug design to obtain isotype-specific inverse agonists of this potential new drug target. Despite a practically filled LBP, the finding that a suitable ligand can induce an opening of the cavity also has broad implications for other orphan nuclear hormone receptors (e.g. the NGFI-B subfamily). | INTRODUCTION |
|---|
|
|
|---|
(ERR
; NR3B1) was the first orphan NR to be identified on the basis of its similarity with estrogen receptor
(ER
; NR3A1) (3). ERR
and its relatives ERR
(NR3B2) and ERR
(NR3B3) form a small family of orphan NRs that are evolutionarily related to the estrogen receptors ER
and ER
. ERRs preferentially bind to DNA sites composed of a single half-site preceded by three nucleotides with the consensus sequence TNAAGGTCA, referred to as an ERR-response element. It has been shown that ERR
also efficiently binds to estrogen-response elements and that these receptors share common target genes (4). This observation was further supported by studies demonstrating cross-talk between the ER and ERR pathways (reviewed in Ref. 5). Several studies have highlighted ERR
as a main player in mitochondrial biogenesis and oxidative metabolism (6–8), suggesting that ERR
could be used for therapeutic intervention in diabetes or metabolic diseases. In addition to this central role in metabolism, ERR
is also now accepted as an emerging target in cancer (9, 10). Finally, few papers have suggested a role of ERR
in bone metabolism (reviewed in Ref. 11). The importance of ERR
as a drug target has been recently reviewed (12), which further re-emphasizes the urge for new synthetic ligands in the ERR subfamily.
Despite their overall sequence similarity with the ERs, ERRs seem to regulate transcription in the absence of known natural agonist ligands. The presence of a phenylalanine residue (Phe328 on helix H3) in the ligand binding pocket (LBP) has been found to be essential for the constitutive activity of ERR
, because its mutation to alanine abolishes constitutive activity (13). The x-ray structure of apoERR
, in complex with a coactivator peptide containing the L3 site of peroxisome proliferator-activated receptor coactivator-1
(PGC-1
), had revealed that the LBP is practically filled with side chains (14). The unoccupied volume was found to be only
100 Å3. It was predicted that a ligand with an equivalent of more than
4–5 carbon atoms could only bind if Phe328 would drastically change its conformation. This would also require a displacement of Phe510 and thus of H12, making the ligand an inverse agonist. The term "inverse agonist" (instead of "antagonist") is used because such a ligand would display an intrinsic activity, namely an inhibition of the constitutively active ERR
. In contrast a neutral antagonist does not have an intrinsic activity by itself, it just counteracts the binding of an agonist.
Indeed, searches for ERR
ligands have so far mostly identified inverse agonists. Few synthetic ligands, classically linked to the ER pathway, have been shown to modulate ERR
activity, namely 4-hydroxytamoxifen (4-OHT) and diethylstilbestrol (DES) (reviewed in Ref. 15). In addition, a new synthetic ligand acting as an inverse agonist has been identified for ERR
. It was shown that it acts at submicromolar activities in a cell-based assay and that this compound inhibits ERR
induction of an ERR
target gene, monoamine oxidase B (16, 17).
For ERR
, the crystal structures for the complexes of its ligand binding domain (LBD) with the inverse agonists DES and 4-OHT (18), a 4-OHT derivative (19), and 4-OHT together with a corepressor peptide (20) have been published. The cavity volume of the LBP for unliganded ERR
was reported as 220 Å3 (21), which is small but still considerably larger than the 100 Å3 described for apoERR
. This difference is mainly explained by the substitution of Phe328 (ERR
) with Ala350 (ERR
). In all reported complexes of ERR
with inverse agonists, the activation helix H12 was found to be completely dissociated from the LBD body (18–20) and also partially disordered (20). The ERR
structures are thus in contrast to selective ER modulator (SERM) complexes for ER
, where binding of raloxifene or 4-hydroxytamoxifen revealed H12 in the coactivator groove. In the latter ER
complexes, accompanying structural adaptations of the C-terminal end of H11 and of the H11/H12 loop (22, 23) were observed. For ER
LBD complexed with the full antagonist ICI 164,384 (ICI), H12 was not visible because of high mobility (24).
For ERR
it was difficult to predict the exact details of how (and whether at all) a ligand (most likely acting as an inverse agonist) would bind in the LBP, because of the multiple conformational changes required to create the necessary space and interactions. Similarly, it was not clear what exact consequences ligand binding would have on H12, except that it would probably be displaced from the agonist position.
Here we report the crystal structure of the human ERR
LBD in complex with the synthetic inverse agonist cyclohexylmethyl-(1-p-tolyl-1H-indol-3-ylmethyl)-amine (Fig. 1C, compound 1a), at 2.3 Å resolution. Compound 1a is a derivative of compound 1b (Fig. 1C), which lacks the methyl group and originally was discovered by high throughput screening. We have introduced the mutation C325S in ERR
in order to reduce biochemical instability problems during purification and crystallization, associated with cysteine oxidation. Isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC) measurements have demonstrated that there are no significant differences in the thermodynamic binding parameters between wild type ERR
LBD and the C325S mutant for compound 1a (or 1b). Also binding studies by NMR spectroscopy of compound 1b to wild type ERR
in solution (with and without a spin label attached to Cys325) revealed no violations of the observed relaxation effects with the crystal structure of the compound 1a complex.
Our results explain the novel molecular mechanism of an inverse agonist for ERR
and provide the basis for rational drug design to obtain selective inverse agonists of this potential new drug target.
| MATERIALS AND METHODS |
|---|
|
|
|---|
LBD, amino acids 290–519, numbering according to Swiss-Prot P11474
[GenBank]
(14)) by site-specific mutagenesis using the QuikChange mutagenesis kit from Stratagene. The correctness of the construction was verified by DNA sequence analysis, and the correct clone was called pXI392b. Generation of recombinant baculovirus by transfection with BacPAK8TM, plaque cloning, and amplification was done as described (14). For production of ERR
(C325S) LBD-expressing cells, 1-liter shake flask cultures with Sf9 cells adapted to SF900 II (Invitrogen) at a density of 1.5 x 106 c/ml were infected with baculovirus at a multiplicity of infection = 1 and were cultured for 3 days at 28 °C and 90 rpm, as described previously (14). Harvest of the cells was carried out as described (14). SDS-PAGE analysis (Coomassie staining and Western blotting with anti-penta-His antibody from Qiagen) revealed good and soluble expression of ERR
(C325S) LBD. The purification of the frozen Sf9 cell pellets was done by three-step chromatography as described previously (14). The purified ERR
LBD protein in 50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 5 mM DTT was concentrated to 10.4 mg/ml prior to crystallization. The resulting protein was estimated to be >95% pure and homogeneous by SDS-PAGE and reverse phase high pressure liquid chromatography coupled to mass spectrometry. The measured molecular mass of 27,025.3 Da of the protein corresponded to the acetylated des-Met form of the ERR
(290–519; C325S) LBD.
Crystallization, Data Collection, and Structure Determination—Compounds 1a and 1b were synthesized in-house. For cocrystallization, a 5:1 stoichiometric excess of compound 1a from a 100 mM stock solution in Me2SO was added to the protein solution. Crystals were obtained at 4 °C by the vapor diffusion method in 1-µl hanging drops containing equal volumes of protein (10.4 mg/ml ERR
(C325S) LBD, 100 mM NaCl, 5 mM DTT, 50 mM Tris-HCl, pH 7.5) and crystallization buffer (0.2 M MgCl2, 12% PEG400, 0.1 M HEPES, pH 7.5). The initial crystallization hit was obtained with Crystal Screen I from Hampton Research. Crystals were directly mounted in cryoloops and flash-frozen in the nitrogen stream. Diffraction data at 100 K were collected at the Swiss Light Source (beamline X10SA), using a Marresearch CCD detector and an incident monochromatic x-ray beam with 0.8 Å wavelength. Raw diffraction data were processed and scaled with the HKL program suite version 1.96.1 (26). The estimated B-factor by Wilson plot analysis is 48.6 Å2. The structure was determined by molecular replacement with MOLREP (27, 28) using as starting model the coordinates of ERR
LBD (Protein Data Bank access code 1XB7
[PDB]
) refined to 2.5 Å resolution (14), with H12 removed. The program REFMAC version 5.0 (28, 29) was used for refinement. Bulk solvent correction, an initial anisotropic B-factor correction, and restrained isotropic atomic B-factor refinement were applied. The refinement target was the maximum likelihood target using amplitudes. No
cut-off was applied on the structure factor amplitudes. Cross-validation was used throughout refinement using a test set comprising 5.1% (1056) of the unique reflections. Water molecules were identified with the program ARP/wARP (28, 30) and selected based on difference peak height (greater than 3.0
) and distance criteria. Water molecules with temperature factors greater than 60 Å2 were rejected. The program O version 7.0 (31) was used for model rebuilding, and the quality of the final refined model was assessed with the programs PROCHECK version 3.3 (32) and REFMAC version 5.0 (28, 29). Crystal data, data collection, and refinement statistics are shown in Table 1.
|
|
LBD (wild type and C325S mutant at 40 and 30 µM) in 50 mM Tris, pH 7.5, 100 mM NaCl, and 0.5 mM TCEP. The syringe was loaded with compound 1b (400 µM) in the same buffer. All solutions were degassed for 10 min. Titrations were performed at 25 °C with injection volumes of 8 µl and a spacing of 240 s. The base line was set to zero assuming that the final injections of each titration represent only the heat of dilution. The data were fit using a one-site binding model available in the Origin ITC data analysis software (version 7.0). DSC analysis was performed on buffer solutions (50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 1 mM TCEP) containing the ERR
LBD-compound 1a complex (1:10 molar ratio). The ERR
LBD concentration was 25 µM (both for wild type or the C325S mutant), and 300 µl of protein complex solution were injected. DSC scans were obtained using a MicroCal VP-capDSC system (MicroCal, LLC, Northampton, MA) at a scan rate of 200 °C/h. The fluorescence energy transfer (FRET) measurements were conducted as time-resolved measurements (TR-FRET) in a miniaturized 1536 well plate format. Typically, 100 nl of compound solution (2 mM in 90% Me2SO) were transferred to the assay plate. 3 µl of assay buffer (50 mM Tris-Cl, pH 7.5, 100 mM NaCl, 0.05% bovine serum albumin, 5 mM DTT, 0.1% Pluronic F-127) containing LANCE EU-W1024-labeled anti-His6 antibody (PerkinElmer Life Sciences; 1 nM final assay concentration) were added. 1 µl of His-ERR
LBD (50 nM final assay concentration) was added, and the plates were left for incubation at 20 °C for 30 min. As FRET acceptor, 1 µl of an N-terminally Cy5-labeled PGC-1
derived peptide (Cy5-RPCSELLKYLTT, 50 nM final assay concentration) was added. After 20 min of incubation the TR-FRET readout was performed on an Envision 2102 multilabel reader (PerkinElmer Life Sciences) with an excitation at 350 nm, first emission (donor) at 615 nm, and second emission (acceptor) at 665 nm (30 flashes, 100-µs delay time). The readout was calculated according to X = donor/acceptor x 1000. Data were fitted with XLfit4 after plate normalization of the derived data.
NMR Spectroscopy—Ligand binding studies were performed with wild type and spin-labeled ERR
LBD. The spin-labeled protein was obtained by addition of a maleimide-substituted tetramethyl pyrrolidine N-oxide (TEMPO) paramagnetic tag (Aldrich) that reacts selectively with the only cysteine residue present (Cys325). Selective reaction was completed within minutes at room temperature prior to the measurement. T1
relaxation and water-LOGSY experiments were performed at 600 MHz using a room temperature probe.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
LBD and the Inverse Agonist Compound 1a Reveals H12 in a Position Similar to PGC-1
—We have solved the crystal structure of ERR
LBD (containing the mutation C325S) in complex with the inverse agonist cyclohexylmethyl-(1-p-tolyl-1H-indol-3-yl-methyl)-amine (compound 1a) (Fig. 1C), at 2.3 Å resolution (space group H3). Compound 1a is a derivative of compound 1b, which originally was discovered by high throughput screening. The asymmetric unit contains a homodimer of ERR
complexes, i.e. ligand binding did not interfere with dimer formation. The results of the crystallographic refinement are summarized in Table 1. In general the electron density for the ERR
LBD is well defined, except for the His tag, residues Pro309–His317 (H2/H3 loop), residues Arg462–Glu470 (H9/H10 loop), and the C-terminal residues Met518–Asp519 (all of which were not modeled). The protein part of the refined model consists of the PreScissionTM site (LEVLFQGP) followed by amino acids Val290–Met308, Leu318–Gly461, and Arg471–Met517 of the ERR
LBD. The refined model also contains 236 water molecules and two compound 1a molecules. The two complexes in the asymmetric unit are very similar, with a root mean square deviation (r.m.s.d.) of 0.13 Å for 210 C
atoms of residues 291–308, 332–461, and 471–498.
A comparison of apoERR
(bound to a PGC-1
peptide (14)) with the inverse agonist complex (Fig. 1) revealed no major conformational changes for the main chain atoms from H1 to the middle of H11 (r.m.s.d. of 0.55 Å for 210 C
atoms as above), with the exception of H3. The N-terminal part of H3 moves away from the LBP, in order to create necessary space for the ligand (Fig. 1). In particular the C
s of Val321 and Leu324 move by 4.3 and 3.2 Å, respectively. Important structural perturbations of the main chain due to compound 1a binding also occur at the C terminus of H11 (unwinding of the last turn of H11, i.e. after Lys499), the H11/H12 loop (reorganization of its conformation, so that Pro505 is now at the beginning of H12), and H12 (new location in the coactivator binding groove) (Fig. 1). In the apoERR
complex with PGC-1
(14), H12 is formed between His507 and Ala516 (i.e. Pro505 and Met506 are in the H11/H12 loop), whereas in the inverse agonist complex H12 already starts at Pro505 and is visible until Met517. In the apoERR
structure, the PGC-1
peptide (which adopts an
-helical conformation) utilizes two Leu side chains (Leu210 and Leu214 of its inverted LXXLL-motif) to make hydrophobic interactions in the coactivator groove. Unexpectedly, similar interactions are made by two Met side chains (Met513 and Met517) of H12 in the compound 1a complex. The hydrogen bond made between Lys340(H3) and the PGC-1
main chain (contributing to the "charge clamp" interaction) is replaced by a hydrogen bond between Lys340 and the carbonyl oxygen of Ala516(H12). H12 in the compound 1a complex is making further stabilizing interactions via a salt bridge between the side chains of Lys508(H12) and Asp329(H3) and via a hydrogen bond between the carbonyl oxygen of Met506(H12) and the side chain of Gln358(H4). Finally, in the inverse agonist complex Leu509(H12) is located in the same hydrophobic pocket (with Trp361 at the bottom) as Met513(H12) in the apoERR
structure. In the crystal lattice, Phe510 makes an aromatic stacking interaction with Phe510 of a neighboring molecule, but these packing interactions seem to be of less importance than the interactions in the coactivator groove itself (described above).
The finding that for the ERR
-compound 1a complex H12 adopts a well defined position in the coactivator groove is in contrast to ERR
, where the binding of diethylstilbestrol or 4-hydroxytamoxifen led to complete dissociation of H12 from the LBD (18–20). On the other hand, the ERR
structure resembles the situation for ER
, where binding of the SERMs raloxifene or 4-hydroxytamoxifen revealed H12 in the coactivator groove, with accompanying structural adaptations of the C-terminal end of H11 and of the H11/H12 loop (22, 23). For ER
LBD complexed with the full antagonist ICI, H12 was not visible because of its high mobility (24).
Compound 1a Binds in the LBP of ERR
by Displacing Phe328(H3) and Phe510(H12); the Novel Mechanism of an ERR
Inverse Agonist—The empty cavity of the LBP in apoERR
has a volume of only
100 Å3 (14), so that multiple structural adaptations are required to enable ligand binding. In apoERR
the aromatic side chains of Phe328(H3), Phe495(H11), and Phe510(H12) cluster (Fig. 2A) and in particular the presence of Phe328 (which is replaced by Ala for ERR
and ERR
) leads to an almost complete occupation of the LBP with side chains. Compound 1a binding induces a new side chain conformation for Phe328(H3) and as a consequence a displacement of Phe510 (to avoid a steric clash with Phe328), which triggers a dramatic movement of H12 (Fig. 2B). Phe495(H11), on the other hand, has almost the same position in both structures. In the compound 1a complex, Phe495 makes an aromatic stacking interaction with Phe328, and they form together a hydrophobic lid on top of the ligand (Fig. 2B). In addition to the reorientation of Phe328, also the N-terminal part of H3 moves away from the LBP, in order to create space for the ligand via removal of Val321 and Leu324 (Fig. 3B). Taken together, Phe510 could be regarded as a sensor that transduces the conformational state of Phe328 (which depends on the ligand occupation of the LBP) into a position of H12 (turning the transcriptional activity of ERR
"on" or "off" (Fig. 1)).
|
-sheet (Phe382, Leu386), H6 (Ala396), the H6/H7 loop (Leu398), H7 (Leu405), H11 (Val491, Phe495, Val498), and the H11/H12 loop (Leu500) (Fig. 3A). All residues in the LBP have well defined electron density. The amino nitrogen of compound 1a makes a direct salt bridge and a water-mediated hydrogen bond interaction with Glu331 (Fig. 3A). The latter water molecule is well ordered (B-factor = 25 Å2) and is further stabilized by a hydrogen bond with the carbonyl oxygen of Phe382. Glu331 of ERR
corresponds to Glu353 of ER
, for which it forms e.g. a hydrogen bond with the hydroxyl of 4-OHT (23). Importantly, the side chain of Glu331 reorients toward the ligand in order to form interactions with the amino nitrogen (Fig. 3B). This reorientation is made possible because of the side chain movement of Phe328 (Fig. 3B), which creates space for Glu331 (and the ligand). The cyclohexyl moiety of compound 1a is deeply embedded inside ERR
in a mainly hydrophobic pocket. There is some space left for substituents on the cyclohexyl pointing in the direction toward Arg372. The indole moiety of compound 1a is nicely fitting in a hydrophobic pocket, which was occupied by Leu324 in the apoERR
structure (Fig. 3B). Leu324 now contributes to one side of the indole pocket. Finally, the tolyl group is the moiety that is closest to the solvent, with the methyl group at the entrance of the LBP. The carbonyl oxygens of Ala396 and Gly397 on the H6/H7 loop are potential targets to obtain additional hydrogen bond interactions. In summary, the LBP has adapted via multiple conformational changes to the presence of compound 1a. The complementarity of fit in the hydrophobic regions is good, and the salt bridge interaction with Glu331 is likely to contribute substantially to the binding affinity.
ITC Measurements Reveal That Binding Is Enthalpically Driven—ITC measurements have revealed that the binding of compound 1b to wild type ERR
LBD is enthalpically driven with
H0 =–12.7 kcal/mol (Fig. 4A). A negative enthalpic contribution is compatible with the formation of an important salt bridge (if desolvation of the partners is not too unfavorable). Of course, because
H0 is composed of many contributions, it is not possible to strictly deduce the formation of a salt bridge. The entropic contribution to binding is unfavorable (–T
S0 = 4.6 kcal/mol). This is evidence for the hypothesis that the apo-form is more flexible than the ligand-bound form and that H12 does not become more disordered upon inverse agonist binding. Indeed the x-ray structure shows a defined position of H12 in the coactivator groove. In addition, the unfavorable entropic term indicates that desolvation (e.g. for the salt bridge partners) does not play a dominant role. Importantly the Kd values, as well as the separate enthalpic and entropic contributions, are very similar for ERR
wild type and the C325S mutant (Fig. 4) (Kd values of 770 and 930 nM, respectively). This indicates that the binding mode has not been modified by the introduction of the mutation C325S. DSC measurements (data not shown), comparing compound 1a binding to ERR
wild type and the C325S mutant, also confirmed very similar affinities (
Tm values of
1.7–1.8 °C in both cases). The IC50 value for compound 1b binding as determined in the FRET assay is 700 nM, which is close to the Kd of 770 nM determined by ITC. Compound 1a has an improved IC50 value of 190 nM, which can be explained by additional van der Waals interactions (with the side chains of Val321, Val498, and Leu500) contributed by its methyl group.
|
|
LBD, the binding was confirmed by T1
relaxation and water-LOGSY experiments (33). The T1
relaxation spectra of the compound alone and in the presence of ERR
(10-fold excess of compound) are shown in Fig. 5, A and B. Apart from strong T1
relaxation effects observed in the presence of the protein, significant line broadening is also observed. This indicates that intermediate exchange processes occur during ligand binding. As we know from the crystal structure presented in this study, significant conformational rearrangements within ERR
are required for ligand binding. The line broadening is compatible with millisecond time scale motions, as is observed by the strong increase of T2 relaxation of the ligand.
Because the expected LBP is close to Cys325 (the only cysteine residue of ERR
LBD), we tried to obtain spatial information utilizing site-specific spin labeling in conjunction with T1
relaxation measurements (34). TEMPO-labeled ERR
was obtained after reaction with maleimide TEMPO, which reacts specifically with freely accessible cysteines. Significant relaxation effects induced by the paramagnetic center attached to Cys325 were measured on the individual resonances of the ligand (Fig. 5, C and D). Strong effects were observed for the resonances of the indole and phenyl groups, and weak effects were observed for the aliphatic ring, indicating that the ligand binds in proximity of Cys325 (with the cyclohexyl ring farthest away). The crystal structure of the compound 1a complex shows that the phenyl and indole groups are about 6–10 Å away from Cys325, whereas the aliphatic ring is about 13 Å away. The paramagnetic relaxation measured at the resonances of compound 1b bound to wild type ERR
in solution thus is compatible with the binding mode observed in the crystal structure of compound 1a in complex with the C325S mutant ERR
. The reactivity of Cys325 with maleimide TEMPO demonstrates the accessibility of this residue. In the crystal structure of apoERR
, however, this residue is not accessible, because of the position of H12 and the H11/H12 loop. The fast reaction indicates that apoERR
is in a dynamic equilibrium between several states in which H12 can move away, which is important for opening up the LBP to enable ligand binding.
|
and ERR
; H3 Behaves Differently and Different Mechanisms of H12 Displacement—For the 4-OHT and DES complexes of ERR
(18–20), the N-terminal part of H3 was not displaced with respect to apoERR
. It also has a very similar position as in apoERR
(14). In particular Leu268(H3) in the inverse agonist ERR
complexes has practically the same position as the corresponding Leu324(H3) for apoERR
. In contrast, for the compound 1a complex of ERR
, Leu324 is displaced by the indole moiety, and Val321 is displaced by the tolyl moiety (Fig. 3B). Compound 1a binding thus induces in ERR
a different shape and location of the LBP (generated by a movement of the N-terminal part of H3), compared with the 4-OHT and DES complexes of ERR
.
For ERR
, DES- and 4-OHT-mediated antagonism is caused by the rotation of the side chain of Phe435(H11), which upon ligand binding flips out of the LBP and sterically interferes with H12 in the agonist position (18–20). In contrast, for ERR
the side chain of Phe495 (which corresponds to Phe435 of ERR
) practically does not move upon ligand binding (Fig. 3B). Compound 1a does not make a steric clash with Phe495, rather there are favorable interactions between this side chain and the ligand. For ERR
, it is the new side chain conformation of Phe328 that acts as a trigger for H12 displacement, removing it from its agonist position (Fig. 2). For ERR
, the corresponding residue Ala272 does not have to move, because its small side chain does not lead to steric clashes with DES or 4-OHT. Docking shows also that compound 1a could fit in the ERR
LBP without movement of Ala272 (the N-terminal part of H3 would have to bend away from the LBP similarly as for ERR
). On the other hand, the hydrophobic lid generated by Phe328 and Phe495 of ERR
on top of compound 1a (Fig. 3A) is not possible for ERR
. This leads to the prediction that compound 1a should bind with a reduced affinity to ERR
(and ERR
), compared with ERR
.
Taken together, a comparison of the inverse agonist complexes of ERR
and ERR
reveals different mechanisms of action leading to H12 displacement. For ERR
the ligand (compound 1a) displaces Phe328(H3) and thus indirectly Phe510(H12), whereas for ERR
the ligands (DES and 4-OHT) displace Phe435(H11) and thus indirectly Leu454 and Phe450(H12). The mechanism of ERR
inactivation by compound 1a is also different from the "active" AF2 antagonists and SERMs of ER
or ER
(where e.g. the bulky extensions of raloxifene, 4-OHT, or ICI directly displace H12) or from the "passive" AF2 antagonists of ER
(18–20, 22–24, 35, 36).
Differences in H12 Positions and Interactions for the Inverse Agonist Complexes of ERR
, ERR
, and ER
—In the crystal structure of the ERR
-compound 1a complex H12 adopts a well defined position in the coactivator groove of its LBD. This is in contrast to ERR
, where the binding of DES or 4-OHT led to complete dissociation of H12 from the LBD (18–20). This dissociation was explained by the absence of an LXXLL "coactivator" motif in H12 of ERR
(18). By comparison, H12 of ER
does contain an LXXLL motif, for which Leu536 and Leu540, i.e. the first and last Leu of the motif, mediate anchoring of H12 in the coactivator cleft for the raloxifene and 4-OHT complexes (22, 23). For ERR
, the corresponding residues Met446 and Phe450 were thought not to allow similar stabilizing interactions of H12 with the coactivator groove (18), and thus to be an explanation for H12 dissociation. In contradiction to this hypothesis, ERR
has the same residues as ERR
in the above motif, namely Met506 and Phe510, but nevertheless displays H12 in the coactivator cleft of its LBD for the compound 1a complex. Actually, there is no sequence difference between ERR
and ERR
in the whole region between the ERR
residues Lys499 (C terminus of H11) and Ala516 (close to C terminus of H12). The first difference occurs at the C terminus of H12, with Met517 of ERR
replacing Lys457 of ERR
. Strikingly, in the ERR
complex with compound 1a, Met517 is well ordered and is practically at the same position as Leu214 of PGC-1
bound to ERR
(14). Lys457 of ERR
cannot make similar hydrophobic interactions in the coactivator groove, so this substitution is a possible explanation for the differences in H12 positions.
In addition to Leu536 and Leu540 of the LXXLL motif, ER
also utilizes Leu544 of H12 for interactions in the coactivator cleft. According to the sequence alignment, these residues would correspond to Met506, Phe510, and Leu514 of ERR
, but interestingly they are not pointing into the groove. Instead, Leu509, Met513, and Met517 of ERR
take up this role, enabled by a lengthening (compared with ER
) of the H11/H12 loop structure by three residues. Knowledge of these unexpected interactions can further help in the design of peptide or nonpeptide antagonists, which bind in the coactivator groove (instead of the LBP), but do not necessarily contain an LXXLL motif. Recently, ERR
-selective peptide antagonists containing an LXXLL motif were described, which were identified by screening combinatorial phage libraries (37).
| CONCLUSION |
|---|
|
|
|---|
LBD in complex with an inverse agonist, cyclohexylmethyl-(1-p-tolyl-1H-indol-3-ylmethyl)-amine (compound 1a), to a resolution of 2.3 Å. Our data reveal for ERR
the multiple conformational changes in the LBP induced by the binding of compound 1a. These changes are the trigger for a novel mechanism of H12 displacement. In addition, they show an unexpected movement of H3 leading to an enlargement of the LBP and a surprising mode of H12 interaction in the coactivator binding groove. The structure also provides the basis for rational drug design to obtain isotype-specific inverse agonists of this potential new drug target. The finding that, despite a practically filled LBP in apoERR
, a suitable ligand can induce an opening of the cavity has broad implications for other orphan NRs. In particular, it is of interest for the recently published LBD crystal structures in the NGFI-B family (38–40), for which the LBP is completely filled with four aromatic residues conserved within the subfamily. Given the dramatic structural modifications that can occur for ERR
, it is tempting to speculate that synthetic ligands (probably acting as inverse agonists) might also be found for members of the NGFI-B subfamily. | FOOTNOTES |
|---|
The atomic coordinates and structure factors (code 2PJL) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/). ![]()
1 To whom correspondence should be addressed. E-mail: Joerg.kallen{at}novartis.com.
2 The abbreviations used are: NR, nuclear hormone receptor; ERR, estrogen-related receptor; ER, estrogen receptor; LBP, ligand binding pocket; PGC-1
, peroxisome proliferator-activated receptor coactivator-1
; 4-OHT, 4-hydroxytamoxifen; DES, diethylstilbestrol; LBD, ligand binding domain; SERM, selective ER modulator; ICI, ICI 164,384; ITC, isothermal titration calorimetry; DSC, differential scanning calorimetry; FRET, fluorescence energy transfer; TR-FRET, time-resolved fluorescence energy transfer; TEMPO, tetramethylpyrrolidine N-oxide; r.m.s.d., root mean square deviation; DTT, dithiothreitol; TCEP, tris(2-carboxyethyl)phosphine; ITC, isothermal titration calorimetry. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Greschik, M. Althage, R. Flaig, Y. Sato, V. Chavant, C. Peluso-Iltis, L. Choulier, P. Cronet, N. Rochel, R. Schule, et al. Communication between the ERR{alpha} Homodimer Interface and the PGC-1{alpha} Binding Surface via the Helix 8-9 Loop J. Biol. Chem., July 18, 2008; 283(29): 20220 - 20230. [Abstract] [Full Text] [PDF] |
||||