Computational analysis of the regulation of EGFR by protein tyrosine phosphatases

Systems Biology of the Cancer Cell

Questions in cancer have engaged systems biologists for decades. During that time, the quantity of molecular data has exploded, but the need for abstractions, formal models, and simplifying insights has remained the same. This review brings together classic breakthroughs and recent findings in the field of cancer systems biology, focusing on cancer cell pathways for tumorigenesis and therapeutic response. Cancer cells mutate and transduce information from their environment to alter gene expression, metabolism, and phenotypic states. Understanding the molecular architectures that make each of these steps possible is a long-term goal of cancer systems biology pursued by iterating between quantitative models and experiments. We argue that such iteration is the best path to deploying targeted therapies intelligently so that each patient receives the maximum benefit for their cancer.

Evaluation of information flows in the RAS-MAPK system using transfer entropy measurements

The RAS-MAPK system plays an important role in regulating various cellular processes, including growth, differentiation, apoptosis, and transformation. Dysregulation of this system has been implicated in genetic diseases and cancers affecting diverse tissues. To better understand the regulation of this system, we employed information flow analysis based on transfer entropy (TE) between the activation dynamics of two key elements in cells stimulated with EGF: SOS, a guanine nucleotide exchanger for the small GTPase RAS, and RAF, a RAS effector serine/threonine kinase. TE analysis allows for model-free assessment of the timing, direction, and strength of the information flow regulating the system response. We detected significant amounts of TE in both directions between SOS and RAF, indicating feedback regulation. Importantly, the amount of TE did not simply follow the input dose or the intensity of the causal reaction, demonstrating the uniqueness of TE. TE analysis proposed regulatory networks containing multiple tracks and feedback loops and revealed temporal switching in the reaction pathway primarily responsible for reaction control. This proposal was confirmed by the effects of an MEK inhibitor on TE. Furthermore, TE analysis identified the functional disorder of a SOS mutation associated with Noonan syndrome, a human genetic disease, of which the pathogenic mechanism has not been precisely known yet. TE assessment holds significant promise as a model-free analysis method of reaction networks in molecular pharmacology and pathology.

Mathematical Modeling and Inference of Epidermal Growth Factor-Induced Mitogen-Activated Protein Kinase Cell Signaling Pathways

The mitogen-activated protein kinase (MAPK) pathway is an important intracellular signaling cascade that plays a key role in various cellular processes. Understanding the regulatory mechanisms of this pathway is essential for developing effective interventions and targeted therapies for related diseases. Recent advances in single-cell proteomic technologies have provided unprecedented opportunities to investigate the heterogeneity and noise within complex, multi-signaling networks across diverse cells and cell types. Mathematical modeling has become a powerful interdisciplinary tool that bridges mathematics and experimental biology, providing valuable insights into these intricate cellular processes. In addition, statistical methods have been developed to infer pathway topologies and estimate unknown parameters within dynamic models. This review presents a comprehensive analysis of how mathematical modeling of the MAPK pathway deepens our understanding of its regulatory mechanisms, enhances the prediction of system behavior, and informs experimental research, with a particular focus on recent advances in modeling and inference using single-cell proteomic data.

An integrated mechanistic and data-driven computational model predicts cell responses to high- and low-affinity EGFR ligands

The biophysical properties of ligand binding heavily influence the ability of receptors to specify cell fates. Understanding the rules by which ligand binding kinetics impact cell phenotype is challenging, however, because of the coupled information transfers that occur from receptors to downstream signaling effectors and from effectors to phenotypes. Here, we address that issue by developing an integrated mechanistic and data-driven computational modeling platform to predict cell responses to different ligands for the epidermal growth factor receptor (EGFR). Experimental data for model training and validation were generated using MCF7 human breast cancer cells treated with the high- and low-affinity ligands epidermal growth factor (EGF) and epiregulin (EREG), respectively. The integrated model captures the unintuitive, concentration-dependent abilities of EGF and EREG to drive signals and phenotypes differently, even at similar levels of receptor occupancy. For example, the model correctly predicts the dominance of EREG over EGF in driving a cell differentiation phenotype through AKT signaling at intermediate and saturating ligand concentrations and the ability of EGF and EREG to drive a broadly concentration-sensitive migration phenotype through cooperative ERK and AKT signaling. Parameter sensitivity analysis identifies EGFR endocytosis, which is differentially regulated by EGF and EREG, as one of the most important determinants of the alternative phenotypes driven by different ligands. The integrated model provides a new platform to predict how phenotypes are controlled by the earliest biophysical rate processes in signal transduction and may eventually be leveraged to understand receptor signaling system performance depends on cell context.

One-sentence summary

Integrated kinetic and data-driven EGFR signaling model identifies the specific signaling mechanisms that dictate cell responses to EGFR activation by different ligands.

Visualization of Concentration Gradients and Colloidal Dynamics under Electrodiffusiophoresis

In this work, we present an experimental study of the dynamics of charged colloids under direct currents and gradients of chemical species (electrodiffusiophoresis). In our approach, we simultaneously visualize the development of concentration polarization and the ensuing dynamics of charged colloids near electrodes. With the aid of confocal microscopy and fluorescent probes, we show that the passage of current through water confined between electrodes, separated about a hundred microns, results in significant pH gradients. Depending on the current density and initial conditions, steep pH gradients develop, thus becoming a significant factor in the behavior of charged colloids. Furthermore, we show that steep pH gradients induce the focusing of charged colloids away from both electrodes. Our results provide the experimental basis for further development of models of electrodiffusiophoresis and the design of non-equilibrium strategies for the fabrication of advanced materials.

Ten Approaches That Improve Immunostaining: A Review of the Latest Advances for the Optimization of Immunofluorescence

Immunostaining has emerged as one of the most common and valuable techniques that allow the localization of proteins at a quantitative level within cells and tissues using antibodies coupled to enzymes, fluorochromes, or colloidal nanogold particles. The application of fluorochromes during immunolabeling is referred to as immunofluorescence, a method coupled to widefield or confocal microscopy and extensively applied in basic research and clinical diagnosis. Notwithstanding, there are still disadvantages associated with the application of this technique due to technical challenges in the process, such as sample fixation, permeabilization, antibody incubation times, and fluid exchange, etc. These disadvantages call for continuous updates and improvements to the protocols extensively described in the literature. This review contributes to protocol optimization, outlining 10 current methods for improving sample processing in different stages of immunofluorescence, including a section with further recommendations. Additionally, we have extended our own antibody signal enhancer method, which was reported to significantly increase antibody signals and is useful for cervical cancer detection, to improve the signals of fluorochrome-conjugated staining reagents in fibrous tissues. In summary, this review is a valuable tool for experienced researchers and beginners when planning or troubleshooting the immunofluorescence assay.

Exendin-4 restores airway mucus homeostasis through the GLP1R-PKA-PPARγ-FOXA2-phosphatase signaling

Goblet cell hyperplasia and metaplasia and excessive mucus are prominent pathologies of chronic airway diseases such as chronic obstructive pulmonary disease (COPD), cystic fibrosis (CF), and chronic bronchitis. Chronic infection by respiratory pathogens, including Pseudomonas aeruginosa, exacerbates cyclical proinflammatory responses and mucus hypersecretion. P. aeruginosa and its virulence factor pyocyanin contribute to these pathologies by inhibiting FOXA2, a key transcriptional regulator of mucus homeostasis, through activation of antagonistic signaling pathways EGFR-AKT/ERK1/2 and IL-4/IL-13-STAT6-SPDEF. However, FOXA2-targeted therapy has not been previously explored. Here, we examined the feasibility of repurposing the incretin mimetic Exendin-4 to restore FOXA2-mediated airway mucus homeostasis. We have found that Exendin-4 restored FOXA2 expression, attenuated mucin production in COPD and CF-diseased airway cells, and reduced mucin and P. aeruginosa burden in mouse lungs. Mechanistically, Exendin-4 activated the GLP1R-PKA-PPAR-γ-dependent phosphatases PTEN and PTP1B, which inhibited key kinases within both EGFR and STAT6 signaling cascades. Our results may lead to the repurposing of Exendin-4 and other incretin mimetics to restore FOXA2 function and ultimately regulate excessive mucus in diseased airways.

Disrupting the transmembrane domain-mediated oligomerization of protein tyrosine phosphatase receptor J inhibits EGFR-driven cancer cell phenotypes

Receptor protein tyrosine phosphatases (RPTPs) play critical regulatory roles in mammalian signal transduction. However, the structural basis for the regulation of their catalytic activity is not fully understood, and RPTPs are generally not therapeutically targetable. This knowledge gap is partially due to the lack of known natural ligands or selective agonists of RPTPs. Contrary to what is known from structure-function studies of receptor tyrosine kinases (RTKs), RPTP activities have been reported to be suppressed by dimerization, which may prevent RPTPs from accessing their RTK substrates. We report here that homodimerization of protein tyrosine phosphatase receptor J (PTPRJ, also known as DEP-1) is regulated by specific transmembrane (TM) residues. We found that disrupting these interactions destabilizes homodimerization of full-length PTPRJ in cells, reduces the phosphorylation of the known PTPRJ substrate epidermal growth factor receptor (EGFR) and of other downstream signaling effectors, antagonizes EGFR-driven cell phenotypes, and promotes substrate access. We demonstrate these observations in human cancer cells using mutational studies and identified a peptide that binds to the PTPRJ TM domain and represents the first example of an allosteric agonist of RPTPs. The results of our study provide fundamental structural and functional insights into how PTPRJ activity is tuned by TM interactions in cells. Our findings also open up opportunities for developing peptide-based agents that could be used as tools to probe RPTPs' signaling mechanisms or to manage cancers driven by RTK signaling.

GP78 Cooperates with Dual-Specificity Phosphatase 1 To Stimulate Epidermal Growth Factor Receptor-Mediated Extracellular Signal-Regulated Kinase Signaling

GP78 is an autocrine motility factor (AMF) receptor (AMFR) with E3 ubiquitin ligase activity that plays a significant role in tumor cell proliferation, motility, and metastasis. Aberrant extracellular signal-regulated kinase (ERK) activation via receptor tyrosine kinases promotes tumor proliferation and invasion. The activation of GP78 leads to ERK activation, but its underlying mechanism is not fully understood. Here, we show that GP78 is required for epidermal growth factor receptor (EGFR)-mediated ERK activation. On one hand, GP78 interacts with and promotes the ubiquitination and subsequent degradation of dual-specificity phosphatase 1 (DUSP1), an endogenous negative regulator of mitogen-activated protein kinases (MAPKs), resulting in ERK activation. On the other hand, GP78 maintains the activation status of EGFR, as evidenced by the fact that EGF fails to induce EGFR phosphorylation in GP78-deficient cells. By the regulation of both EGFR and ERK activation, GP78 promotes cell proliferation, motility, and invasion. Therefore, this study identifies a previously unknown signaling pathway by which GP78 stimulates ERK activation via DUSP1 degradation to mediate EGFR-dependent cancer cell proliferation and invasion.

A conformational sensor based on genetic code expansion reveals an autocatalytic component in EGFR activation

Epidermal growth factor receptor (EGFR) activation by growth factors (GFs) relies on dimerization and allosteric activation of its intrinsic kinase activity, resulting in trans-phosphorylation of tyrosines on its C-terminal tail. While structural and biochemical studies identified this EGF-induced allosteric activation, imaging collective EGFR activation in cells and molecular dynamics simulations pointed at additional catalytic EGFR activation mechanisms. To gain more insight into EGFR activation mechanisms in living cells, we develop a Förster resonance energy transfer (FRET)-based conformational EGFR indicator (CONEGI) using genetic code expansion that reports on conformational transitions in the EGFR activation loop. Comparing conformational transitions, self-association and auto-phosphorylation of CONEGI and its Y845F mutant reveals that Y₈₄₅ phosphorylation induces a catalytically active conformation in EGFR monomers. This conformational transition depends on EGFR kinase activity and auto-phosphorylation on its C-terminal tail, generating a looped causality that leads to autocatalytic amplification of EGFR phosphorylation at low EGF dose.

Upon ligand binding epidermal growth factor receptor (EGFR) dimerizes and activates its intrinsic kinase to auto-phosphorylate EGFR. Here, the authors engineer and image a FRET-based conformational EGFR indicator which reveals that activation loop phosphorylation induces a catalytically active conformation in EGFR monomers.

Phosphatases control PKA-dependent functional microdomains at the outer mitochondrial membrane

The selective phosphorylation of spatially distinct PKA targets is key for the pleiotropy of the cAMP cascade. This characteristic of the pathway is currently attributed to the ability of phosphodiesterases or adenylate cyclases to create subcellular sites (microdomains) where the concentration of cAMP is distinct from that of the surrounding areas. The role of phosphatases in this process has not been tested. Here we show that limited access of phosphatases to the PKA targets present at the outer mitochondrial membrane generates distinct microdomains of PKA phosphorylated proteins despite there being no differences in the local cAMP levels. These results describe an alternative mechanism capable of generating functional cAMP/PKA-dependent microdomains and may be extrapolated to the compartmentalization of other kinase-dependent events.

Evidence supporting the heterogeneity in cAMP and PKA signaling is rapidly accumulating and has been largely attributed to the localization or activity of adenylate cyclases, phosphodiesterases, and A-kinase–anchoring proteins in different cellular subcompartments. However, little attention has been paid to the possibility that, despite homogeneous cAMP levels, a major heterogeneity in cAMP/PKA signaling could be generated by the spatial distribution of the final terminators of this cascade, i.e., the phosphatases. Using FRET-based sensors to monitor cAMP and PKA-dependent phosphorylation in the cytosol and outer mitochondrial membrane (OMM) of primary rat cardiomyocytes, we demonstrate that comparable cAMP increases in these two compartments evoke higher levels of PKA-dependent phosphorylation in the OMM. This difference is most evident for small, physiological increases of cAMP levels and with both OMM-located probes and endogenous OMM proteins. We demonstrate that this disparity depends on differences in the rates of phosphatase-dependent dephosphorylation of PKA targets in the two compartments. Furthermore, we show that the activity of soluble phosphatases attenuates PKA-driven activation of the cAMP response element-binding protein while concurrently enhancing PKA-dependent mitochondrial elongation. We conclude that phosphatases can sculpt functionally distinct cAMP/PKA domains even in the absence of gradients or microdomains of this messenger. We present a model that accounts for these unexpected results in which the degree of PKA-dependent phosphorylation is dictated by both the subcellular distribution of the phosphatases and the different accessibility of membrane-bound and soluble phosphorylated substrates to the cytosolic enzymes.

Integrative meta-modeling identifies endocytic vesicles, late endosome and the nucleus as the cellular compartments primarily directing RTK signaling

Recently, intracellular receptor signaling has been identified as a key component mediating cell responses for various receptor tyrosine kinases (RTKs). However, the extent each endocytic compartment (endocytic vesicle, early endosome, recycling endosome, late endosome, lysosome and nucleus) contributes to receptor signaling has not been quantified. Furthermore, our understanding of endocytosis and receptor signaling is complicated by cell- or receptor-specific endocytosis mechanisms. Therefore, towards understanding the differential endocytic compartment signaling roles, and identifying how to achieve signal transduction control for RTKs, we delineate how endocytosis regulates RTK signaling. We achieve this via a meta-analysis across eight RTKs, integrating computational modeling with experimentally derived cell (compartment volume, trafficking kinetics and pH) and ligand-receptor (ligand/receptor concentration and interaction kinetics) physiology. Our simulations predict the abundance of signaling from eight RTKs, identifying the following hierarchy in RTK signaling: PDGFRβ > IGFR1 > EGFR > PDGFRα > VEGFR1 > VEGFR2 > Tie2 > FGFR1. We find that endocytic vesicles are the primary cell signaling compartment; over 43% of total receptor signaling occurs within the endocytic vesicle compartment for these eight RTKs. Mechanistically, we found that high RTK signaling within endocytic vesicles may be attributed to their low volume (5.3 × 10^-19 L) which facilitates an enriched ligand concentration (3.2 μM per ligand molecule within the endocytic vesicle). Under the analyzed physiological conditions, we identified extracellular ligand concentration as the most sensitive parameter to change; hence the most significant one to modify when regulating absolute compartment signaling. We also found that the late endosome and nucleus compartments are important contributors to receptor signaling, where 26% and 18%, respectively, of average receptor signaling occurs across the eight RTKs. Conversely, we found very low membrane-based receptor signaling, exhibiting <1% of the total receptor signaling for these eight RTKs. Moreover, we found that nuclear translocation, mechanistically, requires late endosomal transport; when we blocked receptor trafficking from late endosomes to the nucleus we found a 57% reduction in nuclear translocation. In summary, our research has elucidated the significance of endocytic vesicles, late endosomes and the nucleus in RTK signal propagation.

Cyclic trans-phosphorylation in a homodimer as the predominant mechanism of EGFRvIII action and regulation

Despite intensive research no therapies targeted against the oncogenic EGFRvIII are present in the clinic. One of the reasons is the elusive nature of the molecular structure and activity of the truncated receptor. The recent publications indicate the EGF-bound wild-type EGFR to trans-phosphorylate the EGFRvIII initiating aberrant signaling cascade. The elevated stability of the mutant receptor contributes towards oncogenic potential, preventing termination of signaling by receptor degradation. Here, we show that inhibition of phosphatases leads to a marked increase in phosphorylation of wild-type EGFR and EGFRvIII, indicating that both undergo cyclic rounds of phosphorylation and dephosphorylation on all investigated tyrosine residues, including Tyr1045. Still, we observe elevated stability of the mutant receptor, suggesting phosphorylation as insufficient to cause degradation. Hyperphosphorylation of EGFRvIII was hindered only by EGFR tyrosine kinase inhibitors. Co-immunoprecipitation as well as semi-native Western blotting structural analyses together with functional investigation of EGFRvIII's phosphorylation following depletion of wild-type EGFR by shRNA or EGF-mediated degradation indicated homodimerization as the predominant quaternary structure of the mutant receptor. Dimers were observed only under non-reducing conditions, suggesting that homodimerization is mediated by covalent bonds. Previous reports indicated cysteine at position 16 to mediate covalent homodimerization. Upon its substitution to serine, we have observed impaired formation of dimers and lower phosphorylation levels of the mutated oncogene. Based on the obtained results we propose that EGFRvIII is predominantly regulated dynamically by phosphatases that counteract the process of trans-phosphorylation occurring within the homodimers.

An engineering design approach to systems biology

Measuring and modeling the integrated behavior of biomolecular-cellular networks is central to systems biology. Over several decades, systems biology has been shaped by quantitative biologists, physicists, mathematicians, and engineers in different ways. However, the basic and applied versions of systems biology are not typically distinguished, which blurs the separate aspirations of the field and its potential for real-world impact. Here, we articulate an engineering approach to systems biology, which applies educational philosophy, engineering design, and predictive models to solve contemporary problems in an age of biomedical Big Data. A concerted effort to train systems bioengineers will provide a versatile workforce capable of tackling the diverse challenges faced by the biotechnological and pharmaceutical sectors in a modern, information-dense economy.

PTPRO represses ERBB2-driven breast oncogenesis by dephosphorylation and endosomal internalization of ERBB2

The plasma membrane-associated tyrosine phosphatase PTPRO is frequently transcriptionally repressed in cancers and signifies poor prognosis of breast cancer patients. In this study, deletion of Ptpro in MMTV-Erbb2 transgenic mice dramatically shortened the mammary tumor latency and accelerated tumor growth due to loss of Ptpro within the breast cancer cells but not in surrounding tissue as confirmed by hetero-transplantation studies. Both in vitro and in vivo data demonstrated that the phosphatase activity was required for the inactivation of ERBB2 and its downstream signaling. PTPRO regulated the phosphorylation status of ERBB2 at Y1248. Co-immunoprecipitation and proximity ligation assay (Duolink) indicated that PTPRO directly physically interacted with ERBB2. Moreover, PTPRO phosphatase activity shortened the half-life of ERBB2 by increasing endocytotic degradation. PTPRO reexpression by demethylation treatment using 5-azacytidine reduced the proliferation and colony formation potential in ERBB2-positive breast cancer cells. Taken together, PTPRO inhibited ERBB2-driven breast cancer through dephosphorylation leading to dual effects of ERBB2 signaling suppression and endosomal internalization of ERBB2, Therefore, reexpression of PTPRO may be a potential therapy for ERBB2-overexpressing breast cancer.

The Dark Side of Cell Signaling: Positive Roles for Negative Regulators

Cell signaling is dominated by analyzing positive responses to stimuli. Signal activation is balanced by negative regulators that are generally considered to terminate signaling. Rather than exerting only negative effects, however, many such regulators play important roles in enhancing cell-signaling control. Considering responses downstream of selected cell-surface receptors, we discuss how receptor internalization affects signaling specificity and how rapid kinase/phosphatase and GTP/GDP cycles increase responsiveness and allow kinetic proofreading in receptor signaling. We highlight the blurring of distinctions between positive and negative signals, recasting signal termination as the response to a switch-like transition into a new cellular state.

Cell signaling regulation by protein phosphorylation: a multivariate, heterogeneous, and context-dependent process

Proper spatiotemporal regulation of protein phosphorylation in cells and tissues is required for normal development and homeostasis, but aberrant protein phosphorylation regulation leads to various diseases. The study of signaling regulation by protein phosphorylation is complicated in part by the sheer scope of the kinome and phosphoproteome, dependence of signaling protein functionality on cellular localization, and the complex multivariate relationships that exist between protein phosphorylation dynamics and the cellular phenotypes they control. Additional complexities arise from the ability of microenvironmental factors to influence phosphorylation-dependent signaling and from the tendency for some signaling processes to occur heterogeneously among cells. These considerations should be taken into account when measuring cell signaling regulation by protein phosphorylation.

Single Particle Tracking Reveals that EGFR Signaling Activity Is Amplified in Clathrin-Coated Pits

Signaling from the epidermal growth factor receptor (EGFR) via phosphorylation on its C-terminal tyrosine residues requires self-association, which depends on the diffusional properties of the receptor and its density in the plasma membrane. Dimerization is a key event for EGFR activation, but the role of higher order clustering is unknown. We employed single particle tracking to relate the mobility and aggregation of EGFR to its signaling activity. EGFR mobility alternates between short-lived free, confined and immobile states. In the immobile state, EGFR tends to aggregate in clathrin-coated pits, which is further enhanced in a phosphorylation-dependent manner and does not require ligand binding. EGFR phosphorylation is further amplified by cross-phosphorylation in clathrin-coated pits. Because phosphorylated receptors can escape from the pits, local gradients of signaling active EGFR are formed. These results show that amplification of EGFR phosphorylation by receptor clustering in clathrin-coated pits supports signal activation at the plasma membrane.

Quantitative analysis reveals how EGFR activation and downregulation are coupled in normal but not in cancer cells

Ubiquitination of the epidermal growth factor receptor (EGFR) that occurs when Cbl and Grb2 bind to three phosphotyrosine residues (pY1045, pY1068 and pY1086) on the receptor displays a sharp threshold effect as a function of EGF concentration. Here we use a simple modelling approach together with experiments to show that the establishment of the threshold requires both the multiplicity of binding sites and cooperative binding of Cbl and Grb2 to the EGFR. While the threshold is remarkably robust, a more sophisticated model predicted that it could be modulated as a function of EGFR levels on the cell surface. We confirmed experimentally that the system has evolved to perform optimally at physiological levels of EGFR. As a consequence, this system displays an intrinsic weakness that causes--at the supraphysiological levels of receptor and/or ligand associated with cancer--uncoupling of the mechanisms leading to signalling through phosphorylation and attenuation through ubiquitination.

The AXL Receptor is a Sensor of Ligand Spatial Heterogeneity

The AXL receptor is a TAM (Tyro3, AXL, MerTK) receptor tyrosine kinase (RTK) important in physiological inflammatory processes such as blood clotting, viral infection, and innate immune-mediated cell clearance. Overexpression of the receptor in a number of solid tumors is increasingly appreciated as a key drug resistance and tumor dissemination mechanism. Although the ligand-receptor (Gas6-AXL) complex structure is known, literature reports on ligand-mediated signaling have provided conflicting conclusions regarding the influence of other factors such as phosphatidylserine binding, and a detailed, mechanistic picture of AXL activation has not emerged. Integrating quantitative experiments with mathematical modeling, we show here that AXL operates to sense local spatial heterogeneity in ligand concentration, a feature consistent with its physiological role in inflammatory cell responses. This effect arises as a result of an intricate reaction-diffusion interaction. Our results demonstrate that AXL functions distinctly from other RTK families, a vital insight for envisioned design of AXL-targeted therapeutic intervention.

Identifying Determinants of EGFR-Targeted Therapeutic Biochemical Efficacy Using Computational Modeling

We modeled cellular epidermal growth factor receptor (EGFR) tyrosine phosphorylation dynamics in the presence of receptor-targeting kinase inhibitors (e.g., gefitinib) or antibodies (e.g., cetuximab) to identify systematically the factors that contribute most to the ability of the therapeutics to antagonize EGFR phosphorylation, an effect we define here as biochemical efficacy. Our model identifies distinct processes as controlling gefitinib or cetuximab biochemical efficacy, suggests biochemical efficacy is favored in the presence of certain EGFR ligands, and suggests new drug design principles. For example, the model predicts that gefitinib biochemical efficacy is preferentially sensitive to perturbations in the activity of tyrosine phosphatases regulating EGFR, but that cetuximab biochemical efficacy is preferentially sensitive to perturbations in ligand binding. Our results highlight numerous other considerations that determine biochemical efficacy beyond those reflected by equilibrium affinities. By integrating these considerations, our model also predicts minimum therapeutic combination concentrations to maximally reduce receptor phosphorylation.

Differential parsing of EGFR endocytic flux among parallel internalization pathways in lung cancer cells with EGFR-activating mutations

Due to the existence of parallel pathways for receptor endocytosis and their complexities, a quantitative understanding of receptor endocytosis in normal and pathological settings requires computational analysis. Here, we develop a mechanistic model of epidermal growth factor receptor (EGFR) endocytosis to determine the relative contributions of three parallel pathways: clathrin-dependent internalization mediated by mitogen-inducible gene 6 (MIG6), an endogenous EGFR kinase inhibitor that links EGFR to endocytic proteins; clathrin-dependent internalization mediated by the ubiquitin ligase CBL, which can be sequestered by the regulatory protein Sprouty2; or alternative pathways that may be non-clathrin mediated. We applied the model to interpret our previous measurements of EGFR endocytosis in lung cancer cells. Interestingly, our results suggest that MIG6 is responsible for at least as much wild-type EGFR internalization as CBL, indicating that a significant fraction of internalizing EGFR may be incapable of driving signaling. Model results also suggest that MIG6's endocytic function is reduced for the kinase-activated and internalization-impaired EGFR mutants found in some lung cancers. Analysis of Sprouty2 knockdown data indicates that Sprouty2 regulates EGFR endocytosis primarily by controlling EGFR expression, rather than by sequestering CBL, and supports the notion that CBL-mediated internalization is impaired for EGFR mutants. We further demonstrate that differences in internalization between wild-type and mutant EGFR cannot explain differences in EGF-mediated EGFR degradation without concomitant changes in EGFR recycling, which we previously quantified. This work provides new quantitative insights into EGFR trafficking in lung cancer and provides a framework for studying parallel endocytosis pathways for other receptors.