ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra

Background: Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agent-based models, with the goal to gain a better understanding of the system. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. Although there exist sophisticated algorithms to determine the dynamics of discrete models, their implementations usually require labor-intensive formatting of the model formulation, and they are oftentimes not accessible to users without programming skills. Efficient analysis methods are needed that are accessible to modelers and easy to use. Method: By converting discrete models into algebraic models, tools from computational algebra can be used to analyze their dynamics. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a long-studied problem in computer algebra. Results: A method for efficiently identifying attractors, and the web-based tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness, i.e., while the number of nodes in a biological network may be quite large, each node is affected only by a small number of other nodes, and robustness, i.e., small number of attractors

esyN: network building, sharing and publishing.

The construction and analysis of networks is increasingly widespread in biological research. We have developed esyN ("easy networks") as a free and open source tool to facilitate the exchange of biological network models between researchers. esyN acts as a searchable database of user-created networks from any field. We have developed a simple companion web tool that enables users to view and edit networks using data from publicly available databases. Both normal interaction networks (graphs) and Petri nets can be created. In addition to its basic tools, esyN contains a number of logical templates that can be used to create models more easily. The ability to use previously published models as building blocks makes esyN a powerful tool for the construction of models and network graphs. Users are able to save their own projects online and share them either publicly or with a list of collaborators. The latter can be given the ability to edit the network themselves, allowing online collaboration on network construction. esyN is designed to facilitate unrestricted exchange of this increasingly important type of biological information. Ultimately, the aim of esyN is to bring the advantages of Open Source software development to the construction of biological networks.This is the final published version of the paper. It's also available from PLOS at http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0106035

Reconstruction of extended Petri nets from time series data and its application to signal transduction and to gene regulatory networks

<p>Abstract</p> <p>Background</p> <p>Network inference methods reconstruct mathematical models of molecular or genetic networks directly from experimental data sets. We have previously reported a mathematical method which is exclusively data-driven, does not involve any heuristic decisions within the reconstruction process, and deliveres all possible alternative minimal networks in terms of simple place/transition Petri nets that are consistent with a given discrete time series data set.</p> <p>Results</p> <p>We fundamentally extended the previously published algorithm to consider catalysis and inhibition of the reactions that occur in the underlying network. The results of the reconstruction algorithm are encoded in the form of an extended Petri net involving control arcs. This allows the consideration of processes involving mass flow and/or regulatory interactions. As a non-trivial test case, the phosphate regulatory network of enterobacteria was reconstructed using <it>in silico</it>-generated time-series data sets on wild-type and <it>in silico </it>mutants.</p> <p>Conclusions</p> <p>The new exact algorithm reconstructs extended Petri nets from time series data sets by finding all alternative minimal networks that are consistent with the data. It suggested alternative molecular mechanisms for certain reactions in the network. The algorithm is useful to combine data from wild-type and mutant cells and may potentially integrate physiological, biochemical, pharmacological, and genetic data in the form of a single model.</p

ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra

Abstract Background Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agent-based models, to gain a better understanding of them. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. There exist software tools to analyze discrete models, but they either lack the algorithmic functionality to analyze complex models deterministically or they are inaccessible to many users as they require understanding the underlying algorithm and implementation, do not have a graphical user interface, or are hard to install. Efficient analysis methods that are accessible to modelers and easy to use are needed. Results We propose a method for efficiently identifying attractors and introduce the web-based tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a long-studied problem in computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness and robustness. For a large set of published complex discrete models, ADAM identified the attractors in less than one second. Conclusions Discrete modeling techniques are a useful tool for analyzing complex biological systems and there is a need in the biological community for accessible efficient analysis tools. ADAM provides analysis methods based on mathematical algorithms as a web-based tool for several different input formats, and it makes analysis of complex models accessible to a larger community, as it is platform independent as a web-service and does not require understanding of the underlying mathematics

Petri Net modelling approach for analysing the behaviour of Wnt/[inline-formula removed] -catenin and Wnt/Ca 2+ signalling pathways in arrhythmogenic right ventricular cardiomyopathy.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart muscle disease that may result in arrhythmia, heart failure and sudden death. The hallmark pathological findings are progressive myocyte loss and fibro fatty replacement, with a predilection for the right ventricle. This study focuses on the adipose tissue formation in cardiomyocyte by considering the signal transduction pathways including Wnt/[inline-formula removed]-catenin and Wnt/Ca2+ regulation system. These pathways are modelled and analysed using stochastic petri nets (SPN) in order to increase our comprehension of ARVC and in turn its treatment regimen. The Wnt/[inline-formula removed]-catenin model predicts that the dysregulation or absence of Wnt signalling, inhibition of dishevelled and elevation of glycogen synthase kinase 3 along with casein kinase I are key cytotoxic events resulting in apoptosis. Moreover, the Wnt/Ca2+ SPN model demonstrates that the Bcl2 gene inhibited by c-Jun N-terminal kinase protein in the event of endoplasmic reticulum stress due to action potential and increased amount of intracellular Ca2+ which recovers the Ca2+homeostasis by phospholipase C, this event positively regulates the Bcl2 to suppress the mitochondrial apoptosis which causes ARVC

Biomodelkit - a framework for modular biomodel-engineering

Otto-von-Guericke-Universität Magdeburg, Fakultät für Naturwissenschaften, Dissertation, 2017von Dipl.-Ing. Mary-Ann BlätkeLiteraturverzeichnis: Seite [177]-18

A Unique Transformation from Ordinary Differential Equations to Reaction Networks

Many models in Systems Biology are described as a system of Ordinary Differential Equations, which allows for transient, steady-state or bifurcation analysis when kinetic information is available. Complementary structure-related qualitative analysis techniques have become increasingly popular in recent years, like qualitative model checking or pathway analysis (elementary modes, invariants, flux balance analysis, graph-based analyses, chemical organization theory, etc.). They do not rely on kinetic information but require a well-defined structure as stochastic analysis techniques equally do. In this article, we look into the structure inference problem for a model described by a system of Ordinary Differential Equations and provide conditions for the uniqueness of its solution. We describe a method to extract a structured reaction network model, represented as a bipartite multigraph, for example, a continuous Petri net (CPN), from a system of Ordinary Differential Equations (ODEs). A CPN uniquely defines an ODE, and each ODE can be transformed into a CPN. However, it is not obvious under which conditions the transformation of an ODE into a CPN is unique, that is, when a given ODE defines exactly one CPN. We provide biochemically relevant sufficient conditions under which the derived structure is unique and counterexamples showing the necessity of each condition. Our method is implemented and available; we illustrate it on some signal transduction models from the BioModels database. A prototype implementation of the method is made available to modellers at http://contraintes.inria.fr/~soliman/ode2pn.html, and the data mentioned in the “Results” section at http://contraintes.inria.fr/~soliman/ode2pn_data/. Our results yield a new recommendation for the import/export feature of tools supporting the SBML exchange format

Engineering a novel self-powering electrochemical biosensor

This paper records the efforts of a multi-disciplinary team of undergraduate students from Glasgow University to collectively design and carry out a 10 week project in Synthetic Biology as part of the international Genetic Engineered Machine competition (iGEM). The aim of the project was to design and build a self-powering electrochemical biosensor called ‘ElectrEcoBlu’. The novelty of this engineered machine lies in coupling a biosensor with a microbial fuel cell to transduce a pollution input into an easily measurable electrical output signal. The device consists of two components; the sensor element which is modular, allowing for customisation to detect a range of input signals as required, and the universal reporter element which is responsible for generating an electrical signal as an output. The genetic components produce pyocyanin, a competitive electron mediator for microbial fuel cells, thus enabling the generation of an electrical current in the presence of target chemical pollutants. The pollutants tested in our implementation were toluene and salicylate. ElectrEcoBlu is expected to drive forward the development of a new generation of biosensors. Our approach exploited a range of state-of-the-art modelling techniques in a unified framework of qualitative, stochastic and continuous approaches to support the design and guide the construction of this novel biological machine. This work shows that integrating engineering techniques with scientific methodologies can provide new insights into genetic regulation and can be considered as a reference framework for the development of biochemical systems in synthetic biology

Modeling Biological Gradient Formation: Combining Partial Differential Equations and Petri Nets

Algorithms and the Foundations of Software technologyComputer Systems, Imagery and Medi