Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published January 20, 2011 | Supplemental Material + Accepted Version
Journal Article Open

Development of a ReaxFF Reactive Force Field for Glycine and Application to Solvent Effect and Tautomerization

Abstract

Tautomerization of amino acids between the neutral form (NF) and the zwitterionic form (ZW) in water has been extensively studied, often using glycine as a model to understand this fundamental process. In spite of many advanced studies, the tautomerization reaction remains poorly understood because of the intrinsic complexities of the system, including multiple accessible reaction pathways, charge transfer, and variations of solvation structure. To establish an accurate model that can be used for molecular dynamics simulations, a ReaxFF reactive force field has been developed for glycine. A training set for the ReaxFF hydrocarbon potential was augmented with several glycine conformers and glycine−water complexes. The force field parameters were optimized to reproduce the quantum mechanically derived energies of the species in the training set. The optimized potential could accurately describe the properties of gas-phase glycine. It was applied to investigate the effect of solvation on the conformational distribution of glycine. Molecular dynamics simulations indicated significant differences in the dominant conformers in the gas phase and in water. This suggests that the tautomerization of glycine occurs through a conformational isomerization followed by the proton transfer event. The direct reaction mechanism of the NF → ZW proton transfer reaction in water, as well as mechanisms mediated by one or two water molecules, were investigated using molecular dynamics simulations. The results suggest that the proton transfer reaction is most likely mediated by a single water molecule. The ReaxFF potential developed in this work provides an accurate description of proton transfer in glycine and thus provides a useful methodology for simulating proton transfer reactions in organic molecules in the aqueous environment.

Additional Information

© 2011 American Chemical Society. Received: September 10, 2010; Revised Manuscript Received: November 18, 2010. Publication Date (Web): December 17, 2010. A.C.T.v.D. and W.A.G. acknowledge NSF-funding DMR 0427177: ITR for the development of the ReaxFF water-potential. A.C.T.v.D. acknowledges funding from KISK startup Grant No. C000032472. D.J.D. and O.R. acknowledge support from NIH under P42 ES010344-06A2.

Attached Files

Accepted Version - nihms259594.pdf

Supplemental Material - jp108642r_si_001.txt

Files

jp108642r_si_001.txt
Files (1.8 MB)
Name Size Download all
md5:aed0fd03f56983b9d9bb8ebf5dc5711d
40.1 kB Preview Download
md5:403247c9cb69df89a5f0fe61004152a1
1.7 MB Preview Download

Additional details

Created:
August 19, 2023
Modified:
October 23, 2023