SimBoolNet--a Cytoscape plugin for dynamic simulation of signaling networks

In silico analysis and experimental validation to exhibit anti-nasopharyngeal carcinoma effects of plumbagin, an anti-cancer compound

BACKGROUND

Nasopharyngeal carcinoma (NPC) is publicly known as a malignant tumor. Our previous study reported that plumbagin exhibits potent anti-cancer actions. Nevertheless, more mechanical details of plumbagin against NPC remain unknown. The present study aimed to unmask the core targets/genes and anti-NPC mechanisms involved in the signaling pathways of plumbagin prior to biochemical validation.

METHODS

A network pharmacology approach was employed to respective identification of mutual and core targets/genes in plumbagin and/treating NPC. Molecular docking determination was used to identify core target proteins for biochemical validation using human and cell line samples.

RESULTS

In total, 60 anti-NPC genes of plumbagin were screened out, and then nine core target genes of plumbagin against NPC were identified accordingly. The enrichment findings revealed detailed biological functions and pharmacological pathways of plumbagin against NPC. Moreover, in silico analysis using molecular docking had determined the core targets for further experimental validation, comprising protein kinase B (AKT1) and sarcoma gene (SRC). In human sample validation, clinical NPC sections showed increased positive expression of AKT1 and SRC. Additionally, plumbagin-treated NPC cells resulted in inactivated protein expression of AKT1 and SRC.

CONCLUSION

The re-identified core targets/genes in the molecular docking report may function as plumbagin-related pharmacological targets for treating NPC via experimental validation. Furthermore, additional anti-NPC molecular mechanisms of plumbagin action were disclosed on the basis of enrichment findings. © 2022 Society of Chemical Industry.

SteadyCellPhenotype: a web-based tool for the modeling of biological networks with ternary logic

SUMMARY

We introduce SteadyCellPhenotype, a browser-based interface for the analysis of ternary biological networks. It includes tools for deterministically finding all steady states of a network, as well as the simulation and visualization of trajectories with publication quality graphics. Simulations allow us to approximate the size of the basin for attractors and deterministic simulations of trajectories nearby specified points allow us to explore the behavior of the system in that neighborhood.

AVAILABILITY AND IMPLEMENTATION

https://github.com/knappa/steadycellphenotype MIT License.

Boolink: a graphical interface for open access Boolean network simulations and use in guard cell CO2 signaling

Signaling networks are at the heart of almost all biological processes. Most of these networks contain large number of components, and often either the connections between these components are not known or the rate equations that govern the dynamics of soluble signaling components are not quantified. This uncertainty in network topology and parameters can make it challenging to formulate detailed mathematical models. Boolean networks, in which all components are either on or off, have emerged as viable alternatives to detailed mathematical models that contain rate constants and other parameters. Therefore, open-source platforms of Boolean models for community use are desirable. Here, we present Boolink, a freely available graphical user interface that allows users to easily construct and analyze existing Boolean networks. Boolink can be applied to any Boolean network. We demonstrate its application using a previously published network for abscisic acid (ABA)-driven stomatal closure in Arabidopsis spp. (Arabidopsis thaliana). We also show how Boolink can be used to generate testable predictions by extending the network to include CO2 regulation of stomatal movements. Predictions of the model were experimentally tested, and the model was iteratively modified based on experiments showing that ABA effectively closes Arabidopsis stomata at near-zero CO2 concentrations (1.5-ppm CO2). Thus, Boolink enables public generation and the use of existing Boolean models, including the prior developed ABA signaling model with added CO2 signaling components.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Capturing dynamic relevance in Boolean networks using graph theoretical measures

Motivation

Interaction graphs are able to describe regulatory dependencies between compounds without capturing dynamics. In contrast, mathematical models that are based on interaction graphs allow to investigate the dynamics of biological systems. However, since dynamic complexity of these models grows exponentially with their size, exhaustive analyses of the dynamics and consequently screening all possible interventions eventually becomes infeasible. Thus, we designed an approach to identify dynamically relevant compounds based on the static network topology.

Results

Here, we present a method only based on static properties to identify dynamically influencing nodes. Coupling vertex betweenness and determinative power, we could capture relevant nodes for changing dynamics with an accuracy of 75% in a set of 35 published logical models. Further analyses of the selected compounds’ connectivity unravelled a new class of not highly connected nodes with high impact on the networks’ dynamics, which we call gatekeepers. We validated our method’s working concept on logical models, which can be readily scaled up to complex interaction networks, where dynamic analyses are not even feasible.

Availability and implementation

Code is freely available at https://github.com/sysbio-bioinf/BNStatic.

Supplementary information

Supplementary data are available at Bioinformatics online.

Concepts in Boolean network modeling: What do they all mean?

Boolean network models are one of the simplest models to study complex dynamic behavior in biological systems. They can be applied to unravel the mechanisms regulating the properties of the system or to identify promising intervention targets. Since its introduction by Stuart Kauffman in 1969 for describing gene regulatory networks, various biologically based networks and tools for their analysis were developed. Here, we summarize and explain the concepts for Boolean network modeling. We also present application examples and guidelines to work with and analyze Boolean network models.

Integrative workflows for network analysis

Due to genetic heterogeneity across patients, the identification of effective disease signatures and therapeutic targets is challenging. Addressing this challenge, we have previously developed a network-based approach, which integrates heterogeneous sources of biological information to identify disease specific core-regulatory networks. In particular, our workflow uses a multi-objective optimization function to calculate a ranking score for network components (e.g. feedback/feedforward loops) based on network properties, biomedical and high-throughput expression data. High ranked network components are merged to identify the core-regulatory network(s) that is then subjected to dynamical analysis using stimulus–response and in silico perturbation experiments for the identification of disease gene signatures and therapeutic targets. In a case study, we implemented our workflow to identify bladder and breast cancer specific core-regulatory networks underlying epithelial–mesenchymal transition from the E2F1 molecular interaction map.

In this study, we review our workflow and described how it has developed over time to understand the mechanisms underlying disease progression and prediction of signatures for clinical decision making.

ViSiBooL-visualization and simulation of Boolean networks with temporal constraints

Summary

Mathematical models and their simulation are increasingly used to gain insights into cellular pathways and regulatory networks. Dynamics of regulatory factors can be modeled using Boolean networks (BNs), among others. Text-based representations of models are precise descriptions, but hard to understand and interpret. ViSiBooL aims at providing a graphical way of modeling and simulating networks. By providing visualizations of static and dynamic network properties simultaneously, it is possible to directly observe the effects of changes in the network structure on the behavior. In order to address the challenges of clear design and a user-friendly graphical user interface (GUI), ViSiBooL implements visual representations of BNs. Additionally temporal extensions of the BNs for the modeling of regulatory time delays are incorporated. The GUI of ViSiBooL allows to model, organize, simulate and visualize BNs as well as corresponding simulation results such as attractors. Attractor searches are performed in parallel to the modeling process. Hence, changes in the network behavior are visualized at the same time.

Availability and Implementation

ViSiBooL (Java 8) is freely available at http://sysbio.uni-ulm.de/?Software:ViSiBooL .

Contact

hans.kestler@uni-ulm.de.

A Network-Based Integrative Workflow to Unravel Mechanisms Underlying Disease Progression

Unraveling mechanisms underlying diseases has motivated the development of systems biology approaches. The key challenges for the development of mathematical models and computational tool are (1) the size of molecular networks, (2) the nonlinear nature of spatio-temporal interactions, and (3) feedback loops in the structure of interaction networks. We here propose an integrative workflow that combines structural analyses of networks, high-throughput data, and mechanistic modeling. As an illustration of the workflow, we use prostate cancer as a case study with the aim of identifying key functional components associated with primary to metastasis transitions. Analysis carried out by the workflow revealed that HOXD10, BCL2, and PGR are the most important factors affected in primary prostate samples, whereas, in the metastatic state, STAT3, JUN, and JUNB are playing a central role. The identified key elements of each network are validated using patient survival analysis. The workflow presented here allows experimentalists to use heterogeneous data sources for the identification of diagnostic and prognostic signatures.

Simulation and visualization of multiple KEGG pathways using BioNSi

Motivation: Many biologists are discouraged from using network simulation tools because these require manual, often tedious network construction. This situation calls for building new tools or extending existing ones with the ability to import biological pathways previously deposited in databases and analyze them, in order to produce novel biological insights at the pathway level.

Results: We have extended a network simulation tool (BioNSi), which now allows merging of multiple pathways from the KEGG pathway database into a single, coherent network, and visualizing its properties. Furthermore, the enhanced tool enables loading experimental expression data into the network and simulating its dynamics under various biological conditions or perturbations. As a proof of concept, we tested two sets of published experimental data, one related to inflammatory bowel disease condition and the other to breast cancer treatment. We predict some of the major observations obtained following these laboratory experiments, and provide new insights that may shed additional light on these results.

Tool requirements: Cytoscape 3.x, JAVA 8

Availability: The tool is freely available at http://bionsi.wix.com/bionsi, where a complete user guide and a step-by-step manual can also be found.

Advanced Boolean modeling of biological networks applied to systems pharmacology

Motivation

Literature on complex diseases is abundant but not always quantitative. Many molecular pathways are qualitatively well described but this information cannot be used in traditional quantitative mathematical models employed in drug development. Tools for analysis of discrete networks are useful to capture the available information in the literature but have not been efficiently integrated by the pharmaceutical industry. We propose an expansion of the usual analysis of discrete networks that facilitates the identification/validation of therapeutic targets.

Results

In this article, we propose a methodology to perform Boolean modeling of Systems Biology/Pharmacology networks by using SPIDDOR (Systems Pharmacology for effIcient Drug Development On R) R package. The resulting models can be used to analyze the dynamics of signaling networks associated to diseases to predict the pathogenesis mechanisms and identify potential therapeutic targets.

Availability and Implementation

The source code is available at https://github.com/SPIDDOR/SPIDDOR .

Contact

itzirurzun@alumni.unav.es , itroconiz@unav.es.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sig2GRN: a software tool linking signaling pathway with gene regulatory network for dynamic simulation

Background

Linking computational models of signaling pathways to predicted cellular responses such as gene expression regulation is a major challenge in computational systems biology. In this work, we present Sig2GRN, a Cytoscape plugin that is able to simulate time-course gene expression data given the user-defined external stimuli to the signaling pathways.

Methods

A generalized logical model is used in modeling the upstream signaling pathways. Then a Boolean model and a thermodynamics-based model are employed to predict the downstream changes in gene expression based on the simulated dynamics of transcription factors in signaling pathways.

Results

Our empirical case studies show that the simulation of Sig2GRN can predict changes in gene expression patterns induced by DNA damage signals and drug treatments.

Conclusions

As a software tool for modeling cellular dynamics, Sig2GRN can facilitate studies in systems biology by hypotheses generation and wet-lab experimental design.

Availability: http://histone.scse.ntu.edu.sg/Sig2GRN/

Power-Law Modeling of Cancer Cell Fates Driven by Signaling Data to Reveal Drug Effects

Extracellular signals are captured and transmitted by signaling proteins inside a cell. An important type of cellular responses to the signals is the cell fate decision, e.g., apoptosis. However, the underlying mechanisms of cell fate regulation are still unclear, thus comprehensive and detailed kinetic models are not yet available. Alternatively, data-driven models are promising to bridge signaling data with the phenotypic measurements of cell fates. The traditional linear model for data-driven modeling of signaling pathways has its limitations because it assumes that the a cell fate is proportional to the activities of signaling proteins, which is unlikely in the complex biological systems. Therefore, we propose a power-law model to relate the activities of all the measured signaling proteins to the probabilities of cell fates. In our experiments, we compared our nonlinear power-law model with the linear model on three cancer datasets with phosphoproteomics and cell fate measurements, which demonstrated that the nonlinear model has superior performance on cell fates prediction. By in silico simulation of virtual protein knock-down, the proposed model is able to reveal drug effects which can complement traditional approaches such as binding affinity analysis. Moreover, our model is able to capture cell line specific information to distinguish one cell line from another in cell fate prediction. Our results show that the power-law data-driven model is able to perform better in cell fate prediction and provide more insights into the signaling pathways for cancer cell fates than the linear model.

Properties of Boolean dynamics by node classification using feedback loops in a network

Background

Biological networks keep their functions robust against perturbations. Many previous studies through simulations or experiments have shown that feedback loop (FBL) structures play an important role in controlling the network robustness without fully explaining how they do it. Hence, there is a pressing need to more rigorously analyze the influence of FBL structures on network robustness.

Results

In this paper, I propose a novel node classification notion based on the FBL structures involved. More specifically, I classify a node as a no-FBL-in-upstream (NFU) or no-FBL-in-downstream (NFD) node if no feedback loop is involved with any upstream or downstream path of the node, respectively. Based on those definitions, I first prove that every NFU node is eventually frozen in Boolean dynamics. Thus, NFU nodes converge to a fixed value determined by the upstream source nodes. Second, I prove that a network is robust against an arbitrary state perturbation subject to a non-source NFD node. This implies that a network state eventually sustains the attractor despite a perturbation subject to a non-source NFD node. Inspired by this result, I further propose a perturbation-sustainable probability that indicates how likely a perturbation effect is to be sustained through propagations. I show that genes with a high perturbation-sustainable probability are likely to be essential, disease, and drug-target genes in large human signaling networks.

Conclusion

Taken together, these results will promote understanding of the effects of FBL on network robustness in a more rigorous manner.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0322-z) contains supplementary material, which is available to authorized users.

BioNSi: A Discrete Biological Network Simulator Tool

Modeling and simulation of biological networks is an effective and widely used research methodology. The Biological Network Simulator (BioNSi) is a tool for modeling biological networks and simulating their discrete-time dynamics, implemented as a Cytoscape App. BioNSi includes a visual representation of the network that enables researchers to construct, set the parameters, and observe network behavior under various conditions. To construct a network instance in BioNSi, only partial, qualitative biological data suffices. The tool is aimed for use by experimental biologists and requires no prior computational or mathematical expertise. BioNSi is freely available at http://bionsi.wix.com/bionsi , where a complete user guide and a step-by-step manual can also be found.

CABeRNET: a Cytoscape app for augmented Boolean models of gene regulatory NETworks

BACKGROUND

Dynamical models of gene regulatory networks (GRNs) are highly effective in describing complex biological phenomena and processes, such as cell differentiation and cancer development. Yet, the topological and functional characterization of real GRNs is often still partial and an exhaustive picture of their functioning is missing.

RESULTS

We here introduce CABERNET, a Cytoscape app for the generation, simulation and analysis of Boolean models of GRNs, specifically focused on their augmentation when a only partial topological and functional characterization of the network is available. By generating large ensembles of networks in which user-defined entities and relations are added to the original core, CABERNET allows to formulate hypotheses on the missing portions of real networks, as well to investigate their generic properties, in the spirit of complexity science.

CONCLUSIONS

CABERNET offers a series of innovative simulation and modeling functions and tools, including (but not being limited to) the dynamical characterization of the gene activation patterns ruling cell types and differentiation fates, and sophisticated robustness assessments, as in the case of gene knockouts. The integration within the widely used Cytoscape framework for the visualization and analysis of biological networks, makes CABERNET a new essential instrument for both the bioinformatician and the computational biologist, as well as a computational support for the experimentalist. An example application concerning the analysis of an augmented T-helper cell GRN is provided.

Generalized logical model based on network topology to capture the dynamical trends of cellular signaling pathways

Background

Cellular responses to extracellular perturbations require signaling pathways to capture and transmit the signals. However, the underlying molecular mechanisms of signal transduction are not yet fully understood, thus detailed and comprehensive models may not be available for all the signaling pathways. In particular, insufficient knowledge of parameters, which is a long-standing hindrance for quantitative kinetic modeling necessitates the use of parameter-free methods for modeling and simulation to capture dynamic properties of signaling pathways.

Results

We present a computational model that is able to simulate the graded responses to degradations, the sigmoidal biological relationships between signaling molecules and the effects of scheduled perturbations to the cells. The simulation results are validated using experimental data of protein phosphorylation, demonstrating that the proposed model is capable of capturing the main trend of protein activities during the process of signal transduction. Compared with existing simulators, our model has better performance on predicting the state transitions of signaling networks.

Conclusion

The proposed simulation tool provides a valuable resource for modeling cellular signaling pathways using a knowledge-based method.

Open source libraries and frameworks for biological data visualisation: a guide for developers

Recent advances in high-throughput experimental techniques have led to an exponential increase in both the size and the complexity of the data sets commonly studied in biology. Data visualisation is increasingly used as the key to unlock this data, going from hypothesis generation to model evaluation and tool implementation. It is becoming more and more the heart of bioinformatics workflows, enabling scientists to reason and communicate more effectively. In parallel, there has been a corresponding trend towards the development of related software, which has triggered the maturation of different visualisation libraries and frameworks. For bioinformaticians, scientific programmers and software developers, the main challenge is to pick out the most fitting one(s) to create clear, meaningful and integrated data visualisation for their particular use cases. In this review, we introduce a collection of open source or free to use libraries and frameworks for creating data visualisation, covering the generation of a wide variety of charts and graphs. We will focus on software written in Java, JavaScript or Python. We truly believe this software offers the potential to turn tedious data into exciting visual stories.

Toward an integrated software platform for systems pharmacology

Understanding complex biological systems requires the extensive support of computational tools. This is particularly true for systems pharmacology, which aims to understand the action of drugs and their interactions in a systems context. Computational models play an important role as they can be viewed as an explicit representation of biological hypotheses to be tested. A series of software and data resources are used for model development, verification and exploration of the possible behaviors of biological systems using the model that may not be possible or not cost effective by experiments. Software platforms play a dominant role in creativity and productivity support and have transformed many industries, techniques that can be applied to biology as well. Establishing an integrated software platform will be the next important step in the field. © 2013 The Authors. Biopharmaceutics & Drug Disposition published by John Wiley & Sons, Ltd.

SBML qualitative models: a model representation format and infrastructure to foster interactions between qualitative modelling formalisms and tools

BACKGROUND

Qualitative frameworks, especially those based on the logical discrete formalism, are increasingly used to model regulatory and signalling networks. A major advantage of these frameworks is that they do not require precise quantitative data, and that they are well-suited for studies of large networks. While numerous groups have developed specific computational tools that provide original methods to analyse qualitative models, a standard format to exchange qualitative models has been missing.

RESULTS

We present the Systems Biology Markup Language (SBML) Qualitative Models Package ("qual"), an extension of the SBML Level 3 standard designed for computer representation of qualitative models of biological networks. We demonstrate the interoperability of models via SBML qual through the analysis of a specific signalling network by three independent software tools. Furthermore, the collective effort to define the SBML qual format paved the way for the development of LogicalModel, an open-source model library, which will facilitate the adoption of the format as well as the collaborative development of algorithms to analyse qualitative models.

CONCLUSIONS

SBML qual allows the exchange of qualitative models among a number of complementary software tools. SBML qual has the potential to promote collaborative work on the development of novel computational approaches, as well as on the specification and the analysis of comprehensive qualitative models of regulatory and signalling networks.

Boolean modeling of biological regulatory networks: a methodology tutorial

Given the complexity and interactive nature of biological systems, constructing informative and coherent network models of these systems and subsequently developing efficient approaches to analyze the assembled networks is of immense importance. The integration of network analysis and dynamic modeling enables one to investigate the behavior of the underlying system as a whole and to make experimentally testable predictions about less-understood aspects of the processes involved. In this paper, we present a tutorial on the fundamental steps of Boolean modeling of biological regulatory networks. We demonstrate how to infer a Boolean network model from the available experimental data, analyze the network using graph-theoretical measures, and convert it into a predictive dynamic model. For each step, the pitfalls one may encounter and possible ways to circumvent them are also discussed. We illustrate these steps on a toy network as well as in the context of the Drosophila melanogaster segment polarity gene network.

Dynamics of influenza virus and human host interactions during infection and replication cycle

The replication and life cycle of the influenza virus is governed by an intricate network of intracellular regulatory events during infection, including interactions with an even more complex system of biochemical interactions of the host cell. Computational modeling and systems biology have been successfully employed to further the understanding of various biological systems, however, computational studies of the complexity of intracellular interactions during influenza infection is lacking. In this work, we present the first large-scale dynamical model of the infection and replication cycle of influenza, as well as some of its interactions with the host's signaling machinery. Specifically, we focus on and visualize the dynamics of the internalization and endocytosis of the virus, replication and translation of its genomic components, as well as the assembly of progeny virions. Simulations and analyses of the models dynamics qualitatively reproduced numerous biological phenomena discovered in the laboratory. Finally, comparisons of the dynamics of existing and proposed drugs, our results suggest that a drug targeting PB1:PA would be more efficient than existing Amantadin/Rimantaine or Zanamivir/Oseltamivir.

Bio-logic builder: a non-technical tool for building dynamical, qualitative models

Computational modeling of biological processes is a promising tool in biomedical research. While a large part of its potential lies in the ability to integrate it with laboratory research, modeling currently generally requires a high degree of training in mathematics and/or computer science. To help address this issue, we have developed a web-based tool, Bio-Logic Builder, that enables laboratory scientists to define mathematical representations (based on a discrete formalism) of biological regulatory mechanisms in a modular and non-technical fashion. As part of the user interface, generalized “bio-logic” modules have been defined to provide users with the building blocks for many biological processes. To build/modify computational models, experimentalists provide purely qualitative information about a particular regulatory mechanisms as is generally found in the laboratory. The Bio-Logic Builder subsequently converts the provided information into a mathematical representation described with Boolean expressions/rules. We used this tool to build a number of dynamical models, including a 130-protein large-scale model of signal transduction with over 800 interactions, influenza A replication cycle with 127 species and 200+ interactions, and mammalian and budding yeast cell cycles. We also show that any and all qualitative regulatory mechanisms can be built using this tool.

Analysis of cell adhesion during early stages of colon cancer based on an extended multi-valued logic approach

Cell adhesion in the normal colon is typically associated with differentiated cells, whereas in cancerous colon it is associated with advanced tumors. For advanced tumors growing evidence supports the existence of stem-like cells that have originated from transdifferentiation. Because stem cells can also be transformed in their own niche, at the base of the Lieberkühn's crypts, we conjectured that cell adhesion can also be critical in early tumorigenesis. To assess this hypothesis we built an annotated, multi-valued logic model addressing cell adhesion of normal and tumorigenic stem cells in the human colon. The model accounts for (i) events involving intercellular adhesion structures, (ii) interactions involving cytoskeleton-related structures, (iii) compartmental distribution of α/β/γ/δ-catenins, and (iv) variations in critical cell adhesion regulators (e.g., ILK, FAK, IQGAP, SNAIL, Caveolin). We developed a method that can deal with graded multiple inhibitions, something which is not possible with conventional logical approaches. The model comprises 315 species (including 26 genes), interconnected by 269 reactions. Simulations of the model covered six scenarios, which considered two types of colonic cells (stem vs. differentiated cells), under three conditions (normal, stressed and tumor). Each condition results from the combination of 92 inputs. We compared our multi-valued logic approach with the conventional Boolean approach for one specific example and validated the predictions against published data. Our analysis suggests that stem cells in their niche synthesize high levels of cytoplasmatic E-cadherin and CdhEP(Ser684,686,692), even under normal-mitogenic stimulus or tumorigenic conditions. Under these conditions, E-cadherin would be incorporated into the plasmatic membrane, but only as a non-adhesive CdhE_β-catenin_IQGAP complex. Under stress conditions, however, this complex could be displaced, yielding adhesive CdhE_β-catenin((cis/trans)) complexes. In the three scenarios tested with stem cells, desmosomes or tight junctions were not assembled. Other model predictions include expected levels of the nuclear complex β-catenin_TCF4 and the anti-apoptotic protein Survivin for both normal and tumorigenic colonic stem cells.

Dynamic network-based epistasis analysis: boolean examples

In this article we focus on how the hierarchical and single-path assumptions of epistasis analysis can bias the inference of gene regulatory networks. Here we emphasize the critical importance of dynamic analyses, and specifically illustrate the use of Boolean network models. Epistasis in a broad sense refers to gene interactions, however, as originally proposed by Bateson, epistasis is defined as the blocking of a particular allelic effect due to the effect of another allele at a different locus (herein, classical epistasis). Classical epistasis analysis has proven powerful and useful, allowing researchers to infer and assign directionality to gene interactions. As larger data sets are becoming available, the analysis of classical epistasis is being complemented with computer science tools and system biology approaches. We show that when the hierarchical and single-path assumptions are not met in classical epistasis analysis, the access to relevant information and the correct inference of gene interaction topologies is hindered, and it becomes necessary to consider the temporal dynamics of gene interactions. The use of dynamical networks can overcome these limitations. We particularly focus on the use of Boolean networks that, like classical epistasis analysis, relies on logical formalisms, and hence can complement classical epistasis analysis and relax its assumptions. We develop a couple of theoretical examples and analyze them from a dynamic Boolean network model perspective. Boolean networks could help to guide additional experiments and discern among alternative regulatory schemes that would be impossible or difficult to infer without the elimination of these assumption from the classical epistasis analysis. We also use examples from the literature to show how a Boolean network-based approach has resolved ambiguities and guided epistasis analysis. Our article complements previous accounts, not only by focusing on the implications of the hierarchical and single-path assumption, but also by demonstrating the importance of considering temporal dynamics, and specifically introducing the usefulness of Boolean network models and also reviewing some key properties of network approaches.

Software for systems biology: from tools to integrated platforms

Understanding complex biological systems requires extensive support from software tools. Such tools are needed at each step of a systems biology computational workflow, which typically consists of data handling, network inference, deep curation, dynamical simulation and model analysis. In addition, there are now efforts to develop integrated software platforms, so that tools that are used at different stages of the workflow and by different researchers can easily be used together. This Review describes the types of software tools that are required at different stages of systems biology research and the current options that are available for systems biology researchers. We also discuss the challenges and prospects for modelling the effects of genetic changes on physiology and the concept of an integrated platform.

NetDS: a Cytoscape plugin to analyze the robustness of dynamics and feedforward/feedback loop structures of biological networks

SUMMARY

NetDS is a novel Cytoscape plugin that conveniently simulates dynamics related to robustness, and examines structural properties with respect to feedforward/feedback loops. It can evaluate how robustly a network sustains a stable state against mutations by employing a Boolean network model. In addition, the plugin can examine all feedforward/feedback loops appearing in a network and determine whether or not a pair of loops is coupled. Random networks can also be generated to evaluate whether or not an interesting finding in real biological networks is significantly random.

AVAILABILITY

NetDS is freely available for non-commercial purposes at http://netds.sourceforge.net/.

CONTACT

kwonyk@ulsan.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Visualization in simulation tools: requirements and a tool specification to support the teaching of dynamic biological processes

Simulation tools are playing an increasingly important role behind advances in the field of systems biology. However, the current generation of biological science students has either little or no experience with such tools. As such, this educational glitch is limiting both the potential use of such tools as well as the potential for tighter cooperation between the designers and users. Although some simulation tool producers encourage their use in teaching, little attempt has hitherto been made to analyze and discuss their suitability as an educational tool for noncomputing science students. In general, today's simulation tools assume that the user has a stronger mathematical and computing background than that which is found in most biological science curricula, thus making the introduction of such tools a considerable pedagogical challenge. This paper provides an evaluation of the pedagogical attributes of existing simulation tools for cell signal transduction based on Cognitive Load theory. Further, design recommendations for an improved educational simulation tool are provided. The study is based on simulation tools for cell signal transduction. However, the discussions are relevant to a broader biological simulation tool set.