Robust perfect adaptation in bacterial chemotaxis through integral feedback control

PD-Based ADRC using time-varying gains: An Application to Microalgal-based bioprocess

Microalgae cultivation has gained significant attention in recent years due to its potential applications in various fields. However, achieving high productivity in these bioprocesses requires efficient process control. The dynamics of growth models typically comprise both known and unknown components due to mismatches between the nonlinear dynamics and their mathematical representations. Additionally, microalgal culture is subject to external disturbances. To address these issues, a classical proportional-derivative (PD) providing the feedback error, assisted by a time-varying gain extended state observer (ESO), maintaining the structure of an Active Disturbance Rejection Control (ADRC), was implemented. The formulation is aided by a time-varying gain extended state observer to avoid high-peaking estimation values. The optimal operating conditions were identified using the GEKKO Python package. The proposed controller was applied to the growth model of the microalga Isochrysis galbana, and numerical results demonstrated the effectiveness of the control strategy in eliminating steady-state error and ensuring asymptotic convergence to the optimal equilibrium despite unknown disturbances. A detailed analysis of the photobioreactor model, including stability under steady-state conditions, was also conducted. The results indicated that the model exhibits one, two, or no stable steady-state solutions when the dilution rate (D(t)) is manipulated.

Lattice ultrasensitivity amplifies signals in E. coli without fine-tuning

The E. coli chemosensory lattice, consisting of receptors, kinases, and adaptor proteins, is an important test case for biochemical signal processing. Kinase output is characterized by precise adaptation to a wide range of background ligand levels and large gain in response to small relative changes in concentration. Existing models of this lattice achieve their gain through allosteric interactions between either receptors or core units of receptors and kinases. Here we introduce a model which operates through an entirely different mechanism in which receptors gate inherently far from equilibrium enzymatic reactions between neighboring kinases. Our lattice model achieves gain through a mechanism more closely related to zero-order ultrasensitivity than to allostery. Thus, we call it lattice ultrasensitivity (LU). Unlike other lattice critical models, the LU model can achieve arbitrarily high gain through time-scale separation, rather than through finetuning. The model also captures qualitative experimental results which are difficult to reconcile with existing models. We discuss possible implementations in the lattice's baseplate where long flexible linkers could potentially mediate interactions between neighboring core units.

A functional approach to homeostatic regulation

In this work, we present a novel modeling framework for understanding the dynamics of homeostatic regulation. Inspired by engineering control theory, this framework incorporates unique features of biological systems. First, biological variables often play physiological roles, and taking this functional context into consideration is essential to fully understand the goals and constraints of homeostatic regulation. Second, biological signals are not abstract variables, but rather material molecules that may undergo complex turnover processes of synthesis and degradation. We suggest that the particular nature of biological signals may condition the type of information they can convey, and their potential role in shaping the dynamics and the ultimate purpose of homeostatic systems. We show that the dynamic interplay between regulated variables and control signals is a key determinant of biological homeostasis, challenging the necessity and the convenience of strictly extrapolating concepts from engineering control theory in modeling the dynamics of homeostatic systems. This work provides a simple, unified framework for studying biological regulation and identifies general principles that transcend molecular details of particular homeostatic mechanisms. We show how this approach can be naturally applied to apparently different regulatory systems, contributing to a deeper understanding of homeostasis as a fundamental process in living systems.

Brain and cognition: The need for a broader biological perspective to overcome old biases

Even with accumulating knowledge, no consensus regarding the understanding of intelligence or cognition exists, and the universal brain bases for these functions remain unclear. Traditionally, our understanding of cognition is based on self-evident principles that appear indisputable when looking only at our species; however, this can distance us from understanding its essence (anthropocentrism, corticocentrism, intellectocentrism, neurocentrism, and idea of orthogenesis of brain evolution). Herein, we use several examples from biology to demonstrate the usefulness of comparative ways of thinking in relativizing these biases. We discuss the relationship between the number of neurons and cognition and draw attention to the highly developed cognitive performance of animals with small brains, to some "tricks" of evolution, to how animals cope with small hardware, to some animals with high-quality brains with an alternative architecture to vertebrates, and to surprising basal cognitive skills in aneural, unicellular organisms. Cognition can be supplemented by the idea of a multicellular organism as a continuum, with many levels of cognitive function, including the possible basal learning in single cells.

Receptor polarization through localized activity and global sensitization

Eukaryotic cells chemosense concentration gradients of extracellular ligands using membrane-bound receptors that polarize their activity. Receptors from several chemosensing families are preferentially degraded after activation and undergo significant lateral diffusion, both of which may blunt their polarization. To explore the combined role of these two seemingly detrimental phenomena on active receptor polarization, we use a reaction/diffusion model. The model elucidates a counterintuitive principle that governs receptor polarization under external gradients: Localized Activity and Global Sensitization (LAGS). In LAGS, receptor activity is localized through receptor degradation or ligand unbinding. In contrast, uniform sensitivity to ligands is maintained over the plasma membrane through lateral receptor diffusion. Surprisingly, increasing preferential degradation of active receptors and increasing lateral diffusion of all receptors both sharpen active receptor polarization. Additionally, when combined with receptor oligomerization, an increase in preferential degradation allows cells to sense relative ligand gradients over a larger range of background ligand concentrations. An analytical model identifies parameter regimes that dictate which processes dominate receptor polarization. A survey of kinetic parameters suggests that receptor polarization in many mammalian pathways can be modeled using LAGS.

Bifurcation enhances temporal information encoding in the olfactory periphery

Living systems continually respond to signals from the surrounding environment. Survival requires that their responses adapt quickly and robustly to the changes in the environment. One particularly challenging example is olfactory navigation in turbulent plumes, where animals experience highly intermittent odor signals while odor concentration varies over many length- and timescales. Here, we show theoretically that Drosophila olfactory receptor neurons (ORNs) can exploit proximity to a bifurcation point of their firing dynamics to reliably extract information about the timing and intensity of fluctuations in the odor signal, which have been shown to be critical for odor-guided navigation. Close to the bifurcation, the system is intrinsically invariant to signal variance, and information about the timing, duration, and intensity of odor fluctuations is transferred efficiently. Importantly, we find that proximity to the bifurcation is maintained by mean adaptation alone and therefore does not require any additional feedback mechanism or fine-tuning. Using a biophysical model with calcium-based feedback, we demonstrate that this mechanism can explain the measured adaptation characteristics of Drosophila ORNs.

Bifurcation enhances temporal information encoding in the olfactory periphery

Living systems continually respond to signals from the surrounding environment. Survival requires that their responses adapt quickly and robustly to the changes in the environment. One particularly challenging example is olfactory navigation in turbulent plumes, where animals experience highly intermittent odor signals while odor concentration varies over many length- and timescales. Here, we show theoretically that Drosophila olfactory receptor neurons (ORNs) can exploit proximity to a bifurcation point of their firing dynamics to reliably extract information about the timing and intensity of fluctuations in the odor signal, which have been shown to be critical for odor-guided navigation. Close to the bifurcation, the system is intrinsically invariant to signal variance, and information about the timing, duration, and intensity of odor fluctuations is transferred efficiently. Importantly, we find that proximity to the bifurcation is maintained by mean adaptation alone and therefore does not require any additional feedback mechanism or fine-tuning. Using a biophysical model with calcium-based feedback, we demonstrate that this mechanism can explain the measured adaptation characteristics of Drosophila ORNs.

CheB localizes to polar receptor arrays during repellent adaptation

Adaptation of the response to stimuli is a fundamental process for all organisms. Here, we show that the adaptation enzyme CheB methylesterase of Escherichia coli assembles to the ON state receptor array after exposure to the repellent l-isoleucine and dissociates from the array after adaptation is complete. The duration of increased CheB localization and the time of highly clockwise-biased flagellar rotation were similar and depended on the strength of the stimulus. The increase in CheB at the receptor array and the decrease in cytoplasmic CheB were both ~100 molecules, which represents 15 to 20% of the total cellular content of CheB. We confirmed that the main binding site for CheB in the ON state array is the P2 domain of phosphorylated CheA, with a second minor site being the carboxyl-terminal pentapeptide of the serine chemoreceptor. Thus, we have been able to quantify the regulation of the signal output of the receptor array by the intracellular dynamics of an adaptation enzyme.

Background compensation revisited: Conserved phase response curves in frequency controlled homeostats with coherent feedback

Background compensation is the ability of a controlled variable to respond to an applied perturbation in an unchanged manner and independent of different but constant background signals which act in parallel to the perturbation. Background compensation occurs by 'coherent feedback' mechanisms where additional control variables feed directly back to the controlled variable. This paper extends a previous study on background compensation to include phase responses in frequency controlled coherent feedback oscillators. While the frequency resetting amplitude in coherent feedback oscillators is found to be dependent on the inflow/outflow perturbation of the controlled variable and thereby become phase dependent, the frequency resetting itself and the corresponding phase response curves are found to be background compensated. It is speculated that this type of background compensation may be an additional way how ambient noise can be 'ignored' by organisms.

Identifying and explaining resilience in ecological networks

Resilient ecological systems are more likely to persist and function in the Anthropocene. Current methods for estimating an ecosystem's resilience rely on accurately parameterized ecosystem models, which is a significant empirical challenge. In this paper, we adapt tools from biochemical kinetics to identify ecological networks that exhibit 'structural resilience', a strong form of resilience that is solely a property of the network structure and is independent of model parameters. We undertake an exhaustive search for structural resilience across all three-species ecological networks, under a generalized Lotka-Volterra modelling framework. Out of 20,000 possible network structures, approximately 2% display structural resilience. The properties of these networks provide important insights into the mechanisms that could promote resilience in ecosystems, provide new theoretical avenues for qualitative modelling approaches and provide a foundation for identifying robust forms of ecological resilience in large, realistic ecological networks.

Revisiting kinetic Monte Carlo algorithms for time-dependent processes: From open-loop control to feedback control

Simulating stochastic systems with feedback control is challenging due to the complex interplay between the system's dynamics and the feedback-dependent control protocols. We present a single-step-trajectory probability analysis to time-dependent stochastic systems. Based on this analysis, we revisit several time-dependent kinetic Monte Carlo (KMC) algorithms designed for systems under open-loop-control protocols. Our analysis provides a unified alternative proof to these algorithms, summarized into a pedagogical tutorial. Moreover, with the trajectory probability analysis, we present a novel feedback-controlled KMC algorithm that accurately captures the dynamics systems controlled by an external signal based on the measurements of the system's state. Our method correctly captures the system dynamics and avoids the artificial Zeno effect that arises from incorrectly applying the direct Gillespie algorithm to feedback-controlled systems. This work provides a unified perspective on existing open-loop-control KMC algorithms and also offers a powerful and accurate tool for simulating stochastic systems with feedback control.

Signal integration and adaptive sensory diversity tuning in Escherichia coli chemotaxis

In uncertain environments, phenotypic diversity can be advantageous for survival. However, as the environmental uncertainty decreases, the relative advantage of having diverse phenotypes decreases. Here, we show how populations of E. coli integrate multiple chemical signals to adjust sensory diversity in response to changes in the prevalence of each ligand in the environment. Measuring kinase activity in single cells, we quantified the sensitivity distribution to various chemoattractants in different mixtures of background stimuli. We found that when ligands bind uncompetitively, the population tunes sensory diversity to each signal independently, decreasing diversity when the signal's ambient concentration increases. However, among competitive ligands, the population can only decrease sensory diversity one ligand at a time. Mathematical modeling suggests that sensory diversity tuning benefits E. coli populations by modulating how many cells are committed to tracking each signal proportionally as their prevalence changes.

Design Principles for Perfect Adaptation in Biological Networks with Nonlinear Dynamics

Establishing a mapping between the emergent biological properties and the repository of network structures has been of great relevance in systems and synthetic biology. Adaptation is one such biological property of paramount importance that promotes regulation in the presence of environmental disturbances. This paper presents a nonlinear systems theory-driven framework to identify the design principles for perfect adaptation with respect to external disturbances of arbitrary magnitude. Based on the prior information about the network, we frame precise mathematical conditions for adaptation using nonlinear systems theory. We first deduce the mathematical conditions for perfect adaptation for constant input disturbances. Subsequently, we translate these conditions to specific necessary structural requirements for adaptation in networks of small size and then extend to argue that there exist only two classes of architectures for a network of any size that can provide local adaptation in the entire state space, namely, incoherent feed-forward (IFF) structure and negative feedback loop with buffer node (NFB). The additional positiveness constraints further narrow the admissible set of network structures. This also aids in establishing the global asymptotic stability for the steady state given a constant input disturbance. The proposed method does not assume any explicit knowledge of the underlying rate kinetics, barring some minimal assumptions. Finally, we also discuss the infeasibility of certain IFF networks in providing adaptation in the presence of downstream connections. Moreover, we propose a generic and novel algorithm based on non-linear systems theory to unravel the design principles for global adaptation. Detailed and extensive simulation studies corroborate the theoretical findings.

Signal Integration and Adaptive Sensory Diversity Tuning in Escherichia coli Chemotaxis

In uncertain environments, phenotypic diversity can be advantageous for survival. However, as the environmental uncertainty decreases, the relative advantage of having diverse phenotypes decreases. Here, we show how populations of E. coli integrate multiple chemical signals to adjust sensory diversity in response to changes in the prevalence of each ligand in the environment. Measuring kinase activity in single cells, we quantified the sensitivity distribution to various chemoattractants in different mixtures of background stimuli. We found that when ligands bind uncompetitively, the population tunes sensory diversity to each signal independently, decreasing diversity when the signal ambient concentration increases. However, amongst competitive ligands the population can only decrease sensory diversity one ligand at a time. Mathematical modeling suggests that sensory diversity tuning benefits E. coli populations by modulating how many cells are committed to tracking each signal proportionally as their prevalence changes.

Strong delayed negative feedback

In this paper, we analyze the strong feedback limit of two negative feedback schemes which have proven to be efficient for many biological processes (protein synthesis, immune responses, breathing disorders). In this limit, the nonlinear delayed feedback function can be reduced to a function with a threshold nonlinearity. This will considerably help analytical and numerical studies of networks exhibiting different topologies. Mathematically, we compare the bifurcation diagrams for both the delayed and non-delayed feedback functions and show that Hopf classical theory needs to be revisited in the strong feedback limit.

Design principles as minimal models

In this essay I suggest that we view design principles in systems biology as minimal models, for a design principle usually exhibits universal behaviors that are common to a whole range of heterogeneous (living and nonliving) systems with different underlying mechanisms. A well-known design principle in systems biology, integral feedback control, is discussed, showing that it satisfies all the conditions for a model to be a minimal model. This approach has significant philosophical implications: it not only accounts for how design principles explain, but also helps clarify one dispute over design principles, e.g., whether design principles provide mechanistic explanations or a distinct kind of explanations called design explanations.

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.

Spatial control of sensory adaptation modulates mechanosensing in Pseudomonas aeruginosa

Sensory signaling pathways use adaptation to dynamically respond to changes in their environment. Here, we report the mechanism of sensory adaptation in the Pil-Chp mechanosensory system, which the important human pathogen Pseudomonas aeruginosa uses to sense mechanical stimuli during surface exploration. Using biochemistry, genetics, and cell biology, we discovered that the enzymes responsible for adaptation, a methyltransferase and a methylesterase, are segregated to opposing cell poles as P. aeruginosa explore surfaces. By coordinating the localization of both enzymes, we found that the Pil-Chp response regulators influence local receptor methylation, the molecular basis of bacterial sensory adaptation. We propose a model in which adaptation during mechanosensing spatially resets local receptor methylation, and thus Pil-Chp signaling, to modulate the pathway outputs, which are involved in P. aeruginosa virulence. Despite decades of bacterial sensory adaptation studies, our work has uncovered an unrecognized mechanism that bacteria use to achieve adaptation to sensory stimuli.

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.

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

Relation between learning process and morphology of transport tube network in plasmodium of Physarum polycephalum

The question of whether a single-celled organism without a brain could have functions such as learning and memory has been the subject of much debate in recent years. The plasmodium of the true slime mold, Physarum polycephalum, is an ideal model organism for such a question. The plasmodium exhibits behaviors that resemble intelligence, including solving mazes, mimicking optimal rail transportation networks, predicting the weather, and solving traveling salesman problems. In addition, the plasmodium has recently been shown to have the simplest form of learning: habituation. In the experiments in which plasmodia were repeatedly allowed to cross bridges containing aversive chemicals, the habituation behavior has been confirmed. It has been shown that the habituation process involves chemicals that are stored internally. However, it is not clear how these chemicals result in change in the behavior of plasmodium during habituation learning. This study focused on the transport tube network formed in plasmodium during the above experiments. Then, the role of the network morphology in the habituation learning process was investigated. The results showed that the network morphology changes from tree to mesh type during habituation learning, and disrupting the learned network reduces habituation behavior. In addition, it was shown that the thickness oscillation frequency depends on the network morphology. The study found that in the plasmodium of P. polycephalum, a primitive organism without a brain, transport tube networks, instead of neuronal networks, play an important role in habituation learning and the resulting decision making.

Cell behaviors underlying Myxococcus xanthus aggregate dispersal

IMPORTANCE

Understanding the processes behind bacterial biofilm formation, maintenance, and dispersal is essential for addressing their effects on health and ecology. Within these multicellular communities, various cues can trigger differentiation into distinct cell types, allowing cells to adapt to their specific local environment. The soil bacterium Myxococcus xanthus forms biofilms in response to starvation, marked by cells aggregating into mounds. Some aggregates persist as spore-filled fruiting bodies, while others disperse after initial formation for unknown reasons. Here, we use a combination of cell tracking analysis and computational simulations to identify behaviors at the cellular level that contribute to aggregate dispersal. Our results suggest that cells in aggregates actively determine whether to disperse or persist and undergo a transition to sporulation based on a self-produced cue related to the aggregate size. Identifying these cues is an important step in understanding and potentially manipulating bacterial cell-fate decisions.

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.

Local negative feedback of Rac activity at the leading edge underlies a pilot pseudopod-like program for amoeboid cell guidance

To migrate efficiently, neutrophils must polarize their cytoskeletal regulators along a single axis of motion. This polarization process is thought to be mediated through local positive feedback that amplifies leading edge signals and global negative feedback that enables sites of positive feedback to compete for dominance. Though this two-component model efficiently establishes cell polarity, it has potential limitations, including a tendency to "lock" onto a particular direction, limiting the ability of cells to reorient. We use spatially defined optogenetic control of a leading edge organizer (PI3K) to probe how neutrophil-like HL-60 cells balance "decisiveness" needed to polarize in a single direction with the flexibility needed to respond to new cues. Underlying this balancing act is a local Rac inhibition process that destabilizes the leading edge to promote exploration. We show that this local inhibition enables cells to process input signal dynamics, linking front stability and orientation to local temporal increases in input signals.

Coherent feedback leads to robust background compensation in oscillatory and non-oscillatory homeostats

When in a reaction kinetic integral controller a step perturbation is applied besides a constant background, the concentration of a controlled variable (described as A) will generally respond with decreased response amplitudes ΔA as backgrounds increase. The controller variable E will at the same time provide the necessary compensatory flux to move A back to its set-point. A typical example of decreased response amplitudes at increased backgrounds is found in retinal light adaptation. Due to remarks in the literature that retinal light adaptation would also involve a compensation of backgrounds we became interested in conditions how background compensation could occur. In this paper we describe novel findings how background influences can be robustly eliminated. When such a background compensation is active, oscillatory controllers will respond to a defined perturbation with always the same (damped or undamped) frequency profile, or in the non-oscillatory case, with the same response amplitude ΔA, irrespective of the background level. To achieve background compensation we found that two conditions need to apply: (i) an additional set of integral controllers (here described as I1 and I2) have to be employed to keep the manipulated variable E at a defined set-point, and (ii), I1 and I2 need to feed back to the A-E signaling axis directly through the controlled variable A. In analogy to a similar feedback applied in quantum control theory, we term these feedback conditions as 'coherent feedback'. When analyzing retinal light adaptations in more detail, we find no evidence of the presence of background compensation mechanisms. Although robust background compensation, as described theoretically here, appears to be an interesting regulatory property, relevant biological or biochemical examples still need to be identified.

A General Approach for the Modelling of Negative Feedback Physiological Control Systems

Mathematical models can improve the understanding of physiological systems behaviour, which is a fundamental topic in the bioengineering field. Having a reliable model enables researchers to carry out in silico experiments, which require less time and resources compared to their in vivo and in vitro counterparts. This work's objective is to capture the characteristics that a nonlinear dynamical mathematical model should exhibit, in order to describe physiological control systems at different scales. The similarities among various negative feedback physiological systems have been investigated and a unique general framework to describe them has been proposed. Within such a framework, both the existence and stability of equilibrium points are investigated. The model here introduced is based on a closed-loop topology, on which the homeostatic process is based. Finally, to validate the model, three paradigmatic examples of physiological control systems are illustrated and discussed: the ultrasensitivity mechanism for achieving homeostasis in biomolecular circuits, the blood glucose regulation, and the neuromuscular reflex arc (also referred to as muscle stretch reflex). The results show that, by a suitable choice of the modelling functions, the dynamic evolution of the systems under study can be described through the proposed general nonlinear model. Furthermore, the analysis of the equilibrium points and dynamics of the above-mentioned systems are consistent with the literature.

Perspectives on adaptive dynamical systems

Adaptivity is a dynamical feature that is omnipresent in nature, socio-economics, and technology. For example, adaptive couplings appear in various real-world systems, such as the power grid, social, and neural networks, and they form the backbone of closed-loop control strategies and machine learning algorithms. In this article, we provide an interdisciplinary perspective on adaptive systems. We reflect on the notion and terminology of adaptivity in different disciplines and discuss which role adaptivity plays for various fields. We highlight common open challenges and give perspectives on future research directions, looking to inspire interdisciplinary approaches.

Evolutionary Stability of Small Molecular Regulatory Networks That Exhibit Near-Perfect Adaptation

Large-scale protein regulatory networks, such as signal transduction systems, contain small-scale modules ('motifs') that carry out specific dynamical functions. Systematic characterization of the properties of small network motifs is therefore of great interest to molecular systems biologists. We simulate a generic model of three-node motifs in search of near-perfect adaptation, the property that a system responds transiently to a change in an environmental signal and then returns near-perfectly to its pre-signal state (even in the continued presence of the signal). Using an evolutionary algorithm, we search the parameter space of these generic motifs for network topologies that score well on a pre-defined measure of near-perfect adaptation. We find many high-scoring parameter sets across a variety of three-node topologies. Of all possibilities, the highest scoring topologies contain incoherent feed-forward loops (IFFLs), and these topologies are evolutionarily stable in the sense that, under 'macro-mutations' that alter the topology of a network, the IFFL motif is consistently maintained. Topologies that rely on negative feedback loops with buffering (NFLBs) are also high-scoring; however, they are not evolutionarily stable in the sense that, under macro-mutations, they tend to evolve an IFFL motif and may-or may not-lose the NFLB motif.

Universal structures for adaptation in biochemical reaction networks

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

Wild animals suppress the spread of socially transmitted misinformation

Understanding the mechanisms by which information and misinformation spread through groups of individual actors is essential to the prediction of phenomena ranging from coordinated group behaviors to misinformation epidemics. Transmission of information through groups depends on the rules that individuals use to transform the perceived actions of others into their own behaviors. Because it is often not possible to directly infer decision-making strategies in situ, most studies of behavioral spread assume that individuals make decisions by pooling or averaging the actions or behavioral states of neighbors. However, whether individuals may instead adopt more sophisticated strategies that exploit socially transmitted information, while remaining robust to misinformation, is unknown. Here, we study the relationship between individual decision-making and misinformation spread in groups of wild coral reef fish, where misinformation occurs in the form of false alarms that can spread contagiously through groups. Using automated visual field reconstruction of wild animals, we infer the precise sequences of socially transmitted visual stimuli perceived by individuals during decision-making. Our analysis reveals a feature of decision-making essential for controlling misinformation spread: dynamic adjustments in sensitivity to socially transmitted cues. This form of dynamic gain control can be achieved by a simple and biologically widespread decision-making circuit, and it renders individual behavior robust to natural fluctuations in misinformation exposure.

Cellular learning: Habituation sans neurons in a unicellular organism

Stentor coeruleus cells stochastically switch between non-responsive (contracted) and responsive (extended) states. Learning is accomplished via habituation, in which the internal model is updated to reflect the current environment by tuning the transition rates according to the time series properties of mechanical stimuli.

Epigenetic inheritance of gene silencing is maintained by a self-tuning mechanism based on resource competition

Biological systems can maintain memories over long timescales, with examples including memories in the brain and immune system. It is unknown how functional properties of memory systems, such as memory persistence, can be established by biological circuits. To address this question, we focus on transgenerational epigenetic inheritance in Caenorhabditis elegans. In response to a trigger, worms silence a target gene for multiple generations, resisting strong dilution due to growth and reproduction. Silencing may also be maintained indefinitely upon selection according to silencing levels. We show that these properties imply the fine-tuning of biochemical rates in which the silencing system is positioned near the transition to bistability. We demonstrate that this behavior is consistent with a generic mechanism based on competition for synthesis resources, which leads to self-organization around a critical state with broad silencing timescales. The theory makes distinct predictions and offers insights into the design principles of long-term memory systems.

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.

Design Principles Underlying Robust Adaptation of Complex Biochemical Networks

Biochemical networks are often characterized by tremendous complexity-both in terms of the sheer number of interacting molecules ("nodes") and in terms of the varied and incompletely understood interactions among these molecules ("interconnections" or "edges"). Strikingly, the vast and intricate networks of interacting proteins that exist within each living cell have the capacity to perform remarkably robustly, and reproducibly, despite significant variations in concentrations of the interacting components from one cell to the next and despite mutability over time of biochemical parameters. Here we consider the ubiquitously observed and fundamentally important signalling response known as robust perfect adaptation (RPA). We have recently shown that all RPA-capable networks, even the most complex ones, must satisfy an extremely rigid set of design principles, and are modular, being decomposable into just two types of network building-blocks-opposer modules and balancer modules. Here we present an overview of the design principles that characterize all RPA-capable network topologies through a detailed examination of a collection of simple examples. We also introduce a diagrammatic method for studying the potential of a network to exhibit RPA, which may be applied without a detailed knowledge of the complex mathematical principles governing RPA.

The Chemoreceptor Sensory Adaptation System Produces Coordinated Reversals of the Flagellar Motors on an Escherichia coli Cell

In isotropic environments, an Escherichia coli cell exhibits coordinated rotational switching of its flagellar motors, produced by fluctuations in the intracellular concentration of phosphorylated CheY (CheY-P) emanating from chemoreceptor signaling arrays. In this study, we show that these CheY-P fluctuations arise through modifications of chemoreceptors by two sensory adaptation enzymes: the methyltransferase CheR and the methylesterase CheB. A cell containing CheR, CheB, and the serine chemoreceptor Tsr exhibited motor synchrony, whereas a cell lacking CheR and CheB or containing enzymatically inactive forms did not. Tsr variants with different combinations of methylation-mimicking Q residues at the adaptation sites also failed to show coordinated motor switching in cells lacking CheR and CheB. Cells containing CheR, CheB, and Tsr [NDND], a variant in which the adaptation site residues are not substrates for CheR or CheB modifications, also lacked motor synchrony. TsrΔNWETF, which lacks a C-terminal pentapeptide-binding site for CheR and CheB, and the ribose-galactose receptor Trg, which natively lacks this motif, failed to produce coordinated motor switching, despite the presence of CheR and CheB. However, addition of the NWETF sequence to Trg enabled Trg-NWETF to produce motor synchrony, as the sole receptor type in cells containing CheR and CheB. Finally, CheBc, the catalytic domain of CheB, supported motor coordination in combination with CheR and Tsr. These results indicate that the coordination of motor switching requires CheR/CheB-mediated changes in receptor modification state. We conclude that the opposing receptor substrate-site preferences of CheR and CheB produce spontaneous blinking of the chemoreceptor array's output activity. IMPORTANCE Under steady-state conditions with no external stimuli, an Escherichia coli cell coordinately switches the rotational direction of its flagellar motors. Here, we demonstrate that the CheR and CheB enzymes of the chemoreceptor sensory adaptation system mediate this coordination. Stochastic fluctuations in receptor adaptation states trigger changes in signal output from the receptor array, and this array blinking generates fluctuations in CheY-P concentration that coordinate directional switching of the flagellar motors. Thus, in the absence of chemoeffector gradients, the sensory adaptation system controls run-tumble swimming of the cell, its optimal foraging strategy.

Feedback-controlled hydrogels with homeostatic oscillations and dissipative signal transduction

Driving systems out of equilibrium under feedback control is characteristic for living systems, where homeostasis and dissipative signal transduction facilitate complex responses. This feature not only inspires dissipative dynamic functionalities in synthetic systems but also poses great challenges in designing novel pathways. Here we report feedback-controlled systems comprising two coupled hydrogels driven by constant light, where the system can be tuned to undergo stable homeostatic self-oscillations or damped steady states of temperature. We demonstrate that stable temperature oscillations can be utilized for dynamic colours and cargo transport, whereas damped steady states enable signal transduction pathways. Here mechanical triggers cause temperature changes that lead to responses such as bending motions inspired by the single-touch mechanoresponse in Mimosa pudica and the frequency-gated snapping motion inspired by the plant arithmetic in the Venus flytrap. The proposed concepts suggest generalizable feedback pathways for dissipative dynamic materials and interactive soft robotics.

Design principles and mechanistic explanation

In this essay I propose that what design principles in systems biology and systems neuroscience do is to present abstract characterizations of mechanisms, and thereby facilitate mechanistic explanation. To show this, one design principle in systems neuroscience, i.e., the multilayer perceptron, is examined. However, Braillard (2010) contends that design principles provide a sort of non-mechanistic explanation due to two related reasons: they are very general and describe non-causal dependence relationships. In response to this, I argue that, on the one hand, all mechanisms are more or less general (or abstract), and on the other, many (if not all) design principles are causal systems.

Universal structural requirements for maximal robust perfect adaptation in biomolecular networks

Adaptation is a running theme in biology. It allows a living system to survive and thrive in the face of unpredictable environments by maintaining key physiological variables at their desired levels through tight regulation. When one such variable is maintained at a certain value at the steady state despite perturbations to a single input, this property is called robust perfect adaptation (RPA). Here we address and solve the fundamental problem of maximal RPA (maxRPA), whereby, for a designated output variable, RPA is achieved with respect to perturbations in virtually all network parameters. In particular, we show that the maxRPA property imposes certain structural constraints on the network. We then prove that these constraints are fully characterized by simple linear algebraic stoichiometric conditions which differ between deterministic and stochastic descriptions of the dynamics. We use our results to derive a new internal model principle (IMP) for biomolecular maxRPA networks, akin to the celebrated IMP in control theory. We exemplify our results through several known biological examples of robustly adapting networks and construct examples of such networks with the aid of our linear algebraic characterization. Our results reveal the universal requirements for maxRPA in all biological systems, and establish a foundation for studying adaptation in general biomolecular networks, with important implications for both systems and synthetic biology.

Vertex results for the robust analysis of uncertain biochemical systems

We consider the problem of assessing the sensitivity of uncertain biochemical systems in the presence of input perturbations (either constant or periodic) around a stable steady state. In particular, we propose approaches for the robust sensitivity analysis of systems with uncertain parameters assumed to take values in a hyper-rectangle. We highlight vertex results, which allow us to check whether a property is satisfied for all parameter choices in the hyper-rectangle by simply checking whether it is satisfied for all parameter choices at the vertices of the hyper-rectangle. We show that, for a vast class of systems, including (bio)chemical reaction networks with mass-action kinetics, the system Jacobian has a totally multiaffine structure (namely, all minors of the Jacobian matrix are multiaffine functions of the uncertain parameters), which can be exploited to obtain several vertex results. We consider different problems: robust non-singularity; robust stability of the steady-state; robust steady-state sensitivity analysis, in the case of constant perturbations; robust frequency-response sensitivity analysis, in the presence of periodic perturbations; and robust adaptation analysis. The developed theory is then applied to gain insight into some examples of uncertain biochemical systems, including the incoherent feed-forward loop, the coherent feed-forward loop, the Brusselator oscillator and the Goldbeter oscillator.

CoRa-A general approach for quantifying biological feedback control

Feedback control is a fundamental underpinning of life, underlying homeostasis of biological processes at every scale of organization, from cells to ecosystems. The ability to evaluate the contribution and limitations of feedback control mechanisms operating in cells is a critical step for understanding and ultimately designing feedback control systems with biological molecules. Here, we introduce CoRa—or Control Ratio—a general framework that quantifies the contribution of a biological feedback control mechanism to adaptation using a mathematically controlled comparison to an identical system that does not contain the feedback. CoRa provides a simple and intuitive metric with broad applicability to biological feedback systems.

A negative feedback loop in the GPCR pathway underlies efficient coding of external stimuli

Efficient navigation based on chemical cues is an essential feature shared by all animals. These cues may be encountered in complex spatiotemporal patterns and with orders of magnitude varying intensities. Nevertheless, sensory neurons accurately extract the relevant information from such perplexing signals. Here, we show how a single sensory neuron in Caenorhabditis elegans animals can cell-autonomously encode complex stimulus patterns composed of instantaneous sharp changes and of slowly changing continuous gradients. This encoding relies on a simple negative feedback in the G-protein-coupled receptor (GPCR) signaling pathway in which TAX-6/Calcineurin plays a key role in mediating the feedback inhibition. This negative feedback supports several important coding features that underlie an efficient navigation strategy, including exact adaptation and adaptation to the magnitude of the gradient's first derivative. A simple mathematical model explains the fine neural dynamics of both wild-type and tax-6 mutant animals, further highlighting how the calcium-dependent activity of TAX-6/Calcineurin dictates GPCR inhibition and response dynamics. As GPCRs are ubiquitously expressed in all sensory neurons, this mechanism may be a general solution for efficient cell-autonomous coding of external stimuli.

Kinetics and mechanisms of catalyzed dual-E (antithetic) controllers

Homeostasis plays a central role in our understanding how cells and organisms are able to oppose environmental disturbances and thereby maintain an internal stability. During the last two decades there has been an increased interest in using control engineering methods, especially integral control, in the analysis and design of homeostatic networks. Several reaction kinetic mechanisms have been discovered which lead to integral control. In two of them integral control is achieved, either by the removal of a single control species E by zero-order kinetics (“single-E controllers”), or by the removal of two control species by second-order kinetics (“antithetic or dual-E control”). In this paper we show results when the control species E₁ and E₂ in antithetic control are removed enzymatically by ping-pong or ternary-complex mechanisms. Our findings show that enzyme-catalyzed dual-E controllers can work in two control modes. In one mode, one of the two control species is active, but requires zero-order kinetics in its removal. In the other mode, both controller species are active and both are removed enzymatically. Conditions for the two control modes are put forward and biochemical examples with the structure of enzyme-catalyzed dual-E controllers are discussed.

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.

A genetic mammalian proportional-integral feedback control circuit for robust and precise gene regulation

To survive in the harsh environments they inhabit, cells have evolved sophisticated regulatory mechanisms that can maintain a steady internal milieu or homeostasis. This robustness, however, does not generally translate to engineered genetic circuits, such as the ones studied by synthetic biology. Here, we introduce an implementation of a minimal and universal gene regulatory motif that produces robust perfect adaptation for mammalian cells, and we improve on it by enhancing the precision of its regulation.

The processes that keep a cell alive are constantly challenged by unpredictable changes in its environment. Cells manage to counteract these changes by employing sophisticated regulatory strategies that maintain a steady internal milieu. Recently, the antithetic integral feedback motif has been demonstrated to be a minimal and universal biological regulatory strategy that can guarantee robust perfect adaptation for noisy gene regulatory networks in Escherichia coli. Here, we present a realization of the antithetic integral feedback motif in a synthetic gene circuit in mammalian cells. We show that the motif robustly maintains the expression of a synthetic transcription factor at tunable levels even when it is perturbed by increased degradation or its interaction network structure is perturbed by a negative feedback loop with an RNA-binding protein. We further demonstrate an improved regulatory strategy by augmenting the antithetic integral motif with additional negative feedback to realize antithetic proportional–integral control. We show that this motif produces robust perfect adaptation while also reducing the variance of the regulated synthetic transcription factor. We demonstrate that the integral and proportional–integral feedback motifs can mitigate the impact of gene expression burden, and we computationally explore their use in cell therapy. We believe that the engineering of precise and robust perfect adaptation will enable substantial advances in industrial biotechnology and cell-based therapeutics.

Programming Molecular Systems To Emulate a Learning Spiking Neuron

Hebbian theory seeks to explain how the neurons in the brain adapt to stimuli to enable learning. An interesting feature of Hebbian learning is that it is an unsupervised method and, as such, does not require feedback, making it suitable in contexts where systems have to learn autonomously. This paper explores how molecular systems can be designed to show such protointelligent behaviors and proposes the first chemical reaction network (CRN) that can exhibit autonomous Hebbian learning across arbitrarily many input channels. The system emulates a spiking neuron, and we demonstrate that it can learn statistical biases of incoming inputs. The basic CRN is a minimal, thermodynamically plausible set of microreversible chemical equations that can be analyzed with respect to their energy requirements. However, to explore how such chemical systems might be engineered de novo, we also propose an extended version based on enzyme-driven compartmentalized reactions. Finally, we show how a purely DNA system, built upon the paradigm of DNA strand displacement, can realize neuronal dynamics. Our analysis provides a compelling blueprint for exploring autonomous learning in biological settings, bringing us closer to realizing real synthetic biological intelligence.

Dichotomous feedback: a signal sequestration-based feedback mechanism for biocontroller design

We introduce a new design framework for implementing negative feedback regulation in synthetic biology, which we term ‘dichotomous feedback’. Our approach is different from current methods, in that it sequesters existing fluxes in the process to be controlled, and in this way takes advantage of the process’s architecture to design the control law. This signal sequestration mechanism appears in many natural biological systems and can potentially be easier to realize than ‘molecular sequestration’ and other comparison motifs that are nowadays common in biomolecular feedback control design. The loop is closed by linking the strength of signal sequestration to the process output. Our feedback regulation mechanism is motivated by two-component signalling systems, where a second response regulator could be competing with the natural response regulator thus sequestering kinase activity. Here, dichotomous feedback is established by increasing the concentration of the second response regulator as the level of the output of the natural process increases. Extensive analysis demonstrates how this type of feedback shapes the signal response, attenuates intrinsic noise while increasing robustness and reducing crosstalk.

A hierarchy of biomolecular proportional-integral-derivative feedback controllers for robust perfect adaptation and dynamic performance

Proportional-Integral-Derivative (PID) feedback controllers are the most widely used controllers in industry. Recently, the design of molecular PID-controllers has been identified as an important goal for synthetic biology and the field of cybergenetics. In this paper, we consider the realization of PID-controllers via biomolecular reactions. We propose an array of topologies offering a compromise between simplicity and high performance. We first demonstrate that different biomolecular PI-controllers exhibit different performance-enhancing capabilities. Next, we introduce several derivative controllers based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop, which can be augmented by PI-components to yield PID-controllers. We demonstrate that PID-controllers can enhance stability and dynamic performance, and can also reduce stochastic noise. Finally, we provide an experimental demonstration using a hybrid setup where in silico PID-controllers regulate a genetic circuit in single yeast cells.

Robust and tunable signal processing in mammalian cells via engineered covalent modification cycles

Engineered signaling networks can impart cells with new functionalities useful for directing differentiation and actuating cellular therapies. For such applications, the engineered networks must be tunable, precisely regulate target gene expression, and be robust to perturbations within the complex context of mammalian cells. Here, we use bacterial two-component signaling proteins to develop synthetic phosphoregulation devices that exhibit these properties in mammalian cells. First, we engineer a synthetic covalent modification cycle based on kinase and phosphatase proteins derived from the bifunctional histidine kinase EnvZ, enabling analog tuning of gene expression via its response regulator OmpR. By regulating phosphatase expression with endogenous miRNAs, we demonstrate cell-type specific signaling responses and a new strategy for accurate cell type classification. Finally, we implement a tunable negative feedback controller via a small molecule-stabilized phosphatase, reducing output expression variance and mitigating the context-dependent effects of off-target regulation and resource competition. Our work lays the foundation for establishing tunable, precise, and robust control over cell behavior with synthetic signaling networks.