The lac operon (lactose operon) is an operon required for the transport and metabolism of lactose in Escherichia coli and many other enteric bacteria. Although glucose is the preferred carbon source for most bacteria, the lac operon allows for the effective digestion of lactose when glucose is not available through the activity of beta-galactosidase. Gene regulation of the lac operon was the first genetic regulatory mechanism to be understood clearly, so it has become a foremost example of prokaryoticgene regulation. It is often discussed in introductory molecular and cellular biology classes for this reason. This lactose metabolism system was used by Jacob and Monod to determine how a cell knows which enzyme to synthesize. Their work on the lac operon won them the Nobel Prize in Physiology in 1965.
Bacterial operons are polycistronic transcripts that are able to produce multiple proteins from one mRNA transcript. In this case, when lactose is required as a sugar source for the bacterium, the three genes of the lac operon can be expressed and their subsequent proteins translated: lacZ, lacY, and lacA. The gene product of lacZ is β-galactosidase which cleaves lactose, a disaccharide, into glucose and galactose. lacY encodes Beta-galactoside permease, a protein which becomes embedded in the cytoplasmic membrane to enable transport of lactose into the cell. Finally, lacA encodes β-galactoside transacetylase.
It would be wasteful to produce the enzymes when there is no lactose available or if there is a more preferable energy source available, such as glucose. The lac operon uses a two-part control mechanism to ensure that the cell expends energy producing the enzymes encoded by the lac operon only when necessary. In the absence of lactose, the lac repressor, lacI, halts production of the enzymes encoded by the lac operon. The lac repressor is always expressed unless a co-inducer binds to it. In other words, it is transcribed only in the presence of small molecule co-inducer. In the presence of glucose, the catabolite activator protein (CAP), required for production of the enzymes, remains inactive, and EIIAGlc shuts down lactose permease to prevent transport of lactose into the cell. This dual control mechanism causes the sequential utilization of glucose and lactose in two distinct growth phases, known as diauxie.
- The lac operon consists of three structural genes, and a promoter, a terminator, regulator, and an operator. The three structural genes are: lacZ, lacY, and lacA.
Only lacZ and lacY appear to be necessary for lactose catabolism.
Three-letter abbreviations are used to describe phenotypes in bacteria including E. coli.
- Lac (the ability to use lactose),
- His (the ability to synthesize the amino acid histidine)
- Mot (swimming motility)
- SmR (resistance to the antibiotic streptomycin)
In the case of Lac, wild type cells are Lac+ and are able to use lactose as a carbon and energy source, while Lac− mutant derivatives cannot use lactose. The same three letters are typically used (lower-case, italicized) to label the genes involved in a particular phenotype, where each different gene is additionally distinguished by an extra letter. The lac genes encoding enzymes are lacZ, lacY, and lacA. The fourth lac gene is lacI, encoding the lactose repressor—"I" stands for inducibility.
One may distinguish between structural genes encoding enzymes, and regulatory genes encoding proteins that affect gene expression. Current usage expands the phenotypic nomenclature to apply to proteins: thus, LacZ is the protein product of the lacZ gene, β-galactosidase. Various short sequences that are not genes also affect gene expression, including the lac promoter, lac p, and the lac operator, lac o. Although it is not strictly standard usage, mutations affecting lac o are referred to as lac oc, for historical reasons.
Specific control of the lac genes depends on the availability of the substrate lactose to the bacterium. The proteins are not produced by the bacterium when lactose is unavailable as a carbon source. The lac genes are organized into an operon; that is, they are oriented in the same direction immediately adjacent on the chromosome and are co-transcribed into a single polycistronic mRNA molecule. Transcription of all genes starts with the binding of the enzyme RNA polymerase (RNAP), a DNA-binding protein, which binds to a specific DNA binding site, the promoter, immediately upstream of the genes. Binding of RNA polymerase to the promoter is aided by the cAMP-bound catabolite activator protein (CAP, also known as the cAMP receptor protein). However, the lacI gene (regulatory gene for lac operon) produces a protein that blocks RNAP from binding to the promoter of the operon. This protein can only be removed when allolactose binds to it, and inactivates it. The protein that is formed by the lacI gene is known as the lac repressor. The type of regulation that the lac operon undergoes is referred to as negative inducible, meaning that the gene is turned off by the regulatory factor (lac repressor) unless some molecule (lactose) is added. Because of the presence of the lac repressor protein, genetic engineers who replace the lacZ gene with another gene will have to grow the experimental bacteria on agar with lactose available on it. If they do not, the gene they are trying to express will not be expressed as the repressor protein is still blocking RNAP from binding to the promoter and transcribing the gene. Once the repressor is removed, RNAP then proceeds to transcribe all three genes (lacZYA) into mRNA. Each of the three genes on the mRNA strand has its own Shine-Dalgarno sequence, so the genes are independently translated. The DNA sequence of the E. coli lac operon, the lacZYAmRNA, and the lacI genes are available from GenBank(view).
The first control mechanism is the regulatory response to lactose, which uses an intracellular regulatory protein called the lactose repressor to hinder production of β-galactosidase in the absence of lactose. The lacI gene coding for the repressor lies nearby the lac operon and is always expressed (constitutive). If lactose is missing from the growth medium, the repressor binds very tightly to a short DNA sequence just downstream of the promoter near the beginning of lacZ called the lac operator. The repressor binding to the operator interferes with binding of RNAP to the promoter, and therefore mRNA encoding LacZ and LacY is only made at very low levels. When cells are grown in the presence of lactose, however, a lactose metabolite called allolactose, made from lactose by the product of the lacZ gene, binds to the repressor, causing an allosteric shift. Thus altered, the repressor is unable to bind to the operator, allowing RNAP to transcribe the lac genes and thereby leading to higher levels of the encoded proteins.
The second control mechanism is a response to glucose, which uses the catabolite activator protein (CAP) homodimer to greatly increase production of β-galactosidase in the absence of glucose. Cyclic adenosine monophosphate (cAMP) is a signal molecule whose prevalence is inversely proportional to that of glucose. It binds to the CAP, which in turn allows the CAP to bind to the CAP binding site (a 16 bp DNA sequence upstream of the promoter on the left in the diagram below, about 60 bp upstream of the transcription start site), which assists the RNAP in binding to the DNA. In the absence of glucose, the cAMP concentration is high and binding of CAP-cAMP to the DNA significantly increases the production of β-galactosidase, enabling the cell to hydrolyse lactose and release galactose and glucose.
More recently inducer exclusion was shown to block expression of the lac operon when glucose is present. Glucose is transported into the cell by the PEP-dependent phosphotransferase system. The phosphate group of phosphoenolpyruvate is transferred via a phosphorylation cascade consisting of the general PTS (phosphotransferase system) proteins HPr and EIA and the glucose-specific PTS proteins EIIAGlc and EIIBGlc, the cytoplasmic domain of the EII glucose transporter. Transport of glucose is accompanied by its phosphorylation by EIIBGlc, draining the phosphate group from the other PTS proteins, including EIIAGlc. The unphosphorylated form of EIIAGlc binds to the lac permease and prevents it from bringing lactose into the cell. Therefore, if both glucose and lactose are present, the transport of glucose blocks the transport of the inducer of the lac operon.
The lac repressor is a four-part protein, a tetramer, with identical subunits. Each subunit contains a helix-turn-helix (HTH) motif capable of binding to DNA. The operator site where repressor binds is a DNA sequence with inverted repeat symmetry. The two DNA half-sites of the operator together bind to two of the subunits of the repressor. Although the other two subunits of repressor are not doing anything in this model, this property was not understood for many years.
Eventually it was discovered that two additional operators are involved in lac regulation. One (O3) lies about -90 bp upstream of O1 in the end of the lacI gene, and the other (O2) is about +410 bp downstream of O1 in the early part of lacZ. These two sites were not found in the early work because they have redundant functions and individual mutations do not affect repression very much. Single mutations to either O2 or O3 have only 2 to 3-fold effects. However, their importance is demonstrated by the fact that a double mutant defective in both O2 and O3 is dramatically de-repressed (by about 70-fold).
In the current model, lac repressor is bound simultaneously to both the main operator O1 and to either O2 or O3. The intervening DNA loops out from the complex. The redundant nature of the two minor operators suggests that it is not a specific looped complex that is important. One idea is that the system works through tethering; if bound repressor releases from O1 momentarily, binding to a minor operator keeps it in the vicinity, so that it may rebind quickly. This would increase the affinity of repressor for O1.
Mechanism of induction
The repressor is an allosteric protein, i.e. it can assume either one of two slightly different shapes, which are in equilibrium with each other. In one form the repressor will bind to the operator DNA with high specificity, and in the other form it has lost its specificity. According to the classical model of induction, binding of the inducer, either allolactose or IPTG, to the repressor affects the distribution of repressor between the two shapes. Thus, repressor with inducer bound is stabilized in the non-DNA-binding conformation. However, this simple model cannot be the whole story, because repressor is bound quite stably to DNA, yet it is released rapidly by addition of inducer. Therefore, it seems clear that inducer can also bind to the repressor when the repressor is already bound to DNA. It is still not entirely known what the exact mechanism of binding is.
Role of non-specific binding
Non-specific binding of the repressor to DNA plays a crucial role in the repression and induction of the Lac-operon. The specific binding site for the Lac-repressor protein is the operator. The non-specific interaction is mediated mainly by charge-charge interactions while binding to the operator is reinforced by hydrophobic interactions. Additionally, there is an abundance of non-specific DNA sequences to which the repressor can bind. Essentially, any sequence that is not the operator, is considered non-specific. Studies have shown, that without the presence of non-specific binding, induction (or unrepression) of the Lac-operon could not occur even with saturated levels of inducer. It had been demonstrated that, without non-specific binding, the basal level of induction is ten thousand times smaller than observed normally. This is because the non-specific DNA acts as sort of a "sink" for the repressor proteins, distracting them from the operator. The non-specific sequences decrease the amount of available repressor in the cell. This in turn reduces the amount of inducer required to unrepress the system.
A number of lactose derivatives or analogs have been described that are useful for work with the lac operon. These compounds are mainly substituted galactosides, where the glucose moiety of lactose is replaced by another chemical group.
- Isopropyl-β-D-thiogalactopyranoside (IPTG) is frequently used as an inducer of the lac operon for physiological work. IPTG binds to repressor and inactivates it, but is not a substrate for β-galactosidase. One advantage of IPTG for in vivo studies is that since it cannot be metabolized by E. coli its concentration remains constant and the rate of expression of lac p/o-controlled genes, is not a variable in the experiment. IPTG intake is dependent on the action of lactose permease in P. fluorescens, but not in E. coli.
- Phenyl-β-D-galactose (phenyl-Gal) is a substrate for β-galactosidase, but does not inactivate repressor and so is not an inducer. Since wild type cells produce very little β-galactosidase, they cannot grow on phenyl-Gal as a carbon and energy source. Mutants lacking repressor are able to grow on phenyl-Gal. Thus, minimal medium containing only phenyl-Gal as a source of carbon and energy is selective for repressor mutants and operator mutants. If 108 cells of a wild type strain are plated on agar plates containing phenyl-Gal, the rare colonies which grow are mainly spontaneous mutants affecting the repressor. The relative distribution of repressor and operator mutants is affected by the target size. Since the lacI gene encoding repressor is about 50 times larger than the operator, repressor mutants predominate in the selection.
- Thiomethyl galactosidase [TMG] is another lactose analog. These inhibit the lacI repressor. At low inducer concentrations, both TMG and IPTG can enter the cell through the lactose permease. However at high inducer concentrations, both analogs can enter the cell independently. TMG can can reduce growth rates at high extracellular concentrations.
- Other compounds serve as colorful indicators of β-galactosidase activity.
- ONPG is cleaved to produce the intensely yellow compound, orthonitrophenol and galactose, and is commonly used as a substrate for assay of β-galactosidase in vitro.
- Colonies that produce β-galactosidase are turned blue by X-gal (5-bromo-4-chloro-3-indolyl-β-D-galactoside) which is an artificial substrate for B-galactosidase whose cleavage results in galactose and 4-Cl,3-Br indigo thus producing a deep blue color.
- Allolactose is an isomer of lactose and is the inducer of the lac operon. Lactose is galactose-(β1->4)-glucose, whereas allolactose is galactose-(β1->6)-glucose. Lactose is converted to allolactose by β-galactosidase in an alternative reaction to the hydrolytic one. A physiological experiment which demonstrates the role of LacZ in production of the "true" inducer in E. coli cells is the observation that a null mutant of lacZ can still produce LacY permease when grown with IPTG but not when grown with lactose. The explanation is that processing of lactose to allolactose (catalyzed by β-galactosidase) is needed to produce the inducer inside the cell.
Development of the classic model
The experimental microorganism used by François Jacob and Jacques Monod was the common laboratory bacterium, E. coli, but many of the basic regulatory concepts that were discovered by Jacob and Monod are fundamental to cellular regulation in all organisms. The key idea is that proteins are not synthesized when they are not needed--- E. coli conserves cellular resources and energy by not making the three Lac proteins when there is no need to metabolize lactose, such as when other sugars like glucose are available. The following section discusses how E. coli controls certain genes in response to metabolic needs.
During World War II, Monod was testing the effects of combinations of sugars as nutrient sources for E. coli and B. subtilis. Monod was following up on similar studies that had been conducted by other scientists with bacteria and yeast. He found that bacteria grown with two different sugars often displayed two phases of growth. For example, if glucose and lactose were both provided, glucose was metabolized first (growth phase I, see Figure 2) and then lactose (growth phase II). Lactose was not metabolized during the first part of the diauxic growth curve because β-galactosidase was not made when both glucose and lactose were present in the medium. Monod named this phenomenon diauxie.
Monod then focused his attention on the induction of β-galactosidase formation that occurred when lactose was the sole sugar in the culture medium.
Classification of regulatory mutants
A conceptual breakthrough of Jacob and Monod was to recognize the distinction between regulatory substances and sites where they act to change gene expression. A former soldier, Jacob used the analogy of a bomber that would release its lethal cargo upon receipt of a special radio transmission or signal. A working system requires both a ground transmitter and a receiver in the airplane. Now, suppose that the usual transmitter is broken. This system can be made to work by introduction of a second, functional transmitter. In contrast, he said, consider a bomber with a defective receiver. The behavior of this bomber cannot be changed by introduction of a second, functional aeroplane.
To analyze regulatory mutants of the lac operon, Jacob developed a system by which a second copy of the lac genes (lacI with its promoter, and lacZYA with promoter and operator) could be introduced into a single cell. A culture of such bacteria, which are diploid for the lac genes but otherwise normal, is then tested for the regulatory phenotype. In particular, it is determined whether LacZ and LacY are made even in the absence of IPTG (due to the lactose repressor produced by the mutant gene being non-functional). This experiment, in which genes or gene clusters are tested pairwise, is called a complementation test.
This test is illustrated in the figure (lacA is omitted for simplicity). First, certain haploid states are shown (i.e. the cell carries only a single copy of the lac genes). Panel (a) shows repression, (b) shows induction by IPTG, and (c) and (d) show the effect of a mutation to the lacI gene or to the operator, respectively. In panel (e) the complementation test for repressor is shown. If one copy of the lac genes carries a mutation in lacI, but the second copy is wild type for lacI, the resulting phenotype is normal---but lacZ is expressed when exposed to inducer IPTG. Mutations affecting repressor are said to be recessive to wild type (and that wild type is dominant), and this is explained by the fact that repressor is a small protein which can diffuse in the cell. The copy of the lac operon adjacent to the defective lacI gene is effectively shut off by protein produced from the second copy of lacI.
If the same experiment is carried out using an operator mutation, a different result is obtained (panel (f)). The phenotype of a cell carrying one mutant and one wild type operator site is that LacZ and LacY are produced even in the absence of the inducer IPTG; because the damaged operator site, does not permit binding of the repressor to inhibit transcription of the structural genes. The operator mutation is dominant. When the operator site where repressor must bind is damaged by mutation, the presence of a second functional site in the same cell makes no difference to expression of genes controlled by the mutant site.
A more sophisticated version of this experiment uses marked operons to distinguish between the two copies of the lac genes and show that the unregulated structural gene(s) is(are) the one(s) next to the mutant operator (panel (g). For example, suppose that one copy is marked by a mutation inactivating lacZ so that it can only produce the LacY protein, while the second copy carries a mutation affecting lacY and can only produce LacZ. In this version, only the copy of the lac operon that is adjacent to the mutant operator is expressed without IPTG. We say that the operator mutation is cis-dominant, it is dominant to wild type but affects only the copy of the operon which is immediately adjacent to it.
This explanation is misleading in an important sense, because it proceeds from a description of the experiment and then explains the results in terms of a model. But in fact, it is often true that the model comes first, and an experiment is fashioned specifically to test the model. Jacob and Monod first imagined that there must be a site in DNA with the properties of the operator, and then designed their complementation tests to show this.
The dominance of operator mutants also suggests a procedure to select them specifically. If regulatory mutants are selected from a culture of wild type using phenyl-Gal, as described above, operator mutations are rare compared to repressor mutants because the target-size is so small. But if instead we start with a strain which carries two copies of the whole lac region (that is diploid for lac), the repressor mutations (which still occur) are not recovered because complementation by the second, wild type lacI gene confers a wild type phenotype. In contrast, mutation of one copy of the operator confers a mutant phenotype because it is dominant to the second, wild type copy.
Regulation by cyclic AMP
Explanation of diauxie depended on the characterization of additional mutations affecting the lac genes other than those explained by the classical model. Two other genes, cya and crp, subsequently were identified that mapped far from lac, and that, when mutated, result in a decreased level of expression in the presence of IPTG and even in strains of the bacterium lacking the repressor or operator. The discovery of cAMP in E. coli led to the demonstration that mutants defective the cya gene but not the crp gene could be restored to full activity by the addition of cAMP to the medium.
The cya gene encodes adenylate cyclase, which produces cAMP. In a cya mutant, the absence of cAMP makes the expression of the lacZYA genes more than ten times lower than normal. Addition of cAMP corrects the low Lac expression characteristic of cya mutants. The second gene, crp, encodes a protein called catabolite activator protein (CAP) or cAMP receptor protein (CRP).
However the lactose metabolism enzymes are made in small quantities in the presence of both glucose and lactose (sometimes called leaky expression) due to the fact that the LacI repressor rapidly associates/dissociates from the DNA rather than tightly binding to it, which can allow time for RNAP to bind and transcribe mRNAs of lacZYA. Leaky expression is necessary in order to allow for metabolism of some lactose after the glucose source is expended, but before lac expression is fully activated.
- When lactose is absent then there is very little Lac enzyme production (the operator has Lac repressor bound to it).
- When lactose is present but a preferred carbon source (like glucose) is also present then a small amount of enzyme is produced (Lac repressor is not bound to the operator).
- When glucose is absent, CAP-cAMP binds to a specific DNA site upstream of the promoter and makes a direct protein-protein interaction with RNAP that facilitates the binding of RNAP to the promoter.
The delay between growth phases reflects the time needed to produce sufficient quantities of lactose-metabolizing enzymes. First, the CAP regulatory protein has to assemble on the lac promoter, resulting in an increase in the production of lacmRNA. More available copies of the lac mRNA results in the production (see translation) of significantly more copies of LacZ (β-galactosidase, for lactose metabolism) and LacY (lactose permease to transport lactose into the cell). After a delay needed to increase the level of the lactose metabolizing enzymes, the bacteria enter into a new rapid phase of cell growth.
Two puzzles of catabolite repression relate to how cAMP levels are coupled to the presence of glucose, and secondly, why the cells should even bother. After lactose is cleaved it actually forms glucose and galactose (easily converted to glucose). In metabolic terms, lactose is just as good a carbon and energy source as glucose. The cAMP level is related not to intracellular glucose concentration but to the rate of glucose transport, which influences the activity of adenylate cyclase. (In addition, glucose transport also leads to direct inhibition of the lactose permease.) As to why E. coli works this way, one can only speculate. All enteric bacteria ferment glucose, which suggests they encounter it frequently. It is possible that a small difference in efficiency of transport or metabolism of glucose v. lactose makes it advantageous for cells to regulate the lac operon in this way.
Use in molecular biology
The lac gene and its derivatives are amenable to use as a reporter gene in a number of bacterial-based selection techniques such as two hybrid analysis, in which the successful binding of a transcriptional activator to a specific promoter sequence must be determined. In LB plates containing X-gal, the colour change from white colonies to a shade of blue corresponds to about 20–100 β-galactosidase units, while tetrazolium lactose and MacConkey lactose media have a range of 100–1000 units, being most sensitive in the high and low parts of this range respectively. Since MacConkey lactose and tetrazolium lactose media both rely on the products of lactose breakdown, they require the presence of both lacZ and lacY genes. The many lac fusion techniques which include only the lacZ gene are thus suited to X-gal plates or ONPG liquid broths.
- ^Griffiths, Anthony J.F.; Wessler, Susan R.; Carroll, Sean B.; Doebley, John (2015). An Introduction to Genetic Analysis (11 ed.). Freeman, W.H. & Company. pp. 400–412. ISBN 9781464109485.
- ^Griffiths, Anthony J.F.; Wessler, Susan R.; Carroll, Sean B.; Doebley, John (2015). An Introduction to Genetic Analysis (11 ed.). Freeman, W.H. & Company. pp. 400–412. ISBN 9781464109485.
- ^"Prokaryotic Gene Expression". www.ndsu.edu.
- ^"Prokaryotic Gene Expression". www.ndsu.edu. Retrieved 2017-05-19.
- ^Busby S., Ebright RH. (2001). "Transcription activation by catabolite activator protein (CAP)". J. Mol. Biol. 293: 199–213. doi:10.1006/jmbi.1999.3161. PMID 10550204.
- ^Kennell, David; Riezman, Howard (July 1977). "Transcription and translation initiation frequencies of the Escherichia coli lac operon". Journal of Molecular Biology. 114 (1): 1–21. doi:10.1016/0022-2836(77)90279-0.
- ^Malan, T. Philip; Kolb, Annie; Buc, Henri; McClure, William (December 1984). "Mechanism of CRP-cAMP Activation of lac Operon Transcription Initiation Activation of the P1 Promoter". J. Mol. Biol. 180: 881–909. doi:10.1016/0022-2836(84)90262-6. PMID 6098691.
- ^Görke B, Stülke J (August 2008). "Carbon catabolite repression in bacteria: many ways to make the most out of nutrients". Nature Reviews. Microbiology. 6 (8): 613–24. doi:10.1038/nrmicro1932. PMID 18628769.
- ^Oehler, S.; Eismann, E. R.; Krämer, H.; Müller-Hill, B. (1990). "The three operators of the lac operon cooperate in repression". The EMBO Journal. 9 (4): 973–979. PMC 551766. PMID 2182324.
- ^Griffiths, Anthony JF; Gelbart, William M.; Miller, Jeffrey H.; Lewontin, Richard C. (1999). "Regulation of the Lactose System".
- ^von Hippel, P.H.; Revzin, A.; Gross, C.A.; Wang, A.C. (December 1974). "Non-specific DNA binding of genome regulating proteins as a biological control mechanism: I. The lac operon: equilibrium aspects". PNAS. 71 (12): 4808–12. doi:10.1073/pnas.71.12.4808. PMC 433986. PMID 4612528.
- ^Hansen LH, Knudsen S, Sørensen SJ (June 1998). "The effect of the lacY gene on the induction of IPTG inducible promoters, studied in Escherichia coli and Pseudomonas fluorescens". Curr. Microbiol. 36 (6): 341–7. doi:10.1007/s002849900320. PMID 9608745. [permanent dead link]
- ^Marbach A, Bettenbrock K (Jan 2012). "lac operon induction in Escherichia coli: Systematic comparison of IPTG and TMG induction and influence of the transacetylase LacA". Science Direct. 157 (1): 82–88. doi:10.1016/j.jbiotec.2011.10.009. PMID 22079752.
- ^"ONPG (β-Galactosidase) test". September 2000.
- ^ abcdJoung J, Ramm E, Pabo C (2000). "A bacterial two-hybrid selection system for studying protein–DNA and protein–protein interactions". Proc Natl Acad Sci USA. 97 (13): 7382–7. doi:10.1073/pnas.110149297. PMC 16554. PMID 10852947.
- ^"Milestone 2 - A visionary pair : Nature Milestones in gene expression". www.nature.com. Retrieved 2015-12-27.
- ^Muller-Hill, Benno (1996). The lac Operon, a Short History of a Genetic Paradigm. Berlin: Walter de Gruyter. pp. 7–10. ISBN 3-11-014830-7.
- ^McKnight, Steven L. (1992). Transcriptional Regulation. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press. pp. 3–24. ISBN 0-87969-410-6.
- ^Jacob F.; Monod J (June 1961). "Genetic regulatory mechanisms in the synthesis of proteins". J Mol Biol. 3 (3): 318–56. doi:10.1016/S0022-2836(61)80072-7. PMID 13718526.
- ^Montminy, M. (1997). "Transcriptional regulation by cyclic AMP". Annual Review of Biochemistry. 66: 807–822. doi:10.1146/annurev.biochem.66.1.807. ISSN 0066-4154. PMID 9242925.
- ^Botsford, J L; Harman, J G (March 1992). "Cyclic AMP in prokaryotes". Microbiological Reviews. 56 (1): 100–122. ISSN 0146-0749. PMC 372856. PMID 1315922.
- ^Vazquez A, Beg QK, Demenezes MA, et al. (2008). "Impact of the solvent capacity constraint on E. coli metabolism". BMC Syst Biol. 2: 7. doi:10.1186/1752-0509-2-7. PMC 2270259. PMID 18215292.
- ^"Induction of the lac operon in E. coli"(PDF). SAPS. Retrieved 29 June 2016.
1: RNA Polymerase, 2: Repressor, 3: Promoter, 4: Operator, 5: Lactose, 6: lacZ, 7: lacY, 8: lacA.
Quantitative approaches to the study of bistability in the lac operon of Escherichia coli
Moisés Santillán1,2,* and Michael C. Mackey2
1Centro de Investigación y Estudios Avanzados del IPN, Unidad Monterrey, Parque de Investigación e Innovación Tecnológica, 66600 Apodaca, NL, Mexico
2Centre for Nonlinear Dynamics in Physiology and Medicine, McGill University, 3655 Promenade Sir William Osler, Montreal, Que., Canada H3G 1Y6
*Author and address for correspondence: Centro de Investigación y Estudios Avanzados del IPN, Unidad Monterrey, Parque de Investigación e Innovación Tecnológica, 66600 Apodaca NL, Mexico (xm.vatsevnic@nallitnasm)
Author information ►Article notes ►Copyright and License information ►
Received 2008 Feb 27; Accepted 2008 Mar 27.
Copyright © 2008 The Royal Society
This article has been cited by other articles in PMC.
In this paper, the history and importance of the lac operon in the development of molecular and systems biology are briefly reviewed. We start by presenting a description of the regulatory mechanisms in this operon, taking into account the most recent discoveries. Then we offer a survey of the history of the lac operon, including the discovery of its main elements and the subsequent influence on the development of molecular and systems biology. Next the bistable behaviour of the operon is discussed, both with respect to its discovery and its molecular origin. A review of the literature in which this bistable phenomenon has been studied from a mathematical modelling viewpoint is then given. We conclude with some brief remarks.
Keywords: lac operon, multistability, gene regulatory networks, systems biology, computational biology, mathematical model
Glucose is the favourite carbon and energy source for Escherichia coli, as well as for many other organisms. Although this bacterium can also feed on other sugars, it only does so when glucose is absent. Thus, if a bacterial culture grows in a medium containing a mixture of glucose and another sugar (such as lactose), it will exclusively feed on the former until it is exhausted, before switching on to the second one. A consequence of this behaviour is that the bacterial growth curve shows two distinctive phases, as can be seen in figure 1. This phenomenon was originally studied by Monod (1941), who described it as diauxic growth. It is worth mentioning at this point that diauxic growth only occurs in batch cultures, and simultaneous usage of sugars is often observed in continuous cultures (Lendenmann et al. 1996).
Typical diauxic growth curve. Note the existence of two different exponential growth phases, separated by a short interval in which the culture does not grow. The first (second) phase corresponds to the bacterial culture feeding on glucose (lactose),...
Molecular level understanding of how an organism sequentially uses different metabolites has been attracting tremendous interest for the past fifty years. Jacob & Monod (1961a,b) tackled this problem and conceptually outlined how bacterial cultures could switch from one mode of growth to another so rapidly and completely. In the process of doing so, they introduced the operon concept (Jacob et al. 1960), which has become a paradigmatic example of gene regulation.
According to the construction of Jacob & Monod, an operon consists of a set of structural genes that are regulated together, depending on the cell metabolic requirements. These structural genes code for a group of enzymes or proteins that are responsible for a specific task or metabolic process, and their regulation is achieved via one or more common regulatory mechanisms. Repression was the first regulatory mechanism to be discovered by them. In it, a repressor molecule binds a specific DNA site (which they termed the operator) located upstream from the structural genes, and inhibits their transcription. The regulation of the structural genes' expression is achieved by varying the number of active repressor molecules. Although Jacob and Monod originally thought of the repressors as RNA molecules, they are now known to be proteins.
The lactose (or simply lac) operon is composed of three structural genes: lacZ, lacY and lacA. These genes, respectively, code for β-galactosidase, lac permease and a transacetylase. β-Galactosidase acts to cleave lactose into galactose and glucose, which is the first step in lactose metabolism; lac permease is a transmembrane protein, which is necessary for lactose uptake; transacetylase transfers an acetyl group from coenzyme A (CoA) to the hydroxyl group of the galactosides. Of these proteins, only β-galactosidase and lac permease play an active role in the regulation of the lac operon. The regulatory gene lacI (in a different operon) codes for the lac repressor which, when active, is capable of inhibiting the transcription of the structural genes by binding an operator. The lac repressor is inactivated when it is bound by allolactose, a by-product of lactose metabolism. Finally, the lac operon genes are also upregulated by an activator that increases the affinity of the mRNA polymerase for the lac promoter, and whose production is controlled by the concentration of extracellular glucose.
The lac operon has been pivotal in the development of molecular biology and is currently having an important impact on the development of systems biology. In §2 we present a more detailed description of the lac operon regulatory mechanisms, in §3 we review the most significant aspects of the development of molecular biology influenced by the lac operon and in §4 we give a comprehensive review of the quantitative studies of the lac operon bistable behaviour.
2. The lac operon regulatory mechanisms in detail
In §1, the lac operon control system was briefly reviewed. Although that brief review gives a good idea of how this system functions, in reality it is far more complex. In the following paragraphs, a more detailed description of all the regulatory mechanisms in the lac operon, including the most recent discoveries, will be given.
The lac operon regulatory elements (pictured in figure 2a) are distributed along the DNA chain as follows (Reznikoff 1992; Müller-Hill 1998): the lac promoter is located between bp −36 (bp stands for base pair, and positions are referred relative to the starting point of gene lacZ, bp +1) and bp −7. Operator O1 is 21 bp long and is centred around bp +11. There are two additional operators, denoted O2 and O3, which are, respectively, located at 401 bp downstream and 92 bp upstream from O1. Finally, the activator (CAP)-binding site spans from bp −72 to bp −50.
(a) Schematic of the regulatory elements located in lac operon DNA. P denotes the promoter, O1, O2 and O3 correspond to the three operators (repressor-binding sites), and C is the binding site for the cAMP–CRP complex. The different ways in which...
The lac repressor is a homotetramer (consisting of two functional homodimers) of lacI polypeptides (Lewis 2005; Wilson et al. 2007). Each functional dimer can bind operators O1, O2 and O3. Furthermore, DNA can also fold in such a way that a single repressor binds two operators simultaneously, one per dimer. Each monomer in the lac repressor can be bound by an allolactose molecule, inhibiting the capability of the corresponding dimer to bind an operator. This means that free repressors can bind one operator (figure 2b) or two of them simultaneously (figure 2c), repressors with three free monomers can bind one but not two operators (figure 2d), repressors with two free monomers can bind one operator, if the bound monomers belong to the same dimer (figure 2e), or none at all, and that repressors with only one free monomer are unable to bind any operator, as are repressors with all four monomers bound by allolactose (Narang 2007).
Deletion experiments have shown that a repressor bound to O1 inhibits transcription initiation, while a repressor bound to either O2 or O3 has almost no effect on the expression of the lac operon structural genes. Nevertheless, O2 and O3 do have an indirect effect because the complex formed by a single repressor simultaneously bound to O1 and either O2 or O3 is far more stable than that of a repressor bound only to O1. The consequence of this is that by interacting with the lac repressor operator O1 is only capable of decreasing the expression of the operon genes 18 times; when it cooperates with O2, the repression level can be as high as 700-fold; when O1 and O3 act together, they can reduce the operon activity up to 440 times; when all three operators are present, the repression intensity can be as high as 1300-fold (Oehler et al. 1990).
The intracellular production of cyclic AMP (cAMP) decreases as the concentration of extracellular glucose increases. cAMP further binds a specific receptor molecule (CRP) to form the so-called CAP complex. Finally, CAP binds a specific DNA site (denoted here as C) upstream from the lac promoter, and by doing so it increases the affinity of the mRNA polymerase for this promoter (Reznikoff 1992). This regulatory mechanism is known as catabolite repression.
A novel source of cooperativity has been recently discovered in the lac operon: when a CAP complex is bound to site C, it bends DNA locally and increases the probability of the complex in which a repressor simultaneously binds operators O1 and O3 (Kuhlman et al. 2007).
The last regulatory mechanism in the lac operon is the so-called inducer exclusion. In it, external glucose decreases the efficiency of lac permease to transport lactose (Reznikoff 1992), and by doing so negatively affects the induction of the operon genes.
These regulatory mechanisms reviewed above are summarized in figure 3. As we have seen, the activity of the lac operon is regulated by extracellular glucose and lactose. While extracellular glucose decreases the operon activity via catabolite repression and inducer exclusion, extracellular lactose increases the operon expression level by deactivating the repressor. Another fact worth noticing is the existence of a positive feedback loop: as more molecules of lac permease and β-galactosidase are produced, there is an elevated lactose uptake flux and an increased lactose metabolism rate; this further increases the production of allolactose and, as a consequence, diminishes the amount of active repressor. This, in turn, increases the operon activity, and thus more lac permease and β-galactosidase molecules are produced.
Schematic of the lac operon regulatory mechanisms. This operon consists of genes lacZ, lacY and lacA. Protein LacY is a permease that transports external lactose into the cell. Protein LacZ polymerizes into a homotetramer named β-galactosidase....
The reader interested in the details of the lac operon regulatory mechanisms is referred to the excellent review by Beckwith (1987) and the references there. A good description of the operon regulatory elements and their location on the DNA chain can be found in Reznikoff (1992). The most recent discoveries regarding the cooperativity between CAP-binding site and operator O3 are reported in Kuhlman et al. (2007).
3. Importance of the lac operon
The operon model, developed by Jacob and Monod in their attempt to explain diauxic growth, depicted how genetic mechanisms can control metabolic events in response to environmental stimuli via the coordinated transcription of a set of genes with related function. It literally became a paradigm for gene regulation in prokaryotes, where many more operons have been discovered, and has also influenced the understanding of gene regulation in eukaryotes. Furthermore, not only the lac operon as a whole, but also its individual components, such as the lac repressor, the three known lac operators, the enzyme β-galactosidase and the protein lac permease, have influenced the development of molecular biology themselves.
3.1 The lac repressor
A year after Jacob and Monod received the Nobel Prize for their contributions to gene regulation, Müller-Hill and Gilbert isolated the lac repressor. This is a protein of 360 amino acids which associates with a homotetramer with 154 520 Da molecular mass (Lewis 2005).
The molecular mechanism of repressing the lac operon requires the repressor to be capable of binding both operator DNA and allolactose (or similar inducers). The possibility of competitive binding by these ligands was eliminated by the demonstration that protease digestion selectively cleaves the repressor into two fragments: a tetrameric ‘core’ (residues 60–360 of the monomer) that retains inducer-binding properties, and a monomeric N-terminal headpiece (amino acids 1–59) capable of binding DNA. To explain the repressor inactivation by allolactose, Monod, Changeux and Jacob proposed that the repressor undergoes a conformational transition in response to bound ligands, and that this alters its ability to bind DNA. They named this phenomenon allostery (Monod et al. 1963).
In the early 1970s, several hundred milligrams of the repressor were purified and used for crystallization. Yet, its three-dimensional architecture remained elusive until the early 1990s (Lewis 2005). The three-dimensional structure of the lac repressor provided insight into how the repressor may function, as well as the three-dimensional framework for interpreting a huge amount of biochemical and genetic information. Most importantly, when the biochemical and genetic data were viewed in the context of the structure, a detailed molecular model could be constructed to provide a physical basis for the allosteric response, as well as a more detailed understanding of the genetic switch in the lac operon.
The allosteric response discovered in the lactose repressor opened a whole new area of research. Allostery has been found in many other proteins and has also been extended to a variety of cellular signalling pathways in all organisms. Notwithstanding, the transcendence of the lactose repressor does not end there. For instance, its monomer has recently been used as a model system for experimental and theoretical explorations of protein-folding mechanisms (Wilson et al. 2005).
Those interested in knowing more about the lactose repressor can refer to the review papers by Lewis (2005) and Wilson et al. (2007).
3.2 The three lac operators
The primary operator site (O1) for the lac operon was sequenced by Gilbert & Maxam (1973) nearly a decade after Jacob and Monod had published their model. In addition to O1, two auxiliary operators (O2 and O3) were identified with sequences similar to those of the primary operator (Reznikoff et al. 1974). We now know that tetrameric lac repressor is ideally suited to bind two operators simultaneously, creating the so-called ‘repression loops’ (Reznikoff 1992). DNA looping enhances the repressor affinity for multi-operator sequences, and supercoiling these DNAs yields complexes with remarkable stability. A number of synthetic operator variants have also been constructed and have proved very useful for understanding the molecular mechanisms of repression (Wilson et al. 2007).
It was thought for decades that all the signals that control the initiation of bacterial gene transcription are clustered at the 5′ ends of operons, as proposed originally in the models of Jacob, Monod and co-workers. This aspect of their pioneering work is now known to be an oversimplification, as initiation control signals have since been found within, downstream and upstream of the genes regulated by them. Although the phenomenon is not as widespread as in higher cells, its study in bacteria can be, in particular, illuminating. Together with phage lambda switch, the three-operator system of the lac operon has been extremely useful (Gralla 1989) in this respect. Moreover, not only have the lac operators been helpful to understand the molecular mechanisms of gene regulation, but they have also been employed for other practical purposes. For instance, a technique for in vivo visualizing chromosome dynamics using lac operator–repressor binding has been proposed (Belmont & Straight 1998).
3.3 The β-galactosidase enzyme
Few genes have a history of study as long and distinguished as lacZ. The lacZ gene encodes an open reading frame of 1024 amino acids and is one of the first large genes to be completely sequenced. In E. coli, the biologically active β-galactosidase protein exists as a tetramer of four identical subunits and has a molecular weight of approximately 480–500 kDa. The primary enzymatic function of β-galactosidase relevant to its role as a biotechnological tool is to cleave the chemical bond between the anomeric carbon and glycosyl oxygen of appropriate substrates (Serebriiskii & Golemis 2000).
lacZ was chosen as the target of a very extensive early analysis, in part owing to specific experimental advantages accompanying work with β-galactosidase. These advantages continue to provide a rationale for using this protein in biotechnological applications today. They are given below (Serebriiskii & Golemis 2000).
Induction of β-galactosidase synthesis occurs over a large dynamic range (up to 10 000-fold over baseline levels with some inducers). This large range is achievable, in part, because the β-galactosidase protein can be tolerated at extremely high levels in E. coli, as well as in many other organisms such as yeasts, Caenorhabditis elegans, Drosophila melanogaster and mammals. Further, the β-galactosidase protein is readily purified by a number of relatively simple techniques, facilitating in vitro analysis of its activity.
A number of substrates (inducers) for β-galactosidase are either naturally available or very easily chemically synthesized. These enhance the development of models for β-galactosidase enzymatic activity and also provide a practical tool to finely modulate the expression or dissect catalytic activity of the β-galactosidase protein product.
β-Galactosidase activity is easily assayed, both in vivo and in vitro. Assays that have achieved prominence involve the use of colorimetric substrates in which the cleavage of specific β-d-galactopyranoside-coupled aglycone moieties releases coloured dyes. More recently, the panel of available colorimetric substrates has been augmented with fluorescent or chemiluminescent alternative substrates, which further expand sensitivity and applications.
The β-galactosidase protein is structurally malleable. It consists of three separable functional domains: alpha (α, amino-terminal), beta (β, central) and omega (ω, carboxy-terminal). Independent coexpression of the separated domains of the β-galactosidase protein successfully reconstitutes the activity of the full enzyme. This ability, as well as the additional capacity of β-galactosidase to function enzymatically when expressed as a translational fusion to a varied group of protein or peptide moieties, enables further applications.
These characteristics have allowed the usage of lacZ and its product (β-galactosidase) in many scientific and technological applications. Reviewing them all is quite beyond the scope of this paper, but the excellent review by Silhavy & Beckwith (1985) can be consulted. Further, those interested in learning more about the applications of the gene lacZ and the protein β-galactosidase can consult Silhavy & Beckwith (1985), Josephy (1996), Serebriiskii & Golemis (2000) and Shuman & Silhavy (2003).
3.4 The lac permease protein
Active transporters (pumps) require a cellular energy source (i.e. ATP hydrolysis) to catalyse the transport of charged components against an electrochemical gradient. Depending on their energy source, active transporters are classified as primary or secondary. Secondary transporters, in particular, use the free energy stored in a given electrochemical ion gradient (Abramson et al. 2004). LacY is a secondary transporter that couples free energy released from downhill translocation of protons to drive the uphill translocation of galactosides against a concentration gradient.
Lactose permease of E. coli (LacY) is composed of 417 amino acid residues and has 12 helices that transverse the membrane in zigzag fashion, connected by relatively hydrophilic loops with both N and C termini on the cytoplasm side. This protein is encoded by lacY, the second structural gene in the lac operon. lacY was the first gene encoding a membrane transport protein to be cloned into a recombinant plasmid, overexpressed and sequenced (see Kaback 2005 and references therein). This success in the early days of molecular biology opened the study of secondary active transport at the molecular level. Thus, LacY was the first protein of its class to be solubilized and purified in a completely functional state, thereby demonstrating that this single gene product is solely responsible for all the translocation reactions catalysed by the galactoside transport system in E. coli. It has also been shown that LacY is both structurally and functionally a monomer in the membrane (Kaback 2005).
Since the discovery of lactose permease, a number of molecular biological, biochemical and biophysical approaches have been used to study this protein. Analysis of this extensive data, and in particular of recent discoveries regarding the LacY structure and the properties of mutants in the irreplaceable residues (Kaback 2005), has led to the formulation of a model for this protein transport mechanism (Abramson et al. 2004). Furthermore, comparison of the structures of LacY and other MFS transporters (such as the Pi/glycerol-3-phosphate antiporter (GlpT)) has yielded valuable information on the functioning of secondary active transporters in general.
Abramson et al. (2004) and Kaback (2005) review the state of knowledge on secondary active transporters, of which the lactose permease is a paradigm.
In summary, we can see from the above considerations that the lac system has been extremely important and continues to advance our molecular understanding of genetic control and the relationship between sequence, structure and function.
4. Quantitative experimental and theoretical approaches
Two different interpretations of the lac operon dynamic performance existed in the beginning. Monod argued that the inducer concentration in the growing environment completely determines the operon induction level. On the other hand, Cohn & Horibata (1959) proposed a more subtle interpretation of their experiments. They suggested that the lac system provides an ‘experimental example of the Delbrück model’ (Delbrück 1949). According to Delbrück, biological systems with identical genotypes may display different behaviours under particular external conditions, due to ‘epigenetic’ differences that can be transmitted in the cell lineage in the absence of genetic modification. This hypothesis corresponds to a very early formulation of the general principle of phenotypic inheritance.
Novick & Weiner (1957) and Cohn & Horibata (1959) discovered the so-called ‘maintenance effect’, according to which a single cell may have two alternative states: induced, in which it can metabolize lactose, or uninduced, in which the corresponding genes are switched off and lactose metabolism does not occur. Their experimental protocol was as follows. First, a large amount of inducer was added to the extracellular medium of a culture of uninduced E. coli cells; then, the culture was split into two parts: U and I. Part U was immediately diluted, and so the cells remained uninduced; part I was diluted after several minutes, allowing the cells in this subculture to become induced. They further observed that, when induced cells were transferred to a medium with an intermediate ‘maintenance concentration’ of inducer, they and their progeny remained induced. Similarly, when uninduced cells were transferred to a medium with a ‘maintenance’ concentration, they and their progeny remained uninduced.
The ‘maintenance effect’ was interpreted as the consequence of a high permease concentration in induced cells, which would also have high inducer pumping efficiency. This would enable these cells to maintain the induced state and to transmit it to their progeny, even if placed in a medium with a low concentration of inducer. This interpretation accounts for the existence of two distinct phenotypes and provides an explanation of why induced cells placed in media with low inducer concentrations remain indefinitely induced, whereas cells that have never been induced stay uninduced. However, it does not explain what makes the cells switch between alternative states. This switching remained a mystery for a long time and it had to wait for the introduction of the concept of multistability to be fully explained.
Griffith (1968) developed a mathematical model (using ordinary differential equations, ODEs) for a single gene controlled by a positive feedback loop. He found that, under certain conditions, two stable states may be accessible for the system simultaneously. However, Griffith did not use his model to explain the maintenance effect of the lac operon. The first models that took into account all the relevant processes to unravel the dynamics of the lac operon were by Babloyantz & Sanglier (1972) and Nicolis & Prigogine (1977). Using a mathematical modelling approach, they interpreted the maintenance effect as the biological facet of the physical process of multistability. This model, involving a nonlinear feedback loop, accounted for the main behavioural features of the lactose–operon bistable transition. However, even though the mathematical description of the model required five differential equations (plus one conservation equation), the model did not take into account the detailed information available concerning molecular interactions between the operon components.
4.1 The Ozbudak et al. experiments
In the last few years, the interest in the bistable behaviour of the lactose operon has been renewed, and this is in part due to a paper recently published by Ozbudak et al. (2004). In it, the authors report the results of a series of ingenious experiments designed to study the bistable lac operon response when induced with lactose and the artificial non-metabolizable inducer thiomethylgalactoside (TMG). Ozbudak et al. incorporated a single copy of the green fluorescent protein gene (gfp) under the control of the lac promoter into the chromosome of E. coli. The cells also contained a plasmid encoding a red fluorescent reporter (HcRed) under the control of the galactitol (gat) promoter. This promoter includes a CRP-binding site, as well as a binding site for the galactitol repressor GatR. However, GatR is absent in E. coli. Therefore, transcription at the gat promoter, measured by red fluorescence, is a direct measure of CRP-cAMP levels. They further measured the response of single cells, initially in a given state of lac expression, to exposure to various combinations of glucose and TMG levels.
Ozbudak et al. report that, for a given concentration of extracellular glucose, the lac operon is uninduced at low TMG concentrations and fully induced at high TMG concentrations regardless of the cell's history. Between these switching thresholds, however, the system response is hysteretic (history dependent). By measuring the fluorescence of single cells, Ozbudak et al. obtained bimodal distributions between the switching thresholds, confirming the existence of bistability. When the experiments were repeated with lactose, instead of TMG, no evidence of bistability was found. The results in this paper not only confirmed bistability in the lac operon when induced with TMG, but also provided new and novel quantitative data that raise questions that may be answered via a modelling approach.
4.2 A minimal model
A number of mathematical models have been developed to investigate the dynamic behaviour of the lac operon. In §4.3, the characteristics of some of these models will be analysed by contrasting them with the minimal model introduced below.
Let M, E and L, respectively, denote the intracellular concentrations of mRNA, LacZ polypeptide and lactose. The differential equations governing the dynamics of these variables are
In the above equations, kM represents the maximum rate of mRNA production per cell, kE is the maximum translation initiation rate per mRNA and kE is the maximum lactose uptake rate per permease. Furthermore, γM, γE and γL, respectively, represent the degradation plus dilution rates of M, E and L. The function 𝒫R(L,Ge) is an increasing function of L (due to the positive regulatory feedback loop described above) and a decreasing function of the extracellular glucose concentration (Ge) owing to catabolite repression. 𝒫R(L,Ge) represents the probability that an mRNA polymerase is bound to the promoter and no repressor is bound to operator O1. The function 𝒫L(Le,Ge) represents the lac permease activity level. It is an increasing function of the external lactose concentration (Le) and a decreasing function of Ge due to inducer exclusion. Finally, ℳ(L) is the rate of lactose metabolism per β-galactosidase.
This simple model structure relies on the assumptions given below.
The processes governing the dynamics of all the chemical species other than M, E and L are assumed to be fast enough to make quasi-steady state approximations to the corresponding dynamic equations.
Since half of the lactose taken up is directly metabolized into glucose and galactose by β-galactosidase, while the rest is turned into allolactose (which is also later metabolized into glucose and galactose), it can be assumed that the intracellular concentrations of lactose and allolactose are very similar (see Santillán et al. 2007 for more details).
The translation and degradation rates of genes lacZ and lacY are assumed to be identical. Thus since β-galactosidase (lac permease) is a tetramer (monomer), its concentration is assumed to be one-quarter of (equal to) E.
It is important to note that the published lac operon models differ in the way the function 𝒫R(L,Ge) is formulated. In some cases, heuristic reasoning is used to propose Hill-type equations for 𝒫R(L,Ge). Some other models take into account, with different levels of detail, the interactions between the mRNA polymerase and the repressor molecules with the DNA chain to model this function. We discuss these differences in our review of various models in §4.3.
4.3 Recent modelling approaches
Wong et al. (1997) developed a 13-dimensional model for the lac operon. Besides the structural genes' mRNA and the intracellular lactose concentrations, the variables they consider are repressor mRNA and protein concentrations; β-galactosidase and permease concentrations (each governed by a different differential equation); the internal concentrations of allolactose, cAMP, glucose and phosphorylated glucose; and the external concentrations of glucose and lactose. Their model includes catabolite repression, inducer exclusion, lactose hydrolysis to glucose and galactose, synthesis and degradation of allolactose, and a variable growth rate. Wong et al. employed their model to study the diauxic growth of E. coli on glucose and lactose. For this, they tested different models for catabolite repression and the phosphorylation of the glucose produced from lactose hydrolysis, and analysed the influence of the model parameters on the two diauxic phases. Besides being quite detailed, this model has the virtue that most of the parameters in it were estimated from experimental data. However, even though Wong et al. considered the existence of the three known operator sites, they ignored their cooperative behaviour and incorrectly assumed that a repressor bound to any operator inhibits transcription initiation. Furthermore, Wong et al. also took into account the effect of the CAP activator by assuming that it must be bound to its specific site in the DNA regulatory region in order for the polymerase to bind the promoter.
Vilar et al. (2003) introduced a simple model of the lac operon to illustrate the applicability and limitations of mathematical modelling of the dynamics of cellular networks. In particular, they study the lac operon induction dynamics and its relation to bistability. Vilar et al. aimed at integrating three different levels of description (molecular, cellular and that of cell population) into a single model, and used it to investigate the system dynamics when an artificial (non-metabolizable) inducer is employed to activate it. In contrast to the minimal model, that of Vilar et al. lumps (through a quasi-steady state assumption) the mRNA dynamics into the equations governing the corresponding protein concentrations. It further accounts for the LacZ and LacY dynamics by means of two differential equations, assumes that LacY can be in either a non-functional or a functional state and includes one more equation for this last chemical species. Finally, this model takes into account neither catabolite repression nor inducer exclusion and accounts for the repression mechanism by means of a Hill-type equation.
Yildirim & Mackey (2003) investigated the bistable behaviour of the lac operon. For this, they introduced a five-dimensional mathematical model. The model of Yildirim and Mackey accounts for the dynamics of β-galactosidase and lac permease by means of two different differential equations, as well as for the dynamics of intracellular lactose and allolactose. This model also takes into account the delays due to transcription and translation processes. The authors paid particular attention to the estimation of the parameters in the model. They further tested their model against two sets of β-galactosidase activity versus time data, and against a set of data on β-galactosidase activity during periodic phosphate feeding. Their analytical and numerical studies indicate that for physiologically realistic values of external lactose and the bacterial growth rate, a regime exists where there may be bistability, and that this corresponds to a cusp bifurcation in the model dynamics. Deficiencies of Yildirim & Mackey's model are that it does not take into account catabolite repression or the inducer exclusion regulatory mechanisms. Furthermore, though they built the repression regulatory function by taking into account the repressor–operator and the polymerase–promoter interactions, they ignored the existence of three operators and considered operator O1 only.
In a later paper, Yildirim et al. (2004) attempted to identify as the origin of bistability one of the mechanisms involved in the regulation of the lac operon. To do this, they simplified the model presented in Yildirim & Mackey (2003) by ignoring permease dynamics and assuming a constant permease concentration. They numerically and analytically analysed the steady states of the reduced model and showed that it may indeed display bistability, depending on the extracellular lactose concentration and growth rate.
Santillán & Mackey (2004) developed a mathematical model of the lac operon, which accounts for all of the known regulatory mechanisms, including catabolite repression and inducer exclusion (both of which depend on external glucose concentrations), as well as the time delays inherent to transcription and translation. With this model, they investigated the influence of catabolite repression and inducer exclusion on the bistable behaviour of the lac operon. The model of Santillán & Mackey is six dimensional and the free variables are the lacZ and lacY mRNA concentrations, the β-galactosidase and lac permease concentrations, and the allolactose and cAMP concentrations. It is important to note that Santillán & Mackey's model considers all three known operators and the cooperativity among them, and that all the parameters in it were estimated from experimental data. In particular, they used a thermodynamic approach to model the interactions between the CAP activator, the repressor and the polymerase with their respective binding sites along the DNA chain, as well as the cooperative behaviour of the three known operators.
Van Hoek & Hogeweg (2006) carried out in silico simulation of the lac operon evolution in bacterial populations. From their results, the parameters that control the expression of the lac operon genes evolve in such a way that the system avoids bistability with respect to lactose, but does exhibit bistability with respect to artificial inducers. Thus, they argue from their computational experiments that the wild-type lac operon, which regulates lactose metabolism, is not a bistable switch under natural conditions. The model used by van Hoek & Hogeweg contains 10 independent differential equations and is based on the model of Wong et al. (1997). This model takes into account all known regulatory mechanisms. However, rather than considering the chemical details, van Hoek & Hogeweg modelled the repressor–DNA and the CAP activator–DNA interactions by means of Hill-type equations.
In a later paper, van Hoek & Hogeweg (2007) modified the lac operon model in Van Hoek & Hogeweg (2006) to incorporate stochasticity and study its effects from an evolutionary point of view. Through a mutation-selection process, they evolved the shape of the promoter function, and thus the effective amount of stochasticity. Van Hoek and Hogeweg concluded from their results that noise values for lactose, the natural inducer, are much lower than those for artificial, non-metabolizable inducers, because these artificial inducers experience a stronger positive feedback. They further showed that a high repression rate and hence high stochasticity increase the delay in lactose uptake in a variable environment. From this, the authors concluded that the lac operon has evolved such that the impact of stochastic gene expression is minor in its natural environment, but happens to respond with much stronger stochasticity when confronted with artificial inducers.
Santillán et al. (2007) investigated the origin of bistability in the lac operon. For this, they developed a mathematical model for the regulatory pathway in this system and compared the model predictions with the experimental results of Ozbudak et al. (2004). Santillán et al. examined the effect of lactose metabolism using this model, and showed that it greatly modifies the bistable region in the external lactose versus external glucose parameter space. The model also predicts that lactose metabolism can cause bistability to disappear for very low external glucose concentrations. The authors further carried out stochastic numerical simulations for several levels of external glucose and lactose and concluded from their results that bistability can help guarantee that E. coli consumes glucose and lactose in the most efficient possible way. Namely, the lac operon is induced only when there is almost no glucose in the growing medium, but if the external lactose is high, the operon induction level increases abruptly when the levels of glucose in the environment decrease to very low values. Finally, they demonstrated that this behaviour could not be obtained without bistability if the stability of the induced and uninduced states is to be preserved.
In a continuation of the work in Santillán et al. (2007), Santillán (2008) improved the mathematical model to account, in a more detailed way, for the interaction of the repressor molecules with the three lac operators. Besides, Santillán includes in the model a recently discovered cooperative interaction between the CAP molecule (an activator of the lactose operon) and operator 3, which influences DNA folding. Finally, this new model also includes the growth rate dependence on bacterial energy input rate in the form of transported glucose molecules and of metabolized lactose molecules. A large number of numerical experiments were carried out with the improved model, and the results are discussed along the same lines as in Santillán et al. (2007), including a detailed examination of the effect of a variable growth rate on the system dynamics. The models in both Santillán et al. (2007) and Santillán (2008) have the same structure as the minimal model above. Furthermore, both models take into account the chemical details of the repressor–DNA and CAP activator–DNA interactions, as well as the cooperativity observed between repressor molecules bound to different operators.
The models here reviewed are summarized in table 1. All of them deal with the bistable behaviour of the lac operon. However, most of the ones published prior to Ozbudak et al. (2004) only consider the use of lactose as inducer, and they predict that the lac operon shows bistability for physiological lactose concentrations. In this sense, the experimental work of Ozbudak et al. provided new data and opened new questions suitable for a mathematical modelling approach. One of these questions is why bistability cannot be observed when the lac operon natural inducer is employed. The most recent quantitative approaches have made use of Ozbudak et al.'s results to develop more accurate models, and two different answers to the above question have been proposed. Van Hoek & Hogeweg (2006, 2007) argue that bistability disappears altogether due to bacteria evolutionary adaptation to a fluctuating environment of glucose and lactose, while Santillán et al. (2007) and Santillán (2008) assert that bistability does not disappear but becomes extremely hard to identify with the experimental setup of Ozbudak et al. They, furthermore, discuss its significance from an evolutionary perspective. New experiments are needed to resolve this discrepancy.
Summary of the mathematical models of the lac operon that we have reviewed here. (The dimension number refers to the number of dependent variables. A model type can be either deterministic (D) or stochastic (S). Finally, the inducer column states whether...
Most of the models reviewed in this subsection involve ODEs. Given that chemical kinetics is the formalism behind ODE models, they are valid only when the molecule count (N) is such that is small enough. However, in the lac operon, the lacZ mRNA degradation rate is so high that the average number of mRNA molecules per bacterium is approximately 0.75, when the operon is fully induced (Santillán 2008). Furthermore, the lac repressor LacI is present in only approximately 10 tetramers per cell (Müller-Hill 1998). It follows from this that an essential aspect of modelling the lac operon is its stochastic nature. As seen in table 1, only the most recent models (Santillán et al. 2007; van Hoek & Hogeweg 2007; Santillán 2008) account for this inherent system stochasticity. Van Hoek & Hogeweg claim that noise has a large effect on the evolution of the lac operon: cells evolve such that noise has little effect on the system dynamical behaviour. By contrast, Santillán et al. do claim that noise has a large effect on the system dynamical behaviour. It is our belief that a huge amount of work remains to be done on this issue.
In this paper, we have presented a description of the regulatory mechanisms in the lac operon, taking into consideration the most recent discoveries. The system history has been surveyed as well, emphasizing the discovery of its main elements and the influence they have had on the development of molecular and systems biology. The operon bistable behaviour has also been analysed, including the discovery and origin of this complex phenomenon.
Multistability (of which bistability is the simplest example) corresponds to a true switch between alternate and coexisting steady states, and so allows a graded signal to be turned into a discontinuous evolution of the system along several different possible pathways. Multistability has certain unique properties not shared by other mechanisms of integrative control. These properties may play an essential role in the dynamics of living cells and organisms. Moreover, multistability has been invoked to explain catastrophic events in ecology (Rietkerk et al. 2004), mitogen-activated protein kinase (MAPK) cascades in animal cells (Ferrell & Machleder 1998; Bagowski & Ferrell 2001; Bhalla et al. 2002), cell cycle regulatory circuits in Xenopus and Saccharomyces cerevisiae (Cross et al. 2002; Pomerening et al. 2003), the generation of switch-like biochemical responses (Ferrell & Machleder 1998; Bagowski & Ferrell 2001; Bagowski et al. 2003), and the establishment of cell cycle oscillations and mutually exclusive cell cycle phases (Pomerening et al. 2003; Sha et al. 2003), among other biological phenomena.
Not only was the lac operon the first system in which bistability was discovered but, as the literature reviewed in this paper demonstrates, it has been and is still one of the ideal model systems to analyse this complex behaviour. According to the quantitative studies we have reviewed, the lac operon has helped to understand the origin, biological implications and subtleness of bistability; it may also help to tackle similar questions in other systems and organisms. For this and other reasons not addressed in this paper, it is our opinion that the lactose operon may be as influential in the development of the nascent field of systems biology as it was in the development of molecular biology.
The different mathematical models of the lac operon here reviewed present a good example of what the philosophy of model making is and how mathematical models can influence the development of a given scientific discipline. To discuss this issue we make extensive use of the excellent essay by Rosenblueth & Wiener (1945), which we shall quote a number of times in the forthcoming paragraphs.
According to Rosenblueth and Wiener, models are a central necessity of scientific procedure because no substantial part of the universe is so simple that it can be grasped and controlled without abstraction, abstraction being replacing the part of the universe under consideration by a model of similar but simpler structure. Rosenblueth and Wiener further classify scientific models as either material or formal. In their view, a material model is the representation of a complex system by a system that is assumed simpler and is also assumed to have some properties similar to those selected for study in the original complex system. By contrast, a formal model is a symbolic assertion in logical terms of an idealized relatively simple situation sharing the structural properties of the original factual system. Here we are concerned with formal models of which the mathematical ones are a subset.
Closed boxes (in which a finite number of output variables are causally related to a finite number of input variables, without the knowledge of the detailed mechanisms inside the box) are often employed in formal models. Indeed, according to Rosenblueth and Wiener, all scientific problems begin as closed-box problems (i.e. only a few of the significant variables are recognized), and scientific progress consists of a progressive opening of those boxes. The successive addition of variables gradually leads to more elaborate theoretical models, hence to a hierarchy in these models, from relatively simple, highly abstract ones to more complex, more concrete theoretical structures. The setting up of a simple model for a closed box is based on the assumption that a number of variables are only loosely coupled with the rest of those belonging to the system. As the successive models become progressively more sophisticated, the number of closed regions usually increases, because the process may be compared to the subdivision of an original single box into several smaller shut compartments. Many of these small compartments may be deliberately left closed, because they are considered only functionally, but not structurally, important.
Early closed-box models, such as that of Griffith (1968), acknowledged that a gene subjected to positive feedback regulation can show bistability. However, since having positive feedback does not guarantee bistability, these models could not predict whether the lac operon would show this behaviour and, if so, what the responsible mechanisms are. More detailed models were needed to address these questions, and they have been developed as the necessary experimental information is available now. As foreseen by Rosenblueth and Wiener more than 60 years ago, these more detailed models are heterogeneous assemblies of elements, some treated in detail (that is, specifically or structurally) and some treated merely with respect to their overall performance (that is, generically or functionally). The models by Santillán et al. and van Hoek and Hogeweg are good examples. For instance, while the Santillán et al. models take into consideration the details of the polymerase–DNA, the repressor–DNA and the activator–DNA interactions, the models of van Hoek and Hogeweg lump them together into a closed box. On the other hand, van Hoek and Hogeweg model with more detail the dynamics of the lactose, allolactose, glucose, cAMP, and ATP intracellular concentrations, and this is reflected in the number of independent equations that the van Hoek and Hogeweg and the Santillán et al. models have. Finally, while van Hoek and Hogeweg use an evolutionary modelling approach to study bistability in the lac operon, Santillán et al. use a more static approach, trying to model the lac operon in full detail. To develop more elaborate models, more accurate quantitative data on the system components and their interactions are required. Ideally, a cycle of modelling and experimental efforts shall continue, with one of the outcomes being models progressively more sophisticated and capable of addressing more specific questions. However, there is no point in carrying this process out until its obvious limit. To explain this, we refer once more to Rosenblueth & Wiener (1945), who asserted that as a model becomes more detailed and accurate, it will tend to become identical with the original system. As a limit, it will become the system itself. The ideal model would then be one which agrees with the system in its full complexity and which leaves no closed boxes. However, any one capable of elaborating and comprehending such a model in its entirety would find it unnecessary, because he/she could then grasp the complete system directly as a whole.
We hope that readers of this paper will appreciate that mathematical modelling is a process that constantly evolves as the predictions of the models are iterated against laboratory data. The results of the past three decades in modelling the dynamics of the lac operon exemplify this. The reader will, no doubt, also realize that each model has its positive and negative aspects. The level of detail of the model depends on the availability and quality of the data and also on the questions we want to address. The more the detail, the more complicated the model will be. A mathematical analysis might then be hard to undertake and the conclusions may only be based on numerical experiments that many, including us, find less than satisfactory. On the other hand, a simple model may be easier to analyse and a mathematical analysis can give more insights into the dynamical properties or the underlying system; but it may oversimplify and fail to capture some important features of the reality.
The issue of model complexity is intimately tied to the issue of the dimensionality of the parameter space, and this is tied directly to one of the quandaries that face every modeller. The more complex the model, the more parameters must be estimated. It is a virtual truism in mathematical biology that one is almost never able to obtain all of the parameters in a model from the same laboratory setting using the same procedures and techniques and subjects. So, as mathematical model construction is something of an art in itself, the same can be said for parameter estimation. The senior author (M.C.M.) with over 45 years of experience in mathematical biology suggests that the hardest part of the modelling exercise is obtaining decent parameter estimations.
This research was partially supported by Consejo Nacional de Ciencia y Tecnología (CONACyT, Mexico) under grant 55228, the Natural Sciences and Engineering Research Council (NSERC, Canada) and the Mathematics of Information Technology and Complex Systems (MITACS, Canada). We thank the anonymous referees whose comments and criticisms greatly helped us to improve this paper.
One contribution of 10 to a Theme Supplement ‘Biological switches and clocks’.
- Abramson J., Iwata S., Kaback H.R. Lactose permease as a paradigm for membrane transport proteins (review) Mol. Membr. Biol. 2004;21:227–236. doi:10.1080/09687680410001716862[PubMed]
- Babloyantz A., Sanglier M. Chemical instabilities of “all-or-none” type in beta-galactosidase induction and active transport. FEBS Lett. 1972;23:364–366. doi:10.1016/0014-5793(72)80317-X[PubMed]
- Bagowski C.P., Ferrell J.E. Bistability in the JNK cascade. Curr. Biol. 2001;11:1176–1182. doi:10.1016/S0960-9822(01)00330-X[PubMed]
- Bagowski C., Besser J., Frey C.R., Ferrell J.E., Jr The JNK cascade as a biochemical switch in mammalian cells: ultrasensitive and all-or-none responses. Curr. Biol. 2003;13:315–320. doi:10.1016/S0960-9822(03)00083-6[PubMed]
- Beckwith, J. 1987 The lactose operon. In Escherichia coli and Salmonella thyphymurium: cellular and molecular biology, vol. 2 (eds F. C. Neidhart, J. L. Ingraham, K.B. Low, B. Magasanik & H. E. Umbarger), pp. 1439–1443. Washington, DC: American Society for Microbiology.
- Belmont A.S., Straight F. In vivo visualization of chromosomes using lac operator–repressor binding. Trends Cell Biol. 1998;8:121–124. doi:10.1016/S0962-8924(97)01211-7[PubMed]
- Bhalla U.S., Ram P.T., Iyengar R. MAP kinase phosphatase as a locus of flexibility in a mitogen-activated protein kinase signaling network. Science. 2002;297:1018–1023. doi:10.1126/science.1068873[PubMed]
- Cohn M., Horibata K. Analysis of the differentiation and of the heterogeneity within a population of Escherichia coli undergoing induced beta-galactosidase synthesis. J. Bacteriol. 1959;78:613–623.[PMC free article][PubMed]
- Cross F.R., Archambault V., Miller M., Klovstad M. Testing a mathematical model of the yeast cell cycle. Mol. Biol. Cell. 2002;13:52–70. doi:10.1091/mbc.01-05-0265[PMC free article][PubMed]
- Delbrück, M. 1949 Discussion. In Unités Biologiques douées de Continuité Génétique, pp. 33–34. Paris, France: Editions du CNRS.
- Ferrell J.E., Jr, Machleder E.M. The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science. 1998;280:895–898. doi:10.1126/science.280.5365.895[PubMed]
- Gilbert W., Maxam A. The nucleotide sequence of the lac repressor. Proc. Natl Acad. Sci. USA. 1973;70:3581–3584. doi:10.1073/pnas.70.12.3581[PMC free article][PubMed]
- Gralla J.D. Bacterial gene regulation from distant DNA sites. Cell. 1989;57:193–195. doi:10.1016/0092-8674(89)90955-0[PubMed]
- Griffith J.S. Mathematics of cellular control processes. II. Positive feedback to one gene. J. Theor. Biol. 1968;20:209–216. doi:10.1016/0022-5193(68)90190-2[PubMed]
- Jacob F., Monod J. Genetic regulatory mechanisms in the synthesis of proteins. J. Mol. Biol. 1961a;3:318–356.[PubMed]
- Jacob F., Monod J. On the regulation of gene activity. Cold Spring Harb. Symp. Quant. Biol. 1961b;26:193–211.
- Jacob F., Perrin D., Sanchez C., Monod J. Operon: a group of genes with the expression coordinated by an operator. C. R. Hebd. Seances Acad. Sci. 1960;250:1727–1729.[PubMed]
- Josephy P.D. The Escherichia coli lacZ reversion mutagenicity assay. Mutat. Res. 1996;455:71–80.[PubMed]
- Kaback H.R. Structure and mechanism of the lactose permease. C. R. Biol. 2005;328:557–567. doi:10.1016/j.crvi.2005.03.008[PubMed]
- Kuhlman T., Zhang Z., Saier M.H., Hwa T. Combinatorial transcriptional control of the lactose operon of Escherichia coli. Proc. Natl Acad. Sci. USA. 2007;104:6043–6048. doi:10.1073/pnas.0606717104[PMC free article][PubMed]
- Lendenmann U., Snozzi M., Egli T. Kinetics of the simultaneous utilization of sugar mixtures by Escherichia coli in continuous culture. Appl. Environ. Microbiol. 1996;62:1493–1499.[PMC free article][PubMed]
- Lewis M. The lac repressor. C. R. Biol. 2005;328:521–548. doi:10.1016/j.crvi.2005.04.004[PubMed]