An improved method to predict the entropy term with the MM/PBSA approach

A method is suggested to calculate improved entropies within the MM/PBSA approach (molecular mechanics combined with Poisson-Boltzmann and surface area calculations) to estimate protein-ligand binding affinities. In the conventional approach, the protein is truncated outside ~8 A from the ligand. This system is freely minimised using a distance-dependent dielectric constant (to simulate the removed protein and solvent). However, this can lead to extensive changes in the molecular geometry, giving rise to a large standard deviation in this term. In our new approach, we introduce a buffer region approximately 4 A outside the truncated protein (including solvent molecules) and keep it fixed during the minimisation. Thereby, we reduce the standard deviation by a factor of 2-4, ensuring that the entropy term no longer limits the precision of the MM/PBSA predictions. The new method is tested for the binding of seven biotin analogues to avidin, eight amidinobenzyl-indole-carboxamide inhibitors to factor Xa, and two substrates to cytochrome P450 3A4 and 2C9. It is shown that it gives more stable results and often improved predictions of the relative binding affinities.