# Modeling biochemical transformation processes and information processing with Narrator

- Johannes J Mandel
^{1, 2}Email author, - Hendrik Fuß
^{2}, - Niall M Palfreyman
^{1}and - Werner Dubitzky
^{2}

**8**:103

**DOI: **10.1186/1471-2105-8-103

© Mandel et al; licensee BioMed Central Ltd. 2007

**Received: **20 January 2007

**Accepted: **27 March 2007

**Published: **27 March 2007

## Abstract

### Background

Software tools that model and simulate the dynamics of biological processes and systems are becoming increasingly important. Some of these tools offer sophisticated graphical user interfaces (GUIs), which greatly enhance their acceptance by users. Such GUIs are based on symbolic or graphical notations used to describe, interact and communicate the developed models. Typically, these graphical notations are geared towards conventional biochemical pathway diagrams. They permit the user to represent the transport and transformation of chemical species and to define inhibitory and stimulatory dependencies. A critical weakness of existing tools is their lack of supporting an integrative representation of transport, transformation as well as biological information processing.

### Results

Narrator is a software tool facilitating the development and simulation of biological systems as *Co-dependence models*. The *Co-dependence Methodology* complements the representation of species transport and transformation together with an explicit mechanism to express biological information processing. Thus, Co-dependence models explicitly capture, for instance, signal processing structures and the influence of exogenous factors or events affecting certain parts of a biological system or process. This combined set of features provides the system biologist with a powerful tool to describe and explore the dynamics of life phenomena. Narrator's GUI is based on an expressive graphical notation which forms an integral part of the Co-dependence Methodology. Behind the user-friendly GUI, Narrator hides a flexible feature which makes it relatively easy to map models defined via the graphical notation to mathematical formalisms and languages such as ordinary differential equations, the Systems Biology Markup Language or Gillespie's direct method. This powerful feature facilitates reuse, interoperability and conceptual model development.

### Conclusion

Narrator is a flexible and intuitive systems biology tool. It is specifically intended for users aiming to construct and simulate dynamic models of biology without recourse to extensive mathematical detail. Its design facilitates mappings to different formal languages and frameworks. The combined set of features makes Narrator unique among tools of its kind. Narrator is implemented as Java software program and available as open-source from http://www.narrator-tool.org.

## Background

### Motivation

One goal of investigating the dynamic behavior of biochemical systems and processes is to understand biological causality in terms of regulation and control mechanisms. Computerized mathematical models of such dynamic behavior have become an important methodology in the pursuit of this goal. Many modeling and simulation tools have been reported in the literature. The Systems Biology Markup Language (SBML) Web site [1] lists more than 100 such tools, all of which adhere to the SBML standard, facilitating the exchange of biochemical reaction network models. In addition to systems modeling and simulation tools, the SBML site posts a large number of tools addressing other tasks, including tools to convert SBML files into other languages or formats (such as R [2], Matlab [3] or Mathematica [4]), to modify and combine different SBML files (see for example SBMLmerge [5]), and to interactively analyze and explore a large number of already developed biological models (e.g., Web-based kinetic modeling [6, 7] and kinetic model databases or modelbases [8]).

Several excellent tools have been developed for modeling and simulating biochemical networks, some of which provide sophisticated, GUI-based design tools for visual model development. Most of these tools are freely available as software packages or as online tools. Examples include Copasi [34], Virtual Cell [9], JDesigner [10], CellWare [12] and CellDesigner [11, 33], which are well developed and widely used in the systems biology community.

A critical weakness of the symbolic notations of these tools is their lack of supporting an integrative representation of biological transformation and transportation processes as well as biological information processing.

Within the context of this work, we define *species transformation* as the alteration or change of species participating in biochemical processes. Examples of species transformation include stoichiometric reactions, enzymatic reactions, synthesis and degradation processes or modification processes such as methylation and phosphorylation. *Biological transport* refers to processes that transport organic, biochemical or inorganic substances between or within cells. Protein translocation across subcellular compartments or material flow through biological membranes are examples of biological transport. *Biological information processing* controls the processes involved in species transformation and species transport. By processing the information provided by one or more information sources and by regulating biological transport and transformation processes accordingly, biological systems can control both the maintenance of a state and also state transitions from one state to the next. Based on this definition information sources can be of different biological and non-biological nature such as a set of transcription factors controlling the expression of a gene, or physical parameters influencing a biological process (e.g. temperature, osmotic pressure or pH-values). Typical examples of biological information processing include gene regulatory processes and cell-cell or cell-environment interactions mediated through receptor proteins. To simultaneously model species transformation and biological information processing is important because only this integrative approach can describe biochemical processes including the influencing factors that control such processes and their interplay. For example, a model of osmotic shock in yeast should capture species transformation and information processing because it is not possible to formulate the causal dependencies of osmotic shrinking and swelling of yeast cells by defining transformation processes exclusively. Instead the processing of internal and external ion concentration levels is fundamental in the response processes initiated by the yeast cells and thus must be integrated into the model to provide contextual causality of the underlying system. While the main focus of most existing tools is placed on the description of species transformation, the modeling of information processing which controls such transformations is an important feature missing in those tools. The main novelty of the Narrator tool is that it facilitates the integrative modeling and simulation of biological transformation, transport and information processing.

### Co-dependence models

*Co-dependence Methodology*[14], which consists of a set of concepts and tools facilitating the modeling and simulation of biological systems and processes. Among others, the development of this methodology was guided by the following needs: The methodology should

- 1.
Enable life scientists to express the full complexity of biological functioning, including species transport, transformation and information processing which guides these processes.

- 2.
Provide a visual or graphical notation suitable for use in informal pen-and-paper discussions and implementation as sophisticated graphical user interface of computerized modeling tools.

- 3.
Be executable on a computer to enable the testing of proposed dynamics

*in silico*.

### Co-dependence Notation

Most notations currently used in systems biology tools facilitate the description of conventional biochemical pathways and are based on mathematical *bipartite graphs* or *directed hypergraphs*. A bipartite graph (or bigraph) is defined by a set of graph vertices (nodes) decomposed into two disjoint sets such that no two graph vertices within the same set are adjacent, i.e., they are not joined by a graph edge (link or arc). A directed hypergraph is a graph in which generalized edges (also known as hyperarcs) may connect more than two nodes. The two types of graphs have been shown to be equivalent when used for representing transport and transformation of chemical species [13, 40]. In bipartite graphs, there are two kinds of node, and no edge can relate nodes of the same kind. This is useful for linking biological species to processes or reactants to products. Hyperarcs of directed hypergraphs can be used to relate *sets* of reactants to *sets* of products, where each set belong to the same kind of node. Additionally, both notations can represent inhibitory and stimulatory dependencies between chemical species and their transformation processes. Petri nets [18], the E-Cell process variable model or graphical notations used in JDesigner or CellWare belong to this class of graph.

A novel aspect in the Co-dependence Notation is the explicit labeling of *information sources* using small, unfilled circles. Interpreted dynamically, information sources are considered as mathematical variables that are integrated into the function or rate law of their associated process. This is illustrated in Figure 3c by the small Co-dependence model and its corresponding equations. Species A and species B, for instance, have an influence on the process transforming species A, B and C.

A thick arc in combination with a small circle indicates that there exists material as well as informational flow between a species and a process node. In this *co-dependent* relationship between a species and a process, the species influences the process that is transforming the species.

Co-dependence models are fully isomorph to systems of ordinary differential equations (ODEs). The rate of change of each species represents an ODE and, as shown in Figure 3c and Figure 3d, it is the visual emphasis of information sources in the Co-dependence Notation (see information source labels) that enables the presentation of the individual factors playing a functional role in the ODEs. Process nodes that are representing the transformation of species are then describing the coupling of the ODEs.

An additional novel aspect of the Co-dependence Notation is the optional use of *information processing nodes* to represent constants or computed logical entities such as rate constants, external stimuli or genetic switches. Information processing nodes also belong to the class of process nodes, yet, as shown in Figure 3d, they are rendered as large circles. Information processing nodes also serve as information sources and can be linked to other process nodes accordingly. Thus, links in Co-dependence models cannot only connect species nodes with process nodes but can also directly relate two or more process nodes with each other. This facilitates the representation of both species transformation and information processing within a single conceptual framework.

From a graph-theoretical point of view the visual representation of information processing is made possible by defining a link or edge type which is able to connect process and species nodes as well as relating two process nodes. Such an informational relation between two process nodes has no direct biochemical interpretation and therefore is not considered in notations derived from biochemical pathway diagrams. Yet, as we shall demonstrate, this extension is extremely useful when incorporating abstractions or additional information into models as is for example required when integrating biochemical and biophysical aspects in one model. Consequently, in Co-dependence models, process nodes describe both the transformation of concentration levels of their participating species as well as the processing and supply of information. In particular, the ability to link processes *informationally* is of considerable value in scenarios illustrated by the following four scenarios.

#### Scenario 1: Incorporating exogenous variables

*SLN1*is directly affected by the

*Turgor pressure*which again depends on the

*Internal*and

*External osmotic pressure*. Thus, we can visually describe which species is contributing as an osmosensor to the regulation of the HOG pathway and how this species depends on various biophysical influences.

#### Scenario 2: Describing detailed steps of information processing

Figure 4 illustrates what we understand by the term *detailed steps of information processing*. Sometimes it is convenient to elaborate a single informational dependency into a sequence of processing steps. For example, in Figure 4, where the simple dependency of process *v1* on the species *glycerol* has been expanded into the more detailed description involving *Experimental input signal*, *External* and *Internal osmotic pressure*, and *Turgor pressure*. These detailed steps of information processing allow us to visually close the loop between the biological system *internal osmotic pressure* and its environments (*external osmotic pressure*).

#### Scenario 3: Reducing diagram complexity with auxiliary processes

In ODE-based modeling, *differential-algebraic equations* are often used to reduce the complexity of the differential equations by simplifying their right-hand expressions with auxiliary variables. Biochemical models frequently make use of this method. However, although SBML's AssignmentRule tag provides support for it, the process nodes in Narrator can improve the readability of equation systems by structuring information processing in a similar way. Differential-algebraic equations thus translate to Co-dependence models without any increase in diagram complexity.

*src_activity*is used to calculate the total enzymatic activity of that kinase. In the model depicted in Figure 5, Src is represented in four different phosphorylation states. Due to allosteric effects each state contributes to kinase activity to a different degree (represented by different weights of informational links A-D). The value of the

*src_activity*node is then reused in processes that represent reactions catalyzed by this kinase (links labeled E, F, G in Figure 5).

*Xenopus laevis*embryos [17] (shown in Figure 6), we also make use of an auxiliary process node to structure the composition and usage of the

*Pool*of all removed forms of

*Cyclin:Cdk2*dimers contained in the model. This variable

*Pool*is defined as auxiliary variable in the differential-algebraic equation system [17] and could be seamlessly translated into the Co-dependence Notation.

Both examples (Figure 5 and Figure 6) illustrate how the proliferation of links can be avoided by introducing auxiliary process nodes.

#### Scenario 4: Emphasizing important process parameters

Quantitative models of biological processes strongly depend on constant parameters such as Hill coefficients, half-saturation constants, rate constants or specific growth and degradation rates. In a Co-dependence model we can describe selected parameters explicitly with process nodes to emphasize their importance to the overall model or to allow direct manipulation of the parameter value via our software tool Narrator.

In all abovementioned examples it is up to the modeler to decide which level of detail is appropriate for a model. For example, if we incorporated only biochemical transformation processes into our model it would look very similar to other bipartite approaches such as Hybrid Object Net [18] or E-Cell [19].

On the other hand, information processing structures resemble the models of the *System Dynamics* methodology. In this method, developed by Forrester [32], dynamic system models are built with stock and flow diagrams which implement the *principle of accumulation* [20]. Here, dynamic behavior is described by flows that accumulate in stocks and a third node type called *converter* can be used to describe graphically the information processing in a system. This graphical notation however is not suitable for modeling stoichiometric relations in metabolic reactions and thus is not applicable to biochemical modeling [21].

As shown in Figure 3e, the Co-dependence Methodology provides a mechanism for decomposing complex systems into simpler elements called *compartments*. Compartments may themselves contain other nested elements [14]. This hierarchical structure allows to model systems of any appreciable complexity and is also implemented in Narrator. This mechanism can be used to encapsulate parts of a model that correspond to one cell compartment but it can be also useful to provide different levels of detail of a model or to integrate a set of models as submodels into one model.

## Implementation

### Architecture and Workflow

Narrator is implemented in Java 5 on the basis of the model-view-controller design pattern, which separates core (application) model functionality from the presentation and control logic that uses this functionality [45].

Figure 7 depicts a typical Narrator workflow. Via the user interface, the modeler designs and manipulates dynamic models according to the Co-dependence Methodology [14]. This involves defining the model structure and assigning initial concentration values and functions to all species and process nodes. Narrator provides a click-drag-drop environment for placing and controlling nodes and links. Narrator contains its own interpreter system which compiles functions from a mathematical expression in string form into an execution tree. The instructions of such trees are executed during simulation to return the runtime function values. A set of classes for different mathematical operators has been created using JavaCC [22] technology. To facilitate simulation, the software tool implements a set of numerical integration methods. Currently these methods include Euler's method, and 2^{nd} and 4^{th} order Runge-Kutta and Runge-Kutta-Fehlberg with adaptive step size (for a detailed description of the integration methods and their implementation see for example [46]). Additionally, Narrator implements Gillespie's direct method to simulate stochastic processes [23]. To examine the resulting dynamics of a model the user can create time-series diagrams, phase portraits [47] and tabular output.

### Narrator models

*species*,

*process*,

*compartment*and

*link*. Compartment objects encapsulate a finite number of Species and Process objects and can be used to model cellular compartments or to decompose a model into constituent elements. Link objects connect Process objects with Species or other Process objects. They model the propagation of change between species and processes (e.g. biochemical transformation) or the transport of information to the processes (e.g. stimulation, inhibition or information processing).

Each model created using Narrator contains one root compartment, which in turn can contain additional compartments, processes and species. The relationship among the classes Link, Process and Species further enables the mapping to and execution of numerical simulations of Narrator models directly within the model structure. Process objects coordinate via Link objects the augmentation and depletion of the concentration levels represented by Species objects. As shown in Figure 8, Link objects can associate Node and Process objects for this purpose.

Similarly, the class structure can be used to instantiate and simulate stochastic processes. Currently, Narrator implements Gillespie's direct method [23] where the runtime function values in the process nodes determine the probability of each reaction to occur in a certain time interval. This time interval is also calculated probabilistically, and when a reaction occurs, all participating Species objects are updated according to their stoichiometry, which can be defined through Link objects.

### Reporters

Narrator provides two reporting components to present the resulting dynamics data of a Narrator model. These are portraits and tables respectively. Portraits are based on the open-source chart library jfreechart [24] and can either be used to generate time-series plots or phase plots. Tables present the time-dependent development of selected variables and can export their information as ASCII-formatted files.

### Animation

### SBML Export and Import

To facilitate a standard-compliant exchange of Co-dependence models, Narrator provides an export and import mechanism for storing Co-dependence models in SBML level 2 version 1 [48]. However, Narrator does only embrace elements of the SBML language that match to elements of the Co-dependence Methodology. Essentially, these are species, reactions, parameters, rules and the assignment of initial values or functions for these components. Narrator does currently not support unit and function definitions and provides no automatic mechanisms to layout imported SBML models. Consequently, Narrator cannot read all SBML models created by other SBML-platforms, but SBML models created with Narrator are SBML compatible.

For transforming mathematical String expressions into MathML we make use of the open-source library jeks [25] and for reading MathML elements we us the class jigcell.sbml2.math.MathMLConvertor which is part of the open-source distribution JigCell [35].

## Results

To evaluate Narrator, we described and simulated a model of the protein kinase C (PKC) pathway taken from Pettinen et al. [26] which used this pathway model as a test case for comparing and evaluating different simulation tools for biochemical networks. In their study Pettinen et al. emphasize the importance of usability and the benefits of exchanging standards, and conclude that simulation tools should support the integration of exogenous variables for the modeling of stimuli since biological systems are not separated from their environment. Another important feature discussed in this article is the automatic estimation of parameters where no information is available.

*CAStimulus*linked to the reaction

*R7*to model the exogenous Ca

^{2+}stimulus for the PKC pathway following a sine wave with time period of 50 seconds.

We used the same model to evaluate the SBML compatibility of Narrator. The PKC Narrator model was exported into SBML and the generated XML file was validated with the SBML Validator available on the SBML Web site [30]. SBML confirmed that the exported file adhered with the SBML format. We also imported this XML file successfully into the SBML platform Copasi after setting the stimulus of CA^{2+} to a constant value.

Feature and requirement summary of Narrator. Summary of Narrator features and requirements according to criteria used in [42], addressing systems requirements, cross-compatibility, diagrammatic user interfaces or analytical capabilities of the simulators.

Type/Name of Criterion | Value of Criterion |
---|---|

Website | |

Hardware requirements | ~10 MB disk space |

Input (model definition) | Diagrammatic |

Software requirements | JRE version 1.5 |

Operating System | OS independent |

Import | SBML |

Export | SBML |

Open source | Yes (Sourceforge) |

Systems Biology Workbench | No |

Graphical representation of multiple compartments | Yes |

Representation of interface reactions | Yes |

Representation of multiple reactant types | No |

Representation of regulatory interactions | Yes |

Automated network layout | No |

Library of predefined rate expressions | No |

Editable list of all parameters/initial conditions | No |

Implementation of non-integer stoichiometries | Yes |

Ordinary differential equations | Yes |

Partial differential equations | No |

Stochastic | Yes |

Automated set-up of compartmentalization | No |

Direct calculations of steady states | No |

Steady-state stability analysis | No |

Steady-state sensitivity analysis | No |

Time-dependent sensitivity analysis | No |

Automated parameter scanning | No |

Time-dependent perturbations | Yes |

Parameter estimation/optimization | No |

Narrator provides a set of example models which are directly accessible via the menu bar. Some models are taken from Klipp et al. [15] and demonstrate Narrator's ability to model and exchange small systems and to simulate their dynamic behavior. The other models are example models developed for this study and in earlier work introducing the Co-dependence Methodology [14].

## Discussion

In general, graphical notations of systems biology tools are designed to capture critical elements that determine the *structure* and govern the *behavior* of a biological system or process. Some of the existing notations and tools emphasize the structural and others the behavioral dimension. Typically, structure-oriented notations use a plethora of graphical symbols and notational conventions to characterize qualitative details of species types and their biochemical transformations. Dynamics-oriented notations, on the other hand, capture information relevant mainly to the dynamic aspects of a system and usually subsume different species and transformation types within a small collection of graphical symbols (see for example bipartite graphs and hypergraphs in Figure 3ab).

Prominent structure-oriented notations include molecular interaction maps (MIMs) [36], process diagrams [33], and the Edinburgh [37] notation. To ensure efficient and correct use of these notations, many have been implemented in graph-based modeling tools (e.g., CellDesigner or Edinburgh Pathway Editor) and are accompanied by simulation environments (e.g. CellDesigner). MIMs however lack of software tools to assist the use of this graphical notation.

The graphical notation of the Narrator software tool emphasizes the development and exploration of *dynamic* models. This is explained by Narrator's main focus on the dynamic structure and properties of a biological system. With regard to SBML, structure-oriented notations correspond to SBML including Systems Biology Ontology terms for the SBML model elements, while dynamics-oriented notations can be understood as the notational counterpart of pure SBML, focusing on the dynamic description of biochemical models. In SBML level 2 version 2, Systems Biology Ontology (SBO) terms [41] can be optionally integrated into SBML model elements. The motivation of the SBO effort is to provide controlled vocabularies tailored to systems biology. In addition to facilitating consistency and sharing, SBO will also support visual rendering of the biochemical reactions.

A recently initiated effort called Systems Biology Graphical Notation aims to develop and standardize a structure-oriented graphical notation for representing biochemical networks [39]. Once available, this standard is likely to make user interfaces, model development, sharing and communication, and system interoperation more effective and efficient [39].

There are many methods listed in Table 1, such as methods for the automatic estimation of parameters values or steady-state analysis, that Narrator currently does not support but which do other modeling tools support. Copasi for instance integrates diverse numerical methods used in systems biology comprising deterministic, stochastic and hybrid numerical integration methods. Theses methods can, for example, divide networks into deterministic and stochastic subnetworks which run significantly faster than pure stochastic networks, or can deal with the numerical integration of deterministic networks in the presence of stiffness [34]. Similarly, Cellware provides an extensive algorithm library also comprising deterministic, stochastic and hypbrid integration methods [51]. The very special feature of Cellware however is that it is grid enabled supporting the distributed simulation of larger biochemical networks [12]. For 2D and 3D simulations Virtual Cell provides integrations methods for solving ODEs and partial differential equations. These methods can then be used for calculating the reaction and diffusion equation on arbitrary surfaces. Copasi, CellWare and Virtual Cell also provide parameter estimation and parameter scan components. Some tools such as JDesigner and CellDesigner have an interface to the Systems Biology Workbench (SBW) [50] which is an open source framework for connecting different software applications for modeling, analyzing or visualizing biochemical network models. Via the SBW interface, JDesigner for example uses the script based biochemical network simulation tool Jarnac [49] as a simulation server.

While ongoing developments implement more features such as bifurcation analysis, the objective of Narrator is not to cover the full spectrum of systems biology methods. Instead, the emphasis of Narrator and one of its major advantages over similar tools lies in the expressiveness and intuitiveness of its graphical notation and the sophisticated and highly usable implementation of a GUI supporting the use of this notation in the context of dynamics modeling and simulation. The focus of Copasi, for example, lies in the many numerical integration, parameter estimation and parameter scan methods. The graphical user interface of Copasi, however, is based on simple dialog boxes providing entry masks to define the chemical equation and rate expression for each reaction. The dynamics-based graphical design tools of JDesigner, Virtual Cell and CellWare, on the other hand, are based on directed hypergraphs or bipartite graphs. As shown in Figure 3 and discussed above (see section Co-dependence Notation), a possibility which is not explicitly allowed in such approaches is the ability to represent non-biological information sources such as temperature or osmotic pressure and their biological information processing accordingly. As also demonstrated above, the unique feature of Narrator, to link process nodes informationally, is a valuable modeling tool in many modeling scenarios and is not supported by the graphical notations of the abovementioned tools since edges in bipartite and hypergraph approaches can only relate biological species to processes or reactants to products. The fundamental structure of the structure-based graphical notation of CellDesigner is similar to directed hypergraphs and thus has the same restrictions in this respect.

## Conclusion

Narrator uses the simple, intuitive graphical notation of Co-dependence models [14] to express the dynamics of a biological system. A Co-dependence model describes dynamical relationships between the components of a system. Narrator simulates the behavior of such a system as an unfolding 'narrative' played out among the components.

Narrator is a stand-alone Java application for the graph-based modeling of dynamic network structures of biochemical systems. It extends the visual description of species transformation with the notion of information processing and therefore introduces a new concept into the multitude of graphical modeling tools used in the field of systems biology. Here, the provided intuitive design tool, which implements the Co-dependence Notation, supports the model building in a way that is accessible without an advanced mathematical background but which is still more oriented to mathematical dynamic systems than to biochemical pathway diagrams. This makes Narrator also applicable to any discipline requiring dynamic modeling.

The primary purpose of Narrator is to facilitate the creation and manipulation of computational models. Therefore, an ongoing effort in the development of Narrator is to maintain and improve the tool's SBML compatibility allowing a further exploration of Narrator models with simulation and analysis methods implemented in other tools.

## Availability and requirements

Narrator is available at [31] under the Lesser GNU Public license. The tool is entirely written in Java and runs with Java 5. It is platform-independent, which we have approved on Windows, Unix and Mac OS X machines. Additional information on installation and usage is also provided at [31]. The source code for Narrator is available at sourceforge http://sourceforge.net/projects/narratortool/.

## Declarations

### Acknowledgements

The authors would like to thank the 2004 to 2006 third- and fourth-year bioinformatics students of Weihenstephan University of Applied Sciences for their help in developing the software tool Narrator.

## Authors’ Affiliations

## References

- Systems Biology Markup Language[http://www.sbml.org]
- SBMLR[http://epbi-radivot.cwru.edu/SBMLR/]
- Schmidt H, Jirstrand M: Systems biology toolbox for MATLAB: A computational platform for research in systems biology. Bioinformatics 2006, 22: 514–515. 10.1093/bioinformatics/bti799View ArticlePubMedGoogle Scholar
- MathSBML[http://sbml.org/software/mathsbml]
- SBMLmerge Project[http://sysbio.molgen.mpg.de/sbmlmerge/index.html]
- Olivier BG, Snoep JL: Web-based kinetic modelling using JWS online. Bioinformatics 2004, 20: 2143–2144. 10.1093/bioinformatics/bth200View ArticlePubMedGoogle Scholar
- Lee D, Yun C, Cho A, Hou BK, Park S, Lee SY: WebCell: A web-based environment for kinetic modeling and dynamic simulation of cellular networks. Bioinformatics 2006, 22: 1150–1151. 10.1093/bioinformatics/btl091View ArticlePubMedGoogle Scholar
- Le Novere N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, Li L, Sauro H, Schilstra M, Shapiro B, Snoep JL, Hucka M: BioModels database: A free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucl Acids Res 2006, 34: D689–691. 10.1093/nar/gkj092PubMed CentralView ArticlePubMedGoogle Scholar
- Schaff J, Loew LM: The virtual cell. Pacific Symposium on Biocomputing 1999, 228–239.Google Scholar
- Sauro HM: An introduction to biochemical modeling using JDesigner. Claremont, CA: Keck Graduate Institute; 2004.Google Scholar
- Funahashi A, Tanimura N, Morohashi M, Kitano H: CellDesigner: A process diagram editor for gene-regulatory and biochemical networks. BIOSILICO 2003, 159–162. 10.1016/S1478-5382(03)02370-9Google Scholar
- Dhar PK, Meng TC, Somani S, Ye L, Sakharkar K, Krishnan A, Ridwan ABM, Wah SHK, Chitre M, Hao Z: Grid cellware: The first grid-enabled tool for modelling and simulating cellular processes. Bioinformatics 2005, 21: 1284–1287. 10.1093/bioinformatics/bti143View ArticlePubMedGoogle Scholar
- Deville Y, Gilbert D, Helden JV, Wodak SJ: An overview of data models for the analysis of biochemical pathways. Briefings in Bioinformatics 2003, 4: 246–259. 10.1093/bib/4.3.246View ArticlePubMedGoogle Scholar
- Mandel JJ, Palfreyman NM, Dubitzky W: Modelling Codependence in Biological Systems. IET Systems Biology 2007, 1(1):18–32. 10.1049/iet-syb:20060002View ArticlePubMedGoogle Scholar
- Klipp E, Herwig R, Kowald A, Wierling C, Lehrach H: Systems Biology in Practice: Concepts, Implementation and Application. Wiley-VCH, Weinheim 2005. ISBN 3–527–31078–9 ISBN 3-527-31078-9Google Scholar
- Fuß H, Dubitzky W, Downes CS, Kurth MJ: Bistable Switching and Excitable Behaviour in the Activation of Src at Mitosis. Bioinformatics 2006, 22(14):e158-e164. 10.1093/bioinformatics/btl201View ArticlePubMedGoogle Scholar
- Ciliberto A, Petrus MJ, Tyson JJ, Sible JC: A kinetic model of the cyclin e/cdk2 developmental timer in xenopus laevis embryos. Biophys Chem 2003, 104: 573–589. 10.1016/S0301-4622(03)00060-7View ArticlePubMedGoogle Scholar
- Drath R: Hybrid object nets: An object oriented concept for modeling complex hybrid systems. Hybrid Dynamical Systems. 3rd International Confernce on Automation of Mixed Processes, ADPM'98:Reims 1998.Google Scholar
- Takahashi K, Ishikawa N, Sadamoto Y, Sasamoto H, Ohta S, Shiozawa A, Miyoshi F, Naito Y, Nakayama Y, Tomita M: E-cell 2: Multi-platform E-cell simulation system. Bioinformatics 2003, 19: 1727–1729. 10.1093/bioinformatics/btg221View ArticlePubMedGoogle Scholar
- Sterman J: Business Dynamics. McGraw-Hill, Inc; 2000.Google Scholar
- Mandel JJ, Palfreyman NM, Lopez JA, Dubitzky W: Representing bioinformatics causality. Brief Bioinform 2004, 5: 270–283. 10.1093/bib/5.3.270View ArticlePubMedGoogle Scholar
- Java Compiler Compiler[https://javacc.dev.java.net/]
- Gillespie D: A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J Comp Phys 1976, 22: 403–434. 10.1016/0021-9991(76)90041-3View ArticleGoogle Scholar
- JFreeChart[http://www.jfree.org/jfreechart/]
- Jeks[http://www.eteks.com/jeks/]
- Pettinen A, Aho T, Smolander O, Manninen T, Saarinen A, Taattola K, Yli-Harja O, Linne M: Simulation tools for biochemical networks: Evaluation of performance and usability. Bioinformatics 2005, 21: 357–363. 10.1093/bioinformatics/bti018View ArticlePubMedGoogle Scholar
- Manninen T, Saarinen A, and Linne ML: Simulation study of differential equation model for protein kinase C signaling. In Report 3. Institute of Signal Processing and Institute of Mathematics, Tampere, Finland; 2004.Google Scholar
- Octave[http://www.gnu.org/software/octave/]
- Radhakrishnan K, Hindmarsh AC: Description and Use of LSODE, the Livermore Solver for Ordinary Differential Equations. Lawrence Livermore Laboratory, Livermore 1993. CAUCRL-ID-113855 CAUCRL-ID-113855Google Scholar
- SBML Validator[http://sbml.org/validator/validate.php]
- Narrator Project[http://www.narrator-tool.org]
- Forrester J: Industrial Dynamics. MIT Press Cambridge; 1961.Google Scholar
- Kitano H, Funahasi A, Matsuoka Y, Oda : Using process diagrams for the graphical representation of biological networks. Nature Biotechnology 2005, 23: 961–966. 10.1038/nbt1111View ArticlePubMedGoogle Scholar
- 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, 22(24):3067–3074. 10.1093/bioinformatics/btl485View ArticlePubMedGoogle Scholar
- Vass M, Allen N, Shaffer CA, Ramakrishnan N, Watson LT, Tyson JJ: The jigcell model builder and run manager. Bioinformatics 2004, 20(18):3680–3681. 10.1093/bioinformatics/bth422View ArticlePubMedGoogle Scholar
- Kohn KW, Aladjem MI, Weinstein JN, Pommier Y: Molecular Interaction Maps of Bioregulatory Networks: A General Rubric for Systems Biology. Mol Biol Cell 2006, 17: 1–13. 10.1091/mbc.E05-09-0824PubMed CentralView ArticlePubMedGoogle Scholar
- Moodie SL, Sorokin A, Groyanin I, Ghazal P: A Graphical Notation to describe the Logical Interactions of Biological Pathways. Journal of Integrative Bioinformatics 2006, 3(2):36. [http://journal.imbio.de/index.php?paper_id=36]Google Scholar
- CellML language[http://www.cellml.org/tutorial/notation/]
- Systems Biology Graphical Notation[http://www.sbgn.org/]
- Gauges R, Rost U, Sahle S, Wegner K: A model diagram layout extension for SBML. Bioinformatics 2006, 22(15):1879–1885. 10.1093/bioinformatics/btl195View ArticlePubMedGoogle Scholar
- Systems Biology Ontology[http://www.ebi.ac.uk/sbo/]
- Alves R, Antunes F, Salvador A: Tools for kinetic modeling of biochemical networks. Nature Biotechnology 2006, 24: 667–672. 10.1038/nbt0606-667View ArticlePubMedGoogle Scholar
- D'haeseleer P: Closing the circle of osmoregulation. Nature Biotechnology 2005, 23: 941–942. 10.1038/nbt0805-941View ArticlePubMedGoogle Scholar
- Gennemark P, Nordlander B, Hohmann S, Wedelin D: A simple mathematical model of adaptation to high osmolarity in yeast. In Silico Biology 2006, 6: 0018.Google Scholar
- Sun Microsystems, Java BluePrints: Model-View-Controller.2002. [http://java.sun.com/blueprints/patterns/MVC-detailed.html]Google Scholar
- Press WH, Teukolsky SA, Vetterling WT, Flannery BP: Numerical Recipes in C: The Art of Scientific Computing. Cambridge University Press New York, USA; 1992.Google Scholar
- Tabor M: Chaos and integrability in nonlinear dynamics: an introduction. Wiley New York; 1989.Google Scholar
- Hucka M, Finney A, Bornstein BJ, Keating SM, Shapiro BE, Matthews J, Kovitz BL, Schilstra MJ, Funahashi A, Doyle JC, Kitano H: Evolving a Lingua Franca and Associated Software Infrastructure for Computational Systems Biology: The Systems Biology Markup Language (SBML) Project. Systems Biology 2004, 1: 41–53. 10.1049/sb:20045008View ArticlePubMedGoogle Scholar
- Sauro HM: JARNAC: a system for interactive metabolic analysis. In Animating the Cellular Map. Edited by: Hofmeyr J-HS, Rohwer JM, Snoep JL. Stellenbosch, University Press Stellenbosch; 2001:221–228.Google Scholar
- Sauro HM, Hucka M, Finney A, Wellock C, Bolouri H, Doyle J, Kitano H: Next generation simulation tools: the Systems BiologyWorkbench and BioSPICE integration. OMICS 2003, 7(4):355–72. 10.1089/153623103322637670View ArticlePubMedGoogle Scholar
- Dhar P, Meng TC, Somani S, Ye L, Sairam A, Chitre M, Hao Z, Sakharkar K: Cellware – a multi-algorithmic software for computational systems biology. Bioinformatics 2004, 20(8):1319–1321. 10.1093/bioinformatics/bth067View ArticlePubMedGoogle Scholar

## 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.