Transforming Boolean models to continuous models: methodology and application to T-cell receptor signaling

BackgroundThe understanding of regulatory and signaling networks has long been a core objective in Systems Biology. Knowledge about these networks is mainly of qualitative nature, which allows the construction of Boolean models, where the state of a component is either 'off' or 'on'. While often able to capture the essential behavior of a network, these models can never reproduce detailed time courses of concentration levels.Nowadays however, experiments yield more and more quantitative data. An obvious question therefore is how qualitative models can be used to explain and predict the outcome of these experiments.ResultsIn this contribution we present a canonical way of transforming Boolean into continuous models, where the use of multivariate polynomial interpolation allows transformation of logic operations into a system of ordinary differential equations (ODE). The method is standardized and can readily be applied to large networks. Other, more limited approaches to this task are briefly reviewed and compared. Moreover, we discuss and generalize existing theoretical results on the relation between Boolean and continuous models. As a test case a logical model is transformed into an extensive continuous ODE model describing the activation of T-cells. We discuss how parameters for this model can be determined such that quantitative experimental results are explained and predicted, including time-courses for multiple ligand concentrations and binding affinities of different ligands. This shows that from the continuous model we may obtain biological insights not evident from the discrete one.ConclusionThe presented approach will facilitate the interaction between modeling and experiments. Moreover, it provides a straightforward way to apply quantitative analysis methods to qualitatively described systems.

[1]  D. Thieffry,et al.  Dynamical behaviour of biological regulatory networks—I. Biological role of feedback loops and practical use of the concept of the loop-characteristic state , 1995 .

[2]  H. Othmer,et al.  The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. , 2003, Journal of theoretical biology.

[3]  Luhua Lai,et al.  Robustness and modular design of the Drosophila segment polarity network , 2006, Molecular systems biology.

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

[5]  Ronald N Germain,et al.  Modeling T Cell Antigen Discrimination Based on Feedback Control of Digital ERK Responses , 2005, PLoS biology.

[6]  Barbara Hausmann,et al.  A motif in the αβ T-cell receptor controls positive selection by modulating ERK activity , 2000, Nature.

[7]  Steffen Klamt,et al.  A methodology for the structural and functional analysis of signaling and regulatory networks , 2006, BMC Bioinformatics.

[8]  Thomas Mestl,et al.  A methodological basis for description and analysis of systems with complex switch-like interactions , 1998, Journal of mathematical biology.

[9]  Denis Thieffry,et al.  Genetic control of flower morphogenesis in Arabidopsis thaliana: a logical analysis , 1999, Bioinform..

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

[11]  Madalena Chaves,et al.  Robustness and fragility of Boolean models for genetic regulatory networks. , 2005, Journal of theoretical biology.

[12]  Chi-Ying F. Huang,et al.  Ultrasensitivity in the mitogen-activated protein kinase cascade. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Barbara Hausmann,et al.  Signaling Life and Death in the Thymus: Timing Is Everything , 2003, Science.

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

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

[16]  A. Hill,et al.  The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves , 1910 .

[17]  Barbara Hausmann,et al.  Thymic selection threshold defined by compartmentalization of Ras/MAPK signalling , 2006, Nature.

[18]  S. Bornholdt,et al.  Boolean Network Model Predicts Cell Cycle Sequence of Fission Yeast , 2007, PloS one.

[19]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[20]  Lindsay B. Nicholson,et al.  Reduced Self-Reactivity of an Autoreactive T Cell After Activation with Cross-reactive Non–Self-Ligand , 2002, The Journal of experimental medicine.

[21]  E. Palmer,et al.  A motif in the alphabeta T-cell receptor controls positive selection by modulating ERK activity. , 2000, Nature.

[22]  Ioannis Xenarios,et al.  A method for the generation of standardized qualitative dynamical systems of regulatory networks , 2005, Theoretical Biology and Medical Modelling.

[23]  Mats Jirstrand,et al.  Systems biology Systems Biology Toolbox for MATLAB : a computational platform for research in systems biology , 2006 .

[24]  T. Sauer,et al.  On the history of multivariate polynomial interpolation , 2000 .

[25]  Steffen Klamt,et al.  Structural and functional analysis of cellular networks with CellNetAnalyzer , 2007, BMC Systems Biology.

[26]  Melissa L. Kemp,et al.  Quantitative Network Signal Combinations Downstream of TCR Activation Can Predict IL-2 Production Response1 , 2007, The Journal of Immunology.

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

[28]  Steffen Klamt,et al.  A Logical Model Provides Insights into T Cell Receptor Signaling , 2007, PLoS Comput. Biol..