Antibody-mediated targeting of TNFR2 activates CD8+ T cells in mice and promotes antitumor immunity.

Tumor necrosis factor receptor 2 (TNFR2) is the alternate receptor for TNF and can mediate both pro- and anti-inflammatory activities of T cells. Although TNFR2 has been linked to enhanced suppressive activity of regulatory T cells (Tregs) in autoimmune diseases, the viability of TNFR2 as a target for cancer immunotherapy has been underappreciated. Here, we show that new murine monoclonal anti-TNFR2 antibodies yield robust antitumor activity and durable protective memory in multiple mouse cancer cell line models. The antibodies mediate potent Fc-dependent T cell costimulation and do not result in significant depletion of Tregs Corresponding human agonistic monoclonal anti-TNFR2 antibodies were identified and also had antitumor effects in humanized mouse models. Anti-TNFR2 antibodies could be developed as a novel treatment option for patients with cancer.

MM-131, a bispecific anti-Met/EpCAM mAb, inhibits HGF-dependent and HGF-independent Met signaling through concurrent binding to EpCAM.

Activation of the Met receptor tyrosine kinase, either by its ligand, hepatocyte growth factor (HGF), or via ligand-independent mechanisms, such as MET amplification or receptor overexpression, has been implicated in driving tumor proliferation, metastasis, and resistance to therapy. Clinical development of Met-targeted antibodies has been challenging, however, as bivalent antibodies exhibit agonistic properties, whereas monovalent antibodies lack potency and the capacity to down-regulate Met. Through computational modeling, we found that the potency of a monovalent antibody targeting Met could be dramatically improved by introducing a second binding site that recognizes an unrelated, highly expressed antigen on the tumor cell surface. Guided by this prediction, we engineered MM-131, a bispecific antibody that is monovalent for both Met and epithelial cell adhesion molecule (EpCAM). MM-131 is a purely antagonistic antibody that blocks ligand-dependent and ligand-independent Met signaling by inhibiting HGF binding to Met and inducing receptor down-regulation. Together, these mechanisms lead to inhibition of proliferation in Met-driven cancer cells, inhibition of HGF-mediated cancer cell migration, and inhibition of tumor growth in HGF-dependent and -independent mouse xenograft models. Consistent with its design, MM-131 is more potent in EpCAM-high cells than in EpCAM-low cells, and its potency decreases when EpCAM levels are reduced by RNAi. Evaluation of Met, EpCAM, and HGF levels in human tumor samples reveals that EpCAM is expressed at high levels in a wide range of Met-positive tumor types, suggesting a broad opportunity for clinical development of MM-131.

Analysis of Single-Cell RNA-Seq Identifies Cell-Cell Communication Associated with Tumor Characteristics.

Tumor ecosystems are composed of multiple cell types that communicate by ligand-receptor interactions. Targeting ligand-receptor interactions (for instance, with immune checkpoint inhibitors) can provide significant benefits for patients. However, our knowledge of which interactions occur in a tumor and how these interactions affect outcome is still limited. We present an approach to characterize communication by ligand-receptor interactions across all cell types in a microenvironment using single-cell RNA sequencing. We apply this approach to identify and compare the ligand-receptor interactions present in six syngeneic mouse tumor models. To identify interactions potentially associated with outcome, we regress interactions against phenotypic measurements of tumor growth rate. In addition, we quantify ligand-receptor interactions between T cell subsets and their relation to immune infiltration using a publicly available human melanoma dataset. Overall, this approach provides a tool for studying cell-cell interactions, their variability across tumors, and their relationship to outcome.

Modeling chemotherapy-induced stress to identify rational combination therapies in the DNA damage response pathway.

Cells respond to DNA damage by activating complex signaling networks that decide cell fate, promoting not only DNA damage repair and survival but also cell death. We have developed a multiscale computational model that quantitatively links chemotherapy-induced DNA damage response signaling to cell fate. The computational model was trained and calibrated on extensive data from U2OS osteosarcoma cells, including the cell cycle distribution of the initial cell population, signaling data measured by Western blotting, and cell fate data in response to chemotherapy treatment measured by time-lapse microscopy. The resulting mechanistic model predicted the cellular responses to chemotherapy alone and in combination with targeted inhibitors of the DNA damage response pathway, which we confirmed experimentally. Computational models such as the one presented here can be used to understand the molecular basis that defines the complex interplay between cell survival and cell death and to rationally identify chemotherapy-potentiating drug combinations.

Estimation of immune cell content in tumour tissue using single-cell RNA-seq data.

As interactions between the immune system and tumour cells are governed by a complex network of cell-cell interactions, knowing the specific immune cell composition of a solid tumour may be essential to predict a patient's response to immunotherapy. Here, we analyse in depth how to derive the cellular composition of a solid tumour from bulk gene expression data by mathematical deconvolution, using indication-specific and cell type-specific reference gene expression profiles (RGEPs) from tumour-derived single-cell RNA sequencing data. We demonstrate that tumour-derived RGEPs are essential for the successful deconvolution and that RGEPs from peripheral blood are insufficient. We distinguish nine major cell types, as well as three T cell subtypes. Using the tumour-derived RGEPs, we can estimate the content of many tumours associated immune and stromal cell types, their therapeutically relevant ratios, as well as an improved gene expression profile of the malignant cells.

Predicting ligand-dependent tumors from multi-dimensional signaling features.

Targeted therapies have shown significant patient benefit in about 5-10% of solid tumors that are addicted to a single oncogene. Here, we explore the idea of ligand addiction as a driver of tumor growth. High ligand levels in tumors have been shown to be associated with impaired patient survival, but targeted therapies have not yet shown great benefit in unselected patient populations. Using an approach of applying Bagged Decision Trees (BDT) to high-dimensional signaling features derived from a computational model, we can predict ligand dependent proliferation across a set of 58 cell lines. This mechanistic, multi-pathway model that features receptor heterodimerization, was trained on seven cancer cell lines and can predict signaling across two independent cell lines by adjusting only the receptor expression levels for each cell line. Interestingly, for patient samples the predicted tumor growth response correlates with high growth factor expression in the tumor microenvironment, which argues for a co-evolution of both factors in vivo.

IL-1ß-induced and p38MAPK-dependent activation of the mitogen-activated protein kinase-activated protein kinase 2 (MK2) in hepatocytes: Signal transduction with robust and concentration-independent signal amplification.

The IL-1ß induced activation of the p38MAPK/MAPK-activated protein kinase 2 (MK2) pathway in hepatocytes is important for control of the acute phase response and regulation of liver regeneration. Many aspects of the regulatory relevance of this pathway have been investigated in immune cells in the context of inflammation. However, very little is known about concentration-dependent activation kinetics and signal propagation in hepatocytes and the role of MK2. We established a mathematical model for IL-1ß-induced activation of the p38MAPK/MK2 pathway in hepatocytes that was calibrated to quantitative data on time- and IL-1ß concentration-dependent phosphorylation of p38MAPK and MK2 in primary mouse hepatocytes. This analysis showed that, in hepatocytes, signal transduction from IL-1ß via p38MAPK to MK2 is characterized by strong signal amplification. Quantification of p38MAPK and MK2 revealed that, in hepatocytes, at maximum, 11.3% of p38MAPK molecules and 36.5% of MK2 molecules are activated in response to IL-1ß. The mathematical model was experimentally validated by employing phosphatase inhibitors and the p38MAPK inhibitor SB203580. Model simulations predicted an IC50 of 1-1.2 µm for SB203580 in hepatocytes. In silico analyses and experimental validation demonstrated that the kinase activity of p38MAPK determines signal amplitude, whereas phosphatase activity affects both signal amplitude and duration. p38MAPK and MK2 concentrations and responsiveness toward IL-1ß were quantitatively compared between hepatocytes and macrophages. In macrophages, the absolute p38MAPK and MK2 concentration was significantly higher. Finally, in line with experimental observations, the mathematical model predicted a significantly higher half-maximal effective concentration for IL-1ß-induced pathway activation in macrophages compared with hepatocytes, underscoring the importance of cell type-specific differences in pathway regulation.

Identification of Cell Type-Specific Differences in Erythropoietin Receptor Signaling in Primary Erythroid and Lung Cancer Cells.

Lung cancer, with its most prevalent form non-small-cell lung carcinoma (NSCLC), is one of the leading causes of cancer-related deaths worldwide, and is commonly treated with chemotherapeutic drugs such as cisplatin. Lung cancer patients frequently suffer from chemotherapy-induced anemia, which can be treated with erythropoietin (EPO). However, studies have indicated that EPO not only promotes erythropoiesis in hematopoietic cells, but may also enhance survival of NSCLC cells. Here, we verified that the NSCLC cell line H838 expresses functional erythropoietin receptors (EPOR) and that treatment with EPO reduces cisplatin-induced apoptosis. To pinpoint differences in EPO-induced survival signaling in erythroid progenitor cells (CFU-E, colony forming unit-erythroid) and H838 cells, we combined mathematical modeling with a method for feature selection, the L1 regularization. Utilizing an example model and simulated data, we demonstrated that this approach enables the accurate identification and quantification of cell type-specific parameters. We applied our strategy to quantitative time-resolved data of EPO-induced JAK/STAT signaling generated by quantitative immunoblotting, mass spectrometry and quantitative real-time PCR (qRT-PCR) in CFU-E and H838 cells as well as H838 cells overexpressing human EPOR (H838-HA-hEPOR). The established parsimonious mathematical model was able to simultaneously describe the data sets of CFU-E, H838 and H838-HA-hEPOR cells. Seven cell type-specific parameters were identified that included for example parameters for nuclear translocation of STAT5 and target gene induction. Cell type-specific differences in target gene induction were experimentally validated by qRT-PCR experiments. The systematic identification of pathway differences and sensitivities of EPOR signaling in CFU-E and H838 cells revealed potential targets for intervention to selectively inhibit EPO-induced signaling in the tumor cells but leave the responses in erythroid progenitor cells unaffected. Thus, the proposed modeling strategy can be employed as a general procedure to identify cell type-specific parameters and to recommend treatment strategies for the selective targeting of specific cell types.

Disentangling the Complexity of HGF Signaling by Combining Qualitative and Quantitative Modeling.

Signaling pathways are characterized by crosstalk, feedback and feedforward mechanisms giving rise to highly complex and cell-context specific signaling networks. Dissecting the underlying relations is crucial to predict the impact of targeted perturbations. However, a major challenge in identifying cell-context specific signaling networks is the enormous number of potentially possible interactions. Here, we report a novel hybrid mathematical modeling strategy to systematically unravel hepatocyte growth factor (HGF) stimulated phosphoinositide-3-kinase (PI3K) and mitogen activated protein kinase (MAPK) signaling, which critically contribute to liver regeneration. By combining time-resolved quantitative experimental data generated in primary mouse hepatocytes with interaction graph and ordinary differential equation modeling, we identify and experimentally validate a network structure that represents the experimental data best and indicates specific crosstalk mechanisms. Whereas the identified network is robust against single perturbations, combinatorial inhibition strategies are predicted that result in strong reduction of Akt and ERK activation. Thus, by capitalizing on the advantages of the two modeling approaches, we reduce the high combinatorial complexity and identify cell-context specific signaling networks.

Deciphering the mechanism behind Fibroblast Growth Factor (FGF) induced biphasic signal-response profiles.

BACKGROUND: The Fibroblast Growth Factor (FGF) pathway is driving various aspects of cellular responses in both normal and malignant cells. One interesting characteristic of this pathway is the biphasic nature of the cellular response to some FGF ligands like FGF2. Specifically, it has been shown that phenotypic behaviors controlled by FGF signaling, like migration and growth, reach maximal levels in response to intermediate concentrations, while high levels of FGF2 elicit weak responses. The mechanisms leading to the observed biphasic response remains unexplained. RESULTS: A combination of experiments and computational modeling was used to understand the mechanism behind the observed biphasic signaling responses. FGF signaling involves a tertiary surface interaction that we captured with a computational model based on Ordinary Differential Equations (ODEs). It accounts for FGF2 binding to FGF receptors (FGFRs) and heparan sulfate glycosaminoglycans (HSGAGs), followed by receptor-phosphorylation, activation of the FRS2 adapter protein and the Ras-Raf signaling cascade. Quantitative protein assays were used to measure the dynamics of phosphorylated ERK (pERK) in response to a wide range of FGF2 ligand concentrations on a fine-grained time scale for the squamous cell lung cancer cell line H1703. We developed a novel approach combining Particle Swarm Optimization (PSO) and feature-based constraints in the objective function to calibrate the computational model to the experimental data. The model is validated using a series of extracellular and intracellular perturbation experiments. We demonstrate that in silico model predictions are in accordance with the observed in vitro results. CONCLUSIONS: Using a combined approach of computational modeling and experiments we found that competition between binding of the ligand FGF2 to HSGAG and FGF receptor leads to the biphasic response. At low to intermediate concentrations of FGF2 there are sufficient free FGF receptors available for the FGF2-HSGAG complex to enable the formation of the trimeric signaling unit. At high ligand concentrations the ligand binding sites of the receptor become saturated and the trimeric signaling unit cannot be formed. This insight into the pathway is an important consideration for the pharmacological inhibition of this pathway.

Comparison of approaches for parameter identifiability analysis of biological systems.

MOTIVATION: Modeling of dynamical systems using ordinary differential equations is a popular approach in the field of Systems Biology. The amount of experimental data that are used to build and calibrate these models is often limited. In this setting, the model parameters may not be uniquely determinable. Structural or a priori identifiability is a property of the system equations that indicates whether, in principle, the unknown model parameters can be determined from the available data. RESULTS: We performed a case study using three current approaches for structural identifiability analysis for an application from cell biology. The approaches are conceptually different and are developed independently. The results of the three approaches are in agreement. We discuss strength and weaknesses of each of them and illustrate how they can be applied to real world problems. AVAILABILITY AND IMPLEMENTATION: For application of the approaches to further applications, code representations (DAISY, Mathematica and MATLAB) for benchmark model and data are provided on the authors webpage. CONTACT: andreas.raue@fdm.uni-freiburg.de.

A Systems Biology Study on NFκB Signaling in Primary Mouse Hepatocytes.

The cytokine tumor necrosis factor-alpha (TNFα) is one of the key factors during the priming phase of liver regeneration as well as in hepatocarcinogenesis. TNFα activates the nuclear factor κ-light-chain-enhancer of activated B cells (NFκB) signaling pathway and contributes to the conversion of quiescent hepatocytes to activated hepatocytes that are able to proliferate in response to growth factor stimulation. Different mathematical models have been previously established for TNFα/NFκB signaling in the context of tumor cells. Combining these mathematical models with time-resolved measurements of expression and phosphorylation of TNFα/NFκB pathway constituents in primary mouse hepatocytes revealed that an additional phosphorylation step of the NFκB isoform p65 has to be considered in the mathematical model in order to sufficiently describe the dynamics of pathway activation in the primary cells. Also, we addressed the role of basal protein turnover by experimentally measuring the degradation rate of pivotal players in the absence of TNFα and including this information in the model. To elucidate the impact of variations in the protein degradation rates on TNFα/NFκB signaling on the overall dynamic behavior we used global sensitivity analysis that accounts for parameter uncertainties and showed that degradation and translation of p65 had a major impact on the amplitude and the integral of p65 phosphorylation. Finally, our mathematical model of TNFα/NFκB signaling was able to predict the time-course of the complex formation of p65 and of the inhibitor of NFκB (IκB) in primary mouse hepatocytes, which was experimentally verified. Hence, we here present a mathematical model for TNFα/NFκB signaling in primary mouse hepatocytes that provides an important basis to quantitatively disentangle the complex interplay of multiple factors in liver regeneration and tumorigenesis.

Dynamic mathematical modeling of IL13-induced signaling in Hodgkin and primary mediastinal B-cell lymphoma allows prediction of therapeutic targets.

Primary mediastinal B-cell lymphoma (PMBL) and classical Hodgkin lymphoma (cHL) share a frequent constitutive activation of JAK (Janus kinase)/STAT signaling pathway. Because of complex, nonlinear relations within the pathway, key dynamic properties remained to be identified to predict possible strategies for intervention. We report the development of dynamic pathway models based on quantitative data collected on signaling components of JAK/STAT pathway in two lymphoma-derived cell lines, MedB-1 and L1236, representative of PMBL and cHL, respectively. We show that the amounts of STAT5 and STAT6 are higher whereas those of SHP1 are lower in the two lymphoma cell lines than in normal B cells. Distinctively, L1236 cells harbor more JAK2 and less SHP1 molecules per cell than MedB-1 or control cells. In both lymphoma cell lines, we observe interleukin-13 (IL13)-induced activation of IL4 receptor α, JAK2, and STAT5, but not of STAT6. Genome-wide, 11 early and 16 sustained genes are upregulated by IL13 in both lymphoma cell lines. Specifically, the known STAT-inducible negative regulators CISH and SOCS3 are upregulated within 2 hours in MedB-1 but not in L1236 cells. On the basis of this detailed quantitative information, we established two mathematical models, MedB-1 and L1236 model, able to describe the respective experimental data. Most of the model parameters are identifiable and therefore the models are predictive. Sensitivity analysis of the model identifies six possible therapeutic targets able to reduce gene expression levels in L1236 cells and three in MedB-1. We experimentally confirm reduction in target gene expression in response to inhibition of STAT5 phosphorylation, thereby validating one of the predicted targets.

Covering a broad dynamic range: information processing at the erythropoietin receptor.

Cell surface receptors convert extracellular cues into receptor activation, thereby triggering intracellular signaling networks and controlling cellular decisions. A major unresolved issue is the identification of receptor properties that critically determine processing of ligand-encoded information. We show by mathematical modeling of quantitative data and experimental validation that rapid ligand depletion and replenishment of the cell surface receptor are characteristic features of the erythropoietin (Epo) receptor (EpoR). The amount of Epo-EpoR complexes and EpoR activation integrated over time corresponds linearly to ligand input; this process is carried out over a broad range of ligand concentrations. This relation depends solely on EpoR turnover independent of ligand binding, which suggests an essential role of large intracellular receptor pools. These receptor properties enable the system to cope with basal and acute demand in the hematopoietic system.