Plato's Cave Algorithm: Inferring Functional Signaling Networks from Early Gene Expression Shadows

Improving the ability to reverse engineer biochemical networks is a major goal of systems biology. Lesions in signaling networks lead to alterations in gene expression, which in principle should allow network reconstruction. However, the information about the activity levels of signaling proteins conveyed in overall gene expression is limited by the complexity of gene expression dynamics and of regulatory network topology. Two observations provide the basis for overcoming this limitation: a. genes induced without de-novo protein synthesis (early genes) show a linear accumulation of product in the first hour after the change in the cell's state; b. The signaling components in the network largely function in the linear range of their stimulus-response curves. Therefore, unlike most genes or most time points, expression profiles of early genes at an early time point provide direct biochemical assays that represent the activity levels of upstream signaling components. Such expression data provide the basis for an efficient algorithm (Plato's Cave algorithm; PLACA) to reverse engineer functional signaling networks. Unlike conventional reverse engineering algorithms that use steady state values, PLACA uses stimulated early gene expression measurements associated with systematic perturbations of signaling components, without measuring the signaling components themselves. Besides the reverse engineered network, PLACA also identifies the genes detecting the functional interaction, thereby facilitating validation of the predicted functional network. Using simulated datasets, the algorithm is shown to be robust to experimental noise. Using experimental data obtained from gonadotropes, PLACA reverse engineered the interaction network of six perturbed signaling components. The network recapitulated many known interactions and identified novel functional interactions that were validated by further experiment. PLACA uses the results of experiments that are feasible for any signaling network to predict the functional topology of the network and to identify novel relationships.

[1]  S. Sealfon,et al.  Transcriptome Fingerprints Distinguish Hallucinogenic and Nonhallucinogenic 5-Hydroxytryptamine 2A Receptor Agonist Effects in Mouse Somatosensory Cortex , 2003, The Journal of Neuroscience.

[2]  J. Hasty,et al.  Reverse engineering gene networks: Integrating genetic perturbations with dynamical modeling , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  R. Seger,et al.  Stimulation of Jun N-terminal kinase (JNK) by gonadotropin-releasing hormone in pituitary alpha T3-1 cell line is mediated by protein kinase C, c-Src, and CDC42. , 1998, Molecular endocrinology.

[4]  Claude E. Shannon,et al.  Prediction and Entropy of Printed English , 1951 .

[5]  S. Sealfon,et al.  Coupling of GnRH concentration and the GnRH receptor-activated gene program. , 2002, Molecular endocrinology.

[6]  P. Gehler,et al.  An introduction to graphical models , 2001 .

[7]  Mark P. Styczynski,et al.  Overview of computational methods for the inference of gene regulatory networks , 2005, Comput. Chem. Eng..

[8]  R. Milo,et al.  Network motifs in integrated cellular networks of transcription-regulation and protein-protein interaction. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[9]  S. Shen-Orr,et al.  Network motifs: simple building blocks of complex networks. , 2002, Science.

[10]  Changbao Wu,et al.  Jackknife, Bootstrap and Other Resampling Methods in Regression Analysis , 1986 .

[11]  D. Lauffenburger,et al.  A Systems Model of Signaling Identifies a Molecular Basis Set for Cytokine-Induced Apoptosis , 2005, Science.

[12]  C. Sander,et al.  Models from experiments: combinatorial drug perturbations of cancer cells , 2008, Molecular systems biology.

[13]  S. Sealfon,et al.  Mixed Analog/Digital Gonadotrope Biosynthetic Response to Gonadotropin-releasing Hormone* , 2006, Journal of Biological Chemistry.

[14]  J. Olefsky,et al.  GnRH activates ERK1/2 leading to the induction of c-fos and LHbeta protein expression in LbetaT2 cells. , 2002, Molecular endocrinology.

[15]  Prahlad T. Ram,et al.  Formation of Regulatory Patterns During Signal Propagation in a Mammalian Cellular Network , 2005, Science.

[16]  S. Shen-Orr,et al.  Networks Network Motifs : Simple Building Blocks of Complex , 2002 .

[17]  J. Olefsky,et al.  GnRH Activates ERK1/2 Leading to the Induction of c-fos and LHβ Protein Expression in LβT2 Cells , 2002 .

[18]  Stuart C. Sealfon,et al.  Hallucinogens Recruit Specific Cortical 5-HT2A Receptor-Mediated Signaling Pathways to Affect Behavior , 2007, Neuron.

[19]  B. Kholodenko,et al.  Quantification of information transfer via cellular signal transduction pathways , 1997, FEBS letters.

[20]  A. Loinger,et al.  Analysis of genetic toggle switch systems encoded on plasmids. , 2009, Physical review letters.

[21]  Martin A. Nowak,et al.  Inferring Cellular Networks Using Probabilistic Graphical Models , 2004 .

[22]  K. Fukunaga,et al.  Regulatory role of autophosphorylation of Ca2+/calmodulin-dependent protein kinase II. , 1990, Advances in second messenger and phosphoprotein research.

[23]  Paul A Clemons,et al.  The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease , 2006, Science.

[24]  H. Margalit,et al.  Stochastic Analysis of the SOS Response in Escherichia coli , 2009, PloS one.

[25]  S. Sealfon,et al.  Monitoring G-protein-coupled receptor signaling with DNA microarrays and real-time polymerase chain reaction. , 2002, Methods in enzymology.

[26]  S. Teichmann,et al.  Gene regulatory network growth by duplication , 2004, Nature Genetics.

[27]  Olaf Wolkenhauer,et al.  A unified framework for unraveling the functional interaction structure of a biomolecular network based on stimulus‐response experimental data , 2005, FEBS letters.

[28]  Florian Markowetz,et al.  How to Understand the Cell by Breaking It: Network Analysis of Gene Perturbation Screens , 2009, PLoS Comput. Biol..

[29]  Rachel B. Brem,et al.  Integrating large-scale functional genomic data to dissect the complexity of yeast regulatory networks , 2008, Nature Genetics.

[30]  Chun-Wei Hsu,et al.  Characterization of microRNA‐regulated protein‐protein interaction network , 2008, Proteomics.

[31]  Nicola J. Rinaldi,et al.  Transcriptional Regulatory Networks in Saccharomyces cerevisiae , 2002, Science.

[32]  Anne-Laure Boulesteix,et al.  Partial least squares: a versatile tool for the analysis of high-dimensional genomic data , 2006, Briefings Bioinform..

[33]  B. Kholodenko,et al.  Modular response analysis of cellular regulatory networks. , 2002, Journal of theoretical biology.

[34]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[35]  P. Rorsman,et al.  Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels. , 2005, Genome research.

[36]  N. Lee,et al.  Computational and experimental approaches for modeling gene regulatory networks. , 2007, Current pharmaceutical design.

[37]  S. Shen-Orr,et al.  Network motifs in the transcriptional regulation network of Escherichia coli , 2002, Nature Genetics.

[38]  R. Lefkowitz,et al.  Regulation of tyrosine kinase cascades by G-protein-coupled receptors. , 1999, Current opinion in cell biology.

[39]  K. Fukunaga,et al.  Serum and growth factors rapidly elicit phosphorylation of the Ca2+/calmodulin-dependent protein kinase II in intact quiescent rat 3Y1 cells. , 1988, The Journal of biological chemistry.

[40]  Ofer Biham,et al.  Stochastic simulations of genetic switch systems. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  Ronan M. T. Fleming,et al.  Genome-Scale Reconstruction of Escherichia coli's Transcriptional and Translational Machinery: A Knowledge Base, Its Mathematical Formulation, and Its Functional Characterization , 2009, PLoS Comput. Biol..

[42]  Nicola J. Rinaldi,et al.  Control of Pancreas and Liver Gene Expression by HNF Transcription Factors , 2004, Science.

[43]  Eduardo D. Sontag,et al.  Inferring dynamic architecture of cellular networks using time series of gene expression, protein and metabolite data , 2004, Bioinform..

[44]  John G. Albeck,et al.  Cue-Signal-Response Analysis of TNF-Induced Apoptosis by Partial Least Squares Regression of Dynamic Multivariate Data , 2004, J. Comput. Biol..

[45]  Nir Friedman,et al.  Inferring Cellular Networks Using Probabilistic Graphical Models , 2004, Science.

[46]  G. Friedlander,et al.  Regulation of gene expression by small non-coding RNAs: a quantitative view , 2007, Molecular systems biology.

[47]  S. Sealfon,et al.  Noise propagation and scaling in regulation of gonadotrope biosynthesis. , 2007, Biophysical journal.

[48]  A. Loinger,et al.  Stochastic simulations of the repressilator circuit. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[49]  Ezgi O. Booth,et al.  Epistasis analysis with global transcriptional phenotypes , 2005, Nature Genetics.

[50]  M. Gerstein,et al.  A Bayesian Networks Approach for Predicting Protein-Protein Interactions from Genomic Data , 2003, Science.

[51]  R. Lefkowitz,et al.  New mechanisms in heptahelical receptor signaling to mitogen activated protein kinase cascades , 2001, Oncogene.

[52]  S. Sealfon,et al.  Gonadotropin-releasing Hormone Receptor-coupled Gene Network Organization * , 2001 .

[53]  Andrew J. Bulpitt,et al.  A Primer on Learning in Bayesian Networks for Computational Biology , 2007, PLoS Comput. Biol..

[54]  J. Collins,et al.  Chemogenomic profiling on a genome-wide scale using reverse-engineered gene networks , 2005, Nature Biotechnology.