Codon usage determines translation rate in Escherichia coli

Molecular considerations in the evolution of bacterial genes

Synonymous and nonsynonymous substitution rates at the loci encoding glyceraldehyde-3-phosphate dehydrogenase (gap) and outer membrane protein 3A (ompA) were examined in 12 species of enteric bacteria. By examining homologous sequences in species of varying degrees of relatedness and of known phylogenetic relationships, we analyzed the patterns of synonymous and nonsynonymous substitutions within and among these genes. Although both loci accumulate synonymous substitutions at reduced rates due to codon usage bias, portions of the gap and ompA reading frames show significant deviation in synonymous substitution rates not attributable to local codon bias. A paucity of synonymous substitutions in portions of the ompA gene may reflect selection for a novel mRNA secondary structure. In addition, these studies allow comparisons of homologous protein-coding sequences (gap) in plants, animals, and bacteria, revealing differences in evolutionary constraints on this glycolytic enzyme in these lineages.

Low-usage codons in Escherichia coli, yeast, fruit fly and primates

Codon usage is compared between four classes of species, with an emphasis on characterization of low-usage codons. The classes of species analyzed include the bacterium Escherichia coli (ECO), the yeast Saccharomyces cerevisiae (YSC), the fruit fly Drosophila melanogaster (DRO), and several species of primates (PRI) (taken as a group; includes eleven species for which nucleotide sequence data have been reported to GenBank, however, greater than 90% of the sequences were from Homo sapiens). The number of protein-coding sequences analyzed were 968 for ECO, 484 for YSC, 244 for DRO, and 1518 for PRI. Three methods have been used to determine low-usage codons in these species. The first and most common way of assessing codon usage is by summing the number of time codons appear in reading frames of the genome in question. The second way is to examine the distribution of usage in different genes by scoring the number of protein reading frames in which a particular codon does not appear. The third way starts with a similar notion, but instead considers combinations of codons that are missing from the maximum number of genes. These three methods give very similar results. Each species has a unique combination of eight least-used codons, but all species contain the arginine codons, CGA and CGG. The agreement between YSC and PRI is particularly striking as they share six low-usage codons. All six carry the dinucleotide sequence, CG. The eight least-used codons in PRI include all codons that contain the CG dinucleotide sequence. Low-usage codons are clearly avoided in genes encoding abundant proteins for ECO, YSC DRO. In all species, proteins containing a high percentage of low-usage codons could be characterized as cases where an excess of the protein could be detrimental. Low codon usage is relatively insensitive to gross base composition. However, dinucleotide usage can sometimes influence codon usage. This is particularly notable in the case of CG dinucleotides in PRI.

Cysteine, even in low concentrations, induces transient amino acid starvation in Escherichia coli

Cysteine, in concentrations down to 0.04 micrograms/ml, induces transient amino acid starvation in Escherichia coli growing in minimal medium. The duration depends on the concentration and is 5 min at 2 micrograms of cysteine per ml. At low cysteine concentrations, threonine and isoleucine almost completely abolish the starvation.

Codon bias and gene expression

The frequencies with which individual synonymous codons are used to code their cognate amino acids is quite variable from genome to genome and within genomes, from gene to gene. One particularly well documented codon bias is that associated with highly expressed genes in bacteria as well as in yeast; this is the so-called major codon bias. Here, it is suggested that the major codon bias is not an arrangement for regulating individual gene expression. Instead, the data suggest that this codon bias, which is correlated with a corresponding bias of tRNA abundance, is a global arrangement for optimizing the growth efficiency of cells. On the practical side, it is suggested that heterologous gene expression is not as sensitive to codon bias as previously thought, but that it is quite sensitive to other characteristics of the heterologous gene.

Expression of an L-alanine dehydrogenase gene in Zymomonas mobilis and excretion of L-alanine

An approach to broaden the product range of the ethanologenic, gram-negative bacterium Zymomonas mobilis by means of genetic engineering is presented. Gene alaD for L-alanine dehydrogenase (EC 1.4.1.1.) from Bacillus sphaericus was cloned and introduced into Z. mobilis. Under the control of the strong promoter of the pyruvate decarboxylase (pdc) gene, the enzyme was expressed up to a specific activity of nearly 1 mu mol . min -1 . mg of protein -1 in recombinant cells. As a results of this high L-alanine dehydrogenase activity, growing cells excreted up to 10 mmol of alanine per 280 mmol of glucose utilized into a mineral salts medium. By the addition of 85 mM NH4+ to the medium, growth of the recombinant cells stopped, and up to 41 mmol alanine was secreted. As alanine dehydrogenase competed with pyruvate decarboxylase (PDC) (EC 4.1.1.1.) for the same substrate (pyruvate), PDC activity was reduced by starvation for the essential PDC cofactor thiamine PPi. A thiamine auxotrophy mutant of Z. mobilis which carried the alaD gene was starved for 40 h in glucose-supplemented mineral salts medium and then shifted to mineral salts medium with 85 mM NH4+ and 280 mmol of glucose. The recombinants excreted up to 84 mmol of alanine (7.5 g/liter) over 25 h. Alanine excretion proceeded at an initial velocity of 238 nmol . min-1 . mg [dry weight]-1. Despite this high activity, the excretion rate seemed to be a limiting factor, as the intracellular concentration of alanine was as high as 260 mM at the beginning of the excretion phase and decreased to 80 to 90 mM over 24 h.

Role of transcription pausing in the control of the pyrE attenuator in Escherichia coli

Expression of the Escherichia coli pyrE gene is regulated by transcription attenuation in the intercistronic orfE-pyrE region and modulated by the distance between the transcribing RNA polymerase and the leading ribosome as a function of the supply of UTP and GTP. In this communication we show that pyrE expression is hyper-repressed in vivo following addition of uracil in strains carrying the nusAcs10 mutation. This phenotype, previously seen in rpsL1204 strains whose ribosomes are pseudodependent on streptomycin and work at suboptimal elongation rate, indicates that RNA polymerase escapes from the ribosomes in the pyrE attenuator region in the nusA mutant. In vitro transcription studies revealed that the build-up of the full-length attenuated orfE transcript occurred more slowly in the presence of the NusA protein than in its absence. Moreover, the NusA protein enhanced several transcription pauses through the orfE gene. These effects were more pronounced when low concentrations of either UTP or GTP were used than at low concentrations of either CTP or ATP. The results indicate that the NusA protein is required for proper regulation of pyrE gene expression and is involved, together with the NTP pools, in maintaining the coupling between transcription and translation in the pyrE attenuator region by inhibiting RNA chain elongation.

Codon context

The analysis of coding sequences reveals nonrandomness in the context of both sense and stop codons. Part of this is related to nucleotide doublet preference, seen also in non-coding sequences and thought to arise from the dependence of mutational events on surrounding sequence. Another nonrandom context element, relating the wobble nucleotides of successive codons, is observed even when doublet preference, codon usage and bias in amino acid doublets are all allowed for. Several phenomena related to protein synthesis have been shown in vivo to be affected by the nucleotide sequence around codons. Thus, nonsense and missense suppression, elongation rate, precision of tRNA selection and polypeptide chain termination are all affected by codon context. At present, it remains unclear how these phenomena may influence the evolution of nonrandomness in the context of codons in natural sequences.

Processivity errors of gene expression in Escherichia coli

Not all ribosomes that initiate translation of an mRNA sequence will successfully complete it and produce a full-length protein product. By comparing the amounts of lacZ monomer and lacZ dimer protein expressed from a plasmid in a strictly controlled assay, we calculate a dimer to monomer ratio of 0.76. We interpret this to mean that ribosomes have a 76% chance of completing the synthesis of a beta-galactosidase polypeptide. The remaining 24% of the initiated chains end in processivity accidents. For the wild-type, premature RNA polymerase termination is found to account for roughly one-third of the processivity accidents. For the hyperaccurate SmP mutant, we observe a processivity of 0.28, but the presence of streptomycin improves this to 0.50. Thus, the hyperaccuracy with respect to missense substitutions for this mutant is accompanied by a reduced processivity. Addition of streptomycin increase the first error class and reduces the second one. This finding is relevant to the optimization of ribosome function and the growth performance of ribosome mutants.

Growth rate dependence of global amino acid composition

The global amino acid composition of bacteria growing in different media has been studied. The data reveal significant changes in the amino acid composition in the growth rate range between 0.5 and 2.1 doublings per hour at 37 degrees C. The changes are consistent with a progressive simplification of the protein population and mRNA pools as the growth rates increase.

Nucleotide sequence, organisation and structural analysis of the products of genes in the nirB-cysG region of the Escherichia coli K-12 chromosome

The DNA sequence and derived amino-acid sequence of a 5618-base region in the 74-min area of the Escherichia coli chromosome has been determined in order to locate the structural gene, nirB, for the NADH-dependent nitrite reductase and a gene, cysG, required for the synthesis of the sirohaem prosthetic group. Three additional open reading frames, nirD, nirE and nirC, were found between nirB and cysG. Potential binding sites on the NirB protein for NADH and FAD, as well as conserved central core and interface domains, were deduced by comparing the derived amino-acid sequence with those of database proteins. A directly repeated sequence, which includes the motif -Cys-Xaa-Xaa-Cys-, is suggested as the binding site for either one [4Fe-4S] or two [2Fe-2S] clusters. The nirD gene potentially encodes a soluble, cytoplasmic protein of unknown function. No significant similarities were found between the derived amino-acid sequence of NirD and either NirB or any other protein in the database. If the nirE open reading frame is translated, it would encode a 33-amino-acid peptide of unknown function which includes 8 phenylalanyl residues. The product of the nirC gene is a highly hydrophobic protein with regions of amino-acid sequence similar to cytochrome oxidase polypeptide 1.

Translation initiation in Escherichia coli: old and new questions

We discuss the features of Escherichia coli mRNAs which determine where and how efficiently translation is initiated. We have shown that DNA fragments comprising 60-80 nucleotides that bracket the initiation codon of real genes generally promote translation when inserted within a foreign mRNA, while those not corresponding to an authentic gene start do not do so even if they include a Shine-Dalgarno-like element followed by AUG or GUG. Therefore, the information that pinpoints the correct start sites, while extending beyond the mere presence of these elements, remains essentially local. The possible nature of this information is discussed. Next, we point out that, in order to remain accessible, translational starts must escape long-range base-pairing within large mRNAs, and we argue that the tight coupling between translation and transcription plays an important role in achieving this. Finally, we discuss two intriguing situations in which the initiation frequency should be dependent upon the rate of translation elongation.

Metabolic growth rate control in Escherichia coli may be a consequence of subsaturation of the macromolecular biosynthetic apparatus with substrates and catalytic components

In this paper, the Escherichia coli cell is considered as a system designed for rapid growth, but limited by the medium. We propose that this very design causes the cell to become subsaturated with precursors and catalytic components at all levels of macromolecular biosynthesis and leads to a molecular sharing economy at a high level of competition inside the cell. Thus, the promoters compete with each other in the binding of a limited amount of free RNA polymerase and the ribosome binding sites on the mRNA chains compete with each other for the free ribosomes. The macromolecular chain elongation reactions sequester a considerable proportion of the total amount of RNA polymerase and ribosomes in the cells. We propose that the degree of subsaturation of the macromolecular biosynthetic apparatus renders a variable fraction of RNA polymerase and ribosomes unavailable for the initiation of new chain synthesis and that this, at least in part, determines the composition of the cell as a function of the growth rate. Thus, at rapid growth, the high speed of the elongation reactions enables the cell to increase the concentrations of free RNA polymerase and ribosomes for initiation purposes. Furthermore, it is proposed that the speed of RNA polymerase movement is adjusted to the performance speed of the ribosomes. Mechanistically, this adjustment of the coupling between transcription and translation involves transcriptional pause sites along the RNA chains, the adjustment of the saturation level of RNA polymerase with the nucleoside triphosphate substrates, and the concentration of ppGpp, which is known to inhibit RNA chain elongation. This model is able to explain the stringent response and the control of stable RNA and of ribosome synthesis in steady states and in shifts, as well as the rate of overall protein synthesis as a function of the growth rate.

Codon preferences in free-living microorganisms

A popular interpretation of the major codon preference is that it reflects the operation of a regulatory device that controls the expression of individual proteins. In this popular model, rapidly translated codons are thought to promote the accumulation of the highly expressed proteins and slowly translated codons are thought to retard the expression of poorly expressed proteins. However, this widely accepted model is not supported by kinetic theory or by experimental results. A less fashionable model in which the major codon preference has nothing to do with the expression level of the individual proteins is forwarded. In this model, the major codon preference is viewed as a global strategy to support the efficient function of the translation system and thereby to maximize the growth rates of cells under favorable conditions.

Identification and comparison of stable and unstable mRNAs in Saccharomyces cerevisiae

We developed a procedure to measure mRNA decay rates in the yeast Saccharomyces cerevisiae and applied it to the determination of half-lives for 20 mRNAs encoded by well-characterized genes. The procedure utilizes Northern (RNA) or dot blotting to quantitate the levels of individual mRNAs after thermal inactivation of RNA polymerase II in an rpb1-1 temperature-sensitive mutant. We compared the results of this procedure with results obtained by two other procedures (approach to steady-state labeling and inhibition of transcription with Thiolutin) and also evaluated whether heat shock alter mRNA decay rates. We found that there are no significant differences in the mRNA decay rates measured in heat-shocked and non-heat-shocked cells and that, for most mRNAs, different procedures yield comparable relative decay rates. Of the 20 mRNAs studied, 11, including those encoded by HIS3, STE2, STE3, and MAT alpha 1, were unstable (t1/2 less than 7 min) and 4, including those encoded by ACT1 and PGK1, were stable (t1/2 greater than 25 min). We have begun to assess the basis and significance of such differences in the decay rates of these two classes of mRNA. Our results indicate that (i) stable and unstable mRNAs do not differ significantly in their poly(A) metabolism; (ii) deadenylation does not destabilize stable mRNAs; (iii) there is no correlation between mRNA decay rate and mRNA size; (iv) the degradation of both stable and unstable mRNAs depends on concomitant translational elongation; and (v) the percentage of rare codons present in most unstable mRNAs is significantly higher than in stable mRNAs.

Identification and comparison of stable and unstable mRNAs in Saccharomyces cerevisiae

We developed a procedure to measure mRNA decay rates in the yeast Saccharomyces cerevisiae and applied it to the determination of half-lives for 20 mRNAs encoded by well-characterized genes. The procedure utilizes Northern (RNA) or dot blotting to quantitate the levels of individual mRNAs after thermal inactivation of RNA polymerase II in an rpb1-1 temperature-sensitive mutant. We compared the results of this procedure with results obtained by two other procedures (approach to steady-state labeling and inhibition of transcription with Thiolutin) and also evaluated whether heat shock alter mRNA decay rates. We found that there are no significant differences in the mRNA decay rates measured in heat-shocked and non-heat-shocked cells and that, for most mRNAs, different procedures yield comparable relative decay rates. Of the 20 mRNAs studied, 11, including those encoded by HIS3, STE2, STE3, and MAT alpha 1, were unstable (t1/2 less than 7 min) and 4, including those encoded by ACT1 and PGK1, were stable (t1/2 greater than 25 min). We have begun to assess the basis and significance of such differences in the decay rates of these two classes of mRNA. Our results indicate that (i) stable and unstable mRNAs do not differ significantly in their poly(A) metabolism; (ii) deadenylation does not destabilize stable mRNAs; (iii) there is no correlation between mRNA decay rate and mRNA size; (iv) the degradation of both stable and unstable mRNAs depends on concomitant translational elongation; and (v) the percentage of rare codons present in most unstable mRNAs is significantly higher than in stable mRNAs.

Analysis of regulation of the ilvGMEDA operon by using leader-attenuator-galK gene fusions

Five of the genes for the biosynthesis of isoleucine and valine form the ilvGMEDA operon of Escherichia coli K-12. Expression of the operon responds to changes in the availability of isoleucine, leucine, and valine (ILV). Addition of an excess of all three amino acids results in reduced expression of the operon, whereas limitation for one of the three amino acids causes an increase in expression. The operon is preceded by a leader-attenuator which clearly regulates the increased expression that occurs due to reduced aminoacylation of tRNA. To assess the factors that result in the reduced expression of this operon upon the addition of ILV, a series of plasmids were constructed in which the ilv regulatory region was fused to galK. In response to addition of the amino acids, expression of the galK gene fused to the leader-attenuator decreased five- to sevenfold, instead of the twofold observed for the chromosomal operon. A deletion analysis with these plasmids indicated that the ILV-specific decrease in expression required an intact leader-attenuator but not ilvGp2 or the DNA that precedes this promoter. This conclusion was supported by both S1 nuclease analysis of transcription initiation and determination of galK mRNA levels by RNA-RNA hybridization.

Suppression of the negative effect of minor arginine codons on gene expression; preferential usage of minor codons within the first 25 codons of the Escherichia coli genes

AGA and AGG codons for arginine are the least used codons in Escherichia coli, which are encoded by a rare tRNA, the product of the dnaY gene. We examined the positions of arginine residues encoded by AGA/AGG codons in 678 E. coli proteins. It was found that AGA/AGG codons appear much more frequently within the first 25 codons. This tendency becomes more significant in those proteins containing only one AGA or AGG codon. Other minor codons such as CUA, UCA, AGU, ACA, GGA, CCC and AUA are also found to be preferentially used within the first 25 codons. The effects of the AGG codon on gene expression were examined by inserting one to five AGG codons after the 10th codon from the initiation codon of the lacZ gene. The production of beta-galactosidase decreased as more AGG codons were inserted. With five AGG codons, the production of beta-galactosidase (Gal-AGG5) completely ceased after a mid-log phase of cell growth. After 22 hr induction of the lacZ gene, the overall production of Gal-AGG5 was 11% of the control production (no insertion of arginine codons). When five CGU codons, the major arginine codon were inserted instead of AGG, the production of beta-galactosidase (Gal-CGU5) continued even after stationary phase and the overall production was 66% of the control. The negative effect of the AGG codons on the Gal-AGG5 production was found to be dependent upon the distance between the site of the AGG codons and the initiation codon. As the distance was increased by inserting extra sequences between the two codons, the production of Gal-AGG5 increased almost linearly up to 8 fold. From these results, we propose that the position of the minor codons in an mRNA plays an important role in the regulation of gene expression possibly by modulating the stability of the initiation complex for protein synthesis.

Translational control of prokaryotic gene expression

Awareness of the importance of post-transcriptional control of gene expression in prokaryotes has grown enormously over the past ten years. In particular, translation features as a step where both control over constitutive rates of gene expression, as well as cis and trans regulation are exercised. Recent research has provided us with new insights into the molecular basis of these phenomena.

Relationship between protein synthesis and concentrations of charged and uncharged tRNATrp in Escherichia coli

We have continuously monitored Trp-tRNA(Trp) concentrations in vivo and, in the same cultures, measured rates of protein synthesis in isogenic stringent and relaxed strains. We have also manipulated cellular charged and uncharged [tRNA(Trp)] by two means: (i) the strain used contains a Trp-tRNA synthetase mutation that increases the Km for Trp; thus, varying exogenous Trp varies cellular Trp-tRNA(Trp); and (ii) we have introduced into the mutant strain a plasmid containing the tRNA(Trp) gene behind an inducible promoter; thus, total [tRNA(Trp)] also can be varied depending on length of induction. The use of these conditions, combined with a previously characterized assay system, has allowed us to demonstrate that (i) the rate of incorporation of Trp into protein is proportional to the fraction of tRNA(Trp) that is charged; for any given total [tRNA(Trp)], this rate is also proportional to the [Trp-tRNA(Trp)]; (ii) uncharged tRNA(Trp) inhibits incorporation of Trp into protein; and (iii) rates of incorporation into protein of at least two other amino acids, Lys and Cys, are also sensitive to [Trp-tRNA(Trp)] and are inhibited by uncharged tRNA(Trp). Our results are consistent with models of translational control that postulate modulating polypeptide chain elongation efficiency by varying concentrations of specific tRNAs.

Expression of tetanus toxin fragment C in E. coli: high level expression by removing rare codons

Tetanus toxin fragment C had been previously expressed in Escherichia coli at 3-4% cell protein. The codon bias for tetanus toxin in Clostridium tetani is very different from that of highly expressed homologous genes in E. coli, resulting in the presence of many rare E. coli codons in the sequence encoding fragment C. We have replaced the coding sequence by sequence optimized for codon usage in E. coli, and show that the expression of fragment C is increased. Although the level of mRNA also increased this appeared to be a secondary consequence of more efficient translation. Complete sequence replacement increased expression to approximately 11-14% cell protein but only after the promoter strength had been improved.