About this Journal Submit a Manuscript Table of Contents
Journal of Amino Acids
Volume 2013 (2013), Article ID 872058, 6 pages
http://dx.doi.org/10.1155/2013/872058
Research Article

In Silico Insights of L-Glutamate: Structural Features in Vacuum and in Complex with Its Receptor

Departamento de Nutrición y Bioquímica, Facultad de Ciencias, Pontificia Universidad Javeriana, Bogotá D.C., Colombia

Received 15 May 2013; Accepted 5 September 2013

Academic Editor: Hari S. Sharma

Copyright © 2013 Janneth Gonzalez and George E. Barreto. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Structural properties of the glutamate in vacuum and in complex with its receptor were analyzed. The analysis was focused on global properties, attempting to characterize features such as overall flexibility and common trends in the conformation set. The glutamate, as other ligands in complex with the receptor, adopts a spatial conformation that corresponds to one of the possible molecular equilibrium states in physiological conditions. The glutamate forms an extended structure for all cases, but the energy of the glutamate round out form is lower than the extended glutamate form. The results showed the glutamate as a flexible molecule, which can easily adapt to different interacting environments, and it can be considered as an approximation to address why glutamate interacts with a great number of molecules.

1. Introduction

The amino acid L-glutamate is considered the main excitatory neurotransmitter in mammals Central Nervous System (CNS) [1, 2]. Most of the CNS synapses use glutamate as a fast neurotransmitter, and at least 60% of the synapses in the human brain are glutamatergic [3]. The glutamate has an important physiological role in many aspects of the normal brain function such as cognition, memory, learning, nervous system development, cellular migration, cellular differentiation, and neural death [4, 5]. After being released from a presynaptic cell, glutamate diffuses across the synaptic cleft and binds to its specific receptors in the cell membrane of a postsynaptic cell. Once across the synaptic cleft, glutamate is recognized by glutamate receptors from a high variety of other molecules. Although the mechanism of molecular recognition has long been considered in a key-keyhole relationship, short-range forces have been considered as the primary cause of such interactions [6].

The glutamate (L-isomer) causes depolarization and excitation of neurons, but it does so by acting on a variety of receptors. The AMPA/kainate receptors respond to the glutamic acid analogues, α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) and kainic acid. Other receptors, N-methyl-D-aspartate (NMDA) receptor, also belong to the ion-channel-linked superfamily, and there is a population of metabotropic receptors, as upon activation they simulate a second messenger transduction system [7].

The behavior of different neurotransmitters has been mainly studied using physiological and pharmacological techniques in both vertebrate and invertebrate brains, and recent techniques in molecular biology have clarified the amino acid sequences of their binding sites [816]. However, it is important to note that natural and synthesized analogues with small differences in their chemical structures exert completely different influences and could even have opposite ones. This fact gives rise to an important question: which structural and chemical characteristics are necessary for the ligand to recognize its receptor? The accepted ideas about the molecular recognition pattern are that it occurs because a ligand and its receptor are geared toward each other in a relationship resembling that between a key and a keyhole, so their geometrical structures are determinant.

The glutamate, as other ligands in complex with the receptor, adopts a spatial conformation or tridimensional structure that corresponds to one of the possible molecular equilibrium states in physiological conditions. The conformations adopted by glutamate in complex with the receptor are the result of ligand-receptor interactions at an atomic level. It is important to consider that the allowed conformations for a molecule could have dramatic effects on the molecule activity or reactivity. Therefore, it is necessary to consider all the possible conformations for molecular properties studies.

To assess the structural properties of the glutamate, we analyzed a conformational ensemble. The analysis was focused on global properties, attempting to characterize features such as overall flexibility and common trends in the conformation set. It is important to have in mind that the different conformational analysis can be applied to any collection of molecular conformations. These may be generated by sampling techniques with theoretical structures but can also have an experimental origin, such as NMR models or different X-ray structures of the same molecule (or analogous molecules). In the present study, we calculated the structural conformations of glutamate (a) in complex with receptor and (b) in vacuum form, as an approach for understanding the mechanism of the recognition of a ligand by its receptor.

2. Materials and Methods

2.1. Experimental Structures

For this study, we used 22 structures of glutamate receptor monomers from different subunits solved by X-ray crystallography and deposited in PDB [17]. The 22 structures were derived from 7 crystals, which 5 of them belong to ionotropic receptors and 2 belong to metabotropic receptors (Table 1).

tab1
Table 1: Selected glutamate receptors solved by X-ray crystallography.

To estimate the uncertainty grade in the atomic positions from the selected glutamate structures, we measured the following internal variables for the heavy atoms: (a) bonds length; (b) bond angles; and (c) dihedral angles. The statistical analyses were made with program Statistica v.6.0 (v. 5.5). To make a comparison among structures, we performed a multiple structural alignment with the software VMD for molecular analysis [18].

2.2. Theoretical Structures

Under neutral pH in solvated conditions, the carboxyl groups and amino groups in glutamate are deprotonated and protonated, respectively. We thus used the ionized models and calculated them in vacuum to obtain the best structural conformation of glutamate. The initial conformation of glutamate was built and optimized with MOE (v. 2007-11) using standard values for the internal 0 variables. The standard nomenclature use by PDB was considered to point out the heavy atoms in the glutamate molecule (Figure 1).

872058.fig.001
Figure 1: Chemical structure of glutamate. The side chain exists in its negatively charged deprotonated carboxylate form at pHs greater than 4.1; therefore, it is also negatively charged at physiological pH ranging from 7.35 to 7.45.

It is important to note that the number of tridimensional structures expands geometrically with the number of possible combinations for the dihedral angles in the structure. Considering that the conformers differ primarily in their torsion angles, in this study we focused on the conformational space. The exploration of the conformational space was made by changing the dihedral angles in the 30 degrees step (Figure 2).

872058.fig.002
Figure 2: Dihedral angles in glutamate. The angles show the main degrees of freedom for the backbone ( and angles) and the side chain ( angles) of glutamate. The figure shows a ball-and-stick representation of glutamate, which has three angles. The fading conformations in the background illustrate a rotation around .

The database of glutamate conformers used in this study was based on minimal energy glutamate conformations, elimination of duplicated structures, and those with an atomic RMSD value higher than 0.1 Å. As a final result, 10 stable structures of glutamate were found. These structures were optimized again by GAMMES program (GAMMES v. Sep 7/2006 (R4)) using RHF/6-31G** approach. Geometry optimization was further performed for the most stable structure of each molecule, and the frequency was then calculated to confirm that the optimized structure was the ground state geometry, using the 6-31G** basis set.

3. Results

3.1. Experimental Structures

It is often assumed that crystal structures have to be obtained at sufficiently high resolution in order to perform macromolecular refinement, and what are considered “acceptable” have been pushed to lower diffraction resolutions. In most cases, the crystal quality is measured according to the resolution value. Although in this study we used two crystals resolved with low resolution, namely, 2CMO (2.65 Å) and 1ISR (4.00 Å), we found no difference when these structures were compared with the high-resolution structures.

The results show that the average of bond lengths values is not statistically different neither between monomers, which come from the same crystal, nor for monomers coming from different crystals. This result shows that at least for the studied structures, the bond length can be viewed at random, and they do not respond to the fact of coming from different crystals or the resolution degree. According to these results, all structures obtained may be used as models for the following analysis. Table 2 shows the descriptive statistics for bond length variable ( ).

tab2
Table 2: Summary of descriptive statistic for bond length between heavy atoms of glutamate conformers. Variation coefficient is in percentages. The values presented in the table come from monomers from different crystal and different experimental resolutions.

The variation coefficient is less than 1.00%, and this is consistent with the fact that the maximum range reported for the CG-CD bonds with a value that corresponds to 3.86% of the all bonds values average [17]. On the other hand, the similarity among bond angles values does not correspond, due to the fact that monomers come from different crystal or experimental resolutions. It can be considered for this variable that all the chosen structures belong to one population with close variations values (Table 3).

tab3
Table 3: Summary of descriptive statistic for bond angle between heavy atoms of glutamate conformers. Variation coefficient is in percentages. The values represent monomers from different crystal and different experimental resolutions.

The variation coefficients higher than 1.00% are the maximum values in which the CA-C-O angles have an average value of 6,554°. The degree of similarity observed for the internal variables discussed previously (bond lengths and bond angles) was not observed for the dihedral angles. These differences are mainly due to the different conformations that glutamate can take in complex with its respective original receptor nearly 22 cases analyzed. This result is a qualitative indicator that the glutamate is a dynamic molecule that shows structural heterogeneity with discrete conformational states that share energetic values in close ranges.

The RMSD value was selected as a quantitative indicator of the structures variability for alternative conformation within the same molecule. The visual analysis of multiple structural alignments allows the identification of three major clusters. These conformations emerged as a result of specific environments of the binding sites interaction with each of the receptors (Figure 3).

872058.fig.003
Figure 3: The highest scoring alignment of glutamate. The clusters formed by different conformations are correctly aligned and emerged as a result of specific environments of the binding sites interaction with each of the receptors.

The first and second groups are composed of 16 and 3 isolated molecules from ionotropic glutamate receptors crystals, and the third group is composed of 3 isolated molecules from metabotropic glutamate receptors crystals. The distribution of clusters based on the RMSD is shown in Table 4.

tab4
Table 4: ID of molecules and receptors from clusters based on pairwise main-chain RMSD. Structures in the first cluster (group 1) represent > two-thirds of the sample with a RMSD of 0.3 Å RMSD. The second cluster (groups 2 and 3) has a RMSD of 2.0 Å.

After conformational analysis, 22 refined structures were grouped into three clusters based on pairwise main-chain RMSD (Figure 3). Structures in the first cluster (cluster 1) represent > two-thirds of the sample with an RMSD of 0.3 Å RMSD from its center. The second group (clusters 2 and 3) has an RMSD of 2.0 Å from the first cluster center. A structural comparison between minimum energy structures of the three clusters is shown in Figure 4. Although it is hard to distinguish clusters by the way that the experimental restraints are satisfied (i.e., by RMSD deviations; see Table 3), we note that the 16 lowest energy structures as evaluated by the MOE scoring function all belong to cluster 1 (structural statistics are given in Table 3).

872058.fig.004
Figure 4: Differences between structural alignments. A backbone overlap according to the same alignment is shown; the backbone RMSD is ≤0.5. The figure shows that the conformations adopted by the glutamate are most stable when rounded out.
3.2. Theoretical Structures

The database obtained included 118 structures, which represent locally stable glutamate conformations; neither bond lengths nor bond angles are statistically different, and the values are in a close range (0.001 Å to 0.02 Å) for the bond angles of the 118 glutamate conformers. The results are shown in Table 5.

tab5
Table 5: The minimum and maximum values for bond angles of the 118  glutamate conformers. The values are in the range from 0.001 Å to 0.02 Å.

From the comparison of the dihedral angles from the 118 selected structures, we identified 10 groups of conformations with high structural similarity. These 10 groups were found based on different arrangements for multiple structural alignments (RMSD among groups ≥ 1.552).

The conformations included in each group share structural similarity ( ), and the graphic representation of the structural alignments between selected conformations is shown in Figure 4.

The minimum energy conformations adopted by the glutamate molecule show that the amino acid is most stable when rounded out, despite the extended form that it adopts when interacting with the receptor. This extended conformation is adopted because the lateral chain is bonded to water molecules through hydrogen bonding.

Since the ionized glutamate model in free state is not subjected to interactions with surrounded molecules, it is hypothesized that the rounded form adopted is due to the attraction between the amino and carboxyl side chain groups(Figure 5); the attraction and charge compensation between groups in the same molecule generate a propensity to form a hydrogen bond between these two groups (distance H-O 1.211 Å).

872058.fig.005
Figure 5: Representation of the glutamate in vacuo shows the propensity to form a hydrogen bond among the amino and carboxyl groups of the principal and side chain, respectively.

According to the study of the preferential glutamate conformations, it is observed that the maximal variability in the relative position of the heavy atoms is found in the side chain. Figure 6 shows the relationship among the variables energy and dihedrals. The glutamate structure in free form is determined for the dihedrals angles combinations for which the molecule energy has the lowest values.

872058.fig.006
Figure 6: Plot of energy versus dihedrals for the selected glutamate conformers. The glutamate structure in free form is determined for the dihedrals angles combinations for which the molecule energy has the lowest values.

Although all the conformers obtained share energetic values in a close range (−82.401 y−82,397 Hartrees), the local minimum differs in 1.5 Å (RMSD) average among structures.

4. Discussion

The analyses of the classical internal variables of the 22 conformers of glutamate allowed us to identify different groups of conformations according to the type of receptor, which can lead to the idea of the ligand conformation as an indicator of the binding site molecular environment. These results showed the glutamate as a flexible molecule, which can easily adapt to different interacting environments, and it can be considered as an approximation to address why glutamate interacts with a great number of molecules [16].

The comparison through the standard method of multiple structures alignment showed that the RMSD values were not dependent on the crystal resolution degree of the diffracted structures. This may be correlated with the standard error of the positional atomic parameters, which can increase if the resolution decreases. Considering that it is not possible to use reliable indicators for the positional standard errors in molecular structures, the use of the RMSD values to compare molecular structures must be considered [19].

The atomic fluctuations of the different conformational states and the structural variability were inherent conditions due to the binding site molecular environment. The molecule flexibility plays a key role in its function [2023]. Since the RMSD values of the groups obtained are closely related to each other, it was possible to use a representative structural model for each group in order to establish relationships type structural between them for subsequent studies.

The analyses of the RMSD value, according to the glutamate atomic position changes, showed that the variability increased with the atomic distance from the CA. Some hypotheses can be generated from our study: (1) the movement of the principal chain is more conservative in comparison with the side chain if molecule flexibility is considered; (2) the atomic positions of the side chain are variables as the RMSD values show. According to glutamate structure in complex with the receptor binding site reported by Tsuchiya et al. [24] and validated with different subunits of the receptor [816], it was found that the glutamate forms an extended structure for all cases, suggesting that there is no number of possibilities to find different geometries when its structure is resolved by X-ray crystallography. The extended form (Figure 3) is different to the round-out form (Figure 4). The energy of the glutamate round out form is lower than the extended glutamate form.

5. Conclusion

Glutamic acid is the key amino acid in living organisms. In order to characterize features such as overall flexibility and common trends in the conformation of glutamate in vacuum and in complex with its receptor, we have undertaken theoretical studies. The results showed the glutamate as a flexible molecule, which can easily adapt to different interacting environments, and it can be considered as an approximation to address why glutamate interacts with a great number of molecules; also our results might provide a possible approach for understanding the mechanism of the recognition of glutamate by its receptor.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgment

This work was supported in part by Grants PUJ 004837, 004247, and 004960.

References

  1. J. C. Watkins and R. H. Evans, “Excitatory amino acid transmitters,” Annual Review of Pharmacology and Toxicology, vol. 21, pp. 165–204, 1981. View at Scopus
  2. P. M. Headley and S. Grillner, “Excitatory amino acids and synaptic transmission: the evidence for a physiological function,” Trends in Pharmacological Sciences, vol. 11, no. 5, pp. 205–211, 1990. View at Scopus
  3. R. A. Nichols, T. S. Sihra, A. J. Czernik, A. C. Nairn, and P. Greengard, “Calcium/calmodulin-dependent protein kinase II increases glutamate and noradrenaline release from synaptosomes,” Nature, vol. 343, no. 6259, pp. 647–651, 1990. View at Publisher · View at Google Scholar · View at Scopus
  4. E. Shimizu, Y.-P. Tang, C. Rampon, and J. Z. Tsien, “NMDA receptor-dependent synaptic reinforcement as a crucial process for memory consolidation,” Science, vol. 290, no. 5494, pp. 1170–1174, 2000. View at Publisher · View at Google Scholar · View at Scopus
  5. D. E. Berman and Y. Dudai, “Memory extinction, learning anew, and learning the new: dissociations in the molecular machinery of learning in cortex,” Science, vol. 291, no. 5512, pp. 2417–2419, 2001. View at Publisher · View at Google Scholar · View at Scopus
  6. B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, and P. Walter, Molecular Biology of the Cell, Garland, New York, NY, USA, 4th edition, 2002.
  7. G. M. Cooper and R. E. Hausman, The Cell: A Molecular Approach, American Society for Microbiology; Washington and Sinauer Associates, Sunderland, Mass, USA, 3rd edition, 2003.
  8. M. Arvola and K. Keinänen, “Characterization of the ligand-binding domains of glutamate receptor (GluR)-B and GluR-D subunits expressed in Escherichia coli as periplasmic proteins,” Journal of Biological Chemistry, vol. 271, no. 26, pp. 15527–15532, 1996. View at Publisher · View at Google Scholar · View at Scopus
  9. K. Keinänen, A. Jouppila, and A. Kuusinen, “Characterization of the kainate-binding domain of the glutamate receptor GluR-6 subunit,” Biochemical Journal, vol. 330, no. 3, pp. 1461–1467, 1998. View at Scopus
  10. H. Furukawa and E. Gouaux, “Mechanisms of activation, inhibition and specificity: crystal structures of the NMDA receptor NR1 ligand-binding core,” EMBO Journal, vol. 22, no. 12, pp. 2873–2885, 2003. View at Publisher · View at Google Scholar · View at Scopus
  11. Y. Paas, M. Eisenstein, F. Medevielle, V. I. Teichberg, and A. Devillers-Thiéry, “Identification of the amino acid subsets accounting for the ligand binding specificity of a glutamate receptor,” Neuron, vol. 17, no. 5, pp. 979–990, 1996. View at Publisher · View at Google Scholar · View at Scopus
  12. H. Takahashi, E. Inagaki, C. Kuroishi, and T. H. Tahirov, “Structure of the Thermus thermophilus putative periplasmic glutamate/glutamine-binding protein,” Acta Crystallographica Section D, vol. 60, no. 10, pp. 1846–1854, 2004. View at Publisher · View at Google Scholar · View at Scopus
  13. M. L. Mayer, “Crystal structures of the GluR5 and GluR6 ligand binding cores: molecular mechanisms underlying kainate receptor selectivity,” Neuron, vol. 45, no. 4, pp. 539–552, 2005. View at Publisher · View at Google Scholar · View at Scopus
  14. P. Naur, B. Vestergaard, L. K. Skov, J. Egebjerg, M. Gajhede, and J. S. Kastrup, “Crystal structure of the kainate receptor GluR5 ligand-binding core in complex with (S)-glutamate,” FEBS Letters, vol. 579, no. 5, pp. 1154–1160, 2005. View at Publisher · View at Google Scholar · View at Scopus
  15. M. L. Mayer, A. Ghosal, N. P. Dolman, and D. E. Jane, “Crystal structures of the kainate receptor GluR5 ligand binding core dimer with novel GluR5-selective antagonists,” Journal of Neuroscience, vol. 26, no. 11, pp. 2852–2861, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. H. Furukawa, S. K. Singh, R. Mancusso, and E. Gouaux, “Subunit arrangement and function in NMDA receptors,” Nature, vol. 438, no. 7065, pp. 185–192, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. Protein Data bank, http://www.rcsb.org/pdb.
  18. M. H. Nanao, T. Green, Y. Stern-Bach, S. F. Heinemann, and S. Choe, “Structure of the kainate receptor subunit GluR6 agonist-binding domain complexed with domoic acid,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 5, pp. 1708–1713, 2005. View at Publisher · View at Google Scholar · View at Scopus
  19. G. Dodson and C. S. Verma, “Protein flexibility: its role in structure and mechanism revealed by molecular simulations,” Cellular and Molecular Life Sciences, vol. 63, no. 2, pp. 207–219, 2006. View at Publisher · View at Google Scholar · View at Scopus
  20. J. Chen, N. Sawyer, and L. Regan, “Protein-protein interactions: general trends in the relationship between binding affinity and interfacial buried surface area,” Protein Science, vol. 22, pp. 510–511, 2013.
  21. A. Q. Zhou, S. C. O'Hern, and L. Regan, “The power of hard-sphere models: explaining side-chain dihedral angle distributions of Thr and Val,” Biophysical Journal, vol. 102, pp. 2345–2352, 2012.
  22. D. A. Liberles, S. A. Teichmann, I. Bahar et al., “The interface of protein structure, protein biophysics, and molecular evolution,” Protein Science, vol. 21, pp. 769–785, 2012.
  23. A. Q. Zhou, S. C. O'Hern, and L. Regan, “Revisiting the Ramachandran plot from a new angle,” Protein Science, vol. 20, pp. 1166–1171, 2011.
  24. D. Tsuchiya, N. Kunishima, N. Kamiya, H. Jingami, and K. Morikawa, “Structural views of the ligand-binding cores of a metabotropic glutamate receptor complexed with an antagonist and both glutamate and Gd3+,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 5, pp. 2660–2665, 2002. View at Publisher · View at Google Scholar · View at Scopus