PUZZLE: a new method for automated protein docking based on surface shape complementarity.

We describe here a novel procedure for automated protein docking, based only on geometric criteria. In our algorithm we project protein surfaces into bi-dimensional matrices; the search for complementary regions is performed by detecting matching sub-matrices. An exhaustive sampling of the rotation space is made in order to analyse all the possible relative orientations of the two proteins, but nevertheless this procedure requires a relatively short processing time (3 h to 24 h cpu time on a SG4D320, depending on the complexity of the input information). When tested with co-crystallized, free components and models of components of known protein-protein complexes, the method gave very satisfactory results. The procedure selects no more than four relative orientations of the molecular components, but the correct orientation is always present among them, ranking either first or second. In more than half the cases the "wrong" solutions nevertheless correctly identify most of the residues involved in the interaction. This is remarkable also in view of the fact that the chosen test complexes (trypsin-trypsin inhibitor and antibody-lysozyme) have a very different geometry of surface complementarity: trypsin inhibitor inserts a long side-chain into the deep specificity pocket of the protease, while the interface between antibody and lysozyme is rather flat and contains buried water molecules (not included in the calculation). In order to simulate a more realistic protein docking problem, we also used a trypsin inhibitor and an anti-lysozyme antibody model in our simulations, again with satisfying results.