Quasi-robust control of biochemical reaction networks via stochastic morphing

Implementation of an Ultrasensitive Biomolecular Controller for Enzymatic Reaction Processes With Delay Using DNA Strand Displacement

In this article, a set of abstract chemical reactions has been employed to construct a novel nonlinear biomolecular controller, i.e, the Brink controller (BC) with direct positive autoregulation (DPAR) (namely BC-DPAR controller). In comparison to dual rail representation-based controllers such as the quasi sliding mode (QSM) controller, the BC-DPAR controller directly reduces the number of CRNs required for realizing an ultrasensitive input-output response because it does not involve the subtraction module, reducing the complexity of DNA implementations. Then, the action mechanism and steady-state condition constraints of two nonlinear controllers, BC-DPAR controller and QSM controller, are investigated further. Considering the mapping relationship between CRNs and DNA implementation, a CRNs-based enzymatic reaction process with delay is constructed, and a DNA strand displacement (DSD) scheme representing time delay is proposed. The BC-DPAR controller, when compared to the QSM controller, can reduce the number of abstract chemical reactions and DSD reactions required by 33.3% and 31.8%, respectively. Finally, an enzymatic reaction scheme with BC-DPAR controller is designed using DSD reactions. According to the findings, the enzymatic reaction process's output substance can approach the target level at a quasi-steady state in both delay-free and non-zero delay conditions, but the target level can only be achieved during a finite-time period, mainly due to the fuel stand depletion.

Stochastic approximations of higher-molecular by bi-molecular reactions

Reactions involving three or more reactants, called higher-molecular reactions, play an important role in mathematical modelling in systems and synthetic biology. In particular, such reactions underpin a variety of important bio-dynamical phenomena, such as multi-stability/multi-modality, oscillations, bifurcations, and noise-induced effects. However, as opposed to reactions involving at most two reactants, called bi-molecular reactions, higher-molecular reactions are biochemically improbable. To bridge the gap, in this paper we put forward an algorithm for systematically approximating arbitrary higher-molecular reactions with bi-molecular ones, while preserving the underlying stochastic dynamics. Properties of the algorithm and convergence are established via singular perturbation theory. The algorithm is applied to a variety of higher-molecular biochemical networks, and is shown to play an important role in synthetic biology.

Transition graph decomposition for complex balanced reaction networks with non-mass-action kinetics

Reaction networks are widely used models to describe biochemical processes. Stochastic fluctuations in the counts of biological macromolecules have amplified consequences due to their small population sizes. This makes it necessary to favor stochastic, discrete population, continuous time models. The stationary distributions provide snapshots of the model behavior at the stationary regime, and as such finding their expression in terms of the model parameters is of great interest. The aim of the present paper is to describe when the stationary distributions of the original model, whose state space is potentially infinite, coincide exactly with the stationary distributions of the process truncated to finite subsets of states, up to a normalizing constant. The finite subsets of states we identify are called copies and are inspired by the modular topology of reaction network models. With such a choice we prove a novel graphical characterization of the concept of complex balancing for stochastic models of reaction networks. The results of the paper hold for the commonly used mass-action kinetics but are not restricted to it, and are in fact stated for more general setting.

Exploiting Noise, Non-Linearity, and Feedback for Differential Control of Multiple Synthetic Cells with a Single Optogenetic Input

Synthetic biology seeks to develop modular biocircuits that combine to produce complex, controllable behaviors. These designs are often subject to noisy fluctuations and uncertainties, and most modern synthetic biology design processes have focused to create robust components to mitigate the noise of gene expression and reduce the heterogeneity of single-cell responses. However, a deeper understanding of noise can achieve control goals that would otherwise be impossible. We explore how an "Optogenetic Maxwell Demon" could selectively amplify noise to control multiple cells using single-input-multiple-output (SIMO) feedback. Using data-constrained stochastic model simulations and theory, we show how an appropriately selected stochastic SIMO controller can drive multiple different cells to different user-specified configurations irrespective of initial conditions. We explore how controllability depends on cells' regulatory structures, the amount of information available to the controller, and the accuracy of the model used. Our results suggest that gene regulation noise, when combined with optogenetic feedback and non-linear biochemical auto-regulation, can achieve synergy to enable precise control of complex stochastic processes.