Role of functionality in two-component signal transduction: a stochastic study

The inverse correlation between robustness and sensitivity to autoregulation in two-component systems

Two-component systems (TCS) are signal transduction systems in bacteria and many other organisms that relay the sensory signal to genetic components. TCS consist of two proteins: a histidine kinase and a response regulator that the histidine kinase activates. This seemingly simple machinery can generate complex regulatory dynamics that enables the level of gene expression that matches the input signal: many TCS response regulators act on their own genes as transcription factors, resulting in a positive autoregulation mechanism. This regulation, in return, modulates the transcription factor activity as a function of the input signal. Positive autoregulation does not necessarily result in positive feedback. Sensitivity to autoregulation is quantified as the output level amplification resulting from the positive autoregulation mechanism. Another structural property of these systems is formally characterized as "robustness": in a robust TCS, the output of the system is solely a function of the input signal. Thus, a robust TCS remains insensitive to fluctuations in the concentrations of its protein components and, this way, maintains the precision in the output transcription factor activity in response to input stimulus. In this paper, we show with a formal model that TCS operate on a spectrum of inverse correlation between robustness and sensitivity to autoregulation. Our model predicts that the modulation by positive autoregulation is a function of loss in TCS robustness, for example, by spontaneous dephosphorylation of the histidine kinase. Consequently, the loss in robustness provides a proportional modulation by positive autoregulation to widen the response range with a scaled amplification of the output. At the other end of the spectrum, in the presence of a strictly robust TCS machinery, amplification of the transcription factor activity by autoregulation is diminished. We show that our results are in agreement with published experimental results. Our results suggest that these TCS evolve to converge at a trade-off between robustness and positive autoregulation.

Uncertainty propagation for deterministic models of biochemical networks using moment equations and the extended Kalman filter

Differential equation models of biochemical networks are frequently associated with a large degree of uncertainty in parameters and/or initial conditions. However, estimating the impact of this uncertainty on model predictions via Monte Carlo simulation is computationally demanding. A more efficient approach could be to track a system of low-order statistical moments of the state. Unfortunately, when the underlying model is nonlinear, the system of moment equations is infinite-dimensional and cannot be solved without a moment closure approximation which may introduce bias in the moment dynamics. Here, we present a new method to study the time evolution of the desired moments for nonlinear systems with polynomial rate laws. Our approach is based on solving a system of low-order moment equations by substituting the higher-order moments with Monte Carlo-based estimates from a small number of simulations, and using an extended Kalman filter to counteract Monte Carlo noise. Our algorithm provides more accurate and robust results compared to traditional Monte Carlo and moment closure techniques, and we expect that it will be widely useful for the quantification of uncertainty in biochemical model predictions.

Trade-offs between effectiveness and cost in bifunctional enzyme circuit with concentration robustness

A fundamental trade-off in biological systems is whether they consume resources to perform biological functions or save resources. Bacteria need to reliably and rapidly respond to input signals by using limited cellular resources. However, excessive resource consumption will become a burden for bacteria growth. To investigate the relationship between functional effectiveness and resource cost, we study the ubiquitous bifunctional enzyme circuit, which is robust to fluctuations in protein concentration and responds quickly to signal changes. We show that trade-off relationships exist between functional effectiveness and protein cost. Expressing more proteins of the circuit increases concentration robustness and response speed but affects bacterial growth. In particular, our study reveals a general relationship between free-energy dissipation rate, response speed, and concentration robustness. The dissipation of free energy plays an important role in the concentration robustness and response speed. High robustness can only be achieved with a large amount of free-energy consumption and protein cost. In addition, the noise of the output increases with increasing protein cost, while the noise of the response time decreases with increasing protein cost. We also calculate the trade-off relationships in the EnvZ-OmpR system and the nitrogen assimilation system, which both have the bifunctional enzyme. Similar results indicate that these relationships are mainly derived from the specific feature of the bifunctional enzyme circuits and are not relevant to the details of the models. According to the trade-off relationships, bacteria take a compromise solution that reliably performs biological functions at a reasonable cost.

Evidence of Robustness in a Two-Component System Using a Synthetic Circuit

Variation in the concentration of biological components is inescapable for any cell. Robustness in any biological circuit acts as a cushion against such variation and enables the cells to produce homogeneous output despite the fluctuation. The two-component system (TCS) with a bifunctional sensor kinase (that possesses both kinase and phosphatase activities) is proposed to be a robust circuit. Few theoretical models explain the robustness of a TCS, although the criteria and extent of robustness by these models differ. Here, we provide experimental evidence to validate the extent of the robustness of a TCS signaling pathway. We have designed a synthetic circuit in Escherichia coli using a representative TCS of Mycobacterium tuberculosis, MprAB, and monitored the in vivo output signal by systematically varying the concentration of either of the components or both. We observed that the output of the TCS is robust if the concentration of MprA is above a threshold value. This observation is further substantiated by two in vitro assays, in which we estimated the phosphorylated MprA pool or MprA-dependent transcription yield by varying either of the components of the TCS. This synthetic circuit could be used as a model system to analyze the relationship among different components of gene regulatory networks.IMPORTANCE Robustness in essential biological circuits is an important feature of the living organism. A few pieces of evidence support the existence of robustness in vivo in the two-component system (TCS) with a bifunctional sensor kinase (SK). The assays were done under physiological conditions in which the SK was much lower than the response regulator (RR). Here, using a synthetic circuit, we varied the concentrations of the SK and RR of a representative TCS to monitor output robustness in vivo. In vitro assays were also performed under conditions where the concentration of the SK was greater than that of the RR. Our results demonstrate the extent of output robustness in the TCS signaling pathway with respect to the concentrations of the two components.

Investigation of chemical noise in multisite phosphorylation chain using linear noise approximation

Quantitative and qualitative nature of chemical noise propagation in biochemical reaction networks depend crucially on the topology of the networks. Multisite reversible phosphorylation-dephosphorylation of target proteins is one such recurrently found topology that regulates host of key functions in living cells. Here we analytically calculated the stochasticity in multistep reversible chemical reactions by determining variance of phosphorylated species at the steady state using linear noise approximation to investigate the effect of mass action and Michaelis-Menten kinetics on the noise of phosphorylated species. We probed the dependence of noise on the number of phosphorylation sites and the equilibrium constants of the reaction equilibria to investigate the chemical noise propagation in the multisite phosphorylation chain.

Interplay of synergy and redundancy in diamond motif

The formalism of partial information decomposition provides a number of independent components which altogether constitute the total information provided by the source variable(s) about the target variable(s). These non-overlapping terms are recognized as unique information, synergistic information, and redundant information. The metric of net synergy conceived as the difference between synergistic and redundant information is capable of detecting effective synergy, effective redundancy, and information independence among stochastic variables. The net synergy can be quantified using appropriate combinations of different Shannon mutual information terms. The utilization of the net synergy in network motifs with the nodes representing different biochemical species, involved in information sharing, uncovers rich store for exciting results. In the current study, we use this formalism to obtain a comprehensive understanding of the relative information processing mechanism in a diamond motif and two of its sub-motifs, namely, bifurcation and integration motif embedded within the diamond motif. The emerging patterns of effective synergy and effective redundancy and their contribution toward ensuring high fidelity information transmission are duly compared in the sub-motifs. Investigation on the metric of net synergy in independent bifurcation and integration motifs are also executed. In all of these computations, the crucial roles played by various systemic time scales, activation coefficients, and signal integration mechanisms at the output of the network topologies are especially emphasized. Following this plan of action, we become confident that the origin of effective synergy and effective redundancy can be architecturally justified by decomposing a diamond motif into bifurcation and integration motif. According to our conjecture, the presence of a common source of fluctuations creates effective redundancy. Our calculations reveal that effective redundancy empowers signal fidelity. Moreover, to achieve this, input signaling species avoids strong interaction with downstream intermediates. This strategy is capable of making the diamond motif noise-tolerant. Apart from the topological features, our study also puts forward the active contribution of additive and multiplicative signal integration mechanisms to nurture effective redundancy and effective synergy.

Transient absolute robustness in stochastic biochemical networks

Absolute robustness allows biochemical networks to sustain a consistent steady-state output in the face of protein concentration variability from cell to cell. This property is structural and can be determined from the topology of the network alone regardless of rate parameters. An important question regarding these systems is the effect of discrete biochemical noise in the dynamical behaviour. In this paper, a variable freezing technique is developed to show that under mild hypotheses the corresponding stochastic system has a transiently robust behaviour. Specifically, after finite time the distribution of the output approximates a Poisson distribution, centred around the deterministic mean. The approximation becomes increasingly accurate, and it holds for increasingly long finite times, as the total protein concentrations grow to infinity. In particular, the stochastic system retains a transient, absolutely robust behaviour corresponding to the deterministic case. This result contrasts with the long-term dynamics of the stochastic system, which eventually must undergo an extinction event that eliminates robustness and is completely different from the deterministic dynamics. The transiently robust behaviour may be sufficient to carry out many forms of robust signal transduction and cellular decision-making in cellular organisms.

Redundancy in information transmission in a two-step cascade

We present a stochastic framework to study signal transmission in a generic two-step cascade S→X→Y. Starting from a set of Langevin equations obeying Gaussian noise processes we calculate the variance and covariance while considering both linear and nonlinear production terms for different biochemical species of the cascade. These quantities are then used to calculate the net synergy within the purview of partial information decomposition. We show that redundancy in information transmission is essentially an important consequence of Markovian property of the two-step cascade motif. We also show that redundancy increases fidelity of the signaling pathway.