Key to the fragment optimization process is the need to accurately capture the changes in affinity that are associated with a given set of chemical modifications. Due to the weakly binding nature of fragments, this has proven to be a challenging task, despite recent advancements in leveraging experimental and computational methods. In this work, we evaluate the use of Absolute Binding Free Energy (ABFE) calculations in guiding fragment optimization decisions, retrospectively calculating binding free energies for 59 ligands across 4 fragment elaboration campaigns. We first demonstrate that ABFEs can be used to accurately rank fragment-sized binders with an overall Spearman’s r of 0.89 and a Kendall τ of 0.67, although often deviating from experiment in absolute free energy values with an RMSE of 2.75 kcal/mol. We then also show that in several cases, retrospective fragment optimization decisions can be supported by the ABFE calculations. Cases that were not supported were often limited by large uncertainties in the free energy estimates, however generally the right direction in ΔΔG is still observed. Comparing against cheaper endpoint methods, namely Nwat-MM/GBSA, we find that ABFEs offer better outcomes in ranking binders, improving correlation metrics, although a similar confidence in retrospective synthetic decisions is achieved. Our results indicate that ABFE calculations are currently at the level of accuracy that can be usefully employed to gauge which fragment elaborations are likely to offer the best gains in affinity.
Full data sets are now released on zenodo and github. A more detailed section on HSP90 discussing the issues of loop motions is also included.
Evaluating the use of Absolute Binding Free Energy in the fragment optimization process