Tyrosine kinase inhibitors were found to be clinically effective for treatment of patients with certain subsets of cancers carrying somatic mutations in receptor tyrosine kinases. However, the duration of clinical response is often limited, and patients ultimately develop drug resistance. Here, we use single-cell RNA sequencing to demonstrate the existence of multiple cancer cell subpopulations within cell lines, xenograft tumors and patient tumors. These subpopulations exhibit epigenetic changes and differential therapeutic sensitivity. Recurrently overrepresented ontologies in genes that are differentially expressed between drug tolerant cell populations and drug sensitive cells include epithelial-to-mesenchymal transition, epithelium development, vesicle mediated transport, drug metabolism and cholesterol homeostasis. We show analysis of identified markers using the LINCS database to predict and functionally validate small molecules that target selected drug tolerant cell populations. In combination with EGFR inhibitors, crizotinib inhibits the emergence of a defined subset of EGFR inhibitor-tolerant clones. In this study, we describe the spectrum of changes associated with drug tolerance and inhibition of specific tolerant cell subpopulations with combination agents.

Hexokinase 2 (HK2), which catalyzes the first committed step in glucose metabolism, is induced in cancer cells. HK2’s role in tumorigenesis has been attributed to its glucose kinase activity. Here, we describe a kinase independent HK2 activity, which contributes to metastasis. HK2 binds and sequesters glycogen synthase kinase 3 (GSK3) and acts as a scaffold forming a ternary complex with the regulatory subunit of protein kinase A (PRKAR1a) and GSK3β to facilitate GSK3β phosphorylation and inhibition by PKA. Thus, HK2 functions as an A-kinase anchoring protein (AKAP). Phosphorylation by GSK3β targets proteins for degradation. Consistently, HK2 increases the level and stability of GSK3 targets, MCL1, NRF2, and particularly SNAIL. In addition to GSK3 inhibition, HK2 kinase activity mediates SNAIL glycosylation, which prohibits its phosphorylation by GSK3. Finally, in mouse models of breast cancer metastasis, HK2 deficiency decreases SNAIL protein levels and inhibits SNAIL-mediated epithelial mesenchymal transition and metastasis.