Quantification of reachable attractors in asynchronous discrete dynamics

Motivation: Models of discrete concurrent systems often lead to huge and complex state transition graphs that represent their dynamics. This makes difficult to analyse dynamical properties. In particular, for logical models of biological regulatory networks, it is of real interest to study attractors and their reachability from specific initial conditions, i.e. to assess the potential asymptotical behaviours of the system. Beyond the identification of the reachable attractors, we propose to quantify this reachability. Results: Relying on the structure of the state transition graph, we estimate the probability of each attractor reachable from a given initial condition or from a portion of the state space. First, we present a quasi-exact solution with an original algorithm called Firefront, based on the exhaustive exploration of the reachable state space. Then, we introduce an adapted version of Monte Carlo simulation algorithm, termed Avatar, better suited to larger models. Firefront and Avatar methods are validated and compared to other related approaches, using as test cases logical models of synthetic and biological networks. Availability: Both algorithms are implemented as Perl scripts that can be freely downloaded from http://compbio.igc.gulbenkian.pt/nmd/node/59 along with Supplementary Material.

[1]  Charles M. Grinstead,et al.  Introduction to probability , 1999, Statistics for the Behavioural Sciences.

[2]  Giovanni De Micheli,et al.  Synchronous versus asynchronous modeling of gene regulatory networks , 2008, Bioinform..

[3]  Réka Albert,et al.  An effective network reduction approach to find the dynamical repertoire of discrete dynamic networks. , 2013, Chaos.

[4]  Assieh Saadatpour,et al.  Discrete dynamic modeling of signal transduction networks. , 2012, Methods in molecular biology.

[5]  Thomas A. Henzinger,et al.  Sliding Window Abstraction for Infinite Markov Chains , 2009, CAV.

[6]  D. Thieffry,et al.  Segmenting the fly embryo: logical analysis of the role of the segment polarity cross-regulatory module. , 2008, The International journal of developmental biology.

[7]  Aurélien Naldi,et al.  Decision Diagrams for the Representation and Analysis of Logical Models of Genetic Networks , 2007, CMSB.

[8]  Hans A. Kestler,et al.  BoolNet - an R package for generation, reconstruction and analysis of Boolean networks , 2010, Bioinform..

[9]  Aurélien Naldi,et al.  Diversity and Plasticity of Th Cell Types Predicted from Regulatory Network Modelling , 2010, PLoS Comput. Biol..

[10]  Aurélien Naldi,et al.  Logical modelling of gene regulatory networks with GINsim. , 2012, Methods in molecular biology.

[11]  Aurélien Naldi,et al.  Dynamical analysis of a generic Boolean model for the control of the mammalian cell cycle , 2006, ISMB.

[12]  Leon Glass,et al.  Logical and symbolic analysis of robust biological dynamics. , 2010, Current opinion in genetics & development.

[13]  Bernard Prum Chaînes de Markov et absorption. Application à l’algorithme de Fu en génomique , 2012 .

[14]  Emmanuel Barillot,et al.  Continuous time boolean modeling for biological signaling: application of Gillespie algorithm , 2012, BMC Systems Biology.

[15]  M. Khammash,et al.  The finite state projection algorithm for the solution of the chemical master equation. , 2006, The Journal of chemical physics.

[16]  John N. Tsitsiklis,et al.  Introduction to Probability , 2002 .

[17]  Robert E. Tarjan,et al.  Depth-First Search and Linear Graph Algorithms , 1972, SIAM J. Comput..

[18]  Elizabeth L. Wilmer,et al.  Markov Chains and Mixing Times , 2008 .