BioModels Database Research Articles

Computing Weakly Reversible Deficiency Zero Network Translations Using Elementary Flux Modes.

We present a computational method for performing structural translation, which has been studied recently in the context of analyzing the steady states and dynamical behavior of mass-action systems derived from biochemical reaction networks. Our procedure involves solving a binary linear programming problem where the decision variables correspond to interactions between the reactions of the original network. We call the resulting network a reaction-to-reaction graph and formalize how such a construction relates to the original reaction network and the structural translation. We demonstrate the efficacy and efficiency of the algorithm by running it on 508 networks from the European Bioinformatics Institutes' BioModels database. We also summarize how this work can be incorporated into recently proposed algorithms for establishing mono- and multistationarity in biochemical reaction systems.

Functional expression and impact of testis-specific gene antigen 10 in breast cancer: a combined in vitro and in silico analysis.

Introduction: Testis-specific gene antigen 10 (TSGA10) is a less-known gene, which is involved in the vague biological paths of different cancers. Here, we investigated the TSGA10 expression using different concentrations of glucose under hypoxia and also its interaction with the hypoxia-inducible factor 1 (HIF-1). Methods: The breast cancer MDA-MB-231 and MCF-7 cells were cultured with different concentrations of glucose (5.5, 11.0 and 25.0 mM) under normoxia/hypoxia for 24, 48, and 72 hours and examined for the HIF-1α expression and cell migration by Western blotting and scratch assays. The qPCR was employed to analyze the expression of TSGA10. Three-dimensional (3D) structure and the energy minimization of the interacting domain of TSGA10 were performed by MODELLER v9.17 and Swiss-PDB viewer v4.1.0/UCSF Chimera v1.11. The UCSF Chimera v1.13.1 and Hex 6.0 were used for the molecular docking simulation. The Cytoscape v3.7.1 and STRING v11.0 were used for protein-protein interaction (PPI) network analysis. The HIF-1a related hypoxia pathways were obtained from BioModels database and reconstructed in CellDesigner v4.4.2. Results: The increased expression of TSGA10 was found to be significantly associated with the reduced metastasis in the MDA-MB-231 cells, while an inverse relationship was seen between the TSGA10 mRNA level and cellular migration but not in the MCF-7 cells. The C-terminal domain of TSGA10 interacted with HIF-1α with high affinity, resulting in PPI network with 10 key nodes (HIF-1α, VEGFA, HSP90AA1, AKT1, ARNT, TP53, TSGA10, VHL, JUN, and EGFR). Conclusions: Collectively, TSGA10 functional expression alters under the hyper-/hypo-glycemia and hypoxia, which indicates its importance as a candidate bio-target for the cancer therapy.

Evolution of computational models in BioModels Database and the Physiome Model Repository

BackgroundA useful model is one that is being (re)used. The development of a successful model does not finish with its publication. During reuse, models are being modified, i.e. expanded, corrected, and refined. Even small changes in the encoding of a model can, however, significantly affect its interpretation. Our motivation for the present study is to identify changes in models and make them transparent and traceable.MethodsWe analysed 13734 models from BioModels Database and the Physiome Model Repository. For each model, we studied the frequencies and types of updates between its first and latest release. To demonstrate the impact of changes, we explored the history of a Repressilator model in BioModels Database.ResultsWe observed continuous updates in the majority of models. Surprisingly, even the early models are still being modified. We furthermore detected that many updates target annotations, which improves the information one can gain from models. To support the analysis of changes in model repositories we developed MoSt, an online tool for visualisations of changes in models. The scripts used to generate the data and figures for this study are available from GitHub github.com/binfalse/BiVeS-StatsGenerator and as a Docker image at hub.docker.com/r/binfalse/bives-statsgenerator. The website most.bio.informatik.uni-rostock.de provides interactive access to model versions and their evolutionary statistics.ConclusionThe reuse of models is still impeded by a lack of trust and documentation. A detailed and transparent documentation of all aspects of the model, including its provenance, will improve this situation. Knowledge about a model’s provenance can avoid the repetition of mistakes that others already faced. More insights are gained into how the system evolves from initial findings to a profound understanding. We argue that it is the responsibility of the maintainers of model repositories to offer transparent model provenance to their users.

MAGPIE: Simplifying access and execution of computational models in the life sciences.

Over the past decades, quantitative methods linking theory and observation became increasingly important in many areas of life science. Subsequently, a large number of mathematical and computational models has been developed. The BioModels database alone lists more than 140,000 Systems Biology Markup Language (SBML) models. However, while the exchange within specific model classes has been supported by standardisation and database efforts, the generic application and especially the re-use of models is still limited by practical issues such as easy and straight forward model execution. MAGPIE, a Modeling and Analysis Generic Platform with Integrated Evaluation, closes this gap by providing a software platform for both, publishing and executing computational models without restrictions on the programming language, thereby combining a maximum on flexibility for programmers with easy handling for non-technical users. MAGPIE goes beyond classical SBML platforms by including all models, independent of the underlying programming language, ranging from simple script models to complex data integration and computations. We demonstrate the versatility of MAGPIE using four prototypic example cases. We also outline the potential of MAGPIE to improve transparency and reproducibility of computational models in life sciences. A demo server is available at magpie.imb.medizin.tu-dresden.de.

A computational approach to extinction events in chemical reaction networks with discrete state spaces

Recent work of Johnston et al. has produced sufficient conditions on the structure of a chemical reaction network which guarantee that the corresponding discrete state space system exhibits an extinction event. The conditions consist of a series of systems of equalities and inequalities on the edges of a modified reaction network called a domination-expanded reaction network. In this paper, we present a computational implementation of these conditions written in Python and apply the program on examples drawn from the biochemical literature. We also run the program on 458 models from the European Bioinformatics Institute’s BioModels Database and report our results.

An improved hybrid of particle swarm optimization and the gravitational search algorithm to produce a kinetic parameter estimation of aspartate biochemical pathways

Mathematical modelling is fundamental to understand the dynamic behavior and regulation of the biochemical metabolisms and pathways that are found in biological systems. Pathways are used to describe complex processes that involve many parameters. It is important to have an accurate and complete set of parameters that describe the characteristics of a given model. However, measuring these parameters is typically difficult and even impossible in some cases. Furthermore, the experimental data are often incomplete and also suffer from experimental noise. These shortcomings make it challenging to identify the best-fit parameters that can represent the actual biological processes involved in biological systems. Computational approaches are required to estimate these parameters. The estimation is converted into multimodal optimization problems that require a global optimization algorithm that can avoid local solutions. These local solutions can lead to a bad fit when calibrating with a model. Although the model itself can potentially match a set of experimental data, a high-performance estimation algorithm is required to improve the quality of the solutions.This paper describes an improved hybrid of particle swarm optimization and the gravitational search algorithm (IPSOGSA) to improve the efficiency of a global optimum (the best set of kinetic parameter values) search. The findings suggest that the proposed algorithm is capable of narrowing down the search space by exploiting the feasible solution areas. Hence, the proposed algorithm is able to achieve a near-optimal set of parameters at a fast convergence speed. The proposed algorithm was tested and evaluated based on two aspartate pathways that were obtained from the BioModels Database. The results show that the proposed algorithm outperformed other standard optimization algorithms in terms of accuracy and near-optimal kinetic parameter estimation. Nevertheless, the proposed algorithm is only expected to work well in small scale systems. In addition, the results of this study can be used to estimate kinetic parameter values in the stage of model selection for different experimental conditions.

Whole-genome metabolic model of Trichoderma reesei built by comparative reconstruction.

Background Trichoderma reesei is one of the main sources of biomass-hydrolyzing enzymes for the biotechnology industry. There is a need for improving its enzyme production efficiency. The use of metabolic modeling for the simulation and prediction of this organism’s metabolism is potentially a valuable tool for improving its capabilities. An accurate metabolic model is needed to perform metabolic modeling analysis.ResultsA whole-genome metabolic model of T. reesei has been reconstructed together with metabolic models of 55 related species using the metabolic model reconstruction algorithm CoReCo. The previously published CoReCo method has been improved to obtain better quality models. The main improvements are the creation of a unified database of reactions and compounds and the use of reaction directions as constraints in the gap-filling step of the algorithm. In addition, the biomass composition of T. reesei has been measured experimentally to build and include a specific biomass equation in the model.ConclusionsThe improvements presented in this work on the CoReCo pipeline for metabolic model reconstruction resulted in higher-quality metabolic models compared with previous versions. A metabolic model of T. reesei has been created and is publicly available in the BIOMODELS database. The model contains a biomass equation, reaction boundaries and uptake/export reactions which make it ready for simulation. To validate the model, we dem1onstrate that the model is able to predict biomass production accurately and no stoichiometrically infeasible yields are detected. The new T. reesei model is ready to be used for simulations of protein production processes.Electronic supplementary materialThe online version of this article (doi:10.1186/s13068-016-0665-0) contains supplementary material, which is available to authorized users.

A rule-based model of insulin signalling pathway.

BackgroundThe insulin signalling pathway (ISP) is an important biochemical pathway, which regulates some fundamental biological functions such as glucose and lipid metabolism, protein synthesis, cell proliferation, cell differentiation and apoptosis. In the last years, different mathematical models based on ordinary differential equations have been proposed in the literature to describe specific features of the ISP, thus providing a description of the behaviour of the system and its emerging properties. However, protein-protein interactions potentially generate a multiplicity of distinct chemical species, an issue referred to as “combinatorial complexity”, which results in defining a high number of state variables equal to the number of possible protein modifications. This often leads to complex, error prone and difficult to handle model definitions.ResultsIn this work, we present a comprehensive model of the ISP, which integrates three models previously available in the literature by using the rule-based modelling (RBM) approach. RBM allows for a simple description of a number of signalling pathway characteristics, such as the phosphorylation of signalling proteins at multiple sites with different effects, the simultaneous interaction of many molecules of the signalling pathways with several binding partners, and the information about subcellular localization where reactions take place. Thanks to its modularity, it also allows an easy integration of different pathways.After RBM specification, we simulated the dynamic behaviour of the ISP model and validated it using experimental data. We the examined the predicted profiles of all the active species and clustered them in four clusters according to their dynamic behaviour. Finally, we used parametric sensitivity analysis to show the role of negative feedback loops in controlling the robustness of the system.ConclusionsThe presented ISP model is a powerful tool for data simulation and can be used in combination with experimental approaches to guide the experimental design. The model is available at http://sysbiobig.dei.unipd.it/ was submitted to Biomodels Database (https://www.ebi.ac.uk/biomodels-main/# MODEL 1604100005).Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-016-0281-4) contains supplementary material, which is available to authorized users.

Challenges in horizontal model integration.

BackgroundSystems Biology has motivated dynamic models of important intracellular processes at the pathway level, for example, in signal transduction and cell cycle control. To answer important biomedical questions, however, one has to go beyond the study of isolated pathways towards the joint study of interacting signaling pathways or the joint study of signal transduction and cell cycle control. Thereby the reuse of established models is preferable, as it will generally reduce the modeling effort and increase the acceptance of the combined model in the field.ResultsObtaining a combined model can be challenging, especially if the submodels are large and/or come from different working groups (as is generally the case, when models stored in established repositories are used). To support this task, we describe a semi-automatic workflow based on established software tools. In particular, two frequent challenges are described: identification of the overlap and subsequent (re)parameterization of the integrated model.ConclusionsThe reparameterization step is crucial, if the goal is to obtain a model that can reproduce the data explained by the individual models. For demonstration purposes we apply our workflow to integrate two signaling pathways (EGF and NGF) from the BioModels Database.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-016-0266-3) contains supplementary material, which is available to authorized users.

Kinetic modeling and the rise of systems pharmacology

Kinetic modeling and the rise of systems pharmacology

An algorithm to detect and communicate the differences in computational models describing biological systems.

Motivation: Repositories support the reuse of models and ensure transparency about results in publications linked to those models. With thousands of models available in repositories, such as the BioModels database or the Physiome Model Repository, a framework to track the differences between models and their versions is essential to compare and combine models. Difference detection not only allows users to study the history of models but also helps in the detection of errors and inconsistencies. Existing repositories lack algorithms to track a model’s development over time.Results: Focusing on SBML and CellML, we present an algorithm to accurately detect and describe differences between coexisting versions of a model with respect to (i) the models’ encoding, (ii) the structure of biological networks and (iii) mathematical expressions. This algorithm is implemented in a comprehensive and open source library called BiVeS. BiVeS helps to identify and characterize changes in computational models and thereby contributes to the documentation of a model’s history. Our work facilitates the reuse and extension of existing models and supports collaborative modelling. Finally, it contributes to better reproducibility of modelling results and to the challenge of model provenance.Availability and implementation: The workflow described in this article is implemented in BiVeS. BiVeS is freely available as source code and binary from sems.uni-rostock.de. The web interface BudHat demonstrates the capabilities of BiVeS at budhat.sems.uni-rostock.de.Contact: martin.scharm@uni-rostock.deSupplementary information: Supplementary data are available at Bioinformatics online.

Finding the positive feedback loops underlying multi-stationarity.

BackgroundBistability is ubiquitous in biological systems. For example, bistability is found in many reaction networks that involve the control and execution of important biological functions, such as signaling processes. Positive feedback loops, composed of species and reactions, are necessary for bistability, and generally for multi-stationarity, to occur. These loops are therefore often used to illustrate and pinpoint the parts of a multi-stationary network that are relevant (‘responsible’) for the observed multi-stationarity. However positive feedback loops are generally abundant in reaction networks but not all of them are important for understanding the network’s dynamics.ResultsWe present an automated procedure to determine the relevant positive feedback loops of a multi-stationary reaction network. The procedure only reports the loops that are relevant for multi-stationarity (that is, when broken multi-stationarity disappears) and not all positive feedback loops of the network. We show that the relevant positive feedback loops must be understood in the context of the network (one loop might be relevant for one network, but cannot create multi-stationarity in another). Finally, we demonstrate the procedure by applying it to several examples of signaling processes, including a ubiquitination and an apoptosis network, and to models extracted from the Biomodels database. The procedure is implemented in Maple.ConclusionsWe have developed and implemented an automated procedure to find relevant positive feedback loops in reaction networks. The results of the procedure are useful for interpretation and summary of the network’s dynamics.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-015-0164-0) contains supplementary material, which is available to authorized users.

A computational approach to persistence, permanence, and endotacticity of biochemical reaction systems.

We introduce a mixed-integer linear programming (MILP) framework capable of determining whether a chemical reaction network possesses the property of being endotactic or strongly endotactic. The network property of being strongly endotactic is known to lead to persistence and permanence of chemical species under genetic kinetic assumptions, while the same result is conjectured but as yet unproved for general endotactic networks. The algorithms we present are the first capable of verifying endotacticity of chemical reaction networks for systems with greater than two constituent species. We implement the algorithms in the open-source online package CoNtRol and apply them to a large sample of networks from the European Bioinformatics Institute's BioModels Database. We use strong endotacticity to establish for the first time the permanence of a well-studied circadian clock mechanism.

Annotation-based feature extraction from sets of SBML models.

BackgroundModel repositories such as BioModels Database provide computational models of biological systems for the scientific community. These models contain rich semantic annotations that link model entities to concepts in well-established bio-ontologies such as Gene Ontology. Consequently, thematically similar models are likely to share similar annotations. Based on this assumption, we argue that semantic annotations are a suitable tool to characterize sets of models. These characteristics improve model classification, allow to identify additional features for model retrieval tasks, and enable the comparison of sets of models.ResultsIn this paper we discuss four methods for annotation-based feature extraction from model sets. We tested all methods on sets of models in SBML format which were composed from BioModels Database. To characterize each of these sets, we analyzed and extracted concepts from three frequently used ontologies, namely Gene Ontology, ChEBI and SBO. We find that three out of the methods are suitable to determine characteristic features for arbitrary sets of models: The selected features vary depending on the underlying model set, and they are also specific to the chosen model set. We show that the identified features map on concepts that are higher up in the hierarchy of the ontologies than the concepts used for model annotations. Our analysis also reveals that the information content of concepts in ontologies and their usage for model annotation do not correlate.ConclusionsAnnotation-based feature extraction enables the comparison of model sets, as opposed to existing methods for model-to-keyword comparison, or model-to-model comparison.Electronic supplementary materialThe online version of this article (doi:10.1186/s13326-015-0014-4) contains supplementary material, which is available to authorized users.

Combining computational models, semantic annotations and simulation experiments in a graph database.

Model repositories such as the BioModels Database, the CellML Model Repository or JWS Online are frequently accessed to retrieve computational models of biological systems. However, their storage concepts support only restricted types of queries and not all data inside the repositories can be retrieved. In this article we present a storage concept that meets this challenge. It grounds on a graph database, reflects the models’ structure, incorporates semantic annotations and simulation descriptions and ultimately connects different types of model-related data. The connections between heterogeneous model-related data and bio-ontologies enable efficient search via biological facts and grant access to new model features. The introduced concept notably improves the access of computational models and associated simulations in a model repository. This has positive effects on tasks such as model search, retrieval, ranking, matching and filtering. Furthermore, our work for the first time enables CellML- and Systems Biology Markup Language-encoded models to be effectively maintained in one database. We show how these models can be linked via annotations and queried.Database URL: https://sems.uni-rostock.de/projects/masymos/

Diagnostics for stochastic genome-scale modeling via model slicing and debugging.

Modeling of biological behavior has evolved from simple gene expression plots represented by mathematical equations to genome-scale systems biology networks. However, due to obstacles in complexity and scalability of creating genome-scale models, several biological modelers have turned to programming or scripting languages and away from modeling fundamentals. In doing so, they have traded the ability to have exchangeable, standardized model representation formats, while those that remain true to standardized model representation are faced with challenges in model complexity and analysis. We have developed a model diagnostic methodology inspired by program slicing and debugging and demonstrate the effectiveness of the methodology on a genome-scale metabolic network model published in the BioModels database. The computer-aided identification revealed specific points of interest such as reversibility of reactions, initialization of species amounts, and parameter estimation that improved a candidate cell's adenosine triphosphate production. We then compared the advantages of our methodology over other modeling techniques such as model checking and model reduction. A software application that implements the methodology is available at http://gel.ym.edu.tw/gcs/.

Effects of small particle numbers on long-term behaviour in discrete biochemical systems.

Motivation: The functioning of many biological processes depends on the appearance of only a small number of a single molecular species. Additionally, the observation of molecular crowding leads to the insight that even a high number of copies of species do not guarantee their interaction. How single particles contribute to stabilizing biological systems is not well understood yet. Hence, we aim at determining the influence of single molecules on the long-term behaviour of biological systems, i.e. whether they can reach a steady state.Results: We provide theoretical considerations and a tool to analyse Systems Biology Markup Language models for the possibility to stabilize because of the described effects. The theory is an extension of chemical organization theory, which we called discrete chemical organization theory. Furthermore we scanned the BioModels Database for the occurrence of discrete chemical organizations. To exemplify our method, we describe an application to the Template model of the mitotic spindle assembly checkpoint mechanism.Availability and implementation: http://www.biosys.uni-jena.de/Services.html.Contact: bashar.ibrahim@uni-jena.de or dittrich@minet.uni-jena.deSupplementary information: Supplementary data are available at Bioinformatics online.

BioModels linked dataset.

BackgroundBioModels Database is a reference repository of mathematical models used in biology. Models are stored as SBML files on a file system and metadata is provided in a relational database. Models can be retrieved through a web interface and programmatically via web services. In addition to those more traditional ways to access information, Linked Data using Semantic Web technologies (such as the Resource Description Framework, RDF), is becoming an increasingly popular means to describe and expose biological relevant data.ResultsWe present the BioModels Linked Dataset, which exposes the models’ content as a dereferencable interlinked dataset. BioModels Linked Dataset makes use of the wealth of annotations available within a large number of manually curated models to link and integrate data and models from other resources.ConclusionsThe BioModels Linked Dataset provides users with a dataset interoperable with other semantic web resources. It supports powerful search queries, some of which were not previously available to users and allow integration of data from multiple resources. This provides a distributed platform to find similar models for comparison, processing and enrichment.

Integrating Systems Biology Markup Language(SBML) with NEURON

The NEURON simulator software is widely used by the computational neuroscience community for neuronal electrophysiology modeling and we have recently extended it to support reaction-diffusion dynamics [1]. This class of dynamics is ideal for studying intracellular dynamics, such as calcium-induced calcium release. NEURON’s original reaction-diffusion support introduced a custom syntax; although this allowed taking advantage of NEURON’s full capabilities, it required existing cell biology models to be rewritten. To eliminate this barrier, we have introduced SBML support for NEURON. SBML [2] is an XML-based representation format used for specifying computational models of biological processes. This format is one of the most widely used, and is supported by over 200 software packages. It has become a standard used by systems biologists for the exchange of information about chemical reaction dynamics. We have introduced SBML support to the reaction diffusion module of the NEURON simulation software. This enables NEURON users to import a large number of previously developed cell biology models and use them in computational neuroscience research. We have added the capability to import and export models in SBML format to NEURON’s reaction-diffusion module using the libSBML library [3]. This allows one to combine models from ModelDB [4] with cell models from the BioModels database [5] in a way that makes it possible to match state variables and to readily change parameters. The procedure is as follows: • Electrophysiology model is loaded from ModelDB or constructed de novo. • NEURON loads SBML data and instantiates appropriate reaction diffusion objects: rxd.Region, rxd.Species, rxd.Reaction. User must interactively or preemptively match state-variable names across combined models to combine the two models into one. • SBML models do not include diffusion so diffusion constants can be added to provide the spatial element of the combined model. • User can adjust the parameters and make simulation runs of the model. We have used these new features to study the interaction of a calcium induced calcium release model from BioModels with electrical activity in a detailed model of a hippocampal CA1 pyramidal neuron.

Model exchange with the NeuroML model database

The Neural Open Markup Language project, NeuroML, is an international, collaborative initiative to develop a language for describing and sharing complex, multiscale neuron and neuronal network models [1]. The project focuses on the key objects that need to be exchanged among software applications used by computational neuroscientists. Examples of these objects include descriptions of neuronal morphology, the dynamics of ion channels and synaptic mechanisms, and the connectivity patterns of networks of model neurons. This modular approach brings additional benefits: not only can entire models be published and exchanged in this format, but each individual object or component, such as a specific calcium channel or excitatory synapse, can be shared and re-implemented in a different model. The use of a standardized description language based on XML also facilitates the development of tools that promote simulator interoperability. The NeuroML Model Database (NeuroML-db.org) is a relational database that provides a means for exchanging multiscale NeuroML model descriptions and their components. There are several existing, curated model databases that include neuroscience models in multiple formats such as ModelDB [2] and the BioModels Database [3], which provide excellent resources for sharing diverse, published neuroscience models in multiple formats. In contrast, NeuroML-db is more focused, including only NeuroML models. This emphasis allows the database design and search to take advantage of this specific format, and in particular, allows for efficient searches over sub-components of models. In addition, the NeuroML Model Database can be used to search over public models developed as part of the Open Source Brain (OSB) initiative, which provides a software infrastructure for the collaborative development and evaluation of models [4]. NeuroML is the preferred model description format for OSB, and OSB models that are described using NeuroML version 2.0 can be flagged for automatic inclusion in NeuroML-db. An overarching design goal for the database is to provide a simple keyword search interface. Because many search terms involve terminology from neuroscience, we complement a metadata-based keyword search with a mechanism that can provide information about the biological meaning of the keywords, where this semantic information is available through an existing ontology for neuroscience called NeuroLex [5]. Through implicit query reformulation, this semantic information is used to expand the search results. For example, if the query is a single brain region, all cell and network models for cells and networks from that region and its sub-regions are provided to the user.