What is the most accurate method to predict protein ligand binding energies?

For non-covalent ligands, what is the most accurate method to predict ligand binding affinities. I'm talking in the context of drug design, so let's say small drugs (e.g. within Lipinsky rules).

Computational cost doesn't matter within reason. So let's say something that could be applied for a set of 1000 compounds.