mRNA Codon Research Articles

Reprogramming the genetic code with flexizymes.

In the canonical genetic code, the 61 sense codons are assigned to the 20 proteinogenic amino acids. Advancements in genetic code manipulation techniques have enabled the ribosomal incorporation of nonproteinogenic amino acids (npAAs). The critical molecule for translating messenger RNA (mRNA) into peptide sequences is aminoacyl-transfer RNA (tRNA), which recognizes the mRNA codon through its anticodon. Because aminoacyl-tRNA synthetases (ARSs) are highly specific for their respective amino acid-tRNA pairs, it is not feasible to use natural ARSs to prepare npAA-tRNAs. However, flexizymes are adaptable aminoacylation ribozymes that can be used to prepare diverse aminoacyl-tRNAs at will using amino acids activated with suitable leaving groups. Regarding recognition elements, flexizymes require only an aromatic ring in either the leaving group or side chain of the activated amino acid, and the conserved 3'-end CCA of the tRNA. Therefore, flexizymes allow virtually any amino acid to be charged onto any tRNA. The flexizyme system can handle not only L-α-amino acids with side chain modifications but also various backbone-modified npAAs. This Review describes the development of flexizyme variants and discusses their structure and mechanism and their applications in genetic code reprogramming for the synthesis of unique peptides and proteins.

HERV-W ENV transcription in B cells predicting symptomatic COVID-19 and risk for long COVID can express a full-length protein despite stop codon in mRNA from chromosome X via a ribosome readthrough.

The human genome comprises 8% of endogenous retroviruses (HERVs). Though HERVS contribute to physiological functions, copies retained pathogenic potential. The HERV-W ENV protein was shown expressed in patients with worse COVID-19 symptoms and post-COVID syndrome. A significant detection of the mRNA encoding HERV-W ENV from patients with COVID-19 in B cells from RNAseq reads obtained from peripheral blond mononuclear cells. This data stratified with increased COVID-19 symptoms or with post-acute sequelae of COVID-19 (long COVID) after 3 months. The HERV-W ENV-U3R RNA was confirmed to display the best alignment with chromosome X ERVWE2 locus. However, a stop codon precluding its translation was re-addressed after recent understandings of ribosome readthrough mechanisms. Experimental results evidenced that this HERV gene can effectively express a full-length protein in the presence of molecules allowing translation via a readthrough mechanism at the ribosome level. Results not only confirm HERV-W ENV RNA origin in these patients but show for the first time how a defective HERV copy can be translated into a complete protein when specific factors make it possible at the ribosome level. The present proof of concept now requires further studies to identify the factors involved in this newly understood mechanism, following SARS-CoV-2 exposure.

N1-Methylpseudouridine and pseudouridine modifications modulate mRNA decoding during translation

The ribosome utilizes hydrogen bonding between mRNA codons and aminoacyl-tRNAs to ensure rapid and accurate protein production. Chemical modification of mRNA nucleobases can adjust the strength and pattern of this hydrogen bonding to alter protein synthesis. We investigate how the N1-methylpseudouridine (m1Ψ) modification, commonly incorporated into therapeutic and vaccine mRNA sequences, influences the speed and fidelity of translation. We find that m1Ψ does not substantially change the rate constants for amino acid addition by cognate tRNAs or termination by release factors. However, we also find that m1Ψ can subtly modulate the fidelity of amino acid incorporation in a codon-position and tRNA dependent manner in vitro and in human cells. Our computational modeling shows that altered energetics of mRNA:tRNA interactions largely account for the context dependence of the low levels of miscoding we observe on Ψ and m1Ψ containing codons. The outcome of translation on modified mRNA bases is thus governed by the sequence context in which they occur.

Quantitative profiling of m6A at single base resolution across the life cycle of rice and Arabidopsis

N6-methyladenosine (m6A) plays critical roles in regulating mRNA metabolism. However, comprehensive m6A methylomes in different plant tissues with single-base precision have yet to be reported. Here, we present transcriptome-wide m6A maps at single-base resolution in different tissues of rice and Arabidopsis using m6A-SAC-seq. Our analysis uncovers a total of 205,691 m6A sites distributed across 22,574 genes in rice, and 188,282 m6A sites across 19,984 genes in Arabidopsis. The evolutionarily conserved m6A sites in rice and Arabidopsis ortholog gene pairs are involved in controlling tissue development, photosynthesis and stress response. We observe an overall mRNA stabilization effect by 3’ UTR m6A sites in certain plant tissues. Like in mammals, a positive correlation between the m6A level and the length of internal exons is also observed in plant mRNA, except for the last exon. Our data suggest an active m6A deposition process occurring near the stop codon in plant mRNA. In addition, the MTA-installed plant mRNA m6A sites correlate with both translation promotion and translation suppression, depicting a more complicated regulatory picture. Our results therefore provide in-depth resources for relating single-base resolution m6A sites with functions in plants and uncover a suppression-activation model controlling m6A biogenesis across species.

GNN Codon Adjacency Tunes Protein Translation.

The central dogma treats the ribosome as a molecular machine that reads one mRNA codon at a time as it adds each amino acid to its growing peptide chain. However, this and previous studies suggest that ribosomes actually perceive pairs of adjacent codons as they take three-nucleotide steps along the mRNA. We examined GNN codons, which we find are surprisingly overrepresented in eukaryote protein-coding open reading frames (ORFs), especially immediately after NNU codons. Ribosome profiling experiments in yeast revealed that ribosomes with NNU at their aminoacyl (A) site have particularly elevated densities when NNU is immediately followed (3') by a GNN codon, indicating slower mRNA threading of the NNU codon from the ribosome's A to peptidyl (P) sites. Moreover, if the assessment was limited to ribosomes that have only recently arrived at the next codon, by examining 21-nucleotide ribosome footprints (21-nt RFPs), elevated densities were observed for multiple codon classes when followed by GNN. This striking translation slowdown at adjacent 5'-NNN GNN codon pairs is likely mediated, in part, by the ribosome's CAR surface, which acts as an extension of the A-site tRNA anticodon during ribosome translocation and interacts through hydrogen bonding and pi stacking with the GNN codon. The functional consequences of 5'-NNN GNN codon adjacency are expected to influence the evolution of protein coding sequences.

Investigation of Phe-tRNA interaction with EF-Tu in GDP/GTP nucleotide bound states: a molecular dynamics simulation study

Elongation factor Tu (EF-Tu) is an important class of translational GTPases (trGTPases) involved in the elongation process. EF-Tu functions by forming the tight complex with the aminoacyl-tRNA (aa-tRNA: transfer RNA carrying an amino acid) in the GTP-bound active state and subsequently transports the aa-tRNA to the A-site (aminoacyl-site) of the ribosome. The correct interaction between mRNA's codon and aa-tRNA anticodon triggers GTP hydrolysis, resulting in the transition of EF-Tu from the active state to the inactive state (EF-Tu bound to GDP). This, in turn, results in the dissociation of EF-Tu:GDP from the tRNA and ribosome. However, the conformational changes and interactions responsible for the tight complex formation between EF-Tu and aa-tRNA in the GTP-bound state and the disruption of the interaction between EF-Tu and aa-tRNA in the GDP-bound state are not well understood. Therefore, to explore the conformational changes in the EF-Tu:Phe-tRNA complex and to get insight into the interaction between Phe-tRNA (tRNA carrying Phenylalanine) and EF-Tu in GDP/GTP nucleotide bound state, 200 ns MD (molecular dynamics) simulation have been carried out. The RMSD, RMSF, cluster, and DSSP analyses suggest that the GTP-bound state attains a more favorable conformation, which may facilitate EF-Tu’s interactions with the tRNA and ribosome. Further investigation of the non-bonded interaction energy calculation revealed the significance of domainII in interaction with the Phe-tRNA in a GTP-bound state. In addition, the H-bonds calculated between the tRNA’s Phe (Phenylalanine attached to the tRNA is considered) and domainII highlights the contribution of Val285 in recognizing tRNA’s Phe by forming an H-bond throughout the simulation time in the GTP bound state. Overall these results provide insight into how GTP nucleotide influences the interaction between EF-Tu and Phe-tRNA.

Information, Coding, and Biological Function: The Dynamics of Life.

In the mid-20th century, two new scientific disciplines emerged forcefully: molecular biology and information-communication theory. At the beginning, cross-fertilization was so deep that the term genetic code was universally accepted for describing the meaning of triplets of mRNA (codons) as amino acids. However, today, such synergy has not taken advantage of the vertiginous advances in the two disciplines and presents more challenges than answers. These challenges not only are of great theoretical relevance but also represent unavoidable milestones for next-generation biology: from personalized genetic therapy and diagnosis to Artificial Life to the production of biologically active proteins. Moreover, the matter is intimately connected to a paradigm shift needed in theoretical biology, pioneered a long time ago, that requires combined contributions from disciplines well beyond the biological realm. The use of information as a conceptual metaphor needs to be turned into quantitative and predictive models that can be tested empirically and integrated in a unified view. Successfully achieving these tasks requires a wide multidisciplinary approach, including Artificial Life researchers, to address such an endeavour.

Tuning gene expression by hairpin elements near the start codon of mRNA in mammalian cells

The hairpin element (HpE) near the start codon in the 5′ untranslated region (UTR) was developed to tune the mRNA translation in mammalian cells. The parameters of HpEs including thermodynamic stability, the content of guanine and cytosine, and distance between HpEs and the 5′ cap were investigated. These parameters influenced the mRNA expression level. In addition, the start codon and the upstream open reading frame sequestered within the structures of HpEs also reduced the translation initiation. In summary, this study shows that the simple engineering HpE structure can be efficiently adopted for gene expression regulation. The predictable controllability of this simple cloning strategy can potentially achieve precise gene expression regulation in different mammalian cell types.

Preclinical Assessment for Gene Addition Strategy in Familial Hemophagocytic Lymphohistiocytosis Related to Munc 13-4 Deficiency

Introduction: Familial Hemophagocytic Lymphohistiocytosis (FHL) comprises a group of rare genetic disorders with defect cytotoxicity in T and NK cells predisposing to severe hyperinflammation with the clinical hallmark of the Hemophagocytic LymphoHistiocytosis (HLH) syndrome. The uncontrolled proliferation of activated T-CD8 lymphocytes and macrophages infiltrates lymphoid tissues, bone marrow, and multiple organs including the brain. Thus, if untreated full-scale FHL is usually fatal. FHL3 is due to mutations in the UNC13D gene coding for Munc 13-4 protein, one of the components of the perforin-dependent cytotoxicity apparatus. This condition accounts for 30% to 35 % of all genetic cases of HLH. Hematopoietic stem and progenitor cell (HSPC) transplantation, which is the only curative treatment for FHL3 to date, remains difficult even when a compatible donor is available because of the patients' critical inflammatory background. In this context, gene therapy could be a promising therapeutic option, especially for those patients without HLA-compatible donor. Method: The feasibility and efficacy of gene addition strategy for FHL3 have been previously demonstrated in mouse models. We then translated this approach to human cells using a third-generation self-inactivated lentivirus vector pseudo-typed with a conventional vesicular stomatitis virus-G (VSVG) expressing a codon-optimized UNC13D gene under the transcriptional control of EF1 alpha (LV-EF1a-UNC13D). HSPC or CD3+ T-cells were transduced with a GMP-like-grade LV-EF1a-UNC13D lentiviral vector. We studied the restoration of cytotoxic function. Results: Our GMP transduction protocol led to efficient transcription of optimized mRNA codon and expression of Munc 13-4 protein as evidenced by RT qPCR and western blot, respectively. Interestingly, we show restored cytotoxicity activity after CD34+ differentiation toward CD3+ T-cells and NK cells. The vector copy number is compatible with a clinical approach. In vivo experiments are ongoing to validate the absence of genotoxicity in this setting. Discussion conclusion: Our results support the feasibility and efficacy of UNCD13D gene transfer into hematopoietic stem and progenitor cells or T cells from Munc 13-4 deficient patients. Alternatively, gene therapy of HSPCs could be combined with transduced autologous T-cells infusion. The first attempts regarding this option are encouraging. UNC13D gene addition strategy could dramatically improve the care of FHL3 patients lacking an HLA-genoidentical donor.

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding

Accurate protein synthesis is determined by the two-subunit ribosome’s capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site’s physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.

Modulation of translational decoding by m6A modification of mRNA

N6-methyladenosine (m6A) is an abundant, dynamic mRNA modification that regulates key steps of cellular mRNA metabolism. m6A in the mRNA coding regions inhibits translation elongation. Here, we show how m6A modulates decoding in the bacterial translation system using a combination of rapid kinetics, smFRET and single-particle cryo-EM. We show that, while the modification does not impair the initial binding of aminoacyl-tRNA to the ribosome, in the presence of m6A fewer ribosomes complete the decoding process due to the lower stability of the complexes and enhanced tRNA drop-off. The mRNA codon adopts a π-stacked codon conformation that is remodeled upon aminoacyl-tRNA binding. m6A does not exclude canonical codon-anticodon geometry, but favors alternative more dynamic conformations that are rejected by the ribosome. These results highlight how modifications outside the Watson-Crick edge can still interfere with codon-anticodon base pairing and complex recognition by the ribosome, thereby modulating the translational efficiency of modified mRNAs.

Direct Nanopore Sequencing for the 17 RNA Modification Types in 36 Locations in the E. coli Ribosome Enables Monitoring of Stress-Dependent Changes.

The bacterium Escherichia coli possesses 16S and 23S rRNA strands that have 36 chemical modification sites with 17 different structures. Nanopore direct RNA sequencing using a protein nanopore sensor and helicase brake, which is also a sensor, was applied to the rRNAs. Nanopore current levels, base calling profile, and helicase dwell times for the modifications relative to unmodified synthetic rRNA controls found signatures for nearly all modifications. Signatures for clustered modifications were determined by selective sequencing of writer knockout E. coli and sequencing of synthetic RNAs utilizing some custom-synthesized nucleotide triphosphates for their preparation. The knowledge of each modification's signature, apart from 5-methylcytidine, was used to determine how metabolic and cold-shock stress impact rRNA modifications. Metabolic stress resulted in either no change or a decrease, and one site increased in modification occupancy, while cold-shock stress led to either no change or a decrease. The double modification m4Cm1402 resides in 16S rRNA, and it decreased with both stressors. Using the helicase dwell time, it was determined that the N4 methyl group is lost during both stressors, and the 2'-OMe group remained. In the ribosome, this modification stabilizes binding to the mRNA codon at the P-site resulting in increased translational fidelity that is lost during stress. The E. coli genome has seven rRNA operons (rrn), and the earlier studies aligned the nanopore reads to a single operon (rrnA). Here, the reads were aligned to all seven operons to identify operon-specific changes in the 11 pseudouridines. This study demonstrates that direct sequencing for >16 different RNA modifications in a strand is achievable.

Analysing the genetic code degeneracy: a consequence towards bacterial staining

As 20 naturally occurring amino acids are coded by 61 mRNA codons out of 64, it is obvious that 61→20 cannot have one-to-one mapping which generates the problem of codon degeneracy. Despite several efforts there is no specific outcome which can describe this well-known enigmatic degeneracy of the codon table. Since, every biological behaviour is regulated by protein which in turn consists of amino acids bearing the inherent characteristics of degeneracy among mRNA codons (Crick F.H.C. The Origin of the Genetic Code. J. Mol. Biol.1968; 38: 367-379), it is worthy to analyse the impact of such degeneracy on biological behaviours. Here, based on mathematical models using the concept of b-type of the nucleotide bases and hamming distances, an effort has been initiated to understand the impact of biasness of genetic code degeneracy on biological behaviours. The proposed models have been utilized to understand the characteristic features of bacterial genes of gram-positive and gram-negative bacteria. To the best of our knowledge, this is the first mathematical model to capture the effect of genetic code degeneracy, showing a paradigm towards understanding the behavioural difference between gram-positive and gram-negative bacteria, and thereby opening a new avenue for revealing differential biological properties. Communicated by Ramaswamy H. Sarma

The Fitness Effects of Codon Composition of the Horizontally Transferred Antibiotic Resistance Genes Intensify at Sub-lethal Antibiotic Levels.

The rampant variability in codon bias existing between bacterial genomes is expected to interfere with horizontal gene transfer (HGT), a phenomenon that drives bacterial adaptation. However, delineating the constraints imposed by codon bias on functional integration of the transferred genes is complicated by multiple genomic and functional barriers controlling HGT, and by the dependence of the evolutionary outcomes of HGT on the host's environment. Here, we designed an experimental system in which codon composition of the transferred genes is the only variable triggering fitness change of the host. We replaced Escherichia coli's chromosomal folA gene encoding dihydrofolate reductase, an essential enzyme that constitutes a target for trimethoprim, with combinatorial libraries of synonymous codons of folA genes from trimethoprim-sensitive Listeria grayi and trimethoprim-resistant Neisseria sicca. The resulting populations underwent selection at a range of trimethoprim concentrations, and the ensuing changes in variant frequencies were used to infer the fitness effects of the individual combinations of codons. We found that when HGT causes overstabilization of the 5'-end mRNA, the fitness contribution of mRNA folding stability dominates over that of codon optimality. The 5'-end overstabilization can also lead to mRNA accumulation outside of the polysome, thus preventing the decay of the foreign transcripts despite the codon composition-driven reduction in translation efficiency. Importantly, the fitness effects of mRNA stability or codon optimality become apparent only at sub-lethal levels of trimethoprim individually tailored for each library, emphasizing the central role of the host's environment in shaping the codon bias compatibility of horizontally transferred genes.

RPL3L-containing ribosomes determine translation elongation dynamics required for cardiac function

Although several ribosomal protein paralogs are expressed in a tissue-specific manner, how these proteins affect translation and why they are required only in certain tissues have remained unclear. Here we show that RPL3L, a paralog of RPL3 specifically expressed in heart and skeletal muscle, influences translation elongation dynamics. Deficiency of RPL3L-containing ribosomes in RPL3L knockout male mice resulted in impaired cardiac contractility. Ribosome occupancy at mRNA codons was found to be altered in the RPL3L-deficient heart, and the changes were negatively correlated with those observed in myoblasts overexpressing RPL3L. RPL3L-containing ribosomes were less prone to collisions compared with RPL3-containing canonical ribosomes. Although the loss of RPL3L-containing ribosomes altered translation elongation dynamics for the entire transcriptome, its effects were most pronounced for transcripts related to cardiac muscle contraction and dilated cardiomyopathy, with the abundance of the encoded proteins being correspondingly decreased. Our results provide further insight into the mechanisms and physiological relevance of tissue-specific translational regulation.

MRNA decoding in human is kinetically and structurally distinct from bacteria

In all species, ribosomes synthesize proteins by faithfully decoding messenger RNA (mRNA) nucleotide sequences using aminoacyl-tRNA substrates. Current knowledge of the decoding mechanism derives principally from studies on bacterial systems1. Although key features are conserved across evolution2, eukaryotes achieve higher-fidelity mRNA decoding than bacteria3. In human, changes in decoding fidelity are linked to ageing and disease and represent a potential point of therapeutic intervention in both viral and cancer treatment4–6. Here we combine single-molecule imaging and cryogenic electron microscopy methods to examine the molecular basis of human ribosome fidelity to reveal that the decoding mechanism is both kinetically and structurally distinct from that of bacteria. Although decoding is globally analogous in both species, the reaction coordinate of aminoacyl-tRNA movement is altered on the human ribosome and the process is an order of magnitude slower. These distinctions arise from eukaryote-specific structural elements in the human ribosome and in the elongation factor eukaryotic elongation factor 1A (eEF1A) that together coordinate faithful tRNA incorporation at each mRNA codon. The distinct nature and timing of conformational changes within the ribosome and eEF1A rationalize how increased decoding fidelity is achieved and potentially regulated in eukaryotic species.

The power of 2,6-diaminopurine in correcting UGA nonsense codons in CFTR mRNA

The power of 2,6-diaminopurine in correcting UGA nonsense codons in CFTR mRNA

Structural basis for reduced ribosomal A-site fidelity in response to P-site codon–anticodon mismatches

Rapid and accurate translation is essential in all organisms to produce properly folded and functional proteins. mRNA codons that define the protein-coding sequences are decoded by tRNAs on the ribosome in the aminoacyl (A) binding site. The mRNA codon and the tRNA anticodon interaction is extensively monitored by the ribosome to ensure accuracy in tRNA selection. While other polymerases that synthesize DNA and RNA can correct for misincorporations, the ribosome is unable to correct mistakes. Instead, when a misincorporation occurs, the mismatched tRNA-mRNA pair moves to the peptidyl (P) site and, from this location, causes a reduction in the fidelity at the A site, triggering post-peptidyl transfer quality control. This reduced fidelity allows for additional incorrect tRNAs to be accepted and for release factor 2 (RF2) to recognize sense codons, leading to hydrolysis of the aberrant peptide. Here, we present crystal structures of the ribosome containing a tRNALys in the P site with a U•U mismatch with the mRNA codon. We find that when the mismatch occurs in the second position of the P-site codon-anticodon interaction, the first nucleotide of the A-site codon flips from the mRNA path to engage highly conserved 16S rRNA nucleotide A1493 in the decoding center. We propose that this mRNA nucleotide mispositioning leads to reduced fidelity at the A site. Further, this state may provide an opportunity for RF2 to initiate premature termination before erroneous nascent chains disrupt the cellular proteome.

Drop-off-reinitiation at the amino termini of nascent peptides and its regulation by IF3, EF-G, and RRF.

In translation initiation in prokaryotes, IF3 recognizes the interaction between the initiator codon of mRNA and the anticodon of fMet-tRNAini and then relocates the fMet-tRNAini to an active position. Here, we have surveyed 328 codon-anticodon combinations for the preference of IF3. At the first and second base of the codon, only Watson-Crick base pairs are tolerated. At the third base, stronger base pairs, for example, Watson-Crick, are more preferred, but other types of base pairs, for example, G/U wobble, are also tolerated; weaker base pairs are excluded by IF3. When the codon-anticodon combinations are unfavorable for IF3 or the concentration of IF3 is too low to recognize any codon-anticodon combinations, IF3 fails to set the P-site fMet-tRNAini at the active position and causes its drop-off from the ribosome. Thereby, translation reinitiation occurs from the second aminoacyl-tRNA at the A site to yield a truncated peptide lacking the amino-terminal fMet. We refer to this event as the amino-terminal drop-off-reinitiation. We also showed that EF-G and RRF are involved in disassembling such an aberrant ribosome complex bearing inactive fMet-tRNAini Thereby EF-G and RRF are able to exclude unfavorable codon-anticodon combinations with weaker base pairs and alleviate the amino-terminal drop-off-reinitiation.

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding.

Accurate protein synthesis is determined by the two-subunit ribosome's capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular and structural basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site's physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.