Our findings demonstrate that traditional flexible docking often produces highly strained conformer geometries. The final hit lists generated using low-energy, physically feasible QM conformers differ significantly from those obtained using traditional flexible ligand docking with the same intermolecular scoring function. The computational costs associated with these new methods have also been discussed in detail.