Push-Pull and Feedback Mechanisms Can Align Signaling System Outputs with Inputs

Negative feedback systems for modelling NF-κB transcription factor oscillatory activity

Low-dimensional negative feedback systems (NFSs) were developed within a signal flow model to describe the oscillatory activities of NF-κB caused by interactions with its inhibitor IκBα. The NFSs were established as 3rd- and 4th-order linear systems containing unperturbed and perturbed negative feedback (NF) loops with constant or time-varying NF strengths and a feed-forward loop. NF-related analytical solutions to the NFSs representing the time courses of NF-κB and IκBα were determined and their exact mathematical relationship was found. The NFS's parameters were determined to fit the experimental time courses of NF-κB in TNF-α-stimulated embryonic fibroblasts, rela-/- embryonic fibroblasts reconstituted with RelA, C9L cells, GFP-p65 knock-in embryonic fibroblasts and embryogenic fibroblasts lacking Iκβ and IκBε, LPS-stimulated IC-21 macrophages treated or not with DCPA, and anti-IgM-stimulated DT40 B-lymphocytes. The unperturbed and perturbed NFSs describing the above biosystems generated isochronous and non-isochronous solutions, depending on a constant or time-varying NF strength, respectively. The oscillation period of the NF-coupled solutions, the phase difference between them and the time delays in the appearance of cytoplasmic IκBα after stimulation of NF-κB were determined. A significant divergence between the IκBα solutions to the NFSs and the IκBα experimental courses led to a rejection of the NF coupling between NF-κB and IκBα in the above biosystems. It was shown that neither the linearity nor the low dimensionality of the NFSs altered the NF relationship and the divergence between the IκBα solutions to the NFS and IκBα experimental time courses. Although the NF relationship between IκBα and NF-κB was not confirmed in all the experimental data analyzed, delayed negative feedback was found in some cases.

Signal integration and integral feedback control with biochemical reaction networks

Biochemical reaction networks perform a variety of signal processing functions, one of which is computing the integrals of signal values. This is often used in integral feedback control, where it enables a system's output to respond to changing inputs, but to then return exactly back to some pre-determined setpoint value afterward. To gain a deeper understanding of how biochemical networks are able to both integrate signals and perform integral feedback control, we investigated these abilities for several simple reaction networks. We found imperfect overlap between these categories, with some networks able to perform both tasks, some able to perform integration but not integral feedback control, and some the other way around. Nevertheless, networks that could either integrate or perform integral feedback control shared key elements. In particular, they included a chemical species that was neutrally stable in the open loop system (no feedback), meaning that this species does not have a unique stable steady-state concentration. Neutral stability could arise from zeroth order decay reactions, binding to a partner that was produced at a constant rate (which occurs in antithetic control), or through a long chain of covalent cycles. Mathematically, it arose from rate equations for the reaction network that were underdetermined when evaluated at steady-state.

Design patterns of biological cells

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

A guide to ERK dynamics, part 2: downstream decoding

Signaling by the extracellular signal-regulated kinase (ERK) pathway controls many cellular processes, including cell division, death, and differentiation. In this second installment of a two-part review, we address the question of how the ERK pathway exerts distinct and context-specific effects on multiple processes. We discuss how the dynamics of ERK activity induce selective changes in gene expression programs, with insights from both experiments and computational models. With a focus on single-cell biosensor-based studies, we summarize four major functional modes for ERK signaling in tissues: adjusting the size of cell populations, gradient-based patterning, wave propagation of morphological changes, and diversification of cellular gene expression states. These modes of operation are disrupted in cancer and other related diseases and represent potential targets for therapeutic intervention. By understanding the dynamic mechanisms involved in ERK signaling, there is potential for pharmacological strategies that not only simply inhibit ERK, but also restore functional activity patterns and improve disease outcomes.

Design principles of improving the dose-response alignment in coupled GTPase switches

"Dose-response alignment" (DoRA), where the downstream response of cellular signaling pathways closely matches the fraction of activated receptor, can improve the fidelity of dose information transmission. The negative feedback has been experimentally identified as a key component for DoRA, but numerical simulations indicate that negative feedback is not sufficient to achieve perfect DoRA, i.e., perfect match of downstream response and receptor activation level. Thus a natural question is whether there exist design principles for signaling motifs within only negative feedback loops to improve DoRA to near-perfect DoRA. Here, we investigated several model formulations of an experimentally validated circuit that couples two molecular switches-mGTPase (monomeric GTPase) and tGTPase (heterotrimeric GTPases) - with negative feedback loops. In the absence of feedback, the low and intermediate mGTPase activation levels benefit DoRA in mass action and Hill-function models, respectively. Adding negative feedback has versatile roles on DoRA: it may impair DoRA in the mass action model with low mGTPase activation level and Hill-function model with intermediate mGTPase activation level; in other cases, i.e., the mass action model with a high mGTPase activation level or the Hill-function model with a non-intermediate mGTPase activation level, it improves DoRA. Furthermore, we found that DoRA in a longer cascade (i.e., tGTPase) can be obtained using Hill-function kinetics under certain conditions. In summary, we show how ranges of activity of mGTPase, reaction kinetics, the negative feedback, and the cascade length affect DoRA. This work provides a framework for improving the DoRA performance in signaling motifs with negative feedback.

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.

Real-Time Analysis of AKT Signaling Activities at Single-Cell Resolution Using Cyclic Peptide-Based Probes

Here we present a protocol for interrogating AKT signaling activities in living single cells, using a pair of cyclic peptide-based fluorescent probes. These probes are encapsulated in liposomes and delivered into cells, where they continuously report on AKT signaling activities through a Föster resonance energy transfer mechanism. We describe the use of a microwell chip to achieve single-cell resolution and demonstrate the procedure for on-chip immunostaining. Finally, we provide a method for data extraction, correction, and processing.

Signaling Pathways in Proton and Non-proton ASIC1a Activation

Acid-sensing ion channels (ASICs) regulate synaptic activities and play important roles in neurodegenerative diseases as well as pain conditions. Classically, ASICs are described as transiently activated by a reduced pH, followed by desensitization; the activation allows sodium influx, and in the case of ASIC1a-composed channels, also calcium to some degree. Several factors are emerging and extensively analyzed as modulators, activating, inhibiting, and potentiating specific channel subunits. However, the signaling pathways triggered by channel activation are only starting to be revealed.The channel has been recently shown to be activated through a mechanism other than proton-mediated. Indeed, the large extracellular loop of these channels opens the possibility that other non-proton ligands might exist. One such molecule discovered was a toxin present in the Texas coral snake venom. The finding was associated with the activation of the channel at neutral pH via the toxin and causing intense and unremitting pain.By using different pharmacological tools, we analyzed the downstream signaling pathway triggered either by the proton and non-proton activation for human, mouse, and rat ASIC1a-composed channels in in vitro models. We show that for all species analyzed, the non-protonic mode of activation determines the activation of the ERK signaling cascade at a higher level and duration compared to the proton mode.This study adds to the growing evidence of the important role ASIC1a channels play in different physiological and pathological conditions and also hints at a possible pathological mechanism for a sustained effect.

Oncogenic mutant RAS signaling activity is rescaled by the ERK/MAPK pathway

Activating mutations in RAS are present in ~ 30% of human tumors, and the resulting aberrations in ERK/MAPK signaling play a central role in oncogenesis. However, the form of these signaling changes is uncertain, with activating RAS mutants linked to both increased and decreased ERK activation in vivo. Rationally targeting the kinase activity of this pathway requires clarification of the quantitative effects of RAS mutations. Here, we use live‐cell imaging in cells expressing only one RAS isoform to quantify ERK activity with a new level of accuracy. We find that despite large differences in their biochemical activity, mutant KRAS isoforms within cells have similar ranges of ERK output. We identify roles for pathway‐level effects, including variation in feedback strength and feedforward modulation of phosphatase activity, that act to rescale pathway sensitivity, ultimately resisting changes in the dynamic range of ERK activity while preserving responsiveness to growth factor stimuli. Our results reconcile seemingly inconsistent reports within the literature and imply that the signaling changes induced by RAS mutations early in oncogenesis are subtle.

Digitonin-facilitated delivery of imaging probes enables single-cell analysis of AKT signalling activities in suspension cells

Analyzing intracellular signalling protein activities in living cells promises a better understanding of the signalling cascade and related biological processes. We have previously developed cyclic peptide-based probes for analyzing intracellular AKT signalling activities, but these peptide probes were not cell-permeable. Implementing fusogenic liposomes as delivery vehicles could circumvent the problem when analyzing adherent cells, but it remained challenging to study suspension cells using similar approaches. Here, we present a method for delivering these imaging probes into suspension cells using digitonin, which could transiently perforate the cell membrane. Using U87, THP-1, and Jurkat cells as model systems representing suspended adherent cells, myeloid cells, and lymphoid cells, we demonstrated that low concentrations of digitonin enabled a sufficient amount of probes to enter the cytosol without affecting cell viability. We further combined this delivery method with a microwell single-cell chip and interrogated the AKT signalling dynamics in THP-1 and Jurkat cells, followed by immunofluorescence-based quantitation of AKT expression levels. We resolved the cellular heterogeneity in AKT signalling activities and showed that the kinetic patterns of AKT signalling and the AKT expression levels were related in THP-1 cells, but decoupled in Jurkat cells. We expect that our approach can be adapted to study other suspension cells.

A precisely adjustable, variation-suppressed eukaryotic transcriptional controller to enable genetic discovery

Conditional expression of genes and observation of phenotype remain central to biological discovery. Current methods enable either on/off or imprecisely controlled graded gene expression. We developed a 'well-tempered' controller, WTC₈₄₆, for precisely adjustable, graded, growth condition independent expression of genes in Saccharomyces cerevisiae. Controlled genes are expressed from a strong semisynthetic promoter repressed by the prokaryotic TetR, which also represses its own synthesis; with basal expression abolished by a second, 'zeroing' repressor. The autorepression loop lowers cell-to-cell variation while enabling precise adjustment of protein expression by a chemical inducer. WTC₈₄₆ allelic strains in which the controller replaced the native promoters recapitulated known null phenotypes (CDC42, TPI1), exhibited novel overexpression phenotypes (IPL1), showed protein dosage-dependent growth rates and morphological phenotypes (CDC28, TOR2, PMA1 and the hitherto uncharacterized PBR1), and enabled cell cycle synchronization (CDC20). WTC₈₄₆ defines an 'expression clamp' allowing protein dosage to be adjusted by the experimenter across the range of cellular protein abundances, with limited variation around the setpoint.

A Circuit for Secretion-coupled Cellular Autonomy in Multicellular Eukaryotes.. [DOI: 10.1101/2021.03.18.436048]

Cancers represent complex autonomous systems, displaying self-sufficiency in growth signaling. Autonomous growth is fueled by a cancer cell’s ability to ‘secrete-and-sense’ growth factors: a poorly understood phenomenon. Using an integrated systems and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giαβγ and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control (CLC), allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimulus-proportionate secretion. Proteomic studies and protein-protein interaction network analyses pinpoint growth factors (e.g., the epidermal growth factor; EGF) as a key stimuli for such self-sustenance. Findings highlight how enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy of growth factors.

SYNOPSIS IMAGE

STANDFIRST TEXT

This work defines the inner workings of a Golgi-localized molecular circuitry comprised of coupled GTPases, which empowers cells to achieve self-sufficiency in growth factor signaling by creating a secrete-and-sense autocrine loop.

HIGHLIGHTS/MAIN FINDINGS

Modeling and experimental approaches were used to dissect a coupled GTPase circuit.

Coupling enables closed loop feedback and mutual control of GTPases.

Coupling generates dose response alignment behavior of sensing and secretion of growth factors.

Coupling is critical for multiscale feedback control to achieve secretion-coupled autonomy.

Molecular switch architecture determines response properties of signaling pathways

Many intracellular signaling pathways are composed of molecular switches, proteins that transition between two states-on and off Typically, signaling is initiated when an external stimulus activates its cognate receptor that, in turn, causes downstream switches to transition from off to on using one of the following mechanisms: activation, in which the transition rate from the off state to the on state increases; derepression, in which the transition rate from the on state to the off state decreases; and concerted, in which activation and derepression operate simultaneously. We use mathematical modeling to compare these signaling mechanisms in terms of their dose-response curves, response times, and abilities to process upstream fluctuations. Our analysis elucidates several operating principles for molecular switches. First, activation increases the sensitivity of the pathway, whereas derepression decreases sensitivity. Second, activation generates response times that decrease with signal strength, whereas derepression causes response times to increase with signal strength. These opposing features allow the concerted mechanism to not only show dose-response alignment, but also to decouple the response time from stimulus strength. However, these potentially beneficial properties come at the expense of increased susceptibility to upstream fluctuations. We demonstrate that these operating principles also hold when the models are extended to include additional features, such as receptor removal, kinetic proofreading, and cascades of switches. In total, we show how the architecture of molecular switches govern their response properties. We also discuss the biological implications of our findings.

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.

Experimental and engineering approaches to intracellular communication

Communication between and within cells is essential for multicellular life. While intracellular signal transduction pathways are often specified in molecular terms, the information content they transmit remains poorly defined. Here, we review research efforts to merge biological experimentation with concepts of communication that emerge from the engineering disciplines of signal processing and control theory. We discuss the challenges of performing experiments that quantitate information transfer at the molecular level, and we highlight recent studies that have advanced toward a clearer definition of the information content carried by signaling molecules. Across these studies, we emphasize a theme of increasingly well-matched experimental and theoretical approaches to decode the data streams directing cellular behavior.

Transferring information without distortion

Despite employing diverse molecular mechanisms, many different cell signaling systems avoid losing information by transmitting it in a linear manner.

Signaling pathways as linear transmitters

One challenge in biology is to make sense of the complexity of biological networks. A good system to approach this is signaling pathways, whose well-characterized molecular details allow us to relate the internal processes of each pathway to their input-output behavior. In this study, we analyzed mathematical models of three metazoan signaling pathways: the canonical Wnt, MAPK/ERK, and Tgfβ pathways. We find an unexpected convergence: the three pathways behave in some physiological contexts as linear signal transmitters. Testing the results experimentally, we present direct measurements of linear input-output behavior in the Wnt and ERK pathways. Analytics from each model further reveal that linearity arises through different means in each pathway, which we tested experimentally in the Wnt and ERK pathways. Linearity is a desired property in engineering where it facilitates fidelity and superposition in signal transmission. Our findings illustrate how cells tune different complex networks to converge on the same behavior.

Single-cell profiling screen identifies microtubule-dependent reduction of variability in signaling

Populations of isogenic cells often respond coherently to signals, despite differences in protein abundance and cell state. Previously, we uncovered processes in the Saccharomyces cerevisiae pheromone response system (PRS) that reduced cell-to-cell variability in signal strength and cellular response. Here, we screened 1,141 non-essential genes to identify 50 "variability genes". Most had distinct, separable effects on strength and variability of the PRS, defining these quantities as genetically distinct "axes" of system behavior. Three genes affected cytoplasmic microtubule function: BIM1, GIM2, and GIM4 We used genetic and chemical perturbations to show that, without microtubules, PRS output is reduced but variability is unaffected, while, when microtubules are present but their function is perturbed, output is sometimes lowered, but its variability is always high. The increased variability caused by microtubule perturbations required the PRS MAP kinase Fus3 and a process at or upstream of Ste5, the membrane-localized scaffold to which Fus3 must bind to be activated. Visualization of Ste5 localization dynamics demonstrated that perturbing microtubules destabilized Ste5 at the membrane signaling site. The fact that such microtubule perturbations cause aberrant fate and polarity decisions in mammals suggests that microtubule-dependent signal stabilization might also operate throughout metazoans.

Analysis of the thresholds for transcriptional activation by the yeast MAP kinases Fus3 and Kss1

Signaling in the pheromone response pathway of budding yeast activates two distinct MAP kinases (MAPKs), Fus3 and Kss1. Either MAPK alone can mediate pheromone-induced transcription, but it has been unclear to what degree each one contributes to transcriptional output in wild-type cells. Here, we report that transcription reflects the ratio of active to inactive MAPK, and not simply the level of active MAPK. For Kss1 the majority of MAPK molecules must be converted to the active form, whereas for Fus3 only a small minority must be activated. These different activation thresholds reflect two opposing effects of each MAPK, in which the inactive forms inhibit transcription, whereas the active forms promote transcription. Moreover, negative feedback from Fus3 limits activation of Kss1 so that it does not meet its required threshold in wild-type cells but does so only when hyperactivated in cells lacking Fus3. The results suggest that the normal transcriptional response involves asymmetric contributions from the two MAPKs, in which pheromone signaling reduces the negative effect of Kss1 while increasing the positive effect of Fus3. These findings reveal new functional distinctions between these MAPKs, and help illuminate how inhibitory functions shape positive pathway outputs in both pheromone and filamentation pathways.

Yeast GPCR signaling reflects the fraction of occupied receptors, not the number

According to receptor theory, the effect of a ligand depends on the amount of agonist-receptor complex. Therefore, changes in receptor abundance should have quantitative effects. However, the response to pheromone in Saccharomyces cerevisiae is robust (unaltered) to increases or reductions in the abundance of the G-protein-coupled receptor (GPCR), Ste2, responding instead to the fraction of occupied receptor. We found experimentally that this robustness originates during G-protein activation. We developed a complete mathematical model of this step, which suggested the ability to compute fractional occupancy depends on the physical interaction between the inhibitory regulator of G-protein signaling (RGS), Sst2, and the receptor. Accordingly, replacing Sst2 by the heterologous hsRGS4, incapable of interacting with the receptor, abolished robustness. Conversely, forcing hsRGS4:Ste2 interaction restored robustness. Taken together with other results of our work, we conclude that this GPCR pathway computes fractional occupancy because ligand-bound GPCR-RGS complexes stimulate signaling while unoccupied complexes actively inhibit it. In eukaryotes, many RGSs bind to specific GPCRs, suggesting these complexes with opposing activities also detect fraction occupancy by a ratiometric measurement. Such complexes operate as push-pull devices, which we have recently described.

Signaling Tug-of-War Delivers the Whole Message

How do cells transmit biochemical signals accurately? It turns out, pushing and pulling can go a long way.