Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Bioinformatics**|**v.11; 2010**|**PMC2873544

Formats

Article sections

- Abstract
- Background
- Implementation
- Results and Discussion
- Conclusions
- Availability and requirements
- Authors' contributions
- References

Authors

Related links

BMC Bioinformatics. 2010; 11: 233.

Published online 2010 May 7. doi: 10.1186/1471-2105-11-233

PMCID: PMC2873544

Reviewed by Jan Krumsiek,^{1} Sebastian Pölsterl,^{1} Dominik M Wittmann,^{1,}^{2} and Fabian J Theis^{}^{1,}^{2}

Jan Krumsiek: ed.nehcneum-ztlohmleh@keismurk.naj; Sebastian Pölsterl: ed.uml.ifi.pic@lretsleop; Dominik M Wittmann: ed.nehcneum-ztlohmleh@nnamttiw.kinimod; Fabian J Theis: ed.nehcneum-ztlohmleh@sieht.naibaf

Received 2009 November 23; Accepted 2010 May 7.

Copyright ©2010 Krumsiek et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Phenomenological information about regulatory interactions is frequently available and can be readily converted to Boolean models. Fully quantitative models, on the other hand, provide detailed insights into the precise dynamics of the underlying system. In order to connect discrete and continuous modeling approaches, methods for the conversion of Boolean systems into systems of ordinary differential equations have been developed recently. As biological interaction networks have steadily grown in size and complexity, a fully automated framework for the conversion process is desirable.

We present *Odefy*, a MATLAB- and Octave-compatible toolbox for the automated transformation of Boolean models into systems of ordinary differential equations. Models can be created from sets of Boolean equations or graph representations of Boolean networks. Alternatively, the user can import Boolean models from the CellNetAnalyzer toolbox, GINSim and the PBN toolbox. The Boolean models are transformed to systems of ordinary differential equations by multivariate polynomial interpolation and optional application of sigmoidal Hill functions. Our toolbox contains basic simulation and visualization functionalities for both, the Boolean as well as the continuous models. For further analyses, models can be exported to SQUAD, GNA, MATLAB script files, the SB toolbox, SBML and R script files. Odefy contains a user-friendly graphical user interface for convenient access to the simulation and exporting functionalities. We illustrate the validity of our transformation approach as well as the usage and benefit of the Odefy toolbox for two biological systems: a mutual inhibitory switch known from stem cell differentiation and a regulatory network giving rise to a specific spatial expression pattern at the mid-hindbrain boundary.

Odefy provides an easy-to-use toolbox for the automatic conversion of Boolean models to systems of ordinary differential equations. It can be efficiently connected to a variety of input and output formats for further analysis and investigations. The toolbox is open-source and can be downloaded at http://cmb.helmholtz-muenchen.de/odefy.

The ultimate goal of the increasingly popular systems biology approach is the integration of many different molecular interactions into extensive computer models that closely reflect real-life behavior of their underlying biological systems. Mathematical implementations of various biological systems have been proposed recently, including cell cycle control in yeast [1] and *Caulobacter crescentus *[2], and circadian rhythms of *Arabidopsis thaliana *[3] to name but just a few. Such studies are primarily designed to match known measurable phenotypes of the respective systems and reveal insights into the precise quantitative evolution of biochemical species over time. With a reasonable *in silico *implementation of a biological system at hand, predictions of knockout and perturbation effects can be performed by the computer.

For most biological systems, however, only qualitative information about regulatory interactions is available, which is not sufficient to implement precise kinetic rate laws for each biochemical reaction. A well-established workaround for this lack of information is the application of discretized modeling approaches. In Boolean methodology, for example, we abstract from actual molecule quantities and assign each player in the system the state *on *or *off *(e.g. active or inactive). Despite the simplicity of Boolean models we still assume them to provide information about the general dynamics and capabilities of the underlying system. Recently proposed Boolean models include developmental processes in *D. melanogaster *[4], the regulation of the mammalian cell cycle [5], the activation of T-cells [6] and EGFR signaling in human hepatocytes [7].

In [8] we described a novel technique called *HillCube *for the automatic transformation of Boolean models into systems of autonomous first-order ordinary differential equations (ODEs). HillCubes are based on multivariate polynomial interpolation and incorporate Hill kinetics (see Implementation), which are known to provide a good generalized approximation of the synergistic dynamics of gene regulation [9,10]. Important properties of the system like steady-state behavior are preserved during the transformation. Our methodology allows to enrich Boolean models built up from coarse information by features of quantitative models, such as intermediate expression levels, continuous transitions and different time-scales. Other approaches for the analysis of purely phenomenological regulatory networks have been developed recently (cf. e.g. [6,11]) but do not employ continuous, quantitative modeling.

Here we present a user-friendly implementation of the HillCube technique suitable for large-scale networks in a MATLAB/Octave toolbox called *Odefy*. This software provides convenient access to different model sources, the conversion process itself and various analysis and export methods (Figure (Figure1).1). Boolean models may be entered as sets of Boolean equations directly or created with the yEd graph editor [12]. The user may build conventional *interaction graphs *with activating and inhibiting edges or use an intuitive hypergraph representation of Boolean models [13]. In addition, models can be imported from the *CellNetAnalyzer *toolbox [6], GINsim [14] and the PBN toolbox [15]. After generating the ODEs, the user can easily adjust model parameters and perform time-course simulations using Odefy's graphical user interface. The ODE systems can be exported to MATLAB script files for further usage in MATLAB programs, to ODE script files for the R computing platform, to the SBML format, or to the well-established MATLAB Systems Biology Toolbox [16]. Due to the nice mathematical properties of the produced ODEs and the integration with state-of-the-art modeling tools, a variety of analysis methods can be immediately applied to the models generated by Odefy, including bifurcation analysis, parameter estimation, parameter sensitivity analysis and so on. For compatibility, we also integrated export options to the discrete model formats of the *Genetic Network Analyzer *[17] and *SQUAD *[18].

In this manuscript we first discuss the mathematical backgrounds and implementation details of the Odefy toolbox, including the different model import sources, analysis methods and export options. In the results section, two examples of quantitative modeling with our toolbox are given, namely a motif from stem cell differentiation and the regulatory network responsible for the establishment and stable maintenance of the mid-hindbrain boundary. We show the ease-of-use of the Odefy toolbox and demonstrate similar dynamical properties between a molecular model of the stem cell motif and the corresponding derived Odefy model. The mid-hindbrain example specifically emphasizes the importance of a fully automated conversion method from discrete to continuous models.

This section outlines the mathematical formulation of our automatic Boolean model conversion technique. For a detailed description of this methodology along with motivations, comparisons to similar approaches and application to a T-cell signaling model, we refer the reader to [8]. A Boolean model consists of *N *species *X*_{1}, *X*_{2}, ..., *X*_{N }each taking a value *x*_{i } {0, 1}. The value of *X*_{i }at time *t *+ 1 depends on the species *X*_{i1}, *X*_{i2}, ..., and is given by the Boolean update function *B*_{i }(*x*_{i1}, *x*_{i2}, ..., ) {0, 1}. In a discrete simulation, time is discretized and the values of *x*_{1}, *x*_{2}, ..., *x*_{N }at time *t *+ 1 are determined by synchronously setting:

The main idea is to convert the above discrete model into a continuous ODE model, where species *X*_{i }is allowed to take values [0, 1], and its temporal development is described by the ordinary differential equation (ODE):

The right hand side of this equation consists of two parts, an activation function describing the production of species *X*_{i }and a first-order decay term. An additional parameter *τ*_{i }is introduced to the system, which can be understood as the life-time of species *X*_{i}. can be considered a continuous homologue of the Boolean update function. The key point is how it can be obtained from *B*_{i }in a computationally efficient manner.

In Odefy, three different methods to transform *B*_{i }into are implemented. They are shortly described in the following. For simplicity of notation, we omit the subscript *i*.

The basis of all three transformation methods are the so-called *BooleCubes*:

which we obtain by multilinearly interpolating the Boolean function *B*, see Figure Figure2A2A.

The functions are affine multilinear. Many molecular interactions, however, are known to show a switch-like behavior, which can be modeled using sigmoid shaped *Hill functions *, see Figure Figure2B.2B. The two parameters *n *and *k *have a clear biological meaning. The Hill coefficient *n *determines the slope of the curve and is a measure of the cooperativity of the interaction. The parameter *k *corresponds to the threshold in the Boolean model, above which one defines the state of a species as *on*. Mathematically speaking, it is the value at which the activation is half maximal, i.e. equal to 0.5. We now introduce a Hill function *f*_{i }with parameters (*n*_{i}, *k*_{i}) for every interaction and define a new continuous function:

which we call *HillCubes*, see Figure Figure2C.2C. We can show that for sufficiently large Hill exponents *n*, there will be a steady-state of the continuous system in the neighborhood of each Boolean steady-state [8].

Note that Hill functions never assume the value 1, but approach it asymptotically. Hence, the HillCubes are not perfect homologues of the Boolean update function *B*. If this is desired a simple solution is to normalize the Hill functions to the unit interval. This yields another continuous (perfect) homologue of the Boolean function *B*:

which we call *normalized HillCube*, see Figure Figure2D2D.

The core functionality of Odefy is accessible through a set of functions for the MATLAB/Octave command line or via a Java-based graphical user interface. Figure Figure11 shows an overview of the complete Odefy tool-box. The following section provides detailed descriptions of the model definition and import process, ODE generation, model simulation and exporting.

An Odefy input model consists of a set of Boolean update rules for the underlying regulatory system. Our toolbox currently supports several possibilities to define such models:

(i) The user may enter a set of symbolic Boolean equations in text-form, allowing for the quick and intuitive generation of model structures (Figure (Figure3A).3A). Boolean equations consist of model variables and the three Boolean operators AND, OR and NOT. For the Odefy import process, we represent these operators by the MATLAB language-specific operators &&, || and ~, respectively. Throughout this manuscript, we stick to the common mathematical notation of ˄ for AND, for OR and ¬ for NOT.

(ii) Models can be derived from directed graphs created in the free yEd graph editing software [12]. The user builds an interaction graph of activating and inhibiting edges, which is then converted to an Odefy Boolean model (Figure (Figure3B).3B). Note that we need to specify how multiple regulatory inputs of a single factor are combined into a Boolean update rule. For this a generic logic of the form *f*(*X*) = (*A*_{1 } *A*_{2 } ... *A*_{m}) ¬(*I*_{1 } *I*_{2 } ... *I*_{n}) defined by three Boolean operators , , {˄, } is used, where *A*_{1}, ..., *A*_{m }is the set of activators and *I*_{1}, ..., *I*_{n }represent all the inhibitors of *X*. The Odefy default setting is to activate the output if at least one activator and no inhibitors are active. In order to create this behavior we choose = , = ˄, = resulting in

The assignment of Boolean operators can be changed during the import process into the Odefy toolbox. In addition to the possibility of inputing interaction graphs, we implemented an intuitive hypergraph-based representation of Boolean models in the *sum of product *form, which is capable of describing any Boolean update function [13] (Figure (Figure3C3C).

(iii) Odefy can be tightly integrated with the well-established CellNetAnalyzer (CNA) toolbox [6]. By a plugin-like menu interface the user can execute Odefy from within CNA and convert existing CNA models into systems of differential equations. Furthermore, parameter settings made in the CNA user interface are directly passed to Odefy and used for simulation and exporting.

(iv) Finally, Boolean models can be directly imported from the GINsim XML format [14] and the Probabilistic Boolean Networks toolbox [15].

The Odefy toolbox can efficiently handle large-scale models containing 50 players and more. One of the largest cellular Boolean model, a T-cell model with 94 nodes and a total of 123 regulatory interactions [19], can be transformed and simulated in less than one second on a standard workstation. Internally, Boolean models are stored as multidimensional arrays (i.e. hypercubes with edge length 2) for rapid element access and Boolean function evaluation. The time complexity of model generation lies in (2^{N}) with *N *being the highest degree of all nodes, yielding an exponential growth of computational runtime. The limiting size of Odefy models is thus not the number of nodes contained, but rather the highest number of incoming edges for any node in the model. For most regulated genes, however, we assume the number of modeled input regulatory factors to be equal to or less than 10, which can be handled on the order of one second per node by Odefy.

To account for systems consisting of multiple cells or, more generally, compartments driven by identical regulatory networks, Odefy contains an automated multicompartment expansion procedure. Given a Boolean model and the assignment of an *intercompartment *flag for a given set of factors in the model, Odefy generates a larger model corresponding to a linear row of connected compartments. Factors flagged as intercompartmental exhibit their influence towards the two neighboring cells and are combined using an OR logic (see also: Mid-hindbrain example below).

After model creation, the resulting ODE systems can be simulated directly by numerical integration algorithms or, alternatively, exported to various external model formats. Note that Boolean models as such are parameter-free, and the dynamical parameters for the ODEs have to be set externally. For convenience, Odefy employs a set of reasonable default values for all parameters in order to allow for a quick analysis of the system. Import, parameter adjustment, simulation and exporting can be accessed by the Odefy command line functionality as well as a graphical user interface (Figure (Figure4).4). These Java Swing-based user dialogs provide a platform-independent look and feel. They use the MATLAB-internal Java engine and therefore do not require an external Java runtime environment. For advanced MATLAB users and users of the Octave environment, we provide functions for convenient parameter access, Boolean state analysis (steady states and state-transition graph) and phase plane visualization of dynamic simulations.

Export formats for Odefy models include the MATLAB/Octave ODE script files, the Systems Biology (SB) Toolbox [16], the SBML format, script files for the R computing platform, the *Genetic Network Analyzer *(GNA) [17] and *SQUAD *[18]. SB Toolbox contains various advanced analysis functions for dynamical systems like parameter sensitivity and bifurcation analysis. The SBML format can be read by various systems biology software tools like COPASI [20] and CellDesigner [21] and thus provides a versatile interchange format. The GNA allows a structural analysis and qualitative simulations of systems of piece-wise linear ODEs. SQUAD analyzes discrete and continuous models using the standardized qualitative dynamical systems approach.

The Odefy toolbox is platform-independent due to the availability of MATLAB and Octave for all major operating systems and the direct integration of the Java Runtime Environment into MATLAB. It was verified to run smoothly on Windows, Linux KDE and GNOME as well as recent versions of Mac OS X. A detailed HTML documentation is included in the download package, which also provides a quick start guide to start working with the toolbox. Odefy is free for non-commerical and academic use. The toolbox including source codes can be downloaded at http://cmb.helmholtz-muenchen.de/odefy.

In the following we demonstrate the use of Odefy for the analysis of a simple regulatory motif. The mutual inhibitory switch (Figure (Figure5A)5A) is a well-known circuit involved in developmental processes and stem cell differentiation, e.g. in the hematopoietic system [22]. Despite its simplicity the circuit displays remarkable dynamic characteristics leading to the fate decision between opposing differentiation lineages. Various theoretical studies have been published recently investigating different aspects and molecular assumptions for this motif [23-25]. We discuss two different ways of formulating the interactions in this network in terms of Boolean equations. Multiple regulatory inputs (in this case self-activation and cross-inhibition) can either be combined using an AND or an OR logic for both factors. Figure Figure5B5B shows the MATLAB code that analyzes the OR logics version of our mutual switch network. After creating the model structure we calculate and output the steady states of the Boolean model (Figure (Figure5C5C).

We demonstrate the actual conversion into an ODE model and subsequent simulation within the Odefy toolbox. A two-dimensional phase plane projection of various initial values is drawn that displays the attractor landscape generated by the dynamical system (Figure (Figure5D,5D, the phase plane visualization for the corresponding AND logics model is shown in Figure Figure5E).5E). Note that this analysis reveals continuous decision boundaries between different attractors not apparent in the discrete model alone. Furthermore, two unstable steady states emerge which mark the switching points from one attractor basin to the other. In stem cell research, the central state is considered to be a pre-differentiation *priming *state whereas the other two states correspond to the regulatory program leading to the commitment to a certain cell lineage [26]. With our continuous mathematical representations we gain insights into the putative switching dynamics of this important differentiation switch in stem cells. After fitting simulated trajectories to observed time series of expression data, we could now determine rate parameters and understand the detailed time dynamics of the system.

We now employ the mutual inhibitory switch model discussed in the preceding section to address an important question for our novel modeling approach, namely whether the quantitative dynamics added to the discrete model are reasonable, or whether spurious, artifical effects are created by the method. In the study by Roeder et al. [23], a mechanistical model of the switch motif was proposed, which is based on actual biochemical reactions like promotor binding, transcription/translation and protein-protein interactions. The system was reduced to a two factor ODE by applying quasi-steady assumptions for the DNA and RNA species in the system. A comparison between simulation trajectories of the Odefy-converted model of the AND-gated switch and the Roeder model is displayed in Figure Figure6.6. Both systems have two non-zero stable steady states at similar positions, and the attractor basins for both states are virtually identical. Furthermore, both systems comprise a third, trivial steady state where both factors are zero. Interestingly, the parameter assignments we made for the simulation of the Odefy model, in order to achieve similarity between the model simulations, is qualitatively comparable with the parameter settings from the Roeder model. More precisely, the Roeder model employs a high *unspecific transcription rate *(we refer to the original publication for more details on the parameters), which effectively reduces the mutual inhibitory influences in relation to the autoregulatory activation of both factors. Accordingly, in our model we set the self-activation threshold to 0.01, which renders both factors strongly sensitive to their own expression levels. Taken together, we can reproduce important dynamical features of the reaction-based system by Roeder et al., inluding multistability, steady state positions, and the general shape of the attractor basins.

Our second example of dynamic modeling using Odefy concentrates on a multicellular biological system. During vertebrate development, the differentiation between mid- and hindbrain is determined by several transcription factors (e.g. Otx2, Gbx2) and secreted factors (e.g. Fgf8, Wnt1). These genes are stably expressed in a well-defined spatial pattern around the boundary between prospective mid- and hindbrain, the so-called mid-hindbrain boundary (Figure (Figure7A).7A). In a recent publication, we have applied both Boolean modeling and the HillCube conversion approach to this system [27]. In the following we will show how to use Odefy for automated model selection, that is the evaluation of an ensemble of regulatory networks with respect to stability of the known expression patterns. Figure Figure7B7B displays the MATLAB code required to fulfill this goal. First, we load a set of 9 regulatory networks known to give rise to the expression pattern along with 1000 random Boolean equations (not shown in the code) as a representative set of arbitrarily chosen regulatory networks. Then we iterate over all equation systems, generate a 6 cell multicompartment version of this model where the species representing the signaling molecules Fgf8 and Wnt1 are flagged as intercompartmental (Figure (Figure7C).7C). The multi-cell system is converted to an ODE system and simulated starting from the known stable expression state with default parameters *n *= 3, *k *= 0.5, *τ *= 1. If the activity of all players, in terms of exceeding the Hill threshold parameter, is still correct after a given amount of time steps, we consider the model to be valid. The results of this experiment show that indeed only 9 networks can give rise to the desired system behavior (Figure (Figure7D).7D). Analyzing these networks we see, in particular, that the maintenance of the boundary requires a mutual inhibition of Otx2 and Gbx2 and that these two transcription factors have antagonistic effects on Fgf8 and Wnt1 expression. Moreover, we find that Fgf8 and Wnt1 require each other for their stable maintenance. This agrees well with results from various loss-and gain-of-function experiments [28]. Note that while the small network in the former sections could still be handled manually, the model selection problem for the mid-hindbrain network demonstrates the absolute necessity for fully automated approaches as implemented in our toolbox. The system contains 6·4 = 24 differential equations with a total of up to 20 kinetic parameters for each compartment. Obviously, a model system of this size with parameter interdependencies due to multicompartmentality cannot reasonably be handled by manual mathematical modeling.

Precise mechanistic details about regulatory interactions required for the quantitative modeling of biological systems are rare. However, more qualitative, phenomenological information like *activation *and *inhibition *is frequently available. With Odefy we created a simple yet useful toolbox to bridge the gap between qualitative and quantitative modeling of regulatory networks. A variety of such discrete models is already available and can immediately be converted into ODE systems by our tool.

Quantitative modeling might reveal features not present in the original Boolean models. For instance, quantitative models allow for the estimation of system robustness with respect to parameter perturbations, even with ad-hoc parameter values. This provides insights into the general capability of the system to withstand external or intracellular fluctuations and has been demonstrated for various biological systems like Drosophila segmentation patterns [29] and the mid-hindbrain specification mentioned in this report. Furthermore, in [8] we determined parameter values by least-square fitting to experimental data in a T-cell signaling model. We could, amongst others, successfully predict relations between binding affinity constants of ligand-receptor interactions, which represent biochemical quantities not capturable in a Boolean framework.

In this report we explained the concepts of automatic conversion from Boolean models to systems of ordinary differential equations. Two example cases were discussed stressing (a) the ease-of-use of the Odefy toolbox as well as (b) the requirement for automated conversion methods for more realistic biological systems like the Mid-hindbrain boundary network. We demonstrated that a discrete model converted to an ODE by Odefy displays similar dynamical properties as a mechanistically derived ODE model of the same system. Here we could show that, even though the identity of dynamical parameters between both modeling approaches is substantially different, qualitatively similar parameter changes show similar results.

The integration of Odefy with other modeling applications through the import and export of models extends the scope of our toolbox. In particular, the SBML export functionality connects our toolbox to a broad variety of systems biology softwares supporting this common interchange format. With its novel modeling technique and its easy usability, Odefy will be a valuable tool for researchers aiming to understand the dynamics of gene regulation.

• **Project name**: Odefy

• **Project home page**: http://cmb.helmholtz-muenchen.de/odefy

• **Operating system(s)**: Platform independent

• **Programming language**: MATLAB/Octave

• **Other requirements**: MATLAB 7.1 or higher (no additional toolboxes required), Octave for non-GUI mode

• **License**: Free for non-commercial purposes

DMW and FJT developed the automatic conversion technique and carried out the Mid-hindbrain study. JK and SP developed the Odefy application. JK and DMW wrote the initial manuscript. All authors revised and approved the final manuscript.

The authors would like to thank Steffen Klamt for valuable feedback during the development of the toolbox and the integration of Odefy into the CellNetAnalyzer package, and Florian Bloechl for stimulating discussions about the manuscript and the methodology. The authors thank the anonymous reviewers for valuable comments and suggestions. This research was partially supported by the Initiative and Networking Fund of the Helmholtz Association within the Helmholtz Alliance on Systems Biology (project CoReNe).

- Chen K, Calzone L, Csikasz-Nagy A, Cross F, Novak B, Tyson J. Integrative analysis of cell cycle control in budding yeast. Molecular Biology of the Cell. 2004;11(8):3841–3862. doi: 10.1091/mbc.E03-11-0794. [PMC free article] [PubMed] [Cross Ref]
- Li S, Brazhnik P, Sobral B, Tyson J. A quantitative study of the division cycle of Caulobacter crescentus stalked cells. PLoS Comput Biol. 2008;11:e9. doi: 10.1371/journal.pcbi.0040009. [PMC free article] [PubMed] [Cross Ref]
- Locke JCW, Southern MM, Kozma-Bognar L, Hibberd V, Brown PE, Turner MS, Millar AJ. Extension of a genetic network model by iterative experimentation and mathematical analysis. Mol Syst Biol. 2005;11 doi: 10.1038/msb4100018. 2005.0013. [PMC free article] [PubMed] [Cross Ref]
- Albert R, Othmer HG. The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. J Theor Biol. 2003;11:1–18. doi: 10.1016/S0022-5193(03)00035-3. [PubMed] [Cross Ref]
- Fauré A, Naldi A, Chaouiya C, Thieffry D. Dynamical analysis of a generic Boolean model for the control of the mammalian cell cycle. Bioinformatics. 2006;11(14):e124–e131. doi: 10.1093/bioinformatics/btl210. [PubMed] [Cross Ref]
- Klamt S, Saez-Rodriguez J, Gilles ED. Structural and functional analysis of cellular networks with CellNetAnalyzer. BMC Syst Biol. 2007;11:2. doi: 10.1186/1752-0509-1-2. [PMC free article] [PubMed] [Cross Ref]
- Samaga R, Saez-Rodriguez J, Alexopoulos LG, Sorger PK, Klamt S. The logic of EGFR/ErbB signaling: theoretical properties and analysis of high-throughput data. PLoS Comput Biol. 2009;11(8):e1000438. doi: 10.1371/journal.pcbi.1000438. [PMC free article] [PubMed] [Cross Ref]
- Wittmann DM, Krumsiek J, Saez-Rodriguez J, Lauffenburger DA, Klamt S, Theis FJ. Transforming Boolean models to continuous models: methodology and application to T-cell receptor signaling. BMC Syst Biol. 2009;11:98. doi: 10.1186/1752-0509-3-98. [PMC free article] [PubMed] [Cross Ref]
- Hill A. The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves. J Physiol. 1910;11:4–7.
- Alon U. An Introduction to Systems Biology: Design Principles of Biological Circuits (Chapman & Hall/Crc Mathematical and Computational Biology Series) Chapman & Hall/CRC; 2006.
- Ruths D, Nakhleh L, Ram PT. Rapidly exploring structural and dynamic properties of signaling networks using PathwayOracle. BMC Syst Biol. 2008;11:76. doi: 10.1186/1752-0509-2-76. [PMC free article] [PubMed] [Cross Ref]
- yED graph editor. http://www.yworks.com/products/yed/
- Klamt S, Haus UU, Theis F. Hypergraphs and cellular networks. PLoS Comput Biol. 2009;11(5):e1000385. doi: 10.1371/journal.pcbi.1000385. [PMC free article] [PubMed] [Cross Ref]
- Gonzalez AG, Naldi A, Sánchez L, Thieffry D, Chaouiya C. GINsim: a software suite for the qualitative modelling, simulation and analysis of regulatory networks. Biosystems. 2006;11(2):91–100. doi: 10.1016/j.biosystems.2005.10.003. [PubMed] [Cross Ref]
- Probabilistic Boolean Networks toolbox. http://personal.systemsbiology.net/ilya/PBN/PBN.htm
- Schmidt H, Jirstrand M. Systems Biology Toolbox for MATLAB: a computational platform for research in systems biology. Bioinformatics. 2006;11(4):514–515. doi: 10.1093/bioinformatics/bti799. [PubMed] [Cross Ref]
- de Jong H, Geiselmann J, Hernandez C, Page M. Genetic Network Analyzer: qualitative simulation of genetic regulatory networks. Bioinformatics. 2003;11(3):336–344. doi: 10.1093/bioinformatics/btf851. [PubMed] [Cross Ref]
- Di Cara A, Garg A, De Micheli G, Xenarios I, Mendoza L. Dynamic simulation of regulatory networks using SQUAD. BMC Bioinformatics. 2007;11:462. doi: 10.1186/1471-2105-8-462. [PMC free article] [PubMed] [Cross Ref]
- Saez-Rodriguez J, Simeoni L, Lindquist JA, Hemenway R, Bommhardt U, Arndt B, Haus UU, Weismantel R, Gilles ED, Klamt S, Schraven B. A Logical Model Provides Insights into T Cell Receptor Signaling. PLoS Comput Biol. 2007;11(8):e163. doi: 10.1371/journal.pcbi.0030163. [PubMed] [Cross Ref]
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U. COPASI-a COmplex PAthway SImulator. Bioinformatics. 2006;11(24):3067–3074. doi: 10.1093/bioinformatics/btl485. [PubMed] [Cross Ref]
- Funahashi A, Matsuoka Y, Jouraku A, Morohashi M, Kikuchi N, Kitano H. CellDesigner 3.5: A Versatile Modeling Tool for Biochemical Networks. Proceedings of the IEEE. 2008;11(8):1254–1265. doi: 10.1109/JPROC.2008.925458. [Cross Ref]
- Cantor AB, Orkin SH. Hematopoietic development: a balancing act. Curr Opin Genet Dev. 2001;11(5):513–519. doi: 10.1016/S0959-437X(00)00226-4. [PubMed] [Cross Ref]
- Roeder I, Glauche I. Towards an understanding of lineage specification in hematopoietic stem cells: a mathematical model for the interaction of transcription factors GATA-1 and PU.1. J Theor Biol. 2006;11(4):852–865. [PubMed]
- Huang S, Guo YP, May G, Enver T. Bifurcation dynamics in lineage-commitment in bipotent progenitor cells. Dev Biol. 2007;11(2):695–713. doi: 10.1016/j.ydbio.2007.02.036. [PubMed] [Cross Ref]
- Chickarmane V, Enver T, Peterson C. Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility. PLoS Comput Biol. 2009;11:e1000268. doi: 10.1371/journal.pcbi.1000268. [PMC free article] [PubMed] [Cross Ref]
- Akashi K, He X, Chen J, Iwasaki H, Niu C, Steenhard B, Zhang J, Haug J, Li L. Transcriptional accessibility for genes of multiple tissues and hematopoietic lineages is hierarchically controlled during early hematopoiesis. Blood. 2003;11(2):383–389. doi: 10.1182/blood-2002-06-1780. [PubMed] [Cross Ref]
- Wittmann DM, Bölchl F, Trümbach D, Wurst W, Prakash N, Theis FJ. Spatial Analysis of Expression Patterns Predicts Genetic Interactions at the Mid-Hindbrain Boundary. PLoS Comput Biol. 2009;11(11):e1000569. doi: 10.1371/journal.pcbi.1000569. [PMC free article] [PubMed] [Cross Ref]
- Prakash N, Wurst W. Specification of midbrain territory. Cell Tissue Res. 2004;11:5–14. doi: 10.1007/s00441-004-0955-x. [PubMed] [Cross Ref]
- von Dassow G, Meir E, Munro EM, Odell GM. The segment polarity network is a robust developmental module. Nature. 2000;11(6792):188–192. doi: 10.1038/35018085. [PubMed] [Cross Ref]

Articles from BMC Bioinformatics are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |