Abstract:A major challenge in computer-aided drug design is predicting relative binding energies of different molecules to a target protein using fast and accurate free-energy calculation methods. Free-energy calculations are primarily computed by utilizing classical molecular dynamics simulations based on all-atom force fields (FF) to model the interactions in the system. The present standard classical all-atom FFs contain fixed partial charges on the atoms, and hence electrostatic interactions are modeled between them. The parametrization process to determine these partial charges usually relies on quantum mechanics or semiempirical calculations of the molecule in the gas phase or homogeneous water surrounding. These present standard parametrization schemes of the partial charges neglect, therefore, polarization effects from the protein surrounding. The absence of protein polarization effects can lead to significant errors in free-energy calculations in proteins. We present a parametrization scheme for the partial charges of ligands, named protein-induced polarization (PIP) charges, which account for the electrostatic polarization due to the protein surrounding. The scheme involves single-point quantum mechanics/molecular mechanics calculations of the ligand charges in the protein/water surrounding. Using PIP ligand partial charges, we have calculated the relative binding free energies (RBFEs) of well-studied protein-ligand systems. We show here that RBFEs computed with PIP charges are either significantly improved or at least comparable to those computed with nonpolarized standard GAFF charges. Overall, we present a simple-to-use parametrization scheme to include protein polarization in any type of binding free-energy calculations. The parametrization scheme increases the accuracy in RBFE calculations, while it does not add significant computation time to standard parametrization procedures.

Impact of force field polarization on the collective motions of proteins

Protein's Native Structure is Dynamically Stabilized by Electronic Polarization

Effect of Polarization on HIV-1protease and Fluoro-Substituted Inhibitors Binding Energies by Large Scale Molecular Dynamics Simulations

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

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

Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields

Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy

Intra-Protein Hydrogen Bonding Is Dynamically Stabilized By Electronic Polarization

Effective Inclusion of Electronic Polarization Improves the Description of Electrostatic Interactions: The prosECCo75 Biomolecular Force Field

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

Polarisable force fields: what do they add in biomolecular simulations?

Molecular Fluctuations as a Ruler of Force-Induced Protein Conformations

Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data

Accounting for electronic polarization in nonpolarizable force fields

RNA-protein complexes and force field polarizability

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field

Accurate biomolecular simulations account for electronic polarization

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

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

Systematic Validation of Protein Force Fields against Experimental Data

Non-Aqueous Ion Pairing Exemplifies the Case for Including Electronic Polarization in Molecular Dynamics Simulations