Design Principles of a Genetic Alarm Clock

Turning genes on and off is a mechanism by which cells and tissues make phenotypic decisions. Gene network motifs capable of supporting two or more steady states and thereby providing cells with a plurality of possible phenotypes are referred to as genetic switches. Modeled on the bases of naturally occurring genetic networks, synthetic biologists have successfully constructed artificial switches, thus opening a door to new possibilities for improvement of the known, but also the design of new synthetic genetic circuits. One of many obstacles to overcome in such efforts is to understand and hence control intrinsic noise which is inherent in all biological systems. For some motifs the noise is negligible; for others, fluctuations in the particle number can be comparable to its average. Due to their slowed dynamics, motifs with positive autoregulation tend to be highly sensitive to fluctuations of their chemical environment and are in general very noisy, especially during transition (switching). In this article we use stochastic simulations (Gillespie algorithm) to model such a system, in particular a simple bistable motif consisting of a single gene with positive autoregulation. Due to cooperativety, the dynamical behavior of this kind of motif is reminiscent of an alarm clock – the gene is (nearly) silent for some time after it is turned on and becomes active very suddenly. We investigate how these sudden transitions are affected by noise and show that under certain conditions accurate timing can be achieved. We also examine how promoter complexity influences the accuracy of this timing mechanism.     

Stochastic Bistability and Bifurcation in a Mesoscopic Signaling System with Autocatalytic Kinase

Bistability is a nonlinear phenomenon widely observed in nature including in biochemical reaction networks. Deterministic chemical kinetics studied in the past has shown that bistability occurs in systems with strong (cubic) nonlinearity. For certain mesoscopic, weakly nonlinear (quadratic) biochemical reaction systems in a small volume, however, stochasticity can induce bistability and bifurcation that have no macroscopic counterpart. We report the simplest yet known reactions involving driven phosphorylation-dephosphorylation cycle kinetics with autocatalytic kinase. We show that the noise-induced phenomenon is correlated with free energy dissipation and thus conforms with the open-chemical system theory. A previous reported noise-induced bistability in futile cycles is found to have originated from the kinase synchronization in a bistable system with slow transitions, as reported here.     

Characterizing multistationarity regimes in biochemical reaction networks

Switch like responses appear as common strategies in the regulation of cellular systems. Here we present a method to characterize bistable regimes in biochemical reaction networks that can be of use to both direct and reverse engineering of biological switches. In the design of a synthetic biological switch, it is important to study the capability for bistability of the underlying biochemical network structure. Chemical Reaction Network Theory (CRNT) may help at this level to decide whether a given network has the capacity for multiple positive equilibria, based on their structural properties. However, in order to build a working switch, we also need to ensure that the bistability property is robust, by studying the conditions leading to the existence of two different steady states. In the reverse engineering of biological switches, knowledge collected about the bistable regimes of the underlying potential model structures can contribute at the model identification stage to a drastic reduction of the feasible region in the parameter space of search. In this work, we make use and extend previous results of the CRNT, aiming not only to discriminate whether a biochemical reaction network can exhibit multiple steady states, but also to determine the regions within the whole space of parameters capable of producing multistationarity. To that purpose we present and justify a condition on the parameters of biochemical networks for the appearance of multistationarity, and propose an efficient and reliable computational method to check its satisfaction through the parameter space.     

Antagonism and bistability in protein interaction networks

A protein interaction network (PIN) is a set of proteins that modulate one another's activities by regulated synthesis and degradation, by reversible binding to form complexes, and by catalytic reactions (e.g., phosphorylation and dephosphorylation). Most PINs are so complex that their dynamical characteristics cannot be deduced accurately by intuitive reasoning alone. To predict the properties of such networks, many research groups have turned to mathematical models (differential equations based on standard biochemical rate laws, e.g., mass-action, Michaelis-Menten, Hill). When using Michaelis-Menten rate expressions to model PINs, care must be exercised to avoid making inconsistent assumptions about enzyme-substrate complexes. We show that an appealingly simple model of a PIN that functions as a bistable switch is compromised by neglecting enzyme-substrate intermediates. When the neglected intermediates are put back into the model, bistability of the switch is lost. The theory of chemical reaction networks predicts that bistability can be recovered by adding specific reaction channels to the molecular mechanism. We explore two very different routes to recover bistability. In both cases, we show how to convert the original 'phenomenological' model into a consistent set of mass-action rate laws that retains the desired bistability properties. Once an equivalent model is formulated in terms of elementary chemical reactions, it can be simulated accurately either by deterministic differential equations or by Gillespie's stochastic simulation algorithm.     

A Fundamental Trade-off in Covalent Switching and Its Circumvention by Enzyme Bifunctionality in Glucose Homeostasis

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

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.     

Bistability: Requirements on Cell-Volume,Protein Diffusion,and Thermodynamics

Bistability is considered wide-spread among bacteria and eukaryotic cells, useful e.g. for enzyme induction, bet hedging, and epigenetic switching. However, this phenomenon has mostly been described with deterministic dynamic or well-mixed stochastic models. Here, we map known biological bistable systems onto the well-characterized biochemical Schlögl model, using analytical calculations and stochastic spatiotemporal simulations. In addition to network architecture and strong thermodynamic driving away from equilibrium, we show that bistability requires fine-tuning towards small cell volumes (or compartments) and fast protein diffusion (well mixing). Bistability is thus fragile and hence may be restricted to small bacteria and eukaryotic nuclei, with switching triggered by volume changes during the cell cycle. For large volumes, single cells generally loose their ability for bistable switching and instead undergo a first-order phase transition.     

Multistable Decision Switches for Flexible Control of Epigenetic Differentiation

It is now recognized that molecular circuits with positive feedback can induce two different gene expression states (bistability) under the very same cellular conditions. Whether, and how, cells make use of the coexistence of a larger number of stable states (multistability) is however largely unknown. Here, we first examine how autoregulation, a common attribute of genetic master regulators, facilitates multistability in two-component circuits. A systematic exploration of these modules'' parameter space reveals two classes of molecular switches, involving transitions in bistable (progression switches) or multistable (decision switches) regimes. We demonstrate the potential of decision switches for multifaceted stimulus processing, including strength, duration, and flexible discrimination. These tasks enhance response specificity, help to store short-term memories of recent signaling events, stabilize transient gene expression, and enable stochastic fate commitment. The relevance of these circuits is further supported by biological data, because we find them in numerous developmental scenarios. Indeed, many of the presented information-processing features of decision switches could ultimately demonstrate a more flexible control of epigenetic differentiation.     

A modular approach for modeling the cell cycle based on functional response curves

Modeling biochemical reactions by means of differential equations often results in systems with a large number of variables and parameters. As this might complicate the interpretation and generalization of the obtained results, it is often desirable to reduce the complexity of the model. One way to accomplish this is by replacing the detailed reaction mechanisms of certain modules in the model by a mathematical expression that qualitatively describes the dynamical behavior of these modules. Such an approach has been widely adopted for ultrasensitive responses, for which underlying reaction mechanisms are often replaced by a single Hill function. Also time delays are usually accounted for by using an explicit delay in delay differential equations. In contrast, however, S-shaped response curves, which by definition have multiple output values for certain input values and are often encountered in bistable systems, are not easily modeled in such an explicit way. Here, we extend the classical Hill function into a mathematical expression that can be used to describe both ultrasensitive and S-shaped responses. We show how three ubiquitous modules (ultrasensitive responses, S-shaped responses and time delays) can be combined in different configurations and explore the dynamics of these systems. As an example, we apply our strategy to set up a model of the cell cycle consisting of multiple bistable switches, which can incorporate events such as DNA damage and coupling to the circadian clock in a phenomenological way.     

Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity.     

Dynamic bistable switches enhance robustness and accuracy of cell cycle transitions

Bistability is a common mechanism to ensure robust and irreversible cell cycle transitions. Whenever biological parameters or external conditions change such that a threshold is crossed, the system abruptly switches between different cell cycle states. Experimental studies have uncovered mechanisms that can make the shape of the bistable response curve change dynamically in time. Here, we show how such a dynamically changing bistable switch can provide a cell with better control over the timing of cell cycle transitions. Moreover, cell cycle oscillations built on bistable switches are more robust when the bistability is modulated in time. Our results are not specific to cell cycle models and may apply to other bistable systems in which the bistable response curve is time-dependent.     

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.     

The Eps1p Protein Disulfide Isomerase Conserves Classic Thioredoxin Superfamily Amino Acid Motifs but Not Their Functional Geometries

The widespread thioredoxin superfamily enzymes typically share the following features: a characteristic α-β fold, the presence of a Cys-X-X-Cys (or Cys-X-X-Ser) redox-active motif, and a proline in the cis configuration abutting the redox-active site in the tertiary structure. The Cys-X-X-Cys motif is at the solvent-exposed amino terminus of an α-helix, allowing the first cysteine to engage in nucleophilic attack on substrates, or substrates to attack the Cys-X-X-Cys disulfide, depending on whether the enzyme functions to reduce, isomerize, or oxidize its targets. We report here the X-ray crystal structure of an enzyme that breaks many of our assumptions regarding the sequence-structure relationship of thioredoxin superfamily proteins. The yeast Protein Disulfide Isomerase family member Eps1p has Cys-X-X-Cys motifs and proline residues at the appropriate primary structural positions in its first two predicted thioredoxin-fold domains. However, crystal structures show that the Cys-X-X-Cys of the second domain is buried and that the adjacent proline is in the trans, rather than the cis isomer. In these configurations, neither the “active-site” disulfide nor the backbone carbonyl preceding the proline is available to interact with substrate. The Eps1p structures thus expand the documented diversity of the PDI oxidoreductase family and demonstrate that conserved sequence motifs in common folds do not guarantee structural or functional conservation.     

Mechanisms of Zero-Lag Synchronization in Cortical Motifs

Zero-lag synchronization between distant cortical areas has been observed in a diversity of experimental data sets and between many different regions of the brain. Several computational mechanisms have been proposed to account for such isochronous synchronization in the presence of long conduction delays: Of these, the phenomenon of “dynamical relaying” – a mechanism that relies on a specific network motif – has proven to be the most robust with respect to parameter mismatch and system noise. Surprisingly, despite a contrary belief in the community, the common driving motif is an unreliable means of establishing zero-lag synchrony. Although dynamical relaying has been validated in empirical and computational studies, the deeper dynamical mechanisms and comparison to dynamics on other motifs is lacking. By systematically comparing synchronization on a variety of small motifs, we establish that the presence of a single reciprocally connected pair – a “resonance pair” – plays a crucial role in disambiguating those motifs that foster zero-lag synchrony in the presence of conduction delays (such as dynamical relaying) from those that do not (such as the common driving triad). Remarkably, minor structural changes to the common driving motif that incorporate a reciprocal pair recover robust zero-lag synchrony. The findings are observed in computational models of spiking neurons, populations of spiking neurons and neural mass models, and arise whether the oscillatory systems are periodic, chaotic, noise-free or driven by stochastic inputs. The influence of the resonance pair is also robust to parameter mismatch and asymmetrical time delays amongst the elements of the motif. We call this manner of facilitating zero-lag synchrony resonance-induced synchronization, outline the conditions for its occurrence, and propose that it may be a general mechanism to promote zero-lag synchrony in the brain.