Search for: All records

A full understanding of cell signaling processes requires knowledge of protein structure–function relationships, protein–protein interactions, and the abilities of pathways to control phenotypes. Computational models offer a valuable framework for integrating that knowledge to predict the effects of system perturbations and interventions in health and disease. Whereas mechanistic models are well suited for understanding the biophysical basis for signal transduction and principles of therapeutic design, data-driven models are particularly suited to distill complex signaling relationships among samples and between multivariate signaling changes and phenotypes. Both approaches have limitations and provide incomplete representations of signaling biology, but their careful implementation and integration can provide new understanding for how manipulating system variables impacts cellular decisions. 

Activation of the epidermal growth factor (EGF) receptor (EGFR) at the cell surface initiates signaling through the RAS-RAF-MAPK/ERK1/2 pathway and receptor endocytosis. Whether this signaling continues from endosomes remains unclear, because RAS is predominantly located on the plasma membrane, and the localization of endogenous RAF kinases, downstream effectors of RAS, is not defined. To examine RAF localization, we labeled endogenous RAF1 with mVenus using gene editing. From 10 to 15% of RAF1-mVenus (<2000 molecules/cell), which was initially entirely cytosolic, transiently translocated to the plasma membrane after EGF stimulation. Following an early burst of translocation, the membrane-associated RAF1-mVenus was undetectable by microscopy or subcellular fractionation, and this pool was estimated to be <200 molecules per cell. In contrast, persistent EGF-dependent translocation of RAF1-mVenus to the plasma membrane was driven by the RAF inhibitor sorafenib, which increases the affinity of Ras-GTP:RAF1 interactions. RAF1-mVenus was not found in EGFR-containing endosomes under any conditions. Computational modeling of RAF1 dynamics revealed that RAF1 membrane abundance is controlled most prominently by association and dissociation rates from RAS-GTP and by RAS-GTP concentration. The model further suggested that the relatively protracted activation of the RAF-MEK1/2-ERK1/2 module, in comparison with RAF1 membrane localization, may involve multiple rounds of cytosolic RAF1 rebinding to active RAS at the membrane.