Stochastic focusing coupled with negative feedback enables robust regulation in biochemical reaction networks

Inherent regulatory asymmetry emanating from network architecture in a prevalent autoregulatory motif

Predicting gene expression from DNA sequence remains a major goal in the field of gene regulation. A challenge to this goal is the connectivity of the network, whose role in altering gene expression remains unclear. Here, we study a common autoregulatory network motif, the negative single-input module, to explore the regulatory properties inherited from the motif. Using stochastic simulations and a synthetic biology approach in E. coli, we find that the TF gene and its target genes have inherent asymmetry in regulation, even when their promoters are identical; the TF gene being more repressed than its targets. The magnitude of asymmetry depends on network features such as network size and TF-binding affinities. Intriguingly, asymmetry disappears when the growth rate is too fast or too slow and is most significant for typical growth conditions. These results highlight the importance of accounting for network architecture in quantitative models of gene expression.

Dynamic disorder in simple enzymatic reactions induces stochastic amplification of substrate

A growing amount of evidence over the last two decades points to the fact that many enzymes exhibit fluctuations in their catalytic activity, which are associated with conformational changes on a broad range of timescales. The experimental study of this phenomenon, termed dynamic disorder, has become possible thanks to advances in single-molecule enzymology measurement techniques, through which the catalytic activity of individual enzyme molecules can be tracked in time. The biological role and importance of these fluctuations in a system with a small number of enzymes, such as a living cell, have only recently started being explored. In this work, we examine a simple stochastic reaction system consisting of an inflowing substrate and an enzyme with a randomly fluctuating catalytic reaction rate that converts the substrate into an outflowing product. To describe analytically the effect of rate fluctuations on the average substrate abundance at steady state, we derive an explicit formula that connects the relative speed of enzymatic fluctuations with the mean substrate level. Under fairly general modelling assumptions, we demonstrate that the relative speed of rate fluctuations can have a dramatic effect on the mean substrate, and lead to large positive deviations from predictions based on the assumption of deterministic enzyme activity. Our results also establish an interesting connection between the amplification effect and the mixing properties of the Markov process describing the enzymatic activity fluctuations, which can be used to easily predict the fluctuation speed above which such deviations become negligible. As the techniques of single-molecule enzymology continuously evolve, it may soon be possible to study the stochastic phenomena due to enzymatic activity fluctuations within living cells. Our work can be used to formulate experimentally testable hypotheses regarding the nature and magnitude of these fluctuations, as well as their phenotypic consequences.

Perfect Adaptation and Optimal Equilibrium Productivity in a Simple Microbial Biofuel Metabolic Pathway Using Dynamic Integral Control

The production of complex biomolecules by genetically engineered organisms is one of the most promising applications of metabolic engineering and synthetic biology. To obtain processes with high productivity, it is therefore crucial to design and implement efficient dynamic in vivo regulation strategies. We consider here the microbial biofuel production model of Dunlop et al. (2010) for which we demonstrate that an antithetic dynamic integral control strategy can achieve robust perfect adaptation for the intracellular biofuel concentration in the presence of poorly known network parameters and implementation errors in certain rate parameters of the controller. We also show that the maximum equilibrium extracellular biofuel productivity is fully defined by some of the network parameters and, in this respect, it can only be achieved when all the corresponding parameters are perfectly known. Since this optimum is a network property, it cannot be improved by the use of any controller that measures the intracellular biofuel concentration and acts on the production of pump proteins. Additional intrinsic fundamental properties for the process are also unveiled, the most important ones being the existence of a conservation relation between the productivity and the toxicity, a low sensitivity of the optimal productivity with respect to a poor implementation of the set-point for the intracellular biofuel, and a strong intrinsic robustness property of the optimal productivity with respect to poorly known parameters. Taken together, these results demonstrate that a high and robust equilibrium rate of production for the extracellular biofuel can be achieved when the parameters of the model are poorly known and those of the controllers are poorly implemented. Finally, several advantages of the proposed dynamic strategy over a static one are also emphasized.

Mechanisms of stochastic focusing and defocusing in biological reaction networks: insight from accurate chemical master equation (ACME) solutions

Stochasticity plays important roles in regulation of biochemical reaction networks when the copy numbers of molecular species are small. Studies based on Stochastic Simulation Algorithm (SSA) has shown that a basic reaction system can display stochastic focusing (SF) by increasing the sensitivity of the network as a result of the signal noise. Although SSA has been widely used to study stochastic networks, it is ineffective in examining rare events and this becomes a significant issue when the tails of probability distributions are relevant as is the case of SF. Here we use the ACME method to solve the exact solution of the discrete Chemical Master Equations and to study a network where SF was reported. We showed that the level of SF depends on the degree of the fluctuations of signal molecule. We discovered that signaling noise under certain conditions in the same reaction network can lead to a decrease in the system sensitivities, thus the network can experience stochastic defocusing. These results highlight the fundamental role of stochasticity in biological reaction networks and the need for exact computation of probability landscape of the molecules in the system.

The dynamic mechanism of noisy signal decoding in gene regulation

Experimental evidence supports that signaling pathways can induce different dynamics of transcription factor (TF) activation, but how an input signal is encoded by such a dynamic, noisy TF and further decoded by downstream genes remains largely unclear. Here, using a system of stochastic transcription with signal regulation, we show that (1) keeping the intensity of the signal noise invariant but prolonging the signal duration can both enhance the mutual information (MI) and reduce the energetic cost (EC); (2) if the signal duration is fixed, the larger MI needs the larger EC, but if the signal period is fixed, there is an optimal time that the signal spends at one lower branch, such that MI reaches the maximum; (3) if both the period and the duration are simultaneously fixed, increasing the input noise can always enhance MI in the case of transcription regulation rather than in the case of degradation regulation. In addition, we find that the input noise can induce stochastic focusing in a regulation-dependent manner. These results reveal not only the dynamic mechanism of noisy signal decoding in gene regulation but also the essential role of external noise in controlling gene expression levels.

Minimal moment equations for stochastic models of biochemical reaction networks with partially finite state space

Many stochastic models of biochemical reaction networks contain some chemical species for which the number of molecules that are present in the system can only be finite (for instance due to conservation laws), but also other species that can be present in arbitrarily large amounts. The prime example of such networks are models of gene expression, which typically contain a small and finite number of possible states for the promoter but an infinite number of possible states for the amount of mRNA and protein. One of the main approaches to analyze such models is through the use of equations for the time evolution of moments of the chemical species. Recently, a new approach based on conditional moments of the species with infinite state space given all the different possible states of the finite species has been proposed. It was argued that this approach allows one to capture more details about the full underlying probability distribution with a smaller number of equations. Here, I show that the result that less moments provide more information can only stem from an unnecessarily complicated description of the system in the classical formulation. The foundation of this argument will be the derivation of moment equations that describe the complete probability distribution over the finite state space but only low-order moments over the infinite state space. I will show that the number of equations that is needed is always less than what was previously claimed and always less than the number of conditional moment equations up to the same order. To support these arguments, a symbolic algorithm is provided that can be used to derive minimal systems of unconditional moment equations for models with partially finite state space.