The structural comparison of protein binding sites is increasingly important in drug design; identifying structurally similar sites can be useful for techniques such as drug repurposing, and also in a polypharmacological approach to deliberately affect multiple targets in a disease pathway, or to explain unwanted off-target effects. Once similar sites are identified, identifying local differences can aid in the design of selectivity. Such an approach moves away from the classical "one target one drug" approach and toward a wider systems biology paradigm. Here, we report a semiautomated approach, called BioGPS, that is based on the software FLAP which combines GRID Molecular Interactions Fields (MIFs) and pharmacophoric fingerprints. BioGPS comprises the automatic preparation of protein structure data, identification of binding sites, and subsequent comparison by aligning the sites and directly comparing the MIFs. Chemometric approaches are included to reduce the complexity of the resulting data on large datasets, enabling focus on the most relevant information. Individual site similarities can be analyzed in terms of their Pharmacophoric Interaction Field (PIF) similarity, and importantly the differences in their PIFs can be extracted. Here we describe the BioGPS approach, and demonstrate its applicability to rationalize off-target effects (ERα and SERCA), to classify protein families and explain polypharmacology (ABL1 kinase and NQO2), and to rationalize selectivity between subfamilies (MAP kinases p38α/ERK2 and PPARδ/PPARγ). The examples shown demonstrate a significant validation of the method and illustrate the effectiveness of the approach.