Next Article in Journal
A Novel Approach of Heat Rate Enhancement in Rectangular Channels with Thin Porous Layer at the Channel Walls
Next Article in Special Issue
Malonaldehyde-like Systems: BeF2 Clusters—A Subtle Balance between Hydrogen Bonds, Beryllium Bonds, and Resonance
Previous Article in Journal
Association between Transcription Factor 7-like-2 Polymorphisms and Type 2 Diabetes Mellitus in a Ghanaian Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations

by
Maria Cristina Caputo
1 and
Patricio Federico Provasi
2,*
1
Departamento de Física, Universidad de Buenos Aires, Instituto de Física de de Buenos Aires (IFIBA)—Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires C1053ABH, Argentina
2
Department of Physics, University of Northeastern, Instituto de Modelado e Innovación Tecnológica (IMIT)—Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Corrientes CP W 3404 AAS, Argentina
*
Author to whom correspondence should be addressed.
Submission received: 27 August 2021 / Revised: 5 October 2021 / Accepted: 8 October 2021 / Published: 9 November 2021
(This article belongs to the Special Issue Define What Is Not Defined: In Chemistry and Beyond)

Abstract

:
Glycine in its neutral form can exist in the gas phase while its zwitterion form is more stable in water solution, but how many waters are actually necessary to stabilize the zwitterionic structure in the gas phase? Are the intramolecular isotropic spin spin coupling constants sensitive enough to accuse the change in the environment? or the conformer observed? These and related questions have been investigated by a computational study at the level of density functional theory employing the B3LYP functional and the 6-31++G**-J basis set. We found that at least two water molecules explicitly accounted for in the super-molecule structure are necessary to stabilize both conformers of glycine within a water polarizable continuum model. At least half of the SSCCs of both conformers are very stable to changes in the environment and at least four of them differ significantly between Neutral and Zwitterion conformation.

1. Introduction

Glycine is the smallest amino-acid conforming living organisms, and therefore of great interest in theoretical studies mainly due to the great relative computational saving its investigation demand. In particular, many of its properties are already known experimentally [1,2,3,4,5,6,7,8,9]. Among the diversity of works on glycine, some are aimed to elucidate the mechanisms of tautomerization. Hence, the proton transfer of the neutral glycine, stable in the gas phase, to zwitterionic conformer, stable in aqueous solution, has an estimated free energy barrier of about 7.3 k c a l m o l [1,2,3,5,6,8,10,11,12,13,14]. Some other recent works aim to study small clusters of Glycine: ( H 2 O ) n complexes in order to elucidate structures and binding energies [15,16,17,18], or the stabilization of Z-Glycine in water and its NMR signature [19] and the election of a method for fast accurate of 1 J ( C , H ) for application in proteins [20].
Császár [21] found the three most stable conformers of neutral glycine and named them Ip, IIp and IIIp), then Godfrey et al. [22] found a slightly different set and named them 1, 3 and 2, later Sauer et al. [23] and Aikens and Gordon [4] made the same finding as Császár’s naming the conformers, respectively, A, C and B and N1, N3 and N2). More recently, Caputo et al. [24] found the structures proposed by Császár.
According to Császár [21] and Ding and Krogh-Jespersen [25], the zwitterion spontaneously tends to structure IIp, 2, C or N3 when it is let in the gas phase, which is obviously the shorter path for the proton. Moreover, the most stable geometry we found when including a continuum polarizable medium (PCM) is precisely this last one, named in the present work by N for neutral. This conformer corresponds to the second most stable structure found by Császár, Sauer et al., Aikens and Gordon and Caputo et al.. Correspondingly, we name the other two structures analyzed in this work with Z (for zwitterion) and TS (for transition state), which were found within the PCM model.
Alonso et al. [16] reported that the conformer of glycine III, which coincides with the neutral conformer chosen for this work, in the complex glycine:water IIIa was not observed in their experiment based in a supersonic jet and they explain the reason as being due to its collisional relaxation to the most stable conformer of glycine (I). However, it is not the purpose of the present work to analyze the collisions of glycine:water.
In this work, we use structure IIp, 2, C or N3 for the neutral glycine and its zwitterion to analyze their stability within the PCM model and including explicit water molecules, up to reaching saturation of the hydrophilic positions.
The number of water molecules necessary to reach saturation vary a bit according to the conditions established for the hydrogen bond (HB) formation, see for instance Ref. [6]. Hence, we only restricted rings of water molecules that do not involve the glycine in them, so the saturation is reached with eight molecules of water. The energies scheme arising from it allows one to analyze the contributions in the energy of the HBs formed in the super-molecule model. The scheme contemplated for adding water molecules is shown in the computational details section. Finally we analyze the NMR indirect spin–spin coupling constants (J) between the intramolecular heavy atoms.

2. Materials and Methods

The calculations were performed with the Gaussian 09 suit of program [26]. The geometry optimization were done at B3LYP [27,28] with 6-31+G(d,p) [29,30,31,32] basis sets. The local minima have passed the frequency test.
The calculations of SSCCs were performed using the density functional theory with the B3LYP functional, and employed a reliable basis set 6-311++G * * -J [33]. Basis functions, for both glycine and water molecules, were taken from Basis Set Exchange [34,35,36].
Calculations include a liquid solvent described by Polarizable Continuum Model (PCM) approximation [37] with the variant (SMD) [38] that differs from the default (IEF−PCM) also by Truhlar and co-workers solvation model [39,40,41,42,43] that do not include the electrostatic terms. Vibrational corrections were not included due to these normally being too expensive to calculate in standard applications to larger molecules.
The geometries obtained in this work are essentially the same as the ones obtained by Caputo et al. in Ref. [24] and are represented in Figure 1, where the transition state (TS) was found with the help of N and Z structures within the PCM approximation mentioned above and QST3 option.
Analyses of the neutral glycine molecule allows one to identify four hydrophilic sites. They are, respectively, the two oxygens (syn- and anti-periplanar to the C-N bond) named in Table 1 as Ac (by Acid) and O T (by trans); the amino group named Am (by amino) and, finally, the possibility of simultaneous HB formation between Ac and Am, which is named AcAm.
Therefore, subsequent water molecules were added following the order and position indicated in Figure 2 and Table 1. Up to eight water molecules effectively joined to glycine were feasible to be included for the three conformers, that is, without having rings of water molecules which do not include the glycine in them. Therefore, the reported SSCCs are single point calculations in all cases.
The theory of indirect nuclear SSCCs, Reference [44] and different computational methods for calculating them has been described extensively in the literature [45,46,47,48,49,50,51]. However, it should be useful, as a reminder, to mention that there are four contributions to the SSCC: the Fermi contact (FC) and the spin-dipolar (SD), which come from the interaction of the nuclear magnetic moments with the spin of the electrons, as well as the diamagnetic spin orbital (DSO) and the paramagnetic spin orbital (PSO), which are due to the interaction of the nuclear spins with the orbital angular momentum of the electrons.

3. Results and Discussion

3.1. Basics Concepts

We have analyzed the hydration of one conformer of glycine, stable in gas phase, its zwitterion counterpart, which is stable in liquid phase, and the transition state obtained from the former two, depicted in Figure 1. The hydration of these conformers was performed in a systematic manner adding one explicit molecule of water each time until saturation was reached and also embedding these complexes in a continuum dielectric constant of water, through the PCM model, which allowed us to have a stable zwitterion conformer. The sequence for incorporating the water molecules explicitly was referred to in Figure 2 and Table 1.
The zwitterion structure needs a minimum of two water molecules in the super-molecule structure to preserve its stability in the gas phase, but within the PCM approximation the three structures are stable by themselves, that is, without any explicit water molecule. As a reminder, the SMD variant is an IEF-PCM that includes the non-electrostatic terms for Truhlar and coworkers’ solvation model. Within the PCM approximation, four water molecules are needed for Z-Gly to become more stable than the whole series of hydrated N-Gly and TS-Gly molecules, see Figure 3.
The energies [in a.u.] for the different super-molecular systems of Gly:( H 2 O ) n (with n = 1 to 8) were obtained as: E G l y + ( 8 n ) · E H 2 O . The E G l y and E H 2 O are the energies of glycine and water, respectively. All of them were obtained within the PCM model and represented in Figure 3 together with the energies of each conformer, N, TS and Z, namely with subindex 0 (zero), i.e., 0 (zero) water molecule. Clearly it can be seen in Figure 3 that two clusters of dots in the structures N-Gly and TS-Gly for n = 1 & 2 (C 1 in the figure) and n = 3 & 4 (C 2 in the figure), which means the presence of a gap between n = 0 and 1 that indicates a decrease in the formation energy due to the formation of an H-bond in the zone labeled AmAc, (1) from Table 1. The other gap, between n = 2 and 3, shows the formation of a stronger H-Bond in the zone labeled O T , (3) from Table 1.
Finally, the gap between n = 7 and 8, which is common for N and Z-Gly structures, again shows a new formation of a stronger H-Bond in the zone labeled O T and the additional H-Bonds formed between water molecules, (8) from Table 1.
Worthy of being mentioned is the fact that the addition of a water molecule in the Ac zone, (1 and 6) of Table 1; and in the Am zone, (2 and 4) of Table 1, have respectively a major and minor effect on the energy decrease in the super-molecular system. While the addition of subsequent water molecules in the O T zone, (3, 5 and 7) of Table 1 has a decreasing influence on the decrease in energy; and the addition of a water molecule in the AmAc zone, (8) of Table 1 produce a large decrease in the energy effect of the supermolecule system for N and Z-Gly.
Structure Z-Gly has a slightly different pattern with at least four gaps (G1 - G4 in the Figure) in its pattern. The first two (G1 and G2) coincide with those just seen for N-Gly and TS-Gly structures, that is, the gap between n = 2 and 3 (G1) and the gap between n = 4 and 5 (G2). Meanwhile, there are two extra gaps, one is between n = 5 and 6 (G3), that appear with the formation of the H-bond between water molecules in the Ac zone, (6) from Table 1 and finally is the gap between n = 7 and 8 (G4) that has to do with the formation of an H-bond in the AmAc zone, (8) from Table 1 and coincides with the corresponding N-Gly. Thus, taking into account the size of the gaps (Table S1 of the Supplementary Material), one sees that the first water molecules added in the O T zone (3), the second in the Ac zone (6) and the AmAc zone (8), contribute more to the downing of the energy of the super-molecular system of the Z-Gly. All the others produce a lesser downing in the super-molecule energy.
Figure 3 also shows that the super-molecular structures of the Z-Gly:W n always have lower energies than the other two. From n = 4, the Z-Gly structures are under the smallest energies for structures N-Gly and TS-Gly. Our calculations showed that adding up to eight explicit water molecules, the complex has at least one H-bond between each water and glycine and that, beyond this number and the addition of more explicit water molecules, new H-bonds are solely formed between water molecules. Hence, the super-molecule structure was limited to eight water molecules [24].
Figure 4 represents the energy difference [in a.u.] between neutral and zwitterionic conformers and their hydrated structures. In all cases the difference of energy, Δ E N Z = E N E Z , is positive, which agrees well with the pattern observed in Figure 3, where the energies of zwitterion are lower than their neutral counterparts. Only for the super-molecule n = 7, Gly:7W, the energy difference, 0.0282 a.u., is slightly smaller than that for n = 6, 0.0277 a.u. This fact, together with the almost imperceptible increase for n = 8, 0.0283 a.u. is directly related to the saturation of the super-molecule systems and the difference between N and Z tends to be constant.

3.2. Spin–Spin Coupling Constants

The intramolecular SSCCs between the heavy atoms, other than hydrogen, from one to three bonds in glycine have been investigated for the neutral and zwitterionic forms in a electrostatic embedding within PCM and adding to the calculations up to eight molecules of water in explicit form.
The isotropic values and their contributions are collected in Figure 5a–f and Tables S2–S4 of the Supplementary Material and summarized in Table 2 and Table 3.
The N- and Z-Gly with just the PCM embedding, 0W, and up to eight explicit water molecules, 8W, are shown in Table 2 (Columns 1 and 2) and compared with the best calculations of Valverde et al. [19] (Columns 3 and 4) and same experiments Ref. [52,53,54]. At first glance can be noted two main characteristics, one is the better agreement with experiment for 1 J ( N , C N ) in the present work and that the values for 1 J ( O A c , C O s ) and 1 J ( O T , C O s ) seems to be exchanged between this work and the work by Valverde et al. [19]. Discarding the possibility of human error, according to private communication, it could arise from differences in the geometries or basis sets, but a further investigation would be needed.
The complexity of the variations along the series are clearly exposed in the Figure 5, where the qualitative general behavior of the intramolecular coupling in glycine can be seen. A summarized pattern can be extracted from Table 3. Where the differences between Gly:4W - Gly:0W, Gly:8W - Gly:4W and the total difference Gly:8W - Gly:0W are exposed. The last one is just the addition of the former two. The FC contributions, which determine, in general, the behavior of the total, are also shown.
The one-bond SSCCs change as much as ∼ 3.0 Hz for Z-Gly, along the series Gly:nW, when incorporating explicitly water molecules, from n = 0 to 8, in the calculations. In particular, the 1 J ( C O s , C N ) and 1 J ( N , C N ) of the Z-Gly, 2.83 Hz and −1.31 Hz are the most sensitive to the presence of explicit water than the corresponding counterparts of N-Gly, −0.77 Hz and 0.42 Hz, see Figure 5a,d. While for Z-Gly 1 J ( C O s , C N ) increases with the incorporation of water molecules 1 J ( N , C N ) decreases, and the opposite happens with N-Gly.
The changes in 1 J ( C O s , C N ) for N-Gly arise from the first round of added water molecules and are reduced a bit by the second round. While for Z-Gly the main variation arises from the second round of water molecules and is enlarged a bit by the first round. The changes in 1 J ( N , C N ) for N- and Z-Gly seems to mirror what occurs in 1 J ( C O s , C N ) .
The other two one-bond SSCCs, 1 J ( O A c , C O s ) and 1 J ( O T , C O s ) , vary much less with the incorporation of explicit water molecules, and for Z-Gly they decrease, −0.59 Hz and −1.26 Hz respectively. While for N-Gly, the decrease is of −0.8570 Hz and −1.4373 Hz.
For N-Gly, the changes observed in 1 J ( O A c , C O s ) arise from the first round of added water molecules and is counteracted by the second round in about 46%; while for Z-Gly the second round act in the same direction as the main variation was introduced by the first round of added water molecules. Again, in 1 J ( O T , C O s ) , the effect of the first and second round of added water molecules for both N- and Z-Gly seems to mirror the previous behavior.
Remarkable is the fact that these last two SSCCs have, along the series, an almost constant difference between N-Gly and Z-Gly of about 6.0 Hz in favor of Z-Gly for the former coupling constant and about 5.0 Hz in favor of Z-Gly for the latter coupling constant. This characteristic makes any of the one-bond oxygen–carbon SSCCs an excellent marker for distinguishing between the conformers.
The one-bond carbon–carbon and nitrogen–carbon SSCCs contrast with the previous remark, since they do not seem to show many differences between conformers, Figure 5a,d. Instead, they vary in different directions according to the conformer and the coupling constant, see Table 2.
It is worth to note that due to the difference in scale of ordinate axes between Figure 5a–f, the variations observed for two- and three-bonds SSCCs are not that much as they look like at first glance, instead they all are rather quite stable.
The general behavior of all gemial SSCCs is that they are roughly stable between n = 0 and 8, with the addition of explicit water molecules with only two exceptions, 2 J ( N , C O s ) and 2 J ( O T , C N ) , which shows changes over 1.0 Hz for N-Gly and Z-Gly, respectively, see Table 3.
Figure 5b,e also allows one to see that 2 J ( N , C O s ) has a steady increment, along the series, for N-Gly but it is almost unchanged for Z-Gly, just 0.1 Hz. In the same manner, 2 J ( O T , O A c ) decreases steadily for the the first round of water molecules and increase a bit for the second round producing a net change of −0.5 Hz for Z-Gly and stays standing for N-Gly. However, the average difference along the series between conformers is roughly constant and is about 4.0 Hz in both cases. Again, this fact would let experiments distinguish between conformers in a quite reliable manner.
For the two-bond SSCCs, the second round of added water molecules is more influential than the first round, with only two exceptions, 2 J ( O A c , C N ) and 2 J ( O T , O A c ) , in both cases for Z-Gly.
The net changes from 0W to 8W in the three-bond SSCCs are between the smallest and are of only a few tenths, Table 3. However, Figure 5c,f shows large variations in the intermediate stages of 3 J ( O A c , N ) for both compounds, but mainly for N-Gly. These large variations can be attributable to the geometrical changes in the conformation of hydrogens in the glycines compounds, N and Z, which affect more to long nuclear interactions, mainly due to their small values.
The 3 J ( O T , N ) coupling constants shows similar trends and magnitudes for both compounds. This fact makes this coupling constant very reliable for experiments but would be very difficult to distinguish between conformers.
As seen previously, at least two of the four stable SSCCs that show a large difference between N- and Z-Gly contain at least one oxygen, 1 J ( N , C N ) , 1 J ( O A c , C O s ) , 2 J ( N , C O s ) and 2 J ( O T , O A c ) , see Figure 5 and Tables S2–S4 of the Supplementary Material.
The 3 J ( O T , N ) SSCC is very stable but shows similar magnitudes for both conformers therefore is unable to discern between them.
Thus, it is worth mentioning that the oxygen, in general, has received less attention than the other elements of the first and second row of the periodic table, mainly due to the small abundance of its only magnetic nucleus, which is 0.037% and the fact that it possess a nuclear spin I = 5/2 that implies a quadrupole moments, responsible of the broadening of the NMR line widths, see for instance, Ref. [55,56,57,58,59,60,61,62] among many others.
Moreover, it is also worth to note that for the geminal coupling constants 2 J ( O T , O A c ) the PSO contribution is dominant by far and determines the sign of the coupling constant. The FC contribution is only 75% of it. The FC and DSO contributions hardly compensate for the dominance of the PSO contribution. For the 1 J ( O A c , C O s ) , 1 J ( O T , C O s ) coupling constants the PSO contributions are about 30 % of the total.

4. Conclusions

In this work, we have systematically studied the energy and the intramolecular SSCCs of the heavy nuclei, other than hydrogens, for neutral and zwitterionic glycine within a dielectric constant of water using the PCM model and including explicit water molecules increasing the number of them from zero to eight where all geometrical parameters are allowed to relax.
Within the PCM embedding, the three structures, N, Z and TS, are stable. The zwitterionic structure needs at least two explicit water molecules in the super-molecule structure to preserve its stability in the gas phase and four explicit water molecules to become more stable than the whole series of hydrated neutral molecules.
From six to eight water molecules, the energy difference between N-Gly and Z-Gly becomes approximately constant.
The analyses of the intramolecular SSCCs’ reflexes, along the series, found some very stable coupling constants for both conformers, like the four one-bond coupling constants, 1 J ( C O s , C N ) , 1 J ( N , C N ) , 1 J ( O A c , C O s ) and 1 J ( O T , C O s ) ; two of the two-bond SSCCs, 2 J ( N , C O s ) and 2 J ( O T , O A c ) ; and one of the three-bond SSCCs, 3 J ( O T , N ) . Three of these coupling constants have very similar magnitudes to each other’s conformers, that is, they do not distinguish either between conformers nor suggest changes in the environment. This is the case of SSCCs 1 J ( C O s , C N ) , 1 J ( N , C N ) and 3 J ( O T , N ) .
Therefore, the other four SSCCs are also very stable along the series; they almost do not suggest changes in the environment, but they exhibit quite different magnitudes for both conformers, over 4.0 Hz, making them excellent markers for experiments.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/sci3040041/s1, Table S1: Energy [in a. u.] for the N-, TS- and Z-Glycine super-molecular systems with 0 to 8 water molecules, Table S2: One- two- and three-bond spin-spin coupling constants (in Hz) obtained at the B3LYP/6-31G* level of calculations using PCM for both neutral (N) and zwitterionic (Z) conformers in function of the number of water molecules, Table S3: One- two- and three-bond spin-spin coupling constants (in Hz) obtained at the B3LYP/6-31G* level of calculations using PCM for both neutral (N) and zwitterionic (Z) conformers in function of the number of water molecules, Table S4: One- two- and three-bond spin-spin coupling constants (in Hz) obtained at the B3LYP/6-31G** level of calculations using PCM for both neutral (N) and zwitterionic (Z) conformers in function of the number of water molecules.

Author Contributions

Conceptualization, M.C.C. and P.F.P.; methodology, M.C.C. and P.F.P.; software, M.C.C. and P.F.P.; validation, M.C.C. and P.F.P.; formal analysis, M.C.C. and P.F.P.; investigation, M.C.C. and P.F.P.; resources, M.C.C. and P.F.P.; data curation, M.C.C. and P.F.P.; writing—original draft preparation, M.C.C. and P.F.P.; writing—review and editing, M.C.C. and P.F.P.; visualization, M.C.C. and P.F.P.; supervision, M.C.C. and P.F.P.; project administration, M.C.C. and P.F.P.; funding acquisition, M.C.C. and P.F.P. All authors have read and agreed to the published version of the manuscript.

Funding

The author acknowledge financial support to the present research from CONICET (PIP: KE3 11220200100467CO).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jensen, J.H.; Gordon, M.S. On the Number of Water Molecules Necessary To Stabilize the Glycine Zwitterion. J. Am. Chem. Soc. 1995, 117, 8159–8170. [Google Scholar] [CrossRef] [Green Version]
  2. Fernandez-Ramos, A.; Smedarchina, Z.; Siebrand, W.; Zgierski, M.Z. A direct-dynamics study of the zwitterion-to-neutral interconversion of glycine in aqueous solution. J. Chem. Phys. 2000, 113, 9714–9721. [Google Scholar] [CrossRef]
  3. Karmacharya, R.; Antoniou, D.; Schwartz, S.D. Nonequilibrium Solvation and the Quantum Kramers Problem: Proton Transfer in Aqueous Glycine. J. Phys. Chem. A 2001, 105, 2563–2567. [Google Scholar] [CrossRef]
  4. Aikens, C.M.; Gordon, M.S. Incremental Solvation of Nonionized and Zwitterionic Glycine. J. Am. Chem. Soc. 2006, 128, 12835–12850. [Google Scholar] [CrossRef] [Green Version]
  5. Leung, K.; Rempe, S.B. Ab initio molecular dynamics study of glycine intramolecular proton transfer in water. J. Chem. Phys. 2005, 122, 184506. [Google Scholar] [CrossRef] [Green Version]
  6. Campo, M.G. Molecular dynamics simulation of glycine zwitterion in aqueous solution. J. Chem. Phys. 2006, 125, 114511. [Google Scholar] [CrossRef] [PubMed]
  7. Balabin, R.M. The First Step in Glycine Solvation: The Glycine-Water Complex. J. Phys. Chem. B 2010, 114, 15075–15078. [Google Scholar] [CrossRef]
  8. Takenaka, N.; Kitamura, Y.; Koyano, Y.; Asada, T.; Nagaoka, M. Reaction path optimization and vibrational frequency analysis via ab initio QM/MM free energy gradient (FEG) method: Application to isomerization process of glycine in aqueous solution. Theor. Chem. Acc. 2011, 130, 215–226. [Google Scholar] [CrossRef]
  9. Kim, J.Y.; Ahn, D.S.; Park, S.W.; Lee, S. Gas phase hydration of amino acids and dipeptides: Effects on the relative stability of zwitterion vs. canonical conformers. RSC Adv. 2014, 4, 16352–16361. [Google Scholar] [CrossRef] [Green Version]
  10. Wada, G.; Tamura, E.; Okina, M.; Nacamura, M. On the ratio of zuitterion form to uncharged form of glycine at equilibrium in various aqueous media. Bull. Chem. Soc. Jpn. 1982, 55, 3064–3067. [Google Scholar] [CrossRef]
  11. Slifkin, M.A.; All, S.M. Thermodinamic parameters of the activation of glycine zwitterion protonation reactions. J. Mol. Liq. 1984, 28, 215–221. [Google Scholar] [CrossRef]
  12. Peteanu, L.A.; Levy, D.H. Spectroscopy of Complexes of Tryptamine and 3-Indolepropionic Acid with Various Solvents. J. Phys. Chem. 1988, 92, 6554–6561. [Google Scholar] [CrossRef]
  13. Xu, S.; Nilles, J.M.; Bowen, K.H. Zwitterion formation in hydrated amino acid, dipole bound anions: How many water molecules are required? J. Chem. Phys. 2003, 119, 10696. [Google Scholar] [CrossRef] [Green Version]
  14. Diken, E.G.; Hammer, N.I.; Johnson, M.A. Preparation and photoelectron spectrum of the glycine molecular anion: Assignment to a dipole-bound electron species with a high-dipole moment, non-zwitterionic form of the neutral core. J. Chem. Phys. 2004, 120, 9899–9902. [Google Scholar] [CrossRef]
  15. Nonose, S.; Iwaoka, S.; Mori, K.; Shibata, Y.; Fuke, K. Structures and reactions of hydrated biomolecular cluster ions. Eur. Phys. J. D 2005, 34, 315–319. [Google Scholar] [CrossRef]
  16. Alonso, J.L.; Cocinero, E.J.; Lesarri, A.; Sanz, M.E.; López, J.C. The Glycine-Water Complex. Angew. Chem. 2006, 118, 3551–3554. [Google Scholar] [CrossRef]
  17. da Silva, A.M.; Ghosh, A.; Chaudhuri, P. Effect of Hydrogen Bond Formation on the NMR Properties of Glycine-HCN Complexes. J. Phys. Chem. A 2013, 117, 10274–10285. [Google Scholar] [CrossRef]
  18. Chaudhuri, P.; Canuto, S.P.P.F. NMR spin–spin coupling constants in hydrogen-bonded glycine clusters. Int. J. Quantum Chem. 2018. [Google Scholar] [CrossRef]
  19. Valverde, D.; da Costa Ludwig, Z.M.; da Costa, C.R.; Ludwig, V.; Georg, H.C. Zwitterionization of glycine in water environment: Stabilization mechanism and NMR spectral signatures. J. Chem. Phys. 2018, 148, 024305. [Google Scholar] [CrossRef]
  20. Arroyuelo, A.; Martin, O.A.; Scheraga, H.A.; Vila, J.A. Assessing the One-Bond Cα-H spin–spin Coupling Constants in Proteins: Pros and Cons of Different Approaches. J. Phys. Chem. B 2020, 124, 735–741. [Google Scholar] [CrossRef]
  21. Császár, A.G. Conformers of Gaseous Glycine. J. Am. Chem. Soc. 1992, 114, 9568–9575. [Google Scholar] [CrossRef]
  22. Godfrey, P.D.; Brown, R.D.; Rodgers, F.M. The missing conformers of glycine and alanine: Relaxation in seeded supersonic jets. J. Mol. Struct. 1996, 376, 65–81. [Google Scholar] [CrossRef]
  23. Sauer, S.P.A.; Oddershede, J.; Sabin, J.R. Directional Dependence of the Mean Excitation Energy and Spectral Moments of the Dipole Oscillator Strength Distribution of Glycine and Its Zwitterion. J. Phys. Chem. A 2006, 110, 8811–8817. [Google Scholar] [CrossRef]
  24. Caputo, M.C.; Provasi, P.F.; Sauer, S.P.A. The role of explicit solvent molecules in the calculation of NMR chemical shifts of glycine in water. Theor. Chem. Acc. 2018, 137, 88. [Google Scholar] [CrossRef] [Green Version]
  25. Ding, Y.; Krogb-Jespersen, K. The glycine zwitterion does not exist in the gas phase: Results from a detailed ab initio electronic structure study. Chem. Phys. Lett. 1992, 199, 261–266. [Google Scholar] [CrossRef]
  26. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09 Revision D.01; Gaussian Inc.: Wallingford, UK, 2009. [Google Scholar]
  27. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  28. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Clark, T.; Chandrasekhar, J.; Spitznagel, G.W.; Schleyer, P.V.R. Efficient diffuse function-augmented basis sets for anion calculations. III. The 3-21+G basis set for first-row elements, Li-F. J. Comput. Chem. 1983, 4, 294–301. [Google Scholar] [CrossRef]
  30. Ditchfield, R.; Hehre, W.J.; Pople, J.A. Self-Consistent Molecular-Orbital Methods. IX. An Extended Gaussian- Type Basis for Molecular-Orbital Studies of Organic Molecules. J. Chem. Phys. 1971, 54, 724–728. [Google Scholar] [CrossRef]
  31. Hariharan, P.C.; Pople, J.A. The influence of polarization functions on molecular orbital hydrogenation energies. Theor. Chim. Acta 1973, 28, 213–222. [Google Scholar] [CrossRef]
  32. Hehre, W.J.; Ditchfield, R.; Pople, J.A. Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257–2261. [Google Scholar] [CrossRef]
  33. Kjær, H.; Sauer, S.P.A. Pople Style Basis Sets for the Calculation of NMR spin–spin Coupling Constants: The 6-31G-J and 6-311G-J Basis Sets. J. Chem. Theory Comput. 2011, 7, 4070–4076. [Google Scholar] [CrossRef] [PubMed]
  34. Pritchard, B.P.; Altarawy, D.; Didier, B.; Gibsom, T.D.; Windus, T.L. A New Basis Set Exchange: An Open, Up-to-date Resource for the Molecular Sciences Community. J. Chem. Inf. Model. 2019, 59, 4814–4820. [Google Scholar] [CrossRef]
  35. Feller, D. The role of databases in support of computational chemistry calculations. J. Comput. Chem. 1996, 17, 1571–1586. [Google Scholar] [CrossRef]
  36. Schuchardt, K.L.; Didier, B.T.; Elsethagen, T.; Lisong Sun, L.; Gurumoorthi, V.; Chase, J.; Li, J.; Windus, T.L. Basis Set Exchange: A Community Database for Computational Sciences. Magn. Reson. Chem. 2007, 47, 1045–1052. [Google Scholar] [CrossRef] [Green Version]
  37. Mennucci, B.; Tomasi, J. Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries. J. Chem. Phys. 1997, 106, 5151–5158. [Google Scholar] [CrossRef]
  38. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
  39. Miertus, S.; Scroco, E.; Tomasi, J. Electrostatic interaction of a solute with a continuum. A direct utilizaion of ab initio molecular potentials for the prevision of solvent effects. Chem. Phys. 1981, 55, 117–129. [Google Scholar] [CrossRef]
  40. Miertus, S.; Tomasi, J. Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes. Chem. Phys. 1982, 65, 239–245. [Google Scholar] [CrossRef]
  41. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999–3093. [Google Scholar] [CrossRef] [PubMed]
  42. Cances, E.; Mennucci, B.; Tomasi, J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J. Chem. Phys. 1997, 107, 3032. [Google Scholar] [CrossRef]
  43. Cances, E.; Mennucci, B. New applications of integral equations methods for solvation continuum models: Ionic solutions and liquid crystals. J. Math. Chem. 1998, 23, 309–326. [Google Scholar] [CrossRef]
  44. Ramsey, N.F. Electron Coupled Interactions between Nuclear Spins in Molecules. Phys. Rev. 1953, 91, 303–307. [Google Scholar] [CrossRef]
  45. Autschbach, J.; Guennic, B.L. Analyzing and interpreting NMR spin–spin coupling constants using molecular orbitalcalculations. J. Chem. Educ. 2007, 84, 156–171. [Google Scholar] [CrossRef]
  46. Sauer, S.P.A. Molecular Electromagnetism: A Computational Chemistry Approach; Oxford University Press: Oxford, UK, 2011. [Google Scholar]
  47. Helgaker, T.; Jaszuński, M.; Ruud, K. Ab Initio Methods for the Calculation of NMR Shielding and Indirect spin–spin Coupling Constants. Chem. Rev. 1999, 99, 293–352. [Google Scholar] [CrossRef]
  48. Krivdin, L.B.; Contreas, R.H. Recent advances in theoretical calculations of indirect spin–spin coupling constants. Annu. Rep. NMR Spectrosc. 2007, 61, 133–245. [Google Scholar]
  49. Vaara, J. Theory and computation of nuclear magnetic resonance parameters. Phys. Chem. Chem. Phys. 2007, 9, 5399–5418. [Google Scholar] [CrossRef] [PubMed]
  50. Helgaker, T.; Jaszunski, M.; Pecul, M. The quantum-chemical calculation of NMR indirect spin–spin coupling constants. Prog. Nucl. Magn. Reson. Spectrosc. 2008, 53, 249–268. [Google Scholar] [CrossRef]
  51. Helgaker, T.; Coriani, S.; Jørgensen, P.; Kristensen, K.; Olsen, J.; Ruud, K. Recent advances in wave function-based methods of molecular-property calculations. Chem. Rev. 2012, 112, 543–631. [Google Scholar]
  52. Sogn, J.A.; Craig, L.C.; Gibbons, W.A. Carbon-13-carbon-13 coupling constants in a series of carbon-13-enriched amino acids. J. Am. Chem. Soc. 1974, 96, 4694–4696. [Google Scholar] [CrossRef]
  53. Fermandjian, S.; Trandinh, S.; Savrda, J.; Sala, E.; Mermetbouvier, R.; Bricas, E.; Fromageot, P. 13C-nuclear magnetic resonance studies of 85% 13C-enriched amino acids and small peptides: pH effects on the chemical shifts, coupling constants, kinetics of cis-trans isomerisation and conformation aspects. Biochim. Biophys. Acta Gen. Subj. 1975, 399, 313–338. [Google Scholar] [CrossRef]
  54. Lichter, R.L.; Roberts, J.D. 15N Nuclear Magnetic Resonance Spectroscopy. XIII. Pyridine-15N1. J. Am. Chem. Soc. 1971, 93, 5218–5224. [Google Scholar] [CrossRef]
  55. Lemaitre, V.; Smith, M.E.; Watts, A. A review of oxygen-17 solid-state NMR of organic materials-towards biological applications. Solid State Nucl. Magn. Reson. 2004, 26, 215–235. [Google Scholar] [CrossRef]
  56. Sergeyev, N.M.; Sergeyeva, N.D.; Strelenko, Y.A.; Raynes, W.T. The 1H-2H, 17O-1H coupling constants and the 16O/18O induced proton isotope shift in water. Chem. Phys. Lett. 1997, 277, 142–146. [Google Scholar] [CrossRef]
  57. Halle, B.; Karlstrom, G. Prototropic charge migration in water. Part 1.-Rate constants in light and heavy water and in salt solution from oxygen-17 spin relaxation. J. Chem. Soc. Faraday Trans. 2 1983, 79, 1031–1046. [Google Scholar] [CrossRef]
  58. Broze, M.; Luz, Z. Oxygen-17 spin–spincoupling with manganese-55 and carbon-13. J. Phys. Chem. 1969, 73, 1600–1602. [Google Scholar] [CrossRef]
  59. Klemperer, W.G. 17O-NMR Spectroscopy as a Structural Probe. Angew. Chem. Int. Ed. Engl. 1978, 17, 246–254. [Google Scholar] [CrossRef]
  60. Amour, T.S.; Fiat, D. Oxygen-17 Magnetic Resonance. Bull. Magn. Reson. 1979, 1, 118–129. [Google Scholar]
  61. Gerothanassis, I.P. Oxygen-17 NMR spectroscopy: Basic principles and applications (Part I). Prog. Nucl. Magn. Reson. Spectrosc. 2010, 56, 95–197. [Google Scholar] [CrossRef]
  62. Gerothanassis, I.P. Oxygen-17 NMR spectroscopy: Basic principles and applications (part II). Prog. Nucl. Magn. Reson. Spectrosc. 2010, 57, 1–110. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Optimized geometries of glycine: N (left) Neutral, TS (centre) Transition State and Z (right) Zwitterionic form: using B3LYP/6-31+G(d,p) within PCM with water dielectric constant. Red represents oxygen, blue represents nitrogen, black represents carbon and gray represents hydrogen.
Figure 1. Optimized geometries of glycine: N (left) Neutral, TS (centre) Transition State and Z (right) Zwitterionic form: using B3LYP/6-31+G(d,p) within PCM with water dielectric constant. Red represents oxygen, blue represents nitrogen, black represents carbon and gray represents hydrogen.
Sci 03 00041 g001
Figure 2. The possible H-bonds between H 2 O and N-Gly (left)/Z-Gly (right).
Figure 2. The possible H-bonds between H 2 O and N-Gly (left)/Z-Gly (right).
Sci 03 00041 g002
Figure 3. Energies [in a.u.] of the three super-molecular systems Gly:( H 2 O ) n (with n = 1 to 8) accounting for the equation: [ E G l y + ( 8 n ) · E H 2 O ].
Figure 3. Energies [in a.u.] of the three super-molecular systems Gly:( H 2 O ) n (with n = 1 to 8) accounting for the equation: [ E G l y + ( 8 n ) · E H 2 O ].
Sci 03 00041 g003
Figure 4. Difference between the energies [in a.u.] of the Glycine, Neutral and Zwitterion conformers, in function of the number of water molecules.
Figure 4. Difference between the energies [in a.u.] of the Glycine, Neutral and Zwitterion conformers, in function of the number of water molecules.
Sci 03 00041 g004
Figure 5. One- (a,b) two- (c,d) and three-bond (e,f) spin–spin coupling constants [in Hz] obtained at the B3LYP/6-311++G * * -J within PCM level of calculations in function of the number of water molecules explicitly incorporated into the calculations.
Figure 5. One- (a,b) two- (c,d) and three-bond (e,f) spin–spin coupling constants [in Hz] obtained at the B3LYP/6-311++G * * -J within PCM level of calculations in function of the number of water molecules explicitly incorporated into the calculations.
Sci 03 00041 g005
Table 1. The possible H-bond sites between H 2 O and N-Gly/Z-Gly.
Table 1. The possible H-bond sites between H 2 O and N-Gly/Z-Gly.
H-Bond SitesN/Z-StructureInsertion Order
AcH-O-H ⋯ (OH)-C1, 6
AmH 2 O ⋯ H 2 N-C2, 4
O T H-O-H ⋯ O=C3, 5, 7
AmAcC-O-H ⋯ (OH)-H ⋯ (NH 2 )-C8
Table 2. One- two- and three-bond of the SSCC’s for N-Gly:8W and Z-Gly:8W (Last column of Tables S2–S4) of the Supplementary Material.
Table 2. One- two- and three-bond of the SSCC’s for N-Gly:8W and Z-Gly:8W (Last column of Tables S2–S4) of the Supplementary Material.
CouplingMolecule0 W8 WPCM a ASEC-FEG a Exp.
1 J ( C O s , C N ) N56.3855.61
Z55.0457.8742.356.553.6 b , c
1 J ( N , C N ) N−5.24−4.82
Z−4.09−5.400.53.36.2 d
1 J ( O A c , C O s ) N29.3828.52
Z35.1134.5136.630.7
1 J ( O T , C O s ) N28.5327.10
Z33.7032.4337.631.4
2 J ( N , C O s ) N−4.94−3.71
Z0.06−0.05
2 J ( O A c , C n ) N−1.86−1.69
Z−6.94−6.70
2 J ( O T , C N ) N−3.32−2.71
Z−10.83−9.22
2 J ( O T , O A c ) N−0.51−0.17
Z−4.50−4.99
3 J ( O A c , N ) N2.231.49
Z0.610.05
3 J ( O T , N ) N0.370.25
Z0.220.20
The values were respectivelly taken from: a Ref. [19], b Ref. [52], c taken from Ref. [53] and d taken from Ref. [54].
Table 3. Difference of the intramolecular one- two- and three-bond spin–spin coupling constants (in Hz) shown in Tables S2–S4 of the Supplementary Material, for the three situation: glycine alone (0 W), glycine + 4 water molecules (4 W) and glycine + 8 water molecules (8 W).
Table 3. Difference of the intramolecular one- two- and three-bond spin–spin coupling constants (in Hz) shown in Tables S2–S4 of the Supplementary Material, for the three situation: glycine alone (0 W), glycine + 4 water molecules (4 W) and glycine + 8 water molecules (8 W).
CouplingMolecule4W–0W8W–4W8W–0W
FCTotalFCTotalFCTotal
1 J ( C O s , C N ) N−0.91−0.940.380.17−0.53−0.77
Z0.600.472.532.363.132.83
1 J ( N , C N ) N−0.45−0.350.770.770.310.42
Z−1.17−1.23−0.15−0.09−1.32−1.31
1 J ( O A c , C O s ) N−1.71−1.601.050.74−0.66−0.86
Z−0.28−0.41−0.14−0.18−0.42−0.59
1 J ( O T , C O s ) N0.350.43−1.81−1.87−1.46−1.44
Z−0.32−0.23−1.35−1.02−1.67−1.26
2 J ( N , C O s ) N0.860.870.380.371.241.24
Z0.110.10−0.22−0.21−0.11−0.11
2 J ( O A c , C N ) N0.030.060.120.110.150.17
Z0.090.180.030.070.120.25
2 J ( O T , C N ) N−0.29−0.310.890.920.610.61
Z0.230.341.251.271.481.61
2 J ( O T , O A c ) N−1.36−1.251.441.580.080.33
Z−0.83−0.780.440.30−0.40−0.48
3 J ( O A c , N ) N1.891.82−2.69−2.56−0.80−0.74
Z0.590.56−1.21−1.12−0.63−0.56
3 J ( O T , N ) N−0.06−0.04−0.05−0.08−0.10−0.12
Z−0.030.000.03−0.020.00−0.02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Caputo, M.C.; Provasi, P.F. Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations. Sci 2021, 3, 41. https://doi.org/10.3390/sci3040041

AMA Style

Caputo MC, Provasi PF. Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations. Sci. 2021; 3(4):41. https://doi.org/10.3390/sci3040041

Chicago/Turabian Style

Caputo, Maria Cristina, and Patricio Federico Provasi. 2021. "Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations" Sci 3, no. 4: 41. https://doi.org/10.3390/sci3040041

APA Style

Caputo, M. C., & Provasi, P. F. (2021). Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations. Sci, 3(4), 41. https://doi.org/10.3390/sci3040041

Article Metrics

Back to TopTop