Skip to main content

Modelling genetic regulatory networks from specified behaviours


Modelling and simulation are needed to understand genetic regulatory networks. But parameters of the models are usually difficult to determine. To deal with this problem we propose a methodology in which the qualitative approach developed by R. Thomas [1] is used. The parameters of the model, which are related to the kinetic parameters of a differential description, may be unknown. We translate the set of possible models into one Symbolic Transition System, and the known behaviours into temporal logic formulas; then the constraints on the parameters corresponding to all models having the specified behaviours can be determined.


In the asynchronous discrete modelling of regulatory networks [1], each variable x represents the concentration of a constituent of the network. In each state, the value of x is an integer bounded by the number of variables that x can regulate. Each state and each variable is associated with a parameter that has an integer value. This parameter is the value toward which the variable tends in the associated state.


Pseudomonas aeruginosa are bacteria that secrete mucus in lungs affected by cystic fibrosis, but not in common environment. As it increases respiratory deficiency, this is a major cause of mortality in this disease. The regulatory network proposed in [2], contains the protein AlgU, and an inhibitor complex anti-AlgU (Fig. 1A). Bacteriophage lambda is a virus whose DNA can integrate into bacterial chromosome. After infection, most of the bacteria display a lytic response, but some display a lysogenic response, i.e. survive and carry lambda genome, becoming immune to infection. The graph of interactions described in [3] involves four genes, cI, cro, cII and N (Fig. 1B). The lytic (resp. lysogenic) response leads to the states where cro (resp. cI) is fully expressed. In these two cases the parameters are unknown.

Figure 1
figure 1

Graphs of interactions: A. Mucus production in P. aeruginosa. B. Immunity control in lambda-phage. Each arrow indicates an interaction; the sign indicates a positive or negative effect, and the integer is the threshold of the interaction; when the regulator is greater than or equal to the threshold, the interaction is effective.


The set of all discrete models associated with a graph of interactions are translated into a Symbolic Transition System. Then we apply symbolic execution techniques [4], to construct a tree of states sequences, such that each path is associated with the constraint on parameters necessary to its existence: this constraint is called path condition. To search a specific path in the symbolic execution tree, we have adapted model-checking techniques for Linear Temporal Logic (LTL): all paths verifying the LTL formula are selected, and the disjunction of the associated path conditions is synthesised. The resulting constraint represents all parameters compatible with the specified behaviour.


It has been observed that mucoid P. aeruginosa can continue to produce mucus isolated from infected lungs. It is commonly thought that the mucoidy is due to a mutation which cancels the inhibition of algU gene; an alternative hypothesis is that it is an epigenetic modification, occurring without mutation [2]. With the method described here it is possible to find the constraints such that the resulting models have two stable behaviours, one mucoid and one non-mucoid: the 8 selected models are compatible with the epigenetic hypothesis. In the case of lambda-phage, there are 2156 coherent models with pathways from initial state to lysis and to lysogeny. But in all these models, there is a common path to lysis, and at least one of two precise paths to lysogeny.


Because of the partial knowledge of systems, even with a qualitative formalism, different models can fit with experimental results. Our method allows manipulating not only one model, but a set of coherent models. Then we can efficiently respond to two kinds of questions: is there any selected model coherent with a hypothetic behaviour (as the epigenetic modification in P. aeruginosa)? Are there common behaviours in selected models (as pathways to lysis or lysogeny in lambda-phage)? Moreover, by keeping this set, new experimental results can be added incrementally to restrict and refine the models.


  1. Thomas R, D'Ari R: Biological feedback. 1990, Boca Raton: CRC Press

    Google Scholar 

  2. Guespin-Michel JF, Bernot G, Comet JP, Merieau A, Richard A, Hulen C, Polack B: Epigenesis and dynamic similarity in two regulatory networks in Pseudomonas aeruginosa. Acta Biotheor. 2004, 52 (4): 379-390. 10.1023/B:ACBI.0000046604.18092.a7

    Article  PubMed  Google Scholar 

  3. Thieffry D, Thomas R: Dynamical behaviour of biological regulatory networks–II. Immunity control in bacteriophage lambda. Bull Math Biol. 1995, 57 (2): 277-297.

    PubMed  CAS  Google Scholar 

  4. Gaston C, Le Gall P, Rapin N, Touil A: Symbolic execution techniques for test purpose definition. TestCom, Lecture Notes in Computer Science. 2006, 3964: 1-18. 10.1007/11754008_1.

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations


Corresponding author

Correspondence to Daniel Mateus.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mateus, D., Comet, JP., Gallois, JP. et al. Modelling genetic regulatory networks from specified behaviours. BMC Syst Biol 1 (Suppl 1), S4 (2007).

Download citation

  • Published:

  • DOI: