To determine the relative binding free energy of a specific ligand for each of its identified binding sites, we use two protocols: ESMACS, 66 and the so-called "direct" binding free energy calculation method. By running the direct protocol, we also derive insights into the reproducibility of expectation values that are computed from "converged" simulations that are multiple microseconds in length.