Identifying Drug Effects via Pathway Alterations using an Integer Linear Programming Optimization Formulation on Phosphoproteomic Data

Understanding the mechanisms of cell function and drug action is a major endeavor in the pharmaceutical industry. Drug effects are governed by the intrinsic properties of the drug (i.e., selectivity and potency) and the specific signaling transduction network of the host (i.e., normal vs. diseased cells). Here, we describe an unbiased, phosphoproteomicbased approach to identify drug effects by monitoring drug-induced topology alterations. With the proposed method, drug effects are investigated under several conditions on a cell-type specific signaling network. First, starting with a generic pathway made of logical gates, we build a cell-type specific map by constraining it to fit 13 key phopshoprotein signals under 55 experimental cases. Fitting is performed via a formulation as an Integer Linear Program (ILP) and solution by standard ILP solvers; a procedure that drastically outperforms previous fitting schemes. Then, knowing the cell topology, we monitor the same key phopshoprotein signals under the presence of drug and cytokines and we re-optimize the specific map to reveal the drug-induced topology alterations. To prove our case, we make a pathway map for the hepatocytic cell line HepG2 and we evaluate the effects of 4 drugs: 3 selective inhibitors for the Epidermal Growth Factor Receptor (EGFR) and a non selective drug. We confirm effects easily predictable from the drugs’ main target (i.e. EGFR inhibitors blocks the EGFR pathway) but we also uncover unanticipated effects due to either drug promiscuity or the cell’s specific topology. An interesting finding is that the selective EGFR inhibitor Gefitinib is able to inhibit signaling downstream the Interleukin-1alpha (IL-1α) pathway; an effect that cannot be extracted from binding affinity based approaches. Our method represents an unbiased approach to identify drug effects on a small to medium size pathways and is scalable to larger topologies with any type of signaling perturbations (small molecules, 3 RNAi etc). The method is a step towards a better picture of drug effects in pathways, the cornerstone in identifying the mechanisms of drug efficacy and toxicity

De Novo Inference of Systems-Level Mechanistic Models of Development from Live-Imaging-Based Phenotype Analysis

SummaryElucidation of complex phenotypes for mechanistic insights presents a significant challenge in systems biology. We report a strategy to automatically infer mechanistic models of cell fate differentiation based on live-imaging data. We use cell lineage tracing and combinations of tissue-specific marker expression to assay progenitor cell fate and detect fate changes upon genetic perturbation. Based on the cellular phenotypes, we further construct a model for how fate differentiation progresses in progenitor cells and predict cell-specific gene modules and cell-to-cell signaling events that regulate the series of fate choices. We validate our approach in C. elegans embryogenesis by perturbing 20 genes in over 300 embryos. The result not only recapitulates current knowledge but also provides insights into gene function and regulated fate choice, including an unexpected self-renewal. Our study provides a powerful approach for automated and quantitative interpretation of complex in vivo information

The how and why of lncRNA function: An innate immune perspective.

Next-generation sequencing has provided a more complete picture of the composition of the human transcriptome indicating that much of the "blueprint" is a vastness of poorly understood non-protein-coding transcripts. This includes a newly identified class of genes called long noncoding RNAs (lncRNAs). The lack of sequence conservation for lncRNAs across species meant that their biological importance was initially met with some skepticism. LncRNAs mediate their functions through interactions with proteins, RNA, DNA, or a combination of these. Their functions can often be dictated by their localization, sequence, and/or secondary structure. Here we provide a review of the approaches typically adopted to study the complexity of these genes with an emphasis on recent discoveries within the innate immune field. Finally, we discuss the challenges, as well as the emergence of new technologies that will continue to move this field forward and provide greater insight into the biological importance of this class of genes. This article is part of a Special Issue entitled: ncRNA in control of gene expression edited by Kotb Abdelmohsen

Algorithms for the analysis of protein interaction networks

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2012.Cataloged from PDF version of thesis.Includes bibliographical references (p. 107-117).In the decade since the human genome project, a major research trend in biology has been towards understanding the cell as a system. This interest has stemmed partly from a deeper appreciation of how important it is to understand the emergent properties of cellular systems (e.g., they seem to be the key to understanding diseases like cancer). It has also been enabled by new high-throughput techniques that have allowed us to collect new types of data at the whole-genome scale. We focus on one sub-domain of systems biology: the understanding of protein interactions. Such understanding is valuable: interactions between proteins are fundamental to many cellular processes. Over the last decade, high-throughput experimental techniques have allowed us to collect a large amount of protein-protein interaction (PPI) data for many species. A popular abstraction for representing this data is the protein interaction network: each node of the network represents a protein and an edge between two nodes represents a physical interaction between the two corresponding proteins. This abstraction has proven to be a powerful tool for understanding the systems aspects of protein interaction. We present some algorithms for the augmentation, cleanup and analysis of such protein interaction networks: 1. In many species, the coverage of known PPI data remains partial. Given two protein sequences, we describe an algorithm to predict if two proteins physically interact, using logistic regression and insights from structural biology. We also describe how our predictions may be further improved by combining with functional-genomic data. 2. We study systematic false positives in a popular experimental protocol, the Yeast 2-Hybrid method. Here, some "promiscuous" proteins may lead to many false positives. We describe a Bayesian approach to modeling and adjusting for this error. 3. Comparative analysis of PPI networks across species can provide valuable insights. We describe IsoRank, an algorithm for global network alignment of multiple PPI networks. The algorithm first constructs an eigenvalue problem that encapsulates the network and sequence similarity constraints. The solution of the problem describes a k-partite graph that is further processed to find the alignment. 4. For a given signaling network, we describe an algorithm that combines RNA-interference data with PPI data to produce hypotheses about the structure of the signaling network. Our algorithm constructs a multi-commodity flow problem that expresses the constraints described by the data and finds a sparse solution to it.by Rohit Singh.Ph.D

Cross-Talk Categorisations in Data-Driven Models of Signalling Networks: A System-Level View

Data-driven models of signalling networks are becoming increasingly important in systems biology in order to reflect the dynamic patterns of signalling activities in a context-specific manner. State-of-the-art approaches for categorising and detecting signalling cross-talks may not be suitable for such models since they rely on static topologies of cell signalling networks and prior biological knowledge. In this chapter, we review state-of-the-art approaches that categorise all possible cross-talks in signalling networks and propose a novel categorisation specific to data-driven network models. Considering such models as undirected networks, we propose two categories of signalling cross-talks between any two given signalling pathways. In a Type-I cross-talk, a signalling link {gi ,gj } connects two signalling pathways, where gi and gj are signalling nodes that belong to two distinct pathways. In a Type-II cross-talk, two signalling links {gi ,gj } and {gj ,gk } meet at the intersection of two signalling pathways at a shared signalling node gj . We compared our categorisation approach with others and found that all the types of cross-talks defined by those approaches can be mapped to Type-I and Type-II cross-talks when underlying signalling activities are considered as non-causal relationships. Next, we provided a simple but intuitive algorithm called XDaMoSiN (cross-talks in data-driven models of signalling networks) to detect both Type-I and Type-II cross-talks between any two given signalling pathways in a data-driven network model. By detecting cross-talks in such network models, our approach can be used to analyse and decipher latent mechanisms of various cell phenotypes, such as cancer or acquired drug resistance, that may evolve due to the highly adaptable and dynamic nature of signal transduction networks

The Local Edge Machine: inference of dynamic models of gene regulation

We present a novel approach, the Local Edge Machine, for the inference of regulatory interactions directly from time-series gene expression data. We demonstrate its performance, robustness, and scalability on in silico datasets with varying behaviors, sizes, and degrees of complexity. Moreover, we demonstrate its ability to incorporate biological prior information and make informative predictions on a well-characterized in vivo system using data from budding yeast that have been synchronized in the cell cycle. Finally, we use an atlas of transcription data in a mammalian circadian system to illustrate how the method can be used for discovery in the context of large complex networks.Department of Applied Mathematic

Discovery of progenitor cell signatures by time-series synexpression analysis during Drosophila embryonic cell immortalization

The use of time series profiling to identify groups of functionally related genes (synexpression groups) is a powerful approach for the discovery of gene function. Here we apply this strategy during RasV12 immortalization of Drosophila embryonic cells, a phenomenon not well characterized. Using high-resolution transcriptional time-series datasets, we generated a gene network based on temporal expression profile similarities. This analysis revealed that common immortalized cells are related to adult muscle precursors (AMPs), a stem cell-like population contributing to adult muscles and sharing properties with vertebrate satellite cells. Remarkably, the immortalized cells retained the capacity for myogenic differentiation when treated with the steroid hormone ecdysone. Further, we validated in vivo the transcription factor CG9650, the ortholog of mammalian Bcl11a/b, as a regulator of AMP proliferation predicted by our analysis. Our study demonstrates the power of time series synexpression analysis to characterize Drosophila embryonic progenitor lines and identify stem/progenitor cell regulators

Logic-Based Models for the Analysis of Cell Signaling Networks

Computational models are increasingly used to analyze the operation of complex biochemical networks, including those involved in cell signaling networks. Here we review recent advances in applying logic-based modeling to mammalian cell biology. Logic-based models represent biomolecular networks in a simple and intuitive manner without describing the detailed biochemistry of each interaction. A brief description of several logic-based modeling methods is followed by six case studies that demonstrate biological questions recently addressed using logic-based models and point to potential advances in model formalisms and training procedures that promise to enhance the utility of logic-based methods for studying the relationship between environmental inputs and phenotypic or signaling state outputs of complex signaling networks.National Institutes of Health (U.S.) (Grant P50- GM68762)National Institutes of Health (U.S.) (Grant U54-CA112967)United States. Dept. of Defense (Institute for Collaborative Biotechnologies

Plasticity of the MAPK signaling network in response to mechanical stress

Cells display versatile responses to mechanical inputs and recent studies have identified the mitogen-activated protein kinase (MAPK) cascades mediating the biological effects observed upon mechanical stimulation. Although, MAPK pathways can act insulated from each other, several mechanisms facilitate the crosstalk between the components of these cascades. Yet, the combinatorial complexity of potential molecular interactions between these elements have prevented the understanding of their concerted functions. To analyze the plasticity of the MAPK signaling network in response to mechanical stress we performed a non-saturating epistatic screen in resting and stretched conditions employing as readout a JNK responsive dJun-FRET biosensor. By knocking down MAPKs, and JNK pathway regulators, singly or in pairs in Drosophila S2R+ cells, we have uncovered unexpected regulatory links between JNK cascade kinases, Rho GTPases, MAPKs and the JNK phosphatase Puc. These relationships have been integrated in a system network model at equilibrium accounting for all experimentally validated interactions. This model allows predicting the global reaction of the network to its modulation in response to mechanical stress. It also highlights its context-dependent sensitivity