Setting the standards for signal transduction research

Multiplexed Single-Cell Plasmonic Immunoassay of Intracellular Signaling Proteins Enables Non-Destructive Monitoring of Cell Fate

It is of significant importance in cancer biology to identify signaling pathways that play key roles in cell fate determination. Dissecting cellular signaling pathways requires the measurement of a large number of signaling proteins. However, tools for simultaneously monitoring multiple signaling pathway components in single living cells remain limited at present. Herein, we describe an approach, termed multiplexed single-cell plasmonic immunosandwich assay (mxscPISA), for simultaneous detection of multiple signaling proteins in individual living cells. This approach enabled simultaneous non-destructive monitoring of multiple (up to five, currently the highest multiplexing capacity in living cells) cytoplasmic and nucleus signaling proteins in individual cells with ultrahigh detection sensitivity. As a proof of principle, the epidermal growth factor receptor (EGFR) pathway, which plays a central role in cell fate determination, was investigated using this approach in this study. We found that there were differential attenuation rate of pro-survival and accumulation rate of pro-death signaling protein of the EGFR pathway in response to EGFR inactivation. These findings implicate that, after EGFR inactivation, a transient imbalance between survival and apoptotic signaling outputs contributed to the final cell fate of death. The mxscPISA approach can be a promising tool to reveal a signaling dynamic pattern at the single-cell level and to identify key components of signaling pathways that contribute to the final cell fate using only a limited number of cells.

Network reconstruction based on proteomic data and prior knowledge of protein connectivity using graph theory

Modeling of signal transduction pathways is instrumental for understanding cells’ function. People have been tackling modeling of signaling pathways in order to accurately represent the signaling events inside cells’ biochemical microenvironment in a way meaningful for scientists in a biological field. In this article, we propose a method to interrogate such pathways in order to produce cell-specific signaling models. We integrate available prior knowledge of protein connectivity, in a form of a Prior Knowledge Network (PKN) with phosphoproteomic data to construct predictive models of the protein connectivity of the interrogated cell type. Several computational methodologies focusing on pathways’ logic modeling using optimization formulations or machine learning algorithms have been published on this front over the past few years. Here, we introduce a light and fast approach that uses a breadth-first traversal of the graph to identify the shortest pathways and score proteins in the PKN, fitting the dependencies extracted from the experimental design. The pathways are then combined through a heuristic formulation to produce a final topology handling inconsistencies between the PKN and the experimental scenarios. Our results show that the algorithm we developed is efficient and accurate for the construction of medium and large scale signaling networks. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGF/TNFA stimulation against made up experimental data. To avoid the possibility of erroneous predictions, we performed a cross-validation analysis. Finally, we validate that the introduced approach generates predictive topologies, comparable to the ILP formulation. Overall, an efficient approach based on graph theory is presented herein to interrogate protein–protein interaction networks and to provide meaningful biological insights.

Modeling of signaling pathways in chondrocytes based on phosphoproteomic and cytokine release data

OBJECTIVE

Chondrocyte signaling is widely identified as a key component in cartilage homeostasis. Dysregulations of the signaling processes in chondrocytes often result in degenerative diseases of the tissue. Traditionally, the literature has focused on the study of major players in chondrocyte signaling, but without considering the cross-talks between them. In this paper, we systematically interrogate the signal transduction pathways in chondrocytes, on both the phosphoproteomic and cytokine release levels.

METHODS

The signaling pathways downstream 78 receptors of interest are interrogated. On the phosphoproteomic level, 17 key phosphoproteins are measured upon stimulation with single treatments of 78 ligands. On the cytokine release level, 55 cytokines are measured in the supernatant upon stimulation with the same treatments. Using an Integer Linear Programming (ILP) formulation, the proteomic data is combined with a priori knowledge of proteins' connectivity to construct a mechanistic model, predictive of signal transduction in chondrocytes.

RESULTS

We were able to validate previous findings regarding major players of cartilage homeostasis and inflammation (e.g., IL1B, TNF, EGF, TGFA, INS, IGF1 and IL6). Moreover, we studied pro-inflammatory mediators (IL1B and TNF) together with pro-growth signals for investigating their role in chondrocytes hypertrophy and highlighted the role of underreported players such as Inhibin beta A (INHBA), Defensin beta 1 (DEFB1), CXCL1 and Flagellin, and uncovered the way they cross-react in the phosphoproteomic level.

CONCLUSIONS

The analysis presented herein, leveraged high throughput proteomic data via an ILP formulation to gain new insight into chondrocytes signaling and the pathophysiology of degenerative diseases in articular cartilage.

Phosphoproteomics in drug discovery

Several important aspects of the drug discovery process, including target identification, mechanism of action determination and biomarker identification as well as drug repositioning, require complete understanding of the effects of drugs on protein phosphorylation in relevant biological systems. Novel high-throughput phosphoproteomic technologies can be employed to measure these phosphorylation events. In this review, we describe the advantages and limitations of state-of-the-art phosphoproteomic approaches such as mass spectrometry and antibody-based technologies in terms of sample and data throughput as well as data quality. We then discuss how datasets from each technology can be analyzed and how the results can be and have been applied to advance different aspects of the drug discovery process.

Detecting and removing inconsistencies between experimental data and signaling network topologies using integer linear programming on interaction graphs

Cross-referencing experimental data with our current knowledge of signaling network topologies is one central goal of mathematical modeling of cellular signal transduction networks. We present a new methodology for data-driven interrogation and training of signaling networks. While most published methods for signaling network inference operate on Bayesian, Boolean, or ODE models, our approach uses integer linear programming (ILP) on interaction graphs to encode constraints on the qualitative behavior of the nodes. These constraints are posed by the network topology and their formulation as ILP allows us to predict the possible qualitative changes (up, down, no effect) of the activation levels of the nodes for a given stimulus. We provide four basic operations to detect and remove inconsistencies between measurements and predicted behavior: (i) find a topology-consistent explanation for responses of signaling nodes measured in a stimulus-response experiment (if none exists, find the closest explanation); (ii) determine a minimal set of nodes that need to be corrected to make an inconsistent scenario consistent; (iii) determine the optimal subgraph of the given network topology which can best reflect measurements from a set of experimental scenarios; (iv) find possibly missing edges that would improve the consistency of the graph with respect to a set of experimental scenarios the most. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGFR/ErbB signaling against a library of high-throughput phosphoproteomic data measured in primary hepatocytes. Our methods detect interactions that are likely to be inactive in hepatocytes and provide suggestions for new interactions that, if included, would significantly improve the goodness of fit. Our framework is highly flexible and the underlying model requires only easily accessible biological knowledge. All related algorithms were implemented in a freely available toolbox SigNetTrainer making it an appealing approach for various applications.

Cellular signal transduction is orchestrated by communication networks of signaling proteins commonly depicted on signaling pathway maps. However, each cell type may have distinct variants of signaling pathways, and wiring diagrams are often altered in disease states. The identification of truly active signaling topologies based on experimental data is therefore one key challenge in systems biology of cellular signaling. We present a new framework for training signaling networks based on interaction graphs (IG). In contrast to complex modeling formalisms, IG capture merely the known positive and negative edges between the components. This basic information, however, already sets hard constraints on the possible qualitative behaviors of the nodes when perturbing the network. Our approach uses Integer Linear Programming to encode these constraints and to predict the possible changes (down, neutral, up) of the activation levels of the involved players for a given experiment. Based on this formulation we developed several algorithms for detecting and removing inconsistencies between measurements and network topology. Demonstrated by EGFR/ErbB signaling in hepatocytes, our approach delivers direct conclusions on edges that are likely inactive or missing relative to canonical pathway maps. Such information drives the further elucidation of signaling network topologies under normal and pathological phenotypes.

Using network biology to bridge pharmacokinetics and pharmacodynamics in oncology

If mathematical modeling is to be used effectively in cancer drug development, future models must take into account both the mechanistic details of cellular signal transduction networks and the pharmacokinetics (PK) of drugs used to inhibit their oncogenic activity. In this perspective, we present an approach to building multiscale models that capture systems-level architectural features of oncogenic signaling networks, and describe how these models can be used to design combination therapies and identify predictive biomarkers in silico.

Measurement and modeling of signaling at the single-cell level

It has long been recognized that a deeper understanding of cell function, with respect to execution of phenotypic behaviors and their regulation by the extracellular environment, is likely to be achieved by analyzing the underlying molecular processes for individual cells selected from across a population, rather than averages of many cells comprising that population. In recent years, experimental and computational methods for undertaking these analyses have advanced rapidly. In this review, we provide a perspective on both measurement and modeling facets of biochemistry at a single-cell level. Our central focus is on receptor-mediated signaling networks that regulate cell phenotypic functions.

State-time spectrum of signal transduction logic models

Despite the current wealth of high-throughput data, our understanding of signal transduction is still incomplete. Mathematical modeling can be a tool to gain an insight into such processes. Detailed biochemical modeling provides deep understanding, but does not scale well above relatively a few proteins. In contrast, logic modeling can be used where the biochemical knowledge of the system is sparse and, because it is parameter free (or, at most, uses relatively a few parameters), it scales well to large networks that can be derived by manual curation or retrieved from public databases. Here, we present an overview of logic modeling formalisms in the context of training logic models to data, and specifically the different approaches to modeling qualitative to quantitative data (state) and dynamics (time) of signal transduction. We use a toy model of signal transduction to illustrate how different logic formalisms (Boolean, fuzzy logic and differential equations) treat state and time. Different formalisms allow for different features of the data to be captured, at the cost of extra requirements in terms of computational power and data quality and quantity. Through this demonstration, the assumptions behind each formalism are discussed, as well as their advantages and disadvantages and possible future developments.

Construction of large signaling pathways using an adaptive perturbation approach with phosphoproteomic data

Construction of large and cell-specific signaling pathways is essential to understand information processing under normal and pathological conditions. On this front, gene-based approaches offer the advantage of large pathway exploration whereas phosphoproteomic approaches offer a more reliable view of pathway activities but are applicable to small pathway sizes. In this paper, we demonstrate an experimentally adaptive approach to construct large signaling pathways from phosphoproteomic data within a 3-day time frame. Our approach--taking advantage of the fast turnaround time of the xMAP technology--is carried out in four steps: (i) screen optimal pathway inducers, (ii) select the responsive ones, (iii) combine them in a combinatorial fashion to construct a phosphoproteomic dataset, and (iv) optimize a reduced generic pathway via an Integer Linear Programming formulation. As a case study, we uncover novel players and their corresponding pathways in primary human hepatocytes by interrogating the signal transduction downstream of 81 receptors of interest and constructing a detailed model for the responsive part of the network comprising 177 species (of which 14 are measured) and 365 interactions.

Predictive mathematical models of cancer signalling pathways

Complex intracellular signalling networks integrate extracellular signals and convert them into cellular responses. In cancer cells, the tightly regulated and fine-tuned dynamics of information processing in signalling networks is altered, leading to uncontrolled cell proliferation, survival and migration. Systems biology combines mathematical modelling with comprehensive, quantitative, time-resolved data and is most advanced in addressing dynamic properties of intracellular signalling networks. Here, we introduce different modelling approaches and their application to medical systems biology, focusing on the identifiability of parameters in ordinary differential equation models and their importance in network modelling to predict cellular decisions. Two related examples are given, which include processing of ligand-encoded information and dual feedback regulation in erythropoietin (Epo) receptor signalling. Finally, we review the current understanding of how systems biology could foster the development of new treatment strategies in the context of lung cancer and anaemia.