PROJECT SUMMARY Tumor heterogeneity is a major contributor to variable response and treatment failure in cancer patients. Usually, heterogeneity in cancer is thought of in terms of resistance-conferring genetic mutations that pre- exist or emerge during treatment. However, recent studies, including our own, increasingly point to non- genetic sources of heterogeneity as critical factors in the early stages of tumor response. Non-genetic mechanisms are known to underlie cellular processes such as stem cell differentiation and epithelial-to- mesenchymal transitions. In bacteria, isogenic cell populations have been shown to diversify in the absence of perturbations (e.g., drugs) into a variety of cellular phenotypes, each with differential fitness to potential stressors. This “bet hedging” strategy increases the odds that a portion of the population will survive a future, unknown challenge. We, and others, have recently hypothesized that cancer cells employ a similar survival strategy to withstand the initial onslaught of anticancer drugs. So-called “drug tolerant” cells may persist within a patient for extended periods of time before acquiring genetic resistance mutations that lead to tumor recurrence. The objective of this proposal is to uncover the molecular factors that control non-genetic heterogeneity in cancer cell populations using a combined computational and experimental approach. In Aim 1, I propose to construct a detailed kinetic model of the biochemical signaling networks that control division and death decisions in individual cancer cells. It is well established that complex biochemical networks can give rise to multiple stable equilibrium states, known as “attractors.” Each attractor corresponds to a cellular phenotype and can be conceptualized as a basin within an “epigenetic landscape.” Cells can transition between phenotypes with rates dependent upon the depths of the basins and the heights of the barriers separating them. Using a dynamical systems analysis approach, I will mathematically solve for the epigenetic landscape of the biochemical division/death model and quantify molecule signatures for all attractors. In Aim 2, using BRAF-mutant melanoma and EGFR- mutant lung cancer as in vitro model systems, I will use clonal and single-cell RNA sequencing and chromatin accessibility sequencing (ATAC-seq) to enumerate the number and molecular signatures of non-genetic phenotypic states. I will also utilize whole-exome sequencing to establish the non-genetic nature of the phenotypes and immunocompromised mouse models to validate model predictions. Differences between the experimental and in silico molecular signatures will lead to model refinement and further experimentation. Quantifying the epigenetic landscapes of cancer cells will lay the groundwork for novel therapies based on rationally modifying the landscape to favor phenotypes with increased drug sensitivity, an approach termed “targeted landscaping.” This would reduce the size ...