Codon usage determines translation rate in Escherichia coli

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

Expression and purification of the full-length N-acetylgalactosaminyltransferase and galactosyltransferase from Campylobacter jejuni in Escherichia coli

The successful enzymatic synthesis of various ganglioside-related oligosaccharides requires many available glycan-processing enzymes. However, the number of available glycan-processing enzymes remains limited. In this study, the full-length CgtA43456 (β-(1→4)-N-acetylgalactosaminyltransferase) and CgtB11168 (β-(1→3)-galactosyltransferase) were successfully produced from Escherichia coli through the optimization of E. coli-preferable codon usage, selection of E. coli strain, and use of the molecular chaperone GroEL-GroES (GroEL/ES). The CgtA43456 enzyme was produced as a soluble form in E. coli C41(DE3) co-expressed with codon-optimized CgtA43456 and GroEL/ES. However, soluble CgtB11168 was well expressed in E. coli C41(DE3) with only the codon-optimized CgtB11168. Rather, when co-expressed with GroEL/ES, total production of CgtB11168 was reduced. Using immobilized-metal affinity chromatography, the CgtA43456 and CgtB11168 proteins were obtained with approximately 75-78 % purity. The purified CgtA43456 showed a specific activity of 21 mU/mg using UDP-N-acetylgalactosamine and GM3 trisaccharide as donor and acceptor, respectively. The purified CgtB11168 catalyzed the transfer of galactose from UDP-Gal to GM2 tetrasaccharide with a specific activity of 16 mU/mg. We propose that they could be used as catalysts for enzymatic synthesis of GM1 ganglioside-related oligosaccharides.

Dissimilation of synonymous codon usage bias in virus-host coevolution due to translational selection

Eighteen of the 20 amino acids are each encoded by more than one synonymous codon. Due to differential transfer RNA supply within the cell, synonymous codons are not used with equal frequency, a phenomenon termed codon usage bias (CUB). Previous studies have demonstrated that CUB of endogenous genes trans-regulates the translational efficiency of other genes. We hypothesized similar effects for CUB of exogenous genes on host translation, and tested it in the case of viral infection, a common form of naturally occurring exogenous gene translation. We analysed public Ribo-Seq datasets from virus-infected yeast and human cells and showed that virus CUB trans-regulated tRNA availability, and therefore the relative decoding time of codons. Manipulative experiments in yeast using 37 synonymous fluorescent proteins confirmed that an exogenous gene with CUB more similar to that of the host would apply decreased translational load on the host per unit of expression, whereas expression of the exogenous gene was elevated. The combination of these two effects was that exogenous genes with CUB overly similar to that of the host severely impeded host translation. Finally, using a manually curated list of viruses and natural and symptomatic hosts, we found that virus CUB tended to be more similar to that of symptomatic hosts than that of natural hosts, supporting a general deleterious effect of excessive CUB similarity between virus and host. Our work revealed repulsion between virus and host CUBs when they are overly similar, a previously unrecognized complexity in the coevolution of virus and host.

Bacterial stress defense: the crucial role of ribosome speed

In nature, bacteria are constantly adapting to various stressful conditions. Timely activation of stress response programs is crucial for bacteria to smoothly survive under stressful conditions. Stress response, demanding the de novo synthesis of many defense proteins, is generally activated at the transcriptional level by specific regulators. However, the effect of the global protein translational status on stress response has been largely overlooked. The translational capacity is limited by the number of translating ribosomes and the translational elongation rate. Recent work has shown that certain environmental stressors (e.g. oxidative stress) could severely compromise the stress response progress of bacteria by causing either slow-down or even complete stalling of the translational elongation process. The maintenance of ribosome elongation rate, being crucial for timely synthesis of stress defense proteins, becomes the physiological bottleneck that limits the survival of bacteria in some stressful conditions. Here, we briefly summarize some recent progress on the translational status of bacteria under two distinct stress conditions, nutrient deprivation and oxidative stress. We further discuss several important open questions on the translational regulation of bacteria during stress. The ribosome translation should be investigated in parallel with traditional transcriptional regulation in order to gain a better understanding on bacterial stress defense.

Cloning and expression of in silico modeled protein enriched with branched chain amino acids in Pichia pastoris

We have earlier in silico designed the 3-dimensional structure of a protein enriched with branched chain amino acids (BCAA, 56.4%), having only α-helical coiled-coil structure. Here, homology modeling was used to improve the in silico designed protein model. The secondary and tertiary structures of improved protein model were predicted, and validated using various online bioinformatics tools. The amino acid sequence of the final predicted Protein Model-51 was EQLTKLEIVIRVLKLLKLIGGLVSLVEWVLTALVTLLGDKVLDDILTDVIMLVKKIL DKVIGIVYVLAILALILSEVLDILWLLEKLVEILEGHHHHHH. The amino acid sequence of the protein model was reverse translated to DNA sequence and codons were optimized using codon optimization tool. The chemically synthesized BCAA51 gene was cloned to pPICZαC vector, and transformed into DH5α E. coli strain. After successful transformation, the protein was expressed in P. pastoris system by inducing with 0.5% methanol, every 24 h for up to 144 h. The expressed protein was purified by His Select Nickel affinity chromatography with an yield of 1.412 mg/L. The recombinant protein was confirmed by SDS-PAGE and western blot analysis, which showed a clear band at the expected molecular weight of ~11 kDa. Thus, here we have shown that the in silico designed protein is successfully cloned and expressed in P. pastoris.

A rare codon-based translational program of cell proliferation

BACKGROUND

The speed of translation elongation is primarily determined by the abundance of tRNAs. Thus, the codon usage influences the rate with which individual mRNAs are translated. As the nature of tRNA pools and modifications can vary across biological conditions, codon elongation rates may also vary, leading to fluctuations in the protein production from individual mRNAs. Although it has been observed that functionally related mRNAs exhibit similar codon usage, presumably to provide an effective way to coordinate expression of multiple proteins, experimental evidence for codon-mediated translation efficiency modulation of functionally related mRNAs in specific conditions is scarce and the associated mechanisms are still debated.

RESULTS

Here, we reveal that mRNAs whose expression increases during cell proliferation are enriched in rare codons, poorly adapted to tRNA pools. Ribosome occupancy profiling and proteomics measurements show that upon increased cell proliferation, transcripts enriched in rare codons undergo a higher translation boost than transcripts with common codons. Re-coding of a fluorescent reporter with rare codons increased protein output by ~ 30% relative to a reporter re-coded with common codons. Although the translation capacity of proliferating cells was higher compared to resting cells, we did not find evidence for the regulation of individual tRNAs. Among the models that were proposed so far to account for codon-mediated translational regulation upon changing conditions, the one that seems most consistent with our data involves a global upregulation of ready-to-translate tRNAs, which we show can lead to a higher increase in the elongation velocity at rare codons compared to common codons.

CONCLUSIONS

We propose that the alleviation of translation bottlenecks in rapidly dividing cells enables preferential upregulation of pro-proliferation proteins, encoded by mRNAs that are enriched in rare codons.

Grad-cryo-EM: Tool to Isolate Translation Initiation Complexes from Rabbit Reticulocyte Lysate Suitable for Structural Studies

Since its development, single-particle cryogenic electron microscopy (cryo-EM) has played a central role in the study at medium resolution of both bacterial and eukaryotic ribosomal complexes. With the advent of the direct electron detectors and new processing software which allow obtaining structures at atomic resolution, formerly obtained only by X-ray crystallography, cryo-EM has become the method of choice for the structural analysis of the translation machinery. In most of the cases, the ribosomal complexes at different stages of the translation process are assembled in vitro from purified components, which limit the analysis to previously well-characterized complexes with known factors composition. The initiation phase of the protein synthesis is a very dynamic process during which several proteins interact with the translation apparatus leading to the formation of a chronological series of initiation complexes (ICs). Here we describe a method to isolate ICs assembled on natural in vitro transcribed mRNA directly from rabbit reticulocyte lysate (RRL) by sucrose density gradient centrifugation . The Grad-cryo-EM approach allows investigating structures and composition of intermediate ribosomal complexes prepared in near-native condition by cryo-EM and mass spectrometry analyses. This is a powerful approach, which could be used to study translation initiation of any mRNAs, including IRES containing ones, and which could be adapted to different cell extracts.

Deep sequencing of tRNA's 3'-termini sheds light on CCA-tail integrity and maturation

The 3'-termini of tRNA are the point of amino acid linkage and thus crucial for their function in delivering amino acids to the ribosome and other enzymes. Therefore, to provide tRNA functionality, cells have to ensure the integrity of the 3'-terminal CCA-tail, which is generated during maturation by the 3'-trailer processing machinery and maintained by the CCA-adding enzyme. We developed a new tRNA sequencing method that is specifically tailored to assess the 3'-termini of E. coli tRNA. Intriguingly, we found a significant fraction of tRNAs with damaged CCA-tails under exponential growth conditions and, surprisingly, this fraction decreased upon transition into stationary phase. Interestingly, tRNAs bearing guanine as a discriminator base are generally unaffected by CCA-tail damage. In addition, we showed tRNA species-specific 3'-trailer processing patterns and reproduced in vitro findings on preferences of the maturation enzyme RNase T in vivo.

Elucidation of Codon Usage Signatures across the Domains of Life

Because of the degeneracy of the genetic code, multiple codons are translated into the same amino acid. Despite being “synonymous,” these codons are not equally used. Selective pressures are thought to drive the choice among synonymous codons within a genome, while GC content, which is typically attributed to mutational drift, is the major determinant of variation across species. Here, we find that in addition to GC content, interspecies codon usage signatures can also be detected. More specifically, we show that a single amino acid, arginine, is the major contributor to codon usage bias differences across domains of life. We then exploit this finding and show that domain-specific codon bias signatures can be used to classify a given sequence into its corresponding domain of life with high accuracy. We then wondered whether the inclusion of codon usage codon autocorrelation patterns, which reflects the nonrandom distribution of codon occurrences throughout a transcript, might improve the classification performance of our algorithm. However, we find that autocorrelation patterns are not domain-specific, and surprisingly, are unrelated to tRNA reusage, in contrast to previous reports. Instead, our results suggest that codon autocorrelation patterns are a by-product of codon optimality throughout a sequence, where highly expressed genes display autocorrelated “optimal” codons, whereas lowly expressed genes display autocorrelated “nonoptimal” codons.

In silicoEvaluation of Substrate Binding Site and Rare Codons in the Structure of CYP152A1

Background:

The Cytochromes P450 (CYPs) have an essential role in the oxidation of endogenous and exogenous molecules. The CYPs are identified in all domains of life, but the CYP152A1 from Bacillus subtilis is specially considered for clinical and industrial applications. The molecular cloning of a new type of CYP from Bacillus subtilis was reported, previously. Here, we describe the hidden layer of biological information of the CYP152A1 enzyme, which can help researchers for better understanding of enzyme application. In this study, four rare codons of enzyme, including Arg63, Arg187, Arg276, and Arg338 were identified and evaluated using the bioinformatics web servers.

Methods:

Through in silico modeling of CYP152A1 via the I-TASSER server, the above-mentioned rare codons were studied in the structure of enzyme that may have an important role in the proper folding of CYP152A1. In the following, the substrate binding site of CYP152A1 was studied by AutoDock Vina, and the heme and palmitic acid were considered as the substrates.

Results:

The results of docking study elucidated the Arg242 in the active site is closely related to the substrate binding site of CYP152A1, which help us to further clarify the mechanism of the enzyme reaction.

Conclusion:

Studies of these hidden information’s can enhance our understanding of CYP152A1 folding and protein expression challenges. Moreover, identification of rare codons can help in the rational design of new and effective drugs.

Computational resources and strategies to assess single-molecule dynamics of the translation process in S. cerevisiae

This work provides a systematic and comprehensive overview of available resources for the molecular-scale modelling of the translation process through agent-based modelling. The case study is the translation in Saccharomyces cerevisiae, one of the most studied yeasts. The data curation workflow encompassed structural information about the yeast (i.e. the simulation environment), and the proteins, ribonucleic acids and other types of molecules involved in the process (i.e. the agents). Moreover, it covers the main process events, such as diffusion (i.e. motion of molecules in the environment) and collision efficiency (i.e. interaction between molecules). Data previously determined by wet-lab techniques were preferred, resorting to computational predictions/extrapolations only when strictly necessary. The computational modelling of the translation processes is of added industrial interest, since it may bring forward knowledge on how to control such phenomena and enhance the production of proteins of interest in a faster and more efficient manner.

Power series method for solving TASEP-based models of mRNA translation

We develop a method for solving mathematical models of messenger RNA (mRNA) translation based on the totally asymmetric simple exclusion process (TASEP). Our main goal is to demonstrate that the method is versatile and applicable to realistic models of translation. To this end we consider the TASEP with codon-dependent elongation rates, premature termination due to ribosome drop-off and translation reinitiation due to circularisation of the mRNA. We apply the method to the model organism Saccharomyces cerevisiae under physiological conditions and find an excellent agreement with the results of stochastic simulations. Our findings suggest that the common view on translation as being rate-limited by initiation is oversimplistic. Instead we find theoretical evidence for ribosome interference and also theoretical support for the ramp hypothesis which argues that codons at the beginning of genes have slower elongation rates in order to reduce ribosome density and jamming.

On the role of mRNA secondary structure in bacterial translation

Messenger RNA (mRNA) is no longer considered as a mere informational molecule whose sole function is to convey the genetic information specified by DNA to the ribosome. Beyond this primary function, mRNA also contains additional instructions that influence the way and the extent to which this message is translated by the ribosome into protein(s). Indeed, owing to its intrinsic propensity to quickly and dynamically fold and form higher order structures, mRNA exhibits a second layer of structural information specified by the sequence itself. Besides influencing transcription and mRNA stability, this additional information also affects translation, and more precisely the frequency of translation initiation, the choice of open reading frame by recoding, the elongation speed, and the folding of the nascent protein. Many studies in bacteria have shown that mRNA secondary structure participates to the rapid adaptation of these versatile organisms to changing environmental conditions by efficiently tuning translation in response to diverse signals, such as the presence of ligands, regulatory proteins, or small RNAs. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Translation Regulation.

A Single Synonymous Variant (c.354G>A [p.P118P]) in ADAMTS13 Confers Enhanced Specific Activity

Synonymous variants within coding regions may influence protein expression and function. We have previously reported increased protein expression levels ex vivo (~120% in comparison to wild-type) from a synonymous polymorphism variant, c.354G>A [p.P118P], of the ADAMTS13 gene, encoding a plasma protease responsible for von Willebrand Factor (VWF) degradation. In the current study, we investigated the potential mechanism(s) behind the increased protein expression levels from this variant and its effect on ADAMTS13 physico-chemical properties. Cell-free assays showed enhanced translation of the c.354G>A variant and the analysis of codon usage characteristics suggested that introduction of the frequently used codon/codon pair(s) may have been potentially responsible for this effect. Limited proteolysis, however, showed no substantial influence of altered translation on protein conformation. Analysis of post-translational modifications also showed no notable differences but identified three previously unreported glycosylation markers. Despite these similarities, p.P118P variant unexpectedly showed higher specific activity. Structural analysis using modeled interactions indicated that subtle conformational changes arising from altered translation kinetics could affect interactions between an exosite of ADAMTS13 and VWF resulting in altered specific activity. This report highlights how a single synonymous nucleotide variation can impact cellular expression and specific activity in the absence of measurable impact on protein structure.

Unraveling the Role of Silent Mutation in the ω-Subunit of Escherichia coli RNA Polymerase: Structure Transition Inhibits Transcription

The bacterial RNA polymerase is a multi-subunit enzyme complex composed of six subunits, α2ββ'σω. The function of this enzyme is to transcribe the DNA base sequence to the RNA intermediate, which is ultimately translated to protein. Though the contribution of each subunit in RNA synthesis has been clearly elucidated, the role of the smallest ω-subunit is still unclear despite several studies. Recently, a study on a dominant negative mutant of rpoZ has been reported in which the mutant was shown to render the RNA polymerase defective in transcription initiation (ω6, N60D) and gave an insight on the function of ω in RNA polymerase. Serendipitously, we also obtained a silent mutant, and the mutant was found to be lethal during the isolation of toxic mutants. The primary focus of this study is to understand the mechanistic details of this lethality. Isolated ω shows a predominantly unstructured circular dichroism profile and becomes α-helical in the enzyme complex. This structural transition is perhaps the reason for this lack of function. Subsequently, we generated several silent mutants of ω to investigate the role of codon bias and the effect of rare codons with respect to their position in rpoZ. Not all silent mutations affect the structure. RNA polymerase when reconstituted with structurally altered silent mutants of ω is transcriptionally inactive. The CodonPlus strain, which has surplus tRNA, was used to assess for the rescue of the phenotype in lethal silent mutants.

Parsing the synonymous mutations in the maize genome: isoaccepting mutations are more advantageous in regions with codon co-occurrence bias

BACKGROUND

Synonymous mutations do not change amino acids but do sometimes change the tRNAs (anticodons) that decode a particular codon. An isoaccepting codon is a synonymous codon that shares the same tRNA. If a mutated codon could base pair with the same anticodon as the original, the mutation is termed an isoaccepting mutation. An interesting but less-studied type of codon bias is codon co-occurrence bias. There is a trend to cluster the isoaccepting codons in the genome. The proposed advantage of codon co-occurrence bias is that the tRNA released from the ribosome E site could be quickly recharged and subsequently decode the following isoaccepting codons. This advantage would enhance translation efficiency. In plant species, whether there are signals of positive selection on isoaccepting mutations in the codon co-occurred regions has not been studied.

RESULTS

We termed polymorphic mutations in coding regions using publicly available RNA-seq data in maize (Zea mays). Next, we classified all synonymous mutations into three categories according to the context, i.e., the relationship between the focal codon and the previous codon, as follows: isoaccepting, nonisoaccepting and nonsynonymous. We observed higher fractions of isoaccepting mutations in the isoaccepting context. If we looked at the minor allele frequency (MAF) spectrum, the isoaccepting mutations have a higher MAF in the isoaccepting context than that in other regions, and accordingly, the nonisoaccepting mutations have a higher MAF in the nonisoaccepting context.

CONCLUSION

Our results indicate that in regions with codon co-occurrence bias, natural selection maintains this pattern by suppressing the nonisoaccepting mutations. However, if the consecutive codons are nonisoaccepting, mutations tend to switch these codons to become isoaccepting. Our study demonstrates that the codon co-occurrence bias in the maize genome is selectively maintained by natural selection and that the advantage of this trend could potentially be the rapid recharging and reuse of tRNAs to increase translation efficiency.

Predicting Functional Effects of Synonymous Variants: A Systematic Review and Perspectives

Recent advances in high-throughput experimentation have put the exploration of genome sequences at the forefront of precision medicine. In an effort to interpret the sequencing data, numerous computational methods have been developed for evaluating the effects of genome variants. Interestingly, despite the fact that every person has as many synonymous (sSNV) as non-synonymous single nucleotide variants, our ability to predict their effects is limited. The paucity of experimentally tested sSNV effects appears to be the limiting factor in development of such methods. Here, we summarize the details and evaluate the performance of nine existing computational methods capable of predicting sSNV effects. We used a set of observed and artificially generated variants to approximate large scale performance expectations of these tools. We note that the distribution of these variants across amino acid and codon types suggests purifying evolutionary selection retaining generated variants out of the observed set; i.e., we expect the generated set to be enriched for deleterious variants. Closer inspection of the relationship between the observed variant frequencies and the associated prediction scores identifies predictor-specific scoring thresholds of reliable effect predictions. Notably, across all predictors, the variants scoring above these thresholds were significantly more often generated than observed. which confirms our assumption that the generated set is enriched for deleterious variants. Finally, we find that while the methods differ in their ability to identify severe sSNV effects, no predictor appears capable of definitively recognizing subtle effects of such variants on a large scale.

Characterizing the heat response of Arabidopsis thaliana from the perspective of codon usage bias and translational regulation

mRNA translation is carefully regulated at both the initiation and the elongation step. Under heat stress, it is known that particular genes change their expression and translation levels to respond to the environment. Attention has been paid to the detailed mechanisms of how a few proteins work, and little is done to analyze whether the global evolutionary patterns affect the translational changes. Determinants like codon usage bias and its related evolutionary features are less studied in heat stress experiments, especially for plants. Utilizing the RNA-seq and Ribo-seq data of normal and heat-stressed Arabidopsis thaliana generated from a previous study, we conducted gene-level (global) and codon-resolution (local) translation analyses. We studied how codon usage bias and other evolutionary features could impact the translation patterns in the heat response of the plant. We found that the evolutionary features including codon usage bias, tAI, nitrogen cost, and conservation (identity) could affect the global and local translation efficiency. Under heat stress, the optimal and conserved codons are more likely to alter their local translation elongation speed to modulate the global translation of host genes. Meanwhile, we also verified the widely accepted notions that the secondary structures and proline codons could largely slow down the translation rate. Our results revealed the effect of codon usage bias and other evolutionary patterns on the translation regulation under heat stress. Unveiling the effect of these features on translational regulation of plants might be helpful in understanding the relationship and interaction between plants and the environment.

Nature and Regulation of Protein Folding on the Ribosome

Co-translational protein folding is an essential process by which cells ensure the safe and efficient production and assembly of new proteins in their functional native states following biosynthesis on the ribosome. In this review, we describe recent progress in probing the changes during protein synthesis of the free energy landscapes that underlie co-translational folding and discuss the critical coupling between these landscapes and the rate of translation that ultimately determines the success or otherwise of the folding process. Recent developments have revealed a variety of mechanisms by which both folding and translation can be modulated or regulated, and we discuss how these effects are utilised by the cell to optimise the outcome of protein biosynthesis.

Evolutionary Changes after Translational Challenges Imposed by Horizontal Gene Transfer

Genes acquired by horizontal gene transfer (HGT) may provide the recipient organism with potentially new functions, but proper expression level and integration of the transferred genes in the novel environment are not granted. Notably, transferred genes can differ from the receiving genome in codon usage preferences, leading to impaired translation and reduced functionality. Here, we characterize the genomic and proteomic changes undergone during experimental evolution of Escherichia coli after HGT of three synonymous versions, presenting very different codon usage preference, of an antibiotic resistance gene. The experimental evolution was conducted with and without the corresponding antibiotic and the mutational patterns and proteomic profiles after 1,000 generations largely depend on the experimental growth conditions (e.g., mutations in antibiotic off-target genes), and on the synonymous gene version transferred (e.g., mutations in genes responsive to translational stress). The transfer of an exogenous gene extensively modifies the whole proteome, and these proteomic changes are different for the different version of the transferred gene. Additionally, we identified conspicuous changes in global regulators and in intermediate metabolism, confirmed the evolutionary ratchet generated by mutations in DNA repair genes and highlighted the plasticity of bacterial genomes accumulating large and occasionally transient duplications. Our results support a central role of HGT in fuelling evolution as a powerful mechanism promoting rapid, often dramatic genotypic and phenotypic changes. The profound reshaping of the pre-existing geno/phenotype allows the recipient bacteria to explore new ways of functioning, far beyond the mere acquisition of a novel function.

Re-introducing non-optimal synonymous codons into codon-optimized constructs enhances soluble recovery of recombinant proteins from Escherichia coli

Gene synthesis services have largely superseded traditional PCR methods for the generation of cDNAs destined for bacterial expression vectors. This, in turn, has increased the application of codon-optimized cDNAs where codons rarely used by Escherchia coli are replaced with common synonymous codons to accelerate translation of the target. A markedly accelerated rate of expression often results in a significant uplift in the levels of target protein but a substantial proportion of the enhanced yield can partition to the insoluble fraction rendering a significant portion of the gains unavailable for native purification. We have assessed several expression attenuation strategies for their utility in the manipulation of the soluble fraction towards higher levels of soluble target recovery from codon optimized systems. Using a set of human small GTPases as a case study, we compare the degeneration of the T7 promoter sequence, the use of alternative translational start codons and the manipulation of synonymous codon usage. Degeneration of both the T7 promoter and the translational start codon merely depressed overall expression and did not increase the percentage of product recovered in native purification of the soluble fraction. However, the selective introduction of rare non-optimal codons back into the codon-optimized sequence resulted in significantly elevated recovery of soluble targets. We propose that slowing the rate of extension during translation using a small number of rare codons allows more time for the co-translational folding of the nascent polypeptide. This increases the proportion of the target recovered in the soluble fraction by immobilized metal affinity chromatography (IMAC). Thus, a "de-optimization" of codon-optimized cDNAs, to attenuate or pause the translation process, may prove a useful strategy for improved recombinant protein production.

Nonsynonymous, synonymous and nonsense mutations in human cancer-related genes undergo stronger purifying selections than expectation

Background

Nonsynonymous mutations change the protein sequences and are frequently subjected to natural selection. The same goes for nonsense mutations that introduce pre-mature stop codons into CDSs (coding sequences). Synonymous mutations, however, are intuitively thought to be functionally silent and evolutionarily neutral. Now researchers know that the optimized synonymous codon usage is advantageous in the speedy mRNA translation process. With the advent of NGS technique, the explosion of NGS data generated from the tumor tissues help researchers identify driver mutations in cancer-related genes, but relatively less attention is paid to the SNP data in healthy human populations when studying cancer.

Methods

Here, we analyzed the publically available human SNPs. We classified these SNPs according to their functional and evolutionary categories. By simply dividing the human genes into cancer-related genes and other genes, we compared the features of nonsynonymous, synonymous and nonsense mutations in these two gene sets from multiple aspects.

Results

We provided lines of evidence that the nonsynonymous, synonymous and nonsense mutations in cancer-related genes undergo stronger purifying selection when compared to the expected pattern in other genes. The lower nonsynonymous to synonymous ratio observed in cancer-related genes suggests the suppression of amino acid substitutions in these genes. The synonymous SNPs, after excluding those in splicing regions, exhibit preferred changes in codon usage and higher codon frequencies in cancer-related genes compared to other genes, indicating the constraint exerted on these mutations. Nonsense mutations are less frequent and located closer to stop codons in cancer-related genes than in other genes, which putatively minimize their deleterious effects.

Conclusion

Our study demonstrated the evolutionary constraint on mutations in CDS of cancer-related genes without the requirement of data from cancer tissues or patients. Our work provides novel perspectives on interpreting the constraint on mutations in cancer-related genes. We reveal extra constraint on synonymous mutations in cancer-related genes which is related to codon usage bias and is in addition to the splicing effect.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5572-x) contains supplementary material, which is available to authorized users.

Novel insights into gene expression regulation during meiosis revealed by translation elongation dynamics

The ability to dynamically control mRNA translation has a great impact on many intracellular processes. Whereas it is believed that translational control in eukaryotes occurs mainly at initiation, the condition-specific changes at the elongation level and their potential regulatory role remain unclear. Using computational approaches applied to ribosome profiling data, we show that elongation rate is dynamic and can change considerably during the yeast meiosis to facilitate the selective translation of stage-specific transcripts. We observed unique elongation changes during meiosis II, including a global inhibition of translation elongation at the onset of anaphase II accompanied by a sharp shift toward increased elongation for genes required at this meiotic stage. We also show that ribosomal proteins counteract the global decreased elongation by maintaining high initiation rates. Our findings provide new insights into gene expression regulation during meiosis and demonstrate that codon usage evolved, among others, to optimize timely translation.

An improved estimation of tRNA expression to better elucidate the coevolution between tRNA abundance and codon usage in bacteria

The degree to which codon usage can be explained by tRNA abundance in bacterial species is often inadequate, partly because differential tRNA abundance is often approximated by tRNA copy numbers. To better understand the coevolution between tRNA abundance and codon usage, we provide a better estimate of tRNA abundance by profiling tRNA mapped reads (tRNA tpm) using publicly available RNA Sequencing data. To emphasize the feasibility of our approach, we demonstrate that tRNA tpm is consistent with tRNA abundances derived from RNA fingerprinting experiments in Escherichia coli, Bacillus subtilis, and Salmonella enterica. Furthermore, we do not observe an appreciable reduction in tRNA sequencing efficiency due to post-transcriptional methylations in the seven bacteria studied. To determine optimal codons, we calculate codon usage in highly and lowly expressed genes determined by protein per transcript. We found that tRNA tpm is sensitive to identify more translationally optimal codons than gene copy number and early tRNA fingerprinting abundances. Additionally, tRNA tpm improves the predictive power of tRNA adaptation index over codon preference. Our results suggest that dependence of codon usage on tRNA availability is not always associated with species growth-rate. Conversely, tRNA availability is better optimized to codon usage in fast-growing than slow-growing species.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

The fitness consequences of synonymous mutations in Escherichia coli: Experimental evidence for a pleiotropic effect of translational selection

Codon usage bias (CUB) is a universal feature of genomes, and in most species CUB of protein coding genes is positively correlated with expression level and degree of evolutionary conservation. There is mounting experimental evidence that CUB is due in part to selection for translational efficiency and/or accuracy, i.e., translational selection. However, there is a paucity of experimental data on whether and how CUB acts in trans - does the usage of preferred codons in a highly expressed gene affect the translation of other genes by freeing up more ribosomes, thereby increasing their availability to translate all mRNA transcripts in the cell? We investigated this question by creating two extreme versions of the highly expressed Escherichia coli β-lactamase (bla) gene, one comprised almost entirely of unpreferred codons, and a second comprised almost entirely of preferred codons. We monitored the fitness effects of these synonymous mutations over hundreds of generations in two selective environments that allowed us to disentangle translational effects acting in cis from those acting in trans. In a selective environment for maximizing translational efficiency in trans of a gene (tetA) encoding a tetracycline resistance protein, unpreferred synonymous mutations had a negative impact on long-term fitness, whereas preferred mutations had a positive impact on long-term fitness, providing strong experimental evidence for a pleiotropic effect of translational selection.

The codon sequences predict protein lifetimes and other parameters of the protein life cycle in the mouse brain

The homeostasis of the proteome depends on the tight regulation of the mRNA and protein abundances, of the translation rates, and of the protein lifetimes. Results from several studies on prokaryotes or eukaryotic cell cultures have suggested that protein homeostasis is connected to, and perhaps regulated by, the protein and the codon sequences. However, this has been little investigated for mammals in vivo. Moreover, the link between the coding sequences and one critical parameter, the protein lifetime, has remained largely unexplored, both in vivo and in vitro. We tested this in the mouse brain, and found that the percentages of amino acids and codons in the sequences could predict all of the homeostasis parameters with a precision approaching experimental measurements. A key predictive element was the wobble nucleotide. G-/C-ending codons correlated with higher protein lifetimes, protein abundances, mRNA abundances and translation rates than A-/U-ending codons. Modifying the proportions of G-/C-ending codons could tune these parameters in cell cultures, in a proof-of-principle experiment. We suggest that the coding sequences are strongly linked to protein homeostasis in vivo, albeit it still remains to be determined whether this relation is causal in nature.

Hepatitis A Virus Codon Usage: Implications for Translation Kinetics and Capsid Folding

Codon usage bias is universal to all genomes. Hepatitis A virus (HAV) codon usage is highly biased and deoptimized with respect to its host. Accordingly, HAV is unable to induce cellular translational shutoff and its internal ribosome entry site (IRES) is inefficient. Codon usage deoptimization may be seen as a hawk (host cell) versus dove (HAV) game strategy for accessing transfer RNA (tRNA). HAV avoids use of abundant host cell codons and thereby eludes competition for the corresponding tRNAs. Instead, codons that are abundant or rare in cellular messenger RNAs (mRNAs) are used relatively rarely in its genome, although intermediately abundant host cell codons are abundant in the viral genome. Rare codons in the capsid coding region slow down the translation elongation rate, and in doing so intrinsically modulate capsid folding, which is critical to the stability of a virus transmitted through the fecal-oral route. HAV is a paradigmatic example of what has been proposed as a codon usage "code" for protein structure.

Evaluation of GFP reporter utility for analysis of transcriptional slippage during gene expression

Background

Epimutations arising from transcriptional slippage seem to have more important role in regulating gene expression than earlier though. Since the level and the fidelity of transcription primarily determine the overall efficiency of gene expression, all factors contributing to their decrease should be identified and optimized.

Results

To examine the influence of A/T homopolymeric sequences on introduction of erroneous nucleotides by slippage mechanism green fluorescence protein (GFP) reporter was chosen. The in- or out-of-frame gfp gene was fused to upstream fragment with variable number of adenine or thymine stretches resulting in several hybrid GFP proteins with diverse amino acids at N-terminus. Here, by using T7 phage expression system we showed that the intensity of GFP fluorescence mainly depends on the number of the retained natural amino acids. While the lack of serine (S₂) residue results in negligible effects, the lack of serine and lysine (S₂K₃) contributed to a significant reduction in fluorescence by 2.7-fold for polyA-based in-frame controls and twofold for polyTs. What is more, N-terminal tails amino acid composition was rather of secondary importance, since the whole-cell fluorescence differed in a range of 9–18% between corresponding polyA- and polyT-based constructs.

Conclusions

Here we present experimental evidence for utility of GFP reporter for accurate estimation of A/T homopolymeric sequence contribution in transcriptional slippage induction. We showed that the intensity of GFP hybrid fluorescence mainly depends on the number of retained natural amino acids, thus fluorescence raw data need to be referred to appropriate positive control. Moreover, only in case of GFP hybrids with relatively short N-terminal tags the fluorescence level solely reflects production yield, what further indicates the impact of an individual slippage sequence. Our results demonstrate that in contrast to the E. coli enzyme, T7 RNA polymerase exhibits extremely high propensity to slippage even on runs as short as 3 adenine or 4 thymine residues.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0999-3) contains supplementary material, which is available to authorized users.

Codon usage influences fitness through RNA toxicity

Many organisms are subject to selective pressure that gives rise to unequal usage of synonymous codons, known as codon bias. To experimentally dissect the mechanisms of selection on synonymous sites, we expressed several hundred synonymous variants of the GFP gene in Escherichia coli, and used quantitative growth and viability assays to estimate bacterial fitness. Unexpectedly, we found many synonymous variants whose expression was toxic to E. coli Unlike previously studied effects of synonymous mutations, the effect that we discovered is independent of translation, but it depends on the production of toxic mRNA molecules. We identified RNA sequence determinants of toxicity and evolved suppressor strains that can tolerate the expression of toxic GFP variants. Genome sequencing of these suppressor strains revealed a cluster of promoter mutations that prevented toxicity by reducing mRNA levels. We conclude that translation-independent RNA toxicity is a previously unrecognized obstacle in bacterial gene expression.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

Selection-driven cost-efficiency optimization of transcripts modulates gene evolutionary rate in bacteria

BACKGROUND

Most amino acids are encoded by multiple synonymous codons. However, synonymous codons are not used equally, and this biased codon use varies between different organisms. It has previously been shown that both selection acting to increase codon translational efficiency and selection acting to decrease codon biosynthetic cost contribute to differences in codon bias. However, it is unknown how these two factors interact or how they affect molecular sequence evolution.

RESULTS

Through analysis of 1320 bacterial genomes, we show that bacterial genes are subject to multi-objective selection-driven optimization of codon use. Here, selection acts to simultaneously decrease transcript biosynthetic cost and increase transcript translational efficiency, with highly expressed genes under the greatest selection. This optimization is not simply a consequence of the more translationally efficient codons being less expensive to synthesize. Instead, we show that transfer RNA gene copy number alters the cost-efficiency trade-off of synonymous codons such that, for many species, selection acting on transcript biosynthetic cost and translational efficiency act in opposition. Finally, we show that genes highly optimized to reduce cost and increase efficiency show reduced rates of synonymous and non-synonymous mutation.

CONCLUSIONS

This analysis provides a simple mechanistic explanation for variation in evolutionary rate between genes that depends on selection-driven cost-efficiency optimization of the transcript. These findings reveal how optimization of resource allocation to messenger RNA synthesis is a critical factor that determines both the evolution and composition of genes.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

The Amount of Nitrogen Used for Photosynthesis Modulates Molecular Evolution in Plants

Genome and transcript sequences are composed of long strings of nucleotide monomers (A, C, G, and T/U) that require different quantities of nitrogen atoms for biosynthesis. Here, it is shown that the strength of selection acting on transcript nitrogen content is influenced by the amount of nitrogen plants require to conduct photosynthesis. Specifically, plants that require more nitrogen to conduct photosynthesis experience stronger selection on transcript sequences to use synonymous codons that cost less nitrogen to biosynthesize. It is further shown that the strength of selection acting on transcript nitrogen cost constrains molecular sequence evolution such that genes experiencing stronger selection evolve at a slower rate. Together these findings reveal that the plant molecular clock is set by photosynthetic efficiency, and provide a mechanistic explanation for changes in plant speciation rates that occur concomitant with improvements in photosynthetic efficiency and changes in the environment such as light, temperature, and atmospheric CO2 concentration.

Expression and purification of biologically active recombinant rabbit monocyte chemoattractant protein1 in Escherichia coli

Monocyte chemoattractant protein 1 (MCP1) with recruiting monocytes is an important factor at the beginning of inflammatory disorders such as atherosclerosis which seems its blocking preclude this process and help improvement of related diseases. To perform clinical research in this field, MCP1 protein is required but firstly, animal studies should be done. As the rabbit is a suitable model for many inflammatory disorders, and Escherichia coli BL21(DE3) (BL21) cell is a high-efficiency host for protein expression, we decided to produce recombinant rabbit MCP1 (rRMCP1) in BL21/pET28a system. After codon usage, a construct containing RMCP1 sequence was synthesized, cloned into the pET28a plasmid, and overexpressed in BL21 cells. Followed that, with changing expression condition such as cell concentration before the induction, time period, temperature, shaking rate and inducer concentration (IPTG), rRMCP1 expression was optimized, and purified by Ni-NTA. The biological activity of the expressed protein was verified using monocyte migration assay. Using this expression system, nearly 28 mg/mL rRMCP1 was produced at 26°C/180 rpm for 24 h in LB broth medium with 1 mM IPTG. Therefore, we were succeeded to express the intermediate level of rRMCP1 with this method. This amount of protein is sufficient for biological researches in the laboratory.

Deciphering mRNA Sequence Determinants of Protein Production Rate

One of the greatest challenges in biophysical models of translation is to identify coding sequence features that affect the rate of translation and therefore the overall protein production in the cell. We propose an analytic method to solve a translation model based on the inhomogeneous totally asymmetric simple exclusion process, which allows us to unveil simple design principles of nucleotide sequences determining protein production rates. Our solution shows an excellent agreement when compared to numerical genome-wide simulations of S. cerevisiae transcript sequences and predicts that the first 10 codons, which is the ribosome footprint length on the mRNA, together with the value of the initiation rate, are the main determinants of protein production rate under physiological conditions. Finally, we interpret the obtained analytic results based on the evolutionary role of the codons' choice for regulating translation rates and ribosome densities.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.

Influence of multiplicative stochastic variation on translational elongation rates

Experimental data indicate that stochastic effects exerted at the level of translation contribute substantially to the variation in abundance of proteins expressed at moderate to high levels. This study analyzes the theoretical consequences of fluctuations in residue-specific elongation rates during translation. A simple analytical framework shows that rate variation during elongation gives rise to protein production rates that consist of sums of products of random variables. Simulations show that because the contribution to total variation of products of random variables greatly exceeds that of sums of random variables, the overall distribution exhibits approximately log-normal behavior. Empirical fits of the data can be satisfied by either sums of log-normal distributions, or sums of log-normal and log-logistic distributions. Elongation rate stochastic variation offers an accounting for a major component of biological variation. The analysis provided here highlights a probability distribution that is a natural extension of the Poisson and has broad applicability to many types of multiplicative noise processes.

Local genic base composition impacts protein production and cellular fitness

The maintenance of a G + C content that is higher than the mutational input to a genome provides support for the view that selection serves to increase G + C contents in bacteria. Recent experimental evidence from Escherichia coli demonstrated that selection for increasing G + C content operates at the level of translation, but the precise mechanism by which this occurs is unknown. To determine the substrate of selection, we asked whether selection on G + C content acts across all sites within a gene or is confined to particular genic regions or nucleotide positions. We systematically altered the G + C contents of the GFP gene and assayed its effects on the fitness of strains harboring each variant. Fitness differences were attributable to the base compositional variation in the terminal portion of the gene, suggesting a connection to the folding of a specific protein feature. Variants containing sequence features that are thought to result in rapid translation, such as low G + C content and high levels of codon adaptation, displayed highly reduced growth rates. Taken together, our results show that purifying selection acting against A and T mutations most likely results from their tendency to increase the rate of translation, which can perturb the dynamics of protein folding.

Unique features of nucleotide and codon usage patterns in mycoplasmas revealed by information entropy

Currently, the comparison between GC usage pattern at the 3rd codon position and codon usage index is commonly used to estimate the roles of evolutionary forces in shaping synonymous codon usages, however, this kind of analysis often losses the information about the role of A/T usage bias in shaping synonymous codon usage bias. To overcome this limitation and better understand the interplay between nucleotide and codon usages for the evolution of bacteria at gene levels, in this study, we employed the information entropy method with some modification to estimate roles of nucleotide compositions in the overall codon usage bias for 18 mycoplasma species in combination with Davies-Bouldin index. At gene levels, the overall nucleotide usage bias represents A content as the highest, followed by T, G and C for mycoplasmas, resulting in a low GC content. This feature is universal across these species derived from different hosts, suggesting that the hosts have the limited impact on nucleotide usage bias of mycoplasmas. Information entropy and Davies-Bouldin index can better reveal that the nucleotide usage bias at the 3rd codon position is essential in shaping the overall nucleotide bias for all given mycoplasmas except M. pneumoniae M129. Davies-Bouldin index revealed that the 1st and 2nd codon position play more important role in synonymous codon usage bias than that of the 3rd position at gene levels. To our knowledge, this is the first comprehensive investigation into cooperation between nucleotide and codon usages for mycoplasma and extends our knowledge of the mechanisms that contribute to codon usage and evolution of this microorganism.

Translation efficiency is maintained at elevated temperature in Escherichia coli

Cellular protein levels are dictated by the balance between gene transcription, mRNA translation, and protein degradation, among other factors. Translation requires the interplay of several RNA hybridization processes, which are expected to be temperature-sensitive. We used ribosome profiling to monitor translation in Escherichia coli at 30 °C and to investigate how this changes after 10-20 min of heat shock at 42 °C. Translation efficiencies are robustly maintained after thermal heat shock and after mimicking the heat-shock response transcriptional program at 30 °C by overexpressing the heat shock σ factor encoded by the rpoH gene. We compared translation efficiency, the ratio of ribosome footprint reads to mRNA reads for each gene, to parameters derived from gene sequences. Genes with stable mRNA structures, non-optimal codon use, and those whose gene product is cotranslationally translocated into the inner membrane are generally less highly translated than other genes. Comparison with other published datasets suggests a role for translational elongation in coupling mRNA structures to translation initiation. Genome-wide calculations of the temperature dependence of mRNA structure predict that relatively few mRNAs show a melting transition between 30 and 42 °C, consistent with the observed lack of changes in translation efficiency. We developed a linear model with six parameters that can predict 38% of the variation in translation efficiency between genes, which may be useful in interpreting transcriptome data.

Amino acid misincorporation in recombinant proteins

Proteins provide the molecular basis for cellular structure, catalytic activity, signal transduction, and molecular transport in biological systems. Recombinant protein expression is widely used to prepare and manufacture novel proteins that serve as the foundation of many biopharmaceutical products. However, protein translation bioprocesses are inherently prone to low-level errors. These sequence variants caused by amino acid misincorporation have been observed in both native and recombinant proteins. Protein sequence variants impact product quality, and their presence can be exacerbated through cellular stress, overexpression, and nutrient starvation. Therefore, the cell line selection process, which is used in the biopharmaceutical industry, is not only directed towards maximizing productivity, but also focuses on selecting clones which yield low sequence variant levels, thereby proactively avoiding potentially inauspicious patient safety and efficacy outcomes. Here, we summarize a number of hallmark studies aimed at understanding the mechanisms of amino acid misincorporation, as well as exacerbating factors, and mitigation strategies. We also describe key advances in analytical technologies in the identification and quantification of sequence variants, and some practical considerations when using LC-MS/MS for detecting sequence variants.

Codon usage regulates protein structure and function by affecting translation elongation speed in Drosophila cells

Codon usage biases are found in all eukaryotic and prokaryotic genomes and have been proposed to regulate different aspects of translation process. Codon optimality has been shown to regulate translation elongation speed in fungal systems, but its effect on translation elongation speed in animal systems is not clear. In this study, we used a Drosophila cell-free translation system to directly compare the velocity of mRNA translation elongation. Our results demonstrate that optimal synonymous codons speed up translation elongation while non-optimal codons slow down translation. In addition, codon usage regulates ribosome movement and stalling on mRNA during translation. Finally, we show that codon usage affects protein structure and function in vitro and in Drosophila cells. Together, these results suggest that the effect of codon usage on translation elongation speed is a conserved mechanism from fungi to animals that can affect protein folding in eukaryotic organisms.

Tuning of Recombinant Protein Expression in Escherichia coli by Manipulating Transcription, Translation Initiation Rates, and Incorporation of Noncanonical Amino Acids

Protein synthesis in cells has been thoroughly investigated and characterized over the past 60 years. However, some fundamental issues remain unresolved, including the reasons for genetic code redundancy and codon bias. In this study, we changed the kinetics of the Eschrichia coli transcription and translation processes by mutating the promoter and ribosome binding domains and by using genetic code expansion. The results expose a counterintuitive phenomenon, whereby an increase in the initiation rates of transcription and translation lead to a decrease in protein expression. This effect can be rescued by introducing slow translating codons into the beginning of the gene, by shortening gene length or by reducing initiation rates. On the basis of the results, we developed a biophysical model, which suggests that the density of co-transcriptional-translation plays a role in bacterial protein synthesis. These findings indicate how cells use codon bias to tune translation speed and protein synthesis.