- Open Access
Rapidly exploring structural and dynamic properties of signaling networks using PathwayOracle
© Ruths et al; licensee BioMed Central Ltd. 2008
- Received: 18 December 2007
- Accepted: 19 August 2008
- Published: 19 August 2008
In systems biology the experimentalist is presented with a selection of software for analyzing dynamic properties of signaling networks. These tools either assume that the network is in steady-state or require highly parameterized models of the network of interest. For biologists interested in assessing how signal propagates through a network under specific conditions, the first class of methods does not provide sufficiently detailed results and the second class requires models which may not be easily and accurately constructed. A tool that is able to characterize the dynamics of a signaling network using an unparameterized model of the network would allow biologists to quickly obtain insights into a signaling network's behavior.
We introduce PathwayOracle, an integrated suite of software tools for computationally inferring and analyzing structural and dynamic properties of a signaling network. The feature which differentiates PathwayOracle from other tools is a method that can predict the response of a signaling network to various experimental conditions and stimuli using only the connectivity of the signaling network. Thus signaling models are relatively easy to build. The method allows for tracking signal flow in a network and comparison of signal flows under different experimental conditions. In addition, PathwayOracle includes tools for the enumeration and visualization of coherent and incoherent signaling paths between proteins, and for experimental analysis – loading and superimposing experimental data, such as microarray intensities, on the network model.
PathwayOracle provides an integrated environment in which both structural and dynamic analysis of a signaling network can be quickly conducted and visualized along side experimental results. By using the signaling network connectivity, analyses and predictions can be performed quickly using relatively easily constructed signaling network models. The application has been developed in Python and is designed to be easily extensible by groups interested in adding new or extending existing features. PathwayOracle is freely available for download and use.
- Signaling Network
- Signal Flow
- Time Block
- Signaling Node
- Transition Firing
Reconstructing cellular signaling networks and understanding how they work are major endeavors in cell biology. The scale and complexity of these networks, however, render their analysis using experimental biology approaches alone very challenging. As a result, computational methods have been developed and combined with experimental biology approaches, producing powerful tools for the analysis of these networks. These tools aid biologists in interpreting existing experimental findings, evaluating hypotheses, enumerating possible biological behaviors, and, ultimately, in quickly designing experiments that maximize the amount of useful information gained. By assisting biologists in maximizing the amount of information obtained from their experiments through improved experimental design and more thorough analysis of results, computational tools increase the pace of scientific discovery.
Biological network analysis can generally be classified as either structural or dynamic . Structural analysis provides insights into global properties of the network, among them decomposition of the network into functional modules (e.g., ), enumeration of signaling paths connecting arbitrary protein pairs (e.g., [3–5]), and the identification of key pathways that determine the behavior of the network (e.g., [2, 6–10]). Dynamic methods, on the other hand, simulate the actual propagation of signals through a network by predicting the changes in the concentration of signaling proteins over time. These predictions will be of varying degrees of resolution and accuracy, depending largely on the accuracy and level of detail of the model from which they are produced.
The prevailing methods for dynamic analysis involve systems of ordinary differential equations (ODEs) [11, 12]. These approaches require kinetic parameters for the individual biochemical reactions involved in the signaling process. This requirement often poses a significant hurdle for researchers as the numerical values of such parameters are difficult to obtain and may be the object of the researcher's project in the first place. In , we presented a novel signaling network simulation method which uses a non-parametric Petri net model of network to predict the signal flow under various experimental conditions. Our simulation method uses a novel technique to approximate the interaction speeds and predicts the qualitative behavior of the signaling network dynamics.
The advantage of our method over ODEs is the wide availability of connectivity-based models of signaling networks, and the relative speed with which they can be constructed. Numerous databases exist which catalog known signaling interactions (e.g., [14–16]). Thus, the existence and type (activating or inhibition) of an interaction can often be inferred directly from literature and/or these databases. This presents a stark contrast to the kinetic parameters required by ODEs, the numerical values for many of which must be determined experimentally for each experimental condition and cell line of interest .
In this paper, we present the software tool PathwayOracle, an integrated environment for connectivity-based structural and dynamic analysis of signaling networks, supporting
• visualization of signaling network connectivity;
• two versions of the simulation method described in  where
– the first allows prediction of signal flow through a given network for a specific experimental condition, and
– the second predicts the difference in signal flow through a given network induced by two different experimental conditions;
• enumeration of the paths connecting arbitrary pairs of nodes in the network; and
• visualization of experimental concentration data on the signaling network display.
In future releases we plan on expanding capabilities in all three areas of analysis – dynamic, structural, and experimental – with a focus on providing effective ways of integrating results from each together.
PathwayOracle has been designed in a modular fashion in order to facilitate extension of existing capabilities and the addition of new features.
Since PathwayOracle's most distinctive analytical capability involves the signaling Petri net simulator, a new dynamic analysis technique for signaling networks, we first provide an overview of the signaling Petri net modeling approach. Then in subsequent sections, we focus on PathwayOracle and explain the architecture and core concepts underlying the tool and then examine the individual features, how they can be used, and how they compare to existing tools.
The Signaling Petri Net Simulator
The signaling Petri net simulator models signal flow as the pattern of token accumulation and dissipation within proteins over time in the Petri net. Through transition firings, the source can influence the marking of (the number of tokens assigned to) the target, modeling the way that signals propagate through protein interactions in cellular signaling networks.
In order to overcome the issue of modeling reaction rates in the network, signaling dynamics are simulated by executing the signaling Petri net (SPN) for a set number of steps (called a run) multiple times, each time beginning at the same initial marking. For each run, the individual signaling rates are simulated via generation of random orders of transition firings (interaction occurrences). When the results of a large enough number of runs are averaged together, we find that the change in distribution of tokens in the network correlate with experimentally measured changes in the activity-levels of individual proteins in the underlying signaling network. In essence, the tokenized activity-levels computed by our method should be taken as abstract quantities whose changes over time correlate to changes that occur in the amounts of active proteins present in the cell. It is worth noting that some of the most widely used experimental techniques for protein quantification – western blots and microarrays – also yield results that are treated as indications, but not exact measurements, of protein activity-levels within the cell. Thus in some respects, the predictions returned by our SPN-based simulator can be interpreted like the results of a western blot or microarray experiment looking at changes relative to "control".
During a simulation run, the simulator imposes a strict ordering on transition firing such that it creates a two-time scale simulation. The smaller time scale is discretized as the firing of a single transition. This unit is referred to as the firing time scale. Firing steps are nested within a larger time scale, called time blocks, within which each transition is fired exactly once. The values returned by the simulator are the averaged token-counts for each protein at each time-block (across all runs).
In the first time block, transition t2 fires first: it reads 2 tokens out of Grb2 and places 2 additional tokens in Ras. Transition t1 fires second, reading 3 tokens out of Grb2. Transition t3 is evaluated last. The final marking for the network, highlighted as the red column in block 1 is used by the simulator as the marking for that block when averaging across runs.
At the conclusion of block 2, compare the values highlighted in red in the Initial column and at the end of both blocks. Note how the distribution of tokens have changed over the course of the simulation. Grb2 has the same number of tokens, implying that its activity-level has remained unchanged – this is consistent with the signaling network since no activating or inhibiting edges affect it in the model. AKT s token-count has risen, consistent with the fact that it is only activated in the signaling network. Ras's token-count has fallen which is one plausible behavior of the system since it is activated by Grb2, but inhibited by AKT.
PathwayOracle is written in Python . The user experience is oriented around visualization of and interaction with three main types of data: the signaling network, markings, and paths. At any given time, one signaling network is open, which is the basis for all analyses. Any simulation or concentration data is loaded and inspected as markings. Currently all static analyses revolve around paths, which are the third data type. In the following subsections, these individual data types and the user interfaces to them are discussed in more detail.
The Signaling Network Model
While the implementation of our methods use the signaling Petri net model discussed in an earlier section of this paper, we provide a simpler and more convenient representation of the network to the user which omits the internal topology of the transitions and allows the user to specify interactions simply as either activating or inhibiting. Thus, for the remainder of this paper we use the following definition of the signaling network which is consistent with the experience the user will have when working with PathwayOracle. The signaling network connectivity is a directed graph G = (V, E) where
• V is the set of nodes, which are signaling proteins and complexes (hereafter referred to collectively as signaling nodes) and
• E is the set of edges, which are signaling interactions. Each edge is of one of two types: u → v for activation and u ⊣ v for inhibition.
Within PathwayOracle, each signaling node has a name, unique within the network. A signaling edge has no properties besides its type and is only defined by its source and target.
u -> v or u - | v
where u is the name of the source signaling node and v is the name of the target signaling node. Each node is taken to represent the active form of the protein it is named for. Thus, from the example above, the interaction PI-3-K→AKT means that the active form of PI-3-K increases the activity-level of AKT whereas the interaction PTEN⊣AKT means that the active form of PTEN decreases the activity-level of AKT. While these types of unparameterized relationships can be represented in SBML, SBML was designed for encoding much more information than just connectivity . As a result, we deemed it appropriate to design a more concise format for our purposes. However, in a future release, PathwayOracle will support loading and saving in the SBML format.
At a given point in time, only one signaling network can be open in PathwayOracle. The main window displays a graphical representation of the network. The layout of the network can be modified by dragging nodes or by shift-clicking on edges to create, remove, or move waypoints. These layouts can be saved with the network and loaded again.
Signaling Network Markings
In signaling networks, signal flow is measured and quantified as the fluctuation of concentrations of various forms of signaling proteins over time. In PathwayOracle, we model concentrations using the concept of a network marking, which was adapted from Petri nets in which it was first used .
In PathwayOracle, a marking, μ is an assignment of real values to the nodes of a signaling network such that every signaling node receives a value. Earlier, the concept of a marking was introduced as the assignment of tokens to protein places in the signaling Petri net. In a signaling Petri net, tokens are discrete. In PathwayOracle, a marking is an average of the markings from many independent simulation runs, which gives rise to the real, rather than integral values, assigned by the marking.
As discussed earlier, the value of the marking of a signaling node, μ(v), can be interpreted as an estimate of the concentration or change in concentration of the active form of the signaling protein v (we call the amount of the active form of the signaling protein its activity-level). The two different versions of the simulator generate markings with these different meanings. The first simulator predicts the signal flow due to an experimental condition and generates markings whose values are taken to represent the actual activity-level of signaling protein present over the assumed basal levels. The second version of the simulator predicts the difference in signaling due to changing experimental conditions. The values assigned by markings produced by this simulator correspond to the change in the activity-level of the protein induced by the change in experimental condition. This will be discussed further in the Results and Discussion section.
In order to model signal flow, a single marking is not enough since it only provides a single snapshot of concentrations throughout the network. A marking series is an sequence of markings, (μ1, μ2,..., μ T ) in which the marking μ t is a snapshot of the concentration distribution at time step t. Thus, it is possible to see how the activity-level of protein v changed by plotting the values μ1(v), μ2(v),..., μ T (v). PathwayOracle provides the ability to do this.
In many experiments, the activity-level of various proteins are sampled at different time points and under different experimental conditions. Since the marking series is not able to represent changes due to different experimental conditions, we introduced the more general concept of a marking group in which each marking can correspond to an arbitrary activity-level distribution. Each marking is given a descriptive label that can be used to identify the conditions under which the activity-level was sampled.
Like the marking series, a marking group is loaded from a .csv file. However, unlike the marking series in which each row corresponds to a time step, in the marking group, each row corresponds to an independent marking (experimental condition). As shown in Figure 5(b), the first row is a header row specifying the molecule names for each column, the first column specifies the names for the individual markings (experimental conditions).
The Marking Manager
PathwayOracle includes a specific user-interface, the Marking Manager, designed to manage the three different types of markings. The Marking Manager provides a central interface within which it is possible to view all markings loaded and inspect them in ways that are relevant to their type (marking, marking series, or marking group). The specific ways in which markings can be inspected will be discussed further in the Results section.
The current structural analysis capabilities available in PathwayOracle allow inspection of signaling paths within the network. A signaling path p is a sequence of nodes, (v1, v2,..., v k ) where v i ∈ V ∀1 ≤ i ≤ k, and (v i , v i + 1) ∈ E ∀1 ≤ i <k. In this case, we say that node v 1 is the source of path p, and node v k is the target of p. Given a path, a variety of statistics may be of interest to the user. Additionally, it may be useful to view the path within the larger network. PathwayOracle provides these capabilites which will be discussed in the Results and Discussion section.
PathwayOracle provides a variety of tools for analyzing the structural and dynamic properties of a signaling network based on its connectivity. While its main differentiating feature is the ability to predict signal flow through a network using only the connectivity of the signaling network, PathwayOracle also provides the ability to visualize the network, analyze its connectivity, and inspect concentration-based experimental data.
Another important distinguishing characteristic of PathwayOracle is the combination of features that it supports. Biological network analysis is a multi-faceted process that may involve structural, dynamic, and data analysis in parallel. Whereas other tools tend to focus on one or two of these general areas of analysis, we considered it important for PathwayOracle to incorporate all three in order to provide the researcher a single environment in which all their analysis could be done. In future releases we plan to increase PathwayOracle's support for all three of these directions of investigation: structural, dynamic, and data analysis.
In the remainder of this section, we discuss the features currently available in PathwayOracle.
As in many other computational analysis tools for signaling networks (e.g., [20, 21]), an interactive graphical representation of the signaling network connectivity is at the center of the PathwayOracle interface. The main window provides a visualization of the signaling network connectivity. This visualization interface allows the user to edit the layout of the network by clicking on and dragging nodes and by shift-clicking on edges to create, remove, or move waypoints. Waypoints are points that lie on an edge. Holding down shift will display all edge waypoints. Existing waypoints can be dragged to change the path that an edge follows. Right-clicking on a waypoint will remove it. Left-clicking on a straight segment of the edge will create a new waypoint.
The network visualization also provides a view onto which path and experimental data analysis may be mapped. As will be discussed in subsequent sections, selected paths may be highlighted in this view and markings from experiments can set the colorings of individual nodes.
Network Signal Flow Simulation
The main feature differentiating PathwayOracle from other tools, such as CellDesigner  and COPASI , is its ability to simulate signal flow using an unparameterized signaling network model. Simulations can be performed in two different ways. In the first (Single Simulation), the simulator predicts the signal flow through the network for a specific experimental condition. In the second (Differential Simulation), the simulator predicts the difference in signal flow due to two different experimental conditions on the same network. These simulation methods themselves are described in . Here we focus on how simulations are configured, run, and analyzed.
Whereas the consensus networks typically represent the connectivity in normal cells, many experiments are conducted on abnormal cells in which oncogenic mutations, gene knockous, and pharmacological inhibitors have altered the behavior of various signaling nodes in the network. In PathwayOracle users can model these cell- and experiment-specific conditions by specifying each signaling node as either High, Low, or Free. The High state models any condition under which a protein's activity-level is held high for the duration of the experiment. This may be due to external stimulation or a known mutation in the protein that makes it constitutively active, for example. Similarly, a Low state models any phenomenon that forces a protein to have a persistently suppressed activity-level. This may be due to mutations that render the protein inactive, gene knockouts, or pharmacological inhibitors that force the activity-level of the protein low. In general, most signaling nodes will be Free, which means that their activity-level is unconstrained throughout the simulation. Only those nodes designated as High or Low will have their activity-level fixed for the duration of the simulation.
• the number of simulation runs to perform and
• the number of time blocks
The number of runs sets the number of independent simulations whose time block markings are averaged together to yield the overall simulation markings. In general, using more runs is a tradeoff between reliability of the results and simulation speed. In practice, the number of runs needed is dependent on the signaling network model and should be selected by observing the reproducability of the simulation results. An appropriate number of iterations will be large enough so that for a given experimental condition, the results are very similar across multiple simulations.
The time block, as discussed earlier, is a fundamental unit of time in the simulator. The appropriate number of time blocks for which to simulate will vary depending on the size of the signaling network and the scale of the network behavior of interest. Generally it should be selected by running simulations for a variety of time block values and determining which yields the most biologically reasonable activity-level changes for a known protein. While this is a manual process in the current version of PathwayOracle, we are investigating automated methods for estimating the number of time blocks by training against experimental time series data.
In PathwayOracle, the setup window for the Single Simulation (see Figure 8(a)) prompts the user for a single experimental condition. The setup window for the Differential Simulation (see Figure 8(b)) prompts the user for two experimental conditions. Both simulators produce a marking series. The tokenized simulation marking series corresponds to the activity-level time series predicted for the specified experimental condition. The differential simulation marking series corresponds to the change in activity-levels over time produced by switching from experimental condition 2 to experimental condition 1. The marking series produced by a simulation can be accessed through the Marking Manager. Choosing to inspect a marking series will present the user with a blank plot. By selecting signaling nodes, the plot is populated by the marking series values for individual nodes over time, as shown in Figure 8(c).
While this plot generation capability exists in many other dynamic simulation tools, the simplicity of the model used for simulation and the speed with which a simulation runs set PathwayOracle apart from other tools which require specification of the numerical values of kinetic parameters for each reaction in the network of interest (e.g., [20, 22]). PathwayOracle, because of its novel approach, does not have such requirements. It is worth noting, however, where PathwayOracle provides approximations of signal flow, an ODE generates the actual concentration changes using extremely detailed and accurate models of the underlying biochemistry. The simulators in PathwayOracle provide an attractive, time- and resource-saving alternative this more exhaustively parameterized techniques. In particular, PathwayOracle's features will benefit researchers interested in quickly assessing characteristics of signal flow in their network.
For some networks, biologists will have partial knowledge of kinetic parameters or of other biological details which the signaling Petri net model does not, at present, consider. By integrating this knowledge into the simulator, it may be possible to improve the simulator's predictions. We identify this as a direction for future investigation. As the signaling Petri net simulator is extended, these new capabilities will be incorporated in future releases of PathwayOracle.
Signaling Path Analysis
Experimental Data Analysis
Our goal is to develop PathwayOracle into an integrated and expansive suite of tools that allow the biologist to extract as much information as possible from models of signaling network connectivity and experimental data relating to those models. We consider future directions for PathwayOracle to fall into several categories: network construction, network augmentation, experimental and computational analysis integration, and architecture.
One of the benefits of working with connectivity models of signaling networks is the abundance of databases and other online resources that publish connectivity-level data. Future versions of PathwayOracle will have support for querying such databases for connectivity components and, ultimately, for automated connectivity construction based on a set of signaling nodes specified by the user.
Analysis of network connectivity and topology is increasingly relevant to biological research. We intend to expand PathwayOracle's structural analysis features to include the ability to search for and identify motifs in the signaling networks.
Network connectivity can also be inferred from experimental data, which provides another direction for research and development. By using experimental results to identify inconsistencies between experimental results and the current network model, it may be possible for PathwayOracle to augment the network with new connectivity based on hints supplied by experimental results. At present only experimental concentration data is supported. However, as experiments produce more information beyond concentration profiles of signaling nodes, we plan to expand the experimental data that PathwayOracle can load, visualize, and use as part of network analyses.
Experimental results can also provide computational analysis methods information that can improve their final predictions or decompositions. Taking advantage of the additional, potentially obfuscated, information present in experimental results to improve the results returned by computational tools is a major goal for future versions of PathwayOracle.
A longer term direction for PathwayOracle is the integration of transcriptional and metabolic network analysis. In the biological systems of interest, the behavior of any one of these networks is dependent on the characteristics of the other two. As a result, developing a complete understanding of signaling, transcriptional regulation, or metabolism depends in part on integrating knowledge from the others. Finally, an ongoing priority in the design of PathwayOracle is its role as an open platform for the development and deployment of new analytical capabilities by other groups. Currently PathwayOracle employes a modular architecture that facilitates easy integration of new functionality. However, in future releases we plan to expose a plugin interface which will make it easier to developers and researchers to develop and deploy tools within PathwayOracle.
PathwayOracle is an integrated software environment in which biologists may conduct structural and dynamic analysis of signaling networks of interest. PathwayOracle is distinguished from other tools in the field of systems biology by its ability to predict the signal flow through a network using a simplified, connectivity-based model of the signaling network. Simulations are fast and, based on a published study, predictors of signal propagation. This novel simulation capability, combined with support for structural analysis of connectivity between pairs of proteins and for analysis of certain kinds of experimental data make PathwayOracle a powerful asset in the experimentalist's endeavor to gain a more complete understanding of the cellular signaling network.
Project name: PathwayOracle
Project home page: http://bioinfo.cs.rice.edu/pathwayoracle
Operating system(s): Platform independent
Programming language: Python
Other requirements: Python 2.4 or higher
License: GNU GPL
Any restrictions to use by non-academics: None
DR and LN are supported in part by a Seed Grant awarded to LN from the Gulf Coast Center for Computational Cancer Research, funded by John and Ann Doerr Fund for Computational Biomedicine. PTR is supported in part by a Department of Defense grant BC044268.
- Papin JA, Hunter T, Palsson BO, Subramaniam S: Reconstruction of cellular signalling networks and analysis of their properties. Nature Reviews Molecular Cell Biology. 2005, 6: 99-111. 10.1038/nrm1570View ArticlePubMedGoogle Scholar
- Sackmann A, Heiner M, Koch I: Application of Petri net based analysis techniques to signal transduction pathways. BMC Bioinformatics. 2006, 7: 482- 10.1186/1471-2105-7-482PubMed CentralView ArticlePubMedGoogle Scholar
- Ruths D, Tseng JT, Nakhleh L, Ram PT: De novo Signaling Pathway Predictions based on Protein-Protein Interaction, Targeted Therapy, and Protein Microarray Analysis. Proceedings of the RECOMB Satellite Workshop on Systems Biology and Proteomics, Lecture Notes in Bioinformatics (LNBI #4466). 2007, 62-72.Google Scholar
- Ruths D, Nakhleh L, Iyengar MS, Reddy SAG, Ram PT: Graph-theoretic Hypothesis Generation in Biological Signaling Networks. Journal of Computational Biology. 2006, 13 (9): 1546-1557. 10.1089/cmb.2006.13.1546View ArticlePubMedGoogle Scholar
- Klamt S, Saez-Rodriguez J, Lindquist JA, Simeoni L, Gilles ED: A methodology for the structural and functional analysis of signaling and regulatory networks. BMC Bioinformatics. 2006, 6: 56-10.1186/1471-2105-7-56.View ArticleGoogle Scholar
- Papin JA, Palsson BO: The JAK-STAT signaling network in the human B-cell: an extreme signaling pathway analysis. Biophysical Journal. 2004, 87: 37-46. 10.1529/biophysj.103.029884PubMed CentralView ArticlePubMedGoogle Scholar
- Papin JA, Price ND, Wiback SJ, Fell DA, Palsson BO: Metabolic pathways in the post-genomic era. TRENDS in Biochemical Sciences. 2003, 28 (5): 250-258. 10.1016/S0968-0004(03)00064-1View ArticlePubMedGoogle Scholar
- Schilling CH, Letscher D, Palsson BO: Theory for the Systemic Definition of Metabolic Pathways and their use in Interpreting Metabolic Function from a Pathway-Oriented Perspective. Journal of Theoretical Biology. 2000, 203: 229-248. 10.1006/jtbi.2000.1073View ArticlePubMedGoogle Scholar
- Peleg M, Rubin D, Altman RB: Using Petri Net Tools to Study Properties and Dynamics of Biological Systems. Journal of the American Medical Informatics Association. 2005, 12 (2): 181-199. 10.1197/jamia.M1637PubMed CentralView ArticlePubMedGoogle Scholar
- Chaouiya C: Petri net modelling of biological networks. Briefings in Bioinformatics. 2007, 8 (4): 210-219. 10.1093/bib/bbm029View ArticlePubMedGoogle Scholar
- Eungdamrong NJ, Iyengar R: Modeling cell signaling networks. Biology of the Cell. 2004, 96 (5): 355-362. 10.1016/j.biolcel.2004.03.004PubMed CentralView ArticlePubMedGoogle Scholar
- Eungdamrong NJ, Iyengar R: Computational Approaches for modeling regulatory cellular networks. Trends Cell Biology. 2004, 14 (12): 661-669. 10.1016/j.tcb.2004.10.007View ArticleGoogle Scholar
- Ruths D, Muller M, Tseng JT, Nakhleh L, Ram PT: The Signaling Petri Net-based Simulator: A Non-Parametric Strategy for Characterizing the Dynamics of Cell-Specific Signaling Networks. PLoS Computational Biology. 2008, 4 (2): e1000005- 10.1371/journal.pcbi.1000005PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research. 2000, 28 (15): 27-30. 10.1093/nar/28.1.27PubMed CentralView ArticlePubMedGoogle Scholar
- The Cancer Cell Map. http://cancer.cellmap.org
- Thomas PD, Campbell MJ, Kejariwal A, Mi J, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A: PANTHER: a library of protein families and subfamilies indexed by function. Genome Research. 2003, 13: 2129-2141. 10.1101/gr.772403PubMed CentralView ArticlePubMedGoogle Scholar
- David R, Alla H: Discrete, Continuous, and Hybrid Petri Nets. 2005, SpringerGoogle Scholar
- Official Website for Python Programming Language. http://www.python.org
- Hucka M, Finney A, Sauro HM, H B, Doyle J, Kitano H, Arkin A, Bornstein B, Bray D, Cornish-Bowden A, Cuellar A, Dronov S, Gilles E, Ginkel M, Gor V, Goryanin I, Hedley W, Hodgman T, Hofmeyr J, Hunter P, Juty N, Kasberger J, Kremling A, Kummer U, Le Novere N, Loew L, Lucio D, Mendes P, Minch E, Mjolsness E, Nakayama Y, Nelson M, Nielsen P, Sakurada T, Schaff J, Shapiro B, Shimizu T, Spence H, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J: The Systems Biology Markup Language (SBML): A medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19 (4): 524-531. 10.1093/bioinformatics/btg015View ArticlePubMedGoogle Scholar
- Funahashi A, Tanimura N, Morohashi M, Kitano H: CellDesigner: a process diagram editor for gene-regulatory and biochemical networks. BIOSILICO. 2003, 1: 159-162. 10.1016/S1478-5382(03)02370-9.View ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Research. 2003, 13 (11): 2498-504. 10.1101/gr.1239303PubMed CentralView 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: 3067-74. 10.1093/bioinformatics/btl485View ArticlePubMedGoogle Scholar
- Alon U: An Introduction to Systems Biology: Design Principles of Biological Circuits. Mathematical and Computational Biology Series. 2007, Chapman & Hall/CRCGoogle Scholar
- Nagasaki M, Doi A, Matsuno H, Miyano S: Genomic Object Net: I. A platform for modelling and simulating biopathways. Applied Bioinformatics. 2003, 2 (3): 181-184.PubMedGoogle Scholar
- Klamt S, Saez-Rodriguez J, Gilles ED: Structural and functional analysis of cellular networks with CellNetAnalyzer. BMC Systems Biology. 2007, 1: 2- 10.1186/1752-0509-1-2PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidt H, Jirstand M: Systems Biology Toolbox for MATLAB: A computational platform for research in Systems Biology. Bioinformatics. 2006, 22 (4): 514-515. 10.1093/bioinformatics/bti799View ArticlePubMedGoogle Scholar
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.