Timo Lubitz

Annotation and merging of SBML models with semanticSBML

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.

Systematic Construction of Kinetic Models from Genome-Scale Metabolic Networks

The quantitative effects of environmental and genetic perturbations on metabolism can be studied in silico using kinetic models. We present a strategy for large-scale model construction based on a logical layering of data such as reaction fluxes, metabolite concentrations, and kinetic constants. The resulting models contain realistic standard rate laws and plausible parameters, adhere to the laws of thermodynamics, and reproduce a predefined steady state. These features have not been simultaneously achieved by previous workflows. We demonstrate the advantages and limitations of the workflow by translating the yeast consensus metabolic network into a kinetic model. Despite crudely selected data, the model shows realistic control behaviour, a stable dynamic, and realistic response to perturbations in extracellular glucose concentrations. The paper concludes by outlining how new data can continuously be fed into the workflow and how iterative model building can assist in directing experiments.

New types of experimental data shape the use of enzyme kinetics for dynamic network modeling

Since the publication of Leonor Michaelis and Maude Mentens paper on the reaction kinetics of the enzyme invertase in 1913, molecular biology has evolved tremendously. New measurement techniques allow in vivo characterization of the whole genome, proteome or transcriptome of cells, whereas the classical enzyme essay only allows determination of the two Michaelis–Menten parameters V and Km. Nevertheless, Michaelis–Menten kinetics are still commonly used, not only in the in vitro context of enzyme characterization but also as a rate law for enzymatic reactions in larger biochemical reaction networks. In this review, we give an overview of the historical development of kinetic rate laws originating from Michaelis–Menten kinetics over the past 100 years. Furthermore, we briefly summarize the experimental techniques used for the characterization of enzymes, and discuss web resources that systematically store kinetic parameters and related information. Finally, describe the novel opportunities that arise from using these data in dynamic mathematical modeling. In this framework, traditional in vitro approaches may be combined with modern genome‐scale measurements to foster thorough understanding of the underlying complex mechanisms.

Parameter Balancing in Kinetic Models of Cell Metabolism

Kinetic modeling of metabolic pathways has become a major field of systems biology. It combines structural information about metabolic pathways with quantitative enzymatic rate laws. Some of the kinetic constants needed for a model could be collected from ever-growing literature and public web resources, but they are often incomplete, incompatible, or simply not available. We address this lack of information by parameter balancing, a method to complete given sets of kinetic constants. Based on Bayesian parameter estimation, it exploits the thermodynamic dependencies among different biochemical quantities to guess realistic model parameters from available kinetic data. Our algorithm accounts for varying measurement conditions in the input data (pH value and temperature). It can process kinetic constants and state-dependent quantities such as metabolite concentrations or chemical potentials, and uses prior distributions and data augmentation to keep the estimated quantities within plausible ranges. An online service and free software for parameter balancing with models provided in SBML format (Systems Biology Markup Language) is accessible at www.semanticsbml.org. We demonstrate its practical use with a small model of the phosphofructokinase reaction and discuss its possible applications and limitations. In the future, parameter balancing could become an important routine step in the kinetic modeling of large metabolic networks.

Glucose de-repression by yeast AMP-activated protein kinase SNF1 is controlled via at least two independent steps

The AMP‐activated protein kinase, AMPK, controls energy homeostasis in eukaryotic cells but little is known about the mechanisms governing the dynamics of its activation/deactivation. The yeast AMPK, SNF1, is activated in response to glucose depletion and mediates glucose de‐repression by inactivating the transcriptional repressor Mig1. Here we show that overexpression of the Snf1‐activating kinase Sak1 results, in the presence of glucose, in constitutive Snf1 activation without alleviating glucose repression. Co‐overexpression of the regulatory subunit Reg1 of the Glc‐Reg1 phosphatase complex partly restores glucose regulation of Snf1. We generated a set of 24 kinetic mathematical models based on dynamic data of Snf1 pathway activation and deactivation. The models that reproduced our experimental observations best featured (a) glucose regulation of both Snf1 phosphorylation and dephosphorylation, (b) determination of the Mig1 phosphorylation status in the absence of glucose by Snf1 activity only and (c) a regulatory step directing active Snf1 to Mig1 under glucose limitation. Hence it appears that glucose de‐repression via Snf1‐Mig1 is regulated by glucose via at least two independent steps: the control of activation of the Snf1 kinase and directing active Snf1 to inactivating its target Mig1.

SBtab: a flexible table format for data exchange in systems biology

Summary: 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

Network reconstruction and validation of the Snf1/AMPK pathway in baker’s yeast based on a comprehensive literature review

Background/Objectives:The SNF1/AMPK protein kinase has a central role in energy homeostasis in eukaryotic cells. It is activated by energy depletion and stimulates processes leading to the production of ATP while it downregulates ATP-consuming processes. The yeast SNF1 complex is best known for its role in glucose derepression.Methods:We performed a network reconstruction of the Snf1 pathway based on a comprehensive literature review. The network was formalised in the rxncon language, and we used the rxncon toolbox for model validation and gap filling.Results:We present a machine-readable network definition that summarises the mechanistic knowledge of the Snf1 pathway. Furthermore, we used the known input/output relationships in the network to identify and fill gaps in the information transfer through the pathway, to produce a functional network model. Finally, we convert the functional network model into a rule-based model as a proof-of-principle.Conclusions:The workflow presented here enables large scale reconstruction, validation and gap filling of signal transduction networks. It is analogous to but distinct from that established for metabolic networks. We demonstrate the workflow capabilities, and the direct link between the reconstruction and dynamic modelling, with the Snf1 network. This network is a distillation of the knowledge from all previous publications on the Snf1/AMPK pathway. The network is a knowledge resource for modellers and experimentalists alike, and a template for similar efforts in higher eukaryotes. Finally, we envisage the workflow as an instrumental tool for reconstruction of large signalling networks across Eukaryota.