PM6/COSMO optimization of the ligand's geometry within the protein prior to scoring offered no improvement to the ∆H of binding (Figure ). We note that, at least in theory, to compute accurate ligand deformation energy for the ∆H of ligand binding, the geometry of the bound state and of the free state must be optimized on the same level of theory, which apparently contradicts our observations in this study.