14,542 research outputs found
Synthesizing and tuning chemical reaction networks with specified behaviours
We consider how to generate chemical reaction networks (CRNs) from functional
specifications. We propose a two-stage approach that combines synthesis by
satisfiability modulo theories and Markov chain Monte Carlo based optimisation.
First, we identify candidate CRNs that have the possibility to produce correct
computations for a given finite set of inputs. We then optimise the reaction
rates of each CRN using a combination of stochastic search techniques applied
to the chemical master equation, simultaneously improving the of correct
behaviour and ruling out spurious solutions. In addition, we use techniques
from continuous time Markov chain theory to study the expected termination time
for each CRN. We illustrate our approach by identifying CRNs for majority
decision-making and division computation, which includes the identification of
both known and unknown networks.Comment: 17 pages, 6 figures, appeared the proceedings of the 21st conference
on DNA Computing and Molecular Programming, 201
Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
Stochasticity is a key characteristic of intracellular processes such as gene
regulation and chemical signalling. Therefore, characterising stochastic
effects in biochemical systems is essential to understand the complex dynamics
of living things. Mathematical idealisations of biochemically reacting systems
must be able to capture stochastic phenomena. While robust theory exists to
describe such stochastic models, the computational challenges in exploring
these models can be a significant burden in practice since realistic models are
analytically intractable. Determining the expected behaviour and variability of
a stochastic biochemical reaction network requires many probabilistic
simulations of its evolution. Using a biochemical reaction network model to
assist in the interpretation of time course data from a biological experiment
is an even greater challenge due to the intractability of the likelihood
function for determining observation probabilities. These computational
challenges have been subjects of active research for over four decades. In this
review, we present an accessible discussion of the major historical
developments and state-of-the-art computational techniques relevant to
simulation and inference problems for stochastic biochemical reaction network
models. Detailed algorithms for particularly important methods are described
and complemented with MATLAB implementations. As a result, this review provides
a practical and accessible introduction to computational methods for stochastic
models within the life sciences community
Formal executable descriptions of biological systems
The similarities between systems of living entities and systems of concurrent processes may support biological experiments in silico. Process calculi offer a formal framework to describe biological systems, as well as to analyse their behaviour, both from a qualitative and a quantitative point of view. A couple of little examples help us in showing how this can be done. We mainly focus our attention on the qualitative and quantitative aspects of the considered biological systems, and briefly illustrate which kinds of analysis are possible. We use a known stochastic calculus for the first example. We then present some statistics collected by repeatedly running the specification, that turn out to agree with those obtained by experiments in vivo. Our second example motivates a richer calculus. Its stochastic extension requires a non trivial machinery to faithfully reflect the real dynamic behaviour of biological systems
Landauer's principle as a special case of Galois connection
It is demonstrated how to construct a Galois connection between two related
systems with entropy. The construction, called the Landauer's connection,
describes coupling between two systems with entropy. It is straightforward and
transfers changes in one system to the other one preserving ordering structure
induced by entropy. The Landauer's connection simplifies the description of the
classical Landauer's principle for computational systems. Categorification and
generalization of the Landauer's principle opens area of modelling of various
systems in presence of entropy in abstract terms.Comment: 24 pages, 3 figure
Gene regulatory networks: a coarse-grained, equation-free approach to multiscale computation
We present computer-assisted methods for analyzing stochastic models of gene
regulatory networks. The main idea that underlies this equation-free analysis
is the design and execution of appropriately-initialized short bursts of
stochastic simulations; the results of these are processed to estimate
coarse-grained quantities of interest, such as mesoscopic transport
coefficients. In particular, using a simple model of a genetic toggle switch,
we illustrate the computation of an effective free energy and of a
state-dependent effective diffusion coefficient that characterize an
unavailable effective Fokker-Planck equation. Additionally we illustrate the
linking of equation-free techniques with continuation methods for performing a
form of stochastic "bifurcation analysis"; estimation of mean switching times
in the case of a bistable switch is also implemented in this equation-free
context. The accuracy of our methods is tested by direct comparison with
long-time stochastic simulations. This type of equation-free analysis appears
to be a promising approach to computing features of the long-time,
coarse-grained behavior of certain classes of complex stochastic models of gene
regulatory networks, circumventing the need for long Monte Carlo simulations.Comment: 33 pages, submitted to The Journal of Chemical Physic
Variable-free exploration of stochastic models: a gene regulatory network example
Finding coarse-grained, low-dimensional descriptions is an important task in
the analysis of complex, stochastic models of gene regulatory networks. This
task involves (a) identifying observables that best describe the state of these
complex systems and (b) characterizing the dynamics of the observables. In a
previous paper [13], we assumed that good observables were known a priori, and
presented an equation-free approach to approximate coarse-grained quantities
(i.e, effective drift and diffusion coefficients) that characterize the
long-time behavior of the observables. Here we use diffusion maps [9] to
extract appropriate observables ("reduction coordinates") in an automated
fashion; these involve the leading eigenvectors of a weighted Laplacian on a
graph constructed from network simulation data. We present lifting and
restriction procedures for translating between physical variables and these
data-based observables. These procedures allow us to perform equation-free
coarse-grained, computations characterizing the long-term dynamics through the
design and processing of short bursts of stochastic simulation initialized at
appropriate values of the data-based observables.Comment: 26 pages, 9 figure
The Energetic Costs of Cellular Computation
Cells often perform computations in response to environmental cues. A simple
example is the classic problem, first considered by Berg and Purcell, of
determining the concentration of a chemical ligand in the surrounding media. On
general theoretical grounds (Landuer's Principle), it is expected that such
computations require cells to consume energy. Here, we explicitly calculate the
energetic costs of computing ligand concentration for a simple two-component
cellular network that implements a noisy version of the Berg-Purcell strategy.
We show that learning about external concentrations necessitates the breaking
of detailed balance and consumption of energy, with greater learning requiring
more energy. Our calculations suggest that the energetic costs of cellular
computation may be an important constraint on networks designed to function in
resource poor environments such as the spore germination networks of bacteria.Comment: 9 Pages (including Appendix); 4 Figures; v3 corrects even more typo
Analysis of signalling pathways using the prism model checker
We describe a new modelling and analysis approach for signal
transduction networks in the presence of incomplete data. We illustrate
the approach with an example, the RKIP inhibited ERK pathway
[1]. Our models are based on high level descriptions of continuous time
Markov chains: reactions are modelled as synchronous processes and concentrations
are modelled by discrete, abstract quantities. The main advantage
of our approach is that using a (continuous time) stochastic logic
and the PRISM model checker, we can perform quantitative analysis of
queries such as if a concentration reaches a certain level, will it remain at
that level thereafter? We also perform standard simulations and compare
our results with a traditional ordinary differential equation model. An
interesting result is that for the example pathway, only a small number
of discrete data values is required to render the simulations practically
indistinguishable
- …