Oversampling Free Energy Perturbation Simulation in Determination of the Ligand-Binding Free Energy

Determination of the ligand-binding affinity is an extremely interesting problem. Normally, the free energy perturbation (FEP) method provides an appropriate result. However, it is of great interest to improve the accuracy and precision of this method. In this context, temperature replica exchange molecular dynamics implementation of the FEP computational approach, which we call replica exchange free energy perturbation (REP) was proposed. In particular, during REP simulations, the system can easily escape from being trapped in local minima by exchanging configurations with high temperatures, resulting in significant improvement in the accuracy and precision of protein–ligand binding affinity calculations. The distribution of the decoupling free energy was enlarged, and its mean values were decreased. This results in changes in the magnitude of the calculated binding free energies as well as in alteration in the binding mechanism. Moreover, the REP correlation coefficient with respect to experiment ( R^REP = 0.85 ± 0.15 ) is significantly boosted in comparison with the FEP one ( R^FEP = 0.64 ± 0.30 ). Furthermore, the root-mean-square error (RMSE) of REP is also smaller than FEP, RMSE^REP = 4.28 ± 0.69 versus RMSE^FEP = 5.80 ± 1.11 kcal/mol, respectively. © 2019 Wiley Periodicals, Inc.