Case-specific performance of MM-PBSA, MM-GBSA, and SIE in virtual screening

In drug discovery the reliable prediction of binding free energies is of crucial importance. Methods that combine molecular mechanics force fields with continuum solvent models have become popular because of their high accuracy and relatively good computational efficiency. In this research we studied the performance of molecular mechanics generalized Born surface area (MM-GBSA), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and solvated interaction energy (SIE) both in their virtual screening efficiency and their ability to predict experimentally determined binding affinities for five different protein targets. The protein-ligand complexes were derived with two different approaches important in virtual screening: molecular docking and ligand-based similarity search methods. The results show significant differences between the different binding energy calculation methods. However, the length of the molecular dynamics simulation was not of crucial importance for accuracy of results.