An unbiased sensitivity analysis reveals important parameters controlling periodicity of circadian clock

To assess the importance of model parameters in kinetic models, sensitivity analysis is generally employed to provide key measures. However, it is quite often that no information is available for a significant number of parameters in biochemical models. Therefore, the results of sensitivity analysis that heavily rely on the accuracy of parameters are largely ambiguous. In this study, we propose a computational approach to determine the relative importance of parameters controlling the performance of the circadian clock in Drosophila. While previous attempts to sensitivity analysis largely depend on the knowledge of model parameters which are generally unknown, our study depicts a consistent picture of sensitivity assessment for a large number of parameters, even when the values of these parameters are not available in vivo. The resulting parametric sensitivity analysis suggests that PER/TIM negative loop is critical to maintain the stable periodicity of the circadian clock, which is consistent to the previously experimental and computational findings. Furthermore, our analysis generates a rich hypothesis of important parameters in the circadian clock that can be further tested experimentally. This approach can also be extended to assess the sensitivity of parameters in any biochemical system where a large number of parameters have unknown values. Biotechnol. Bioeng. 2010; 105: 250–259. © 2009 Wiley Periodicals, Inc.

[1]  P. Hardin,et al.  The Circadian Timekeeping System of Drosophila , 2005, Current Biology.

[2]  M. Nitabach Circadian rhythms: Clock coordination , 2005, Nature.

[3]  Paul Smolen,et al.  Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK. , 2004, Biophysical journal.

[4]  D. Kulasiri,et al.  Modelling Circadian Rhythms in Drosophila and Investigation of VRI and PDP1 Feedback Loops Using a New Mathematical Model , 2008 .

[5]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[6]  Jeffrey C. Hall,et al.  An Ultrashort Clock Mutation at the period Locus of Drosophila melanogaster That Reveals Some New Features of the Fly's Circadian System , 1994, Journal of biological rhythms.

[7]  S. Kay,et al.  Time zones: a comparative genetics of circadian clocks , 2001, Nature Reviews Genetics.

[8]  J. Tyson,et al.  A simple model of circadian rhythms based on dimerization and proteolysis of PER and TIM. , 1999, Biophysical journal.

[9]  Michael W. Young,et al.  Rhythmic Expression of timeless: A Basis for Promoting Circadian Cycles in period Gene Autoregulation , 1995, Science.

[10]  William G. Boissonnault Chapter 7 – Review of Systems , 2005 .

[11]  D. A. Baxter,et al.  Modeling Circadian Oscillations with Interlocking Positive and Negative Feedback Loops , 2001, The Journal of Neuroscience.

[12]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[13]  A. Goldbeter,et al.  Modeling the molecular regulatory mechanism of circadian rhythms in Drosophila. , 2000, BioEssays : news and reviews in molecular, cellular and developmental biology.

[14]  D. Allen,et al.  A circadian enhancer mediates PER-dependent mRNA cycling in Drosophila melanogaster , 1997, Molecular and cellular biology.

[15]  A. Goldbeter Computational approaches to cellular rhythms , 2002, Nature.

[16]  G. Hornberger,et al.  Approach to the preliminary analysis of environmental systems , 1981 .

[17]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[18]  J L Hargrove,et al.  The kinetics of mammalian gene expression. , 1991, BioEssays : news and reviews in molecular, cellular and developmental biology.

[19]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[20]  D. Kulasiri,et al.  Modelling of circadian rhythms in Drosophila incorporating the interlocked PER/TIM and VRI/PDP1 feedback loops. , 2007, Journal of theoretical biology.

[21]  Hua Wu,et al.  Parametric sensitivity in chemical systems , 1999 .

[22]  J. Stelling,et al.  Robustness properties of circadian clock architectures. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[23]  M. Stephens Use of the Kolmogorov-Smirnov, Cramer-Von Mises and Related Statistics without Extensive Tables , 1970 .

[24]  Kenzo Hirose,et al.  Intercellular coupling mechanism for synchronized and noise-resistant circadian oscillators. , 2002, Journal of theoretical biology.

[25]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[26]  Erik D Herzog,et al.  Circadian Rhythms: In the Loop at Last , 2003, Science.

[27]  M. W. Young,et al.  Isolation and analysis of six timeless alleles that cause short- or long-period circadian rhythms in Drosophila. , 2000, Genetics.

[28]  J. Dunlap Molecular Bases for Circadian Clocks , 1999, Cell.

[29]  A. Goldbeter,et al.  Toward a detailed computational model for the mammalian circadian clock , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Oliver T. Fackler,et al.  The Nef Protein of Human Immunodeficiency Virus Establishes Superinfection Immunity by a Dual Strategy to Downregulate Cell-Surface CCR5 and CD4 , 2005, Current Biology.

[31]  Paul E. Brown,et al.  Extension of a genetic network model by iterative experimentation and mathematical analysis , 2005, Molecular systems biology.

[32]  Douglas B. Kell,et al.  Non-linear optimization of biochemical pathways: applications to metabolic engineering and parameter estimation , 1998, Bioinform..

[33]  S. Leibler,et al.  Mechanisms of noise-resistance in genetic oscillators , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[34]  D A Rand,et al.  Mapping global sensitivity of cellular network dynamics: sensitivity heat maps and a global summation law , 2008, Journal of The Royal Society Interface.

[35]  S. Samarasinghe,et al.  A Review of Systems Biology Perspective on Genetic Regulatory Networks with Examples , 2008 .

[36]  François Rouyer,et al.  The F-box protein Slimb controls the levels of clock proteins Period and Timeless , 2002, Nature.