A novel hybrid optimisation method, called quantum stochastic tunnelling, is introduced. This method combines path integral Monte Carlo and stochastic tunnelling approaches. We have implemented this method within a flexible ligand docking program and found it to be successful in exploring the complex potential energy surfaces of protein-ligand complexes and finding the correct crystallographic ligand binding modes. (C) 2003 Elsevier Science B.V. All rights reserved.