AbstractIdentifying protein binding sites provides important clues to the function of a protein. Experimental methods to identify the binding sites such as determining the crystal structures of protein complexes are extremely laborious and expensive. Here, we present a computational technique called spatial aggregation propensity (SAP) based on molecular simulations to predict protein binding sites. We apply this technique to two model proteins, an IgG1 antibody and epidermal growth factor receptor (EGFR) and demonstrate that SAP predicts protein binding regions with very good accuracy. In the case of the IgG1 antibody, SAP accurately predicts binding regions with the Fc‐receptor, protein‐A, and protein‐G. For EGFR, SAP accurately predicts binding regions with EGF, TGFα, and with another EGFR. The resolution of SAP is varied to obtain a detailed picture of these binding sites. We also show that some of these binding sites overlap with protein self‐aggregation prone regions. We demonstrate how SAP analysis can be used to engineer the protein to remove unfavorable aggregation prone regions without disturbing protein binding regions. The SAP technique could be also used to predict the yet unknown binding sites of numerous proteins, thereby providing clues to their function. Proteins 2011. © 2010 Wiley‐Liss, Inc.