Abstract:Entropy effects play an important role in drug-target interactions, but the entropic contribution to ligand-binding affinity is often neglected by end-point binding free energy calculation methods, such as MM/GBSA and MM/PBSA, due to the expensive computational cost of normal mode analysis (NMA). Here, we systematically investigated entropy effects on the prediction power of MM/GBSA and MM/PBSA using >1500 protein-ligand systems and six representative AMBER force fields. Two computationally efficient methods, including NMA based on truncated structures and the interaction entropy approach, were used to estimate the entropic contributions to ligand-target binding free energies. In terms of the overall accuracy, we found that, for the minimized structures, in most cases the inclusion of the conformational entropies predicted by truncated NMA (enthalpynmode_min_9Å) compromises the overall accuracy of MM/GBSA and MM/PBSA compared with the enthalpies calculated based on the minimized structures (enthalpymin). However, for the MD trajectories, the binding free energies can be improved by the inclusion of the conformation entropies predicted by either truncated-NMA for a relatively high dielectric constant (εin = 4) or the interaction entropy method for εin = 1-4. In terms of reproducing the absolute binding free energies, the binding free energies estimated by including the truncated-NMA entropies based on the MD trajectories (ΔGnmode_md_9Å) give the lowest average absolute deviations against the experimental data among all the tested strategies for both MM/GBSA and MM/PBSA. Although the inclusion of the truncated NMA based on the MD trajectories (ΔGnmode_md_9Å) for a relatively high dielectric constant gave the overall best result and the lowest average absolute deviations against the experimental data (for the ff03 force field), it needs too much computational time. Alternatively, considering that the interaction entropy method does not incur any additional computational cost and can give comparable (at high dielectric constant, εin = 4) or even better (at low dielectric constant, εin = 1-2) results than the truncated-NMA entropy (ΔGnmode_md_9Å), the interaction entropy approach is recommended to estimate the entropic component for MM/GBSA and MM/PBSA based on MD trajectories, especially for a diverse dataset. Furthermore, we compared the predictions of MM/GBSA with six different AMBER force fields. The results show that the ff03 force field (ff03 for proteins and gaff with AM1-BCC charges for ligands) performs the best, but the predictions given by the tested force fields are comparable, implying that the MM/GBSA predictions are not very sensitive to force fields.

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

An Accurate Free Energy Estimator: Based on MM/PBSA Combined with Interaction Entropy for Protein-Ligand Binding Affinity.

Interaction Entropy: A New Paradigm for Highly Efficient and Reliable Computation of Protein-Ligand Binding Free Energy.

Interaction Entropy for Computational Alanine Scanning in Protein–protein Binding

Interaction Entropy for Computational Alanine Scanning.

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Improving the Performance of MM/PBSA in Protein-Protein Interactions Via the Screening Electrostatic Energy.

Interaction Entropy for Protein-Protein Binding.

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

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

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

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

Achieving Accurate Standard Protein-Protein Binding Free Energy Calculations through the Geometrical Route and Ergodic Sampling

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

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

Assessing the Performance of MM/PBSA and MM/GBSA Methods. 7. Entropy Effects on the Performance of End-Point Binding Free Energy Calculation Approaches.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Quantifying Protein-Protein Binding Energy And Entropy Using Molecular Dynamics Simulations

How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations?

Assessing the performance of the MM/PBSA and MM/GBSA methods. 10. Impacts of enhanced sampling and variable dielectric model on protein-protein Interactions