ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra
 Franziska Hinkelmann^{1, 2},
 Madison Brandon†^{3},
 Bonny Guang†^{4},
 Rustin McNeill†^{5},
 Grigoriy Blekherman^{1},
 Alan VelizCuba^{1, 2} and
 Reinhard Laubenbacher^{1, 2}Email author
DOI: 10.1186/1471210512295
© Hinkelmann et al; licensee BioMed Central Ltd. 2011
Received: 14 February 2011
Accepted: 20 July 2011
Published: 20 July 2011
Abstract
Background
Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agentbased models, to gain a better understanding of them. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. There exist software tools to analyze discrete models, but they either lack the algorithmic functionality to analyze complex models deterministically or they are inaccessible to many users as they require understanding the underlying algorithm and implementation, do not have a graphical user interface, or are hard to install. Efficient analysis methods that are accessible to modelers and easy to use are needed.
Results
We propose a method for efficiently identifying attractors and introduce the webbased tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a longstudied problem in computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness and robustness. For a large set of published complex discrete models, ADAM identified the attractors in less than one second.
Conclusions
Discrete modeling techniques are a useful tool for analyzing complex biological systems and there is a need in the biological community for accessible efficient analysis tools. ADAM provides analysis methods based on mathematical algorithms as a webbased tool for several different input formats, and it makes analysis of complex models accessible to a larger community, as it is platform independent as a webservice and does not require understanding of the underlying mathematics.
Background
Mathematical modeling is a crucial tool in understanding the dynamic behavior of complex biological systems. In addition to the popular ordinary differential equations (ODE) models, discrete models are now increasingly used for this purpose [1–3]. Model types include (probabilistic) Boolean networks, logical networks, Petri nets, cellular automata, and agentbased (individualbased) models, to name the most commonly found ones [4–9]. While discrete models tend to be more intuitive than those based on differential equations, they do not have the broad range of mathematical analysis tools available that have been developed for ODE models. For small models, exhaustive enumeration of all possible state transitions of the model is the method of choice. But since the size of the state space grows exponentially in the number of model variables, this method is very limited in its applicability. For larger models sampling methods can be used to get some information about model dynamics. There are several existing sophisticated software tools available that allow users to analyze and simulate discrete networks, focused on a particular model type. These tools use a variety of computational and analytical tools for analysis purposes, with a range of different user interfaces; see, e.g., [10–16]. They will be discussed in detail in a later section.
The software tool introduced in this paper, Analysis of Dynamic Algebraic Models (ADAM) complements existing software packages in several ways. By translating models into the rich mathematical framework of polynomial dynamical systems over a finite number system, we can bring to bear a variety of theoretical results, computational algorithms, and available software packages from computer algebra and computational algebraic geometry on the analysis of model dynamics. For this purpose we provide implemented algorithms that import models created in with other packages, so that the user does not need to learn a new mathematical framework [17, 18]. The basic computational workhorse underlying our software tool is the (symbolic) solution of systems of (nonlinear) polynomial equations over a finite number system. This is a wellstudied problem in computer algebra and sophisticated algorithms are implemented for this purpose, which we make use of. An efficient computational implementation results in the ability to analyze model dynamics for quite large discrete models without having to resort to heuristic algorithms. We offer ADAM as a web service, avoiding the problems associated with software downloads and different computational platforms.
Results and Discussion
In this manuscript, we present the webbased tool ADAM, Analysis of Dynamic Algebraic Models [19], a tool to study the dynamics of a wide range of discrete models. ADAM provides efficient analysis methods based on mathematical algorithms as a webbased tool for several different input formats, and it makes analysis of complex models accessible to a larger community, as it is platform independent as a webservice and does not require understanding of the underlying mathematics. ADAM is the successor to DVD, Discrete Visualizer of Dynamics [20], a tool to visualize the temporal evolution of small polynomial dynamical systems.
As the underlying computational approach, we propose a novel method to identify attractors of a discrete model. This method relies on the fact that many discrete models can be translated into the algebraic framework of polynomial dynamical systems. Using these polynomials, one can construct a system of polynomial equations, such that its solutions correspond to fixed points or limit cycles. Thus, the problem of identifying attractors becomes equivalent to solving a system of polynomial equations over a finite field. This is a longstudied problem in computer algebra, and can usually be solved efficiently by using Gröbner basis methods [21]. We emphasize that this method is not a new mathematical algorithm to solve polynomial equations, but a novel approach to the analysis of discrete dynamical systems that uses a novel encoding of the periodic points of such a system as the solutions of polynomial systems derived from the model when expressed in the algebraic framework. ADAM allows users unfamiliar with polynomial dynamical systems or Gröbner bases to benefit from this efficient algorithm. We tested the method on several examples and had an average run time of less than one second, comparable to the performance of other software tools; and we were able to identify limit cycles of systems with more than 32 variables in less than one second.
In addition to providing access to mathematical theory for efficient analysis, algebraic models are a unifying framework and systematic approach for several model types. This allows for an effective comparison of heterogeneous models, such as a Boolean network model and an agentbased model. For community integration in the biological sciences, ADAM contains a model repository of previously published models available in ADAM specific format [22]. This allows new users to familiarize themselves quickly with ADAM and to validate and experiment with existing models. In the following section, we discuss general features of ADAM briefly and explain new features in more detail.
General Features of ADAM
ADAM is a tool for analyzing different types of discrete models. It automatically converts discrete models into polynomial dynamical systems, that is, time and state discrete dynamical systems described by polynomials over a finite field (see Appendix A.1 for definition and example). The dynamics of the models is then analyzed by using various computational algebra techniques. Even for large systems, ADAM computes key dynamic features, such as steady states, in a matter of seconds. ADAM is available online and free of charge. It is platform independent and does not require the installation of software or a computer algebra system.
ADAM translates the following inputs into (probabilistic) polynomial dynamical systems and can then analyze them.
ADAM also translates Petri nets generated with Snoopy and we plan to implement analysis methods for Petri nets in future versions.
ADAM's main application is the analysis of the dynamic features of a model, which includes the identification of stable attractors. These are either steady states, i.e., timeinvariant states, or limit cycles, i.e., timeinvariant sets of states. ADAM is capable of identifying all steady states and limit cycles of length up to a userspecified length m. The process of finding long limit cycles is quite slow for large models. However, in biological models limit cycles are likely to be short, so that m can be chosen to be small in general, i.e., less than 10.
The temporal evolution of the model can be visualized by the state transition graph, the directed graph of all possible states and edges indicating their transitions, also called the state space. For small enough models, i.e., less than eleven variables, ADAM generates a graph of the complete state space; for larger models, ADAM uses algebraic algorithms to determine dynamic properties. Independent of network size, ADAM generates a wiring diagram. The wiring diagram, also known as dependency graph, shows the static relationship between the variables. All edges in ADAM's wiring diagrams are functional edges, that is, there exists at least one state such that a change in the input variable causes a change in the output variable (see Appendix A.2 for more details). This means that ADAM determines all nonfunctional edges, which is oftentimes of interest.
With ADAM, one can also study the temporal evolution of userspecified initial states. The trajectory of a state describes the state's evolution, and it can be computed by repeatedly applying the transition function until an attractor is reached.
All of these features can be computed assuming synchronous updates or sequential updates according to an updateschedule specified by the user. Note that the steady states are the same independent of the update schedule. This is due to the fact that updating any variable at a steady state does not change its value. It is irrelevant for a steady state analysis whether updates are considered to happen sequentially or simultaneously.
For probabilistic networks, i.e., models in which each variable has several choices of local update rules, ADAM can generate a graph of all possible updates. This means that states in the phase space can have outdegree greater than one, since different transitions are possible. ADAM can find all true steady states, in the context of probabilistic networks, meaning all states that are timeinvariant independent of the choice of update function. For further information of probabilistic networks, see [6].
For Boolean networks, ADAM calculates all functional circuits (see Appendix A.2). Positive functional circuits are a necessary condition for multistationarity. For a certain class of Boolean networks, namely conjunctive/disjunctive networks, ADAM computes a complete description of the phase space as described in [23]. For further details on conjunctive networks, see Appendix B.2.
In summary, ADAM can generate the following outputs.

wiring diagram

phase space for small models

steady states (for deterministic and probabilistic systems)

limit cycles of specified length m

trajectories originating from a given initial state until a stable attractor is found

dynamics for synchronous or sequential updates

functional circuits for Boolean networks

a complete description of the phase space for conjunctive/disjunctive networks.
Applications
We show how to use ADAM on a wellunderstood model of the expression pattern of the segment polarity genes in Drosophila melanogaster. Albert and Othmer developed a model for embryonic pattern formation in the fruit fly Drosophila melanogaster [24]. Their Boolean model consists of 60 variables, resulting in a state space with more than 10^{18} states. They analyze the model for steady states by manually solving a system of Boolean equations. They also analyze the temporal evolution of a specific initial state corresponding to the wild type expression pattern by repeatedly applying the Boolean update rules until a steady state is found. The update schedule of the model is synchronous with the exception of activation of SMO and the binding of PTC to HH (activation of PH), which are assumed to happen instantaneously. This can be accounted for by substituting the equations for SMO and PH into the update rules for other genes and proteins, rather than using SMO and PH themselves.
Correspondence of genes and variable names
Cell 1  SLP x1  wg x2  WG x3  en x4  EN x5  hh x6  HH x7  ptc x8  PTC x9  PH x10  SMO x11  ci x12  CI x13  CIA x14  CIR x15 

Cell 2  SLP x16  wg x17  WG x18  en x19  EN x20  hh x21  HH x22  ptc x23  PTC x24  PH x25  SMO x26  ci x27  CI x28  CIA x29  CIR x30 
Cell 3  SLP x31  wg x32  WG x33  en x34  EN x35  hh x36  HH x37  ptc x38  PTC x39  PH x40  SMO x41  ci x42  CI x43  CIA x44  CIR x45 
Cell 4  SLP x46  wg x47  WG x48  en x49  EN x50  hh x51  HH x52  ptc x53  PTC x54  PH x55  SMO x56  ci x57  CI x58  CIA x59  CIR x60 
Genes and proteins present in steady state
compartment 1  en, EN, hh, HH, SMO 

compartment 2  ptc, PTC, PH, SMO, ci, CI, CIA 
compartment 3  SLP, PTC, ci, CI, CIR 
compartment 4  SLP, wg, WG, ptc, PTC, PH, SMO, ci, CI, CIA 
Genes and proteins present in initial state
compartment 1  en, hh 

compartment 2  ptc, ci 
compartment 3  SLP, ptc, ci 
compartment 4  SLP, wg, ptc, ci 
Furthermore, we used ADAM to verify that there are no limit cycles of length two or three. The model has not been analyzed previously for limit cycles. The absence of two and threecycles strengthens confidence in the model, since oscillatory behavior has not been observed experimentally. Computations for limit cycles of length greater than three have not been conducted, as composing the system several times with itself is computationally complex. The model file in ADAM format can be accessed at [22].
Benchmark Calculations
In addition to the published models in [25], we analyzed randomly generated networks that have the same structure that we expect from biological systems, namely sparse, i.e., while the number of nodes in a biological network may be quite large, each node is affected only by a small number of other nodes, and robust, i.e., small number of attractors. We tested a total of 50 networks with 50150 nodes (10^{15}  10^{45} states) and an average of average indegrees of 1.6848. The steady state calculations took less than half a second for each network on a 2.7 GHz computer.
Comparison to Other Systems
Software Comparison
Steady State Analysis  Limit Cycle Analysis  Input Format  System Requirements  

ADAM  Yes^{‡}  Yes^{◊}  Boolean (or polynomial) functions Logical Models (GINsim)  None, web based 
GINsim  Yes^{‡}  For small models  Parameters (nonzero truth tables) Logical Model  Java virtual machine^{○} 
BoolNet R package  For small ^{†} models  For small ^{†} models  Boolean functions  R statistics software 
DDLab  For small models  For small models  Logical tables  ○ 
BN/PBN Matlab Toolbox  For small models  For small models  Logical tables  Matlab 
GINsim (Gene Interaction Network simulation) is a package designed for the analysis of gene regulatory networks [10]. As input, it accepts logical models. Logical models are an extension of Boolean models; they consist of similar switchlike rules, but allow for a finer data discretization with more than two states per variable, e.g., low, medium, and high. Logical models can be updated synchronously or asynchronously. For the latter, the temporal evolution of a logical model is nondeterministic because the variables are updated randomly in an asynchronous fashion. In either case, updates of every variable are continuous, meaning that no variable changes its value by more than one unit in one timestep, see section Remarks about Logical Models for a detailed discussion.
GINsim provides algorithms that use binary decision diagrams (BDD) for the determination of steady states [7]. Analysis of limit cycles is executed by simulating every trajectory, i.e., generating the complete state space, called state transition graph in GINsim, and therefore limited by network size. We tested GINsim on logical models with up to 72 variables; determining the steady states took less than one second. More complex logical networks were not available to us.
BoolNet R package provides methods for inference and analysis of synchronous, asynchronous, and probabilistic Boolean networks [11]. It is a package for the free statistics software R, and it is run via the R commandline. It is helpful, if the user is already familiar with R. Steady state analysis is implemented as exhaustive search of the state space, heuristic search, random walk, or Markov Chain analysis [6].
Nonheuristic analysis is limited to networks with 29 variables. For larger networks, steady states can be inferred heuristically, which does not guarantee that all steady states are identified.
DDLab is an interactive graphics software for discrete models, including cellular automata, Boolean and multivalued networks [13]. As it is mainly a visualization tool, analysis is based on exhaustive enumeration of the state space, and model size is limited to 31 variables.
BN/PBN Toolbox is a toolbox written in Matlab[14]. It uses the state transition matrix to compute attractors. Statistics for networks with more than 27 variables cannot be computed ("Maximum variable size allowed by the program is exceeded"). In addition to analyzing deterministic Boolean networks, the toolbox can analyze probabilistic Boolean networks and calculate statistics such as numbers and sizes of attractors, basins, transient lengths, Derrida curves, percolation on 2D lattices, and influence matrices.
Remarks about Logical Models
In this manuscript, we distinguish between three different update types: synchronous, sequential according to an update schedule, and asynchronous. ADAM allows for synchronous or sequential updates according to a given update schedule. In models with synchronous updates, all variables are updated simultaneously at every time step. In models with sequential updates according to an update schedule, all variables are updated at every timestep in the order given by the schedule. Both these model types are deterministic.
In models with asynchronous updates, as is common for logical models, one variable is updated at random at every time step, which results in a nondeterministic model. Models with sequential updates according to an update schedule produce dynamics that are different from that of models with asynchronous updates, i.e., logical models.
In GINsim, all models are continuous in the sense that at each timestep, each variable increases or decreases by at most one unit. Though logical models are discrete, there are no jumps skipping intermediate states. For example, in a model with three states, low, medium, and high, no variable can drop from high to low in a single update step. This interpretation is different from the common meaning of continuous, which usually refers to models of ordinary or partial differential equations. The parameters entered in GINsim specify the target value towards which the variable changes, i.e., the value increases by one, decreases by one, or remains constant if the target value is larger, smaller, or equal than the initial value, respectively. The state transition graph generated with ADAM might differ from the state transition graph generated in GINsim. To obtain the exact same state transition graph, every variable in the logical model must contain an explicit selfloop, and all parameters must be entered such that the target value differs by at most one from the value of the variable to be updated. Any logical model can be specified in this way without changing its state transition graph. Boolean models are always continuous.
Multivalued models
next state of x_{2}  low x_{2}  medium x_{2}  high x_{2} 

x_{1} absent  low x_{2}  medium x_{2}  high x_{2} 
x_{1} present  medium x_{2}  high x_{2}  high x_{2} 
extension x_{1} present  medium x_{2}  high x_{2}  high x_{2} 
Architecture
ADAM is available as a webbased tool that does not require any software installation. ADAM's user interface is implemented in HTML. We use JavaScript to generate a dynamic website that adapts as the user makes various choices. This simplifies the process of entering a model. For example, after defining the model type, i.e., Polynomial Dynamical System, Probabilistic Network, Petri net, and Logical Model the next line changes to the number of states, kbound, or nothing, appropriately. Input can be entered directly into the text area, or uploaded as a text document.
All mathematical algorithms are programmed in Macaulay2 [26]. Macaulay2 is a powerful computer algebra system. The routines for which fast execution is crucial are implemented in C/C++ as part of the Macaulay2 core. Logical Models and Petri nets in XML format are parsed using Ruby's XmlSimple library. The interplay between HTML and Macaulay2 is also programmed in Ruby.
Output graphs are generated with Graphviz's dot command. When Simulation is chosen as analysis method, Graphviz's ccomps  connected components filter for graphs is used to count the connected components. A Perl script directs the execution of the Graphviz commands.
Model Repository
A model repository is part of the ADAM website [22]. The repository consists of a collection of several previously published models in ADAM format. The models are extracted from publications, and rewritten in ADAM specific format, i.e., all variables are renamed to x_{ i } and the update rules from the original publication are reformulated as Boolean rules or polynomials. The central repository with models in a unified format allows for quick verification and experimentation with published models. By changing parameters or initial states, users can gain a better understanding of the models.
New users can also use the repository to quickly familiarize themselves with the main functionalities of ADAM. In addition to the model itself, the database entries contain a short summary of the biological system and relevant graphs, together with an analysis of dynamic features determined by ADAM and their biological explanation. The repository is work in progress by researchers from several institutions generating more entries for the repository. We invite all interested researchers to submit their models. Because of their intuitive nature, discrete models are an excellent introduction to mathematical modeling for students of the life sciences. ADAM's model repository is a great starting point to familiarize students with the abstraction of discrete models such as Boolean networks.
Conclusions
Discrete modeling techniques are a useful tool for analyzing complex biological systems and there is a need in the biological community for easy to use analysis software. ADAM provides efficient methods as a webbased tool and will allow a larger community to use complex modeling techniques, as it is platform independent and does not require the user to understand the underlying mathematics. Upon translating discrete models, such as logical networks, Petri nets, or agentbased models into algebraic models, rich mathematical theory becomes available for model analysis, e.g., for steady state and limit cycle analysis.
After extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for sparse systems with few attractors, a structure maintained by most biological systems. All algorithms have been included in the software package ADAM[19], which is userfriendly and available as a free webbased tool. ADAM is highly suitable to be used in a classroom as a first introduction to discrete models because students can use it without going through an installation process.
ADAM provides methods to analyze the key dynamic features, such as steady states and limit cycles, for largescale (probabilistic) Boolean networks and logical models. ADAM unifies different modeling types by providing analysis methods for all of them and thus can be used by a larger community.
We hope to expand ADAM to a more comprehensive Discrete Toolkit which incorporates new analysis methods, better visualization, and automatic conversion for more model types. We also hope to analyze controlled algebraic models and expand theory to stochastic systems.
Methods
Logical models, Petri nets, and Boolean networks are converted automatically into the corresponding polynomial dynamical system as described in [17], so that algorithms from computational algebra can be used to analyze the dynamics. In polynomial dynamical systems over a finite field, states of a variable are assigned to values in the field, and the update (or transition) rule for each variable is given as a polynomial rather than a Boolean or logical expression. For more details, see Appendix A.1. Using these polynomials, one can construct systems of polynomial equations, such that their solutions correspond to fixed points or limit cycles. Thus, the problem of identifying attractors becomes equivalent to solving a system of polynomial equations over a finite field. This is a longstudied problem in computer algebra, and can usually be solved efficiently by Gröbner basis methods.
Gröbner basis calculation is for polynomial systems what GaussJordan elimination is for linear systems: a structured way to transform the original system to triangular shape without changing its solution space. The triangular shape of the resulting systems allows for stepwise retrieval of the solutions of the system. For a more in depth discussion of Gröbner bases, see for example [21].
In the worst case, computing Gröbner bases for a set of polynomials has complexity doubly exponential in the number of solutions to the system. However, in practice, Gröbner bases are computable in a reasonable time. It has been suggested, that in robust gene regulatory networks genes are regulated by only a handful of regulators [27]. Thus, the polynomial dynamical systems representing such biological networks are sparse, i.e., each function depends only on a small subset of the model variables. From our experience, a Gröbner basis calculation for sparse systems with few attractors, a structure common for biological systems, is actually quite fast.
A Mathematical Background
A.1 Polynomial Dynamical Systems
with coordinate functions f_{ i } ∈ k[x_{1},..., x_{ n } ], the ring of polynomials in the variables x_{ i } , with coefficients in k. Iteration of f results in a timediscrete dynamical system. A PDS can be used to describe the dynamic behavior of a biological system: every variable x_{ i } corresponds to a biological substrate, for example a protein or gene, and the polynomials f_{ i } describe the evolution of x_{ i } depending on the previous state of the variables x_{1},..., x_{ n } .
PDS have several dynamic features of biological relevance. These include the number of components, component sizes, steady states, limit cycles, and limit cycle lengths.
together with a probability distribution for every coordinate that assigns the probability that a specific function is chosen to update that coordinate. The coordinate functions f_{ i, j } are elements in k[x_{1},..., x_{ n } ]. Probabilistic PDS, specifically Boolean probabilistic networks (PBN), have been studied extensively in [6]. ADAM analyzes probabilistic PDS. It can simulate the complete state transition graph for sufficently small models, by generating every possible transition and labeling the edge with its probability according to the distribution. If no distribution is given, ADAM assumes a uniform distribution on all functions. For large networks, ADAM's Algorithm choice computes steady states of probabilistic networks.
A.2 Functional Edges
An edge in the wiring diagram from x_{ i } to x_{ j } is considered functional, if there exists a state such that , where a and b are values for x_{ i } , in other words, if there is at least one state, such that changing only x_{ i } but keeping all other values fixed, changes the next state of x_{ j } . In ADAM, all edges in the wiring diagram are functional. For Boolean networks, ADAM identifies all functional elementary circuits. An elementary circuit is a finite closed path in the wiring diagram in which all the nodes are distinct. The existence of functional circuits is a necessary condition for multistationarity and limit cycles. For a further discussion of functional circuits, see [7]. For multivalued networks, circuit analysis has not yet been implemented.
B Algorithms
B.1 Analysis of stable attractors
Every attractor in a PDS is either a steady state or a limit cycle. For small models, ADAM determines the complete phase space by enumeration, for large models, ADAM computes steady states and limit cycles of a given length. A state is a steady state, if it transitions to itself after one update of the system. A state is part of a limit cycle of length m, if, after m updates, it results in itself. Any steady state of a PDS satisfies the equation f(x) = x, as no coordinate of x is changing as it is updated. Similarly, states of a limit cycle of length m satisfy the equation f ^{ m } (x) = x. ADAM computes all steady states by solving the system f_{ i } (x)  x_{ i } = 0 for i ∈ {1,..., n} simultaneously. To efficiently solve the resulting systems of polynomial equations, we first compute the Gröbner basis in lexicographic order for the ideal generated by the equations. Choosing a lexicographic order allows to easily obtain the solutions [21]. We use the Gröbner basis algorithms distributed with Macaulay2, version 1.3.1.1, and found that for quotient rings over a finite field the implementation 'Sugarless' is more efficient than the default algorithm with 'Sugar' [26, 29]. For limit cycles of length m, the solutions of f^{ m } (x) = x are found and then grouped into cycles, by applying f to each of the solutions.
results in the solutions (0, 0, 0), (0, 1, 0), (0, 1, 1), (1, 1, 1). (0, 0, 0) is a steady state, and (0, 1, 0), (0, 1, 1), (1, 1, 1) are elements of a limit cycle of length 3. For all m > 3, f ^{ m } (x) = x has no solutions, that means the system f has exactly two attractors, a steady state a a limit cycle of length 3.
B.2 Conjunctive/Disjunctive Networks
Some classes of networks have a certain structure that can be exploited to achieve faster calculations. Jarrah et al. show that for conjunctive (disjunctive) networks key dynamic features can be found with almost no computational effort [23]. Conjunctive (respectively disjunctive) networks consist of functions using only the AND (respectively OR) operator. ADAM comes with an implementation of this algorithm to analyze dynamics in the case of conjunctive (disjunctive) networks. Currently, this option is only implemented for networks with strongly connected dependency graphs.
Declarations
Acknowledgements
Dimitrova, Clemson University; J. Adeyeye, WinstonSalem State University; B. Stigler, Southern Methodist University; R. Isokpehi, Jackson State University are currently expanding ADAM's Model Repository. Funding for this work was provided through U.S. Army Research Office Grant Nr. W911NF0910538, National Science Foundation Grant Nr. CMMI0908201, and National Science Foundation Grant Nr. 0755322.
Authors’ Affiliations
References
 Sahin O, Frohlich H, Lobke C, Korf U, Burmester S, Majety M, Mattern J, Schupp I, Chaouiya C, Thieffry D, Poustka A, Wiemann S, Beissbarth T, Arlt D: Modeling ERBB receptorregulated G1/S transition to find novel targets for de novo trastuzumab resistance. BMC Systems Biology 2009, 3: 1. [http://www.biomedcentral.com/1752–0509/3/1] 10.1186/1752050931PubMed CentralView ArticlePubMed
 Heiser L, Wang N, Talcott C, Laderoute K, Knapp M, Guan Y, Hu Z, Ziyad S, Weber B, Laquerre S, Jackson J, Wooster R, Kuo L, Gray J, Spellman P: Integrated analysis of breast cancer cell lines reveals unique signaling pathways. Genome Biology 2009, 10: 1–17. 10.1186/gb2009101r1View Article
 Calzone L, Tournier L, Fourquet S, Thieffry D, Zhivotovsky B, Barillot E, Zinovyev A: Mathematical Modelling of CellFate Decision in Response to Death Receptor Engagement. PLoS Comput Biol 2010, 6(3):e1000702. 10.1371/journal.pcbi.1000702PubMed CentralView ArticlePubMed
 Steggles LJ, Banks R, Shaw O, Wipat A: Qualitatively modelling and analysing genetic regulatory networks: a Petri net approach. Bioinformatics 2007, 23: 336–343. 10.1093/bioinformatics/btl596View ArticlePubMed
 Sackmann A, Heiner M, Koch I: Application of Petri net based analysis techniques to signal transduction pathways. BMC Bioinformatics 2006, 7: 482. [http://www.biomedcentral.com/1471–2105/7/482] 10.1186/147121057482PubMed CentralView ArticlePubMed
 Shmulevich I, Dougherty ER, Kim S, Zhang W: Probabilistic Boolean networks: a rulebased uncertainty model for gene regulatory networks. Bioinformatics 2002, 18(2):261–274. 10.1093/bioinformatics/18.2.261View ArticlePubMed
 Naldi A, Thieffry D, Chaouiya C: Decision diagrams for the representation and analysis of logical models of genetic networks. In CMSB'07: Proceedings of the 2007 international conference on Computational methods in systems biology. Berlin, Heidelberg: SpringerVerlag; 2007:233–247.
 Bonchev D, Thomas S, Apte A, Kier LB: Cellular automata modelling of biomolecular networks dynamics. SAR and QSAR in Environmental Research 2010, 21: 77–102. [http://www.informaworld.com/10.1080/10629360903568580] 10.1080/10629360903568580View ArticlePubMed
 Laubenbacher R, Jarrah AS, Mortveit H, Ravi SS: Encyclopedia of Complexity and System Science. New York: Springer Verlag 2009 chap. A mathematical foundation for agentbased computer simulation;
 Naldi A, Berenguier D, Fauré A, Lopez F, Thieffry D, Chaouiya C: Logical modelling of regulatory networks with GINsim 2.3. Biosystems 2009, 97(2):134–139. 10.1016/j.biosystems.2009.04.008View ArticlePubMed
 Müssel C, Hopfensitz M, Kestler HA: BoolNetan R package for generation, reconstruction and analysis of Boolean networks. Bioinformatics 2010, 26(10):1378–1380. 10.1093/bioinformatics/btq124View ArticlePubMed
 Rohr C, Marwan W, Heiner M: Snoopya unifying Petri net framework to investigate biomolecular networks. Bioinformatics 2010, 26(7):974–975. 10.1093/bioinformatics/btq050View ArticlePubMed
 Wuensche A: DDLab  Discrete Dynamics Lab.2010. [http://www.ddlab.com/]
 Shmulevich I, Lähdesmäki H: BN/PBN MATLAB Toolbox.2005. [http://personal.systemsbiology.net/ilya/PBN/PBN.htm]
 Dubrova E, Teslenko M: A SATBased Algorithm for Finding Attractors in Synchronous Boolean Networks. IEEE/ACM Transactions on Computational Biology and Bioinformatics 2010, in press.
 Ruths D, Muller M, Tseng JTT, Nakhleh L, Ram PT: The signaling petri netbased simulator: a nonparametric strategy for characterizing the dynamics of cellspecific signaling networks. PLoS computational biology 2008, 4(2):e1000005+. 10.1371/journal.pcbi.1000005PubMed CentralView ArticlePubMed
 VelizCuba A, Jarrah AS, Laubenbacher R: Polynomial algebra of discrete models in systems biology. Bioinformatics 2010, 26(13):1637–1643. 10.1093/bioinformatics/btq240View ArticlePubMed
 Hinkelmann F, Murrugarra D, Jarrah AS, Laubenbacher R: A Mathematical Framework for Agent Based Models of Complex Biological Networks. Bulletin of Mathematical Biology 2010, 1–20.
 Hinkelmann F, Brandon M, Guang B, McNeill R, VelizCuba A, Laubenbacher R: ADAM, Analysis of Dynamic Algebraic Models.[http://adam.vbi.vt.edu]
 Laubenbacher R: DVD  Discrete Visualizer of Dynamics.[http://dvd.vbi.vt.edu/]
 Cox DA, Little J, O'Shea D: Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, 3/e (Undergraduate Texts in Mathematics). Secaucus, NJ, USA: SpringerVerlag New York, Inc; 2007.View Article
 Hinkelmann F: Model Repository of Discrete Biological Systems.[http://dvd.vbi.vt.edu/repository.html]
 Jarrah A, Laubenbacher R, VelizCuba A: The Dynamics of Conjunctive and Disjunctive Boolean Network Models. Bulletin of Mathematical Biology 2010, 72: 1425–1447. 10.1007/s115380109501zView ArticlePubMed
 Albert R, Othmer HG: The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster . Journal of Theoretical Biology 2003, 223: 1–18. 10.1016/S00225193(03)000353View ArticlePubMed
 Naldi A, Thieffry D, Chaouiya C: GINsim  model repository.[http://gin.univmrs.fr/GINsim/model_repository.html]
 Grayson DR, Stillman ME: Macaulay2, a software system for research in algebraic geometry.1992. [http://www.math.uiuc.edu/Macaulay2/]
 Leclerc RobertD: Survival of the sparsest: robust gene networks are parsimonious. Mol Syst Biol 2008., 4: 10.1038/msb.2008.52
 Jarrah AS, Laubenbacher R, Stigler B, Stillman ME: Reverseengineering of polynomial dynamical systems. Adv Appl Math 2007, 39: 477–489. 10.1016/j.aam.2006.08.004View Article
 Giovini A, Mora T, Niesi G, Robbiano L, Traverso C: "One sugar cube, please" or selection strategies in the Buchberger algorithm. In ISSAC '91: Proceedings of the 1991 international symposium on Symbolic and algebraic computation. New York, NY, USA: ACM; 1991:49–54.View Article
Copyright
This article is published under license to 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.