Transcription attenuation in Salmonella typhimurium: the significance of rare leucine codons in the leu leader

Mining regulatory 5'UTRs from cDNA deep sequencing datasets

Regulatory 5' untranslated regions (r5'UTRs) of mRNAs such as riboswitches modulate the expression of genes involved in varied biological processes in both bacteria and eukaryotes. New high-throughput sequencing technologies could provide powerful tools for discovery of novel r5'UTRs, but the size and complexity of the datasets generated by these technologies makes it difficult to differentiate r5'UTRs from the multitude of other types of RNAs detected. Here, we developed and implemented a bioinformatic approach to identify putative r5'UTRs from within large datasets of RNAs recently identified by pyrosequencing of the Vibrio cholerae small transcriptome. This screen yielded only approximately 1% of all non-overlapping RNAs along with 75% of previously annotated r5'UTRs and 69 candidate V. cholerae r5'UTRs. These candidates include several putative functional homologues of diverse r5'UTRs characterized in other species as well as numerous candidates upstream of genes involved in pathways not known to be regulated by r5'UTRs, such as fatty acid oxidation and peptidoglycan catabolism. Two of these novel r5'UTRs were experimentally validated using a GFP reporter-based approach. Our findings suggest that the number and diversity of pathways regulated by r5'UTRs has been underestimated and that deep sequencing-based transcriptomics will be extremely valuable in the search for novel r5'UTRs.

Inhibition of translation by consecutive rare leucine codons in E. coli: absence of effect of varying mRNA stability

Consecutive homologous codons that are rarely used in E. coli are known to inhibit translation to varying degrees. As few as two consecutive rare arginine codons exhibit a profound inhibition of translation when they are located in the 5' portion of a gene in E. coli. We have previously shown that nine consecutive rare CUA leucine codons cause almost complete inhibition of translation when they are placed after the 13th codon of a test message (although they do not inhibit translation when they are placed in the middle of the message). In the present work, we report that five consecutive rare CUA leucine codons exhibit approximately a threefold inhibition of translation when they are similarly placed after the 13th codon of a test message, compared to five consecutive common CUG leucine codons, in a T7 RNA polymerase-driven system. Further, by removing RNase III processing sites at the 3' ends of the mRNAs, we have manipulated the stability of the mRNAs encoding the test and control messages to see if decreasing mRNA stability might have an effect on the extent of translation inhibition by the rare leucine codons. However, the inhibition with the less stable mRNAs was similar to that with the stable mRNAs, approximately 3.4-fold, indicating that mRNA stability per se does not have a major influence on the effects of rare codons in this system.

What makes ribosome-mediated transcriptional attenuation sensitive to amino acid limitation?

Ribosome-mediated transcriptional attenuation mechanisms are commonly used to control amino acid biosynthetic operons in bacteria. The mRNA leader of such an operon contains an open reading frame with “regulatory” codons, cognate to the amino acid that is synthesized by the enzymes encoded by the operon. When the amino acid is in short supply, translation of the regulatory codons is slow, which allows transcription to continue into the structural genes of the operon. When amino acid supply is in excess, translation of regulatory codons is rapid, which leads to termination of transcription. We use a discrete master equation approach to formulate a probabilistic model for the positioning of the RNA polymerase and the ribosome in the attenuator leader sequence. The model describes how the current rate of amino acid supply compared to the demand in protein synthesis (signal) determines the expression of the amino acid biosynthetic operon (response). The focus of our analysis is on the sensitivity of operon expression to a change in the amino acid supply. We show that attenuation of transcription can be hyper-sensitive for two main reasons. The first is that its response depends on the outcome of a race between two multi-step mechanisms with synchronized starts: transcription of the leader of the operon, and translation of its regulatory codons. The relative change in the probability that transcription is aborted (attenuated) can therefore be much larger than the relative change in the time it takes for the ribosome to read a regulatory codon. The second is that the general usage frequencies of codons of the type used in attenuation control are small. A small percentage decrease in the rate of supply of the controlled amino acid can therefore lead to a much larger percentage decrease in the rate of reading a regulatory codon. We show that high sensitivity further requires a particular choice of regulatory codon among several synonymous codons for the same amino acid. We demonstrate the importance of a high fraction of regulatory codons in the control region. Finally, our integrated model explains how differences in leader sequence design of the trp and his operons of Escherichia coli and Salmonella typhimurium lead to high basal expression and low sensitivity in the former case, and to large dynamic range and high sensitivity in the latter. The model clarifies how mechanistic and systems biological aspects of the attenuation mechanism contribute to its overall sensitivity. It also explains structural differences between the leader sequences of the trp and his operons in terms of their different functions.

When cells grow and divide, they must continually construct new proteins from the 20 amino acid building blocks according to the instructions of the genetic code. Proteins are made by large macromolecular complexes, ribosomes, where information encoded as base triplets (codons) in messenger RNA sequences, transcribed from the DNA sequences of the genes, is translated into amino acid sequences that determine the functions of all proteins. Rapid growth of cells requires that the supply of each free amino acid is balanced to the demand for it in protein synthesis.

The present work mathematically models a common control mechanism in bacteria, which regulates synthesis of amino acids to eliminate deviations from balanced supply and demand. The mechanism “measures” the speed by which the ribosome translates the codons of a regulated amino acid. When supply is less than demand, translation of these “control” codons is slow, which is sensed by the mechanism and used to increase synthesis of the amino acid.

This paper explains why the mechanism is “hyper-sensitive” to relative changes in supply and demand, and why it is differently designed for control of the enzymes that synthesize the amino acids histidine and tryptophan.

Transcription attenuation

In this review, we describe a variety of mechanisms that bacteria use to regulate transcription elongation in order to control gene expression in response to changes in their environment. Together, these mechanisms are known as attenuation and antitermination, and both involve controlling the formation of a transcription terminator structure in the RNA transcript prior to a structural gene or operon. We examine attenuation and antitermination from the point of view of the different biomolecules that are used to influence the RNA structure. Attenuation of many amino acid biosynthetic operons, particularly in enteric bacteria, is controlled by ribosomes translating leader peptides. RNA-binding proteins regulate attenuation, particularly in gram-positive bacteria such as Bacillus subtilis. Transfer RNA is also used to bind to leader RNAs and influence transcription antitermination in a large number of amino acyl tRNA synthetase genes and several biosynthetic genes in gram-positive bacteria. Finally, antisense RNA is involved in mediating transcription attenuation to control copy number of several plasmids.

Transcriptional and translational regulation of alpha-acetolactate decarboxylase of Lactococcus lactis subsp. lactis

The alpha-acetolactate decarboxylase (ALDC) gene, aldB, is the penultimate gene of the leu-ilv-ald operon, which encodes the three branched-chain amino acid (BCAA) biosynthesis genes in Lactococcus lactis. Its product plays a dual role in the cell: (i) it catalyzes the second step of the acetoin pathway, and (ii) it controls the pool of alpha-acetolactate during leucine and valine synthesis. It can be transcribed from the two promoters present upstream of the leu and ilv genes (P1 and P2) or independently under the control of its own promoter (P3). In this paper we show that the production of ALDC is limited by two mechanisms. First, the strength of P3 decreases greatly during starvation for BCAAs and under other conditions that generally provoke the stringent response. Second, although aldB is actively transcribed from P1 and P2 during BCAA starvation, ALDC is not significantly produced from these transcripts. The aldB ribosome binding site (RBS) appears to be entrapped in a stem-loop, which is itself part of a more complex RNA folding structure. The function of the structure was studied by mutagenesis, using translational fusions with luciferase genes to assess its activity. The presence of the single stem-loop entrapping the aldB RBS was responsible for a 100-fold decrease in the level of aldB translation. The presence of a supplementary secondary structure upstream of the stem-loop led to an additional fivefold decrease of aldB translation. Finally, the translation of the ilvA gene terminating in the latter structure decreased the level of translation of aldB fivefold more, leading to the complete extinction of the reporter gene activity. Since three leucines and one valine are present among the last six amino acids of the ilvA product, we propose that pausing of the ribosomes during translation could modulate the folding of the messenger, as a function of BCAA availability. The purpose of the structure-dependent regulation could be to ensure the minimal production of ALDC required for the control of the acetolactate pool during BCAA synthesis but to avoid its overproduction, which would dissipate acetolactate. Large amounts of ALDC, necessary for operation of the acetoin pathway, could be produced under favorable conditions from the P3 transcripts, which do not contain the secondary structures.

Protein secondary structural types are differentially coded on messenger RNA

Tricodon regions on messenger RNAs corresponding to a set of proteins from Escherichia coli were scrutinized for their translation speed. The fractional frequency values of the individual codons as they occur in mRNAs of highly expressed genes from Escherichia coli were taken as an indicative measure of the translation speed. The tricodons were classified by the sum of the frequency values of the constituent codons. Examination of the conformation of the encoded amino acid residues in the corresponding protein tertiary structures revealed a correlation between codon usage in mRNA and topological features of the encoded proteins. Alpha helices on proteins tend to be preferentially coded by translationally fast mRNA regions while the slow segments often code for beta strands and coil regions. Fast regions correspondingly avoid coding for beta strands and coil regions while the slow regions similarly move away from encoding alpha helices. Structural and mechanistic aspects of the ribosome peptide channel support the relevance of sequence fragment translation and subsequent conformation. A discussion is presented relating the observation to the reported kinetic data on the formation and stabilization of protein secondary structural types during protein folding. The observed absence of such strong positive selection for codons in non-highly expressed genes is compatible with existing theories that mutation pressure may well dominate codon selection in non-highly expressed genes.

Ribosome-mediated translational pause and protein domain organization

Because regions on the messenger ribonucleic acid differ in the rate at which they are translated by the ribosome and because proteins can fold cotranslationally on the ribosome, a question arises as to whether the kinetics of translation influence the folding events in the growing nascent polypeptide chain. Translationally slow regions were identified on mRNAs for a set of 37 multidomain proteins from Escherichia coli with known three-dimensional structures. The frequencies of individual codons in mRNAs of highly expressed genes from E. coli were taken as a measure of codon translation speed. Analysis of codon usage in slow regions showed a consistency with the experimentally determined translation rates of codons; abundant codons that are translated with faster speeds compared with their synonymous codons were found to be avoided; rare codons that are translated at an unexpectedly higher rate were also found to be avoided in slow regions. The statistical significance of the occurrence of such slow regions on mRNA spans corresponding to the oligopeptide domain termini and linking regions on the encoded proteins was assessed. The amino acid type and the solvent accessibility of the residues coded by such slow regions were also examined. The results indicated that protein domain boundaries that mark higher-order structural organization are largely coded by translationally slow regions on the RNA and are composed of such amino acids that are stickier to the ribosome channel through which the synthesized polypeptide chain emerges into the cytoplasm. The translationally slow nucleotide regions on mRNA possess the potential to form hairpin secondary structures and such structures could further slow the movement of ribosome. The results point to an intriguing correlation between protein synthesis machinery and in vivo protein folding. Examination of available mutagenic data indicated that the effects of some of the reported mutations were consistent with our hypothesis.

Decoding with the A:I wobble pair is inefficient

tRNAs with inosine (I) in the first position read three codons ending in U, C and A. However, A-ending codons read with I are rarely used. In Escherichia coli, CGA/U/C are all read solely by tRNAICGArg. CGU and CGC are very common codons, but CGA is very rare. Three independent in vivo assays show that translation of CGA is relatively inefficient. In the first, nine tandem CGA cause a strong rho-mediated polar effect on expression of a lacZ reporter gene. The inhibition is made more extreme by a mutation in ribosomal protein S12 (rpsL), which indicates that ribosomal binding by tRNAICGArg is slow and/or unstable in the CGA cluster. The second assay, in which codons are substituted for the regulatory UGA of the RF2 frameshift, confirms that aa-tRNA selection is slow and/or unstable at CGA. In the third assay, CGA is found to be a poor 5' context for amber suppression, which suggests that an A:I base pair in the P site can interfere with translation of a codon in the A site. Two possible errors, frameshifting and premature termination by RF2, are not significant causes for inefficiency at CGA. It is concluded that the A:I pair destabilizes codon:anticodon complexes during two successive ribosomal cycles, and it is suggested that these properties contribute to the rare usage of codons read with the A:I base pair.

Organization and regulation of genes for amino acid biosynthesis in lactic acid bacteria

The recent description of large clusters of biosynthetic genes in the chromosome of Lactococcus lactis and, to a lesser extent, of Lactobacillus, has brought some information on gene organization and control of gene expression in these organisms. The genes involved in a given amino acid biosynthetic pathway are clustered at a single chromosomal location and form an operon. Additional genes which are not required for the biosynthesis are present within some operons. Genetic signals are, in general, similar to those found in other prokaryotes. Several systems controlling gene expression have been identified and transcription attenuation seems frequent. Among the attenuation mechanisms identified, one resembles that controlling amino acid biosynthesis in many bacteria by ribosome stalling at codons corresponding to limiting amino acid. The others are different and might be related to a new class of attenuation mechanism. Preliminary evidence for a new type of regulatory mechanism, involving a metabolic shunt, is also reviewed.

Unusual codon bias occurring within insertion sequences in Escherichia coli

The large open reading frames of insertion sequences from Escherichia coli were examined for their spatial pattern of codon usage bias and distribution of rarely used codons. There is a bias in codon usage that is generally lower toward the terminal ends of the coding regions, which is reflected in the occurrence of an excess of nonpreferred codons in the 3' portions of the coding regions as compared with the 5' portions. In contrast, typical chromosomal genes have a lower codon usage bias toward the 5' ends of the coding regions. These results imply that the selective forces reflected in codon usage bias may differ according to position within the coding sequence. In addition, these constraints apparently differ in important ways between genes contained in insertion sequences and those in the chromosome.

Transcription attenuation-mediated control of leu operon expression: influence of the number of Leu control codons

Four adjacent Leu codons within the leu leader RNA are critically important in transcription attenuation-mediated control of leu operon expression in Salmonella typhimurium and Escherichia coli (P. W. Carter, D. L. Weiss, H. L. Weith, and J. M. Calvo, J. Bacteriol. 162:943-949, 1985). The leader region from S. typhimurium was altered by site-directed mutagenesis to produce constructs having between one and seven adjacent Leu codons, all CUA. leu operon expression was measured in strains containing six of these constructs, each integrated into the chromosome in a single copy. Operon expression was sufficiently high that all strains grew in minimal medium unsupplemented by leucine. Expression of the operon was measured in strains cultured in such a way that their growth was limited by the intracellular concentration of either leucine or of leucyl-tRNA. In general, the leu operon for each construct responded similarly to the parent construct in terms of the degree of expression as a function of the degree of limitation. However, a strain containing (CUA)1 and, to a certain extent, a strain having (CUA)2 responded somewhat more sluggishly and strains containing (CUA)6 and (CUA)7 responded more sensitively to limitations than did the parent construct. In addition, DNA fragments containing the leu promoter and leader region were used as templates in in vitro transcription reactions employing purified RNA polymerase. With nucleoside triphosphate concentrations of 200 microM, RNA polymerase paused during transcription of the leu leader region at a site about 95 bp downstream from the site of transcription initiation. The halftimes of the pause were 1 min at 37 degrees C and 3 min at 22 degrees C. The pause was lengthened substantially when the GTP concentration was lowered to 20 micromoles. Our results are interpreted most easily in terms of an all-or-none model. Given two Leu control codons, the operon responds with nearly maximum output over a wide range of leucine limitation, and that outcome does not change much with increasing numbers of control codons.

A regulatory element within a gene of a ribosomal protein operon of Escherichia coli negatively controls expression by decreasing the translational efficiency

The trmD operon of Escherichia coli consists of the genes for the ribosomal protein (r-protein) S16, a 21 kDa protein (21K) of unknown function, the tRNA(m1G37)methyltransferase (TrmD), and r-protein L19, in this order. Previously we have shown that the steady-state amount of the two r-proteins exceeds that of the 21K and TrmD proteins 12- and 40-fold, respectively, and that this differential expression is solely explained by translational regulation. Here we have constructed translational gene fusions of the trmD operon and lacZ. The expression of a lacZ fusion containing the first 18 codons of the 21K protein gene is 15-fold higher than the expression of fusions containing 49 or 72 codons of the gene. This suggests that sequences between the 18th and the 49th codon may act as a negative element controlling the expression of the 21K protein gene. Evidence is presented which demonstrates that this regulation is achieved by reducing the efficiency of translation.

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

The late operon of bacteriophage lambda contains the genes encoding the morphogenetic proteins of the phage. These genes are transcribed equally from the single late promoter. Although the functional half-lives of the mRNA for the various genes of this operon vary less than 2-fold, their relative rates of expression have been shown to vary by nearly 1000-fold. This variation could result from differing rates of translation initiation, from overlapping upstream translation, or from differential elongation rates due to the presence of codons for which the corresponding tRNAs are rare. To distinguish between these possibilities, we have cloned sequences surrounding the initiator codons of several of these genes and measured their ability to drive synthesis of hybrid lambda-beta-galactosidase proteins. The rates of expression of the hybrid genes thus produced correlate very well with the natural rates of expression of the corresponding phage genes, suggesting that the rate of initiation of translation controls the relative expression rates of these genes.

Codon choice and gene expression: synonymous codons differ in their ability to direct aminoacylated-transfer RNA binding to ribosomes in vitro

Phe-tRNA (anticodon GAA)--polypeptide-chain elongation factor Tu-GTP ternary complexes react faster with ribosomes programmed with UUC codons than with ribosomes programmed with UUU codons. A similar preference is shown by Leu-tRNA2 (anticodon GAG) complexes, which react faster with ribosomes programmed with CUC than with those programmed with CUU. The difference is seen in the rate of ternary-complex binding to the ribosome; no differences are seen in peptide-bond formation. Highly expressed mRNAs in Escherichia coli favor codons terminating in cytosine rather than uracil when both codons are read by a single tRNA with an anticodon beginning with guanine. The results suggest that intrinsic differences between the efficiencies of synonymous codons play an important role in modulating gene expression in E. coli.

The puf operon region of Rhodobacter sphaeroides

The puf operon of the purple nonsulfur photosynthetic bacterium, Rhodobacter sphaeroides, contains structural gene information for at least two functionally distinct bacteriochlorophyll-protein complexes (light harvesting and reaction center) which are present in a fixed ratio within the photosynthetic intracytoplasmic membrane. Two proximal genes (pufBA) specify subunits of a long wavelength absorbing (i.e., 875 nm) light harvesting complex which are present in the photosynthetic membrane in ≃15 fold excess relative to the reaction center subunits which are encoded by the pufLM genes. This review summarizes recent studies aimed at determining how expression of the R. sphaeroides puf operon region relates to the ratio of individual bacteriochlorophyll-protein complexes found within the photosynthetic membrane. These experiments indicate that puf operon expression may be regulated at the transcriptional, post-transcriptional, translation and post-translational levels. In addition, this review discusses the possible role(s) of newly identified loci upstream of pufB which may be involved in regulating either synthesis or assembly of individual bacteriochrlorophyll-protein complexes as well as the pufX gene, the most distal genetic element within the puf operon whose function is still unknown.

Role of codon choice in the leader region of the ilvGMEDA operon of Serratia marcescens

Leucine participates in multivalent repression of the Serratia marcescens ilvGMEDA operon by attenuation (J.-H. Hsu, E. Harms, and H.E. Umbarger, J. Bacteriol. 164:217-222, 1985), although there is only one single leucine codon that could be involved in this type of control. This leucine codon is the rarely used CUA. The contribution of this leucine codon to the control of transcription by attenuation was examined by replacing it with the commonly used leucine codon CUG and with a nonregulatory proline codon, CCG. These changes left intact the proposed secondary structure of the leader. The effects of the codon changes were assessed by placing the mutant leader regions upstream of the ilvGME structural genes or the cat gene and measuring acetohydroxy acid synthase II, transaminase B, or chloramphenicol acetyltransferase activities in cells grown under limiting and repressing conditions. The presence of the common leucine codon in place of the rare leucine codon reduced derepression by about 70%. Eliminating the leucine codon by converting it to proline abolished leucine control. Furthermore, a possible context effect of the adjacent upstream serine codon on leucine control was examined by changing it into a glycine codon.

Stage-specific selection of alternative transcriptional initiation sites from the 5C actin gene of Drosophila melanogaster

The transcription unit of the 5C actin gene exhibits a complex organization that is unique among the six actin genes of Drosophila melanogaster. Three different mRNA size classes showing distinct patterns of accumulation throughout development are detected on Northern blots. We have determined the structure of the various 5C actin transcripts by exon mapping using strand-specific RNA probes, primer extension analysis, and DNA sequences analysis of both cDNA and genomic clones. All the transcripts share a single protein-coding nucleotide sequence but are heterogeneous in the 5' and 3' untranslated regions. The 5' untranslated region of each transcript consists of either one of two small exons (exon 1 and exon 2) which are alternatively spliced to a single acceptor site 8 bp upstream from the translation initiation codon in exon 3. Results from primer extension analysis suggest that transcription can initiate from either exon 1 or exon 2, and also from a third site within exon 2. We detect an increase in the relative abundance of exon 1-containing transcripts at larval and pupal stages, as well as a change in the proportion of transcripts that initiate at either of the two exon 2 sites. Five polyadenylation sites have been found within three termination/processing regions that define the three size classes of polyadenylated transcripts. The results of our experiments indicate the existence in vivo of all possible combinations of 5' exon with 3' polyadenylation site. However, particular combinations of 5' initiation site and 3' polyadenylation site are preferred at certain developmental stages.

Expression of leucine genes from an extremely thermophilic bacterium in Escherichia coli

The organisation of the leucine genes in Thermus thermophilus HB8 was analysed by examining the ability of recombinant DNAs to complement Escherichia coli mutations. The arrangement of the genes is different from that in the mesophilic bacteria E. coli and Salmonella typhimurium. The promoter responsible for the expression of the leuB, leuC and leuD genes of Thermus HB8 in E. coli was identified. The sequence of Thermus DNA containing this promoter revealed structural similarities to the promoter and attenuator regions of the E. coli leucine operon.