SUMMARY: The human body harbours a plethora of microbes that play a fundamental role in the well-being of the host. Still, the contribution of many microorganisms to human health remains undiscovered. To understand the composition of their communities, the accurate genome-scale metabolic network models of participating microorganisms are integrated to construct a community that mimics the normal bacterial flora of humans. So far, tools for modelling the communities have transformed the community into various optimization problems and model compositions. Therefore, any knockout or modification of each submodel (each species) necessitates the up-to-date creation of the community to incorporate rebuildings. To solve this complexity, we refer to the context of SBML in a hierarchical model composition, wherein each species's genome-scale metabolic model is imported as a submodel in another model. Hence, the community is a model composed of submodels defined in separate files. We combine all these files upon parsing to a so-called 'flattened' model, i.e., a comprehensive and valid SBML file of the entire community that COBRApy can parse for further processing. The hierarchical model facilitates the analysis of the whole community irrespective of any changes in the individual submodels. AVAILABILITY AND IMPLEMENTATION: The module is freely available at https://github.com/manuelgloeckler/ncmw.

MOTIVATION: The number and size of computational models in biology have drastically increased over the past years and continue to grow. Modeled networks are becoming more complex, and reconstructing them from the beginning in an exchangeable and reproducible manner is challenging. Using precisely defined ontologies enables the encoding of field-specific knowledge and the association of disparate data types. In computational modeling, the medium for representing domain knowledge is the set of orthogonal structured controlled vocabularies named Systems Biology Ontology (SBO). The SBO terms enable modelers to explicitly define and describe model entities, including their roles and characteristics. RESULTS: Here, we present the first standalone tool that automatically assigns SBO terms to multiple entities of a given SBML model, named the SBOannotator. The main focus lies on the reactions, as the correct assignment of precise SBO annotations requires their extensive classification. Our implementation does not consider only top-level terms but examines the functionality of the underlying enzymes to allocate precise and highly specific ontology terms to biochemical reactions. Transport reactions are examined separately and are classified based on the mechanism of molecule transport. Pseudo-reactions that serve modeling purposes are given reasonable terms to distinguish between biomass production and the import or export of metabolites. Finally, other model entities, such as metabolites and genes, are annotated with appropriate terms. Including SBO annotations in the models will enhance the reproducibility, usability, and analysis of biochemical networks. AVAILABILITY AND IMPLEMENTATION: SBOannotator is freely available from https://github.com/draeger-lab/SBOannotator/.

Abstract Summary The Computational Modelling of Systems Biology (SysMod) Community of Special Interest (COSI) convenes annually at the Intelligent Systems for Molecular Biology (ISMB) conference to facilitate knowledge dissemination and exchange of research findings on systems modelling from interdisciplinary domains. The SysMod meeting 2022 was held in a hybrid mode in Madison, Wisconsin, spanning a 1-day duration centred on modelling techniques, applications, and single-cell technology implementations. The meeting showcased innovative approaches to modelling biological systems using cell-specific and multiscale modelling, multiomics data integration, and novel tools to develop systems models using single-cell and multiomics technology. The meeting also recognized outstanding research by awarding the three best posters. This report summarizes the key highlights and outcomes of the meeting. Availability and implementation All resources and further information are freely accessible at https://sysmod.info.

COVID-19 is one of the deadliest respiratory diseases, and its emergence caught the pharmaceutical industry off guard. While vaccines have been rapidly developed, treatment options for infected people remain scarce, and COVID-19 poses a substantial global threat. This study presents a novel workflow to predict robust druggable targets against emerging RNA viruses using metabolic networks and information of the viral structure and its genome sequence. For this purpose, we implemented pymCADRE and PREDICATE to create tissue-specific metabolic models, construct viral biomass functions and predict host-based antiviral targets from more than one genome. We observed that pymCADRE reduces the computational time of flux variability analysis for internal optimizations. We applied these tools to create a new metabolic network of primary bronchial epithelial cells infected with SARS-CoV-2 and identified enzymatic reactions with inhibitory effects. The most promising reported targets were from the purine metabolism, while targeting the pyrimidine and carbohydrate metabolisms seemed to be promising approaches to enhance viral inhibition. Finally, we computationally tested the robustness of our targets in all known variants of concern, verifying our targets' inhibitory effects. Since laboratory tests are time-consuming and involve complex readouts to track processes, our workflow focuses on metabolic fluxes within infected cells and is applicable for rapid hypothesis-driven identification of potentially exploitable antivirals concerning various viruses and host cell types.

Abstract Summary The Community of Special Interest (COSI) in Computational Modelling of Biological Systems (SysMod) brings together interdisciplinary scientists interested in combining data-driven computational modelling, multi-scale mechanistic frameworks, large-scale -omics data and bioinformatics. SysMod’s main activity is an annual meeting at the Intelligent Systems for Molecular Biology (ISMB) conference, a meeting for computer scientists, biologists, mathematicians, engineers and computational and systems biologists. The 2021 SysMod meeting was conducted virtually due to the ongoing COVID-19 pandemic (coronavirus disease 2019). During the 2-day meeting, the development of computational tools, approaches and predictive models was discussed, along with their application to biological systems, emphasizing disease mechanisms. This report summarizes the meeting. Availability and implementation All resources and further information are freely accessible at https://sysmod.info.

SUMMARY: Studying biological systems generally relies on computational modeling and simulation, e.g., model-driven discovery and hypothesis testing. Progress in standardization efforts led to the development of interrelated file formats to exchange and reuse models in systems biology, such as SBML, the Simulation Experiment Description Markup Language (SED-ML) or the Open Modeling EXchange format. Conducting simulation experiments based on these formats requires efficient and reusable implementations to make them accessible to the broader scientific community and to ensure the reproducibility of the results. The Systems Biology Simulation Core Library (SBSCL) provides interpreters and solvers for these standards as a versatile open-source API in JavaTM. The library simulates even complex bio-models and supports deterministic Ordinary Differential Equations; Stochastic Differential Equations; constraint-based analyses; recent SBML and SED-ML versions; exchange of results, and visualization of in silico experiments; open modeling exchange formats (COMBINE archives); hierarchically structured models; and compatibility with standard testing systems, including the Systems Biology Test Suite and published models from the BioModels and BiGG databases. AVAILABILITY AND IMPLEMENTATION: SBSCL is freely available at https://draeger-lab.github.io/SBSCL/ and via Maven Central. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational models of biological systems can exploit a broad range of rapidly developing approaches, including novel experimental approaches, bioinformatics data analysis, emerging modelling paradigms, data standards and algorithms. A discussion about the most recent advances among experts from various domains is crucial to foster data-driven computational modelling and its growing use in assessing and predicting the behaviour of biological systems. Intending to encourage the development of tools, approaches and predictive models, and to deepen our understanding of biological systems, the Community of Special Interest (COSI) was launched in Computational Modelling of Biological Systems (SysMod) in 2016. SysMod's main activity is an annual meeting at the Intelligent Systems for Molecular Biology (ISMB) conference, which brings together computer scientists, biologists, mathematicians, engineers, computational and systems biologists. In the five years since its inception, SysMod has evolved into a dynamic and expanding community, as the increasing number of contributions and participants illustrate. SysMod maintains several online resources to facilitate interaction among the community members, including an online forum, a calendar of relevant meetings and a YouTube channel with talks and lectures of interest for the modelling community. For more than half a decade, the growing interest in computational systems modelling and multi-scale data integration has inspired and supported the SysMod community. Its members get progressively more involved and actively contribute to the annual COSI meeting and several related community workshops and meetings, focusing on specific topics, including particular techniques for computational modelling or standardisation efforts.

MOTIVATION: The novel coronavirus (SARS-CoV-2) currently spreads worldwide, causing the disease COVID-19. The number of infections increases daily, without any approved antiviral therapy. The recently released viral nucleotide sequence enables the identification of therapeutic targets, e.g. by analyzing integrated human-virus metabolic models. Investigations of changed metabolic processes after virus infections and the effect of knock-outs on the host and the virus can reveal new potential targets. RESULTS: We generated an integrated host-virus genome-scale metabolic model of human alveolar macrophages and SARS-CoV-2. Analyses of stoichiometric and metabolic changes between uninfected and infected host cells using flux balance analysis (FBA) highlighted the different requirements of host and virus. Consequently, alterations in the metabolism can have different effects on host and virus, leading to potential antiviral targets. One of these potential targets is guanylate kinase (GK1). In FBA analyses, the knock-out of the GK1 decreased the growth of the virus to zero, while not affecting the host. As GK1 inhibitors are described in the literature, its potential therapeutic effect for SARS-CoV-2 infections needs to be verified in in-vitro experiments. AVAILABILITY AND IMPLEMENTATION: The computational model is accessible at https://identifiers.org/biomodels.db/MODEL2003020001.

BACKGROUND: New technologies have given rise to an abundance of -omics data, particularly metabolomic data. The scale of these data introduces new challenges for the interpretation and extraction of knowledge, requiring the development of innovative computational visualization methodologies. Here, we present GEM-Vis, an original method for the visualization of time-course metabolomic data within the context of metabolic network maps. We demonstrate the utility of the GEM-Vis method by examining previously published data for two cellular systems-the human platelet and erythrocyte under cold storage for use in transfusion medicine. RESULTS: The results comprise two animated videos that allow for new insights into the metabolic state of both cell types. In the case study of the platelet metabolome during storage, the new visualization technique elucidates a nicotinamide accumulation that mirrors that of hypoxanthine and might, therefore, reflect similar pathway usage. This visual analysis provides a possible explanation for why the salvage reactions in purine metabolism exhibit lower activity during the first few days of the storage period. The second case study displays drastic changes in specific erythrocyte metabolite pools at different times during storage at different temperatures. CONCLUSIONS: The new visualization technique GEM-Vis introduced in this article constitutes a well-suitable approach for large-scale network exploration and advances hypothesis generation. This method can be applied to any system with data and a metabolic map to promote visualization and understand physiology at the network level. More broadly, we hope that our approach will provide the blueprints for new visualizations of other longitudinal -omics data types. The supplement includes a comprehensive user's guide and links to a series of tutorial videos that explain how to prepare model and data files, and how to use the software SBMLsimulator in combination with further tools to create similar animations as highlighted in the case studies.

Life science researchers use computational models to articulate and test hypotheses about the behavior of biological systems. Semantic annotation is a critical component for enhancing the interoperability and reusability of such models as well as for the integration of the data needed for model parameterization and validation. Encoded as machine-readable links to knowledge resource terms, semantic annotations describe the computational or biological meaning of what models and data represent. These annotations help researchers find and repurpose models, accelerate model composition and enable knowledge integration across model repositories and experimental data stores. However, realizing the potential benefits of semantic annotation requires the development of model annotation standards that adhere to a community-based annotation protocol. Without such standards, tool developers must account for a variety of annotation formats and approaches, a situation that can become prohibitively cumbersome and which can defeat the purpose of linking model elements to controlled knowledge resource terms. Currently, no consensus protocol for semantic annotation exists among the larger biological modeling community. Here, we report on the landscape of current annotation practices among the COmputational Modeling in BIology NEtwork community and provide a set of recommendations for building a consensus approach to semantic annotation.

The 17th International NETTAB workshop was held in Palermo, Italy, on October 16-18, 2017. The special topic for the meeting was "Methods, tools and platforms for Personalised Medicine in the Big Data Era", but the traditional topics of the meeting series were also included in the event. About 40 scientific contributions were presented, including four keynote lectures, five guest lectures, and many oral communications and posters. Also, three tutorials were organised before and after the workshop. Full papers from some of the best works presented in Palermo were submitted for this Supplement of BMC Bioinformatics. Here, we provide an overview of meeting aims and scope. We also shortly introduce selected papers that have been accepted for publication in this Supplement, for a complete presentation of the outcomes of the meeting.

During various inflammatory processes circulating cytokines including IL-6, IL-1β, and TNFα elicit a broad and clinically relevant impairment of hepatic detoxification that is based on the simultaneous downregulation of many drug metabolizing enzymes and transporter genes. To address the question whether a common mechanism is involved we treated human primary hepatocytes with IL-6, the major mediator of the acute phase response in liver, and characterized acute phase and detoxification responses in quantitative gene expression and (phospho-)proteomics data sets. Selective inhibitors were used to disentangle the roles of JAK/STAT, MAPK, and PI3K signaling pathways. A prior knowledge-based fuzzy logic model comprising signal transduction and gene regulation was established and trained with perturbation-derived gene expression data from five hepatocyte donors. Our model suggests a greater role of MAPK/PI3K compared to JAK/STAT with the orphan nuclear receptor RXRα playing a central role in mediating transcriptional downregulation. Validation experiments revealed a striking similarity of RXRα gene silencing versus IL-6 induced negative gene regulation (rs = 0.79; P<0.0001). These results concur with RXRα functioning as obligatory heterodimerization partner for several nuclear receptors that regulate drug and lipid metabolism.

UNLABELLED: JSBML, the official pure Java programming library for the Systems Biology Markup Language (SBML) format, has evolved with the advent of different modeling formalisms in systems biology and their ability to be exchanged and represented via extensions of SBML. JSBML has matured into a major, active open-source project with contributions from a growing, international team of developers who not only maintain compatibility with SBML, but also drive steady improvements to the Java interface and promote ease-of-use with end users. AVAILABILITY AND IMPLEMENTATION: Source code, binaries and documentation for JSBML can be freely obtained under the terms of the LGPL 2.1 from the website http://sbml.org/Software/JSBML. More information about JSBML can be found in the user guide at http://sbml.org/Software/JSBML/docs/. CONTACT: [email protected] or [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Escher is a web application for visualizing data on biological pathways. Three key features make Escher a uniquely effective tool for pathway visualization. First, users can rapidly design new pathway maps. Escher provides pathway suggestions based on user data and genome-scale models, so users can draw pathways in a semi-automated way. Second, users can visualize data related to genes or proteins on the associated reactions and pathways, using rules that define which enzymes catalyze each reaction. Thus, users can identify trends in common genomic data types (e.g. RNA-Seq, proteomics, ChIP)--in conjunction with metabolite- and reaction-oriented data types (e.g. metabolomics, fluxomics). Third, Escher harnesses the strengths of web technologies (SVG, D3, developer tools) so that visualizations can be rapidly adapted, extended, shared, and embedded. This paper provides examples of each of these features and explains how the development approach used for Escher can be used to guide the development of future visualization tools.

UNLABELLED: GRN2SBML automatically encodes gene regulatory networks derived from several inference tools in systems biology markup language. Providing a graphical user interface, the networks can be annotated via the simple object access protocol (SOAP)-based application programming interface of BioMart Central Portal and minimum information required in the annotation of models registry. Additionally, we provide an R-package, which processes the output of supported inference algorithms and automatically passes all required parameters to GRN2SBML. Therefore, GRN2SBML closes a gap in the processing pipeline between the inference of gene regulatory networks and their subsequent analysis, visualization and storage. AVAILABILITY: GRN2SBML is freely available under the GNU Public License version 3 and can be downloaded from http://www.hki-jena.de/index.php/0/2/490. SUPPLEMENTARY INFORMATION: General information on GRN2SBML, examples and tutorials are available at the tool's web page.

MOTIVATION: The biological pathway exchange language (BioPAX) and the systems biology markup language (SBML) belong to the most popular modeling and data exchange languages in systems biology. The focus of SBML is quantitative modeling and dynamic simulation of models, whereas the BioPAX specification concentrates mainly on visualization and qualitative analysis of pathway maps. BioPAX describes reactions and relations. In contrast, SBML core exclusively describes quantitative processes such as reactions. With the SBML qualitative models extension (qual), it has recently also become possible to describe relations in SBML. Before the development of SBML qual, relations could not be properly translated into SBML. Until now, there exists no BioPAX to SBML converter that is fully capable of translating both reactions and relations. RESULTS: The entire nature pathway interaction database has been converted from BioPAX (Level 2 and Level 3) into SBML (Level 3 Version 1) including both reactions and relations by using the new qual extension package. Additionally, we present the new webtool BioPAX2SBML for further BioPAX to SBML conversions. Compared with previous conversion tools, BioPAX2SBML is more comprehensive, more robust and more exact. AVAILABILITY: BioPAX2SBML is freely available at http://webservices.cs.uni-tuebingen.de/ and the complete collection of the PID models is available at http://www.cogsys.cs.uni-tuebingen.de/downloads/Qualitative-Models/.

SUMMARY: CySBML is a plugin designed to work with Systems Biology Markup Language (SBML) in Cytoscape having the following features: SBML import, support of the SBML layout and qualitative model packages, navigation in network layouts based on SBML structure, access to MIRIAM and SBO-based annotations and SBML validation. CySBML includes an importer for BioModels to load SBML from standard repositories. AVAILABILITY AND IMPLEMENTATION: Freely available for non-commercial purposes through the Cytoscape plugin manager or for download at http://sourceforge.net/projects/cysbml/. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Tutorial, usage guide, installation instructions and additional figures are available for download at http://www.charite.de/sysbio/people/koenig/software/cysbml/.

Higher organisms are able to respond to continuously changing external conditions by transducing cellular signals into specific regulatory programs, which control gene expression states of thousands of different genes. One of the central problems in understanding gene regulation is to decipher how combinations of transcription factors control sets of co expressed genes under specific experimental conditions. Existing methods in this field mainly focus on sequence aspects and pattern recognition, e.g., by detecting cis-regulatory modules (CRMs) based on gene expression profiling data. We propose a novel approach by combining experimental data with a priori knowledge of respective experimental conditions. These various sources of evidence are likewise considered using multi-objective evolutionary optimization. In this work, we present three objective functions that are especially designed for stimulus-response experiments and can be used to integrate a priori knowledge into the detection of gene regulatory modules. This method was tested and evaluated on whole-genome microarray measurements of drug-response in human hepatocytes.

SUMMARY: The specifications of the Systems Biology Markup Language (SBML) define standards for storing and exchanging computer models of biological processes in text files. In order to perform model simulations, graphical visualizations and other software manipulations, an in-memory representation of SBML is required. We developed JSBML for this purpose. In contrast to prior implementations of SBML APIs, JSBML has been designed from the ground up for the Java programming language, and can therefore be used on all platforms supported by a Java Runtime Environment. This offers important benefits for Java users, including the ability to distribute software as Java Web Start applications. JSBML supports all SBML Levels and Versions through Level 3 Version 1, and we have strived to maintain the highest possible degree of compatibility with the popular library libSBML. JSBML also supports modules that can facilitate the development of plugins for end user applications, as well as ease migration from a libSBML-based backend. AVAILABILITY: Source code, binaries and documentation for JSBML can be freely obtained under the terms of the LGPL 2.1 from the website http://sbml.org/Software/JSBML.

MOTIVATION: Statins are the most widely used cholesterol-lowering drugs. The primary target of statins is HMG-CoA reductase, a key enzyme in cholesterol synthesis. However, statins elicit pleitropic responses including beneficial as well as adverse effects in the liver or other organs. Today, the regulatory mechanisms that cause these pleiotropic effects are not sufficiently understood. RESULTS: In this work, genome-wide RNA expression changes in primary human hepatocytes of six individuals were measured at up to six time points upon atorvastatin treatment. A computational analysis workflow was applied to reconstruct regulatory mechanisms based on these drug-response data and available knowledge about transcription factor (TF) binding specificities and protein-drug interactions. Several previously unknown TFs were predicted to be involved in atorvastatin-responsive gene expression. The novel relationships of nuclear receptors NR2C2 and PPARA on CYP3A4 were successfully validated in wet-lab experiments. AVAILABILITY: Microarray data are available at the Gene Expression Omnibus (GEO) database at www.ncbi.nlm.nih.gov/geo/, under accession number GSE29868. CONTACT: [email protected]; [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SUMMARY: The KEGG PATHWAY database provides a widely used service for metabolic and nonmetabolic pathways. It contains manually drawn pathway maps with information about the genes, reactions and relations contained therein. To store these pathways, KEGG uses KGML, a proprietary XML-format. Parsers and translators are needed to process the pathway maps for usage in other applications and algorithms. We have developed KEGGtranslator, an easy-to-use stand-alone application that can visualize and convert KGML formatted XML-files into multiple output formats. Unlike other translators, KEGGtranslator supports a plethora of output formats, is able to augment the information in translated documents (e.g. MIRIAM annotations) beyond the scope of the KGML document, and amends missing components to fragmentary reactions within the pathway to allow simulations on those. AVAILABILITY: KEGGtranslator is freely available as a Java(™) Web Start application and for download at http://www.cogsys.cs.uni-tuebingen.de/software/KEGGtranslator/. KGML files can be downloaded from within the application. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SUMMARY: The XML-based Systems Biology Markup Language (SBML) has emerged as a standard for storage, communication and interchange of models in systems biology. As a machine-readable format XML is difficult for humans to read and understand. Many tools are available that visualize the reaction pathways stored in SBML files, but many components, e.g. unit declarations, complex kinetic equations or links to MIRIAM resources, are often not made visible in these diagrams. For a broader understanding of the models, support in scientific writing and error detection, a human-readable report of the complete model is needed. We present SBML2L(A)T(E)X, a Java-based stand-alone program to fill this gap. A convenient web service allows users to directly convert SBML to various formats, including DVI, L(A)T(E)X and PDF, and provides many settings for customization. AVAILABILITY: Source code, documentation and a web service are freely available at (http://www.ra.cs.uni-tuebingen.de/software/SBML2LaTeX).

SUMMARY: BioJava is a mature open-source project that provides a framework for processing of biological data. BioJava contains powerful analysis and statistical routines, tools for parsing common file formats and packages for manipulating sequences and 3D structures. It enables rapid bioinformatics application development in the Java programming language. AVAILABILITY: BioJava is an open-source project distributed under the Lesser GPL (LGPL). BioJava can be downloaded from the BioJava website (http://www.biojava.org). BioJava requires Java 1.5 or higher. All queries should be directed to the BioJava mailing lists. Details are available at http://biojava.org/wiki/BioJava:MailingLists.

An important problem in systems biology is parameter estimation for biochemical system models. Our work concentrates on the metabolic subnetwork of the valine and leucine biosynthesis in Corynebacterium glutamicum, an anaerobic actinobacterium of high biotechnological importance. Using data of an in vivo experiment measuring 13 metabolites during a glucose stimulus-response experiment we investigate the performance of various evolutionary algorithms on the parameter inference problem in biochemical modeling. Due to the inconclusive information on the reversibility of the reactions in the pathway, we develop both a reversible and an irreversible differential equation model based on the recent convenience kinetics approach. As the reversible model allows better approximation on the whole, we use it to analyze the impact of different settings on four especially promising EAs. We show that Particle Swarm Optimization as well as Differential Evolution are useful methods for parameter estimation on convenience kinetics models outperforming Genetic Algorithm and Evolution Strategy approaches and nearly reaching the quality of independent spline approximations on the raw data.

Parameter estimation for biochemical model systems has become an important problem in systems biology. Here we focus on the metabolic subnetwork of the valine and leucine biosynthesis in C. glutamicum. Due to the lack of indisputable information regarding reversibility of the reactions in the pathway we derived two alternative ordinary differential equation models based on the formalisms of the generalized mass-action rate law. We introduced two alternative modeling approaches for feedback inhibition and evaluated the applicability of six optimization procedures (multi start Hill Climber, binary and real valued Genetic Algorithm, standard and covariance matrix adaption Evolution Strategy as well as Simulated Annealing) to the problem of parameter fitting. The model considering irreversible reactions performed worse and was therefore rejected from further analysis. We benchmarked the impact of different mutation and crossover operators as well as the influence of the population size on the remaining system and the two best optimization procedures namely binary Genetic Algorithm and the Evolution Strategy. The GA performed best on average and found the best total result based on the relative squared error.