Synonymous Codons Direct Cotranslational Folding toward Different Protein Conformations

Synonymous GATA2 mutations result in selective loss of mutated RNA and are common in patients with GATA2 deficiency

Deficiency of the transcription factor GATA2 is a highly penetrant genetic disorder predisposing to myelodysplastic syndromes (MDS) and immunodeficiency. It has been recognized as the most common cause underlying primary MDS in children. Triggered by the discovery of a recurrent synonymous GATA2 variant, we systematically investigated 911 patients with phenotype of pediatric MDS or cellular deficiencies for the presence of synonymous alterations in GATA2. In total, we identified nine individuals with five heterozygous synonymous mutations: c.351C>G, p.T117T (N = 4); c.649C>T, p.L217L; c.981G>A, p.G327G; c.1023C>T, p.A341A; and c.1416G>A, p.P472P (N = 2). They accounted for 8.2% (9/110) of cases with GATA2 deficiency in our cohort and resulted in selective loss of mutant RNA. While for the hotspot mutation (c.351C>G) a splicing error leading to RNA and protein reduction was identified, severe, likely late stage RNA loss without splicing disruption was found for other mutations. Finally, the synonymous mutations did not alter protein function or stability. In summary, synonymous GATA2 substitutions are a new common cause of GATA2 deficiency. These findings have broad implications for genetic counseling and pathogenic variant discovery in Mendelian disorders.

Human γS-Crystallin-Copper Binding Helps Buffer against Aggregation Caused by Oxidative Damage

Divalent metal cations can play a role in protein aggregation diseases, including cataract. Here we compare the aggregation of human γS-crystallin, a key structural protein of the eye lens, via mutagenesis, ultraviolet light damage, and the addition of metal ions. All three aggregation pathways result in globular, amorphous-looking structures that do not elongate into fibers. We also investigate the molecular mechanism underlying copper(II)-induced aggregation. This work was motivated by the observation that zinc(II)-induced aggregation of γS-crystallin is driven by intermolecular bridging of solvent-accessible cysteine residues, while in contrast, copper(II)-induced aggregation of this protein is exacerbated by the removal of solvent-accessible cysteines via mutation. Here we find that copper(II)-induced aggregation results from a complex mechanism involving multiple interactions with the protein. The initial protein-metal interactions result in the reduction of Cu(II) to Cu(I) with concomitant oxidation of γS-crystallin. In addition to the intermolecular disulfides that represent a starting point for aggregation, intramolecular disulfides also occur in the cysteine loop, a region of the N-terminal domain that was previously found to mediate the early stages of cataract formation. This previously unobserved ability of γS-crystallin to transfer disulfides intramolecularly suggests that it may serve as an oxidation sink for the lens after glutathione levels have become depleted during aging. γS-Crystallin thus serves as the last line of defense against oxidation in the eye lens, a result that underscores the chemical functionality of this protein, which is generally considered to play a purely structural role.

Analysis of computational codon usage models and their association with translationally slow codons

Improved computational modeling of protein translation rates, including better prediction of where translational slowdowns along an mRNA sequence may occur, is critical for understanding co-translational folding. Because codons within a synonymous codon group are translated at different rates, many computational translation models rely on analyzing synonymous codons. Some models rely on genome-wide codon usage bias (CUB), believing that globally rare and common codons are the most informative of slow and fast translation, respectively. Others use the CUB observed only in highly expressed genes, which should be under selective pressure to be translated efficiently (and whose CUB may therefore be more indicative of translation rates). No prior work has analyzed these models for their ability to predict translational slowdowns. Here, we evaluate five models for their association with slowly translated positions as denoted by two independent ribosome footprint (RFP) count experiments from S. cerevisiae, because RFP data is often considered as a “ground truth” for translation rates across mRNA sequences. We show that all five considered models strongly associate with the RFP data and therefore have potential for estimating translational slowdowns. However, we also show that there is a weak correlation between RFP counts for the same genes originating from independent experiments, even when their experimental conditions are similar. This raises concerns about the efficacy of using current RFP experimental data for estimating translation rates and highlights a potential advantage of using computational models to understand translation rates instead.

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.

Codon Harmonization of a Kir3.1-KirBac1.3 Chimera for Structural Study Optimization

The expression of functional, folded, and isotopically enriched membrane proteins is an enduring bottleneck for nuclear magnetic resonance (NMR) studies. Indeed, historically, protein yield optimization has been insufficient to allow NMR analysis of many complex Eukaryotic membrane proteins. However, recent work has found that manipulation of plasmid codons improves the odds of successful NMR-friendly protein production. In the last decade, numerous studies showed that matching codon usage patterns in recombinant gene sequences to those in the native sequence is positively correlated with increased protein yield. This phenomenon, dubbed codon harmonization, may be a powerful tool in optimizing recombinant expression of difficult-to-produce membrane proteins for structural studies. Here, we apply this technique to an inward rectifier K⁺ Channel (Kir) 3.1-KirBac1.3 chimera. Kir3.1 falls within the G protein-coupled inward rectifier K⁺ (GIRK) channel family, thus NMR studies may inform on the nuances of GIRK gating action in the presence and absence of its G Protein, lipid, and small molecule ligands. In our hands, harmonized plasmids increase protein yield nearly two-fold compared to the traditional ‘fully codon optimized’ construct. We then employ a fluorescence-based functional assay and solid-state NMR correlation spectroscopy to show the final protein product is folded and functional.

Impact of In Vivo Protein Folding Probability on Local Fitness Landscapes

It is incompletely understood how biophysical properties like protein stability impact molecular evolution and epistasis. Epistasis is defined as specific when a mutation exclusively influences the phenotypic effect of another mutation, often at physically interacting residues. In contrast, nonspecific epistasis results when a mutation is influenced by a large number of nonlocal mutations. As most mutations are pleiotropic, the in vivo folding probability-governed by basal protein stability-is thought to determine activity-enhancing mutational tolerance, implying that nonspecific epistasis is dominant. However, evidence exists for both specific and nonspecific epistasis as the prevalent factor, with limited comprehensive data sets to support either claim. Here, we use deep mutational scanning to probe how in vivo enzyme folding probability impacts local fitness landscapes. We computationally designed two different variants of the amidase AmiE with statistically indistinguishable catalytic efficiencies but lower probabilities of folding in vivo compared with wild-type. Local fitness landscapes show slight alterations among variants, with essentially the same global distribution of fitness effects. However, specific epistasis was predominant for the subset of mutations exhibiting positive sign epistasis. These mutations mapped to spatially distinct locations on AmiE near the initial mutation or proximal to the active site. Intriguingly, the majority of specific epistatic mutations were codon dependent, with different synonymous codons resulting in fitness sign reversals. Together, these results offer a nuanced view of how protein folding probability impacts local fitness landscapes and suggest that transcriptional-translational effects are as important as stability in determining evolutionary outcomes.

Molecular mechanism of translational stalling by inhibitory codon combinations and poly(A) tracts

Inhibitory codon pairs and poly(A) tracts within the translated mRNA cause ribosome stalling and reduce protein output. The molecular mechanisms that drive these stalling events, however, are still unknown. Here, we use a combination of in vitro biochemistry, ribosome profiling, and cryo-EM to define molecular mechanisms that lead to these ribosome stalls. First, we use an in vitro reconstituted yeast translation system to demonstrate that inhibitory codon pairs slow elongation rates which are partially rescued by increased tRNA concentration or by an artificial tRNA not dependent on wobble base-pairing. Ribosome profiling data extend these observations by revealing that paused ribosomes with empty A sites are enriched on these sequences. Cryo-EM structures of stalled ribosomes provide a structural explanation for the observed effects by showing decoding-incompatible conformations of mRNA in the A sites of all studied stall- and collision-inducing sequences. Interestingly, in the case of poly(A) tracts, the inhibitory conformation of the mRNA in the A site involves a nucleotide stacking array. Together, these data demonstrate a novel mRNA-induced mechanisms of translational stalling in eukaryotic ribosomes.

Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness

In the cell, proteins are synthesized from N to C terminus and begin to fold during translation. Cotranslational folding mechanisms are therefore linked to elongation rate, which varies as a function of synonymous codon usage. However, synonymous codon substitutions can affect many distinct cellular processes, which has complicated attempts to deconvolve the extent to which synonymous codon usage can promote or frustrate proper protein folding in vivo. Although previous studies have shown that some synonymous changes can lead to different final structures, other substitutions will likely be more subtle, perturbing predominantly the protein folding pathway without radically altering the final structure. Here we show that synonymous codon substitutions encoding a single essential enzyme lead to dramatically slower cell growth. These mutations do not prevent active enzyme formation; instead, they predominantly alter the protein folding mechanism, leading to enhanced degradation in vivo. These results support a model in which synonymous codon substitutions can impair cell fitness by significantly perturbing cotranslational protein folding mechanisms, despite the chaperoning provided by the cellular protein homeostasis network.

TissueCoCoPUTs: Novel Human Tissue-Specific Codon and Codon-Pair Usage Tables Based on Differential Tissue Gene Expression

Protein expression in multicellular organisms varies widely across tissues. Codon usage in the transcriptome of each tissue is derived from genomic codon usage and the relative expression level of each gene. We created a comprehensive computational resource that houses tissue-specific codon, codon-pair, and dinucleotide usage data for 51 Homo sapiens tissues (TissueCoCoPUTs: https://hive.biochemistry.gwu.edu/review/tissue_codon), using transcriptome data from the Broad Institute Genotype-Tissue Expression (GTEx) portal. Distances between tissue-specific codon and codon-pair frequencies were used to generate a dendrogram based on the unique patterns of codon and codon-pair usage in each tissue that are clearly distinct from the genomic distribution. This novel resource may be useful in unraveling the relationship between codon usage and tRNA abundance, which could be critical in determining translation kinetics and efficiency across tissues. Areas of investigation such as biotherapeutic development, tissue-specific genetic engineering, and genetic disease prediction will greatly benefit from this resource.

Cotranslational folding allows misfolding-prone proteins to circumvent deep kinetic traps

Many large proteins suffer from slow or inefficient folding in vitro. It has long been known that this problem can be alleviated in vivo if proteins start folding cotranslationally. However, the molecular mechanisms underlying this improvement have not been well established. To address this question, we use an all-atom simulation-based algorithm to compute the folding properties of various large protein domains as a function of nascent chain length. We find that for certain proteins, there exists a narrow window of lengths that confers both thermodynamic stability and fast folding kinetics. Beyond these lengths, folding is drastically slowed by nonnative interactions involving C-terminal residues. Thus, cotranslational folding is predicted to be beneficial because it allows proteins to take advantage of this optimal window of lengths and thus avoid kinetic traps. Interestingly, many of these proteins' sequences contain conserved rare codons that may slow down synthesis at this optimal window, suggesting that synthesis rates may be evolutionarily tuned to optimize folding. Using kinetic modeling, we show that under certain conditions, such a slowdown indeed improves cotranslational folding efficiency by giving these nascent chains more time to fold. In contrast, other proteins are predicted not to benefit from cotranslational folding due to a lack of significant nonnative interactions, and indeed these proteins' sequences lack conserved C-terminal rare codons. Together, these results shed light on the factors that promote proper protein folding in the cell and how biomolecular self-assembly may be optimized evolutionarily.

Cotranslational Folding of Proteins on the Ribosome

Many proteins in the cell fold cotranslationally within the restricted space of the polypeptide exit tunnel or at the surface of the ribosome. A growing body of evidence suggests that the ribosome can alter the folding trajectory in many different ways. In this review, we summarize the recent examples of how translation affects folding of single-domain, multiple-domain and oligomeric proteins. The vectorial nature of translation, the spatial constraints of the exit tunnel, and the electrostatic properties of the ribosome-nascent peptide complex define the onset of early folding events. The ribosome can facilitate protein compaction, induce the formation of intermediates that are not observed in solution, or delay the onset of folding. Examples of single-domain proteins suggest that early compaction events can define the folding pathway for some types of domain structures. Folding of multi-domain proteins proceeds in a domain-wise fashion, with each domain having its role in stabilizing or destabilizing neighboring domains. Finally, the assembly of protein complexes can also begin cotranslationally. In all these cases, the ribosome helps the nascent protein to attain a native fold and avoid the kinetic traps of misfolding.

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.

Towards next generation therapies for cystic fibrosis: Folding, function and pharmacology of CFTR

The treatment of cystic fibrosis (CF) has been transformed by orally-bioavailable small molecule modulators of the cystic fibrosis transmembrane conductance regulator (CFTR), which restore function to CF mutants. However, CFTR modulators are not available to all people with CF and better modulators are required to prevent disease progression. Here, we review selectively recent advances in CFTR folding, function and pharmacology. We highlight ensemble and single-molecule studies of CFTR folding, which provide new insight into CFTR assembly, its perturbation by CF mutations and rescue by CFTR modulators. We discuss species-dependent differences in the action of the F508del-CFTR mutation on CFTR expression, stability and function, which might influence pharmacological studies of CFTR modulators in CF animal models. Finally, we illuminate the identification of combinations of two CFTR potentiators (termed co-potentiators), which restore therapeutically-relevant levels of CFTR activity to rare CF mutations. Thus, mechanistic studies of CFTR folding, function and pharmacology inform the development of highly effective CFTR modulators.

Studying Plant MIF/D-DT-Like Genes and Proteins (MDLs)

Human macrophage migration inhibitory factor (MIF) is an inflammatory cytokine with chemokine-like characteristics and an upstream regulator of host innate immunity. It is a critical mediator of a variety of human diseases, such as acute and chronic inflammatory diseases, autoimmunity, atherosclerosis, and cancer. MIF is an atypical chemokine that not only signals through its cognate receptor CD74, but also interacts with the classical chemokine receptors CXCR2 and CXCR4. MIF and its homolog D-dopachrome tautomerase (D-DT)/MIF-2 are structurally unique proteins that are conserved across kingdoms and that share a remarkable homology with bacterial tautomerases/isomerases, albeit the relevance of the tautomerase activity in mammalian systems has remained unclear. Intriguingly, in silico analysis also predicts MIF orthologs in plants such as in the model plant Arabidopsis thaliana. There are three predicted MIF orthologs in A. thaliana, which have been termed A. thaliana MIF/D-DT-like proteins (AtMDLs). Anticipating that there will be a future research interest in studying AtMDLs or other plant MDLs, here we describe methods how to clone, recombinantly express and purify AtMDL proteins, taking into account codon usage differences between plant and mammalian cell systems.

Fail-safe genetic codes designed to intrinsically contain engineered organisms

One challenge in engineering organisms is taking responsibility for their behavior over many generations. Spontaneous mutations arising before or during use can impact heterologous genetic functions, disrupt system integration, or change organism phenotype. Here, we propose restructuring the genetic code itself such that point mutations in protein-coding sequences are selected against. Synthetic genetic systems so-encoded should fail more safely in response to most spontaneous mutations. We designed fail-safe codes and simulated their expected effects on the evolution of so-encoded proteins. We predict fail-safe codes supporting expression of 20 or 15 amino acids could slow protein evolution to ∼30% or 0% the rate of standard-encoded proteins, respectively. We also designed quadruplet-codon codes that should ensure all single point mutations in protein-coding sequences are selected against while maintaining expression of 20 or more amino acids. We demonstrate experimentally that a reduced set of 21 tRNAs is capable of expressing a protein encoded by only 20 sense codons, whereas a standard 64-codon encoding is not expressed. Our work suggests that biological systems using rationally depleted but otherwise natural translation systems should evolve more slowly and that such hypoevolvable organisms may be less likely to invade new niches or outcompete native populations.

eRF1 mediates codon usage effects on mRNA translation efficiency through premature termination at rare codons

Codon usage bias is a universal feature of eukaryotic and prokaryotic genomes and plays an important role in regulating gene expression levels. A major role of codon usage is thought to regulate protein expression levels by affecting mRNA translation efficiency, but the underlying mechanism is unclear. By analyzing ribosome profiling results, here we showed that codon usage regulates translation elongation rate and that rare codons are decoded more slowly than common codons in all codon families in Neurospora. Rare codons resulted in ribosome stalling in manners both dependent and independent of protein sequence context and caused premature translation termination. This mechanism was shown to be conserved in Drosophila cells. In both Neurospora and Drosophila cells, codon usage plays an important role in regulating mRNA translation efficiency. We found that the rare codon-dependent premature termination is mediated by the translation termination factor eRF1, which recognizes ribosomes stalled on rare sense codons. Silencing of eRF1 expression resulted in codon usage-dependent changes in protein expression. Together, these results establish a mechanism for how codon usage regulates mRNA translation efficiency.

Modulation of HIV-1 Gag/Gag-Pol frameshifting by tRNA abundance

A hallmark of translation in human immunodeficiency virus type 1 (HIV-1) is a –1 programmed ribosome frameshifting event that produces the Gag-Pol fusion polyprotein. The constant Gag to Gag-Pol ratio is essential for the virion structure and infectivity. Here we show that the frameshifting efficiency is modulated by Leu-tRNA^Leu that reads the UUA codon at the mRNA slippery site. This tRNA^Leu isoacceptor is particularly rare in human cell lines derived from T-lymphocytes, the cells that are targeted by HIV-1. When UUA decoding is delayed, the frameshifting follows an alternative route, which maintains the Gag to Gag-Pol ratio constant. A second potential slippery site downstream of the first one is normally inefficient but can also support –1-frameshifting when altered by a compensatory resistance mutation in response to current antiviral drug therapy. Together these different regimes allow the virus to maintain a constant –1-frameshifting efficiency to ensure successful virus propagation.

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.

[Synonymous Codon Usage-a Guide for Co-Translational Protein Folding in the Cell]

In the cell, protein folding begins during protein synthesis/translation and thus is a co-translational process. Co-translational protein folding is tightly linked to translation elongation, which is not a uniform process. While there are many reasons for translation non-uniformity, it is generally believed that non-uniform synonymous codon usage is one of the key factors modulating translation elongation rates. Frequent/optimal codons as a rule are translated more rapidly than infrequently used ones and vice versa. Over 30 years ago, it was hypothesized that changes in synonymous codon usage affecting translation elongation rates could impinge on co-translation protein folding and that many synonymous codons are strategically placed within mRNA to ensure a particular translation kinetics facilitating productive step-by-step co-translational folding of proteins. It was suggested that this particular translation kinetics (and, specifically, translation pause sites) may define the window of opportunity for the protein parts to fold locally, particularly at the critical points where folding is far from equilibrium. It was thus hypothesized that synonymous codons may provide a secondary code for protein folding in the cell. Although, mostly accepted now, this hypothesis appeared to be difficult to prove and many convincing results were obtained only relatively recently. Here, I review the progress in the field and explain, why this simple idea appeared to be so challenging to prove.

Effects of codon optimization on coagulation factor IX translation and structure: Implications for protein and gene therapies

Synonymous codons occur with different frequencies in different organisms, a phenomenon termed codon usage bias. Codon optimization, a common term for a variety of approaches used widely by the biopharmaceutical industry, involves synonymous substitutions to increase protein expression. It had long been presumed that synonymous variants, which, by definition, do not alter the primary amino acid sequence, have no effect on protein structure and function. However, a critical mass of reports suggests that synonymous codon variations may impact protein conformation. To investigate the impact of synonymous codons usage on protein expression and function, we designed an optimized coagulation factor IX (FIX) variant and used multiple methods to compare its properties to the wild-type FIX upon expression in HEK293T cells. We found that the two variants differ in their conformation, even when controlling for the difference in expression levels. Using ribosome profiling, we identified robust changes in the translational kinetics of the two variants and were able to identify a region in the gene that may have a role in altering the conformation of the protein. Our data have direct implications for codon optimization strategies, for production of recombinant proteins and gene therapies.

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.

Conservation of location of several specific inhibitory codon pairs in the Saccharomyces sensu stricto yeasts reveals translational selection

Synonymous codons provide redundancy in the genetic code that influences translation rates in many organisms, in which overall codon use is driven by selection for optimal codons. It is unresolved if or to what extent translational selection drives use of suboptimal codons or codon pairs. In Saccharomyces cerevisiae, 17 specific inhibitory codon pairs, each comprised of adjacent suboptimal codons, inhibit translation efficiency in a manner distinct from their constituent codons, and many are translated slowly in native genes. We show here that selection operates within Saccharomyces sensu stricto yeasts to conserve nine of these codon pairs at defined positions in genes. Conservation of these inhibitory codon pairs is significantly greater than expected, relative to conservation of their constituent codons, with seven pairs more highly conserved than any other synonymous pair. Conservation is strongly correlated with slow translation of the pairs. Conservation of suboptimal codon pairs extends to two related Candida species, fungi that diverged from Saccharomyces ∼270 million years ago, with an enrichment for codons decoded by I•A and U•G wobble in both Candida and Saccharomyces. Thus, conservation of inhibitory codon pairs strongly implies selection for slow translation at particular gene locations, executed by suboptimal codon pairs.

Structural Basis of the Subcellular Topology Landscape of Escherichia coli

Cellular proteomes are distributed in multiple compartments: on DNA, ribosomes, on and inside membranes, or they become secreted. Structural properties that allow polypeptides to occupy subcellular niches, particularly to after crossing membranes, remain unclear. We compared intrinsic and extrinsic features in cytoplasmic and secreted polypeptides of the Escherichia coli K-12 proteome. Structural features between the cytoplasmome and secretome are sharply distinct, such that a signal peptide-agnostic machine learning tool distinguishes cytoplasmic from secreted proteins with 95.5% success. Cytoplasmic polypeptides are enriched in aliphatic, aromatic, charged and hydrophobic residues, unique folds and higher early folding propensities. Secretory polypeptides are enriched in polar/small amino acids, β folds, have higher backbone dynamics, higher disorder and contact order and are more often intrinsically disordered. These non-random distributions and experimental evidence imply that evolutionary pressure selected enhanced secretome flexibility, slow folding and looser structures, placing the secretome in a distinct protein class. These adaptations protect the secretome from premature folding during its cytoplasmic transit, optimize its lipid bilayer crossing and allowed it to acquire cell envelope specific chemistries. The latter may favor promiscuous multi-ligand binding, sensing of stress and cell envelope structure changes. In conclusion, enhanced flexibility, slow folding, looser structures and unique folds differentiate the secretome from the cytoplasmome. These findings have wide implications on the structural diversity and evolution of modern proteomes and the protein folding problem.

Assessing cell-specific effects of genetic variations using tRNA microarrays

Background

By definition, effect of synonymous single-nucleotide variants (SNVs) on protein folding and function are neutral, as they alter the codon and not the encoded amino acid. Recent examples indicate tissue-specific and transfer RNA (tRNA)-dependent effects of some genetic variations arguing against neutrality of synonymous SNVs for protein biogenesis.

Results

We performed systematic analysis of tRNA abunandance across in various models used in cystic fibrosis (CF) research and drug development, including Fischer rat thyroid (FRT) cells, patient-derived primary human bronchial epithelia (HBE) from lung biopsies, primary human nasal epithelia (HNE) from nasal curettage, intestinal organoids, and airway progenitor-directed differentiation of human induced pluripotent stem cells (iPSCs). These were compared to an immortalized CF bronchial cell model (CFBE41o⁻) and two widely used laboratory cell lines, HeLa and HEK293. We discovered that specific synonymous SNVs exhibited differential effects which correlated with variable concentrations of cognate tRNAs.

Conclusions

Our results highlight ways in which the presence of synonymous SNVs may alter local kinetics of mRNA translation; and thus, impact protein biogenesis and function. This effect is likely to influence results from mechansistic analysis and/or drug screeining efforts, and establishes importance of cereful model system selection based on genetic variation profile.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5864-1) contains supplementary material, which is available to authorized users.

Structural insight into co-translational membrane protein folding

Membrane protein folding studies lag behind those of water-soluble proteins due to immense difficulties of experimental study, resulting from the need to provide a hydrophobic lipid-bilayer environment when investigated in vitro. A sound understanding of folding mechanisms is important for membrane proteins as they contribute to a third of the proteome and are frequently associated with disease when mutated and/or misfolded. Membrane proteins largely consist of α-helical, hydrophobic transmembrane domains, which insert into the membrane, often using the SecYEG/Sec61 translocase system. This mini-review highlights recent advances in techniques that can further our understanding of co-translational folding and notably, the structure and insertion of nascent chains as they emerge from translating ribosomes. This article is part of a Special Issue entitled: Molecular biophysics of membranes and membrane proteins.

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.

Variation and selection on codon usage bias across an entire subphylum

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

A pan-cancer analysis of synonymous mutations

Synonymous mutations have been viewed as silent mutations, since they only affect the DNA and mRNA, but not the amino acid sequence of the resulting protein. Nonetheless, recent studies suggest their significant impact on splicing, RNA stability, RNA folding, translation or co-translational protein folding. Hence, we compile 659194 synonymous mutations found in human cancer and characterize their properties. We provide the user-friendly, comprehensive resource for synonymous mutations in cancer, SynMICdb (http://SynMICdb.dkfz.de), which also contains orthogonal information about gene annotation, recurrence, mutation loads, cancer association, conservation, alternative events, impact on mRNA structure and a SynMICdb score. Notably, synonymous and missense mutations are depleted at the 5'-end of the coding sequence as well as at the ends of internal exons independent of mutational signatures. For patient-derived synonymous mutations in the oncogene KRAS, we indicate that single point mutations can have a relevant impact on expression as well as on mRNA secondary structure.

Synonymous mutations do not alter amino acid sequence but may exert oncogenic effects in other ways. Here, the authors present a catalogue of synonymous mutations in cancer and characterise their properties.

Codon and Codon-Pair Usage Tables (CoCoPUTs): Facilitating Genetic Variation Analyses and Recombinant Gene Design

Usage of sequential codon-pairs is non-random and unique to each species. Codon-pair bias is related to but clearly distinct from individual codon usage bias. Codon-pair bias is thought to affect translational fidelity and efficiency and is presumed to be under the selective pressure. It was suggested that changes in codon-pair utilization may affect human disease more significantly than changes in single codons. Although recombinant gene technologies often take codon-pair usage bias into account, codon-pair usage data/tables are not readily available, thus potentially impeding research efforts. The present computational resource (https://hive.biochemistry.gwu.edu/review/codon2) systematically addresses this issue. Building on our recent HIVE-Codon Usage Tables, we constructed a new database to include genomic codon-pair and dinucleotide statistics of all organisms with sequenced genome, available in the GenBank. We believe that the growing understanding of the importance of codon-pair usage will make this resource an invaluable tool to many researchers in academia and pharmaceutical industry.

"CodonWizard" - An intuitive software tool with graphical user interface for customizable codon optimization in protein expression efforts

Optimization of coding sequences to maximize protein expression yield is often outsourced to external service providers during commercial gene synthesis and thus unfortunately remains a black box for many researchers. The presented software program "CodonWizard" offers scientists a powerful but easy-to-use tool for customizable codon optimization: The intuitive graphical user interface empowers even scientists inexperienced in the art to straightforward design, modify, test and save complex codon optimization strategies and to publicly share successful otimization strategies among the scientific community. "Codon Wizard" provides highly flexible features for sequence analysis and completely customizable modification/optimization of codon usage of any given input sequence data (DNA/RNA/peptide) using freely combinable algorithms, allowing for implementation of contemporary, well-established optimization strategies as well as novel, proprietary ones alike. Contrary to comparable tools, "Codon Wizard" thus finally opens up ways for an empirical approach to codon optimization and may also >be used completely offline to protect resulting intellectual property. As a benchmark, the reliability, intuitiveness and utility of the application could be demonstrated by increasing the yield of recombinant TEV-protease expressed in E. coli by several orders of magnitude after codon optimization using "CodonWizard" - Permanently available for download on the web at http://schwalbe.org.chemie.uni-frankfurt.de/node/3324.

A New Look at Codon Usage and Protein Expression

%MinMax, a model of intra-gene translational elongation rate, relies on codon usage frequencies. Historically, %MinMax has used tables that measure codon usage bias for all genes in an organism, such as those found at HIVE-CUT. In this paper, we provide evidence that codon usage bias based on all genes is insufficient to accurately measure absolute translation rate. We show that alternative "High-ϕ" codon usage tables, generated by another model (ROC-SEMPPR), are a promising alternative. By creating a hybrid model, future codon usage analyses and their applications (e.g., codon harmonization) are likely to more accurately measure the "tempo" of translation elongation. We also suggest a High-ϕ alternative to the Codon Adaptation Index (CAI), a classic metric of codon usage bias based on highly expressed genes. Significantly, our new alternative is equally well correlated with empirical data as traditional CAI without using experimentally determined expression counts as input.

The Ribosome Cooperates with a Chaperone to Guide Multi-domain Protein Folding

Multi-domain proteins, containing several structural units within a single polypeptide, constitute a large fraction of all proteomes. Co-translational folding is assumed to simplify the conformational search problem for large proteins, but the events leading to correctly folded, functional structures remain poorly characterized. Similarly, how the ribosome and molecular chaperones promote efficient folding remains obscure. Using optical tweezers, we have dissected early folding events of nascent elongation factor G, a multi-domain protein that requires chaperones for folding. The ribosome and the chaperone trigger factor reduce inter-domain misfolding, permitting folding of the N-terminal G-domain. Successful completion of this step is a crucial prerequisite for folding of the next domain. Unexpectedly, co-translational folding does not proceed unidirectionally; emerging unfolded polypeptide can denature an already-folded domain. Trigger factor, but not the ribosome, protects against denaturation. The chaperone thus serves a previously unappreciated function, helping multi-domain proteins overcome inherent challenges during co-translational folding.

Unravelling Ribosome Function Through Structural Studies

Ribosomes are biological nanomachine that synthesise all proteins within a cell. It took decades to reveal the architecture of this essential cellular component. To understand the structure -function relationship of this nanomachine needed the utilisisation of different biochemical, biophysical and structural techniques. Structural studies combined with mutagenesis of the different ribosomal complexes comprising various RNAs and proteins enabled us to understand how this machine works inside a cell. Nowadays quite a number of ribosomal structures were published that confirmed biochemical studies on particular steps of protein synthesis by the ribosome . Four major steps were identified: initiation , elongation, termination and recycling. These steps lead us to the important question how the ribosome function can be regulated. Advances in technology for cryo electron microscopy: sample preparations, image recording, developments in algorithms for image analysis and processing significantly helped in revelation of structural details of the ribosome . We now have a library of ribosome structures from prokaryotes to eukaryotes that enable us to understand the complex mechanics of this nanomachine. As this structural library continues to grow, we gradually improve our understanding of this process and how it can be regulated and how the specific ribosomes can be stalled or activated, or completely disabled. This article provides a comprehensive overview of ribosomal structures that represent structural snapshots of the ribosome at its different functional states. Better understanding rises more particular questions that have to be addressed by determination structures of more complexes.Synopsis: Structural biology of the ribosome.

The new strategies to overcome challenges in protein production in bacteria

Recombinant proteins are essential for biotechnology. Here we review some of the key points for improving the production of heterologous proteins, and what can be the future of the field.

The stop-and-go traffic regulating protein biogenesis: How translation kinetics controls proteostasis

Generating a functional proteome requires the ribosome to carefully regulate disparate co-translational processes that determine the fate of nascent polypeptides. With protein synthesis being energetically expensive, the ribosome must balance the costs of efficiently making a protein with those of properly folding it. Emerging as a primary means of regulating this trade-off is the nonuniform rate of translation elongation that defines translation kinetics. The varying speeds with which the ribosome progresses along a transcript have been implicated in several aspects of protein biogenesis, including co-translational protein folding and translational fidelity, as well as gene expression by mediating mRNA decay and protein quality control pathways. The optimal translation kinetics required to efficiently execute these processes can be distinct. Thus, the ribosome is tasked with tightly regulating translation kinetics to balance these processes while maintaining adaptability for changing cellular conditions. In this review, we first discuss the regulatory role of translation elongation in protein biogenesis and what factors influence elongation kinetics. We then describe how changes in translation kinetics signal downstream pathways that dictate the fate of nascent polypeptides. By regulating these pathways, the kinetics of translation elongation has emerged as a critical tool for driving gene expression and maintaining proteostasis through varied mechanisms, including nascent chain folding and binding different ribosome-associated machinery. Indeed, a growing number of examples demonstrate the important role of local changes in elongation kinetics in modulating the pathophysiology of human disease.

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.

Folding up and Moving on-Nascent Protein Folding on the Ribosome

All cellular proteins are synthesized by the ribosome, an intricate molecular machine that translates the information of protein coding genes into the amino acid alphabet. The linear polypeptides synthesized by the ribosome must generally fold into specific three-dimensional structures to become biologically active. Folding has long been recognized to begin before synthesis is complete. Recently, biochemical and biophysical studies have shed light onto how the ribosome shapes the folding pathways of nascent proteins. Here, we discuss recent progress that is beginning to define the role of the ribosome in the folding of newly synthesized polypeptides.

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3 kcal mol⁻¹. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.

•

Adjacent N and C terminal strands are most structured region in the transition state.

•

Two circular permutants retain the same folding pathway as wild-type MATH.

•

On the ribosome, early emergence of terminal strands does not promote earlier folding.

•

Formation of both β-sheets is energetically critical for folding on the ribosome.

•

Folding pathway minimizes formation of partly structured states prone to mis-folding.

Codon usage regulates human KRAS expression at both transcriptional and translational levels

KRAS and HRAS are highly homologous oncogenic Ras GTPase family members that are mutated in a wide spectrum of human cancers. Despite having high amino acid identity, KRAS and HRAS have very different codon usage biases: the HRAS gene contains many common codons, and KRAS is enriched for rare codons. Rare codons in KRAS suppress its protein expression, which has been shown to affect both normal and cancer biology in mammals. Here, using HRAS or KRAS expression in different human cell lines and in vitro transcription and translation assays, we show that KRAS rare codons inhibit both translation efficiency and transcription and that the contribution of these two processes varies among different cell lines. We observed that codon usage regulates mRNA translation efficiency such that WT KRAS mRNA is poorly translated. On the other hand, common codons increased transcriptional rates by promoting activating histone modifications and recruitment of transcriptional coactivators. Finally, we found that codon usage also influences KRAS protein conformation, likely because of its effect on co-translational protein folding. Together, our results reveal that codon usage has multidimensional effects on protein expression, ranging from effects on transcription to protein folding in human cells.

Peripheral Protein Quality Control as a Novel Drug Target for CFTR Stabilizer

Conformationally defective cystic fibrosis transmembrane conductance regulator (CFTR) including rescued ΔF508-CFTR is rapidly eliminated from the plasma membrane (PM) even in the presence of a CFTR corrector and potentiator, limiting the therapeutic effort of the combination therapy. CFTR elimination from the PM is determined by the conformation-dependent ubiquitination as a part of the peripheral quality control (PQC) mechanism. Recently, the molecular machineries responsible for the CFTR PQC mechanism which includes molecular chaperones and ubiquitination enzymes have been revealed. This review summarizes the molecular mechanism of the CFTR PQC and discusses the possibility that the peripheral ubiquitination mechanism becomes a novel drug target to develop the CFTR stabilizer as a novel class of CFTR modulator.

Translation in Prokaryotes

This review summarizes our current understanding of translation in prokaryotes, focusing on the mechanistic and structural aspects of each phase of translation: initiation, elongation, termination, and ribosome recycling. The assembly of the initiation complex provides multiple checkpoints for messenger RNA (mRNA) and start-site selection. Correct codon-anticodon interaction during the decoding phase of elongation results in major conformational changes of the small ribosomal subunit and shapes the reaction pathway of guanosine triphosphate (GTP) hydrolysis. The ribosome orchestrates proton transfer during peptide bond formation, but requires the help of elongation factor P (EF-P) when two or more consecutive Pro residues are to be incorporated. Understanding the choreography of transfer RNA (tRNA) and mRNA movements during translocation helps to place the available structures of translocation intermediates onto the time axis of the reaction pathway. The nascent protein begins to fold cotranslationally, in the constrained space of the polypeptide exit tunnel of the ribosome. When a stop codon is reached at the end of the coding sequence, the ribosome, assisted by termination factors, hydrolyzes the ester bond of the peptidyl-tRNA, thereby releasing the nascent protein. Following termination, the ribosome is dissociated into subunits and recycled into another round of initiation. At each step of translation, the ribosome undergoes dynamic fluctuations between different conformation states. The aim of this article is to show the link between ribosome structure, dynamics, and function.

Synonymous mutations make dramatic contributions to fitness when growth is limited by a weak-link enzyme

Synonymous mutations do not alter the specified amino acid but may alter the structure or function of an mRNA in ways that impact fitness. There are few examples in the literature, however, in which the effects of synonymous mutations on microbial growth rates have been measured, and even fewer for which the underlying mechanism is understood. We evolved four populations of a strain of Salmonella enterica in which a promiscuous enzyme has been recruited to replace an essential enzyme. A previously identified point mutation increases the enzyme’s ability to catalyze the newly needed reaction (required for arginine biosynthesis) but decreases its ability to catalyze its native reaction (required for proline biosynthesis). The poor performance of this enzyme limits growth rate on glucose. After 260 generations, we identified two synonymous mutations in the first six codons of the gene encoding the weak-link enzyme that increase growth rate by 41 and 67%. We introduced all possible synonymous mutations into the first six codons and found substantial effects on growth rate; one doubles growth rate, and another completely abolishes growth. Computational analyses suggest that these mutations affect either the stability of a stem-loop structure that sequesters the start codon or the accessibility of the region between the Shine-Dalgarno sequence and the start codon. Thus, these mutations would be predicted to affect translational efficiency and thereby indirectly affect mRNA stability because translating ribosomes protect mRNA from degradation. Experimental data support these hypotheses. We conclude that the effects of the synonymous mutations are due to a combination of effects on mRNA stability and translation efficiency that alter levels of the weak-link enzyme. These findings suggest that synonymous mutations can have profound effects on fitness under strong selection and that their importance in evolution may be under-appreciated.

When a new enzyme is needed, microbes often recruit a pre-existing enzyme with a promiscuous activity corresponding to the newly needed activity. Such enzymes are often the “weak-link” in metabolism because they have not evolved to efficiently catalyze the new reaction. Under these circumstances, increasing the level of the weak-link enzyme can improve fitness. We evolved a strain of S. enterica in which a weak-link enzyme–E383A ProA–serves essential functions in synthesis of proline and arginine for 260 generations and then sequenced the genomes of several evolved strains. A mutation in the promoter of the operon encoding E383A ProA increased growth rate 9-fold. More surprisingly, a mutation upstream of the start codon and two synonymous mutations within the first six codons also increased growth rate by up to 68%. Introduction of all possible synonymous mutations in the first six codons showed that some doubled growth rate, while others slowed or even prevented growth. Computational and experimental data suggest that these effects were due to enhanced translational efficiency of the weak-link enzyme. These results show that synonymous mutations, once assumed to be selectively neutral, can have strong impacts on fitness when growth rate is limited by a weak-link enzyme.

Genetic code expansion via integration of redundant amino acid assignment by finely tuning tRNA pools

In all translation systems, the genetic code assigns codons to amino acids as building blocks of polypeptides, defining their chemical, structural and physiological properties. The canonical genetic code, however, utilizes only 20 proteinogenic amino acids redundantly encoded in 61 codons. In order to expand the building block repertoire, this redundancy was reduced by tuning composition of the transfer RNA (tRNA) mixture in vitro. Depletion of particular tRNAs from the total tRNA mixture or its reconstitution with in vitro-transcribed tRNA_SNNs (S = C or G, N = U, C, A or G) divided a codon box to encode two amino acids, expanding the repertoire to 23. The expanded genetic codes may benefit analysis of cellular regulatory pathways and drug screening.