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.

Accuracy Of Alternate Nonpolarizable Force Fields For The Determination Of Protein-Ligand Binding Affinities Dominated By Cation-Pi Interactions

Accurate Calculation of Relative Binding Free Energies Between Ligands with Different Net Charges

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Effect of Protein-Polarized Ligand Charges on Relative Protein Ligand Binding Affinities

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

Accurate calculation of absolute protein-ligand binding free energies

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

Implementing Protein-Polarized Ligand Charges into Relative Protein Ligand Binding Affinities Calculations

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

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

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant

Fast, Accurate, and Reliable Protocols for Routine Calculations of Protein–Ligand Binding Affinities in Drug Design Projects Using AMBER GPU-TI with ff14SB/GAFF

Chelator-Based Parameterization of the 12-6-4 Lennard-Jones Molecular Mechanics Potential for More Realistic Metal Ion–Protein Interactions

Fluctuations in a two-way interferometer

Rigorous Free Energy Perturbation Approach to Estimating Relative Binding Affinities Between Ligands with Multiple Protonation and Tautomeric States.

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Efficient Calculation of Protein-Ligand Binding Free Energy Using GFN Methods: the Power of the Cluster Model.

Improving the accuracy of predicting protein-ligand binding-free energy with semiempirical quantum chemistry charge.

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations