Sensitivity summation theorems for stochastic biochemical reaction systems

Synthetic biology and regulatory networks: where metabolic systems biology meets control engineering

Metabolic pathways can be engineered to maximize the synthesis of various products of interest. With the advent of computational systems biology, this endeavour is usually carried out through in silico theoretical studies with the aim to guide and complement further in vitro and in vivo experimental efforts. Clearly, what counts is the result in vivo, not only in terms of maximal productivity but also robustness against environmental perturbations. Engineering an organism towards an increased production flux, however, often compromises that robustness. In this contribution, we review and investigate how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed by systems and control engineering. While trade-offs between production optimality and cellular robustness have already been studied diagnostically and statically, the dynamics also matter. Integration of the dynamic design aspects of control engineering with the more diagnostic aspects of metabolic, hierarchical control and regulation analysis is leading to the new, conceptual and operational framework required for the design of robust and productive dynamic pathways.

Ranking of Reactions Based on Sensitivity of Protein Noise Depends on the Choice of Noise Measure

Gene expression is a stochastic process. Identification of the step maximally affecting noise in the protein level is an important aspect of investigation of gene product distribution. There are numerous experimental and theoretical studies that seek to identify this important step. However, these studies have used two different measures of noise, viz. coefficient of variation and Fano factor, and have compared different processes leading to contradictory observations regarding the important step. In this study, we performed systematic global and local sensitivity analysis on two models of gene expression to investigate relative contribution of reaction rate parameters to steady state noise in the protein level using both the measures of noise. We analytically and computationally showed that the ranking of parameters based on the sensitivity of the noise to variation in a given parameter is a strong function of the choice of the noise measure. If the Fano factor is used as the noise measure, translation is the important step whereas for coefficient of variation, transcription is the important step. We derived an analytical expression for local sensitivity and used it to explain the distinct contributions of each reaction parameter to the two measures of noise. We extended the analysis to a generic linear catalysis reaction system and observed that the reaction network topology was an important factor influencing the local sensitivity of the two measures of noise. Our study suggested that, for the analysis of contributions of reactions to the noise, consideration of both the measures of noise is important.

Controlling E. coli Gene Expression Noise

Intracellular protein copy numbers show significant cell-to-cell variability within an isogenic population due to the random nature of biological reactions. Here we show how the variability in copy number can be controlled by perturbing gene expression. Depending on the genetic network and host, different perturbations can be applied to control variability. To understand more fully how noise propagates and behaves in biochemical networks we developed stochastic control analysis (SCA) which is a sensitivity-based analysis framework for the study of noise control. Here we apply SCA to synthetic gene expression systems encoded on plasmids that are transformed into Escherichia coli. We show that (1) dual control of transcription and translation efficiencies provides the most efficient way of noise-versus-mean control. (2) The expressed proteins follow the gamma distribution function as found in chromosomal proteins. (3) One of the major sources of noise, leading to the cell-to-cell variability in protein copy numbers, is related to bursty translation. (4) By taking into account stochastic fluctuations in autofluorescence, the correct scaling relationship between the noise and mean levels of the protein copy numbers was recovered for the case of weak fluorescence signals.

Stochastic flux analysis of chemical reaction networks

Background

Chemical reaction networks provide an abstraction scheme for a broad range of models in biology and ecology. The two common means for simulating these networks are the deterministic and the stochastic approaches. The traditional deterministic approach, based on differential equations, enjoys a rich set of analysis techniques, including a treatment of reaction fluxes. However, the discrete stochastic simulations, which provide advantages in some cases, lack a quantitative treatment of network fluxes.

Results

We describe a method for flux analysis of chemical reaction networks, where flux is given by the flow of species between reactions in stochastic simulations of the network. Extending discrete event simulation algorithms, our method constructs several data structures, and thereby reveals a variety of statistics about resource creation and consumption during the simulation. We use these structures to quantify the causal interdependence and relative importance of the reactions at arbitrary time intervals with respect to the network fluxes. This allows us to construct reduced networks that have the same flux-behavior, and compare these networks, also with respect to their time series. We demonstrate our approach on an extended example based on a published ODE model of the same network, that is, Rho GTP-binding proteins, and on other models from biology and ecology.

Conclusions

We provide a fully stochastic treatment of flux analysis. As in deterministic analysis, our method delivers the network behavior in terms of species transformations. Moreover, our stochastic analysis can be applied, not only at steady state, but at arbitrary time intervals, and used to identify the flow of specific species between specific reactions. Our cases study of Rho GTP-binding proteins reveals the role played by the cyclic reverse fluxes in tuning the behavior of this network.

Adaptation for protein synthesis efficiency in a naturally occurring self-regulating operon

The korAB operon in RK2 plasmids is a beautiful natural example of a negatively and cooperatively self-regulating operon. It has been particularly well characterized both experimentally and with mathematical models. We have carried out a detailed investigation of the role of the regulatory mechanism using a biologically grounded mechanistic multi-scale stochastic model that includes plasmid gene regulation and replication in the context of host growth and cell division. We use the model to compare four hypotheses for the action of the regulatory mechanism: increased robustness to extrinsic factors, decreased protein fluctuations, faster response-time of the operon and reduced host burden through improved efficiency of protein production. We find that the strongest impact of all elements of the regulatory architecture is on improving the efficiency of protein synthesis by reduction in the number of mRNA molecules needed to be produced, leading to a greater than ten-fold reduction in host energy required to express these plasmid proteins. A smaller but still significant role is seen for speeding response times, but this is not materially improved by the cooperativity. The self-regulating mechanisms have the least impact on protein fluctuations and robustness. While reduction of host burden is evident in a plasmid context, negative self-regulation is a widely seen motif for chromosomal genes. We propose that an important evolutionary driver for negatively self-regulated genes is to improve the efficiency of protein synthesis.

Adjusting phenotypes by noise control

Genetically identical cells can show phenotypic variability. This is often caused by stochastic events that originate from randomness in biochemical processes involving in gene expression and other extrinsic cellular processes. From an engineering perspective, there have been efforts focused on theory and experiments to control noise levels by perturbing and replacing gene network components. However, systematic methods for noise control are lacking mainly due to the intractable mathematical structure of noise propagation through reaction networks. Here, we provide a numerical analysis method by quantifying the parametric sensitivity of noise characteristics at the level of the linear noise approximation. Our analysis is readily applicable to various types of noise control and to different types of system; for example, we can orthogonally control the mean and noise levels and can control system dynamics such as noisy oscillations. As an illustration we applied our method to HIV and yeast gene expression systems and metabolic networks. The oscillatory signal control was applied to p53 oscillations from DNA damage. Furthermore, we showed that the efficiency of orthogonal control can be enhanced by applying extrinsic noise and feedback. Our noise control analysis can be applied to any stochastic model belonging to continuous time Markovian systems such as biological and chemical reaction systems, and even computer and social networks. We anticipate the proposed analysis to be a useful tool for designing and controlling synthetic gene networks.

Stochastic gene expression at the single cell level can lead to significant phenotypic variation at the population level. To obtain a desired phenotype, the noise levels of intracellular protein concentrations may need to be tuned and controlled. Noise levels often decrease in relative amount as the mean values increase. This implies that the noise levels can be passively controlled by changing the mean values. In an engineering perspective, the noise levels can be further controlled while the mean values can be simultaneously adjusted to desired values. Here, systematic schemes for such simultaneous control are described by identifying where and by how much the system needs to be perturbed. The schemes can be applied to the design process of a potential therapeutic HIV-drug that targets a certain set of reactions that are identified by the proposed analysis, to prevent stochastic transition to the lytic state. In some cases, the simultaneous control cannot be performed efficiently, when the noise levels strongly change with the mean values. This problem is shown to be resolved by applying extra noise and feedback.