The software I have used in the past did not allow us to include moveable non-bound substrate ligands and attach them (pseudo-covalently) to the end of some adjacent residue. Then the problem is coming up with potentials (scoring functions) between this attached ligand and other side-chain types. Is there modern software tools for better doing this kind of searching?