BackgroundTumor heterogeneity is a main contributor of resistance to anti-cancer targeted agents though it has proven difficult to study. Unfortunately, model systems to functionally characterize and mechanistically study dynamic responses to treatment across coexisting subpopulations of cancer cells remain a missing need in oncology.MethodsUsing single cell cloning and expansion techniques, we established monoclonal cell subpopulations (MCPs) from a commercially available epidermal growth factor receptor (EGFR)-mutant non-small cell lung cancer cell line. We then used this model sensitivity to the EGFR inhibitor osimertinib across coexisting cell populations within the same tumor. Pathway-centered signaling dynamics associated with response to treatment and morphological characteristics of the MCPs were assessed using Reverse Phase Protein Microarray. Signaling nodes differentially activated in MCPs less sensitive to treatment were then pharmacologically inhibited to identify target signaling proteins putatively implicated in promoting drug resistance.ResultsMCPs demonstrated highly heterogeneous sensitivities to osimertinib. Cell viability after treatment increased > 20% compared to the parental line in selected MCPs, whereas viability decreased by 75% in other MCPs. Reduced treatment response was detected in MCPs with higher proliferation rates, EGFR L858R expression, activation of EGFR binding partners and downstream signaling molecules, and expression of epithelial-to-mesenchymal transition markers. Levels of activation of EGFR binding partners and MCPs’ proliferation rates were also associated with response to c-MET and IGFR inhibitors.ConclusionsMCPs represent a suitable model system to characterize heterogeneous biomolecular behaviors in preclinical studies and identify and functionally test biological mechanisms associated with resistance to targeted therapeutics.