Molecular dynamics simulations of the ligand-binding domain of the ionotropic glutamate receptor GluR2

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.

Comparative analysis of different competitive antagonists interaction with NR2A and NR2B subunits of N-methyl-D-aspartate (NMDA) ionotropic glutamate receptor

The antagonist-bound conformation of the NR2A and NR2B subunits of N-methyl-D-aspartate (NMDA) ionotropic glutamate receptor are modeled using the crystal structure of the DCKA (5,7-dichloro-kynurenic acid)-bound form of the NR1 subunit ligand-binding core (S1S2). Five different competitive NMDA receptor antagonists [(1) DL-AP5; (2) DL-AP7; (3) CGP-37847; (4) CGP 39551; (5) (RS)-CPP] have been docked into both NR2A and NR2B subunits. Experimental studies report NR2A and NR2B subunits having dissimilar interactions and affinities towards the antagonists. However, the molecular mechanism of this difference remains unexplored. The distinctive features in the antagonist's interaction with these two different but closely related (approximately 80% sequence identity at this region) subunits are analyzed from the patterns of their hydrogen bonding. The regions directly involved in the antagonist binding have been classified into seven different interaction sites. Two conserved hydrophilic pockets located at both the S1 and S2 domains are found to be crucial for antagonist binding. The positively charged (Lys) residues present at the second interaction site and the invariant residue (Arg) located at the fourth interaction site are seen to influence ligand binding. The geometry of the binding pockets of NR2A and NR2B subunits have been determined from the distance between the C-alpha atoms in the residues interacting with the ligands. The binding pockets are found to be different for NR2A and NR2B. There are gross dissimilarities in competitive antagonist binding between these two subunits. The binding pocket geometry identified in this study may have the potential for future development of selective antagonists for the NR2A or NR2B subunit.

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.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Structural dynamics of an ionotropic glutamate receptor

Ionotropic glutamate receptors (iGluRs) are postsynaptic ion channels involved in excitatory neurotransmission. iGluRs play important roles in development and in forms of synaptic plasticity that underlie higher order processes such as learning and memory. Neurobiological and biochemical studies have long characterized iGluRs in detail. However, the structural basis for the function of iGluRs has not yet been investigated, because there is insufficient information about their three-dimensional structures. In 1998, a crystal structure called S1S2 lobes was first solved for the extracellular bilobed ligand-binding domain of the GluR2 subunit. Since then, the crystal structures for the S1S2 lobes both in the apo and in various liganded states have been reported, and recent biophysical studies have further elucidated the dynamic aspects of the structure of the S1S2 lobes. In this review, the dynamic structures of the S1S2 lobes and their ligands are summarized, and the importance of their structural flexibility and fluctuation is discussed in light of the mechanisms of ligand recognition, activation, and desensitization of the receptor.

Nucleotide-dependent conformational changes in HisP: molecular dynamics simulations of an ABC transporter nucleotide-binding domain

ATP-binding cassette (ABC) transporters mediate the movement of molecules across cell membranes in both prokaryotes and eukaryotes. In ABC transporters, solute translocation occurs after ATP is either bound or hydrolyzed at the intracellular nucleotide-binding domains (NBDs). Molecular dynamics (MD) simulations have been employed to study the interactions of nucleotide with NBD. The results of extended (approximately 20 ns) MD simulations of HisP (total simulation time approximately 80 ns), the NBD of the histidine transporter HisQMP2J from Salmonella typhimurium, are presented. Analysis of the MD trajectories reveals conformational changes within HisP that are dependent on the presence of ATP in the binding pocket of the protein, and are sensitive to the presence/absence of Mg ions bound to the ATP. These changes are predominantly confined to the alpha-helical subdomain of HisP. Specifically there is a rotation of three alpha-helices within the subdomain, and a movement of the signature sequence toward the bound nucleotide. In addition, there is considerable conformational flexibility in a conserved glutamine-containing loop, which is situated at the interface between the alpha-helical subdomain and the F1-like subdomain. These results support the mechanism for ATP-induced conformational transitions derived from the crystal structures of other NBDs.

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.

Towards a three-alpha-helix bundle protein that binds volatile general anesthetics

The general anesthetics halothane and chloroform are capable of binding to synthetic water-soluble four-alpha-helix bundles, which model the putative in vivo receptors. In this study, we investigate the binding of these anesthetics to synthetic water-soluble three-alpha-helix bundles. A series of variants containing up to four X-to-Ala and up to four X-to-Met substitutions was made; and the effect of these substitutions on structure, stability and anesthetic binding affinity was examined. Generally, the amount of alpha-helix and the stability of the three-alpha-helix bundles decreased as the number of X-to-Ala substitutions increased. A concomitant red-shift in tryptophan fluorescence lambdamax was seen, suggesting an increased flexibility of the native structure. Up to four X-to-Met substitutions had little effect on the amount of alpha-helix, but an increase in tryptophan lambdamax was seen for the variants with three and four methionine substitutions. The exceptions were a) a variant with a clustering of alanine and methionine residues at one end of the three-alpha-helix bundle, suggesting a gate structure that can admit ligand molecules; and b) a variant with a single Leu35Ala substitution, suggesting that at select positions, the size of the side chain is important for defining anesthetic binding affinity.

Evolutionary trace analysis of ionotropic glutamate receptor sequences and modeling the interactions of agonists with different NMDA receptor subunits

The ionotropic N-methyl- d-aspartate (NMDA) receptor is of importance in neuronal development, functioning, and degeneration in the mammalian central nervous system. The functional NMDA receptor is a heterotetramer comprising two NR1 and two NR2 or NR3 subunits. We have carried out evolutionary trace (ET) analysis of forty ionotropic glutamate receptor (IGRs) sequences to identify and characterize the residues forming the binding socket. We have also modeled the ligand binding core (S1S2) of NMDA receptor subunits using the recently available crystal structure of NR1 subunit ligand binding core which shares approximately 40% homology with other NMDA receptor subunits. A short molecular dynamics simulation of the glycine-bound form of wild-type and double-mutated (D481N; K483Q) NR1 subunit structure shows considerable RMSD at the hinge region of S1S2 segment, where pore forming transmembrane helices are located in the native receptor. It is suggested that the disruption of domain closure could affect ion-channel activation and thereby lead to perturbations in normal animal behavior. In conclusion, we identified the amino acids that form the ligand-binding pocket in many ionotropic glutamate receptors and studied their hydrogen bonded and nonbonded interaction patterns. Finally, the disruption in the S1S2 domain conformation (of NR1 subunit- crystal structure) has been studied with a short molecular dynamics simulation and correlated with some experimental observations. [figure: see text]. The figure shows the binding mechanism of glutamate with NR2B subunit of the NMDA receptor. Glutamate is shown in cpk, hydrogen bonds in dotted lines and amino acids in blue. The amino acids shown here are within a 4-A radius of the ligand (glutamate).

Ligand-binding residues integrate affinity and efficacy in the NMDA receptor

The interaction of an agonist with its receptor can be characterized by two fundamental properties, affinity and efficacy. Affinity defines how tightly the agonist associates with its receptor, and efficacy measures the ability of the bound ligand to activate the receptor. Although affinity and efficacy are independent properties, the binding and activation processes that they describe are tightly coupled. This strong coupling has complicated the interpretation of concentration-response phenotypes caused by receptor mutations. We present an approach that quantifies the role of individual amino acids in defining affinity and efficacy. This method, which employs partial agonists and covalent modification of introduced cysteines, was applied to the ligand-binding sites of the NMDA receptor. Recent crystallographic structures for glutamate receptor ligand-binding cores allowed identification of residues that are either known or are predicted to be critical for ligand binding in the NR1 and NR2A subunit, respectively. Mutation of amino acids whose sidechains would directly coordinate bound ligands affected both agonist affinity and efficacy. In contrast, positions predicted to stabilize the closed-cleft conformation contributed only to agonist efficacy. The results provide a molecular basis for the tight coupling of agonist binding and receptor activation.

Emerging structural explanations of ionotropic glutamate receptor function

High-resolution studies of ionotropic glutamate receptor (iGluR) extracellular domains are beginning to bridge the gap between structure and function. Crystal structures have defined the ligand binding pocket well beyond what was suggested by mutational analysis and homology models alone, providing initial suggestions about the mechanisms of channel gating and desensitization. NMR-derived backbone dynamics and molecular dynamics simulations have added further insights into the role of protein dynamics in receptor function. As a whole, the current knowledge of iGluR structure in conjunction with new advances in the understanding of K+ channels provides a vastly improved understanding of iGluR function. This review focuses on structural and dynamic studies of the extracellular ligand binding domain of iGluRs and the pore region of K+ channels that have contributed to mechanistic insights into the processes of iGluR gating and desensitization

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Ionotropic Glutamate Receptor Recognition and Activation

Ionotropic glutamate receptors are the major excitatory neurotransmitters in mammalian brain but are found throughout the animal kingdom as well as in plants and bacteria. A great deal of progress in understanding the structure of these essential neurotransmitter receptors has been made since the first examples were cloned and sequenced in 1989. The atomic structure of the ligand-binding domain of several ionotropic glutamate receptors has been determined, and a great deal of progress has been made in relating the structural properties of the binding site to the function of the intact receptor. In addition, the identification of glutamate receptors from a wide variety of organisms ranging from several types of bacteria to Arabidopsis to a range of animal species has made glutamate receptors a molecular laboratory for studying the evolution of proteins. The fact that glutamate receptors are a particularly ancient intercellular signaling molecule suggests a potential role in the transition from single celled to multicellular organisms. This review focuses on the structure and dynamics of ionotropic glutamate receptors and their relation to the function and evolution of these proteins.

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.

A molecular dynamics investigation of mono and dimeric states of the outer membrane enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca(2+)) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca(2+)-bound dimer, and of the Ca(2+)-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.

The structural basis for kainoid selectivity at AMPA receptors revealed by low-mode docking calculations

The kainoids are a class of excitatory and excitotoxic pyrrolidine dicarboxylates that act at ionotropic glutamate receptors. The kainoids bind kainate receptors with high affinity and, while binding affinity is lower at AMPA (alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptors, they are active in functional assays at this receptor subtype as well. However, kainoids are only partial agonists at AMPA receptors. Currents evoked by kainoids have been described as either slowly desensitizing, partially desensitizing, or non-desensitizing. Recently acquired X-ray crystal structures of the ligand binding domain of the iGluR2, AMPA sensitive receptor suggest that differences in ligand-receptor interactions may influence functional properties of an agonist. In an effort to identify important ligand-receptor interactions of various kainoids, we have conducted a series of low-mode docking searches of AMPA agonists in the iGluR2 binding domain. Kainic acid exhibited alternate low-lying geometries, with loss of hydrogen bonds to domain 2, which may represent a dissociation route not available to other kainoids. The most potent of the kainoids are capable of forming hydrogen bonding interactions that span the two domains of the receptor. In particular, a hydrogen bond between the domoic acid C6' carboxylic acid and Ser652 may prevent a peptide bond rotation that is associated with the desensitized state of the receptor.

Stereostructure-activity studies on agonists at the AMPA and kainate subtypes of ionotropic glutamate receptors

(S)-Glutamic acid (Glu), the major excitatory neurotransmitter in the central nervous system, operates through ionotropic as well as metabotropic receptors and is considered to be involved in certain neurological disorders and degenerative brain diseases that are currently without any satisfactory therapeutic treatment. Until recently, development of selective Glu receptor agonists had mainly been based on lead compounds, which were frequently naturally occurring excitants structurally related to Glu. These Glu receptor agonists generally contain heterocyclic acidic moieties, which has stimulated the use of bioisosteric replacement approaches for the design of subtype-selective agonists. Furthermore, most of these leads are conformationally restricted and stereochemically well-defined Glu analogs. Crystallization of the agonist binding domain of the GluR2 subunit of the (RS)-2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid (AMPA) receptor subtype of ionotropic Glu receptors in the presence or absence of an agonist has provided important information about ligand-receptor interaction mechanisms. The availability of these binding domain crystal structures has formed the basis for rational design of ligands, especially for the AMPA and kainate subtypes of ionotropic Glu receptors. This mini-review will focus on structure-activity relationships on AMPA and kainate receptor agonists with special emphasis on stereochemical and three-dimensional aspects.

Molecular dynamics simulation of the ligand binding domain of mGluR1 in response to agonist and antagonist binding

The interdomain movements of the ligand binding domain (LBD) of mGluR1 in response to agonist or antagonist binding are studied by 2 ns molecular dynamics (MD) simulations. Our results indicate that MD is able to reproduce many of the experimentally determined features of the open and closed conformations of LBD. Analysis of the ligand behavior over time allows to delineate some of the molecular determinants responsible for the agonist-induced or antagonist-blocked LBD responses.

Glutamate receptors: desensitizing dimers

Recent structural studies show, not only how the desensitization of a ligand-gated ion channel with bound agonist can be rationalized in terms of subunit-subunit instability, but also how a previously unknown mode of interaction may provide clues into how the receptor is tetramerically assembled in vivo.

Structural mobility of the extracellular ligand-binding core of an ionotropic glutamate receptor. Analysis of NMR relaxation dynamics

Ionotropic glutamate receptors play important roles in a variety of neuronal processes and have been implicated in multiple neurodegenerative diseases. The extracellular ligand-binding (S1S2) core of the GluR2 subtype can be expressed in bacteria as a soluble, monomeric protein with binding properties essentially identical to those of the intact receptor. The crystal structure of this protein has been determined in the presence and absence of various agonists and antagonists [Armstrong, N., Sun, Y., Chen, G. Q., and Gouaux, E. (1998) Nature 395, 913-917; Armstrong, N., and Gouaux, E. (2000) Neuron 28, 165-181]. The protein consists of two lobes, with the S1 segment composing the majority of lobe 1 and the S2 segment composing most of lobe 2. A domain closure upon ligand binding has been postulated, but details of intradomain motions have not been investigated. In this paper, the backbone motions of the ligand-binding core of GluR2 bound to glutamate were studied using (15)N longitudinal (T1) and transverse (T2) relaxation measurements as well as [1H]-15N nuclear Overhauser effects at 500 and 600 MHz. Residues in the agonist-binding pocket exhibited two main classes of motion. Those contacting the alpha-substituents of the ligand glutamate exhibited minimal internal motion, while those contacting the gamma-constituents exhibited exchange dynamics, indicating two dynamically distinct portions of the binding pocket. Also, two residues in transdomain linkers between lobes 1 and 2 show exchange, lending new insight into the previously proposed domain closure hypothesis. Finally, concerted motion of helix F suggests a pathway for ligand dissociation without the necessity of domain reopening.