Revised AMBER Parameters for Bioorganic Phosphates
Citations Over TimeTop 10% of 2012 papers
Abstract
We report AMBER force field parameters for biological simulations involving phosphorylation of serine, threonine or tyrosine. The initial parameters used RESP fitting for the atomic partial charges and standard values for all other parameters such as Lennard-Jones coefficients. These were refined with the aid of a thermodynamic cycle consisting of experimentally determined pKa values, solvation energies from molecular dynamics free energy simulations, and gas phase basicities from QM calculations. A polarization energy term was included to account for the charge density change between the gas-phase and solution, and solvation free energies were determined using thermodynamic integration. Parameter adjustment is required to obtain consistent thermodynamic results with better balanced electrostatic interactions between water and the phosphate oxygens. To achieve this we modified the phosphate oxygen radii. A thermodynamically consistent parameter set can be derived for monoanions and requires an increase of the van der Waals phosphate oxygen radii of approximately 0.09 Å. Larger, residue-specific radii appear to be needed for dianions. The revised parameters developed here should be of particular interest for environments where simulations of multiple protonation states may be of interest.
Related Papers
- → Calculation of the free energy of solvation for neutral analogs of amino acid side chains(2002)189 cited
- → Optimizing Protein–Protein van der Waals Interactions for the AMBER ff9x/ff12 Force Field(2013)27 cited
- → On the sensitivity of MD trajectories to changes in water-protein interaction parameters: The potato carboxypeptidase inhibitor in water as a test case for the GROMOS force field(1996)40 cited
- → A Solvation/Desolvation Indicator Based on van der Waals Interactions between Solvents and Porphyrins(2018)9 cited
- → Contribution of Solvation Energy in Protein‐Peptide Recognition Systems(2001)3 cited