SMBioNet homepage
Selection of Models of Biological Networks

Adrien RICHARD, Jean-Paul COMET, Gilles BERNOT


arichard@lami.univ-evry.fr, comet@lami.univ-evry.fr

SMBioNet is a software, implemented in Java, which allows us to make an exhaustive qualitative modeling of the dynamics of biological networks using the generalized logical analysis of R. Thomas and coworkers, the feedback circuit functionality and the model checking (CTL/NuSMV).




Introduction Download Example Releted works



Overview



Using an extension of the so-called generalized logical formalism of R. Thomas and coworkers, SMBioNet allows the selection of all models (values of parameters) associated to a given biological regulatory network which lead to a qualitative dynamics (a state transition graph) satisfying the biological knowledges (or hypothesis) on the behavior of studied system. These knowledges on the behavior, called temporal properties, can be mainly expressed as:
  • a list of steady states (regular or singular),
  • a list of functional feedback circuits,
  • a CTL formula (checked with the NuSMV model checker).
R. Thomas' formalism extracts the essential qualitative features of the dynamics described by piecewise-linear ordinary differential equations systems often use to model biological regulatory networks (for example, all the steady states of this continuous modeling can be detected in the corresponding discrete dynamics). Contrary to the quantitative approach, this discrete description used logical parameters wich can take a finite number of possibles values. Thus, even if the values of logical parameters are often unknown (as for the kinetic ones), SMBioNet allows a restriction of the number of parameterizations (models) to consider according to experimental knowlegdes of (or hypothesis on) its behavior:

SMBioNet generates all the possible model of a given network and selects those leading to a qualitative behavior coherent with the specified temporal properties (steady states, functional circuits, CTL formula...).

If not any models are selected, it means that the given biological regulatory network is not compatible with the specified properties. One of the two kind of input (static/dynamic) has to be called into question. Otherwise, the set coherent models can be used to formally test hypothesis (is there, among the selected model, a model satisfying our hypothesis?) make predictions, built new experiments in the aim to refine the modeling (to decrease the number of coherent models)...





SMBioNet DownLoad



Download:



Example of modeling with SMBioNet



The context of this example is given in :

J. Guespin, G. Bernot, J.-P. Comet, A. Mérieau, A. Richard, C. Hulen and B. Polack. Epigenesis and dynamic similarity in two regulatory networks in Pseudomonas aeruginosa. Acta Biotheoretica, in press. (2004).

The following regulatory networks are possible simplified networks controlling the mucoidy and the cytotoxicity of Pseudomonas aeruginosa :

networks 1

networks 2

Acquisition of both phenotype (mucoidy and cytotoxicity) may be due to an epigenetic modification. This epigenetic hypothesis can be translated into the following CTL formula :

((x=0)->AG(x‹2)) & (x=2)->AX(AF(x=2))

When (x=2) the bacteria products mucus and is cytotoxic. When (x=0) it has not both phenotypes. Thus the formula expresses that :
  • if the bacteria has not the phenotypes, in all futures (AG) it has not the phenotypes.
  • if the bacteria has the phenotypes, in all strict futures there is a time (AX AF) where it has the phenotypes.
SMBioNet shows that for both networks 4 models (values of logical parameters) fulfill the epigenetic hypothesis (to obtain this result, we use the Snoussi constraints and we set Kx,{}=0 and Ky,{}=0).

Models satisfying the epigenetic hypothesis

Parameter table of networks 1 Parameter table of networks 2
xy KxKy
00Kx,{y}Ky,{}
01Kx,{}Ky,{}
10Kx,{y}Ky,{x}
11Kx,{}Ky,{x}
20Kx,{xy}Ky,{x}
21Kx,{x}Ky,{x}
xy KxKy
00Kx,{y}Ky,{}
01Kx,{}Ky,{}
10Kx,{xy}Ky,{}
11Kx,{x}Ky,{}
20Kx,{xy}Ky,{x}
21Kx,{x}Ky,{x}

Models of network 1


Models of network 2

Model 1
Model 1
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
002200
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
001200
State graph of model 1 State graph of model 1


Model 2


Model 2
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
012200
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
002200
State graph of model 2 State graph of model 2


Model 3


Model 3
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
002201
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
001201
State graph of model 3 State graph of model 3


Model 4


Model 4
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
012201
Kx,{}Kx,{y}Kx,{x}Kx,{xy}Ky,{}Ky,{x}
002201
State graph of model 4 State graph of model 4


Among this 8 models, SMBioNet automatically shows (with the observability constraints) that 2 models, the model 4 of network 1 and the model 3 of networks 2, are such that all interactions of the networks are active : the two corresponding state graphs cannot be associated to simplest networks (with less interactions).



SMBioNet Related Works



  1. G. Bernot, J.-P. Comet, A. Richard, J. Guespin, "A Fruitful Application of Formal Methods to Biological Regulatory Networks: Extending Thomas' Asynchronous Logical Approach with Temporal Logic", J.T.B., 229(3):339-347, 2004.

  2. J. Guespin, G. Bernot, J.-P. Comet, A. Mérieau, A. Richard, C. Hulen and B. Polack. Epigenesis and dynamic similarity in two regulatory networks in Pseudomonas aeruginosa. Acta Biotheoretica. 52 (4): 379-390, 2004.

    The result of the modeling, obtained with SMBioNet and not given in this paper, is presented in the example of modeling page.

  3. G. Bernot, J. Guespin-Michel, J.-P. Comet, P. Amar, A. Zemirline, F. Delaplace, P. Ballet and A. Richard. Modelling, observability and experiment: a case study. In Proc. of the Dieppe Spring school on Modelling and simulation of biological processes in the context of genomics (eds. P. Amar, F. Képés, V. Norris and P. Tracqui) pp. 49-55, Publisher Frontier group, ISBN : 2 84704 036 6, 2003.






Last update : 07/02/2005