Structural insights into ligand-binding pocket formation in Nurr1 by molecular dynamics simulations

Abstract

The nuclear receptor Nurr1 (NR4A2) has been identified as a potential target for the treatment of Parkinson’s disease. In contrast to most other nuclear receptors, the X-ray crystal structure of the Nurr1 ligand-binding domain (LBD) lacks any ligand-binding pocket (LBP). However, NMR spectroscopy measurements have revealed that the known Nurr1 agonist docosahexaenoic acid (DHA) binds to a region within the LBD that corresponds to the classical NR ligand-binding pocket (LBP). In order to investigate the structural dynamics of the Nurr1 LBD and to study potential LBP formation, the conformational space of the receptor was sampled using a molecular dynamics (MD) simulation. Docking of DHA into 50,000 LBD structures extracted from the simulation revealed the existence of a transient LBP that is capable to fully harbor the compound. The location of the identified pocket overlaps with the ligand-binding site suggested by NMR experiments. Structural analysis of the protein-ligand complex showed that only modest structural rearrangements within the Nurr1 LBD are required for LBP formation. These findings may support structure-based drug discovery campaigns for the development of receptor-specific agonists.     

Hydrogen/deuterium exchange reveals distinct agonist/partial agonist receptor dynamics within vitamin D receptor/retinoid X receptor heterodimer

Regulation of nuclear receptor (NR) activity is driven by alterations in the conformational dynamics of the receptor upon ligand binding. Previously, we demonstrated that hydrogen/deuterium exchange (HDX) can be applied to determine novel mechanism of action of PPARγ ligands and in predicting tissue specificity of selective estrogen receptor modulators. Here, we applied HDX to probe the conformational dynamics of the ligand binding domain (LBD) of the vitamin D receptor (VDR) upon binding its natural ligand 1α,25-dihydroxyvitamin D3 (1,25D3), and two analogs, alfacalcidol and ED-71. Comparison of HDX profiles from ligands in complex with the LBD with full-length receptor bound to its cognate receptor retinoid X receptor (RXR) revealed unique receptor dynamics that could not be inferred from static crystal structures. These results demonstrate that ligands modulate the dynamics of the heterodimer interface as well as provide insight into the role of AF-2 dynamics in the action of VDR partial agonists.     

核雌激素受体配体结合域的晶体结构特性与功能

雌激素或类雌激素活性物质通过细胞核雌激素受体（nuclear estrogen receptor, nER）通路发挥相应的生理性作用。当这些配体被nER的配体结合域（ligand binding domain, LBD）识别后进入疏水性配体结合空腔内并引起受体构象发生改变,使得原先处于高度活动性的helix 12(H12)被固定从而进一步稳定空腔结构|同时nER也能通过招募一系列辅助调节因子及其他共调节蛋白质,最终调控基因转录。但是,由于不同的配体和受体结合形成的晶体结构并不完全相同,导致这些复合体具有不同的性质,从而影响基因的转录活性。本文综述了nER配体结合域及结合配体后形成的相应晶体结构与活性以及不同配体对受体结构和基因转录的影响。     

Hormone selectivity in thyroid hormone receptors

Separate genes encode thyroid hormone receptor subtypes TRalpha (NR1A1) and TRbeta (NR1A2). Products from each of these contribute to hormone action, but the subtypes differ in tissue distribution and physiological response. Compounds that discriminate between these subtypes in vivo may be useful in treating important medical problems such as obesity and hypercholesterolemia. We previously determined the crystal structure of the rat (r) TRalpha ligand-binding domain (LBD). In the present study, we determined the crystal structure of the rTRalpha LBD in a complex with an additional ligand, Triac (3,5, 3'-triiodothyroacetic acid), and two crystal structures of the human (h) TRbeta receptor LBD in a complex with either Triac or a TRbeta-selective compound, GC-1 [3,5-dimethyl-4-(4'-hydroy-3'-isopropylbenzyl)-phenoxy acetic acid]. The rTRalpha and hTRbeta LBDs show close structural similarity. However, the hTRbeta structures extend into the DNA-binding domain and allow definition of a structural "hinge" region of only three amino acids. The two TR subtypes differ in the loop between helices 1 and 3, which could affect both ligand recognition and the effects of ligand in binding coactivators and corepressors. The two subtypes also differ in a single amino acid residue in the hormone-binding pocket, Asn (TRbeta) for Ser (TRalpha). Studies here with TRs in which the subtype-specific residue is exchanged suggest that most of the selectivity in binding derives from this amino acid difference. The flexibility of the polar region in the TRbeta receptor, combined with differential recognition of the chemical group at the 1-carbon position, seems to stabilize the complex with GC-1 and contribute to its beta-selectivity. These results suggest a strategy for development of subtype-specific compounds involving modifications of the ligand at the 1-position.     

Dynamic stabilization of nuclear receptor ligand binding domains by hormone or corepressor binding

We have developed a novel assembly assay to examine structural changes in the ligand binding domain (LBD) of the thyroid hormone receptor (TR). Fragments including the first helix of the TR LBD interact only weakly with the remainder of the LBD in the absence of hormone, but this interaction is strongly enhanced by the addition of either hormone or the corepressor NCoR. Since neither the ligand nor the corepressor shows direct interaction with this helix, we propose that both exert their effects by stabilizing the overall structure of the LBD. Current models of activation of nuclear hormone receptors focus on a ligand-induced allosteric shift in the position of the C-terminal helix 12 that generates the coactivator binding site. Our results suggest that ligand binding also has more global effects that dynamically alter the structure of the receptor LBD.     

New insights into receptor ligand binding domains from a novel assembly assay

Previous studies have demonstrated that hormone binding stabilizes the ligand binding domain (LBD) of the nuclear hormone receptors against proteolysis. We have confirmed and extended this observation using a newly developed assembly assay. In this assay, the LBD is divided into two parts, of which one includes the first helix of this domain and the other corresponds to the remainder of the LBD. Several independent criteria demonstrate that these two fragments can assemble into a functional LBD in the presence of a ligand, but not in its absence, and that this is a reflection of the stabilizing effect of ligand. We have also used this assay to demonstrate that binding of the nuclear receptor corepressor NCoR can directly stabilize the LBD. Overall, these results highlight the dynamic nature of the LBD and suggest that current models for activation based solely on allosteric effects on the C-terminal helix may be too limited.     

Evolutionary selection across the nuclear hormone receptor superfamily with a focus on the NR1I subfamily (vitamin D,pregnane X,and constitutive androstane receptors)

Background

The nuclear hormone receptor (NR) superfamily complement in humans is composed of 48 genes with diverse roles in metabolic homeostasis, development, and detoxification. In general, NRs are strongly conserved between vertebrate species, and few examples of molecular adaptation (positive selection) within this superfamily have been demonstrated. Previous studies utilizing two-species comparisons reveal strong purifying (negative) selection of most NR genes, with two possible exceptions being the ligand-binding domains (LBDs) of the pregnane X receptor (PXR, NR1I2) and the constitutive androstane receptor (CAR, NR1I3), two proteins involved in the regulation of toxic compound metabolism and elimination. The aim of this study was to apply detailed phylogenetic analysis using maximum likelihood methods to the entire complement of genes in the vertebrate NR superfamily. Analyses were carried out both across all vertebrates and limited to mammals and also separately for the two major domains of NRs, the DNA-binding domain (DBD) and LBD, in addition to the full-length sequences. Additional functional data is also reported for activation of PXR and the vitamin D receptor (VDR; NR1I1) to gain further insight into the evolution of the NR1I subfamily.

Results

The NR genes appear to be subject to strong purifying selection, particularly in the DBDs. Estimates of the ratio of the non-synonymous to synonymous nucleotide substitution rates (the ω ratio) revealed that only the PXR LBD had a sub-population of codons with an estimated ω ratio greater than 1. CAR was also unusual in showing high relative ω ratios in both the DBD and LBD, a finding that may relate to the recent appearance of the CAR gene (presumably by duplication of a pre-mammalian PXR gene) just prior to the evolution of mammals. Functional analyses of the NR1I subfamily show that human and zebrafish PXRs show similar activation by steroid hormones and early bile salts, properties not shared by sea lamprey, mouse, or human VDRs, or by Xenopus laevis PXRs.

Conclusion

NR genes generally show strong sequence conservation and little evidence for positive selection. The main exceptions are PXR and CAR, genes that may have adapted to cross-species differences in toxic compound exposure.

     

Dynamics and allostery of the ionotropic glutamate receptors and the ligand binding domain

The dynamics of the ligand‐binding domain (LBD) and the intact ionotropic glutamate receptor (iGluR) were studied using Gaussian Network Model (GNM) analysis. The dynamics of LBDs with various allosteric modulators is compared using a novel method of multiple alignment of GNM modes of motion. The analysis reveals that allosteric effectors change the dynamics of amino acids at the upper lobe interface of the LBD dimer as well as at the hinge region between the upper‐ and lower‐ lobes. For the intact glutamate receptor the analysis show that the clamshell‐like movement of the LBD upper and lower lobes is coupled to the bending of the trans‐membrane domain (TMD) helices which may open the channel pore. The results offer a new insight on the mechanism of action of allosteric modulators on the iGluR and support the notion of TMD helices bending as a possible mechanism for channel opening. In addition, the study validates the methodology of multiple GNM modes alignment as a useful tool to study allosteric effect and its relation to proteins dynamics. Proteins 2016; 84:267–277. © 2015 Wiley Periodicals, Inc.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.