AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine

Abstract
We report a consistent set of AMBER force-field parameters for the most common phosphorylated amino acids, phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine in different protonation states. The calculation of atomic charges followed the original restrained electrostatic potential fitting procedure used to determine the charges for the parm94/99 parameter set, taking α-helical and β-strand conformations of the corresponding ACE-/NME-capped model peptide backbone into account. Missing force-field parameters were taken directly from the general AMBER force field (gaff) and the parm99 data set with minor modifications, or were newly generated based on ab initio calculations for model systems. Final parameters were validated by geometry optimizations and molecular-dynamics simulations. Template libraries for the phosphorylated amino acids in Leap format and corresponding frcmod parameter files are made available.