|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 18, 13648-13655, May 4, 2007
Crystal Structure of the T877A Human Androgen Receptor Ligand-binding Domain Complexed to Cyproterone Acetate Provides Insight for Ligand-induced Conformational Changes and Structure-based Drug Design*![]() ![]() ![]() 1
From the
Received for publication, December 21, 2006 , and in revised form, February 20, 2007.
Cyproterone acetate (CPA) is a steroidal antiandrogen used clinically in the treatment of prostate cancer. Compared with steroidal agonists for the androgen receptor (AR) (e.g. dihydrotestosterone, R1881), CPA is bulkier in structure and therefore seemingly incompatible with the binding pockets observed in currently available x-ray crystal structures of the AR ligand-binding domain (LBD). We solved the x-ray crystal structure of the human AR LBD bound to CPA at 1.8Å in the T877A variant, a mutation known to increase the agonist activity of CPA and therefore facilitate purification and crystal formation of the receptor·drug complex. The structure demonstrates that bulk from the 17 -acetate group of CPA induces movement of the Leu-701 side chain, which results in partial unfolding of the C-terminal end of helix 11 and displacement of the loop between helices 11 and 12 in comparison to all other AR LBD crystal structures published to date. This structural alteration leads to an expansion of the AR binding cavity to include an additional pocket bordered by Leu-701, Leu-704, Ser-778, Met-780, Phe-876, and Leu-880. Further, we found that CPA invokes transcriptional activation in the L701A AR at low nanomolar concentrations similar to the T877A mutant. Analogous mutations in the glucocorticoid receptor (GR) and progesterone receptor were constructed, and we found that CPA was also converted into a potent agonist in the M560A GR. Altogether, these data offer information for structure-based drug design, elucidate flexible regions of the AR LBD, and provide insight as to how CPA antagonizes the AR and GR.
The AR2 is a ligand-inducible steroid hormone receptor involved in regulation of prostate growth, spermatogenesis, and bone and muscle mass. Ligands that block the actions of androgens are used clinically for the treatment of prostate cancer and include the steroidal antiandrogen, CPA (CyprostatTM), and the nonsteroidal antiandrogens flutamide (EulexinTM) and bicalutamide (CasodexTM). Mutations in the AR often occur in patients being treated with these compounds that increase their agonist activity and thus result in drug resistance (1-7). Paradoxically, discontinuation of drug therapy at this point may lead to tumor regression known as antiandrogen withdrawal syndrome. A number of antiandrogen withdrawal syndrome mutations occurs in regions of direct contact with the ligand that compensate for the increased bulkiness of antagonist compared with agonist compounds. Such mutations allow antagonist ligands to invoke the same AR LBD conformation as the agonist-bound AR LBD (1, 8, 9).
Molecular modeling studies have offered a number of potential docking solutions to explain how drugs more bulky than steroidal agonists can be accommodated within the AR (10-14). Our recent report of the crystal structure of the AR LBD complexed to S-3-(4-fluorophenoxy)-2-hydroxy-2-methyl-N-[4-nitro-3-(trifluoromethyl)phenyl] propionamide (hereafter referred to as S-1 as in Ref. 15), an aryl propionamide selective androgen receptor modulator that acts as an agonist for the AR, demonstrated that movement of Trp-741 provides the additional space necessary for the compound to bind (8). Similarly, CPA is seemingly too large a molecule to fit within the AR binding pocket due to its bulky substituents at the 17-position distinguishing it from steroidal AR agonists. The LNCaP prostate cancer cell line contains an AR with the T877A point mutation, which allows CPA and hydroxyflutamide (HF) to act as agonists for the AR, whereas bicalutamide remains an antagonist. Agonist activity of CPA observed in the PR may be attributed to the presence of a cysteine (Cys-891) instead of a threonine at this position. However, CPA still does not appear compatible with the binding pocket observed in the T877A AR (16) or the PR. The estrogen receptor and GR have been shown to assume different conformations when liganded to agonists and antagonists (9, 17) due to the ability of antagonists to partially unfold the LBD of the receptor. X-ray crystal structures demonstrating changes in protein folding of the AR LBD when complexed to different ligands however are not currently available. Here we report structural evidence that CPA binding to the T877A AR induces unraveling of the C terminus of helix 11 and displacement of the loop between helices 11 and 12. Based on these findings, we show that the AR mutation, L701A, confers agonist activity to CPA similar to T877A and that a similar mutation in the GR converts CPA into a potent agonist. Furthermore, identification of the expanded binding pocket occupied by CPA provides insight for exploiting receptor interactions using structure-based drug design.
Cloning, Expression, and PurificationAn AR-LBD-(663-919) was obtained by PCR amplification from a full-length AR expression construct (pCMVhAR; generously provided by Dr. Donald J. Tindall, Mayo Clinic and Mayo Foundation, Rochester, MN) with primers containing flanking restriction sites and inserted into the pGEX6P-1 plasmid vector (Amersham Biosciences). Mutations were created in the pGEX6P1-AR-(663-919), pCMVhAR, pCMVhGR, and pCMVhPR via the Stratagene QuikChange mutagenesis kit according to the manufacturer's instructions. AR LBD expression and purification were performed essentially as previously described (1). The AR LBD was expressed as a glutathione S-transferase (GST) fusion protein in Escherichia coli BL21(DE3) at 15 °C for 16 h by induction with 30 µM isopropyl 1-thio- -D-galactopyranoside. Cells were lysed in a buffer containing 150 mM NaCl, 50 mM Tris, pH 8.0, 5 mM EDTA, 10% glycerol, 1 mg/ml lysozyme, 10 units/ml DNase I, 10 mM MgCl2, 10 mM DTT, 0.5% CHAPS, 100 µM CPA (Sigma), and 100 µM phenylmethylsulfonyl fluoride by 3 cycles of freeze thaw. The supernatant from ultracentrifugation was incubated for 1 h at 4 °C with glutathione-Sepharose (Amersham Biosciences) and washed with 150 mM NaCl, 50 mM Tris, pH 8.0, 5 mM EDTA, 10% glycerol, 10 mM ATP, 100 µM CPA, 0.1% n-octyl- -glucoside, and 1 mM DTT. The GST-LBD fusion protein was cleaved in a buffer containing 150 mM NaCl, 50 mM Tris, pH 7.0, 10% glycerol, 100 µM CPA, 0.1% n-octyl- -glucoside, 1 mM DTT, and 5 units/mg of protein Pre-Scission protease (Amersham Biosciences) at 4 °C overnight releasing the AR LBD from the glutathione-Sepharose resin. The supernatant was then diluted 3-fold in 10 mM Hepes, pH 7.2, 10% glycerol, 100 µM CPA, 0.1% n-octyl- -glucoside, and 1 mM DTT and loaded onto an HP SP cation exchange-column (Amersham Biosciences). Protein was eluted with a gradient of 50 to 500 mM NaCl in the same dilution buffer. The buffer was exchanged in a Millipore 10-kDa cut-off concentrator to a buffer containing 150 mM Li2SO4, 50 mM Hepes, pH 7.2, 10% glycerol, 100 µM CPA, 0.1% n-octyl- -glucoside, and 10 mM DTT, and protein was concentrated to 8 mg/ml. Crystallization, Data Collection, and Structure DeterminationAR LBD crystals formed in 1-2 days using the hanging drop vapor diffusion method in reservoirs containing 0.1 M Hepes, pH 7.5, and 0.5 M to 0.8 M sodium citrate in drops of equal amounts of protein and reservoir solution. Crystals grew to maximum size within 2 weeks. Prior to flash freezing in liquid nitrogen, crystals were transferred to a solution consisting of 0.1 M Hepes, pH 7.5, 0.7 M sodium citrate, and 20% ethylene glycol. Diffraction data were collected using a Rigaku RU300 rotating anode generator and an R-axis IV++ image plate (Rigaku, The Woodlands, TX), and processed with Crystal Clear software (Molecular Structure Corp., The Woodlands, TX). The x-ray crystal structure of the WT AR LBD in complex with the nonsteroidal propionamide selective androgen receptor modulator, S-1 (8) (PDB code 2AXA), was used as a starting structure for refinement using the Crystallography & NMR System (18). After an initial round of refinement, electron density maps allowed for accurate fitting of the ligand. Model building and water molecules were added using the program O (19), and further rounds of refinement were performed using rigid body, torsion angle simulated annealing, and individual temperature factor modules of the Crystallography & NMR System. The figures were prepared with PyMOL.3
Cotransfection AssayCV-1 cells (green monkey kidney cells) were grown in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum, 2 mM glutamine, and 1% streptomycin and penicillin to confluency in 75-cm2 tissue culture flasks. Cells were then transfected with 3 µg of the CMVhAR, CMVhGR, or CMVhPR expression vectors (WT or mutant), 30 µg of a luciferase reporter construct (pMMTV-luc, generously provided by Dr. Ronald Evans at The Salk Institute, San Diego, CA), and 30 µg of a Docking StudiesStructures were prepared for docking using Sybyl 7.1 (Tripos Inc., St. Louis, MO) by removing the ligand and water molecules, adding polar hydrogens, and assigning Kollman united atom charges. Ligands were extracted from PDB files and modified with the Sketch Molecule tool in Sybyl, and Gasteiger-Marsili charges were assigned. AutoDock3.05 (20) was used to perform docking simulations using default parameters with a grid centered around the x, y, and z coordinates of 29, 0, and 4 Å in the WT·DHT (PDB code 1I37 [PDB] ), T877A·DHT (PDB code 1I38 [PDB] ), and T877A·CPA complexes. The solutions with the lowest free energy of binding of 10 clusters were reported.
T877A-AR LBD·CPA ComplexThe final model of the T877A AR LBD·CPA crystal structure at 1.8-Å resolution consists of a single monomer of AR residues 671-919, which includes 2027 protein atoms in the asymmetric unit. Electron density maps for AR residues 886-890 in the T877A·CPA were weak and significantly displaced from the starting model (S-1-bound WT AR LBD). The backbone of these amino acids was fit using the 2Fo - Fc map at the 0.6 level, but side chains were omitted in the model. The model also includes 134 water molecules. CPA was fit unambiguously in a single conformation using the 2Fo - Fc map contoured at the 1.0 level. The refined temperature factors for the atoms of CPA were similar to those of surrounding protein atoms indicating full occupancy of the ligand. Crystallography statistics are listed in Table 1.
Binding Conformation of CPA to the T877A ARA hydrogen bond network is present in the x-ray crystal structure between the 3-keto group of CPA with Arg-752, Gln-711, and a water molecule (see Fig. 2a). The chlorine atom on the C6-position is accommodated in a groove surrounded by hydrophobic residues of helices 4, 7, and 11. The oxygen atom of the -acetyl group branched from the 17-position of CPA is oriented for a hydrogen bond with the nitrogen 2 of Asn-705, whereas the methyl of this acetyl group denoted C-21 (Fig. 1) on CPA is located only 3.2 Å from Ala-877. No hydrophilic interactions are available for the 17 -acetate group of CPA, which binds in a cavity bordered by Leu-701, Leu-704, Met-780, Phe-876, and Leu-880. The backbone oxygen of Ser-778 is located <4 Å from the C-24 methyl group of the CPA acetate, but the side chain of this residue faces away from the ligand to form a hydrogen bond with the backbone nitrogen of Met-780.
Comparison of T877A·CPA and WT·DHT AR LBDThe steroid scaffold of CPA binds in the same plane as DHT (Fig. 2b). The increased bulkiness of CPA in comparison to DHT is a result of a number of structural modifications, including a three-member ring incorporating the 1- and 2-positions of the steroid core structure, a chlorine group at the 6-position, an acetyl group substituted in place of the 17
Transcriptional Activation of CPA in the WT, T877A, and L701A ARStructural evidence that Leu-701 is displaced by CPA upon AR binding and causes a significant effect on the packing of the loop between helices 11 and 12 led us to the hypothesis that interaction with Leu-701 is important for the antiandrogenic properties of CPA similar to Thr-877. The L701A AR mutant was therefore constructed to investigate whether this mutation would also increase the ability of CPA to stimulate transcriptional activation at low concentrations. In the WT AR, DHT elicited a maximum response by 0.1 nM (Fig. 3a). CPA elicited only a small percentage of transcriptional activation at 100 nM and slightly more at 1000 nM relative to DHT, whereas HF and R-bicalutamide showed nearly no stimulation even at 1000 nM. In the T877A mutant, both CPA and HF elicited a substantial response at 10 nM (Fig. 3b). DHT maintained its potent activity in the T877A AR, and R-bicalutamide showed no response. CPA also induced a drastic increase in transcriptional activation in the L701A AR compared with the WT (Fig. 3c), yet no activity was observed for HF or R-bicalutamide in this mutant. Furthermore, a loss of activity was seen with DHT in the L701A AR inciting a near equal response as CPA in this mutant.
Transcriptional Activation of CPA in the WT and M560A GRGiven the high degree of homology, we aligned the AR with the GR and found that the Met-560 corresponds to the Leu-701 residue (Fig. 2d). Based on the position of the 17 Transcriptional Activation of CPA in the WT and L715A PRSimilarly, we constructed the L715A mutant in the PR, because the Leu-715 corresponds to Leu-701 in the AR. CPA is known to act as an agonist for the WT PR similar to progesterone and has therefore clinical use in hormonal contraception. We found that CPA in fact does elicit transcriptional activation in the murine mammary tumor virus-reporter gene assay at a concentration as low as 1 nM similar to progesterone (Fig. 4c). However, the maximal effect seen with CPA was substantially less than that observed with progesterone. The decrease in activity seen in the L715A mutant with progesterone was extreme with complete loss of transcription up to 100 nM and only a minor effect at 1000 nM. CPA, on the other hand, retained its activity with this mutant (Fig. 4d). Comparison of Binding Pockets Observed in Various AR LBD StructuresRecently, we reported crystal structures of the AR LBD bound to S-1 in the WT, T877A, and W741L mutants (1, 8). Each of these AR variants accommodated an identical binding conformation of S-1, which included an expansion of the binding pocket as compared with the AR when complexed to steroidal agonists such as DHT and R1881 (16, 23). Although the binding conformation that S-1 and structurally related compounds adopted in the AR was poorly predicted through molecular modeling (11, 13), the x-ray crystal structure of the AR·S-1 complex demonstrated the induced fit of S-1 through displacement of Trp-741 (8). The repositioning of this residue resulted in expansion of the binding pocket (Fig. 5, a and b) for direct ligand interactions with residues in the ligand-dependent activation function (AF-2), including Gln-738 and Ile-898 and a basis for under-standing the structure-activity relationships observed with structural analogs of S-1 (15, 24-28). Interestingly, the crystal structure of the T877A·CPA complex demonstrated expansion of the AR binding pocket in the opposite direction as compared with S-1. The induced fit of CPA in the AR from displacement of Leu-701 thus opens another region for which drug interactions can be exploited. In particular, the backbone carbonyl of Ser-778 appears available for hydrogen bonding (Fig. 5c).
AR Docking StudiesWe used AutoDock3.05 to examine docking solutions obtained to various AR LBD crystal structures. DHT was removed and docked into both the WT (PDB code 1I37 [PDB] ) and T877A (PDB code 1I38 [PDB] ) AR LBD with low r.m.s.d. values from the bound conformation of DHT observed in the crystal structures. Hydrogen bonds to Arg-752, Asn-705, and Thr-877 in the WT and Arg-752 and Asn-705 in the T877A mutant were also present with DHT, further validating the consistency of our docking methods with crystallographic data. In addition, these solutions exhibited low free binding energies (Table 2) with an appropriate increase observed in the T877A docking run due to loss of this hydrogen bond partner. Conversely, the binding energies obtained from docking CPA to the WT AR (PDB code 1I37 [PDB] ) and T877A (PDB code 1I38 [PDB] ) structures that were co-crystallized with DHT (16) were much higher likely due to the poor ability of these structures to accommodate CPA. The solution with the most favorable free energy of binding showed conservation of the hydrogen bonds to Arg-752 and Asn-705 with the 3-keto group and 17 -acetyl group, respectively, but CPA was rotated out of the steroidal plane relative to the x-ray structure to accommodate the acetate between Leu-701 and Phe-876. The docking solution with the lowest free energy of binding to both the WT and T877A coordinate files demonstrated the same conformation for CPA, yet a substantial improvement in binding energy to the T877A structure was observed likely due to the lack of a steric interaction with the Thr-877 side chain. CPA was docked to the T877A AR structure reported herein with a significantly lower predicted binding energy as compared with the DHT-complexed structures, and an orientation almost identical to that observed in the x-ray crystal structure as demonstrated by the r.m.s.d. of only 0.45 Å. The T877A·CPA structure is thus useful for predicting binding conformations of ligands that occupy this expanded cavity and studying interactions in this otherwise occluded region of the AR in silico.
Herein we presented the crystal structure of CPA bound to the T877A AR LBD. Based on the structural findings observed in the T877A·CPA complex, we predicted that the L701A AR mutant would increase agonist activity of CPA due to the relief of the strain on helix 3 and recovery of the packing of the loop between helices 11 and 12 into the groove observed in other AR LBD structures. As predicted, CPA stimulated transcriptional activation in the L701A AR similar to T877A AR. Conversely, HF exhibited no response in the L701A AR, demonstrating the ability of this mutation to selectively confer agonist activity to CPA. This result further suggests that the gain in AR function observed in the L701A AR is due to the interaction with CPA, because HF does not alter the position of Leu-701 in the T877A AR LBD (8) relative to the published WT·DHT structures (16, 21). Electron density maps for AR residues 886-890 in the T877A·CPA structure demonstrated that this region is not well ordered. The flexibility of this loop, however, provides insight to understanding the dynamics of the AR LBD. Unlike other published AR LBD complexes, the T877A·CPA structure elucidates a gateway for the ligand to enter the binding pocket. Crystal structures of the estrogen receptor LBD have shown that the dynamic properties of helix 12 act as a lid that provides an entrance for ligands to reach the binding cavity (9). However, the C terminus of the AR contains a -sheet that may restrict the mobility of helix 12 relative to the crystal structures of the estrogen receptor , which were determined with a truncated C terminus. Furthermore, evidence of the alteration in the folding of C terminus of helix 11 and adjacent loop region in the T877A·CPA structure suggests that this dynamic portion of the AR may provide ligand accessibility to the binding pocket. As seen from Fig. 2 (b and c), it can be anticipated that CPA binding in the presence of the Thr-877 side chain would push helix 11 away from the binding pocket and thus intensify the unfolding of the AR LBD that already occurs in this region. Such an alteration to the AR LBD structure would resemble the change induced by the ligand, GC-24, in the human thyroid hormone receptor (29). Transcriptional activation data suggest that both the T877A and L701A AR relieve a sufficient amount of strain when complexed to CPA to recover functional AF-2 formation as seen from the ability of CPA to induce transcriptional activation at a low concentration (i.e. 10 nM). Conversely, bicalutamide likely unfolds the AR LBD at the N terminus of helix 12 from steric interaction with Met-895 as shown from our previous report in which W741L and M895T convert it to an agonist by reducing the bulk in this region of the AR (8).
Although CPA acts as an antagonist in the GR similar to the AR, it is an agonist ligand for the PR. The corresponding residue to Leu-701 in the GR and PR were mutated to decrease bulk in the region presumably occupied by the 17 The expanded binding pocket seen in the AR T877A·CPA complex demonstrates an expanded binding pocket not seen in other reported AR LBD crystal structures to date. This finding explains the discrepancy observed in docking experiments and the incompatibility of CPA with the ligand-binding pocket in the AR from the WT·DHT and T877A·DHT complexes (16). The receptor interactions in this cavity provide information useful for structure-based drug design. Ligands that hydrogen bond to Ser-778 will likely offer a means of optimizing binding affinities of new compounds that target the AR. Unfortunately, receptor flexibility is difficult to predict with molecular modeling. Docking software such as FlexX and AutoDock employ rigid protein structures, which limits the ability to accurately predict binding conformations of AR ligands (10, 15). Here we demonstrate the use of ligand-receptor co-crystallography to determine flexible regions of the AR and identify new insight for structure-based drug design and receptor function.
The atomic coordinates and structure factors (code 2Oz7) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).
* This work was supported by National Institutes of Health Grants R01 DK59800 and R01 DK065227. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. 1 To whom correspondence should be addressed: Dept. of Pharmaceutics, Ohio State University, 500 W. 12th Ave., Columbus, OH 43210. Tel.: 614-688-3797; Fax: 614-292-7766; E-mail: dalton.1{at}osu.edu.
2 The abbreviations used are: AR, androgen receptor; CPA, cyproterone acetate; DHT, dihydrotestosterone; GR, estrogen receptor; HF, hydroxyflutamide; LBD, ligand-binding domain; PR, progesterone receptor; r.m.s.d., root mean square distance; S-1, S-3-(4-fluorophenoxy)-2-hydroxy-2-methyl-N-[4-nitro-3-(trifluoromethyl)phenyl] propionamide; GST, glutathione S-transferase; DTT, dithiothreitol; CHAPS, 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonic acid; WT, wild type; CMV, cytomegalovirus.
3 W. L. DeLano (2002) The PyMOL Molecular Graphics System, DeLano Scientific, San Carlos, CA.
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||