Abstract Identification of protein-ligand interaction networks on a proteome scale is crucial to address a wide range of biological problems such as correlating molecular functions to physiological processes and designing safe and efficient therapeutics. We have developed a novel computational strategy to identify ligand binding profiles of proteins across gene families and applied it to predicting protein functions, elucidating molecular mechanisms of drug adverse effects, and repositioning safe pharmaceuticals to treat different diseases.