Stochastic Simulation of Process Calculi for Biology

Biological systems typically involve large numbers of components with complex, highly parallel interactions and intrinsic stochasticity. To model this complexity, numerous programming languages based on process calculi have been developed, many of which are expressive enough to generate unbounded numbers of molecular species and reactions. As a result of this expressiveness, such calculi cannot rely on standard reaction-based simulation methods, which require fixed numbers of species and reactions. Rather than implementing custom stochastic simulation algorithms for each process calculus, we propose to use a generic abstract machine that can be instantiated to a range of process calculi and a range of reaction-based simulation algorithms. The abstract machine functions as a just-in-time compiler, which dynamically updates the set of possible reactions and chooses the next reaction in an iterative cycle. In this short paper we give a brief summary of the generic abstract machine, and show how it can be instantiated with the stochastic simulation algorithm known as Gillespie's Direct Method. We also discuss the wider implications of such an abstract machine, and outline how it can be used to simulate multiple calculi simultaneously within a common framework.Comment: In Proceedings MeCBIC 2010, arXiv:1011.005

Abstract Interpretation for Probabilistic Termination of Biological Systems

In a previous paper the authors applied the Abstract Interpretation approach for approximating the probabilistic semantics of biological systems, modeled specifically using the Chemical Ground Form calculus. The methodology is based on the idea of representing a set of experiments, which differ only for the initial concentrations, by abstracting the multiplicity of reagents present in a solution, using intervals. In this paper, we refine the approach in order to address probabilistic termination properties. More in details, we introduce a refinement of the abstract LTS semantics and we abstract the probabilistic semantics using a variant of Interval Markov Chains. The abstract probabilistic model safely approximates a set of concrete experiments and reports conservative lower and upper bounds for probabilistic termination

External Control in Process Algebra for Systems Biology

A critical aspect in the modeling of biological systems is the description view point. On the one hand, the Stochastic \u3c0-calculus formalism provides an intuitive and compact representation from an internal perspective. On the other hand, other proposed languages such as Hybrid Automata and Stochastic Concurrent Constraint Programming introduce in the system description an external control and provide more structured models. This work aims at bridging the above discussed gap. In particular, we propose a different approach for the encoding of biological systems in Stochastic \u3c0-calculus in the direction of introducing an external control and comparing different formalisms. We show the effectiveness of our method on some examples

Theoretical basis of the community effect in development

Peer reviewedPublisher PD

Multi-Level Kinetic Model of mRNA Delivery via Transfection of Lipoplexes

Recent work on the use of mRNA lipoplexes for gene delivery demonstrates the need for a mathematical model that simulates and predicts kinetics and transfection efficiency. The small copy numbers involved make it necessary to use stochastic models and include statistical analysis of the variation observed in the experimental data. The modeling requirements are further complicated by the multi-level nature of the problem, where mRNA molecules are contained in lipoplexes, which are in turn contained in endosomes, where each of these entities displays a behavior of its own. We have created a mathematical model that reproduces both the time courses and the statistical variance observed in recent experiments using single-cell tracking of GFP expression after transfection. By applying a few key simplifications and assumptions, we have limited the number of free parameters to five, which we optimize to match five experimental determinants by means of a simulated annealing algorithm. The models demonstrate the need for modeling of nested species in order to reproduce the shape of the dose-response and expression-level curves

Multi-Level Kinetic Model of mRNA Delivery via Transfection of Lipoplexes

Recent work on the use of mRNA lipoplexes for gene delivery demonstrates the need for a mathematical model that simulates and predicts kinetics and transfection efficiency. The small copy numbers involved make it necessary to use stochastic models and include statistical analysis of the variation observed in the experimental data. The modeling requirements are further complicated by the multi-level nature of the problem, where mRNA molecules are contained in lipoplexes, which are in turn contained in endosomes, where each of these entities displays a behavior of its own. We have created a mathematical model that reproduces both the time courses and the statistical variance observed in recent experiments using single-cell tracking of GFP expression after transfection. By applying a few key simplifications and assumptions, we have limited the number of free parameters to five, which we optimize to match five experimental determinants by means of a simulated annealing algorithm. The models demonstrate the need for modeling of nested species in order to reproduce the shape of the dose-response and expression-level curves

Stochastic models for the in silico simulation of synaptic processes

Background: Research in life sciences is benefiting from a large availability of formal description techniques and analysis methodologies. These allow both the phenomena investigated to be precisely modeled and virtual experiments to be performed in silico. Such experiments may result in easier, faster, and satisfying approximations of their in vitro/vivo counterparts. A promising approach is represented by the study of biological phenomena as a collection of interactive entities through process calculi equipped with stochastic semantics. These exploit formal grounds developed in the theory of concurrency in computer science, account for the not continuous, nor discrete, nature of many phenomena, enjoy nice compositional properties and allow for simulations that have been demonstrated to be coherent with data in literature. Results: Motivated by the need to address some aspects of the functioning of neural synapses, we have developed one such model for synaptic processes in the calyx of Held, which is a glutamatergic synapse in the auditory pathway of the mammalia. We have developed such a stochastic model starting from existing kinetic models based on ODEs of some sub-components of the synapse, integrating other data from literature and making some assumptions about non-fully understood processes. Experiments have confirmed the coherence of our model with known biological data, also validating the assumptions made. Our model overcomes some limitations of the kinetic ones and, to our knowledge, represents the first model of synaptic processes based on process calculi. The compositionality of the approach has permitted us to independently focus on tuning the models of the pre- and post- synaptic traits, and then to naturally connect them, by dealing with “interface” issues. Furthermore, we have improved the expressiveness of the model, e.g. by embedding easy control of element concentration time courses. Sensitivity analysis over several parameters of the model has provided results that may help clarify the dynamics of synaptic transmission, while experiments with the model of the complete synapse seem worth explaining short-term plasticity mechanisms. Conclusions: Specific presynaptic and postsynaptic mechanisms can be further analysed under various conditions, for instance by studying the presynaptic behaviour under repeated activations. The level of details of the description can be refined, for instance by further specifying the neurotransmitter generation and release steps. Taking advantage of the compositionality of the approach, an enhanced model could then be composed with other neural models, designed within the same framework, in order to obtain a more detailed and comprehensive model. In the long term, we are interested, in particular, in addressing models of synaptic plasticity, i.e. activity dependent mechanisms, which are the bases of memory and learning processes. More on the computer science side, we plan to follow some directions to improve the underlying computational model and the linguistic primitives it provides as suggested by the experiments carried out, e.g. by introducing a suitable notion of (spatial) locality

On Designing Multicore-Aware Simulators for Systems Biology Endowed with OnLine Statistics

The paper arguments are on enabling methodologies for the design of a fully parallel, online, interactive tool aiming to support the bioinformatics scientists .In particular, the features of these methodologies, supported by the FastFlow parallel programming framework, are shown on a simulation tool to perform the modeling, the tuning, and the sensitivity analysis of stochastic biological models. A stochastic simulation needs thousands of independent simulation trajectories turning into big data that should be analysed by statistic and data mining tools. In the considered approach the two stages are pipelined in such a way that the simulation stage streams out the partial results of all simulation trajectories to the analysis stage that immediately produces a partial result. The simulation-analysis workflow is validated for performance and effectiveness of the online analysis in capturing biological systems behavior on a multicore platform and representative proof-of-concept biological systems. The exploited methodologies include pattern-based parallel programming and data streaming that provide key features to the software designers such as performance portability and efficient in-memory (big) data management and movement. Two paradigmatic classes of biological systems exhibiting multistable and oscillatory behavior are used as a testbed