Parameter estimation for gene regulatory networks from microarray data: cold shock response in Saccharomyces cerevisiae.

*(English)*Zbl 1336.92026Summary: We investigated the dynamics of a gene regulatory network controlling the cold shock response in budding yeast, Saccharomyces cerevisiae. The medium-scale network, derived from published genome-wide location data, consists of 21 transcription factors that regulate one another through 31 directed edges. The expression levels of the individual transcription factors were modeled using mass balance ordinary differential equations with a sigmoidal production function. Each equation includes a production rate, a degradation rate, weights that denote the magnitude and type of influence of the connected transcription factors (activation or repression), and a threshold of expression. The inverse problem of determining model parameters from observed data is our primary interest. We fit the differential equation model to published microarray data using a penalized nonlinear least squares approach. Model predictions fit the experimental data well, within the 95 % confidence interval. Tests of the
model using randomized initial guesses and model-generated data also lend confidence to the fit. The results have revealed activation and repression relationships between the transcription factors. Sensitivity analysis indicates that the model is most sensitive to changes in the production rate parameters, weights, and thresholds of Yap1, Rox1, and Yap6, which form a densely connected core in the network. The modeling results newly suggest that Rap1, Fhl1, Msn4, Rph1, and Hsf1 play an important role in regulating the early response to cold shock in yeast. Our results demonstrate that estimation for a large number of parameters can be successfully performed for nonlinear dynamic gene regulatory networks using sparse, noisy microarray data.

##### MSC:

92C40 | Biochemistry, molecular biology |

92C42 | Systems biology, networks |

92D10 | Genetics and epigenetics |

##### Software:

nlmdl
PDF
BibTeX
XML
Cite

\textit{K. D. Dahlquist} et al., Bull. Math. Biol. 77, No. 8, 1457--1492 (2015; Zbl 1336.92026)

Full Text:
DOI

##### References:

[1] | Ackleh, AS; Fitzpatrick, BG; Scribner, R; Simonsen, N; Thibodeaux, JJ, Ecosystem modeling of college drinking: parameter estimation and comparing models to data, Math Comput Model, 50, 481-497, (2009) · Zbl 1185.91141 |

[2] | Aguilera, J; Randez-Gil, F; Prieto, JA, Cold response in saccharomyces cerevisiae: new functions for old mechanisms, FEMS Microbiol Rev, 31, 327-341, (2007) |

[3] | Al-Fageeh, MB; Smales, CM, Control and regulation of the cellular responses to cold shock: the responses in yeast and Mammalian systems, Biochem J, 397, 247-259, (2006) |

[4] | Alon U (2007) An introduction to systems biology: design principles of biological circuits. Chapman & Hall/CRC, Boca Raton · Zbl 1141.92002 |

[5] | Angeli, D; Hirsch, MW; Sontag, ED, Attractors in coherent systems of differential equations, J Differ Equ, 246, 3058-3076, (2009) · Zbl 1172.34036 |

[6] | Bailey, KR; Fitzpatrick, BG, Estimation of groundwater flow parameters using least squares, Math Comput Model, 26, 117-127, (1997) · Zbl 1185.76821 |

[7] | Banks, HT; Fitzpatrick, BG, Statistical methods for model comparison in parameter estimation problems for distributed systems, J Math Biol, 28, 501-527, (1990) · Zbl 0732.62061 |

[8] | Belle, A; Tanay, A; Bitincka, L; Shamir, R; O’Shea, EK, Quantification of protein half-lives in the budding yeast proteome, Proc Natl Acad Sci USA, 103, 13004-13009, (2006) · Zbl 1355.70003 |

[9] | Berger JO (1993) Statistical decision theory and Bayesian analysis, 2nd edn. Springer, New York |

[10] | Blossey, R; Cardelli, L; Phillips, A, Compositionality, stochasticity, and cooperativity in dynamic models of gene regulation, HFSP J, 2, 17-28, (2008) |

[11] | Cao, J; Zhao, H, Estimating dynamic models for gene regulation networks, Bioinformatics, 24, 1619-1624, (2008) |

[12] | Causton, HC; Ren, B; Koh, SS; Harbison, CT; Kanin, E; Jennings, EG; Lee, TI; True, HL; Lander, ES; Young, RA, Remodeling of yeast genome expression in response to environmental changes, Mol Biol Cell, 12, 323-337, (2001) |

[13] | Chen, K-C; Wang, T-Y; Tseng, H-H; Huang, C-YF; Kao, C-Y, A stochastic differential equation model for quantifying transcriptional regulatory network in saccharomyces cerevisiae, Bioinformatics, 21, 2883-2890, (2005) |

[14] | Chen, T; He, HL; Church, GM, Modeling gene expression with differential equations, Pac Symp Biocomput, 4, 4, (1999) |

[15] | Climescu-Haulica, A; Quirk, MD, A stochastic differential equation model for transcriptional regulatory networks, BMC Bioinf, 8, s4, (2007) |

[16] | Dawes, IW; Dickinson, JR (ed.); Schweizer, M (ed.), Stress responses, 376-438, (2004), Boca Raton |

[17] | Fan M, Kuwahara H, Wang X, Wang S, Gao X (2015) Parameter estimation methods for gene circuit modeling from time-series mRNA data: a comparative study. Brief Bioinf bbv015. doi:10.1093/bib/bbv015 |

[18] | Fitzpatrick, BG, Parameter estimation in conservation laws, J Math Syst Est Control, 3, 413-425, (1993) · Zbl 0785.93027 |

[19] | Fitzpatrick, BG, Statistical considerations and techniques for understanding physiological data, modeling, and treatments, Cardiovasc Eng, 8, 135-143, (2008) |

[20] | Fitzpatrick, BG; Keeling, SL, On approximation in total variation penalization for image reconstruction and inverse problems, Numer Func Anal Opt, 18, 941-958, (1997) · Zbl 0904.65131 |

[21] | Gallant AR (1987) Nonlinear statistical models. Wiley, New York · Zbl 0611.62071 |

[22] | Gasch, AP; Spellman, PT; Kao, CM; Carmel-Harel, O; Eisen, MB; Storz, G; Botstein, D; Brown, PO, Genomic expression programs in the response of yeast cells to environmental changes, Mol Biol Cell, 11, 4241-4257, (2000) |

[23] | Golub, GH; Heath, M; Wahba, G, Generalized cross-validation as a method for choosing a good ridge parameter, Technometrics, 21, 215-223, (1979) · Zbl 0461.62059 |

[24] | Hansen, PC; O’Leary, DP, The use of the L-curve in the regularization of discrete ill-posed problems, SIAM J Sci Comput, 14, 1487-1503, (1993) · Zbl 0789.65030 |

[25] | Harbison, CT; Gordon, DB; Lee, TI; Rinaldi, NJ; Macisaac, KD; Danford, TW; Hannett, NM; Tagne, JB; Reynolds, DB; Yoo, J; Jennings, EG; Zeitlinger, J; Pokholok, DK; Kellis, M; Rolfe, PA; Takusagawa, KT; Lander, ES; Gifford, DK; Fraenkel, E; Young, RA, Transcriptional regulatory code of a eukaryotic genome, Nature, 431, 99-104, (2004) |

[26] | Hecker, M; Lambeck, S; Toepfer, S; Someren, E; Guthke, R, Gene regulatory network inference: data integration in dynamic models—a review, Biosystems, 96, 86-103, (2009) |

[27] | Huet S, Bouvier A, Poursat M-A, Jolivet E (2004) Statistical tools for nonlinear regression: a practical guide with S-PLUS and R examples, 2nd edn. Springer, New York · Zbl 1041.62053 |

[28] | Kandror, O; Bretschneider, N; Kreydin, E; Cavalieri, D; Goldberg, AL, Yeast adapt to near-freezing temperatures by STRE/msn2,4-dependent induction of trehalose synthesis and certain molecular chaperones, Mol Cell, 13, 771-781, (2004) |

[29] | Karlebach, G; Shamir, R, Modelling and analysis of gene regulatory networks, Nat Rev Mol Cell Biol, 9, 770-780, (2008) |

[30] | Kauffman, KJ; Prakash, P; Edwards, JS, Advances in flux balance analysis, Curr Opin Biotechnol, 14, 491-496, (2003) |

[31] | Kuwahara, H; Fan, M; Wang, S; Gao, X, A framework for scalable parameter estimation of gene circuit models using structural information, Bioinformatics, 29, i98-i107, (2013) |

[32] | Lee, TI; Rinaldi, NJ; Robert, F; Odom, DT; Bar-Joseph, Z; Gerber, GK; Hannett, NM; Harbison, CT; Thompson, CM; Simon, I; Zeitlinger, J; Jennings, EG; Murray, HL; Gordon, DB; Ren, B; Wyrick, JJ; Tagne, JB; Volkert, TL; Fraenkel, E; Gifford, DK; Young, RA, Transcriptional regulatory networks in saccharomyces cerevisiae, Science, 298, 799-804, (2002) |

[33] | Lillacci, G; Khammash, M, Parameter estimation and model selection in computational biology, PLoS Comput Biol, 6, e1000696, (2010) |

[34] | Mendoza, L; Xenarios, I, A method for the generation of standardized qualitative dynamical systems of regulatory networks, Theor Biol Med Model, 3, 13, (2006) |

[35] | Morano, KA; Grant, CM; Moye-Rowley, WS, The response to heat shock and oxidative stress in saccharomyces cerevisiae, Genetics, 190, 1157-1195, (2012) |

[36] | Murata, Y; Homma, T; Kitagawa, E; Momose, Y; Sato, MS; Odani, M; Shimizu, H; Hasegawa-Mizusawa, M; Matsumoto, R; Mizukami, S; Fujita, K; Parveen, M; Komatsu, Y; Iwahashi, H, Genome-wide expression analysis of yeast response during exposure to 4 degrees C, Extremophiles, 10, 117-128, (2006) |

[37] | Orzechowski Westholm, J; Tronnersjö, S; Nordberg, N; Olsson, I; Komorowski, J; Ronne, H, Gis1 and rph1 regulate glycerol and acetate metabolism in glucose depleted yeast cells, PLoS ONE, 7, e31577, (2012) |

[38] | Sahara, T; Goda, T; Ohgiya, S, Comprehensive expression analysis of time-dependent genetic responses in yeast cells to low temperature, J Biol Chem, 277, 50015-50021, (2002) |

[39] | Schade, B; Jansen, G; Whiteway, M; Entian, KD; Thomas, DY, Cold adaptation in budding yeast, Mol Biol Cell, 15, 5492-5502, (2004) |

[40] | Schlitt, T; Brazma, A, Current approaches to gene regulatory network modelling, BMC Bioinf, 8, s9, (2007) |

[41] | Shea, MA; Ackers, GK, The OR control system of bacteriophage lambda. A physical-chemical model for gene regulation, J Mol Biol, 181, 211-230, (1985) |

[42] | Shore, D; Nasmyth, K, Purification and cloning of a DNA binding protein from yeast that binds to both silencer and activator elements, Cell, 51, 721-732, (1987) |

[43] | Smolen, P; Baxter, DA; Byrne, JH, Modeling transcriptional control in gene networks-methods, recent results, and future directions, Bull Math Biol, 62, 247-292, (2000) · Zbl 1323.92088 |

[44] | Sontag, ED, Monotone and near-monotone biochemical networks, Syst Synth Biol, 1, 59-87, (2007) |

[45] | Stekel D (2003) Microarray bioinformatics. Cambridge University Press, Cambridge |

[46] | Tai, SL; Daran-Lapujade, P; Walsh, MC; Pronk, JT; Daran, J-M, Acclimation of saccharomyces cerevisiae to low temperature: a chemostat-based transcriptome analysis, Mol Biol Cell, 18, 5100-5112, (2007) |

[47] | Tang, L; Liu, X; Clarke, ND, Inferring direct regulatory targets from expression and genome location analyses: a comparison of transcription factor deletion and overexpression, BMC Genom, 7, 215, (2006) |

[48] | Thieringer HA, Jones PG, Inouye M (1998) Cold shock and adaptation. Bioessays 20:49-57. doi: 10.1002/(SICI)1521-1878(199801)20:1\(<\)3.0.CO;2-N · Zbl 1185.91141 |

[49] | Vohradský, J, Neural network model of gene expression, FASEB J, 15, 846-854, (2001) |

[50] | Vu, TT; Vohradsky, J, Nonlinear differential equation model for quantification of transcriptional regulation applied to microarray data of saccharomyces cerevisiae, Nucleic Acids Res, 35, 279-287, (2007) |

[51] | Wilkinson DJ (2006) Stochastic modelling for systems biology. Taylor & Francis, Boca Raton · Zbl 1099.92004 |

[52] | Winzeler, EA; Shoemaker, DD; Astromoff, A; Liang, H; Anderson, K; Andre, B; Bangham, R; Benito, R; Boeke, JD; Bussey, H; Chu, AM; Connelly, C; Davis, K; Dietrich, F; Dow, SW; Bakkoury, M; Foury, F; Friend, SH; Gentalen, E; Giaever, G; Hegemann, JH; Jones, T; Laub, M; Liao, H; Liebundguth, N; Lockhart, DJ; Lucau-Danila, A; Lussier, M; M’Rabet, N; Menard, P; Mittmann, M; Pai, C; Rebischung, C; Revuelta, JL; Riles, L; Roberts, CJ; Ross-MacDonald, P; Scherens, B; Snyder, M; Sookhai-Mahadeo, S; Storms, RK; Véronneau, S; Voet, M; Volckaert, G; Ward, TR; Wysocki, R; Yen, GS; Yu, K; Zimmermann, K; Philippsen, P; Johnston, M; Davis, RW, Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis, Science, 285, 901-906, (1999) |

[53] | Xiao, L; Grove, A, Coordination of ribosomal protein and ribosomal RNA gene expression in response to TOR signaling, Curr Genomics, 10, 198-205, (2009) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.