51
|
Erickson KE, Rukhlenko OS, Shahinuzzaman M, Slavkova KP, Lin YT, Suderman R, Stites EC, Anghel M, Posner RG, Barua D, Kholodenko BN, Hlavacek WS. Modeling cell line-specific recruitment of signaling proteins to the insulin-like growth factor 1 receptor. PLoS Comput Biol 2019; 15:e1006706. [PMID: 30653502 PMCID: PMC6353226 DOI: 10.1371/journal.pcbi.1006706] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Revised: 01/30/2019] [Accepted: 12/09/2018] [Indexed: 12/27/2022] Open
Abstract
Receptor tyrosine kinases (RTKs) typically contain multiple autophosphorylation sites in their cytoplasmic domains. Once activated, these autophosphorylation sites can recruit downstream signaling proteins containing Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains, which recognize phosphotyrosine-containing short linear motifs (SLiMs). These domains and SLiMs have polyspecific or promiscuous binding activities. Thus, multiple signaling proteins may compete for binding to a common SLiM and vice versa. To investigate the effects of competition on RTK signaling, we used a rule-based modeling approach to develop and analyze models for ligand-induced recruitment of SH2/PTB domain-containing proteins to autophosphorylation sites in the insulin-like growth factor 1 (IGF1) receptor (IGF1R). Models were parameterized using published datasets reporting protein copy numbers and site-specific binding affinities. Simulations were facilitated by a novel application of model restructuration, to reduce redundancy in rule-derived equations. We compare predictions obtained via numerical simulation of the model to those obtained through simple prediction methods, such as through an analytical approximation, or ranking by copy number and/or KD value, and find that the simple methods are unable to recapitulate the predictions of numerical simulations. We created 45 cell line-specific models that demonstrate how early events in IGF1R signaling depend on the protein abundance profile of a cell. Simulations, facilitated by model restructuration, identified pairs of IGF1R binding partners that are recruited in anti-correlated and correlated fashions, despite no inclusion of cooperativity in our models. This work shows that the outcome of competition depends on the physicochemical parameters that characterize pairwise interactions, as well as network properties, including network connectivity and the relative abundances of competitors. Cells rely on networks of interacting biomolecules to sense and respond to environmental perturbations and signals. However, it is unclear how information is processed to generate appropriate and specific responses to signals, especially given that these networks tend to share many components. For example, receptors that detect distinct ligands and regulate distinct cellular activities commonly interact with overlapping sets of downstream signaling proteins. Here, to investigate the downstream signaling of a well-studied receptor tyrosine kinase (RTK), the insulin-like growth factor 1 (IGF1) receptor (IGF1R), we formulated and analyzed 45 cell line-specific mathematical models, which account for recruitment of 18 different binding partners to six sites of receptor autophosphorylation in IGF1R. The models were parameterized using available protein copy number and site-specific affinity measurements, and restructured to allow for network generation. We find that recruitment is influenced by the protein abundance profile of a cell, with different patterns of recruitment in different cell lines. Furthermore, in a given cell line, we find that pairs of IGF1R binding partners may be recruited in a correlated or anti-correlated fashion. We demonstrate that the simulations of the model have greater predictive power than protein copy number and/or binding affinity data, and that even a simple analytical model cannot reproduce the predicted recruitment ranking obtained via simulations. These findings represent testable predictions and indicate that the outputs of IGF1R signaling depend on cell line-specific properties in addition to the properties that are intrinsic to the biomolecules involved.
Collapse
Affiliation(s)
- Keesha E. Erickson
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | | | - Md Shahinuzzaman
- Department of Chemical and Biochemical Engineering, University of Missouri Science and Technology, Rolla, Missouri, United States of America
| | - Kalina P. Slavkova
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Yen Ting Lin
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Ryan Suderman
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Edward C. Stites
- The Salk Institute for Biological Studies, La Jolla, California, United States of America
| | - Marian Anghel
- Information Sciences Group, Computer, Computational and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Richard G. Posner
- Department of Biological Sciences, Northern Arizona University, Flagstaff, Arizona, United States of America
| | - Dipak Barua
- Department of Chemical and Biochemical Engineering, University of Missouri Science and Technology, Rolla, Missouri, United States of America
| | - Boris N. Kholodenko
- Systems Biology Ireland, University College Dublin, Belfield, Dublin, Ireland
- School of Medicine and Medical Science and Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin, Ireland
| | - William S. Hlavacek
- Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
- * E-mail:
| |
Collapse
|
52
|
|
53
|
Molecular pathway activation – New type of biomarkers for tumor morphology and personalized selection of target drugs. Semin Cancer Biol 2018; 53:110-124. [DOI: 10.1016/j.semcancer.2018.06.003] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2018] [Revised: 06/19/2018] [Accepted: 06/19/2018] [Indexed: 02/06/2023]
|
54
|
Zhou W, Chen Z, Liu Z, Wang Y. Stochasticity and robustness analysis of microRNA-mediated ERK signaling network. Comput Biol Chem 2018; 76:318-326. [PMID: 30144727 DOI: 10.1016/j.compbiolchem.2018.07.018] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Revised: 07/27/2018] [Accepted: 07/30/2018] [Indexed: 11/30/2022]
Abstract
MicroRNAs (miRNAs) play a critical role in regulating the signaling network such as the extracellular signal regulated kinase (ERK) pathway. However, the mechanisms of miRNA-mediated posttranscriptional regulations with regard to their impacts on signaling require further investigation from a systematic level. Therefore, we develop a mathematical model by analyzing the precise dynamic mechanisms with a bevy of miRNAs involved in the Ras/Raf/MEK/ERK pathway. A systems-based analysis approach is introduced into this model and the dynamics have been implemented deterministically and stochastically. Our analysis reveals that miRNAs are key participants regulating the gene expression of ERK network, and the cooperative actions of miRNAs are important in keeping the normal biological characteristics and amplification effects of the system. Meanwhile, the appearance of system disorder in the absence of miRNAs suggests that miRNAs may play a role in the pathological processes, such as tumor and inflammation. The sensitivity analysis and the kinetic parameter perturbation show that the binding of receptor (EGFR) and adaptor protein (Shc, Grb2 and Sos) under normal physiological conditions is crucial for the robustness of the whole pathway. In addition, the stochastic dynamic patterns are in a good agreement with the deterministic results, further demonstrating that the variability of the system due to the presence of some stochastic noise is low. All these will be helpful for a deeper understanding of the dynamic mechanism of miRNA-mediated ERK signal network, which might present a rich area of future research with the relevant regulatory roles of miRNAs in cell signaling pathway.
Collapse
Affiliation(s)
- Wei Zhou
- Department of Allergy, The Third Affiliated Hospital of Shenzhen University, Shenzhen, 518020, China.
| | - Ziyi Chen
- Musculoskeletal Research Laboratory, Department of Orthopaedics & Traumatology, Faculty of Medicine, The Chinese University of Hong Kong, Shatin, 999077, Hong Kong Special Administrative Region
| | - Zhigang Liu
- Department of Allergy, The Third Affiliated Hospital of Shenzhen University, Shenzhen, 518020, China
| | - Yonghua Wang
- College of Life Sciences, Northwest University, Shaanxi, 712100, China.
| |
Collapse
|
55
|
Dissecting RAF Inhibitor Resistance by Structure-based Modeling Reveals Ways to Overcome Oncogenic RAS Signaling. Cell Syst 2018; 7:161-179.e14. [PMID: 30007540 DOI: 10.1016/j.cels.2018.06.002] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2017] [Revised: 03/09/2018] [Accepted: 06/04/2018] [Indexed: 12/19/2022]
Abstract
Clinically used RAF inhibitors are ineffective in RAS mutant tumors because they enhance homo- and heterodimerization of RAF kinases, leading to paradoxical activation of ERK signaling. Overcoming enhanced RAF dimerization and the resulting resistance is a challenge for drug design. Combining multiple inhibitors could be more effective, but it is unclear how the best combinations can be chosen. We built a next-generation mechanistic dynamic model to analyze combinations of structurally different RAF inhibitors, which can efficiently suppress MEK/ERK signaling. This rule-based model of the RAS/ERK pathway integrates thermodynamics and kinetics of drug-protein interactions, structural elements, posttranslational modifications, and cell mutational status as model rules to predict RAF inhibitor combinations for inhibiting ERK activity in oncogenic RAS and/or BRAFV600E backgrounds. Predicted synergistic inhibition of ERK signaling was corroborated by experiments in mutant NRAS, HRAS, and BRAFV600E cells, and inhibition of oncogenic RAS signaling was associated with reduced cell proliferation and colony formation.
Collapse
|
56
|
In silico prediction of ErbB signal activation from receptor expression profiles through a data analytics pipeline. J Biosci 2018. [DOI: 10.1007/s12038-018-9747-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
57
|
Huang X, Post JN, Zhong L, Leijten J, Larsson S, Karperien M, Struglics A. Dickkopf-related protein 1 and gremlin 1 show different response than frizzled-related protein in human synovial fluid following knee injury and in patients with osteoarthritis. Osteoarthritis Cartilage 2018. [PMID: 29526783 DOI: 10.1016/j.joca.2018.02.904] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
OBJECTIVE To explore the involvement of the wingless-type MMTV integration site (WNT) and bone morphogenetic protein (BMP) antagonists dickkopf-related protein 1 (DKK1), frizzled-related protein (FRZB) and gremlin 1 (GREM1) in knee injury and osteoarthritis (OA). DESIGN The antagonists were immunoassayed in synovial fluid from a cross-sectional cohort of nine knee healthy reference subjects, patients with recent (0-77 days, n = 158) or old (1-37 years, n = 50) knee injuries, and OA (n = 22). Cartilage (ARGS-aggrecan, cartilage oligomeric matrix protein and C2C type II collagen) and other biomarkers were assessed in synovial fluid in a subset of samples. Statistical analysis was by Kendall's tau (τ) correlation, Mann-Whitney U test, and linear regression analysis. RESULTS Compared to references, median concentration of GREM1 (but not DKK1 and FRZB) was elevated 1.5-fold immediately after injury, and FRZB was reduced 1000-folds in OA. All three antagonists decreased with increasing time after injury as well as with increasing age, but the temporal change after injury was less accentuated for FRZB (peaked 8-22 days after injury) compared to that of DKK1 and GREM1 (peaked immediately after injury). In the recent injury group, there was a correlation between GREM1 and DKK1 (τ = 0.172); FRZB concentrations correlated with concentrations of cartilage biomarkers (τ between 0.257 and 0.369), while DKK1 and GREM1 were inversely correlated (τ between -0.177 and -0.217) with these markers. CONCLUSIONS Our results indicate separate roles for the antagonists, where DKK1 and GREM1 had similarities in response to injury and in OA, with a different response for FRZB.
Collapse
Affiliation(s)
- X Huang
- Department of Developmental BioEngineering, MIRA Institute for Biomedical Technology and Technical Medicine, University of Twente, Enschede, The Netherlands.
| | - J N Post
- Department of Developmental BioEngineering, MIRA Institute for Biomedical Technology and Technical Medicine, University of Twente, Enschede, The Netherlands.
| | - L Zhong
- Department of Developmental BioEngineering, MIRA Institute for Biomedical Technology and Technical Medicine, University of Twente, Enschede, The Netherlands.
| | - J Leijten
- Department of Developmental BioEngineering, MIRA Institute for Biomedical Technology and Technical Medicine, University of Twente, Enschede, The Netherlands.
| | - S Larsson
- Lund University, Faculty of Medicine, Department of Clinical Sciences Lund, Orthopaedics, Lund, Sweden.
| | - M Karperien
- Department of Developmental BioEngineering, MIRA Institute for Biomedical Technology and Technical Medicine, University of Twente, Enschede, The Netherlands.
| | - A Struglics
- Lund University, Faculty of Medicine, Department of Clinical Sciences Lund, Orthopaedics, Lund, Sweden.
| |
Collapse
|
58
|
Abstract
In quantitative analyses of biological processes, one may use many different scales of models (e.g. spatial or non-spatial, deterministic or stochastic, time-varying or at steady-state) or many different approaches to match models to experimental data (e.g. model fitting or parameter uncertainty/sloppiness quantification with different experiment designs). These different analyses can lead to surprisingly different results, even when applied to the same data and the same model. We use a simplified gene regulation model to illustrate many of these concerns, especially for ODE analyses of deterministic processes, chemical master equation and finite state projection analyses of heterogeneous processes, and stochastic simulations. For each analysis, we employ Matlab and Python software to consider a time-dependent input signal (e.g. a kinase nuclear translocation) and several model hypotheses, along with simulated single-cell data. We illustrate different approaches (e.g. deterministic and stochastic) to identify the mechanisms and parameters of the same model from the same simulated data. For each approach, we explore how uncertainty in parameter space varies with respect to the chosen analysis approach or specific experiment design. We conclude with a discussion of how our simulated results relate to the integration of experimental and computational investigations to explore signal-activated gene expression models in yeast (Neuert et al 2013 Science 339 584-7) and human cells (Senecal et al 2014 Cell Rep. 8 75-83)5.
Collapse
Affiliation(s)
- Lisa Weber
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, CO
| | - William Raymond
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, CO
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO
| | - Brian Munsky
- Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, CO
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO
| |
Collapse
|
59
|
Inferring a nonlinear biochemical network model from a heterogeneous single-cell time course data. Sci Rep 2018; 8:6790. [PMID: 29717206 PMCID: PMC5931614 DOI: 10.1038/s41598-018-25064-w] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2017] [Accepted: 04/09/2018] [Indexed: 12/30/2022] Open
Abstract
Mathematical modeling and analysis of biochemical reaction networks are key routines in computational systems biology and biophysics; however, it remains difficult to choose the most valid model. Here, we propose a computational framework for data-driven and systematic inference of a nonlinear biochemical network model. The framework is based on the expectation-maximization algorithm combined with particle smoother and sparse regularization techniques. In this method, a “redundant” model consisting of an excessive number of nodes and regulatory paths is iteratively updated by eliminating unnecessary paths, resulting in an inference of the most likely model. Using artificial single-cell time-course data showing heterogeneous oscillatory behaviors, we demonstrated that this algorithm successfully inferred the true network without any prior knowledge of network topology or parameter values. Furthermore, we showed that both the regulatory paths among nodes and the optimal number of nodes in the network could be systematically determined. The method presented in this study provides a general framework for inferring a nonlinear biochemical network model from heterogeneous single-cell time-course data.
Collapse
|
60
|
Hastings JF, Skhinas JN, Fey D, Croucher DR, Cox TR. The extracellular matrix as a key regulator of intracellular signalling networks. Br J Pharmacol 2018; 176:82-92. [PMID: 29510460 DOI: 10.1111/bph.14195] [Citation(s) in RCA: 129] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Revised: 02/06/2018] [Accepted: 02/13/2018] [Indexed: 12/11/2022] Open
Abstract
The extracellular matrix (ECM) is a salient feature of all solid tissues within the body. This complex, acellular entity is composed of hundreds of individual molecules whose assembly, architecture and biomechanical properties are critical to controlling the behaviour and phenotype of the different cell types residing within tissues. Cells are the basic unit of life and the core building block of tissues and organs. At their simplest, they follow a set of rules, governed by their genetic code and effected through the complex protein signalling networks that these genes encode. These signalling networks assimilate and process the information received by the cell to control cellular decisions that govern cell fate. The ECM is the biggest provider of external stimuli to cells and as such is responsible for influencing intracellular signalling dynamics. In this review, we discuss the inclusion of ECM as a central regulatory signalling sub-network in computational models of cellular decision making, with a focus on its role in diseases such as cancer. LINKED ARTICLES: This article is part of a themed section on Translating the Matrix. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.1/issuetoc.
Collapse
Affiliation(s)
- Jordan F Hastings
- The Garvan Institute of Medical Research and The Kinghorn Cancer Centre, Cancer Division, Darlinghurst, NSW, 2010, Australia
| | - Joanna N Skhinas
- The Garvan Institute of Medical Research and The Kinghorn Cancer Centre, Cancer Division, Darlinghurst, NSW, 2010, Australia
| | - Dirk Fey
- Systems Biology Ireland, University College Dublin, Dublin 4, Ireland
| | - David R Croucher
- The Garvan Institute of Medical Research and The Kinghorn Cancer Centre, Cancer Division, Darlinghurst, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW Sydney, Kensington, NSW, 2010, Australia.,School of Medicine and Medical Science, University College Dublin, Dublin 4, Ireland
| | - Thomas R Cox
- The Garvan Institute of Medical Research and The Kinghorn Cancer Centre, Cancer Division, Darlinghurst, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW Sydney, Kensington, NSW, 2010, Australia
| |
Collapse
|
61
|
Arkun Y, Yasemi M. Dynamics and control of the ERK signaling pathway: Sensitivity, bistability, and oscillations. PLoS One 2018; 13:e0195513. [PMID: 29630631 PMCID: PMC5891012 DOI: 10.1371/journal.pone.0195513] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2018] [Accepted: 03/22/2018] [Indexed: 02/06/2023] Open
Abstract
Cell signaling is the process by which extracellular information is transmitted into the cell to perform useful biological functions. The ERK (extracellular-signal-regulated kinase) signaling controls several cellular processes such as cell growth, proliferation, differentiation and apoptosis. The ERK signaling pathway considered in this work starts with an extracellular stimulus and ends with activated (double phosphorylated) ERK which gets translocated into the nucleus. We model and analyze this complex pathway by decomposing it into three functional subsystems. The first subsystem spans the initial part of the pathway from the extracellular growth factor to the formation of the SOS complex, ShC-Grb2-SOS. The second subsystem includes the activation of Ras which is mediated by the SOS complex. This is followed by the MAPK subsystem (or the Raf-MEK-ERK pathway) which produces the double phosphorylated ERK upon being activated by Ras. Although separate models exist in the literature at the subsystems level, a comprehensive model for the complete system including the important regulatory feedback loops is missing. Our dynamic model combines the existing subsystem models and studies their steady-state and dynamic interactions under feedback. We establish conditions under which bistability and oscillations exist for this important pathway. In particular, we show how the negative and positive feedback loops affect the dynamic characteristics that determine the cellular outcome.
Collapse
Affiliation(s)
- Yaman Arkun
- Department of Chemical and Biological Engineering, Koc University, Rumeli Feneri Yolu, Sariyer, Istanbul, Turkey
- * E-mail:
| | - Mohammadreza Yasemi
- Department of Chemical and Biological Engineering, Koc University, Rumeli Feneri Yolu, Sariyer, Istanbul, Turkey
| |
Collapse
|
62
|
Bouhaddou M, Barrette AM, Stern AD, Koch RJ, DiStefano MS, Riesel EA, Santos LC, Tan AL, Mertz AE, Birtwistle MR. A mechanistic pan-cancer pathway model informed by multi-omics data interprets stochastic cell fate responses to drugs and mitogens. PLoS Comput Biol 2018; 14:e1005985. [PMID: 29579036 PMCID: PMC5886578 DOI: 10.1371/journal.pcbi.1005985] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Revised: 04/05/2018] [Accepted: 01/16/2018] [Indexed: 01/02/2023] Open
Abstract
Most cancer cells harbor multiple drivers whose epistasis and interactions with expression context clouds drug and drug combination sensitivity prediction. We constructed a mechanistic computational model that is context-tailored by omics data to capture regulation of stochastic proliferation and death by pan-cancer driver pathways. Simulations and experiments explore how the coordinated dynamics of RAF/MEK/ERK and PI-3K/AKT kinase activities in response to synergistic mitogen or drug combinations control cell fate in a specific cellular context. In this MCF10A cell context, simulations suggest that synergistic ERK and AKT inhibitor-induced death is likely mediated by BIM rather than BAD, which is supported by prior experimental studies. AKT dynamics explain S-phase entry synergy between EGF and insulin, but simulations suggest that stochastic ERK, and not AKT, dynamics seem to drive cell-to-cell proliferation variability, which in simulations is predictable from pre-stimulus fluctuations in C-Raf/B-Raf levels. Simulations suggest MEK alteration negligibly influences transformation, consistent with clinical data. Tailoring the model to an alternate cell expression and mutation context, a glioma cell line, allows prediction of increased sensitivity of cell death to AKT inhibition. Our model mechanistically interprets context-specific landscapes between driver pathways and cell fates, providing a framework for designing more rational cancer combination therapy. Cancer is a complex and diverse disease. Two people with the same cancer type often respond differently to the same treatment. These differences are primarily driven by the fact that two type-matched tumors can possess distinct sets of mutations and gene expression profiles, provoking differential sensitivity to drugs. Over the past few decades, we have seen a shift away from more broadly cytotoxic drugs to more targeted molecules therapies; but how to match a patient with a specific drug or drug cocktail remains a difficult problem. Here, we build a mechanistic ordinary differential equation model describing the interactions between commonly mutated pan-cancer signaling pathways—receptor tyrosine kinases, Ras/RAF/ERK, PI3K/AKT, mTOR, cell cycle, DNA damage, and apoptosis. We develop methods for how to tailor the model to multi-omics data from a specific biological context, devise a novel stochastic algorithm to induce non-genetic cell-to-cell fluctuations in mRNA and protein quantities over time, and train the model against a wealth of biochemical and cell fate data to gain insight into the systems-level, context-specific control of proliferation and death. One day, we hope models of this kind could be tailored to patient-derived tumor mRNA sequencing data and used to prioritize patient-specific drug regimens.
Collapse
Affiliation(s)
- Mehdi Bouhaddou
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Anne Marie Barrette
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Alan D. Stern
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Rick J. Koch
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Matthew S. DiStefano
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Eric A. Riesel
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Luis C. Santos
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Annie L. Tan
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Alex E. Mertz
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
| | - Marc R. Birtwistle
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States of America
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, SC, United States of America
- * E-mail:
| |
Collapse
|
63
|
Brogden KA, Parashar D, Hallier AR, Braun T, Qian F, Rizvi NA, Bossler AD, Milhem MM, Chan TA, Abbasi T, Vali S. Genomics of NSCLC patients both affirm PD-L1 expression and predict their clinical responses to anti-PD-1 immunotherapy. BMC Cancer 2018; 18:225. [PMID: 29486723 PMCID: PMC5897943 DOI: 10.1186/s12885-018-4134-y] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2017] [Accepted: 02/15/2018] [Indexed: 12/21/2022] Open
Abstract
Background Programmed Death Ligand 1 (PD-L1) is a co-stimulatory and immune checkpoint protein. PD-L1 expression in non-small cell lung cancers (NSCLC) is a hallmark of adaptive resistance and its expression is often used to predict the outcome of Programmed Death 1 (PD-1) and PD-L1 immunotherapy treatments. However, clinical benefits do not occur in all patients and new approaches are needed to assist in selecting patients for PD-1 or PD-L1 immunotherapies. Here, we hypothesized that patient tumor cell genomics influenced cell signaling and expression of PD-L1, chemokines, and immunosuppressive molecules and these profiles could be used to predict patient clinical responses. Methods We used a recent dataset from NSCLC patients treated with pembrolizumab. Deleterious gene mutational profiles in patient exomes were identified and annotated into a cancer network to create NSCLC patient-specific predictive computational simulation models. Validation checks were performed on the cancer network, simulation model predictions, and PD-1 match rates between patient-specific predicted and clinical responses. Results Expression profiles of these 24 chemokines and immunosuppressive molecules were used to identify patients who would or would not respond to PD-1 immunotherapy. PD-L1 expression alone was not sufficient to predict which patients would or would not respond to PD-1 immunotherapy. Adding chemokine and immunosuppressive molecule expression profiles allowed patient models to achieve a greater than 85.0% predictive correlation among predicted and reported patient clinical responses. Conclusions Our results suggested that chemokine and immunosuppressive molecule expression profiles can be used to accurately predict clinical responses thus differentiating among patients who would and would not benefit from PD-1 or PD-L1 immunotherapies. Electronic supplementary material The online version of this article (10.1186/s12885-018-4134-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Kim A Brogden
- Iowa Institute for Oral Health Research, College of Dentistry, The University of Iowa, 801 Newton Road, Iowa City, IA, 52242, USA.
| | - Deepak Parashar
- Cellworks Research India Ltd., Whitefield, Bangalore, 560066, India
| | - Andrea R Hallier
- Biomedical Engineering, The University of Iowa, 5318 SC, Iowa City, IA, 52242, USA
| | - Terry Braun
- Biomedical Engineering, The University of Iowa, 5318 SC, Iowa City, IA, 52242, USA
| | - Fang Qian
- Iowa Institute for Oral Health Research, College of Dentistry, The University of Iowa, 801 Newton Road, Iowa City, IA, 52242, USA.,Division of Biostatistics and Research Design, College of Dentistry, The University of Iowa, 801 Newton Road, Iowa City, IA, 52242, USA
| | - Naiyer A Rizvi
- Division of Hematology/Oncology, Columbia University Medical Center, 177 Fort Washington Avenue, New York, NY, 10032, USA
| | - Aaron D Bossler
- Molecular Pathology Laboratory, Department of Pathology, University of Iowa Hospitals and Clinics, 200 Hawkins Dr., C606GH, Iowa City, IA, 52242, USA
| | - Mohammed M Milhem
- Clinical Services, Experimental Therapeutics, Melanoma and Sarcoma Program, Holden Comprehensive Cancer Center, The University of Iowa, Iowa City, IA, 52242, USA
| | - Timothy A Chan
- Department of Radiation Oncology, Human Oncology and Pathogenesis Program, Immunogenomics and Precision Oncology Platform, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
| | - Taher Abbasi
- Cellworks Group, Inc., 2033 Gateway Place Suite 500, San Jose, CA, 95110, USA
| | - Shireen Vali
- Cellworks Group, Inc., 2033 Gateway Place Suite 500, San Jose, CA, 95110, USA
| |
Collapse
|
64
|
Varusai TM, Nguyen LK. Dynamic modelling of the mTOR signalling network reveals complex emergent behaviours conferred by DEPTOR. Sci Rep 2018; 8:643. [PMID: 29330362 PMCID: PMC5766521 DOI: 10.1038/s41598-017-18400-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Accepted: 12/01/2017] [Indexed: 01/10/2023] Open
Abstract
The mechanistic Target of Rapamycin (mTOR) signalling network is an evolutionarily conserved network that controls key cellular processes, including cell growth and metabolism. Consisting of the major kinase complexes mTOR Complex 1 and 2 (mTORC1/2), the mTOR network harbours complex interactions and feedback loops. The DEP domain-containing mTOR-interacting protein (DEPTOR) was recently identified as an endogenous inhibitor of both mTORC1 and 2 through direct interactions, and is in turn degraded by mTORC1/2, adding an extra layer of complexity to the mTOR network. Yet, the dynamic properties of the DEPTOR-mTOR network and the roles of DEPTOR in coordinating mTORC1/2 activation dynamics have not been characterised. Using computational modelling, systems analysis and dynamic simulations we show that DEPTOR confers remarkably rich and complex dynamic behaviours to mTOR signalling, including abrupt, bistable switches, oscillations and co-existing bistable/oscillatory responses. Transitions between these distinct modes of behaviour are enabled by modulating DEPTOR expression alone. We characterise the governing conditions for the observed dynamics by elucidating the network in its vast multi-dimensional parameter space, and develop strategies to identify core network design motifs underlying these dynamics. Our findings provide new systems-level insights into the complexity of mTOR signalling contributed by DEPTOR.
Collapse
Affiliation(s)
- Thawfeek M Varusai
- European Bioinformatics Institute, EMBL-EBI, Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK.,Systems Biology Ireland, Conway Institute, University College Dublin, Belfield, Dublin 4, Ireland
| | - Lan K Nguyen
- Department of Biochemistry and Molecular Biology, Monash University, Melbourne, Victoria, 3800, Australia. .,Biomedicine Discovery Institute, Monash University, Melbourne, Victoria, 3800, Australia. .,Systems Biology Ireland, Conway Institute, University College Dublin, Belfield, Dublin 4, Ireland.
| |
Collapse
|
65
|
Caldieri G, Malabarba MG, Di Fiore PP, Sigismund S. EGFR Trafficking in Physiology and Cancer. PROGRESS IN MOLECULAR AND SUBCELLULAR BIOLOGY 2018; 57:235-272. [PMID: 30097778 DOI: 10.1007/978-3-319-96704-2_9] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Signaling from the epidermal growth factor receptor (EGFR) elicits multiple biological responses, including cell proliferation, migration, and survival. Receptor endocytosis and trafficking are critical physiological processes that control the strength, duration, diversification, and spatial restriction of EGFR signaling through multiple mechanisms, which we review in this chapter. These mechanisms include: (i) regulation of receptor density and activation at the cell surface; (ii) concentration of receptors into distinct nascent endocytic structures; (iii) commitment of the receptor to different endocytic routes; (iv) endosomal sorting and postendocytic trafficking of the receptor through distinct pathways, and (v) recycling to restricted regions of the cell surface. We also highlight how communication between organelles controls EGFR activity along the endocytic route. Finally, we illustrate how abnormal trafficking of EGFR oncogenic mutants, as well as alterations of the endocytic machinery, contributes to aberrant EGFR signaling in cancer.
Collapse
Affiliation(s)
- Giusi Caldieri
- Dipartimento di Oncologia ed Emato-oncologia, Università degli Studi di Milano, Via Santa Sofia 9/1, 20122, Milan, Italy
- Istituto Europeo di Oncologia, Via Ripamonti 435, 20141, Milan, Italy
| | - Maria Grazia Malabarba
- Dipartimento di Oncologia ed Emato-oncologia, Università degli Studi di Milano, Via Santa Sofia 9/1, 20122, Milan, Italy
- Istituto Europeo di Oncologia, Via Ripamonti 435, 20141, Milan, Italy
| | - Pier Paolo Di Fiore
- Dipartimento di Oncologia ed Emato-oncologia, Università degli Studi di Milano, Via Santa Sofia 9/1, 20122, Milan, Italy
- Istituto Europeo di Oncologia, Via Ripamonti 435, 20141, Milan, Italy
| | - Sara Sigismund
- Dipartimento di Oncologia ed Emato-oncologia, Università degli Studi di Milano, Via Santa Sofia 9/1, 20122, Milan, Italy.
- Istituto Europeo di Oncologia, Via Ripamonti 435, 20141, Milan, Italy.
| |
Collapse
|
66
|
Behinaein B, Rudie K, Sangrar W. Petri net siphon analysis and graph theoretic measures for identifying combination therapies in cancer. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:231-243. [PMID: 28113516 DOI: 10.1109/tcbb.2016.2614301] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Epidermal Growth Factor Receptor (EGFR) signaling to the Ras-MAPK pathway is implicated in the development and progression of cancer and is a major focus of targeted combination therapies. Physiochemical models have been used for identifying and testing the signal-inhibiting potential of targeted therapies, however, their application to larger multi-pathway networks is limited by the availability of experimentally-determined rate and concentration parameters. An alternate strategy for identifying and evaluating drug-targetable nodes is proposed. A physiochemical model of EGFR-Ras-MAPK signaling is implemented and calibrated to experimental data. Essential topological features of the model are converted into a Petri net and nodes that behave as siphons-a structural property of Petri nets-are identified. Siphons represent potential drug-targets since they are unrecoverable if their values fall below a threshold. Centrality measures are then used to prioritize siphons identified as candidate drug-targets. Single and multiple drug-target combinations are identified which correspond to clinically relevant drug targets and exhibit inhibition synergy in physiochemical simulations of EGF-induced EGFR-Ras-MAPK signaling. Taken together, these studies suggest that siphons and centrality analyses are a promising computational strategy to identify and rank drug-targetable nodes in larger networks as they do not require knowledge of the dynamics of the system, but rely solely on topology.
Collapse
|
67
|
Rocha HL, Almeida RC, Lima EABF, Resende ACM, Oden JT, Yankeelov TE. A HYBRID THREE-SCALE MODEL OF TUMOR GROWTH. MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES : M3AS 2018; 28:61-93. [PMID: 29353950 PMCID: PMC5773147 DOI: 10.1142/s0218202518500021] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Cancer results from a complex interplay of different biological, chemical, and physical phenomena that span a wide range of time and length scales. Computational modeling may help to unfold the role of multiple evolving factors that exist and interact in the tumor microenvironment. Understanding these complex multiscale interactions is a crucial step towards predicting cancer growth and in developing effective therapies. We integrate different modeling approaches in a multiscale, avascular, hybrid tumor growth model encompassing tissue, cell, and sub-cell scales. At the tissue level, we consider the dispersion of nutrients and growth factors in the tumor microenvironment, which are modeled through reaction-diffusion equations. At the cell level, we use an agent based model (ABM) to describe normal and tumor cell dynamics, with normal cells kept in homeostasis and cancer cells differentiated apoptotic, hypoxic, and necrotic states. Cell movement is driven by the balance of a variety of forces according to Newton's second law, including those related to growth-induced stresses. Phenotypic transitions are defined by specific rule of behaviors that depend on microenvironment stimuli. We integrate in each cell/agent a branch of the epidermal growth factor receptor (EGFR) pathway. This pathway is modeled by a system of coupled nonlinear differential equations involving the mass laws of 20 molecules. The rates of change in the concentration of some key molecules trigger proliferation or migration advantage response. The bridge between cell and tissue scales is built through the reaction and source terms of the partial differential equations. Our hybrid model is built in a modular way, enabling the investigation of the role of different mechanisms at multiple scales on tumor progression. This strategy allows representating both the collective behavior due to cell assembly as well as microscopic intracellular phenomena described by signal transduction pathways. Here, we investigate the impact of some mechanisms associated with sustained proliferation on cancer progression. Specifically, we focus on the intracellular proliferation/migration-advantage-response driven by the EGFR pathway and on proliferation inhibition due to accumulation of growth-induced stresses. Simulations demonstrate that the model can adequately describe some complex mechanisms of tumor dynamics, including growth arrest in avascular tumors. Both the sub-cell model and growth-induced stresses give rise to heterogeneity in the tumor expansion and a rich variety of tumor behaviors.
Collapse
Affiliation(s)
- H L Rocha
- National Laboratory for Scientific Computing (LNCC), Av. Getúlio Vargas, 333, Quitandinha, Petrópolis, Rio de Janeiro, 25651-075, Brazil
| | - R C Almeida
- National Laboratory for Scientific Computing (LNCC), Av. Getúlio Vargas, 333, Quitandinha, Petrópolis, Rio de Janeiro, 25651-075, Brazil
| | - E A B F Lima
- Center of Computational Oncology, Institute for Computational Engineering and Sciences, The University of Texas at Austin, 201 East 24th St, Austin, TX, 78712-1229, USA
| | - A C M Resende
- National Laboratory for Scientific Computing (LNCC), Av. Getúlio Vargas, 333, Quitandinha, Petrópolis, Rio de Janeiro, 25651-075, Brazil
| | - J T Oden
- Center of Computational Oncology, Institute for Computational Engineering and Sciences, The University of Texas at Austin, 201 East 24th St, Austin, TX, 78712-1229, USA
| | - T E Yankeelov
- Center of Computational Oncology, Institute for Computational Engineering and Sciences, The University of Texas at Austin, 201 East 24th St, Austin, TX, 78712-1229, USA
- Department of Biomedical Engineering, The University of Texas at Austin, 107 W. Dean Keeton, Austin, TX, 78712, USA
- Department of Internal Medicine, Livestrong Cancer Institutes, Dell Medical School, The University of Texas at Austin
| |
Collapse
|
68
|
Weddell JC, Chen S, Imoukhuede PI. VEGFR1 promotes cell migration and proliferation through PLCγ and PI3K pathways. NPJ Syst Biol Appl 2017; 4:1. [PMID: 29263797 PMCID: PMC5736688 DOI: 10.1038/s41540-017-0037-9] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 11/08/2017] [Accepted: 11/21/2017] [Indexed: 12/16/2022] Open
Abstract
The ability to control vascular endothelial growth factor (VEGF) signaling offers promising therapeutic potential for vascular diseases and cancer. Despite this promise, VEGF-targeted therapies are not clinically effective for many pathologies, such as breast cancer. VEGFR1 has recently emerged as a predictive biomarker for anti-VEGF efficacy, implying a functional VEGFR1 role beyond its classically defined decoy receptor status. Here we introduce a computational approach that accurately predicts cellular responses elicited via VEGFR1 signaling. Aligned with our model prediction, we show empirically that VEGFR1 promotes macrophage migration through PLCγ and PI3K pathways and promotes macrophage proliferation through a PLCγ pathway. These results provide new insight into the basic function of VEGFR1 signaling while offering a computational platform to quantify signaling of any receptor.
Collapse
Affiliation(s)
- Jared C. Weddell
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| | - Si Chen
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| | - P. I. Imoukhuede
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| |
Collapse
|
69
|
Bahrami A, Miraie-Ashtiani SR, Sadeghi M, Najafi A, Ranjbar R. Dynamic modeling of folliculogenesis signaling pathways in the presence of miRNAs expression. J Ovarian Res 2017; 10:76. [PMID: 29258623 PMCID: PMC5735818 DOI: 10.1186/s13048-017-0371-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2017] [Accepted: 11/23/2017] [Indexed: 11/10/2022] Open
Abstract
Background TEK signaling plays a very important role in folliculogenesis. It activates Ras/ERK/MYC, PI3K/AKT/mTORC1 and ovarian steroidogenesis activation pathways. These are the main pathways for cell growth, differentiation, migration, adhesion, proliferation, survival and protein synthesis. Results TEK signaling on each of the two important pathways where levels of pERK, pMYC, pAkt, pMCL1 and pEIF4EBP1 are increased in dominant follicles and pMYC is decreased in dominant follicles. Over activation of ERK and MYC which are the main cell growth and proliferation and over activation of Akt, MCl1, mTORC1 and EIF4EBP1 which are the main cell survival and protein synthesis factors act as promoting factors for folliculogenesis. In case of over expression of hsa-miR-30d-3p and hsa-miR-451a, MYC activity level is considerably increased in subordinate follicles. Our simulation results show that in the presence of has-miR-548v and bta-miR-22-3p, downstream factors of pathways are inhibited. Conclusions Our work offers insight into the design of natural biological procedures and makes predictions that can guide further experimental studies on folliculogenesis pathways. Moreover, it defines a simple signal processing unit that may be useful for engineering synthetic biology and genes circuits to carry out cell-based computation. Electronic supplementary material The online version of this article (doi:10.1186/s13048-017-0371-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Abolfazl Bahrami
- Department of Animal Science, University college of Agriculture and Natural Resources, University of Tehran, Karaj, Iran.
| | - Seyed Reza Miraie-Ashtiani
- Department of Animal Science, University college of Agriculture and Natural Resources, University of Tehran, Karaj, Iran.
| | - Mostafa Sadeghi
- Department of Animal Science, University college of Agriculture and Natural Resources, University of Tehran, Karaj, Iran
| | - Ali Najafi
- Molecular Biology Research Center, Baqiyatallah University of Medical Sciences, Tehran, Iran
| | - Reza Ranjbar
- Molecular Biology Research Center, Baqiyatallah University of Medical Sciences, Tehran, Iran.
| |
Collapse
|
70
|
Wang YM, Gu ML, Meng FS, Jiao WR, Zhou XX, Yao HP, Ji F. Histone acetyltransferase p300/CBP inhibitor C646 blocks the survival and invasion pathways of gastric cancer cell lines. Int J Oncol 2017; 51:1860-1868. [PMID: 29075795 DOI: 10.3892/ijo.2017.4176] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2017] [Accepted: 09/15/2017] [Indexed: 12/31/2022] Open
Abstract
The histone acetyltransferases (HATs) adenovirus E1A-associated protein (p300) and CREB binding protein (CBP) serve as coactivators during a diverse assortment of cellular processes. In the present study, p300 and CBP were highly expressed in 5 gastric cancer (GC) cell lines (SGC‑7901, MKN45, MGC-803, BGC-823 and KATO III) compared with human normal gastric epithelial cell line (GES-1). C646, a selective inhibitor of p300 and CBP, inhibited cell viability and cell cycle and promoted cell apoptosis in all 5 GC cell lines. In addition, C646 suppressed the migration and invasion capability of the GC cell lines, except for the middle-differentiated SGC-7901 cell line. Furthermore, we detected the differential expression of corresponding oncogenic signalling molecules, such as c-Met, Akt, Bcl-2, Bax, cyclin D1, MMP7 and MMP9, in GC cells following C646 treatment. In conclusion, our results suggest that C646 inhibits the acetylation of histone H3 via inactivation of p300 and CBP, resulting in antineoplastic effects toward GC cells. Thus, the selective HAT inhibitor C646 could be a promising antitumour reagent for GC treatment.
Collapse
Affiliation(s)
- Ya-Mei Wang
- Department of Gastroenterology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| | - Meng-Li Gu
- Department of Gastroenterology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| | - Fan-Sheng Meng
- Department of Gastroenterology, Linyi People's Hospital, Linyi, Shandong 276000, P.R. China
| | - Wen-Rui Jiao
- Department of Gastroenterology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| | - Xin-Xin Zhou
- Department of Gastroenterology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| | - Hang-Ping Yao
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| | - Feng Ji
- Department of Gastroenterology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang 310009, P.R. China
| |
Collapse
|
71
|
Freed DM, Bessman NJ, Kiyatkin A, Salazar-Cavazos E, Byrne PO, Moore JO, Valley CC, Ferguson KM, Leahy DJ, Lidke DS, Lemmon MA. EGFR Ligands Differentially Stabilize Receptor Dimers to Specify Signaling Kinetics. Cell 2017; 171:683-695.e18. [PMID: 28988771 DOI: 10.1016/j.cell.2017.09.017] [Citation(s) in RCA: 230] [Impact Index Per Article: 32.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2017] [Revised: 07/31/2017] [Accepted: 09/12/2017] [Indexed: 12/21/2022]
Abstract
Epidermal growth factor receptor (EGFR) regulates many crucial cellular programs, with seven different activating ligands shaping cell signaling in distinct ways. Using crystallography and other approaches, we show how the EGFR ligands epiregulin (EREG) and epigen (EPGN) stabilize different dimeric conformations of the EGFR extracellular region. As a consequence, EREG or EPGN induce less stable EGFR dimers than EGF-making them partial agonists of EGFR dimerization. Unexpectedly, this weakened dimerization elicits more sustained EGFR signaling than seen with EGF, provoking responses in breast cancer cells associated with differentiation rather than proliferation. Our results reveal how responses to different EGFR ligands are defined by receptor dimerization strength and signaling dynamics. These findings have broad implications for understanding receptor tyrosine kinase (RTK) signaling specificity. Our results also suggest parallels between partial and/or biased agonism in RTKs and G-protein-coupled receptors, as well as new therapeutic opportunities for correcting RTK signaling output.
Collapse
Affiliation(s)
- Daniel M Freed
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA; Yale Cancer Biology Institute, Yale University, West Haven, CT 06516, USA
| | - Nicholas J Bessman
- Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, PA 19104-6059, USA
| | - Anatoly Kiyatkin
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA; Yale Cancer Biology Institute, Yale University, West Haven, CT 06516, USA
| | - Emanuel Salazar-Cavazos
- Department of Pathology and UNM Comprehensive Cancer Center, University of New Mexico Health Science Center, Albuquerque, NM 87131, USA
| | - Patrick O Byrne
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712, USA
| | - Jason O Moore
- Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, PA 19104-6059, USA
| | - Christopher C Valley
- Department of Pathology and UNM Comprehensive Cancer Center, University of New Mexico Health Science Center, Albuquerque, NM 87131, USA
| | - Kathryn M Ferguson
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA; Yale Cancer Biology Institute, Yale University, West Haven, CT 06516, USA
| | - Daniel J Leahy
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712, USA
| | - Diane S Lidke
- Department of Pathology and UNM Comprehensive Cancer Center, University of New Mexico Health Science Center, Albuquerque, NM 87131, USA
| | - Mark A Lemmon
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA; Yale Cancer Biology Institute, Yale University, West Haven, CT 06516, USA; Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, PA 19104-6059, USA.
| |
Collapse
|
72
|
Deisboeck TS, Zhang L, Martin S. Advancing Cancer Systems Biology: Introducing the Center for the Development of a Virtual Tumor, CViT. Cancer Inform 2017. [DOI: 10.1177/117693510700500001] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Integrative cancer biology research relies on a variety of data-driven computational modeling and simulation methods and techniques geared towards gaining new insights into the complexity of biological processes that are of critical importance for cancer research. These include the dynamics of gene-protein interaction networks, the percolation of sub-cellular perturbations across scales and the impact they may have on tumorigenesis in both experiments and clinics. Such innovative ‘systems’ research will greatly benefit from enabling Information Technology that is currently under development, including an online collaborative environment, a Semantic Web based computing platform that hosts data and model repositories as well as high-performance computing access. Here, we present one of the National Cancer Institute's recently established Integrative Cancer Biology Programs, i.e. the Center for the Development of a Virtual Tumor, CViT, which is charged with building a cancer modeling community, developing the aforementioned enabling technologies and fostering multi-scale cancer modeling and simulation.
Collapse
Affiliation(s)
- Thomas S. Deisboeck
- Complex Biosystems Modeling Laboratory, Harvard-MIT (HST) Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA 02129
| | - Le Zhang
- Complex Biosystems Modeling Laboratory, Harvard-MIT (HST) Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA 02129
| | - Sean Martin
- IBM Advanced Internet Technology, Cambridge, MA 02142
| |
Collapse
|
73
|
Zangooei MH, Habibi J. Hybrid multiscale modeling and prediction of cancer cell behavior. PLoS One 2017; 12:e0183810. [PMID: 28846712 PMCID: PMC5573302 DOI: 10.1371/journal.pone.0183810] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Accepted: 08/13/2017] [Indexed: 12/03/2022] Open
Abstract
Background Understanding cancer development crossing several spatial-temporal scales is of great practical significance to better understand and treat cancers. It is difficult to tackle this challenge with pure biological means. Moreover, hybrid modeling techniques have been proposed that combine the advantages of the continuum and the discrete methods to model multiscale problems. Methods In light of these problems, we have proposed a new hybrid vascular model to facilitate the multiscale modeling and simulation of cancer development with respect to the agent-based, cellular automata and machine learning methods. The purpose of this simulation is to create a dataset that can be used for prediction of cell phenotypes. By using a proposed Q-learning based on SVR-NSGA-II method, the cells have the capability to predict their phenotypes autonomously that is, to act on its own without external direction in response to situations it encounters. Results Computational simulations of the model were performed in order to analyze its performance. The most striking feature of our results is that each cell can select its phenotype at each time step according to its condition. We provide evidence that the prediction of cell phenotypes is reliable. Conclusion Our proposed model, which we term a hybrid multiscale modeling of cancer cell behavior, has the potential to combine the best features of both continuum and discrete models. The in silico results indicate that the 3D model can represent key features of cancer growth, angiogenesis, and its related micro-environment and show that the findings are in good agreement with biological tumor behavior. To the best of our knowledge, this paper is the first hybrid vascular multiscale modeling of cancer cell behavior that has the capability to predict cell phenotypes individually by a self-generated dataset.
Collapse
Affiliation(s)
| | - Jafar Habibi
- Department of Computer Engineering, Sharif University of Technology, Tehran, Iran
- * E-mail:
| |
Collapse
|
74
|
Borisov N, Suntsova M, Sorokin M, Garazha A, Kovalchuk O, Aliper A, Ilnitskaya E, Lezhnina K, Korzinkin M, Tkachev V, Saenko V, Saenko Y, Sokov DG, Gaifullin NM, Kashintsev K, Shirokorad V, Shabalina I, Zhavoronkov A, Mishra B, Cantor CR, Buzdin A. Data aggregation at the level of molecular pathways improves stability of experimental transcriptomic and proteomic data. Cell Cycle 2017; 16:1810-1823. [PMID: 28825872 DOI: 10.1080/15384101.2017.1361068] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
High throughput technologies opened a new era in biomedicine by enabling massive analysis of gene expression at both RNA and protein levels. Unfortunately, expression data obtained in different experiments are often poorly compatible, even for the same biologic samples. Here, using experimental and bioinformatic investigation of major experimental platforms, we show that aggregation of gene expression data at the level of molecular pathways helps to diminish cross- and intra-platform bias otherwise clearly seen at the level of individual genes. We created a mathematical model of cumulative suppression of data variation that predicts the ideal parameters and the optimal size of a molecular pathway. We compared the abilities to aggregate experimental molecular data for the 5 alternative methods, also evaluated by their capacity to retain meaningful features of biologic samples. The bioinformatic method OncoFinder showed optimal performance in both tests and should be very useful for future cross-platform data analyses.
Collapse
Affiliation(s)
- Nicolas Borisov
- a Centre for Convergence of Nano-, Bio-, Information and Cognitive Sciences and Technologies, National Research Centre "Kurchatov Institute" , Moscow , Russia.,b Department of R&D, First Oncology Research and Advisory Center , Moscow , Russia
| | - Maria Suntsova
- c Department of R&D, Center for Biogerontology and Regenerative Medicine , Moscow , Russia.,d Laboratory of Bioinformatics, D. Rogachyov Federal Research Center of Pediatric Hematology, Oncology and Immunology , Moscow , Russia
| | - Maxim Sorokin
- a Centre for Convergence of Nano-, Bio-, Information and Cognitive Sciences and Technologies, National Research Centre "Kurchatov Institute" , Moscow , Russia.,e Group for Genomic Regulation of Cell Signaling Systems, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry , Moscow , Russia
| | - Andrew Garazha
- c Department of R&D, Center for Biogerontology and Regenerative Medicine , Moscow , Russia.,f Department of R&D, OmicsWay Corporation , Walnut , CA , USA
| | - Olga Kovalchuk
- g Department of Biological Sciences , University of Lethbridge , Lethbridge , AB , Canada
| | - Alexander Aliper
- d Laboratory of Bioinformatics, D. Rogachyov Federal Research Center of Pediatric Hematology, Oncology and Immunology , Moscow , Russia
| | - Elena Ilnitskaya
- c Department of R&D, Center for Biogerontology and Regenerative Medicine , Moscow , Russia
| | - Ksenia Lezhnina
- b Department of R&D, First Oncology Research and Advisory Center , Moscow , Russia
| | - Mikhail Korzinkin
- c Department of R&D, Center for Biogerontology and Regenerative Medicine , Moscow , Russia
| | - Victor Tkachev
- f Department of R&D, OmicsWay Corporation , Walnut , CA , USA
| | - Vyacheslav Saenko
- h Technological Research Institute S.P. Kapitsa , Ulyanovsk State University , Ulyanovsk , Russia
| | - Yury Saenko
- h Technological Research Institute S.P. Kapitsa , Ulyanovsk State University , Ulyanovsk , Russia
| | - Dmitry G Sokov
- i Chemotherapy Department, Moscow 1st Oncological Hospital , Moscow , Russia
| | - Nurshat M Gaifullin
- j Faculty of Fundamental Medicine , Lomonosov Moscow State University , Moscow , Russia.,k Department of Oncology, Russian Medical Postgraduate Academy , Moscow , Russia
| | - Kirill Kashintsev
- l Chemotherapy Department, Moscow Oncological Hospital 62 , Stepanovskoye , Russia
| | - Valery Shirokorad
- l Chemotherapy Department, Moscow Oncological Hospital 62 , Stepanovskoye , Russia
| | - Irina Shabalina
- m Faculty of Mathematics and Information Technologies , Petrozavodsk State University , Petrozavodsk , Russia
| | - Alex Zhavoronkov
- d Laboratory of Bioinformatics, D. Rogachyov Federal Research Center of Pediatric Hematology, Oncology and Immunology , Moscow , Russia
| | | | - Charles R Cantor
- o Department of Biomedical Engineering , Boston University , Boston , MA , USA
| | - Anton Buzdin
- a Centre for Convergence of Nano-, Bio-, Information and Cognitive Sciences and Technologies, National Research Centre "Kurchatov Institute" , Moscow , Russia.,b Department of R&D, First Oncology Research and Advisory Center , Moscow , Russia.,e Group for Genomic Regulation of Cell Signaling Systems, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry , Moscow , Russia.,f Department of R&D, OmicsWay Corporation , Walnut , CA , USA
| |
Collapse
|
75
|
Getz MC, Nirody JA, Rangamani P. Stability analysis in spatial modeling of cell signaling. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2017; 10. [PMID: 28787545 DOI: 10.1002/wsbm.1395] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Revised: 06/06/2017] [Accepted: 06/11/2017] [Indexed: 11/08/2022]
Abstract
Advances in high-resolution microscopy and other techniques have emphasized the spatio-temporal nature of information transfer through signal transduction pathways. The compartmentalization of signaling molecules and the existence of microdomains are now widely acknowledged as key features in biochemical signaling. To complement experimental observations of spatio-temporal dynamics, mathematical modeling has emerged as a powerful tool. Using modeling, one can not only recapitulate experimentally observed dynamics of signaling molecules, but also gain an understanding of the underlying mechanisms in order to generate experimentally testable predictions. Reaction-diffusion systems are commonly used to this end; however, the analysis of coupled nonlinear systems of partial differential equations, generated by considering large reaction networks is often challenging. Here, we aim to provide an introductory tutorial for the application of reaction-diffusion models to the spatio-temporal dynamics of signaling pathways. In particular, we outline the steps for stability analysis of such models, with a focus on biochemical signal transduction. WIREs Syst Biol Med 2018, 10:e1395. doi: 10.1002/wsbm.1395 This article is categorized under: Biological Mechanisms > Cell Signaling Analytical and Computational Methods > Dynamical Methods Models of Systems Properties and Processes > Mechanistic Models.
Collapse
Affiliation(s)
- Michael C Getz
- Chemical Engineering Graduate Program, University of California San Diego, La Jolla, CA, USA
| | - Jasmine A Nirody
- Biophysics Graduate Program, University of California, Berkeley, CA, USA
| | - Padmini Rangamani
- Department of Mechanical and Aerospace Engineering, University of California San Diego, La Jolla, CA, USA
| |
Collapse
|
76
|
Bravi B, Sollich P. Statistical physics approaches to subnetwork dynamics in biochemical systems. Phys Biol 2017; 14:045010. [PMID: 28510539 DOI: 10.1088/1478-3975/aa7363] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We apply a Gaussian variational approximation to model reduction in large biochemical networks of unary and binary reactions. We focus on a small subset of variables (subnetwork) of interest, e.g. because they are accessible experimentally, embedded in a larger network (bulk). The key goal is to write dynamical equations reduced to the subnetwork but still retaining the effects of the bulk. As a result, the subnetwork-reduced dynamics contains a memory term and an extrinsic noise term with non-trivial temporal correlations. We first derive expressions for this memory and noise in the linearized (Gaussian) dynamics and then use a perturbative power expansion to obtain first order nonlinear corrections. For the case of vanishing intrinsic noise, our description is explicitly shown to be equivalent to projection methods up to quadratic terms, but it is applicable also in the presence of stochastic fluctuations in the original dynamics. An example from the epidermal growth factor receptor signalling pathway is provided to probe the increased prediction accuracy and computational efficiency of our method.
Collapse
Affiliation(s)
- B Bravi
- Current affiliation: Institute of Theoretical Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
77
|
Hochman G, Halevi-Tobias K, Kogan Y, Agur Z. Extracellular inhibitors can attenuate tumorigenic Wnt pathway activity in adenomatous polyposis coli mutants: Predictions of a validated mathematical model. PLoS One 2017; 12:e0179888. [PMID: 28708837 PMCID: PMC5510801 DOI: 10.1371/journal.pone.0179888] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2017] [Accepted: 06/06/2017] [Indexed: 12/21/2022] Open
Abstract
Background Despite considerable investigational efforts, no method to overcome the pathogenesis caused by loss of function (LoF) mutations in tumor suppressor genes has been successfully translated to the clinic. The most frequent LoF mutation in human cancers is Adenomatous polyposis coli (APC), causing aberrant activation of the Wnt pathway. In nearly all colon cancer tumors, the APC protein is truncated, but still retains partial binding abilities. Objective & methods Here, we tested the hypothesis that extracellular inhibitors of the Wnt pathway, although acting upstream of the APC mutation, can restore normal levels of pathway activity in colon cancer cells. To this end, we developed and simulated a mathematical model for the Wnt pathway in different APC mutants, with or without the effects of the extracellular inhibitors, Secreted Frizzled-Related Protein1 (sFRP1) and Dickhopf1 (Dkk1). We compared our model predictions to experimental data in the literature. Results Our model accurately predicts T-cell factor (TCF) activity in mutant cells that vary in APC mutation. Model simulations suggest that both sFRP1 and DKK1 can reduce TCF activity in APC1638N/1572T and Apcmin/min mutants, but restoration of normal activity levels is possible only in the former. When applied in combination, synergism between the two inhibitors can reduce their effective doses to one-fourth of the doses required under single inhibitor application. Overall, re-establishment of normal Wnt pathway activity is predicted for every APC mutant in whom TCF activity is increased by up to 11 fold. Conclusions Our work suggests that extracellular inhibitors can effectively restore normal Wnt pathway activity in APC-truncated cancer cells, even though these LoF mutations occur downstream of the inhibitory action. The insufficient activity of the truncated APC can be quantitatively balanced by the upstream intervention. This new concept of upstream intervention to control the effects of downstream mutations may be considered also for other partial LoF mutations in other signaling pathways.
Collapse
Affiliation(s)
- Gili Hochman
- Institute for Medical BioMathematics, Bene Ataroth, Israel
| | | | - Yuri Kogan
- Institute for Medical BioMathematics, Bene Ataroth, Israel
| | - Zvia Agur
- Institute for Medical BioMathematics, Bene Ataroth, Israel
- * E-mail:
| |
Collapse
|
78
|
Kirouac DC, Schaefer G, Chan J, Merchant M, Orr C, Huang SMA, Moffat J, Liu L, Gadkar K, Ramanujan S. Clinical responses to ERK inhibition in BRAFV600E-mutant colorectal cancer predicted using a computational model. NPJ Syst Biol Appl 2017; 3:14. [PMID: 28649441 PMCID: PMC5460205 DOI: 10.1038/s41540-017-0016-1] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2017] [Revised: 04/18/2017] [Accepted: 05/04/2017] [Indexed: 12/11/2022] Open
Abstract
Approximately 10% of colorectal cancers harbor BRAFV600E mutations, which constitutively activate the MAPK signaling pathway. We sought to determine whether ERK inhibitor (GDC-0994)-containing regimens may be of clinical benefit to these patients based on data from in vitro (cell line) and in vivo (cell- and patient-derived xenograft) studies of cetuximab (EGFR), vemurafenib (BRAF), cobimetinib (MEK), and GDC-0994 (ERK) combinations. Preclinical data was used to develop a mechanism-based computational model linking cell surface receptor (EGFR) activation, the MAPK signaling pathway, and tumor growth. Clinical predictions of anti-tumor activity were enabled by the use of tumor response data from three Phase 1 clinical trials testing combinations of EGFR, BRAF, and MEK inhibitors. Simulated responses to GDC-0994 monotherapy (overall response rate = 17%) accurately predicted results from a Phase 1 clinical trial regarding the number of responding patients (2/18) and the distribution of tumor size changes ("waterfall plot"). Prospective simulations were then used to evaluate potential drug combinations and predictive biomarkers for increasing responsiveness to MEK/ERK inhibitors in these patients.
Collapse
Affiliation(s)
- Daniel C. Kirouac
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Gabriele Schaefer
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Jocelyn Chan
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Mark Merchant
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Christine Orr
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Shih-Min A. Huang
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - John Moffat
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Lichuan Liu
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Kapil Gadkar
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| | - Saroja Ramanujan
- Genentech Research & Early Development, 1 DNA Way, South San Francisco, CA 94080 USA
| |
Collapse
|
79
|
Masunaga H, Sugimoto Y, Magi S, Itasaki R, Okada-Hatakeyama M, Kurata H. Robustness analysis of the detailed kinetic model of an ErbB signaling network by using dynamic sensitivity. PLoS One 2017; 12:e0178250. [PMID: 28542548 PMCID: PMC5443533 DOI: 10.1371/journal.pone.0178250] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2016] [Accepted: 05/02/2017] [Indexed: 11/19/2022] Open
Abstract
The ErbB receptor signaling pathway plays an important role in the regulation of cellular proliferation, survival and differentiation, and dysregulation of the pathway is linked to various types of human cancer. Mathematical models have been developed as a practical complementary approach to deciphering the complexity of ErbB receptor signaling and elucidating how the pathways discriminate between ligands to induce different cell fates. In this study, we developed a simulator to accurately calculate the dynamic sensitivity of extracellular-signal-regulated kinase (ERK) activity (ERK*) and Akt activity (Akt*), downstream of the ErbB receptors stimulated with epidermal growth factor (EGF) and heregulin (HRG). To demonstrate the feasibility of this simulator, we estimated how the reactions critically responsible for ERK* and Akt* change with time and in response to different doses of EGF and HRG, and predicted that only a small number of reactions determine ERK* and Akt*. ERK* increased steeply with increasing HRG dose until saturation, while showing a gently rising response to EGF. Akt* had a gradual wide-range response to HRG and a blunt response to EGF. Akt* was sensitive to perturbations of intracellular kinetics, while ERK* was more robust due to multiple, negative feedback loops. Overall, the simulator predicted reactions that were critically responsible for ERK* and Akt* in response to the dose of EGF and HRG, illustrated the response characteristics of ERK* and Akt*, and estimated mechanisms for generating robustness in the ErbB signaling network.
Collapse
Affiliation(s)
- Hiroyuki Masunaga
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka, Japan
| | - Yurie Sugimoto
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka, Japan
| | - Shigeyuki Magi
- Laboratory for Integrated Cellular Systems, RIKEN Center for Integrative Medical Sciences (IMS), Tsurumi-ku, Yokohama, Kanagawa, Japan
| | - Ryunosuke Itasaki
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka, Japan
| | - Mariko Okada-Hatakeyama
- Laboratory for Integrated Cellular Systems, RIKEN Center for Integrative Medical Sciences (IMS), Tsurumi-ku, Yokohama, Kanagawa, Japan
- Laboratory of Cell Systems, Institute for Protein Research, Osaka University, Suita-shi, Osaka, Japan
| | - Hiroyuki Kurata
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka, Japan
- Biomedical Informatics R&D Center, Kyushu Institute of Technology, Iizuka, Fukuoka, Japan
- * E-mail:
| |
Collapse
|
80
|
Bates AM, Lanzel EA, Qian F, Abbasi T, Vali S, Brogden KA. Cell genomics and immunosuppressive biomarker expression influence PD-L1 immunotherapy treatment responses in HNSCC-a computational study. Oral Surg Oral Med Oral Pathol Oral Radiol 2017; 124:157-164. [PMID: 28756882 DOI: 10.1016/j.oooo.2017.05.474] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2017] [Revised: 05/09/2017] [Accepted: 05/09/2017] [Indexed: 12/31/2022]
Abstract
OBJECTIVES Programmed death-ligand 1 (PD-L1) expression is correlated with objective response rates to PD-1 and PD-L1 immunotherapies. However, both immunotherapies have only demonstrated 12%-24.8% objective response rates in patients with head and neck squamous cell carcinoma (HNSCC), demonstrating a need for a more accurate method to identify those who will respond before their therapy. Immunohistochemistry to detect PD-L1 reactivity in tumors can be challenging, and additional methods are needed to predict and confirm PD-L1 expression. Here, we hypothesized that HNSCC tumor cell genomics influences cell signaling and downstream effects on immunosuppressive biomarkers and that these profiles can predict patient clinical responses. STUDY DESIGN We identified deleterious gene mutations in SCC4, SCC15, and SCC25 and created cell line-specific predictive computational simulation models. The expression of 24 immunosuppressive biomarkers were then predicted and used to sort cell lines into those that would respond to PD-L1 immunotherapy and those that would not. RESULTS SCC15 and SCC25 were identified as cell lines that would respond to PD-L1 immunotherapy treatment and SCC4 was identified as a cell line that would not likely respond to PD-L1 immunotherapy treatment. CONCLUSIONS This approach, when applied to HNSCC cells, has the ability to predict PD-L1 expression and predict PD-1- or PD-L1-targeted treatment responses in these patients.
Collapse
Affiliation(s)
- Amber M Bates
- Iowa Institute for Oral Health Research, College of Dentistry, University of Iowa, Iowa City, IA, USA
| | - Emily A Lanzel
- Department of Oral Pathology, Radiology and Medicine, College of Dentistry, University of Iowa, Iowa City, IA, USA
| | - Fang Qian
- Iowa Institute for Oral Health Research, College of Dentistry, University of Iowa, Iowa City, IA, USA; Division of Biostatistics and Research Design, College of Dentistry, University of Iowa, Iowa City, IA, USA
| | | | | | - Kim A Brogden
- Iowa Institute for Oral Health Research, College of Dentistry, University of Iowa, Iowa City, IA, USA.
| |
Collapse
|
81
|
Varga A, Ehrenreiter K, Aschenbrenner B, Kocieniewski P, Kochanczyk M, Lipniacki T, Baccarini M. RAF1/BRAF dimerization integrates the signal from RAS to ERK and ROKα. Sci Signal 2017; 10:eaai8482. [PMID: 28270557 DOI: 10.1126/scisignal.aai8482] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Downstream of growth factor receptors and of the guanine triphosphatase (GTPase) RAS, heterodimers of the serine/threonine kinases BRAF and RAF1 are critical upstream kinases and activators of the mitogen-activated protein kinase (MAPK) module containing the mitogen-activated and extracellular signal-regulated kinase kinase (MEK) and their targets, the extracellular signal-regulated kinase (ERK) family. Either direct or scaffold protein-mediated interactions among the components of the ERK module (the MAPKKKs BRAF and RAF1, MEK, and ERK) facilitate signal transmission. RAF1 also has essential functions in the control of tumorigenesis and migration that are mediated through its interaction with the kinase ROKα, an effector of the GTPase RHO and regulator of cytoskeletal rearrangements. We combined mutational and kinetic analysis with mathematical modeling to show that the interaction of RAF1 with ROKα is coordinated with the role of RAF1 in the ERK pathway. We found that the phosphorylated form of RAF1 that interacted with and inhibited ROKα was generated during the interaction of RAF1 with the ERK module. This mechanism adds plasticity to the ERK pathway, enabling signal diversification at the level of both ERK and RAF. Furthermore, by connecting ERK activation with the regulation of ROKα and cytoskeletal rearrangements by RAF1, this mechanism has the potential to precisely coordinate the proper timing of proliferation with changes in cell shape, adhesion, or motility.
Collapse
Affiliation(s)
- Andrea Varga
- Department of Microbiology, Immunology and Genetics, Max F. Perutz Laboratories, University of Vienna, Vienna, Austria
| | - Karin Ehrenreiter
- Department of Microbiology, Immunology and Genetics, Max F. Perutz Laboratories, University of Vienna, Vienna, Austria
| | - Bertram Aschenbrenner
- Department of Microbiology, Immunology and Genetics, Max F. Perutz Laboratories, University of Vienna, Vienna, Austria
| | - Pawel Kocieniewski
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Marek Kochanczyk
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Tomasz Lipniacki
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Manuela Baccarini
- Department of Microbiology, Immunology and Genetics, Max F. Perutz Laboratories, University of Vienna, Vienna, Austria.
| |
Collapse
|
82
|
Das AA, Ajayakumar Darsana T, Jacob E. Agent-based re-engineering of ErbB signaling: a modeling pipeline for integrative systems biology. Bioinformatics 2017; 33:726-732. [PMID: 27998938 DOI: 10.1093/bioinformatics/btw709] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2016] [Accepted: 11/08/2016] [Indexed: 11/14/2022] Open
Abstract
Motivation Experiments in systems biology are generally supported by a computational model which quantitatively estimates the parameters of the system by finding the best fit to the experiment. Mathematical models have proved to be successful in reverse engineering the system. The data generated is interpreted to understand the dynamics of the underlying phenomena. The question we have sought to answer is that - is it possible to use an agent-based approach to re-engineer a biological process, making use of the available knowledge from experimental and modelling efforts? Can the bottom-up approach benefit from the top-down exercise so as to create an integrated modelling formalism for systems biology? We propose a modelling pipeline that learns from the data given by reverse engineering, and uses it for re-engineering the system, to carry out in-silico experiments. Results A mathematical model that quantitatively predicts co-expression of EGFR-HER2 receptors in activation and trafficking has been taken for this study. The pipeline architecture takes cues from the population model that gives the rates of biochemical reactions, to formulate knowledge-based rules for the particle model. Agent-based simulations using these rules, support the existing facts on EGFR-HER2 dynamics. We conclude that, re-engineering models, built using the results of reverse engineering, opens up the possibility of harnessing the power pack of data which now lies scattered in literature. Virtual experiments could then become more realistic when empowered with the findings of empirical cell biology and modelling studies. Availability and Implementation Implemented on the Agent Modelling Framework developed in-house. C ++ code templates available in Supplementary material . Contact liz.csir@gmail.com. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
|
83
|
Computational investigation of sphingosine kinase 1 (SphK1) and calcium dependent ERK1/2 activation downstream of VEGFR2 in endothelial cells. PLoS Comput Biol 2017; 13:e1005332. [PMID: 28178265 PMCID: PMC5298229 DOI: 10.1371/journal.pcbi.1005332] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2016] [Accepted: 12/23/2016] [Indexed: 01/14/2023] Open
Abstract
Vascular endothelial growth factor (VEGF) is a powerful regulator of neovascularization. VEGF binding to its cognate receptor, VEGFR2, activates a number of signaling pathways including ERK1/2. Activation of ERK1/2 is experimentally shown to involve sphingosine kinase 1 (SphK1) activation and its calcium-dependent translocation downstream of ERK1/2. Here we construct a rule-based computational model of signaling downstream of VEGFR2, by including SphK1 and calcium positive feedback mechanisms, and investigate their consequences on ERK1/2 activation. The model predicts the existence of VEGF threshold in ERK1/2 activation that can be continuously tuned by cellular concentrations of SphK1 and sphingosine 1 phosphate (S1P). The computer model also predicts powerful effects of perturbations in plasma and ER calcium pump rates and the current through the CRAC channels on ERK1/2 activation dynamics, highlighting the critical role of intracellular calcium in shaping the pERK1/2 signal. The model is then utilized to simulate anti-angiogenic therapeutic interventions targeting VEGFR2-ERK1/2 axis. Simulations indicate that monotherapies that exclusively target VEGFR2 phosphorylation, VEGF, or VEGFR2 are ineffective in shutting down signaling to ERK1/2. By simulating therapeutic strategies that target multiple nodes of the pathway such as Raf and SphK1, we conclude that combination therapy should be much more effective in blocking VEGF signaling to EKR1/2. The model has important implications for interventions that target signaling pathways in angiogenesis relevant to cancer, vascular diseases, and wound healing. Vascular endothelial growth factor (VEGF) signaling is a potent regulator of angiogenesis, the growth and development of new vessels out of a preexisting vascular network. Angiogenesis requires enhanced survival, proliferation, and motility of the vascular endothelial cells. Crucial signaling endpoints in VEGF-mediated angiogenic response include elevation in intracellular calcium and the activation of the proteins ERK1 and 2 (ERK1/2). In this study, we have developed a novel computer model for the activation of ERK1/2 and calcium downstream of VEGF receptor type 2 (VEGFR2). Our model is the first of its kind to incorporate and investigate the consequences of calcium elevation and the role of a cellular lipid modifier known as sphingosine kinase 1 (SphK1). We also utilize the model to simulate therapeutic strategies targeting VEGF signaling to ERK1/2 indicating inefficiency of single therapies known as tyrosine kinase inhibitors (TKI) that target receptor phosphorylation. Computer simulations indicate that combination therapy is essential for effective blockade of this important pathway. Our results have important implications for human diseases such as cancer where plethora of anti-VEGF therapies are currently employed. Overall, our computer model sheds new light on a complex feedback involving SphK1 and calcium that radically alters the response of cells to VEGF.
Collapse
|
84
|
Huang L, Jiang Y, Chen Y. Predicting Drug Combination Index and Simulating the Network-Regulation Dynamics by Mathematical Modeling of Drug-Targeted EGFR-ERK Signaling Pathway. Sci Rep 2017; 7:40752. [PMID: 28102344 PMCID: PMC5244366 DOI: 10.1038/srep40752] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2016] [Accepted: 12/06/2016] [Indexed: 02/05/2023] Open
Abstract
Synergistic drug combinations enable enhanced therapeutics. Their discovery typically involves the measurement and assessment of drug combination index (CI), which can be facilitated by the development and applications of in-silico CI predictive tools. In this work, we developed and tested the ability of a mathematical model of drug-targeted EGFR-ERK pathway in predicting CIs and in analyzing multiple synergistic drug combinations against observations. Our mathematical model was validated against the literature reported signaling, drug response dynamics, and EGFR-MEK drug combination effect. The predicted CIs and combination therapeutic effects of the EGFR-BRaf, BRaf-MEK, FTI-MEK, and FTI-BRaf inhibitor combinations showed consistent synergism. Our results suggest that existing pathway models may be potentially extended for developing drug-targeted pathway models to predict drug combination CI values, isobolograms, and drug-response surfaces as well as to analyze the dynamics of individual and combinations of drugs. With our model, the efficacy of potential drug combinations can be predicted. Our method complements the developed in-silico methods (e.g. the chemogenomic profile and the statistically-inferenced network models) by predicting drug combination effects from the perspectives of pathway dynamics using experimental or validated molecular kinetic constants, thereby facilitating the collective prediction of drug combination effects in diverse ranges of disease systems.
Collapse
Affiliation(s)
- Lu Huang
- The Ministry-Province Jointly Constructed Base for State Key Lab and Shenzhen Technology and Engineering Lab for Personalized Cancer Diagnostics and Therapeutics Tsinghua University Shenzhen Graduate School, and Shenzhen Kivita Innovative Drug Discovery Institute, Shenzhen, 518055, P.R. China
- Institute of Molecular Biology (IMB), Ackermannweg 4, 55128 Mainz, Germany
- Department of Pharmacy, and Center for Computational Science and Engineering, National University of Singapore, Blk S16, Level 8, 3 Science Drive 2, 117543 Singapore
| | - Yuyang Jiang
- The Ministry-Province Jointly Constructed Base for State Key Lab and Shenzhen Technology and Engineering Lab for Personalized Cancer Diagnostics and Therapeutics Tsinghua University Shenzhen Graduate School, and Shenzhen Kivita Innovative Drug Discovery Institute, Shenzhen, 518055, P.R. China
| | - Yuzong Chen
- Department of Pharmacy, and Center for Computational Science and Engineering, National University of Singapore, Blk S16, Level 8, 3 Science Drive 2, 117543 Singapore
- State Key Laboratory of Biotherapy, West China Hospital, West China School of Medicine, Sichuan University, Chengdu, China
| |
Collapse
|
85
|
Modeling of Receptor Tyrosine Kinase Signaling: Computational and Experimental Protocols. Methods Mol Biol 2017; 1636:417-453. [PMID: 28730495 DOI: 10.1007/978-1-4939-7154-1_27] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The advent of systems biology has convincingly demonstrated that the integration of experiments and dynamic modelling is a powerful approach to understand the cellular network biology. Here we present experimental and computational protocols that are necessary for applying this integrative approach to the quantitative studies of receptor tyrosine kinase (RTK) signaling networks. Signaling by RTKs controls multiple cellular processes, including the regulation of cell survival, motility, proliferation, differentiation, glucose metabolism, and apoptosis. We describe methods of model building and training on experimentally obtained quantitative datasets, as well as experimental methods of obtaining quantitative dose-response and temporal dependencies of protein phosphorylation and activities. The presented methods make possible (1) both the fine-grained modeling of complex signaling dynamics and identification of salient, course-grained network structures (such as feedback loops) that bring about intricate dynamics, and (2) experimental validation of dynamic models.
Collapse
|
86
|
Aliper AM, Korzinkin MB, Kuzmina NB, Zenin AA, Venkova LS, Smirnov PY, Zhavoronkov AA, Buzdin AA, Borisov NM. Mathematical Justification of Expression-Based Pathway Activation Scoring (PAS). Methods Mol Biol 2017; 1613:31-51. [PMID: 28849557 DOI: 10.1007/978-1-4939-7027-8_3] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Although modeling of activation kinetics for various cell signaling pathways has reached a high grade of sophistication and thoroughness, most such kinetic models still remain of rather limited practical value for biomedicine. Nevertheless, recent advancements have been made in application of signaling pathway science for real needs of prescription of the most effective drugs for individual patients. The methods for such prescription evaluate the degree of pathological changes in the signaling machinery based on two types of data: first, on the results of high-throughput gene expression profiling, and second, on the molecular pathway graphs that reflect interactions between the pathway members. For example, our algorithm OncoFinder evaluates the activation of molecular pathways on the basis of gene/protein expression data in the objects of the interest.Yet, the question of assessment of the relative importance for each gene product in a molecular pathway remains unclear unless one call for the methods of parameter sensitivity /stiffness analysis in the interactomic kinetic models of signaling pathway activation in terms of total concentrations of each gene product.Here we show two principal points: 1. First, the importance coefficients for each gene in pathways that were obtained using the extremely time- and labor-consuming stiffness analysis of full-scaled kinetic models generally differ from much easier-to-calculate expression-based pathway activation score (PAS) not more than by 30%, so the concept of PAS is kinetically justified. 2. Second, the use of pathway-based approach instead of distinct gene analysis, due to the law of large numbers, allows restoring the correlation between the similar samples that were examined using different transcriptome investigation techniques.
Collapse
Affiliation(s)
- Alexander M Aliper
- Drug Research and Design Department, Pathway Pharmaceuticals, Wan Chai, Hong Kong, Hong Kong SAR
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
| | - Michael B Korzinkin
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
| | - Natalia B Kuzmina
- Laboratory of Systems Biology, A.I. Burnazyan Federal Medical Biophysical Center, Moscow, 123182, Russia
| | - Alexander A Zenin
- Laboratory of Systems Biology, A.I. Burnazyan Federal Medical Biophysical Center, Moscow, 123182, Russia
| | - Larisa S Venkova
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
| | - Philip Yu Smirnov
- Laboratory of Systems Biology, A.I. Burnazyan Federal Medical Biophysical Center, Moscow, 123182, Russia
| | - Alex A Zhavoronkov
- Drug Research and Design Department, Pathway Pharmaceuticals, Wan Chai, Hong Kong, Hong Kong SAR
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
| | - Anton A Buzdin
- Drug Research and Design Department, Pathway Pharmaceuticals, Wan Chai, Hong Kong, Hong Kong SAR
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia
- Group for Genomic Regulation of Cell Signaling Systems, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia
- National Research Centre "Kurchatov Institute", Centre for Convergence of Nano-, Bio-, Information and Cognitive Sciences and Technologies, Moscow, Russia
| | - Nikolay M Borisov
- Department of Personalized Medicine, First Oncology Research and Advisory Center, Moscow, Russia.
- Laboratory of Bioinformatics, D. Rogachev Federal Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia.
- National Research Centre "Kurchatov Institute", Centre for Convergence of Nano-, Bio-, Information and Cognitive Sciences and Technologies, Moscow, Russia.
| |
Collapse
|
87
|
Shin SY, Nguyen LK. Dissecting Cell-Fate Determination Through Integrated Mathematical Modeling of the ERK/MAPK Signaling Pathway. Methods Mol Biol 2017; 1487:409-432. [PMID: 27924583 DOI: 10.1007/978-1-4939-6424-6_29] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The past three decades have witnessed an enormous progress in the elucidation of the ERK/MAPK signaling pathway and its involvement in various cellular processes. Because of its importance and complex wiring, the ERK pathway has been an intensive subject for mathematical modeling, which facilitates the unraveling of key dynamic properties and behaviors of the pathway. Recently, however, it became evident that the pathway does not act in isolation but closely interacts with many other pathways to coordinate various cellular outcomes under different pathophysiological contexts. This has led to an increasing number of integrated, large-scale models that link the ERK pathway to other functionally important pathways. In this chapter, we first discuss the essential steps in model development and notable models of the ERK pathway. We then use three examples of integrated, multipathway models to investigate how crosstalk of ERK signaling with other pathways regulates cell-fate decision-making in various physiological and disease contexts. Specifically, we focus on ERK interactions with the phosphoinositide-3 kinase (PI3K), c-Jun N-terminal kinase (JNK), and β-adrenergic receptor (β-AR) signaling pathways. We conclude that integrated modeling in combination with wet-lab experimentation have been and will be instrumental in gaining an in-depth understanding of ERK signaling in multiple biological contexts.
Collapse
Affiliation(s)
- Sung-Young Shin
- Department of Biochemistry and Molecular Biology, School of Biomedical Sciences, Monash University, Clayton, VIC, 3800, Australia.,Biomedicine Discovery Institute, Monash University, Clayton, VIC, 3800, Australia
| | - Lan K Nguyen
- Department of Biochemistry and Molecular Biology, School of Biomedical Sciences, Monash University, Clayton, VIC, 3800, Australia. .,Biomedicine Discovery Institute, Monash University, Clayton, VIC, 3800, Australia.
| |
Collapse
|
88
|
Vanhaelen Q, Aliper AM, Zhavoronkov A. A comparative review of computational methods for pathway perturbation analysis: dynamical and topological perspectives. MOLECULAR BIOSYSTEMS 2017; 13:1692-1704. [DOI: 10.1039/c7mb00170c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Stem cells offer great promise within the field of regenerative medicine but despite encouraging results, the large scale use of stem cells for therapeutic applications still faces challenges when it comes to controlling signaling pathway responses with respect to environmental perturbations.
Collapse
Affiliation(s)
- Q. Vanhaelen
- Insilico Medicine Inc
- Johns Hopkins University
- ETC
- USA
| | - A. M. Aliper
- Insilico Medicine Inc
- Johns Hopkins University
- ETC
- USA
| | | |
Collapse
|
89
|
Halasz M, Kholodenko BN, Kolch W, Santra T. Integrating network reconstruction with mechanistic modeling to predict cancer therapies. Sci Signal 2016; 9:ra114. [PMID: 27879396 DOI: 10.1126/scisignal.aae0535] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Signal transduction networks are often rewired in cancer cells. Identifying these alterations will enable more effective cancer treatment. We developed a computational framework that can identify, reconstruct, and mechanistically model these rewired networks from noisy and incomplete perturbation response data and then predict potential targets for intervention. As a proof of principle, we analyzed a perturbation data set targeting epidermal growth factor receptor (EGFR) and insulin-like growth factor 1 receptor (IGF1R) pathways in a panel of colorectal cancer cells. Our computational approach predicted cell line-specific network rewiring. In particular, feedback inhibition of insulin receptor substrate 1 (IRS1) by the kinase p70S6K was predicted to confer resistance to EGFR inhibition, suggesting that disrupting this feedback may restore sensitivity to EGFR inhibitors in colorectal cancer cells. We experimentally validated this prediction with colorectal cancer cell lines in culture and in a zebrafish (Danio rerio) xenograft model.
Collapse
Affiliation(s)
- Melinda Halasz
- Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Ireland. .,School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
| | - Boris N Kholodenko
- Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Ireland.,School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland.,Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin 4, Ireland
| | - Walter Kolch
- Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Ireland. .,School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland.,Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin 4, Ireland
| | - Tapesh Santra
- Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Ireland.
| |
Collapse
|
90
|
Unveiling Hidden Dynamics of Hippo Signalling: A Systems Analysis. Genes (Basel) 2016; 7:genes7080044. [PMID: 27527217 PMCID: PMC4999832 DOI: 10.3390/genes7080044] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2016] [Revised: 07/27/2016] [Accepted: 07/29/2016] [Indexed: 12/29/2022] Open
Abstract
The Hippo signalling pathway has recently emerged as an important regulator of cell apoptosis and proliferation with significant implications in human diseases. In mammals, the pathway contains the core kinases MST1/2, which phosphorylate and activate LATS1/2 kinases. The pro-apoptotic function of the MST/LATS signalling axis was previously linked to the Akt and ERK MAPK pathways, demonstrating that the Hippo pathway does not act alone but crosstalks with other signalling pathways to coordinate network dynamics and cellular outcomes. These crosstalks were characterised by a multitude of complex regulatory mechanisms involving competitive protein-protein interactions and phosphorylation mediated feedback loops. However, how these different mechanisms interplay in different cellular contexts to drive the context-specific network dynamics of Hippo-ERK signalling remains elusive. Using mathematical modelling and computational analysis, we uncovered that the Hippo-ERK network can generate highly diverse dynamical profiles that can be clustered into distinct dose-response patterns. For each pattern, we offered mechanistic explanation that defines when and how the observed phenomenon can arise. We demonstrated that Akt displays opposing, dose-dependent functions towards ERK, which are mediated by the balance between the Raf-1/MST2 protein interaction module and the LATS1 mediated feedback regulation. Moreover, Ras displays a multi-functional role and drives biphasic responses of both MST2 and ERK activities; which are critically governed by the competitive protein interaction between MST2 and Raf-1. Our study represents the first in-depth and systematic analysis of the Hippo-ERK network dynamics and provides a concrete foundation for future studies.
Collapse
|
91
|
Huang H, Haar Petersen M, Ibañez-Vea M, Lassen PS, Larsen MR, Palmisano G. Simultaneous Enrichment of Cysteine-containing Peptides and Phosphopeptides Using a Cysteine-specific Phosphonate Adaptable Tag (CysPAT) in Combination with titanium dioxide (TiO2) Chromatography. Mol Cell Proteomics 2016; 15:3282-3296. [PMID: 27281782 PMCID: PMC5054350 DOI: 10.1074/mcp.m115.054551] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2015] [Indexed: 12/13/2022] Open
Abstract
Cysteine is a rare and conserved amino acid involved in most cellular functions. The thiol group of cysteine can be subjected to diverse oxidative modifications that regulate many physio-pathological states. In the present work, a Cysteine-specific Phosphonate Adaptable Tag (CysPAT) was synthesized to selectively label cysteine-containing peptides (Cys peptides) followed by their enrichment with titanium dioxide (TiO2) and subsequent mass spectrometric analysis. The CysPAT strategy was developed using a synthetic peptide, a standard protein and subsequently the strategy was applied to protein lysates from Hela cells, achieving high specificity and enrichment efficiency. In particular, for Cys proteome analysis, the method led to the identification of 7509 unique Cys peptides from 500 μg of HeLa cell lysate starting material. Furthermore, the method was developed to simultaneously enrich Cys peptides and phosphorylated peptides. This strategy was applied to SILAC labeled Hela cells subjected to 5 min epidermal growth factor (EGF) stimulation. In total, 10440 unique reversibly modified Cys peptides (3855 proteins) and 7339 unique phosphopeptides (2234 proteins) were simultaneously identified from 250 μg starting material. Significant regulation was observed in both phosphorylation and reversible Cys modification of proteins involved in EGFR signaling. Our data indicates that EGF stimulation can activate the well-known phosphorylation of EGFR and downstream signaling molecules, such as mitogen-activated protein kinases (MAPK1 and MAPK3), however, it also leads to substantial modulation of reversible cysteine modifications in numerous proteins. Several protein tyrosine phosphatases (PTPs) showed a reduction of the catalytic Cys site in the conserved putative phosphatase HC(X)5R motif indicating an activation and subsequent de-phosphorylation of proteins involved in the EGF signaling pathway. Overall, the CysPAT strategy is a straight forward, easy and promising method for studying redox proteomics and the simultaneous enrichment strategy offers an excellent solution for characterization of cross-talk between phosphorylation and redox induced reversible cysteine modifications.
Collapse
Affiliation(s)
- Honggang Huang
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark; §The Danish Diabetes Academy, Odense, Denmark
| | - Martin Haar Petersen
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark; ¶Institute of Molecular Medicine, Cancer & Inflammation Research, University of Southern Denmark
| | - Maria Ibañez-Vea
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark
| | - Pernille S Lassen
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark
| | - Martin R Larsen
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark
| | - Giuseppe Palmisano
- From the ‡Department of Biochemistry and Molecular Biology, University of Southern Denmark; ‖Department of Parasitology, ICB, University of São Paulo, Brazil
| |
Collapse
|
92
|
Morrison ES, Badyaev AV. The Landscape of Evolution: Reconciling Structural and Dynamic Properties of Metabolic Networks in Adaptive Diversifications. Integr Comp Biol 2016; 56:235-46. [PMID: 27252203 DOI: 10.1093/icb/icw026] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
The network of the interactions among genes, proteins, and metabolites delineates a range of potential phenotypic diversifications in a lineage, and realized phenotypic changes are the result of differences in the dynamics of the expression of the elements and interactions in this deterministic network. Regulatory mechanisms, such as hormones, mediate the relationship between the structural and dynamic properties of networks by determining how and when the elements are expressed and form a functional unit or state. Changes in regulatory mechanisms lead to variable expression of functional states of a network within and among generations. Functional properties of network elements, and the magnitude and direction of evolutionary change they determine, depend on their location within a network. Here, we examine the relationship between network structure and the dynamic mechanisms that regulate flux through a metabolic network. We review the mechanisms that control metabolic flux in enzymatic reactions and examine structural properties of the network locations that are targets of flux control. We aim to establish a predictive framework to test the contributions of structural and dynamic properties of deterministic networks to evolutionary diversifications.
Collapse
Affiliation(s)
- Erin S Morrison
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721-0001, USA
| | - Alexander V Badyaev
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721-0001, USA
| |
Collapse
|
93
|
Hathcock D, Sheehy J, Weisenberger C, Ilker E, Hinczewski M. Noise Filtering and Prediction in Biological Signaling Networks. ACTA ACUST UNITED AC 2016. [DOI: 10.1109/tmbmc.2016.2633269] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
94
|
Abstract
To understand the behaviour of complex systems, it is often necessary to use models that describe the dynamics of subnetworks. It has previously been established using projection methods that such subnetwork dynamics generically involves memory of the past and that the memory functions can be calculated explicitly for biochemical reaction networks made up of unary and binary reactions. However, many established network models involve also Michaelis-Menten kinetics, to describe, e.g., enzymatic reactions. We show that the projection approach to subnetwork dynamics can be extended to such networks, thus significantly broadening its range of applicability. To derive the extension, we construct a larger network that represents enzymes and enzyme complexes explicitly, obtain the projected equations, and finally take the limit of fast enzyme reactions that gives back Michaelis-Menten kinetics. The crucial point is that this limit can be taken in closed form. The outcome is a simple procedure that allows one to obtain a description of subnetwork dynamics, including memory functions, starting directly from any given network of unary, binary, and Michaelis-Menten reactions. Numerical tests show that this closed form enzyme elimination gives a much more accurate description of the subnetwork dynamics than the simpler method that represents enzymes explicitly and is also more efficient computationally.
Collapse
Affiliation(s)
- Katy J Rubin
- Department of Mathematics, King's College London, Strand, London WC2R 2LS, United Kingdom
| | - Peter Sollich
- Department of Mathematics, King's College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
95
|
Arampatzis G, Katsoulakis MA, Rey-Bellet L. Efficient estimators for likelihood ratio sensitivity indices of complex stochastic dynamics. J Chem Phys 2016; 144:104107. [DOI: 10.1063/1.4943388] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Georgios Arampatzis
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Markos A. Katsoulakis
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Luc Rey-Bellet
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| |
Collapse
|
96
|
Maya-Bernal JL, Ramírez-Santiago G. Spatio-temporal dynamics of a cell signal pathway with negative feedbacks: the MAPK/ERK pathway. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2016; 39:28. [PMID: 26987732 DOI: 10.1140/epje/i2016-16028-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2015] [Accepted: 01/25/2016] [Indexed: 06/05/2023]
Abstract
We studied the spatio-temporal dynamics of a cell signal cascade with negative feedback that quantitatively emulates the regulative process that occurs in the Mitogen Activated Protein Kinase/Extracellular Regulated Kinase (MAPK/ERK) pathway. The model consists of a set of six coupled reaction-diffusion equations that describes the dynamics of the six-module pathway. In the basic module the active form of the protein transmits the signal to the next pathway’s module. As suggested by experiments, the model considers that the fifth module's kinase down-regulates the first and third modules. The feedback parameter is defined as, μ(r)( j)= k(kin)5/k(kin)(j), (j = 1, 3). We analysed the pathway's dynamics for μ(r)( j) = 0.10, 1.0, and 10 in the kinetic regimes: i) saturation of both kinases and phosphatases, ii) saturation of the phosphatases and iii) saturation of the kinases. For a regulated pathway the Total Activated Protein Profiles (TAPPs) as a function of time develop a maximum during the transient stage in the three kinetic regimes. These maxima become higher and their positions shift to longer times downstream. This scenario also applies to the TAPP's regulatory kinase that sums up its inhibitory action to that of the phosphatases leading to a maximum. Nevertheless, when μ(r)(j)= 1.0 , the TAPPs develop two maxima, with the second maximum being almost imperceptible. These results are in qualitative agreement with experimental data obtained from NIH 3T3 mouse fibroblasts. In addition, analyses of the stationary states as a function of position indicate that in the kinetic regime i) which is of physiological interest, signal transduction occurs with a relatively large propagation length for the three values of the regulative parameter. However, for μ(r)(j)= 0.10 , the sixth module concentration profile is transmitted with approximately 45% of its full value. The results obtained for μ(r)(j) = 10 , indicate that the first five concentration profiles are small with a short propagation length; nonetheless, the last concentration profile, c6, attains more than 90% of its full value with a relatively large propagation length as an indication of signal transduction. Signal transduction also occurred favourably in the kinetic regimes ii) and iii), but the signal was longer-ranged in the regime ii).
Collapse
Affiliation(s)
- José Luis Maya-Bernal
- Posgrado en Ciencias Bioquímicas, Universidad Nacional Autónoma de México, Coyoacán D.F., Mexico
| | - Guillermo Ramírez-Santiago
- Instituto de Matemáticas, Universidad Nacional Autónoma de México, C.P. 76230, Juriquilla Querétaro, Mexico.
| |
Collapse
|
97
|
Abstract
Despite extensive study of the EGF receptor (EGFR) signaling network, the immediate posttranslational changes that occur in response to growth factor stimulation remain poorly characterized; as a result, the biological mechanisms underlying signaling initiation remain obscured. To address this deficiency, we have used a mass spectrometry-based approach to measure system-wide phosphorylation changes throughout the network with 10-s resolution in the 80 s after stimulation in response to a range of eight growth factor concentrations. Significant changes were observed on proteins far downstream in the network as early as 10 s after stimulation, indicating a system capable of transmitting information quickly. Meanwhile, canonical members of the EGFR signaling network fall into clusters with distinct activation patterns. Src homology 2 domain containing transforming protein (Shc) and phosphoinositol 3-kinase (PI3K) phosphorylation levels increase rapidly, but equilibrate within 20 s, whereas proteins such as Grb2-associated binder-1 (Gab1) and SH2-containing tyrosine phosphatase (SHP2) show slower, sustained increases. Proximity ligation assays reveal that Shc and Gab1 phosphorylation patterns are representative of separate timescales for physical association with the receptor. Inhibition of phosphatases with vanadate reveals site-specific regulatory mechanisms and also uncovers primed activating components in the network, including Src family kinases, whose inhibition affects only a subset of proteins within the network. The results presented highlight the complexity of signaling initiation and provide a window into exploring mechanistic hypotheses about receptor tyrosine kinase (RTK) biology.
Collapse
|
98
|
Stites EC, Aziz M, Creamer MS, Von Hoff DD, Posner RG, Hlavacek WS. Use of mechanistic models to integrate and analyze multiple proteomic datasets. Biophys J 2016; 108:1819-1829. [PMID: 25863072 DOI: 10.1016/j.bpj.2015.02.030] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2014] [Revised: 02/18/2015] [Accepted: 02/24/2015] [Indexed: 11/30/2022] Open
Abstract
Proteins in cell signaling networks tend to interact promiscuously through low-affinity interactions. Consequently, evaluating the physiological importance of mapped interactions can be difficult. Attempts to do so have tended to focus on single, measurable physicochemical factors, such as affinity or abundance. For example, interaction importance has been assessed on the basis of the relative affinities of binding partners for a protein of interest, such as a receptor. However, multiple factors can be expected to simultaneously influence the recruitment of proteins to a receptor (and the potential of these proteins to contribute to receptor signaling), including affinity, abundance, and competition, which is a network property. Here, we demonstrate that measurements of protein copy numbers and binding affinities can be integrated within the framework of a mechanistic, computational model that accounts for mass action and competition. We use cell line-specific models to rank the relative importance of protein-protein interactions in the epidermal growth factor receptor (EGFR) signaling network for 11 different cell lines. Each model accounts for experimentally characterized interactions of six autophosphorylation sites in EGFR with proteins containing a Src homology 2 and/or phosphotyrosine-binding domain. We measure importance as the predicted maximal extent of recruitment of a protein to EGFR following ligand-stimulated activation of EGFR signaling. We find that interactions ranked highly by this metric include experimentally detected interactions. Proteins with high importance rank in multiple cell lines include proteins with recognized, well-characterized roles in EGFR signaling, such as GRB2 and SHC1, as well as a protein with a less well-defined role, YES1. Our results reveal potential cell line-specific differences in recruitment.
Collapse
Affiliation(s)
- Edward C Stites
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona; Department of Pathology & Immunology, Washington University School of Medicine, St. Louis, Missouri.
| | - Meraj Aziz
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona
| | - Matthew S Creamer
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona; Department of Molecular, Cellular, and Developmental Biology, Yale University, New Haven, Connecticut
| | - Daniel D Von Hoff
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona
| | - Richard G Posner
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona; Department of Biological Sciences, Northern Arizona University, Flagstaff, Arizona.
| | - William S Hlavacek
- Clinical Translational Research Division, Translational Genomics Research Institute, Phoenix, Arizona; Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico.
| |
Collapse
|
99
|
Chapter Six - The Ubiquitin Network in the Control of EGFR Endocytosis and Signaling. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2016; 141:225-76. [DOI: 10.1016/bs.pmbts.2016.03.002] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
|
100
|
Rowland MA, Harrison B, Deeds EJ. Phosphatase specificity and pathway insulation in signaling networks. Biophys J 2015; 108:986-996. [PMID: 25692603 DOI: 10.1016/j.bpj.2014.12.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2014] [Revised: 11/13/2014] [Accepted: 12/05/2014] [Indexed: 12/31/2022] Open
Abstract
Phosphatases play an important role in cellular signaling networks by regulating the phosphorylation state of proteins. Phosphatases are classically considered to be promiscuous, acting on tens to hundreds of different substrates. We recently demonstrated that a shared phosphatase can couple the responses of two proteins to incoming signals, even if those two substrates are from otherwise isolated areas of the network. This finding raises a potential paradox: if phosphatases are indeed highly promiscuous, how do cells insulate themselves against unwanted crosstalk? Here, we use mathematical models to explore three possible insulation mechanisms. One approach involves evolving phosphatase KM values that are large enough to prevent saturation by the phosphatase's substrates. Although this is an effective method for generating isolation, the phosphatase becomes a highly inefficient enzyme, which prevents the system from achieving switch-like responses and can result in slow response kinetics. We also explore the idea that substrate degradation can serve as an effective phosphatase. Assuming that degradation is unsaturatable, this mechanism could insulate substrates from crosstalk, but it would also preclude ultrasensitive responses and would require very high substrate turnover to achieve rapid dephosphorylation kinetics. Finally, we show that adaptor subunits, such as those found on phosphatases like PP2A, can provide effective insulation against phosphatase crosstalk, but only if their binding to substrates is uncoupled from their binding to the catalytic core. Analysis of the interaction network of PP2A's adaptor domains reveals that although its adaptors may isolate subsets of targets from one another, there is still a strong potential for phosphatase crosstalk within those subsets. Understanding how phosphatase crosstalk and the insulation mechanisms described here impact the function and evolution of signaling networks represents a major challenge for experimental and computational systems biology.
Collapse
Affiliation(s)
- Michael A Rowland
- Center for Computational Biology, University of Kansas, Lawrence, Kansas
| | - Brian Harrison
- Center for Computational Biology, University of Kansas, Lawrence, Kansas
| | - Eric J Deeds
- Center for Computational Biology, University of Kansas, Lawrence, Kansas; Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas; Santa Fe Institute, Santa Fe, New Mexico.
| |
Collapse
|