How Sensitive Are Nanosecond Molecular Dynamics Simulations of Proteins to Changes in the Force Field?

Citations Over TimeTop 22% of 2007 papers

Abstract

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

Related Papers

• → Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D(2021)44 cited
• → Molecular Dynamics Simulations of Proteins: Can the Explicit Water Model Be Varied?(2007)71 cited
• → The quest for the best nonpolarizable water model from the adaptive force matching method(2010)59 cited
• → Self-similarity and protein compactness(2009)18 cited
• → Why the OPLS-AA Force Field Cannot Produce the β-Hairpin Structure of H1 Peptide in Solution When Comparing with the GROMOS 43A1 Force Field?(2011)11 cited