Computational Analysis of Cholangiocarcinoma Phosphoproteomes Identifies Patient-Specific Drug Targets

Cholangiocarcinoma is a form of hepatobiliary cancer with an abysmal prognosis. Despite advances in our understanding of cholangiocarcinoma pathophysiology and its genomic landscape, targeted therapies have not yet made a significant impact on its clinical management. The low response rates of targeted therapies in cholangiocarcinoma suggest that patient heterogeneity contributes to poor clinical outcome. Here we used mass spectrometry–based phosphoproteomics and computational methods to identify patient-specific drug targets in patient tumors and cholangiocarcinoma-derived cell lines. We analyzed 13 primary tumors of patients with cholangiocarcinoma with matched nonmalignant tissue and 7 different cholangiocarcinoma cell lines, leading to the identification and quantification of more than 13,000 phosphorylation sites. The phosphoproteomes of cholangiocarcinoma cell lines and patient tumors were significantly correlated. MEK1, KIT, ERK1/2, and several cyclin-dependent kinases were among the protein kinases most frequently showing increased activity in cholangiocarcinoma relative to nonmalignant tissue. Application of the Drug Ranking Using Machine Learning (DRUML) algorithm selected inhibitors of histone deacetylase (HDAC; belinostat and CAY10603) and PI3K pathway members as high-ranking therapies to use in primary cholangiocarcinoma. The accuracy of the computational drug rankings based on predicted responses was confirmed in cell-line models of cholangiocarcinoma. Together, this study uncovers frequently activated biochemical pathways in cholangiocarcinoma and provides a proof of concept for the application of computational methodology to rank drugs based on efficacy in individual patients. SIGNIFICANCE: Phosphoproteomic and computational analyses identify patient-specific drug targets in cholangiocarcinoma, supporting the potential of a machine learning method to predict personalized therapies.