Issue No.01 - January-February (2011 vol.8)
pp: 166-181
Delphine Ropers , INRIA Grenoble-Rhône-Alpes, Saint Ismier Cedex
Valentina Baldazzi , INRIA Grenoble-Rhône-Alpes, Saint Ismier Cedex
Hidde de Jong , INRIA Grenoble-Rhône-Alpes, Saint Ismier Cedex
The adaptation of the bacterium Escherichia coli to carbon starvation is controlled by a large network of biochemical reactions involving genes, mRNAs, proteins, and signalling molecules. The dynamics of these networks is difficult to analyze, notably due to a lack of quantitative information on parameter values. To overcome these limitations, model reduction approaches based on quasi-steady-state (QSS) and piecewise-linear (PL) approximations have been proposed, resulting in models that are easier to handle mathematically and computationally. These approximations are not supposed to affect the capability of the model to account for essential dynamical properties of the system, but the validity of this assumption has not been systematically tested. In this paper, we carry out such a study by evaluating a large and complex PL model of the carbon starvation response in E. coli using an ensemble approach. The results show that, in comparison with conventional nonlinear models, the PL approximations generally preserve the dynamics of the carbon starvation response network, although with some deviations concerning notably the quantitative precision of the model predictions. This encourages the application of PL models to the qualitative analysis of bacterial regulatory networks, in situations where the reference time scale is that of protein synthesis and degradation.
Bacterial regulatory networks, piecewise-linear differential equations, ordinary differential equations, model reduction, Escherichia coli, stress response, quasi-steady-state approximation.
Delphine Ropers, Valentina Baldazzi, Hidde de Jong, "Model Reduction Using Piecewise-Linear Approximations Preserves Dynamic Properties of the Carbon Starvation Response in Escherichia coli", IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol.8, no. 1, pp. 166-181, January-February 2011, doi:10.1109/TCBB.2009.49
[1] R. Gutierrez-Ríos, J. Freyre-Gonzalez, O. Resendis, J. Collado Vides, M. Saier, and G. Gosset, "Identification of Regulatory Network Topological Units Coordinating the Genome-Wide Transcriptional Response to Glucose in Escherichia coli," BMC Microbiology, vol. 7, p. 53, 2007.
[2] T. Hardiman, K. Lemuth, M. Keller, M. Reuss, and M. Siemann-Herzberg, "Topology of the Global Regulatory Network of Carbon Limitation in Escherichia coli," J. Biotechnology, vol. 132, pp. 359-374, 2007.
[3] H. de Jong and D. Ropers, "Strategies for Dealing with Incomplete Information in the Modeling of Molecular Interaction Networks," Briefings in Bioinformatics, vol. 7, no. 4, pp. 354-363, 2006.
[4] J. Heijnen, "Approximative Kinetic Formats Used in Metabolic Network Modeling," Biotechnology and Bioeng., vol. 91, no. 5, pp. 534-545, 2005.
[5] R. Heinrich and S. Schuster, The Regulation of Cellular Systems. Chapman & Hall, 1996.
[6] T. Millat, E. Bullinger, J. Rohwer, and O. Wolkenhauer, "Approximations and Their Consequences for Dynamic Modelling of Signal Transduction Pathways," Math. Biosciences, vol. 207, no. 1, pp. 40-57, 2007.
[7] M. Roussel and S. Fraser, "Invariant Manifold Methods for Metabolic Model Reduction," Chaos, vol. 11, no. 1, pp. 196-206, 2001.
[8] M. Okino and M. Mavrovouniotis, "Simplification of Mathematical Models of Chemical Reaction Systems," Chemical Rev., vol. 98, no. 2, pp. 391-408, 1998.
[9] J. Papin, J. Stelling, N. Price, S. Klamt, S. Schuster, and B. Palsson, "Comparison of Network-Based Pathway Analysis Methods," Trends in Biotechnology, vol. 22, no. 8, pp. 400-405, 2004.
[10] E. Pecou, "Splitting the Dynamics of Large Biochemical Interaction Networks," J. Theoretical Biology, vol. 232, no. 3, pp. 375-384, 2005.
[11] M. Savageau, "Design Principles for Elementary Gene Circuits: Elements, Methods, and Examples," Chaos, vol. 11, no. 1, pp. 142-159, 2001.
[12] R. Thomas and M. Kaufman, "Multistationarity, the Basis of Cell Differentiation and Memory: II. Logical Analysis of Regulatory Networks in Terms of Feedback Circuits," Chaos, vol. 11, no. 1, pp. 180-195, 2001.
[13] D. Ropers, H. de Jong, M. Page, D. Schneider, and J. Geiselmann, "Qualitative Simulation of the Carbon Starvation Response in Escherichia coli," Biosystems, vol. 84, no. 2, pp. 124-152, 2006.
[14] L. Glass and S. Kauffman, "The Logical Analysis of Continuous Non-Linear Biochemical Control Networks," J. Theoretical Biology, vol. 39, no. 1, pp. 103-129, 1973.
[15] G. Batt, H. de Jong, M. Page, and J. Geiselmann, "Symbolic Reachability Analysis of Genetic Regulatory Networks Using Discrete Abstractions," Automatica, vol. 44, no. 4, pp. 982-989, 2008.
[16] H. de Jong, J.-L. Gouzé, C. Hernandez, M. Page, T. Sari, and J. Geiselmann, "Qualitative Simulation of Genetic Regulatory Networks Using Piecewise-Linear Models," Bull. Math. Biology, vol. 66, no. 2, pp. 301-340, 2004.
[17] A.U. Viretta and M. Fussenegger, "Modeling the Quorum Sensing Regulatory Network of Human-Pathogenic Pseudomonas aeruginosa," Biotechnology Progress, vol. 20, no. 3, pp. 670-678, 2004.
[18] H. de Jong, J. Geiselmann, G. Batt, C. Hernandez, and M. Page, "Qualitative Simulation of the Initiation of Sporulation in B. subtilis," Bull. Math. Biology, vol. 66, no. 2, pp. 261-299, 2004.
[19] J.-A. Sepulchre, S. Reverchon, and W. Nasser, "Modeling the Onset of Virulence in a Pectinolytic Bacterium," J. Theoretical Biology, vol. 44, no. 2, pp. 239-257, 2007.
[20] M. Chaves, E. Sontag, and R. Albert, "Methods of Robustness Analysis for Boolean Models of Gene Control Networks," IET Systems Biology, vol. 153, no. 4, pp. 154-167, 2006.
[21] A. Halász, V. Kumar, M. Imielinski, C. Belta, O. Sokolsky, S. Pathak, and H. Rubin, "Analysis of Lactose Metabolism in E. coli Using Reachability Analysis of Hybrid Systems," IET Systems Biology, vol. 1, no. 2, pp. 130-148, 2007.
[22] G. von Dassow, E. Meir, E. Munro, and G. Odell, "The Segment Polarity Network is a Robust Developmental Module," Nature, vol. 406, no. 6792, pp. 188-192, 2000.
[23] J. Stelling, U. Sauer, Z. Szallasi, F. Doyle, and J. Doyle, "Robustness of Cellular Functions," Cell, vol. 118, no. 6, pp. 675-686, 2004.
[24] T. Erneux and A. Goldbeter, "Rescue of the Quasi-Steady-State Approximation in a Model for Oscillations in an Enzymatic Cascade," SIAM J. Applied Math., vol. 67, no. 2, pp. 305-320, 2007.
[25] R. Krüger and R. Heinrich, "Model Reduction and Analysis of Robustness for the Wnt/$\beta$ -Catenin Signal Transduction Pathway," Genome Informatics, vol. 15, no. 1, pp. 138-148, 2004.
[26] O. Radulescu, A. Gorban, A. Zinovyev, and A. Lilienbaum, "Robust Simplifications of Multiscale Biochemical Networks," BMC Systems Biology, vol. 2, p. 86, 2008.
[27] A. Ciliberto, F. Capuani, and J. Tyson, "Modeling Networks of Coupled Enzymatic Reactions Using the Total Quasi-Steady State Approximation," PLoS Computational Biology, vol. 3, no. 3, p. e45, 2007.
[28] L. Segel and M. Slemrod, "The Quasi-Steady-State Assumption: A Case-Study in Perturbation," SIAM Rev., vol. 31, no. 3, pp. 446-477, 1989.
[29] S. Fraser, "The Steady State and Equilibrium Approximations: A Geometrical Picture," J. Chemical Physics, vol. 88, no. 8, pp. 4732-4738, 1988.
[30] T. Mestl, E. Plahte, and S. Omholt, "A Mathematical Framework for Describing and Analysing Gene Regulatory Networks," J. Theoretical Biology, vol. 176, no. 2, pp. 291-300, 1995.
[31] M. Ptashne, A Genetic Switch: Phage $\lambda$ Revisited, third ed., Cold Spring Harbor Laboratory Press, 2004.
[32] D. Battogtokh, D. Asch, M. Case, J. Arnold, and H. Schüttler, "An Ensemble Method for Identifying Regulatory Circuits with Special Reference to the qa Gene Cluster of Neurospora crassa," Proc. Nat'l Academy of Sciences USA, vol. 99, no. 26, pp. 16904-16909, 2002.
[33] L. Kuepfer, M. Peter, U. Sauer, and J. Stelling, "Ensemble Modeling for Analysis of Cell Signaling Dynamics," Nature Biotechnology, vol. 25, no. 9, pp. 1001-1006, 2007.
[34] R. Alves and M. Savageau, "Comparing Systemic Properties of Ensembles of Biological Networks by Graphical and Statistical Methods," Bioinformatics, vol. 16, no. 6, pp. 527-533, 2000.
[35] R. Alves and M. Savageau, "Systemic Properties of Ensembles of Metabolic Networks: Application of Graphical and Statistical Methods to Simple Unbranched Pathways," Bioinformatics, vol. 16, no. 6, pp. 534-547, 2000.
[36] H. Kwakernaak and R. Sivan, Modern Signals and Systems. Prentice-Hall, 1991.
[37] R. Casey, H. de Jong, and J.-L. Gouzé, "Piecewise-Linear Models of Genetic Regulatory Networks: Equilibria and Their Stability," J. Math. Biology, vol. 52, no. 1, pp. 27-56, 2006.
[38] H. de Jong and M. Page, "Search for Steady States of Piecewise-Linear Differential Equation Models of Genetic Regulatory Networks," IEEE/ACM Trans. Computational Biology and Bioinformatics, vol. 5, no. 2, pp. 208-222, Apr.-June 2008.
[39] K. Bettenbrock, S. Fischer, A. Kremling, K. Jahreis, T. Sauter, and E. Gilles, "A Quantitative Approach to Catabolite Repression in Escherichia coli," J. Biological Chemistry, vol. 281, no. 5, pp. 2578-2584, 2006.
[40] J. Rohwer, N. Meadow, S. Roseman, H. Westerhoff, and P. Postma, "Understanding Glucose Transport by the Bacterial Phosphoenolpyruvate: Glycose Phosphotransferase System on the Basis of Kinetic Measurements in vitro," J. Biological Chemistry, vol. 275, no. 45, pp. 34909-34921, 2000.
[41] C. Chassagnole, N. Noisommit-Rizzi, J. Schmid, K. Mauch, and M. Reuss, "Dynamic Modeling of the Central Carbon Metabolism of Escherichia coli," Biotechnology and Bioeng., vol. 79, no. 1, pp. 53-73, 2002.
[42] P. Dennis, M. Ehrenberg, and H. Bremer, "Control of rRNA Synthesis in Escherichia coli: A Systems Biology Approach," Microbiology and Molecular Biology Rev., vol. 68, no. 4, pp. 639-668, 2004.
[43] G. Hatfield and C. Benham, "DNA Topology-Mediated Control of Global Gene Expression in Escherichia coli," Ann. Rev. Genetics, vol. 36, pp. 175-203, 2002.
[44] R. Hengge-Aronis, "The General Stress Response in Escherichia coli," Bacterial Stress Responses, G. Storz and R. Hengge-Aronis, eds., ASM Press, pp. 161-177, 2000.
[45] M.H. Saier, T.M. Ramseier, and J. Reizer, "Regulation of Carbon Utilization," Escherichia coli and Salmonella: Cellular and Molecular Biology, F. Neidhardt, R. Curtiss III, J. Ingraham, E. Lin, K. Low, B. Magasanik, W. Reznikoff, M. Riley, M. Schaechter, and H. Umbarger, eds., ASM Press, pp. 1325-1343, 1996.
[46] J. Barthelmes, C. Ebeling, A. Chang, I. Schomburg, and D. Schomburg, "BRENDA, AMENDA and FRENDA: The Enzyme Information System in 2007," Nucleic Acids Research, vol. 35, pp. D511-D514, 2007.
[47] P. Reddy, A. Peterkofsky, and K. McKenney, "Hyperexpression and Purification of Escherichia coli Adenylate Cyclase Using a Vector Designed for Expression of Lethal Gene Products," Nucleic Acids Research, vol. 17, pp. 10473-10488, 1989.
[48] J. Yang and W. Epstein, "Purification and Characterization of Adenylate Cyclase from Escherichia coli K12," J. Biological Chemistry, vol. 258, no. 6, pp. 3750-3758, 1983.
[49] W. Liebermeister and E. Klipp, "Biochemical Networks with Uncertain Parameters," IET Systems Biology, vol. 152, no. 3, pp. 97-107, 2005.
[50] A.G. Gonzalez, A. Naldi, L. Sánchez, D. Thieffry, and C. Chaouiya, "GINsim: a Software Suit for the Qualitative Modelling, Simulation and Analysis of Regulatory Networks," Biosystems, vol. 84, no. 2, pp. 91-100, 2006.
[51] R. Thomas and R. d'Ari, Biological Feedback. CRC Press, 1990.
[52] E. Plahte and S. Kjóglum, "Analysis and Generic Properties of Gene Regulatory Networks with Graded Response Functions," Physica D, vol. 201, no. 1, pp. 150-176, 2005.
[53] M. Chaves, R. Albert, and E. Sontag, "Robustness and Fragility of Boolean Models for Genetic Regulatory Networks," J. Theoretical Biology, vol. 235, no. 3, pp. 431-449, 2005.
[54] M. Davidich and S. Bornholdt, "The Transition from Differential Equations to Boolean Networks: A Case Study in Simplifying a Regulatory Network Model," J. Theoretical Biology, vol. 255, pp. 269-277, 2008.
[55] A. Andronov, E. Leontovich, I. Gordon, and A. Maier, Qualitative Theory of Second-Order Dynamic Systems. Wiley & Sons, 1973.
[56] C. Belta and L. Habets, "Controlling a Class of Nonlinear Systems on Rectangles," IEEE Trans. Automatic Control, vol. 51, no. 11, pp. 1749-1759, Nov. 2006.
[57] M. Musters, H. de Jong, P. van den Bosch, and N. van Riel, "Qualitative Analysis of Nonlinear Biochemical Networks with Piecewise-Affine Functions," Proc. Int'l Workshop Hybrid Systems: Computation and Control (HSCC '07), A. Bemporad, A. Bicchi, and G. Buttazzo, eds., pp. 727-730, 2007.
[58] M. Quach, N. Brunel, and F. d'Alché Buc, "Estimating Parameters and Hidden Variables in Non-Linear State-Space Models Based on ODEs for Biological Networks Inference," Bioinformatics, vol. 23, no. 23, pp. 3209-3216, 2007.
[59] J. Zwolak, J. Tyson, and L. Watson, "Parameter Estimation for a Mathematical Model of the Cell Cycle in Frog Eggs," J. Computational Biology, vol. 12, no. 1, pp. 48-63, 2005.
[60] N. van Riel and E. Sontag, "Parameter Estimation in Models Combining Signal Transduction and Metabolic Pathways: The Dependent Input Approach," IET Systems Biology, vol. 153, no. 4, pp. 263-274, 2006.
[61] T. Gardner, D. di Bernardo, D. Lorenz, and J. Collins, "Inferring Genetic Networks and Identifying Compound Mode of Action via Expression Profiling," Science, vol. 301, no. 5629, pp. 102-105, 2003.
[62] C. Moles, P. Mendes, and J. Banga, "Parameter Estimation in Biochemical Pathways: A Comparison of Global Optimization Methods," Genome Research, vol. 13, no. 11, pp. 2467-2474, 2003.
[63] M. Vilela, I. Chou, S. Vinga, A. Vasconcelos, E. Voit, and J. Almeida, "Parameter Optimization in S-System Models," BMC Systems Biology, vol. 2, p. 35, 2008.
[64] R. Porreca, S. Drulhe, H. de Jong, and G. Ferrari-Trecate, "Structural Identification of Piecewise-Linear Models of Genetic Regulatory Networks," J. Computational Biology, vol. 15, no. 10, pp. 1365-1380, 2008.
[65] R. Westra, G. Hollanders, G. Bex, M. Gyssens, and K. Tuyls, "The Identification of Dynamic Gene-Protein Networks," Proc. First Int'l Workshop Knowledge Discovery and Emergent Complexity in Bioinformatics (KDECB), K. Tuyls, R. Westra, Y. Saeys, and A. Nowé, eds., pp. 157-170, 2007.
[66] M. Ronen, R. Rosenberg, B. Shraiman, and U. Alon, "Assigning Numbers to the Arrows: Parameterizing a Gene Regulation Network by Using Accurate Expression Kinetics," Proc. Nat'l Academy of Sciences USA vol. 99, no. 16, pp. 10555-10560, 2002.