To aid in the problem of binding mode determination and binding mode sampling, here, we developed a new simulation technique based on MD and NCMC to sample binding modes of flexible ligands, and incorporated this into our BLUES package for binding mode sampling. In our approach, we propagate the system using standard MD with regular NCMC move proposals in between to overcome energy barriers.