Estimating confidence intervals in predicted responses for oscillatory biological models

BMC Systems Biology

Francis J Doyle¹  Peter C StJohn¹ 

[1] Department of Chemical Engineering, University of California, Santa Barbara, CA 93106-5080, USA

关键词: Parameter estimation;    Sensitivity analysis;    Circadian rhythms;    Oscillatory models;    Identifiability;    Bootstrap;   

Others  :  1142468
DOI  :  10.1186/1752-0509-7-71

 received in 2013-01-05, accepted in 2013-07-15,  发布年份 2013

【 摘 要 】

Background

The dynamics of gene regulation play a crucial role in a cellular control: allowing the cell to express the right proteins to meet changing needs. Some needs, such as correctly anticipating the day-night cycle, require complicated oscillatory features. In the analysis of gene regulatory networks, mathematical models are frequently used to understand how a network’s structure enables it to respond appropriately to external inputs. These models typically consist of a set of ordinary differential equations, describing a network of biochemical reactions, and unknown kinetic parameters, chosen such that the model best captures experimental data. However, since a model’s parameter values are uncertain, and since dynamic responses to inputs are highly parameter-dependent, it is difficult to assess the confidence associated with these in silico predictions. In particular, models with complex dynamics - such as oscillations - must be fit with computationally expensive global optimization routines, and cannot take advantage of existing measures of identifiability. Despite their difficulty to model mathematically, limit cycle oscillations play a key role in many biological processes, including cell cycling, metabolism, neuron firing, and circadian rhythms.

Results

In this study, we employ an efficient parameter estimation technique to enable a bootstrap uncertainty analysis for limit cycle models. Since the primary role of systems biology models is the insight they provide on responses to rate perturbations, we extend our uncertainty analysis to include first order sensitivity coefficients. Using a literature model of circadian rhythms, we show how predictive precision is degraded with decreasing sample points and increasing relative error. Additionally, we show how this method can be used for model discrimination by comparing the output identifiability of two candidate model structures to published literature data.

Conclusions

Our method permits modellers of oscillatory systems to confidently show that a model’s dynamic characteristics follow directly from experimental data and model structure, relaxing assumptions on the particular parameters chosen. Ultimately, this work highlights the importance of continued collection of high-resolution data on gene and protein activity levels, as they allow the development of predictive mathematical models.

【 授权许可】

   
2013 St. John and Doyle; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

【 参考文献 】

• [1]Gutenkunst RN, Waterfall JJ, Casey FP, Brown KS, Myers CR, Sethna JP: Universally sloppy parameter sensitivities in systems biology models. PLoS Comput Biol 2007, 3(10):1871-1878.
• [2]Knowles JD, Kell DB: The role of modeling in systems biology. In System Modeling in, Cellular Biology: From Concepts to Nuts and Bolts. Edited by Szallasi Z, Stelling J, Periwal V. Cambridge: MIT Press; 2006:3-18.
• [3]Raue A, Kreutz C, Maiwald T, Bachmann J, Schilling M, Klingmüller U, Timmer J: Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood. Bioinformatics 2009, 25(15):1923-1929.
• [4]Nihtilä M, Virkkunen J: Practical identifiability of growth and substrate consumption models. Biotechnol Bioeng 1977, 19(12):1831-1850.
• [5]Jiménez-Hornero JE, Santos-Dueñas IM, García-García I: Optimization of biotechnological processes. The acetic acid fermentation. Part II: Practical identifiability analysis and parameter estimation. Biochem Eng J 2009, 45:7-21.
• [6]Holmberg A: On the practical identifiability of microbial growth models incorporating Michaelis-Menten type nonlinearities. Math Biosci 1982, 62:23-43.
• [7]Joshi M, Seidel-Morgenstern A, Kremling A: Exploiting the bootstrap method for quantifying parameter confidence intervals in dynamical systems. Metab Eng 2006, 8(5):447-455.
• [8]Goldbeter A, Keizer J: Biochemical Oscillations and Cellular Rhythms: The Molecular Bases of Periodic and Chaotic Behaviour Volume 51. Cambridge: Cambridge University Press; 1996.
• [9]Forger DB, Peskin CS: A detailed predictive model of the mammalian circadian clock. Proc Natl Acad Sci USA 2003, 100(25):14806-14811.
• [10]Leloup JC, Goldbeter A: Toward a detailed computational model for the mammalian circadian clock. Proc Natl Acad Sci USA 2003, 100(12):7051-7056.
• [11]Mirsky HP, Liu AC, Welsh DK, Kay SA, Doyle III FJ: A model of the cell-autonomous mammalian circadian clock. Proc Natl Acad Sci USA 2009, 106(27):11107-11112.
• [12]Wilkins AK, Tidor B, White J, Barton PI: Sensitivity analysis for oscillating dynamical systems. SIAM J Sci Comput 2009, 31(4):2706-2732.
• [13]Hirota T, Lee JW, St. John PC, Sawa M, Iwaisako K, Noguchi T, Pongsawakul PY, Sonntag T, Welsh DK, Brenner DA, Doyle III FJ, Schultz PG, Kay SA: Identification of small molecule activators of Cryptochrome. Science 2012, 337(6098):1094-1097.
• [14]Ko CH, Takahashi JS: Molecular components of the mammalian circadian clock. Hum Mol Genet 2006, 15 Spec No 2:R271-R277.
• [15]Doyle III FJ, Gunawan R, Bagheri N, Mirsky HP, To TL: Circadian rhythm: A natural, robust, multi-scale control system. Comput Chem Eng 2006, 30(10–12):1700-1711.
• [16]Herzog ED: Neurons and networks in daily rhythms. Nat Rev Neurosci 2007, 8(10):790-802.
• [17]Hughes ME, DiTacchio L, Hayes KR, Vollmers C, Pulivarthy S, Baggs JE, Panda S, Hogenesch JB: Harmonics of circadian gene transcription in mammals. PLoS Genet 2009, 5(4):e1000442.
• [18]Zhang EE, Liu AC, Hirota T, Miraglia LJ, Welch G, Pongsawakul PY, Liu X, Atwood A, Huss JW, Janes J, Su AI, Hogenesch JB, Kay SA: A genome-wide RNAi screen for modifiers of the circadian clock in human cells. Cell 2009, 139:199-210.
• [19]Hirota T, Lee JW, Lewis WG, Zhang EE, Breton G, Liu X, Garcia M, Peters EC, Etchegaray JP, Traver D, Schultz PG, Kay SA: High-throughput chemical screen identifies a novel potent modulator of cellular circadian rhythms and reveals CKIα, as a clock regulatory kinase. PLoS Biol 2010, 8(12):e1000559.
• [20]Hirota T, Lewis WG, Liu AC, Lee JW, Schultz PG, Kay SA: A chemical biology approach reveals period shortening of the mammalian circadian clock by specific inhibition of GSK-3beta. Proc Natl Acad Sci USA 2008, 105(52):20746-20751.
• [21]Biegler LT: Nonlinear Programming: Concepts, Algorithms, and Applications to Chemical Processes. Philidelphia: SIAM; 2010.
• [22]Bock HG, Kostina E, Schloder JP: Numerical methods for parameter estimation in nonlinear differential algebraic equations. GAMM-âĂR̆Mitteilungen 2007, 408(2):376-408.
• [23]Bure EG, Rosenwasser E: The study of the sensitivity of oscillatory systems. Autom Rem Contr 1974, 7:1045-1052.
• [24]Kramer M, Rabitz H, Calo J: Sensitivity analysis of oscillatory systems. Appl Math Modell 1984, 8(5):328-350.
• [25]Zak DE, Gonye GE, Schwaber JS, Doyle III FJ: Importance of input perturbations and stochastic gene expression in the reverse engineering of genetic regulatory networks: insights from an identifiability analysis of an in silico network. Genome Res 2003, 13(11):2396-2405.
• [26]Lee C, Etchegaray JP, Cagampang FR, Loudon ASI, Reppert SM: Posttranslational mechanisms regulate the mammalian circadian clock. Cell 2001, 107(7):855-867.
• [27]Wachter A, Biegler LT: On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math Program 2006, 106:25-57.
• [28]Andersson J, Åkesson J, Diehl M: CasADi – A symbolic package for automatic differentiation and optimal control. In Recent Advances in Algorithmic Differentiation, Lecture Notes in Computational Science and Engineering. Edited by Forth S, Hovland P, Phipps E, Utke J, Walther A. Berlin: Springer; 2012.
• [29]Hindmarsh AC, Brown PN, Grant KE, Lee SL, Serban R, Shumaker DANE, Woodward CS: SUNDIALS: Suite of nonlinear and differential / algebraic equation solvers. ACM T Math Softw 2005, 31(3):363-396.