Structural Mechanism of N-Methyl-D-Aspartate Receptor Type 1 Partial Agonism

N-methyl-D-aspartate (NMDA) receptors belong to a family of ionotropic glutamate receptors that contribute to the signal transmission in the central nervous system. NMDA receptors are heterotetramers that usually consist of two GluN1 and GluN2 monomers. The extracellular ligand-binding domain (LBD) of a monomer is comprised of discontinuous segments that form the functional domains D1 and D2. While the binding of a full agonist glycine to LBD of GluN1 is linked to cleft closure and subsequent ion-channel opening, partial agonists are known to activate the receptor only sub-maximally. Although the crystal structures of the LBD of related GluA2 receptor explain the mechanism for the partial agonism, structures of GluN1-LBD cannot distinguish the difference between full and partial agonists. It is, however, probable that the partial agonists of GluN1 alter the structure of the LBD in order to result in a different pharmacological response than seen with full agonists. In this study, we used molecular dynamics simulations to reveal an intermediate closure-stage for GluN1, which is unseen in crystal structures. According to our calculations, this intermediate closure is not a transient stage but an energetically stable conformation. Our results demonstrate that the partial agonist cannot exert firm GluN1-LBD closure, especially if there is even a small force that disrupts the LBD closure. Accordingly, this result suggests the importance of forces from the ion channel for the relationship between pharmacological response and the structure of the LBD of members of this receptor family.     

The loss of an electrostatic contact unique to AMPA receptor ligand binding domain 2 slows channel activation

Ligand-gated ion channels undergo conformational changes that transfer the energy of agonist binding to channel opening. Within ionotropic glutamate receptor (iGluR) subunits, this process is initiated in their bilobate ligand binding domain (LBD) where agonist binding to lobe 1 favors closure of lobe 2 around the agonist and allows formation of interlobe hydrogen bonds. AMPA receptors (GluAs) differ from other iGluRs because glutamate binding causes an aspartate-serine peptide bond in a flexible part of lobe 2 to rotate 180° (flipped conformation), allowing these residues to form cross-cleft H-bonds with tyrosine and glycine in lobe 1. This aspartate also contacts the side chain of a lysine residue in the hydrophobic core of lobe 2 by a salt bridge. We investigated how the peptide flip and electrostatic contact (D655-K660) in GluA3 contribute to receptor function by examining pharmacological and structural properties with an antagonist (CNQX), a partial agonist (kainate), and two full agonists (glutamate and quisqualate) in the wildtype and two mutant receptors. Alanine substitution decreased the agonist potency of GluA3(i)-D655A and GluA3(i)-K660A receptor channels expressed in HEK293 cells and differentially affected agonist binding affinity for isolated LBDs without changing CNQX affinity. Correlations observed in the crystal structures of the mutant LBDs included the loss of the D655-K660 electrostatic contact, agonist-dependent differences in lobe 1 and lobe 2 closure, and unflipped D(A)655-S656 bonds. Glutamate-stimulated activation was slower for both mutants, suggesting that efficient energy transfer of agonist binding within the LBD of AMPA receptors requires an intact tether between the flexible peptide flip domain and the rigid hydrophobic core of lobe 2.     

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.     

Microsecond-to-millisecond conformational dynamics demarcate the GluR2 glutamate receptor bound to agonists glutamate, quisqualate, and AMPA

Chemical shift changes and internal motions on microsecond-to-millisecond time scales of the S1S2 ligand-binding domain of the GluR2 ionotropic glutamate receptor have been studied by NMR spectroscopy in the presence of the agonists glutamic acid (glutamate), quisqualic acid (quisqualate), and alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA). Although the crystal structures of the three agonist-bound forms of GluR2 S1S2 ligand-binding domain are very similar, chemical shift changes imply that AMPA-bound GluR2 S1S2 is conformationally distinct from glutamate- and quisqualate-bound forms of GluR2 S1S2. NMR spin relaxation measurements for backbone amide (15)N nuclei reveal that GluR2 S1S2 exhibits reduced chemical exchange line broadening, resulting from microsecond-to-millisecond conformational dynamics, in AMPA-bound compared to glutamate- and quisqualate-bound states. The largest changes in line broadening are observed for two regions of GluR2 S1S2: Val683 and the segment around Lys716-Cys718. The differences in binding affinity of these agonists do not explain the differences in microsecond-to-millisecond conformational dynamics because quisqualate and AMPA bind with similar affinities that are 10-fold greater than the affinity of glutamate. Differences in conformational mobility may reflect differences in the binding mode of AMPA in the GluR2 S1S2 active site compared to the other two ligands. The sites of conformational mobility in GluR2 S1S2 imply that subtle differences exist between the agonists glutamate, quisqualate, and AMPA in modulating glutamate receptor function.     

Luminescence resonance energy transfer investigation of conformational changes in the ligand binding domain of a kainate receptor

The apo state structure of the isolated ligand binding domain of the GluR6 subunit and the conformational changes induced by agonist binding to this protein have been investigated by luminescence resonance energy transfer (LRET) measurements. The LRET-based distances show that agonist binding induces cleft closure, and the extent of cleft closure is proportional to the extent of activation over a wide range of activations, thus establishing that the cleft closure conformational change is one of the mechanisms by which the agonist mediates receptor activation. The LRET distances also provide insight into the apo state structure, for which there is currently no crystal structure available. The distance change between the glutamate-bound state and the apo state is similar to that observed between the glutamate-bound and antagonist UBP-310-bound form of the GluR5 ligand binding domain, indicating that the cleft for the apo state of the GluR6 ligand binding domain should be similar to the UBP-310-bound form of GluR5. This observation implies that te apo state of GluR6 undergoes a cleft closure of 29-30 degrees upon binding full agonists, one of the largest observed in the glutamate receptor family.     

Dynamics of Cleft Closure of the GluA2 Ligand-binding Domain in the Presence of Full and Partial Agonists Revealed by Hydrogen-Deuterium Exchange

The majority of excitatory neurotransmission in the CNS is mediated by tetrameric AMPA receptors. Channel activation begins with a series of interactions with an agonist that binds to the cleft between the two lobes of the ligand-binding domain of each subunit. Binding leads to a series of conformational transitions, including the closure of the two lobes of the binding domain around the ligand, culminating in ion channel opening. Although a great deal has been learned from crystal structures, determining the molecular details of channel activation, deactivation, and desensitization requires measures of dynamics and stabilities of hydrogen bonds that stabilize cleft closure. The use of hydrogen-deuterium exchange at low pH provides a measure of the variation of stability of specific hydrogen bonds among agonists of different efficacy. Here, we used NMR measurements of hydrogen-deuterium exchange to determine the stability of hydrogen bonds in the GluA2 (AMPA receptor) ligand-binding domain in the presence of several full and partial agonists. The results suggest that the stabilization of hydrogen bonds between the two lobes of the binding domain is weaker for partial than for full agonists, and efficacy is correlated with the stability of these hydrogen bonds. The closure of the lobes around the agonists leads to a destabilization of the hydrogen bonding in another portion of the lobe interface, and removing an electrostatic interaction in Lobe 2 can relieve the strain. These results provide new details of transitions in the binding domain that are associated with channel activation and desensitization.     

Medium chain fatty acids are selective peroxisome proliferator activated receptor (PPAR) γ activators and pan-PPAR partial agonists

Thiazolidinediones (TZDs) act through peroxisome proliferator activated receptor (PPAR) γ to increase insulin sensitivity in type 2 diabetes (T2DM), but deleterious effects of these ligands mean that selective modulators with improved clinical profiles are needed. We obtained a crystal structure of PPARγ ligand binding domain (LBD) and found that the ligand binding pocket (LBP) is occupied by bacterial medium chain fatty acids (MCFAs). We verified that MCFAs (C8-C10) bind the PPARγ LBD in vitro and showed that they are low-potency partial agonists that display assay-specific actions relative to TZDs; they act as very weak partial agonists in transfections with PPARγ LBD, stronger partial agonists with full length PPARγ and exhibit full blockade of PPARγ phosphorylation by cyclin-dependent kinase 5 (cdk5), linked to reversal of adipose tissue insulin resistance. MCFAs that bind PPARγ also antagonize TZD-dependent adipogenesis in vitro. X-ray structure B-factor analysis and molecular dynamics (MD) simulations suggest that MCFAs weakly stabilize C-terminal activation helix (H) 12 relative to TZDs and this effect is highly dependent on chain length. By contrast, MCFAs preferentially stabilize the H2-H3/β-sheet region and the helix (H) 11-H12 loop relative to TZDs and we propose that MCFA assay-specific actions are linked to their unique binding mode and suggest that it may be possible to identify selective PPARγ modulators with useful clinical profiles among natural products.     

Solution X-ray scattering evidence for agonist- and antagonist-induced modulation of cleft closure in a glutamate receptor ligand-binding domain

Agonist-induced conformational changes in the ligand-binding domains (LBD) of glutamate receptor ion channels provide the driving force for molecular rearrangements that mediate channel opening and subsequent desensitization. The resulting regulated transmembrane ion fluxes form the basis for most excitatory neuronal signaling in the brain. Crystallographic analysis of the GluR2 LBD core has revealed a ligand-binding cleft located between two lobes. Channel antagonists stabilize an open cleft, whereas agonists stabilize a closed cleft. The crystal structure of the apo form is similar to the antagonist-bound, open state. To understand the conformational behavior of the LBD in the absence of crystal lattice constraints, and thus better to appreciate the thermodynamic constraints on ligand binding, we have undertaken a solution x-ray scattering study using two different constructs encoding either the core or an extended LBD. In agreement with the GluR2 crystal structures, the LBD is more compact in the presence of agonist than it is in the presence of antagonist. However, the time-averaged conformation of the ligand-free core in solution is intermediate between the open, antagonist-bound state and the closed, agonist-bound state, suggesting a conformational equilibrium. Addition of peptide moieties that connect the core domain to the other functional domains in each channel subunit appears to constrain the conformational equilibrium in favor of the open state.     

Conformational changes at the agonist binding domain of the N-methyl-D-aspartic acid receptor

The conformational changes in the agonist binding domain of the glycine-binding GluN1 and glutamate-binding GluN2A subunits of the N-methyl D-aspartic acid receptor upon binding agonists of varying efficacy have been investigated by luminescence resonance energy transfer (LRET) measurements. The LRET-based distances indicate a cleft closure conformational change at the GluN1 subunit upon binding agonists; however, no significant changes in the cleft closure are observed between partial and full agonists. This is consistent with the previously reported crystal structures for the isolated agonist binding domain of this receptor. Additionally, the LRET-based distances show that the agonist binding domain of the glutamate-binding GluN2A subunit exhibits a graded cleft closure with the extent of cleft closure being proportional to the extent of activation, indicating that the mechanism of activation in this subunit is similar to that of the glutamate binding α-amino-5-methyl-3-hydroxy-4-isoxazole propionate and kainate subtypes of the ionotropic glutamate receptors.     

Mechanism of AMPA receptor activation by partial agonists: disulfide trapping of closed lobe conformations

The mechanism by which agonist binding to an ionotropic glutamate receptor leads to channel opening is a central issue in molecular neurobiology. Partial agonists are useful tools for studying the activation mechanism because they produce full channel activation with lower probability than full agonists. Structural transitions that determine the efficacy of partial agonists can provide information on the trigger that begins the channel-opening process. The ligand-binding domain of AMPA receptors is a bilobed structure, and the closure of the lobes is associated with channel activation. One possibility is that partial agonists sterically block full lobe closure but that partial degrees of closure trigger the channel with a lower probability. Alternatively, full lobe closure may be required for activation, and the stability of the fully closed state could determine efficacy with the fully closed state having a lower stability when bound to partial relative to full agonists. Disulfide-trapping experiments demonstrated that even extremely low efficacy ligands such as 6-cyano-7-nitroquinoxaline-2,3-dione can produce a full lobe closure, presumably with low probability. The results are consistent the hypothesis that the efficacy is determined at least in part by the stability of the state in which the lobes are fully closed.     

Mechanisms for activation and antagonism of an AMPA-sensitive glutamate receptor: crystal structures of the GluR2 ligand binding core

Crystal structures of the GluR2 ligand binding core (S1S2) have been determined in the apo state and in the presence of the antagonist DNQX, the partial agonist kainate, and the full agonists AMPA and glutamate. The domains of the S1S2 ligand binding core are expanded in the apo state and contract upon ligand binding with the extent of domain separation decreasing in the order of apo > DNQX > kainate > glutamate approximately equal to AMPA. These results suggest that agonist-induced domain closure gates the transmembrane channel and the extent of receptor activation depends upon the degree of domain closure. AMPA and glutamate also promote a 180 degrees flip of a trans peptide bond in the ligand binding site. The crystal packing of the ligand binding cores suggests modes for subunit-subunit contact in the intact receptor and mechanisms by which allosteric effectors modulate receptor activity.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

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.     

Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs

G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.     

Structural Dynamics of the Glycine-binding Domain of the N-Methyl-d-Aspartate Receptor

N-Methyl-d-aspartate receptors mediate the slow component of excitatory neurotransmission in the central nervous system. These receptors are obligate heteromers containing glycine- and glutamate-binding subunits. The ligands bind to a bilobed agonist-binding domain of the receptor. Previous x-ray structures of the glycine-binding domain of NMDA receptors showed no significant changes between the partial and full agonist-bound structures. Here we have used single molecule fluorescence resonance energy transfer (smFRET) to investigate the cleft closure conformational states that the glycine-binding domain of the receptor adopts in the presence of the antagonist 5,7-dichlorokynurenic acid (DCKA), the partial agonists 1-amino-1-cyclobutanecarboxylic acid (ACBC) and l-alanine, and full agonists glycine and d-serine. For these studies, we have incorporated the unnatural amino acid p-acetyl-l-phenylalanine for specific labeling of the protein with hydrazide derivatives of fluorophores. The single molecule fluorescence resonance energy transfer data show that the agonist-binding domain can adopt a wide range of cleft closure states with significant overlap in the states occupied by ligands of varying efficacy. The difference lies in the fraction of the protein in a more closed-cleft form, with full agonists having a larger fraction in the closed-cleft form, suggesting that the ability of ligands to select for these states could dictate the extent of activation.     

Ligand based conformational space studies of the μ-opioid receptor

BackgroundG protein-coupled receptors (GPCRs) comprise a family of membrane proteins that can be activated by a variety of external factors. The μ-opioid receptor (MOR), a class A GPCR, is the main target of morphine. Recently, enhanced sampling molecular dynamics simulations of a constitutively active mutant of MOR in its apo form allowed us to capture the novel intermediate states of activation, as well as the active state. This prompted us to apply the same techniques to wild type MOR in complex with ligands, in order to explore their contributions to the receptor conformational changes in the activation process.MethodsMOR was modeled in complex with agonists (morphine, BU72), a partial agonist (naloxone benzoylhydrazone) and an antagonist (naloxone). Replica exchange with solute tempering (REST2) molecular dynamics simulations were carried out for all systems. Trajectory frames were clustered, and the activation state of each cluster was assessed by two different methods.ResultsCluster sizes and activation indices show that while agonists stabilized structures in a higher activation state, the antagonist behaved oppositely. Morphine tends to drive the receptor towards increasing R165-T279 distances, while naloxone tends to increase the NPxxYA motif conformational change.ConclusionsDespite not observing a full transition between inactive and active states, an important conformational change of transmembrane helix 5 was observed and associated with a ligand-driven step of the process.General significanceThe activation process of GPCRs is widely studied but still not fully understood. Here we carried out a step forward in the direction of gaining more details of this process.     

Stereochemistry of glutamate receptor agonist efficacy: engineering a dual-specificity AMPA/kainate receptor

Upon agonist binding, the bilobate ligand-binding domains of the ionotropic glutamate receptors (iGluR) undergo a cleft closure whose magnitude correlates broadly with the efficacy of the agonist. AMPA (alpha-amino-5-methyl-3-hydroxy-4-isoxazolepropionic acid) and kainate are nonphysiological agonists that distinguish between subsets of iGluR. Kainate acts with low efficacy at AMPA receptors. Here we report that the structure-based mutation L651V converts the GluR4 AMPA receptor into a dual-specificity AMPA/kainate receptor fully activated by both agonists. To probe the stereochemical basis of partial agonism, we have also investigated the correlation between agonist efficacy and a series of vibrational and fluorescence spectroscopic signals of agonist binding to the corresponding wild-type and mutant GluR4 ligand-binding domains. Two signals track the extent of channel activation: the maximal change in intrinsic tryptophan fluorescence and the environment of the single non-disulfide bonded C426, which appears to probe the strength of interactions with the ligand alpha-amino group. Both of these signals arise from functional groups that are poised to detect changes in the extent of channel cleft closure and thus provide additional information about the coupling between conformational changes in the ligand-binding domain and activation of the intact receptor.     

Population Shift Mechanism for Partial Agonism of AMPA Receptor

α-amino-3-hydroxy-5-methyl-4-isoaxazolepropionic acid (AMPA) ionotropic glutamate receptors mediate fast excitatory neurotransmission in the central nervous system, and their dysfunction is associated with neurological diseases. Glutamate binding to ligand-binding domains (LBDs) of AMPA receptors induces channel opening in the transmembrane domains of the receptors. The T686A mutation reduces glutamate efficacy so that the glutamate behaves as a partial agonist. The crystal structures of wild-type and mutant LBDs are very similar and cannot account for the observed behavior. To elucidate the molecular mechanism inducing partial agonism of the T686A mutant, we computed the free-energy landscapes governing GluA2 LBD closure using replica-exchange umbrella sampling simulations. A semiclosed state, not observed in crystal structures, appears in the mutant during simulation. In this state, the LBD cleft opens slightly because of breaking of interlobe hydrogen bonds, reducing the efficiency of channel opening. The energy difference between the LBD closed and semiclosed states is small, and transitions between the two states would occur by thermal fluctuations. Evidently, glutamate binding to the T686A mutant induces a population shift from a closed to a semiclosed state, explaining the partial agonism in the AMPA receptor.