As discussed in this section, there does not exist a general solution to rapidly sample ligand binding modes and obtain correct populations or free energies, even in cases where binding modes differ primarily only by rotation of one or several rotatable bonds. To address this problem, we developed a hybrid MD and nonequilibrium candidate Monte Carlo based method in this work, which could potentially provide a more general approach.