Operating regimes of signaling cycles: statics, dynamics, and noise filtering

Generating synthetic signaling networks for in silico modeling studies

Predictive models of signaling pathways have proven to be difficult to develop. Traditional approaches to developing mechanistic models rely on collecting experimental data and fitting a single model to that data. This approach works for simple systems but has proven unreliable for complex systems such as biological signaling networks. Thus, there is a need to develop new approaches to create predictive mechanistic models of complex systems. To meet this need, we developed a method for generating artificial signaling networks that were reasonably realistic and thus could be treated as ground truth models. These synthetic models could then be used to generate synthetic data for developing and testing algorithms designed to recover the underlying network topology and associated parameters. We defined the reaction degree and reaction distance to measure the topology of reaction networks, especially to consider enzymes. To determine whether our generated signaling networks displayed meaningful behavior, we compared them with signaling networks from the BioModels Database. This comparison indicated that our generated signaling networks had high topological similarities with BioModels signaling networks with respect to the reaction degree and distance distributions. In addition, our synthetic signaling networks had similar behavioral dynamics with respect to both steady states and oscillations, suggesting that our method generated synthetic signaling networks comparable with BioModels and thus could be useful for building network evaluation tools.

First-order ultrasensitivity in phosphorylation cycles

Cellular signal transduction takes place through a network of phosphorylation cycles. These pathways take the form of a multi-layered cascade of cycles. This work focuses on the sensitivity of single, double and n length cycles. Cycles that operate in the zero-order regime can become sensitive to changes in signal, resulting in zero-order ultrasensitivity (ZOU). Using frequency analysis, we confirm previous efforts that cascades can act as noise filters by computing the bandwidth. We show that n length cycles display what we term first-order ultrasensitivity which occurs even when the cycles are not operating in the zero-order regime. The magnitude of the sensitivity, however, has an upper bound equal to the number of cycles. It is known that ZOU can be significantly reduced in the presence of retroactivity. We show that the first-order ultrasensitivity is immune to retroactivity and that the ZOU and first-order ultrasensitivity can be blended to create systems with constant sensitivity over a wider range of signal. We show that the ZOU in a double cycle is only modestly higher compared with a single cycle. We therefore speculate that the double cycle has evolved to enable amplification even in the face of retroactivity.

Implementation of an Ultrasensitive Biomolecular Controller for Enzymatic Reaction Processes With Delay Using DNA Strand Displacement

In this article, a set of abstract chemical reactions has been employed to construct a novel nonlinear biomolecular controller, i.e, the Brink controller (BC) with direct positive autoregulation (DPAR) (namely BC-DPAR controller). In comparison to dual rail representation-based controllers such as the quasi sliding mode (QSM) controller, the BC-DPAR controller directly reduces the number of CRNs required for realizing an ultrasensitive input-output response because it does not involve the subtraction module, reducing the complexity of DNA implementations. Then, the action mechanism and steady-state condition constraints of two nonlinear controllers, BC-DPAR controller and QSM controller, are investigated further. Considering the mapping relationship between CRNs and DNA implementation, a CRNs-based enzymatic reaction process with delay is constructed, and a DNA strand displacement (DSD) scheme representing time delay is proposed. The BC-DPAR controller, when compared to the QSM controller, can reduce the number of abstract chemical reactions and DSD reactions required by 33.3% and 31.8%, respectively. Finally, an enzymatic reaction scheme with BC-DPAR controller is designed using DSD reactions. According to the findings, the enzymatic reaction process's output substance can approach the target level at a quasi-steady state in both delay-free and non-zero delay conditions, but the target level can only be achieved during a finite-time period, mainly due to the fuel stand depletion.

Sensitivity and Frequency Response of Biochemical Cascades

Signal transduction from a cell's surface to cytoplasmic and nuclear targets takes place through a complex network of interconnected pathways. Phosphorylation cycles are common components of many pathways and may take the form of a multi-layered cascade of cycles or incorporate species with multiple phosphorylation sites that effectively create a sequence of cycles with increasing states of phosphorylation. This work focuses on the frequency response and sensitivity of such systems, two properties that have not been thoroughly examined. Starting with a singularly phosphorylated single-cycle system, we compare the sensitivity to perturbation at steady-state across a range of input signal strengths. This is followed by a frequency response analysis focusing on the gain and associated bandwidth. Next, we consider a two-layer cascade of single phosphorylation cycles and focus on how the two cycles interact to produce various effects on the bandwidth and damping properties. Then we consider the (ultra)sensitivity of a doubly phosphorylated system, where we describe in detail first-order ultrasensitivity, a unique property of these systems, which can be blended with zero-order ultrasensitivity to create systems with relatively constant gain over a range of signal input. Finally, we give an in-depth analysis of the sensitivity of an n-phosphorylated system.

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.

Dynamics and Sensitivity of Signaling Pathways

Purpose of Review

Signaling pathways serve to communicate information about extracellular conditions into the cell, to both the nucleus and cytoplasmic processes to control cell responses. Genetic mutations in signaling network components are frequently associated with cancer and can result in cells acquiring an ability to divide and grow uncontrollably. Because signaling pathways play such a significant role in cancer initiation and advancement, their constituent proteins are attractive therapeutic targets. In this review, we discuss how signaling pathway modeling can assist with identifying effective drugs for treating diseases, such as cancer. An achievement that would facilitate the use of such models is their ability to identify controlling biochemical parameters in signaling pathways, such as molecular abundances and chemical reaction rates, because this would help determine effective points of attack by therapeutics.

Recent Findings

We summarize the current state of understanding the sensitivity of phosphorylation cycles with and without sequestration. We also describe some basic properties of regulatory motifs including feedback and feedforward regulation.

Summary

Although much recent work has focused on understanding the dynamics and particularly the sensitivity of signaling networks in eukaryotic systems, there is still an urgent need to build more scalable models of signaling networks that can appropriately represent their complexity across different cell types and tumors.

Frequency-preference response in covalent modification cycles under substrate sequestration conditions

Covalent modification cycles (CMCs) are basic units of signaling systems and their properties are well understood. However, their behavior has been mostly characterized in situations where the substrate is in excess over the modifying enzymes. Experimental data on protein abundance suggest that the enzymes and their target proteins are present in comparable concentrations, leading to substrate sequestration by the enzymes. In this enzyme-in-excess regime, CMCs have been shown to exhibit signal termination, the ability of the product to return to a stationary value lower than its peak in response to constant stimulation, while this stimulation is still active, with possible implications for the ability of systems to adapt to environmental inputs. We characterize the conditions leading to signal termination in CMCs in the enzyme-in-excess regime. We also demonstrate that this behavior leads to a preferred frequency response (band-pass filters) when the cycle is subjected to periodic stimulation, whereas the literature reports that CMCs investigated so far behave as low-pass filters. We characterize the relationship between signal termination and the preferred frequency response to periodic inputs and we explore the dynamic mechanism underlying these phenomena. Finally, we describe how the behavior of CMCs is reflected in similar types of responses in the cascades of which they are part. Evidence of protein abundance in vivo shows that enzymes and substrates are present in comparable concentrations, thus suggesting that signal termination and frequency-preference response to periodic inputs are also important dynamic features of cell signaling systems, which have been overlooked.

Retroactivity induced operating regime transition in an enzymatic futile cycle

Activated phosphorylation-dephosphorylation biochemical reaction cycles are a class of enzymatic futile cycles. A futile cycle such as a single MAPK cascade governed by two underlying enzymatic reactions permits Hyperbolic (H), Signal transducing (ST), Threshold-hyperbolic (TH) and Ultrasensitive (U) operating regimes that characterize input-output behaviour. Retroactive signalling caused by load due to sequestration of phosphorylated or unphosphorylated form of the substrate in a single enzymatic cascade without explicit feedback can introduce two-way communication, a feature not possible otherwise. We systematically characterize the operating regimes of a futile cycle subject to retroactivity in either of the substrate forms. We demonstrate that increasing retroactivity strength, which quantifies the downstream load, can trigger five possible regime transitions. Retroactivity strength is a reflection of the fraction of the substrate sequestered by its downstream target. Remarkably, the minimum required retroactivity strength to evidence any sequestration triggered regime transition demands 23% of the substrate bound to its downstream target. This minimum retroactivity strength corresponds to the transition of the dose-response curve from ST to H regime. We show that modulation of the saturation and unsaturation levels of the enzymatic reactions by retroactivity is the fundamental mechanism governing operating regime transition.

Crosstalk and ultrasensitivity in protein degradation pathways

Protein turnover is vital to cellular homeostasis. Many proteins are degraded efficiently only after they have been post-translationally “tagged” with a polyubiquitin chain. Ubiquitylation is a form of Post-Translational Modification (PTM): addition of a ubiquitin to the chain is catalyzed by E3 ligases, and removal of ubiquitin is catalyzed by a De-UBiquitylating enzyme (DUB). Nearly four decades ago, Goldbeter and Koshland discovered that reversible PTM cycles function like on-off switches when the substrates are at saturating concentrations. Although this finding has had profound implications for the understanding of switch-like behavior in biochemical networks, the general behavior of PTM cycles subject to synthesis and degradation has not been studied. Using a mathematical modeling approach, we found that simply introducing protein turnover to a standard modification cycle has profound effects, including significantly reducing the switch-like nature of the response. Our findings suggest that many classic results on PTM cycles may not hold in vivo where protein turnover is ubiquitous. We also found that proteins sharing an E3 ligase can have closely related changes in their expression levels. These results imply that it may be difficult to interpret experimental results obtained from either overexpressing or knocking down protein levels, since changes in protein expression can be coupled via E3 ligase crosstalk. Understanding crosstalk and competition for E3 ligases will be key in ultimately developing a global picture of protein homeostasis.

Previous work has shown that substrates of Post-Translational Modification (PTM) cycles can have coupled responses if those substrates share enzymes. This implies that modifications leading to substrate degradation (e.g. ubiquitylation by an E3 ligase) could introduce coupling in concentrations of substrates sharing a ligase. Using mathematical models, we found adding protein turnover to a PTM cycle diminishes both sensitivity and ultrasensitivity, particularly in models admitting long ubiquitin chains. We also found that proteins sharing an E3 ligase can indeed have coupled changes in both expression and sensitivity to signals. These results imply that accounting for crosstalk in protein degradation networks is crucial for the interpretation of results from a wide variety of common experimental perturbations to living systems.

Misuse of the Michaelis-Menten rate law for protein interaction networks and its remedy

For over a century, the Michaelis-Menten (MM) rate law has been used to describe the rates of enzyme-catalyzed reactions and gene expression. Despite the ubiquity of the MM rate law, it accurately captures the dynamics of underlying biochemical reactions only so long as it is applied under the right condition, namely, that the substrate is in large excess over the enzyme-substrate complex. Unfortunately, in circumstances where its validity condition is not satisfied, especially so in protein interaction networks, the MM rate law has frequently been misused. In this review, we illustrate how inappropriate use of the MM rate law distorts the dynamics of the system, provides mistaken estimates of parameter values, and makes false predictions of dynamical features such as ultrasensitivity, bistability, and oscillations. We describe how these problems can be resolved with a slightly modified form of the MM rate law, based on the total quasi-steady state approximation (tQSSA). Furthermore, we show that the tQSSA can be used for accurate stochastic simulations at a lower computational cost than using the full set of mass-action rate laws. This review describes how to use quasi-steady state approximations in the right context, to prevent drawing erroneous conclusions from in silico simulations.

Bayesian modeling reveals metabolite-dependent ultrasensitivity in the cyanobacterial circadian clock

Mathematical models can enable a predictive understanding of mechanism in cell biology by quantitatively describing complex networks of interactions, but such models are often poorly constrained by available data. Owing to its relative biochemical simplicity, the core circadian oscillator in Synechococcus elongatus has become a prototypical system for studying how collective dynamics emerge from molecular interactions. The oscillator consists of only three proteins, KaiA, KaiB, and KaiC, and near-24-h cycles of KaiC phosphorylation can be reconstituted in vitro. Here, we formulate a molecularly detailed but mechanistically naive model of the KaiA-KaiC subsystem and fit it directly to experimental data within a Bayesian parameter estimation framework. Analysis of the fits consistently reveals an ultrasensitive response for KaiC phosphorylation as a function of KaiA concentration, which we confirm experimentally. This ultrasensitivity primarily results from the differential affinity of KaiA for competing nucleotide-bound states of KaiC. We argue that the ultrasensitive stimulus-response relation likely plays an important role in metabolic compensation by suppressing premature phosphorylation at nighttime.

Nonlinear Regulation of Enzyme-Free DNA Circuitry with Ultrasensitive Switches

DNA is used to construct synthetic chemical reaction networks (CRNs), such as inorganic oscillators and gene regulatory networks. Nonlinear regulation with a simpler molecular mechanism is particularly important in large-scale CRNs with complex dynamics, such as bistability, adaptation, and oscillation of cellular functions. Here we introduce a new approach based on ultrasensitive switches as modular regulatory elements to nonlinearly regulate DNA-based CRNs. The nonlinear behavior of the systems can be finely tuned by programmable regulation of the linker length and the ligand binding sites, of which the Hill coefficients (n_H) are in the range of 1.00-2.32. By integrating two different strand displacement reactions with low-order nonlinearities (n_H ≈ 1.44 and 1.54), we could construct CRNs exhibiting high-order nonlinearities with Hill coefficients of up to ∼2.70. In addition, this could provide an efficient approach for designing CRNs at will with complex chemical dynamics by incorporating our design with previously developed enzyme-free DNA circuits.

The Design Principles of Biochemical Timers: Circuits that Discriminate between Transient and Sustained Stimulation

Many cellular responses for which timing is critical display temporal filtering-the ability to suppress response until stimulated for longer than a given minimal time. To identify biochemical circuits capable of kinetic filtering, we comprehensively searched the space of three-node enzymatic networks. We define a metric of "temporal ultrasensitivity," the steepness of activation as a function of stimulus duration. We identified five classes of core network motifs capable of temporal filtering, each with distinct functional properties such as rejecting high-frequency noise, committing to response (bistability), and distinguishing between long stimuli. Combinations of the two most robust motifs, double inhibition (DI) and positive feedback with AND logic (PF_AND), underlie several natural timer circuits involved in processes such as cell cycle transitions, T cell activation, and departure from the pluripotent state. The biochemical network motifs described in this study form a basis for understanding common ways cells make dynamic decisions.

Operating regimes in a single enzymatic cascade at ensemble-level

Single enzymatic cascade, ubiquitously found in cellular signaling networks, is a phosphorylation-dephosphorylation reaction cycle causing a transition between inactive and active states of a protein catalysed by kinase and phosphatase, respectively. Steady-state information processing ability of such a cycle (e.g., MAPK cascade) has been classified into four qualitatively different operating regimes, viz., hyperbolic (H), signal-transducing (ST), threshold-hyperbolic (TH) and ultrasensitive (U). These four regimes represent qualitatively different dose-response curves, that is, relationship between concentrations of input kinase (e.g., pMEK) and response activated protein (e.g., pERK). Regimes were identified using a deterministic model accounting for population-averaged behavior only. Operating regimes can be strongly influenced by the inherently present cell-to-cell variability in an ensemble of cells which is captured in the form of pMEK and pERK distributions using reporter-based single-cell experimentation. In this study, we show that such experimentally acquired snapshot pMEK and pERK distribution data of a single MAPK cascade can be directly used to infer the underlying operating regime even in the absence of a dose-response curve. This deduction is possible primarily due to the presence of a monotonic relationship between experimental observables R_IQR, ratio of the inter-quartile range of the pERK and pMEK distribution pairs and R_M, ratio of the medians of the distribution pair. We demonstrate this relationship by systematic analysis of a quasi-steady state approximated model superimposed with an input gamma distribution constrained by the stimulus strength specific pMEK distribution measured on Jurkat-T cells stimulated with PMA. As a first, we show that introduction of cell-to-cell variability only in the upstream kinase achieved by superimposition of an appropriate input pMEK distribution on the dose-response curve can predict bimodal response pERK distribution in ST regime. Implementation of the proposed method on the input-response distribution pair obtained in stimulated Jurkat-T cells revealed that while low-dosage PMA stimulation preserves the H regime observed in resting cells, high-dosage causes H to ST regime transition.

Monostationarity and Multistationarity in Tree Networks of Goldbeter-Koshland Loops

A major challenge in systems biology is to elicit general properties in the face of molecular complexity. Here, we introduce a class of enzyme-catalysed biochemical networks and examine how the existence of a single positive steady state (monostationarity) depends on the network structure, enzyme mechanisms, kinetic rate laws and parameter values. We consider Goldbeter-Koshland (GK) covalent modification loops arranged in a tree network, so that a substrate form in one loop can be an enzyme in another loop. GK loops are a canonical motif in cell signalling and trees offer a generalisation of linear cascades which accommodate network complexity while remaining mathematically tractable. In particular, they permit a modular, recursive proof strategy which may be more widely applicable. We show that if each enzyme follows its own complex reaction mechanism under mass action kinetics, then any network is monostationary for all appropriate parameter values. If the kinetics is non-mass action with a plausible monotonicity requirement, and each enzyme follows the Michaelis-Menten mechanism, then monostationarity is preserved. Surprisingly, a single GK loop with a complex enzyme mechanism under non-mass action monotone kinetics can have more than one positive steady state (multistationarity). The broader interplay between network structure, enzyme mechanism and kinetics remains an intriguing open problem.

Protein Serine/Threonine Phosphatases: Keys to Unlocking Regulators and Substrates

Protein serine/threonine phosphatases (PPPs) are ancient enzymes, with distinct types conserved across eukaryotic evolution. PPPs are segregated into types primarily on the basis of the unique interactions of PPP catalytic subunits with regulatory proteins. The resulting holoenzymes dock substrates distal to the active site to enhance specificity. This review focuses on the subunit and substrate interactions for PPP that depend on short linear motifs. Insights about these motifs from structures of holoenzymes open new opportunities for computational biology approaches to elucidate PPP networks. There is an expanding knowledge base of posttranslational modifications of PPP catalytic and regulatory subunits, as well as of their substrates, including phosphorylation, acetylation, and ubiquitination. Cross talk between these posttranslational modifications creates PPP-based signaling. Knowledge of PPP complexes, signaling clusters, as well as how PPPs communicate with each other in response to cellular signals should unlock the doors to PPP networks and signaling "clouds" that orchestrate and coordinate different aspects of cell physiology.

Properties of cell signaling pathways and gene expression systems operating far from steady-state

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are basic units of signaling systems and their steady-state properties are well understood. However, the behavior of such systems before steady-state is poorly characterized. Here, we analyzed the properties of input-output curves for each of these systems as they approach steady-state. In ligand-receptor systems, the EC₅₀ (concentration of the ligand that occupies 50% of the receptors) is higher before the system reaches steady-state. Based on this behavior, we have previously defined PRESS (for pre-equilibrium sensing and signaling), a general “systems level” mechanism cells may use to overcome input saturation. Originally, we showed that, given a step stimulation, PRESS operates when the kinetics of ligand-receptor binding are slower than the downstream signaling steps. Now, we show that, provided the input increases slowly, it is not essential for the ligand binding reaction itself to be slow. In addition, we demonstrate that covalent modification cycles and gene expression systems may also operate in PRESS mode. Thus, nearly all biochemical processes may operate in PRESS mode, suggesting that this mechanism may be ubiquitous in cell signaling systems.

Intrinsic limits of information transmission in biochemical signalling motifs

All living things have evolved to sense changes in their environment in order to respond in adaptive ways. At the cellular level, these sensing systems generally involve receptor molecules at the cell surface, which detect changes outside the cell and relay those changes to the appropriate response elements downstream. With the advent of experimental technologies that can track signalling at the single-cell level, it has become clear that many signalling systems exhibit significant levels of 'noise,' manifesting as differential responses of otherwise identical cells to the same environment. This noise has a large impact on the capacity of cell signalling networks to transmit information from the environment. Application of information theory to experimental data has found that all systems studied to date encode less than 2.5 bits of information, with the majority transmitting significantly less than 1 bit. Given the growing interest in applying information theory to biological data, it is crucial to understand whether the low values observed to date represent some sort of intrinsic limit on information flow given the inherently stochastic nature of biochemical signalling events. In this work, we used a series of computational models to explore how much information a variety of common 'signalling motifs' can encode. We found that the majority of these motifs, which serve as the basic building blocks of cell signalling networks, can encode far more information (4-6 bits) than has ever been observed experimentally. In addition to providing a consistent framework for estimating information-theoretic quantities from experimental data, our findings suggest that the low levels of information flow observed so far in living system are not necessarily due to intrinsic limitations. Further experimental work will be needed to understand whether certain cell signalling systems actually can approach the intrinsic limits described here, and to understand the sources and purpose of the variation that reduces information flow in living cells.

When More Is Less: Dual Phosphorylation Protects Signaling Off State against Overexpression

Kinases in signaling pathways are commonly activated by multisite phosphorylation. For example, the mitogen-activated protein kinase Erk is activated by its kinase Mek by two consecutive phosphorylations within its activation loop. In this article, we use kinetic models to study how the activation of Erk is coupled to its abundance. Intuitively, Erk activity should rise with increasing amounts of Erk protein. However, a mathematical model shows that the signaling off state is robust to increasing amounts of Erk, and Erk activity may even decline with increasing amounts of Erk. This counterintuitive, bell-shaped response of Erk activity to increasing amounts of Erk arises from the competition of the unmodified and single phosphorylated form of Erk for access to its kinase Mek. This shows that phosphorylation cycles can contain an intrinsic robustness mechanism that protects signaling from aberrant activation e.g., by gene expression noise or kinase overexpression after gene duplication events in diseases like cancer.

Computing Mathematical Functions using DNA via Fractional Coding

This paper discusses the implementation of mathematical functions such as exponentials, trigonometric functions, the sigmoid function and the perceptron function with molecular reactions in general, and DNA strand displacement reactions in particular. The molecular constructs for these functions are predicated on a novel representation for input and output values: a fractional encoding, in which values are represented by the relative concentrations of two molecular types, denoted as type-1 and type-0. This representation is inspired by a technique from digital electronic design, termed stochastic logic, in which values are represented by the probability of 1's in a stream of randomly generated 0's and 1's. Research in the electronic realm has shown that a variety of complex functions can be computed with remarkably simple circuitry with this stochastic approach. This paper demonstrates how stochastic electronic designs can be translated to molecular circuits. It presents molecular implementations of mathematical functions that are considerably more complex than any shown to date. All designs are validated using mass-action simulations of the chemical kinetics of DNA strand displacement reactions.

Exploring the Function of Dynamic Phosphorylation-Dephosphorylation Cycles

Protein phosphorylation is a dynamic post-translational modification critical for biological responses. At the level of individual molecules, phosphorylation dynamics can have important functional implications, but this information is rarely quantified. We discuss how rapid phosphorylation-dephosphorylation cycles could underlie important signaling properties, including the ability to rapidly bind and release proteins.

Analysis of network motifs in cellular regulation: Structural similarities, input-output relations and signal integration

Much of the complexity of regulatory networks derives from the necessity to integrate multiple signals and to avoid malfunction due to cross-talk or harmful perturbations. Hence, one may expect that the input-output behavior of larger networks is not necessarily more complex than that of smaller network motifs which suggests that both can, under certain conditions, be described by similar equations. In this review, we illustrate this approach by discussing the similarities that exist in the steady state descriptions of a simple bimolecular reaction, covalent modification cycles and bacterial two-component systems. Interestingly, in all three systems fundamental input-output characteristics such as thresholds, ultrasensitivity or concentration robustness are described by structurally similar equations. Depending on the system the meaning of the parameters can differ ranging from protein concentrations and affinity constants to complex parameter combinations which allows for a quantitative understanding of signal integration in these systems. We argue that this approach may also be extended to larger regulatory networks.

Signaling Architectures that Transmit Unidirectional Information Despite Retroactivity

A signaling pathway transmits information from an upstream system to downstream systems, ideally in a unidirectional fashion. A key obstacle to unidirectional transmission is retroactivity, the additional reaction flux that affects a system once its species interact with those of downstream systems. This raises the fundamental question of whether signaling pathways have developed specialized architectures that overcome retroactivity and transmit unidirectional signals. Here, we propose a general procedure based on mathematical analysis that provides an answer to this question. Using this procedure, we analyze the ability of a variety of signaling architectures to transmit one-way (from upstream to downstream) signals, as key biological parameters are tuned. We find that single stage phosphorylation and phosphotransfer systems that transmit signals from a kinase show a stringent design tradeoff that hampers their ability to overcome retroactivity. Interestingly, cascades of these architectures, which are highly represented in nature, can overcome this tradeoff and thus enable unidirectional transmission. By contrast, phosphotransfer systems, and single and double phosphorylation cycles that transmit signals from a substrate, are unable to mitigate retroactivity effects, even when cascaded, and hence are not well suited for unidirectional information transmission. These results are largely independent of the specific reaction-rate constant values, and depend on the topology of the architectures. Our results therefore identify signaling architectures that, allowing unidirectional transmission of signals, embody modular processes that conserve their input/output behavior across multiple contexts. These findings can be used to decompose natural signal transduction networks into modules, and at the same time, they establish a library of devices that can be used in synthetic biology to facilitate modular circuit design.

Operating regimes of covalent modification cycles at high enzyme concentrations

The Goldbeter-Koshland model has been a paradigm for ultrasensitivity in biological networks for more than 30 years. Despite its simplicity the validity of this model is restricted to conditions when the substrate is in excess over the converter enzymes - a condition that is easy to satisfy in vitro, but which is rarely satisfied in vivo. Here, we analyze the Goldbeter-Koshland model by means of the total quasi-steady state approximation which yields a comprehensive classification of the steady state operating regimes under conditions when the enzyme concentrations are comparable to or larger than that of the substrate. Where possible we derive simple expressions characterizing the input-output behavior of the system. Our analysis suggests that enhanced sensitivity occurs if the concentration of at least one of the converter enzymes is smaller (but not necessarily much smaller) than that of the substrate and if that enzyme is saturated. Conversely, if both enzymes are saturated and at least one of the enzyme concentrations exceeds that of the substrate the system exhibits concentration robustness with respect to changes in that enzyme concentration. Also, depending on the enzyme's saturation degrees and the ratio between their maximal reaction rates the total fraction of phosphorylated substrate may increase, decrease or change nonmonotonically as a function of the total substrate concentration. The latter finding may aid the interpretation of experiments involving genetic perturbations of enzyme and substrate abundances.

Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations

Ultrasensitive response motifs, capable of converting graded stimuli into binary responses, are well-conserved in signal transduction networks. Although it has been shown that a cascade arrangement of multiple ultrasensitive modules can enhance the system's ultrasensitivity, how a given combination of layers affects a cascade's ultrasensitivity remains an open question for the general case. Here, we introduce a methodology that allows us to determine the presence of sequestration effects and to quantify the relative contribution of each module to the overall cascade's ultrasensitivity. The proposed analysis framework provides a natural link between global and local ultrasensitivity descriptors and it is particularly well-suited to characterize and understand mathematical models used to study real biological systems. As a case study, we have considered three mathematical models introduced by O'Shaughnessy et al. to study a tunable synthetic MAPK cascade, and we show how our methodology can help modelers better understand alternative models.

Biologically inspired design of feedback control systems implemented using DNA strand displacement reactions

The use of abstract chemical reaction networks (CRNs) as a modelling and design framework for the implementation of computing and control circuits using enzyme-free, entropy driven DNA strand displacement (DSD) reactions is starting to garner widespread attention in the area of synthetic biology. Previous work in this area has demonstrated the theoretical plausibility of using this approach to design biomolecular feedback control systems based on classical proportional-integral (PI) controllers, which may be constructed from CRNs implementing gain, summation and integrator operators. Here, we propose an alternative design approach that utilises the abstract chemical reactions involved in cellular signalling cycles to implement a biomolecular controller - termed a signalling-cycle (SC) controller. We compare the performance of the PI and SC controllers in closed-loop with a nonlinear second-order chemical process. Our results show that the SC controller outperforms the PI controller in terms of both performance and robustness, and also requires fewer abstract chemical reactions to implement, highlighting its potential usefulness in the construction of biomolecular control circuits.

On the role of topology in regulating transcriptional cascades

Topology of interactions in a transcriptional cascade determines the behavior of its signal-response profile and the activation states of genes.

Exploiting the dynamic properties of covalent modification cycle for the design of synthetic analog biomolecular circuitry

Background

Cycles of covalent modification are ubiquitous motifs in cellular signalling. Although such signalling cycles are implemented via a highly concise set of chemical reactions, they have been shown to be capable of producing multiple distinct input-output mapping behaviours – ultrasensitive, hyperbolic, signal-transducing and threshold-hyperbolic.

Results

In this paper, we show how the set of chemical reactions underlying covalent modification cycles can be exploited for the design of synthetic analog biomolecular circuitry. We show that biomolecular circuits based on the dynamics of covalent modification cycles allow (a) the computation of nonlinear operators using far fewer chemical reactions than purely abstract designs based on chemical reaction network theory, and (b) the design of nonlinear feedback controllers with strong performance and robustness properties.

Conclusions

Our designs provide a more efficient route for translation of complex circuits and systems from chemical reactions to DNA strand displacement-based chemistry, thus facilitating their experimental implementation in future Synthetic Biology applications.

Signaling cascades transmit information downstream and upstream but unlikely simultaneously

Background

Signal transduction is the process through which cells communicate with the external environment, interpret stimuli and respond to them. This mechanism is controlled by signaling cascades, which play the role of intracellular transmitter, being able to transmit biochemical information between cell membrane and nucleus. In theory as well as in practice, it has been shown that a perturbation can propagate upstream (and not only downstream) a cascade, by a mechanism known as retroactivity. This study aims to compare the conditions on biochemical parameters which favor one or the other direction of signaling in such a cascade.

Results

From a mathematical point of view, we show that the steady states of a cascade of arbitrary length n are described by an iterative map of second order, meaning that the cascade tiers are actually coupled three-by-three. We study the influence of the biochemical parameters in the control of the direction of transmission – upstream and/or downstream – along a signaling cascade. A numerical and statistical approach, based on the random scan of parameters describing a 3-tier signaling cascade, provides complementary findings to the analytical study. In particular, computing the likelihood of parameters with respect to various signaling regimes, we identify conditions on biochemical parameters which enhance a specific direction of propagation corresponding to forward or retro-signaling regimes. A compact graphical representation is designed to relay the gist of these conditions.

Conclusions

The values of biochemical parameters such as kinetic rates, Michaelis-Menten constants, total concentrations of kinases and of phosphatases, determine the propensity of a cascade to favor or impede downstream or upstream signal transmission. We found that generally there is an opposition between parameter sets favoring forward and retro-signaling regimes. Therefore, on one hand our study supports the idea that in most cases, retroactive effects can be neglected when a cascade which is efficient in forward signaling, is perturbed by an external ligand inhibiting the activation at some tier of the cascade. This result is relevant for therapeutic methodologies based on kinase inhibition. On the other hand, our study highlights a less-known part of the parameter space where, although the forward signaling is inefficient, the cascade can interestingly act as a retro-signaling device.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0303-2) contains supplementary material, which is available to authorized users.

Ultrasensitive Negative Feedback Control: A Natural Approach for the Design of Synthetic Controllers

Many of the most important potential applications of Synthetic Biology will require the ability to design and implement high performance feedback control systems that can accurately regulate the dynamics of multiple molecular species within the cell. Here, we argue that the use of design strategies based on combining ultrasensitive response dynamics with negative feedback represents a natural approach to this problem that fully exploits the strongly nonlinear nature of cellular information processing. We propose that such feedback mechanisms can explain the adaptive responses observed in one of the most widely studied biomolecular feedback systems—the yeast osmoregulatory response network. Based on our analysis of such system, we identify strong links with a well-known branch of mathematical systems theory from the field of Control Engineering, known as Sliding Mode Control. These insights allow us to develop design guidelines that can inform the construction of feedback controllers for synthetic biological systems.

The Dark Side of Cell Signaling: Positive Roles for Negative Regulators

Cell signaling is dominated by analyzing positive responses to stimuli. Signal activation is balanced by negative regulators that are generally considered to terminate signaling. Rather than exerting only negative effects, however, many such regulators play important roles in enhancing cell-signaling control. Considering responses downstream of selected cell-surface receptors, we discuss how receptor internalization affects signaling specificity and how rapid kinase/phosphatase and GTP/GDP cycles increase responsiveness and allow kinetic proofreading in receptor signaling. We highlight the blurring of distinctions between positive and negative signals, recasting signal termination as the response to a switch-like transition into a new cellular state.

Implementing Nonlinear Feedback Controllers Using DNA Strand Displacement Reactions

We show how an important class of nonlinear feedback controllers can be designed using idealized abstract chemical reactions and implemented via DNA strand displacement (DSD) reactions. Exploiting chemical reaction networks (CRNs) as a programming language for the design of complex circuits and networks, we show how a set of unimolecular and bimolecular reactions can be used to realize input-output dynamics that produce a nonlinear quasi sliding mode (QSM) feedback controller. The kinetics of the required chemical reactions can then be implemented as enzyme-free, enthalpy/entropy driven DNA reactions using a toehold mediated strand displacement mechanism via Watson-Crick base pairing and branch migration. We demonstrate that the closed loop response of the nonlinear QSM controller outperforms a traditional linear controller by facilitating much faster tracking response dynamics without introducing overshoots in the transient response. The resulting controller is highly modular and is less affected by retroactivity effects than standard linear designs.

Phosphatase specificity and pathway insulation in signaling networks

Phosphatases play an important role in cellular signaling networks by regulating the phosphorylation state of proteins. Phosphatases are classically considered to be promiscuous, acting on tens to hundreds of different substrates. We recently demonstrated that a shared phosphatase can couple the responses of two proteins to incoming signals, even if those two substrates are from otherwise isolated areas of the network. This finding raises a potential paradox: if phosphatases are indeed highly promiscuous, how do cells insulate themselves against unwanted crosstalk? Here, we use mathematical models to explore three possible insulation mechanisms. One approach involves evolving phosphatase KM values that are large enough to prevent saturation by the phosphatase's substrates. Although this is an effective method for generating isolation, the phosphatase becomes a highly inefficient enzyme, which prevents the system from achieving switch-like responses and can result in slow response kinetics. We also explore the idea that substrate degradation can serve as an effective phosphatase. Assuming that degradation is unsaturatable, this mechanism could insulate substrates from crosstalk, but it would also preclude ultrasensitive responses and would require very high substrate turnover to achieve rapid dephosphorylation kinetics. Finally, we show that adaptor subunits, such as those found on phosphatases like PP2A, can provide effective insulation against phosphatase crosstalk, but only if their binding to substrates is uncoupled from their binding to the catalytic core. Analysis of the interaction network of PP2A's adaptor domains reveals that although its adaptors may isolate subsets of targets from one another, there is still a strong potential for phosphatase crosstalk within those subsets. Understanding how phosphatase crosstalk and the insulation mechanisms described here impact the function and evolution of signaling networks represents a major challenge for experimental and computational systems biology.

Simple biochemical pathways far from steady state can provide switchlike and integrated responses

Covalent modification cycles (systems in which the activity of a substrate is regulated by the action of two opposing enzymes) and ligand/receptor interactions are ubiquitous in signaling systems and their steady-state properties are well understood. However, the behavior of such systems far from steady state remains unclear. Here, we analyze the properties of covalent modification cycles and ligand/receptor interactions driven by the accumulation of the activating enzyme and the ligand, respectively. We show that for a large range of parameters both systems produce sharp switchlike response and yet allow for temporal integration of the signal, two desirable signaling properties. Ultrasensitivity is observed also in a region of parameters where the steady-state response is hyperbolic. The temporal integration properties are tunable by regulating the levels of the deactivating enzyme and receptor, as well as by adjusting the rate of accumulation of the activating enzyme and ligand. We propose that this tunability is used to generate precise responses in signaling systems.

Digital microfluidic immunocytochemistry in single cells

We report a new technique called Digital microfluidic Immunocytochemistry in Single Cells (DISC). DISC automates protocols for cell culture, stimulation and immunocytochemistry, enabling the interrogation of protein phosphorylation on pulsing with stimulus for as little as 3 s. DISC was used to probe the phosphorylation states of platelet-derived growth factor receptor (PDGFR) and the downstream signalling protein, Akt, to evaluate concentration- and time-dependent effects of stimulation. The high time resolution of the technique allowed for surprising new observations-for example, a 10 s pulse stimulus of a low concentration of PDGF is sufficient to cause >30% of adherent fibroblasts to commit to Akt activation. With the ability to quantitatively probe signalling events with high time resolution at the single-cell level, we propose that DISC may be an important new technique for a wide range of applications, especially for screening signalling responses of a heterogeneous cell population.

Tandem phosphorylation within an intrinsically disordered region regulates ACTN4 function

Phosphorylated residues occur preferentially in the intrinsically disordered regions of eukaryotic proteins. In the disordered amino-terminal region of human α-actinin-4 (ACTN4), Tyr(4) and Tyr(31) are phosphorylated in cells stimulated with epidermal growth factor (EGF), and a mutant with phosphorylation-mimicking mutations of both tyrosines exhibits reduced interaction with actin in vitro. Cleavage of ACTN4 by m-calpain, a protease that in motile cells is predominantly activated at the rear, removes the Tyr(4) site. We found that introducing a phosphomimetic mutation at only Tyr(31) was sufficient to inhibit the interaction with actin in vitro. However, molecular dynamics simulations predicted that Tyr(31) is mostly buried and that phosphorylation of Tyr(4) would increase the solvent exposure and thus kinase accessibility of Tyr(31). In fibroblast cells, EGF stimulation increased tyrosine phosphorylation of a mutant form of ACTN4 with a phosphorylation-mimicking residue at Tyr(4), whereas a truncated mutant representing the product of m-calpain cleavage exhibited EGF-stimulated tyrosine phosphorylation at a background amount similar to that observed for a double phosphomimetic mutant of Tyr(4) and Tyr(31). We also found that inhibition of the receptor tyrosine kinases of the TAM family, such as AXL, blocked EGF-stimulated tyrosine phosphorylation of ACTN4. Mathematical modeling predicted that the kinetics of phosphorylation at Tyr(31) can be dictated by the kinase affinity for Tyr(4). This study suggests that tandem-site phosphorylation within intrinsically disordered regions provides a mechanism for a site to function as a switch to reveal a nearby function-regulating site.

Quantitative Analysis of Robustness of Dynamic Response and Signal Transfer in Insulin mediated PI3K/AKT Pathway

Robustness is a critical feature of signaling pathways ensuring signal propagation with high fidelity in the event of perturbations. Here we present a detailed quantitative analysis of robustness in insulin mediated PI3K/AKT pathway, a critical signaling pathway maintaining self-renewal in human embryonic stem cells. Using global sensitivity analysis, we identified robustness promoting mechanisms that ensure (1) maintenance of a first order or overshoot dynamics of self-renewal molecule, p-AKT and (2) robust transfer of signals from oscillatory insulin stimulus to p-AKT in the presence of noise. Our results indicate that negative feedback controls the robustness to most perturbations. Faithful transfer of signal from the stimulating ligand to p-AKT occurs even in the presence of noise, albeit with signal attenuation and high frequency cut-off. Negative feedback contributes to signal attenuation, while positive regulators upstream of PIP₃ contribute to signal amplification. These results establish precise mechanisms to modulate self-renewal molecules like p-AKT.

Three level signal transduction cascades lead to reliably timed switches

Signaling cascades proliferate signals received on the cell membrane to the nucleus. While noise filtering, ultra-sensitive switches, and signal amplification have all been shown to be features of such signaling cascades, it is not understood why cascades typically show three or four layers. Using singular perturbation theory, Michaelis-Menten type equations are derived for open enzymatic systems. Cascading these equations we demonstrate that the output signal as a function of time becomes sigmoidal with the addition of more layers. Furthermore, it is shown that the activation time will speed up to a point, after which more layers become superfluous. It is shown that three layers create a reliable sigmoidal response progress curve from a wide variety of time-dependent signaling inputs arriving at the cell membrane, suggesting the evolutionary benefit of the observed cascades.

Impact of upstream and downstream constraints on a signaling module's ultrasensitivity

Much work has been done on the study of the biochemical mechanisms that result in ultrasensitive behavior of simple biochemical modules. However, in a living cell, such modules are embedded in a bigger network that constrains the range of inputs that the module will receive as well as the range of the module's outputs that network will be able to detect. Here, we studied how the effective ultrasensitivity of a modular system is affected by these restrictions. We use a simple setup to explore to what extent the dynamic range spanned by upstream and downstream components of an ultrasensitive module impact on the effective sensitivity of the system. Interestingly, we found for some ultrasensitive motifs that dynamic range limitations imposed by downstream components can produce effective sensitivities much larger than that of the original module when considered in isolation.

Ultrasensitivity part I: Michaelian responses and zero-order ultrasensitivity

Quantitative studies of signal transduction systems have shown that ultrasensitive responses - switch-like, sigmoidal input/output relationships - are commonplace in cell signaling. Ultrasensitivity is important for various complex signaling systems, including signaling cascades, bistable switches, and oscillators. In this first installment of a series on ultrasensitivity we survey the occurrence of ultrasensitive responses in signaling systems. We review why the simplest mass action systems exhibit Michaelian responses, and then move on to zero-order ultrasensitivity, a phenomenon that occurs when signaling proteins are operating near saturation. We also discuss the physiological relevance of zero-order ultrasensitivity to cellular regulation.

Reciprocal regulation as a source of ultrasensitivity in two-component systems with a bifunctional sensor kinase

Two-component signal transduction systems, where the phosphorylation state of a regulator protein is modulated by a sensor kinase, are common in bacteria and other microbes. In many of these systems, the sensor kinase is bifunctional catalyzing both, the phosphorylation and the dephosphorylation of the regulator protein in response to input signals. Previous studies have shown that systems with a bifunctional enzyme can adjust the phosphorylation level of the regulator protein independently of the total protein concentrations – a property known as concentration robustness. Here, I argue that two-component systems with a bifunctional enzyme may also exhibit ultrasensitivity if the input signal reciprocally affects multiple activities of the sensor kinase. To this end, I consider the case where an allosteric effector inhibits autophosphorylation and, concomitantly, activates the enzyme's phosphatase activity, as observed experimentally in the PhoQ/PhoP and NRII/NRI systems. A theoretical analysis reveals two operating regimes under steady state conditions depending on the effector affinity: If the affinity is low the system produces a graded response with respect to input signals and exhibits stimulus-dependent concentration robustness – consistent with previous experiments. In contrast, a high-affinity effector may generate ultrasensitivity by a similar mechanism as phosphorylation-dephosphorylation cycles with distinct converter enzymes. The occurrence of ultrasensitivity requires saturation of the sensor kinase's phosphatase activity, but is restricted to low effector concentrations, which suggests that this mode of operation might be employed for the detection and amplification of low abundant input signals. Interestingly, the same mechanism also applies to covalent modification cycles with a bifunctional converter enzyme, which suggests that reciprocal regulation, as a mechanism to generate ultrasensitivity, is not restricted to two-component systems, but may apply more generally to bifunctional enzyme systems.

Bacteria often use two-component systems to sense and respond to environmental changes, which involves autophosphorylation of a sensor kinase and phosphotransfer to a cognate response regulator. However, despite conservation of this ‘classical’ scheme there exist substantial variations in the mechanism of phosphotransfer among systems. Also, many sensor kinases exhibit phosphatase activity raising the question whether such a bifunctional architecture enables special regulatory properties in the response behavior to input signals. According to previous studies, classical two-component systems are unlikely to produce sigmoidal response curves (ultrasensitivity) if the sensor protein is bifunctional. Here, I argue that this is not necessarily true if the input stimulus (allosteric effector) reciprocally affects multiple activities of the sensor kinase, as it seems to be common for bifunctional enzymes. To this end, I propose and analyze an extension of the experimentally well-supported Batchelor-Goulian model which shows that ultrasensitivity requires a high-affinity effector and saturation of the phosphatase activity. The underlying mechanism involves sequestration of the effector by the sensor kinase which restricts the occurrence of ultrasensitivity to sufficiently low effector concentrations. Hence, this operating regime might be useful to sense effector limitations or to amplify weak input signals.

A fundamental trade-off in covalent switching and its circumvention by enzyme bifunctionality in glucose homeostasis

Covalent modification provides a mechanism for modulating molecular state and regulating physiology. A cycle of competing enzymes that add and remove a single modification can act as a molecular switch between "on" and "off" and has been widely studied as a core motif in systems biology. Here, we exploit the recently developed "linear framework" for time scale separation to determine the general principles of such switches. These methods are not limited to Michaelis-Menten assumptions, and our conclusions hold for enzymes whose mechanisms may be arbitrarily complicated. We show that switching efficiency improves with increasing irreversibility of the enzymes and that the on/off transition occurs when the ratio of enzyme levels reaches a value that depends only on the rate constants. Fluctuations in enzyme levels, which habitually occur due to cellular heterogeneity, can cause flipping back and forth between on and off, leading to incoherent mosaic behavior in tissues, that worsens as switching becomes sharper. This trade-off can be circumvented if enzyme levels are correlated. In particular, if the competing catalytic domains are on the same protein but do not influence each other, the resulting bifunctional enzyme can switch sharply while remaining coherent. In the mammalian liver, the switch between glycolysis and gluconeogenesis is regulated by the bifunctional 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2). We suggest that bifunctionality of PFK-2/FBPase-2 complements the metabolic zonation of the liver by ensuring coherent switching in response to insulin and glucagon.

Tunable signal processing through a kinase control cycle: the IKK signaling node

The transcription factor NFκB, a key component of the immune system, shows intricate stimulus-specific temporal dynamics. Those dynamics are thought to play a role in controlling the physiological response to cytokines and pathogens. Biochemical evidence suggests that the NFκB inducing kinase, IKK, a signaling hub onto which many signaling pathways converge, is regulated via a regulatory cycle comprising a poised, an active, and an inactive state. We hypothesize that it operates as a modulator of signal dynamics, actively reshaping the signals generated at the receptor proximal level. Here we show that a regulatory cycle can function in at least three dynamical regimes, tunable by regulating a single kinetic parameter. In particular, the simplest three-state regulatory cycle can generate signals with two well-defined phases, each with distinct coding capabilities in terms of the information they can carry about the stimulus. We also demonstrate that such a kinase cycle can function as a signal categorizer classifying diverse incoming signals into outputs with a limited set of temporal activity profiles. Finally, we discuss the extension of the results to other regulatory motifs that could be understood in terms of the regimes of the three-state cycle.

Zero-order ultrasensitivity: a study of criticality and fluctuations under the total quasi-steady state approximation in the linear noise regime

Zero-order ultrasensitivity (ZOU) is a long known and interesting phenomenon in enzyme networks. Here, a substrate is reversibly modified by two antagonistic enzymes (a 'push-pull' system) and the fraction in modified state undergoes a sharp switching from near-zero to near-unity at a critical value of the ratio of the enzyme concentrations, under saturation conditions. ZOU and its extensions have been studied for several decades now, ever since the seminal paper of Goldbeter and Koshland (1981); however, a complete probabilistic treatment, important for the study of fluctuations in finite populations, is still lacking. In this paper, we study ZOU using a modular approach, akin to the total quasi-steady state approximation (tQSSA). This approach leads to a set of Fokker-Planck (drift-diffusion) equations for the probability distributions of the intermediate enzyme-bound complexes, as well as the modified/unmodified fractions of substrate molecules. We obtain explicit expressions for various average fractions and their fluctuations in the linear noise approximation (LNA). The emergence of a 'critical point' for the switching transition is rigorously established. New analytical results are derived for the average and variance of the fractional substrate concentration in various chemical states in the near-critical regime. For the total fraction in the modified state, the variance is shown to be a maximum near the critical point and decays algebraically away from it, similar to a second-order phase transition. The new analytical results are compared with existing ones as well as detailed numerical simulations using a Gillespie algorithm.

Characterization of tradeoffs in biomolecular signaling

Systems-level tradeoffs are fundamental in engineering, and recent work has highlighted an analogous role for them in biology. However, the extent of validity of these tradeoffs, especially for biomolecular systems, is generally unclear. Here, we address this issue for signaling tradeoffs that can constrain, for a fixed concentration of the signaling protein, a simultaneous enhancement of the gain and range of an amplifier or of the gain and threshold of a switch. We find that these gain-related tradeoffs persist in mathematical models of biomolecular reaction mechanisms that are at the core of large classes of signaling systems. Further, we find that these tradeoffs are also prevalent in the parametric functional forms commonly used to describe input-output curves in experimental analyses. Finally, we find that these tradeoffs can persist even in the presence of transcriptional feedback mechanisms that can change the concentration of the signaling protein. These results present a systematic characterization of these tradeoffs in biomolecular signaling systems.