BioSimulators: a central registry of simulation engines and services for recommending specific tools.

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.

eQuilibrator 3.0: a database solution for thermodynamic constant estimation.

eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

Long-Term Effects of a Web-Based Low-FODMAP Diet Versus Probiotic Treatment for Irritable Bowel Syndrome, Including Shotgun Analyses of Microbiota: Randomized, Double-Crossover Clinical Trial.

BACKGROUND: The long-term management of irritable bowel syndrome (IBS) poses many challenges. In short-term studies, eHealth interventions have been demonstrated to be safe and practical for at-home monitoring of the effects of probiotic treatments and a diet low in fermentable oligosaccharides, disaccharides, monosaccharides, and polyols (FODMAPs). IBS has been linked to alterations in the microbiota. OBJECTIVE: The aim of this study was to determine whether a web-based low-FODMAP diet (LFD) intervention and probiotic treatment were equally good at reducing IBS symptoms, and whether the response to treatments could be explained by patients' microbiota. METHODS: Adult IBS patients were enrolled in an open-label, randomized crossover trial (for nonresponders) with 1 year of follow-up using the web application IBS Constant Care (IBS CC). Patients were recruited from the outpatient clinic at the Department of Gastroenterology, North Zealand University Hospital, Denmark. Patients received either VSL#3 for 4 weeks (2 × 450 billion colony-forming units per day) or were placed on an LFD for 4 weeks. Patients responding to the LFD were reintroduced to foods high in FODMAPs, and probiotic responders received treatments whenever they experienced a flare-up of symptoms. Treatment response and symptom flare-ups were defined as a reduction or increase, respectively, of at least 50 points on the IBS Severity Scoring System (IBS-SSS). Web-based ward rounds were performed daily by the study investigator. Fecal microbiota were analyzed by shotgun metagenomic sequencing (at least 10 million 2 × 100 bp paired-end sequencing reads per sample). RESULTS: A total of 34 IBS patients without comorbidities and 6 healthy controls were enrolled in the study. Taken from participating subjects, 180 fecal samples were analyzed for their microbiota composition. Out of 21 IBS patients, 12 (57%) responded to the LFD and 8 (38%) completed the reintroduction of FODMAPs. Out of 21 patients, 13 (62%) responded to their first treatment of VSL#3 and 7 (33%) responded to multiple VSL#3 treatments. A median of 3 (IQR 2.25-3.75) probiotic treatments were needed for sustained symptom control. LFD responders were reintroduced to a median of 14.50 (IQR 7.25-21.75) high-FODMAP items. No significant difference in the median reduction of IBS-SSS for LFD versus probiotic responders was observed, where for LFD it was -126.50 (IQR -196.75 to -76.75) and for VSL#3 it was -130.00 (IQR -211.00 to -70.50; P>.99). Responses to either of the two treatments were not able to be predicted using patients' microbiota. CONCLUSIONS: The web-based LFD intervention and probiotic treatment were equally efficacious in managing IBS symptoms. The response to treatments could not be explained by the composition of the microbiota. The IBS CC web application was shown to be practical, safe, and useful for clinical decision making in the long-term management of IBS. Although this study was underpowered, findings from this study warrant further research in a larger sample of patients with IBS to confirm these long-term outcomes. TRIAL REGISTRATION: ClinicalTrials.gov NCT03586622; https://clinicaltrials.gov/ct2/show/NCT03586622.

multiTFA: a Python package for multi-variate thermodynamics-based flux analysis.

MOTIVATION: We achieve a significant improvement in thermodynamic-based flux analysis (TFA) by introducing multivariate treatment of thermodynamic variables and leveraging component contribution, the state-of-the-art implementation of the group contribution methodology. Overall, the method greatly reduces the uncertainty of thermodynamic variables. RESULTS: We present multiTFA, a Python implementation of our framework. We evaluated our application using the core Escherichia coli model and achieved a median reduction of 6.8 kJ/mol in reaction Gibbs free energy ranges, while three out of 12 reactions in glycolysis changed from reversible to irreversible. AVAILABILITY AND IMPLEMENTATION: Our framework along with documentation is available on https://github.com/biosustain/multitfa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes.

For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical 'Rosetta Stone' to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org/biochem and KBase.

Author Correction: A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Community standards to facilitate development and address challenges in metabolic modeling.

Standardization of data and models facilitates effective communication, especially in computational systems biology. However, both the development and consistent use of standards and resources remain challenging. As a result, the amount, quality, and format of the information contained within systems biology models are not consistent and therefore present challenges for widespread use and communication. Here, we focused on these standards, resources, and challenges in the field of constraint-based metabolic modeling by conducting a community-wide survey. We used this feedback to (i) outline the major challenges that our field faces and to propose solutions and (ii) identify a set of features that defines what a "gold standard" metabolic network reconstruction looks like concerning content, annotation, and simulation capabilities. We anticipate that this community-driven outline will help the long-term development of community-inspired resources as well as produce high-quality, accessible models within our field. More broadly, we hope that these efforts can serve as blueprints for other computational modeling communities to ensure the continued development of both practical, usable standards and reproducible, knowledge-rich models.

Publisher Correction: MEMOTE for standardized genome-scale metabolic model testing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism.

Genome-scale metabolic models (GEMs) represent extensive knowledgebases that provide a platform for model simulations and integrative analysis of omics data. This study introduces Yeast8 and an associated ecosystem of models that represent a comprehensive computational resource for performing simulations of the metabolism of Saccharomyces cerevisiae--an important model organism and widely used cell-factory. Yeast8 tracks community development with version control, setting a standard for how GEMs can be continuously updated in a simple and reproducible way. We use Yeast8 to develop the derived models panYeast8 and coreYeast8, which in turn enable the reconstruction of GEMs for 1,011 different yeast strains. Through integration with enzyme constraints (ecYeast8) and protein 3D structures (proYeast8DB), Yeast8 further facilitates the exploration of yeast metabolism at a multi-scale level, enabling prediction of how single nucleotide variations translate to phenotypic traits.

Cameo: A Python Library for Computer Aided Metabolic Engineering and Optimization of Cell Factories.

Computational systems biology methods enable rational design of cell factories on a genome-scale and thus accelerate the engineering of cells for the production of valuable chemicals and proteins. Unfortunately, the majority of these methods' implementations are either not published, rely on proprietary software, or do not provide documented interfaces, which has precluded their mainstream adoption in the field. In this work we present cameo, a platform-independent software that enables in silico design of cell factories and targets both experienced modelers as well as users new to the field. It is written in Python and implements state-of-the-art methods for enumerating and prioritizing knockout, knock-in, overexpression, and down-regulation strategies and combinations thereof. Cameo is an open source software project and is freely available under the Apache License 2.0. A dedicated Web site including documentation, examples, and installation instructions can be found at http://cameo.bio . Users can also give cameo a try at http://try.cameo.bio .

Effect of database drift on network topology and enrichment analyses: a case study for RegulonDB.

RegulonDB is a database storing the biological information behind the transcriptional regulatory network (TRN) of the bacterium Escherichia coli. It is one of the key bioinformatics resources for Systems Biology investigations of bacterial gene regulation. Like most biological databases, the content drifts with time, both due to the accumulation of new information and due to refinements in the underlying biological concepts. Conclusions based on previous database versions may no longer hold. Here, we study the change of some topological properties of the TRN of E. coli, as provided by RegulonDB across 16 versions, as well as a simple index, digital control strength, quantifying the match between gene expression profiles and the transcriptional regulatory networks. While many of network characteristics change dramatically across the different versions, the digital control strength remains rather robust and in tune with previous results for this index. Our study shows that: (i) results derived from network topology should, when possible, be studied across a range of database versions, before detailed biological conclusions are derived, and (ii) resorting to simple indices, when interpreting high-throughput data from a network perspective, may help achieving a robustness of the findings against variation of the underlying biological information. Database URL: www.regulondb.ccg.unam.mx.

Artefacts in statistical analyses of network motifs: general framework and application to metabolic networks.

Few-node subgraphs are the smallest collective units in a network that can be investigated. They are beyond the scale of individual nodes but more local than, for example, communities. When statistically over- or under-represented, they are called network motifs. Network motifs have been interpreted as building blocks that shape the dynamic behaviour of networks. It is this promise of potentially explaining emergent properties of complex systems with relatively simple structures that led to an interest in network motifs in an ever-growing number of studies and across disciplines. Here, we discuss artefacts in the analysis of network motifs arising from discrepancies between the network under investigation and the pool of random graphs serving as a null model. Our aim was to provide a clear and accessible catalogue of such incongruities and their effect on the motif signature. As a case study, we explore the metabolic network of Escherichia coli and show that only by excluding ever more artefacts from the motif signature a strong and plausible correlation with the essentiality profile of metabolic reactions emerges.