Abstract:The ff15ipq protein force field is a fixed charge model built by automated tools based on the two charge sets of the implicitly polarized charge method: one set (appropriate for vacuum) for deriving bonded parameters and the other (appropriate for aqueous solution) for running simulations. The duality is intended to treat water-induced electronic polarization with an understanding that fitting data for bonded parameters will come from quantum mechanical calculations in the gas phase. In this study, we compare ff15ipq to two alternatives produced with the same fitting software and a further expanded data set but following more conventional methods for tailoring bonded parameters (harmonic angle terms and torsion potentials) to the charge model. First, ff15ipq-Qsolv derives bonded parameters in the context of the ff15ipq solution phase charge set. Second, ff15ipq-Vac takes ff15ipq's bonded parameters and runs simulations with the vacuum phase charge set used to derive those parameters. The IPolQ charge model and associated protocol for deriving bonded parameters are shown to be an incremental improvement over protocols that do not account for the material phases of each source of their fitting data. Both force fields incorporating the polarized charge set depict stable globular proteins and have varying degrees of success modeling the metastability of short (5-19 residues) peptides. In this particular case, ff15ipq-Qsolv increases stability in a number of α-helices, correctly obtaining 70% helical character in the K19 system at 275 K and showing appropriately diminishing content up to 325 K, but overestimating the helical fraction of AAQAA3 by 50% or more, forming long-lived α-helices in simulations of a β-hairpin, and increasing the likelihood that the disordered p53 N-terminal peptide will also form a helix. This may indicate a systematic bias imparted by the ff15ipq-Qsolv parameter development strategy, which has the hallmarks of strategies used to develop other popular force fields, and may explain some of the need for manual corrections in this force fields' evolution. In contrast, ff15ipq-Vac incorrectly depicts globular protein unfolding in numerous systems tested, including Trp cage, villin, lysozyme, and GB3, and does not perform any better than ff15ipq or ff15ipq-Qsolv in tests on short peptides. We analyze the free energy surfaces of individual amino acid dipeptides and the electrostatic potential energy surfaces of each charge model to explain the differences.

Ab initio folding of Chignolins by replica-exchange molecular dynamics simulations

A Test of AMBER Force Fields in Predicting the Secondary Structure of Α-Helical and Β-Hairpin Peptides

Electronic Polarization Stabilizes Tertiary Structure Prediction of HP-36.

Direct Folding Simulation of Helical Proteins Using an Effective Polarizable Bond Force Field

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains.

Solvation and Desolvation Effects in Protein Folding: Native Flexibility, Kinetic Cooperativity and Enthalpic Barriers under Isostability Conditions

All-atom Replica Exchange Molecular Simulation of Protein BBL

The Folding Transition State of Protein L is Extensive with Nonnative Interactions (and Not Small and Polarized)

Protein folding simulations: from coarse-grained model to all-atom model.

An Evaluation of Force Field Accuracy for the Mini-Protein Chignolin using Markov State Models

Protein Folding Simulations: from Coarse-Grained Model to All-Atom Model

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

A Test of Implicit Solvent Models on the Folding Simulation of the GB1 Peptide.

Folding Kinetics and Volume Variation of the β-Hairpin Peptide Chignolin upon Ultrafast pH-Jumps

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

Effects of Turn Stability and Side-Chain Hydrophobicity on the Folding of Β-Structures

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

Links between the charge model and bonded parameter force constants in biomolecular force fields

The role of hydrophobic interactions in folding of $β$-sheets

Folding Simulations of an Α-Helical Hairpin Motif Αtα with Residue-Specific Force Fields.