![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Pharmacology, University of Medicine & Dentistry of New Jersey (UMDNJ), Robert Wood Johnson Medical School and UMDNJ Informatics Institute, Piscataway, New Jersey
Received January 11, 2006; accepted February 1, 2006
| Abstract |
|---|
|
|
|---|
AF or AF-2 helix in other nuclear receptors is responsible for transcription activation by recruiting coactivators through conformational change. In the absence of ligands, PXR can also repress gene expression by interacting with transcriptional corepressors, such as the silencing mediator for retinoid and thyroid hormone receptor (SMRT). We first constructed homology models of the complete LBD with two SMRT nuclear receptor (NR)-interacting domains (ID1 and ID2), respectively. We then performed energy minimization and molecular dynamics simulations on these systems to study the specific interactions between the interacting domains and LBD. Further experimental results supported and validated the observed preference of SMRT toward ID2 over ID1. Our modeling results revealed the key interactions that account for the binding preference. Here, we propose structural models of the PXR-LBD/SMRT-ID1 and PXR-LBD/SMRT-ID2 complexes to understand their molecular interactions and potential inhibitory mechanism.
AF is responsible for transcription activation by recruiting coactivators through conformational rearrangement (Leo and Chen, 2000
SMRT is a general corepressor for many nuclear receptors, such as thyroid receptor and peroxisome proliferators-activated receptor
(PPAR
). It has two nuclear receptor (NR)-interacting domains, named ID1 and ID2 (Li et al., 1997
). Each ID has a signature, CoRNR box or LxxxIxxxI/L motif, a corepressor motif responsible for interactions with NRs (Hu and Lazar, 1999
; Nagy et al., 1999
; Perissi et al., 1999
). The two IDs can interact with a single NR or with a NR dimer, respectively, and different NRs have different binding preferences. In a recent report, Johnson et al. (2006
) have demonstrated that PXR binds preferably to ID2 over ID1. Interpretation of the molecular modeling results in this report provides clues to the key residues involved in the binding preference.
|
| Materials and Methods |
|---|
|
|
|---|
Homology Modeling. Numerous X-ray crystal structures of PXR are available for use as a template (PDB codes 1ILG
[PDB]
, 1ILH, 1M13, 1NRL, 1SKX) to build our homology models; however, none of them contain an intact LBD. In particular, residues 178 to 209 form a highly flexible loop, which is unresolved in these crystal structures. This loop is purported to be involved in PXRs promiscuity with respect to the diversity of its ligands. The crystal structure of PXR in apo-form (PDB code 1ILG
[PDB]
) was chosen to be the initial template, because this structure represents the most suitable conformation of PXR LBD in complex with its corepressor. The flexible loop of 15 amino residues, although not resolved in the X-ray crystal structure, was generated by using Loopy module of Jackal (http://trantor.bioc.columbia.edu/programs/jackal), a protein structure modeling package. The final model for PXR/SMRT-ID2 complex was built from the X-ray crystal structure of SMRT bound to PPAR
and GW6471 (PDB code 1KKQ) (Xu et al., 2002
), and the intact LBD of PXR. The PPAR
/SMRT-ID2 structure was first aligned and superimposed onto the PXR-LBD structure to minimize the root-mean-square deviation (RMSD) of backbone atoms in the core of the LBD using Sybyl 6.9 (http://www.tripos.com). This superimposition minimized the RMSD between the core PXR helices and the corresponding PPAR
helices and effectively superimposed their corepressor binding sites. The coordinates and amino acids of SMRT-ID2 peptide were transferred from the PPAR
structure to the PXR-LBD, and then examined visually. A few residue side chains of SMRT were rotated to improve its fit with the
AF helix of PXR using the program IN-SIGHT II (http://www.accelrys.com). The final PXR/SMRT-ID2 structural model was then subjected to further EM and MD simulations. For the PXR/SMRT-ID1 complex, we mutated the ID2 sequence to the corresponding ID1 sequence while preserving the remaining structural components. EM and MD simulations were employed for each construct to explore protein-protein interactions in the binary complex (PXR-LBD/SMRT) under normal physiological conditions. Two structural models were generated and refined: PXR-LBD/SMRT-ID2 and PXR-LBD/SMRT-ID1 (Fig. 1, b and c).
Energy Minimization and Molecular Dynamics Simulations. Refinement of the PXR/SMRT-ID1 and -ID2 structural models was conducted by EM and MD calculations using AMBER7 (http://amber.scripps.edu/). After EM, PXR-LBD in complex with SMRT-ID1 was solvated by a 9-Å radius shell of transferable intermolecular potentials (TIP3) water molecules (Ryckaert et al., 1977
). The solvated system was energy minimized in two steps: first, 1000 iterations of constrained steepest descent (SD) whereby only the water molecules were free to move to eliminate improper steric interactions; second, 5000 interations of SD and 5000 iterations of conjugated gradient minimization was conducted on the entire system. MD simulations were then performed on the resulting system using the standard force-field parameter set parm99 in AMBER7 with dielectric constant
= 1 and cutoff distance = 9.0 Å applied for both electrostatic and van der Waals interactions. The SHAKE algorithm (Ryckaert et al., 1977
) was implemented for bonds involving hydrogen atoms and the time step = 2 fs. The system was then coupled to a Berendsen bath at 300 K by using a coupling constant
T = 2 ps (Berendsen et al., 1984
). The system was gradually increased in temperature from 100 to 300 K over 10 ps of simulation time with the volume held constant (ensemble NVT). When the system approached T
300 K and the density
1.0 g/ml, constant pressure and temperature controls (NPT) were applied to the simulation was conducted for 1 ns. The pressure of the system was raised to 1 atm (ensemble NPT). The same force-field parameters were then employed for EM and MD simulations of the PXR/SMRT-ID2 system using the same protocol as for PXR/SMRT-ID1. Structures retrieved from the equilibrated (post-500 ps) system were averaged and energy-minimized using 1000 iterations of SD and 5000 iterations of conjugated gradient. The final refined structures were analyzed with PROCHECK (Morris et al., 1992
; Laskowski et al., 1993
).
Experimental Validation Using GST Pull-Down Assay. The human pregnane X receptor (GenBank accession no. NR1I2), human retinoid acid receptor-
(NR1B1), and human thyroid receptor-
(NR1A2) were created by standard molecular cloning techniques, including restriction digestion with endonucleases (New England BioLabs), and PCR using Pfu polymerase (Stratagene, La Jolla, CA), followed by end-joining ligation with T4 DNA ligase (Roche Molecular Biochemicals, Indianapolis, IN). Derivatives of SMRT (Chen and Evans, 1995
; Park et al., 1999
) were created similarly. All constructs used for in vitro transcription/translation were in the pCMX vector (Umesono et al., 1991
). Glutathione S-transferase (GST) fusion protein purification was conducted as described previously (Frangioni and Neel, 1993
) with glutathione-conjugated agarose beads (Sigma). GST pull-down assay was conducted as described previously (Ghosh et al., 2002
).
|
|
| Results and Discussion |
|---|
|
|
|---|
AF helix of PXR, SMRT-ID1, SMRT-ID2, and the AF-2 of human retinoid acid (RAR-AF2) (Fig. 1a). We observed a similarity between the PXR
AF helix and the corepressor LxxIxxxI/L motif. In contrast, the RAR AF-2 helix is more related to the coactivator LxxLL motif (Fig. 1a). In general, a nuclear receptor recruits a coactivator or corepressor through a combination of hydrophobic interactions by the canonical motifs and other protein-protein interactions by the "charged clamps" and/or hydrogen bonds. The sequence similarity of the PXR
AF helix with the corepressor motif is further supported by the fact that the helix contains a two-turn
-helix, similar to the proposed structure of the corepressor motif. The primary amino acid sequences and the secondary structure of the helix support the ability of the
AF helix to mediate or inhibit PXR binding to SMRT. To further assess the ability of PXR to recruit SMRT, we first constructed homology models of PXR with various the SMRT NR-interacting domains (ID1 and ID2). We subsequently performed EM and MD simulations to study the specific interactions between the two IDs and the PXR LBD. We further validated the structural models and the binding of ID1 and ID2 to PXR using GST pull-down assays.
Numerous X-ray crystal structures for use as a structural template are not intact and lack of a highly flexible loop of residues 178 to 209, which are unresolved in these crystal structures. The crystal structure of PXR in apo-form (PDB ID = 1ILG [PDB] ) was chosen to be the initial template. This unresolved loop of 15 amino residues was then generated by using Loopy module of Jackal (http://trantor.bioc.columbia.edu/programs/jackal), a protein structure-modeling package. The structural model for the PXR-LBD/SMRT-ID2 complex was built from the X-ray crystal structure of SMRT bound to PPAR and GW6471 (PDB code 1KKQ), and the X-ray crystal structure of the apo-form PXR (PDB code 1ILG [PDB] ). For the PXR-LBD/SMRT-ID1 complex, we mutated the ID2 sequence to ID1 sequence and kept the rest of structure complex intact. EM and MD simulations were employed for each construct. Two structural models were generated and refined: PXR-LBD/SMRT-ID1 and PXR-LBD/SMRT-ID2.
MD simulations on the PXR-LBD/SMRT-ID1 structure confirm the stability of three-dimensional structure and the equilibrium of binary complex, as evidenced by the rapid convergence of the RMSD of the
-carbon backbone of the sampled conformation compared with the initial structure (Fig. 2, straight line). Likewise, RMSD analysis of the MD trajectory (Fig. 2, dashed line) of the PXR-LBD/SMRT-ID2 structural model indicated a stable 3D structure. PRO-CHECK assessment also displayed acceptable quality, aside from a few amino acid residues falling outside the "allowed" regions associated with the flexible loop mentioned above (residues 178-209).
Based on the assumption that the corepressor site overlaps with the coactivator site, SMRT-ID2 and the PXR LBD were shown to form three key interactions: one charge clamp and two salt bridges (Fig. 3b). The essential charge clamp typically formed by Lys259 of helix 3 (Watkins et al., 2003
) and coactivator was conserved for both ID1 and ID2 domains of SMRT. In particular, Lys259 formed a charge clamp with Ile2143 of ID1 and Leu2350 of ID2, respectively (Fig. 3). Another essential charge clamp typically formed by Glu427 of the
AF helix (Watkins et al., 2003
) was not found in either domain. The preference of SMRT for ID2 over ID1 can be attributed to two additional key interactions: Glu427 of
AF with Lys2348 of ID2 and Glu270 of helix 4 with Arg2347 of ID2 (Fig. 3b). These two additional interactions are salt bridges and, as such, are stronger than a standard hydrogen bond. Our structural model suggested that Pro423 and Lys2348 of ID2 formed another hydrogen bond; however, this interaction is presumably much weaker than the two salt bridges.
|
AF helix. Furthermore, ID1 is anchored with the
AF helix through only two hydrogen bonds: one between Glu427 of
AF and Gln2137 of ID1, and the other between Pro423 of
AF and His2129 of ID1 (Fig. 3a). The preceding discussion underscores the difference in ascribing the differential binding among ID1 and ID2 to a solitary factor. Indeed, ID1 and ID2 bear great similarity in their primary amino acid sequence and secondary structure, and both share the canonical corepressor motif. Our computational approaches enabled us to elucidate putative key interactions that account for the observed binding difference in binding of these two similar interacting domains to SMRT. The differences in the type and number of nonbonded interactions between the two domains of SMRT and PXR LBD revealed by our structural analysis may explain the specific preference for the recruitment of SMRT-ID2 by PXR.
The results from the MD studies are in accordance with the experimental results. Biological experimental evidence using GST-pull down assay showed the differential preferences of three nuclear receptors toward individual SMRT domains (Fig. 4). Human RAR
in general binds preferentially to ID1, whereas human TR
and human PXR bind preferentially to SMRT-ID2. Further biological studies, including immunofluorescence colocalization and luciferase reporter assays, are detailed in the recent publication by Johnson et al. (2006
) Taken together, our computational and experimental results provide new insights into the assembly of the PXR-LBD/SMRT-ID2 complex by proposing the key interactions that account for the binding preferences and possible molecular mechanism of
AF helix-mediated inhibition of SMRT binding. In future work, we will perform molecular simulations and biological experiments to assess the binding of PXR-
AF to NR-domains (ID1 and ID2) of SMRT. These future studies will complement ongoing experimental mutational analysis of the specific residues (e.g., Phe423 and Glu427) identified in the
AF helix to validate the present modeling observations.
| Conclusions |
|---|
|
|
|---|
In this report, we have provided a potential molecular model of interactions that determine the assembly of the PXR/SMRT-ID1 and -ID2 complexes. Our modeling results reveal the structural basis for the binding preference of SMRT-ID2 over SMRT-ID1 to the PXR
AF helix. Because there are very few known antagonists of PXR, the selectivity of NR-interacting domains can be further employed to regulate the gene expression mediated by PXR. These studies provide insights into the molecular interactions between PXR and SMRT corepressor, which may have important implications for understanding the role of corepressors in regulating the manifold biological activities of PXR.
| Acknowledgements |
|---|
| Footnotes |
|---|
ABBREVIATIONS: PXR, pregnane X receptor; SXR, steroid and xenobiotic receptor; SMRT, silencing mediator for retinoid and thyroid hormone receptor; NR, nuclear receptor; EM, energy minimization; MD, molecular dynamics; LBD, ligand binding domain; RMSD, root-mean-square deviation; PPAR, peroxisome proliferator-activated receptor; SD, steepest descent; GST, glutathione S-transferase; ID, nuclear receptor interacting domain; RAR, retinoic acid receptor; GW6471, N-((2S)-2-((1Z)-1-methyl-3-oxo-3-(4-(trifluoromethyl)phenyl)prop-1-enyl)amino)-3-(4-(2-(5-methyl-2-phenyl-1,3-oxazol-4-yl)ethoxy)phenyl)propyl)propanamide;
AF, activation function 2; PDB, Protein Data Bank.
Address correspondence to: Dr. William J. Welsh, Department of Pharmacology, University of Medicine and Dentistry of New Jersey, Robert Wood Johnson Medical School, 675 Hoes Lane, Piscataway, NJ 08854. E-mail: welshwj{at}umdnj.edu
| References |
|---|
|
|
|---|
Blumberg B, Sabbagh W, Juguilon H, Bolado JJ, van Meter CM, Ong ES, and Evans RM (1998) SXR, a novel steroid and xenobiotic sensing nuclear receptor. Genes Dev 12: 3195-3205.
Chen JD and Evans RM (1995) A transcriptional co-repressor that interacts with nuclear hormone receptors. Nature (Lond) 377: 454-457.[CrossRef][Medline]
Frangioni JV and Neel BG (1993) Solubilization and purification of enzymatically active glutathione S-transferase (pGEX) fusion proteins. Anal Biochem 210: 179-187.[CrossRef][Medline]
Ghosh JC, Yang X, Zhang A, Lambert MH, Li H, Xu HE, and Chen JD (2002) Interactions that determine the assembly of a retinoid X receptor/corepressor complex. Proc Natl Acad Sci USA 99: 5842-5847.
Horlein AJ, Naar AM, Heinzel T, Torchia J, Gloss B, Kurokawa R, Ryan A, Kamei Y, Soderstrom M, Glass CK, et al. (1995) Ligand-independent repression by the thyroid hormone receptor mediated by a nuclear receptor co-repressor. Nature (Lond) 377: 397-404.[CrossRef][Medline]
Hu X and Lazar MA (1999) The CoRNR motif controls the recruitment of corepressors by nuclear hormone receptors. Nature (Lond) 402: 93-96.[CrossRef][Medline]
Johnson DR, Li C-W, Chen L-Y, Ghosh JC, and Chen JD (2006) Regulation and binding of pregnane X receptor by nuclear receptor corepressor silencing mediator of retinoid and thyroid hormone receptors (SMRT). Mol Pharmacol 69: 99-108.
Kliewer SA, Moore JT, Wade L, Staudinger JL, Watson MA, Jones SA, McKee DD, Oliver BB, Willson TM, and Zetterstrom RH (1998) An orphan nuclear receptor activated by pregnanes defines a novel steroid signaling pathway. Cell 92: 73-82.[CrossRef][Medline]
Laskowski RA, MacArthur MW, Moss DS, and Thornton JM (1993) PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crys 26: 283-291.[CrossRef]
Leo C and Chen JD ((2000)) The SRC family of nuclear receptor coactivators. Gene 245: 1-11.[CrossRef][Medline]
Li H, Leo C, Schroen DJ, and Chen JD (1997) Characterization of Receptor Interaction and Transcriptional Repression by the Corepressor SMRT. Mol Endocrinol 11: 2025-2037.
Morris AL, MacArthur MW, Hutchinson EG, and Thornton JM (1992) Stereochemical quality of protein structure coordinates. Proteins 12: 345-364.[CrossRef][Medline]
Nagy L, Kao H-Y, Love JD, Li C, Banayo E, Gooch JT, Krishna V, Chatterjee K, Evans RM, and Schwabe JWR (1999) Mechanism of corepressor binding and release from nuclear hormone receptors. Genes Dev 13: 3209-3216.
Ordentlich P, Downes M, Xie W, Genin A, Spinner NB, and Evans RM (1999) Unique forms of human and mouse nuclear receptor corepressor SMRT. Proc Natl Acad Sci USA 96: 2639-2644.
Park EJ, Schroen DJ, Yang M, Li H, Li L, and Chen JD (1999) SMRTe, a silencing mediator for retinoid and thyroid hormone receptors-extended isoform that is more related to the nuclear receptor corepressor. Proc Natl Acad Sci USA 96: 2639-2644.
Perissi V, Staszewski LM, McInerney EM, Kurokawa R, Krones A, Rose DW, Lambert MH, Milburn MV, Glass CK, and Rosenfeld MG (1999) Molecular determinants of nuclear receptor-corepressor interaction. Genes Dev 13: 3198-3208.
Ryckaert JP, Ciccotti G, and Berendsen HJC (1977) Numerical integration of Cartesian equation of motion of a system with constraints: molecular dynamics of N-alkanes. J Comput Phys 23: 327-341.[CrossRef]
Synold TW, Dussault I, and Forman BM (2001) The orphan nuclear receptor SXR coordinately regulates drug metabolism and efflux. Nat Med 7: 584-590.[CrossRef][Medline]
Tabb MM, VK, Grun F, Zhou C, Welsh WJ, and Blumberg B (2004) Highly chlorinated PCBs inhibit the human xenobiotic response mediated by the steroid and xenobiotic receptor (SXR). Environ Health Perspect 112: 163-169.[Medline]
Umesono K, Murakami KK, Thompson CC, and Evans RM (1991) Direct repeats as selective response elements for the thyroid hormone, retinoic acid and vitamin D3 receptors. Cell 65: 1255-1266.[CrossRef][Medline]
Watkins RE, Davis-Searles PR, Lambert MH, and Redinbo MR (2003) Coactivator binding promotes the specific interaction between ligand and the pregnane X receptor. J Mol Biol 331: 815-828.[CrossRef][Medline]
Watkins RE, Wisely GB, Moore LB, Collins JL, Lambert MH, Williams SP, Wilson TM, Kliewer SA, and Redinbo MR (2001) The human nuclear xenobiotic receptor PXR: Structural determinants of directed promiscuity. Science (Wash DC) 292: 2329-2333.
Wilson TM and Kliewer SA (2002) PXR, CAR and metabolism. Nat Rev Drug Discov 1: 259-266.[CrossRef][Medline]
Xu HE, Stanley TB, Montana VG, Lambert MH, Shearer BG, Cobb JE, McKee DD, Galardi CM, Plunket KD, Nolte RT, et al. (2002) Structural basis for antagonist-mediated recruitment of nuclear co-repressors by PPAR[alpha]. Nature (Lond) 415: 813-817.[Medline]
This article has been cited by other articles:
![]() |
S. Ekins, V. Kholodovych, N. Ai, M. Sinz, J. Gal, L. Gera, W. J. Welsh, K. Bachmann, and S. Mani Computational Discovery of Novel Low Micromolar Human Pregnane X Receptor Antagonists Mol. Pharmacol., September 1, 2008; 74(3): 662 - 672. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ekins, C. Chang, S. Mani, M. D. Krasowski, E. J. Reschly, M. Iyer, V. Kholodovych, N. Ai, W. J. Welsh, M. Sinz, et al. Human Pregnane X Receptor Antagonists and Agonists Define Molecular Requirements for Different Binding Sites Mol. Pharmacol., September 1, 2007; 72(3): 592 - 603. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||