SBMLsqueezer 2: context-sensitive creation of kinetic equations in biochemical networks

BioNexusSentinel: a visual tool for bioregulatory network and cytohistological RNA-seq genetic expression profiling within the context of multicellular simulation research using ChatGPT-augmented software engineering

Summary

Motivated by the need to parameterize ongoing multicellular simulation research, this paper documents the culmination of a ChatGPT augmented software engineering cycle resulting in an integrated visual platform for efficient cytohistological RNA-seq and bioregulatory network exploration. As contrasted to other systems and synthetic biology tools, BioNexusSentinel was developed de novo to uniquely combine these features. Reactome served as the primary source of remotely accessible biological models, accessible using BioNexusSentinel's novel search engine and REST API requests. The innovative, feature-rich gene expression profiler component was developed to enhance the exploratory experience for the researcher, culminating in the cytohistological RNA-seq explorer based on Human Protein Atlas data. A novel cytohistological classifier would be integrated via pre-processed analysis of the RNA-seq data via R statistical language, providing for useful analytical functionality and good performance for the end-user. Implications of the work span prospects for model orthogonality evaluations, gap identification in network modelling, prototyped automatic kinetics parameterization, and downstream simulation and cellular biological state analysis. This unique computational biology software engineering collaboration with generative natural language processing artificial intelligence was shown to enhance worker productivity, with evident benefits in terms of accelerating coding and machine-human intelligence transfer.

Availability and implementation

BioNexusSentinel project releases, with corresponding data and installation instructions, are available at https://github.com/RichardMatzko/BioNexusSentinel.

An automated model annotation system (AMAS) for SBML models

MOTIVATION

Annotations of biochemical models provide details of chemical species, documentation of chemical reactions, and other essential information. Unfortunately, the vast majority of biochemical models have few, if any, annotations, or the annotations provide insufficient detail to understand the limitations of the model. The quality and quantity of annotations can be improved by developing tools that recommend annotations. For example, recommender tools have been developed for annotations of genes. Although annotating genes is conceptually similar to annotating biochemical models, there are important technical differences that make it difficult to directly apply this prior work.

RESULTS

We present AMAS, a system that predicts annotations for elements of models represented in the Systems Biology Markup Language (SBML) community standard. We provide a general framework for predicting model annotations for a query element based on a database of annotated reference elements and a match score function that calculates the similarity between the query element and reference elements. The framework is instantiated to specific element types (e.g. species, reactions) by specifying the reference database (e.g. ChEBI for species) and the match score function (e.g. string similarity). We analyze the computational efficiency and prediction quality of AMAS for species and reactions in BiGG and BioModels and find that it has subsecond response times and accuracy between 80% and 95% depending on specifics of what is predicted. We have incorporated AMAS into an open-source, pip-installable Python package that can run as a command-line tool that predicts and adds annotations to species and reactions to an SBML model.

AVAILABILITY AND IMPLEMENTATION

Our project is hosted at https://github.com/sys-bio/AMAS, where we provide examples, documentation, and source code files. Our source code is licensed under the MIT open-source license.

An Automated Model Annotation System (AMAS) for SBML Models

Motivation

Annotations of biochemical models provide details of chemical species, documentation of chemical reactions, and other essential information. Unfortunately, the vast majority of biochemical models have few, if any, annotations, or the annotations provide insufficient detail to understand the limitations of the model. The quality and quantity of annotations can be improved by developing tools that recommend annotations. For example, recommender tools have been developed for annotations of genes. Although annotating genes is conceptually similar to annotating biochemical models, there are important technical differences that make it difficult to directly apply this prior work.

Results

We present AMAS, a system that predicts annotations for elements of models represented in the Systems Biology Markup Language (SBML) community standard. We provide a general framework for predicting model annotations for a query element based on a database of annotated reference elements and a match score function that calculates the similarity between the query element and reference elements. The framework is instantiated to specific element types (e.g., species, reactions) by specifying the reference database (e.g., ChEBI for species) and the match score function (e.g., string similarity). We analyze the computational efficiency and prediction quality of AMAS for species and reactions in BiGG and BioModels and find that it has sub-second response times and accuracy between 80% and 95% depending on specifics of what is predicted. We have incorporated AMAS into an open-source, pip-installable Python package that can run as a command-line tool that predicts and adds annotations to species and reactions to an SBML model.

Availability

Our project is hosted at https://github.com/sys-bio/AMAS, where we provide examples, documentation, and source code files. Our source code is licensed under the MIT open-source license.

SBMLKinetics: a tool for annotation-independent classification of reaction kinetics for SBML models

BACKGROUND

Reaction networks are widely used as mechanistic models in systems biology to reveal principles of biological systems. Reactions are governed by kinetic laws that describe reaction rates. Selecting the appropriate kinetic laws is difficult for many modelers. There exist tools that attempt to find the correct kinetic laws based on annotations. Here, I developed annotation-independent technologies that assist modelers by focusing on finding kinetic laws commonly used for similar reactions.

RESULTS

Recommending kinetic laws and other analyses of reaction networks can be viewed as a classification problem. Existing approaches to determining similar reactions rely heavily on having good annotations, a condition that is often unsatisfied in model repositories such as BioModels. I developed an annotation-independent approach to find similar reactions via reaction classifications. I proposed a two-dimensional kinetics classification scheme (2DK) that analyzed reactions along the dimensions of kinetics type (K type) and reaction type (R type). I identified approximately ten mutually exclusive K types, including zeroth order, mass action, Michaelis-Menten, Hill kinetics, and others. R types were organized by the number of distinct reactants and the number of distinct products in reactions. I constructed a tool, SBMLKinetics, that inputted a collection of SBML models and then calculated reaction classifications as the probability of each 2DK class. The effectiveness of 2DK was evaluated on BioModels, and the scheme classified over 95% of the reactions.

CONCLUSIONS

2DK had many applications. It provided a data-driven annotation-independent approach to recommending kinetic laws by using type common for the kind of models in combination with the R type of the reactions. Alternatively, 2DK could also be used to alert users that a kinetic law was unusual for the K type and R type. Last, 2DK provided a way to analyze groups of models to compare their kinetic laws. I applied 2DK to BioModels to compare the kinetics of signaling networks with the kinetics of metabolic networks and found significant differences in K type distributions.

System biology analysis of endosulfan biodegradation in bacteria and its effect in other living systems: modeling and simulation studies

Endosulfan is a broadly applied cyclodiene insecticide which has been in use across 80 countries since last 5 decades. Owing to its recalcitrant nature, endosulfan residues have been reported from air, water and soil causing toxicity to various non-target organisms. Microbial decontamination of endosulfan has been reported previously by several authors. In the current study, we have evaluated the pathways of endosulfan degradation and its hazardous impact on other living beings including insects, humans, plants, aquatic life and environment by in-silico methods. For establishment of the endosulfan metabolism in different ecosystems, cell designer was employed. The established model was thereafter assessed and simulated to understand the biochemical and physiological metabolism of the endosulfan in various systems of the network. Topological investigation analysis of the endosulfan metabolism validated the presence of 207 nodes and 274 edges in the network. We have concluded that biomagnification of the endosulfan generally occurs in the various elements of the ecosystem. Dynamics study of endosulfan degrading enzymes suggested the important role of monooxygenase I, II and hydrolase in endosulfan bioremediation. Endosulfan shows toxicity in human beings, fishes and plants, however it is biodegraded by the microbes. To date, there are no reports of in- silico analysis of bioremediation of endosulfan and its hazardous effects on the environment. Thus, this report can be important in terms of modelling and simulation of biodegradation network of endosulfan and similar compounds and their impact on several other systems.Communicated by Ramaswamy H. Sarma.

Symbolic kinetic models in python (SKiMpy): intuitive modeling of large-scale biological kinetic models

MOTIVATION

Large-scale kinetic models are an invaluable tool to understand the dynamic and adaptive responses of biological systems. The development and application of these models have been limited by the availability of computational tools to build and analyze large-scale models efficiently. The toolbox presented here provides the means to implement, parameterize and analyze large-scale kinetic models intuitively and efficiently.

RESULTS

We present a Python package (SKiMpy) bridging this gap by implementing an efficient kinetic modeling toolbox for the semiautomatic generation and analysis of large-scale kinetic models for various biological domains such as signaling, gene expression and metabolism. Furthermore, we demonstrate how this toolbox is used to parameterize kinetic models around a steady-state reference efficiently. Finally, we show how SKiMpy can implement multispecies bioreactor simulations to assess biotechnological processes.

AVAILABILITY AND IMPLEMENTATION

The software is available as a Python 3 package on GitHub: https://github.com/EPFL-LCSB/SKiMpy, along with adequate documentation.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modeling and simulation of atrazine biodegradation in bacteria and its effect in other living systems

Atrazine is the most commonly used herbicide worldwide in the agricultural system. The increased environmental concentration of the atrazine showed the toxic effects on the non-target living species. Biodegradation of the atrazine is possible with the bacterial systems. The present study investigated biodegradation potential of atrazine degrading bacteria and the impact of atrazine on environmental systems. Model of atrazine fate in ecological systems constructed using the cell designer. The used model further analyzed and simulated to know the biochemistry and physiology of the atrazine in different cellular networks. Topological analysis of the atrazine degradation confirmed the 289 nodes and 300 edges. Our results showed that the overall biomagnification of the atrazine in the different environmental systems. Atrazine is showing toxic effects on humans and plants, whereas degraded by the bacterial systems. To date, no one has analyzed the complete degradation and poisonous effects of the atrazine in the environment. Therefore, this study is useful for overall system biology based modeling and simulation analysis of atrazine in living systems.Communicated by Ramaswamy H. Sarma.

Insights into zinc-sensing metalloregulator 'Zur' deciphering mechanism of zinc transportation in Bacillus spp. by modeling, simulation and molecular docking

To comprehend the molecular mechanism of zinc transportation by bacteria tends to be a very complicated and time-consuming method. To date, fragmented and scanty studies are available about the mechanism of zinc transportation at molecular level. So, the present study scrutinizes in silico pathways of zinc fractions transportation, specifically in Bacillus spp. stimulating dynamic performance of zinc. For this, the constructed model reveals Zur to be the prime regulatory transport protein maintaining bacterial survivability at fluctuation in zinc concentrations, thereby attaining zinc homeostasis. Topology for hub nodes displays appropriate evidence of the molecular basis of bacterial zinc imports and exports. Further, the molecular docking reveals interaction of Zur protein with the zinc ligands (ZnCO₃ and ZnSO₄). By validation of binding affinity, binding energy and docking score via Autodock Vina and X-Score, the ZnSO₄ compound was found to possess excellent stability in the active pocket site of Zur, stating Zur-ZnSO₄ complex to be the most potential. Owing to which, the Zur-ZnSO₄ complex was selected and subjected to molecular dynamics simulation, revealing RMSD, RG, RMSF, SASA and interaction energy for 20 ns trajectory period. Henceforth,the study provides novel insight into revealing the unrecognized Zur protein pathway, assisting zinc transportation, besides retaining best interaction with ZnSO₄ ligand. This is the first system biology where molecular docking and molecular dynamics simulation-based investigation decipher the role of Zur transport protein system and interaction of its amino acids with zinc ligands in a simpler and economical form via in silico techniques.Communicated by Ramaswamy H. Sarma.

SABIO-RK: an updated resource for manually curated biochemical reaction kinetics

SABIO-RK (http://sabiork.h-its.org/) is a manually curated database containing data about biochemical reactions and their reaction kinetics. The data are primarily extracted from scientific literature and stored in a relational database. The content comprises both naturally occurring and alternatively measured biochemical reactions and is not restricted to any organism class. The data are made available to the public by a web-based search interface and by web services for programmatic access. In this update we describe major improvements and extensions of SABIO-RK since our last publication in the database issue of Nucleic Acid Research (2012). (i) The website has been completely revised and (ii) allows now also free text search for kinetics data. (iii) Additional interlinkages with other databases in our field have been established; this enables users to gain directly comprehensive knowledge about the properties of enzymes and kinetics beyond SABIO-RK. (iv) Vice versa, direct access to SABIO-RK data has been implemented in several systems biology tools and workflows. (v) On request of our experimental users, the data can be exported now additionally in spreadsheet formats. (vi) The newly established SABIO-RK Curation Service allows to respond to specific data requirements.

Insights into Dynamic Network States Using Metabolomic Data

Metabolomic data is the youngest of the high-throughput data types; however, it is potentially one of the most informative, as it provides a direct, quantitative biochemical phenotype. There are a number of ways in which metabolomic data can be analyzed in systems biology; however, the thermodynamic and kinetic relevance of these data cannot be overstated. Genome-scale metabolic network reconstructions provide a natural context to incorporate metabolomic data in order to provide insight into the condition-specific kinetic characteristics of metabolic networks. Herein we discuss how metabolomic data can be incorporated into constraint-based models in a flexible framework that enables scaling from small pathways to cell-scale models, while being able to accommodate coarse-grained to more detailed, allosteric interactions, all using the well-known principle of mass action.

DMPy: a Python package for automated mathematical model construction of large-scale metabolic systems

Background

Models of metabolism are often used in biotechnology and pharmaceutical research to identify drug targets or increase the direct production of valuable compounds. Due to the complexity of large metabolic systems, a number of conclusions have been drawn using mathematical methods with simplifying assumptions. For example, constraint-based models describe changes of internal concentrations that occur much quicker than alterations in cell physiology. Thus, metabolite concentrations and reaction fluxes are fixed to constant values. This greatly reduces the mathematical complexity, while providing a reasonably good description of the system in steady state. However, without a large number of constraints, many different flux sets can describe the optimal model and we obtain no information on how metabolite levels dynamically change. Thus, to accurately determine what is taking place within the cell, finer quality data and more detailed models need to be constructed.

Results

In this paper we present a computational framework, DMPy, that uses a network scheme as input to automatically search for kinetic rates and produce a mathematical model that describes temporal changes of metabolite fluxes. The parameter search utilises several online databases to find measured reaction parameters. From this, we take advantage of previous modelling efforts, such as Parameter Balancing, to produce an initial mathematical model of a metabolic pathway. We analyse the effect of parameter uncertainty on model dynamics and test how recent flux-based model reduction techniques alter system properties. To our knowledge this is the first time such analysis has been performed on large models of metabolism. Our results highlight that good estimates of at least 80% of the reaction rates are required to accurately model metabolic systems. Furthermore, reducing the size of the model by grouping reactions together based on fluxes alters the resulting system dynamics.

Conclusion

The presented pipeline automates the modelling process for large metabolic networks. From this, users can simulate their pathway of interest and obtain a better understanding of how altering conditions influences cellular dynamics. By testing the effects of different parameterisations we are also able to provide suggestions to help construct more accurate models of complete metabolic systems in the future.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0584-8) contains supplementary material, which is available to authorized users.

Modeling of the jasmonate signaling pathway in Arabidopsis thaliana with respect to pathophysiology of Alternaria blight in Brassica

The productivity of Oilseed Brassica, one of the economically important crops of India, is seriously affected by the disease, Alternaria blight. The disease is mainly caused by two major necrotrophic fungi, Alternaria brassicae and Alternaria brassicicola which are responsible for significant yield losses. Till date, no resistant source is available against Alternaria blight, hence plant breeding methods can not be used to develop disease resistant varieties. Jasmonate mediated signalling pathway, which is known to play crucial role during defense response against necrotrophs, could be strengthened in Brassica plants to combat the disease. Since scanty information is available in Brassica-Alternaria pathosystems at molecular level therefore, in the present study efforts have been made to model jasmonic acid pathway in Arabidopsis thaliana to simulate the dynamic behaviour of molecular species in the model. Besides, the developed model was also analyzed topologically for investigation of the hubs node. COI1 is identified as one of the promising candidate genes in response to Alternaria and other linked components of plant defense mechanisms against the pathogens. The findings from present study are therefore informative for understanding the molecular basis of pathophysiology and rational management of Alternaria blight for securing food and nutritional security.

From word models to executable models of signaling networks using automated assembly

Word models (natural language descriptions of molecular mechanisms) are a common currency in spoken and written communication in biomedicine but are of limited use in predicting the behavior of complex biological networks. We present an approach to building computational models directly from natural language using automated assembly. Molecular mechanisms described in simple English are read by natural language processing algorithms, converted into an intermediate representation, and assembled into executable or network models. We have implemented this approach in the Integrated Network and Dynamical Reasoning Assembler (INDRA), which draws on existing natural language processing systems as well as pathway information in Pathway Commons and other online resources. We demonstrate the use of INDRA and natural language to model three biological processes of increasing scope: (i) p53 dynamics in response to DNA damage, (ii) adaptive drug resistance in BRAF‐V600E‐mutant melanomas, and (iii) the RAS signaling pathway. The use of natural language makes the task of developing a model more efficient and it increases model transparency, thereby promoting collaboration with the broader biology community.

MAMMOTh: A new database for curated mathematical models of biomolecular systems

MOTIVATION

Living systems have a complex hierarchical organization that can be viewed as a set of dynamically interacting subsystems. Thus, to simulate the internal nature and dynamics of the entire biological system, we should use the iterative way for a model reconstruction, which is a consistent composition and combination of its elementary subsystems. In accordance with this bottom-up approach, we have developed the MAthematical Models of bioMOlecular sysTems (MAMMOTh) tool that consists of the database containing manually curated MAMMOTh fitted to the experimental data and a software tool that provides their further integration.

RESULTS

The MAMMOTh database entries are organized as building blocks in a way that the model parts can be used in different combinations to describe systems with higher organizational level (metabolic pathways and/or transcription regulatory networks). The tool supports export of a single model or their combinations in SBML or Mathematica standards. The database currently contains 110 mathematical sub-models for Escherichia coli elementary subsystems (enzymatic reactions and gene expression regulatory processes) that can be combined in at least 5100 complex/sophisticated models concerning more complex biological processes as de novo nucleotide biosynthesis, aerobic/anaerobic respiration and nitrate/nitrite utilization in E. coli. All models are functionally interconnected and sufficiently complement public model resources.

AVAILABILITY

http://mammoth.biomodelsgroup.ru.

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Mechanistic mathematical modeling of biochemical reaction networks using ordinary differential equation (ODE) models has improved our understanding of small- and medium-scale biological processes. While the same should in principle hold for large- and genome-scale processes, the computational methods for the analysis of ODE models which describe hundreds or thousands of biochemical species and reactions are missing so far. While individual simulations are feasible, the inference of the model parameters from experimental data is computationally too intensive. In this manuscript, we evaluate adjoint sensitivity analysis for parameter estimation in large scale biochemical reaction networks. We present the approach for time-discrete measurement and compare it to state-of-the-art methods used in systems and computational biology. Our comparison reveals a significantly improved computational efficiency and a superior scalability of adjoint sensitivity analysis. The computational complexity is effectively independent of the number of parameters, enabling the analysis of large- and genome-scale models. Our study of a comprehensive kinetic model of ErbB signaling shows that parameter estimation using adjoint sensitivity analysis requires a fraction of the computation time of established methods. The proposed method will facilitate mechanistic modeling of genome-scale cellular processes, as required in the age of omics.

In this manuscript, we introduce a scalable method for parameter estimation for genome-scale biochemical reaction networks. Mechanistic models for genome-scale biochemical reaction networks describe the behavior of thousands of chemical species using thousands of parameters. Standard methods for parameter estimation are usually computationally intractable at these scales. Adjoint sensitivity based approaches have been suggested to have superior scalability but any rigorous evaluation is lacking. We implement a toolbox for adjoint sensitivity analysis for biochemical reaction network which also supports the import of SBML models. We show by means of a set of benchmark models that adjoint sensitivity based approaches unequivocally outperform standard approaches for large-scale models and that the achieved speedup increases with respect to both the number of parameters and the number of chemical species in the model. This demonstrates the applicability of adjoint sensitivity based approaches to parameter estimation for genome-scale mechanistic model. The MATLAB toolbox implementing the developed methods is available from http://ICB-DCM.github.io/AMICI/.

Kinetic modeling of cell metabolism for microbial production

Kinetic models of cellular metabolism are important tools for the rational design of metabolic engineering strategies and to explain properties of complex biological systems. The recent developments in high-throughput experimental data are leading to new computational approaches for building kinetic models of metabolism. Herein, we briefly survey the available databases, standards and software tools that can be applied for kinetic models of metabolism. In addition, we give an overview about recently developed ordinary differential equations (ODE)-based kinetic models of metabolism and some of the main applications of such models are illustrated in guiding metabolic engineering design. Finally, we review the kinetic modeling approaches of large-scale networks that are emerging, discussing their main advantages, challenges and limitations.