An individual tumor harbors multiple molecular alterations that promote cell proliferation and prevent apoptosis and differentiation. Drugs that target specific molecular alterations have been introduced into personalized cancer medicine, but their effects can be modulated by the activities of other genes or molecules. Previous studies aiming to identify multiple molecular alterations for combination therapies are limited by available data. Given the recent large scale of available pharmacogenomic data, it is possible to systematically identify multiple biomarkers that contribute jointly to drug sensitivity, and to identify combination therapies for personalized cancer medicine. In this study, we used pharmacogenomic profiling data provided from two independent cohorts in a systematic in silico investigation of perturbed genes cooperatively associated with drug sensitivity. Our study predicted many pairs of molecular biomarkers that may benefit from the use of combination therapies. One of our predicted biomarker pairs, a mutation in the BRAF gene and upregulated expression of the PIM1 gene, was experimentally validated to benefit from a therapy combining BRAF inhibitor and PIM1 inhibitor in lung cancer. This study demonstrates how pharmacogenomic data can be used to systematically identify potentially cooperative genes and provide novel insights to combination therapies in personalized cancer medicine.