For the relative binding free energy (∆∆G) values we used the raw values without any cycle closure correction as they reflect better potential shortcomings of force fields. For the pmx results, we calculated the ∆∆G values as averages over three repeats and the standard deviation across the repeats was used as an error estimate.