Design of a Synthetic Integral Feedback Circuit: Dynamic Analysis and DNA Implementation

Multi-Layer Autocatalytic Feedback Enables Integral Control Amidst Resource Competition and Across Scales

Integral feedback control strategies have proven effective in regulating protein expression in unpredictable cellular environments. These strategies, grounded in model-based designs and control theory, have advanced synthetic biology applications. Autocatalytic integral feedback controllers, utilizing positive autoregulation for integral action, are one class of simplest architectures to design integrators. This class of controllers offers unique features, such as robustness against dilution effects and cellular growth, as well as the potential for synthetic realizations across different biological scales, owing to their similarity to self-regenerative behaviors widely observed in nature. Despite this, their potential has not yet been fully exploited. One key reason, we discuss, is that their effectiveness is often hindered by resource competition and context-dependent couplings. This study addresses these challenges using a multilayer feedback strategy. Our designs enabled population-level integral feedback and multicellular integrators, where the control function emerges as a property of coordinated interactions distributed across different cell populations coexisting in a multicellular consortium. We provide a generalized mathematical framework for modeling resource competition in complex genetic networks, supporting the design of intracellular control circuits. The use of our proposed multilayer autocatalytic controllers is examined in two typical control tasks that pose significant relevance to synthetic biology applications: concentration regulation and ratiometric control. We define a ratiometric control task and solve it using a variant of our controller. The effectiveness of our controller motifs is demonstrated through a range of application examples, from precise regulation of gene expression and gene ratios in embedded designs to population growth and coculture composition control in multicellular designs within engineered microbial ecosystems. These findings offer a versatile approach to achieving robust adaptation and homeostasis from subcellular to multicellular scales.

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.

Analysis and design of antithetic proportional-integral-derivative biocontrol-systems with species dilution

The nonlinearity and non-separability of the antithetic PID (aPID) controller have provided greater flexibility in the design of biochemical reaction networks (BCRNs), resulting in significant impacts on biocontrol-systems. Nevertheless, the dilution of control species is disregarded in designs of aPID controllers, which would lead to the failure of inhibition mechanism in the controller and loss of robust perfect adaptation (RPA)-the biological counterpart of robust steady-state tracking. Here, the impact of dilution processes on the structure of aPID is investigated in this study. It is discovered that the proportional and low-pass filters are altered when the dilution processes is present in control species, which increases the coupling between the controller parameters. Moreover, additional integrations for the reference signal and control output generated by control species dilution further leads to the loss of RPA. Subsequently, a novel aPID controller represented by BCRNs, termed quasi-aPID, has been designed to eliminate the detrimental effects of the dilution processes. In an effort to ameliorate the interdependencies among controller parameters, a degradation inhibition mechanism is employed within this controller. Furthermore, this work establishes the limiting relationship between the controller's reaction rates in order to guarantee RPA, while abstaining from the introduction of supplementary species and biochemical reactions. By using the quasi-aPID controller in both the Escherichia coli gene expression model and the whole-body cholesterol metabolism model, its effectiveness is confirmed. Simulation results demonstrate that, the quasi-aPID exhibits a smaller absolute steady-state error in both models and guarantees the RPA property.

Noise properties of adaptation-conferring biochemical control modules

A key goal of synthetic biology is to develop functional biochemical modules with network-independent properties. Antithetic integral feedback (AIF) is a recently developed control module in which two control species perfectly annihilate each other's biological activity. The AIF module confers robust perfect adaptation to the steady-state average level of a controlled intracellular component when subjected to sustained perturbations. Recent work has suggested that such robustness comes at the unavoidable price of increased stochastic fluctuations around average levels. We present theoretical results that support and quantify this trade-off for the commonly analyzed AIF variant in the idealized limit with perfect annihilation. However, we also show that this trade-off is a singular limit of the control module: Even minute deviations from perfect adaptation allow systems to achieve effective noise suppression as long as cells can pay the corresponding energetic cost. We further show that a variant of the AIF control module can achieve significant noise suppression even in the idealized limit with perfect adaptation. This atypical configuration may thus be preferable in synthetic biology applications.

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.

Biomolecular feedback controllers: from theory to applications

Billions of years of evolution have led to the creation of sophisticated genetic regulatory mechanisms that control various biological processes in a timely and precise fashion, despite their uncertain and noisy environments. Understanding such naturally existing mechanisms and even designing novel ones will have direct implications in various fields such as biotechnology, medicine, and synthetic biology. In particular, many studies have revealed that feedback-based control mechanisms inside the living cells endow the overall system with multiple attractive features, including homeostasis, noise reduction, and high dynamic performance. The remarkable interdisciplinary nature of these studies has brought together disparate disciplines such as systems/synthetic biology and control theory in an effort to design and build more powerful and reliable biomolecular control systems. Here, we review various biomolecular feedback controllers, highlight their characteristics, and point out their promising impact.

Homeostasis at different backgrounds: The roles of overlayed feedback structures in vertebrate photoadaptation

We have studied the resetting behavior of eight basic integral controller motifs with respect to different but constant backgrounds. We found that the controllers split symmetrically into two classes: one class, based on derepression of the compensatory flux, leads to more rapid resetting kinetics as backgrounds increase. The other class, which directly activates the compensatory flux, shows a slowing down in the resetting at increased backgrounds. We found a striking analogy between the resetting kinetics of vertebrate photoreceptors and controllers based on derepression, i.e. vertebrate rod or cone cells show decreased sensitivities and accelerated response kinetics as background illuminations increase. The central molecular model of vertebrate photoadaptation consists of an overlay of three negative feedback loops with cytosolic calcium ([Formula: see text]), cyclic guanosine monophosphate (cGMP) and cyclic nucleotide-gated (CNG) channels as components. While in one of the feedback loops the extrusion of [Formula: see text] by potassium-dependent sodium-calcium exchangers (NCKX) can lead to integral control with cGMP as the controlled variable, the expected robust perfect adaptation of cGMP is lost, because of the two other feedback loops. They avoid that [Formula: see text] levels become too high and toxic. Looking at psychophysical laws, we found that in all of the above mentioned basic controllers Weber's law is followed when a "just noticeable difference" (threshold) of 1% of the controlled variable's set-point was considered. Applying comparable threshold pulses or steps to the photoadaptation model we find, in agreement with experimental results, that Weber's law is followed for relatively high backgrounds, while Stephens' power law gives a better description when backgrounds are low. Limitations of our photoadaption model, in particular with respect to potassium/sodium homeostasis, are discussed. Finally, we discuss possible implication of background perturbations in biological controllers when compensatory fluxes are based on activation.

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.

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.

Frequency spectra and the color of cellular noise

The invention of the Fourier integral in the 19th century laid the foundation for modern spectral analysis methods. This integral decomposes a temporal signal into its frequency components, providing deep insights into its generating process. While this idea has precipitated several scientific and technological advances, its impact has been fairly limited in cell biology, largely due to the difficulties in connecting the underlying noisy intracellular networks to the frequency content of observed single-cell trajectories. Here we develop a spectral theory and computational methodologies tailored specifically to the computation and analysis of frequency spectra of noisy intracellular networks. Specifically, we develop a method to compute the frequency spectrum for general nonlinear networks, and for linear networks we present a decomposition that expresses the frequency spectrum in terms of its sources. Several examples are presented to illustrate how our results provide frequency-based methods for the design and analysis of noisy intracellular networks.

Platforms for Optogenetic Stimulation and Feedback Control

Harnessing the potential of optogenetics in biology requires methodologies from different disciplines ranging from biology, to mechatronics engineering, to control engineering. Light stimulation of a synthetic optogenetic construct in a given biological species can only be achieved via a suitable light stimulation platform. Emerging optogenetic applications entail a consistent, reproducible, and regulated delivery of light adapted to the application requirement. In this review, we explore the evolution of light-induction hardware-software platforms from simple illumination set-ups to sophisticated microscopy, microtiter plate and bioreactor designs, and discuss their respective advantages and disadvantages. Here, we examine design approaches followed in performing optogenetic experiments spanning different cell types and culture volumes, with induction capabilities ranging from single cell stimulation to entire cell culture illumination. The development of automated measurement and stimulation schemes on these platforms has enabled researchers to implement various in silico feedback control strategies to achieve computer-controlled living systems—a theme we briefly discuss in the last part of this review.

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.

Synthetic robust perfect adaptation achieved by negative feedback coupling with linear weak positive feedback

Unlike their natural counterparts, synthetic genetic circuits are usually fragile in the face of environmental perturbations and genetic mutations. Several theoretical robust genetic circuits have been designed, but their performance under real-world conditions has not yet been carefully evaluated. Here, we designed and synthesized a new robust perfect adaptation circuit composed of two-node negative feedback coupling with linear positive feedback on the buffer node. As a key feature, the linear positive feedback was fine-tuned to evaluate its necessity. We found that the desired function was robustly achieved when genetic parameters were varied by systematically perturbing all interacting parts within the topology, and the necessity of the completeness of the topological structures was evaluated by destroying key circuit features. Furthermore, different environmental perturbances were imposed onto the circuit by changing growth rates, carbon metabolic strategies and even chassis cells, and the designed perfect adaptation function was still achieved under all conditions. The successful design of a robust perfect adaptation circuit indicated that the top-down design strategy is capable of predictably guiding bottom-up engineering for robust genetic circuits. This robust adaptation circuit could be integrated as a motif into more complex circuits to robustly implement more sophisticated and critical biological functions.

DOPA Homeostasis by Dopamine: A Control-Theoretic View

Dopamine (DA) is an important signal mediator in the brain as well as in the periphery. The term "dopamine homeostasis" occasionally found in the literature refers to the fact that abnormal DA levels can be associated with a variety of neuropsychiatric disorders. An analysis of the negative feedback inhibition of tyrosine hydroxylase (TH) by DA indicates, with support from the experimental data, that the TH-DA negative feedback loop has developed to exhibit 3,4-dihydroxyphenylalanine (DOPA) homeostasis by using DA as a derepression regulator. DA levels generally decline when DOPA is removed, for example, by increased oxidative stress. Robust DOPA regulation by DA further implies that maximum vesicular DA levels are established, which appear necessary for a reliable translation of neural activity into a corresponding chemical transmitter signal. An uncontrolled continuous rise (windup) in DA occurs when Levodopa treatment exceeds a critical dose. Increased oxidative stress leads to the successive breakdown of DOPA homeostasis and to a corresponding reduction in DA levels. To keep DOPA regulation robust, the vesicular DA loading requires close to zero-order kinetics combined with a sufficiently high compensatory flux provided by TH. The protection of DOPA and DA due to a channeling complex is discussed.

Rapid prototyping and design of cybergenetic single-cell controllers

The design and implementation of synthetic circuits that operate robustly in the cellular context is fundamental for the advancement of synthetic biology. However, their practical implementation presents challenges due to low predictability of synthetic circuit design and time-intensive troubleshooting. Here, we present the Cyberloop, a testing framework to accelerate the design process and implementation of biomolecular controllers. Cellular fluorescence measurements are sent in real-time to a computer simulating candidate stochastic controllers, which in turn compute the control inputs and feed them back to the controlled cells via light stimulation. Applying this framework to yeast cells engineered with optogenetic tools, we examine and characterize different biomolecular controllers, test the impact of non-ideal circuit behaviors such as dilution on their operation, and qualitatively demonstrate improvements in controller function with certain network modifications. From this analysis, we derive conditions for desirable biomolecular controller performance, thereby avoiding pitfalls during its biological implementation.

Noise suppression in stochastic genetic circuits using PID controllers

Inside individual cells, protein population counts are subject to molecular noise due to low copy numbers and the inherent probabilistic nature of biochemical processes. We investigate the effectiveness of proportional, integral and derivative (PID) based feedback controllers to suppress protein count fluctuations originating from two noise sources: bursty expression of the protein, and external disturbance in protein synthesis. Designs of biochemical reactions that function as PID controllers are discussed, with particular focus on individual controllers separately, and the corresponding closed-loop system is analyzed for stochastic controller realizations. Our results show that proportional controllers are effective in buffering protein copy number fluctuations from both noise sources, but this noise suppression comes at the cost of reduced static sensitivity of the output to the input signal. In contrast, integral feedback has no effect on the protein noise level from stochastic expression, but significantly minimizes the impact of external disturbances, particularly when the disturbance comes at low frequencies. Counter-intuitively, integral feedback is found to amplify external disturbances at intermediate frequencies. Next, we discuss the design of a coupled feedforward-feedback biochemical circuit that approximately functions as a derivate controller. Analysis using both analytical methods and Monte Carlo simulations reveals that this derivative controller effectively buffers output fluctuations from bursty stochastic expression, while maintaining the static input-output sensitivity of the open-loop system. In summary, this study provides a systematic stochastic analysis of biochemical controllers, and paves the way for their synthetic design and implementation to minimize deleterious fluctuations in gene product levels.

In the noisy cellular environment, biochemical species such as genes, RNAs and proteins that often occur at low molecular counts, are subject to considerable stochastic fluctuations in copy numbers over time. How cellular biochemical processes function reliably in the face of such randomness is an intriguing fundamental problem. Increasing evidence suggests that random fluctuations (noise) in protein copy numbers play important functional roles, such as driving genetically identical cells to different cell fates. Moreover, many disease states have been attributed to elevated noise levels in specific proteins. Here we systematically investigate design of biochemical systems that function as proportional, integral and derivative-based feedback controllers to suppress protein count fluctuations arising from bursty expression of the protein and external disturbance in protein synthesis. Our results show that different controllers are effective in buffering different noise components, and identify ranges of feedback gain for minimizing deleterious fluctuations in protein levels.

Stability Based on PI Control of Three-Dimensional Chaotic Oscillatory System via DNA Chemical Reaction Networks

The classical proportional integral (PI) controller of SISO linear system is realized by DNA chemical reaction networks (CRNs) in the previous work. Up to now, few works have been done to realize PI controller of chaotic system through DNA CRNs. In this paper, a three-dimensional chaotic oscillatory system and a PI controller of three-dimensional chaotic oscillatory system are proposed by DNA CRNs. The CRNs of chaotic oscillatory system are made up of catalysis modules, degradation module and annihilation module then chemical reaction equations can be compiled into three-dimensional chaotic oscillatory system by the law of mass action to generate chaotic oscillatory signals. The CRNs of PI controller are designed by an integral module, a proportion module and an addition module, which can be compiled into PI controller for stabilizing chaotic oscillatory signals. The simulations of Matlab and Visual DSD are given to show our design achieving the PI control of a three-variable chaotic oscillatory system.

Perfect adaptation in biology

A distinctive feature of many biological systems is their ability to adapt to persistent stimuli or disturbances that would otherwise drive them away from a desirable steady state. The resulting stasis enables organisms to function reliably while being subjected to very different external environments. This perspective concerns a stringent type of biological adaptation, robust perfect adaptation (RPA), that is resilient to certain network and parameter perturbations. As in engineered control systems, RPA requires that the regulating network satisfy certain structural constraints that cannot be avoided. We elucidate these ideas using biological examples from systems and synthetic biology. We then argue that understanding the structural constraints underlying RPA allows us to look past implementation details and offers a compelling means to unravel regulatory biological complexity.

An amplified derepression controller with multisite inhibition and positive feedback

How organisms are able to maintain robust homeostasis has in recent years received increased attention by the use of combined control engineering and kinetic concepts, which led to the discovery of robust controller motifs. While these motifs employ kinetic conditions showing integral feedback and homeostasis for step-wise perturbations, the motifs’ performance differ significantly when exposing them to time dependent perturbations. One type of controller motifs which are able to handle exponentially and even hyperbolically growing perturbations are based on derepression. In these controllers the compensatory reaction, which neutralizes the perturbation, is derepressed, i.e. its reaction rate is increased by the decrease of an inhibitor acting on the compensatory flux. While controllers in this category can deal well with different time-dependent perturbations they have the disadvantage that they break down once the concentration of the regulatory inhibitor becomes too low and the compensatory flux has gained its maximum value. We wondered whether it would be possible to bypass this restriction, while still keeping the advantages of derepression kinetics. In this paper we show how the inclusion of multisite inhibition and the presence of positive feedback loops lead to an amplified controller which is still based on derepression kinetics but without showing the breakdown due to low inhibitor concentrations. By searching for the amplified feedback motif in natural systems, we found it as a part of the plant circadian clock where it is highly interlocked with other feedback loops.

Development of a neuron model based on DNAzyme regulation

Neural networks based on DNA molecular circuits play an important role in molecular information processing and artificial intelligence systems. In fact, some DNA molecular systems can become dynamic units with the assistance of DNAzymes. The complex DNA circuits can spontaneously induce corresponding feedback behaviors when their inputs changed. However, most of the reported DNA neural networks have been implemented by the toehold-mediated strand displacement (TMSD) method. Therefore, it was important to develop a method to build a neural network utilizing the TMSD mechanism and adding a mechanism to account for modulation by DNAzymes. In this study, we designed a model of a DNA neuron controlled by DNAzymes. We proposed an approach based on the DNAzyme modulation of neuronal function, combing two reaction mechanisms: DNAzyme digestion and TMSD. Using the DNAzyme adjustment, each component simulating the characteristics of neurons was constructed. By altering the input and weight of the neuron model, we verified the correctness of the computational function of the neurons. Furthermore, in order to verify the application potential of the neurons in specific functions, a voting machine was successfully implemented. The proposed neuron model regulated by DNAzymes was simple to construct and possesses strong scalability, having great potential for use in the construction of large neural networks.

Ultrasensitive molecular controllers for quasi-integral feedback

Feedback control has enabled the success of automated technologies by mitigating the effects of variability, unknown disturbances, and noise. While it is known that biological feedback loops reduce the impact of noise and help shape kinetic responses, many questions remain about how to design molecular integral controllers. Here, we propose a modular strategy to build molecular quasi-integral feedback controllers, which involves following two design principles. The first principle is to utilize an ultrasensitive response, which determines the gain of the controller and influences the steady-state error. The second is to use a tunable threshold of the ultrasensitive response, which determines the equilibrium point of the system. We describe a reaction network, named brink controller, that satisfies these conditions by combining molecular sequestration and an activation/deactivation cycle. With computational models, we examine potential biological implementations of brink controllers, and we illustrate different example applications.

Frequency switching between oscillatory homeostats and the regulation of p53

Homeostasis is an essential concept to understand the stability of organisms and their adaptive behaviors when coping with external and internal assaults. Many hormones that take part in homeostatic control come in antagonistic pairs, such as glucagon and insulin reflecting the inflow and outflow compensatory mechanisms to control a certain internal variable, such as blood sugar levels. By including negative feedback loops homeostatic controllers can exhibit oscillations with characteristic frequencies. In this paper we demonstrate the associated frequency changes in homeostatic systems when individual controllers -in a set of interlocked feedback loops- gain control in response to environmental changes. Taking p53 as an example, we show how Per2, ATM and Mdm2 feedback loops -interlocked with p53- gain individual control in dependence to the level of DNA damage, and how each of these controllers provide certain functionalities in their regulation of p53. In unstressed cells, the circadian regulator Per2 ensures a basic p53 level to allow its rapid up-regulation in case of DNA damage. When DNA damage occurs the ATM controller increases the level of p53 and defends it towards uncontrolled degradation, which despite DNA damage, would drive p53 to lower values and p53 dysfunction. Mdm2 on its side keeps p53 at a high but sub-apoptotic level to avoid premature apoptosis. However, with on-going DNA damage the Mdm2 set-point is increased by HSP90 and other p53 stabilizers leading finally to apoptosis. An emergent aspect of p53 upregulation during cell stress is the coordinated inhibition of ubiquitin-independent and ubiquitin-dependent degradation reactions. Whether oscillations serve a function or are merely a by-product of the controllers are discussed in view of the finding that homeostatic control of p53, as indicated above, does in principle not require oscillatory homeostats.

In vitro implementation of robust gene regulation in a synthetic biomolecular integral controller

Feedback mechanisms play a critical role in the maintenance of cell homeostasis in the presence of disturbances and uncertainties. Motivated by the need to tune the dynamics and improve the robustness of gene circuits, biological engineers have proposed various designs that mimic natural molecular feedback control mechanisms. However, practical and predictable implementations have proved challenging because of the complexity of synthesis and analysis of complex biomolecular networks. Here, we analyze and experimentally validate a synthetic biomolecular controller executed in vitro. The controller ensures that gene expression rate tracks an externally imposed reference level, and achieves this goal even in the presence of certain kinds of disturbances. Our design relies upon an analog of the well-known principle of integral feedback in control theory. We implement the controller in an Escherichia coli cell-free transcription-translation system, which allows rapid prototyping and implementation. Modeling and theory guide experimental implementation with well-defined operational predictability.

Controlling cell-to-cell variability with synthetic gene circuits

Cell-to-cell variability originating, for example, from the intrinsic stochasticity of gene expression, presents challenges for designing synthetic gene circuits that perform robustly. Conversely, synthetic biology approaches are instrumental in uncovering mechanisms underlying variability in natural systems. With a focus on reducing noise in individual genes, the field has established a broad synthetic toolset. This includes noise control by engineering of transcription and translation mechanisms either individually, or in combination to achieve independent regulation of mean expression and its variability. Synthetic feedback circuits use these components to establish more robust operation in closed-loop, either by drawing on, but also by extending traditional engineering concepts. In this perspective, we argue that major conceptual advances will require new theory of control adapted to biology, extensions from single genes to networks, more systematic considerations of origins of variability other than intrinsic noise, and an exploration of how noise shaping, instead of noise reduction, could establish new synthetic functions or help understanding natural functions.

Homeostatic controllers compensating for growth and perturbations

Cells and organisms have developed homeostatic mechanisms which protect them against a changing environment. How growth and homeostasis interact is still not well understood, but of increasing interest to the molecular and synthetic biology community to recognize and design control circuits which can oppose the diluting effects of cell growth. In this paper we describe the performance of selected negative feedback controllers in response to different applied growth laws and time dependent outflow perturbations of a controlled variable. The approach taken here is based on deterministic mass action kinetics assuming that cell content is instantaneously mixed. All controllers behave ideal in the sense that they for step-wise perturbations in volume and a controlled compound A are able to drive A precisely back to the controllers' theoretical set-points. The applied growth kinetics reflect experimentally observed growth laws, which range from surface to volume ratio growth to linear and exponential growth. Our results show that the kinetic implementation of integral control and the structure of the negative feedback loop are two properties which affect controller performance. Best performance is observed for controllers based on derepression kinetics and controllers with an autocatalytic implementation of integral control. Both are able to defend exponential growth and perturbations, although the autocatalytic controller shows an offset from its theoretical set-point. Controllers with activating signaling using zero-order or bimolecular (antithetic) kinetics for integral control behave very similar but less well. Their performance can be improved by implementing negative feedback structures having repression/derepression steps or by increasing controller aggressiveness. Our results provide a guide what type of feedback structures and integral control kinetics are suitable to oppose the dilution effects by different growth laws and time dependent perturbations on a deterministic level.

Stem cell bioengineering: building from stem cell biology

New fundamental discoveries in stem cell biology have yielded potentially transformative regenerative therapeutics. However, widespread implementation of stem-cell-derived therapeutics remains sporadic. Barriers that impede the development of these therapeutics can be linked to our incomplete understanding of how the regulatory networks that encode stem cell fate govern the development of the complex tissues and organs that are ultimately required for restorative function. Bioengineering tools, strategies and design principles represent core components of the stem cell bioengineering toolbox. Applied to the different layers of complexity present in stem-cell-derived systems - from gene regulatory networks in single stem cells to the systemic interactions of stem-cell-derived organs and tissues - stem cell bioengineering can address existing challenges and advance regenerative medicine and cellular therapies.

Synthetic control systems for high performance gene expression in mammalian cells

Tunable induction of gene expression is an essential tool in biology and biotechnology. In spite of that, current induction systems often exhibit unpredictable behavior and performance shortcomings, including high sensitivity to transactivator dosage and plasmid take-up variation, and excessive consumption of cellular resources. To mitigate these limitations, we introduce here a novel family of gene expression control systems of varying complexity with significantly enhanced performance. These include: (i) an incoherent feedforward circuit that exhibits output tunability and robustness to plasmid take-up variation; (ii) a negative feedback circuit that reduces burden and provides robustness to transactivator dosage variability; and (iii) a new hybrid circuit integrating negative feedback and incoherent feedforward that combines the benefits of both. As with endogenous circuits, the complexity of our genetic controllers is not gratuitous, but is the necessary outcome of more stringent performance requirements. We demonstrate the benefits of these controllers in two applications. In a culture of CHO cells for protein manufacturing, the circuits result in up to a 2.6-fold yield improvement over a standard system. In human-induced pluripotent stem cells they enable precisely regulated expression of an otherwise poorly tolerated gene of interest, resulting in a significant increase in the viability of the transfected cells.

Realizing 'integral control' in living cells: how to overcome leaky integration due to dilution?

A major problem in the design of synthetic genetic circuits is robustness to perturbations and uncertainty. Because of this, there have been significant efforts in recent years in finding approaches to implement integral control in genetic circuits. Integral controllers have the unique ability to make the output of a process adapt perfectly to disturbances. However, implementing an integral controller is challenging in living cells. This is because a key aspect of any integral controller is a 'memory' element that stores the accumulation (integral) of the error between the output and its desired set-point. The ability to realize such a memory element in living cells is fundamentally challenged by the fact that all biomolecules dilute as cells grow, resulting in a 'leaky' memory that gradually fades away. As a consequence, the adaptation property is lost. Here, we propose a general principle for designing integral controllers such that the performance is practically unaffected by dilution. In particular, we mathematically prove that if the reactions implementing the integral controller are all much faster than dilution, then the adaptation error due to integration leakiness becomes negligible. We exemplify this design principle with two synthetic genetic circuits aimed at reaching adaptation of gene expression to fluctuations in cellular resources. Our results provide concrete guidance on the biomolecular processes that are most appropriate for implementing integral controllers in living cells.

A universal biomolecular integral feedback controller for robust perfect adaptation

Homeostasis is a recurring theme in biology that ensures that regulated variables robustly-and in some systems, completely-adapt to environmental perturbations. This robust perfect adaptation feature is achieved in natural circuits by using integral control, a negative feedback strategy that performs mathematical integration to achieve structurally robust regulation^1,2. Despite its benefits, the synthetic realization of integral feedback in living cells has remained elusive owing to the complexity of the required biological computations. Here we prove mathematically that there is a single fundamental biomolecular controller topology³ that realizes integral feedback and achieves robust perfect adaptation in arbitrary intracellular networks with noisy dynamics. This adaptation property is guaranteed both for the population-average and for the time-average of single cells. On the basis of this concept, we genetically engineer a synthetic integral feedback controller in living cells⁴ and demonstrate its tunability and adaptation properties. A growth-rate control application in Escherichia coli shows the intrinsic capacity of our integral controller to deliver robustness and highlights its potential use as a versatile controller for regulation of biological variables in uncertain networks. Our results provide conceptual and practical tools in the area of cybergenetics^3,5, for engineering synthetic controllers that steer the dynamics of living systems^3-9.

Adaptive information processing of network modules to dynamic and spatial stimuli

BACKGROUND

Adaptation and homeostasis are basic features of information processing in cells and seen in a broad range of contexts. Much of the current understanding of adaptation in network modules/motifs is based on their response to simple stimuli. Recently, there have also been studies of adaptation in dynamic stimuli. However a broader synthesis of how different circuits of adaptation function, and which circuits enable a broader adaptive behaviour in classes of more complex and spatial stimuli is largely missing.

RESULTS

We study the response of a variety of adaptive circuits to time-varying stimuli such as ramps, periodic stimuli and static and dynamic spatial stimuli. We find that a variety of responses can be seen in ramp stimuli, making this a basis for discriminating between even similar circuits. We also find that a number of circuits adapt exactly to ramp stimuli, and dissect these circuits to pinpoint what characteristics (architecture, feedback, biochemical aspects, information processing ingredients) allow for this. These circuits include incoherent feedforward motifs, inflow-outflow motifs and transcritical circuits. We find that changes in location in such circuits where a signal acts can result in non-adaptive behaviour in ramps, even though the location was associated with exact adaptation in step stimuli. We also demonstrate that certain augmentations of basic inflow-outflow motifs can alter the behaviour of the circuit from exact adaptation to non-adaptive behaviour. When subject to periodic stimuli, some circuits (inflow-outflow motifs and transcritical circuits) are able to maintain an average output independent of the characteristics of the input. We build on this to examine the response of adaptive circuits to static and dynamic spatial stimuli. We demonstrate how certain circuits can exhibit a graded response in spatial static stimuli with an exact maintenance of the spatial mean-value. Distinct features which emerge from the consideration of dynamic spatial stimuli are also discussed. Finally, we also build on these results to show how different circuits which show any combination of presence or absence of exact adaptation in ramps, exact mainenance of time average output in periodic stimuli and exact maintenance of spatial average of output in static spatial stimuli may be realized.

CONCLUSIONS

By studying a range of network circuits/motifs on one hand and a range of stimuli on the other, we isolate characteristics of these circuits (structural) which enable different degrees of exact adaptive and homeostatic behaviour in such stimuli, how they may be combined, and also identify cases associated with non-homeostatic behaviour. We also reveal constraints associated with locations where signals may act to enable homeostatic behaviour and constraints associated with augmentations of circuits. This consideration of multiple experimentally/naturally relevant stimuli along with circuits of adaptation of relevance in natural and engineered biology, provides a platform for deepening our understanding of adaptive and homeostatic behaviour in natural systems, bridging the gap between models of adaptation and experiments and in engineering homeostatic synthetic circuits.

A Framework for Engineering Stress Resilient Plants Using Genetic Feedback Control and Regulatory Network Rewiring

Crop disease leads to significant waste worldwide, both pre- and postharvest, with subsequent economic and sustainability consequences. Disease outcome is determined both by the plants' response to the pathogen and by the ability of the pathogen to suppress defense responses and manipulate the plant to enhance colonization. The defense response of a plant is characterized by significant transcriptional reprogramming mediated by underlying gene regulatory networks, and components of these networks are often targeted by attacking pathogens. Here, using gene expression data from Botrytis cinerea-infected Arabidopsis plants, we develop a systematic approach for mitigating the effects of pathogen-induced network perturbations, using the tools of synthetic biology. We employ network inference and system identification techniques to build an accurate model of an Arabidopsis defense subnetwork that contains key genes determining susceptibility of the plant to the pathogen attack. Once validated against time-series data, we use this model to design and test perturbation mitigation strategies based on the use of genetic feedback control. We show how a synthetic feedback controller can be designed to attenuate the effect of external perturbations on the transcription factor CHE in our subnetwork. We investigate and compare two approaches for implementing such a controller biologically-direct implementation of the genetic feedback controller, and rewiring the regulatory regions of multiple genes-to achieve the network motif required to implement the controller. Our results highlight the potential of combining feedback control theory with synthetic biology for engineering plants with enhanced resilience to environmental stress.

Antithetic proportional-integral feedback for reduced variance and improved control performance of stochastic reaction networks

The ability of a cell to regulate and adapt its internal state in response to unpredictable environmental changes is called homeostasis and this ability is crucial for the cell's survival and proper functioning. Understanding how cells can achieve homeostasis, despite the intrinsic noise or randomness in their dynamics, is fundamentally important for both systems and synthetic biology. In this context, a significant development is the proposed antithetic integral feedback (AIF) motif, which is found in natural systems, and is known to ensure robust perfect adaptation for the mean dynamics of a given molecular species involved in a complex stochastic biomolecular reaction network. From the standpoint of applications, one drawback of this motif is that it often leads to an increased cell-to-cell heterogeneity or variance when compared to a constitutive (i.e. open-loop) control strategy. Our goal in this paper is to show that this performance deterioration can be countered by combining the AIF motif and a negative feedback strategy. Using a tailored moment closure method, we derive approximate expressions for the stationary variance for the controlled network that demonstrate that increasing the strength of the negative feedback can indeed decrease the variance, sometimes even below its constitutive level. Numerical results verify the accuracy of these results and we illustrate them by considering three biomolecular networks with two types of negative feedback strategies. Our computational analysis indicates that there is a trade-off between the speed of the settling-time of the mean trajectories and the stationary variance of the controlled species; i.e. smaller variance is associated with larger settling-time.

Mathematical Modeling of RNA-Based Architectures for Closed Loop Control of Gene Expression

Feedback allows biological systems to control gene expression precisely and reliably, even in the presence of uncertainty, by sensing and processing environmental changes. Taking inspiration from natural architectures, synthetic biologists have engineered feedback loops to tune the dynamics and improve the robustness and predictability of gene expression. However, experimental implementations of biomolecular control systems are still far from satisfying performance specifications typically achieved by electrical or mechanical control systems. To address this gap, we present mathematical models of biomolecular controllers that enable reference tracking, disturbance rejection, and tuning of the temporal response of gene expression. These controllers employ RNA transcriptional regulators to achieve closed loop control where feedback is introduced via molecular sequestration. Sensitivity analysis of the models allows us to identify which parameters influence the transient and steady state response of a target gene expression process, as well as which biologically plausible parameter values enable perfect reference tracking. We quantify performance using typical control theory metrics to characterize response properties and provide clear selection guidelines for practical applications. Our results indicate that RNA regulators are well-suited for building robust and precise feedback controllers for gene expression. Additionally, our approach illustrates several quantitative methods useful for assessing the performance of biomolecular feedback control systems.

Complexity-aware simple modeling

Mathematical models continue to be essential for deepening our understanding of biology. On one extreme, simple or small-scale models help delineate general biological principles. However, the parsimony of detail in these models as well as their assumption of modularity and insulation make them inaccurate for describing quantitative features. On the other extreme, large-scale and detailed models can quantitatively recapitulate a phenotype of interest, but have to rely on many unknown parameters, making them often difficult to parse mechanistically and to use for extracting general principles. We discuss some examples of a new approach-complexity-aware simple modeling-that can bridge the gap between the small-scale and large-scale approaches.

Perfect Adaptation and Optimal Equilibrium Productivity in a Simple Microbial Biofuel Metabolic Pathway Using Dynamic Integral Control

The production of complex biomolecules by genetically engineered organisms is one of the most promising applications of metabolic engineering and synthetic biology. To obtain processes with high productivity, it is therefore crucial to design and implement efficient dynamic in vivo regulation strategies. We consider here the microbial biofuel production model of Dunlop et al. (2010) for which we demonstrate that an antithetic dynamic integral control strategy can achieve robust perfect adaptation for the intracellular biofuel concentration in the presence of poorly known network parameters and implementation errors in certain rate parameters of the controller. We also show that the maximum equilibrium extracellular biofuel productivity is fully defined by some of the network parameters and, in this respect, it can only be achieved when all the corresponding parameters are perfectly known. Since this optimum is a network property, it cannot be improved by the use of any controller that measures the intracellular biofuel concentration and acts on the production of pump proteins. Additional intrinsic fundamental properties for the process are also unveiled, the most important ones being the existence of a conservation relation between the productivity and the toxicity, a low sensitivity of the optimal productivity with respect to a poor implementation of the set-point for the intracellular biofuel, and a strong intrinsic robustness property of the optimal productivity with respect to poorly known parameters. Taken together, these results demonstrate that a high and robust equilibrium rate of production for the extracellular biofuel can be achieved when the parameters of the model are poorly known and those of the controllers are poorly implemented. Finally, several advantages of the proposed dynamic strategy over a static one are also emphasized.

Realizing 'integral control' in living cells: how to overcome leaky integration due to dilution?

A major problem in the design of synthetic genetic circuits is robustness to perturbations and uncertainty. Because of this, there have been significant efforts in recent years in finding approaches to implement integral control in genetic circuits. Integral controllers have the unique ability to make the output of a process adapt perfectly to disturbances. However, implementing an integral controller is challenging in living cells. This is because a key aspect of any integral controller is a 'memory' element that stores the accumulation (integral) of the error between the output and its desired set-point. The ability to realize such a memory element in living cells is fundamentally challenged by the fact that all biomolecules dilute as cells grow, resulting in a 'leaky' memory that gradually fades away. As a consequence, the adaptation property is lost. Here, we propose a general principle for designing integral controllers such that the performance is practically unaffected by dilution. In particular, we mathematically prove that if the reactions implementing the integral controller are all much faster than dilution, then the adaptation error due to integration leakiness becomes negligible. We exemplify this design principle with two synthetic genetic circuits aimed at reaching adaptation of gene expression to fluctuations in cellular resources. Our results provide concrete guidance on the biomolecular processes that are most appropriate for implementing integral controllers in living cells.

Design of an embedded inverse-feedforward biomolecular tracking controller for enzymatic reaction processes

Feedback control is widely used in chemical engineering to improve the performance and robustness of chemical processes. Feedback controllers require a 'subtractor' that is able to compute the error between the process output and the reference signal. In the case of embedded biomolecular control circuits, subtractors designed using standard chemical reaction network theory can only realise one-sided subtraction, rendering standard controller design approaches inadequate. Here, we show how a biomolecular controller that allows tracking of required changes in the outputs of enzymatic reaction processes can be designed and implemented within the framework of chemical reaction network theory. The controller architecture employs an inversion-based feedforward controller that compensates for the limitations of the one-sided subtractor that generates the error signals for a feedback controller. The proposed approach requires significantly fewer chemical reactions to implement than alternative designs, and should have wide applicability throughout the fields of synthetic biology and biological engineering.

A Blueprint for a Synthetic Genetic Feedback Controller to Reprogram Cell Fate

To artificially reprogram cell fate, experimentalists manipulate the gene regulatory networks (GRNs) that maintain a cell's phenotype. In practice, reprogramming is often performed by constant overexpression of specific transcription factors (TFs). This process can be unreliable and inefficient. Here, we address this problem by introducing a new approach to reprogramming based on mathematical analysis. We demonstrate that reprogramming GRNs using constant overexpression may not succeed in general. Instead, we propose an alternative reprogramming strategy: a synthetic genetic feedback controller that dynamically steers the concentration of a GRN's key TFs to any desired value. The controller works by adjusting TF expression based on the discrepancy between desired and actual TF concentrations. Theory predicts that this reprogramming strategy is guaranteed to succeed, and its performance is independent of the GRN's structure and parameters, provided that feedback gain is sufficiently high. As a case study, we apply the controller to a model of induced pluripotency in stem cells.