Analysis of the alchemical simulations was achieved using the alchemlyb v0.3.0 library 78 and the MBAR estimator as implemented in pyMBAR v3.0.3 . The first 1 ns of each production window was discarded as extra equilibration time and samples were decorrelated based on the derivative of the potential with respect to λ up to a maximum frequency of 100 ps -1 .