Exploring the conformational and reactive dynamics of biomolecules in solution using an extended version of the glycine reactive force field

In order to describe possible reaction mechanisms involving amino acids, and the evolution of the protonation state of amino acid side chains in solution, a reactive force field (ReaxFF-based description) for peptide and protein simulations has been developed as an expansion of the previously reported glycine parameters. This expansion consists of adding to the training set more than five hundred molecular systems, including all the amino acids and some short peptide structures, which have been investigated by means of quantum mechanical calculations. The performance of this ReaxFF protein force field on a relatively short time scale (500 ps) is validated by comparison with classical non-reactive simulations and experimental data of well characterized test cases, comprising capped amino acids, peptides, and small proteins, and reaction mechanisms connected to the pharmaceutical sector. A good agreement of ReaxFF predicted conformations and kinetics with reference data is obtained.