Avoiding transcription factor competition at promoter level increases the chances of obtaining oscillation
© Munteanu et al; licensee BioMed Central Ltd. 2010
Received: 10 September 2009
Accepted: 17 May 2010
Published: 17 May 2010
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.
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.
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.
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 , the next step for synthetic biology was the construction of biological oscillators. The first successful implementation  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 . 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 . 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  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 involve a certain number of genes that regulate one another's expression . 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  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  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 . 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.  for a review). More precisely, it was proved  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 ).
Table associated to Figure 1 representing known genetic oscillators that employ transcriptional regulation alone.
Reference of the study and remarks
Experimental: Swinburne et al. (2008) 
Analytic: Widder et al. (2007) . The minimum Hill exponent to make oscillations possible is n = 3.
Experimental: Atkinson et al. (2003) : The notation 2 + 2 refers to the fact that the DNA loops operate when 2 DNA-bound dimers form a tetramer. It is still not clear if the design functions with competition or not, and this is the reason for employing the 2 + 2 notation instead of the arrows.
Numeric: Scott et al. (2006)  on the Atkinson oscillator.
Numeric: Guantes&Poyatos (2006)  associate their Design I to the Atkinson experimental context. The oscillators they propose due their oscillations to the time-scale difference between activator and repressor life-time. Otherwise, higher multimers than dimers are needed.
Numeric: Hasty et al. (2002) 
Experimental: Stricker et a. (2008) 
Numeric: Smolen et al. (1998) . The activator needs to be at least a dimer for the existence of the oscillations.
Experimental and numeric: Elowitz&Liebler (2000) : n = 2 from Figure (1) above makes reference to the value employed by Elowitz&Leibler in their model. See also the discussion in the Supporting Information from Buchler et al. (2005) .
Analytic: Mueller et al. (2006) : a general case.
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 . 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 . 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 . Nevertheless, the relation between delay, protein half-life and cooperativity required for oscillatory solutions (see page 192 of Ref.  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" ; Figure 1H) and its generalization has been studied theoretically . The clarification of the specific parameter ranges for which oscillations occur is crucial for the design of a synthetic circuit, such as the repressilator . 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 : 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 . 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 . 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].
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. . Guantes&Poyatos  claim that Design I corresponds to Atkinson model. And indeed Atkinson et al.  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 . 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 . 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  and Conrad et al. . Guantes&Poyatos  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 . In a nutshell, Guantes&Poyatos  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.  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?
The list of symbols employed in the ordinary differential equations associated to the designs from Figure 2.
Definition & Comments
concentration of the activator (nM)
concentration of the repressor (nM)
model time (hours)
DNA-binding constant of the activator A
DNA-binding constant of the repressor R
dimerization constant of the activator A to A2
dimerization constant of the repressor R to R2
multimerization constant of the activator A to A n
dimerization constant of the repressor R to R m
degradation rate of the activator A
degradation rate of the repressor R
δ A /δR,non-dimensional parameter of the model
, the nondimensional variable associated to the activator, when the activator is
, the nondimensional variable associated to the repressor, when the repressor is a dimer
tδ R , the nondimensional time variable of the model
the transcriptional synergy conferred by the activator bound to DNA
, the nondimensional parameter of the transcriptional response from the activator promoter, with γ A , the translation rate of the activator A; β A , the transcription rate of the mRNA of the activator; α A , activator mRNA degradation; , total activator promoter number. The root order n denotes the multimerization degree of the TF binding to DNA (i.e. n = 2, dimer; n = 4, tetramer). The subscript from , the equilibrium constant for multimer formation, also denotes its dependence on the multimerization of the binding factor.
, the nondimensional parameter of the transcriptional response from the repressor promoter, similar to β. See also . Notice also the multimerization degree appears here as m with similar consequences as for the β parameter.
, the binding ratio of the repressor compared to the activator.
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 . 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 , the current study is the first to inquire on the differences in the dynamical behavior between the two scenarios.
We have chosen fixed values for Δ, α and σ, and varied β and γ, in order to compare to existing results . 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 , 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) , 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.
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. ). 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 ). 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.
Returning to Figure 4, notice also how the behavior of Design III from panel C matches the behavior of Design II from (their Fig. 3b): depending on the parameter to be varied, a design can present both types of bifurcations. Guantes&Poyatos  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. , 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. ) 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  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 . Nevertheless, the clarification of this issue needs further experimental studies.
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 σx2y2 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  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 , noise-induced oscillations . Moreover, there is clear evidence that the transcription and translation level are indeed noise-prone processes . 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. , 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  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.
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.
The model from eqs. (1-2) is a simplified version of the Atkinson model as described by Scott et al. . 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 . Nevertheless, the simplified model is a more tractable tool for studying the general characteristics of the system.
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 . 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. , 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 .
The Atkinson oscillator introduced and experimentally implemented by  has been modeled both through competition  and non-competition models . 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 . 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.
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.
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.
- Haseloff J, Ajioka J: Synthetic biology: history, challenges and prospects. J R Soc Interface. 2009, 6: S389-S391. 10.1098/rsif.2009.0176.focusPubMed CentralView ArticlePubMedGoogle Scholar
- Benner SA, Sismour AM: Synthetic biology. Nat Rev Genet. 2005, 6 (7): 533-543. 10.1038/nrg1637View ArticlePubMedGoogle Scholar
- Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403 (6767): 339-342. 10.1038/35002131View ArticlePubMedGoogle Scholar
- Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403 (6767): 335-338. 10.1038/35002125View ArticlePubMedGoogle Scholar
- Glass L, Mackey MC: From clocks to chaos The rythms of life. 1988, Princeton,Google Scholar
- 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-5View ArticlePubMedGoogle Scholar
- 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/nature07389View ArticlePubMedGoogle Scholar
- Tigges M, Marquez-Lago TT, Stelling J, Fussenegger M: A tunable synthetic mammalian oscillator. Nature. 2009, 457 (7227): 309-312. 10.1038/nature07616View ArticlePubMedGoogle Scholar
- Isalan M, Morrison M: This title is false. Nature. 2009, 458 (7241): 969- 10.1038/458969aView ArticlePubMedGoogle Scholar
- 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.0155View ArticlePubMedGoogle Scholar
- de Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol. 2002, 9: 67-103. 10.1089/10665270252833208View ArticlePubMedGoogle Scholar
- 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.0020439PubMed CentralView ArticlePubMedGoogle Scholar
- Goldbeter A: Computational biology: a propagating wave of interest. Curr Biol. 2004, 14 (15): R601-602. 10.1016/j.cub.2004.07.036View ArticlePubMedGoogle Scholar
- 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.004View ArticlePubMedGoogle Scholar
- 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/nature04473View ArticlePubMedGoogle Scholar
- 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.1536PubMed CentralView ArticlePubMedGoogle Scholar
- Gertz J, Siggia ED, Cohen BA: Analysis of combinatorial cis-regulation in synthetic and genomic promoters. Nature. 2009, 457 (7226): 215-218. 10.1038/nature07521PubMed CentralView ArticlePubMedGoogle Scholar
- Guantes R, Poyatos JF: Dynamical principles of two-component genetic oscillators. PLoS Comput Biol. 2006, 2 (3): e30- 10.1371/journal.pcbi.0020030PubMed CentralView ArticlePubMedGoogle Scholar
- 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.focusPubMed CentralView ArticlePubMedGoogle Scholar
- Lewin B: Genes IX. 2007, Jones & Bartlett,Google Scholar
- Storz G, Altuvia S, Wassarman KM: An abundance of RNA regulators. Annu Rev Biochem. 2005, 74: 199-217. 10.1146/annurev.biochem.74.082803.133136View ArticlePubMedGoogle Scholar
- Davidson EA, Ellington AD: Synthetic RNA circuits. Nat Chem Biol. 2007, 3: 23-28. 10.1038/nchembio846View ArticlePubMedGoogle Scholar
- 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.0703961104PubMed CentralView ArticlePubMedGoogle Scholar
- 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.1345702View ArticlePubMedGoogle Scholar
- Mukherji S, van Oudenaarden A: Synthetic biology: understanding biological design from synthetic circuits. Nat Rev Genet. 2009, 10 (12): 859-871.PubMed CentralPubMedGoogle Scholar
- 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.focusPubMed CentralView ArticlePubMedGoogle Scholar
- Bhalerao KD: Synthetic gene networks: the next wave in biotechnology?. Trends Biotechnol. 2009, 27 (6): 368-374. 10.1016/j.tibtech.2009.03.003View ArticlePubMedGoogle Scholar
- Malphettes L, Fussenegger M: Impact of RNA interference on gene networks. Metab Eng. 2006, 8 (6): 672-683. 10.1016/j.ymben.2006.07.005View ArticlePubMedGoogle Scholar
- 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.092133899PubMed CentralView ArticlePubMedGoogle Scholar
- 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-5View ArticlePubMedGoogle Scholar
- 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.0153View ArticlePubMedGoogle Scholar
- Novák B, Tyson JJ: Design principles of biochemical oscillators. Nat Rev Mol Cell Biol. 2008, 9: 981-991. 10.1038/nrm2530PubMed CentralView ArticlePubMedGoogle Scholar
- Rössler OE: Chaos in abstract kinetics. Two prototypes. Bull Math Biol. 1977, 39: 275-289.View ArticlePubMedGoogle Scholar
- Scott SK: Chemical chaos. 1993, Oxford University Press,Google Scholar
- 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-6View ArticlePubMedGoogle Scholar
- 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.focusPubMed CentralView ArticlePubMedGoogle Scholar
- 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/nature03508View ArticlePubMedGoogle Scholar
- 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.0307095101PubMed CentralView ArticlePubMedGoogle Scholar
- Danino T, Mondragón-Palomino O, Tsimring L, Hasty J: A synchronized quorum of genetic clocks. Nature. 2010, 463 (7279): 326-330. 10.1038/nature08753PubMed CentralView ArticlePubMedGoogle Scholar
- Pearson JE: Complex Patterns in a Simple System. Science. 1993, 261 (5118): 189-192. 10.1126/science.261.5118.189View ArticlePubMedGoogle Scholar
- 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-7View ArticlePubMedGoogle Scholar
- 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.1696108PubMed CentralView ArticlePubMedGoogle Scholar
- 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.004View ArticlePubMedGoogle Scholar
- 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.2211787View ArticlePubMedGoogle Scholar
- 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.148101View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.0409553102PubMed CentralView ArticlePubMedGoogle Scholar
- 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-9View ArticlePubMedGoogle Scholar
- 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-CView ArticlePubMedGoogle Scholar
- Thieffry D: Dynamical roles of biological regulatory circuits. Brief Bioinform. 2007, 8 (4): 220-225. 10.1093/bib/bbm028View ArticlePubMedGoogle Scholar
- Smith H: Oscillations and multiple steady states in a cyclic gene model with repression. J Math Biol. 1987, 25 (2): 169-190.View ArticlePubMedGoogle Scholar
- 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-2View ArticlePubMedGoogle Scholar
- Smolen P, Baxter DA, Byrne JH: Mathematical modeling of gene networks. Neuron. 2000, 26 (3): 567-580. 10.1016/S0896-6273(00)81194-0View ArticlePubMedGoogle Scholar
- Barkai N, Leibler S: Circadian clocks limited by noise. Nature. 2000, 403 (6767): 267-268.PubMedGoogle Scholar
- Mahaffy J: Models of genetic control by repression with time delays and spatial effects. J Math Biol. 1984, 20: 39-57. 10.1007/BF00275860View ArticlePubMedGoogle Scholar
- 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.0503858102PubMed CentralView ArticlePubMedGoogle Scholar
- 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/R01View ArticlePubMedGoogle Scholar
- 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.0304532101PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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/btm362View ArticlePubMedGoogle Scholar
- 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-66PubMed CentralView ArticlePubMedGoogle Scholar
- 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.5998PubMed CentralView ArticlePubMedGoogle Scholar
- Veitia RA: A sigmoidal transcriptional response: cooperativity, synergy and dosage effects. Biol Rev Camb Philos Soc. 2003, 78: 149-170. 10.1017/S1464793102006036View ArticlePubMedGoogle Scholar
- 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.002View ArticlePubMedGoogle Scholar
- Strogatz SH: Nonlinear Dynamics and Chaos: With Applications in Physics, Biology, Chemistry, and Engineering. 1994, Addison Wesley,Google Scholar
- Nicolis G: Introduction to Nonlinear Science. 1995, Cambridge University Press,View ArticleGoogle Scholar
- Izhikevich EM: Dynamical Systems in Neuroscience. The Geometry of Excitability and Bursting. 2006, MIT Press,Google Scholar
- 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/nrg1615View ArticlePubMedGoogle Scholar
- Raser JM, O'Shea EK: Noise in gene expression: origins, consequences, and control. Science. 2005, 309 (5743): 2010-2013. 10.1126/science.1105891PubMed CentralView ArticlePubMedGoogle Scholar
- Pedraza JM, van Oudenaarden A: Noise propagation in gene networks. Science. 2005, 307 (5717): 1965-1969. 10.1126/science.1109090View ArticlePubMedGoogle Scholar
- 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.0806349105PubMed CentralView ArticlePubMedGoogle Scholar
- Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks. Proc Natl Acad Sci USA. 2001, 98 (15): 8614-8619. 10.1073/pnas.151588598PubMed CentralView ArticlePubMedGoogle Scholar
- 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.1156951PubMed CentralView ArticlePubMedGoogle Scholar
- 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/35066056View ArticlePubMedGoogle Scholar