+ Site Statistics
+ Search Articles
+ Subscribe to Site Feeds
Most Shared
PDF Full Text
+ PDF Full Text
Request PDF Full Text
+ Follow Us
Follow on Facebook
Follow on Twitter
Follow on LinkedIn
+ Translate
+ Recently Requested

Continuous time Boolean modeling for biological signaling: application of Gillespie algorithm

Continuous time Boolean modeling for biological signaling: application of Gillespie algorithm

Bmc Systems Biology 6: 116

Mathematical modeling is used as a Systems Biology tool to answer biological questions, and more precisely, to validate a network that describes biological observations and predict the effect of perturbations. This article presents an algorithm for modeling biological networks in a discrete framework with continuous time. There exist two major types of mathematical modeling approaches: (1) quantitative modeling, representing various chemical species concentrations by real numbers, mainly based on differential equations and chemical kinetics formalism; (2) and qualitative modeling, representing chemical species concentrations or activities by a finite set of discrete values. Both approaches answer particular (and often different) biological questions. Qualitative modeling approach permits a simple and less detailed description of the biological systems, efficiently describes stable state identification but remains inconvenient in describing the transient kinetics leading to these states. In this context, time is represented by discrete steps. Quantitative modeling, on the other hand, can describe more accurately the dynamical behavior of biological processes as it follows the evolution of concentration or activities of chemical species as a function of time, but requires an important amount of information on the parameters difficult to find in the literature. Here, we propose a modeling framework based on a qualitative approach that is intrinsically continuous in time. The algorithm presented in this article fills the gap between qualitative and quantitative modeling. It is based on continuous time Markov process applied on a Boolean state space. In order to describe the temporal evolution of the biological process we wish to model, we explicitly specify the transition rates for each node. For that purpose, we built a language that can be seen as a generalization of Boolean equations. Mathematically, this approach can be translated in a set of ordinary differential equations on probability distributions. We developed a C++ software, MaBoSS, that is able to simulate such a system by applying Kinetic Monte-Carlo (or Gillespie algorithm) on the Boolean state space. This software, parallelized and optimized, computes the temporal evolution of probability distributions and estimates stationary distributions. Applications of the Boolean Kinetic Monte-Carlo are demonstrated for three qualitative models: a toy model, a published model of p53/Mdm2 interaction and a published model of the mammalian cell cycle. Our approach allows to describe kinetic phenomena which were difficult to handle in the original models. In particular, transient effects are represented by time dependent probability distributions, interpretable in terms of cell populations.

(PDF emailed within 0-6 h: $19.90)

Accession: 036558332

Download citation: RISBibTeXText

PMID: 22932419

DOI: 10.1186/1752-0509-6-116

Related references

Polynomial-time algorithm for controllability test of a class of boolean biological networks. Eurasip Journal on Bioinformatics and Systems Biology 2010, 2010

An adjustable aperiodic model class of genomic interactions using continuous time Boolean networks (Boolean delay equations). Chaos 13(4): 1167-1174, 2003

Modeling Langmuir isotherms with the Gillespie stochastic algorithm. Journal of Chromatography. a 1380: 81-87, 2015

Modeling ant battles by means of a diffusion-limited Gillespie algorithm. Theoretical Biology Forum 107(1-2): 57-76, 2015

Transforming Boolean models to continuous models: methodology and application to T-cell receptor signaling. Bmc Systems Biology 3: 98, 2009

An implementation of the Gillespie algorithm for RNA kinetics with logarithmic time update. Nucleic Acids Research 43(12): 5708-5715, 2015

Temporal Gillespie Algorithm: Fast Simulation of Contagion Processes on Time-Varying Networks. Plos Computational Biology 11(10): E1004579, 2016

Simulating the cellular passive transport of glucose using a time-dependent extension of Gillespie algorithm for stochastic pi-calculus. International Journal of Data Mining and Bioinformatics 1(4): 315-336, 2008

A New Type of Discrete Self-Decomposability and Its Application to Continuous-Time Markov Processes for Modeling Count Data Time Series. Stochastic Models 19(2): 235-254, 2003

A Boolean approach to linear prediction for signaling network modeling. Plos One 5(9), 2011

A Boolean network inference from time-series gene expression data using a genetic algorithm. Bioinformatics 34(17): I927-I933, 2018

Modeling mechanistic biological networks: an advanced Boolean approach. Bioinformatics 28(4): 557-563, 2012

Boolean modeling of biological regulatory networks: a methodology tutorial. Methods 62(1): 3-12, 2014

Modeling signaling-dependent pluripotency with Boolean logic to predict cell fate transitions. Molecular Systems Biology 14(1): E7952, 2018

Application of Burg algorithm in time-frequency analysis of Doppler blood flow signal based on AR modeling. Sheng Wu Yi Xue Gong Cheng Xue Za Zhi 22(3): 481-485, 2005