AbstractMolecular Dynamics (MD) simulations have been performed on a set of rigid‐body docking poses, carried out over 25 protein–protein complexes. The results show that fully flexible relaxation increases the fraction of native contacts (NC) by up to 70% for certain docking poses. The largest increase in the fraction of NC is observed for docking poses where anchor residues are able to sample their bound conformation. For each MD simulation, structural snap‐shots were clustered and the centre of each cluster used as the MD‐relaxed docking pose. A comparison between two energy‐based scoring schemes, the first calculated for the MD‐relaxed poses, the second for energy minimized poses, shows that the former are better in ranking complexes with large hydrophobic interfaces. Furthermore, complexes with large interfaces are generally ranked well, regardless of the type of relaxation method chosen, whereas complexes with small hydrophobic interfaces remain difficult to rank. In general, the results indicate that current force‐fields are able to correctly describe direct intermolecular interactions between receptor and ligand molecules. However, these force‐fields still fail in cases where protein–protein complexes are stabilized by subtle energy contributions. Proteins 2007. © 2007 Wiley‐Liss, Inc.