Isostericity and tautomerism of base pairs in nucleic acids

The Expanding Universe of Extensions and Tails: Ribosomal Proteins and Histones in RNA and DNA Complex Signaling and Dynamics

This short review bridges two biological fields: ribosomes and nucleosomes-two nucleoprotein assemblies that, along with many viruses, share proteins featuring long filamentous segments at their N- or C-termini. A central hypothesis is that these extensions and tails perform analogous functions in both systems. The evolution of these structures appears closely tied to the emergence of regulatory networks and signaling pathways, facilitating increasingly complex roles for ribosomes and nucleosome alike. This review begins by summarizing the structures and functions of ribosomes and nucleosomes, followed by a detailed comparison highlighting their similarities and differences, particularly in light of recent findings on the roles of ribosomal proteins in signaling and ribosome dynamics. The analysis seeks to uncover whether these systems operate based on shared principles and mechanisms. The nucleosome-ribosome analogy may offer valuable insights into unresolved questions in both fields. For instance, new structural insights from ribosomes might shed light on potential motifs formed by histone tails. From an evolutionary perspective, this study revisits the origins of signaling and regulation in ancient nucleoprotein assemblies, suggesting that tails and extensions may represent remnants of the earliest network systems governing signaling and dynamic control.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

yaaJ, the tRNA-Specific Adenosine Deaminase, Is Dispensable in Bacillus subtilis

Post-transcriptional modifications of tRNA are crucial for their core function. The inosine (I; 6-deaminated adenosine) at the first position in the anticodon of tRNAArg(ICG) modulates the decoding capability and is generally considered essential for reading CGU, CGC, and CGA codons in eubacteria. We report here that the Bacillus subtilis yaaJ gene encodes tRNA-specific adenosine deaminase and is non-essential for viability. A β-galactosidase reporter assay revealed that the translational activity of CGN codons was not impaired in the yaaJ-deletion mutant. Furthermore, tRNAArg(CCG) responsible for decoding the CGG codon was dispensable, even in the presence or absence of yaaJ. These results strongly suggest that tRNAArg with either the anticodon ICG or ACG has an intrinsic ability to recognize all four CGN codons, providing a fundamental concept of non-canonical wobbling mediated by adenosine and inosine nucleotides in the anticodon. This is the first example of the four-way wobbling by inosine nucleotide in bacterial cells. On the other hand, the absence of inosine modification induced +1 frameshifting, especially at the CGA codon. Additionally, the yaaJ deletion affected growth and competency. Therefore, the inosine modification is beneficial for translational fidelity and proper growth-phase control, and that is why yaaJ has been actually conserved in B. subtilis.

Anionic G•U pairs in bacterial ribosomal rRNAs

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Conformation-based refinement of 18-mer DNA structures

Nine new crystal structures of CG-rich DNA 18-mers with the sequence 5'-GGTGGGGGC-XZ-GCCCCACC-3', which are related to the bacterial repetitive extragenic palindromes, are reported. 18-mer oligonucleotides with the central XZ dinucleotide systematically mutated to all 16 sequences show complex behavior in solution, but all ten so far successfully crystallized 18-mers crystallized as A-form duplexes. The refinement protocol benefited from the recurrent use of geometries of the dinucleotide conformer (NtC) classes as refinement restraints in regions of poor electron density. The restraints are automatically generated at the dnatco.datmos.org web service and are available for download. This NtC-driven protocol significantly helped to stabilize the structure refinement. The NtC-driven refinement protocol can be adapted to other low-resolution data such as cryo-EM maps. To test the quality of the final structural models, a novel validation method based on comparison of the electron density and conformational similarity to the NtC classes was employed.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

U5 snRNA Interactions With Exons Ensure Splicing Precision

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5-the exon junction positioned at U5Loop1 C₃₉|C₃₈; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition-guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson-Crick pairs with the 5'exon in the absence of +5G at the intron start. In 5'exon positions -3 and -5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson-Crick pairs with 3'exon position +1 coincide with substitutions of the conserved -3C at the intron 3'end. Based on mutation and X-ray evidence, we propose that -3C pairs with U2 G₃₁ juxtaposing the branchpoint and the 3'intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3'intron end ensure that the 3'exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson-Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson-Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Can modified DNA base pairs with chalcogen bonding expand the genetic alphabet? A combined quantum chemical and molecular dynamics simulation study

A comprehensive (DFT and MD) computational study is presented with the goal to design and analyze model chalcogen-bonded modified nucleobase pairs that replace one (i.e., A_XY:T, G:C_XY, G_XY:C) or two (G_XY:C_X'Y', X/X' = S, Se and Y/Y' = F, Cl, Br) Watson-Crick (WC) hydrogen bonds of the canonical A:T or G:C pair with chalcogen bond(s). DFT calculations on 18 base pair combinations that replace one WC hydrogen bond with a chalcogen bond reveal that the bases favorably interact in the gas phase (binding strengths up to -140 kJ mol^-1) and water (up to -85 kJ mol^-1). Although the remaining hydrogen bond(s) exhibits similar characteristics to those in the canonical base pairs, the structural features of the (Y-XO) chalcogen bond(s) change significantly with the identity of X and Y. The 36 doubly-substituted (G_XY:C_X'Y') base pairs have structural deviations from canonical G:C similar to those of the singly-substituted modifications (G:C_XY or G_XY:C). Furthermore, despite the replacement of two strong hydrogen bonds with chalcogen bonds, some G_XY:C_X'Y' pairs possess comparable binding energies (up to -132 kJ mol^-1 in the gas phase and up to -92 kJ mol^-1 in water) to the most stable G:C_XY or G_XY:C pairs, as well as canonical G:C. More importantly, G:C-modified pairs containing X = Se (high polarizability) and Y = F (high electronegativity) are the most stable, with comparable or slightly larger (by up to 13 kJ mol^-1) binding energies than G:C. Further characterization of the chalcogen bonding in all modified base pairs (AIM, NBO and NCI analyses) reveals that the differences in the binding energies of modified base pairs are mainly dictated by the differences in the strengths of their chalcogen bonds. Finally, MD simulations on DNA oligonucleotides containing the most stable chalcogen-bonded base pair from each of the four classifications (A_XY:T, G:C_XY, G_XY:C and G_XY:C_X'Y') reveal that the singly-modified G:C pairs best retain the local helical structure and pairing stability to a greater extent than the modified A:T pair. Overall, our study identifies two (G:C_SeF and G_SeF:C) promising pairs that retain chalcogen bonding in DNA and should be synthesized and further explored in terms of their potential to expand the genetic alphabet.

Probing conformational transitions towards mutagenic Watson-Crick-like G·T mismatches using off-resonance sugar carbon R1ρ relaxation dispersion

NMR off-resonance R_1ρ relaxation dispersion measurements on base carbon and nitrogen nuclei have revealed that wobble G·T/U mismatches in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-abundance, and mutagenic Watson-Crick-like conformations. As Watson-Crick-like G·T mismatches have base pairing geometries similar to Watson-Crick base pairs, we hypothesized that they would mimic Watson-Crick base pairs with respect to the sugar-backbone conformation as well. Using off-resonance R_1ρ measurements targeting the sugar C3' and C4' nuclei, a structure survey, and molecular dynamics simulations, we show that wobble G·T mismatches adopt sugar-backbone conformations that deviate from the canonical Watson-Crick conformation and that transitions toward tautomeric and anionic Watson-Crick-like G·T mismatches restore the canonical Watson-Crick sugar-backbone. These measurements also reveal kinetic isotope effects for tautomerization in D₂O versus H₂O, which provide experimental evidence in support of a transition state involving proton transfer. The results provide additional evidence in support of mutagenic Watson-Crick-like G·T mismatches, help rule out alternative inverted wobble conformations in the case of anionic G·T^-, and also establish sugar carbons as new non-exchangeable probes of this exchange process.

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC⁺ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.

Pseudouridines or how to draw on weak energy differences

In many RNA molecules, pseudouridines occur at conserved positions in functional sites. A great diversity of pseudouridine synthases guarantees the insertion of the modified base at precise locations. The accepted structural role of pseudouridines is a reduction of the RNA flexibility around the modification site. However, experiments rarely yield clear-cut evidence. The article "Dynamic stacking of an expected branch point adenosine in duplexes containing pseudouridine-modified or unmodified U2 snRNA sites" published in 2019 in Biochemical and Biophysical Research Communication by Kennedy et al. constitute a provocative case [1]. This example illustrates how a definite conformational state can be selected through small energy differences in a constrained environment.

Symmetry in Nucleic-Acid Double Helices

In nature and in the test tube, nucleic acids occur in many different forms. Apart from single-stranded, coiled molecules, DNA and RNA prefer to form helical arrangements, in which the bases are stacked to shield their hydrophobic surfaces and expose their polar edges. Focusing on double helices, we describe the crucial role played by symmetry in shaping DNA and RNA structure. The base pairs in nucleic-acid double helices display rotational pseudo-symmetry. In the Watson–Crick base pairs found in naturally occurring DNA and RNA duplexes, the symmetry axis lies in the base-pair plane, giving rise to two different helical grooves. In contrast, anti-Watson–Crick base pairs have a dyad axis perpendicular to the base-pair plane and identical grooves. In combination with the base-pair symmetry, the syn/anti conformation of paired nucleotides determines the parallel or antiparallel strand orientation of double helices. DNA and RNA duplexes in nature are exclusively antiparallel. Watson–Crick base-paired DNA or RNA helices display either right-handed or left-handed helical (pseudo-) symmetry. Genomic DNA is usually in the right-handed B-form, and RNA double helices adopt the right-handed A-conformation. Finally, there is a higher level of helical symmetry in superhelical DNA in which B-form double strands are intertwined in a right- or left-handed sense.

Codon Directional Asymmetry Suggests Swapped Prebiotic 1st and 2nd Codon Positions

Background: Codon directional asymmetry (CDA) classifies the 64 codons into palindromes (XYX, CDA = 0), and 5'- and 3'-dominant (YXX and XXY, CDA < 0 and CDA > 0, respectively). Previously, CDA was defined by the purine/pyrimidine divide (A,G/C,T), where X is either a purine or a pyrimidine. For the remaining codons with undefined CDA, CDA was defined by the 5' or 3' nucleotide complementary to Y. This CDA correlates with cognate amino acid tRNA synthetase classes, antiparallel beta sheet conformation index and the evolutionary order defined by the self-referential genetic code evolution model (CDA < 0: class I, high beta sheet index, late genetic code inclusion). Methods: We explore associations of CDAs defined by nucleotide classifications according to complementarity strengths (A:T, weak; C:G, strong) and keto-enol/amino-imino groupings (G,T/A,C), also after swapping 1st and 2nd codon positions with amino acid physicochemical and structural properties. Results: Here, analyses show that for the eight codons whose purine/pyrimidine-based CDA requires using the rule of complementarity with the midposition, using weak interactions to define CDA instead of complementarity increases associations with tRNA synthetase classes, antiparallel beta sheet index and genetic code evolutionary order. CDA defined by keto-enol/amino-imino groups, 1st and 2nd codon positions swapped, correlates with amino acid parallel beta sheet formation indices and Doolittle's hydropathicities. Conclusions: Results suggest (a) prebiotic swaps from N2N1N3 to N1N2N3 codon structures, (b) that tRNA-mediated translation replaced direct codon-amino acid interactions, and (c) links between codon structures and cognate amino acid properties.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

A tRNA-mimic Strategy to Explore the Role of G34 of tRNAGly in Translation and Codon Frameshifting

Decoding of the 61 sense codons of the genetic code requires a variable number of tRNAs that establish codon-anticodon interactions. Thanks to the wobble base pairing at the third codon position, less than 61 different tRNA isoacceptors are needed to decode the whole set of codons. On the tRNA, a subtle distribution of nucleoside modifications shapes the anticodon loop structure and participates to accurate decoding and reading frame maintenance. Interestingly, although the 61 anticodons should exist in tRNAs, a strict absence of some tRNAs decoders is found in several codon families. For instance, in Eukaryotes, G34-containing tRNAs translating 3-, 4- and 6-codon boxes are absent. This includes tRNA specific for Ala, Arg, Ile, Leu, Pro, Ser, Thr, and Val. tRNA^Gly is the only exception for which in the three kingdoms, a G34-containing tRNA exists to decode C3 and U3-ending codons. To understand why G34-tRNA^Gly exists, we analysed at the genome wide level the codon distribution in codon +1 relative to the four GGN Gly codons. When considering codon GGU, a bias was found towards an unusual high usage of codons starting with a G whatever the amino acid at +1 codon. It is expected that GGU codons are decoded by G34-containing tRNA^Gly, decoding also GGC codons. Translation studies revealed that the presence of a G at the first position of the downstream codon reduces the +1 frameshift by stabilizing the G34•U3 wobble interaction. This result partially explains why G34-containing tRNA^Gly exists in Eukaryotes whereas all the other G34-containing tRNAs for multiple codon boxes are absent.

Unique tRNA gene profile suggests paucity of nucleotide modifications in anticodons of a deep-sea symbiotic Spiroplasma

The position 34 of a tRNA is always modified for efficient recognition of codons and accurate integration of amino acids by the translation machinery. Here, we report genomics features of a deep-sea gut symbiotic Spiroplasma, which suggests that the organism does not require tRNA(34) anticodon modifications. In the genome, there is a novel set of tRNA genes composed of 32 species for recognition of the 20 amino acids. Among the anticodons of the tRNAs, we witnessed the presence of both U34- and C34-containing tRNAs required to decode NNR (A/G) 2:2 codons as countermeasure of probable loss of anticodon modification genes. In the tRNA fragments detected in the gut transcriptome, mismatches expected to be caused by some tRNA modifications were not shown in their alignments with the corresponding genes. However, the probable paucity of modified anticodons did not fundamentally change the codon usage pattern of the Spiroplasma. The tRNA gene profile that probably resulted from the paucity of tRNA(34) modifications was not observed in other symbionts and deep-sea bacteria, indicating that this phenomenon was an evolutionary dead-end. This study provides insights on co-evolution of translation machine and tRNA genes and steric constraints of codon-anticodon interactions in deep-sea extreme environment.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

The multiple flavors of GoU pairs in RNA

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.

NMR analyses on N-hydroxymethylated nucleobases - implications for formaldehyde toxicity and nucleic acid demethylases

Formaldehyde is produced in cells by enzyme-catalysed demethylation reactions, including those occurring on N-methylated nucleic acids. Formaldehyde reacts with nucleobases to form N-hydroxymethylated adducts that may contribute to its toxicity/carcinogenicity when added exogenously, but the chemistry of these reactions has been incompletely defined. We report NMR studies on the reactions of formaldehyde with canonical/modified nucleobases. The results reveal that hydroxymethyl hemiaminals on endocyclic nitrogens, as observed with thymidine and uridine monophosphates, are faster to form than equivalent hemiaminals on exocyclic nitrogens; however, the exocyclic adducts, as formed with adenine, guanine and cytosine, are more stable in solution. Nucleic acid demethylase (FTO)-catalysed hydroxylation of (6-methyl)adenosine results in (6-hydroxymethyl)adenosine as the major observed product; by contrast no evidence for a stable 3-hydroxymethyl adduct was accrued with FTO-catalysed oxidation of (3-methyl)thymidine. Collectively, our results imply N-hydroxymethyled adducts of nucleic acid bases, formed either by reactions with formaldehyde or via demethylase catalysis, have substantially different stabilities, with some being sufficiently stable to have functional roles in disease or the regulation of nucleic acid/nucleobase activity.

Deciphering the reading of the genetic code by near-cognate tRNA

Protein translation is a key cellular process in which each codon of mRNAs has to be accurately and efficiently recognized by cognate tRNAs of a large repertoire of noncognate tRNAs. A successful decoding process is largely dependent on the presence of modified nucleotides within the anticodon loop, especially of tRNAs having to read A/U-rich codons. In this latter case, their roles appear to stabilize the codon–anticodon interaction, allowing them to reach an optimal energetic value close to that of other interacting tRNAs involving G/C-rich anticodons. In this work we demonstrate that, while helping an efficient translation of A/U-rich codons, modified nucleotides also allow certain unconventional base pairing to occur, as evidenced in the case of stop codon suppression.

Some codons of the genetic code can be read not only by cognate, but also by near-cognate tRNAs. This flexibility is thought to be conferred mainly by a mismatch between the third base of the codon and the first of the anticodon (the so-called “wobble” position). However, this simplistic explanation underestimates the importance of nucleotide modifications in the decoding process. Using a system in which only near-cognate tRNAs can decode a specific codon, we investigated the role of six modifications of the anticodon, or adjacent nucleotides, of the tRNAs specific for Tyr, Gln, Lys, Trp, Cys, and Arg in Saccharomyces cerevisiae. Modifications almost systematically rendered these tRNAs able to act as near-cognate tRNAs at stop codons, even though they involve noncanonical base pairs, without markedly affecting their ability to decode cognate or near-cognate sense codons. These findings reveal an important effect of modifications to tRNA decoding with implications for understanding the flexibility of the genetic code.

FT-IR and FT-Raman spectra of 5-chlorocytosine: Solid state simulation and tautomerism. Effect of the chlorine substitution in the Watson-Crick base pair 5-chlorodeoxycytidine-deoxyguanosine

The laser Raman and IR spectra of 5-chlorocytosine have been recorded and accurately assigned in the solid state using Density functional calculations (DFT) together with the linear scaling equation procedure (LSE) and the solid state simulation of the crystal unit cell through a tetramer form. These results remarkably improve those reported previously by other authors. Several new scaling equations were proposed to be used in related molecules. The six main tautomers of the biomolecule 5-chlorocytosine were determined and optimized at the MP2 and CCSD levels, using different basis sets. The relative stabilities were compared with those obtained in cytosine and their 5-halo derivatives. Several relationships between energies, geometric parameters and NBO atomic charges were established. The effect of the chlorine substitution in the fifth position was evaluated through the stability of the Watson-Crick (WC) base pair of 5-chlorodeoxycytidine with deoxyguanosine, and through their vibrational spectra.

Amino acid misincorporation in recombinant proteins

Proteins provide the molecular basis for cellular structure, catalytic activity, signal transduction, and molecular transport in biological systems. Recombinant protein expression is widely used to prepare and manufacture novel proteins that serve as the foundation of many biopharmaceutical products. However, protein translation bioprocesses are inherently prone to low-level errors. These sequence variants caused by amino acid misincorporation have been observed in both native and recombinant proteins. Protein sequence variants impact product quality, and their presence can be exacerbated through cellular stress, overexpression, and nutrient starvation. Therefore, the cell line selection process, which is used in the biopharmaceutical industry, is not only directed towards maximizing productivity, but also focuses on selecting clones which yield low sequence variant levels, thereby proactively avoiding potentially inauspicious patient safety and efficacy outcomes. Here, we summarize a number of hallmark studies aimed at understanding the mechanisms of amino acid misincorporation, as well as exacerbating factors, and mitigation strategies. We also describe key advances in analytical technologies in the identification and quantification of sequence variants, and some practical considerations when using LC-MS/MS for detecting sequence variants.

C5-substituents of uridines and 2-thiouridines present at the wobble position of tRNA determine the formation of their keto-enol or zwitterionic forms - a factor important for accuracy of reading of guanosine at the 3΄-end of the mRNA codons

Modified nucleosides present in the wobble position of the tRNA anticodons regulate protein translation through tuning the reading of mRNA codons. Among 40 of such nucleosides, there are modified uridines containing either a sulfur atom at the C2 position and/or a substituent at the C5 position of the nucleobase ring. It is already evidenced that tRNAs with 2-thiouridines at the wobble position preferentially read NNA codons, while the reading mode of the NNG codons by R5U/R5S2U-containing anticodons is still elusive. For a series of 18 modified uridines and 2-thiouridines, we determined the pKa values and demonstrated that both modifying elements alter the electron density of the uracil ring and modulate the acidity of their N3H proton. In aqueous solutions at physiological pH the 2-thiouridines containing aminoalkyl C5-substituents are ionized in ca. 50%. The results, confirmed also by theoretical calculations, indicate that the preferential binding of the modified units bearing non-ionizable 5-substituents to guanosine in the NNG codons may obey the alternative C-G-like (Watson–Crick) mode, while binding of those bearing aminoalkyl C5-substituents (protonated under physiological conditions) and especially those with a sulfur atom at the C2 position, adopt a zwitterionic form and interact with guanosine via a ‘new wobble’ pattern.

Ensemble cryo-EM elucidates the mechanism of translation fidelity

Gene translation depends on accurate decoding of mRNA, the structural mechanism of which remains poorly understood. Ribosomes decode mRNA codons by selecting cognate aminoacyl-tRNAs delivered by elongation factor Tu (EF-Tu). Here we present high-resolution structural ensembles of ribosomes with cognate or near-cognate aminoacyl-tRNAs delivered by EF-Tu. Both cognate and near-cognate tRNA anticodons explore the aminoacyl-tRNA-binding site (A site) of an open 30S subunit, while inactive EF-Tu is separated from the 50S subunit. A transient conformation of decoding-centre nucleotide G530 stabilizes the cognate codon-anticodon helix, initiating step-wise 'latching' of the decoding centre. The resulting closure of the 30S subunit docks EF-Tu at the sarcin-ricin loop of the 50S subunit, activating EF-Tu for GTP hydrolysis and enabling accommodation of the aminoacyl-tRNA. By contrast, near-cognate complexes fail to induce the G530 latch, thus favouring open 30S pre-accommodation intermediates with inactive EF-Tu. This work reveals long-sought structural differences between the pre-accommodation of cognate and near-cognate tRNAs that elucidate the mechanism of accurate decoding.

Copper(II)-Thymine Coordination Polymer Nanoribbons as Potential Oligonucleotide Nanocarriers

The direct reaction between copper nitrate, thymine-1-acetic acid, and 4,4'-bipyridine in water leads to the formation of a blue colloid comprising uniform crystalline nanoribbons (length >1 μm; width ca. 150-185 nm; diameter ca. 15-60 nm) of a coordination polymer. The polymer displays a thymine-based structure freely available for supramolecular interactions. These nanostructures show significant selective interaction with single-stranded oligonucleotides based on adenine. Remarkably, they present low cell toxicity in three cell lines-despite the copper(II) content-and can be used as nanocarriers of oligonucleotides. These results suggest the potential of these types of nanostructures in several biological applications.

New Structural Insights into Translational Miscoding

The fidelity of translation depends strongly on the selection of the correct aminoacyl-tRNA that is complementary to the mRNA codon present in the ribosomal decoding center. The ribosome occasionally makes mistakes by selecting the wrong substrate from the pool of aminoacyl-tRNAs. Here, we summarize recent structural advances that may help to clarify the origin of missense errors that occur during decoding. These developments suggest that discrimination between tRNAs is based primarily on steric complementarity and shape acceptance rather than on the number of hydrogen bonds between the molding of the decoding center and the codon-anticodon duplex. They strengthen the hypothesis that spatial mimicry, due either to base tautomerism or ionization, drives infidelity in ribosomal translation.

Tautomers of a Fluorescent G Surrogate and Their Distinct Photophysics Provide Additional Information Channels

Thienoguanosine ((th) G) is an isomorphic nucleoside analogue acting as a faithful fluorescent substitute of G, with respectable quantum yield in oligonucleotides. Photophysical analysis of (th) G reveals the existence of two ground-state tautomers with significantly shifted absorption and emission wavelengths, and high quantum yield in buffer. Using (TD)-DFT calculations, the tautomers were identified as the H1 and H3 keto-amino tautomers. When incorporated into the loop of (-)PBS, the (-)DNA copy of the HIV-1 primer binding site, both tautomers are observed and show differential sensitivity to protein binding. The red-shifted H1 tautomer is strongly favored in matched (-)/(+)PBS duplexes, while the relative emission of the H3 tautomer can be used to detect single nucleotide polymorphisms. These tautomers and their distinct environmental sensitivity provide unprecedented information channels for analyzing G residues in oligonucleotides and their complexes.

The ribosome prohibits the G•U wobble geometry at the first position of the codon-anticodon helix

Precise conversion of genetic information into proteins is essential to cellular health. However, a margin of error exists and is at its highest on the stage of translation of mRNA by the ribosome. Here we present three crystal structures of 70S ribosome complexes with messenger RNA and transfer RNAs and show that when a G•U base pair is at the first position of the codon-anticodon helix a conventional wobble pair cannot form because of inescapable steric clash between the guanosine of the A codon and the key nucleotide of decoding center adenosine 1493 of 16S rRNA. In our structure the rigid ribosomal decoding center, which is identically shaped for cognate or near-cognate tRNAs, forces this pair to adopt a geometry close to that of a canonical G•C pair. We further strengthen our hypothesis that spatial mimicry due either to base tautomerism or ionization dominates the translation infidelity mechanism.

Watson-Crick-like pairs in CCUG repeats: evidence for tautomeric shifts or protonation

RNA transcripts that include expanded CCUG repeats are associated with myotonic dystrophy type 2. Crystal structures of two CCUG-containing oligomers show that the RNA strands associate into slipped duplexes that contain noncanonical C-U pairs that have apparently undergone tautomeric transition or protonation resulting in an unusual Watson-Crick-like pairing. The overhanging ends of the duplexes interact forming U-U pairs, which also show tautomerism. Duplexes consisting of CCUG repeats are thermodynamically less stable than the trinucleotide repeats involved in the TRED genetic disorders, but introducing LNA residues increases their stability and raises the melting temperature of the studied oligomers by ∼10°C, allowing detailed crystallographic studies. Quantum mechanical calculations were performed to test the possibility of the tautomeric transitions or protonation within the noncanonical pairs. The results indicate that tautomeric or ionic shifts of nucleobases can manifest themselves in biological systems, supplementing the canonical "rules of engagement."

Convergent evolution of AUA decoding in bacteria and archaea

Deciphering AUA codons is a difficult task for organisms, because AUA and AUG specify isoleucine (Ile) and methionine (Met), separately. Each of the other purine-ending sense co-don sets (NNR) specifies a single amino acid in the universal genetic code. In bacteria and archaea, the cytidine derivatives, 2-lysylcytidine (L or lysidine) and 2-agmatinylcytidine (agm(2)C or agmatidine), respectively, are found at the first letter of the anticodon of tRNA(Ile) responsible for AUA codons. These modifications prevent base pairing with G of the third letter of AUG codon, and enable tRNA(Ile) to decipher AUA codon specifically. In addition, these modifications confer a charging ability of tRNA(Ile) with Ile. Despite their similar chemical structures, L and agm(2)C are synthesized by distinctive mechanisms and catalyzed by different classes of enzymes, implying that the analogous decoding systems for AUA codons were established by convergent evolution after the phylogenic split between bacteria and archaea-eukaryotes lineages following divergence from the last universal common ancestor (LUCA).

Construction of RNA-Quantum Dot Chimera for Nanoscale Resistive Biomemory Application

RNA nanotechnology offers advantages to construct thermally and chemically stable nanoparticles with well-defined shape and structure. Here we report the development of an RNA-QD (quantum dot) chimera for resistive biomolecular memory application. Each QD holds two copies of the pRNA three-way junction (pRNA-3WJ) of the bacteriophage phi29 DNA packaging motor. The fixed quantity of two RNAs per QD was achieved by immobilizing the pRNA-3WJ with a Sephadex aptamer for resin binding. Two thiolated pRNA-3WJ serve as two feet of the chimera that stand on the gold plate. The RNA nanostructure served as both an insulator and a mediator to provide defined distance between the QD and gold. Immobilization of the chimera nanoparticle was confirmed with scanning tunneling microscopy. As revealed by scanning tunneling spectroscopy, the conjugated pRNA-3WJ-QD chimera exhibited an excellent electrical bistability signal for biomolecular memory function, demonstrating great potential for the development of resistive biomolecular memory and a nano-bio-inspired electronic device for information processing and computing.

Diverse mutants of HIV RRE IIB recognize wild-type Rev ARM or Rev ARM R35G-N40V

The binding of human immunodeficiency virus Rev protein via its arginine-rich motif (ARM) to an internal loop in the Rev-response element region IIB (RRE IIB) is necessary for viral replication. Many variant RNAs and ARMs that bind Rev and RRE IIB have been found. Despite the essential role of Rev asparagine 40 in recognition, the Rev ARM double-mutant R35G-N40V functions well in a Rev-RRE IIB reporter assay, indicating R35G-N40V uses a distinct recognition strategy. To examine how RRE IIB may evolve specificity to wild-type Rev ARM and R35G-N40V, 10 RRE IIB libraries, each completely randomized in overlapping regions, were screened with wild-type Rev ARM and R35G-N40V using a reporter system based on bacteriophage λ N antitermination. Consistent with previous studies, a core element of RRE IIB did not vary, and substitutions occurred at conserved residues only in the presence of other substitutions. Notably, the groove-widening, non-canonical base-pair G48:G71 was mutable to U48:G71 without strong loss of binding to wild-type Rev ARM, suggesting U48:G71 performs the same role by adopting the nearly isosteric, reverse wobble base pair. Originating from RRE IIB, as few as one or two substitutions are sufficient to confer specificity to wild-type Rev or Rev R35G-N40. The diversity of RRE IIB mutants that maintain binding to wild-type Rev ARM and R35G-N40V supports neutral theories of evolution and illustrates paths by which viral RNA-protein interactions can evolve new specificities. Rev-RRE offers an excellent model with which to study the fine structure of how specificity evolves.

Structural basis for a novel mechanism of DNA bridging and alignment in eukaryotic DSB DNA repair

Eukaryotic DNA polymerase mu of the PolX family can promote the association of the two 3'-protruding ends of a DNA double-strand break (DSB) being repaired (DNA synapsis) even in the absence of the core non-homologous end-joining (NHEJ) machinery. Here, we show that terminal deoxynucleotidyltransferase (TdT), a closely related PolX involved in V(D)J recombination, has the same property. We solved its crystal structure with an annealed DNA synapsis containing one micro-homology (MH) base pair and one nascent base pair. This structure reveals how the N-terminal domain and Loop 1 of Tdt cooperate for bridging the two DNA ends, providing a templating base in trans and limiting the MH search region to only two base pairs. A network of ordered water molecules is proposed to assist the incorporation of any nucleotide independently of the in trans templating base. These data are consistent with a recent model that explains the statistics of sequences synthesized in vivo by Tdt based solely on this dinucleotide step. Site-directed mutagenesis and functional tests suggest that this structural model is also valid for Pol mu during NHEJ.

Structural basis for molecular discrimination by a 3',3'-cGAMP sensing riboswitch

Cyclic dinucleotides are second messengers that target the adaptor STING and stimulate the innate immune response in mammals. Besides protein receptors, there are bacterial riboswitches that selectively recognize cyclic dinucleotides. We recently discovered a natural riboswitch that targets 3',3'-cGAMP, which is distinguished from the endogenous mammalian signal 2',3'-cGAMP by its backbone connectivity. Here, we report on structures of the aptamer domain of the 3',3'-cGAMP riboswitch from Geobacter in the 3',3'-cGAMP and c-di-GMP bound states. The riboswitch adopts a tuning fork-like architecture with a junctional ligand-binding pocket and different orientations of the arms are correlated with the identity of the bound cyclic dinucleotide. Subsequent biochemical experiments revealed that specificity of ligand recognition can be affected by point mutations outside of the binding pocket, which has implications for both the assignment and reengineering of riboswitches in this structural class.

2-Thiouracil deprived of thiocarbonyl function preferentially base pairs with guanine rather than adenine in RNA and DNA duplexes

2-Thiouracil-containing nucleosides are essential modified units of natural and synthetic nucleic acids. In particular, the 5-substituted-2-thiouridines (S2Us) present in tRNA play an important role in tuning the translation process through codon-anticodon interactions. The enhanced thermodynamic stability of S2U-containing RNA duplexes and the preferred S2U-A versus S2U-G base pairing are appreciated characteristics of S2U-modified molecular probes. Recently, we have demonstrated that 2-thiouridine (alone or within an RNA chain) is predominantly transformed under oxidative stress conditions to 4-pyrimidinone riboside (H2U) and not to uridine. Due to the important biological functions and various biotechnological applications for sulfur-containing nucleic acids, we compared the thermodynamic stabilities of duplexes containing desulfured products with those of 2-thiouracil-modified RNA and DNA duplexes. Differential scanning calorimetry experiments and theoretical calculations demonstrate that upon 2-thiouracil desulfuration to 4-pyrimidinone, the preferred base pairing of S2U with adenosine is lost, with preferred base pairing with guanosine observed instead. Therefore, biological processes and in vitro assays in which oxidative desulfuration of 2-thiouracil-containing components occurs may be altered. Moreover, we propose that the H2U-G base pair is a suitable model for investigation of the preferred recognition of 3'-G-ending versus A-ending codons by tRNA wobble nucleosides, which may adopt a 4-pyrimidinone-type structural motif.

New structural hypostases of the A·T and G·C Watson–Crick DNA base pairs caused by their mutagenic tautomerisation in a wobble manner: a QM/QTAIM prediction

Our investigation reveals the hitherto unknown ability of the canonical Watson–Crick DNA base pairs to switch into wobble mismatches with mutagenic tautomers, clarifying the nature of genome instability.