Translation initiation controls the relative rates of expression of the bacteriophage lambda late genes

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Genetic toggle switch controlled by bacterial growth rate

Background

In favorable conditions bacterial doubling time is less than 20 min, shorter than DNA replication time. In E. coli a single round of genome replication lasts about 40 min and it must be accomplished about 20 min before cell division. To achieve such fast growth rates bacteria perform multiple replication rounds simultaneously. As a result, when the division time is as short as 20 min E. coli has about 8 copies of origin of replication (ori) and the average copy number of the genes situated close to ori can be 4 times larger than those near the terminus of replication (ter). It implies that shortening of cell cycle may influence dynamics of regulatory pathways involving genes placed at distant loci.

Results

We analyze this effect in a model of a genetic toggle switch, i.e. a system of two mutually repressing genes, one localized in the vicinity of ori and the other localized in the vicinity of ter. Using a stochastic model that accounts for cell growth and divisions we demonstrate that shortening of the cell cycle can induce switching of the toggle to the state in which expression of the gene placed near ter is suppressed. The toggle bistability causes that the ratio of expression of the competing genes changes more than two orders of magnitude for a two-fold change of the doubling time. The increasing stability of the two toggle states enhances system sensitivity but also its reaction time.

Conclusions

By fusing the competing genes with fluorescent tags this mechanism could be tested and employed to create an indicator of the doubling time. By manipulating copy numbers of the competing genes and locus of the gene situated near ter, one can obtain equal average expression of both genes for any doubling time T between 20 and 120 min. Such a toggle would accurately report departures of the doubling time from T.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0483-4) contains supplementary material, which is available to authorized users.

Stability of bacterial toggle switches is enhanced by cell-cycle lengthening by several orders of magnitude

Bistable regulatory elements are important for nongenetic inheritance, increase of cell-to-cell heterogeneity allowing adaptation, and robust responses at the population level. Here, we study computationally the bistable genetic toggle switch-a small regulatory network consisting of a pair of mutual repressors-in growing and dividing bacteria. We show that as cells with an inhibited growth exhibit high stability of toggle states, cell growth and divisions lead to a dramatic increase of toggling rates. The toggling rates were found to increase with rate of cell growth, and can be up to six orders of magnitude larger for fast growing cells than for cells with the inhibited growth. The effect is caused mainly by the increase of protein and mRNA burst sizes associated with the faster growth. The observation that fast growth dramatically destabilizes toggle states implies that rapidly growing cells may vigorously explore the epigenetic landscape enabling nongenetic evolution, while cells with inhibited growth adhere to the local optima. This can be a clever population strategy that allows the slow growing (but stress resistant) cells to survive long periods of unfavorable conditions. Simultaneously, at favorable conditions, this stress resistant (but slowly growing-or not growing) subpopulation may be replenished due to a high switching rate from the fast growing population.

A balanced ratio of proteins from gene G and frameshift-extended gene GT is required for phage lambda tail assembly

In bacteriophage λ, the overlapping open reading frames G and T are expressed by a programmed translational frameshift similar to that of the gag-pol genes of many retroviruses to produce the proteins gpG and gpGT. An analogous frameshift is widely conserved among other dsDNA tailed phages in their corresponding "G" and "GT" tail genes even in the absence of detectable sequence homology. The longer protein gpGT is known to be essential for tail assembly, but the requirement for the shorter gpG remained unclear because mutations in gene G affect both proteins. A plasmid system that can direct the efficient synthesis of tails was created and used to show that gpG and gpGT are both essential for correct tail assembly. Phage complementation assays under conditions where levels of plasmid-expressed gpG or gpGT could be altered independently revealed that the correct molar ratio of these two related proteins, normally determined by the efficiency of the frameshift, is also crucial for efficient assembly of functional tails. Finally, the physical connection between the G and T domains of gpGT, a consequence of the frameshift mechanism of protein expression, appears to be important for efficient tail assembly.

High-resolution view of bacteriophage lambda gene expression by ribosome profiling

Bacteriophage lambda is one of the most extensively studied organisms and has been a primary model for understanding basic modes of genetic regulation. Here, we examine the progress of lambda gene expression during phage development by ribosome profiling and, thereby, provide a very-high-resolution view of lambda gene expression. The known genes are expressed in a predictable fashion, authenticating the analysis. However, many previously unappreciated potential open reading frames become apparent in the expression analysis, revealing an unexpected complexity in the pattern of lambda gene function.

Toggle switch: noise determines the winning gene

Bistable regulatory elements enhance heterogeneity in cell populations and, in multicellular organisms, allow cells to specialize and specify their fate. Our study demonstrates that in a system of bistable genetic switch, the noise characteristics control in which of the two epigenetic attractors the cell population will settle. We focus on two types of noise: the gene switching noise and protein dimerization noise. We found that the change of magnitudes of these noise components for one of the two competing genes introduces a large asymmetry of the protein stationary probability distribution and changes the relative probability of individual gene activation. Interestingly, an increase of noise associated with a given gene can either promote or suppress the activation of the gene, depending on the type of noise. Namely, each gene is repressed by an increase of its gene switching noise and activated by an increase of its protein-product dimerization noise. The observed effect was found robust to the large, up to fivefold deviations of the model parameters. In summary, we demonstrated that noise itself may determine the relative strength of the epigenetic attractors, which may provide a unique mode of control of cell fate decisions.

Type of noise defines global attractors in bistable molecular regulatory systems

The aim of this study is to demonstrate that in molecular dynamical systems with the underlying bi- or multistability, the type of noise determines the most strongly attracting steady state or stochastic attractor. As an example we consider a simple stochastic model of autoregulatory gene with a nonlinear positive feedback, which in the deterministic approximation has two stable steady state solutions. Three types of noise are considered: transcriptional and translational - due to the small number of gene product molecules and the gene switching noise - due to gene activation and inactivation transitions. We demonstrate that the type of noise in addition to the noise magnitude dictates the allocation of probability mass between the two stable steady states. In particular, we found that when the gene switching noise dominates over the transcriptional and translational noise (which is characteristic of eukaryotes), the gene preferentially activates, while in the opposite case, when the transcriptional noise dominates (which is characteristic of prokaryotes) the gene preferentially remains inactive. Moreover, even in the zero-noise limit, when the probability mass generically concentrates in the vicinity of one of two steady states, the choice of the most strongly attracting steady state is noise type-dependent. Although the epigenetic attractors are defined with the aid of the deterministic approximation of the stochastic regulatory process, their relative attractivity is controlled by the type of noise, in addition to noise magnitude. Since noise characteristics vary during the cell cycle and development, such mode of regulation can be potentially employed by cells to switch between alternative epigenetic attractors.

Assembly mechanism is the key determinant of the dosage sensitivity of a phage structural protein

Altering the expression level of proteins that are subunits of complexes has been proposed to be particularly detrimental because the resulting stoichiometric imbalance among components would lead to misassembly of the complex. Here we show that assembly of the phage HK97 connector complex is severely inhibited by the overexpression of one of its component proteins, gp6. However, this effect is a result of the unusual mechanism by which the oligomerization and assembly of gp6 are controlled. Alteration of this mechanism by single amino acid substitutions leads to a reversal of the response to gp6 overexpression. Surprisingly, the binding partner of gp6 within the phage particle is expressed at a 500-fold higher concentration despite their identical stoichiometry within the complex. Our data emphasize that a generalized prediction of the effects of changes in the expression level of protein complex subunits is very difficult because these effects are dependent upon assembly mechanism.

Conserved translational frameshift in dsDNA bacteriophage tail assembly genes

A programmed translational frameshift similar to frameshifts in retroviral gag-pol genes and bacterial insertion elements was found to be strongly conserved in tail assembly genes of dsDNA phages and to be independent of sequence similarities. In bacteriophage lambda, this frameshift controls production of two proteins with overlapping sequences, gpG and gpGT, that are required for tail assembly. We developed bioinformatic approaches to identify analogous -1 frameshifting sites and experimentally confirmed our predictions for five additional phages. Clear evidence was also found for an unusual but analogous -2 frameshift in phage Mu. Frameshifting sites could be identified for most phages with contractile or noncontractile tails whose length is controlled by a tape measure protein. Phages from a broad spectrum of hosts spanning Eubacteria and Archaea appear to conserve this frameshift as a fundamental component of their tail assembly mechanisms, supporting the idea that their tail genes share a common, distant ancestry.

Complete genomic sequence of the virulent Salmonella bacteriophage SP6

We report the complete genome sequence of enterobacteriophage SP6, which infects Salmonella enterica serovar Typhimurium. The genome contains 43,769 bp, including a 174-bp direct terminal repeat. The gene content and organization clearly place SP6 in the coliphage T7 group of phages, but there is approximately 5 kb at the right end of the genome that is not present in other members of the group, and the homologues of T7 genes 1.3 through 3 appear to have undergone an unusual reorganization. Sequence analysis identified 10 putative promoters for the SP6-encoded RNA polymerase and seven putative rho-independent terminators. The terminator following the gene encoding the major capsid subunit has a termination efficiency of about 50% with the SP6-encoded RNA polymerase. Phylogenetic analysis of phages related to SP6 provided clear evidence for horizontal exchange of sequences in the ancestry of these phages and clearly demarcated exchange boundaries; one of the recombination joints lies within the coding region for a phage exonuclease. Bioinformatic analysis of the SP6 sequence strongly suggested that DNA replication occurs in large part through a bidirectional mechanism, possibly with circular intermediates.

Bacteriophage lambda lysis gene product modified and inserted into Escherichia coli outer membrane: Rz1 lipoprotein

Lysis proteins of bacteriophage lambda were localized in different parts of the host envelope: S in the inner membrane,36 Rz in the membrane adhesion sites,14 and Rz1 in the outer membrane. The R gene product, the transglycosylase destroying bacterial murein, is a soluble protein. Computer-assisted analysis of the Rz1 protein amino acids sequence revealed that its N-terminal part contained the site 15VVVG [symbol: see text] C20, which could be recognizable for the SPase II and cleaved leaving lipid modified C20 as the N-terminal amino acid of the mature protein. Microsequencing of the Rz1 protein isolated from the expression products of E. coli [pSB54] carrying the Rz1 gene showed that the N-terminal part of the protein was cleaved as predicted. Lipid labeling with [3H]palmitate confirmed the expectation that Rz1 was a lipoprotein. E. coli [pSB54] treated with globomycin accumulated prolipoprotein, the Rz1 precursor, which was detectable by the anti-Rz1 serum on electropherograms as the 6.5-kDa protein, larger than mature protein. Physiological function of the Rz1 protein remains to be discovered, but as a first hint we noticed that it evokes increase of the fraction of adhesion sites of outer and inner membranes when overproduced from pSB54. The same effect was observed in induced E. coli (lambda) just before the lysis onset, however, one should be cautious in interpreting the results obtained in conditions of the overproduction of the Rz1 lipoprotein.

The RNA chain elongation rate of the lambda late mRNA is unaffected by high levels of ppGpp in the absence of amino acid starvation

In this study, the effects of high levels of guanosine tetraphosphate (ppGpp) on the decay and RNA chain elongation kinetics of the bacteriophage lambda late transcript in Escherichia coli were examined in the absence of amino acid starvation. The accumulation, mRNA decay kinetics, and RNA chain elongation rate of the lambda late mRNA were determined after heat induction of lambdacI857 lysogens in the presence of high levels of ppGpp induced from a RelAalpha fragment-overproducing plasmid. The accumulation kinetics and elongation rate determinations of the late mRNA were made at long times after induction to allow a new steady state of transcriptional activities under conditions of elevated intracellular levels of ppGpp. The results indicate no prolonged or significant effect on either mRNA decay or the RNA chain elongation rate of the late mRNA as a result of elevated ppGpp levels. Surprisingly, the RNA chain elongation rate determinations indicate an RNA polymerase processivity of approximately 90-100 nucleotides/s for the lambda late transcript despite the presence of high levels of ppGpp. The results are discussed in terms of various models for regulation of stable and messenger RNA synthesis in E. coli.

Translational initiation at the coat-protein gene of phage MS2: native upstream RNA relieves inhibition by local secondary structure

Maximal translation of the coat-protein gene from RNA bacteriophage MS2 requires a contiguous stretch of native MS2 RNA that extends hundreds of nucleotides upstream from the translational start site. Deletion of these upstream sequences from MS2 cDNA plasmids results in a 30-fold reduction of translational efficiency. By site-directed mutagenesis, we show that this low level of expression is caused by a hairpin structure centred around the initiation codon. When this hairpin is destabilized by the introduction of mismatches, expression from the truncated messenger increases 20-fold to almost the level of the full-length construct. Thus, the translational effect of hundreds of upstream nucleotides can be mimicked by a single substitution that destabilizes the structure. The same hairpin is also present in full-length MS2 RNA, but there it does not impair ribosome binding. Apparently, the upstream RNA somehow reduces the inhibitory effect of the structure on translational initiation. The upstream MS2 sequence does not stimulate translation when cloned in front of another gene, nor can unrelated RNA segments activate the coat-protein gene. Several possible mechanisms for the activation are discussed and a function in gene regulation of the phage is suggested.

Bacteriophage lysis: mechanism and regulation

Bacteriophage lysis involves at least two fundamentally different strategies. Most phages elaborate at least two proteins, one of which is a murein hydrolase, or lysin, and the other is a membrane protein, which is given the designation holin in this review. The function of the holin is to create a lesion in the cytoplasmic membrane through which the murein hydrolase passes to gain access to the murein layer. This is necessary because phage-encoded lysins never have secretory signal sequences and are thus incapable of unassisted escape from the cytoplasm. The holins, whose prototype is the lambda S protein, share a common organization in terms of the arrangement of charged and hydrophobic residues, and they may all contain at least two transmembrane helical domains. The available evidence suggests that holins oligomerize to form nonspecific holes and that this hole-forming step is the regulated step in phage lysis. The correct scheduling of the lysis event is as much an essential feature of holin function as is the hole formation itself. In the second strategy of lysis, used by the small single-stranded DNA phage phi X174 and the single-stranded RNA phage MS2, no murein hydrolase activity is synthesized. Instead, there is a single species of small membrane protein, unlike the holins in primary structure, which somehow causes disruption of the envelope. These lysis proteins function by activation of cellular autolysins. A host locus is required for the lytic function of the phi X174 lysis gene E.

In the Escherichia coli lacZ gene the spacing between the translating ribosomes is insensitive to the efficiency of translation initiation

We have constructed a series of 44 Escherichia coli strains in which the chromosomal region corresponding to the Ribosome Binding Site (RBS) of the lacZ gene, has been replaced by small DNA fragments harboring either RBSs from other genes, or artificial RBSs. The beta-galactosidase expression from these strains ranges from 1 to 130 per cent of that of the parental strain. Using this collection, we demonstrate here that strain-to-strain variations in expression are paralleled by nearly equivalent variations in lacZ mRNA content. We propose that, in this system, polarity and mRNA stability are tightly coupled to translation initiation, so that changes in RBS efficiency are detected mainly as changes in mRNA concentration rather than in the spacing between translating ribosomes. In addition, we show that the mRNA sequence immediately downstream from the initiator codon influences per se the lifetime of the lacZ mRNA. We discuss the mechanism of the interdependence between translation, transcription and degradation in this gene, and speculate about the general role of this interdependence in determining the expression of bacterial genes.

Differential gene expression from the Escherichia coli atp operon mediated by segmental differences in mRNA stability

The atp operon of Escherichia coli directs synthesis rates of protein subunits that are well matched to the requirements of assembly of the membrane-bound H(+)-ATPase (alpha 3 beta 3 gamma 1 delta 1 epsilon 1a1b2c10-15). Segmental differences in mRNA stability are shown to contribute to the differential control of atp gene expression. The first two genes of the operon, atpl and atpB, are rapidly inactivated at the mRNA level. The remaining seven genes are more stable. It has previously been established that the translational efficiencies of the atp genes vary greatly. Thus differential expression from this operon is achieved via post-transcriptional control exerted at two levels. Neither enhancement of translational efficiency nor insertion of repetitive extragenic palindromic (REP) sequences into the atplB intercistronic region stabilized atpl. We discuss the implications of these results in terms of the pathway of mRNA degradation and of the role of mRNA stability in the control of gene expression.

Interference with phage lambda development by the small subunit of the phage 21 terminase, gp1

Bacteriophage lambda development is blocked in cells carrying a plasmid that expresses the terminase genes of phage 21. The interference is caused by the small subunit of phage 21 terminase, gp1. Mutants of lambda able to form plaques in the presence of gp1 include sti mutants. One such mutation, sti30, is an A. T-to-G.C transition mutation at base pair 184 on the lambda chromosome. The sti30 mutation extends the length of the ribosome-binding sequence of the Nul gene that is complementary to the 3' end of the 16S rRNA from GGA to GGAG. The sti30 mutation causes a approximately 50-fold increase in the level of expression of a Nul-lacZ reporter gene, indicating that the sti30 mutation overcomes the gp1 inhibition by increasing the level of expression of gpNul. Although the Nul and A genes of lambda overlap, the sti30 mutation has little effect on the level of gpA expression, indicating that translational coupling does not occur.

Ribosome-messenger recognition: mRNA target sites for ribosomal protein S1

Ribosomal protein S1 is known to play an important role in translational initiation, being directly involved in recognition and binding of mRNAs by 30S ribosomal particles. Using a specially developed procedure based on efficient crosslinking of S1 to mRNA induced by UV irradiation, we have identified S1 binding sites on several phage RNAs in preinitiation complexes. Targets for S1 on Q beta and fr RNAs are localized upstream from the coat protein gene and contain oligo(U)-sequences. In the case of Q beta RNA, this S1 binding site overlaps the S-site for Q beta replicase and the site for S1 binding within a binary complex. It is reasonable that similar U-rich sequences represent S1 binding sites on bacterial mRNAs. To test this idea we have used E. coli ssb mRNA prepared in vitro with the T7 promoter/RNA polymerase system. By the methods of toeprinting, enzymatic footprinting, and UV crosslinking we have shown that binding of the ssb mRNA to 30S ribosomes is S1-dependent. The oligo(U)-sequence preceding the SD domain was found to be the target for S1. We propose that S1 binding sites, represented by pyrimidine-rich sequences upstream from the SD region, serve as determinants involved in recognition of mRNA by the ribosome.

Ribosome-messenger recognition: mRNA target sites for ribosomal protein S1

Ribosomal protein S1 is known to play an important role in translational initiation, being directly involved in recognition and binding of mRNAs by 30S ribosomal particles. Using a specially developed procedure based on efficient crosslinking of S1 to mRNA induced by UV irradiation, we have identified S1 binding sites on several phage RNAs in preinitiation complexes. Targets for S1 on Q beta and fr RNAs are localized upstream from the coat protein gene and contain oligo(U)-sequences. In the case of Q beta RNA, this S1 binding site overlaps the S-site for Q beta replicase and the site for S1 binding within a binary complex. It is reasonable that similar U-rich sequences represent S1 binding sites on bacterial mRNAs. To test this idea we have used E. coli ssb mRNA prepared in vitro with the T7 promoter/RNA polymerase system. By the methods of toeprinting, enzymatic footprinting, and UV crosslinking we have shown that binding of the ssb mRNA to 30S ribosomes is S1-dependent. The oligo(U)-sequence preceding the SD domain was found to be the target for S1. We propose that S1 binding sites, represented by pyrimidine-rich sequences upstream from the SD region, serve as determinants involved in recognition of mRNA by the ribosome.

Overexpression of N antitermination proteins of bacteriophages lambda, 21, and P22: loss of N protein specificity

The N protein of bacteriophage lambda (N lambda) modifies Escherichia coli RNA polymerase in such a way that it transcribes through termination signals, a process called antitermination. N antitermination normally occurs only if the template contains a specific utilization or nut site upstream of the terminators and only in the presence of host-encoded Nus proteins. The lambda-related phages 21 and P22 produce N analogs, N21 and N22, but these require different nut sites and show a different pattern of functional interaction with one of the Nus factors, NusA, according to whether this protein is of E. coli or Salmonella origin (NusAEc or NusASal). We report the overproduction of N lambda, N21, or N22, each of which was induced by isopropyl-beta-D-thiogalactopyranoside at 37 degrees C from its cloned position downstream from ptac on a high-expression plasmid, each in a host that provided NusAEc or NusASal. Overproduction of each of these N proteins resulted in relaxed specificity for nut, which was shown by the ability to complement N mutants of heterologous phages; NusA specificity was determined by the N type that was present in these complementation tests. We also observed that excess N was able to suppress transcriptional polarity in the particular case of cloned 'trpA, the last gene of the tryptophan operon, although there was no effect on polarity within chromosomal trpE. Such polarity is attributed to the presence of cryptic intragenic terminators that become exposed in the absence of translation. Because there is no known nut site cis to 'trpA, we suggest that the 'trpA segment itself fortuitously contains a nut sequence that is able to function with excess N of any of the types tested and with either NusAEc or NusASal. We also found that excess N of any specificity, or even inactive N with missense mutation, could cause an increase in the level of NusAEc or NusASal, possibly because interaction between N and NusA, but independent of nut, whether functional or not, interferes with the autoregulation of NusA synthesis. These observations highlight the importance of protein concentration for the specificity of interactions both with other proteins and with nucleic acids. They also indicate that the interaction between N and NusA requires nut participation both for specificity and functionality.