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

Citations Over TimeTop 10% of 2010 papers

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.

Related Papers

• → Parameterization of a reactive force field using a Monte Carlo algorithm(2013)100 cited
• → Global optimization of parameters in the reactive force field ReaxFF for SiOH(2013)103 cited
• → Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach(2019)37 cited
• → A test on reactive force fields for the study of silica dimerization reactions(2015)23 cited
• ReaxFF Force Field in Molecular Dynamic Simulation of Reactive Systems(2020)