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

Avital Sharir-Ivry,Suliman Adam,Itamar Kass,Dana Krepel Zussman,Gal Masarati,Dorit Shemesh

DOI: https://doi.org/10.26434/chemrxiv-2023-h0s0k-v2

2023-12-07

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. 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 semi-empirical calculations of the molecule in the gas phase or homogenous 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 (QM/MM) 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 non-polarized 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 accuracy in RBFE calculations while it does not add significant computation time to standard parametrization procedures.

Chemistry

What problem does this paper attempt to address?

This paper mainly discusses the challenges encountered in computer-aided drug design when predicting the relative binding energies between proteins and different molecules, and how to achieve this goal through fast and accurate free energy calculation methods. The current standard classical all-atom force fields (FF) usually fix partial charges on atoms and ignore the polarization effect caused by the protein environment, which may lead to significant errors in free energy calculations. A new parameterization scheme called Protein-Induced Polarization (PIP) charge is proposed in the study to consider the polarization effect of the protein environment on ligand charges. This method determines the polarized partial charges of the ligand through quantum mechanics/molecular mechanics (QM/MM) calculations, where the protein and water molecules are modeled as point charges. Using these polarized ligand charges for the calculation of relative binding free energies (RBFEs), the results show that compared to the non-polarized standard GAFF charges, the RBFEs calculated using PIP charges are either significantly improved or at least equivalent. In the paper, the researchers used three extensively studied protein-ligand systems for experiments, comparing the performance of the PIP method with GAFF and the state-of-the-art FEP+ method. The results show that using protein polarization charges can improve the accuracy of calculating relative binding energies, especially in certain proteins, without significantly increasing computation time. In conclusion, the paper attempts to address the issue of how to more accurately simulate the impact of charge polarization effects in the protein environment on ligand binding affinity, in order to improve the accuracy and efficiency of computer prediction of the binding affinity between drug molecules and target proteins.

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

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

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

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

Pbsa_e: A Pbsa-Based Free Energy Estimator for Protein-Ligand Binding Affinity

Developing an Effective Polarizable Bond Method for Small Molecules with Application to Optimized Molecular Docking

Improving the Scoring of Protein-Ligand Binding Affinity by Including the Effects of Structural Water and Electronic Polarization.

Some practical approaches to treating electrostatic polarization of proteins.

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

Binding Free Energy Estimation for Protein-Ligand Complex Based on MM-PBSA with Various Partial Charge Models.

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

Accessing the Applicability of Polarized Protein-Specific Charge in Linear Interaction Energy Analysis

Evaluating the functional importance of conformer-dependent atomic partial charge assignment

Accurate Estimation Of Protein-Ligand Binding Free Energies Based On Geometric Restraints

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

Accurate Description of Cation-π Interactions in Proteins with a Nonpolarizable Force Field at No Additional Cost.

Calculation of Protein–ligand Binding Affinities Based on a Fragment Quantum Mechanical Method

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

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate calculation of absolute protein-ligand binding free energies

Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants.