Multi-range ERK responses shape the proliferative trajectory of single cells following oncogene induction.

Oncogene-induced senescence is a phenomenon in which aberrant oncogene expression causes non-transformed cells to enter a non-proliferative state. Cells undergoing oncogenic induction display phenotypic heterogeneity, with some cells senescing and others remaining proliferative. The causes of heterogeneity remain unclear. We studied the sources of heterogeneity in the responses of human epithelial cells to oncogenic BRAFV600E expression. We found that a narrow expression range of BRAFV600E generated a wide range of activities of its downstream effector ERK. In population-level and single-cell assays, ERK activity displayed a non-monotonic relationship to proliferation, with intermediate ERK activities leading to maximal proliferation. We profiled gene expression across a range of ERK activities over time and characterized four distinct ERK response classes, which we propose act in concert to generate the ERK-proliferation response. Altogether, our studies map the input-output relationships between ERK activity and proliferation, elucidating how heterogeneity can be generated during oncogene induction.

Mechanistic model of MAPK signaling reveals how allostery and rewiring contribute to drug resistance.

BRAF is prototypical of oncogenes that can be targeted therapeutically and the treatment of BRAFV600E melanomas with RAF and MEK inhibitors results in rapid tumor regression. However, drug-induced rewiring generates a drug adapted state thought to be involved in acquired resistance and disease recurrence. In this article, we study mechanisms of adaptive rewiring in BRAFV600E melanoma cells using an energy-based implementation of ordinary differential equation (ODE) modeling in combination with proteomic, transcriptomic and imaging data. We develop a method for causal tracing of ODE models and identify two parallel MAPK reaction channels that are differentially sensitive to RAF and MEK inhibitors due to differences in protein oligomerization and drug binding. We describe how these channels, and timescale separation between immediate-early signaling and transcriptional feedback, create a state in which the RAS-regulated MAPK channel can be activated by growth factors under conditions in which the BRAFV600E -driven channel is fully inhibited. Further development of the approaches in this article is expected to yield a unified model of adaptive drug resistance in melanoma.

Multiplexed and reproducible high content screening of live and fixed cells using Dye Drop.

High-throughput measurement of cells perturbed using libraries of small molecules, gene knockouts, or different microenvironmental factors is a key step in functional genomics and pre-clinical drug discovery. However, it remains difficult to perform accurate single-cell assays in 384-well plates, limiting many studies to well-average measurements (e.g., CellTiter-Glo®). Here we describe a public domain Dye Drop method that uses sequential density displacement and microscopy to perform multi-step assays on living cells. We use Dye Drop cell viability and DNA replication assays followed by immunofluorescence imaging to collect single-cell dose-response data for 67 investigational and clinical-grade small molecules in 58 breast cancer cell lines. By separating the cytostatic and cytotoxic effects of drugs computationally, we uncover unexpected relationships between the two. Dye Drop is rapid, reproducible, customizable, and compatible with manual or automated laboratory equipment. Dye Drop improves the tradeoff between data content and cost, enabling the collection of information-rich perturbagen-response datasets.

Structural basis for SHOC2 modulation of RAS signalling.

The RAS-RAF pathway is one of the most commonly dysregulated in human cancers1-3. Despite decades of study, understanding of the molecular mechanisms underlying dimerization and activation4 of the kinase RAF remains limited. Recent structures of inactive RAF monomer5 and active RAF dimer5-8 bound to 14-3-39,10 have revealed the mechanisms by which 14-3-3 stabilizes both RAF conformations via specific phosphoserine residues. Prior to RAF dimerization, the protein phosphatase 1 catalytic subunit (PP1C) must dephosphorylate the N-terminal phosphoserine (NTpS) of RAF11 to relieve inhibition by 14-3-3, although PP1C in isolation lacks intrinsic substrate selectivity. SHOC2 is as an essential scaffolding protein that engages both PP1C and RAS to dephosphorylate RAF NTpS11-13, but the structure of SHOC2 and the architecture of the presumptive SHOC2-PP1C-RAS complex remain unknown. Here we present a cryo-electron microscopy structure of the SHOC2-PP1C-MRAS complex to an overall resolution of 3 Å, revealing a tripartite molecular architecture in which a crescent-shaped SHOC2 acts as a cradle and brings together PP1C and MRAS. Our work demonstrates the GTP dependence of multiple RAS isoforms for complex formation, delineates the RAS-isoform preference for complex assembly, and uncovers how the SHOC2 scaffold and RAS collectively drive specificity of PP1C for RAF NTpS. Our data indicate that disease-relevant mutations affect complex assembly, reveal the simultaneous requirement of two RAS molecules for RAF activation, and establish rational avenues for discovery of new classes of inhibitors to target this pathway.

Overcoming differential tumor penetration of BRAF inhibitors using computationally guided combination therapy.

BRAF-targeted kinase inhibitors (KIs) are used to treat malignancies including BRAF-mutant non-small cell lung cancer, colorectal cancer, anaplastic thyroid cancer, and, most prominently, melanoma. However, KI selection criteria in patients remain unclear, as are pharmacokinetic/pharmacodynamic (PK/PD) mechanisms that may limit context-dependent efficacy and differentiate related drugs. To address this issue, we imaged mouse models of BRAF-mutant cancers, fluorescent KI tracers, and unlabeled drug to calibrate in silico spatial PK/PD models. Results indicated that drug lipophilicity, plasma clearance, faster target dissociation, and, in particular, high albumin binding could limit dabrafenib action in visceral metastases compared to other KIs. This correlated with retrospective clinical observations. Computational modeling identified a timed strategy for combining dabrafenib and encorafenib to better sustain BRAF inhibition, which showed enhanced efficacy in mice. This study thus offers principles of spatial drug action that may help guide drug development, KI selection, and combination.

Receptor-Driven ERK Pulses Reconfigure MAPK Signaling and Enable Persistence of Drug-Adapted BRAF-Mutant Melanoma Cells.

Targeted inhibition of oncogenic pathways can be highly effective in halting the rapid growth of tumors but often leads to the emergence of slowly dividing persister cells, which constitute a reservoir for the selection of drug-resistant clones. In BRAFV600E melanomas, RAF and MEK inhibitors efficiently block oncogenic signaling, but persister cells emerge. Here, we show that persister cells escape drug-induced cell-cycle arrest via brief, sporadic ERK pulses generated by transmembrane receptors and growth factors operating in an autocrine/paracrine manner. Quantitative proteomics and computational modeling show that ERK pulsing is enabled by rewiring of mitogen-activated protein kinase (MAPK) signaling: from an oncogenic BRAFV600E monomer-driven configuration that is drug sensitive to a receptor-driven configuration that involves Ras-GTP and RAF dimers and is highly resistant to RAF and MEK inhibitors. Altogether, this work shows that pulsatile MAPK activation by factors in the microenvironment generates a persistent population of melanoma cells that rewires MAPK signaling to sustain non-genetic drug resistance.

Real-time metabolome profiling of the metabolic switch between starvation and growth.

Metabolic systems are often the first networks to respond to environmental changes, and the ability to monitor metabolite dynamics is key for understanding these cellular responses. Because monitoring metabolome changes is experimentally tedious and demanding, dynamic data on time scales from seconds to hours are scarce. Here we describe real-time metabolome profiling by direct injection of living bacteria, yeast or mammalian cells into a high-resolution mass spectrometer, which enables automated monitoring of about 300 compounds in 15-30-s cycles over several hours. We observed accumulation of energetically costly biomass metabolites in Escherichia coli in carbon starvation-induced stationary phase, as well as the rapid use of these metabolites upon growth resumption. By combining real-time metabolome profiling with modeling and inhibitor experiments, we obtained evidence for switch-like feedback inhibition in amino acid biosynthesis and for control of substrate availability through the preferential use of the metabolically cheaper one-step salvaging pathway over costly ten-step de novo purine biosynthesis during growth resumption.