Effect of electrostatic polarization and bridging water on CDK2–ligand binding affinities calculated using a highly efficient interaction entropy method

Abstract

A new highly efficient interaction entropy (IE) method combined with the polarized protein-specific charge (PPC) force field is employed to investigate the interaction mechanism of CDK2-ligand binding and the effect of the bridging water. Our result shows that the computed binding free energies for five CDK2-ligand complexes using the IE method have a significantly linear correlation with the experimentally measured values with a correlation coefficient of 0.98 in consideration of the bridging water under the PPC force field. And the correlation coefficient is found to be slightly weaker with a value of 0.95 using the traditional normal mode (Nmode) method for calculation of entropy change. Importantly, the rank of the predicted binding free energies is significantly consistent with the experimental rank based on the IE method calculated entropy change using the PPC force field. However, without including the bridging water under PPC simulation, the correlation coefficient is below 0.83. For comparison, the result obtained from the simulation using the nonpolarized AMBER force field gives a much weaker correlation with the correlation coefficients of 0.44 and 0.45 using the Nmode method and IE method, due to the lack of electrostatic polarization. Furthermore, hydrogen bond analysis indicates that the bridging water makes a significant contribution to mediating the hydrogen bond network of protein-ligand binding and stabilizing the complex structure. The current study demonstrates that the new IE method is superior to the standard Nmode method in computing the binding free energy. And our results also emphasize the importance of electronic polarization and bridging water in MD simulations and free energy calculations.

Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV‐1 protease by enhanced conformational sampling

“…Molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) and molecular mechanics generalized Born surface area (MM‐GBSA) methods are fast methods to predict binding free energies of inhibitors to proteins. Based on the systematic comparison from Hou's group, MM‐GBSA method was selected to compute binding free energies of three inhibitors to the WT and MDR PRs using the following equation in which ΔH and TΔ S represent the changes in entropy and enthalpy induced by inhibitor bindings.…”

Insight into selective mechanism of class of I‐BRD9 inhibitors toward BRD9 based on molecular dynamics simulations

“…To probe binding modes of inhibitors to BRDs, we carried out detailed insight into binding selectivity of I‐BRD9 toward BRD9 and BRD7. Meanwhile, three derivations N1D, TVU, and 5V2 of I‐BRD9 were selected to clarify their selective mechanism toward BRD9 and BRD4 by using molecular dynamic (MD) simulation and binding free energy calculations . Nine complexes composed of these four inhibitors and three proteins are H1B‐BRD9/BRD7/BRD4, N1D‐BRD9/BRD4, TVU‐BRD9/BRD4, and 5V2‐BRD9/BRD4.…”

Exploring the Reasons for Decrease in Binding Affinity of HIV-2 Against HIV-1 Protease Complex Using Interaction Entropy Under Polarized Force Field

“…So, it is promising to further improve the accuracy of MM/PBSA methods by improving the accuracy of the calculation of entropy change using a highly efficient and reliable method. In this paper, a novel and effective method, namely interaction entropy (IE) (Duan et al, 2016b, 2017), is used to calculate the entropy change. The method is more theoretically rigorous and time-saving than traditional Nmode method (Cong et al, 2017; Yan et al, 2017; Liu et al, 2018; Qiu et al, 2018; Song et al, 2018).…”