Quantification of noise in the bifunctionality-induced post-translational modification

We present a generic analytical scheme for the quantification of fluctuations due to bifunctionality-induced signal transduction within the members of bacterial two-component system. The proposed model takes into account post-translational modifications in terms of elementary phosphotransfer kinetics. Sources of fluctuations due to autophosphorylation, kinase and phosphatase activity of the sensor kinase have been considered in the model via Langevin equations, which are then solved within the framework of linear noise approximation. The resultant analytical expression of phosphorylated response regulators are then used to quantify the noise profile of biologically motivated single and branched pathways. Enhancement and reduction of noise in terms of extra phosphate outflux and influx, respectively, have been analyzed for the branched system. Furthermore, role of fluctuations of the network output in the regulation of a promoter with random activation/deactivation dynamics has been analyzed.

💡 Research Summary

The manuscript presents a comprehensive analytical framework for quantifying stochastic fluctuations that arise from bifunctionality‑driven signal transduction in bacterial two‑component systems (TCS). In a typical TCS, a membrane‑bound sensor kinase (SK) undergoes autophosphorylation, transfers the phosphate to a cytoplasmic response regulator (RR), and can also act as a phosphatase, removing the phosphate from RR. These three enzymatic activities—autophosphorylation, kinase, and phosphatase—constitute a bifunctional module that is intrinsically noisy because the participating molecules often exist in low copy numbers.

Model Construction
The authors first write down the elementary mass‑action reactions:

1. SK + ATP → SK~P + ADP (autophosphorylation, rate k_aut)
2. SKP + RR → SK + RRP (kinase, rate k_kin)
3. SK + RR~P → SK + RR (phosphatase, rate k_phos)

To capture stochasticity, each reaction is represented by a Langevin equation that adds a Gaussian white‑noise term ξ_i(t) to the deterministic rate equation. The noise intensity D_i is taken proportional to the mean flux of the corresponding reaction, which is a standard assumption in chemical Langevin formulations. The system of stochastic differential equations is then linearized around the deterministic steady state, yielding a set of linear Langevin equations amenable to the Linear Noise Approximation (LNA). By solving the associated Lyapunov equation, the authors obtain closed‑form expressions for the covariance matrix of the molecular species, in particular the variance of phosphorylated RR (RR~P).

Key Analytical Results
The mean concentration of RR~P follows the familiar balance between kinase and phosphatase fluxes, but the variance contains contributions from all three noise sources. The analytical variance expression shows that fluctuations in the autophosphorylation step dominate the overall noise budget, because this step is the primary source of phosphate into the system. The variance also depends on the total protein copy numbers, the kinetic constants, and the network topology.

Single vs. Branched Pathways
The paper extends the analysis to a branched architecture where a single SK phosphorylates two distinct RRs (RR1 and RR2). In this configuration, the phosphate “outflux” to the second RR acts as an additional sink for SKP, effectively increasing the noise in both RRP species. The authors demonstrate analytically that the variance of each RRP grows with the branching ratio, reflecting a redistribution of stochastic phosphate fluxes. Conversely, introducing an external phosphate influx (e.g., a constitutively active SKP source) reduces the relative contribution of intrinsic autophosphorylation noise, leading to a net noise attenuation—a phenomenon the authors term “buffering by influx.”

Coupling to Gene Expression
To assess the functional impact of RRP fluctuations, the authors couple the stochastic RRP dynamics to a downstream promoter that switches randomly between active and inactive states with rates k_on and k_off. The promoter activation probability depends on the instantaneous RRP concentration, creating a two‑stage stochastic cascade. By propagating the covariance through this cascade, the authors derive an expression for the variance of the output protein (e.g., GFP). Their analysis reveals a non‑linear transfer function: in regimes where k_on is comparable to the RRP fluctuation timescale, noise is amplified, whereas when k_on is either much faster or much slower, the promoter acts as a low‑pass filter, attenuating upstream noise. This result highlights how cells can tune promoter kinetics to either exploit or suppress signaling noise.

Biological Implications and Design Principles
The study provides several insights of broad relevance:

1. Noise as a Design Parameter – Bifunctional SKs can modulate cellular noise by adjusting the balance between kinase and phosphatase activities.
2. Branching as a Noise Amplifier – Pathway branching, common in complex bacterial signaling networks, inherently raises stochastic variability, which may be beneficial for bet‑hedging strategies.
3. Inflow‑Induced Noise Suppression – External phosphate sources or constitutively active SK~P can serve as “noise buffers,” a principle that could be harnessed in synthetic circuits.
4. Promoter Kinetics as a Noise Filter – By engineering promoter switching rates, synthetic biologists can control the degree to which upstream signaling noise propagates to gene expression outputs.

Limitations and Future Directions
The Linear Noise Approximation assumes small fluctuations around a stable steady state; thus, its accuracy may diminish in regimes of strong non‑linearity or very low molecule numbers. Moreover, the model neglects feedback loops where RRP can regulate SK activity—a common motif in real TCSs. Extending the framework to include such feedback, as well as spatial heterogeneity and cell‑to‑cell coupling, would provide a more complete picture. Experimental validation using single‑cell fluorescence microscopy to simultaneously monitor SKP, RR~P, and reporter expression would be a logical next step.

Conclusion
Overall, the paper delivers a rigorous, analytically tractable description of stochastic dynamics in bifunctional two‑component signaling systems. By integrating Langevin formalism with linear noise analysis, it quantifies how enzymatic activities, pathway architecture, and downstream promoter kinetics shape the noise landscape. These findings not only deepen our understanding of bacterial signal transduction but also furnish actionable design rules for constructing synthetic circuits with predictable noise characteristics.