Asymptotic Robustness in Biochemical Systems

TL;DR

The paper addresses how approximate concentration robustness can emerge in biochemical networks without exact ACR motifs or finely tuned parameters. It introduces asymptotic ACR (aACR) and proves that aACR arises generically from network structure in systems with positive conservation laws, using an algebraic-geometry–based framework. Through analyses of the EnvZ-OmpR signaling system and archetypal networks, it demonstrates that aACR is widespread and often more robust than classical ACR, persisting under modest structural changes. The work provides a practical toolkit—rooted in polynomial constraints and Gröbner-basis–style methods—to identify aACR relationships and fully characterize steady-state responses, with implications for understanding biological robustness and guiding synthetic circuit design.

Abstract

Living systems maintain stable internal states despite environmental fluctuations. Absolute concentration robustness (ACR) is a striking homeostatic phenomenon in which the steady-state concentration of a species remains invariant despite changes in total supply. Although experimental studies have reported approximate-but not exact-robustness in steady-state concentrations, such behavior has often been attributed to exact ACR motifs perturbed by measurement noise or minor reactions, rather than recognized as a structural property of a network itself. Here, we introduce a previously underappreciated phenomenon, namely asymptotic ACR (aACR): approximate robustness can emerge solely from the network structure, without requiring exact ACR motifs or negligible parameters. We find that aACR is more pervasive than classical ACR, as demonstrated in systems such as the Escherichia coli EnvZ-OmpR osmoregulation system and a futile cycle. Furthermore, we prove that such ubiquity stems solely from network structure without fine-tuning of kinetic parameters. The notion of aACR provides a rigorous and practical tool to analyze robust responses in broad biochemical systems.

Figures (3)

• Figure 1: Illustration of an experimental setup for assessing concentration robustness. (a) A series of experiments is conducted by varying the concentration of the dose (i.e., input species). (b) The corresponding steady-state concentrations of the response (i.e., the output species) are measured. (c) The input–output pairs are used to construct a dose–response curve. If the output concentrations remain similar despite changes in the input, the system is considered to exhibit concentration robustness of the output species with respect to the input species.
• Figure 1: Phase planes and dose-response curves of the reaction networks in \ref{['main_arXiv_251031:crn:archetypal']} and \ref{['main_arXiv_251031:crn:archetypal-modified']}. (a) The reaction network diagram given in \ref{['main_arXiv_251031:crn:archetypal']}. (b) The phase plane describing the dynamics of the reaction network in \ref{['main_arXiv_251031:crn:archetypal']}. Orange lines represent the set of steady states. (c) The dose-response curve of the steady-state concentration of $X$, denoted by $X^{\text{SS}}$, with respect to the conserved quantity in the reaction network, $T=X(0)+Y(0)$. (d) A reaction network modified from (a) by adding one more reaction, $X \to Y$, given in \ref{['main_arXiv_251031:crn:archetypal-modified']}. (e) The phase-plane of the modified network. (f) The dose-response curve the modified network. Although $X^{\text{SS}}$ no longer becomes a constant function, it converges to the same value. Notably, the two dose-response curves in (c) and (f) are nearly indistinguishable.
• Figure 2: Biologically relevant networks showing ACR and aACR phenomena. (a) The EnvZ-OmpR signaling network, a well-studied biochemical system, serves as a key example of a network with ACR in species $Y_\text{P}$. (b) Dose-response curves for $Y_\text{P}$ show that its steady-state concentration reaches a fixed value (horizontal line) when varying the initial concentration of $X$ (orange curve) or $Y$ (blue curve), demonstrating ACR. Here, $Y_\text{P}^{\text{SS}}$ denotes the steady-state concentration of $Y_\text{P}$. Note that these two curves overlap as they reach the same $Y_\text{P}^{\text{SS}}$ value. (c) In contrast, $X_\text{P}Y$ exhibits aACR, denoted by dashed lines. This is demonstrated by the dose-response curve of its steady-state concentration, denoted by $X_\text{P}Y^{\text{SS}}$, approaching but never precisely reaching a fixed value as the initial concentrations of $X$ or $Y$ are increasing. (d) A network modified from the network in (a) by making one irreversible reaction, $XT \to X_\text{P}$, reversible. (e) In this modified network, $Y_\text{P}$ no longer exhibits ACR: its steady-state concentration, $Y_\text{P}^{\text{SS}}$, approaches zero when varying $X(0)$ and shows aACR when varying $Y(0)$. (f) Despite this modification, $X_\text{P}Y$ retains aACR, highlighting the robustness of this property compared to ACR.

Theorems & Definitions (16)

• Definition 3.1
• Definition 3.2
• Definition 3.3
• Definition 5.1
• Definition 5.2
• Theorem 5.3
• Lemma 5.4
• Proof 1: Proof of Lemma \ref{['main_arXiv_251031:lem:LimitExists']}
• Remark 5.5
• Theorem 5.6