Modern methods in genomics and high-throughput crystallography have ensured that we have access to a large and rapidly increasing, number of X-ray structures of protein-ligand complexes. A structure-based approach to drug design aims to exploit this information, but current methods are not suited to the examination of the large numbers of complexes available. We present computational tools that analyse and display multiple protein-ligand interactions and their properties in a simplified way. We illustrate how a novel binding-mode similarity metric is able to cluster 20 ligands complexed to HIV-1 reverse transcriptase into distinct groups. The properties of each cluster are then projected onto a group surface as a series of color gradients. Analysis of these surfaces reveals fundamental similarities and differences in the binding modes of these diverse compounds. In addition, the simplicity of the surface representations facilitates the transfer of information between the crystallographer, computational chemist and the chemist. We also show how two- and three-dimensional (2- and 3-D) similarities can be combined to provide enhanced understanding of 33 factor Xa inhibitor complexes. This methodology has enabled us to identify pharmaceutically relevant relationships between ligands and their binding modes that had previously been hidden in a wealth of data.