To initialize the non-equilibrium alchemical transitions, the first 2.256 ns of all production simulations were discarded and the equilibrium conformations were sampled every 67 ps yielding 165 conformations for each system. The extracted conformers were used to construct an equilibrium ensemble of the decoupled ligand in the protein's binding site and generate proteinligand restraints as described in the Section 1.