Search for Steady States of Piecewise-Linear Differential Equation Models of Genetic Regulatory Networks

The analysis of the attractors of a genetic regulatory network gives a good indication of the possible functional modes of the system. In this paper, we are concerned with the problem of finding all steady states of genetic regulatory networks described by piecewise-linear differential equation (PLDE) models. We show that the problem is NP-hard and translate it into the problem of finding all valuations of a propositional satisfiability (SAT) formula. This allows the use of existing, efficient SAT solvers and has enabled the development of a steady state search module of the computer tool genetic network analyzer (GNA). The practical use of this module is demonstrated by means of the analysis of a number of relatively small bacterial regulatory networks, as well as randomly generated networks of several hundred genes.

[1]  T. Mestl,et al.  A mathematical framework for describing and analysing gene regulatory networks. , 1995, Journal of theoretical biology.

[2]  Jun Gu,et al.  Algorithms for the satisfiability (SAT) problem: A survey , 1996, Satisfiability Problem: Theory and Applications.

[3]  Radu Mateescu,et al.  Validation of qualitative models of genetic regulatory networks by model checking: analysis of the nutritional stress response in Escherichia coli , 2005, ISMB.

[4]  Christopher Edwards,et al.  Sliding mode control : theory and applications , 1998 .

[5]  David S. Johnson,et al.  Computers and Intractability: A Guide to the Theory of NP-Completeness , 1978 .

[6]  R. Thomas,et al.  Boolean formalization of genetic control circuits. , 1973, Journal of theoretical biology.

[7]  J. Stelling,et al.  Robustness of Cellular Functions , 2004, Cell.

[8]  G. Casari,et al.  From molecular networks to qualitative cell behavior , 2005, FEBS letters.

[9]  M. Ptashne.  A genetic switch : phage λ and higher organisms , 1992 .

[10]  Matthew L. Ginsberg,et al.  Generalizing Boolean Satisfiability I: Background and Survey of Existing Work , 2011, J. Artif. Intell. Res..

[11]  Katy C. Kao,et al.  Global Expression Profiling of Acetate-grown Escherichia coli * , 2002, The Journal of Biological Chemistry.

[12]  Vipul Periwal,et al.  System Modeling in Cellular Biology: From Concepts to Nuts and Bolts , 2006 .

[13]  H. Kitano.  Systems Biology: A Brief Overview , 2002, Science.

[14]  Hidde de Jong,et al.  Genetic Network Analyzer: qualitative simulation of genetic regulatory networks , 2003, Bioinform..

[15]  El Houssine Snoussi.  Qualitative dynamics of piecewise-linear differential equations: a discrete mapping approach , 1989 .

[16]  T. Zolezzi,et al.  Differential Inclusions and Sliding Mode Control , 2002 .

[17]  D. Schneider,et al.  Qualitative simulation of the carbon starvation response in Escherichia coli. , 2006, Bio Systems.

[18]  T. Mestl,et al.  Periodic solutions in systems of piecewise- linear differential equations , 1995 .

[19]  D. A. Baxter,et al.  Modeling transcriptional control in gene networks—methods, recent results, and future directions , 2000, Bulletin of mathematical biology.

[20]  R. Edwards.  Analysis of continuous-time switching networks , 2000 .

[21]  L. Glass,et al.  Chaos in high-dimensional neural and gene networks , 1996 .

[22]  E. Davidson,et al.  Modeling transcriptional regulatory networks. , 2002, BioEssays : news and reviews in molecular, cellular and developmental biology.

[23]  Martin Fussenegger,et al.  Modeling the Quorum Sensing Regulatory Network of Human‐Pathogenic Pseudomonas aeruginosa , 2004, Biotechnology progress.

[24]  John Taylor Stallings,et al.  The Search For Satisfaction , 1935 .

[25]  Toby Walsh,et al.  The Search for Satisfaction , 1999 .

[26]  J. Gouzé,et al.  A class of piecewise linear differential equations arising in biological models , 2002 .

[27]  S. Kauffman.  Metabolic stability and epigenesis in randomly constructed genetic nets. , 1969, Journal of theoretical biology.

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

[29]  H. D. Jong,et al.  Qualitative simulation of the initiation of sporulation in Bacillus subtilis , 2004, Bulletin of mathematical biology.

[30]  Martine Labbé,et al.  Identification of all steady states in large networks by logical analysis , 2003, Bulletin of mathematical biology.

[31]  El Houssine Snoussi,et al.  Logical identification of all steady states: The concept of feedback loop characteristic states , 1993 .

[32]  Chen Yang,et al.  Analysis of Gene Expression in Escherichia coli in Response to Changes of Growth-Limiting Nutrient in Chemostat Cultures , 2004, Applied and Environmental Microbiology.

[33]  Farren J. Isaacs,et al.  Computational studies of gene regulatory networks: in numero molecular biology , 2001, Nature Reviews Genetics.

[34]  Erik Plahte,et al.  Analysis and generic properties of gene regulatory networks with graded response functions , 2005 .

[35]  Andrew W. Murray,et al.  The Ups and Downs of Modeling the Cell Cycle , 2004, Current Biology.

[36]  L. Glass,et al.  Stable oscillations in mathematical models of biological control systems , 1978 .

[37]  Hidde de Jong,et al.  Modeling and Simulation of Genetic Regulatory Systems: A Literature Review , 2002, J. Comput. Biol..

[38]  H. D. Jong,et al.  Qualitative simulation of genetic regulatory networks using piecewise-linear models , 2004, Bulletin of mathematical biology.

[39]  H. D. Jong,et al.  Piecewise-linear Models of Genetic Regulatory Networks: Equilibria and their Stability , 2006, Journal of mathematical biology.

[40]  M. Ptashne.  A Genetic Switch , 1986 .

[41]  J. Sepulchre,et al.  Modeling the onset of virulence in a pectinolytic bacterium. , 2007, Journal of theoretical biology.

[42]  Aleksej F. Filippov,et al.  Differential Equations with Discontinuous Righthand Sides , 1988, Mathematics and Its Applications.