Topological analysis of MAPK cascade for kinetic ErbB signaling

Knowledge-based mechanistic modeling accurately predicts disease progression with gefitinib in EGFR-mutant lung adenocarcinoma

Lung adenocarcinoma (LUAD) is associated with a low survival rate at advanced stages. Although the development of targeted therapies has improved outcomes in LUAD patients with identified and specific genetic alterations, such as activating mutations on the epidermal growth factor receptor gene (EGFR), the emergence of tumor resistance eventually occurs in all patients and this is driving the development of new therapies. In this paper, we present the In Silico EGFR-mutant LUAD (ISELA) model that links LUAD patients' individual characteristics, including tumor genetic heterogeneity, to tumor size evolution and tumor progression over time under first generation EGFR tyrosine kinase inhibitor gefitinib. This translational mechanistic model gathers extensive knowledge on LUAD and was calibrated on multiple scales, including in vitro, human tumor xenograft mouse and human, reproducing more than 90% of the experimental data identified. Moreover, with 98.5% coverage and 99.4% negative logrank tests, the model accurately reproduced the time to progression from the Lux-Lung 7 clinical trial, which was unused in calibration, thus supporting the model high predictive value. This knowledge-based mechanistic model could be a valuable tool in the development of new therapies targeting EGFR-mutant LUAD as a foundation for the generation of synthetic control arms.

p52Shc regulates the sustainability of ERK activation in a RAF-independent manner

p52SHC (SHC) and GRB2 are adaptor proteins involved in the RAS/MAPK (ERK) pathway mediating signals from cell-surface receptors to various cytoplasmic proteins. To further examine their roles in signal transduction, we studied the translocation of fluorescently labeled SHC and GRB2 to the cell surface, caused by the activation of ERBB receptors by heregulin (HRG). We simultaneously evaluated activated ERK translocation to the nucleus. Unexpectedly, the translocation dynamics of SHC were sustained when those of GRB2 were transient. The sustained localization of SHC positively correlated with the sustained nuclear localization of ERK, which became more transient after SHC knockdown. SHC-mediated PI3K activation was required to maintain the sustainability of the ERK translocation regulating MEK but not RAF. In cells overexpressing ERBB1, SHC translocation became transient, and the HRG-induced cell fate shifted from a differentiation to a proliferation bias. Our results indicate that SHC and GRB2 functions are not redundant but that SHC plays the critical role in the temporal regulation of ERK activation.

ErbB Family Signalling: A Paradigm for Oncogene Addiction and Personalized Oncology

ErbB family members represent important biomarkers and drug targets for modern precision therapy. They have gained considerable importance as paradigms for oncoprotein addiction and personalized medicine. This review summarizes the current understanding of ErbB proteins in cell signalling and cancer and describes the molecular rationale of prominent cases of ErbB oncoprotein addiction in different cancer types. In addition, we have highlighted experimental technologies for the development of innovative cancer cell models that accurately predicted clinical ErbB drug efficacies. In the future, such cancer models might facilitate the identification and validation of physiologically relevant novel forms of oncoprotein and non-oncoprotein addiction or synthetic lethality. The identification of genotype-drug response relationships will further advance personalized oncology and improve drug efficacy in the clinic. Finally, we review the most important drugs targeting ErbB family members that are under investigation in clinical trials or that made their way already into clinical routine. Taken together, the functional characterization of ErbB oncoproteins have significantly increased our knowledge on predictive biomarkers, oncoprotein addiction and patient stratification and treatment.

Helicobacter pylori Induces Serine Phosphorylation of EGFR via Novel TAK1-p38 Activation Pathway in an HB-EGF-Independent Manner

BACKGROUND

The interaction of Helicobacter pylori with gastric epithelial cells can result in the activation of transcription factor NF-κB via TGF-β-activated kinase 1 (TAK1). In this study, we have demonstrated the role of H. pylori in the activation of EGFR via TAK1-mediated phosphorylation of p38.

MATERIALS AND METHODS

Gastric epithelial AGS or MKN-45 cells were co-cultured with wild-type or cagA(-) H. pylori strains. H. pylori was added to the cells, and the activation of EGFR, p65 (NF-κB) subunit, p38, ERK, and TAK1 was examined by Western blotting. Infected cells were pretreated with or without ligands, chemical inhibitors, anti-HB-EGF antibody, and siRNAs to evaluate the effects on phosphorylation of various EGFR residues. Fluorescence microscopy and flow cytometry were performed to detect the internalization of EGFR.

RESULTS

Incubating cells with wild-type and CagA(-) H. pylori strains resulted in the rapid and transient phosphorylation of serine residues of EGFR. RNAi experiments using siRNA against TAK1 and p38 pathways blocked the phosphorylation of serine residue. Immunofluorescence and flow cytometry revealed that EGFR was internalized in H. pylori-infected cells after EGFR phosphorylation in a p38-dependent manner. In contrast, pretreatment with gefitinib and anti-HB-EGF antibody did not block both the phosphorylation and internalization of EGFR.

CONCLUSION

Helicobacter pylori induces internalization of EGFR via novel TAK1-p38-serine activation pathway which is independent of HB-EGF. The interaction between TAK1 and EGFR in H. pylori-infected cells might open new dimensions in understanding H. pylori-associated gastric carcinogenesis.

Inverse correlation between Thr-669 and constitutive tyrosine phosphorylation in the asymmetric epidermal growth factor receptor dimer conformation

We have recently identified tumor necrosis factor (TNF)-α-induced phosphorylation of epidermal growth factor receptor (EGFR) at Thr-669 and Ser-1046/1047 via ERK and p38 pathways, respectively. In the present study, we investigated the roles of ligand-induced phosphorylation of serine and threonine residues in EGFR-overexpressing MDA-MB-468 breast cancer cells. Epidermal growth factor and heregulin, an ErbB3 ligand, induced the phosphorylation of Thr-669 and Ser-1046/1047. Inversely, constitutive tyrosine phosphorylation of the C-terminal domain, including Tyr-1068, was significantly downregulated on ligand stimulation. Inhibition of the ERK pathway by U0126 blocked ligand-induced Thr-669 phosphorylation as well as Tyr-1068 dephosphorylation. Downregulation of constitutive tyrosine phosphorylation of EGFR in HEK293 cells stably expressing the wild type was abolished by substitution of Thr-669 for Ala. In an asymmetric EGFR homodimer structure, one Thr-669 in the receiver kinase of the dimer was involved in downregulation. Similarly, Thr-669 in an EGFR-ErbB3 heterodimer also participated in tyrosine dephosphorylation. These results indicate that ERK-mediated Thr-669 phosphorylation suppresses constitutive tyrosine phosphosphorylation in the homo- and heterodimer asymmetric conformations of the EGFR.

Random parameter sampling of a generic three-tier MAPK cascade model reveals major factors affecting its versatile dynamics

The mitogen-activated protein kinase (MAPK) pathway is considered to be a central block in many biological signaling networks. Despite the common core cascade structure, the activation of MAPK in different biological systems can exhibit different types of dynamic behaviors. Computer modeling may help to reveal the mechanisms underlying such variations. However, so far most computational models of the MAPK cascade have been system-specific, or to reflect a particular type among the wide spectrum of possible dynamics. To obtain a general and integrated view of the relationship between the dynamics of MAPK activation and the structures and parameters of the MAPK cascade, we constructed a generic model by comparing previous models covering different specific biological systems. We assumed that reliable qualitative results could be predicted through a qualitative model with pseudo parameters. We used randomly sampled parameters instead of a specific set of “best-fit” parameters to avoid biases towards any particular systems. A range of dynamics behaviors for MAPK activation, including ultrasensitivity, bistability, transient activation and oscillation, were successfully predicted by the generic model. The results indicated that the steady state dynamics (ultrasensitivity and bistability) was jointly determined by the three-tiered structure of the MAPK cascade and the competitive substrate binding in the dual-phosphorylation processes of the central components, while the temporal dynamics (transient activation and oscillation) was mainly affected by the upstream signaling pathway and feedbacks. Moreover, MAPK kinase (MAPKK) played more important roles than the other two components in determining the dynamics of MAPK activation. We hypothesize that this is an important and advantageous property for the regulation and for the functional diversity of MAPK pathways in real cells. Finally, to assist developing generic models for signaling motifs through model comparisons, we proposed a reaction-based database to make the model data more flexible and interoperable.

A systematic modeling study on the pathogenic role of p38 MAPK activation in myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are a group of clonal hematopoietic stem cell diseases. In addition to intrinsic genetic alterations, the effects of the extrinsic microenvironment also play a pathological role in MDS development. The presence of increased inflammatory cytokines, such as tumor necrosis factor-alpha (TNF-α), in marrow and abnormal activation of the p38 mitogen-activated protein kinase (MAPK) signaling pathway in hematopoietic cells are associated with the ineffective hematopoiesis in MDS. However, the molecular mechanism of p38 MAPK activation triggered by microenvironment cytokines remains poorly understood. To address this question, we combined computational modeling analysis and molecular biology studies to perform a systematic investigation of signaling events regulated by microenvironment cytokines in hematopoietic cells from MDS patients. We examined dynamic changes of key signaling events, including the p38 MAPK and the c-Jun N-terminal kinase (JNK) pathway in bone marrow mononuclear cells from MDS patients or normal donors in response to TNF-α stimulation using reverse phase protein array technology. The results were analyzed by a novel computational model and preliminarily validated by immunohistochemistry analysis of the bone marrow tissues from twelve MDS patients and normal donors. Our systematic model revealed that the dynamic response patterns of p38 MAPK and JNK to TNF-α stimulation in MDS were different from that observed in normal marrow cells. Particularly, B-cell lymphoma-X (BCL-XL) protein degradation was regulated by the JNK pathway in normal cells, but by p38 MAPK in MDS cells. By immunohistochemistry, BCL-XL was highly expressed in hematopoietic cells from normal marrow, but was minimally expressed in MDS marrow. Additionally, immunostaining for phosphorylated p38 MAPKα showed much higher p38 MAPK activation in MDS marrows, supporting over-activation of p38 MAPK-enhanced degradation of BCL-XL in MDS. The degradation of BCL-XL triggered by p38 MAPK over-activation may contribute to the increasing apoptosis of marrow cells, a phenomenon commonly observed in MDS, and lead to ineffective hematopoiesis. Our study suggests that the combination of molecular biological studies and systematic modeling is a powerful tool for comprehensive investigation of the complex cellular mechanisms involved in MDS pathogenesis.

Model-based global sensitivity analysis as applied to identification of anti-cancer drug targets and biomarkers of drug resistance in the ErbB2/3 network

High levels of variability in cancer-related cellular signalling networks and a lack of parameter identifiability in large-scale network models hamper translation of the results of modelling studies into the process of anti-cancer drug development. Recently global sensitivity analysis (GSA) has been recognised as a useful technique, capable of addressing the uncertainty of the model parameters and generating valid predictions on parametric sensitivities. Here we propose a novel implementation of model-based GSA specially designed to explore how multi-parametric network perturbations affect signal propagation through cancer-related networks. We use area-under-the-curve for time course of changes in phosphorylation of proteins as a characteristic for sensitivity analysis and rank network parameters with regard to their impact on the level of key cancer-related outputs, separating strong inhibitory from stimulatory effects. This allows interpretation of the results in terms which can incorporate the effects of potential anti-cancer drugs on targets and the associated biological markers of cancer. To illustrate the method we applied it to an ErbB signalling network model and explored the sensitivity profile of its key model readout, phosphorylated Akt, in the absence and presence of the ErbB2 inhibitor pertuzumab. The method successfully identified the parameters associated with elevation or suppression of Akt phosphorylation in the ErbB2/3 network. From analysis and comparison of the sensitivity profiles of pAkt in the absence and presence of targeted drugs we derived predictions of drug targets, cancer-related biomarkers and generated hypotheses for combinatorial therapy. Several key predictions have been confirmed in experiments using human ovarian carcinoma cell lines. We also compared GSA-derived predictions with the results of local sensitivity analysis and discuss the applicability of both methods. We propose that the developed GSA procedure can serve as a refining tool in combinatorial anti-cancer drug discovery.

Epidermal growth factor receptor mutation in combination with expression of MIG6 alters gefitinib sensitivity

BACKGROUND

Epidermal growth factor receptor (EGFR) signaling plays an important role in the regulation of cell proliferation, survival, metastasis, and invasion in various tumors. Earlier studies showed that the EGFR is frequently overexpressed in non-small-cell lung cancer (NSCLC) and EGFR mutations at specific amino acid residues in the kinase domain induce altered responsiveness to gefitinib, a small molecule EGFR tyrosine kinase inhibitor. However, the mechanism underlying the drug response modulated by EGFR mutation is still largely unknown. To elucidate drug response in EGFR signal transduction pathway in which complex dynamics of multiple molecules involved, a systematic approach is necessary. In this paper, we performed experimental and computational analyses to clarify the underlying mechanism of EGFR signaling and cell-specific gefitinib responsiveness in three H1299-derived NSCLC cell lines; H1299 wild type (H1299WT), H1299 with an overexpressed wild type EGFR (H1299EGFR-WT), and H1299 with an overexpressed mutant EGFR L858R (H1299L858R; gefitinib sensitive mutant).

RESULTS

We predicted and experimentally verified that Mig6, which is a known negative regulator of EGFR and specifically expressed in H1299L858R cells, synergized with gefitinib to suppress cellular growth. Computational analyses indicated that this inhibitory effect is amplified at the phosphorylation/dephosphorylation steps of MEK and ERK.

CONCLUSIONS

Thus, we showed that L858R receptor mutation in combination with expression of its negative regulator, Mig6, alters signaling outcomes and results in variable drug sensitivity.

TAK1-mediated serine/threonine phosphorylation of epidermal growth factor receptor via p38/extracellular signal-regulated kinase: NF-{kappa}B-independent survival pathways in tumor necrosis factor alpha signaling

The kinase TAK1, a mitogen-activated protein kinase kinase kinase (MAP3K), has been widely accepted as a key kinase activating NF-kappaB and MAPKs in tumor necrosis factor alpha (TNF-alpha) signaling. We have recently reported that TAK1 regulates the transient phosphorylation and endocytosis of epidermal growth factor receptor (EGFR) in a tyrosine kinase activity-independent manner. In the present study, we found that Thr-669 in the juxtamembrane domain and Ser-1046/1047 in the carboxyl-terminal regulatory domain were transiently phosphorylated in response to TNF-alpha. Experiments using chemical inhibitors and small interfering RNA demonstrated that TNF-alpha-mediated phosphorylation of Thr-669 and Ser-1046/7 were differently regulated via TAK1-extracellular signal-regulated kinase (ERK) and TAK1-p38 pathways, respectively. In addition, p38, but not ERK, was involved in the endocytosis of EGFR. Surprisingly, modified EGFR was essential to prevent apoptotic cellular responses; however, the EGFR pathway was independent of the NF-kappaB antiapoptotic pathway. These results demonstrated that TAK1 controls two different signaling pathways, IkappaB kinase-NF-kappaB and MAPK-EGFR, leading to the survival of cells exposed to the death signal from the TNF-alpha receptor.

Positive and negative symptoms in schizophrenia: the NMDA receptor hypofunction hypothesis, neuregulin/ErbB4 and synapse regression

Carlsson has put forward the hypothesis that the positive and negative symptoms of schizophrenia are due to failure of mesolimbic and mesocortical projections consequent on hypofunction of the glutamate N-methyl-d-aspartate (NMDA) receptor. The hypothesis has been recently emphasized in this Journal that the loss of synaptic spines with NMDA receptors, which can be precipitated by stress, can explain the emergence of positive symptoms such as hallucinations and that this synapse regression involves molecules such as neuregulin and its receptor ErbB4 that have been implicated in schizophrenia. In this essay these two hypotheses are brought together in a single scheme in which emphasis is placed on the molecular pathways from neuregulin/ErbB4, to modulation of the NMDA receptors, subsequent changes in the synaptic spine's cytoskeletal apparatus and so regression of the spines. It is suggested that identification of the molecular constituents of this pathway will allow synthesis of suitable substances for removing the hypofunction of NMDA receptors and so the phenotypic consequences that flow from this hypofunction.

Systems pharmacology and genome medicine: a future perspective

Genome medicine uses genomic information in the diagnosis of disease and in prescribing treatment. This transdisciplinary field brings together knowledge on the relationships between genetics, pathophysiology and pharmacology. Systems pharmacology aims to understand the actions and adverse effects of drugs by considering targets in the context of the biological networks in which they exist. Genome medicine forms the base on which systems pharmacology can develop. Experimental and computational approaches enable systems pharmacology to obtain holistic, mechanistic information on disease networks and drug responses, and to identify new drug targets and specific drug combinations. Network analyses of interactions involved in pathophysiology and drug response across various scales of organization, from molecular to organismal, will allow the integration of the systems-level understanding of drug action with genome medicine. The interface of the two fields will enable drug discovery for personalized medicine. Here we provide a perspective on the questions and approaches that drive the development of these new interrelated fields.