Synonymous Codons Direct Cotranslational Folding toward Different Protein Conformations

Modulation of protein activity by small RNA base pairing internal to coding sequences

Most characterized interactions between bacterial small RNAs (sRNAs) and their target mRNAs occur near ribosome binding sites, resulting in changes in translation initiation or target mRNA decay. To understand the consequences of sRNA pairing internal to coding sequences detected by global RNA-RNA interactome approaches, we examined the impact of sRNA overexpression on seven target proteins. Overexpression of the sRNA led to decreased target protein levels for two pairs, but there were no differences for the others. By further examining ArcZ-ligA and ArcZ-hemK, we discovered that ArcZ pairing with the mRNAs leads to translation pausing and increased protein activity. A ligA point mutation that eliminates sRNA pairing resulted in increased sensitivity to DNA damage, revealing the physiological consequences of the regulation. Thus, regulatory RNA pairing in coding sequences can locally slow translation elongation, likely impacting co-translational protein folding and allowing improved incorporation of co-factors or more optimal folding under specific conditions.

Cotranslational Protein Folding Through Non-Native Structural Intermediates

Cotranslational protein folding follows a distinct pathway shaped by the vectorial emergence of the peptide and spatial constraints of the ribosome exit tunnel. Variations in translation rhythm can cause misfolding linked to disease; however, predicting cotranslational folding pathways remains challenging. Here we computationally predict and experimentally validate a vectorial hierarchy of folding resolved at the atomistic level, where early intermediates are stabilized through non-native hydrophobic interactions before rearranging into the native-like fold. Disrupting these interactions destabilizes intermediates and impairs folding. The chaperone Trigger Factor alters the cotranslational folding pathway by keeping the nascent peptide dynamic until the full domain emerges. Our results highlight an unexpected role of surface-exposed residues in protein folding on the ribosome and provide tools to improve folding prediction and protein design.

DNA Sequence Changes Resulting from Codon Optimization Affect Gene Expression in Pichia pastoris by Altering Chromatin Accessibility

Codon optimization is a widely employed strategy to enhance protein expression. However, it occasionally leads to unexpected transcriptional repression despite preserving amino acid sequences. This study investigates the mechanistic basis of such transcriptional attenuation by analyzing two gene candidates (0432 and Fluc) in the common expression chassis P. pastoris. Both genes experienced severe mRNA reduction following codon optimization. Evidenced by histone H3 chromatin immunoprecipitation (ChIP) and a DNase I hypersensitivity assay, gene sequences with transcriptional repression displayed elevated nucleosome occupancy and reduced chromatin accessibility. The above change was caused by an ORF sequence change independent of the promoter, since transcriptional attenuation and compromised chromatin accessibility were still observed after replacing the strong promoter PGAP with Ppor1 or Prps8b. Our findings challenge the conventional view of codon optimization as solely translation-centric, revealing its capacity to preemptively modulate transcription through chromatin accessibility. This work underscores the necessity of integrating chromatin-level considerations into synthetic gene design to avoid unintended transcriptional silencing and optimize expression outcomes.

Variation in the fitness impact of translationally optimal codons among animals

Early studies in invertebrate model organisms (fruit flies, nematodes) showed that their synonymous codon usage is under selective pressure to optimize translation efficiency in highly expressed genes (a process called translational selection). In contrast, mammals show little evidence of selection for translationally optimal codons. To understand this difference, we examined the use of synonymous codons in 223 metazoan species, covering a wide range of animal clades. For each species, we predicted the set of optimal codons based on the pool of tRNA genes present in its genome, and we analyzed how the frequency of optimal codons correlates with gene expression to quantify the intensity of translational selection (S). We observed that few metazoans show clear signs of translational selection. As predicted by the nearly neutral theory, the highest values of S are observed in species with large effective population sizes (N e). Overall, however, N e appears to be a poor predictor of the intensity of translational selection, suggesting important differences in the fitness effect of synonymous codon usage across taxa. We propose that the few animal taxa that are clearly affected by translational selection correspond to organisms with strong constraints for a very rapid growth rate.

Interpreting mammalian synonymous site conservation in light of the unwanted transcript hypothesis

Mammalian genomes are biased towards GC bases at third codon positions, likely due to a GC-biased ancestral genome and the selectively neutral recombination-related process of GC-biased gene conversion. The unwanted transcript hypothesis posits that this high GC content at synonymous sites may be beneficial for protecting against spurious transcripts, particularly in species with low effective population sizes. Utilising a 240 placental mammal genome alignment and single-base resolution conservation scores, we interpret sequence conservation at mammalian four-fold degenerate sites in this context and find evidence in support of the unwanted transcript hypothesis, including a strong GC bias, high conservation at sites relating to exon splicing, less human genetic variation at conserved four-fold degenerate sites, and conservation of sites important for epigenetic regulation of developmental genes. Additionally, we show that high conservation of four-fold degenerate sites in essential developmental genes, including homeobox genes, likely relates to the low mutation rates experienced by these genes.

Role of silent mutations in KRAS -mutant tumors

ABSTRACT

Silent mutations within the RAS  gene have garnered increasing attention for their potential roles in tumorigenesis and therapeutic strategies. Kirsten-RAS ( KRAS ) mutations, predominantly oncogenic, are pivotal drivers in various cancers. While extensive research has elucidated the molecular mechanisms and biological consequences of active KRAS  mutations, the functional significance of silent mutations remains relatively understudied. This review synthesizes current knowledge on KRAS  silent mutations, highlighting their impact on cancer development. Silent mutations, which do not alter protein sequences but can affect RNA stability and translational efficiency, pose intriguing questions regarding their contribution to tumor biology. Understanding these mutations is crucial for comprehensively unraveling KRAS -driven oncogenesis and exploring novel therapeutic avenues. Moreover, investigations into the clinical implications of silent mutations in KRAS -mutant tumors suggest potential diagnostic and therapeutic strategies. Despite being in early stages, research on KRAS  silent mutations holds promise for uncovering novel insights that could inform personalized cancer treatments. In conclusion, this review underscores the evolving landscape of KRAS  silent mutations, advocating for further exploration to bridge fundamental biology with clinical applications in oncology.

Comprehensive analysis of 111 Pleuronectiformes mitochondrial genomes: insights into structure, conservation, variation and evolution

BACKGROUND

Pleuronectiformes, also known as flatfish, are important model and economic animals. However, a comprehensive genome survey of their important organelles, mitochondria, has been limited. Therefore, we aim to analyze the genomic structure, codon preference, nucleotide diversity, selective pressure and repeat sequences, as well as reconstruct the phylogenetic relationship using the mitochondrial genomes of 111 flatfish species.

RESULTS

Our analysis revealed a conserved gene content of protein-coding genes and rRNA genes, but varying numbers of tRNA genes and control regions across species. Various gene rearrangements were found in flatfish species, especially for the rearrangement of nad5-nad6-cytb block in Samaridae family, the swapping rearrangement of nad6 and cytb gene in Bothidae family, as well as the control region translocation and tRNA-Gln gene inversion in the subfamily Cynoglossinae, suggesting their unique evolutionary history and/or functional benefit. Codon usage showed obvious biases, with adenine being the most frequent nucleotide at the third codon position. Nucleotide diversity and selective pressure analysis suggested that different protein-coding genes underwent varying degrees of evolutionary pressure, with cytb and cox genes being the most conserved ones. Phylogenetic analysis using both whole mitogenome information and concatenated independently aligned protein-coding genes largely mirrored the taxonomic classification of the species, but showed different phylogeny. The identification of simple sequence repeats and various long repetitive sequences provided additional complexity of genome organization and offered markers for evolutionary studies and breeding practices.

CONCLUSIONS

This study represents a significant step forward in our comprehension of the flatfish mitochondrial genomes, providing valuable insights into the structure, conservation and variation within flatfish mitogenomes, with implications for understanding their evolutionary history, functional genomics and fisheries management. Future research can delve deeper into conservation biology, evolutionary biology and functional usages of variations.

Predicting gene sequences with AI to study codon usage patterns

Selective pressure acts on the codon use, optimizing multiple, overlapping signals that are only partially understood. We trained AI models to predict codons given their amino acid sequence in the eukaryotes Saccharomyces cerevisiae and Schizosaccharomyces pombe and the bacteria Escherichia coli and Bacillus subtilis to study the extent to which we can learn patterns in naturally occurring codons to improve predictions. We trained our models on a subset of the proteins and evaluated their predictions on large, separate sets of proteins of varying lengths and expression levels. Our models significantly outperformed naïve frequency-based approaches, demonstrating that there are learnable dependencies in evolutionary-selected codon usage. The prediction accuracy advantage of our models is greater for highly expressed genes and is greater in bacteria than eukaryotes, supporting the hypothesis that there is a monotonic relationship between selective pressure for complex codon patterns and effective population size. In S. cerevisiae and bacteria, our models were more accurate for longer proteins, suggesting that the learned patterns may be related to cotranslational folding. Gene functionality and conservation were also important determinants that affect the performance of our models. Finally, we showed that using information encoded in homologous proteins has only a minor effect on prediction accuracy, perhaps due to complex codon-usage codes in genes undergoing rapid evolution. Our study employing contemporary AI methods offers a unique perspective and a deep-learning-based prediction tool for evolutionary-selected codons. We hope that these can be useful to optimize codon usage in endogenous and heterologous proteins.

Yeast Genome Mutagenesis With Multi-Mismatch PCR: A Rapid and Efficient Strategy for Site-Directed Mutagenesis in Saccharomyces cerevisiae

Saccharomyces cerevisiae (S. cerevisiae) provides an array of cost-effective and time-efficient methods for diverse genome modifications. Among these techniques, site-directed mutagenesis of target genes is a powerful strategy to elucidate intricate structure-function relationships and create specific mutations. While various PCR-based and CRISPR/Cas9-based methods have been developed for introducing point mutations into the S. cerevisiae genome, they often involve multiple steps. In this study, we presented a rapid and effective site-directed mutagenesis strategy using one-step multi-mismatch PCR, termed Yeast Genome Mutagenesis with Multi-mismatch PCR (YGMMP). YGMMP incorporated multiple synonymous mutations proximal to the target point mutations, along with a selection marker cassette and flanking homologous sequences, into the gene segment spanning from the desired mutation to the gene's terminus through overlap PCR. The resulting PCR product was introduced into yeast cells to facilitate the selection of target variants. As a proof of concept, we applied YGMMP to generate an ADE2 mutant. The results demonstrated that the introduction of five and nine synonymous mutations, in addition to the desired single-point mutation, yielded mutagenesis efficiencies of approximately 20% and 30%, respectively. This rapid, straightforward, and efficient method has the potential to greatly simplify site-specific modifications within the S. cerevisiae genome.

mRNA vaccine sequence and structure design and optimization: Advances and challenges

Messenger RNA (mRNA) vaccines have emerged as a powerful tool against communicable diseases and cancers, as demonstrated by their huge success during the coronavirus disease 2019 (COVID-19) pandemic. Despite the outstanding achievements, mRNA vaccines still face challenges such as stringent storage requirements, insufficient antigen expression, and unexpected immune responses. Since the intrinsic properties of mRNA molecules significantly impact vaccine performance, optimizing mRNA design is crucial in preclinical development. In this review, we outline four key principles for optimal mRNA sequence design: enhancing ribosome loading and translation efficiency through untranslated region (UTR) optimization, improving translation efficiency via codon optimization, increasing structural stability by refining global RNA sequence and extending in-cell lifetime and expression fidelity by adjusting local RNA structures. We also explore recent advancements in computational models for designing and optimizing mRNA vaccine sequences following these principles. By integrating current mRNA knowledge, addressing challenges, and examining advanced computational methods, this review aims to promote the application of computational approaches in mRNA vaccine development and inspire novel solutions to existing obstacles.

RNA encodes physical information

Most amino acids are encoded by multiple codons, making the genetic code degenerate. Synonymous mutations affect protein translation and folding, but their impact on RNA itself is often neglected. We developed a genetic algorithm that introduces synonymous mutations to control the diversity of structures sampled by an mRNA. The behavior of the designed mRNAs reveals a physical code layered in the genetic code. We find that mRNA conformational heterogeneity directs physical properties and functional outputs of RNA-protein complexes and biomolecular condensates. The role of structure and disorder of proteins in biomolecular condensates is well appreciated, but we find that RNA conformational heterogeneity is equally important. This feature of RNA enables both evolution and engineers to build cellular structures with specific material and responsive properties.

Influence of Stereochemistry in a Local Approach for Calculating Protein Conformations

Protein structure prediction is generally based on the use of local conformational information coupled with long-range distance restraints. Such restraints can be derived from the knowledge of a template structure or the analysis of protein sequence alignment in the framework of models arising from the physics of disordered systems. The accuracy of approaches based on sequence alignment, however, is limited in the case where the number of aligned sequences is small. Here, we derive protein conformations using only local conformations knowledge by means of the interval Branch-and-Prune algorithm. The computation efficiency is directly related to the knowledge of stereochemistry (bond angle and ω values) along the protein sequence and, in particular, to the variations of the torsion angle ω. The impact of stereochemistry variations is particularly strong in the case of protein topologies defined from numerous long-range restraints, as in the case of protein of β secondary structures. The systematic enumeration of the conformations improves the efficiency of the calculations. The analysis of DNA codons permits to connect the variations of torsion angle ω to the positions of rare DNA codons.

Genome Galaxy Identified by the Circular Code Theory

The genome galaxy identified in bacteria is studied by expressing the reading frame retrieval (RFR) function according to the YZ-content (GC-, AG- and GT-content) of bacterial codons. We have developed a simple probabilistic model for ambiguous sequences in order to show that the RFR function is a measure of the gene reading frame retrieval. Indeed, the RFR function increases with the ratio of ambiguous sequences and the ratio of ambiguous sequences decreases when the codon usage dispersion increases. The classical GC-content is the best parameter for characterizing the upper arm, which is related to bacterial genes with a low GC-content, and the lower arm, which is related to bacterial genes with a high GC-content. The galaxy center has a GC-content around 0.5. Then, these results are confirmed by expressing the GC-content of bacterial codons as a function of the codon usage dispersion. Finally, the bacterial genome galaxy is better described with the GC3-content in the 3rd codon site compared to the GC1-content and GC2-content in the 1st and 2nd codons sites, respectively. Whereas the codon usage is used extensively by biologists, its dispersion, which is an important parameter to reveal this genome galaxy, is surprisingly little known and unused. Therefore, we have developed a mathematical theory of codon usage dispersion by deriving several formulæ. It shows three important parameters in codon usage: the minimum and maximum codon probabilities and the number of codons with high frequency, i.e. with a probability at least 1/64. By applying this theory to the evolution of the genetic code, we see that bacteria have optimised the number of codons with high frequency to maximise the codon dispersion, thus maximising the capacity to retrieve the reading frame in genes. The derived formulæ of dispersion can be easily extended to any weighted code over a finite alphabet.

Further Evidence for Strong Nonneutrality of Yeast Synonymous Mutations

Although synonymous mutations are commonly assumed neutral or nearly so, recent years have seen reports of fitness effects of synonymous mutations detected under laboratory conditions. In a previous study, we used genome editing to construct thousands of yeast mutants each carrying a synonymous or nonsynonymous mutation in one of 21 genes, and discovered that most synonymous and most nonsynonymous mutations are deleterious. A concern was raised that this observation could be caused by the fitness effects of potential CRISPR/Cas9 off-target edits and/or secondary mutations, and an experiment that would be refractory to such effects was proposed. Using genome sequencing, we here show that no CRISPR/Cas9 off-target editing occurred, although some mutants did carry secondary mutations. Analysis of mutants with negligible effects from secondary mutations and new data collected from the proposed experiment confirms the original conclusion. These findings, along with other reports of fitness effects of synonymous mutations from both case and systematic studies, necessitate a paradigm shift from assuming (near) neutrality of synonymous mutations.

Kinetic mechanism and determinants of EF-P recruitment to translating ribosomes

EF-P is a translation factor that facilitates the formation of peptide bonds between consecutive prolines. Using FRET between EF-P and ribosomal protein bL33, we studied dynamics and specificity of EF-P binding to the ribosome. Our findings reveal that EF-P rapidly scans for a free E site and can bind to any ribosome containing a P-site tRNA, regardless of the ribosome's functional state. The interaction with uL1 is essential for EF-P binding, while the β-Lys modification of EF-P doubles the association rate. Specific interactions with the D-loop of tRNAPro or tRNAfMet and via the β-Lys group with the tRNA in the peptidyl transferase center reduce the rate of EF-P dissociation from the ribosome, providing the specificity for complexes that need help in catalyzing peptide bond formation. The nature of the E-site codon has little effect on EF-P binding kinetics. Although EF-P dissociation is reduced upon recognizing its correct tRNA substrate, it remains sufficiently rapid compared to tRNA translocation and does not affect the translocation rate. These results highlight the importance of EF-P's scanning-engagement mechanism for dynamic substrate recognition during rapid translation.

Codon-Optimized and de novo-Synthesized E-Selectin/AAV2 Dose-Response Study for Vascular Regeneration Gene Therapy

OBJECTIVE

This study focuses on dose-response investigation using a codon-optimized and de novo-synthesized E-Selectin/AAV2 (E-Sel/AAV2) vector in preparation for Investigational New Drug enabling of subsequent clinical studies.

BACKGROUND

Gene therapy is a potential solution for patients suffering from chronic limb-threatening ischemia. Understanding the dose for effective gene delivery is crucial for future Investigational New Drug-enabling studies.

METHODS

Expression of the codon-optimized E-Selectin gene was assessed by flow cytometry following in vitro cell transfection assay and RT-qPCR for murine limbs injected in vivo with AAV-m-E-Selectin (E-Sel/AAV2). Dose-response studies involved 3 cohorts of FVB/NJ mice (n=6/group) with escalating log doses of E-Selectin/AAV2 injected intramuscularly in divided aliquots, ranging from 2 × 10 9 VG to 2 × 10 11 VG, into ischemic limbs created by left femoral artery/vein ligation/excision and administration of nitric oxide synthase inhibitor, L-NAME. Limb perfusion, extent of gangrene free limb, functional limb recovery, and therapeutic angiogenesis were assessed.

RESULTS

Codon-optimized E-Sel/AAV2 gene therapy exhibits a superior expression level than WT E-Sel/AAV2 gene therapy both in vitro and in vivo. Mice treated with a high dose (2 × 10 11 VG) of E-Sel/AAV2 showed significantly improved perfusion indices, lower Faber scores, increased running stamina, and neovascularization compared with lower doses tested with control groups, indicating a distinct dose-dependent response. No toxicity was detected in any of the animal groups studied.

CONCLUSIONS

E-Sel/AAV2 Vascular Regeneration Gene Therapy holds promise for enhancing the recovery of ischemic hindlimb perfusion and function, with the effective dose identified in this study as 2 × 10 11 VG aliquots injected intramuscularly.

Synonymous codon substitutions modulate transcription and translation of a divergent upstream gene by modulating antisense RNA production

Synonymous codons were originally viewed as interchangeable, with no phenotypic consequences. However, substantial evidence has now demonstrated that synonymous substitutions can perturb a variety of gene expression and protein homeostasis mechanisms, including translational efficiency, translational fidelity, and cotranslational folding of the encoded protein. To date, most studies of synonymous codon-derived perturbations have focused on effects within a single gene. Here, we show that synonymous codon substitutions made far within the coding sequence of Escherichia coli plasmid-encoded chloramphenicol acetyltransferase (cat) can significantly increase expression of the divergent upstream tetracycline resistance gene, tetR. In four out of nine synonymously recoded cat sequences tested, expression of the upstream tetR gene was significantly elevated due to transcription of a long antisense RNA (asRNA) originating from a transcription start site within cat. Surprisingly, transcription of this asRNA readily bypassed the native tet transcriptional repression mechanism. Even more surprisingly, accumulation of the TetR protein correlated with the level of asRNA, rather than total tetR RNA. These effects of synonymous codon substitutions on transcription and translation of a neighboring gene suggest that synonymous codon usage in bacteria may be under selection to both preserve the amino acid sequence of the encoded gene and avoid DNA sequence elements that can significantly perturb expression of neighboring genes. Avoiding such sequences may be especially important in plasmids and prokaryotic genomes, where genes and regulatory elements are often densely packed. Similar considerations may apply to the design of genetic circuits for synthetic biology applications.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

Translation Rates and Protein Folding

The mRNA coding sequence defines not only the amino acid sequence of the protein, but also the speed at which the ribosomes move along the mRNA while making the protein. The non-uniform local kinetics - denoted as translational rhythm - is similar among mRNAs coding for related protein folds. Deviations from this conserved rhythm can result in protein misfolding. In this review we summarize the experimental evidence demonstrating how local translation rates affect cotranslational protein folding, with the focus on the synonymous codons and patches of charged residues in the nascent peptide as best-studied examples. Alterations in nascent protein conformations due to disturbed translational rhythm can persist off the ribosome, as demonstrated by the effects of synonymous codon variants of several disease-related proteins. Charged amino acid patches in nascent chains also modulate translation and cotranslational protein folding, and can abrogate translation when placed at the N-terminus of the nascent peptide. During cotranslational folding, incomplete nascent chains navigate through a unique conformational landscape in which earlier intermediate states become inaccessible as the nascent peptide grows. Precisely tuned local translation rates, as well as interactions with the ribosome, guide the folding pathway towards the native structure, whereas deviations from the natural translation rhythm may favor pathways leading to trapped misfolded states. Deciphering the 'folding code' of the mRNA will contribute to understanding the diseases caused by protein misfolding and to rational protein design.

The Effects of Codon Usage on Protein Structure and Folding

The rate of protein synthesis is slower than many folding reactions and varies depending on the synonymous codons encoding the protein sequence. Synonymous codon substitutions thus have the potential to regulate cotranslational protein folding mechanisms, and a growing number of proteins have been identified with folding mechanisms sensitive to codon usage. Typically, these proteins have complex folding pathways and kinetically stable native structures. Kinetically stable proteins may fold only once over their lifetime, and thus, codon-mediated regulation of the pioneer round of protein folding can have a lasting impact. Supporting an important role for codon usage in folding, conserved patterns of codon usage appear in homologous gene families, hinting at selection. Despite these exciting developments, there remains few experimental methods capable of quantifying translation elongation rates and cotranslational folding mechanisms in the cell, which challenges the development of a predictive understanding of how biology uses codons to regulate protein folding.

The effectiveness of selection in a species affects the direction of amino acid frequency evolution

Nearly neutral theory predicts that species with higher effective population size (N e ) are better able to purge slightly deleterious mutations. We compare evolution in high-N e vs. low-N e vertebrates to reveal which amino acid frequencies are subject to subtle selective preferences. We take three complementary approaches, two measuring flux and one measuring outcomes. First, we fit non-stationary substitution models of amino acid flux using maximum likelihood, comparing the high-N e clade of rodents and lagomorphs to its low-N e sister clade of primates and colugos. Second, we compare evolutionary outcomes across a wider range of vertebrates, via correlations between amino acid frequencies and N e . Third, we dissect the details of flux in human, chimpanzee, mouse, and rat, as scored by parsimony - this also enables comparison to a historical paper. All three methods agree on which amino acids are preferred under more effective selection. Preferred amino acids tend to be smaller, less costly to synthesize, and to promote intrinsic structural disorder. Parsimony-induced bias in the historical study produces an apparent reduction in structural disorder, perhaps driven by slightly deleterious substitutions. Within highly exchangeable pairs of amino acids, arginine is strongly preferred over lysine, and valine over isoleucine, consistent with more effective selection preferring a marginally larger free energy of folding. These two preferences match differences between thermophiles and mesophilic relatives. These results reveal the biophysical consequences of mutation-selection-drift balance, and demonstrate the utility of nearly neutral theory for understanding protein evolution.

Progression of monoclonal gammopathy of undetermined significance to multiple myeloma is associated with enhanced translational quality control and overall loss of surface antigens

BACKGROUND

Despite significant advancements in treatment strategies, multiple myeloma remains incurable. Additionally, there is a distinct lack of reliable biomarkers that can guide initial treatment decisions and help determine suitable replacement or adjuvant therapies when relapse ensues due to acquired drug resistance.

METHODS

To define specific proteins and pathways involved in the progression of monoclonal gammopathy of undetermined significance (MGUS) to multiple myeloma (MM), we have applied super-SILAC quantitative proteomic analysis to CD138 + plasma cells from 9 individuals with MGUS and 37 with MM.

RESULTS

Unsupervised hierarchical clustering defined three groups: MGUS, MM, and MM with an MGUS-like proteome profile (ML) that may represent a group that has recently transformed to MM. Statistical analysis identified 866 differentially expressed proteins between MM and MGUS, and 189 between MM and ML, 177 of which were common between MGUS and ML. Progression from MGUS to MM is accompanied by upregulated EIF2 signaling, DNA repair, and proteins involved in translational quality control, whereas integrin- and actin cytoskeletal signaling and cell surface markers are downregulated.

CONCLUSION

Compared to the premalignant plasma cells in MGUS, malignant MM cells apparently have mobilized several pathways that collectively contribute to ensure translational fidelity and to avoid proteotoxic stress, especially in the ER. The overall reduced expression of immunoglobulins and surface antigens contribute to this and may additionally mediate evasion from recognition by the immune apparatus. Our analyses identified a range of novel biomarkers with potential prognostic and therapeutic value, which will undergo further evaluation to determine their clinical significance.

Navigating the complexities of multi-domain protein folding

Proteome complexity has expanded tremendously over evolutionary time, enabling biological diversification. Much of this complexity is achieved by combining a limited set of structural units into long polypeptides. This widely used evolutionary strategy poses challenges for folding of the resulting multi-domain proteins. As a consequence, their folding differs from that of small single-domain proteins, which generally fold quickly and reversibly. Co-translational processes and chaperone interactions are important aspects of multi-domain protein folding. In this review, we discuss some of the recent experimental progress toward understanding these processes.

Ser/Leu-swapped cell-free translation system constructed with natural/in vitro transcribed-hybrid tRNA set

The Ser/Leu-swapped genetic code can act as a genetic firewall, mitigating biohazard risks arising from horizontal gene transfer in genetically modified organisms. Our prior work demonstrated the orthogonality of this swapped code to the standard genetic code using a cell-free translation system comprised of 21 in vitro transcribed tRNAs. In this study, to advance this system for protein engineering, we introduce a natural/in vitro transcribed-hybrid tRNA set. This set combines natural tRNAs from Escherichia coli (excluding Ser, Leu, and Tyr) and in vitro transcribed tRNAs, encompassing anticodon-swapped tRNASerGAG and tRNALeuGGA. This approach reduces the number of in vitro transcribed tRNAs required from 21 to only 4. In this optimized system, the production of a model protein, superfolder green fluorescent protein, increases to 3.5-fold. With this hybrid tRNA set, the Ser/Leu-swapped cell-free translation system will stand as a potent tool for protein production with reduced biohazard concerns in future biological endeavors.

Molecular Evolution of Protein Sequences and Codon Usage in Monkeypox Viruses

The monkeypox virus (mpox virus, MPXV) epidemic in 2022 has posed a significant public health risk. Yet, the evolutionary principles of MPXV remain largely unknown. Here, we examined the evolutionary patterns of protein sequences and codon usage in MPXV. We first demonstrated the signal of positive selection in OPG027, specifically in the Clade I lineage of MPXV. Subsequently, we discovered accelerated protein sequence evolution over time in the variants responsible for the 2022 outbreak. Furthermore, we showed strong epistasis between amino acid substitutions located in different genes. The codon adaptation index (CAI) analysis revealed that MPXV genes tended to use more non-preferred codons compared to human genes, and the CAI decreased over time and diverged between clades, with Clade I > IIa and IIb-A > IIb-B. While the decrease in fatality rate among the three groups aligned with the CAI pattern, it remains unclear whether this correlation was coincidental or if the deoptimization of codon usage in MPXV led to a reduction in fatality rates. This study sheds new light on the mechanisms that govern the evolution of MPXV in human populations.

Two synonymous single-nucleotide polymorphisms promoting fluoroquinolone resistance of Escherichia coli in the environment

Single-nucleotide polymorphism (SNP) is one of the core mechanisms that respond to antibiotic resistance of Escherichia coli (E. coli), which is a major issue in environmental pollution. A specific type of SNPs, synonymous SNPs, have been generally considered as the "silent" SNPs since they do not change the encoded amino acid. However, the impact of synonymous SNPs on mRNA splicing, nucleo-cytoplasmic export, stability, and translation was gradually discovered in the last decades. Figuring out the mechanism of synonymous SNPs in regulating antibiotic resistance is critical to improve antimicrobial therapy strategies in clinics and biological treatment strategies of antibiotic-resistant E. coli-polluted materials. With our newly designed antibiotic resistant SNPs prediction algorithm, Multilocus Sequence Type based Identification for Phenotype-single nucleotide polymorphism Analysis (MIPHA), and in vivo validation, we identified 2 important synonymous SNPs 522 G>A and 972 C>T, located at hisD gene, which was previously predicted as a fluoroquinolone resistance-related gene without a detailed mechanism in the E. coli samples with environmental backgrounds. We first discovered that hisD causes gyrA mutation via the upregulation of sbmC and its downstream gene umuD. Moreover, those 2 synonymous SNPs of hisD cause its own translational slowdown and further reduce the expression levels of sbmC and its downstream gene umuD, making the fluoroquinolone resistance determining region of gyrA remains unmutated, ultimately causing the bacteria to lose their ability to resist drugs. This study provided valuable insight into the role of synonymous SNPs in mediating antibiotic resistance of bacteria and a new perspective for the treatment of environmental pollution caused by drug-resistant bacteria.

Chloroplast genome analysis and evolutionary insights in the versatile medicinal plant Calendula officinalis L

Calendula officinalis L.is a versatile medicinal plant with numerous applications in various fields. However, its chloroplast genome structure, features, phylogeny, and patterns of evolution and mutation remain largely unexplored. This study examines the chloroplast genome, phylogeny, codon usage bias, and divergence time of C. officinalis, enhancing our understanding of its evolution and adaptation. The chloroplast genome of C. officinalis is a 150,465 bp circular molecule with a G + C content of 37.75% and comprises 131 genes. Phylogenetic analysis revealed a close relationship between C. officinalis, C. arvensis, and Osteospermum ecklonis. A key finding is the similarity in codon usage bias among these species, which, coupled with the divergence time analysis, supports their close phylogenetic proximity. This similarity in codon preference and divergence times underscores a parallel evolutionary adaptation journey for these species, highlighting the intricate interplay between genetic evolution and environmental adaptation in the Asteraceae family. Moreover unique evolutionary features in C. officinalis, possibly associated with certain genes were identified, laying a foundation for future research into the genetic diversity and medicinal value of C. officinalis.

Quantifying negative selection on synonymous variants

Widespread adoption of DNA sequencing has resulted in large numbers of genetic variants, whose contribution to disease is not easily determined. Although many types of variation are known to disrupt cellular processes in predictable ways, for some categories of variants, the effects may not be directly detectable. A particular example is synonymous variants, that is, those single-nucleotide variants that create a codon substitution, such that the produced amino acid sequence is unaffected. Contrary to the original theory suggesting that synonymous variants are benign, there is a growing volume of research showing that, despite their "silent" mechanism of action, some synonymous variation may be deleterious. Here, we studied the extent of the negative selective pressure acting on different classes of synonymous variants by analyzing the relative enrichment of synonymous singleton variants in the human exomes provided by gnomAD. Using a modification of the mutability-adjusted proportion of singletons (MAPS) metric as a measure of purifying selection, we found that some classes of synonymous variants are subject to stronger negative selection than others. For instance, variants that reduce codon optimality undergo stronger selection than optimality-increasing variants. Besides, selection affects synonymous variants implicated in splice-site-loss or splice-site-gain events. To understand what drives this negative selection, we tested a number of predictors in the aim to explain the variability in the selection scores. Our findings provide insights into the effects of synonymous variants at the population level, highlighting the specifics of the role that these variants play in health and disease.

Codon-optimization in gene therapy: promises, prospects and challenges

Codon optimization has evolved to enhance protein expression efficiency by exploiting the genetic code's redundancy, allowing for multiple codon options for a single amino acid. Initially observed in E. coli, optimal codon usage correlates with high gene expression, which has propelled applications expanding from basic research to biopharmaceuticals and vaccine development. The method is especially valuable for adjusting immune responses in gene therapies and has the potenial to create tissue-specific therapies. However, challenges persist, such as the risk of unintended effects on protein function and the complexity of evaluating optimization effectiveness. Despite these issues, codon optimization is crucial in advancing gene therapeutics. This study provides a comprehensive review of the current metrics for codon-optimization, and its practical usage in research and clinical applications, in the context of gene therapy.

Essential factors, advanced strategies, challenges, and approaches involved for efficient expression of recombinant proteins in Escherichia coli

Producing recombinant proteins is a major accomplishment of biotechnology in the past century. Heterologous hosts, either eukaryotic or prokaryotic, are used for the production of these proteins. The utilization of microbial host systems continues to dominate as the most efficient and affordable method for biotherapeutics and food industry productions. Hence, it is crucial to analyze the limitations and advantages of microbial hosts to enhance the efficient production of recombinant proteins on a large scale. E. coli is widely used as a host for the production of recombinant proteins. Researchers have identified certain obstacles with this host, and given the growing demand for recombinant protein production, there is an immediate requirement to enhance this host. The following review discusses the elements contributing to the manifestation of recombinant protein. Subsequently, it sheds light on innovative approaches aimed at improving the expression of recombinant protein. Lastly, it delves into the obstacles and optimization methods associated with translation, mentioning both cis-optimization and trans-optimization, producing soluble recombinant protein, and engineering the metal ion transportation. In this context, a comprehensive description of the distinct features will be provided, and this knowledge could potentially enhance the expression of recombinant proteins in E. coli.

Computational design of mRNA vaccines

mRNA technology has emerged as a successful vaccine platform that offered a swift response to the COVID-19 pandemic. Accumulating evidence shows that vaccine efficacy, thermostability, and other important properties, are largely impacted by intrinsic properties of the mRNA molecule, such as RNA sequence and structure, both of which can be optimized. Designing mRNA sequence for vaccines presents a combinatorial problem due to an extremely large selection space. For instance, due to the degeneracy of the genetic code, there are over 10632 possible mRNA sequences that could encode the spike protein, the COVID-19 vaccines' target. Moreover, designing different elements of the mRNA sequence simultaneously against multiple objectives such as translational efficiency, reduced reactogenicity, and improved stability requires an efficient and sophisticated optimization strategy. Recently, there has been a growing interest in utilizing computational tools to redesign mRNA sequences to improve vaccine characteristics and expedite discovery timelines. In this review, we explore important biophysical features of mRNA to be considered for vaccine design and discuss how computational approaches can be applied to rapidly design mRNA sequences with desirable characteristics.

Engineering Saccharomyces cerevisiae for efficient production of recombinant proteins

Saccharomyces cerevisiae is an excellent microbial cell factory for producing valuable recombinant proteins because of its fast growth rate, robustness, biosafety, ease of operability via mature genomic modification technologies, and the presence of a conserved post-translational modification pathway among eukaryotic organisms. However, meeting industrial and market requirements with the current low microbial production of recombinant proteins can be challenging. To address this issue, numerous efforts have been made to enhance the ability of yeast cell factories to efficiently produce proteins. In this review, we provide an overview of recent advances in S. cerevisiae engineering to improve recombinant protein production. This review focuses on the strategies that enhance protein production by regulating transcription through promoter engineering, codon optimization, and expression system optimization. Additionally, we describe modifications to the secretory pathway, including engineered protein translocation, protein folding, glycosylation modification, and vesicle trafficking. Furthermore, we discuss global metabolic pathway optimization and other relevant strategies, such as the disruption of protein degradation, cell wall engineering, and random mutagenesis. Finally, we provide an outlook on the developmental trends in this field, offering insights into future directions for improving recombinant protein production in S. cerevisiae.

An easy tool to monitor the elemental steps of in vitro translation via gel electrophoresis of fluorescently labeled small peptides

Several methods are available to visualize and assess the kinetics and efficiency of elemental steps of protein biosynthesis. However, each of these methods has its own limitations. Here, we present a novel, simple and convenient tool for monitoring stepwise in vitro translation initiated by BODIPY-Met-tRNA. Synthesis and release of very short, 1-7 amino acids, BODIPY-labeled peptides, can be monitored using urea-polyacrylamide gel electrophoresis. Very short BODIPY-labeled oligopeptides might be resolved this way, in contrast to widely used Tris-tricine gel electrophoresis, which is suitable to separate peptides larger than 1 kDa. The method described in this manuscript allows one to monitor the steps of translation initiation, peptide transfer, translocation, and termination as well as their inhibition at an unprecedented single amino acid resolution.

Proteomics revealed an association between ribosome-associated proteins and amyloid beta deposition in Alzheimer's disease

Most scholars believe that amyloid-beta (Aβ) has the potential to induce apoptosis, stimulate an inflammatory cascade, promote oxidative stress and exacerbate the pathological progression of Alzheimer's disease (AD). Therefore, it is crucial to investigate the deposition of Aβ in AD. At approximately 6 months of age, APP/PS1 double transgenic mice gradually exhibit the development of plaques, as well as spatial and learning impairment. Notably, the hippocampus is specifically affected in the course of AD. Herein, 6-month-old APP/PS1 double transgenic mice were utilized, and the differentially expressed (DE) proteins in the hippocampus were identified and analyzed using 4D label-free quantitative proteomics technology and parallel reaction monitoring (PRM). Compared to wild-type mice, 29 proteins were upregulated and 25 proteins were downregulated in the AD group. Gene Ontology (GO) enrichment analysis of biological processes (BP) indicated that the DE proteins were mainly involved in 'ribosomal large subunit biogenesis'. Molecular function (MF) analysis results were primarily associated with '5.8S rRNA binding' and 'structural constituent of ribosome'. In terms of cellular components (CC), the DE proteins were mainly found in 'polysomal ribosome', 'cytosolic large ribosomal subunit', 'cytosolic ribosome', and 'large ribosomal subunit', among others. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis demonstrated that the results were mainly enriched in the 'Ribosome signaling pathway'. The key target proteins identified were ribosomal protein (Rp)l18, Rpl17, Rpl19, Rpl24, Rpl35, and Rpl6. The PRM verification results were consistent with the findings of the 4D label-free quantitative proteomics analysis. Overall, these findings suggest that Rpl18, Rpl17, Rpl19, Rpl24, Rpl35, and Rpl6 may have potential therapeutic value for the treatment of AD by targeting Aβ.

How the ribosome shapes cotranslational protein folding

During protein synthesis, the growing nascent peptide chain moves inside the polypeptide exit tunnel of the ribosome from the peptidyl transferase center towards the exit port where it emerges into the cytoplasm. The ribosome defines the unique energy landscape of the pioneering round of protein folding. The spatial confinement and the interactions of the nascent peptide with the tunnel walls facilitate formation of secondary structures, such as α-helices. The vectorial nature of protein folding inside the tunnel favors local intra- and inter-molecular interactions, thereby inducing cotranslational folding intermediates that do not form upon protein refolding in solution. Tertiary structures start to fold in the lower part of the tunnel, where interactions with the ribosome destabilize native protein folds. The present review summarizes the recent progress in understanding the driving forces of nascent protein folding inside the tunnel and at the surface of the ribosome.

Broken silence: 22,841 predicted deleterious synonymous variants identified in the human exome through computational analysis

Synonymous single nucleotide variants (sSNVs) do not alter the primary structure of a protein, thus it was previously accepted that they were neutral. Recently, several studies demonstrated their significance to a range of diseases. Still, variant prioritization strategies lack focus on sSNVs. Here, we identified 22,841 deleterious synonymous variants in 125,748 human exomes using two in silico predictors (SilVA and CADD). While 98.2% of synonymous variants are classified as neutral, 1.8% are predicted to be deleterious, yielding an average of 9.82 neutral and 0.18 deleterious sSNVs per exome. Further investigation of prediction features via Heterogeneous Ensemble Feature Selection revealed that impact on amino acid sequence and conservation carry the most weight for a deleterious prediction. Thirty nine detrimental sSNVs are not rare and are located on disease associated genes. Ten distinct putatively non-deleterious sSNVs are likely to be under positive selection in the North-Western European and East Asian populations. Taken together our analysis gives voice to the so-called silent mutations as we propose a robust framework for evaluating the deleteriousness of sSNVs in variant prioritization studies.

Synonymous codon usage regulates translation initiation

Nonoptimal synonymous codons repress gene expression, but the underlying mechanisms are poorly understood. We and others have previously shown that nonoptimal codons slow translation elongation speeds and thereby trigger messenger RNA (mRNA) degradation. Nevertheless, transcript levels are often insufficient to explain protein levels, suggesting additional mechanisms by which codon usage regulates gene expression. Using reporters in human and Drosophila cells, we find that transcript levels account for less than half of the variation in protein abundance due to codon usage. This discrepancy is explained by translational differences whereby nonoptimal codons repress translation initiation. Nonoptimal transcripts are also less bound by the translation initiation factors eIF4E and eIF4G1, providing a mechanistic explanation for their reduced initiation rates. Importantly, translational repression can occur without mRNA decay and deadenylation, and it does not depend on the known nonoptimality sensor, CNOT3. Our results reveal a potent mechanism of regulation by codon usage where nonoptimal codons repress further rounds of translation.

Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene

BACKGROUND

Single synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed.

RESULTS

Synonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo.

CONCLUSION

In summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

Tissue-specific temperature dependence of RNA editing levels in zebrafish

BACKGROUND

RNA editing by adenosine deaminase acting on RNA (ADAR) occurs in all metazoans and fulfils several functions. Here, we examined effects of acclimation temperature (27 °C, 18 °C,13 °C) on editing patterns in six tissues of zebrafish (Danio rerio).

RESULTS

Sites and total amounts of editing differed among tissues. Brain showed the highest levels, followed by gill and skin. In these highly edited tissues, decreases in temperatures led to large increases in total amounts of editing and changes in specific edited sites. Gene ontology analysis showed both similarities (e.g., endoplasmic reticulum stress response) and differences in editing among tissues. The majority of edited sites were in transcripts of transposable elements and the 3'UTR regions of protein coding genes. By experimental validation, translation efficiency was directly related to extent of editing of the 3'UTR region of an mRNA.

CONCLUSIONS

RNA editing increases 3'UTR polymorphism and affects efficiency of translation. Such editing may lead to temperature-adaptive changes in the proteome through altering relative amounts of synthesis of different proteins.

Combination of hydrophobicity and codon usage bias determines sorting of model K+ channel protein to either mitochondria or endoplasmic reticulum

When the K+ channel-like protein Kesv from Ectocarpus siliculosus virus 1 is heterologously expressed in mammalian cells, it is sorted to the mitochondria. This targeting can be redirected to the endoplasmic reticulum (ER) by altering the codon usage in distinct regions of the gene or by inserting a triplet of hydrophobic amino acids (AAs) into the protein's C-terminal transmembrane domain (ct-TMD). Systematic variations in the flavor of the inserted AAs and/or its codon usage show that a positive charge in the inserted AA triplet alone serves as strong signal for mitochondria sorting. In cases of neutral AA triplets, mitochondria sorting are favored by a combination of hydrophilic AAs and rarely used codons; sorting to the ER exhibits the inverse dependency. This propensity for ER sorting is particularly high when a common codon follows a rarer one in the AA triplet; mitochondria sorting in contrast is supported by codon uniformity. Since parameters like positive charge, hydrophobic AAs, and common codons are known to facilitate elongation of nascent proteins in the ribosome the data suggest a mechanism in which local changes in elongation velocity and co-translational folding in the ct-TMD influence intracellular protein sorting.

VeloPro: A pipeline integrating Ribo-seq and AlphaFold deciphers association patterns between translation velocity and protein structure features

VeloPro integrates Ribo-seq data and AlphaFold2-predicted 3D protein structure information for characterization of the association patterns between translation velocity and many protein structure features in prokaryotic and eukaryotic organisms across different taxonomical clades such as bacteria, fungi, protozoa, nematode, plants, insect, and mammals. We illustrated that association patterns between translation velocity and protein structure features differ across organisms, partially reflecting their taxonomical relationship.

microRNAs slow translating ribosomes to prevent protein misfolding in eukaryotes

Slower translation rates reduce protein misfolding. Such reductions in speed can be mediated by the presence of non-optimal codons, which allow time for proper folding to occur. Although this phenomenon is conserved from bacteria to humans, it is not known whether there are additional eukaryote-specific mechanisms which act in the same way. MicroRNAs (miRNAs), not present in prokaryotes, target both coding sequences (CDS) and 3' untranslated regions (UTR). Given their low suppressive efficiency, it has been unclear why miRNAs are equally likely to bind to a CDS. Here, we show that miRNAs transiently stall translating ribosomes, preventing protein misfolding with little negative effect on protein abundance. We first analyzed ribosome profiles and miRNA binding sites to examine whether miRNAs stall ribosomes. Furthermore, either global or specific miRNA deficiency accelerated ribosomes and induced aggregation of a misfolding-prone polypeptide reporter. These defects were rescued by slowing ribosomes using non-cleaving shRNAs as miRNA mimics. We finally show that proinsulin misfolding, associated with type II diabetes, was resolved by non-cleaving shRNAs. Our findings provide a eukaryote-specific mechanism of co-translational protein folding and a previously unknown mechanism of action to target protein misfolding diseases.

Active role of the protein translation machinery in protecting against stress tolerance in Synechococcus elongatus PCC7942

In vivo protein synthesis is crucial for all domains of life. It is accomplished through translational machinery, and a key step is the translocation of tRNA-mRNA by elongation factor G (EF-G). Genome-based analysis revealed two EF-G encoding genes (S0885 and S2082) in the freshwater cyanobacterium model Synechococcus elongatus PCC7942. S0885 is the essential EF-G gene for photosynthesis. We generated a strain of S. elongatus PCC7942 that overexpressed S0885 (OX-S0885) to identify EF-G functionality. RT-PCR and Western blot analyses revealed increased transcriptional and translational levels in OX-S0885 at 10.5-13.5 and 2.0-3.0 fold, respectively. Overexpression of S0885 led to an increase in specific growth rate. Additionally, polysome-to-monosome ratio (P/M) and RNA-to-protein ratio (R/P) were elevated in OX-S0885 compared with the empty vector. Interestingly, R/P in OX-S0885 was retained at more than 70% under oxidative stress while R/P in the empty vector was severely depleted, suggesting the maintenance of translation. Thus, S0885 appeared to be the important target of oxidative stress because it was protected by the stress response system to maintain its function. These results suggest that cyanobacterial EF-G has a primary function in translation and an unrelated activity during stress conditions. These findings support the substantial role of EF-G in the formation and maintenance of cellular protein formation, and in the protection of the global translational mechanism under oxidative stress condition.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Future prospects in mRNA vaccine development

The recent advancements in messenger ribonucleic acid (mRNA) vaccine development have vastly enhanced their use as alternatives to conventional vaccines in the prevention of various infectious diseases and treatment of several types of cancers. This is mainly due to their remarkable ability to stimulate specific immune responses with minimal clinical side effects. This review gives a detailed overview of mRNA vaccines currently in use or at various stages of development, the recent advancements in mRNA vaccine development, and the challenges encountered in their development. Future perspectives on this technology are also discussed.

Developmental Selection and the Perception of Mutation Bias

The notion that mutations are random relative to their fitness effects is central to the Neo-Darwinian view of evolution. However, a recent interpretation of the patterns of mutation accumulation in the genome of Arabidopsis thaliana has challenged this notion, arguing for the presence of a targeted DNA repair mechanism that causes a nonrandom association of mutation rates and fitness effects. Specifically, this mechanism was suggested to cause a reduction in the rates of mutations on essential genes, thus lowering the rates of deleterious mutations. Central to this argument were attempts to rule out selection at the population level. Here, we offer an alternative and parsimonious interpretation of the patterns of mutation accumulation previously attributed to mutation bias, showing how they can instead or additionally be caused by developmental selection, that is selection occurring at the cellular level during the development of a multicellular organism. Thus, the depletion of deleterious mutations in A. thaliana may indeed be the result of a selective process, rather than a bias in mutation. More broadly, our work highlights the importance of considering development in the interpretation of population-genetic analyses of multicellular organisms, and it emphasizes that efforts to identify mechanisms involved in mutational biases should explicitly account for developmental selection.

Optimization and Deoptimization of Codons in SARS-CoV-2 and Related Implications for Vaccine Development

The spread of coronavirus disease 2019 (COVID-19), caused by severe respiratory syndrome coronavirus 2 (SARS-CoV-2), has progressed into a global pandemic. To date, thousands of genetic variants have been identified among SARS-CoV-2 isolates collected from patients. Sequence analysis reveals that the codon adaptation index (CAI) values of viral sequences have decreased over time but with occasional fluctuations. Through evolution modeling, it is found that this phenomenon may result from the virus's mutation preference during transmission. Using dual-luciferase assays, it is further discovered that the deoptimization of codons in the viral sequence may weaken protein expression during virus evolution, indicating that codon usage may play an important role in virus fitness. Finally, given the importance of codon usage in protein expression and particularly for mRNA vaccines, it is designed several codon-optimized Omicron BA.2.12.1, BA.4/5, and XBB.1.5 spike mRNA vaccine candidates and experimentally validated their high levels of expression. This study highlights the importance of codon usage in virus evolution and provides guidelines for codon optimization in mRNA and DNA vaccine development.

Is Posttranslational Folding More Efficient Than Refolding from a Denatured State: A Computational Study

The folding of proteins into their native conformation is a complex process that has been extensively studied over the past half-century. The ribosome, the molecular machine responsible for protein synthesis, is known to interact with nascent proteins, adding further complexity to the protein folding landscape. Consequently, it is unclear whether the folding pathways of proteins are conserved on and off the ribosome. The main question remains: to what extent does the ribosome help proteins fold? To address this question, we used coarse-grained molecular dynamics simulations to compare the mechanisms by which the proteins dihydrofolate reductase, type III chloramphenicol acetyltransferase, and d-alanine-d-alanine ligase B fold during and after vectorial synthesis on the ribosome to folding from the full-length unfolded state in bulk solution. Our results reveal that the influence of the ribosome on protein folding mechanisms varies depending on the size and complexity of the protein. Specifically, for a small protein with a simple fold, the ribosome facilitates efficient folding by helping the nascent protein avoid misfolded conformations. However, for larger and more complex proteins, the ribosome does not promote folding and may contribute to the formation of intermediate misfolded states cotranslationally. These misfolded states persist posttranslationally and do not convert to the native state during the 6 μs runtime of our coarse-grain simulations. Overall, our study highlights the complex interplay between the ribosome and protein folding and provides insight into the mechanisms of protein folding on and off the ribosome.

Non-synonymous, synonymous, and non-coding nucleotide variants contribute to recurrently altered biological processes during retinoblastoma progression

Retinoblastomas form in response to biallelic RB1 mutations or MYCN amplification and progress to more aggressive and therapy-resistant phenotypes through accumulation of secondary genomic changes. Progression-related changes include recurrent somatic copy number alterations and typically non-recurrent nucleotide variants, including synonymous and non-coding variants, whose significance has been unclear. To determine if nucleotide variants recurrently affect specific biological processes, we identified altered genes and over-represented variant gene ontologies in 168 exome or whole-genome-sequenced retinoblastomas and 12 tumor-matched cell lines. In addition to RB1 mutations, MYCN amplification, and established retinoblastoma somatic copy number alterations, the analyses revealed enrichment of variant genes related to diverse biological processes including histone monoubiquitination, mRNA processing (P) body assembly, and mitotic sister chromatid segregation and cytokinesis. Importantly, non-coding and synonymous variants increased the enrichment significance of each over-represented biological process term. To assess the effects of such mutations, we examined the consequences of a 3' UTR variant of PCGF3 (a BCOR-binding component of Polycomb repressive complex I), dual 3' UTR variants of CDC14B (a regulator of sister chromatid segregation), and a synonymous variant of DYNC1H1 (a regulator of P-body assembly). One PCGF3 and one of two CDC14B 3' UTR variants impaired gene expression whereas a base-edited DYNC1H1 synonymous variant altered protease sensitivity and stability. Retinoblastoma cell lines retained only ~50% of variants detected in tumors and enriched for new variants affecting p53 signaling. These findings reveal potentially important differences in retinoblastoma cell lines and tumors and implicate synonymous and non-coding variants, along with non-synonymous variants, in retinoblastoma oncogenesis.