Search

Motivation: Efficient resource allocation can contribute to an organism's fitness and can improve evolutionary success. Resource Balance Analysis (RBA) is a computational framework that models an organism's growth-optimal proteome configurations in various environments. RBA software enables the construction of RBA models on genome scale and the calculation of medium-specific, growth-optimal cell states including metabolic fluxes and the abundance of macromolecular machines. However, existing software lacks a simple programming interface for non-expert users, easy to use and interoperable with other software., Results: The python package RBAtools provides convenient access to RBA models. As a flexible programming interface, it enables the implementation of custom workflows and the modification of existing genome-scale RBA models. Its high-level functions comprise simulation, model fitting, parameter screens, sensitivity analysis, variability analysis and the construction of Pareto fronts. Models and data are represented as structured tables and can be exported to common data formats for fluxomics and proteomics visualization., Availability and Implementation: RBAtools documentation, installation instructions and tutorials are available at https://sysbioinra.github.io/rbatools/. General information about RBA and related software can be found at rba.inrae.fr., (© The Author(s) 2023. Published by Oxford University Press.)

Computational models have great potential to accelerate bioscience, bioengineering, and medicine. However, it remains challenging to reproduce and reuse simulations, in part, because the numerous formats and methods for simulating various subsystems and scales remain siloed by different software tools. For example, each tool must be executed through a distinct interface. To help investigators find and use simulation tools, we developed BioSimulators (https://biosimulators.org), a central registry of the capabilities of simulation tools and consistent Python, command-line and containerized interfaces to each version of each tool. The foundation of BioSimulators is standards, such as CellML, SBML, SED-ML and the COMBINE archive format, and validation tools for simulation projects and simulation tools that ensure these standards are used consistently. To help modelers find tools for particular projects, we have also used the registry to develop recommendation services. We anticipate that BioSimulators will help modelers exchange, reproduce, and combine simulations., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Summary: Measured kinetic constants are key input data for metabolic models, but they are often uncertain, inconsistent and incomplete. Parameter balancing translates such data into complete and consistent parameter sets while accounting for predefined ranges and physical constraints. Based on Bayesian regression, it determines a most plausible parameter set as well as uncertainty ranges for all model parameters. Our tools for parameter balancing support standard model and data formats and enable an easy customization of prior distributions and constraints for biochemical constants. Modellers can balance kinetic constants, thermodynamic data and metabolomic data to obtain thermodynamically consistent metabolic states that comply with user-defined flux directions., Availability and Implementation: An online tool for parameter balancing, a stand-alone Python command line tool, a Python package and a Matlab toolbox (which uses the CPLEX solver) are freely available at www.parameterbalancing.net., (© The Author(s) 2019. Published by Oxford University Press.)

Systems biology models are rapidly increasing in complexity, size and numbers. When building large models, researchers rely on software tools for the retrieval, comparison, combination and merging of models, as well as for version control. These tools need to be able to quantify the differences and similarities between computational models. However, depending on the specific application, the notion of 'similarity' may greatly vary. A general notion of model similarity, applicable to various types of models, is still missing. Here we survey existing methods for the comparison of models, introduce quantitative measures for model similarity, and discuss potential applications of combined similarity measures. To frame model comparison as a general problem, we describe a theoretical approach to defining and computing similarities based on a combination of different model aspects. The six aspects that we define as potentially relevant for similarity are underlying encoding, references to biological entities, quantitative behaviour, qualitative behaviour, mathematical equations and parameters and network structure. We argue that future similarity measures will benefit from combining these model aspects in flexible, problem-specific ways to mimic users' intuition about model similarity, and to support complex model searches in databases., (© The Author 2016. Published by Oxford University Press.)

Unlabelled: : SBtab is a table-based data format for Systems Biology, designed to support automated data integration and model building. It uses the structure of spreadsheets and defines conventions for table structure, controlled vocabularies and semantic annotations. The format comes with predefined table types for experimental data and SBML-compliant model structures and can easily be customized to cover new types of data., Availability and Implementation: SBtab documents can be created and edited with any text editor or spreadsheet tool. The website www.sbtab.net provides online tools for syntax validation and conversion to SBML and HTML, as well as software for using SBtab in MS Excel, MATLAB and R. The stand-alone Python code contains functions for file parsing, validation, conversion to SBML and HTML and an interface to SQLite databases, to be integrated into Systems Biology workflows. A detailed specification of SBtab, including examples and descriptions of table types and available tools, can be found at www.sbtab.net, Contact: : wolfram.liebermeister@gmail.com., (© The Author 2016. Published by Oxford University Press.)

Protein levels are a dominant factor shaping natural and synthetic biological systems. Although proper functioning of metabolic pathways relies on precise control of enzyme levels, the experimental ability to balance the levels of many genes in parallel is a major outstanding challenge. Here, we introduce a rapid and modular method to span the expression space of several proteins in parallel. By combinatorially pairing genes with a compact set of ribosome-binding sites, we modulate protein abundance by several orders of magnitude. We demonstrate our strategy by using a synthetic operon containing fluorescent proteins to span a 3D color space. Using the same approach, we modulate a recombinant carotenoid biosynthesis pathway in Escherichia coli to reveal a diversity of phenotypes, each characterized by a distinct carotenoid accumulation profile. In a single combinatorial assembly, we achieve a yield of the industrially valuable compound astaxanthin 4-fold higher than previously reported. The methodology presented here provides an efficient tool for exploring a high-dimensional expression space to locate desirable phenotypes.

Motivation: Standard rate laws are a key requisite for systematically turning metabolic networks into kinetic models. They should provide simple, general and biochemically plausible formulae for reaction velocities and reaction elasticities. At the same time, they need to respect thermodynamic relations between the kinetic constants and the metabolic fluxes and concentrations., Results: We present a family of reversible rate laws for reactions with arbitrary stoichiometries and various types of regulation, including mass-action, Michaelis-Menten and uni-uni reversible Hill kinetics as special cases. With a thermodynamically safe parameterization of these rate laws, parameter sets obtained by model fitting, sampling or optimization are guaranteed to lead to consistent chemical equilibrium states. A reformulation using saturation values yields simple formulae for rates and elasticities, which can be easily adjusted to the given stationary flux distributions. Furthermore, this formulation highlights the role of chemical potential differences as thermodynamic driving forces. We compare the modular rate laws to the thermodynamic-kinetic modelling formalism and discuss a simplified rate law in which the reaction rate directly depends on the reaction affinity. For automatic handling of modular rate laws, we propose a standard syntax and semantic annotations for the Systems Biology Markup Language., Availability: An online tool for inserting the rate laws into SBML models is freely available at www.semanticsbml.org., Supplementary Information: Supplementary data are available at Bioinformatics online.

Motivation: The availability of modern sequencing techniques has led to a rapid increase in the amount of reconstructed metabolic networks. Using these models as a platform for the analysis of high throughput transcriptomic, proteomic and metabolomic data can provide valuable insight into conditional changes in the metabolic activity of an organism. While transcriptomics and proteomics provide important insights into the hierarchical regulation of metabolic flux, metabolomics shed light on the actual enzyme activity through metabolic regulation and mass action effects. Here we introduce a new method, termed integrative omics-metabolic analysis (IOMA) that quantitatively integrates proteomic and metabolomic data with genome-scale metabolic models, to more accurately predict metabolic flux distributions. The method is formulated as a quadratic programming (QP) problem that seeks a steady-state flux distribution in which flux through reactions with measured proteomic and metabolomic data, is as consistent as possible with kinetically derived flux estimations., Results: IOMA is shown to successfully predict the metabolic state of human erythrocytes (compared to kinetic model simulations), showing a significant advantage over the commonly used methods flux balance analysis and minimization of metabolic adjustment. Thereafter, IOMA is shown to correctly predict metabolic fluxes in Escherichia coli under different gene knockouts for which both metabolomic and proteomic data is available, achieving higher prediction accuracy over the extant methods. Considering the lack of high-throughput flux measurements, while high-throughput metabolomic and proteomic data are becoming readily available, we expect IOMA to significantly contribute to future research of cellular metabolism.

Summary: Systems Biology Markup Language (SBML) is the leading exchange format for mathematical models in Systems Biology. Semantic annotations link model elements with external knowledge via unique database identifiers and ontology terms, enabling software to check and process models by their biochemical meaning. Such information is essential for model merging, one of the key steps towards the construction of large kinetic models. SemanticSBML is a tool that helps users to check and edit MIRIAM annotations and SBO terms in SBML models. Using a large collection of biochemical names and database identifiers, it supports modellers in finding the right annotations and in merging existing models. Initially, an element matching is derived from the MIRIAM annotations and conflicting element attributes are categorized and highlighted. Conflicts can then be resolved automatically or manually, allowing the user to control the merging process in detail., Availability: SemanticSBML comes as a free software written in Python and released under the GPL 3. A Debian package, a source package for other Linux distributions, a Windows installer and an online version of semanticSBML with limited functionality are available at http://www.semanticsbml.org. A preinstalled version can be found on the Linux live DVD SB.OS, available at http://www.sbos.eu., Supplementary Information: Supplementary data are available at Bioinformatics online.

Motivation: The expression of genes is controlled by specific combinations of cellular variables. We applied Independent Component Analysis (ICA) to gene expression data, deriving a linear model based on hidden variables, which we term 'expression modes'. The expression of each gene is a linear function of the expression modes, where, according to the ICA model, the linear influences of different modes show a minimal statistical dependence, and their distributions deviate sharply from the normal distribution., Results: Studying cell cycle-related gene expression in yeast, we found that the dominant expression modes could be related to distinct biological functions, such as phases of the cell cycle or the mating response. Analysis of human lymphocytes revealed modes that were related to characteristic differences between cell types. With both data sets, the linear influences of the dominant modes showed distributions with large tails, indicating the existence of specifically up- and downregulated target genes. The expression modes and their influences can be used to visualize the samples and genes in low-dimensional spaces. A projection to expression modes helps to highlight particular biological functions, to reduce noise, and to compress the data in a biologically sensible way.