| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Protein Engineering, Vol. 14, No. 8, 565-571,
August 2001
© 2001 Oxford University Press
Homology modelling of the ligand-binding domain of glucocorticoid receptor: binding site interactions with cortisol and corticosterone
1 Department of Physics and 3 Department of Chemistry, University of Calcutta, 92 A.P.C. Road, Calcutta-700009, India
| Abstract |
|---|
|
|
|---|
Glucocorticoids are involved in the growth, development and homeostasis of a number of tissues. The physiological effects of this class of lipophilic steroids are mediated by ligand-inducible nuclear transcription factor, the glucocorticoid receptor/mineralocorticoid receptor, a member of the steroid/nuclear receptor superfamily. The glucocorticoid receptor interacts specifically with glucocorticoids, whereas the mineralocorticoid receptor interacts with both glucocorticoids and mineralocorticoids. The molecular structure of progesterone complexed to its receptor obtained from X-ray crystal structure analysis is used to build up a homology model of mouse glucocorticoid receptor ligand-binding domain (mGR LBD). The secondary structure of mGR LBD contains 11 helices, nine turns and four sheets. The mGR LBD contains a long helix, H9, with 30 residues, and exhibits slight deformation when the receptor protein binds with its cognate ligands. The mGR LBD has a 12-residue C-terminal extension (residues 772783) that is essential for hormone binding. This extension is tightly fixed in position by an antiparallel ß-sheet interaction between amino acids 680682 (S3) and 775777 (S4). The three-dimensional model reveals two polar sites located at the extremities of the elongated hydrophobic ligand-binding pocket. Cortisol and corticosterone are docked to this ligand-binding pocket. The difference accessible surface area study revealed the steroid-binding region of mGR LBD.
Keywords: glucocorticoid receptor/homology modelling/ligand-binding domain/steroid receptor complex/transactivation
| Introduction |
|---|
|
|
|---|
The ability of nuclear receptor to activate specific gene transcription requires the binding of cognate ligands to their ligand-binding domains (LBDs) (Wurtz et al., 1996
The DBDs and LBDs for almost all the members of the steroid receptor superfamily are relatively well conserved, whereas the NTDs are less conserved (Hollenberg et al., 1985
; Arriza et al., 1987
). The receptors have overlapping steroid-binding specificity as they bind to the same hormone-responsive elements in target genes (Arriza et al., 1987
; Pearce and Yamamoto, 1993
; Rupprecht et al., 1993b
). Mineralocorticoid receptor binds mineralocorticoids and most glucocorticoids with high affinity, whereas glucocorticoid receptor binds only glucocorticoids with high affinity (Arriza et al., 1987
, 1998; Rupprecht et al., 1993a
,b
). There is a functional similarity in both these classes of steroids, but glucocorticoids specifically regulate the development of specific tissues, glucose metabolism and the immune response (Miller and Tyrrell, 1995
).
A remarkable feature of the glucocorticoids is that they exert pronounced effects on a variety of tissues. In certain sites they stimulate the synthesis of macromolecular components, including nucleic acids and proteins. Sometimes the action of glucocorticoids is catabolic, resulting ultimately in cellular destruction (King, 1974
). This diversity of action is reflected in the contemporary search for specific glucocorticoid-binding proteins in a wide range of experimental systems. Binding to the specific hormone receptors is acutely dependent on the presence of the 3-one A-ring and/or 11ß-hydroxyl group, which is the characteristic structural feature of all biologically active glucocorticoids. The 3-one A-ring of many progestins and corticoids is the area of greatest conformational flexibility in these molecules. This conformational flexibility exhibited by these naturally occurring hormones facilitates their highly specific interaction with a variety of proteins including synthesizing enzymes, transport proteins and target proteins (Duax and Norton, 1975
).
The theoretical model of mGR LBD was developed templating the X-ray structure of progesterone receptor (Williams and Siglar, 1998
). The model receptor protein is then used to study the binding site interactions during complexation with cortisol and corticosterone by molecular modelling. The difference in accessible surface area (DASA) (Lee and Richards, 1971
) between mGR LBD and the ligand-bound protein is then calculated for both the steroidreceptor complexes. The interaction zone of mGR LBD with the steroid was revealed from this study. Binding of steroids to mGR LBD does not necessitate any structural change to the hydrophobic ligand-binding pocket. However, the modelling study did indicate a slight conformational change in the secondary structure on complexation of mGR LBD with both cortisol and corticosterone. This minor conformational change was found in the C-terminal site of the long helix H9 without any gross alteration of the ligand-binding pocket itself. This suggests that ligand binding may trigger a conformational modification which could account for the effect of ligand binding on transactivation by AF-2 (Bourguet et al., 1995
). Perhaps the ability of corticosteroids to stimulate the transactivation function of mGR depends on the stability of the steroidreceptor complexes (Hellal-Levy et al., 1999
).
| Materials and methods |
|---|
|
|
|---|
Starting conformation and sequence alignment
The refined crystal structures of progesterone-bound LBD of the human progesterone receptor (hPR) (Williams and Siglar, 1998
), LBD of estrogen receptor (ER) in complex with the endrogenous estrogen (Brzozowski et al., 1997
) and the LBD of human androgen receptor (hAR) in complex with metribolone (R1881) (Matias et al., 2000
) were taken from the Brookhaven Protein Data Bank (PDB entries 1A28, 1ERRand 1e3g, respectively) as starting materials. The amino acid sequences of three LBDs were extracted from these three crystal structures. The amino acid sequence of house mouse glucocorticoid receptor (mGR) (Danielsen et al., 1986
; Nohno et al., 1989
) was obtained from SWISSPROT Sequence Data Bank and was compared with the sequences of the crystal structures separately by pairwise sequence alignment using the software GAP (Needleman and Wunsch, 1970
) of the GCG package. The best similarity and identity of mGR with hPR LBD encouraged us to select the crystal structure of progesterone-bound LBD of hPR to develop a theoretical model of mGR LBD. A multiple sequence alignment among mGR, hPR, hAR and hER was done using the PILE UP (Feng and Doolittle, 1987
) program of the GCG package.
Coordinate assignment and minimization
The coordinates of mGR LBD were assigned by templating the X-ray structure of hPR LBD after aligning the two sequences as was observed in the GCG output. This coordinate assignment was done using the HOMOLOGY module (Biosym Technologies) of the InsightII program package. The model of mGR LBD was then put through energy minimization for 10 000 steps of the steepest descent method using the DISCOVER module of InsightII. The model was further subjected to energy minimization for 500 steps of the conjugate gradient technique that led to a refined structure of mGR LBD with an r.m.s. deviation of <0.001. The secondary structure prediction of mGR LBD was performed using the program DSSP (Kabsch and Sander, 1983
), July 1995 version. The DISCOVER simulation package (Biosym Technologies) with the consistent valence force-field (Hagler et al., 1985
; Dauber-Osguthorpe et al., 1988
) was employed for minimization calculations.
Superposition and ligand docking
Using the Biopolymer module (Biosym Technologies) of InsightII, the refined model structure of mGR LBD was superposed on the crystal structure of progesterone-bound hPR LBD. The molecular structure of cortisol (Roberts et al., 1973
) was then superposed on progesterone such that their steroid nuclei were practically coincident. The allocated cortisol could then be easily associated with the ligand-binding domain of the superposed glucocorticoid receptor. The hydroxyl group at C-17 of the cortisol in the mGR LBDcortisol complex thus prepared was deleted to mimic an mGR LBDcorticosterone complex. These two complexes were then subjected to energy minimization for 1000 steps of the conjugate gradient technique with the DISCOVER simulation package.
Solvent accessibility
The values of the accessible surface area for both the native protein (mGR LBD) and ligand-bound proteins were calculated using the HOMOLOGY module of InsightII. The differences in accessible surface areas between mGR LBD and ligand-bound protein were calculated for every residue. This DASA study traced the steroid protein interaction regions for the complexes of mGR LBD with both cortisol and corticosterone.
| Results and discussion |
|---|
|
|
|---|
Table I
|
|
|
The majority of the residues of mGR LBD occupy the most favored regions of the Ramachandran plot and the other residues occupy additional allowed regions as defined in Procheck (Laskowski et al., 1993
The postulated model of mGR LBD consists of 250 residues (Leu534Lys783), which folds into a hydrophobic ligand-binding pocket. The secondary structure of this model protein contains 11 helices, nine turns and four sheets. The pairs of helices H3, H4; H5, H6; and H10, H11 are almost contiguous pairwise. mGR LBD contains a relatively longer helix H9 (it contains 30 residues, whereas the PR contains 15, the RAR contains eight, the TR contains six and the ER contains nine). A shortening of the length of the helix H9 is observed in both the complexes of mGR LBD with cortisol and corticosterone. On complexation the helix H9 (718747) of the native protein (mGR LBD) becomes deformed and is shortened by one residue to become helix H9 (718746). This deformation at the C-terminal end of helix H9 is shown in Figure 3
, where the trace of mGR LBDcortisol complex is superposed on the native protein. This superposition (Figure 3
) shows a perfect coincidence everywhere except from Ser740 to Asp748. The r.m.s. deviation in the aligned position is 0.147. The DASA study between mGR LBD and its complex with cortisol is represented by the bar graph in Figure 4(a)
. Figure 4(b)
represents the DASA study between mGR LBD and its corticosterone complex. The positive DASA values reveal the interaction zone of the receptor with its cognate ligands. The negative DASA values are found to be around the region of deformation arises from complexation. Figure 5(a)
represents the C
trace of mGR LBD, where H and S indicate
-helices and ß-sheets, respectively. Bound cortisol is shown in a space filling model. Figure 5(b)
represents the ribbon diagram of mGR LBD where the bound cortisol is shown in the space filling model.
|
|
|
The polar neutral residues Asn570 and Gln576 of mGR LBD are directly involved in binding with cortisol through hydrogen bonding. The C-3-ketone group of the A-ring of cortisol forms a hydrogen bond with the amido NH2 group of Gln576 of helix H2 at site I of the ligand-binding pocket. The required position and orientation of Gln576 are maintained by the supporting role of Arg617 and Phe629, underscored by their conservation at the corresponding sequence position in all steroid receptors. The amido NH2 group of the polar residue Asn570 of helix H2 serves as a donor to O-21 of cortisol and the keto-oxygen of the same residue serves as an acceptor from the hydroxyl group at C-11 of cortisol. The residue Asn570 thereby held stably by forming two hydrogen bonds with the steroid at site II of the ligand-binding pocket. The remaining residues in the binding cavity participate in a number of hydrophobic interactions with the ligand. The mineralocorticoid compound aldosterone bearing an identical 17ß-substituent as in cortisol shows similar binding interactions with MR in site II of the ligand-binding pocket (Fagart et al., 1998
The hydrogen-bonding parameters associated with cortisol and corticosterone in complexation with mGR LBD are given in Table II
. Table III
presents the binding energies of cortisol and corticosterone in their complexes with mGR LBD. Figure 6(a) and (b)![]()
show a stereoscopic view of the binding site interactions of cortisol and corticosterone, respectively, with mGR LBD in a stick diagram. Final energies before complexation of mGR LBD and hPR LBD (X-ray structure) with steroids are given in Table IV
. Table V
represents the final energies of the receptors complexed to steroids. Tables IV and V![]()
clearly reveal a high degree of agreement between the actual crystal structure and model structure, thereby confirming the reliability of the modelled steroidreceptor complex as well as the model of the uncomplexed mGR LBD.
|
|
|
|
|
The overall architechture of this modeled mGR LBD is similar to that seen in the crystal structures of the LBDs of other nuclear receptors such as ER, TR, PR, RXR, AR and RAR (Bourguet et al., 1995
and ß orientational side of the steroid nucleus, respectively. However, the helix H2 is disposed at an angle with the steroid nucleus. The disposition of these three helices around cortisol in complex with mGR LBD is shown in Figure 7
|
| Notes |
|---|
2 To whom correspondence should be addressed. Email: prcxray{at}cubmb.ernet.in and draja63{at}yahoo.co.in
| Acknowledgments |
|---|
The authors thank Professor Asok Banerjee, Bose Institute, Calcutta for providing laboratory facilities. Thanks are due to D.I.C., Bose Institute, Calcutta for providing the computational facilities for the modelling work.
| References |
|---|
|
|
|---|
Alen,P., Claessens,F., Verhoeven,G., Rombauts,W. and Peeters,B. (1999) Mol. Cell. Biol., 19, 60856097.
Arriza,J.L., Weinberger,C., Cerelli,G., Glaser,T.M., Handelin,B.L., Housman,D.E. and Evans,R.M. (1987) Science, 237, 268275.
Arriza,J.L., Simerly,R.B., Swanson,L.W. and Evans,R.M. (1988) Neuron, 1, 887900.[Web of Science][Medline]
Bourguet,W., Ruff,M., Chambon,P., Gronemeyer,H. and Moras,D. (1995) Nature, 375, 377382.[Medline]
Brzozowski,A.M., Pike,A.C., Dauter,Z., Hubbard,R.E., Bonn,T., Engstrom,O., Ohman,L., Greene,G.L., Gustafsson,J.-A. and Carlquist,M. (1997) Nature, 389, 753758.[Medline]
Carlstedt-Duke,J., Stromstedt,P.-E., Wrange,O., Bergman,T., Gustafsson,J.-A. and Jornvall,H. (1987) Proc. Natl Acad. Sci. USA, 84, 44374440.
Danielsen,M., Northrop,J.P. and Ringold,G.M. (1986) EMBO J., 5, 25132522.[Web of Science][Medline]
Dauber-Osguthorpe,P., Roberts,V.A., Osguthorpe,D.J., Wolff,J., Genest,M. and Hagler,A.T. (1988) Proteins: Struct. Funct. Genet., 4, 3147.[Web of Science][Medline]
Duax,W.L. and Norton,D.A. (1975) Atlas of Steroid Structure, Vol. 1. Plenum Press, New York, p. 29.
Evans,R.M. (1988) Science, 240, 889895.
Fagart,J., Wurtz,J.M., Souque,A., Hellal-Levy,C., Moras,D. and Rafestin-Oblin,M.E. (1998) EMBO J., 17, 33173325.[Web of Science][Medline]
Feng,D.F. and Doolittle,R.F. (1987) J. Mol. Evol., 25, 351360.[Web of Science][Medline]
Glass,C.K. (1994) Endocrinol. Rev., 15, 391407.
Hagler,A.T. (1985) In Hurby,V.J. and Meienhofer,J. (eds), The Peptides. Academic Press, New York, pp. 213299.
He,B., Kemppainen,J.A., Gronemeyer,H. and Wilson,E.M. (1999) J. Biol. Chem., 274, 3721937225.
Hellal-Levy,C., Couette,B., Fagart,J., Souque,A., Gomez-Sanchez,C. and Rafestin-Oblin,M.E. (1999) FEBS Lett., 464, 913.[Web of Science][Medline]
Hollenberg,S.M., Weinberger,C., Ong,E.S., Cerelli,G., Oro,A., Lebo,R., Thompson,E.B., Rosenfeld,M.G. and Evans,R.M. (1985) Nature, 318, 635641.[Medline]
Jenster,G. et al. (1991) Mol. Endocrinol., 5, 13961404.
Kabsch,W. and Sander,C. (1983) Biopolymers, 22, 25772637.[Web of Science][Medline]
King,R.J.B. (1974) SteroidCell Interaction. Butterworths, London, p. 109.
Lanz,R. and Rusconi,S. (1996) Endocrinology, 135, 21832195.[Abstract]
Laskowski,R., MacArthur,M., Moss,D. and Thornton,J. (1993) J. Appl. Crystallogr., 26, 9197.
Laudet,V., Hanni,C., Coll,J., Catzeflis,F. and Stehelin,D. (1992) EMBO J., 11, 10031013.[Web of Science][Medline]
Lee,B. and Richards,F.M. (1971) J. Mol. Biol., 55, 379400.[Web of Science][Medline]
Lind,U., Greenidge,P., Gustafsson,J.-A., Wright,A.P.H. and Carlstedt-Duke,J. (1999) J. Biol. Chem., 274, 1851518523.
Matias,P.M. et al. (2000) J. Biol. Chem., 275, 2616426171.
Miller,W.L. and Tyrrell,J.B. (1995) In Felig,P., Baxter,J.D. and Frohman,L.D. (eds), Endocrinology and Metabolism, 3rd edn., McGraw-Hill, New York.
Needleman and Wunsch (1970) J. Mol. Biol., 48, 443453.[Web of Science][Medline]
Nohno,T., Kasai,Y. and Saito,T. (1989) Nucleic Acids Res., 17, 445.
Onate,S.A., Tsai,S.Y., Tsai,M.J. and O'Malley,B.W. (1995) Science, 270, 13541357.
Parker,M.J. and White,R. (1996) Nature Struct. Biol., 3, 113115.[Web of Science][Medline]
Pearce,D. and Yamamoto,K.R. (1993) Science, 259, 11611165.
Renaud,J.-P., Rochel,N., Ruff,M., Vivat,V., Chambon,P., Gronemeyer,H. and Moras,D. (1995) Nature, 378, 681689.[Medline]
Ribeiro,R.J.C., Kushner,P.J. and Baxter,J.D. (1995) Annu. Rev. Med., 46, 443453.[Web of Science][Medline]
Roberts,R.J., Coppola,J.C., Isaacs,N.W. and Kennard,O. (1973) J. Chem., Soc., Perkin Trans. 2, 774781.
Rupprecht,R., Arriza,J.L., Spengler,D., Reul,J.M., Evans,R.M., Holsboer,F. and Damm,K. (1993a) Mol. Endocrinol., 7, 597603.
Rupprecht,R., Reul,J.M., van Steensel,B., Spengler,D., Soder,M., Berning,B., Holsboer,F. and Damm,K. (1993b) Eur. J. Pharmacol., 247, 145154.[Web of Science][Medline]
Torchia,J., Glass,C. and Rosenfeld,M.G. (1998) Curr. Opin. Cell. Biol., 10, 373383.[Web of Science][Medline]
Tsai,M.-J. and O'Malley,B.W. (1994) Annu. Rev. Biochem., 63, 451486.[Web of Science][Medline]
Voegel,J.J., Heine,M.J.S., Tini,M., Vivat,V., Chambon,P. and Gronemeyer,H. (1998) EMBO J., 17, 507519.[Web of Science][Medline]
Wagner,R.L., Apriletti,J.W., Mc Grath,M.E., West,B.L., Baxter,J.D. and Fletterick,R.J. (1995) Nature, 378, 690697.[Medline]
Williams,S.P. and Siglar,P.B. (1998) Nature, 393, 392396.[Medline]
Wurtz,J.M., Bourguet,W., Renaud,J.P., Vivat,V., Chambon,P., Moras,D. and Gronemeyer,H. (1996) Nature Struct. Biol., 3, 8794.[Web of Science][Medline]
Xu,J., Nawaz,Z., Tsai,S.Y., Tsai,M.J. and O'Malley,B.W. (1996) Proc. Natl Acad. Sci. USA, 93, 1219512199.
Zhang,S., Liang,X. and Danielsen,M. (1996) Mol. Endocrinol., 10, 2434.
Received September 29, 2000; revised April 27, 2001; accepted April 30, 2001.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
K. V. Kepple, N. Patel, P. Salamon, and A. M. Segall Interactions between branched DNAs and peptide inhibitors of DNA repair Nucleic Acids Res., September 1, 2008; 36(16): 5319 - 5334. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Yoshikawa, K. Yamamoto, N. Shimizu, S. Yamada, C. Morimoto, and H. Tanaka The Distinct Agonistic Properties of the Phenylpyrazolosteroid Cortivazol Reveal Interdomain Communication within the Glucocorticoid Receptor Mol. Endocrinol., May 1, 2005; 19(5): 1110 - 1124. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Zhang, J. Simisky, F. T. F. Tsai, and D. S. Geller A critical role of helix 3-helix 5 interaction in steroid hormone receptor function PNAS, February 22, 2005; 102(8): 2707 - 2712. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Kumar and E. B. Thompson Transactivation Functions of the N-Terminal Domains of Nuclear Hormone Receptors: Protein Folding and Coactivator Interactions Mol. Endocrinol., January 1, 2003; 17(1): 1 - 10. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Yoshikawa, Y. Makino, K. Okamoto, C. Morimoto, I. Makino, and H. Tanaka Distinct Interaction of Cortivazol with the Ligand Binding Domain Confers Glucocorticoid Receptor Specificity. CORTIVAZOL IS A SPECIFIC LIGAND FOR THE GLUCOCORTICOID RECEPTOR J. Biol. Chem., February 8, 2002; 277(7): 5529 - 5540. [Abstract] [Full Text] [PDF] |
||||













