Design and Development of Software Tools for Bio-PEPA

This paper surveys the design of software tools for the Bio-PEPA process algebra. Bio-PEPA is a high-level language for modelling biological systems such as metabolic pathways and other biochemical reaction networks. Through providing tools for this modelling language we hope to allow easier use of a range of simulators and model-checkers thereby freeing the modeller from the responsibility of developing a custom simulator for the problem of interest. Further, by providing mappings to a range of different analysis tools the Bio-PEPA language allows modellers to compare analysis results which have been computed using independent numerical analysers, which enhances the reliability and robustness of the results computed.

Process Calculi Abstractions for Biology

Several approaches have been proposed to model biological systems by means of the formal techniques and tools available in computer science. To mention just a few of them, some representations are inspired by Petri Nets theory, and some other by stochastic processes. A most recent approach consists in interpreting the living entities as terms of process calculi where the behavior of the represented systems can be inferred by applying syntax-driven rules. A comprehensive picture of the state of the art of the process calculi approach to biological modeling is still missing. This paper goes in the direction of providing such a picture by presenting a comparative survey of the process calculi that have been used and proposed to describe the behavior of living entities. This is the preliminary version of a paper that was published in Algorithmic Bioprocesses. The original publication is available at http://www.springer.com/computer/foundations/book/978-3-540-88868-

Simulation of non-Markovian Processes in BlenX

BlenX is a programming language explicitly designed for modeling biological processes inspired by Beta-binders. The actual framework assumes biochemical interactions being exponentially distributed, i.e., an underlying Markov process is associated with BlenX programs. In this paper we relax this condition by providing formal tools for managing non-Markovian processes within BlenX

Curve classes on irreducible holomorphic symplectic varieties

We prove that the integral Hodge conjecture holds for 1-cycles on irreducible holomorphic symplectic varieties of K3 type and of Generalized Kummer type. As an application, we give a new proof of the integral Hodge conjecture for cubic fourfolds.Comment: 15 page

Narrative-based computational modelling of the Gp130/JAK/STAT signalling pathway.

BACKGROUND: Appropriately formulated quantitative computational models can support researchers in understanding the dynamic behaviour of biological pathways and support hypothesis formulation and selection by "in silico" experimentation. An obstacle to widespread adoption of this approach is the requirement to formulate a biological pathway as machine executable computer code. We have recently proposed a novel, biologically intuitive, narrative-style modelling language for biologists to formulate the pathway which is then automatically translated into an executable format and is, thus, usable for analysis via existing simulation techniques. RESULTS: Here we use a high-level narrative language in designing a computational model of the gp130/JAK/STAT signalling pathway and show that the model reproduces the dynamic behaviour of the pathway derived by biological observation. We then "experiment" on the model by simulation and sensitivity analysis to define those parameters which dominate the dynamic behaviour of the pathway. The model predicts that nuclear compartmentalisation and phosphorylation status of STAT are key determinants of the pathway and that alternative mechanisms of signal attenuation exert their influence on different timescales. CONCLUSION: The described narrative model of the gp130/JAK/STAT pathway represents an interesting case study showing how, by using this approach, researchers can model biological systems without explicitly dealing with formal notations and mathematical expressions (typically used for biochemical modelling), nevertheless being able to obtain simulation and analysis results. We present the model and the sensitivity analysis results we have obtained, that allow us to identify the parameters which are most sensitive to perturbations. The results, which are shown to be in agreement with existing mathematical models of the gp130/JAK/STAT pathway, serve us as a form of validation of the model and of the approach itself

Snazer: the simulations and networks analyzer

<p>Abstract</p> <p>Background</p> <p>Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.</p> <p>From a pure dynamical perspective, simulation represents a relevant <it>way</it>-<it>out</it>. Many simulator tools capitalized on the "noisy" behavior of some systems and used formal models to represent cellular activities as temporal trajectories. Statistical methods have been applied to a fairly large number of replicated trajectories in order to infer knowledge.</p> <p>A tool which both graphically manipulates reactive models and deals with sets of simulation time-course data by aggregation, interpretation and statistical analysis is missing and could add value to simulators.</p> <p>Results</p> <p>We designed and implemented <it>Snazer</it>, the simulations and networks analyzer. Its goal is to aid the processes of visualizing and manipulating reactive models, as well as to share and interpret time-course data produced by stochastic simulators or by any other means.</p> <p>Conclusions</p> <p><it>Snazer </it>is a solid prototype that integrates biological network and simulation time-course data analysis techniques.</p

Verification of Spatial and Temporal Modalities in Biochemical Systems

AbstractBiochemical systems such as metabolic and signaling pathways tend to be arranged in a physical space: the product of one reaction must be in the right place to become the reactant for the subsequent reaction in the pathway. Moreover, in some cases, the behavior of the systems can depend on both, the location of the reactants as well as on the time needed for the reaction to occur. We address the problem of specifying and verifying properties of biochemical systems that exhibit both temporal and spatial modalities at the same time. For that, we use as specification language a fragment of intuitionistic linear logic with subexponentials (SELL). The subexponential signature allows us to capture the spatial relations among the different components of the system and the timed constraints for reactions to occur. We show that our framework is general enough to give a declarative semantics to P-Systems and we show that such logical characterization has a strong level of adequacy. Hence, derivations in SELL follow exactly the behavior of the modeled system

A Stochastic Pi Calculus for Concurrent Objects

International audienceWe present SpiCO, a new modeling and simulation language for system biology, based on the stochastic pi-calculus. SpiCO supports higher level modeling via multi-profile concurrent objects with static inheritance. We present a semantics for SpiCO in terms of continuous time Markov chains, and show how to compile SpiCO back into the biochemical stochastic pi-calculus while preserving semantics

Investigating modularity in the analysis of process algebra models of biochemical systems

Compositionality is a key feature of process algebras which is often cited as one of their advantages as a modelling technique. It is certainly true that in biochemical systems, as in many other systems, model construction is made easier in a formalism which allows the problem to be tackled compositionally. In this paper we consider the extent to which the compositional structure which is inherent in process algebra models of biochemical systems can be exploited during model solution. In essence this means using the compositional structure to guide decomposed solution and analysis. Unfortunately the dynamic behaviour of biochemical systems exhibits strong interdependencies between the components of the model making decomposed solution a difficult task. Nevertheless we believe that if such decomposition based on process algebras could be established it would demonstrate substantial benefits for systems biology modelling. In this paper we present our preliminary investigations based on a case study of the pheromone pathway in yeast, modelling in the stochastic process algebra Bio-PEPA

Probabilistic model checking of complex biological pathways

Probabilistic model checking is a formal verification technique that has been successfully applied to the analysis of systems from a broad range of domains, including security and communication protocols, distributed algorithms and power management. In this paper we illustrate its applicability to a complex biological system: the FGF (Fibroblast Growth Factor) signalling pathway. We give a detailed description of how this case study can be modelled in the probabilistic model checker PRISM, discussing some of the issues that arise in doing so, and show how we can thus examine a rich selection of quantitative properties of this model. We present experimental results for the case study under several different scenarios and provide a detailed analysis, illustrating how this approach can be used to yield a better understanding of the dynamics of the pathway