Conformers of gaseous protonated glycine

Ab initio geometry optimizations were performed on gaseous protonated glycine using the second-order Møller–Plesset perturbation theory with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** basis sets. Eight energy minima and 12 saddle points in the low-energy region of the electronic potential energy surface were characterized. The global minimum was an amino N-protonated conformer containing an ionic H bond between the (SINGLE BOND)NH₃⁺ and O(DOUBLE BOND)C(DIAGONAL BOND)(DIAGONAL BOND) groups. The lowest energy O-protonated conformer was stabilized by a conjugative attraction between the nitrogen lone-pair electrons and the positively charged planar fragment (SINGLE BOND)C(OH)₂⁺. Relative electronic energies of the nine N- and 11 O-protonated species fall in the ranges of 0–10 and 30–40 kcal mol⁻¹. At room temperature the equilibrium distribution contained the most stable N-protonated conformer almost exclusively. Additional subjects for investigation include the effects of basis set and electron correlation on the predicted structures, nonbonded interactions that influence the relative stability of protonated conformers, conformational interconversions based on intrinsic reaction coordinate calculations, and kinetic pathways for protonation and associated changes in Gibbs free energy. The work provides geometric, energetic, and thermodynamic data pertinent to the study of gas-phase ion chemistry of amino acids and peptides. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1862–1876, 1998