Abstract:Classical protein force fields were reported with too weak protein–water interactions relative to protein–protein interactions, leading to more compact structures and artificial protein aggregation. Here we investigated the impacts of scaled Lennard-Jones (LJ) interactions on the hydration of amino acids and the simulation of folded and intrinsically disordered proteins (IDPs). The obtained optimal scaling parameters reproduce accurately hydration free energies of neutral amino acid side chain analogues and do not affect the compactness and structural stability of folded proteins significantly. The scaling leads to less compact IDPs and varies from case to case. Strengthening the interactions between protein and water oxygen or hydrogen atoms by increasing the interacting LJ well depth (ε) appears more effective than weakening protein–protein interactions by reducing the interacting dispersion coefficients (C6). We demonstrate that weakening water–water interactions is a solution as well to obtaining more favorable protein–water interactions in an indirect way, although modern force fields like Amber ff19SB and a99SB-disp tend to use water models with strong water–water interactions. This is likely a compromise between strong protein–protein interactions and strong water–water interactions. Independent optimization of protein force fields and water models is therefore needed to make both interactions more close to reality, leading to good accuracy without bias or scaling.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.1c00339?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.1c00339</a>.Scaling parameters tested for the four strategies (Table S1), system setup for protein simulation (Table S2), LJ parameters for <a class="internalNav" href="#fig3">Figure </a><a class="internalNav" href="#fig3">3</a>b plotting (Table S3), LJ parameters for water oxygen and intermolecular energies (Ei) of water molecules (Table S4), Ei as a function of simulation time for TIP3P and TIPS3P (Figure S1), the radius of gyration (Rg) as a function of simulation time for folded proteins (Figure S2), simulated Rg with error estimates by block averaging (Figure S3), RMSD of protein backbone from crystal structure (Figure S4), and Rg vs simulation time for the IDPs (Figure S5) (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.1c00339/suppl_file/ci1c00339_si_001.pdf">PDF</a>)Scaled nonbonded parameters and simulation MDP files (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.1c00339/suppl_file/ci1c00339_si_002.zip">ZIP</a>)This article has not yet been cited by other publications.

Pace Force Field for Protein Simulations. 1. Full Parameterization of Version 1 and Verification

Pace Force Field for Protein Simulations. 2. Folding Simulations of Peptides

Parameterization of Pace Force Field for Membrane Environment and Simulation of Helical Peptides and Helix-Helix Association

Residue-Speci fi c Force Field Based on the Protein Coil Library

Residue-specific Force Field Based on the Protein Coil Library. RSFF1: Modification of OPLS-AA/L.

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

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

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Coupled Two-Dimensional Main-Chain Torsional Potential for Protein Dynamics II: Performance and Validation.

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Significant Refinement of Protein Structure Models Using a Residue-Specific Force Field.

The Intrinsic Conformational Features of Amino Acids from a Protein Coil Library and Their Applications in Force Field Development.

Systematic Validation of Protein Force Fields against Experimental Data

A New Quantum Calibrated Force Field for Zinc-Protein Complex

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Developments and Applications of Coil-Library-Based Residue-Specific Force Fields for Molecular Dynamics Simulations of Peptides and Proteins

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

PMFF: Development of a Physics-Based Molecular Force Field for Protein Simulation and Ligand Docking

Force Field Benchmark of Amino Acids. 3. Hydration with Scaled Lennard-Jones Interactions

Optimizations of force-field parameters for protein systems with the secondary-structure stability and instability