An Integrated Framework to Model Cellular Phenotype as a Component of Biochemical Networks

Identification of regulatory molecules in signaling pathways is critical for understanding cellular behavior. Given the complexity of the transcriptional gene network, the relationship between molecular expression and phenotype is difficult to determine using reductionist experimental methods. Computational models provide the means to characterize regulatory mechanisms and predict phenotype in the context of gene networks. Integrating gene expression data with phenotypic data in transcriptional network models enables systematic identification of critical molecules in a biological network. We developed an approach based on fuzzy logic to model cell budding in Saccharomyces cerevisiae using time series expression microarray data of the cell cycle. Cell budding is a phenotype of viable cells undergoing division. Predicted interactions between gene expression and phenotype reflected known biological relationships. Dynamic simulation analysis reproduced the behavior of the yeast cell cycle and accurately identified genes and interactions which are essential for cell viability.

[1]  E. Queralt,et al.  Functional Connection Between the Clb5 Cyclin, the Protein Kinase C Pathway and the Swi4 Transcription Factor in Saccharomyces cerevisiae , 2005, Genetics.

[2]  Cheng Li,et al.  Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application , 2001, Genome Biology.

[3]  K. Nasmyth,et al.  Yeast G1 cyclins CLN1 and CLN2 and a GAP‐like protein have a role in bud formation. , 1993, The EMBO journal.

[4]  A. Nordheim,et al.  Mcm1 is required to coordinate G2-specific transcription in Saccharomyces cerevisiae , 1995, Molecular and cellular biology.

[5]  Katherine C. Chen,et al.  Kinetic analysis of a molecular model of the budding yeast cell cycle. , 2000, Molecular biology of the cell.

[6]  Alan Wells,et al.  Modeling of signal-response cascades using decision tree analysis , 2005, Bioinform..

[7]  James E. Andrews,et al.  Combinatorial rule explosion eliminated by a fuzzy rule configuration , 1998, IEEE Trans. Fuzzy Syst..

[8]  Markus J. Herrgård,et al.  Integrating high-throughput and computational data elucidates bacterial networks , 2004, Nature.

[9]  K Nasmyth,et al.  Control of the yeast cell cycle by the Cdc28 protein kinase. , 1993, Current opinion in cell biology.

[10]  M. Kuo,et al.  A library of yeast genomic MCM1 binding sites contains genes involved in cell cycle control, cell wall and membrane structure, and metabolism , 1994, Molecular and cellular biology.

[11]  M. Mendenhall,et al.  An inhibitor of p34CDC28 protein kinase activity from Saccharomyces cerevisiae. , 1993, Science.

[12]  K. Nasmyth,et al.  A new role for MCM1 in yeast: cell cycle regulation of SW15 transcription. , 1991, Genes & development.

[13]  L. Glass,et al.  The logical analysis of continuous, non-linear biochemical control networks. , 1973, Journal of theoretical biology.

[14]  L. Hartwell.  Saccharomyces cerevisiae cell cycle. , 1974, Bacteriological reviews.

[15]  L. Breeden,et al.  A novel Mcm1-dependent element in the SWI4, CLN3, CDC6, and CDC47 promoters activates M/G1-specific transcription. , 1997, Genes & development.

[16]  Andrew A. Quong,et al.  Linear fuzzy gene network models obtained from microarray data by exhaustive search , 2004, BMC Bioinformatics.

[17]  K Nasmyth,et al.  Cdc6 is an unstable protein whose de novo synthesis in G1 is important for the onset of S phase and for preventing a ‘reductional’ anaphase in the budding yeast Saccharomyces cerevisiae. , 1995, The EMBO journal.

[18]  Q. Ouyang,et al.  The yeast cell-cycle network is robustly designed. , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Ronald W. Davis,et al.  Functional profiling of the Saccharomyces cerevisiae genome , 2002, Nature.

[20]  J. Mendel.  Fuzzy logic systems for engineering: a tutorial , 1995, Proc. IEEE.

[21]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[22]  K Nasmyth,et al.  The transcription factor Swi5 regulates expression of the cyclin kinase inhibitor p40SIC1 , 1996, Molecular and cellular biology.

[23]  Katherine C. Chen,et al.  Integrative analysis of cell cycle control in budding yeast. , 2004, Molecular biology of the cell.

[24]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[25]  Shoshana J. Wodak,et al.  CYGD: the Comprehensive Yeast Genome Database , 2004, Nucleic Acids Res..

[26]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[27]  Terence P. Speed,et al.  A benchmark for Affymetrix GeneChip expression measures , 2004, Bioinform..

[28]  Leif H. Finkel,et al.  BIOENGINEERING MODELS OF CELL SIGNALING , 2007 .

[29]  Kara Dolinski,et al.  Saccharomyces Genome Database provides mutant phenotype data , 2009, Nucleic Acids Res..

[30]  D. Botstein,et al.  Genomic binding sites of the yeast cell-cycle transcription factors SBF and MBF , 2001, Nature.

[31]  Suman Datta,et al.  Accelerated search for biomolecular network models to interpret high-throughput experimental data , 2006, BMC Bioinformatics.

[32]  R. Tibshirani,et al.  Significance analysis of microarrays applied to the ionizing radiation response , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[34]  D Kindelberger,et al.  Cell cycle-regulated transcription of the CLB2 gene is dependent on Mcm1 and a ternary complex factor , 1995, Molecular and cellular biology.

[35]  Michael Ruogu Zhang,et al.  Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. , 1998, Molecular biology of the cell.

[36]  Julio Saez-Rodriguez,et al.  Fuzzy Logic Analysis of Kinase Pathway Crosstalk in TNF/EGF/Insulin-Induced Signaling , 2007, PLoS Comput. Biol..

[37]  John J. Wyrick,et al.  Genome-wide location and function of DNA binding proteins. , 2000, Science.

[38]  Shan Wu,et al.  Decision tree modeling predicts effects of inhibiting contractility signaling on cell motility , 2007, BMC Systems Biology.

[39]  M. Mendenhall,et al.  Regulation of Cdc28 Cyclin-Dependent Protein Kinase Activity during the Cell Cycle of the Yeast Saccharomyces cerevisiae , 1998, Microbiology and Molecular Biology Reviews.

[40]  Jason E. Stewart,et al.  Minimum information about a microarray experiment (MIAME)—toward standards for microarray data , 2001, Nature Genetics.

[41]  F. Cross,et al.  Testing a mathematical model of the yeast cell cycle. , 2002, Molecular biology of the cell.

[42]  B. Andrews,et al.  The cyclin family of budding yeast: abundant use of a good idea. , 1998, Trends in genetics : TIG.

[43]  C. Li,et al.  Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. , 2001, Proceedings of the National Academy of Sciences of the United States of America.