Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

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

Andreea Munteanu1*, Marco Constante2, Mark Isalan2 and Ricard V Solé13

Author Affiliations

1 ICREA-Complex Systems Lab, Universitat Pompeu Fabra (PRBB-GRIB), Dr Aiguader 88, 08003 Barcelona, Spain

2 EMBL-CRG Systems Biology Research Unit, Centre for Genomic Regulation (CRG), UPF, Dr Aiguader 88, 08003 Barcelona, Spain

3 Santa Fe Institute, 1399 Hyde Park Road, Santa Fe NM 87501, USA

For all author emails, please log on.

BMC Systems Biology 2010, 4:66  doi:10.1186/1752-0509-4-66


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/4/66


Received:10 September 2009
Accepted:17 May 2010
Published:17 May 2010

© 2010 Munteanu et al; licensee 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.

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.

Additional file 1. Transcriptional regulation and model description. A PDF file containing a brief introduction to transcriptional regulation and the terminology employed in the main text. The details necessary for obtaining the equations employed in this study are also included here.

Format: PDF Size: 139KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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.

thumbnailFigure 1. Transcriptional genetic circuits from the literature. Figure illustrating the circuit architectures associated to Table 1. In the graphical representation, we have made the distinction between non-competition and competition models by solid and dashed arrowed lines, respectively. In the notation n1, n2, ... characterizing the interactions in the schematic figures, the comma separates various binding sites and the degree of the multimer binding to that site. The notation n1 + n2 is employed only for subfigure (c) - see discussion in the Table 1.

Table 1. Table associated to Figure 1 representing known genetic oscillators that employ transcriptional regulation alone.

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.

thumbnailFigure 2. Distinguishing competition versus non-competition of transcription factors at promoter level in a chosen circuit architecture. Three different implementations compatible with the scheme in (A): Design I - inhibition through competition at the activator promoter (B), Design II -inhibition through active degradation of the activator by the repressor (C), and Design III - specific inhibition (d). The cases (B) and (C) are discussed in [18], while there is controversy whether the Atkinson oscillator [6] is based on case (B) [18] or (D) [44].

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/, 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.

Table 2. The list of symbols employed in the ordinary differential equations associated to the designs from Figure 2.

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.

Additional file 2. The C++ codes developed for the numerical simulations, as well as the instructions of their use.

Format: ZIP Size: 45KB Download fileOpen Data

thumbnailFigure 3. Oscillatory features of the competition and non-competition designs associated to the studied circuit architecture. The parameter space (β, γ) and the values of the oscillation period for (A) Design I and (B) Design III. The color code represents the period of oscillations depending on the (β, γ) pair, with α = 50, Δ = 10, σ = 1. Dark blue implies no sustained oscillations. As defined in the main text, the period in the color bar is expressed in units of δR, the degradation rate of the repressor.

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.

Additional file 3. Oscillatory region for an extended parameter space for Design III.

Format: PDF Size: 178KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4. Comparison between Design I and Design III for Δ = 20.

Format: PDF Size: 309KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 5. The behavior of Design II from Guantes&Poyatos[18].

Format: PDF Size: 179KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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.

thumbnailFigure 4. The bifurcation diagram with β as control parameter. The bifurcation diagram obtained by varying β, with α = 50 and Δ = 10, (A) and (C) γ = 0.01 and (B) and (D) γ = 0.05. It corresponds to the crossing of Figures 3A and 3B at constant γ. In black, we plot the stable (solid) and unstable (dashed) fixed points, xo, and in red and blue, the maximum and minimum, respectively, of the stable limit cycle. Notice the change in the way oscillations are triggered depending on the fixed parameters. Here Hopf(b) and Hopf(p) indicate subcritical and supercritical bifurcation, respectively. Gray areas represent damped oscillations.

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).

Additional file 6. Example of a subcritical Hopf bifurcation.

Format: PDF Size: 48KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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.

Additional file 7. The trajectory behavior prior to a Hopf and to a SNIC bifurcation.

Format: PDF Size: 153KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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.

thumbnailFigure 5. The bifurcation diagram with Δ as control parameter. Similar to Figure 4, but having Δ as control parameter. The fixed parameters are α = 50, and β = 1.58, γ = 0.079 for (A) and (C), and α = 50, and β = 0.2, γ = 0.04 for (B) and (D).

thumbnailFigure 6. An example of time series for the two studied designs. Time series for Design I and Design III for the case α = 50, Δ = 11, β = 1.58, γ = 0.079. Notice the pulse-like oscillation of Design I compared to a more sinusoidal one for Design III, as a consequence of the type of bifurcation they originated from (Figure 5A and 5C).

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.

Additional file 8. Robustness to parameter changes.

Format: PDF Size: 796KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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 [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 yx and yy 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 yx and yy, 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 yx 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 yx , and an increase in Δ or γ implies upwards displacement of yy. 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.

Additional file 9. Dynamical differences between the designs: discussion of the nullclines.

Format: PDF Size: 142KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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.

Additional file 10. The influence of internal noise.

Format: PDF Size: 119KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

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:

(1)

(2)

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:

(3)

(4)

thumbnailFigure 7. Oscillatory features of the simplified Atkinson model. The behavior of the simplified model of the Atkinson oscillator from eqs. (1-2) for α = 50, Δ = 1, n = 3 and (A) m = 1,(B) m = 2 and (C) m = 3.

thumbnailFigure 8. The scheme of the experimental design from [6].

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.

thumbnailFigure 9. Oscillatory features of the Atkinson model. The numerical simulations for the corresponding model in eqs. (3-4) and Figure 8 with Δ = 10 and α = 50. When compared to Figure 3A, it is visible that the oscillatory region is significantly reduced.

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):

(5)

(6)

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.

thumbnailFigure 10. Oscillatory features of the circuit introduced by Smolen et al. [46]. The behavior of the oscillator from eqs. (5-6) introduced in [46] with α = 50, Δ = 4, σ = 1. The scale differs from that of Figures 3 and 7, but nevertheless notice the difference in the location and extension of the oscillatory region when compared to these figures.

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.

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.

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.

References

  1. Haseloff J, Ajioka J: Synthetic biology: history, challenges and prospects.

    J R Soc Interface 2009, 6:S389-S391. PubMed Abstract | Publisher Full Text OpenURL

  2. Benner SA, Sismour AM: Synthetic biology.

    Nat Rev Genet 2005, 6(7):533-543. PubMed Abstract | Publisher Full Text OpenURL

  3. Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli.

    Nature 2000, 403(6767):339-342. PubMed Abstract | Publisher Full Text OpenURL

  4. Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators.

    Nature 2000, 403(6767):335-338. PubMed Abstract | Publisher Full Text OpenURL

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

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  8. Tigges M, Marquez-Lago TT, Stelling J, Fussenegger M: A tunable synthetic mammalian oscillator.

    Nature 2009, 457(7227):309-312. PubMed Abstract | Publisher Full Text OpenURL

  9. Isalan M, Morrison M: This title is false.

    Nature 2009, 458(7241):969. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  11. de Jong H: Modeling and simulation of genetic regulatory systems: a literature review.

    J Comput Biol 2002, 9:67-103. PubMed Abstract | Publisher Full Text OpenURL

  12. Cohen JE: Mathematics is biology's next microscope, only better; biology is mathematics' next physics, only better.

    PLoS Biol 2004, 2(12):e439. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Goldbeter A: Computational biology: a propagating wave of interest.

    Curr Biol 2004, 14(15):R601-602. PubMed Abstract | Publisher Full Text OpenURL

  14. Mogilner A, Wollman R, Marshall WF: Quantitative modeling in cell biology: what is it good for?

    Dev Cell 2006, 11(3):279-287. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Gertz J, Siggia ED, Cohen BA: Analysis of combinatorial cis-regulation in synthetic and genomic promoters.

    Nature 2009, 457(7226):215-218. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Guantes R, Poyatos JF: Dynamical principles of two-component genetic oscillators.

    PLoS Comput Biol 2006, 2(3):e30. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Lewin B: Genes IX. Jones & Bartlett; 2007. OpenURL

  21. Storz G, Altuvia S, Wassarman KM: An abundance of RNA regulators.

    Annu Rev Biochem 2005, 74:199-217. PubMed Abstract | Publisher Full Text OpenURL

  22. Davidson EA, Ellington AD: Synthetic RNA circuits.

    Nat Chem Biol 2007, 3:23-28. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Hasty J, Isaacs F, Dolnik M, McMillen D, Collins JJ: Designer gene networks: Towards fundamental cellular control.

    Chaos 2001, 11:207-220. PubMed Abstract | Publisher Full Text OpenURL

  25. Mukherji S, van Oudenaarden A: Synthetic biology: understanding biological design from synthetic circuits.

    Nat Rev Genet 2009, 10(12):859-871. PubMed Abstract | Publisher Full Text OpenURL

  26. Boyle PM, Silver PA: Harnessing nature's toolbox: regulatory elements for synthetic biology.

    J R Soc Interface 2009, 6:S535-S546. PubMed Abstract | Publisher Full Text OpenURL

  27. Bhalerao KD: Synthetic gene networks: the next wave in biotechnology?

    Trends Biotechnol 2009, 27(6):368-374. PubMed Abstract | Publisher Full Text OpenURL

  28. Malphettes L, Fussenegger M: Impact of RNA interference on gene networks.

    Metab Eng 2006, 8(6):672-683. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Painter PR, Bliss RD: Reconsideration of the theory of oscillatory repression.

    J Theor Biol 1981, 90(2):293-298. PubMed Abstract | Publisher Full Text OpenURL

  31. Goldbeter A: A model for circadian oscillations in the Drosophila period protein (PER).

    Proc Biol Sci 1995, 261(1362):319-324. PubMed Abstract | Publisher Full Text OpenURL

  32. Novák B, Tyson JJ: Design principles of biochemical oscillators.

    Nat Rev Mol Cell Biol 2008, 9:981-991. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Rössler OE: Chaos in abstract kinetics. Two prototypes.

    Bull Math Biol 1977, 39:275-289. PubMed Abstract OpenURL

  34. Scott SK: Chemical chaos. Oxford University Press; 1993. OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Danino T, Mondragón-Palomino O, Tsimring L, Hasty J: A synchronized quorum of genetic clocks.

    Nature 2010, 463(7279):326-330. PubMed Abstract | Publisher Full Text OpenURL

  40. Pearson JE: Complex Patterns in a Simple System.

    Science 1993, 261(5118):189-192. PubMed Abstract | Publisher Full Text OpenURL

  41. Lewis J: Autoinhibition with transcriptional delay: a simple mechanism for the zebrafish somitogenesis oscillator.

    Curr Biol 2003, 13(16):1398-1408. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  44. Scott M, Ingalls B, Kaern M: Estimations of intrinsic and extrinsic noise in models of nonlinear genetic networks.

    Chaos 2006, 16(2):026107. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. [http://dx.doi.org/10.1016/0092-8240(94)00036-C] webcite

    Bull Math Biol 1995, 57(2):247-276. PubMed Abstract OpenURL

  50. Thieffry D: Dynamical roles of biological regulatory circuits.

    Brief Bioinform 2007, 8(4):220-225. PubMed Abstract | Publisher Full Text OpenURL

  51. Smith H: Oscillations and multiple steady states in a cyclic gene model with repression.

    J Math Biol 1987, 25(2):169-190. PubMed Abstract OpenURL

  52. Griffith JS: Mathematics of cellular control processes. II. Positive feedback to one gene.

    J Theor Biol 1968, 20(2):209-216. PubMed Abstract | Publisher Full Text OpenURL

  53. Smolen P, Baxter DA, Byrne JH: Mathematical modeling of gene networks.

    Neuron 2000, 26(3):567-580. PubMed Abstract | Publisher Full Text OpenURL

  54. Barkai N, Leibler S: Circadian clocks limited by noise.

    Nature 2000, 403(6767):267-268. PubMed Abstract | Publisher Full Text OpenURL

  55. Mahaffy J: Models of genetic control by repression with time delays and spatial effects.

    J Math Biol 1984, 20:39-57. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Tiana G, Krishna S, Jensen MH, Sneppen K: Oscillations and temporal signalling in cells.

    Phys Biol 2007, 4:R1-17. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract OpenURL

  60. Batt G, Yordanov B, Weiss R, Belta C: Robustness analysis and tuning of synthetic gene networks.

    Bioinformatics 2007, 23(18):2415-2422. PubMed Abstract | Publisher Full Text OpenURL

  61. Chen BS, Wu CH: A systematic design method for robust synthetic biology to satisfy design specifications.

    BMC Syst Biol 2009, 3:66. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Veitia RA: A sigmoidal transcriptional response: cooperativity, synergy and dosage effects.

    Biol Rev Camb Philos Soc 2003, 78:149-170. PubMed Abstract | Publisher Full Text OpenURL

  64. Amouyal M: Gene regulation at-a-distance in E. coli: new insights.

    C R Biol 2005, 328:1-9. PubMed Abstract | Publisher Full Text OpenURL

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

  66. Nicolis G: Introduction to Nonlinear Science. Cambridge University Press; 1995. OpenURL

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

  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. PubMed Abstract | Publisher Full Text OpenURL

  69. Raser JM, O'Shea EK: Noise in gene expression: origins, consequences, and control.

    Science 2005, 309(5743):2010-2013. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Pedraza JM, van Oudenaarden A: Noise propagation in gene networks.

    Science 2005, 307(5717):1965-1969. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks.

    Proc Natl Acad Sci USA 2001, 98(15):8614-8619. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL