Synonymous Codons Direct Cotranslational Folding toward Different Protein Conformations

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Folding of heterologous proteins in bacterial cell factories: Cellular mechanisms and engineering strategies

The expression of correctly folded and functional heterologous proteins is important in many biotechnological production processes, whether it is enzymes, biopharmaceuticals or biosynthetic pathways for production of sustainable chemicals. For industrial applications, bacterial platform organisms, such as E. coli, are still broadly used due to the availability of tools and proven suitability at industrial scale. However, expression of heterologous proteins in these organisms can result in protein aggregation and low amounts of functional protein. This review provides an overview of the cellular mechanisms that can influence protein folding and expression, such as co-translational folding and assembly, chaperone binding, as well as protein quality control, across different model organisms. The knowledge of these mechanisms is then linked to different experimental methods that have been applied in order to improve functional heterologous protein folding, such as codon optimization, fusion tagging, chaperone co-production, as well as strain and protein engineering strategies.

How synonymous mutations alter enzyme structure and function over long timescales

The specific activity of enzymes can be altered over long timescales in cells by synonymous mutations that alter a messenger RNA molecule's sequence but not the encoded protein's primary structure. How this happens at the molecular level is unknown. Here, we use multiscale modelling of three Escherichia coli enzymes (type III chloramphenicol acetyltransferase, D-alanine-D-alanine ligase B and dihydrofolate reductase) to understand experimentally measured changes in specific activity due to synonymous mutations. The modelling involves coarse-grained simulations of protein synthesis and post-translational behaviour, all-atom simulations to test robustness and quantum mechanics/molecular mechanics calculations to characterize enzymatic function. We show that changes in codon translation rates induced by synonymous mutations cause shifts in co-translational and post-translational folding pathways that kinetically partition molecules into subpopulations that very slowly interconvert to the native, functional state. Structurally, these states resemble the native state, with localized misfolding near the active sites of the enzymes. These long-lived states exhibit reduced catalytic activity, as shown by their increased activation energies for the reactions they catalyse.

Reading Frame Retrieval of Genes: A New Parameter of Codon Usage Based on the Circular Code Theory

Based on the circular code theory, we define a new function f that quantifies the property of reading frame retrieval (RFR) of genes from their codon usage. This RFR function f is computed on a massive scale in genes of genomes of bacteria, eukaryotes and archaea. By expressing f as a function of the mean number [Formula: see text] of codons per gene, a "universal" property is identified, whatever the kingdom: the reading frame retrieval is enhanced in large genes. By investigating this property according to the theory developed, a Spearman's rank correlation with a strong negative coefficient is observed between the codon usage dispersion d (from the uniform codon distribution [Formula: see text]) and the RFR function f, whatever the kingdom (p-values [Formula: see text] in bacteria, [Formula: see text] in eukaryotes and [Formula: see text] in archaea). Thus, the reading frame retrieval is enhanced with the codon usage dispersion. Furthermore, this approach identifies a "genome centre" from which emerge two distinct "genome arms": an upper arm and a lower arm, respectively, above and below the linear regression. The RFR function by itself or combined with classical methods (alignment, phylogeny) could also be a new approach to classify the genomes in the future.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies

BACKGROUND

Since the onset of the SARS-CoV-2 pandemic, bioinformatic analyses have been performed to understand the nucleotide and synonymous codon usage features and mutational patterns of the virus. However, comparatively few have attempted to perform such analyses on a considerably large cohort of viral genomes while organizing the plethora of available sequence data for a month-by-month analysis to observe changes over time. Here, we aimed to perform sequence composition and mutation analysis of SARS-CoV-2, separating sequences by gene, clade, and timepoints, and contrast the mutational profile of SARS-CoV-2 to other comparable RNA viruses.

METHODS

Using a cleaned, filtered, and pre-aligned dataset of over 3.5 million sequences downloaded from the GISAID database, we computed nucleotide and codon usage statistics, including calculation of relative synonymous codon usage values. We then calculated codon adaptation index (CAI) changes and a nonsynonymous/synonymous mutation ratio (dN/dS) over time for our dataset. Finally, we compiled information on the types of mutations occurring for SARS-CoV-2 and other comparable RNA viruses, and generated heatmaps showing codon and nucleotide composition at high entropy positions along the Spike sequence.

RESULTS

We show that nucleotide and codon usage metrics remain relatively consistent over the 32-month span, though there are significant differences between clades within each gene at various timepoints. CAI and dN/dS values vary substantially between different timepoints and different genes, with Spike gene on average showing both the highest CAI and dN/dS values. Mutational analysis showed that SARS-CoV-2 Spike has a higher proportion of nonsynonymous mutations than analogous genes in other RNA viruses, with nonsynonymous mutations outnumbering synonymous ones by up to 20:1. However, at several specific positions, synonymous mutations were overwhelmingly predominant.

CONCLUSIONS

Our multifaceted analysis covering both the composition and mutation signature of SARS-CoV-2 gives valuable insight into the nucleotide frequency and codon usage heterogeneity of SARS-CoV-2 over time, and its unique mutational profile compared to other RNA viruses.

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.

The roles of mutated SPINK1 gene in prostate cancer cells

SPINK1-positive prostate cancer (PCa) has been identified as an aggressive PCa subtype. However, there is a lack of definite studies to elucidate the underlying mechanism of the loss of SPINK1 expression in most PCa cells except 22Rv1 cells, which are derived from a human prostatic carcinoma xenograft, CWR22R. The aim of this study was to investigate the mechanisms of SPINK1 protein positive/negative expression and its biological roles in PCa cell lines. SPINK1 mRNA was highly expressed in 22Rv1 cells compared with LNCaP, C4-2B, DU145, and PC-3 cells, and the protein was only detected in 22Rv1 cells. Among these cell lines, the wild-type SPINK1 coding sequence was only found in 22Rv1 cells, and two mutation sites, the c.194G>A missense mutation and the c.210T>C synonymous mutation, were found in other cell lines. Our further research showed that the mutations were associated with a reduction in SPINK1 mRNA and protein levels. Functional experiments indicated that SPINK1 promoted PC-3 cell proliferation, migration, and invasion, while knockdown of SPINK1 attenuated 22Rv1 cell proliferation, migration, and invasion. The wild-type SPINK1 gene can promote the malignant behaviors of cells more than the mutated ones. Cell cycle analysis by flow cytometry showed that SPINK1 decreased the percentage of cells in the G0/G1 phase and increased the percentage of S phase cells. We demonstrated that the c.194G>A and c.210T>C mutations in the SPINK1 gene decreased the mRNA and protein levels. The wild-type SPINK1 gene is related to aggressive biological behaviors of PCa cells and may be a potential therapeutic target for PCa.

A male germ-cell-specific ribosome controls male fertility

Ribosomes are highly sophisticated translation machines that have been demonstrated to be heterogeneous in the regulation of protein synthesis^1,2. Male germ cell development involves complex translational regulation during sperm formation³. However, it remains unclear whether translation during sperm formation is performed by a specific ribosome. Here we report a ribosome with a specialized nascent polypeptide exit tunnel, Ribosome^ST, that is assembled with the male germ-cell-specific protein RPL39L, the paralogue of core ribosome (Ribosome^Core) protein RPL39. Deletion of Ribosome^ST in mice causes defective sperm formation, resulting in substantially reduced fertility. Our comparison of single-particle cryo-electron microscopy structures of ribosomes from mouse kidneys and testes indicates that Ribosome^ST features a ribosomal polypeptide exit tunnel of distinct size and charge states compared with Ribosome^Core. Ribosome^ST predominantly cotranslationally regulates the folding of a subset of male germ-cell-specific proteins that are essential for the formation of sperm. Moreover, we found that specialized functions of Ribosome^ST were not replaceable by Ribosome^Core. Taken together, identification of this sperm-specific ribosome should greatly expand our understanding of ribosome function and tissue-specific regulation of protein expression pattern in mammals.

Bacterial Hsp90 predominantly buffers but does not potentiate the phenotypic effects of deleterious mutations during fluorescent protein evolution

Chaperones facilitate the folding of other ("client") proteins and can thus affect the adaptive evolution of these clients. Specifically, chaperones affect the phenotype of proteins via two opposing mechanisms. On the one hand, they can buffer the effects of mutations in proteins and thus help preserve an ancestral, premutation phenotype. On the other hand, they can potentiate the effects of mutations and thus enhance the phenotypic changes caused by a mutation. We study that how the bacterial Hsp90 chaperone (HtpG) affects the evolution of green fluorescent protein. To this end, we performed directed evolution of green fluorescent protein under low and high cellular concentrations of Hsp90. Specifically, we evolved green fluorescent protein under both stabilizing selection for its ancestral (green) phenotype and directional selection toward a new (cyan) phenotype. While Hsp90 did only affect the rate of adaptive evolution transiently, it did affect the phenotypic effects of mutations that occurred during adaptive evolution. Specifically, Hsp90 allowed strongly deleterious mutations to accumulate in evolving populations by buffering their effects. Our observations show that the role of a chaperone for adaptive evolution depends on the organism and the trait being studied.

Oxidative stress strongly restricts the effect of codon choice on the efficiency of protein synthesis in Escherichia coli

INTRODUCTION

The response of enterobacteria to oxidative stress is usually considered to be regulated by transcription factors such as OxyR and SoxR. Nevertheless, several reports have shown that under oxidative stress the levels, modification and aminoacylation of tRNAs may be altered suggesting a role of codon bias in regulation of gene expression under this condition.

METHODS

In order to characterize the effects of oxidative stress on translation elongation we constructed a library of 61 plasmids, each coding for the green fluorescent protein (GFP) translationally fused to a different set of four identical codons.

RESULTS

Using these reporters, we observed that GFP production levels vary widely (~15 fold) when Escherichia coli K-12 is cultured in minimal media as a consequence of codon choice variations. When bacteria are cultured under oxidative stress caused by paraquat the levels of GFP produced by most clones is reduced and, in contrast to control conditions, the range of GFP levels is restricted to a ~2 fold range. Restricting elongation of particular sequences does not increase the range of GFP production under oxidative stress, but altering translation initiation rates leads to an increase in this range.

DISCUSSION

Altogether, our results suggest that under normal conditions the speed of translation elongation is in the range of the speed of initiation and, consequently, codon choice impacts the speed of protein synthesis. In contrast, under oxidative stress translation initiation becomes much slower than elongation, limiting the speed of translation such that codon choice has at most only subtle effects on the overall output of translation.

The deubiquitinase OTUD1 regulates immunoglobulin production and proteasome inhibitor sensitivity in multiple myeloma

Serum monoclonal immunoglobulin (Ig) is the main diagnostic factor for patients with multiple myeloma (MM), however its prognostic potential remains unclear. On a large MM patient cohort (n = 4146), we observe no correlation between serum Ig levels and patient survival, while amount of intracellular Ig has a strong predictive effect. Focused CRISPR screen, transcriptional and proteomic analysis identify deubiquitinase OTUD1 as a critical mediator of Ig synthesis, proteasome inhibitor sensitivity and tumor burden in MM. Mechanistically, OTUD1 deubiquitinates peroxiredoxin 4 (PRDX4), protecting it from endoplasmic reticulum (ER)-associated degradation. In turn, PRDX4 facilitates Ig production which coincides with the accumulation of unfolded proteins and higher ER stress. The elevated load on proteasome ultimately potentiates myeloma response to proteasome inhibitors providing a window for a rational therapy. Collectively, our findings support the significance of the Ig production machinery as a biomarker and target in the combinatory treatment of MM patients.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi

Zinc binuclear cluster proteins (ZBC; Zn(II)2Cys6) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)2Cys6 coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)2Cys6 genes. The perseverance of GC-rich codons of Zn(II)2Cys6 in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

The High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

Biosynthesis of plant hemostatic dencichine in Escherichia coli

Dencichine is a plant-derived nature product that has found various pharmacological applications. Currently, its natural biosynthetic pathway is still elusive, posing challenge to its heterologous biosynthesis. In this work, we design artificial pathways through retro-biosynthesis approaches and achieve de novo production of dencichine. First, biosynthesis of the two direct precursors L-2, 3-diaminopropionate and oxalyl-CoA is achieved by screening and integrating microbial enzymes. Second, the solubility of dencichine synthase, which is the last and only plant-derived pathway enzyme, is significantly improved by introducing 28 synonymous rare codons into the codon-optimized gene to slow down its translation rate. Last, the metabolic network is systematically engineered to direct the carbon flux to dencichine production, and the final titer reaches 1.29 g L-1 with a yield of 0.28 g g-1 glycerol. This work lays the foundation for sustainable production of dencichine and represents an example of how synthetic biology can be harnessed to generate unnatural pathways to produce a desired molecule.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Lipid nanoparticle-based mRNA vaccines in cancers: Current advances and future prospects

Messenger RNA (mRNA) vaccines constitute an emerging therapeutic method with the advantages of high safety and efficiency as well as easy synthesis; thus, they have been widely used in various human diseases, especially in malignant cancers. However, the mRNA vaccine technology has some limitations, such as instability and low transitive efficiency in vivo, which greatly restrict its application. The development of nanotechnology in the biomedical field offers new strategies and prospects for the early diagnosis and treatment of human cancers. Recent studies have demonstrated that Lipid nanoparticle (LNP)-based mRNA vaccines can address the poor preservation and targeted inaccuracy of mRNA vaccines. As an emerging cancer therapy, mRNA vaccines potentially have broad future applications. Unlike other treatments, cancer mRNA vaccines provide specific, safe, and tolerable treatments. Preclinical studies have used personalized vaccines to demonstrate the anti-tumor effect of mRNA vaccines in the treatment of various solid tumors, including colorectal and lung cancer, using these in a new era of therapeutic cancer vaccines. In this review, we have summarized the latest applications and progress of LNP-based mRNA vaccines in cancers, and discussed the prospects and limitations of these fields, thereby providing novel strategies for the targeted therapy of cancers.

Oligodendrocyte differentiation alters tRNA modifications and codon optimality-mediated mRNA decay

Oligodendrocytes are specialized cells that confer neuronal myelination in the central nervous system. Leukodystrophies associated with oligodendrocyte deficits and hypomyelination are known to result when a number of tRNA metabolism genes are mutated. Thus, for unknown reasons, oligodendrocytes may be hypersensitive to perturbations in tRNA biology. In this study, we survey the tRNA transcriptome in the murine oligodendrocyte cell lineage and find that specific tRNAs are hypomodified in oligodendrocytes within or near the anticodon compared to oligodendrocyte progenitor cells (OPCs). This hypomodified state may be the result of differential expression of key modification enzymes during oligodendrocyte differentiation. Moreover, we observe a concomitant relationship between tRNA hypomodification and tRNA decoding potential; observing oligodendrocyte specific alterations in codon optimality-mediated mRNA decay and ribosome transit. Our results reveal that oligodendrocytes naturally maintain a delicate, hypersensitized tRNA/mRNA axis. We suggest this axis is a potential mediator of pathology in leukodystrophies and white matter disease when further insult to tRNA metabolism is introduced.

When a Synonymous Variant Is Nonsynonymous

Term synonymous variation is widely used, but frequently in a wrong or misleading meaning and context. Twenty three point eight % of possible nucleotide substitution types in the universal genetic code are for synonymous amino acid changes, but when these variants have a phenotype and functional effect, they are very seldom synonymous. Such variants may manifest changes at DNA, RNA and/or protein levels. Large numbers of variations are erroneously annotated as synonymous, which causes problems e.g., in clinical genetics and diagnosis of diseases. To facilitate precise communication, novel systematics and nomenclature are introduced for variants that when looking only at the genetic code seem like synonymous, but which have phenotypes. A new term, unsense variant is defined as a substitution in the mRNA coding region that affects gene expression and protein production without introducing a stop codon in the variation site. Such variants are common and need to be correctly annotated. Proper naming and annotation are important also to increase awareness of these variants and their consequences.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Ribosome selectivity and nascent chain context in modulating the incorporation of fluorescent non-canonical amino acid into proteins

Fluorescence reporter groups are important tools to study the structure and dynamics of proteins. Genetic code reprogramming allows for cotranslational incorporation of non-canonical amino acids at any desired position. However, cotranslational incorporation of bulky fluorescence reporter groups is technically challenging and usually inefficient. Here we analyze the bottlenecks for the cotranslational incorporation of NBD-, BodipyFL- and Atto520-labeled Cys-tRNA^Cys into a model protein using a reconstituted in-vitro translation system. We show that the modified Cys-tRNA^Cys can be rejected during decoding due to the reduced ribosome selectivity for the modified aa-tRNA and the competition with native near-cognate aminoacyl-tRNAs. Accommodation of the modified Cys-tRNA^Cys in the A site of the ribosome is also impaired, but can be rescued by one or several Gly residues at the positions −1 to −4 upstream of the incorporation site. The incorporation yield depends on the steric properties of the downstream residue and decreases with the distance from the protein N-terminus to the incorporation site. In addition to the full-length translation product, we find protein fragments corresponding to the truncated N-terminal peptide and the C-terminal fragment starting with a fluorescence-labeled Cys arising from a StopGo-like event due to a defect in peptide bond formation. The results are important for understanding the reasons for inefficient cotranslational protein labeling with bulky reporter groups and for designing new approaches to improve the yield of fluorescence-labeled protein.

Recent Advances in the Molecular Design and Delivery Technology of mRNA for Vaccination Against Infectious Diseases

Vaccines can prevent many millions of illnesses against infectious diseases and save numerous lives every year. However, traditional vaccines such as inactivated viral and live attenuated vaccines cannot adapt to emerging pandemics due to their time-consuming development. With the global outbreak of the COVID-19 epidemic, the virus continues to evolve and mutate, producing mutants with enhanced transmissibility and virulence; the rapid development of vaccines against such emerging global pandemics becomes more and more critical. In recent years, mRNA vaccines have been of significant interest in combating emerging infectious diseases due to their rapid development and large-scale production advantages. However, their development still suffers from many hurdles such as their safety, cellular delivery, uptake, and response to their manufacturing, logistics, and storage. More efforts are still required to optimize the molecular designs of mRNA molecules with increased protein expression and enhanced structural stability. In addition, a variety of delivery systems are also needed to achieve effective delivery of vaccines. In this review, we highlight the advances in mRNA vaccines against various infectious diseases and discuss the molecular design principles and delivery systems of associated mRNA vaccines. The current state of the clinical application of mRNA vaccine pipelines against various infectious diseases and the challenge, safety, and protective effect of associated vaccines are also discussed.

Structural, functional, and immunogenicity implications of F9 gene recoding

Hemophilia B is a blood clotting disorder caused by deficient activity of coagulation factor IX (FIX). Multiple recombinant FIX proteins are currently approved to treat hemophilia B, and several gene therapy products are currently being developed. Codon optimization is a frequently used technique in the pharmaceutical industry to improve recombinant protein expression by recoding a coding sequence using multiple synonymous codon substitutions. The underlying assumption of this gene recoding is that synonymous substitutions do not alter protein characteristics because the primary sequence of the protein remains unchanged. However, a critical body of evidence shows that synonymous variants can affect cotranslational folding and protein function. Gene recoding could potentially alter the structure, function, and in vivo immunogenicity of recoded therapeutic proteins. Here, we evaluated multiple recoded variants of F9 designed to further explore the effects of codon usage bias on protein properties. The detailed evaluation of these constructs showed altered conformations, and assessment of translation kinetics by ribosome profiling revealed differences in local translation kinetics. Assessment of wild-type and recoded constructs using a major histocompatibility complex (MHC)-associated peptide proteomics assay showed distinct presentation of FIX-derived peptides bound to MHC class II molecules, suggesting that despite identical amino acid sequence, recoded proteins could exhibit different immunogenicity risks. Posttranslational modification analysis indicated that overexpression from gene recoding results in suboptimal posttranslational processing. Overall, our results highlight potential functional and immunogenicity concerns associated with gene-recoded F9 products. These findings have general applicability and implications for other gene-recoded recombinant proteins.

Preferred synonymous codons are translated more accurately: Proteomic evidence, among-species variation, and mechanistic basis

A commonly stated cause of unequal uses of synonymous codons is their differential translational accuracies. A key prediction of this long-standing translational accuracy hypothesis (TAH) of codon usage bias is higher translational accuracies of more frequently used synonymous codons, which, however, has had no direct evidence beyond case studies. Analyzing proteomic data from Escherichia coli, we present direct, global evidence for this prediction. The experimentally measured codon-specific translational accuracies validate a sequence-based proxy; this proxy provides support for the TAH from the vast majority of over 1000 taxa surveyed in all domains of life. We find that the relative translational accuracies of synonymous codons vary substantially among taxa and are strongly correlated with the amounts of cognate transfer RNAs (tRNAs) relative to those of near-cognate tRNAs. These and other observations suggest a model in which selections for translational efficiency and accuracy drive codon usage bias and its coevolution with the tRNA pool.

Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less functional

Some misfolded protein conformations can bypass proteostasis machinery and remain soluble in vivo. This is an unexpected observation, as cellular quality control mechanisms should remove misfolded proteins. Three questions, then, are: how do long-lived, soluble, misfolded proteins bypass proteostasis? How widespread are such misfolded states? And how long do they persist? We address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of proteins exhibit misfolded subpopulations that bypass molecular chaperones, avoid aggregation, and will not be rapidly degraded, with some misfolded states persisting for months or longer. The surface properties of these misfolded states are native-like, suggesting they will remain soluble, while self-entanglements make them long-lived kinetic traps. In terms of function, we predict that one-third of proteins can misfold into soluble less-functional states. For the heavily entangled protein glycerol-3-phosphate dehydrogenase, limited-proteolysis mass spectrometry experiments interrogating misfolded conformations of the protein are consistent with the structural changes predicted by our simulations. These results therefore provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass proteostasis, and indicate, unexpectedly, this may be a wide-spread phenomenon.

Synonymous mutations in representative yeast genes are mostly strongly non-neutral

Synonymous mutations in protein-coding genes do not alter protein sequences and are thus generally presumed to be neutral or nearly neutral^1-5. Here, to experimentally verify this presumption, we constructed 8,341 yeast mutants each carrying a synonymous, nonsynonymous or nonsense mutation in one of 21 endogenous genes with diverse functions and expression levels and measured their fitness relative to the wild type in a rich medium. Three-quarters of synonymous mutations resulted in a significant reduction in fitness, and the distribution of fitness effects was overall similar-albeit nonidentical-between synonymous and nonsynonymous mutations. Both synonymous and nonsynonymous mutations frequently disturbed the level of mRNA expression of the mutated gene, and the extent of the disturbance partially predicted the fitness effect. Investigations in additional environments revealed greater across-environment fitness variations for nonsynonymous mutants than for synonymous mutants despite their similar fitness distributions in each environment, suggesting that a smaller proportion of nonsynonymous mutants than synonymous mutants are always non-deleterious in a changing environment to permit fixation, potentially explaining the common observation of substantially lower nonsynonymous than synonymous substitution rates. The strong non-neutrality of most synonymous mutations, if it holds true for other genes and in other organisms, would require re-examination of numerous biological conclusions about mutation, selection, effective population size, divergence time and disease mechanisms that rely on the assumption that synoymous mutations are neutral.

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Background

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

Results

Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to prevous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding.

Conclusions

Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08635-0).

Codon-specific Ramachandran plots show amino acid backbone conformation depends on identity of the translated codon

Synonymous codons translate into chemically identical amino acids. Once considered inconsequential to the formation of the protein product, there is evidence to suggest that codon usage affects co-translational protein folding and the final structure of the expressed protein. Here we develop a method for computing and comparing codon-specific Ramachandran plots and demonstrate that the backbone dihedral angle distributions of some synonymous codons are distinguishable with statistical significance for some secondary structures. This shows that there exists a dependence between codon identity and backbone torsion of the translated amino acid. Although these findings cannot pinpoint the causal direction of this dependence, we discuss the vast biological implications should coding be shown to directly shape protein conformation and demonstrate the usefulness of this method as a tool for probing associations between codon usage and protein structure. Finally, we urge for the inclusion of exact genetic information into structural databases.

Analysis of codon usage patterns in open reading frame 4 of hepatitis E viruses

BACKGROUND

Hepatitis E virus (HEV) is a member of the family Hepeviridae and causes acute HEV infections resulting in thousands of deaths worldwide. The zoonotic nature of HEV in addition to its tendency from human to human transmission has led scientists across the globe to work on its different aspects. HEV also accounts for about 30% mortality rates in case of pregnant women. The genome of HEV is organized into three open reading frames (ORFs): ORF1 ORF2 and ORF3. A reading frame encoded protein ORF4 has recently been discovered which is exclusive to GT 1 isolates of HEV. The ORF4 is suggested to play crucial role in pregnancy-associated pathology and enhanced replication. Though studies have documented the ORF4's importance, the genetic features of ORF4 protein genes in terms of compositional patterns have not been elucidated. As codon usage performs critical role in establishment of the host-pathogen relationship, therefore, the present study reports the codon usage analysis (based on nucleotide sequences of HEV ORF4 available in the public database) in three hosts along with the factors influencing the codon usage patterns of the protein genes of ORF4 of HEV.

RESULTS

The nucleotide composition analysis indicated that ORF4 protein genes showed overrepresentation of C nucleotide and while A nucleotide was the least-represented, with random distribution of G and T(U) nucleotides. The relative synonymous codon usage (RSCU) analysis revealed biasness toward C/G-ended codons (over U/A) in all three natural HEV-hosts (human, rat and ferret). It was observed that all the ORF4 genes were richly endowed with GC content. Further, our results showed the occurrence of both coincidence and antagonistic codon usage patterns among HEV-hosts. The findings further emphasized that both mutational and selection forces influenced the codon usage patterns of ORF4 protein genes.

CONCLUSIONS

To the best of our knowledge, this is first bioinformatics study evaluating codon usage patterns in HEV ORF4 protein genes. The findings from this study are expected to increase our understanding toward significant factors involved in evolutionary changes of ORF4.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s43088-022-00244-w.

Oral mRNA Vaccines Against Infectious Diseases- A Bacterial Perspective [Invited]

The mRNA vaccines from Pfizer/BioNTech and Moderna were granted emergency approval in record time in the history of vaccinology and played an instrumental role in limiting the pandemic caused by SARS-CoV-2. The success of these vaccines resulted from over 3 decades of research from many scientists. However, the development of orally administrable mRNA vaccine development is surprisingly underexplored. Our group specializing in Salmonella-based vaccines explored the possibility of oral mRNA vaccine development. Oral delivery was made possible by the exploitation of the Semliki Forest viral replicon and Salmonella vehicle for transgene amplification and gene delivery, respectively. Herein we highlight the prospect of developing oral replicon-based mRNA vaccines against infectious diseases based on our recent primary studies on SARS-CoV-2. Further, we discuss the potential advantages and limitations of bacterial gene delivery.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

Synonymous Variants: Necessary Nuance in our Understanding of Cancer Drivers and Treatment Outcomes

Once called "silent mutations" and assumed to have no effect on protein structure and function, synonymous variants are now recognized to be drivers for some cancers. There have been significant advances in our understanding of the numerous mechanisms by which synonymous single nucleotide variants (sSNVs) can affect protein structure and function by affecting pre-mRNA splicing, mRNA expression, stability, folding, miRNA binding, translation kinetics, and co-translational folding. This review highlights the need for considering sSNVs in cancer biology to gain a better understanding of the genetic determinants of human cancers and to improve their diagnosis and treatment. We surveyed the literature for reports of sSNVs in cancer and found numerous studies on the consequences of sSNVs on gene function with supporting in vitro evidence. We also found reports of sSNVs that have statistically significant associations with specific cancer types but for which in vitro studies are lacking to support the reported associations. Additionally, we found reports of germline and somatic sSNVs that were observed in numerous clinical studies and for which in silico analysis predicts possible effects on gene function. We provide a review of these investigations and discuss necessary future studies to elucidate the mechanisms by which sSNVs disrupt protein function and are play a role in tumorigeneses, cancer progression, and treatment efficacy. As splicing dysregulation is one of the most well recognized mechanisms by which sSNVs impact protein function, we also include our own in silico analysis for predicting which sSNVs may disrupt pre-mRNA splicing.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Thermodynamics of co-translational folding and ribosome-nascent chain interactions

Proteins can begin the conformational search for their native structure in parallel with biosynthesis on the ribosome, in a process termed co-translational folding. In contrast to the reversible folding of isolated domains, as a nascent chain emerges from the ribosome exit tunnel during translation the free energy landscape it explores also evolves as a function of chain length. While this presents a substantially more complex measurement problem, this review will outline the progress that has been made recently in understanding, quantitatively, the process by which a nascent chain attains its full native stability, as well as the mechanisms through which interactions with the nearby ribosome surface can perturb or modulate this process.

Advances in COVID-19 mRNA vaccine development

To date, the coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has determined 399,600,607 cases and 5,757,562 deaths worldwide. COVID-19 is a serious threat to human health globally. The World Health Organization (WHO) has declared COVID-19 pandemic a major public health emergency. Vaccination is the most effective and economical intervention for controlling the spread of epidemics, and consequently saving lives and protecting the health of the population. Various techniques have been employed in the development of COVID-19 vaccines. Among these, the COVID-19 messenger RNA (mRNA) vaccine has been drawing increasing attention owing to its great application prospects and advantages, which include short development cycle, easy industrialization, simple production process, flexibility to respond to new variants, and the capacity to induce better immune response. This review summarizes current knowledge on the structural characteristics, antigen design strategies, delivery systems, industrialization potential, quality control, latest clinical trials and real-world data of COVID-19 mRNA vaccines as well as mRNA technology. Current challenges and future directions in the development of preventive mRNA vaccines for major infectious diseases are also discussed.

Bias at the third nucleotide of codon pairs in virus and host genomes

Genomes of different sizes and complexity can be compared using common features. Most genomes contain open reading frames, and most genomes use the same genetic code. Redundancy in the genetic code means that different biases in the third nucleotide position of a codon exist in different genomes. However, the nucleotide composition of viruses can be quite different from host nucleotide composition making it difficult to assess the relevance of these biases. Here we show that grouping codons of a codon-pair according to the GC content of the first two nucleotide positions of each codon reveals patterns in nucleotide usage at the third position of the 1st codon. Differences between the observed and expected biases occur predominantly when the first two nucleotides of the 2nd codon are both S (strong, G or C) or both W (weak, A or T), not a mixture of strong and weak. The data indicates that some codon pairs are preferred because of the strength of the interactions between the codon and anticodon, the adjacent tRNAs and the ribosome. Using base-pairing strength and third position bias facilitates the comparison of genomes of different size and nucleotide composition and reveals patterns not previously described.

Immunogenicity mechanism of mRNA vaccines and their limitations in promoting adaptive protection against SARS-CoV-2

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19) in late 2019, hundreds of millions of people have been infected worldwide. There have been unprecedented efforts in acquiring effective vaccines to confer protection against the disease. mRNA vaccines have emerged as promising alternatives to conventional vaccines due to their high potency with the capacity for rapid development and low manufacturing costs. In this review, we summarize the currently available vaccines against SARS-CoV-2 in development, with the focus on the concepts of mRNA vaccines, their antigen selection, delivery and optimization to increase the immunostimulatory capability of mRNA as well as its stability and translatability. We also discuss the host immune responses to the SARS-CoV-2 infection and expound in detail, the adaptive immune response upon immunization with mRNA vaccines, in which high levels of spike-specific IgG and neutralizing antibodies were detected after two-dose vaccination. mRNA vaccines have been shown to induce a robust CD8⁺T cell response, with a balanced CD4⁺ T_H1/T_H2 response. We further discuss the challenges and limitations of COVID-19 mRNA vaccines, where newly emerging variants of SARS-CoV-2 may render currently deployed vaccines less effective. Imbalanced and inappropriate inflammatory responses, resulting from hyper-activation of pro-inflammatory cytokines, which may lead to vaccine-associated enhanced respiratory disease (VAERD) and rare cases of myocarditis and pericarditis also are discussed.

GroEL/S overexpression helps to purge deleterious mutations and reduce genetic diversity during adaptive protein evolution

Chaperones are proteins that help other proteins fold. They also affect the adaptive evolution of their client proteins by buffering the effect of deleterious mutations and increasing the genetic diversity of evolving proteins. We study how the bacterial chaperone GroE (GroEL + GroES) affects the evolution of green fluorescent protein (GFP). To this end we subjected GFP to multiple rounds of mutation and selection for its color phenotype in four replicate E. coli populations, and studied its evolutionary dynamics through high-throughput sequencing and mutant engineering. We evolved GFP both under stabilizing selection for its ancestral (green) phenotype, and to directional selection for a new (cyan) phenotype. We did so both under low and high expression of the chaperone GroE. In contrast to previous work, we observe that GroE does not just buffer but also helps purge deleterious (fluorescence reducing) mutations from evolving populations. In doing so, GroE helps reduce the genetic diversity of evolving populations. In addition, it causes phenotypic heterogeneity in mutants with the same genotype, helping to enhance their fluorescence in some cells, and reducing it in others. Our observations show that chaperones can affect adaptive evolution in more than one way.

Ribosome slowdown triggers codon-mediated mRNA decay independently of ribosome quality control

The control of mRNA stability plays a central role in regulating gene expression patterns. Recent studies have revealed that codon composition in the open reading frame determines mRNA stability in multiple organisms. Based on genome-wide correlation approaches, this previously unrecognized role for the genetic code is attributable to the kinetics of the codon-decoding process by the ribosome. However, complementary experimental analyses are required to clarify the codon effects on mRNA stability and the related cotranslational mRNA decay pathways, for example, those triggered by aberrant ribosome stalling. In the current study, we performed a set of reporter-based analyses to define codon-mediated mRNA decay and ribosome stall-dependent mRNA decay in zebrafish embryos. Our analysis showed that the effect of codons on mRNA stability stems from the decoding process, independent of the ribosome quality control factor Znf598 and stalling-dependent mRNA decay. We propose that codon-mediated mRNA decay is rather triggered by transiently slowed ribosomes engaging in a productive translation cycle in zebrafish embryos.

Role of cotranslational folding for β-sheet-enriched proteins: A perspective from molecular dynamics simulations

The formations of correct three-dimensional structures of proteins are essential to their functions. Cotranslational folding is vital for proteins to form correct structures in vivo. Although some experiments have shown that cotranslational folding can improve the efficiency of folding, its microscopic mechanism is not yet clear. Previously, we built a model of the ribosomal exit tunnel and investigated the cotranslational folding of a three-helix protein by using all-atom molecular dynamics simulations. Here we study the cotranslational folding of three β-sheet-enriched proteins using the same method. The results show that cotranslational folding can enhance the helical population in most cases and reduce non-native long-range contacts before emerging from the ribosomal exit tunnel. After exiting the tunnel, all proteins fall into local minimal states and the structural ensembles of cotranslational folding show more helical conformations than those of free folding. In particular, for one of the three proteins, the GTT WW domain, we find that one local minimum state of the cotranslational folding is the known folding intermediate, which is not found in free folding. This result suggests that the cotranslational folding may increase the folding efficiency by accelerating the sampling more than by avoiding the misfolded state, which is presently a mainstream viewpoint.

Transcription-wide mapping of dihydrouridine reveals that mRNA dihydrouridylation is required for meiotic chromosome segregation

The epitranscriptome has emerged as a new fundamental layer of control of gene expression. Nevertheless, the determination of the transcriptome-wide occupancy and function of RNA modifications remains challenging. Here we have developed Rho-seq, an integrated pipeline detecting a range of modifications through differential modification-dependent rhodamine labeling. Using Rho-seq, we confirm that the reduction of uridine to dihydrouridine (D) by the Dus reductase enzymes targets tRNAs in E. coli and fission yeast. We find that the D modification is also present on fission yeast mRNAs, particularly those encoding cytoskeleton-related proteins, which is supported by large-scale proteome analyses and ribosome profiling. We show that the α-tubulin encoding mRNA nda2 undergoes Dus3-dependent dihydrouridylation, which affects its translation. The absence of the modification on nda2 mRNA strongly impacts meiotic chromosome segregation, resulting in low gamete viability. Applying Rho-seq to human cells revealed that tubulin mRNA dihydrouridylation is evolutionarily conserved.

Ageing exacerbates ribosome pausing to disrupt cotranslational proteostasis

Ageing is accompanied by a decline in cellular proteostasis, which underlies many age-related protein misfolding diseases^1,2. Yet, how ageing impairs proteostasis remains unclear. As nascent polypeptides represent a substantial burden on the proteostasis network³, we hypothesized that altered translational efficiency during ageing could help to drive the collapse of proteostasis. Here we show that ageing alters the kinetics of translation elongation in both Caenorhabditis elegans and Saccharomyces cerevisiae. Ribosome pausing was exacerbated at specific positions in aged yeast and worms, including polybasic stretches, leading to increased ribosome collisions known to trigger ribosome-associated quality control (RQC)^4-6. Notably, aged yeast cells exhibited impaired clearance and increased aggregation of RQC substrates, indicating that ageing overwhelms this pathway. Indeed, long-lived yeast mutants reduced age-dependent ribosome pausing, and extended lifespan correlated with greater flux through the RQC pathway. Further linking altered translation to proteostasis collapse, we found that nascent polypeptides exhibiting age-dependent ribosome pausing in C. elegans were strongly enriched among age-dependent protein aggregates. Notably, ageing increased the pausing and aggregation of many components of proteostasis, which could initiate a cycle of proteostasis collapse. We propose that increased ribosome pausing, leading to RQC overload and nascent polypeptide aggregation, critically contributes to proteostasis impairment and systemic decline during ageing.

A switch from α-helical to β-strand conformation during co-translational protein folding

Cellular proteins begin to fold as they emerge from the ribosome. The folding landscape of nascent chains is not only shaped by their amino acid sequence but also by the interactions with the ribosome. Here, we combine biophysical methods with cryo‐EM structure determination to show that folding of a β‐barrel protein begins with formation of a dynamic α‐helix inside the ribosome. As the growing peptide reaches the end of the tunnel, the N‐terminal part of the nascent chain refolds to a β‐hairpin structure that remains dynamic until its release from the ribosome. Contacts with the ribosome and structure of the peptidyl transferase center depend on nascent chain conformation. These results indicate that proteins may start out as α‐helices inside the tunnel and switch into their native folds only as they emerge from the ribosome. Moreover, the correlation of nascent chain conformations with reorientation of key residues of the ribosomal peptidyl‐transferase center suggest that protein folding could modulate ribosome activity.

CHARMING: Harmonizing synonymous codon usage to replicate a desired codon usage pattern

There is a growing appreciation that synonymous codon usage, although historically regarded as phenotypically silent, can instead alter a wide range of mechanisms related to functional protein production, a term we use here to describe the net effect of transcription (mRNA synthesis), mRNA half-life, translation (protein synthesis) and the probability of a protein folding correctly to its active, functional structure. In particular, recent discoveries have highlighted the important role that sub-optimal codons can play in modifying co-translational protein folding. These results have drawn increased attention to the patterns of synonymous codon usage within coding sequences, particularly in light of the discovery that these patterns can be conserved across evolution for homologous proteins. Because synonymous codon usage differs between organisms, for heterologous gene expression it can be desirable to make synonymous codon substitutions to match the codon usage pattern from the original organism in the heterologous expression host. Here we present CHARMING (for Codon HARMonizING), a robust and versatile algorithm to design mRNA sequences for heterologous gene expression and other related codon harmonization tasks. CHARMING can be run as a downloadable Python script or via a web portal at http://www.codons.org.

Clinical and molecular characterization of a patient with MBTPS1 related spondyloepiphyseal dysplasia: Evidence of pathogenicity for a synonymous variant

BACKGROUND

A novel autosomal recessive skeletal dysplasia resulting from pathogenic variants in membrane-bound transcription factor peptidase, site 1 (MBTPS1) has been recently delineated. To date, only three patients have been reported.

METHODS

In this study, we reported the clinical and molecular features of a Chinese boy who was diagnosed with spondyloepiphyseal dysplasia. The effects of variants on mRNA splicing were analyzed through transcript analysis in vivo and minigene splice assay in vitro.

RESULTS

The proband mainly showed short stature, special facial features, cataract, hernias, and serious sleep apnea syndrome. Growth hormone stimulation tests suggested the boy had growth hormone deficiency. Imaging examinations suggested abnormal thoracolumbar vertebrae and severely decreased bone mineral density. Genetic analysis of MBTPS1 gene revealed two novel heterozygous variants, a nonsense mutation c.2656C > T (p.Q886*, 167) in exon 20 and a synonymous variant c.774C > T (p.A258=) in exon 6. The transcript analysis in vivo exhibited that the synonymous variant c.774C > T caused exon 6 skipping. The minigene splice assay in vitro confirmed the alteration of MBTPS1 mRNA splicing and the exon skipping was partially restored by an antisense oligonucleotide (ASO) treatment.

CONCLUSION

Notably, we report a Chinese rare case of spondyloepiphyseal dysplasia and validate its pathogenic synonymous variant in the MBTPS1 gene.

Membrane-dependent relief of translation elongation arrest on pseudouridine- and N1-methyl-pseudouridine-modified mRNAs

Expression of therapeutically important proteins has benefited dramatically from the advent of chemically modified mRNAs that feature decreased lability and immunogenicity. This had a momentous effect on the rapid development of COVID-19 mRNA vaccines. Incorporation of the naturally occurring pseudouridine (Ψ) or N1-methyl-pseudouridine (N1mΨ) into in vitro transcribed mRNAs prevents the activation of unwanted immune responses by blocking eIF2α phosphorylation, which inhibits translation. Here, we report that Ψs in luciferase (Luc) mRNA exacerbate translation pausing in nuclease-untreated rabbit reticulocyte lysate (uRRL) and promote the formation of high-order-ribosome structures. The major deceleration of elongation occurs at the Ψ-rich nucleotides 1294-1326 of Ψ-Luc mRNA and results in premature termination of translation. The impairment of translation is mainly due to the shortage of membranous components. Supplementing uRRL with canine microsomal membranes (CMMs) relaxes the impediments to ribosome movement, resolves collided ribosomes, and greatly enhances full-size luciferase production. CMMs also strongly stimulated an extremely inefficient translation of N1mΨ-Luc mRNA in uRRL. Evidence is presented that translational pausing can promote membrane recruitment of polysomes with nascent polypeptides that lack a signal sequence. Our results highlight an underappreciated role of membrane binding to polysomes in the prevention of ribosome collision and premature release of nascent polypeptides.

Identification of conserved slow codons that are important for protein expression and function

ABSTRASTDue to the redundancy of the genetic code most amino acids are encoded by several 'synonymous' codons. These codons are used unevenly, and each organism demonstrates its own unique codon usage bias, where the 'preferred' codons are associated with tRNAs that are found in high concentrations. Therefore, for decades, the prevailing view had been that preferred and non-preferred codons are linked to high or slow translation rates, respectively.However, this simplified view is contrasted by the frequent failures of codon-optimization efforts and by evidence of non-preferred (i.e. 'slow') codons having specific roles important for efficient production of functional proteins. One such specific role of slower codons is the regulation of co-translational protein folding, a complex biophysical process that is very challenging to model or to measure.Here, we combined a genome-wide approach with experiments to investigate the role of slow codons in protein production and co-translational folding. We analysed homologous gene groups from divergent bacteria and identified positions of inter-species conservation of bias towards slow codons. We then generated mutants where the conserved slow codons are substituted with 'fast' ones, and experimentally studied the effects of these codon substitutions. Using cellular and biochemical approaches we find that at certain locations, slow-to-fast codon substitutions reduce protein expression, increase protein aggregation, and impair protein function.This report provides an approach for identifying functionally relevant regions with slower codons and demonstrates that such codons are important for protein expression and function.

Tracing protein and proteome history with chronologies and networks: folding recapitulates evolution

INTRODUCTION

While the origin and evolution of proteins remain mysterious, advances in evolutionary genomics and systems biology are facilitating the historical exploration of the structure, function and organization of proteins and proteomes. Molecular chronologies are series of time events describing the history of biological systems and subsystems and the rise of biological innovations. Together with time-varying networks, these chronologies provide a window into the past.

AREAS COVERED

Here, we review molecular chronologies and networks built with modern methods of phylogeny reconstruction. We discuss how chronologies of structural domain families uncover the explosive emergence of metabolism, the late rise of translation, the co-evolution of ribosomal proteins and rRNA, and the late development of the ribosomal exit tunnel; events that coincided with a tendency to shorten folding time. Evolving networks described the early emergence of domains and a late 'big bang' of domain combinations.

EXPERT OPINION

Two processes, folding and recruitment appear central to the evolutionary progression. The former increases protein persistence. The later fosters diversity. Chronologically, protein evolution mirrors folding by combining supersecondary structures into domains, developing translation machinery to facilitate folding speed and stability, and enhancing structural complexity by establishing long-distance interactions in novel structural and architectural designs.