An effect of codon context on the mistranslation of UGU codons in vitro

Codon pair optimization (CPO): a software tool for synthetic gene design based on codon pair bias to improve the expression of recombinant proteins in Pichia pastoris

Background

Codon optimization is a common method to improve protein expression levels in Pichia pastoris and the current strategy is to replace rare codons with preferred codons to match the codon usage bias. However, codon-pair contexts have a profound effect on translation efficiency by influencing both translational elongation rates and accuracy. Until now, it remains untested whether optimized genes based on codon pair bias results in higher protein expression levels compared to codon usage bias.

Results

In this study, an algorithm based on dynamic programming was introduced to develop codon pair optimization (CPO) which is a software tool to provide simple and efficient codon pair optimization for synthetic gene design in Pichia pastoris. Two reporters (MT1-MMP E2C6 and ADAM17 A9B8 scFvs) were employed to test the effects of codon pair bias and CPO optimization on their protein expression levels. Four variants of MT1-MMP E2C6 and ADAM17 A9B8 for each were generated, one variant with the best codon-pair context, one with the worst codon-pair context, one with unbiased codon-pair context, and another optimized based on codon usage. The expression levels of variants with the worst codon-pair context were almost undetectable by Western blot and the variants with the best codon-pair context were expressed well. The expression levels on MT1-MMP E2C6 and ADAM17 A9B8 were more than five times and seven times higher in the optimized sequences based on codon-pair context compared to that based on codon usage, respectively. The results indicated that the codon-pair context-based codon optimization is more effective in enhancing expression of protein in Pichia pastoris.

Conclusions

Codon-pair context plays an important role on the protein expression in Pichia pastoris. The codon pair optimization (CPO) software developed in this study efficiently improved the protein expression levels of exogenous genes in Pichia pastoris, suggesting gene design based on codon pair bias is an alternative strategy for high expression of recombinant proteins in Pichia pastoris.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01696-y.

"Lost in translation: Seeing the forest by focusing on the trees"

A complex process translates messenger RNA (mRNA) base sequence into protein amino acid sequence. Transfer RNAs must recognize 3-base codons in the mRNA to insert the correct amino acids into the growing protein. Codon degeneracy makes decoding complicated in that multiple (synonymous) triplets can encode a single amino acid and multiple tRNAs can have the same anticodon. Over the last twenty years, new developments in structural biology, genome sequencing and bioinformatics has elucidated the intricacies of the ribosome structure and the details of the translation process. High throughput analyses of sequence information support the idea that mRNA folding has a major effect on expression for codons at the 5'-end of mRNA (N-terminal region of a polypeptide). Despite a forest of sequence data, significant details of the complex translation process can escape detection. However, a sensitive translation assay has allowed a single tree in this forest to be revealing.

Gene Designer: a synthetic biology tool for constructing artificial DNA segments

Background

Direct synthesis of genes is rapidly becoming the most efficient way to make functional genetic constructs and enables applications such as codon optimization, RNAi resistant genes and protein engineering. Here we introduce a software tool that drastically facilitates the design of synthetic genes.

Results

Gene Designer is a stand-alone software for fast and easy design of synthetic DNA segments. Users can easily add, edit and combine genetic elements such as promoters, open reading frames and tags through an intuitive drag-and-drop graphic interface and a hierarchical DNA/Protein object map. Using advanced optimization algorithms, open reading frames within the DNA construct can readily be codon optimized for protein expression in any host organism. Gene Designer also includes features such as a real-time sliding calculator of oligonucleotide annealing temperatures, sequencing primer generator, tools for avoidance or inclusion of restriction sites, and options to maximize or minimize sequence identity to a reference.

Conclusion

Gene Designer is an expandable Synthetic Biology workbench suitable for molecular biologists interested in the de novo creation of genetic constructs.

Comparative context analysis of codon pairs on an ORFeome scale

We have developed a system for comparative codon context analysis of open reading frames in whole genomes, providing insights into the rules that govern the evolution of codon-pair context.

Codon context is an important feature of gene primary structure that modulates mRNA decoding accuracy. We have developed an analytical software package and a graphical interface for comparative codon context analysis of all the open reading frames in a genome (the ORFeome). Using the complete ORFeome sequences of Saccharomyces cerevisiae, Schizosaccharomyces pombe, Candida albicans and Escherichia coli, we show that this methodology permits large-scale codon context comparisons and provides new insight on the rules that govern the evolution of codon-pair context.

Gradients in nucleotide and codon usage along Escherichia coli genes

The usage of codons and nucleotide combinations varies along genes and systematic variation causes gradients in usage. We have studied such gradients of nucleotides and nucleotide combinations and their immediate context in Escherichia coli. To distinguish mutational and selectional effects, the genes were subdivided into three groups with different codon usage bias and the gradients of nucleotide usage were studied in each group. Some combinations that can be associated with a propensity for processivity errors show strong negative gradients that become weaker in genes with low codon bias, consistent with a selection on translational efficiency. One of the strongest gradients is for third position G, which shows a pervasive positive gradient in usage in most contexts of surrounding bases.

Transient conformational states of aminoacyl-tRNA during ribosome binding catalyzed by elongation factor Tu

Conformational transitions of Phe-tRNA(Phe) that take place during elongation factor Tu (EF-Tu)-dependent binding to the A site of Escherichia coli ribosomes were followed by transient fluorescence measurements. The fluorescence signal of proflavin replacing dihydrouracil at position 16 or 17 in yeast tRNA(Phe) was utilized to monitor changes of the conformation of the D loop. The ternary complex EF-Tu.GTP.Phe-TRNA(Phe)(Pf16/17) was purified by gel filtration. Upon binding of the complex to the A site of poly(U)-programmed, P-site-blocked ribosomes, the fluorescence changes in several steps. First, the rapid formation of an initial complex gives rise to a small fluorescence increase. Subsequent codon-anticodon recognition leads to a conformational rearrangement of the D loop of the tRNA that is reflected in a major fluorescence increase. Fluorescence-quenching data indicate an unfolding of the D loop in this state. The latter conformational state is short-lived, and the aminoacyl-tRNA refolds during the following rearrangement that occurs after GTP hydrolysis and accompanies the release of the aminoacyl-tRNA from EF-Tu.GDP and/or its accommodation in the A site. Further experiments show that the status of the P site influences the binding to the A site in that the two rearrangement steps are slowed down when the P site is unoccupied and even more so when it is occupied with the near-cognate tRNA(Leu2). In contrast, the occupancy of the E site has no influence on A-site binding, and vice versa, thus excluding any coupling between the two sites.

Absolute in vivo translation rates of individual codons in Escherichia coli. The two glutamic acid codons GAA and GAG are translated with a threefold difference in rate

We have determined the absolute translation rates for four individual codons in Escherichia coli. We used our previously described system for direct measurements of in vivo translation rates using small, in-frame inserts in the lacZ gene. The inserts consisted of multiple synthetic 30 base-pair DNA oligomers with high densities of the four individual codons, GAA (Glu), GAG (Glu), CCG (Pro) and CGA (Arg). Our method is independent of expression level, of mRNA half-life and of transcription rate. Codon GAA was found to be translated with a rate of 21.6 codons/second whereas codon GAG was translated 3.4-fold slower (6.4 codons/s). These two codons are read by the same tRNA species. Codon CCG and CGA are both read by abundant tRNA species but nevertheless we found them to be translated slowly with rates of 5.8 and 4.2 codons/second, respectively. The context of these codons were varied, but we found no significant influence of context on their translation rates and we suggest a mechanism for why context may not affect translation rates. One insert with a low translation rate gave results that most readily can be explained by assuming queue formation of ribosomes on the insert. Such a queue was found to reduce the expression level by approximately 35%. Our experiments allowed us to estimate the average distance between ribosomes and thereby the translation initiation frequency on the wild-type lacZ mRNA. This was found to be one per three seconds.

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.

Ribosomal frameshifting in the yeast retrotransposon Ty: tRNAs induce slippage on a 7 nucleotide minimal site

Ribosomal frameshifting regulates expression of the TYB gene of yeast Ty retrotransposons. We previously demonstrated that a 14 nucleotide sequence conserved between two families of Ty elements was necessary and sufficient to support ribosomal frameshifting. This work demonstrates that only 7 of these 14 nucleotides are needed for normal levels of frameshifting. Any change to the sequence CUU-AGG-C drastically reduces frameshifting; this suggests that two specific tRNAs, tRNA(UAGLeu) and tRNA(CCUArg), are involved in the event. Our tRNA overproduction data suggest that a leucyl-tRNA, probably tRNA(UAGLeu), an unusual leucine isoacceptor that recognizes all six leucine codons, slips from CUU-Leu onto UUA-Leu (in the +1 reading frame) during a translational pause at the AGG-Arg codon induced by the low availability of tRNA(CCUArg), encoded by a single-copy essential gene. Frameshifting is also directional and reading frame specific. Interestingly, frameshifting is inhibited when the "slip" CUU codon is located three codons downstream, but not four or more codons downstream, of the translational initiation codon.

Nonrandom utilization of codon pairs in Escherichia coli

We have analyzed protein-coding sequences of Escherichia coli and find that codon-pair utilization is highly biased, reflecting overrepresentation or underrepresentation of many pairs compared with their random expectations. This effect is over and above that contributed by nonrandomness in the use of amino acid pairs, which itself is highly evident; it is much weaker when nonadjacent codon pairs are examined and virtually disappears when pairs separated by two or three intervening codons are evaluated. There appears to be a high degree of directionality in this bias: any codon that participates in many nonrandom pairs tends to make both over- and underrepresented pairs, but preferentially as a left- or right-hand member. We show a relationship between codon-pair utilization patterns and levels of gene expression: genes encoding proteins expressed at high levels tend to contain more abundant, but more highly underrepresented, codon pairs, relative to genes expressed at low levels. The nonrandom utilization of codon pairs may be a consequence of their effects on translational efficiency, which in turn may be related to the compatibility of adjacent aminoacyl-tRNA isoacceptors at the A and P sites of a translating ribosome.

Expression of peptide chain release factor 2 requires high-efficiency frameshift

Peptide chain release factors are soluble proteins that participate in the stop codon-dependent termination of polypeptide biosynthesis. In Escherichia coli, two release factors are necessary for peptide chain termination: release factor 1 (RF1) specifies UAG- and UAA-dependent termination whereas release factor 2 (RF2) specifies UGA- and UAA-dependent termination. Release factors are found in low concentrations relative to other translation factors, suggesting that their expression is tightly regulated and, accordingly, making the study of their structure-function relationship difficult. RF1 and RF2 exhibit significant sequence homology, probably reflecting their similar functions and perhaps a common evolutionary origin. DNA and peptide sequencing have suggested the existence of a unique mechanism for the autogenous regulation of RF2 in which an in-frame UGA stop codon requires an obligatory +1 frameshift within the coding region of the RF2 gene. In this report we present in vitro experimental results consistent with the autogenous regulation of RF2. Additionally, we used RF2-lacZ gene fusions to demonstrate that autogenous regulation occurs, at least in part, by premature termination at the in-frame stop codon, since deletion of this stop codon leads to overproduction of the RF2-LacZ fusion protein. Frameshifting at this premature termination codon occurs at the remarkably high rate of 50%.

Constraints on codon context in Escherichia coli genes. Their possible role in modulating the efficiency of translation

The constraints on nucleotide sequences of highly and weakly expressed genes from Escherichia coli have been analysed and compared. Differences in synonymous codon spectra in highly and weakly expressed genes lead to different frequencies of nucleotides (in the first and third codon positions) and dinucleotides in the two groups of genes. It has been found that the choice of synonymous codons in highly expressed genes depends on the nucleotides adjacent to the codon. For example, lysine is preferably encoded by the AAA codon if guanosine is 3' to the lysine codon (AAA-G, P less than 10(-9)). And, on the contrary, AAG is used more often than AAA (P less than 0.001) if cytidine is 3' adjacent to lysine. Guanosine occurs more frequently than adenosine 5' to all the lysine codons (AAR, P less than 10(-5), i.e. NNG codons are preferred over the synonymous NNA codons 5' to the positions of lysine in the genes. The context effect was observed in nonsense and missense suppression experiments. Therefore, a hypothesis has been suggested that the efficiency of translation of some codons (for which the constraints on the adjacent nucleotides were found) can be modulated by the codon context. The rules for preferable synonymous codon choice in highly expressed genes depending on the nucleotides surrounding the codon are presented. These rules can be used in the chemical synthesis of genes designed for expression in E. coli.

Fine tuning of ribosomal accuracy

If the rate constant for peptide bond formation were high just after an amino acid incorporation and occasionally switched to a lower value afterwards, then the ribosome could compensate for tRNA imbalance specifically at hungry codons. A rigorous analysis of the scheme proves its effectiveness. For instance, a 10-fold reduction in cognate tRNA concentration may increase the error rate by only a factor of two.

Sense codons are found in specific contexts

The sequence environment of codons in structural genes has been investigated statistically, using computer methods. A set of Escherichia coli genes with abundant products was compared with a set having low gene product levels, in order to detect potential differences associated with expression. The results show striking non-randomness in the nucleotides occurring near codons. These effects are, unexpectedly, very much larger and more homogeneous among the genes with rare products. The intensity of effects in weakly expressed genes suggests that such non-random sequence environments decrease expression. In the weakly expressed set of genes, the 5' neighbor of a codon, and all positions of the 3' neighbor codon are biased. In the highly expressed genes, the first nucleotide of the next codon is a uniquely affected site. The distribution of non-randomness in weakly expressed genes suggests that sequence bias is primarily due to a constraint acting directly on the secondary or tertiary structure of the codon/anticodon. In highly expressed genes, the observed bias suggests an interaction between the codon/anticodon and a site outside the codon/anticodon. Much of the tendency to non-random near-neighbor sequences in weakly expressed genes can be ascribed to a correlation between nearby nucleotides and the wobble nucleotide of the codon, despite the fact that selection of such correlations will alter the amino acid sequence. The favored pattern, in genes expressed at low level, is R YYR or Y RRY. R indicates purine, Y indicates pyrimidine; the space is the boundary between codons. It seems likely that this preference for nearby sequences is the physical basis of the genetic context effect. Under this assumption such sequence biases will affect expression. On this basis, we predict new sites for contextual mutations which decrease expression, and suggest strategy for the design of messages having optimal translational activity.

Does UGA suppressor tRNATrp from Escherichia coli have a unique CCA anticodon sequence?

The properties of the UGA-suppressor tRNATrp (anticodon CCA) from Escherichia coli has greatly influenced ideas about the specificity of codon-anticodon interactions showing that the anticodon sequence is not the sole determinant. However, a recent hypothesis for the mechanism of suppression by this tRNA proposes that the base change in position 24, in the dihydrouridine stem, leads to a change in translational specificity of the tRNA by increasing post-transcriptional modification of cytidine 34, in the anticodon wobble position [M. Yarus (1982) Science (Wash. DC) 218, 646-652]. The enzyme postulated to do this normally modifies C34 of a minor isoleucine isoacceptor specific for AUA codons. This modification should reduce reading of G in the third codon position: affinity chromatography on columns containing immobilised tRNAPro (anticodon VGG) has therefore been employed to isolate an enriched population of the putative suppressor species, if such a sub-population exists. The results obtained are difficult to reconcile with the presence of a subfraction, modified post-transcriptionally in C34, responsible for the suppression of UGA. This argues in favour of the previously advanced hypothesis for the mechanism of suppression, which depends on events outside the codon-anticodon interaction itself.

Codon context effects in missense suppression

After our first observation of codon context effects in missense suppression ( Murgola & Pagel , 1983), we measured the suppression of missense mutations at two positions in trpA in Escherichia coli. The suppressible codons in the trpA messenger RNA were the lysine codons, AAA and AAG, and the glutamic acid codons, GAA and GAG. The mRNA sites of the codons correspond to amino acids 211 and 234 of the trpA polypeptide, positions at which glycine is the wild-type amino acid. Our data demonstrated codon context effects with both pairs of codons. The results indicate that suppression of AAA and AAG by mutant lysine transfer RNAs was more efficient at 211 than at 234, whereas suppression of GAA and GAG by two different mutant glycine tRNAs was more efficient at 234 than at 211. In general, the context effects were more pronounced with GAG and AAG than with GAA and AAA. (In some instances it appeared that suppression of GAA or AAA at a given position was more effective than suppression of GAG or AAG.) By contrast, no context effects were observed with a glyT suppressor of AAA and AAG, a glyT GAA/G-suppressor, and a glyU suppressor of GAG. Our observation of this phenomenon in missense suppression demonstrates that codon context can affect polypeptide elongation and that the effects can be different depending on the codons and tRNAs examined. It is suggested that tRNA-tRNA interaction on the ribosome is involved in the observed context effects.