A hidden integral structure endows absolute concentration robust systems with resilience to dynamical concentration disturbances

Robust homeostasis of cellular cholesterol is a consequence of endogenous antithetic integral control

Although cholesterol is essential for cellular viability and proliferation, it is highly toxic in excess. The concentration of cellular cholesterol must therefore be maintained within tight tolerances, and is thought to be subject to a stringent form of homeostasis known as Robust Perfect Adaptation (RPA). While much is known about the cellular signalling interactions involved in cholesterol regulation, the specific chemical reaction network structures that might be responsible for the robust homeostatic regulation of cellular cholesterol have been entirely unclear until now. In particular, the molecular mechanisms responsible for sensing excess whole-cell cholesterol levels have not been identified previously, and no mathematical models to date have been able to capture an integral control implementation that could impose RPA on cellular cholesterol. Here we provide a detailed mathematical description of cholesterol regulation pathways in terms of biochemical reactions, based on an extensive review of experimental and clinical literature. We are able to decompose the associated chemical reaction network structures into several independent subnetworks, one of which is responsible for conferring RPA on several intracellular forms of cholesterol. Remarkably, our analysis reveals that RPA in the cholesterol concentration in the endoplasmic reticulum (ER) is almost certainly due to a well-characterised control strategy known as antithetic integral control which, in this case, involves the high-affinity binding of a multi-molecular transcription factor complex with cholesterol molecules that are excluded from the ER membrane. Our model provides a detailed framework for exploring the necessary biochemical conditions for robust homeostatic control of essential and tightly regulated cellular molecules such as cholesterol.

Universal structures for adaptation in biochemical reaction networks

At the molecular level, the evolution of life is driven by the generation and diversification of adaptation mechanisms. A universal description of adaptation-capable chemical reaction network (CRN) structures has remained elusive until now, since currently-known criteria for adaptation apply only to a tiny subset of possible CRNs. Here we identify the definitive structural requirements that characterize all adaptation-capable collections of interacting molecules, however large or complex. We show that these network structures implement a form of integral control in which multiple independent integrals can collaborate to confer the capacity for adaptation on specific molecules. Using an algebraic algorithm informed by these findings, we demonstrate the existence of embedded integrals in a variety of biologically important CRNs that have eluded previous methods, and for which adaptation has been observed experimentally. This definitive picture of biological adaptation at the level of intermolecular interactions represents a blueprint for adaptation-capable signaling networks across all domains of life, and for the design of synthetic biosystems.

Protein-protein complexes can undermine ultrasensitivity-dependent biological adaptation

Robust perfect adaptation (RPA) is a ubiquitously observed signalling response across all scales of biological organization. A major class of network architectures that drive RPA in complex networks is the Opposer module-a feedback-regulated network into which specialized integral-computing 'opposer node(s)' are embedded. Although ultrasensitivity-generating chemical reactions have long been considered a possible mechanism for such adaptation-conferring opposer nodes, this hypothesis has relied on simplified Michaelian models, which neglect the presence of protein-protein complexes. Here we develop complex-complete models of interlinked covalent-modification cycles with embedded ultrasensitivity, explicitly capturing all molecular interactions and protein complexes. Strikingly, we demonstrate that the presence of protein-protein complexes thwarts the network's capacity for RPA in any 'free' active protein form, conferring RPA capacity instead on the concentration of a larger protein pool consisting of two distinct forms of a single protein. We further show that the presence of enzyme-substrate complexes, even at comparatively low concentrations, play a crucial and previously unrecognized role in controlling the RPA response-significantly reducing the range of network inputs for which RPA can obtain, and imposing greater parametric requirements on the RPA response. These surprising results raise fundamental new questions as to the biochemical requirements for adaptation-conferring Opposer modules within complex cellular networks.

Foundations of static and dynamic absolute concentration robustness

Absolute Concentration Robustness (ACR) was introduced by Shinar and Feinberg (Science 327:1389-1391, 2010) as robustness of equilibrium species concentration in a mass action dynamical system. Their aim was to devise a mathematical condition that will ensure robustness in the function of the biological system being modeled. The robustness of function rests on what we refer to as empirical robustness-the concentration of a species remains unvarying, when measured in the long run, across arbitrary initial conditions. Even simple examples show that the ACR notion introduced in Shinar and Feinberg (Science 327:1389-1391, 2010) (here referred to as static ACR) is neither necessary nor sufficient for empirical robustness. To make a stronger connection with empirical robustness, we define dynamic ACR, a property related to long-term, global dynamics, rather than only to equilibrium behavior. We discuss general dynamical systems with dynamic ACR properties as well as parametrized families of dynamical systems related to reaction networks. We find necessary and sufficient conditions for dynamic ACR in complex balanced reaction networks, a class of networks that is central to the theory of reaction networks.

Transition graph decomposition for complex balanced reaction networks with non-mass-action kinetics

Reaction networks are widely used models to describe biochemical processes. Stochastic fluctuations in the counts of biological macromolecules have amplified consequences due to their small population sizes. This makes it necessary to favor stochastic, discrete population, continuous time models. The stationary distributions provide snapshots of the model behavior at the stationary regime, and as such finding their expression in terms of the model parameters is of great interest. The aim of the present paper is to describe when the stationary distributions of the original model, whose state space is potentially infinite, coincide exactly with the stationary distributions of the process truncated to finite subsets of states, up to a normalizing constant. The finite subsets of states we identify are called copies and are inspired by the modular topology of reaction network models. With such a choice we prove a novel graphical characterization of the concept of complex balancing for stochastic models of reaction networks. The results of the paper hold for the commonly used mass-action kinetics but are not restricted to it, and are in fact stated for more general setting.

RNA Compensation: A Positive Feedback Insulation Strategy for RNA-Based Transcription Networks

The lack of signaling modularity of biomolecular systems poses major challenges toward engineering complex networks. Directional signaling between an upstream and a downstream circuit requires the presence of binding events, which result in the consumption of regulatory molecules and can compromise the operation of the upstream circuit. This issue has been previously addressed by introducing insulation strategies that include high-gain negative feedback and activation-deactivation reaction cycles. In this paper, we focus on RNA-based circuits and propose a new positive-feedback strategy to mitigate signal consumption that we propose occurs for each regulatory event due to irreversible binding of the RNA input to the RNA target. To mitigate this, an extra RNA input is added in tandem with transcription output to compensate the RNA consumption, leading to concentration robustness of the input RNA molecule regardless of the amount of downstream modules. We term this strategy RNA compensation, and it can be applied to systems that have a stringent input-output gain, such as Small Transcription Activating RNAs (STARs). Our theoretical analysis shows that RNA compensation not only eliminates the signaling consumption in individual STAR-based regulators, but also improves the composability of STAR cascades and the modularity of RNA bistable systems.