- Research article
- Open Access
- Published:

# Avoiding transcription factor competition at promoter level increases the chances of obtaining oscillation

*BMC Systems Biology***volume 4**, Article number: 66 (2010)

## Abstract

### Background

The ultimate goal of synthetic biology is the conception and construction of genetic circuits that are reliable with respect to their designed function (e.g. oscillators, switches). This task remains still to be attained due to the inherent synergy of the biological building blocks and to an insufficient feedback between experiments and mathematical models. Nevertheless, the progress in these directions has been substantial.

### Results

It has been emphasized in the literature that the architecture of a genetic oscillator must include positive (activating) and negative (inhibiting) genetic interactions in order to yield robust oscillations. Our results point out that the oscillatory capacity is not only affected by the interaction polarity but by how it is implemented at promoter level. For a chosen oscillator architecture, we show by means of numerical simulations that the existence or lack of competition between activator and inhibitor at promoter level affects the probability of producing oscillations and also leaves characteristic fingerprints on the associated period/amplitude features.

### Conclusions

In comparison with non-competitive binding at promoters, competition drastically reduces the region of the parameters space characterized by oscillatory solutions. Moreover, while competition leads to pulse-like oscillations with long-tail distribution in period and amplitude for various parameters or noisy conditions, the non-competitive scenario shows a characteristic frequency and confined amplitude values. Our study also situates the competition mechanism in the context of existing genetic oscillators, with emphasis on the Atkinson oscillator.

## Background

In the relatively young field of synthetic biology [1, 2], there is increasing interest in the conception and construction of genetic circuits that are reliable with respect to their designed function. Having given the first step with the implementation of biological switches [3], the next step for synthetic biology was the construction of biological oscillators. The first successful implementation [4] constituted the onset of the quest for oscillators of tunable amplitude and/or period. Oscillators are important in biology for many reasons, since they are involved in the cell cycle, cell motion, embryonic development [5]. In some cases, the genetic machinery associated with the oscillatory behavior is rather small.

In the context of genetic designs, two- and three-element networks have been shown to be enough for implementing oscillatory behavior either of transient or sustained type [4, 6]. Several theoretical and experimental implementations of suggested robust tunable oscillators exist [7, 8]. However, predictable robustness of genetic networks is still a difficult task. The unpredictability of synthetic designs originates from the more-than-additive effect of assembling multiple building blocks, and has raised philosophical conundrums [9]. Intuition alone cannot grasp the effects of multiple regulatory interactions and thus mathematical models are particularly well suited for unraveling the implications of the underlying nonlinear interactions [10–14].

The way towards understanding and reliably predicting the information encoded by combinations of cis-regulatory sites lies in the coupling of such mathematical models with the experimental synthesis of regulatory promoter libraries [15–17]. Recent studies address the impact of some particular designs on the expected patterns of genetic switches and clocks [18, 19]. The authors of these studies have shown that the interplay between rate constants and circuit structure leads to nontrivial outcomes. Moreover, these studies suggest different views of how to explain the observed behavior. But several aspects have not yet been properly addressed, in particular the importance of the architecture of promoters and the resulting behavior.

In this paper, we show how the dynamics of a simple two-component activator-inhibitor oscillator is drastically affected by the architecture of the promoters and the nature of binding by the transcription factors. We show that competitive DNA-binding of several transcription factors acting at the level of a promoter leads to a more fragile oscillator than with non-competitive binding. The type of binding mechanism and the circuit architecture together determine the features of the resultant oscillations. The intention of the present work is to emphasize the importance of choosing the appropriate mathematical modeling in ensuring the experimental achievement of the desired biological function.

We also argue that detailed studies of these DNA-binding mechanisms could provide new insights into old experimental and theoretical studies of genetic circuits. A clear example is the experimental implementation of the Atkinson oscillator [6] for which several distinct mathematical models (and thus distinct dynamical behaviors) exist in the literature. Which among these models would be most suitable with respect to Atkinson's experimental implementation still remains to be determined, as the details of the biological binding mechanism are not sufficiently understood.

This paper is organized as follows: *Background* section continues with an extended introductory section on transcriptional regulation and genetic circuits that provides the framework of discussion of the existing models and experimental implementations. *Results and discussion* section details the main results on the competitive and non-competitive transcriptional binding, especially the mathematical characteristics of the route towards sustained oscillations. In *Conclusions* section we emphasize the implications of the DNA-binding mechanism on the reliability of the genetic circuit in terms of its desired function.

### Genetic circuits

Genetic circuits involve a certain number of genes that regulate one another's expression [20]. This regulation is achieved by means of the RNA or protein they encode for, that serve as regulators for other genes' expression. As major regulators, the protein factors influence the rate of transcription, translation or post-translational modification. RNAs, too, can have a variety of regulatory functions (i.e translation inhibition by microRNAs). The regulatory mechanisms are still in the process of being elucidated [21] and have been already addressed in synthetic biology [22, 23].

Among the reasons for constructing synthetic biological circuits, there is the goal of understanding the fundamental building blocks and regulation mechanisms of biology [24, 25], and subsequently manipulating and monitoring biological processes at the DNA level [26, 27]. The majority of the existing works in systems and synthetic biology have focused on transcriptional regulation, as the interactions are relatively modular and can be altered at will. Additional File 1 constitutes an appendix on the basics of this regulation mechanism and the definition of the Hill function, which describes transcriptional interactions. Comparatively few studies have focused on other types of regulation, such as the post-transcriptional regulation employed in Fussenegger's Lab [8, 28]. A combination of transcriptional regulation and post-translational modification is common in the literature, where the latter refers to a nonlinear degradation of a transcription factor [18, 29]. In general, this nonlinear degradation involves a saturated response [30–32], sometimes referred to as enzymatic control. The term of enzymatic control was inherited from the first studies on chemical circuits [33-36] which could be considered as the *in silico* progenitors of today's genetic circuits. Moreover, the coupling of transcriptional regulation with metabolic flux [37] has opened yet another door to innovation and control of biological circuits. Even though much work is needed in terms of reliability and control of any individual oscillator design, there are already advances towards the next step: coupling of genetic oscillators [38, 39]. This step allows accessibility to new behaviors that are not possible at the single-oscillator level, similar to the introduction of spatial features into homogeneous systems [40]. Given all these distinct mechanisms and types of control, one should intuitively expect that not only the "*sign*" of genetic interaction (positive or negative), but also the *type* of interaction (transcriptional, post-transcriptional, post-translational) will affect the behavior of the circuit. These dynamic consequences have only begun to be explored and are far from being clarified.

When referring to an interaction network, one considers two types of interactions or loosely called feedbacks: positive (protein A promotes the production of protein B) and negative (A inhibits the production of B). Beyond the details of the interaction itself, the theory of circuits associates an overall sign to the interaction network itself, with the circuit being positive or negative depending only on the parity of the number of negative interactions in the network. A network is positive if it contains either positive interactions alone, or an even number of negative ones, while a negative network has an odd number of negative feedbacks [49, 50]. The general theory of circuits has demonstrated that a positive network is a necessary (though not sufficient) condition for multistability (and thus differentiation), while a negative one, for homeostasis (see Ref. [10] for a review). More precisely, it was proved [51] that there is a qualitative difference in the behavior of circuits containing even or odd negative feedbacks, respectively, with the latter showing the possibility of sustained oscillations. A negative feedback alone produces a stable state, unless there exists a sufficiently long time delay, in which case oscillations are possible. In the case of a genetic circuit with a positive feedback, sustained oscillations are not possible even with time delay (not proved though, but reasonably argued [52]).

The situation might be different in non-genetic feedback loops, for example enzymatic circuits. This is the case commented by Smolen et al. [53] regarding the enzyme phosphofructokinase that shows oscillations in a non-delay (non-genetic) positive loop with limiting rate-supply of precursors. As mentioned above, we are mainly interested in studying the transcriptional interaction, thus we shall not discuss the known circuits that combine transcriptional control with enzymatic one [29, 30, 32, 54], nor shall we comment for the moment on the mRNA-level (post-transcriptional) regulation circuits [8]. Thus, the models discussed in the current study employ only transcriptional regulation. The circuits studied by previous works are schematically represented in Figure 1 and discussed in Table 1.

#### Oscillations from negative feedbacks alone

The simplest genetic circuit capable of producing sustained oscillations is a circuit of one gene with delayed self-inhibiting regulation (Figure 1A). The theoretical analysis of such systems has been performed long before the experimental advances contemplated the measurements and implications of transcriptional and transport delays [55]. More recently, stochastic models [41, 56] have added to these deterministic ones, showing that noise promotes the oscillatory behavior in a delayed negative circuit. Due to its simplicity, the case of self-inhibitory feedback with delay has been addressed in synthetic biology as well [42]. Beyond synthetic circuits, transcriptional and transport delay are expected to be more relevant in eukaryotes than in prokaryotes, and thus can and have been assigned as the causal factor of oscillatory behavior in only a few specific mammalian genetic networks, all composed of short negative feedback loops [57]. Nevertheless, the relation between delay, protein half-life and cooperativity required for oscillatory solutions (see page 192 of Ref. [5] for the analytical relation) is not easily fulfilled under normal conditions, and this is also the reason for not considering these delays in the general-purpose modeling of genetic circuits.

As discussed above, a circuit of odd-number negative feedbacks is equivalent to a positive-and-negative feedback circuit, and thus oscillations are theoretically possible. For example, a three-gene synthetic oscillator based on (non-explicitly delayed) negative interactions alone has been already verified experimentally (the "repressilator" [4]; Figure 1H) and its generalization has been studied theoretically [48]. The clarification of the specific parameter ranges for which oscillations occur is crucial for the design of a synthetic circuit, such as the repressilator [4]. However, it is an issue rarely considered. And the subsequent experimental challenge is to achieve these parameter ranges in the experimental circuit. For this purpose, the repressilator implementation required extreme solutions [47]: very strong promoters, very high degradation rates of the proteins and multiple binding sites to increase the nonlinearity of the interactions (that is, increase the Hill coefficient - see Additional File 1).

## Results and Discussion

### On two-component genetic oscillators

To achieve the goal of reliable synthetic circuits, several strategies have been employed by the community. One strategy is the *in silico* evolution of biochemical and genetic networks based on simplified principles of biological interactions and where the fitness function rewards circuits characterized by the function required to be implemented [58, 59]. Another strategy is to bring together experimental data and modeling framework in search of designs that are tunable and robust to fluctuations [16, 60, 61]. Yet another strategy of a more systems-biology flavor points towards understanding the dynamics of minimal fundamental modules [18, 32, 43], and the current work belongs to this approach.

Sufficient conditions for the occurrence of sustained oscillations are still an open issue in nonlinear dynamics theory, and naturally they are continuing to be delineated in feedback-based gene networks [10]. It is clear though that the combination of both negative and positive feedback is a necessary (but not sufficient) condition in order to obtain oscillations for (non-explicitly delayed) two-gene circuits [43]. Analytically, the simplest genetic oscillator of two genes needs two feedbacks: a positive and a negative one (see Figure 1B). However, this might not be the easiest system from the experimental point of view, due to the high value of cooperativity (Hill exponent *n* ≥ 3) required for sustained oscillations (i.e. non-damped). For a multimer binding to a single binding site, the Hill coefficient coincides with the multimer's degree (i.e. *n* = 2 if the TF is a dimer; *n* = 4, if it is a tetramer). From the experimental point of view, one hopes to obtain a high Hill exponent by introducing several binding sites for the transcription factor (TF). In this way and depending on the binding affinities of the TF, binding cooperativity *may* occur, leading to a higher exponent than the one corresponding to the degree of the multimer employed. For example, a Hill coefficient of *n* = 3 as in the example above may be obtained with dimers binding to two sites if some degree of cooperativity exists [62, 63].

Figure 1C, D, E refer to a genetic *relaxation-based* oscillator containing three interactions: two positive and one negative. When the expression of a TF, say protein **A**, is controlled by two TFs, say **A** and **R**, through a positive and a negative feedback, respectively, one can imagine several transcriptional scenarios. In Figure 2 we illustrate several implementations compatible with the schematic representation of the circuit shown in (*A*). For the competition case shown in (*B*), the activator **A** (circle) and the repressor **R** (square) compete for the binding on the promoter in the activator's module. In this case, the transcription of the activator occurs from two possible states of the activator promoter: occupied with **A** or occupied with **R**. For the design in (*C*), the inhibition occurs through post-translational control: the degradation of the activator is catalyzed by the repressor. These two designs have been studied by Guantes&Poyatos (2006) [18] and we follow here their terminology: Design I and Design II, respectively.

We introduce in Figure 2D Design III which confers to both TFs their binding sites. This implies that the transcription of **A** results from several combinations of the A promoter: *AX*, *XB*, *AB*, with *X* implying free site. Compared to Design I, the mathematical form of Design III is closer to the model associated to the Atkinson oscillator introduced by Scott et al. [44]. Guantes&Poyatos [18] claim that Design I corresponds to Atkinson model. And indeed Atkinson et al. [6] envision that, even though the experimental implementation is similar to Design III from Figure 2, the resultant activation DNA loop and the repression DNA loop would be mutually exclusive [64]. This exclusion might thus imply a competition process, leaning towards effective similarity to Design I. From the experimental perspective, the data has revealed *damped* oscillations rather than the sustained ones expected from the model [6]. There is still no clear identification of the cause of damped oscillations, and no definite study on the extent of the exclusion between the two DNA loops.

Below we shall assess the implications of binding competition by comparing Design I with Design III. Before detailing these results, we wish to establish a framework for their comparison in the light of two existing works: Guantes&Poyatos [18] and Conrad et al. [19]. Guantes&Poyatos [18] argue that there is a fundamental difference, in terms of dynamic behavior, between Design I and Design II. The number of steady states in the two designs and their stability properties completely determine the onset of *sustained* oscillations. The way these equilibria lose their stability, that is the underlying bifurcation type, determines the type of oscillators. More precisely, when a parameter of the system is varied, a *bifurcation* is said to occur at a certain parameter value (bifurcation point) for which a qualitative change in the behavior of the system appears, such as the birth of sustained oscillations [65, 66].

We suggest that circuit architecture AND DNA-binding processes determine the type of bifurcations leading to oscillatory solutions, types that confer specific features to the oscillations. Previous results show that Design I induces oscillations with large activator amplitudes and arbitrarily small frequencies, and acts as an "integrator" of external stimuli, while Design II shows emergence of oscillations with finite, and less variable, frequencies and smaller amplitudes, and detects better frequency-encoded signals ("resonator"). This classification originates from neuronal oscillations where the type of neuron (type I or type II) determines entrainment properties, phase response curves and robustness to noise [67]. In a nutshell, Guantes&Poyatos [18] suggest that the underlying biochemical mechanism of repression may determine whether the genetic circuit behaves as a resonator or as an integrator. This statement is essentially correct, as the mathematical difference between the two models makes the dynamical behavior possible. On the other hand, Conrad et al. [19] argue that the rate constants rather than the biochemical mechanism have the capacity to determine which type of behavior is observed. They show that for a certain parameter interval, Design II passes initially through an integrator-originating bifurcation and subsequently (higher parameter values) through a resonator-type bifurcation. Even though the parameters are responsible for the occurrence of these bifurcations, it is nevertheless the underlying mechanism (how the architecture is biologically implemented) that makes them possible. In our comparative study, we show that indeed by changing different parameters, we allow different types of oscillations-originating routes, but we also show that the biological regulation mechanisms make them possible.

#### To compete or not to compete?

In terms of equations, the systems describing the implementations from Figure 2, following the notations from Guantes&Poyatos [18], are:

With ≡ *dx*/*dτ*, and the rest of the variables and parameters being defined in Table 2. As the details on how Design I and II were obtained and the parameters' definition are extensively presented in the Supplementary File of [18], we have chosen to include in the Additional File 1 only the necessary details for recovering the equations associated to Design III. Nevertheless, we comment briefly on the meaning of the parameters. All parameters are non-dimensional and defined identically in all three models, allowing a direct (and correct) comparison between the results. The factor *α* gives a measure of how much stronger the transcription from the activator-bound promoter is, compared to the basal transcription from the free one. The parameters *β* and *γ* contain the total non-dimensional strength of the transcription and translation for the activator and repressor, respectively. The *σ* parameter contains the DNA-binding properties of the repressor multimer scaled to the activator's. For Design II, the parameter *σ*' differs in expression from that of *σ* in the other models, justifying the different notation (see [18]). For the current study and in order to restrict to fewer parameters, we shall consider *σ* = 1, as [18]. A parameter that plays a crucial role in the existence of oscillations in the models defined above is the quotient of degradation rates Δ = *δ*_{
A
}/*δ*_{
R
}, with *δ*_{
A
}, the activator's degradation rate. By scaling the time to the repressor's degradation rate, *δ*_{
R
}, the period of the oscillations resulting from these models is expressed here in units of *δ*_{
R
}, a detail that needs to be remembered when comparing the resulting periods of the different designs.

In addition to the above discussion of the parameters and their meaning, let us comment briefly on the similarities and differences on the equations describing the three designs. All designs share the same rate law for the repressor, while the definitory features concern the rate law for the activator. All rate-law terms have a form of the type: production term minus degradation term. Even from the equations, one can notice that Design II is different from the other two by its particular degradation term, the post-transcriptional regulation. On the other hand, the production terms for Designs I and III differ in a subtle manner: while the transcription factors appear summed in the denominator of Design I, they appear multiplied for Design III. The *ad hoc* general rule cited in the literature for the case of more proteins affecting a gene is that the functions appear multiplied [59]. By our reasoning, the multiplication is appropriate only in cases of non-competitive transcriptional regulation. Even though cases of explicit competitive regulation are mentioned in the literature [10], the current study is the first to inquire on the differences in the dynamical behavior between the two scenarios.

We developed numerical algorithms to simulate the equations characterizing Design I and Design III from Figure 2 and these codes are included in Additional File 2. From the resultant time series, we determined the oscillation period (if any) and represented it in Figure 3. The figures plot the value of the period for different values of *β* (abscissa) and *γ* (ordinate), where the period value is represented by a color from the color code in the right-hand-side of the figure. The dark-blue color is associated to cases lacking *sustained* oscillations.

We have chosen fixed values for Δ, *α* and *σ*, and varied *β* and *γ*, in order to compare to existing results [18]. By this study, it results that Design III is characterized by a more extended oscillatory region in the parameter space than Design I. Moreover, for Design III, the oscillatory region continues beyond the parameter limits from Figure 3B (see Additional File 3). For completion, we illustrate in Additional File 4 the change in the position of the oscillatory region when Δ (the quotient of degradation rates) is increased. For the same purpose, Additional File 5 contains the oscillatory region for Design II from [18], directly comparable with Figure 3.

By fixing Δ and allowing *β* and *γ* to vary, we may relate to the experimental context. More precisely, the binding constants and could in principle be experimentally-tunable as in the implementation of the oscillator from Stricker et al. (2008) [7], allowing *β* and *γ* to take different values within a certain interval. These parameters could also be varied by varying the total promoter concentration, and and in the expression of *β* and *γ*, respectively (Table 2). The degradation rates of the proteins could also be varied to a certain measure, and thus Δ parameter could be used as control parameter too.

In addition to the existence of the sustained oscillations, we also identified the bifurcation types that allowed this behavior. We illustrate in Figure 4 how the sustained oscillations originate when crossing the parameters space in Figures 3A and 3B from left to right (from lower to higher *β*) for two fixed values of *γ*: *γ* = 0:01 (panels A and C) and *γ* = 0:05 (panels B and D). As already identified by the studies mentioned above, two main bifurcation types are observed: the Hopf bifurcation and the saddle-node on an invariant cycle, or SNIC, bifurcation [65]. Below, we shall comment on the defining features of these two distinct types of bifurcations, without entering in excessive technical details.

As with many other bifurcation types, the Hopf bifurcation comes in two flavors: the *supercritical* (or soft) and the *subcritical* (or hard) (defined below). In both cases, the real part of the complex eigenvalues becomes positive, and the unique stable state of the system loses stability. For the supercritical type (referred to as Hopf(p) in the figures), the oscillations begin with low amplitude near the bifurcation, which increases with the change in the parameter. It is the case of Figure 4B. Even though it has the same eigenvalue behavior, the subcritical case (referred to as Hopf(b) in the figures) is more dramatic: the *sustained* oscillations start or end with high amplitude at the bifurcation point. It is a consequence of two limit cycles, a stable and an unstable one, colliding at the bifurcation point, leaving an unstable fixed point surrounded by a stable limit cycle (Figure 4C; Additional File 6).

We wish to especially stress the fact that systems undergoing Hopf bifurcations, both subcritical and supercritical, always exhibit *damped* oscillations prior to the bifurcation point (where sustained oscillations begin), whereas systems near saddle-node bifurcations do not (page 13. of Ref. [67]). One would observe the following: at a value of the control parameter that leads to negative real eigenvalues, the system stabilizes to a fixed value, without any oscillations; by changing the control parameter until the eigenvalues become complex but still with negative real part, the system would still stabilize to a fixed value, but through damped oscillatory trajectories; finally, a bifurcation is said to occur when, by further changing the control parameter, the system does not stabilize to a fixed value, but continues to oscillate indefinitely, and thus a new qualitative behavior appears. Additional File 7 illustrates an example of how sustained oscillations exist for a certain value of the control parameter Δ after the Hopf bifurcation (Δ = 4.5), and how these are preceded (that is, for values of the control parameter prior to the birth of the stable periodic orbit, Δ = 4.3) by damped oscillations.

In Figure 4, the gray areas denote ranges of *β* for which damped oscillations exist. In these ranges, the eigenvalues are complex, but with negative real part. The existence of an imaginary part of the eigenvalues leads to the oscillatory behavior, but due to the real part being negative, these oscillations are damped.

For example, damped oscillations exist for *β* > *β*_{
c
}= 0.225 in Figure 4C, with a higher damping rate the farther away from the bifurcation point.

In the SNIC bifurcation, a single real eigenvalue changes sign; geometrically, a stable node and a saddle point meet, annihilate each other, and leave in their wake a limit cycle. When crossing a SNIC bifurcation, the oscillations start with very long period, and tend to stabilize to a constant value as the parameter is changed further (see also [18]). This behavior is contrasted by the Hopf bifurcation for which the oscillations have always a characteristic period. With these comments in mind, the types of bifurcations can be easily identified in Figure 3 through the color code. For Design III, the frontier around *β* ≈ 0.1 is characterized by a SNIC bifurcation as *β* is increased, as the region shows long periods (Figure 4C). Crossing the oscillatory region for higher *γ*, the period appears to be constant, as both entering and leaving the oscillatory band occur through Hopf bifurcations (Figure 4D). On the other hand, for Design I, the exit from the oscillatory band occurs always through a SNIC bifurcation.

As discussed by Conrad et al. [19], the parameters determine the type of bifurcations the system experiences, and we illustrate this remark by varying parameters other than *β* and *γ*. In Figure 5 we show the types of bifurcations occurring when Δ is varied for fixed *β*, *γ* . This analysis directly compares with the previous studies [18, 19]. As commented above, the spiking-like oscillations are a consequence of SNIC bifurcation, while Hopf one produces more sinusoidal oscillations. An example of these two types of oscillations is illustrated in Figure 6.

Returning to Figure 4, notice also how the behavior of Design III from panel *C* matches the behavior of Design II from [19](their Fig. 3b): depending on the parameter to be varied, a design can present both types of bifurcations. Guantes&Poyatos [18] stated that Design I acts as a signal integrator (due to SNIC bifurcation) while Design II resonates with specific frequencies of the stimulus as a resonator (due to the Hopf bifurcation). Conrad et al. [19], on the other hand, highlight the fact that Design II can suffer both types of bifurcations depending on the parameters. This is the reason for using two distinct sets of parameters in Figure 4 and Figure 5, emphasizing their statement. According to Figure 4C, Design III behaves very similarly to Design II (see study of Design II in Ref. [18]) in spite of very different underlying biochemical mechanisms that produce the negative feedback. Thus we would suggest an intermediate position that states that indeed the parameters determine the type of bifurcations that occur, but only when the *biochemical mechanism* (or the resulting mathematical formulation) allows these bifurcations to occur at all. For example, Design I cannot pass from an integrator to a resonator by changing one parameter as happens for Design II and Design III.

As a final remark and returning to the discussion on the Atkinson oscillator, it is interesting to notice that the experimental implementation showed damped oscillations. We have mentioned that damped oscillations occur prior to a Hopf bifurcation, and not to SNIC bifurcation. Judging from the discussion of the two designs (Figure 4), it is reasonable to infer that the experimental implementation from [6] might be closer to Design III than to Design I. In other words, the activation and repression DNA loops might not be completely exclusive, as Atkinson et al. claim. This would be in accordance to the model studied by [44]. Nevertheless, the clarification of this issue needs further experimental studies.

#### Dynamical considerations

Apart from the computation of the parameters leading to oscillations, let us now present a dynamical perspective on the differences between Design I and III. For a better understanding of the oscillatory features of the systems, we have studied in more detail the stability properties of the fixed points associated to the three designs. More specifically, we have calculated the eigenvalues for the fixed points associated to the three designs for various parameters spaces. The Additional File 8 includes the results of this study of which we mention here only the main conclusions. We present the stability properties of the fixed point of the system for the parameters space (γ/β, σ) for fixed Δ (Figure S8.1), (γ/β, Δ) for fixed σ (Figure S8.2) and also, for comparison with Figure 3, the space (γ, β) for fixed Δ and σ (Figure S8.4). A first conclusion is that Design III shows oscillations, *both* damped and sustained, for a larger fraction of the parameter space (Figure S8.1 and S8.2), further reinforcing Figure 3. From this perspective, one could say that Design III is substantially more robust than Design I to changes in the parameters' values. Secondly, while Design III appears to be an oscillator of a specific frequency for a wide parameters range (Figure S8.5), Design I presents a wide distribution of possible oscillation periods depending on the parameters' values, due to the SNIC bifurcation. In this aspect, Design II comes as an intermediate step between the other two designs.

It is interesting to mention that, while the repressor's amplitude for the three models is similar, the distribution of activator's amplitude is specific for each design. It ranges from the queue of high amplitudes for the spiky oscillations of Design I, to an almost uniform and narrow distribution for Design III (Figure S8.5). The difference in this aspect between Design I and III is somehow expected, as the amplitude of the latter is further restrained by the term *σx*^{2}*y*^{2} at the denominator. Both Design II and III show an interesting cut-off in activator's amplitude limiting it at low values. From these features can be seen that the biochemical mechanism leaves a characteristic fingerprint on the nature of the resultant oscillations. Having observed these features in amplitude and period, we have inquired on the dynamical causes responsible for these limitations. As already illustrated by [18] in their comparison between Design I and II, the answer to this question could be found by studying the nullclines associated to the systems. The nullclines denoted as *y*_{
x
}and *y*_{
y
}represent the function *y*(*x*) obtained from the conditions = 0 and = 0, respectively. The fixed points - stable or unstable - of the system are found at the intersection of *y*_{
x
}and *y*_{
y
}, and the form of these functions also gives clues on the general dynamics of the system beyond the steady states. The definition of the nullclines and the extended results of this analysis are included in Additional File 9. In spite of very different function for the three systems, it is surprising to find that *y*_{
x
}for Design II and III are tightly related, with the first being the square of the other. This "similarity" however does not imply necessarily a similarity in dynamic behavior. The analysis of the nullclines as detailed in Additional File 9 shows that a change in *β* does not imply a change in the shape of the nullclines, as the shape is mainly controlled by *α*. For this reason, these two designs present a characteristic frequency and a limitation in *x* amplitude. Moreover, from the equations can be seen that, starting from an oscillation-producing parameter case, an increase in *β* implies upwards displacement of *y*_{
x ,
}and an increase in Δ or *γ* implies upwards displacement of *y*_{
y
}. For this reason, by increasing both *β* and *γ* simultaneously, their crossing is maintained and thus the oscillations too. This does not occur for Design I for which a change in *β* implies not only a displacement, but a change in shape, leading to a three-crossing case: no oscillations. Besides explaining the loss of oscillation capability for this design beyond a certain threshold value in *β*, the shape of the nullclines, as discussed in the additional file, is responsible for the high amplitudes in both *x* and *y*, as can be see in Figure S8.5.

#### Influence of internal noise

There has been increasing acknowledgment of the fact that internal noise at cellular level - low number of molecules [68–70] - can in principle change qualitatively the behavior of the system with respect to the one predicted from noise-free (deterministic) modeling: noise-induced stabilization [71], noise-induced oscillations [29]. Moreover, there is clear evidence that the transcription and translation level are indeed noise-prone processes [72]. For these reasons, we consider that the influence of internal noise in the three systems studied here constitutes a work on its own. Nevertheless, we have chosen to mention here how noise affects the behavior of the three oscillators close to the bifurcation regime. It is observed that, for values of the parameters leading to stable states in the deterministic cases close to the bifurcation, the noise induces oscillations in all three designs, but with different characteristics (Additional File 10). While Design I and II show period distribution typical for noise-induced oscillations - long tail distribution of periods-, Design III still maintains its characteristic period. These results as well as the deterministic ones from Figure S8.5 (amplitude distribution) show that while Design I has high amplitudes, it is unreliable at period-level (long tail when noise exists), Design III has low amplitude, but it is reliable in terms of period constancy. We remind the reader that Design III is characterized also by a large fraction of parameters space leading to oscillations. All these reasons considered, Design III appears as a more reliable oscillator.

#### Other two-component oscillators

Once more, we wish to emphasize an important ingredient that makes oscillations possible in the above models: the time scale difference between the degradation rates of the two TFs given by Δ (the quotient of degradation rates). We have seen that for all these models, no oscillations are possible unless a substantial difference between the two degradation rates exists (Δ ≫ 1). The requirement thus says that *δ*_{
A
}≫*δ*_{
R
}. As commented also by Atkinson et al. [6], strategies for generating oscillations include increasing *δ*_{
A
}and/or decreasing *δ*_{
R
}. However, decreasing *δ*_{
R
}is not a good strategy as we have seen that the resulting oscillation period depends on *δ*_{
R
}. The shorter *δ*_{
R
}, the longer the period. For example, a value of *δ*_{
R
}= 0.02 *h*^{-1} typical value of *δ*_{
A
}= 1 *h*^{-1} would result in a good Δ for the occurrence of oscillations, but in a period of about 100 hours for Design III. Such a long period, much longer than the cell-cycle period, could be altered by the intrinsic cellular division period. On the other hand, Elowitz&Leibler [4] follow the strategy of increasing *δ*_{
A
}, that is synthetically design proteins that have a short life time (of the order of minutes). This strategy appears to be widespread nowadays.

Another strategy to facilitate oscillatory conditions is a more pronounced nonlinearity of the regulatory function, that is consider a Hill function of higher exponent. The mathematical model for Design III considered dimers as active TFs and one binding site for each feedback. For example, let us assume the general case for Design III. By the term *general* we refer to a non-competitive design similar to Design III, but having several bindings sites for activators and repressors. In this case thus, an approximate model would employ the general form of Hill function:

However, keep in mind that the parameters *β* and *γ* depend on the degree of multimerization of the TFs. If *n* ≠ 2 and *m* ≠ 2, the parameters *β* and *γ* differ in expression from those employed in the previous section, and thus a direct comparison with the previous results is difficult. Nevertheless, the interesting remark to make here is that oscillations are possible even for Δ = 1 (no difference in degradation rates) when *n* = 3, as seen in Figure 7. The figures illustrate the general behavior of the model to changes in the multimerization degree of the repressor. It is apparent from these figures that *n*, the nonlinearity of the positive feedback, is responsible for the *existence* of the oscillations, while *m* relates to the extension of the parameter space leading to these oscillations: the higher the *m* (degree of repressor's multimerization), the wider the oscillatory region.

Again concerning the nonlinearity of the feedback, let us return to the experimental implementation of the Atkinson oscillator. The experimental design [6] appears in Figure 8 (and also schematically in Figure 1C), including two binding sites for each TF. The full-competition scenario occurs under the hypothesis that when both sites of the activator (or of the repressor) are occupied by the respective dimers, a DNA loop is formed, blocking the formation of the repressor's loop (or of the activator's loop, respectively). The mathematical model describing this scenario would be based on Design I but characterized by *n* = *m* = 4:

and where Δ, *β* and *γ* have the same expression as in Design I and thus are directly comparable. Notice in Figure 9 the reduced sustained-oscillations region for the model in eqs. (3-4) compared to Figure 3A. It is interesting that, in this case, the increase in the nonlinearity of the feedbacks leads to a decrease of the oscillatory region.

The model from eqs. (1-2) is a simplified version of the Atkinson model as described by Scott et al. [44]. There, the model included the approximation of the Hill functions as above, but considered explicitly the dynamics of the messenger RNAs, resulting in a system of 4 variables instead of 2 variables as employed here. The approximation that allowed us to consider only the protein variables, as in the above equations, is based on the fact that translation is much faster than transcription, and thus as soon as the mRNA is formed, one can consider that the subsequent protein is formed too. It is a drastic approximation and it has already been proved to have relevant consequences: explicitly modeling the mRNAs yields a better fit to the experimental data than the simplified model [7]. Nevertheless, the simplified model is a more tractable tool for studying the general characteristics of the system.

Still in two-gene systems, another strategy to facilitate oscillations is introducing another positive feedback for Design I, more precisely an autocatalysis on the repressor. This eliminates the necessity of large Δ, while maintaining dimers and one binding site. Of course, this leads to another architecture than the one shown in Figure 2A. This architecture is contemplated by [46] (see Table 1 and Figure 1G):

Notice in Figure 10 how the oscillatory region considerably changes in this architecture. The parameters are directly comparable with the Designs I, II and III. Let us notice that the activator must be a multimer (at least a dimer) for the occurrence of oscillations, but also a significant value of Δ is also needed, as Δ = 1 is not sufficient for oscillations to occur. As already discussed for the Atkinson model above, the degree of multimerization of the repressor is not as critical as that of the activator. Also Smolen et al. [46] introduce the above architecture considering the repressor as monomer, and emphasizing the necessity of a dimer activator for the oscillations to occur. Once again, the nonlinearity of the transcriptional response is the driving factor in the instability required for the oscillations.

## Conclusions

In the present work, we have analyzed two DNA-binding mechanisms in composed promoter sites and how they affect the existence of oscillations in genetic circuits. There is general consensus that both positive and negative genetic interactions are needed to obtain a robust oscillator [73]. Given a number of such interactions, there exist several mathematical and experimental proposals for oscillators. Even though there is sufficient proof that two positive interactions and a negative one would produce a genetic oscillator, there are various biochemical mechanisms compatible with such an architecture (Figure 2). The differences between these mechanisms and their consequences at the level of robustness and reliability of the resultant oscillator are a new direction of study [18, 47].

From our results, we have emphasized the fact that, when these interactions, either positive or negative, are accomplished through transcriptional regulation, the interactions between various binding sites on the promoter need to be taken into consideration. We showed the non-competition scenario to present robustness to parameter changes defined by the existence of oscillatory solutions for a much wider parameters range than for the competition one. Moreover, the two scenarios imply different routes toward the occurrence of oscillations, that is different bifurcation types. While the competition model is mainly associated to the SNIC bifurcation, the non-competition one relies more on the Hopf bifurcation. While the former bifurcation type leads to more spiking-like oscillations, the latter shows more sinusoidal ones. In addition, a characteristic of the Hopf bifurcation is that sustained oscillations are preceded by damping oscillations. This may prove to be an important feature to consider when interpreting experimental data. Many other interesting properties arise as a consequence of the type of the oscillation-generating bifurcation, an issue that has been extensively analyzed in previous works [18, 19, 67] and thus has not been discussed here. However, consistent with previous works, we found that which parameter is changed - and in which range it is changed - determines the bifurcations that occur. Nevertheless, not only the parameters are responsible for the resultant dynamics of the system, as argued by Conrad et al. [19], but most important of all, it is the underlying logic of the biochemical mechanism that makes the dynamics possible at all. In addition, we conclude from our results that also very different designs (or biochemical mechanisms) can produce very similar behaviors. The non-competition scenario analyzed here produces an overall dynamic behavior similar to the post-transcriptional control scenario described by Guantes&Poyatos [18].

The Atkinson oscillator introduced and experimentally implemented by [6] has been modeled both through competition [18] and non-competition models [44]. As the experimental implementation failed to provide sustained oscillations, it is interesting to assess the implications of the two scenarios explored in the current study in the context of this experimental implementation. The goal of the current study was to pinpoint the consequences of the biochemical mechanisms and the mathematical approximations taken into consideration when constructing the associated theoretical model, and thus not to present a detailed study of the already-implemented experimental design [6]. Nevertheless, we consider that our work sheds new light on the behavior of this genetic oscillator.

Here we have used simplified models of genetic circuits. In this simplification, the modeling of genetic circuits includes only the temporal evolution of the protein concentration [45, 74]. This allows a straight-forward mathematical analysis of the system and its properties. Moreover, we have chosen here to employ the simplified approach as our goal was not a characterization of a specific experimental design, but to emphasize the importance of taking into consideration the underlying biochemical mechanisms and the promoter architecture at the moment of constructing a mathematical model of a genetic circuit. In summary, to have the largest chances of obtaining oscillations, synthetic biologists should build circuits of the form of Design III, avoiding competition between activator and repressor. Additionally, faster degradation rates for activators than repressors are preferred (or adding an autocatalysis for the repressor). Finally, achieving non-linear transcription responses (multimeric factors, at least dimers, and/or multiple binding sites) also favor oscillatory behaviour. Models that are more robust to parameter changes should be easier to construct experimentally. We therefore hope that this work will help bring closer the theoretical and experimental communities in the fields of systems and synthetic biology.

## Methods

For the numerical simulations, we developed C++ codes employing routines from the free-software GNU Scientific Library (Additional File 2). In addition, the bifurcation diagrams employ results obtained with the use of the AUTO software within XPPAUT program.

## References

- 1.
Haseloff J, Ajioka J: Synthetic biology: history, challenges and prospects. J R Soc Interface. 2009, 6: S389-S391. 10.1098/rsif.2009.0176.focus

- 2.
Benner SA, Sismour AM: Synthetic biology. Nat Rev Genet. 2005, 6 (7): 533-543. 10.1038/nrg1637

- 3.
Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403 (6767): 339-342. 10.1038/35002131

- 4.
Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403 (6767): 335-338. 10.1038/35002125

- 5.
Glass L, Mackey MC: From clocks to chaos The rythms of life. 1988, Princeton,

- 6.
Atkinson MR, Savageau MA, Myers JT, Ninfa AJ: Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli. Cell. 2003, 113 (5): 597-607. 10.1016/S0092-8674(03)00346-5

- 7.
Stricker J, Cookson S, Bennett MR, Mather WH, Tsimring LS, Hasty J: A fast, robust and tunable synthetic gene oscillator. Nature. 2008, 456: 516-519. 10.1038/nature07389

- 8.
Tigges M, Marquez-Lago TT, Stelling J, Fussenegger M: A tunable synthetic mammalian oscillator. Nature. 2009, 457 (7227): 309-312. 10.1038/nature07616

- 9.
Isalan M, Morrison M: This title is false. Nature. 2009, 458 (7241): 969- 10.1038/458969a

- 10.
Smolen P, Baxter DA, Byrne JH: Modeling transcriptional control in gene networks-methods, recent results, and future directions. Bull Math Biol. 2000, 62 (2): 247-292. 10.1006/bulm.1999.0155

- 11.
de Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol. 2002, 9: 67-103. 10.1089/10665270252833208

- 12.
Cohen JE: Mathematics is biology's next microscope, only better; biology is mathematics' next physics, only better. PLoS Biol. 2004, 2 (12): e439- 10.1371/journal.pbio.0020439

- 13.
Goldbeter A: Computational biology: a propagating wave of interest. Curr Biol. 2004, 14 (15): R601-602. 10.1016/j.cub.2004.07.036

- 14.
Mogilner A, Wollman R, Marshall WF: Quantitative modeling in cell biology: what is it good for?. Dev Cell. 2006, 11 (3): 279-287. 10.1016/j.devcel.2006.08.004

- 15.
Guido NJ, Wang X, Adalsteinsson D, McMillen D, Hasty J, Cantor CR, Elston TC, Collins JJ: A bottom-up approach to gene regulation. Nature. 2006, 439 (7078): 856-860. 10.1038/nature04473

- 16.
Ellis T, Wang X, Collins JJ: Diversity-based, model-guided construction of synthetic gene networks with predicted functions. Nat Biotechnol. 2009, 27: 465-471. 10.1038/nbt.1536

- 17.
Gertz J, Siggia ED, Cohen BA: Analysis of combinatorial cis-regulation in synthetic and genomic promoters. Nature. 2009, 457 (7226): 215-218. 10.1038/nature07521

- 18.
Guantes R, Poyatos JF: Dynamical principles of two-component genetic oscillators. PLoS Comput Biol. 2006, 2 (3): e30- 10.1371/journal.pcbi.0020030

- 19.
Conrad E, Mayo AE, Ninfa AJ, Forger DB: Rate constants rather than biochemical mechanism determine behaviour of genetic clocks. J R Soc Interface. 2008, 5 (Suppl 1): S9-15. 10.1098/rsif.2008.0046.focus

- 20.
Lewin B: Genes IX. 2007, Jones & Bartlett,

- 21.
Storz G, Altuvia S, Wassarman KM: An abundance of RNA regulators. Annu Rev Biochem. 2005, 74: 199-217. 10.1146/annurev.biochem.74.082803.133136

- 22.
Davidson EA, Ellington AD: Synthetic RNA circuits. Nat Chem Biol. 2007, 3: 23-28. 10.1038/nchembio846

- 23.
Win MN, Smolke CD: A modular and extensible RNA-based gene-regulatory platform for engineering cellular function. Proc Natl Acad Sci USA. 2007, 104 (36): 14283-14288. 10.1073/pnas.0703961104

- 24.
Hasty J, Isaacs F, Dolnik M, McMillen D, Collins JJ: Designer gene networks: Towards fundamental cellular control. Chaos. 2001, 11: 207-220. 10.1063/1.1345702

- 25.
Mukherji S, van Oudenaarden A: Synthetic biology: understanding biological design from synthetic circuits. Nat Rev Genet. 2009, 10 (12): 859-871.

- 26.
Boyle PM, Silver PA: Harnessing nature's toolbox: regulatory elements for synthetic biology. J R Soc Interface. 2009, 6: S535-S546. 10.1098/rsif.2008.0521.focus

- 27.
Bhalerao KD: Synthetic gene networks: the next wave in biotechnology?. Trends Biotechnol. 2009, 27 (6): 368-374. 10.1016/j.tibtech.2009.03.003

- 28.
Malphettes L, Fussenegger M: Impact of RNA interference on gene networks. Metab Eng. 2006, 8 (6): 672-683. 10.1016/j.ymben.2006.07.005

- 29.
Vilar JMG, Kueh HY, Barkai N, Leibler S: Mechanisms of noise-resistance in genetic oscillators. Proc Natl Acad Sci USA. 2002, 99 (9): 5988-5992. 10.1073/pnas.092133899

- 30.
Painter PR, Bliss RD: Reconsideration of the theory of oscillatory repression. J Theor Biol. 1981, 90 (2): 293-298. 10.1016/0022-5193(81)90048-5

- 31.
Goldbeter A: A model for circadian oscillations in the Drosophila period protein (PER). Proc Biol Sci. 1995, 261 (1362): 319-324. 10.1098/rspb.1995.0153

- 32.
Novák B, Tyson JJ: Design principles of biochemical oscillators. Nat Rev Mol Cell Biol. 2008, 9: 981-991. 10.1038/nrm2530

- 33.
Rössler OE: Chaos in abstract kinetics. Two prototypes. Bull Math Biol. 1977, 39: 275-289.

- 34.
Scott SK: Chemical chaos. 1993, Oxford University Press,

- 35.
Tyson JJ, Chen KC, Novak B: Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell. Curr Opin Cell Biol. 2003, 15 (2): 221-31. 10.1016/S0955-0674(03)00017-6

- 36.
Tyson JJ, Albert R, Goldbeter A, Ruoff P, Sible J: Biological switches and clocks. J R Soc Interface. 2008, 5 (Suppl 1): S1-8. 10.1098/rsif.2008.0179.focus

- 37.
Fung E, Wong WW, Suen JK, Bulter T, gu Lee S, Liao JC: A synthetic gene-metabolic oscillator. Nature. 2005, 435 (7038): 118-122. 10.1038/nature03508

- 38.
Garcia-Ojalvo J, Elowitz MB, Strogatz SH: Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing. Proc Natl Acad Sci USA. 2004, 101 (30): 10955-10960. 10.1073/pnas.0307095101

- 39.
Danino T, Mondragón-Palomino O, Tsimring L, Hasty J: A synchronized quorum of genetic clocks. Nature. 2010, 463 (7279): 326-330. 10.1038/nature08753

- 40.
Pearson JE: Complex Patterns in a Simple System. Science. 1993, 261 (5118): 189-192. 10.1126/science.261.5118.189

- 41.
Lewis J: Autoinhibition with transcriptional delay: a simple mechanism for the zebrafish somitogenesis oscillator. Curr Biol. 2003, 13 (16): 1398-1408. 10.1016/S0960-9822(03)00534-7

- 42.
Swinburne IA, Miguez DG, Landgraf D, Silver PA: Intron length increases oscillatory periods of gene expression in animal cells. Genes Dev. 2008, 22 (17): 2342-2346. 10.1101/gad.1696108

- 43.
Widder S, Schicho J, Schuster P: Dynamic patterns of gene regulation I: simple two-gene systems. J Theor Biol. 2007, 246 (3): 395-419. 10.1016/j.jtbi.2007.01.004

- 44.
Scott M, Ingalls B, Kaern M: Estimations of intrinsic and extrinsic noise in models of nonlinear genetic networks. Chaos. 2006, 16 (2): 026107- 10.1063/1.2211787

- 45.
Hasty J, Dolnik M, Rottschäfer V, Collins JJ: Synthetic gene network for entraining and amplifying cellular oscillations. Phys Rev Lett. 2002, 88 (14): 148101- 10.1103/PhysRevLett.88.148101

- 46.
Smolen P, Baxter DA, Byrne JH: Frequency selectivity, multistability, and oscillations emerge from models of genetic regulatory systems. Am J Physiol. 1998, 274: C531-542.

- 47.
Buchler NE, Gerland U, Hwa T: Nonlinear protein degradation and the function of genetic circuits. Proc Natl Acad Sci USA. 2005, 102 (27): 9559-9564. 10.1073/pnas.0409553102

- 48.
Müller S, Hofbauer J, Endler L, Flamm C, Widder S, Schuster P: A generalized model of the repressilator. J Math Biol. 2006, 53 (6): 905-937. 10.1007/s00285-006-0035-9

- 49.
Thomas R, Thieffry D, Kaufman M: Dynamical behaviour of biological regulatory networks-I. Biological role of feedback loops and practical use of the concept of the loop-characteristic state. Bull Math Biol. 1995, 57 (2): 247-276. http://dx.doi.org/10.1016/0092-8240(94)00036-C

- 50.
Thieffry D: Dynamical roles of biological regulatory circuits. Brief Bioinform. 2007, 8 (4): 220-225. 10.1093/bib/bbm028

- 51.
Smith H: Oscillations and multiple steady states in a cyclic gene model with repression. J Math Biol. 1987, 25 (2): 169-190.

- 52.
Griffith JS: Mathematics of cellular control processes. II. Positive feedback to one gene. J Theor Biol. 1968, 20 (2): 209-216. 10.1016/0022-5193(68)90190-2

- 53.
Smolen P, Baxter DA, Byrne JH: Mathematical modeling of gene networks. Neuron. 2000, 26 (3): 567-580. 10.1016/S0896-6273(00)81194-0

- 54.
Barkai N, Leibler S: Circadian clocks limited by noise. Nature. 2000, 403 (6767): 267-268.

- 55.
Mahaffy J: Models of genetic control by repression with time delays and spatial effects. J Math Biol. 1984, 20: 39-57. 10.1007/BF00275860

- 56.
Bratsun D, Volfson D, Tsimring LS, Hasty J: Delay-induced stochastic oscillations in gene regulation. Proc Natl Acad Sci USA. 2005, 102 (41): 14593-14598. 10.1073/pnas.0503858102

- 57.
Tiana G, Krishna S, Jensen MH, Sneppen K: Oscillations and temporal signalling in cells. Phys Biol. 2007, 4: R1-17. 10.1088/1478-3975/4/2/R01

- 58.
François P, Hakim V: Design of genetic networks with specified functions by evolution in silico. Proc Natl Acad Sci USA. 2004, 101 (2): 580-585. 10.1073/pnas.0304532101

- 59.
Paladugu SR, Chickarmane V, Deckard A, Frumkin JP, McCormack M, Sauro HM: In silico evolution of functional modules in biochemical networks. Syst Biol (Stevenage). 2006, 153 (4): 223-235.

- 60.
Batt G, Yordanov B, Weiss R, Belta C: Robustness analysis and tuning of synthetic gene networks. Bioinformatics. 2007, 23 (18): 2415-2422. 10.1093/bioinformatics/btm362

- 61.
Chen BS, Wu CH: A systematic design method for robust synthetic biology to satisfy design specifications. BMC Syst Biol. 2009, 3: 66- 10.1186/1752-0509-3-66

- 62.
Burz DS, Rivera-Pomar R, Jäckle H, Hanes SD: Cooperative DNA-binding by Bicoid provides a mechanism for threshold-dependent gene activation in the Drosophila embryo. EMBO J. 1998, 17 (20): 5998-6009. 10.1093/emboj/17.20.5998

- 63.
Veitia RA: A sigmoidal transcriptional response: cooperativity, synergy and dosage effects. Biol Rev Camb Philos Soc. 2003, 78: 149-170. 10.1017/S1464793102006036

- 64.
Amouyal M: Gene regulation at-a-distance in E. coli: new insights. C R Biol. 2005, 328: 1-9. 10.1016/j.crvi.2004.10.002

- 65.
Strogatz SH: Nonlinear Dynamics and Chaos: With Applications in Physics, Biology, Chemistry, and Engineering. 1994, Addison Wesley,

- 66.
Nicolis G: Introduction to Nonlinear Science. 1995, Cambridge University Press,

- 67.
Izhikevich EM: Dynamical Systems in Neuroscience. The Geometry of Excitability and Bursting. 2006, MIT Press,

- 68.
Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: from theories to phenotypes. Nat Rev Genet. 2005, 6 (6): 451-464. 10.1038/nrg1615

- 69.
Raser JM, O'Shea EK: Noise in gene expression: origins, consequences, and control. Science. 2005, 309 (5743): 2010-2013. 10.1126/science.1105891

- 70.
Pedraza JM, van Oudenaarden A: Noise propagation in gene networks. Science. 2005, 307 (5717): 1965-1969. 10.1126/science.1109090

- 71.
Turcotte M, Garcia-Ojalvo J, Süel GM: A genetic timer through noise-induced stabilization of an unstable state. Proc Natl Acad Sci USA. 2008, 105 (41): 15732-15737. 10.1073/pnas.0806349105

- 72.
Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks. Proc Natl Acad Sci USA. 2001, 98 (15): 8614-8619. 10.1073/pnas.151588598

- 73.
Tsai TYC, Choi YS, Ma W, Pomerening JR, Tang C, Ferrell JE: Robust, tunable biological oscillations from interlinked positive and negative feedback loops. Science. 2008, 321 (5885): 126-129. 10.1126/science.1156951

- 74.
Hasty J, McMillen D, Isaacs F, Collins JJ: Computational studies of gene regulatory networks: in numero molecular biology. Nat Rev Genet. 2001, 2 (4): 268-279. 10.1038/35066056

## Acknowledgements

This work was supported by the EU Grant FP6 SYNLET (NEST-043312). MI is funded by EU Grants: FP7 ERC Zinc-hubs-201249, FP6 Integra-29025; Spanish Ministry of Science (MEC) Grant-SAF2007-61422; The MEC-EMBL Agreement. MC is funded by the grant SFRH/BD/39100/2007 from Fundaçao para a Ciencia e Tecnologia. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Author information

## Additional information

### Authors' contributions

All authors contributed to the design and coordination of the study. AM performed the computational implementations and prepared the original draft, which was revised by RVS, and subsequently by MC and MI. All authors read and approved the final manuscript.

## Electronic supplementary material

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Hopf Bifurcation
- Bifurcation Point
- Synthetic Biology
- Biochemical Mechanism
- Experimental Implementation