DFT analysis of nitrogen isotopic reduction partition function ratios in proton exchange equilibria of glutamic acid species in water and its application to the estimation of nitrogen isotope fractionation

L-Glutamic acid (Glu) plays a pivotal role in amino acid metabolism in living organisms. Theoretical knowledge of how metabolism propagates the nitrogen isotope composition (δ15N) is essential for elucidating the mechanism of amino acid metabolism in vivo. In this paper, to estimate the nitrogen isotope fractionation involving Glu and its protonated/deprotonated species, their nitrogen isotopic reduction partition function ratios (RPFRs) were calculated at the apfd/6-311+G(2d,p) level of theory by the SCRF method. The results show that the RPFR values of the eight possible Glu species at 25 °C are: cation > zwitterion with α-COOH and γ-COO- > zwitterion with α-COO- and γ-COOH > anion with -NH3+ ≫ uncharged molecule > anion with -NH2 and γ-COO- > anion with -NH2 and α-COO- > di-anion. Several correlations between RPFR and bond distances were found. Most importantly, it was found that the shorter the distance of the (C-)N⋯H(-O) hydrogen bond (HB), the greater the RPFR value for that species. This correlation indicates that the tri-coordinate nitrogen atom of the amino group may take on the character of the tetra-coordinate nitrogen atom to some extent by forming a HB of this type. The nitrogen isotope exchange equilibrium between Glu and glycine/serine was evaluated at a typical physiological pH of 7.4. Calculations suggest that the equilibrium constant for the former is an increasing function of temperature from 0 to 100 °C, while the latter is a decreasing function of temperature over the same temperature range.

Conformational analysis of glutamic acid: a density functional approach using implicit continuum solvent model

Amino acids are constituents of proteins and enzymes which take part almost in all metabolic reactions. Glutamic acid, with an ability to form a negatively charged side chain, plays a major role in intra and intermolecular interactions of proteins, peptides, and enzymes. An exhaustive conformational analysis has been performed for all eight possible forms at B3LYP/cc-pVTZ level. All possible neutral, zwitterionic, protonated, and deprotonated forms of glutamic acid structures have been investigated in solution by using polarizable continuum model mimicking water as the solvent. Nine families based on the dihedral angles have been classified for eight glutamic acid forms. The electrostatic effects included in the solvent model usually stabilize the charged forms more. However, the stability of the zwitterionic form has been underestimated due to the lack of hydrogen bonding between the solute and solvent; therefore, it is observed that compact neutral glutamic acid structures are more stable in solution than they are in vacuum. Our calculations have shown that among all eight possible forms, some are not stable in solution and are immediately converted to other more stable forms. Comparison of isoelectronic glutamic acid forms indicated that one of the structures among possible zwitterionic and anionic forms may dominate over the other possible forms. Additional investigations using explicit solvent models are necessary to determine the stability of charged forms of glutamic acid in solution as our results clearly indicate that hydrogen bonding and its type have a major role in the structure and energy of conformers.