Influence of the Number, Nature and Position of Methyl Posttranscriptional Modifications on Nucleobase Stacking in RNA

Exploring the Impact of mRNA Modifications on Translation Efficiency and Immune Tolerance to Self-Antigens

Therapeutic modified mRNAs are being developed for a broad range of human diseases. However, the impact of potential miscoding of modified mRNAs on self-tolerance remains unknown. Additionally, more studies are needed to explore the effects of nucleoside alkylation on translation. While all six tested modifications are tolerated as substrates by T7 RNA polymerase and inhibited mRNA immunogenicity, the translation efficiency varied significantly depending on the type of modification. In contrast to methylation, ethylation at the N1 position of pseudouridine (Ψ) hindered translation, suggesting that the C5-C1' glycosidic bond alone is not a critical element for high translation. Inhibition of mRNA translation was also observed with 5-methoxyuridine modification. However, this inhibition was partially alleviated through the optimization of mRNA coding sequences. BALB/c mice immunized with syngeneic ψ-modified mRNA encoding for Wilms' tumor antigen-1 (WT1) developed a low but significant level of anti-WT1 IgG antibodies compared to those immunized with either unmodified or N1-methyl ψ-modified mRNA. Overall, the data indicate that adding a simple ethyl group (-CH2CH3) at the N1 position of ψ has a major negative effect on translation despite its reduced immunogenicity. Additionally, mRNA containing Ψ may alter translation fidelity at certain codons, which could lead to a breakdown of immune tolerance to self-antigens. This concern should be taken into account during gene replacement therapies, although it could benefit mRNA-based vaccines by generating a diverse repertoire of antigens.

An Atlas of the base inter-RNA stacks involved in bacterial translation

Nucleobase-specific noncovalent interactions play a crucial role in translation. Herein, we provide a comprehensive analysis of the stacks between different RNA components in the crystal structures of the bacterial ribosome caught at different translation stages. Analysis of tRNA||rRNA stacks reveals distinct behaviour; both the A-and E-site tRNAs exhibit unique stacking patterns with 23S rRNA bases, while P-site tRNAs stack with 16S rRNA bases. Furthermore, E-site stacks exhibit diverse face orientations and ring topologies-rare for inter-chain RNA interactions-with higher average interaction energies than A or P-site stacks. This suggests that stacking may be essential for stabilizing tRNA progression through the E-site. Additionally, mRNA||rRNA stacks reveal other geometries, which depend on the tRNA binding site, whereas 16S rRNA||23S rRNA stacks highlight the importance of specific bases in maintaining the integrity of the translational complex by linking the two rRNAs. Furthermore, tRNA||mRNA stacks exhibit distinct geometries and energetics at the E-site, indicating their significance during tRNA translocation and elimination. Overall, both A and E-sites display a more diverse distribution of inter-RNA stacks compared to the P-site. Stacking interactions in the active ribosome are not simply accidental byproducts of biochemistry but are likely invoked to compensate and support the integrity and dynamics of translation.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

Structural Properties of Hachimoji Nucleic Acids and Their Building Blocks: Comparison of Genetic Systems with Four, Six and Eight Alphabets

Expansion of the genetic alphabet is an ambitious goal. A recent breakthrough has led to the eight-base (hachimoji) genetics having canonical and unnatural bases. However, very little is known on the molecular-level features that facilitate the candidature of unnatural bases as genetic alphabets. Here we amalgamated DFT calculations and MD simulations to analyse the properties of the constituents of hachimoji DNA and RNA. DFT reveals the dominant syn conformation for isolated unnatural deoxyribonucleosides and at the 5'-end of oligonucleotides, although an anti/syn mixture is predicted at the nonterminal and 3'-terminal positions. However, isolated ribonucleotides prefer an anti/syn mixture, but mostly prefer anti conformation at the nonterminal positions. Further, the canonical base pairing combinations reveals significant strength, which may facilitate replication of hachimoji DNA. We also identify noncanonical base pairs that can better tolerate the substitution of unnatural pairs in RNA. Stacking strengths of 51 dimers reveals higher average stacking stabilization of dimers of hachimoji bases than canonical bases, which provides clues for choosing energetically stable sequences. A total of 14.4 μs MD simulations reveal the influence of solvent on the properties of hachimoji oligonucleotides and point to the likely fidelity of replication of hachimoji DNA. Our results pinpoint the features that explain the experimentally observed stability of hachimoji DNA.

Role of Main RNA Methylation in Hepatocellular Carcinoma: N6-Methyladenosine, 5-Methylcytosine, and N1-Methyladenosine

RNA methylation is considered a significant epigenetic modification, a process that does not alter gene sequence but may play a necessary role in multiple biological processes, such as gene expression, genome editing, and cellular differentiation. With advances in RNA detection, various forms of RNA methylation can be found, including N6-methyladenosine (m6A), N1-methyladenosine (m1A), and 5-methylcytosine (m5C). Emerging reports confirm that dysregulation of RNA methylation gives rise to a variety of human diseases, particularly hepatocellular carcinoma. We will summarize essential regulators of RNA methylation and biological functions of these modifications in coding and noncoding RNAs. In conclusion, we highlight complex molecular mechanisms of m6A, m5C, and m1A associated with hepatocellular carcinoma and hope this review might provide therapeutic potent of RNA methylation to clinical research.

Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands

Guanidinium ion is a toxic cellular metabolite. The ykkC-III riboswitch, an mRNA stretch, regulates the gene expression by undergoing a conformational change in response to the binding of a free guanidinium ion and thereby plays a potentially important role in alleviating guanidinium toxicity in cells. An experimental crystal structure of the guanidinium-bound aptamer domain of the riboswitch from Thermobifida Fusca revealed the overall RNA architecture and mapped the specific noncovalent interactions that stabilize the ligand within the binding pocket aptamer. However, details of how the aptamer domain discriminates the cognate ligand from its closest structurally analogous physiological metabolites (arginine and urea), and how the binding of cognate ligand arrays information from the aptamer domain to the expression platform for regulating the gene expression, are not well understood. To fill this void, we perform a cumulative of 2 μs all-atom explicit-solvent molecular dynamics (MD) simulations on the full aptamer domain, augmented with quantum-chemical calculations on the ligand-binding pocket, to compare the structural and dynamical details of the guanidinium-bound state with the arginine or urea bound states, as well as the unbound (open) state. Analysis of the ligand-binding pocket reveals that due to unfavorable interactions with the binding-pocket residues, urea cannot bind the aptamer domain and thereby cannot alter the gene expression. Although interaction of the guanidyl moiety of arginine within the binding pocket is either comparable or stronger than the guanidinium ion, additional non-native hydrogen-bonding networks, as well as differences in the dynamical details of the arginine-bound state, explain why arginine cannot transmit the information from the aptamer domain to the expression platform. Based on our simulations, we propose a mechanism of how the aptamer domain communicates with the expression platform. Overall, our work provides interesting insights into the ligand recognition by a specific class of riboswitches and may hopefully inspire future studies to further understand the gene regulation by riboswitches.