Gene regulation at the right operator (OR) of bacteriophage lambda. III. lambda repressor directly activates gene transcription

RNA polymerase efficiently transcribes DNA-scaffolded, cooperative bacteriophage repressor complexes

DNA can act as a scaffold for the cooperative binding of protein oligomers. For example, the phage 186 CI repressor forms a wheel of seven dimers wrapped in DNA with specific binding sites, while phage λ CI repressor dimers bind to two well-separated sets of operators, forming a DNA loop. Atomic force microscopy was used to measure transcription elongation by E. coli RNA polymerase through these protein complexes. 186 CI, or λ CI, bound along unlooped DNA negligibly interfered with transcription by RNAP. Wrapped and looped topologies induced by these scaffolded, cooperatively bound repressor oligomers did not form significantly better roadblocks to transcription. Thus, despite binding with high affinity, these repressors are not effective roadblocks to transcription.

Bistable Expression of a Toxin-Antitoxin System Located in a Cryptic Prophage of Escherichia coli O157:H7

Type II toxin-antitoxin (TA) systems are classically composed of two genes that encode a toxic protein and a cognate antitoxin protein. Both genes are organized in an operon whose expression is autoregulated at the level of transcription by the antitoxin-toxin complex, which binds operator DNA through the antitoxin’s DNA-binding domain. Here, we investigated the transcriptional regulation of a particular TA system located in the immunity region of a cryptic lambdoid prophage in the Escherichia coli O157:H7 EDL933 strain. This noncanonical paaA2-parE2 TA operon contains a third gene, paaR2, that encodes a transcriptional regulator that was previously shown to control expression of the TA. We provide direct evidence that the PaaR2 is a transcriptional regulator which shares functional similarities to the lambda CI repressor. Expression of the paaA2-parE2 TA operon is regulated by two other transcriptional regulators, YdaS and YdaT, encoded within the same region. We argue that YdaS and YdaT are analogous to lambda Cro and CII and that they do not constitute a TA system, as previously debated. We show that PaaR2 primarily represses the expression of YdaS and YdaT, which in turn controls the expression of paaR2-paaA2-parE2 operon. Overall, our results show that the paaA2-parE2 TA is embedded in an intricate lambdoid prophage-like regulation network. Using single-cell analysis, we observed that the entire locus exhibits bistability, which generates diversity of expression in the population. Moreover, we confirmed that paaA2-parE2 is addictive and propose that it could limit genomic rearrangements within the immunity region of the CP-933P cryptic prophage.

A precisely adjustable, variation-suppressed eukaryotic transcriptional controller to enable genetic discovery

Conditional expression of genes and observation of phenotype remain central to biological discovery. Current methods enable either on/off or imprecisely controlled graded gene expression. We developed a 'well-tempered' controller, WTC₈₄₆, for precisely adjustable, graded, growth condition independent expression of genes in Saccharomyces cerevisiae. Controlled genes are expressed from a strong semisynthetic promoter repressed by the prokaryotic TetR, which also represses its own synthesis; with basal expression abolished by a second, 'zeroing' repressor. The autorepression loop lowers cell-to-cell variation while enabling precise adjustment of protein expression by a chemical inducer. WTC₈₄₆ allelic strains in which the controller replaced the native promoters recapitulated known null phenotypes (CDC42, TPI1), exhibited novel overexpression phenotypes (IPL1), showed protein dosage-dependent growth rates and morphological phenotypes (CDC28, TOR2, PMA1 and the hitherto uncharacterized PBR1), and enabled cell cycle synchronization (CDC20). WTC₈₄₆ defines an 'expression clamp' allowing protein dosage to be adjusted by the experimenter across the range of cellular protein abundances, with limited variation around the setpoint.

Thermoinducible expression system for producing recombinant proteins in Escherichia coli: advances and insights

Recombinant protein (RP) production from Escherichia coli has been extensively studied to find strategies for increasing product yields. The thermoinducible expression system is commonly employed at the industrial level to produce various RPs which avoids the addition of chemical inducers, thus minimizing contamination risks. Multiple aspects of the molecular origin and biotechnological uses of its regulatory elements (pL/pR promoters and cI857 thermolabile repressor) derived from bacteriophage λ provide knowledge to improve the bioprocesses using this system. Here, we discuss the main aspects of the potential use of the λpL/pR-cI857 thermoinducible system for RP production in E. coli, focusing on the approaches of investigations that have contributed to the advancement of this expression system. Metabolic and physiological changes that occur in the host cells caused by heat stress and by RP overproduction are also described. Therefore, the current scenario and the future applications of systems that use heat to induce RP production is discussed to understand the relationship between the activation of the bacterial heat shock response, RP accumulation, and its possible aggregation to form inclusion bodies.

Biophysical ambiguities prevent accurate genetic prediction

A goal of biology is to predict how mutations combine to alter phenotypes, fitness and disease. It is often assumed that mutations combine additively or with interactions that can be predicted. Here, we show using simulations that, even for the simple example of the lambda phage transcription factor CI repressing a gene, this assumption is incorrect and that perfect measurements of the effects of mutations on a trait and mechanistic understanding can be insufficient to predict what happens when two mutations are combined. This apparent paradox arises because mutations can have different biophysical effects to cause the same change in a phenotype and the outcome in a double mutant depends upon what these hidden biophysical changes actually are. Pleiotropy and non-monotonic functions further confound prediction of how mutations interact. Accurate prediction of phenotypes and disease will sometimes not be possible unless these biophysical ambiguities can be resolved using additional measurements.

Elements in the λ immunity region regulate phage development: beyond the 'Genetic Switch'

Genetic elements in the bacteriophage λ immunity region contribute to stable maintenance and synchronous induction of the integrated Escherichia coli prophage. There is a bistable switch between lysogenic and lytic growth that is orchestrated by the CI and Cro repressors acting on the lytic (PL and PR ) and lysogenic (PRM ) promoters, referred to as the Genetic Switch. Other less well-characterized elements in the phage immunity region include the PLIT promoter and the immunity terminator, TIMM . The PLIT promoter is repressed by the bacterial LexA protein in λ lysogens. LexA repressor, like the λ CI repressor, is inactivated during the SOS response to DNA damage, and this regulation ensures that the PLIT promoter and the lytic PL and PR promoters are synchronously activated. Proper RexA and RexB protein levels are critical for the switch from lysogeny to lytic growth. Mutation of PLIT reduces RexB levels relative to RexA, compromising cellular energetics and causing a 10-fold reduction in lytic phage yield. The RexA and RexB proteins interact with themselves and each other in a bacterial two-hybrid system. We also find that the transcription terminator, TIMM , is a Rho-independent, intrinsic terminator. Inactivation of TIMM has minimal effect on λ lysogenization or prophage induction.

Sex and death: from cell fate specification to dynamic control of X-chromosome structure and gene expression

Determining sex is a binary developmental decision that most metazoans must make. Like many organisms, Caenorhabditis elegans specifies sex (XO male or XX hermaphrodite) by tallying X-chromosome number. We dissected this precise counting mechanism to determine how tiny differences in concentrations of signals are translated into dramatically different developmental fates. Determining sex by counting chromosomes solved one problem but created another-an imbalance in X gene products. We found that nematodes compensate for the difference in X-chromosome dose between sexes by reducing transcription from both hermaphrodite X chromosomes. In a surprising feat of evolution, X-chromosome regulation is functionally related to a structural problem of all mitotic and meiotic chromosomes: achieving ordered compaction of chromosomes before segregation. We showed the dosage compensation complex is a condensin complex that imposes a specific three--dimensional architecture onto hermaphrodite X chromosomes. It also triggers enrichment of histone modification H4K20me1. We discovered the machinery and mechanism underlying H4K20me1 enrichment and demonstrated its pivotal role in regulating higher-order X-chromosome structure and gene expression.

The Developmental Switch in Bacteriophage λ: A Critical Role of the Cro Protein

Bacteriophage λ of Escherichia coli has two alternative life cycles after infection-host survival with lysogen formation, or host lysis and phage production. In a lysogen, CI represses the two lytic promoters, pR and pL, and activates its own transcription from the auto-regulated pRM promoter. During induction from the lysogenic to lytic state, CI is inactivated, and the two lytic promoters are de-repressed allowing for expression of Cro from pR. Cro is known to repress transcription of CI from pRM to prevent lysogeny. We show here that when Cro and CI are both present but at low levels, the low level of Cro initially stimulates the lytic promoters while CI repressor is still present, stimulating the level of Cro to a concentration required for pRM repression. Cro has no stimulatory effect without the presence of CI. We propose that this early auto-activating role of Cro at lower concentrations is essential in the developmental switch to lytic growth, whereas pRM repression by Cro at relatively higher concentrations avoids restoring lysogeny.

New Insights into the Phage Genetic Switch: Effects of Bacteriophage Lambda Operator Mutations on DNA Looping and Regulation of PR, PL, and PRM

One of the best understood systems in genetic regulatory biology is the so-called "genetic switch" that determines the choice the phage-encoded CI repressor binds cooperatively to tripartite operators, O_L and O_R, in a defined pattern, thus blocking the transcription at two lytic promoters, P_L and P_R, and auto-regulating the promoter, P_RM, which directs CI synthesis by the prophage. Fine-tuning of the maintenance of lysogeny is facilitated by interactions between CI dimers bound to O_R and O_L through the formation of a loop by the intervening DNA segment. By using a purified in vitro transcription system, we have genetically dissected the roles of individual operator sites in the formation of the DNA loop and thus have gained several new and unexpected insights into the system. First, although both O_R and O_L are tripartite, the presence of only a single active CI binding site in one of the two operators is sufficient for DNA loop formation. Second, in P_L, unlike in P_R, the promoter distal operator site, O_L3, is sufficient to directly repress P_L. Third, DNA looping mediated by the formation of CI octamers arising through the interaction of pairs of dimers bound to adjacent operator sites in O_R and O_L does not require O_R and O_L to be aligned "in register", that is, CI bound to "out-of-register" sub-operators, for example, O_L1~O_l2 and O_R2~O_R3, can also mediate loop formation. Finally, based on an examination of the mechanism of activation of P_RM when only O_R1 or O_R2 are wild type, we hypothesize that RNA polymerase bound at P_R interferes with DNA loop formation. Thus, the formation of DNA loops involves potential interactions between proteins bound at numerous cis-acting sites, which therefore very subtly contribute to the regulation of the "switch".

Long-range transcriptional interference in E. coli used to construct a dual positive selection system for genetic switches

We have investigated transcriptional interference between convergent genes in E. coli and demonstrate substantial interference for inter-promoter distances of as far as 3 kb. Interference can be elicited by both strong σ⁷⁰ dependent and T7 promoters. In the presented design, a strong promoter driving gene expression of a ‘forward’ gene interferes with the expression of a ‘reverse’ gene by a weak promoter. This arrangement allows inversely correlated gene expression without requiring further regulatory components. Thus, modulation of the activity of the strong promoter alters expression of both the forward and the reverse gene. We used this design to develop a dual selection system for conditional operator site binding, allowing positive selection both for binding and for non-binding to DNA. This study demonstrates the utility of this novel system using the Lac repressor as a model protein for conditional DNA binding, and spectinomycin and chloramphenicol resistance genes as positive selection markers in liquid culture. Randomized LacI libraries were created and subjected to subsequent dual selection, but mispairing IPTG and selection cues in respect to the wild-type LacI response, allowing the isolation of a LacI variant with a reversed IPTG response within three rounds of library generation and dual selection.

Yet another way that phage λ manipulates its Escherichia coli host: λrexB is involved in the lysogenic-lytic switch

The life cycle of phage λ has been studied extensively. Of particular interest has been the process leading to the decision of the phage to switch from lysogenic to lytic cycle. The principal participant in this process is the λcI repressor, which is cleaved under conditions of DNA damage. Cleaved λcI no longer acts as a repressor, allowing phage λ to switch from its lysogenic to lytic cycle. The well-known mechanism responsible for λcI cleavage is the SOS response. We have recently reported that the Escherichia coli toxin-antitoxin mazEF pathway inhibits the SOS response; in fact, the SOS response is permitted only in E. coli strains deficient in the expression of the mazEF pathway. Moreover, in strains lysogenic for prophage λ, the SOS response is enabled by the presence of λrexB. λRexB had previously been found to inhibit the degradation of the antitoxin MazE, thereby preventing the toxic action of MazF. Thus, phage λ rexB gene not only safeguards the prophage state by preventing death of its E. coli host but is also indirectly involved in the lysogenic-lytic switch.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

Natural biased coin encoded in the genome determines cell strategy

Decision making at a cellular level determines different fates for isogenic cells. However, it is not yet clear how rational decisions are encoded in the genome, how they are transmitted to their offspring, and whether they evolve and become optimized throughout generations. In this paper, we use a game theoretic approach to explain how rational decisions are made in the presence of cooperators and competitors. Our results suggest the existence of an internal switch that operates as a biased coin. The biased coin is, in fact, a biochemical bistable network of interacting genes that can flip to one of its stable states in response to different environmental stimuli. We present a framework to describe how the positions of attractors in such a gene regulatory network correspond to the behavior of a rational player in a competing environment. We evaluate our model by considering lysis/lysogeny decision making of bacteriophage lambda in E. coli.

Effect of supercoiling on the λ switch

The lysogenic state of the λ switch is exceptionally stable, still, it is capable of responding to DNA-damage and rapidly enter the lytic state. We invented an assay where PNA mediated tethering of a plasmid allowed for single molecule investigations of the effect of supercoiling on the efficiency of the epigenetic λ switch. Compared with non-supercoiled DNA, the presence of supercoils enhances the CI-mediated DNA looping probability and renders the transition between the looped and unlooped states steeper, thus increasing the Hill coefficient. Interestingly, the transition occurs exactly at the CI concentration corresponding to the minimum number of CI molecules capable of maintaining the pRM-repressed state. Based on these results we propose that supercoiling maintains the pRM-repressible state as CI concentration decline during induction and thus prevent autoregulation of cI from interfering with induction.

Transcription-factor-mediated DNA looping probed by high-resolution, single-molecule imaging in live E. coli cells

DNA looping mediated by transcription factors plays critical roles in prokaryotic gene regulation. The "genetic switch" of bacteriophage λ determines whether a prophage stays incorporated in the E. coli chromosome or enters the lytic cycle of phage propagation and cell lysis. Past studies have shown that long-range DNA interactions between the operator sequences O(R) and O(L) (separated by 2.3 kb), mediated by the λ repressor CI (accession number P03034), play key roles in regulating the λ switch. In vitro, it was demonstrated that DNA segments harboring the operator sequences formed loops in the presence of CI, but CI-mediated DNA looping has not been directly visualized in vivo, hindering a deep understanding of the corresponding dynamics in realistic cellular environments. We report a high-resolution, single-molecule imaging method to probe CI-mediated DNA looping in live E. coli cells. We labeled two DNA loci with differently colored fluorescent fusion proteins and tracked their separations in real time with ∼40 nm accuracy, enabling the first direct analysis of transcription-factor-mediated DNA looping in live cells. Combining looping measurements with measurements of CI expression levels in different operator mutants, we show quantitatively that DNA looping activates transcription and enhances repression. Further, we estimated the upper bound of the rate of conformational change from the unlooped to the looped state, and discuss how chromosome compaction may impact looping kinetics. Our results provide insights into transcription-factor-mediated DNA looping in a variety of operator and CI mutant backgrounds in vivo, and our methodology can be applied to a broad range of questions regarding chromosome conformations in prokaryotes and higher organisms.

Multilevel autoregulation of λ repressor protein CI by DNA looping in vitro

The prophage state of bacteriophage λ is extremely stable and is maintained by a highly regulated level of λ repressor protein, CI, which represses lytic functions. CI regulates its own synthesis in a lysogen by activating and repressing its promoter, P(RM). CI participates in long-range interactions involving two regions of widely separated operator sites by generating a loop in the intervening DNA. We investigated the roles of each individual site under conditions that permitted DNA loop formation by using in vitro transcription assays for the first time on supercoiled DNA that mimics in vivo situation. We confirmed that DNA loops generated by oligomerization of CI bound to its operators influence the autoactivation and autorepression of P(RM) regulation. We additionally report that different configurations of DNA loops are central to this regulation--one configuration further enhances autoactivation and another is essential for autorepression of P(RM).

The lysis-lysogeny decision of bacteriophage 933W: a 933W repressor-mediated long-distance loop has no role in regulating 933W P(RM) activity

Our data show that unlike bacteriophage λ, repressor bound at O(L) of bacteriophage 933W has no role in regulation of 933W repressor occupancy of 933W O(R)3 or the transcriptional activity of 933W P(RM). This finding suggests that a cooperative long-range loop between repressor tetramers bound at O(R) and O(L) does not form in bacteriophage 933W. Nonetheless, 933W forms lysogens, and 933W prophage display a threshold response to UV induction similar to related lambdoid phages. Hence, the long-range loop thought to be important for constructing a threshold response in lambdoid bacteriophages is dispensable. The lack of a loop requires bacteriophage 933W to use a novel strategy in regulating its lysis-lysogeny decisions. As part of this strategy, the difference between the repressor concentrations needed to bind O(R)2 and activate 933W P(RM) transcription or bind O(R)3 and repress transcription from P(RM) is <2-fold. Consequently, P(RM) is never fully activated, reaching only ∼25% of the maximum possible level of repressor-dependent activation before repressor-mediated repression occurs. The 933W repressor also apparently does not bind cooperatively to the individual sites in O(R) and O(L). This scenario explains how, in the absence of DNA looping, bacteriophage 933W displays a threshold effect in response to DNA damage and suggests how 933W lysogens behave as "hair triggers" with spontaneous induction occurring to a greater extent in this phage than in other lambdoid phages.

Defining cooperativity in gene regulation locally through intrinsic noise

Regulatory networks in cells may comprise a variety of types of molecular interactions. The most basic are pairwise interactions, in which one species controls the behaviour of another (e.g. a transcription factor activates or represses a gene). Higher-order interactions, while more subtle, may be important for determining the function of networks. Here, the authors systematically expand a simple master equation model for a gene to derive an approach for robustly assessing the cooperativity (effective copy number) with which a transcription factor acts. The essential idea is that moments of a joint distribution of protein copy numbers determine the Hill coefficient of a cis-regulatory input function without non-linear fitting. The authors show that this method prescribes a definition of cooperativity that is meaningful even in highly complex situations in which the regulation does not conform to a simple Hill function. To illustrate the utility of the method, the authors measure the cooperativity of the transcription factor CI in simulations of phage- and show how the cooperativity accurately reflects the behaviour of the system. The authors numerically assess the effects of deviations from ideality, as well as possible sources of error. The relationship to other definitions of cooperativity and issues for experimentally realising the procedure are discussed.

Probability landscape of heritable and robust epigenetic state of lysogeny in phage lambda

Computational studies of biological networks can help to identify components and wirings responsible for observed phenotypes. However, studying stochastic networks controlling many biological processes is challenging. Similar to Schrödinger's equation in quantum mechanics, the chemical master equation (CME) provides a basic framework for understanding stochastic networks. However, except for simple problems, the CME cannot be solved analytically. Here we use a method called discrete chemical master equation (dCME) to compute directly the full steady-state probability landscape of the lysogeny maintenance network in phage lambda from its CME. Results show that wild-type phage lambda can maintain a constant level of repressor over a wide range of repressor degradation rate and is stable against UV irradiation, ensuring heritability of the lysogenic state. Furthermore, it can switch efficiently to the lytic state once repressor degradation increases past a high threshold by a small amount. We find that beyond bistability and nonlinear dimerization, cooperativity between repressors bound to O(R)1 and O(R)2 is required for stable and heritable epigenetic state of lysogeny that can switch efficiently. Mutants of phage lambda lack stability and do not possess a high threshold. Instead, they are leaky and respond to gradual changes in degradation rate. Our computation faithfully reproduces the hair triggers for UV-induced lysis observed in mutants and the limitation in robustness against mutations. The landscape approach computed from dCME is general and can be applied to study broad issues in systems biology.

The bacteriophage lambda CI protein finds an asymmetric solution

The CI protein of bacteriophage lambda (lambdaCI) is both a repressor and activator of transcription that has served as a model for understanding how gene regulatory proteins work. A dimeric DNA-binding protein, lambdaCI also forms higher-order oligomers that allow it to bind cooperatively to both adjacent and nonadjacent operator sites within the phage genome. The ability of phage lambda to transition efficiently from one program of gene expression to another depends upon the formation of these higher-order protein-DNA complexes. A recently determined crystal structure of a DNA-bound lambdaCI dimer reveals that the two subunits of the dimer adopt different conformations. This unexpected asymmetry helps explain how these higher-order complexes are assembled.

Sequence-dependent upstream DNA-RNA polymerase interactions in the open complex with lambdaPR and lambdaPRM promoters and implications for the mechanism of promoter interference

Upstream interactions of Escherichia coli RNA polymerase (RNAP) in an open promoter complex (RPo) formed at the P(R) and P(RM) promoters of bacteriophage lambda have been studied by atomic force microscopy. We demonstrate that the previously described 30-nm DNA compaction observed upon RPo formation at P(R) [Rivetti, C., Guthold, M. & Bustamante, C. (1999). Wrapping of DNA around the E. coli RNA polymerase open promoter complex. EMBO J., 18, 4464-4475.] is a consequence of the specific interaction of the RNAP with two AT-rich sequence determinants positioned from -36 to -59 and from -80 to -100. Likewise, RPos formed at P(RM) showed a specific contact between RNAP and the upstream DNA sequence. We further demonstrate that this interaction, which results in DNA wrapping against the polymerase surface, is mediated by the C-terminal domains of alpha-subunits (carboxy-terminal domain). Substitution of these AT-rich sequences with heterologous DNA reduces DNA wrapping but has only a small effect on the activity of the P(R) promoter. We find, however, that the frequency of DNA templates with both P(R) and P(RM) occupied by an RNAP significantly increases upon loss of DNA wrapping. These results suggest that alpha carboxy-terminal domain interactions with upstream DNA can also play a role in regulating the expression of closely spaced promoters. Finally, a model for a possible mechanism of promoter interference between P(R) and P(RM) is proposed.

DNA looping can enhance lysogenic CI transcription in phage lambda

The lysogenic state of bacteriophage lambda is maintained by CI repressor, which negatively regulates two promoters to block lytic gene expression. Expression of CI is itself controlled by positive and negative feedback as CI binds to O(R) to regulate the P(RM) promoter. In addition to direct interactions with operator DNA, CI tetramers bound at O(L) and O(R) can come together to form an octamer, looping the DNA that lies between them and allowing O(L) to assist with negative regulation of P(RM). We used a fluorescent reporter protein to measure the CI concentration for a set of constructs that differ in their ability to assume various forms of the looped structure. Based on the observed steady-state fluorescence for these constructs, the presence of O(L) increases P(RM) activation unless both operators can be fully occupied. By calculating the probabilities for the underlying operator configurations present in each construct, two different models for the mechanism of enhanced activation allow us to predict that when the DNA is looped, P(RM) activation can be 2- to 4-fold higher than is possible for unlooped DNA. Based on our results, transcriptional regulation for lambda's lysogenic/lytic switch includes both activation and repression due to DNA looping.

Positive autoregulation of cI is a dispensable feature of the phage lambda gene regulatory circuitry

Complex gene regulatory circuits contain many features that are likely to contribute to their operation. It is unclear, however, whether all these features are necessary for proper circuit behavior or whether certain ones are refinements that make the circuit work better but are dispensable for qualitatively normal behavior. We have addressed this question using the phage lambda regulatory circuit, which can persist in two stable states, the lytic state and the lysogenic state. In the lysogenic state, the CI repressor positively regulates its own expression by stimulating transcription from the P(RM) promoter. We tested whether this feature is an essential part of the regulatory circuitry. Several phages with a cI mutation preventing positive autoregulation and an up mutation in the P(RM) promoter showed near-normal behavior. We conclude that positive autoregulation is not necessary for proper operation of the lambda circuitry and speculate that it serves a partially redundant function of stabilizing a bistable circuit, a form of redundancy we term "circuit-level redundancy." We discuss our findings in the context of a two-stage model for evolution and elaboration of regulatory circuits from simpler to more complex forms.

Single-cell analysis of lambda immunity regulation

We have examined expression of the lambdacI operon in single cells via a rex Colon, two colons gfp substitution. Although average fluorescence agreed with expectations for expression of lambda-repressor, fluorescence fluctuated greatly from cell-to-cell. Fluctuations in repressor concentration are not predicted by previous models and are tolerated in part by a regulatory response to DNA damage.

Protein-protein and protein-DNA interactions of sigma70 region 4 involved in transcription activation by lambdacI

The cI protein of bacteriophage lambda (lambdacI) activates transcription from promoter P(RM) through an acidic patch on the surface of its DNA-binding domain. Genetic evidence suggests that this acidic patch stimulates transcription from P(RM) through contact with the C-terminal domain (region 4) of the sigma(70) subunit of Escherichia coli RNA polymerase. Here, we identify two basic residues in region 4 of sigma(70) that are critical for lambdacI-mediated activation of transcription from P(RM). On the basis of structural modeling, we propose that one of these sigma(70) residues, K593, facilitates the interaction between lambdacI and region 4 of sigma(70) by inducing a bend in the DNA upstream of the -35 element, whereas the other, R588, interacts directly with a critical acidic residue within the activating patch of lambdacI. Residue R588 of sigma(70) has been shown to play an important role in promoter recognition; our findings suggest that the R588 side-chain has a dual function at P(RM), facilitating the interaction of region 4 with the promoter -35 element and participating directly in the protein-protein interaction with lambdacI.

Genetic screen to dissect protein-protein interactions: ribonuclease inhibitor-ribonuclease A as a model system

Protein-protein interactions are critical for the function of biological systems. Here, we describe a means to dissect a protein-protein interaction. Our method is based on the in vivo interaction between a target protein and the peptide epitopes derived from its partner. This interaction is detected by using hybrid proteins in which the target protein and peptide epitopes are fused to the DNA-binding domain of the lambda repressor protein. An interaction prevents the transcription of a reporter gene. The efficacy of this approach is demonstrated with the ribonuclease inhibitor protein and ribonuclease A, which form a complex with an equilibrium dissociation constant in the femtomolar range. Our method can enable the identification of residues important in a designated protein-protein interaction and the development of antagonists for that interaction.

The lambda switch: cI closes the gap in autoregulation

The bacteriophage lambda genetic switch is still yielding surprises. A recent study reveals that a long-range interaction involving proteins bound 2.4 kilobases away from one another on the phage genome mediates negative autoregulation, solving a long-standing puzzle concerning the regulation of lysogeny.

Crystal structure of the lambda repressor C-terminal domain octamer

The three-dimensional structure of the lambda repressor C-terminal domain (CTD) has been determined at atomic resolution. In the crystal, the CTD forms a 2-fold symmetric tetramer that mediates cooperative binding of two repressor dimers to pairs of operator sites. Based upon this structure, a model was proposed for the structure of an octameric repressor that forms both in the presence and absence of DNA. Here, we have determined the structure of the lambda repressor CTD in three new crystal forms, under a wide variety of conditions. All crystals have essentially the same tetramer, confirming the results of the earlier study. One crystal form has two tetramers bound to form an octamer, which has the same overall architecture as the previously proposed model. An unexpected feature of the octamer in the crystal structure is a unique interaction at the tetramer-tetramer interface, formed by residues Gln209, Tyr210 and Pro211, which contact symmetry-equivalent residues from other subunits of the octamer. Interestingly, these residues are also located at the dimer-dimer interface, where the specific interactions are different. The structures thus indicate specific amino acid residues that, at least in principle, when altered could result in repressors that form tetramers but not octamers.

Octamerization of lambda CI repressor is needed for effective repression of P(RM) and efficient switching from lysogeny

The CI repressor of bacteriophage lambda is a model for the role of cooperativity in the efficient functioning of genetic switches. Pairs of CI dimers interact to cooperatively occupy adjacent operator sites at O(R) and at O(L). These CI tetramers repress the lytic promoters and activate transcription of the cI gene from P(RM). CI is also able to octamerize, forming a large DNA loop between O(R) and O(L), but the physiological role of this is unclear. Another puzzle is that, although a dimer of CI is able to repress P(RM) by binding to the third operator at O(R), O(R)3, this binding seems too weak to affect CI production in the lysogenic state. Here we show that repression of P(RM) at lysogenic CI concentrations is absolutely dependent on O(L), in this case 3.8 kb away. A mutant defective in this CI negative autoregulation forms a lysogen with elevated CI levels that cannot efficiently switch from lysogeny to lytic development. Our results invalidate previous evidence that Cro binding to O(R)3 is important in prophage induction. We propose the octameric CI:O(R)-O(L) complex increases the affinity of CI for O(R)3 by allowing a CI tetramer to link O(R)3 and the third operator at O(L), O(L)3.

Factor recruitment and TIF2/GRIP1 corepressor activity at a collagenase-3 response element that mediates regulation by phorbol esters and hormones

To investigate determinants of specific transcriptional regulation, we measured factor occupancy and function at a response element, col3A, associated with the collagenase-3 gene in human U2OS osteosarcoma cells; col3A confers activation by phorbol esters, and repression by glucocorticoid and thyroid hormones. The subunit composition and activity of AP-1, which binds col3A, paralleled the intracellular level of cFos, which is modulated by phorbol esters and glucocorticoids. In contrast, a similar AP-1 site at the collagenase-1 gene, not inducible in U2OS cells, was not bound by AP-1. The glucocorticoid receptor (GR) associated with col3A through protein-protein interactions with AP-1, regardless of AP-1 subunit composition, and repressed transcription. TIF2/GRIP1, reportedly a coactivator for GR and the thyroid hormone receptor (TR), was recruited to col3A and potentiated GR-mediated repression in the presence of a GR agonist but not antagonist. GRIP1 mutants deficient in GR binding and coactivator functions were also defective for corepression, and a GRIP1 fragment containing the GR-interacting region functioned as a dominant-negative for repression. In contrast, repression by TR was unaffected by GRIP1. Thus, the composition of regulatory complexes, and the biological activities of the bound factors, are dynamic and dependent on cell and response element contexts. Cofactors such as GRIP1 probably contain distinct surfaces for activation and repression that function in a context-dependent manner.

A small protein-protein interaction domain common to KlcB and global regulators KorA and TrbA of promiscuous IncP plasmids

The kor regulon of broad host-range, incompatibility group P (IncP) plasmids uses the KorA, KorB, and KorC repressors to regulate expression of genes for replication, conjugation, segregation, and host range. One operon, kilC, encodes the KorC repressor and two genes of unknown function (klcA and klcB). The predicted sequences of the 51.1 kDa KlcB protein, the 11.3 kDa KorA repressor, and another small (13.5 kDa) regulatory protein, TrbA, show a highly related 35 amino acid residue segment (V-L-P domain). We found that induction of the klcB gene is toxic to Escherichia coli host cells harboring an IncP plasmid. We confirmed a model in which the V-L-P domain of KlcB interacts directly with the V-L-P domain of KorA to derepress KorA-regulated operons, thereby allowing expression of toxic genes. First, a lacZ reporter fused to the kleA promoter, which is regulated by KorA and KorC, revealed that klcB induction specifically releases KorA-repression but has no effect on KorC repression. Second, induced expression of the V-L-P domains from KorA or KlcB is sufficient to release KorA-repression at the kleA promoter. Third, purified GST-KlcB fusion protein interacts specifically with His-tagged KorA. Fourth, fusion of the V-L-P domains of KorA and TrbA and full-length KlcB polypeptide to the DNA-binding domain of bacteriophage lambda repressor leads to the formation of functional, dimeric repressors, and mutations that alter conserved residues of the V-L-P domain adversely affect dimerization. Fifth, crosslinking experiments demonstrated that the V-L-P domain of KorA is able to dimerize in solution and form heterodimers in mixtures with full-length KorA polypeptide. These findings show that the V-L-P domain is a protein-protein interaction module that is likely to be responsible for dimerization of KorA and TrbA, and important for KlcB dimerization. We speculate on the possible significance of KlcB-KorA heterodimers in IncP plasmid maintenance.

Symbolic dynamics and computation in model gene networks

We analyze a class of ordinary differential equations representing a simplified model of a genetic network. In this network, the model genes control the production rates of other genes by a logical function. The dynamics in these equations are represented by a directed graph on an n-dimensional hypercube (n-cube) in which each edge is directed in a unique orientation. The vertices of the n-cube correspond to orthants of state space, and the edges correspond to boundaries between adjacent orthants. The dynamics in these equations can be represented symbolically. Starting from a point on the boundary between neighboring orthants, the equation is integrated until the boundary is crossed for a second time. Each different cycle, corresponding to a different sequence of orthants that are traversed during the integration of the equation always starting on a boundary and ending the first time that same boundary is reached, generates a different letter of the alphabet. A word consists of a sequence of letters corresponding to a possible sequence of orthants that arise from integration of the equation starting and ending on the same boundary. The union of the words defines the language. Letters and words correspond to analytically computable Poincare maps of the equation. This formalism allows us to define bifurcations of chaotic dynamics of the differential equation that correspond to changes in the associated language. Qualitative knowledge about the dynamics found by integrating the equation can be used to help solve the inverse problem of determining the underlying network generating the dynamics. This work places the study of dynamics in genetic networks in a context comprising both nonlinear dynamics and the theory of computation. (c) 2001 American Institute of Physics.

Crystal structure of the lambda repressor C-terminal domain provides a model for cooperative operator binding

Interactions between transcription factors bound to separate operator sites commonly play an important role in gene regulation by mediating cooperative binding to the DNA. However, few detailed structural models for understanding the molecular basis of such cooperativity are available. The c1 repressor of bacteriophage lambda is a classic example of a protein that binds to its operator sites cooperatively. The C-terminal domain of the repressor mediates dimerization as well as a dimer-dimer interaction that results in the cooperative binding of two repressor dimers to adjacent operator sites. Here, we present the x-ray crystal structure of the lambda repressor C-terminal domain determined by multiwavelength anomalous diffraction. Remarkably, the interactions that mediate cooperativity are captured in the crystal, where two dimers associate about a 2-fold axis of symmetry. Based on the structure and previous genetic and biochemical data, we present a model for the cooperative binding of two lambda repressor dimers at adjacent operator sites.

Transcriptional activation by artificial recruitment in yeast is influenced by promoter architecture and downstream sequences

The idea that recruitment of the transcriptional machinery to a promoter suffices for gene activation is based partly on the results of "artificial recruitment" experiments performed in vivo. Artificial recruitment can be effected by a "nonclassical" activator comprising a DNA-binding domain fused to a component of the transcriptional machinery. Here we show that activation by artificial recruitment in yeast can be sensitive to any of three factors: position of the activator-binding elements, sequence of the promoter, and coding sequences downstream of the promoter. In contrast, classical activators worked efficiently at all promoters tested. In all cases the "artificial recruitment" fusions synergized well with classical activators. A classical activator evidently differs from a nonclassical activator in that the former can touch multiple sites on the transcriptional machinery, and we propose that that difference accounts for the broader spectrum of activity of the typical classical activator. A similar conclusion is reached from studies in mammalian cells in the accompanying paper [Nevado, J., Gaudreau, L., Adam, M. & Ptashne, M. (1999) Proc. Natl. Acad. Sci. USA 96, 2674-2677].

Imposing specificity by localization: mechanism and evolvability

Cells detect extracellular signals by allostery and then give those signals meaning by 'regulated localization'. We suggest that this formulation applies to many biological processes and is particularly well illustrated by the mechanisms of gene regulation. Analysis of these mechanisms reveals that regulated localization requires simple molecular interactions that are readily used combinatorially. This system of regulation is highly 'evolvable', and its use accounts, at least in part, for the nature of the complexities observed in biological systems.

Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells

Fluctuations in rates of gene expression can produce highly erratic time patterns of protein production in individual cells and wide diversity in instantaneous protein concentrations across cell populations. When two independently produced regulatory proteins acting at low cellular concentrations competitively control a switch point in a pathway, stochastic variations in their concentrations can produce probabilistic pathway selection, so that an initially homogeneous cell population partitions into distinct phenotypic subpopulations. Many pathogenic organisms, for example, use this mechanism to randomly switch surface features to evade host responses. This coupling between molecular-level fluctuations and macroscopic phenotype selection is analyzed using the phage lambda lysis-lysogeny decision circuit as a model system. The fraction of infected cells selecting the lysogenic pathway at different phage:cell ratios, predicted using a molecular-level stochastic kinetic model of the genetic regulatory circuit, is consistent with experimental observations. The kinetic model of the decision circuit uses the stochastic formulation of chemical kinetics, stochastic mechanisms of gene expression, and a statistical-thermodynamic model of promoter regulation. Conventional deterministic kinetics cannot be used to predict statistics of regulatory systems that produce probabilistic outcomes. Rather, a stochastic kinetic analysis must be used to predict statistics of regulatory outcomes for such stochastically regulated systems.

Simulation of prokaryotic genetic circuits

Biochemical and genetic approaches have identified the molecular mechanisms of many genetic reactions, particularly in bacteria. Now a comparably detailed understanding is needed of how groupings of genes and related protein reactions interact to orchestrate cellular functions over the cell cycle, to implement preprogrammed cellular development, or to dynamically change a cell's processes and structures in response to environmental signals. Simulations using realistic, molecular-level models of genetic mechanisms and of signal transduction networks are needed to analyze dynamic behavior of multigene systems, to predict behavior of mutant circuits, and to identify the design principles applicable to design of genetic regulatory circuits. When the underlying design rules for regulatory circuits are understood, it will be far easier to recognize common circuit motifs, to identify functions of individual proteins in regulation, and to redesign circuits for altered functions.

Activation of transcription in vitro by recruitment of the yeast RNA polymerase II holoenzyme

It has been argued that many transcriptional activators work by "recruitment," that is, by helping the transcriptional machinery bind stably to DNA. We demonstrate here a realization of a strong prediction of this idea in an in vitro transcription reaction performed with purified yeast RNA polymerase II holoenzyme and a classical transcriptional activator. We show that the level of transcription reached by the activator working on low concentrations of holoenzyme can also be reached in the absence of activator by raising the holoenzyme concentration, and that under that condition the activator has no further stimulatory effect. We also show, in agreement with another prediction of the recruitment model, that in a reaction using a holoenzyme purified from cells bearing the "P" mutation, transcription is stimulated by a DNA-tethered peptide that binds the mutant holoenzyme component Gal11P but that lacks a classical activating region.

Transcriptional activation by recruitment

The recruitment model for gene activation stipulates that an activator works by bringing the transcriptional machinery to the DNA. Recent experiments in bacteria and yeast indicate that many genes can be activated by this mechanism. These findings have implications for our understanding of the nature of activating regions and their targets, and for the role of histones in gene regulation.

Intracellular receptors use a common mechanism to interpret signaling information at response elements

The glucocorticoid receptor (GR) activates transcription in certain glucocorticoid response element (GRE) contexts, and represses or displays no activity in others. We isolated point mutations in one GRE, plfG, at which GR activated transcription under conditions in which the wild-type element was inactive or conferred repression, implying that GREs may carry signals that are interpreted by bound receptors. Consistent with this notion, we identified a mutant rat GR, K461A, which activated transcription in all GRE contexts tested, implying that this residue is important in interpretation of GRE signals. In a yeast screen of 60,000 GR mutants for strong activation from plfG, all 13 mutants isolated contained substitutions at K461. This lysine residue is highly conserved in the zinc-binding region (ZBR) of the intracellular receptor (IR) superfamily; when it was mutated in MR and RARbeta, the resulting receptors similarly activated transcription at response elements that their wild-type counterparts repressed or were inactive. We suggest that IR response elements serve in part as signaling components, and that a critical lysine residue serves as an allosteric "lock" that restricts IRs to inactive or repressing configurations except in response element contexts that signal their conversion to transcriptional activators. Therefore, mutation of this residue produces altered receptors that activate in many or all response element contexts.

Circuit simulation of genetic networks

Genetic networks with tens to hundreds of genes are difficult to analyze with currently available techniques. Because of the many parallels in the function of these biochemically based genetic circuits and electrical circuits, a hybrid modeling approach is proposed that integrates conventional biochemical kinetic modeling within the framework of a circuit simulation. The circuit diagram of the bacteriophage lambda lysislysogeny decision circuit represents connectivity in signal paths of the biochemical components. A key feature of the lambda genetic circuit is that operons function as active integrated logic components and introduce signal time delays essential for the in vivo behavior of phage lambda.

Synergistic activation of transcription by bacteriophage lambda cI protein and E. coli cAMP receptor protein

Two heterologous prokaryotic activators, the bacteriophage lambda cI protein (lambda cI) and the Escherichia coli cyclic AMP receptor protein (CRP), were shown to activate transcription synergistically from an artificial promoter bearing binding sites for both proteins. The synergy depends on a functional activation (positive control) surface on each activator. These results imply that both proteins interact directly with RNA polymerase and thus suggest a precise mechanism for transcriptional synergy: the interaction of two activators with two distinct surfaces of RNA polymerase.

Target of the transcriptional activation function of phage lambda cI protein

Activation of transcription initiation by the cI protein of phage lambda is thought to be mediated by a direct interaction between cl and RNA polymerase at the PRM promoter. Two negatively charged amino acid residues in the DNA binding domain of cI play a key role in activation, suggesting that these residues contact RNA polymerase. The subunit of RNA polymerase involved was identified by selecting polymerase mutants that restored the activation function of a mutant form of cI protein. Although previous studies suggest that several activators interact with the alpha subunit of RNA polymerase, the results here suggest that cI interacts with the sigma subunit. An arginine to histidine change near the carboxyl terminus of sigma specifically suppresses an aspartic acid to asparagine change in the activation region of cI. This finding supports the direct-contact model and suggests that a cluster of positively charged residues near the carboxyl terminus of sigma is the target of the negatively charged activation region of cI.

The Cro-like Apl repressor of coliphage 186 is required for prophage excision and binds near the phage attachment site

The Apl protein of the temperature coliphage 186 represses transcription of the immunity repressor gene and down-regulates lytic transcription. It is shown here that an apl- mutant is competent for lytic development and establishes lysogeny normally but is defective in excision of the prophage. The Apl protein binds between the lytic and lysogenic promoters and also near the phage attachment site, suggesting that its role in excision is direct. Apl thus appears to act as an excisionase as well as a repressor. The pattern of Apl-induced DNase I enhancements indicates that the DNA is bent by Apl. Potential Apl recognition sequences are identified; these sequences are directly repeated several times across each binding region and are spaced 10 or 11 bases apart, suggesting that Apl binds to one face of the DNA helix.