3,989 research outputs found
PID Control of Biochemical Reaction Networks
Principles of feedback control have been shown to naturally arise in
biological systems and successfully applied to build synthetic circuits. In
this work we consider Biochemical Reaction Networks (CRNs) as a paradigm for
modelling biochemical systems and provide the first implementation of a
derivative component in CRNs. That is, given an input signal represented by the
concentration level of some species, we build a CRN that produces as output the
concentration of two species whose difference is the derivative of the input
signal. By relying on this component, we present a CRN implementation of a
feedback control loop with Proportional-Integral-Derivative (PID) controller
and apply the resulting control architecture to regulate the protein expression
in a microRNA regulated gene expression model.Comment: 8 Pages, 4 figures, Submitted to CDC 201
PID control of biochemical reaction networks
Principles of feedback control have been shown to naturally arise in biological systems and have been applied with success to build synthetic circuits. Here, we present an implementation of a proportional–integral–derivative (PID) controller as a chemical reaction network with mass-action kinetics. This makes the controller synthesizable in vitro using DNA strand displacement technology, owing to its demonstrated capability of realizing arbitrary reaction-network designs as interacting DNA molecules. Previous related work has studied biological PID architectures using linearizations of nonlinear dynamics arising in both the controller components and in the plant. In this article, we present a proof of correctness of our nonlinear design in closed loop using arguments from singular perturbation theory. As an application to show the effectiveness of our controller, we provide numerical simulations on a genetic model to perform PID feedback control of protein expression
On Estimating Derivatives of Input Signals in Biochemistry
The online estimation of the derivative of an input signal is widespread in
control theory and engineering. In the realm of chemical reaction networks
(CRN), this raises however a number of specific issues on the different ways to
achieve it. A CRN pattern for implementing a derivative block has already been
proposed for the PID control of biochemical processes, and proved correct using
Tikhonov's limit theorem. In this paper, we give a detailed mathematical
analysis of that CRN, thus clarifying the computed quantity and quantifying the
error done as a function of the reaction kinetic parameters. In a synthetic
biology perspective, we show how this can be used to design error correcting
terms to compute online functions involving derivatives with CRNs. In the
systems biology perspective, we give the list of models in BioModels containing
(in the sense of subgraph epimorphisms) the core derivative CRN, most of which
being models of oscillators and control systems in the cell, and discuss in
detail two such examples: one model of the circadian clock and one model of a
bistable switch
Finite time distributions of stochastically modeled chemical systems with absolute concentration robustness
Recent research in both the experimental and mathematical communities has
focused on biochemical interaction systems that satisfy an "absolute
concentration robustness" (ACR) property. The ACR property was first discovered
experimentally when, in a number of different systems, the concentrations of
key system components at equilibrium were observed to be robust to the total
concentration levels of the system. Followup mathematical work focused on
deterministic models of biochemical systems and demonstrated how chemical
reaction network theory can be utilized to explain this robustness. Later
mathematical work focused on the behavior of this same class of reaction
networks, though under the assumption that the dynamics were stochastic. Under
the stochastic assumption, it was proven that the system will undergo an
extinction event with a probability of one so long as the system is
conservative, showing starkly different long-time behavior than in the
deterministic setting. Here we consider a general class of stochastic models
that intersects with the class of ACR systems studied previously. We consider a
specific system scaling over compact time intervals and prove that in a limit
of this scaling the distribution of the abundances of the ACR species converges
to a certain product-form Poisson distribution whose mean is the ACR value of
the deterministic model. This result is in agreement with recent conjectures
pertaining to the behavior of ACR networks endowed with stochastic kinetics,
and helps to resolve the conflicting theoretical results pertaining to
deterministic and stochastic models in this setting
The Two Regime method for optimizing stochastic reaction-diffusion simulations
The computer simulation of stochastic reaction-diffusion processes in biology is often done using either compartment-based (spatially discretized) simulations or molecular-based (Brownian dynamics) approaches. Compartment-based approaches can yield quick and accurate mesoscopic results but lack the level of detail that is characteristic of the more computationally intensive molecular-based models. Often microscopic detail is only required in a small region but currently the best way to achieve this detail is to use a resource intensive model over the whole domain. We introduce the Two Regime Method (TRM) in which a molecular-based algorithm is used in part of the computational domain and a compartment-based approach is used elsewhere in the computational domain. We apply the TRM to two test problems including a model from developmental biology. We thereby show that the TRM is accurate and subsequently may be used to inspect both mesoscopic and microscopic detail of reaction diffusion simulations according to the demands of the modeller
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
- …