Modelling periodic oscillation in gene regulatory networks by cyclic feedback systems

In this paper, we develop a new methodology to analyze and design periodic oscillators of biological networks, in particular gene regulatory networks with multiple genes, proteins and time delays, by using negative cyclic feedback systems. We show that negative cyclic feedback networks have no stable equilibria but stable periodic orbits when certain conditions are satisfied. Specifically, we first prove the basic properties of the biological networks composed of cyclic feedback loops, and then extend our results to general cyclic feedback network with less restriction, thereby making our theoretical analysis and design of oscillators easy to implement, even for large-scale systems. Finally, we use one circadian network formed by a period protein (PER) and per mRNA, and one biologically plausible synthetic gene network, to demonstrate the theoretical results. Since there is less restriction on the network structure, the results of this paper can be expected to apply to a wide variety of areas on modelling, analyzing and designing of biological systems.     

Tuning genetic clocks employing DNA binding sites

Periodic oscillations play a key role in cell physiology from the cell cycle to circadian clocks. The interplay of positive and negative feedback loops among genes and proteins is ubiquitous in these networks. Often, delays in a negative feedback loop and/or degradation rates are a crucial mechanism to obtain sustained oscillations. How does nature control delays and kinetic rates in feedback networks? Known mechanisms include proper selection of the number of steps composing a feedback loop and alteration of protease activity, respectively. Here, we show that a remarkably simple means to control both delays and effective kinetic rates is the employment of DNA binding sites. We illustrate this design principle on a widely studied activator-repressor clock motif, which is ubiquitous in natural systems. By suitably employing DNA target sites for the activator and/or the repressor, one can switch the clock “on” and “off” and precisely tune its period to a desired value. Our study reveals a design principle to engineer dynamic behavior in biomolecular networks, which may be largely exploited by natural systems and employed for the rational design of synthetic circuits.     

Coupling oscillations and switches in genetic networks

Switches (bistability) and oscillations (limit cycle) are omnipresent in biological networks. Synthetic genetic networks producing bistability and oscillations have been designed and constructed experimentally. However, in real biological systems, regulatory circuits are usually interconnected and the dynamics of those complex networks is often richer than the dynamics of simple modules. Here we couple the genetic Toggle switch and the Repressilator, two prototypic systems exhibiting bistability and oscillations, respectively. We study two types of coupling. In the first type, the bistable switch is under the control of the oscillator. Numerical simulation of this system allows us to determine the conditions under which a periodic switch between the two stable steady states of the Toggle switch occurs. In addition we show how birhythmicity characterized by the coexistence of two stable small-amplitude limit cycles, can easily be obtained in the system. In the second type of coupling, the oscillator is placed under the control of the Toggleswitch. Numerical simulation of this system shows that this construction could for example be exploited to generate a permanent transition from a stable steady state to self-sustained oscillations (and vice versa) after a transient external perturbation. Those results thus describe qualitative dynamical behaviors that can be generated through the coupling of two simple network modules. These results differ from the dynamical properties resulting from interlocked feedback loops systems in which a given variable is involved at the same time in both positive and negative feedbacks. Finally the models described here may be of interest in synthetic biology, as they give hints on how the coupling should be designed to get the required properties.     

Robust,tunable genetic memory from protein sequestration combined with positive feedback

Natural regulatory networks contain many interacting components that allow for fine-tuning of switching and memory properties. Building simple bistable switches, synthetic biologists have learned the design principles of complex natural regulatory networks. However, most switches constructed so far are so simple (e.g. comprising two regulators) that they are functional only within a limited parameter range. Here, we report the construction of robust, tunable bistable switches in Escherichia coli using three heterologous protein regulators (ExsADC) that are sequestered into an inactive complex through a partner swapping mechanism. On the basis of mathematical modeling, we accurately predict and experimentally verify that the hysteretic region can be fine-tuned by controlling the interactions of the ExsADC regulatory cascade using the third member ExsC as a tuning knob. Additionally, we confirm that a dual-positive feedback switch can markedly increase the hysteretic region, compared to its single-positive feedback counterpart. The dual-positive feedback switch displays bistability over a 10⁶-fold range of inducer concentrations, to our knowledge, the largest range reported so far. This work demonstrates the successful interlocking of sequestration-based ultrasensitivity and positive feedback, a design principle that can be applied to the construction of robust, tunable, and predictable genetic programs to achieve increasingly sophisticated biological behaviors.     

A computational method for the investigation of multistable systems and its application to genetic switches

Background

Genetic switches exhibit multistability, form the basis of epigenetic memory, and are found in natural decision making systems, such as cell fate determination in developmental pathways. Synthetic genetic switches can be used for recording the presence of different environmental signals, for changing phenotype using synthetic inputs and as building blocks for higher-level sequential logic circuits. Understanding how multistable switches can be constructed and how they function within larger biological systems is therefore key to synthetic biology.

Results

Here we present a new computational tool, called StabilityFinder, that takes advantage of sequential Monte Carlo methods to identify regions of parameter space capable of producing multistable behaviour, while handling uncertainty in biochemical rate constants and initial conditions. The algorithm works by clustering trajectories in phase space, and iteratively minimizing a distance metric. Here we examine a collection of models of genetic switches, ranging from the deterministic Gardner toggle switch to stochastic models containing different positive feedback connections. We uncover the design principles behind making bistable, tristable and quadristable switches, and find that rate of gene expression is a key parameter. We demonstrate the ability of the framework to examine more complex systems and examine the design principles of a three gene switch. Our framework allows us to relax the assumptions that are often used in genetic switch models and we show that more complex abstractions are still capable of multistable behaviour.

Conclusions

Our results suggest many ways in which genetic switches can be enhanced and offer designs for the construction of novel switches. Our analysis also highlights subtle changes in correlation of experimentally tunable parameters that can lead to bifurcations in deterministic and stochastic systems. Overall we demonstrate that StabilityFinder will be a valuable tool in the future design and construction of novel gene networks.

     

Synthetic conversion of a graded receptor signal into a tunable,reversible switch

The ability to engineer an all‐or‐none cellular response to a given signaling ligand is important in applications ranging from biosensing to tissue engineering. However, synthetic gene network ‘switches’ have been limited in their applicability and tunability due to their reliance on specific components to function. Here, we present a strategy for reversible switch design that instead relies only on a robust, easily constructed network topology with two positive feedback loops and we apply the method to create highly ultrasensitive (n_H>20), bistable cellular responses to a synthetic ligand/receptor complex. Independent modulation of the two feedback strengths enables rational tuning and some decoupling of steady‐state (ultrasensitivity, signal amplitude, switching threshold, and bistability) and kinetic (rates of system activation and deactivation) response properties. Our integrated computational and synthetic biology approach elucidates design rules for building cellular switches with desired properties, which may be of utility in engineering signal‐transduction pathways.     

The regulatory circuits for hysteretic switching in cellular signal transduction pathways

Hysteresis can be found in many physical systems, and a hysteretic switch has been used for various mechanical and electrical systems. Such a hysteretic switch can be created by using a single positive feedback loop, as often used in engineering systems. It is, however, intriguing that various cellular signaling systems use coupled positive feedback loops to implement the hysteretic switch. A question then arises about the advantage of using coupled positive feedback loops instead of simple isolated positive feedback for an apparently equivalent hysteretic switch. Through mathematical simulations, we determined that cellular systems with coupled positive feedback loops show enhanced hysteretic switching, and can thereby make a more reliable decision under conditions of noisy signaling. As most intracellular processes are accompanied by intrinsic noise, important cellular decisions such as differentiation and apoptosis need to be highly robust to such noises. The coupled positive feedback loops might have been evolutionarily acquired to enable correct cell fate decisions to be made through enhanced hysteretic switching in noisy cellular environments.     

Dual selection of a genetic switch by a single selection marker

Forward engineering of synthetic genetic circuits in living cells is expected to deliver various applications in biotechnology and medicine and to provide valuable insights into the design principles of natural gene networks. However, lack of biochemical data and complexity of biological environment complicate rational design of such circuits based on quantitative simulation. Previously, we have shown that directed evolution can complement our weakness in designing genetic circuits by screening or selecting functional circuits from a large pool of nonfunctional ones. Here we describe a dual selection strategy that allows selection of both ON and OFF states of genetic circuits using tetA as a single selection marker. We also describe a successful demonstration of a genetic switch selection from a 2000-fold excess background of nonfunctional switches in three rounds of iterative selection. The dual selection system is more robust than the previously reported selection system employing three genes, with no observed false positive mutants during the simulated selections.     

Modelling periodic oscillation of biological systems with multiple timescale networks

In this paper, we aim to develop a new methodology to model and design periodic oscillators of biological networks, in particular gene regulatory networks with multiple genes, proteins and time delays, by using multiple timescale networks (MTN). Fast reactions constitute a positive feedback-loop network (PFN), while slow reactions consist of a cyclic feedback-loop network (CFN), in MTN. Multiple timescales are exploited to simplify models according to singular perturbation theory. We show that a MTN has no stable equilibrium but stable periodic orbits when certain conditions are satisfied. Specifically, we first prove the basic properties of MTNs with only one PFN, and then generalise the result to MTNs with multiple PFNs. Finally, we design a biologically plausible gene regulatory network by the cI and Lac genes, to demonstrate the theoretical results. Since there is less restriction on the network structure of a MTN, it can be expected to apply to a wide variety of areas on the modelling, analysing and designing of biological systems.     

Synthetic genetic circuits for programmable biological functionalities

Living organisms evolve complex genetic networks to interact with the environment. Due to the rapid development of synthetic biology, various modularized genetic parts and units have been identified from these networks. They have been employed to construct synthetic genetic circuits, including toggle switches, oscillators, feedback loops and Boolean logic gates. Building on these circuits, complex genetic machines with capabilities in programmable decision-making could be created. Consequently, these accomplishments have led to novel applications, such as dynamic and autonomous modulation of metabolic networks, directed evolution of biological units, remote and targeted diagnostics and therapies, as well as biological containment methods to prevent release of engineered microorganisms and genetic materials. Herein, we outline the principles in genetic circuit design that have initiated a new chapter in transforming concepts to realistic applications. The features of modularized building blocks and circuit architecture that facilitate realization of circuits for a variety of novel applications are discussed. Furthermore, recent advances and challenges in employing genetic circuits to impart microorganisms with distinct and programmable functionalities are highlighted. We envision that this review gives new insights into the design of synthetic genetic circuits and offers a guideline for the implementation of different circuits in various aspects of biotechnology and bioengineering.     

Coupled positive and negative feedback circuits form an essential building block of cellular signaling pathways

Cellular circuits have positive and negative feedback loops that allow them to respond properly to noisy external stimuli. It is intriguing that such feedback loops exist in many cases in a particular form of coupled positive and negative feedback loops with different time delays. As a result of our mathematical simulations and investigations into various experimental evidences, we found that such coupled feedback circuits can rapidly turn on a reaction to a proper stimulus, robustly maintain its status, and immediately turn off the reaction when the stimulus disappears. In other words, coupled feedback loops enable cellular systems to produce perfect responses to noisy stimuli with respect to signal duration and amplitude. This suggests that coupled positive and negative feedback loops form essential signal transduction motifs in cellular signaling systems.     

Distinctive topologies of partner-switching signaling networks correlate with their physiological roles

Regulatory networks controlling bacterial gene expression often evolve from common origins and share homologous proteins and similar network motifs. However, when functioning in different physiological contexts, these motifs may be re-arranged with different topologies that significantly affect network performance. Here we analyze two related signaling networks in the bacterium Bacillus subtilis in order to assess the consequences of their different topologies, with the aim of formulating design principles applicable to other systems. These two networks control the activities of the general stress response factor sigma(B) and the first sporulation-specific factor sigma(F). Both networks have at their core a "partner-switching" mechanism, in which an anti-sigma factor forms alternate complexes either with the sigma factor, holding it inactive, or with an anti-anti-sigma factor, thereby freeing sigma. However, clear differences in network structure are apparent: the anti-sigma factor for sigma(F) forms a long-lived, "dead-end" complex with its anti-anti-sigma factor and ADP, whereas the genes encoding sigma(B) and its network partners lie in a sigma(B)-controlled operon, resulting in positive and negative feedback loops. We constructed mathematical models of both networks and examined which features were critical for the performance of each design. The sigma(F) model predicts that the self-enhancing formation of the dead-end complex transforms the network into a largely irreversible hysteretic switch; the simulations reported here also demonstrate that hysteresis and slow turn off kinetics are the only two system properties associated with this complex formation. By contrast, the sigma(B) model predicts that the positive and negative feedback loops produce graded, reversible behavior with high regulatory capacity and fast response time. Our models demonstrate how alterations in network design result in different system properties that correlate with regulatory demands. These design principles agree with the known or suspected roles of similar networks in diverse bacteria.     

How robust are switches in intracellular signaling cascades?

Since all-or-none decisions of the cell are controlled by extracellular signals, cells have biochemical switches within their intracellular signaling networks. Central elements of these switches are multisite phosphorylation, enzymic saturation, and amplification by cascades. Moreover, positive feedback can contribute to switch-like behavior termed also ultrasensitivity. Here we analyse the robustness of these mechanisms exemplified by models of the three-molecule MAPK-cascade and the single-molecule Goldbeter-Koshland switch. We show that the ultrasensitivity in the MAPK-cascades is more robust against changes of the kinetic parameters than the Goldbeter-Koshland switch. If multiple parameters are changed randomly, the effects of parameter changes can compensate each other in the cascade leading to a remarkable robustness of the switch-like behavior. The different degrees of robustness can be traced back to the different mechanisms of generating ultrasensitivity. While in the Goldbeter-Koshland switch the saturation of the enzymes are crucial, in the MAPK-cascade the adjustment of working ranges determines the ultrasensitivity. Our results indicate that amplification of ultrasensitivity in cascades and multisite phosphorylation might be a design principle to achieve robust switches.