DNA looping can enhance lysogenic CI transcription in phage lambda

Bacterial resistance to temperate phage is influenced by the frequency of lysogenic establishment

Temperate phages can shape bacterial community dynamics and evolution through lytic and lysogenic life cycles. In response, bacteria that resist phage infection can emerge. This study explores phage-based factors that influence bacterial resistance using a model system of temperate P22 phage and Salmonella both inside and outside the mammalian host. Phages that remained functional despite gene deletions had minimal impact on lysogeny and phage resistance except for deletions in the immI region that substantially reduced lysogeny and increased phage resistance to levels comparable to that observed with an obligately lytic P22. This immI deletion does not make the lysogen less competitive but instead increases the frequency of bacterial lysis. Thus, subtle changes in the balance between lysis and lysogeny during the initial stages of infection can significantly influence the extent of phage resistance in the bacterial population. Our work highlights the complex nature of the phage-bacteria-mammalian host triad.

Improving cooperativity of transcription activators by oligomerization domains in mammalian cells

Cooperative activation is critical for the applications of synthetic biology in mammalian cells. In this study, we have developed cooperative transcription factor by fusing oligomerization domain in mammalian cells. Firstly, we demonstrated that two oligomerized domains (CI434 and CI) successfully improved transcription factor cooperativity in bacterial cells but failed to increase cooperativity in mammalian cells, possibly because the additional mammalian activation domain disrupted their oligomerization capability. Therefore, we chose a different type of oligomerized domain (CarHC), whose ability to oligomerize is not dependent on its C-terminal domains, to fuse with a transcription factor (RpaR) and activation domain (VTR3), forming a potential cooperative transcription activator RpaR-CarH-VTR3 for mammalian regulatory systems. Compared with RpaR-VTR3, the cooperativity of RpaR-CarH-VTR3 was significantly improved with higher Hill coefficient and a narrower input range in the inducible switch system in mammalian cells. Moreover, a mathematical model based on statistical mechanics model was developed and the simulation results supported the hypothesis that the tetramer of the CarH domain in mammalian cells was the reason for the cooperative capacity of RpaR-CarH-VTR3.

Robustness of DNA looping across multiple cell divisions in individual bacteria

DNA looping has emerged as a central paradigm of transcriptional regulation, as it is shared across many living systems. One core property of DNA looping-based regulation is its ability to greatly enhance repression or activation of genes with only a few copies of transcriptional regulators. However, this property based on a small number of proteins raises the question of the robustness of such a mechanism with respect to the large intracellular perturbations taking place during growth and division of the cell. Here we address the issue of sensitivity to variations of intracellular parameters of gene regulation by DNA looping. We use the lac system as a prototype to experimentally identify the key features of the robustness of DNA looping in growing Escherichia coli cells. Surprisingly, we observe time intervals of tight repression spanning across division events, which can sometimes exceed 10 generations. Remarkably, the distribution of such long time intervals exhibits memoryless statistics that is mostly insensitive to repressor concentration, cell division events, and the number of distinct loops accessible to the system. By contrast, gene regulation becomes highly sensitive to these perturbations when DNA looping is absent. Using stochastic simulations, we propose that the observed robustness to division emerges from the competition between fast, multiple rebinding events of repressors and slow initiation rate of the RNA polymerase. We argue that fast rebinding events are a direct consequence of DNA looping that ensures robust gene repression across a range of intracellular perturbations.

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.

Ageing and sources of transcriptional heterogeneity

Cellular heterogeneity is an important contributor to biological function and is employed by cells, tissues and organisms to adapt, compensate, respond, defend and/or regulate specific processes. Research over the last decades has revealed that transcriptional noise is a major driver for cell-to-cell variability. In this review we will discuss sources of transcriptional variability, in particular bursting of gene expression and how it could contribute to cellular states and fate decisions. We will highlight recent developments in single cell sequencing technologies that make it possible to address cellular heterogeneity in unprecedented detail. Finally, we will review recent literature, in which these new technologies are harnessed to address pressing questions in the field of ageing research, such as transcriptional noise and cellular heterogeneity in the course of ageing.

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.

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.

Cell fate decisions emerge as phages cooperate or compete inside their host

The system of the bacterium Escherichia coli and its virus, bacteriophage lambda, is paradigmatic for gene regulation in cell-fate development, yet insight about its mechanisms and complexities are limited due to insufficient resolution of study. Here we develop a 4-colour fluorescence reporter system at the single-virus level, combined with computational models to unravel both the interactions between phages and how individual phages determine cellular fates. We find that phages cooperate during lysogenization, compete among each other during lysis, and that confusion between the two pathways occasionally occurs. Additionally, we observe that phage DNAs have fluctuating cellular arrival times and vie for resources to replicate, enabling the interplay during different developmental paths, where each phage genome may make an individual decision. These varied strategies could separate the selection for replication-optimizing beneficial mutations during lysis from sequence diversification during lysogeny, allowing rapid adaptation of phage populations for various environments.

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

Supercoiling biases the formation of loops involved in gene regulation

The function of DNA as a repository of genetic information is well-known. The post-genomic effort is to understand how this information-containing filament is chaperoned to manage its compaction and topological states. Indeed, the activities of enzymes that transcribe, replicate, or repair DNA are regulated to a large degree by access. Proteins that act at a distance along the filament by binding at one site and contacting another site, perhaps as part of a bigger complex, create loops that constitute topological domains and influence regulation. DNA loops and plectonemes are not necessarily spontaneous, especially large loops under tension for which high energy is required to bring their ends together, or small loops that require accessory proteins to facilitate DNA bending. However, the torsion in stiff filaments such as DNA dramatically modulates the topology, driving it from extended and genetically accessible to more looped and compact, genetically secured forms. Furthermore, there are accessory factors that bias the response of the DNA filament to supercoiling. For example, small molecules like polyamines, which neutralize the negative charge repulsions along the phosphate backbone, enhance flexibility and promote writhe over twist in response to torsion. Such increased flexibility likely pushes the topological equilibrium from twist toward writhe at tensions thought to exist in vivo. A predictable corollary is that stiffening DNA antagonizes looping and bending. Certain sequences are known to be more or less flexible or to exhibit curvature, and this may affect interactions with binding proteins. In vivo all of these factors operate simultaneously on DNA that is generally negatively supercoiled to some degree. Therefore, in order to better understand gene regulation that involves protein-mediated DNA loops, it is critical to understand the thermodynamics and kinetics of looping in DNA that is under tension, negatively supercoiled, and perhaps exposed to molecules that alter elasticity. Recent experiments quantitatively reveal how much negatively supercoiling DNA lowers the free energy of looping, possibly biasing the operation of genetic switches.

Effect of Interaction between Chromatin Loops on Cell-to-Cell Variability in Gene Expression

According to recent experimental evidence, the interaction between chromatin loops, which can be characterized by three factors-connection pattern, distance between regulatory elements, and communication form, play an important role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effect that addresses the question of how changes in these factors affect variability at the expression level in a systematic rather than case-by-case fashion. Here we make such an effort, based on a mechanic model that maps three fundamental patterns for two interacting DNA loops into a 4-state model of stochastic transcription. We first show that in contrast to side-by-side loops, nested loops enhance mRNA expression and reduce expression noise whereas alternating loops have just opposite effects. Then, we compare effects of facilitated tracking and direct looping on gene expression. We find that the former performs better than the latter in controlling mean expression and in tuning expression noise, but this control or tuning is distance-dependent, remarkable for moderate loop lengths, and there is a limit loop length such that the difference in effect between two communication forms almost disappears. Our analysis and results justify the facilitated chromatin-looping hypothesis.

Measurement of gene regulation in individual cells reveals rapid switching between promoter states

In vivo mapping of transcription-factor binding to the transcriptional output of the regulated gene is hindered by probabilistic promoter occupancy, the presence of multiple gene copies, and cell-to-cell variability. We demonstrate how to overcome these obstacles in the lysogeny maintenance promoter of bacteriophage lambda, P(RM). We simultaneously measured the concentration of the lambda repressor CI and the number of messenger RNAs (mRNAs) from P(RM) in individual Escherichia coli cells, and used a theoretical model to identify the stochastic activity corresponding to different CI binding configurations. We found that switching between promoter configurations is faster than mRNA lifetime and that individual gene copies within the same cell act independently. The simultaneous quantification of transcription factor and promoter activity, followed by stochastic theoretical analysis, provides a tool that can be applied to other genetic circuits.

Biological Sources of Intrinsic and Extrinsic Noise in cI Expression of Lysogenic Phage Lambda

Genetically identical cells exposed to homogeneous environment can show remarkable phenotypic difference. To predict how phenotype is shaped, understanding of how each factor contributes is required. During gene expression processes, noise could arise either intrinsically in biochemical processes of gene expression or extrinsically from other cellular processes such as cell growth. In this work, important noise sources in gene expression of phage λ lysogen are quantified using models described by stochastic differential equations (SDEs). Results show that DNA looping has sophisticated impacts on gene expression noise: When DNA looping provides autorepression, like in wild type, it reduces noise in the system; When the autorepression is defected as it is in certain mutants, DNA looping increases expression noise. We also study how each gene operator affects the expression noise by changing the binding affinity between the gene and the transcription factor systematically. We find that the system shows extraordinarily large noise when the binding affinity is in certain range, which changes the system from monostable to bistable. In addition, we find that cell growth causes non-negligible noise, which increases with gene expression level. Quantification of noise and identification of new noise sources will provide deeper understanding on how stochasticity impacts phenotype.

Cooperativity leads to temporally-correlated fluctuations in the bacteriophage lambda genetic switch

Cooperative interactions are widespread in biochemical networks, providing the nonlinear response that underlies behavior such as ultrasensitivity and robust switching. We introduce a temporal correlation function-the conditional activity-to study the behavior of these phenomena. Applying it to the bistable genetic switch in bacteriophage lambda, we find that cooperative binding between binding sites on the prophage DNA lead to non-Markovian behavior, as quantified by the conditional activity. Previously, the conditional activity has been used to predict allosteric pathways in proteins; here, we show that it identifies the rare unbinding events which underlie induction from lysogeny to lysis.

Revisiting bistability in the lysis/lysogeny circuit of bacteriophage lambda

The lysis/lysogeny switch of bacteriophage lambda serves as a paradigm for binary cell fate decision, long-term maintenance of cellular state and stimulus-triggered switching between states. In the literature, the system is often referred to as “bistable.” However, it remains unclear whether this term provides an accurate description or is instead a misnomer. Here we address this question directly. We first quantify transcriptional regulation governing lysogenic maintenance using a single-cell fluorescence reporter. We then use the single-cell data to derive a stochastic theoretical model for the underlying regulatory network. We use the model to predict the steady states of the system and then validate these predictions experimentally. Specifically, a regime of bistability, and the resulting hysteretic behavior, are observed. Beyond the steady states, the theoretical model successfully predicts the kinetics of switching from lysogeny to lysis. Our results show how the physics-inspired concept of bistability can be reliably used to describe cellular phenotype, and how an experimentally-calibrated theoretical model can have accurate predictive power for cell-state switching.

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.

DNA supercoiling enhances cooperativity and efficiency of an epigenetic switch

Bacteriophage λ stably maintains its dormant prophage state but efficiently enters lytic development in response to DNA damage. The mediator of these processes is the λ repressor protein, CI, and its interactions with λ operator DNA. This λ switch is a model on the basis of which epigenetic switch regulation is understood. Using single molecule analysis, we directly examined the stability of the CI-operator structure in its natural, supercoiled state. We marked positions adjacent to the λ operators with peptide nucleic acids and monitored their movement by tethered particle tracking. Compared with relaxed DNA, the presence of supercoils greatly enhances juxtaposition probability. Also, the efficiency and cooperativity of the λ switch is significantly increased in the supercoiled system compared with a linear assay, increasing the Hill coefficient.

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.

Role of cis-acting sites in stimulation of the phage λ P(RM) promoter by CI-mediated looping

The lysogenic state of phage λ is maintained by the CI repressor. CI binds to three operators each in the right operator (O(R)) and left operator (O(L)) regions, which lie 2.4 kb apart. At moderate CI levels, the predominant binding pattern is two dimers of CI bound cooperatively at each regulatory region. The resulting tetramers can then interact, forming an octamer and a loop of the intervening DNA. CI is expressed from the P(RM) promoter, which lies in the O(R) region and is subjected to multiple regulatory controls. Of these, the most recently discovered is stimulation by loop formation. In this work, we have investigated the mechanism by which looping stimulates P(RM). We find that two cis-acting sites lying in the O(L) region are involved. One site, an UP element, is required for stimulation. Based on the behavior of other promoters with UP elements located upstream of the -35 region, we suggest that a subunit of RNA polymerase (RNAP) bound at P(RM) binds to the UP element located in the O(L) region. In addition, adjacent to the UP element lies a binding site for integration host factor (IHF); this site plays a less critical role but is required for stimulation of the weak prm240 allele. A loop with CI at the O(L)2 and O(L)3 operators does not stimulate P(RM), while one with CI only at O(L)2 provides some stimulation. We discuss possible mechanisms for stimulation.

Are gene loops the cause of transcriptional noise?

Expression levels of the same mRNA or protein vary significantly among the cells of an otherwise identical population. Such biological noise has great functional implications and is largely due to transcriptional bursting, the episodic production of mRNAs in short, intense bursts, interspersed by periods of transcriptional inactivity. Bursting has been demonstrated in a wide range of pro- and eukaryotic species, attesting to its universal importance. However, the mechanistic origins of bursting remain elusive. A different type of phenomenon, which has also been suggested to be widespread, is the physical interaction between the promoter and 3' end of a gene. Several functional roles have been proposed for such gene loops, including the facilitation of transcriptional reinitiation. Here, I discuss the most recent findings related to these subjects and argue that gene loops are a likely cause of transcriptional bursting and, thus, biological noise.

Enhancer-like long-range transcriptional activation by λ CI-mediated DNA looping

How distant enhancer elements regulate the assembly of a transcription complex at a promoter remains poorly understood. Here, we use long-range gene regulation by the bacteriophage λ CI protein as a powerful system to examine this process in vivo. A 2.3-kb DNA loop, formed by CI bridging its binding sites at OR and OL, is known already to enhance repression at the lysogenic promoter PRM, located at OR. Here, we show that CI looping also activates PRM by allowing the C-terminal domain of the α subunit of the RNA polymerase bound at PRM to contact a DNA site adjacent to the distal CI sites at OL. Our results establish OL as a multifaceted enhancer element, able to activate transcription from long distances independently of orientation and position. We develop a physicochemical model of our in vivo data and use it to show that the observed activation is consistent with a simple recruitment mechanism, where the α-C-terminal domain to DNA contact need only provide ∼2.7 kcal/mol of additional binding energy for RNA polymerase. Structural modeling of this complete enhancer-promoter complex reveals how the contact is achieved and regulated, and suggests that distal enhancer elements, once appropriately positioned at the promoter, can function in essentially the same way as proximal promoter elements.

The effect of nonspecific binding of lambda repressor on DNA looping dynamics

The λ repressor (CI) protein-induced DNA loop maintains stable lysogeny, yet allows efficient switching to lysis. Herein, the kinetics of loop formation and breakdown has been characterized at various concentrations of protein using tethered particle microscopy and a novel, to our knowledge, method of analysis. Our results show that a broad distribution of rate constants and complex kinetics underlie loop formation and breakdown. In addition, comparison of the kinetics of looping in wild-type DNA and DNA with mutated o3 operators showed that these sites may trigger nucleation of nonspecific binding at the closure of the loop. The average activation energy calculated from the rate constant distribution is consistent with a model in which nonspecific binding of CI between the operators shortens their effective separation, thereby lowering the energy barrier for loop formation and broadening the rate constant distribution for looping. Similarly, nonspecific binding affects the kinetics of loop breakdown by increasing the number of loop-securing protein interactions, and broadens the rate constant distribution for this reaction. Therefore, simultaneous increase of the rate constant for loop formation and reduction of that for loop breakdown stabilizes lysogeny. Given these simultaneous changes, the frequency of transitions between the looped and the unlooped state remains nearly constant. Although the loop becomes more stable thermodynamically with increasing CI concentration, it still opens periodically, conferring sensitivity to environmental changes, which may require switching to lytic conditions.

Stretching DNA to quantify nonspecific protein binding

Nonspecific binding of regulatory proteins to DNA can be an important mechanism for target search and storage. This seems to be the case for the lambda repressor protein (CI), which maintains lysogeny after infection of E. coli. CI binds specifically at two distant regions along the viral genome and induces the formation of a repressive DNA loop. However, single-molecule imaging as well as thermodynamic and kinetic measurements of CI-mediated looping show that CI also binds to DNA nonspecifically and that this mode of binding may play an important role in maintaining lysogeny. This paper presents a robust phenomenological approach using a recently developed method based on the partition function, which allows calculation of the number of proteins bound nonspecific to DNA from measurements of the DNA extension as a function of applied force. This approach was used to analyze several cycles of extension and relaxation of λ DNA performed at several CI concentrations to measure the dissociation constant for nonspecific binding of CI (~100 nM), and to obtain a measurement of the induced DNA compaction (~10%) by CI.

Phage λ--new insights into regulatory circuits

Bacteriophage λ, rediscovered in the early 1950s, has served as a model in molecular biology studies for decades. Although currently more complex organisms and more complicated biological systems can be studied, this phage is still an excellent model to investigate principles of biological processes occurring at the molecular level. In fact, very few other biological models provide possibilities to examine regulations of biological mechanisms as detailed as performed with λ. In this chapter, recent advances in our understanding of mechanisms of bacteriophage λ development are summarized and discussed. Particularly, studies on (i) phage DNA injection, (ii) molecular bases of the lysis-versus-lysogenization decision and the lysogenization process itself, (iii) prophage maintenance and induction, (iv), λ DNA replication, (v) phage-encoded recombination systems, (vi) transcription antitermination, (vii) formation of the virion structure, and (viii) lysis of the host cell, as published during several past years, will be presented.

Decision making in living cells: lessons from a simple system

The life cycle of bacteriophage lambda serves as a simplified paradigm for cell-fate decisions. The ongoing quantitative, high-resolution experimental investigation of this life cycle has produced some important insights in recent years. These insights have to do with the way cells choose among alternative fates, how they maintain long-term memory of their gene-expression state, and how they switch from one stable state to another. The recent studies have highlighted the role of spatiotemporal effects in cellular processes and the importance of distinguishing chemical stochasticity from possible hidden variables in cellular decision making.

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

Lysogen stability is determined by the frequency of activity bursts from the fate-determining gene

Bacterial lysogeny serves as a simple paradigm for cell differentiation.

We characterize the activity of the fate-determining genes, cI and cro, with single-molecule resolution.

Stability of the lysogenic state is found to depend in a simple manner on the frequency of activity bursts from cI.

The ability of living cells to maintain an inheritable memory of their gene-expression state is key to cellular differentiation. Bacterial lysogeny serves as a simple paradigm for long-term cellular memory. In this study, we address the following question: in the absence of external perturbation, how long will a cell stay in the lysogenic state before spontaneously switching away from that state? We show by direct measurement that lysogen stability exhibits a simple exponential dependence on the frequency of activity bursts from the fate-determining gene, cI. We quantify these gene-activity bursts using single-molecule-resolution mRNA measurements in individual cells, analyzed using a stochastic mathematical model of the gene-network kinetics. The quantitative relation between stability and gene activity is independent of the fine details of gene regulation, suggesting that a quantitative prediction of cell-state stability may also be possible in more complex systems.

Evolution of complex gene regulatory circuits by addition of refinements

How do complex gene regulatory circuits evolve? These circuits involve many interacting components, which work together to specify patterns of gene expression. They typically include many subtle mechanistic features, but in most cases it is unclear whether these features are essential for the circuit to work at all, or if instead they make a functional circuit work better. In the latter case, such a feature is here termed 'dispensable', and it is plausible that the feature has been added at a late stage in the evolution of the circuit. This review describes experimental tests of this question, using the phage λ gene regulatory circuit. Several features of this circuit are found to be dispensable, in the sense that the circuitry works without these features, though not as well as the wild type. In some cases, second-site suppressor mutations are needed to confer near-normal behavior in the absence of such a feature. These findings are discussed here in the context of a two-stage model for evolution of gene regulatory circuits. In this model, a circuit evolves by assembly of a primitive or basic form, followed by adjustment of parameters and addition of qualitatively new features. Pathways are suggested for the addition of such features to a more basic form. Selected examples in other systems are described. Some of the dispensable features of phage λ may be evolutionary refinements. Finding that a feature is dispensable, however, does not prove that it is a late addition - it is possible that it was essential early in evolution, and became dispensable as the circuit evolved. Conversely, a late addition might have become essential. As ongoing work provides additional examples of dispensable features, it may become clearer how often they represent refinements.

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.

Stability and instability in the lysogenic state of phage lambda

Complex gene regulatory circuits exhibit emergent properties that are difficult to predict from the behavior of the components. One such property is the stability of regulatory states. Here we analyze the stability of the lysogenic state of phage λ. In this state, the virus maintains a stable association with the host, and the lytic functions of the virus are repressed by the viral CI repressor. This state readily switches to the lytic pathway when the host SOS system is induced. A low level of SOS-dependent switching occurs without an overt stimulus. We found that the intrinsic rate of switching to the lytic pathway, measured in a host lacking the SOS response, was almost undetectably low, probably less than 10(-8)/generation. We surmise that this low rate has not been selected directly during evolution but results from optimizing the rate of switching in a wild-type host over the natural range of SOS-inducing conditions. We also analyzed a mutant, λprm240, in which the promoter controlling CI expression was weakened, rendering lysogens unstable. Strikingly, the intrinsic stability of λprm240 lysogens depended markedly on the growth conditions; lysogens grown in minimal medium were nearly stable but switched at high rates when grown in rich medium. These effects on stability likely reflect corresponding effects on the strength of the prm240 promoter, measured in an uncoupled assay system. Several derivatives of λprm240 with altered stabilities were characterized. This mutant and its derivatives afford a model system for further analysis of stability.

Determination of the number of proteins bound non-specifically to DNA

We have determined the change in the number of proteins bound non-specifically to DNA as a function of applied force using force-extension measurements on tethered DNA. Using magnetic tweezers, single molecules of λ DNA were repeatedly stretched and relaxed in the absence and presence of 170 nM λ repressor protein (CI). CI binds to six specific sites of λ DNA with nanomolar affinity and also binds non-specifically with micromolar affinity. The force versus extension data were analyzed using a recently developed theoretical framework for quantitative determination of protein binding to the DNA. The results indicate that the non-specific binding of CI changes the force-extension relation significantly in comparison to that of naked DNA. The DNA tether used in our experiment would have about 640 bound repressors, if it was completely saturated with bound proteins. We find that as the pulling force on DNA is reduced from 4.81 to 0.13 pN, approximately 138 proteins bind to DNA, which is about 22% of the length of the tethered DNA. Our results show that 0.13 pN is not low enough to cause saturation of DNA by repressor and 4.81 pN is also not high enough to eliminate all the repressors bound to DNA. This demonstrates that the force-extension relation provides an effective approach for estimating the number of proteins bound non-specifically to a DNA molecule.

Predicting gene-regulation functions: lessons from temperate bacteriophages

Gene-regulation functions (GRF) provide a unique characteristic of a cis-regulatory module (CRM), relating the concentrations of transcription factors (input) to the promoter activities (output). The challenge is to predict GRFs from the sequence. Here we systematically consider the lysogeny-lysis CRMs of different temperate bacteriophages such as the Lactobacillus casei phage A2, Escherichia coli phages lambda, and 186 and Lactococcal phage TP901-1. This study allowed explaining a recent experimental puzzle on the role of Cro protein in the lambda switch. Several general conclusions have been drawn: 1), long-range interactions, multilayer assembly and DNA looping may lead to complex GRFs that cannot be described by linear functions of binding site occupancies; 2), in general, GRFs cannot be described by the Boolean logic, whereas a three-state non-Boolean logic suffices for the studied examples; 3), studied CRMs of the intact phages seemed to have a similar GRF topology (the number of plateaus and peaks corresponding to different expression regimes); we hypothesize that functionally equivalent CRMs might have topologically equivalent GRFs for a larger class of genetic systems; and 4) within a given GRF class, a set of mechanistic-to-mathematical transformations has been identified, which allows shaping the GRF before carrying out a system-level analysis.

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.

Single-molecule approaches to probe the structure, kinetics, and thermodynamics of nucleoprotein complexes that regulate transcription

Single-molecule experimentation has contributed significantly to our understanding of the mechanics of nucleoprotein complexes that regulate epigenetic switches. In this minireview, we will discuss the application of the tethered-particle motion technique, magnetic tweezers, and atomic force microscopy to (i) directly visualize and thermodynamically characterize DNA loops induced by the lac, gal, and lambda repressors and (ii) understand the mechanistic role of DNA-supercoiling and DNA-bending cofactors in both prokaryotic and eukaryotic systems.

Quantitative analysis of DNA-looping kinetics from tethered particle motion experiments

In this chapter we show the application of a maximum-likelihood-based method to the reconstruction of DNA-looping single-molecule time traces from tethered particle motion experiments. The method does not require time filtering of the data and improves the time resolution by an order of magnitude with respect to the threshold-crossing approach. Moreover, it is not based on presumed kinetic models, overcoming the limitations of other approaches proposed previously, and allowing its applications to mechanisms with complex kinetic schemes. Numerical simulations have been used to test the performances of this analysis over a wide range of time scales. We have then applied this method to determine the looping kinetics of a well-known DNA-looping protein, the lambda-repressor.

Qualitatively distinct phenotypes in the design space of biochemical systems

Although characterization of the genotype has undergone revolutionary advances as a result of the successful genome projects, the chasm between our understanding of a fully characterized gene sequence and the phenotypic repertoire of the organism is as broad and deep as it was in the pre-genomic era. There are two fundamental unsolved problems that provide the context for the challenges in relating genotype to phenotype. We address one of these and describe a generic method for constructing a system design space in which qualitatively distinct phenotypes can be identified and counted, their relative fitness analyzed and compared, and their tolerance to change measured.

Cytoplasmic membrane protonmotive force energizes periplasmic interactions between ExbD and TonB

The TonB system of Escherichia coli (TonB/ExbB/ExbD) transduces the protonmotive force (pmf) of the cytoplasmic membrane to drive active transport by high-affinity outer membrane transporters. In this study, chromosomally encoded ExbD formed formaldehyde-linked complexes with TonB, ExbB and itself (homodimers) in vivo. Pmf was required for detectable cross-linking between TonB-ExbD periplasmic domains. Consistent with that observation, the presence of inactivating transmembrane domain mutations ExbD(D25N) or TonB(H20A) also prevented efficient formaldehyde cross-linking between ExbD and TonB. A specific site of periplasmic interaction occurred between ExbD(A92C) and TonB(A150C) and required functional transmembrane domains in both proteins. Conversely, neither TonB, ExbB nor pmf were required for ExbD dimer formation. These data suggest two possible models where either dynamic complex formation occurred through transmembrane domains or the transmembrane domains of ExbD and TonB configure their respective periplasmic domains. Analysis of T7-tagged ExbD with anti-ExbD antibodies revealed that a T7 tag was responsible both for our previous failure to detect T7-ExbD-ExbB and T7-ExbD-TonB formaldehyde-linked complexes and for the concomitant artefactual appearance of T7-ExbD trimers.

A simplified model for lysogenic regulation through DNA looping

The lysogenic state of bacteriophage lambda has been a key model for understanding gene regulation. A single protein, CI repressor, maintains this state by blocking gene expression from two promoters separated by approximately 2.3 kbp of DNA. CI controls its own expression by positive and negative feedback by binding to OR to regulate the PRM promoter. Not only does CI interact directly with operator DNA, but CI tetramers bound to OL and OR can form an octamer, looping the DNA that lies between them. Previous studies show that OL can assist with negative regulation of PRM, and we recently showed that DNA looping can also enhance looping activation. In this paper we present a new interpretation of our recent experimental data based a new crystal structure of CI. The simpler model suggested by the new structural information predicts that the most common forms of the DNA loop have similar activation behavior, and are enhanced approximately 2.2-fold relative to unlooped activation.

DNA looping provides stability and robustness to the bacteriophage lambda switch

The bistable gene regulatory switch controlling the transition from lysogeny to lysis in bacteriophage lambda presents a unique challenge to quantitative modeling. Despite extensive characterization of this regulatory network, the origin of the extreme stability of the lysogenic state remains unclear. We have constructed a stochastic model for this switch. Using Forward Flux Sampling simulations, we show that this model predicts an extremely low rate of spontaneous prophage induction in a recA mutant, in agreement with experimental observations. In our model, the DNA loop formed by octamerization of CI bound to the O(L) and O(R) operator regions is crucial for stability, allowing the lysogenic state to remain stable even when a large fraction of the total CI is depleted by nonspecific binding to genomic DNA. DNA looping also ensures that the switch is robust to mutations in the order of the O(R) binding sites. Our results suggest that DNA looping can provide a mechanism to maintain a stable lysogenic state in the face of a range of challenges including noisy gene expression, nonspecific DNA binding, and operator site mutations.

Direct demonstration and quantification of long-range DNA looping by the lambda bacteriophage repressor

Recently, it was proposed that DNA looping by the λ repressor (CI protein) strengthens repression of lytic genes during lysogeny and simultaneously ensures efficient switching to lysis. To investigate this hypothesis, tethered particle motion experiments were performed and dynamic CI-mediated looping of single DNA molecules containing the λ repressor binding sites separated by 2317 bp (the wild-type distance) was quantitatively analyzed. DNA containing all three intact operators or with mutated o3 operators were compared. Modeling the thermodynamic data established the free energy of CI octamer-mediated loop formation as 1.7 kcal/mol, which decreased to –0.7 kcal/mol when supplemented by a tetramer (octamer+tetramer-mediated loop). These results support the idea that loops secured by an octamer of CI bound at oL1, oL2, oR1 and oR2 operators must be augmented by a tetramer of CI bound at the oL3 and oR3 to be spontaneous and stable. Thus the o3 sites are critical for loops secured by the CI protein that attenuate cI expression.

Quantitative transcription factor binding kinetics at the single-molecule level

We investigated the binding interaction between the bacteriophage lambda-repressor CI and its target DNA using total internal reflection fluorescence microscopy. Large stepwise changes in the intensity of the red fluorescent protein fused to CI were observed as it associated with and dissociated from individually labeled single-molecule DNA targets. The stochastic association and dissociation were characterized by Poisson statistics. Dark and bright intervals were measured for thousands of individual events. The exponential distribution of the intervals allowed direct determination of the association and dissociation rate constants (k(a) and k(d), respectively). We resolved in detail how k(a) and k(d) varied as a function of three control parameters: the DNA length L, the CI dimer concentration, and the binding affinity. Our results show that although interactions with nonoperator DNA sequences are observable, CI binding to the operator site is not dependent on the length of flanking nonoperator DNA.

Signatures of combinatorial regulation in intrinsic biological noise

Gene expression is controlled by the action of transcription factors that bind to DNA and influence the rate at which a gene is transcribed. The quantitative mapping between the regulator concentrations and the output of the gene is known as the cis-regulatory input function (CRIF). Here, we show how the CRIF shapes the form of the joint probability distribution of molecular copy numbers of the regulators and the product of a gene. Namely, we derive a class of fluctuation-based relations that relate the moments of the distribution to the derivatives of the CRIF. These relations are useful because they enable statistics of naturally arising cell-to-cell variations in molecular copy numbers to substitute for traditional manipulations for probing regulatory mechanisms. We demonstrate that these relations can distinguish super- and subadditive gene regulatory scenarios (molecular analogs of AND and OR logic operations) in simulations that faithfully represent bacterial gene expression. Applications and extensions to other regulatory scenarios are discussed.