Negative auto-regulation increases the input dynamic-range of the arabinose system of Escherichia coli
© Madar et al; licensee BioMed Central Ltd. 2011
Received: 28 April 2011
Accepted: 12 July 2011
Published: 12 July 2011
Gene regulation networks are made of recurring regulatory patterns, called network motifs. One of the most common network motifs is negative auto-regulation, in which a transcription factor represses its own production. Negative auto-regulation has several potential functions: it can shorten the response time (time to reach halfway to steady-state), stabilize expression against noise, and linearize the gene's input-output response curve. This latter function of negative auto-regulation, which increases the range of input signals over which downstream genes respond, has been studied by theory and synthetic gene circuits. Here we ask whether negative auto-regulation preserves this function also in the context of a natural system, where it is embedded within many additional interactions. To address this, we studied the negative auto-regulation motif in the arabinose utilization system of Escherichia coli, in which negative auto-regulation is part of a complex regulatory network.
We find that when negative auto-regulation is disrupted by placing the regulator araC under constitutive expression, the input dynamic range of the arabinose system is reduced by 10-fold. The apparent Hill coefficient of the induction curve changes from about n = 1 with negative auto-regulation, to about n = 2 when it is disrupted. We present a mathematical model that describes how negative auto-regulation can increase input dynamic-range, by coupling the transcription factor protein level to the input signal.
Here we demonstrate that the negative auto-regulation motif in the native arabinose system of Escherichia coli increases the range of arabinose signals over which the system can respond. In this way, negative auto-regulation may help to increase the input dynamic-range while maintaining the specificity of cooperative regulatory systems. This function may contribute to explaining the common occurrence of negative auto-regulation in biological systems.
NAR has been suggested experimentally and theoretically to have several functions. The first is increased homeostasis or buffering of the auto-regulated gene product concentration against stochastic noise [12–14]. Because protein levels can vary from cell to cell by tens of percents [15–17], such a noise buffering mechanism is useful in cases where precision in TF levels is needed . Low frequency noise in TF production rates tends to be buffered by NAR because negative feedback reduces TF levels if they are too high, and increases them if they are too low, making TF levels more uniform across cells.
A second feature of NAR is its ability to speed the response time of gene circuits [6, 19, 20]. Response time is defined as the time it takes to reach half of the total change in a dynamic process. Theoretical comparison between NAR and a simply regulated promoter with no NAR, with parameters in which both reach the same steady-state level expression, shows that the response time is shorter when the TF is negatively auto-regulated. This speed up is achieved by the use of a strong promoter allowing a rapid initial rise in TF levels, up to its auto-repression threshold, followed by reduction in production rate due to NAR . This speedup feature was observed in a synthetic NAR circuit  as well as in the native SOS DNA repair system of E. coli. Speedup offered by NAR may be advantageous in dynamic environments where rapid responses improve fitness.
Recently, it was shown by Nevozhay et al. that NAR can also linearize dose responses . In this study the response of synthetic, TetR-based transcriptional circuits with and without NAR was studied in S. cerevisiae as a function of inducer (anhydrotetracycline, aTc) levels. NAR was found to transform a sigmoid induction curve into a more linear curve (see also ). This feature was also suggested in theoretical studies [6, 22–25]. This role of NAR can be interpreted as an increase in the input dynamic range - the range of input signals over which the system can respond.
Such theoretical and synthetic-circuit studies are a powerful approach because one can study the function of circuits such as NAR without of interfering effects. In natural systems, however, this motif is embedded inside a large regulatory network with many other interactions. These additional interactions might in principle modify its function. Therefore, to fully test the function of a motif requires, in addition to theory and synthetic circuits, experiments on the motif in its natural context, wired into the full interaction networks of the cell.
Here, we study the function of the NAR motif in a natural system. We chose one of the best-studied gene regulation systems, the arabinose utilization system of E. coli. This system has been characterized over the past decades by Schleif and colleagues ([26–29] for reviews). The arabinose-responsive TF, called AraC, is negatively autoregulated (Figure 1b). We asked whether NAR increases the input dynamic range in this system.
The arabinose system is composed of 9 genes arranged in 5 operons: araC- the system-specific TF; araE, araFGH, araJ- the arabinose transporters [30–32]; and araBAD- arabinose catabolic enzymes. Two operons, araC and araBAD, are divergent and share the same regulatory region (Figure 1d). The system is regulated by cAMP Receptor Protein (CRP) and AraC (Figure 1e), which are activated by cAMP and L-Arabinose respectively [26–29]. AraC represses its own promoter, creating a NAR motif. It both activates and represses the arabinose utilization operon araBAD by means of a DNA looping mechanism [33, 34]. AraC undergoes a conformational change when it binds L-Arabinose, leading to expression of the ara genes. The system includes several interactions and feedback loops, in which metabolic enzymes and transporters downstream of araC affect the level of intracellular arabinose, the inducer that activates AraC (Figure 1e). In a study of the input functions of E. coli sugar systems, it was recently found that promoters from the arabinose system respond to their inducer with a wider input dynamic range compared to other sugar systems (eg. the maltose system) in which the TF is not negatively auto-regulated .
To test the role of NAR, we compared the wild type ara system (Figure 1b) to a variant in which NAR is disrupted by placing the regulator AraC under a constitutive promoter (Figure 1c). We find, using high-temporal resolution measurement of promoter activity, that disrupting NAR in the arabinose system increases the steepness of the sigmoidal response curve. It reduces the input dynamic range by about an order of magnitude. Thus, NAR increases input dynamic range in the context of the natural ara system. We also analyze this mathematically, suggesting that the increase in input dynamic range is due to the increase of AraC protein level with increasing arabinose due to the NAR.
The native input function of the araBAD promoter has a broad input dynamic range
In order to determine the input dynamic range of E. coli promoters we used a fluorescent reporter automated assay [35, 39], with strains from the comprehensive E. coli transcription reporter library . Each strain bears a low-copy plasmid with a green fluorescent protein gene under the control of a full length copy of the promoter of interest. In this study we used reporters for the araBAD and araC promoters in E. coli strain MG1655 (see Methods).
Reporter strains were grown on glucose minimal medium containing saturating amount of cAMP (30 mM, to fully activate CRP) and increasing amounts of L-Arabinose [35, 41]. Promoter activity (PA) was defined as the rate of GFP production per OD (optical density) unit, PA = dGFP/dt/OD (see Methods). The input functions were derived from the promoter activity averaged over a window that spans 1-2 cell generations in exponential phase (5-7 hours after initial 1:600 inoculation). Over this time window, promoter activity was constant to a good approximation (see Additional File 1 for fluorescence and growth curves, p. 2-3, Figure S1 and S2 respectively).
The araC gene is induced by arabinose
Disruption of negative auto-regulation of araC reduces the input dynamic range of araBAD
To study the role of the negative auto-regulation of araC on the input dynamic range of its downstream genes, we decoupled araC expression from its negative auto-regulation (Figure. 1c). For this purpose we deleted the araC open reading frame from the chromosome of the wild-type strain MG1655 and re-introduced araC on a plasmid (pZE11) which provides constitutive expression (strain U426, see Methods). The plasmid has a tetR controlled promoter, repressed by a chromosomal tetR gene. Without induction, this plasmid produces levels of AraC that are comparable to the induced wild-type AraC level, as assessed from the maximal promoter activity of the araBAD reporter. It should be noted that the parental strain in this study (U424) also contains chromosomal tetR as well as an emptly pZE11 vector, in order to preserve genotypic identity between the two strains.
We find that in the absence of NAR, the arabinose-dependent input function of araBAD is significantly steeper than the parental input function (Figure 3b,c), with an apparent Hill coefficient of n = 1.9 ± 0.4 (s.e.), and halfway induction point of K = 42 ± 0.6 mM (s.e.). The measured input dynamic range spans between 14 mM - 135 mM, and thus has R = 10 ± 3 (s.e.), in comparison to R = 100 ± 40 (s.e.) in the parental strain. Thus, decoupling araC from its negative auto-regulation reduces the input dynamic-range of its downstream genes by about an order of magnitude.
A model of NAR and increased input dynamic range
What is the main effect at play that allows negative auto-regulation to increase input dynamic range? To understand this, we analyzed a mathematical model of the NAR motif. We sought to make the model as simple as possible, in order to be able to understand it intuitively, and at the same time not too simple so as not to lose the essence of the problem. A more comprehensive model, based on mass-action kinetics, which includes a dual transcription factor that acts as both a repressor and an activator, is given in Additional File 1 (p. 5-7).
Where K z is the dissociation constant of X * from the promoter of Z, β Z is the maximal production rate of Z, and α is its degradation/dilution rate .
It is at this point that negative auto-regulation has an important effect: instead of a constitutive level of X, negative auto-regulation allows the signal s to modify the concentration of X, an effect termed direct coupling . With negative auto-regulation of the type found in the ara system, the promoter that encodes X is repressed by free X, (denoted X f ) a repression which is relieved when X is bound to the signal.
Where A 2 = K x β x /α. In other words, the transcription factor (X) levels increases with the signal (s) levels (see the relationship between AraC and L-Arabinose in Figure 4).
Where . This function has an input dynamic range of R = 81, which is 9 fold wider than that of Eq. (3). Thus, NAR increases the input dynamic range.
Note that the assumption K x <<X f is not crucial for the increased input dynamic range, and was used only for the sake of simplicity. In Additional File 1 we present a full mass-action model, without these assumptions, and show that the present considerations apply as well.
Furthermore, the model explains how the change in regulator levels caused by NAR can cause a shift in the halfway induction point K of downstream genes (relative to no NAR). The direction and size of the shift depends on the mode of regulation. For repressors, K generally increases with regulator levels, whereas activators show the converse dependence. Since AraC both activates and represses araBAD, the detailed model in Additional File 1 explains the observed increase in K shown in Figure 3 (Additional File 1, p.8-9, Figure S3).
To summarize the conclusions of this analysis, negative auto-regulation causes regulator levels to increase with inducer level. This enhances the input dynamic range by extending the range of inputs that can affect the downstream genes.
This study supports a role for negative auto-regulation in increasing the input dynamic range of downstream genes. Previous studies suggested this role theoretically [6, 22–24] and demonstrated it using synthetic circuits [18, 21]. Here we tested NAR in the context of a natural system, the arabinose system of E. coli, in which NAR is embedded within multiple feedback loops and regulatory interactions. Disruption of the NAR in the arabinose system reduced the input dynamic range by an order of magnitude.
What is the intuitive explanation for the enhancement of the input dynamic range provided by negative auto-regulation? Negative auto-regulation in the arabinose system allows the transcription factor concentration to be modulated by its own input signal. As the concentration of input signal increases, the concentration of transcription factor also increases. This extends the response range of downstream promoters, which would otherwise reach maximal activity when the transcription factor becomes saturated with input signal.
A related but distinct feature was studied by M. Savageau [6, 42], in which proper coupling of inducer levels and transcription factor levels can increase the output (as opposed to input) dynamic range of genes: the ratio of their maximal to minimal expression level.
Use of NAR to increase input dynamic range might be especially useful for regulators that bind the promoter cooperatively. Such cooperative binding is thought to increase specificity . However, a well-known feature of cooperative binding (high Hill coefficient) is a narrow input dynamic range . NAR is a simple way to provide wide input dynamic range, while maintaining cooperativity at the level of regulator binding. The combination of cooperativity and negative auto-regulation might thus provide a response across several decades of input strength and at the same time remains specific.
The present study adds to our understanding of the functions of negative auto-regulation network motif, showing that it can increase the input dynamic range of the response, even when embedded in a relatively complicated native gene circuit. Integration of negative auto-regulation within a system with high cooperativity (high specificity and steep activation curve), enables the system to respond to a wide range of input signal (making the activation curve wide) while maintaining the system's specificity to the signal. This function can be experimentally tested in the numerous additional gene systems which bear this network motif across organisms. Because the negative auto-regulation motif is not limited to transcription networks this feature might also apply to other biological systems including protein-level interactions.
Plasmids (see Table 1)
Plasmids and strains used in this study
sc101 ori, promoterless version of the GFP reporter plasmid, (kanR).
P araBAD reporter
GFP reporter plasmid for the araBAD promoter (P araBAD in pUA66), (kanR).
P araC reporter
GFP reporter plasmid for the araC promoter, (P araC in pUA66), (kanR).
Control plasmid: colE1 ori, P LtetO-1 promoter, (ampR) .
araC controlled by the tet promoter on pZE11, (ampR).
MG1655z1: MG1655 (F- lambda- ilvG- rfb-50 rph-1) with chromosomal tetR, (specR).
U423 +pZE11 +P araBAD reporter plasmid (specR, ampR, kanR).
U423 with ΔaraC chromosomal deletion (specR).
U425 +pZE11-araC +P araBAD reporter plasmid (specR, ampR, kanR).
U423 +pZE11 +pUA66 (specR, ampR, kanR).
MG1655 (F- lambda- ilvG- rfb-50 rph-1) + P araC reporter plasmid (kanR).
MG1655 (F- lambda- ilvG- rfb-50 rph-1) + P araBAD reporter plasmid (kanR).
MG1655 (F- lambda- ilvG- rfb-50 rph-1) + pUA66 (kanR).
GFP reporter plasmids (pUA66 based , sc101 ori, kanR, with gfpmut2) for the araC (coordinates 69973-> 70469) and araBAD (coordinates 70469-> 69973) promoters are from the fluorescent reporter library given in detail in . In short: the intergenic region between araBAD and araC (Figure 1d), with more than 100 bps of both flanking regions, was incorporated twice into the GFP reporter plasmid: once in the plus strand orientation (araC promoter) and once in the minus strand orientation (araBAD promoter).
araC was decoupled from its native regulation by cloning it into the pZE11 plasmid (colE1 ori, ampR, P LtetO-1 ) using the Kpn I and Hind III restriction enzymes. The araC gene (the entire coding region) was PCR amplified from MG1655 genomic DNA with the following start and end coordinates: 70387-71299 (positive strand), by using the following primers: 5' ggcggtaccatggctgaagcgcaaaatgatcc for the 5' end and 5' ggcaagcttccgtcaagccgtcaattgtctg for the 3' end. The PCR product and the pZE11 plasmid were digested with Kpn I and Hind III, and then were ligated, yielding pZE11-araC. A self-ligated pZE11 was generated as well to serve as a control plasmid.
Strains (see Table 1)
In order to achieve maximal genotypic identity between the wild-type (with NAR) and the mutant (without NAR) strains, a modified wild-type strain was constructed. tetR gene (z1, specR) from DH5αZ1 was P1 transduced into the wild-type MG1655 chromosome (K12 strain MG1655: F- lambda- ilvG- rfb-50 rph-1), yielding strain U423 (MG1655z1). The pZE11 plasmid and the araBAD reporter plasmid were transformed into U423, yielding strain U424.
An isogenic ΔaraC strain was obtained by deleting araC from the MG1655z1 chromosome (coordinates 70391-> 71244) using the phage λ Red recombination system [46, 47], yielding strain U425 (MGz1ΔaraC). Primers 5' ggacaattggtttcttctctgaatggtgggagtatgaaaagtatggtgtaggctggagctgcttc 3' (for the 5 prime end) and 5' gccgtcaattgtctgattcgttaccaattatgacaacttgacggctaccatatgaatatcctccttag 3' (for the 3 prime end) were used to amplify the kanamycin resistance gene from the pKD4 plasmid with extensions homologous to the 5' and 3' ends of the araC gene, to allow recombination. Kanamycin resistance was removed from the deleted strain using FLP recombinase, as described . U425 did not grow on L-Arabinose as a sole carbon source. The araC deletion and the integrity of the araC\araBAD divergent chromosomal promoter were verified using PCR and sequencing of the scar region. araBAD reporter and pZE11-araC plasmids were transformed into U425, yielding strain U426. Transformation of pZE11-araC into the U425 strain restored its ability to grow on L-Arabinose as a sole carbon source. This strain produced AraC levels, similar to that of the wild-type strain (assessed from the promoter activity of the araBAD reporter at maximal induction, which was about 70% of that of the wild-type strain).
MG1655z1 with empty-pZE11 and pUA66 promoterless reporter plasmid (strain U427) was used for fluorescence background subtraction for U424 & U426.
Growth conditions and measurements
Strains were grown over-night in M9 minimal medium containing 0.4% glucose 0.05% casamino acids, 50 μg/ml kanamycin and 100 μg/ml ampicillin (dictated by the plasmids in each strain) at 37°C. No aTc was used to induce pZE11-araC, since its basal expression level was found to be close to the wild type AraC level. Using a robotic liquid handler (Freedom Evo, Tecan), flat bottom 96-well plates (Nunc) were prepared with 150 μl of M9 minimal medium containing 0.2% glucose 0.05% casamino acids, 30 mM cAMP, 50 μg/ml ampicillin and 25 μg/ml. L-Arabinose, in increasing concentrations was added. The wells were inoculated with the reporter strain at a 1:600 dilution from the overnight culture. Wells were then covered with 100 μl of mineral oil (Sigma) to prevent evaporation (a step which we previously found not to significantly affect aeration or growth [48, 49], and transferred into an automated incubator. Cells were grown in an incubator with shaking (6 hz) at 30°C for about 20 hr. Every 8 minutes the plate was transferred by the robotic arm into a multi-well fluorimeter (Infinite F200, Tecan) that read OD (600 nm) and GFP fluorescence (535 nm).
Promoter activity for each well was calculated from the OD and GFP measurements after subtracting the OD and GFP backgrounds. GFP background was obtained for each well from the promoterless control strains U427 (for strains U424 and U426) and U66 (for strain U428) (Additional File 1, Figure S1). Promoter activity was calculated by computing the rate of accumulation of GFP per unit time divided by the OD (dGFP/dt/OD) as described .
cyclic adenosine mono phosphate
cAMP receptor protein
green fluorescent protein
We thank all of our group members for fruitful comments and discussions. DM thanks R. Milo for his help. This work was supported by the European Research Council, and the Israel Science Foundation.
- Alon U: An Introduction to Systems Biology: Design Principles of Biological Circuits. 2006, Chapman and Hall/CRC, 1,
- Shen-Orr SS, Milo R, Mangan S, Alon U: Network motifs in the transcriptional regulation network of Escherichia coli. Nat Genet. 2002, 31: 64-68. 10.1038/ng881View ArticlePubMed
- Alon U: Network motifs: theory and experimental approaches. Nat Rev Genet. 2007, 8: 450-461. 10.1038/nrg2102View ArticlePubMed
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298: 824-827. 10.1126/science.298.5594.824View ArticlePubMed
- Xiong W, Ferrell JE: A positive-feedback-based bistable/`memory module/' that governs a cell fate decision. Nature. 2003, 426: 460-465. 10.1038/nature02089View ArticlePubMed
- Savageau MA: Comparison of classical and autogenous systems of regulation in inducible operons. Nature. 1974, 252: 546-549. 10.1038/252546a0View ArticlePubMed
- Thieffry D, Huerta AM, Pérez-Rueda E, Collado-Vides J: From specific gene regulation to genomic networks: a global analysis of transcriptional regulation in Escherichia coli. Bioessays. 1998, 20: 433-440. 10.1002/(SICI)1521-1878(199805)20:5<433::AID-BIES10>3.0.CO;2-2View ArticlePubMed
- Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, Zeitlinger J, Jennings EG, Murray HL, Gordon DB, Ren B, Wyrick JJ, Tagne JB, Volkert TL, Fraenkel E, Gifford DK, Young RA: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298: 799-804. 10.1126/science.1075090View ArticlePubMed
- Kim TH, Kim J, Heslop-Harrison P, Cho KH: Evolutionary design principles and functional characteristics based on kingdom-specific network motifs. Bioinformatics. 2011, 27: 245-251. 10.1093/bioinformatics/btq633View ArticlePubMed
- Graham TGW, Tabei SMA, Dinner AR, Rebay I: Modeling bistable cell-fate choices in the Drosophila eye: qualitative and quantitative perspectives. Development. 2010, 137: 2265-2278. 10.1242/dev.044826PubMed CentralView ArticlePubMed
- Danielli A, Scarlato V: Regulatory circuits in Helicobacter pylori°: network motifs and regulators involved in metal-dependent responses. FEMS Microbiology Reviews. 2010, 34: 738-752.View ArticlePubMed
- Becskei A, Serrano L: Engineering stability in gene networks by autoregulation. Nature. 2000, 405: 590-593. 10.1038/35014651View ArticlePubMed
- Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks. Proc Natl Acad Sci USA. 2001, 98: 8614-8619. 10.1073/pnas.151588598PubMed CentralView ArticlePubMed
- Austin DW, Allen MS, McCollum JM, Dar RD, Wilgus JR, Sayler GS, Samatova NF, Cox CD, Simpson ML: Gene network shaping of inherent noise spectra. Nature. 2006, 439: 608-611. 10.1038/nature04194View ArticlePubMed
- Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene expression: from theories to phenotypes. Nat Rev Genet. 2005, 6: 451-464. 10.1038/nrg1615View ArticlePubMed
- McAdams HH, Arkin A: Stochastic mechanisms in gene expression. Proc Natl Acad Sci USA. 1997, 94: 814-819. 10.1073/pnas.94.3.814PubMed CentralView ArticlePubMed
- Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A: Regulation of noise in the expression of a single gene. Nat Genet. 2002, 31: 69-73. 10.1038/ng869View ArticlePubMed
- Dublanche Y, Michalodimitrakis K, Kümmerer N, Foglierini M, Serrano L: Noise in transcription negative feedback loops: simulation and experimental analysis. Mol Syst Biol. 2006, 2: 41-PubMed CentralView ArticlePubMed
- Rosenfeld N, Elowitz MB, Alon U: Negative Autoregulation Speeds the Response Times of Transcription Networks. Journal of Molecular Biology. 2002, 323: 785-793. 10.1016/S0022-2836(02)00994-4View ArticlePubMed
- Camas FM, Blázquez J, Poyatos JF: Autogenous and nonautogenous control of response in a genetic network. Proceedings of the National Academy of Sciences. 2006, 103: 12718-12723. 10.1073/pnas.0602119103.View Article
- Nevozhay D, Adams RM, Murphy KF, Josić K, Balázsi G: Negative autoregulation linearizes the dose-response and suppresses the heterogeneity of gene expression. Proceedings of the National Academy of Sciences. 2009, 106: 5123-5128. 10.1073/pnas.0809901106.View Article
- Wall ME, Hlavacek WS, Savageau MA: Design of gene circuits: lessons from bacteria. Nat Rev Genet. 2004, 5: 34-42. 10.1038/nrg1244View ArticlePubMed
- Wall ME, Hlavacek WS, Savageau MA: Design Principles for Regulator Gene Expression in a Repressible Gene Circuit. Journal of Molecular Biology. 2003, 332: 861-876. 10.1016/S0022-2836(03)00948-3View ArticlePubMed
- BÁRÁNY E: Negative feedback regulation of intracellular Ca2+ and the slope of the dose-response curve. Acta Physiologica Scandinavica. 1987, 130: 77-81. 10.1111/j.1748-1716.1987.tb08114.xView ArticlePubMed
- Verma M, Rawool S, Bhat PJ, Venkatesh KV: Biological significance of autoregulation through steady state analysis of genetic networks. Biosystems. 2006, 84: 39-48. 10.1016/j.biosystems.2005.10.001View ArticlePubMed
- Schleif R: Two positively regulated systems ara and mal. Escherichia coli and Salmonella: cellular and molecular biology. Edited by: Ingraham JL, Neidhardt FC. 1996, 1: 1300-1309. Washington: ASM Press, 2,
- Schleif R: Regulation of the L-arabinose operon of Escherichia coli. Trends Genet. 2000, 16: 559-565. 10.1016/S0168-9525(00)02153-3View ArticlePubMed
- Schleif R: AraC protein, regulation of the l-arabinose operon in Escherichia coli, and the light switch mechanism of AraC action. FEMS Microbiology Reviews. 2010, 34: 779-796.View ArticlePubMed
- Schleif R: AraC protein: A love-hate relationship. BioEssays. 2003, 25: 274-282. 10.1002/bies.10237View ArticlePubMed
- Novotny CP, Englesberg E: The L-arabinose permease system in Escherichia coli B/r. Biochim Biophys Acta. 1966, 117: 217-230.View ArticlePubMed
- Schleif R: An L-arabinose binding protein and arabinose permeation in Escherichia coli. J Mol Biol. 1969, 46: 185-196. 10.1016/0022-2836(69)90065-5View ArticlePubMed
- Doyle ME, Brown C, Hogg RW, Helling RB: Induction of the ara operon of Escherichia coli B-r. J Bacteriol. 1972, 110: 56-65.PubMed CentralPubMed
- Lobell RB, Schleif RF: DNA looping and unlooping by AraC protein. Science. 1990, 250: 528-532. 10.1126/science.2237403View ArticlePubMed
- Seabold RR, Schleif RF: Apo-AraC actively seeks to loop. Journal of Molecular Biology. 1998, 278: 529-538. 10.1006/jmbi.1998.1713View ArticlePubMed
- Kaplan S, Bren A, Zaslaver A, Dekel E, Alon U: Diverse two-dimensional input functions control bacterial sugar genes. Mol Cell. 2008, 29: 786-792. 10.1016/j.molcel.2008.01.021PubMed CentralView ArticlePubMed
- Goldbeter A, Koshland DE: Ultrasensitivity in biochemical systems controlled by covalent modification. Interplay between zero-order and multistep effects. Journal of Biological Chemistry. 1984, 259: 14441-14447.PubMed
- Goldbeter A, Koshland DE: Sensitivity amplification in biochemical systems. Q Rev Biophys. 1982, 15: 555-591. 10.1017/S0033583500003449View ArticlePubMed
- Goldbeter A, Koshland DE: An amplified sensitivity arising from covalent modification in biological systems. Proc Natl Acad Sci USA. 1981, 78: 6840-6844. 10.1073/pnas.78.11.6840PubMed CentralView ArticlePubMed
- Kaplan S, Bren A, Dekel E, Alon U: The incoherent feed-forward loop can generate non-monotonic input functions for genes. Mol Syst Biol. 2008, 4 (203):
- Zaslaver A, Bren A, Ronen M, Itzkovitz S, Kikoin I, Shavit S, Liebermeister W, Surette MG, Alon U: A comprehensive library of fluorescent transcriptional reporters for Escherichia coli. Nat Meth. 2006, 3: 623-628. 10.1038/nmeth895.View Article
- Setty Y, Mayo AE, Surette MG, Alon U: Detailed map of a cis-regulatory input function. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100: 7702-7707. 10.1073/pnas.1230759100PubMed CentralView ArticlePubMed
- Savageau MA: Design principles for elementary gene circuits: Elements, methods, and examples. Chaos. 2001, 11: 142- 10.1063/1.1349892View ArticlePubMed
- Ptashne M, Gann A: Genes & Signals. 2001, Cold Spring Harbor Laboratory Press, 1,
- Cormack BP, Valdivia RH, Falkow S: FACS-optimized mutants of the green fluorescent protein (GFP). Gene. 1996, 173: 33-38. 10.1016/0378-1119(95)00685-0View ArticlePubMed
- Lutz R, Bujard H: Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements. Nucleic Acids Res. 1997, 25: 1203-1210. 10.1093/nar/25.6.1203PubMed CentralView ArticlePubMed
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006, 2:
- Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proceedings of the National Academy of Sciences. 2000, 97: 6640-6645. 10.1073/pnas.120163297.View Article
- Zaslaver A, Mayo AE, Rosenberg R, Bashkin P, Sberro H, Tsalyuk M, Surette MG, Alon U: Just-in-time transcription program in metabolic pathways. Nat Genet. 2004, 36: 486-491. 10.1038/ng1348View ArticlePubMed
- Ronen M, Rosenberg R, Shraiman BI, Alon U: Assigning numbers to the arrows: Parameterizing a gene regulation network by using accurate expression kinetics. Proceedings of the National Academy of Sciences. 2002, 99: 10555-10560. 10.1073/pnas.152046799.View Article
- Keseler IM, Bonavides-Martinez C, Collado-Vides J, Gama-Castro S, Gunsalus RP, Johnson DA, Krummenacker M, Nolan LM, Paley S, Paulsen IT, Peralta-Gil M, Santos-Zavaleta A, Shearer AG, Karp PD: EcoCyc: A comprehensive view of Escherichia coli biology. Nucleic Acids Research. 2009, 37: D464-D470. 10.1093/nar/gkn751PubMed CentralView ArticlePubMed
- Ecocyc-E.coli K12 MG1655 database. http://www.ecocyc.com/
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.