Conservation analysis of large biochemical networks

MOTIVATION: Large biochemical networks pose a unique challenge from the point of view of evaluating conservation laws. The computational problem in most cases exceeds the capability of available software tools, often resulting in inaccurate computation of the number and form of conserved cycles. Such errors have profound effects on subsequent calculations, particularly in the evaluation of the Jacobian which is a critical quantity in many other calculations. The goal of this paper is to outline a new algorithm that is computationally efficient and robust at extracting the correct conservation laws for very large biochemical networks. RESULTS: We show that our algorithm can perform the conservation analysis of large biochemical networks, and can evaluate the correct conserved cycles when compared with other similar software tools. Biochemical simulators such as Jarnac and COPASI are successful at extracting only a subset of the conservation laws that our algorithm can. This is illustrated with examples for some large networks which show the advantages of our method.     

New computational approaches for analysis of cis-regulatory networks

The investigation and modeling of gene regulatory networks requires computational tools specific to the task. We present several locally developed software tools that have been used in support of our ongoing research into the embryogenesis of the sea urchin. These tools are especially well suited to iterative refinement of models through experimental and computational investigation. They include: BioArray, a macroarray spot processing program; SUGAR, a system to display and correlate large-BAC sequence analyses; SeqComp and FamilyRelations, programs for comparative sequence analysis; and NetBuilder, an environment for creating and analyzing models of gene networks. We also present an overview of the process used to build our model of the Strongylocentrotus purpuratus endomesoderm gene network. Several of the tools discussed in this paper are still in active development and some are available as open source.     

A computational tool for quantitative analysis of vascular networks

Angiogenesis is the generation of mature vascular networks from pre-existing vessels. Angiogenesis is crucial during the organism' development, for wound healing and for the female reproductive cycle. Several murine experimental systems are well suited for studying developmental and pathological angiogenesis. They include the embryonic hindbrain, the post-natal retina and allantois explants. In these systems vascular networks are visualised by appropriate staining procedures followed by microscopical analysis. Nevertheless, quantitative assessment of angiogenesis is hampered by the lack of readily available, standardized metrics and software analysis tools. Non-automated protocols are being used widely and they are, in general, time--and labour intensive, prone to human error and do not permit computation of complex spatial metrics. We have developed a light-weight, user friendly software, AngioTool, which allows for quick, hands-off and reproducible quantification of vascular networks in microscopic images. AngioTool computes several morphological and spatial parameters including the area covered by a vascular network, the number of vessels, vessel length, vascular density and lacunarity. In addition, AngioTool calculates the so-called "branching index" (branch points/unit area), providing a measurement of the sprouting activity of a specimen of interest. We have validated AngioTool using images of embryonic murine hindbrains, post-natal retinas and allantois explants. AngioTool is open source and can be downloaded free of charge.     

CONTROL: software for the analysis of the control of metabolic networks

The program CONTROL is based on metabolic control theory anduses the method developed by Reder (1988). In this theory, twosets of parameters are defined in the vicinity of a steady-state:the elasticity coefficients which describe the local behaviourof the isolated enzymes, and the control coefficients whichexpress the response of the whole metabolic network to perturbationsat a given step. The theory shows that relationships exist betweenthe control coefficients (summation relationships or structuralrelationships) and also between the two types of coefficients(control and elasticity coefficients: connectivity relationships).The program CONTROL is divided into two parts (sub-menus). Thefirst one calculates all the control coefficients (flux andconcentrations) of a metabolic network from the elasticity coefficients.Using the second menu, the symbolic relationships are obtainedbetween the control coefficients (summation relationships) andbetween the control coefficients and the elasticity coefficients(connectivity relationships). These two sub-menus can be appliedindependently to any metabolic network (to date limited to 19steps and 19 metabolites).     

A new computational method to split large biochemical networks into coherent subnets

Background  

Compared to more general networks, biochemical networks have some special features: while generally sparse, there are a small number of highly connected metabolite nodes; and metabolite nodes can also be divided into two classes: internal nodes with associated mass balance constraints and external ones without. Based on these features, reclassifying selected internal nodes (separators) to external ones can be used to divide a large complex metabolic network into simpler subnetworks. Selection of separators based on node connectivity is commonly used but affords little detailed control and tends to produce excessive fragmentation.     

Condor-COPASI: high-throughput computing for biochemical networks

ABSTRACT: BACKGROUND: Mathematical modelling has become a standard technique to improve our understanding of complex biological systems. As models become larger and more complex, simulations and analyses require increasing amounts of computational power. Clusters of computers in a high-throughput computing environment can help to provide the resources required for computationally expensive model analysis. However, exploiting such a system can be difficult for users without the necessary expertise. RESULTS: We present Condor-COPASI, a server-based software tool that integrates COPASI, a biological pathway simulation tool, with Condor, a high-throughput computing environment. Condor-COPASI provides a web-based interface, which makes it extremely easy for a user to run a number of model simulation and analysis tasks in parallel. Tasks are transparently split into smaller parts, and and submitted for execution on a Condor pool. Result output is presented to the user in a number of formats, including tables and interactive graphical displays. CONCLUSIONS: Condor-COPASI can effectively use a Condor high-throughput computing environment to provide significant gains in performance for a number of model simulation and analysis tasks. Condor-COPASI is free, open source software, released under the Artistic License 2.0, and is suitable for use by any institution with access to a Condor pool. Source code is freely available for download at http://code.google.com/p/condor-copasi/, along with full instructions on deployment and usage.     

RMBNToolbox: random models for biochemical networks

Background  

There is an increasing interest to model biochemical and cell biological networks, as well as to the computational analysis of these models. The development of analysis methodologies and related software is rapid in the field. However, the number of available models is still relatively small and the model sizes remain limited. The lack of kinetic information is usually the limiting factor for the construction of detailed simulation models.     

Timescale analysis of rule-based biochemical reaction networks

The flow of information within a cell is governed by a series of protein–protein interactions that can be described as a reaction network. Mathematical models of biochemical reaction networks can be constructed by repetitively applying specific rules that define how reactants interact and what new species are formed on reaction. To aid in understanding the underlying biochemistry, timescale analysis is one method developed to prune the size of the reaction network. In this work, we extend the methods associated with timescale analysis to reaction rules instead of the species contained within the network. To illustrate this approach, we applied timescale analysis to a simple receptor–ligand binding model and a rule‐based model of interleukin‐12 (IL‐12) signaling in naïve CD4+ T cells. The IL‐12 signaling pathway includes multiple protein–protein interactions that collectively transmit information; however, the level of mechanistic detail sufficient to capture the observed dynamics has not been justified based on the available data. The analysis correctly predicted that reactions associated with Janus Kinase 2 and Tyrosine Kinase 2 binding to their corresponding receptor exist at a pseudo‐equilibrium. By contrast, reactions associated with ligand binding and receptor turnover regulate cellular response to IL‐12. An empirical Bayesian approach was used to estimate the uncertainty in the timescales. This approach complements existing rank‐ and flux‐based methods that can be used to interrogate complex reaction networks. Ultimately, timescale analysis of rule‐based models is a computational tool that can be used to reveal the biochemical steps that regulate signaling dynamics. © 2011 American Institute of Chemical Engineers Biotechnol. Prog., 2012     

Cellware--a multi-algorithmic software for computational systems biology

The intracellular environment of a cell hosts a wide variety of enzymatic reactions, diffusion events, molecular binding, polymerization and metabolic channeling. To transform these biological events into a computational framework, distinct modeling strategies are required. While currently no tool is capable of capturing all these events, progress is being made to create an integrated environment for the modeling community. To address this niche requirement, Cellware has been developed to offer a multi-algorithmic environment for modeling and simulating both deterministic and stochastic events in the cell. AVAILABILITY: The software is available for free and can be downloaded from http://www.bii.a-star.edu.sg/sbg/cellware     

Quasi-multiparameter sensitivity measure for robustness analysis of complex biochemical networks

Robustness is the ability to resume reliable operation in the face of different types of perturbations. Analysis of how network structure achieves robustness enables one to understand and design cellular systems. It is typically true that all parameters simultaneously differ from their nominal values in vivo, but there have been few intelligible measures to estimate the robustness of a system's function to the uncertainty of all parameters.We propose a numerical and fast measure of a robust property to the uncertainty of all kinetic parameters, named quasi-multiparameter sensitivity (QMPS), which is defined as the sum of the squared magnitudes of single-parameter sensitivities. Despite its plain idea, it has hardly been employed in analysis of biological models. While QMPS is theoretically derived as a linear model, QMPS can be consistent with the expected variance simulated by the widely used Monte Carlo method in nonlinear biological models, when relatively small perturbations are given. To demonstrate the feasibility of QMPS, it is employed for numerical comparison to analyze the mechanism of how specific regulations generate robustness in typical biological models.QMPS characterizes the robustness much faster than the Monte Carlo method, thereby enabling the extensive search of a large parameter space to perform the numerical comparison between alternative or competing models. It provides a theoretical or quantitative insight to an understanding of how specific network structures are related to robustness. In circadian oscillators, a negative feedback loop with multiple phosphorylations is demonstrated to play a critical role in generating robust cycles to the uncertainty of multiple parameters.     

Bioconductor: open software development for computational biology and bioinformatics

The Bioconductor project is an initiative for the collaborative creation of extensible software for computational biology and bioinformatics. The goals of the project include: fostering collaborative development and widespread use of innovative software, reducing barriers to entry into interdisciplinary scientific research, and promoting the achievement of remote reproducibility of research results. We describe details of our aims and methods, identify current challenges, compare Bioconductor to other open bioinformatics projects, and provide working examples.     

Top-down analysis of temporal hierarchy in biochemical reaction networks

The study of dynamic functions of large-scale biological networks has intensified in recent years. A critical component in developing an understanding of such dynamics involves the study of their hierarchical organization. We investigate the temporal hierarchy in biochemical reaction networks focusing on: (1) the elucidation of the existence of "pools" (i.e., aggregate variables) formed from component concentrations and (2) the determination of their composition and interactions over different time scales. To date the identification of such pools without prior knowledge of their composition has been a challenge. A new approach is developed for the algorithmic identification of pool formation using correlations between elements of the modal matrix that correspond to a pair of concentrations and how such correlations form over the hierarchy of time scales. The analysis elucidates a temporal hierarchy of events that range from chemical equilibration events to the formation of physiologically meaningful pools, culminating in a network-scale (dynamic) structure-(physiological) function relationship. This method is validated on a model of human red blood cell metabolism and further applied to kinetic models of yeast glycolysis and human folate metabolism, enabling the simplification of these models. The understanding of temporal hierarchy and the formation of dynamic aggregates on different time scales is foundational to the study of network dynamics and has relevance in multiple areas ranging from bacterial strain design and metabolic engineering to the understanding of disease processes in humans.     

Graph analysis of functional brain networks: practical issues in translational neuroscience

The brain can be regarded as a network: a connected system where nodes, or units, represent different specialized regions and links, or connections, represent communication pathways. From a functional perspective, communication is coded by temporal dependence between the activities of different brain areas. In the last decade, the abstract representation of the brain as a graph has allowed to visualize functional brain networks and describe their non-trivial topological properties in a compact and objective way. Nowadays, the use of graph analysis in translational neuroscience has become essential to quantify brain dysfunctions in terms of aberrant reconfiguration of functional brain networks. Despite its evident impact, graph analysis of functional brain networks is not a simple toolbox that can be blindly applied to brain signals. On the one hand, it requires the know-how of all the methodological steps of the pipeline that manipulate the input brain signals and extract the functional network properties. On the other hand, knowledge of the neural phenomenon under study is required to perform physiologically relevant analysis. The aim of this review is to provide practical indications to make sense of brain network analysis and contrast counterproductive attitudes.     

The Virtual Cell: a software environment for computational cell biology

The newly emerging field of computational cell biology requires software tools that address the needs of a broad community of scientists. Cell biological processes are controlled by an interacting set of biochemical and electrophysiological events that are distributed within complex cellular structures. Computational modeling is familiar to researchers in fields such as molecular structure, neurobiology and metabolic pathway engineering, and is rapidly emerging in the area of gene expression. Although some of these established modeling approaches can be adapted to address problems of interest to cell biologists, relatively few software development efforts have been directed at the field as a whole. The Virtual Cell is a computational environment designed for cell biologists as well as for mathematical biologists and bioengineers. It serves to aid the construction of cell biological models and the generation of simulations from them. The system enables the formulation of both compartmental and spatial models, the latter with either idealized or experimentally derived geometries of one, two or three dimensions.     

Mass spectrometry software for biochemical analysis in electrospray and fast atom bombardment modes.

New mass spectrometry techniques, such as electrospray ionization (ESI), allow the study of large biomolecules and peptide mixtures. The data produced are complex and interpretation can be a long and tedious process. A new suite of data-processing software is described which allows many of these operations to be carried out in a rapid, automated way. Software is described for the deconvolution of the spectra of multiply charged ions, for both pure compounds and mixtures. The rapid peptide mapping of protein digests from h.p.l.c.-m.s. data and peptide sequence confirmation from multiple-stage (MS)-m.s. data using tandem quadrupole m.s. are also described. In addition preliminary results are presented on first principle sequencing of unknown peptides from MS-m.s. experiments.     

Identifying fragilities in biochemical networks: robust performance analysis of Fas signaling-induced apoptosis

Proper control of apoptotic signaling is critical to immune response and development in multicellular organisms. Two tools from control engineering are applied to a mathematical model of Fas ligand signaling-induced apoptosis. Structured singular value analysis determines the volume in parameter space within which the system parameters may exist and still maintain efficacious signaling, but is limited to linear behaviors. Sensitivity analysis can be applied to nonlinear systems but is difficult to relate to performance criteria. Thus, structured singular value analysis is used to quantify performance during apoptosis rejection, ensuring that the system remains sensitive but not overly so to apoptotic stimuli. Sensitivity analysis is applied when the system has switched to the death-inducing, apoptotic steady state to determine parameters significant to maintaining the bistability. The analyses reveal that the magnitude of the death signal is fragile to perturbations in degradation parameters (failures in the ubiquitin/proteasome mechanism) while the timing of signal expression can be tuned by manipulating local parameters. Simultaneous parameter uncertainty highlights apoptotic fragility to disturbances in the ubiquitin/proteasome system. Sensitivity analysis reveals that the robust signaling characteristics of the apoptotic network is due to network architecture, and the apoptotic signaling threshold is best manipulated by interactions upstream of the apoptosome.     

Extended CADLIVE: a novel graphical notation for design of biochemical network maps and computational pathway analysis

Biochemical network maps are helpful for understanding the mechanism of how a collection of biochemical reactions generate particular functions within a cell. We developed a new and computationally feasible notation that enables drawing a wide resolution map from the domain-level reactions to phenomenological events and implemented it as the extended GUI network constructor of CADLIVE (Computer-Aided Design of LIVing systEms). The new notation presents ‘Domain expansion’ for proteins and RNAs, ‘Virtual reaction and nodes’ that are responsible for illustrating domain-based interaction and ‘InnerLink’ that links real complex nodes to virtual nodes to illustrate the exact components of the real complex. A modular box is also presented that packs related reactions as a module or a subnetwork, which gives CADLIVE a capability to draw biochemical maps in a hierarchical modular architecture. Furthermore, we developed a pathway search module for virtual knockout mutants as a built-in application of CADLIVE. This module analyzes gene function in the same way as molecular genetics, which simulates a change in mutant phenotypes or confirms the validity of the network map. The extended CADLIVE with the newly proposed notation is demonstrated to be feasible for computational simulation and analysis.