Molecular mechanism of ligand recognition by NR3 subtype glutamate receptors

NR3 subtype glutamate receptors have a unique developmental expression profile, but are the least well-characterized members of the NMDA receptor gene family, which have key roles in synaptic plasticity and brain development. Using ligand binding assays, crystallographic analysis, and all atom MD simulations, we investigate mechanisms underlying the binding by NR3A and NR3B of glycine and D-serine, which are candidate neurotransmitters for NMDA receptors containing NR3 subunits. The ligand binding domains of both NR3 subunits adopt a similar extent of domain closure as found in the corresponding NR1 complexes, but have a unique loop 1 structure distinct from that in all other glutamate receptor ion channels. Within their ligand binding pockets, NR3A and NR3B have strikingly different hydrogen bonding networks and solvent structures from those found in NR1, and fail to undergo a conformational rearrangement observed in NR1 upon binding the partial agonist ACPC. MD simulations revealed numerous interdomain contacts, which stabilize the agonist-bound closed-cleft conformation, and a novel twisting motion for the loop 1 helix that is unique in NR3 subunits.     

Molecular dynamics simulations reveal multiple pathways of ligand dissociation from thyroid hormone receptors

Nuclear receptor (NR) ligands occupy a pocket that lies within the core of the NR ligand-binding domain (LBD), and most NR LBDs lack obvious entry/exit routes upon the protein surface. Thus, significant NR conformational rearrangements must accompany ligand binding and release. The precise nature of these processes, however, remains poorly understood. Here, we utilize locally enhanced sampling (LES) molecular dynamics computer simulations to predict molecular motions of x-ray structures of thyroid hormone receptor (TR) LBDs and determine events that permit ligand escape. We find that the natural ligand 3,5,3'-triiodo-L-thyronine (T(3)) dissociates from the TRalpha1 LBD along three competing pathways generated through i), opening of helix (H) 12; ii), separation of H8 and H11 and the Omega-loop between H2 and H3; and iii), opening of H2 and H3, and the intervening beta-strand. Similar pathways are involved in dissociation of T(3) and the TRbeta-selective ligand GC24 from TRbeta; the TR agonist IH5 from the alpha- and beta-TR forms; and Triac from two natural human TRbeta mutants, A317T and A234T, but are detected with different frequencies in simulations performed with the different structures. Path I was previously suggested to represent a major pathway for NR ligand dissociation. We propose here that Paths II and III are also likely ligand escape routes for TRs and other NRs. We also propose that different escape paths are preferred in different situations, implying that it will be possible to design NR ligands that only associate stably with their cognate receptors in specific cellular contexts.     

Structural and Functional Similarity of Amphibian Constitutive Androstane Receptor with Mammalian Pregnane X Receptor

The nuclear receptors and xenosensors constitutive androstane receptor (CAR, NR1I3) and pregnane X receptor (PXR, NR1I2) induce the expression of xenobiotic metabolizing enzymes and transporters, which also affects various endobiotics. While human and mouse CAR feature a high basal activity and low induction upon ligand exposure, we recently identified two constitutive androstane receptors in Xenopus laevis (xlCARα and β) that possess PXR-like characteristics such as low basal activity and activation in response to structurally diverse compounds. Using a set of complementary computational and biochemical approaches we provide evidence for xlCARα being the structural and functional counterpart of mammalian PXR. A three-dimensional model of the xlCARα ligand-binding domain (LBD) reveals a human PXR-like L-shaped ligand binding pocket with a larger volume than the binding pockets in human and murine CAR. The shape and amino acid composition of the ligand-binding pocket of xlCAR suggests PXR-like binding of chemically diverse ligands which was confirmed by biochemical methods. Similarly to PXR, xlCARα possesses a flexible helix 11’. Modest increase in the recruitment of coactivator PGC-1α may contribute to the enhanced basal activity of three gain-of-function xlCARα mutants humanizing key LBD amino acid residues. xlCARα and PXR appear to constitute an example of convergent evolution.     

In Silico Adoption of an Orphan Nuclear Receptor NR4A1

A 4.1μs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the ‘traditional’ nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 μs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.     

Structural Consequences of D481N/K483Q Mutation at Glycine Binding Site of NMDA Ionotropic Glutamate Receptors: A Molecular Dynamics Study

Abstract

N-Methyl-D-Aspartate (NMDA) receptors are the ligand gated as well as voltage sensitive ionotropic glutamate receptors, widely distributed in the vertebrate central nervous system and they play critical role in the pathogenesis of schizophrenia. Molecular dynamics simulations have been carried out on high resolution crystal structure of NR1 subunit of NMDA receptor ligand binding core (S1S2) in four different conformations. We have investigated consequence of D481N/K483Q double mutation of NR1 subunit from simulation results of (a) glycine bound form (WG), (b) unbound (closed-apo) form (WOG), (c) a double mutated form (DM), and (d) the antagonist (5,7-dichlorokynuric acid) bound form (DCKA). The MD simulations and simulated annealing for 4ns show a distinct conformation for the double mutated conformation that neither follows the antagonist nor apo conformation. There are two distinct sites, loop1 and loop2 where the double mutated structure in its glycine bound form shows significant RMSD deviations as compared to the wild-type. The interactions of glycine with the receptor remain theoretically unchanged in the double mutated structure and there is no detachment of S1S2 domains. The results suggest that separation of S1 and S2 domains may not be essential for channel inactivation. Therefore, it is hypothesized that hypoactivation of NMDA receptor channels may arise out of the conformational changes at non-conserved Loop1 and Loop2 regions observed in the mutated structure. The Loop1 and Loop2 regions responsible for inter-subunit interactions in a functional NMDA receptor, may therefore, render the ligand bound form defunct. This may account for behavioral anomalies due to receptor inactivation seen in grin1 mutated mice.     

Structural consequences of D481N/K483Q mutation at glycine binding site of NMDA ionotropic glutamate receptors: a molecular dynamics study

N-Methyl-D-Aspartate (NMDA) receptors are the ligand gated as well as voltage sensitive ionotropic glutamate receptors, widely distributed in the vertebrate central nervous system and they play critical role in the pathogenesis of schizophrenia. Molecular dynamics simulations have been carried out on high resolution crystal structure of NR1 subunit of NMDA receptor ligand binding core (S1S2) in four different conformations. We have investigated consequence of D481N/K483Q double mutation of NR1 subunit from simulation results of (a) glycine bound form (WG), (b) unbound (closed-apo) form (WOG), (c) a double mutated form (DM), and (d) the antagonist (5,7-dichlorokynuric acid) bound form (DCKA). The MD simulations and simulated annealing for 4ns show a distinct conformation for the double mutated conformation that neither follows the antagonist nor apo conformation. There are two distinct sites, loop1 and loop2 where the double mutated structure in its glycine bound form shows significant RMSD deviations as compared to the wild-type. The interactions of glycine with the receptor remain theoretically unchanged in the double mutated structure and there is no detachment of S1S2 domains. The results suggest that separation of S1 and S2 domains may not be essential for channel inactivation. Therefore, it is hypothesized that hypoactivation of NMDA receptor channels may arise out of the conformational changes at non-conserved Loop1 and Loop2 regions observed in the mutated structure. The Loop1 and Loop2 regions responsible for inter-subunit interactions in a functional NMDA receptor, may therefore, render the ligand bound form defunct. This may account for behavioral anomalies due to receptor inactivation seen in grin1 mutated mice.     

Modeling of peptides connecting the ligand‐binding and transmembrane domains in the GluR2 glutamate receptor

Ligand‐gated Glutamate receptors (GluR) mediate synaptic signals in the nervous system. Ionotropic GluRs of AMPA type, the subject of this study, are tetrameric assemblies of monomer subunits, each of which is constructed in a modular fashion from functional subdomains. The extracellular ligand‐binding domain (LBD) changes its conformation upon binding of an agonist ligand followed by opening of a transmembrane (TM) ion channel. Peptides connecting the LBD and TM domains facilitate gating of the channel, and their structure and composition are important for the receptor functioning. In this study, we used replica exchange molecular dynamics (REMD) simulations to model S1M1 and S2M3 connecting peptides of the GluR2 receptor in two implicit solvents, water and interfacial water/lipid medium characterized by lower polarity. Propensity of these peptides to form helical structures was analyzed using helicity measure derived from the free energy of the simulated ensembles of structures. The S1M1 and S2M3 connecting peptides were not helical in our simulations in both dielectric environments in the absence of the rest of the protein. The structures of the LBD fragment with known high‐resolution α‐helical structure and of the TM3 helix were successfully predicted in the simulations, which in part validate our results. The S2M3 peptide, which is important in gating, formed a well‐defined coil structure and salt‐bridges with the S2 domain. The S1M1 peptide formed a loop structure via formation of internal salt‐bridges. Potential implications of these structures on function of the receptor are discussed. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Helix 1/8 interactions influence the activity of nuclear receptor ligand-binding domains

The ligand-binding domain (LBD) of apo-nuclear receptors in solution is thought to be a very dynamic structure with many possible conformations. Upon ligand binding, the structure is stabilized to a more rigid conformation. The dynamic stabilization assay is a LBD reassembly assay that takes advantage of the high specificity of the intramolecular interactions that comprise the ligand-bound LBD. Here, we demonstrate dynamic stabilization for the nuclear receptors peroxisome proliferator-activated receptor (PPAR)gamma and nerve growth factor inducible (NGFIB)beta and identify residues important for stabilization of the intramolecular interactions induced by PPARgamma ligands. Site-directed mutagenesis studies identified residues in helices 1 and 8 required for LBD reassembly. Further, disrupting the helix 1/8 interaction in the context of the holo-LBD alters the response of the receptor in a compound-specific manner, suggesting that residues far from the ligand-binding pocket can influence the stability of the ligand-bound receptor. Thus, these results support and extend models of the apo-LBD of PPARgamma as a dynamic structure.     

Achieving Peptide Binding Specificity and Promiscuity by Loops: Case of the Forkhead-Associated Domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.     

Complete integrin headpiece opening in eight steps

Carefully soaking crystals with Arg-Gly-Asp (RGD) peptides, we captured eight distinct RGD-bound conformations of the α_IIbβ₃ integrin headpiece. Starting from the closed βI domain conformation, we saw six intermediate βI conformations and finally the fully open βI with the hybrid domain swung out in the crystal lattice. The β1-α1 backbone that hydrogen bonds to the Asp side chain of RGD was the first element to move followed by adjacent to metal ion-dependent adhesion site Ca²⁺, α1 helix, α1’ helix, β6-α7 loop, α7 helix, and hybrid domain. We define in atomic detail how conformational change was transmitted over long distances in integrins, 40 Å from the ligand binding site to the opposite end of the βI domain and 80 Å to the far end of the hybrid domain. During these movements, RGD slid in its binding groove toward α_IIb, and its Arg side chain became ordered. RGD concentration requirements in soaking suggested a >200-fold higher affinity after opening. The thermodynamic cycle shows how higher affinity pays the energetic cost of opening.