Interestingly, Baumann et al. identify that the presence of waters in the binding site worsened free energy estimates for their HSP90 test case, although it offered improved convergence in results. It is therefore possible that our attempts to optimise binding site waters through MC/MD in fully interacting protein-ligand complexes may have influenced the deviations in absolute values seen in our results.