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

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

Exploring the bioactivity of pentacyclic triterpenoids as potential antimycobacterial nutraceutics: Insights through comparative biomolecular modelling

A group of bioactive compounds known as triterpenoids, which are often found in plant materials, have been tested to possess nutritional and pharmaceutical activity. These plant components are referred to as nutraceuticals, and are used as therapeutic agents. In this study, we explore the interactions of betulinic acid (BA), oleanolic acid (OA), ursolic acid (UA), and maslinic acid (MA) against FadA5. Studies have identified FadA5, a trifunctional enzyme-like thiolase, as a target towards Mycobacterium tuberculosis inhibition. The investigation involves molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) applications. Analyses of the four pentacyclic triterpenoids binding to FadA5 showed appreciable bioactivity against FadA5. The application of two or more theoretical models to unravel ligand-enzyme binding energies can pave the way for accurate binding affinity prediction and validation.

Effects of tosyl-l-arginine methyl ester (TAME) on the APC/c subunits: An in silico investigation for inhibiting cell cycle

The anaphase-promoting complex/cyclosome (APC/c) is requisite for controlling mitosis, which is activated by Cdh1 and Cdc20 activators. Dysregulation of APC/c is observed in many cancers and is known as a targeted drug particularly in cancer drug resistance. It was shown that tosyl-l-arginine methyl ester (TAME), via mimicking isoleucine-arginine (IR) tail of co-activators, inhibits APC/c functions. However, structure details and interaction of TAME with APC/c are poorly defined. In the current study, a well-established set of computational methods was used to identify the best binding pocket in order to inhibit APC activity. Therefore, the interaction of IR tail and Cbox of co-activators, as well as TAME as an inhibitor, as an inhibitor, with APC3 and APC8 subunits of APC/c were analyzed, regarding structure, molecular docking, molecular dynamics, and free binding energy. The results indicated that TAME bound to APC3 with a higher binding affinity (∼-7.3 kcal/mol) than APC8 (∼-5.7 kcal/mol). Also, the binding free energy value obtained for the APC3-TAME was -22.25 ± 1.12 kcal/mol. According to binding free energies, van der Waals energy was the major favorable contributor to the ligand binding. These results offer that TAME had more affinity to interact with the APC3 subunit, at the IR binding pocket than the APC8 subunit at the Cbox binding pocket. In conclusion, IR binding pocket can serve as an appropriate potential target for TAME as an inhibitor of APC/c.

Structure–Activity relationship studies of two dietary flavonoids and their Nitric Oxide Synthase inhibition activity by spectroscopic and quantum/classical computational techniques

The Nitric Oxide Synthase inhibitory activity of two dietary phytocompounds naringenin and quercetin that belong to the family of flavonoids was studied by spectroscopic and computational methods. The determining role played by the global reactivity parameters that were calculated using spectroscopic and computational methods were correlated with the Nitric Oxide Synthase inhibitory activity by a Quantitative Structure–Activity Relationship (QSAR) study. The inter and intramolecular charge–transfer interactions responsible for the biological activity of naringenin and quercetin were studied using Natural Bond Orbital Analysis. The reactive sites of the title compounds were studied using their molecular electrostatic map. The likelihood of naringenin and quercetin to be a small drug molecule was determined by Lipinski’s rule. The molecular dynamics and docking studies were carried out to test the energetic and structural favorability of the inhibition of Nitric Oxide Synthase by small drug molecules naringenin and quercetin. Experimental cell culture-based in vitro assays were done to test theoretical predictions. A QSAR study revealed that for Nitric Oxide Synthase inhibitory activity to be better, the compounds belonging to the family of flavonoids ought to have a more negative ionization energy and chemical softness while bandgap and electrophilicity index ought to be more positive.

Computing long time scale biomolecular dynamics using quasi-stationary distribution kinetic Monte Carlo (QSD-KMC)

It is a challenge to obtain an accurate model of the state-to-state dynamics of a complex biological system from molecular dynamics (MD) simulations. In recent years, Markov state models have gained immense popularity for computing state-to-state dynamics from a pool of short MD simulations. However, the assumption that the underlying dynamics on the reduced space is Markovian induces a systematic bias in the model, especially in biomolecular systems with complicated energy landscapes. To address this problem, we have devised a new approach we call quasistationary distribution kinetic Monte Carlo (QSD-KMC) that gives accurate long time state-to-state evolution while retaining the entire time resolution even when the dynamics is highly non-Markovian. The proposed method is a kinetic Monte Carlo approach that takes advantage of two concepts: (i) the quasistationary distribution, the distribution that results when a trajectory remains in one state for a long time (the dephasing time), such that the next escape is Markovian, and (ii) dynamical corrections theory, which properly accounts for the correlated events that occur as a trajectory passes from state to state before it settles again. In practice, this is achieved by specifying, for each escape, the intermediate states and the final state that has resulted from the escape. Implementation of QSD-KMC imposes stricter requirements on the lengths of the trajectories than in a Markov state model approach as the trajectories must be long enough to dephase. However, the QSD-KMC model produces state-to-state trajectories that are statistically indistinguishable from an MD trajectory mapped onto the discrete set of states for an arbitrary choice of state decomposition. Furthermore, the aforementioned concepts can be used to construct a Monte Carlo approach to optimize the state boundaries regardless of the initial choice of states. We demonstrate the QSD-KMC method on two one-dimensional model systems, one of which is a driven nonequilibrium system, and on two well-characterized biomolecular systems.

Unraveling the Impact of High-Order Silk Structures on Molecular Drug Binding and Release Behaviors

Silk continues to amaze: over the past decade, new research threads have emerged that include the use of silk fibroin for advanced pharmaceutics, including its suitability for drug delivery. Despite this ongoing interest, the details of silk fibroin structures and their subsequent drug interactions at the molecular level remain elusive, primarily because of the difficulties encountered in modeling the silk fibroin molecule. Here, we generated an atomistic silk model containing amorphous and crystalline regions. We then exploited advanced well-tempered metadynamics simulations to generate molecular conformations that we subsequently exposed to classical molecular dynamics simulations to monitor both drug binding and release. Overall, this study demonstrated the importance of the silk fibroin primary sequence, electrostatic interactions, hydrogen bonding, and higher-order conformation in the processes of drug binding and release.

Activation and desensitization of ionotropic glutamate receptors by selectively triggering pre-existing motions

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that are key players in synaptic transmission and plasticity. They are composed of four subunits, each containing four functional domains, the quaternary packing and collective structural dynamics of which are important determinants of their molecular mechanism of function. With the explosion of structural studies on different members of the family, including the structures of activated open channels, the mechanisms of action of these central signaling machines are now being elucidated. We review the current state of computational studies on two major members of the family, AMPA and NMDA receptors, with focus on molecular simulations and elastic network model analyses that have provided insights into the coupled movements of extracellular and transmembrane domains. We describe the newly emerging mechanisms of activation, allosteric signaling and desensitization, as mainly a selective triggering of pre-existing soft motions, as deduced from computational models and analyses that leverage structural data on intact AMPA and NMDA receptors in different states.

Insights from a combination of theoretical and experimental methods for probing the biomolecular interactions between human serum albumin and clomiphene

In this study, the interaction of clomiphene (CLO), a non-steroidal and ovulatory stimulant drug employed in the treatment of infertility, with human serum albumin (HSA), the most abundant plasma transport protein, was investigated using spectrofluorometric, FT-IR, UV-Vis, and molecular modeling methods. The obtained results indicated that the binding of CLO to HSA led to intense fluorescence quenching of HSA via a static quenching mechanism, and that the process of CLO binding to HSA was enthalpy driven. By using experimental and theoretical methods, it was confirmed that as a result of binding CLO, slight conformational changes in HSA occurred. Also, the negative ΔH of interaction indicated that the binding of CLO with HSA was mainly enthalpy driven. The experimental and computational results suggested that hydrogen bonds and van der Waals interactions played a major role in the binding, with overall binding constants of K = 3.67 × 10⁹ M⁻¹ at 286 K and 6.52 × 10⁵ mol L⁻¹ at 310 K. Moreover, the results of molecular modeling showed that Asp234, Phe228, Leu327, and Arg209 in HSA had the highest interaction energies with the ligand.

In this study, the interaction of clomiphene with human serum albumin (HSA), the most abundant plasma transport protein, was investigated using spectrofluorometric, FT-IR, UV-Vis, and molecular modeling methods.

Neurotransmitter Funneling Optimizes Glutamate Receptor Kinetics

Ionotropic glutamate receptors (iGluRs) mediate neurotransmission at the majority of excitatory synapses in the brain. Little is known, however, about how glutamate reaches the recessed binding pocket in iGluR ligand-binding domains (LBDs). Here we report the process of glutamate binding to a prototypical iGluR, GluA2, in atomistic detail using unbiased molecular simulations. Charged residues on the LBD surface form pathways that facilitate glutamate binding by effectively reducing a three-dimensional diffusion process to a spatially constrained, two-dimensional one. Free energy calculations identify residues that metastably bind glutamate and help guide it into the binding pocket. These simulations also reveal that glutamate can bind in an inverted conformation and also reorient while in its pocket. Electrophysiological recordings demonstrate that eliminating these transient binding sites slows activation and deactivation, consistent with slower glutamate binding and unbinding. These results suggest that binding pathways have evolved to optimize rapid responses of AMPA-type iGluRs at synapses.

An improved quantum biochemistry description of the glutamate–GluA2 receptor binding within an inhomogeneous dielectric function framework

A new methodology to define the inhomogeneous dielectric constant of protein residues, to apply to the calculation of protein–ligand properties such as the electrostatic interaction.

Structural Heterogeneity of CNGA1 Channels Revealed by Electrophysiology and Single-Molecule Force Spectroscopy

The determination at atomic resolution of the three-dimensional molecular structure of membrane proteins such as receptors and several ion channels has been a major breakthrough in structural biology. The molecular structure of several members of the superfamily of voltage-gated ionic channels such as K⁺ and Na⁺ is now available. However, despite several attempts, the molecular structure at atomic resolution of the full cyclic nucleotide-gated (CNG) ion channel, although a member of the same superfamily of voltage-gated ion channels, has not been obtained yet, neither by X-ray crystallography nor by electron cryomicroscopy (cryo-EM). It is possible that CNG channels have a high structural heterogeneity, making difficult crystallization and single-particle analysis. To address this issue, we have combined single-molecule force spectroscopy (SMFS) and electrophysiological experiments to characterize the structural heterogeneity of CNGA1 channels expressed in Xenopus laevis oocytes. The unfolding of the cytoplasmic domain had force peaks, occurring with a probability from 0.2 to 0.96. Force peaks during the unfolding of the transmembrane domain had a probability close to 1, but the distribution of the increase in contour length between two successive force peaks had multiple maxima differing by tens of nanometers. Concomitant electrophysiological experiments showed that the rundown in mutant channels S399C is highly variable and that the effect of thiol reagents when specific residues were mutated was consistent with a dynamic structural heterogeneity. These results show that CNGA1 channels have a wide spectrum of native conformations that are difficult to detect with X-ray crystallography and cryo-EM.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Agonist and antagonist binding in human glycine receptors

The human glycine receptor (hGlyR) is an anion-permeable ligand-gated channel that is part of a larger superfamily of receptors called the Cys-loop family. hGlyRs are particularly amenable to single-channel recordings, thus making them a model experimental system for understanding the Cys-loop receptor family in general. Understanding the relationship between agonist binding and efficacy in Cys-loop receptors should improve our future prospects for making specific agonists or antagonists. However, at present, there is no high-resolution structure for the complete hGlyR, and thus, modeling is needed to provide a physical framework on which to interpret single-channel data. The structure of the glutamate-gated chloride channel from Caenorhabditis elegans shows a much higher level of sequence identity to human hGlyR than previous templates such as AChBP or the bacterial channels, GLIC and ELIC. Thus, we constructed a model of the hGlyR and validated it against previously reported mutagenesis data. We used molecular dynamics to refine the model and to explore binding of both an agonist (glycine) and an antagonist (strychnine). The model shows excellent agreement with previous data but also suggests some unique features: (i) a water molecule that forms part of the binding site and allows us to account for some previous results that were difficult to reconcile, (ii) an interaction of the glycine agonist with S129, and (iii) an effect from E211, both of which we confirmed with new site-directed mutagenesis and patch clamp recordings. Finally, examination of the simulations suggests that strychnine binding induces movement to a conformational state distinct from the glycine-bound or apo state, not only within the ligand-binding domain but also in the transmembrane domain.

Molecular dynamics simulations of the adenosine A2a receptor: structural stability, sampling, and convergence

Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.

Computational study of synthetic agonist ligands of ionotropic glutamate receptors

Neurological glutamate receptors are among the most important and intensely studied protein ligand binding systems in humans. They are crucial for the functioning of the central nervous system and involved in a variety of pathologies. Apart from the neurotransmitter glutamate, several artificial, agonistic and antagonistic ligands are known. Of particular interest here are novel photoswitchable agonists that would open the field of optogenetics to glutamate receptors. The receptor proteins are complex, membrane-bound multidomain oligomers that undergo large scale functional conformational changes, making detailed studies of their atomic structure challenging. Therefore, a thorough understanding of the microscopic details of ligand binding and receptor activation remains elusive in many cases. This topic has been successfully addressed by theoretical studies in the past and in this paper, we present extensive molecular dynamics simulation and free energy calculation results on the binding of AMPA and an AMPA derivative, which is the basis for designing light-sensitive ligands. We provide a two-step model for ligand binding domain activation and predict binding free energies for novel compounds in good agreement to experimental observations.

Atomistic mechanism for the activation and desensitization of an AMPA-subtype glutamate receptor

Ionotropic glutamate receptors (iGluRs) mediate fast excitatory synaptic transmission in the central nervous system. Upon agonist binding, an iGluR opens to allow the flow of cations and subsequently enters into a desensitized state. It remains unclear how agonist binding to the ligand-binding domain (LBD) is transmitted to the transmembrane domain (TMD) for channel activation and desensitization. Here we report molecular dynamics simulations of an AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid)-subtype iGluR in explicit water and membrane. Channel opening and closing were observed in simulations of the activation and desensitization processes, respectively. The motions of the LBD-TMD linkers along the central axis of the receptor and in the lateral plane contributed cooperatively to channel opening and closing. The detailed mechanism of channel activation and desensitization suggested by the simulations here is consistent with existing data and may serve as a guide for new experiments and for the design of pharmacological agents.

Quantifying water-mediated protein-ligand interactions in a glutamate receptor: a DFT study

It is becoming increasingly clear that careful treatment of water molecules in ligand-protein interactions is required in many cases if the correct binding pose is to be identified in molecular docking. Water can form complex bridging networks and can play a critical role in dictating the binding mode of ligands. A particularly striking example of this can be found in the ionotropic glutamate receptors. Despite possessing similar chemical moieties, crystal structures of glutamate and α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) in complex with the ligand-binding core of the GluA2 ionotropic glutamate receptor revealed, contrary to all expectation, two distinct modes of binding. The difference appears to be related to the position of water molecules within the binding pocket. However, it is unclear exactly what governs the preference for water molecules to occupy a particular site in any one binding mode. In this work we use density functional theory (DFT) calculations to investigate the interaction energies and polarization effects of the various components of the binding pocket. Our results show (i) the energetics of a key water molecule are more favorable for the site found in the glutamate-bound mode compared to the alternative site observed in the AMPA-bound mode, (ii) polarization effects are important for glutamate but less so for AMPA, (iii) ligand-system interaction energies alone can predict the correct binding mode for glutamate, but for AMPA alternative modes of binding have similar interaction energies, and (iv) the internal energy is a significant factor for AMPA but not for glutamate. We discuss the results within the broader context of rational drug-design.

Full and partial agonism of ionotropic glutamate receptors indicated by molecular dynamics simulations

Ionotropic glutamate receptors (iGluRs) are synaptic proteins that facilitate signal transmission in the central nervous system. Extracellular iGluR cleft closure is linked to receptor activation; however, the mechanism underlying partial agonism is not entirely understood. Full agonists close the bilobed ligand-binding domain (LBD), while antagonists prevent closure; the transmembrane ion channel either opens or stays closed, respectively. Although some bulky partial agonists produce intermediate iGluR-LBD closure, the available crystal structures also imply that the cleft can be shut with certain partial agonists. Recently, we have shown that the iGluR-LBD closure stage can be recreated by inserting a ligand into the closed cleft and simulating the ligand-receptor complex with molecular dynamics. Our simulations indicate that partial agonist binding does not necessarily prevent full receptor cleft closure; instead, it destabilizes cleft closure. Interdomain hydrogen bonds were studied thoroughly, and one hydrogen bond, in particular, was consistently disrupted by bound partial agonists. Accordingly, the simulation protocol presented here can be used to categorize compounds in silico as partial or full agonists for iGluRs.

Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics

Periplasmic binding proteins (PBPs) are a large family of molecular transporters that play a key role in nutrient uptake and chemotaxis in Gram-negative bacteria. All PBPs have characteristic two-domain architecture with a central interdomain ligand-binding cleft. Upon binding to their respective ligands, PBPs undergo a large conformational change that effectively closes the binding cleft. This conformational change is traditionally viewed as a ligand induced-fit process; however, the intrinsic dynamics of the protein may also be crucial for ligand recognition. Recent NMR paramagnetic relaxation enhancement (PRE) experiments have shown that the maltose binding protein (MBP) - a prototypical member of the PBP superfamily - exists in a rapidly exchanging (ns to µs regime) mixture comprising an open state (approx 95%), and a minor partially closed state (approx 5%). Here we describe accelerated MD simulations that provide a detailed picture of the transition between the open and partially closed states, and confirm the existence of a dynamical equilibrium between these two states in apo MBP. We find that a flexible part of the protein called the balancing interface motif (residues 175-184) is displaced during the transformation. Continuum electrostatic calculations indicate that the repacking of non-polar residues near the hinge region plays an important role in driving the conformational change. Oscillations between open and partially closed states create variations in the shape and size of the binding site. The study provides a detailed description of the conformational space available to ligand-free MBP, and has implications for understanding ligand recognition and allostery in related proteins.

A comparative analysis of the role of water in the binding pockets of ionotropic glutamate receptors

The binding pockets within proteins often contain water molecules. The ligand-binding core of ionotropic glutamate receptors represents an example where the binding pocket has many crystallographically reported waters, but the precise role remains unclear. It is also unclear to what extent the dynamic properties of these waters are conserved across the different receptor subtypes. In order to shed some light on these aspects we have performed multiple molecular dynamics simulations of the ligand binding core of four glutamate bound iGluR structures (GluA2, GluK1, GluK2, and GluN2A) and one apo structure (GluA2). We find that the water positions are reproduced from the simulations, but they also reveal that all but one water molecule in the binding site can be rearranged or replaced with water molecules from the bulk that enter the binding site through transient water channels. This one exception is not reported in the apo crystal structure but within 15 ns of simulation, a water molecule enters the site from the bulk suggesting that it is a favoured position regardless of the state of the protein. Further calculations demonstrate that whilst it is not needed in order to be able to predict the correct binding pose, it does contribute a large favourable interaction energy. We also find that one conserved water has a much stronger interaction with the protein in GluA2, GluK1 and GluK2 compared to the GluN2A receptor. The position of this water molecule is such that it can influence the dynamics of the proposed switch in the GluA2 and GluK1/2 receptors.

Glutamate receptor ion channels: structure, regulation, and function

The mammalian ionotropic glutamate receptor family encodes 18 gene products that coassemble to form ligand-gated ion channels containing an agonist recognition site, a transmembrane ion permeation pathway, and gating elements that couple agonist-induced conformational changes to the opening or closing of the permeation pore. Glutamate receptors mediate fast excitatory synaptic transmission in the central nervous system and are localized on neuronal and non-neuronal cells. These receptors regulate a broad spectrum of processes in the brain, spinal cord, retina, and peripheral nervous system. Glutamate receptors are postulated to play important roles in numerous neurological diseases and have attracted intense scrutiny. The description of glutamate receptor structure, including its transmembrane elements, reveals a complex assembly of multiple semiautonomous extracellular domains linked to a pore-forming element with striking resemblance to an inverted potassium channel. In this review we discuss International Union of Basic and Clinical Pharmacology glutamate receptor nomenclature, structure, assembly, accessory subunits, interacting proteins, gene expression and translation, post-translational modifications, agonist and antagonist pharmacology, allosteric modulation, mechanisms of gating and permeation, roles in normal physiological function, as well as the potential therapeutic use of pharmacological agents acting at glutamate receptors.

Collective dynamics of periplasmic glutamine binding protein upon domain closure

The glutamine binding protein is a vital component of the associated ATP binding cassette transport systems responsible for the uptake of glutamine into the cell. We have investigated the global movements of this protein by molecular dynamics simulations and principal component analysis (PCA). We confirm that the most dominant mode corresponds to the biological function of the protein, i.e., a hinge-type motion upon ligand binding. The closure itself was directly observed from two independent trajectories whereby PCA was used to elucidate the nature of this closing reaction. Two intermediary states are identified and described in detail. The ligand binding induces the structural change of the hinge regions from a discontinuous beta-sheet to a continuous one, which also enhances softness of the hinge and modifies the direction of hinge motion to enable closing. We also investigated the convergence behavior of PCA modes, which were found to converge rather quickly when the associated magnitudes of the eigenvalues are well separated.

Exploring kainate receptor pharmacology using molecular dynamics simulations

Ionotropic glutamate receptors (iGluRs) are enticing targets for pharmaceutical research; however, the search for selective ligands is a laborious experimental process. Here we introduce a purely computational procedure as an approach to evaluate ligand-iGluR pharmacology. The ligands are docked into the closed ligand-binding domain and during the molecular dynamics (MD) simulation the bi-lobed interface either opens (partial agonist/antagonist) or stays closed (agonist) according to the properties of the ligand. The procedure is tested with closely related set of analogs of the marine toxin dysiherbaine bound to GluK1 kainate receptor. The modeling is set against the abundant binding data and electrophysiological analyses to test reproducibility and predictive value of the procedure. The MD simulations produce detailed binding modes for analogs, which in turn are used to define structure-activity relationships. The simulations suggest correctly that majority of the analogs induce full domain closure (agonists) but also distinguish exceptions generated by partial agonists and antagonists. Moreover, we report ligand-induced opening of the GluK1 ligand-binding domain in free MD simulations. The strong correlation between in silico analysis and the experimental data imply that MD simulations can be utilized as a predictive tool for iGluR pharmacology and functional classification of ligands.

Selectivity and cooperativity of modulatory ions in a neurotransmitter receptor

Ions play a modulatory role in many proteins. Kainate receptors, members of the ionotropic glutamate receptor family, require both monovalent anions and cations in the extracellular milieu for normal channel activity. Molecular dynamics simulations and extensive relative binding free energy calculations using thermodynamic integration were performed to elucidate the rank order of binding of monovalent cations, using x-ray crystal structures of the GluR5 kainate receptor dimers with bound cations from the alkali metal family. The simulations show good agreement with experiments and reveal that the underlying backbone structure of the binding site is one of the most rigid regions of the protein. A simplified model where the partial charge of coordinating oxygens was varied suggests that selectivity arises from the presence of two carboxylate groups. Furthermore, using a potential of mean force derived from umbrella sampling, we show that the presence of cations lower the energy barrier for anion approach and binding in the buried anion binding cavity.

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.

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.

Mechanism of partial agonism at the GluR2 AMPA receptor: Measurements of lobe orientation in solution

The mechanism by which the binding of a neurotransmitter to a receptor leads to channel opening is a central issue in molecular neurobiology. The structure of the agonist binding domain of ionotropic glutamate receptors has led to an improved understanding of the changes in structure that accompany agonist binding and have provided important clues about the link between these structural changes and channel activation and desensitization. However, because the binding domain has exhibited different structures under different crystallization conditions, understanding the structure in the absence of crystal packing is of considerable importance. The orientation of the two lobes of the binding domain in the presence of a full agonist, an antagonist, and several partial agonists was measured using NMR spectroscopy by employing residual dipolar couplings. For some partial agonists, the solution conformation differs from that observed in the crystal. A model of channel activation based on the results is discussed.

The NR1 M3 domain mediates allosteric coupling in the N-methyl-D-aspartate receptor

N-Methyl-D-aspartate (NMDA) receptors play a critical role in both development of the central nervous system and adult neuroplasticity. However, although the NMDA receptor presents a valuable therapeutic target, the relationship between its structure and functional properties has yet to be fully elucidated. To further explore the mechanism of receptor activation, we characterized two gain-of-function mutations within the NR1 M3 segment, a transmembrane domain proposed to couple ligand binding and channel opening. Both mutants (A7Q and A7Y) displayed significant glycine-independent currents, indicating that their M3 domains may preferentially adopt a more activated conformation. Substituted cysteine modification experiments revealed that the glycine binding clefts of both A7Q and A7Y are inaccessible to modifying reagents and resistant to competitive antagonism. These data suggest that perturbation of M3 can stabilize the ligand binding domain in a closed cleft conformation, even in the absence of agonist. Both mutants also displayed significant glutamate-independent current and insensitivity to glutamate-site antagonism, indicating partial activation by either glycine or glutamate alone. Furthermore, A7Q and A7Y increased accessibility of the NR2 M3 domain, providing evidence for intersubunit coupling at the transmembrane level and suggesting that these NR1 mutations dominate the properties of the intact heteromeric receptor. The equivalent mutations in NR2 did not exhibit comparable phenotypes, indicating that the NR1 and NR2 M3 domains may play different functional roles. In summary, our data demonstrate that the NR1 M3 segment is functionally coupled to key structural domains in both the NR1 and NR2 subunits.

Rigid body essential X-ray crystallography: distinguishing the bend and twist of glutamate receptor ligand binding domains

The ligand-binding domain (LBD) from the ionotropic glutamate receptor subtype 2 (GluR2) has been shown to adopt a range of ligand-dependent conformational states. These states have been described in terms of the rotation required to fit subdomain (lobe) 2 following superposition of subdomain (lobe) 1. The LBD has a closed-cleft conformation for full agonists, but partial agonists induce a range of closure, which in turn controls the open probability of discrete subconductance states in the full-length receptor. Although this description is useful, it may not account for all physiologically important motions that the receptor undergoes. We have used an approach that combines the methods of essential dynamics and rigid-body dynamics to analyze 124 monomer domains from 55 crystal structures of the GluR2 LBD. We are able to show that partial agonists also induce a significant amount of twist that would not be anticipated using one rotational descriptor between apo and full-agonist-bound states. Furthermore, one of the crystal structures (chain B from 1P1U, the GluR2 L650T-AMPA complex), which has been suggested to represent an agonist-bound inactive form of the receptor, lies at the extreme of this twist motion. We suggest that partial agonists not only prevent full closure but also move the receptor closer to this inactive state. We demonstrate additionally how the method can be used to compare the results of molecular dynamics simulations with the crystallographic data and the extent to which the conformational space explored by both overlaps.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) because of difficulties associated with expression and crystallization. As such, it is one area in which computational studies, particularly molecular dynamics (MD), can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, although not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Constitutive activation of the N-methyl-D-aspartate receptor via cleft-spanning disulfide bonds

Although the N-methyl-D-aspartate (NMDA) receptor plays a critical role in the central nervous system, many questions remain regarding the relationship between its structure and functional properties. In particular, the involvement of ligand-binding domain closure in determining agonist efficacy, which has been reported in other glutamate receptor subtypes, remains unresolved. To address this question, we designed dual cysteine point mutations spanning the NR1 and NR2 ligand-binding clefts, aiming to stabilize these domains in closed cleft conformations. Two mutants, E522C/I691C in NR1 (EI) and K487C/N687C in NR2 (KN) were found to exhibit significant glycine- and glutamate-independent activation, respectively, and co-expression of the two subunits produced a constitutively active channel. However, both individual mutants could be activated above constitutive levels in a concentration-dependent manner, indicating that cleft closure does not completely prevent agonist association. Interestingly, whereas the NR2 KN disulfide was found to potentiate channel gating and M3 accessibility, NR1 EI exhibited the opposite phenotype, suggesting that the EI disulfide may trap the NR1 ligand-binding domain in a lower efficacy conformation. Furthermore, both mutants affected agonist sensitivity at the opposing subunit, suggesting that closed cleft stabilization may contribute to coupling between the subunits. These results support a correlation between cleft stability and receptor activation, providing compelling evidence for the Venus flytrap mechanism of glutamate receptor domain closure.

Characterization of protein conformational states by normal-mode frequencies

Conformational change in polymers including proteins is central to many molecular processes. Defining conformational states, however, remains a difficult and increasingly common problem, with many existing methods based on arbitrary or potentially unrepresentative measures. Furthermore, the expanding length of molecular dynamics simulations and direct observation of transitions between different energy basins suggest that this issue will only become evermore important. Methods commonly used to characterize conformational states include principal component analysis, root-mean-square deviation-based clustering, and geometric measurements such as hinge angles and distances. Here we present a method where the eigenvector frequencies derived from a Gaussian network model (Bahar, I.; Atilgan, A. R.; Erman, B. Folding Des. 1997, 2, 173-181) of a trajectory of structures from a molecular dynamics simulation are used to describe the state of the protein at each time point. We apply the method to three proteins that share the same fold as the type II periplasmic binding proteins: The lysine-arginine-ornithine-binding protein, the glutamine-binding protein, and the ligand-binding domain from the NR1 N-methyl-D-aspartate receptor. We find that the method can distinguish different states in good agreement with a variety of previous analyses and additionally provides information on the dynamic properties of that system at a given time point.

Conformational Change in an MFS Protein: MD Simulations of LacY

Molecular dynamics simulations of lactose permease (LacY) in a phospholipid bilayer reveal the conformational dynamics of the protein. In inhibitor-bound simulations (i.e., those closest to the X-ray structure) the protein was stable, showing little conformational change over a 50 ns timescale. Movement of the bound inhibitor, TDG, to an alternative binding mode was observed, so that it interacted predominantly with the N-terminal domain and with residue E269 from the C-terminal domain. In multiple ligand-free simulations, a degree of domain closure occurred. This switched LacY to a state with a central cavity closed at both the intracellular and periplasmic ends. This may resemble a possible intermediate in the transport mechanism. Domain closure occurs by a combination of rigid-body movements of domains and of intradomain motions of helices, especially TM4, TM5, TM10, and TM11. A degree of intrahelix flexibility appears to be important in the conformational change.

Dynamics and function in a bacterial ABC transporter: simulation studies of the BtuCDF system and its components

ABC transporters are integral membrane proteins which couple the energy of ATP hydrolysis to the translocation of solutes across cell membranes. BtuCD is a approximately 1100-residue protein found in the inner membrane of Gram-negative bacteria which transports vitamin B12. Vitamin B12 is bound in the periplasm by BtuF, which delivers the solute to the periplasmic entrance of the transporter protein complex BtuCD. Molecular dynamics simulations of the BtuCD and BtuCDF complexes (in a lipid bilayer) and of the isolated BtuD and BtuF proteins (in water) have been used to explore the conformational dynamics of this complex transport system. Overall, seven simulations have been performed, with and without bound ATP, corresponding to a total simulation time of 0.1 micros. Binding of ATP drives closure of the nucleotide-binding domains (NBDs) in BtuD in a symmetrical fashion, but not in BtuCD. It seems that ATP constrains the flexibility of the NBDs in BtuCD such that their closure may only occur upon binding of BtuF to the complex. Upon introduction of BtuF, and concomitant with NBD association, one ATP-binding site displays a closure, while the opposite site remains relatively unchanged. This asymmetry may reflect an initial step in the "alternating hydrolysis" mechanism and is consistent with measurements of nucleotide-binding stoichiometries. Principal components analysis of the simulation of BtuCD reveals motions that are comparable to those suggested in current transport models.

In silico mutation of cysteine residues in the ligand-binding domain of an N-methyl-D-aspartate receptor

The precise nature of redox modulation of N-methyl-d-aspartate (NMDA) receptors is still unclear, although it is thought to be related to the formation and breaking of disulfide bonds. Recent structural data demonstrated the way in which disulfide bonds in the ligand-binding core of the NR1 subunit are arranged. However, the structures were not able to reconcile existing experimental data that examined the effects of mutating these cysteine residues. We have used molecular dynamics (MD) simulations of a series of in silico mutations to try and address this in terms of the current structure of the NR1 ligand-binding domain. A double mutation that removes the disulfide bridge between C744 and C798 gives rise to greater interlobe mobility which was predicted from the crystal structure information but, unexpectedly, also appears to predispose the receptor toward greater flexibility in the hinge region. Removal of the disulfide bond between C454 and C420 did not show any appreciable difference from the "wild-type" simulation, suggesting that removal of this would not change receptor properties, which is in agreement with experimental findings. Furthermore, the position of the C454 side chain could be characterized into discrete rotamers, which may reflect the observation of alternative density in the crystal structure for this residue. Simulations in which two of the disulfide bonds are removed via mutations to alanine (C420A and C436A) resulted in a tendency of the protein to adopt a partially closed conformation.

Computer simulations of ABC transporter componentsThis paper is one of a selection of papers published in this Special Issue, entitled CSBMCB — Membrane Proteins in Health and Disease

Current computer simulation techniques provide robust tools for studying the detailed structure and functional dynamics of proteins, as well as their interaction with each other and with other biomolecules. In this minireview, we provide an illustration of recent progress and future challenges in computer modeling by discussing computational studies of ATP-binding cassette (ABC) transporters. ABC transporters have multiple components that work in a well coordinated fashion to enable active transport across membranes. The mechanism by which members of this superfamily execute transport remains largely unknown. Molecular dynamics simulations initiated from high-resolution crystal structures of several ABC transporters have proven to be useful in the investigation of the nature of conformational coupling events that may drive transport. In addition, fruitful efforts have been made to predict unknown structures of medically relevant ABC transporters, such as P-glycoprotein, using homology-based computational methods. The various techniques described here are also applicable to gaining an atomically detailed understanding of the functional mechanisms of proteins in general.

Opening and Closing Motions in the Periplasmic Vitamin B12Binding Protein BtuF†

BtuF is the periplasmic binding protein (PBP) in the vitamin B(12) uptake system in Escherichia coli where it is associated with the ABC transporter BtuCD. When the ligand binds, PBPs generally display large conformational changes, commonly termed the Venus flytrap mechanism. BtuF belongs to a group of PBPs that, on the basis of crystal structures, does not appear to display such behavior. Using 480 ns multicopy molecular dynamics simulations of apo and holo forms of the protein, we investigate the dynamics of BtuF. We find BtuF to be more flexible than previously assumed, displaying clear opening and closing motions which are more pronounced in the apo form. The protein behavior is compatible with a PBP functional model that postulates a closed conformation for the ligand-bound state, whereas the empty form fluctuates between open and closed conformations. Elastic network normal-mode analysis suggests that all BtuF-like PBPs are capable of similar opening and closing motions. It also makes the typical Venus flytrap domain motions a likely common means of how PBP-ABC transporter interaction could occur.

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.

Binding site flexibility: molecular simulation of partial and full agonists within a glutamate receptor

Ionotropic glutamate receptors mediate fast synaptic transmission in the mammalian central nervous system and play an important role in many different functions, including memory and learning. They have also been implicated in a variety of neuropathologies and as such have generated widespread interest in their structure and function. Molecular Dynamics simulations (5 x 20 ns) of the ligand-binding core of the GluR2 glutamate receptor were performed. Through simulations of both wild type and the L650T mutant, we show that the degree of protein flexibility can be correlated with the extent to which the binding cleft is open. In agreement with recent experiments, the simulations of kainate with the wild-type construct show a slight increase in beta-sheet content that we are able to localize to two specific regions. During one simulation, the protein made a transition from an open-cleft conformation to a closed-cleft conformation. This closed cleft conformation closely resembles the closed-cleft crystal structure, thus indicating a potential pathway for conformational change associated with receptor activation. Analysis of the binding pocket suggests that partial agonists possess a greater degree of flexibility within the pocket that may help to explain why they are less efficient at opening the channel than full agonists. Examination of water molecules surrounding the ligands reveals that mobility in distinct subsites can be a discriminator between full and partial agonism and will be an important consideration in the design of drugs against these receptors.

Tweaking agonist efficacy at N-methyl-D-aspartate receptors by site-directed mutagenesis

The structural basis for partial agonism at N-methyl-D-aspartate (NMDA) receptors is currently unresolved. We have characterized several partial agonists at the NR1/NR2B receptor and investigated the mechanisms underlying their reduced efficacy by introducing mutations in the glutamate binding site. Key residues were selected for mutation based on ligand-protein docking studies using a homology model of NR2B-S1S2 built from the X-ray structure of NR1-S1S2 in complex with glycine. Wild-type and mutant forms of NR2B were coexpressed with NR1 in Xenopus laevis oocytes and characterized by two-electrode voltage-clamp electrophysiology. By combining mutagenesis of residues His486 or Val686 with activation by differently substituted partial agonists, we introduce varying degrees of steric clash between the ligand and the two binding domains S1 and S2. In cases where ligand-protein docking predicts increased steric clashes between agonists and the residues forming the S1-S2 interface, the agonists clearly show decreased relative efficacy. Furthermore, we demonstrate that the mutation S690A affects both potency and efficacy in an agonist-specific manner. The results indicate that essential residues in the ligand binding pocket of NR2B may adopt different conformations depending on the agonist bound. Together, these data indicate that agonist efficacy at the NR2B subunit can be controlled by the extent of steric clashes between the agonist and the ligand binding domains and by ligand-dependent arrangements of residues within the binding pocket.

A salt-bridge motif involved in ligand binding and large-scale domain motions of the maltose-binding protein

The uptake of nutrients is essential for the survival of bacterial cells. Many specialized systems have evolved, such as the maltose-dependent ABC transport system that transfers oligosaccharides through the cytoplasmic membrane. The maltose/maltodextrin-binding protein (MBP) serves as an initial high-affinity binding component in the periplasm that delivers the bound sugar into the cognate ABC transporter MalFGK(2). We have investigated the domain motions induced by the binding of the ligand maltotriose into the binding cleft using molecular dynamics simulations. We find that MBP is predominantly in the open state without ligand and in the closed state with ligand bound. Oligosaccharide binding induces a closure motion (30.0 degrees rotation), whereas ligand removal leads to domain opening (32.6 degrees rotation) around a well-defined hinge affecting key areas relevant for chemotaxis and transport. Our simulations suggest that a "hook-and-eye" motif is involved in the binding. A salt bridge between Glu-111 and Lys-15 forms that effectively locks the protein-ligand complex in a semiclosed conformation inhibiting any further opening and promoting complete closure. This previously unrecognized feature seems to secure the ligand in the binding site and keeps MBP in the closed conformation and suggests a role in the initial steps of substrate transport.

On the Binding Determinants of the Glutamate Agonist with the Glutamate Receptor Ligand Binding Domain

Ionotropic glutamate receptors (GluRs) are ligand-gated membrane channel proteins found in the central neural system that mediate a fast excitatory response of neurons. In this paper, we report theoretical analysis of the ligand-protein interactions in the binding pocket of the S1S2 (ligand binding) domain of the GluR2 receptor in the closed conformation. By utilizing several theoretical methods ranging from continuum electrostatics to all-atom molecular dynamics simulations and quantum chemical calculations, we were able to characterize in detail glutamate agonist binding to the wild-type and E705D mutant proteins. A theoretical model of the protein-ligand interactions is validated via direct comparison of theoretical and Fourier transform infrared spectroscopy (FTIR) measured frequency shifts of the ligand's carboxylate group vibrations [Jayaraman et al. (2000) Biochemistry 39, 8693-8697; Cheng et al. (2002) Biochemistry 41, 1602-1608]. A detailed picture of the interactions in the binding site is inferred by analyzing contributions to vibrational frequencies produced by protein residues forming the ligand-binding pocket. The role of mobility and hydrogen-bonding network of water in the ligand-binding pocket and the contribution of protein residues exposed in the binding pocket to the binding and selectivity of the ligand are discussed. It is demonstrated that the molecular surface of the protein in the ligand-free state has mainly positive electrostatic potential attractive to the negatively charged ligand, and the potential produced by the protein in the ligand-binding pocket in the closed state is complementary to the distribution of the electrostatic potential produced by the ligand itself. Such charge complementarity ensures specificity to the unique charge distribution of the ligand.