Building a stable RNA U-turn with a protonated cytidine

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Structural and dynamic effects of pseudouridine modifications on noncanonical interactions in RNA

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Towards a comprehensive understanding of RNA deamination: synthesis and properties of xanthosine-modified RNA

Nucleobase deamination, such as A-to-I editing, represents an important posttranscriptional modification of RNA. When deamination affects guanosines, a xanthosine (X) containing RNA is generated. However, the biological significance and chemical consequences on RNA are poorly understood. We present a comprehensive study on the preparation and biophysical properties of X-modified RNA. Thermodynamic analyses revealed that base pairing strength is reduced to a level similar to that observed for a G•U replacement. Applying NMR spectroscopy and X-ray crystallography, we demonstrate that X can form distinct wobble geometries with uridine depending on the sequence context. In contrast, X pairing with cytidine occurs either through wobble geometry involving protonated C or in Watson–Crick-like arrangement. This indicates that the different pairing modes are of comparable stability separated by low energetic barriers for switching. Furthermore, we demonstrate that the flexible pairing properties directly affect the recognition of X-modified RNA by reverse transcription enzymes. Primer extension assays and PCR-based sequencing analysis reveal that X is preferentially read as G or A and that the ratio depends on the type of reverse transcriptase. Taken together, our results elucidate important properties of X-modified RNA paving the way for future studies on its biological significance.

Mutations of N1 Riboswitch Affect its Dynamics and Recognition by Neomycin Through Conformational Selection

Short, structured fragments of non-coding mRNA may act as molecular switches upon binding specific ligands, regulating the translation of proteins encoded downstream this mRNA sequence. One switch, called riboswitch N1, is regulated by aminoglycosides such as neomycin. Nucleobase mutations in the apical loop, although distant from the binding pocket, significantly affect neomycin affinity and riboswitch regulatory efficiency. To explain this influence, we conducted molecular dynamics simulations using generalized replica exchange with solute tempering (gREST). Translation assay of a reporter protein in a yeast system shows that mutating A17 to G in the riboswitch apical loop reduces 6-fold the translation regulation efficiency of the mutant. Indeed, simulations of the unbound riboswitch show that G17 frequently stacks with base 7, while base 8 is stabilized towards the binding site in a way that it may interfere with the conformational selection mechanism and decrease riboswitch regulatory activity. In the riboswitch complexes, this single-point A to G mutation disrupts a strong hydrogen bond between nucleotides 5 and 17 and, instead, a new hydrogen bond between residue 17 and neomycin is created. This change forces neomycin to occupy a slightly shifted position in the binding pocket, which increases neomycin flexibility. Our simulations of the U14C mutation suggest that the riboswitch complex with neomycin is more stable if cytosine 14 is protonated. A hydrogen bond between the RNA phosphate and protonated cytosine appears as the stabilizing factor. Also, based on the cell-free translation assay and isothermal titration calorimetry experiments, mutations of nucleotides 14 and 15 affect only slightly the riboswitch ability to bind the ligand and its activity. Indeed, the simulation of the unbound U15A mutant suggests conformations preformed for ligand binding, which may explain slightly higher regulatory activity of this mutant. Overall, our results corroborate the in vivo and in vitro experiments on the N1 riboswitch-neomycin system, detail the relationship between nucleobase mutations and RNA dynamics, and reveal the conformations playing the major role in the conformational selection mechanism.

Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations, and NMR Experiments

Phosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single nonbridging oxygen of the phosphate is replaced with a sulfur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as a tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present molecular dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn: an H-bond, an anion-π interaction, and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulfur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases the dispersion component of the anion-π stacking. The sulfur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to the unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulfur parameters and provide suggestions for improving the force-field performance.

Molecular mechanisms for dynamic regulation of N1 riboswitch by aminoglycosides

A synthetic riboswitch N1, inserted into the 5'-untranslated mRNA region of yeast, regulates gene expression upon binding ribostamycin and neomycin. Interestingly, a similar aminoglycoside, paromomycin, differing from neomycin by only one substituent (amino versus hydroxyl), also binds to the N1 riboswitch, but without affecting gene expression, despite NMR evidence that the N1 riboswitch binds all aminoglycosides in a similar way. Here, to explore the details of structural dynamics of the aminoglycoside-N1 riboswitch complexes, we applied all-atom molecular dynamics (MD) and temperature replica-exchange MD simulations in explicit solvent. Indeed, we found that ribostamycin and neomycin affect riboswitch dynamics similarly but paromomycin allows for more flexibility because its complex lacks the contact between the distinctive 6' hydroxyl group and the G9 phosphate. Instead, a transient hydrogen bond of 6'-OH with A17 is formed, which partially diminishes interactions between the bulge and apical loop of the riboswitch, likely contributing to riboswitch inactivity. In many ways, the paromomycin complex mimics the conformations, interactions, and Na+ distribution of the free riboswitch. The MD-derived interaction network helps understand why riboswitch activity depends on aminoglycoside type, whereas for another aminoglycoside-binding site, aminoacyl-tRNA site in 16S rRNA, activity is not discriminatory.

Coarse-grained dynamic RNA titration simulations

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK _a values for each nucleotide.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pK_a shifts, real pK_a shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pK_as of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into ¹H,³¹P-correlation experiments as well as a BEST-HNP experiment exploiting ^3hJ_N,P rather than ^2hJ_H,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

Fast coarse-grained model for RNA titration

A new numerical scheme for RNA (ribonucleic acid) titration based on the Debye-Hückel framework for the salt description is proposed in an effort to reduce the computational costs for further applications to study protein-RNA systems. By means of different sets of Monte Carlo simulations, we demonstrated that this new scheme is able to correctly reproduce the experimental titration behavior and salt pK_a shifts. In comparison with other theoretical approaches, similar or even better outcomes are achieved at much lower computational costs. The model was tested on the lead-dependent ribozyme, the branch-point helix, and the domain 5 from Azotobacter vinelandii Intron 5.

The crystal structure of the 5΄ functional domain of the transcription riboregulator 7SK

In vertebrates, the 7SK RNA forms the scaffold of a complex, which regulates transcription pausing of RNA-polymerase II. By binding to the HEXIM protein, the complex comprising proteins LARP7 and MePCE captures the positive transcription elongation factor P-TEFb and prevents phosphorylation of pausing factors. The HEXIM-binding site embedded in the 5΄-hairpin of 7SK (HP1) encompasses a short signature sequence, a GAUC repeat framed by single-stranded uridines. The present crystal structure of HP1 shows a remarkably straight helical stack involving several unexpected triples formed at a central region. Surprisingly, two uridines of the signature sequence make triple interactions in the major groove of the (GAUC)₂. The third uridine is turned outwards or inward, wedging between the other uridines, thus filling the major groove. A molecular dynamics simulation indicates that these two conformations of the signature sequence represent stable alternatives. Analyses of the interaction with the HEXIM protein confirm the importance of the triple interactions at the signature sequence. Altogether, the present structural analysis of 7SK HP1 highlights an original mechanism of swapping bases, which could represent a possible ‘7SK signature’ and provides new insight into the functional importance of the plasticity of RNA.

Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

General acid-base catalysis is a key mechanistic strategy in protein and RNA enzymes. Ribozymes use hydrated metal ions, nucleobases, and organic cofactors to carry this out. In most small ribozymes, a guanosine is positioned to participate in proton transfer with the nucleophilic 2'-OH. The unshifted pK_a values for nucleobases and solvated metal ions are far from neutrality, however, and thus nonideal for general acid-base catalysis. Herein, evidence is provided for cooperative interaction in the hammerhead ribozyme among the guanine that interacts with the nucleophilic 2'-OH, G12, the -1 nucleobase C17, and Mg²⁺ ions. We introduce global fitting for analyzing ribozyme rate-pH data parametric in Mg²⁺ concentration and benchmark this method on data from the hepatitis delta virus ribozyme. We then apply global fitting to new rate-pH data for the hammerhead ribozyme using a minimal three-dimensional, four-channel cooperative model. The value for the pK_a of G12 that we obtain is channel-dependent and varies from 8.1 to 9.9, shifting closest toward neutrality in the presence of two cationic species: C17H⁺ and a Mg²⁺ ion. The value for the pK_a of the -1 nucleotide, C17, is increased a remarkable 3.5-5 pK_a units toward neutrality. Shifting of the pK_a of C17 appears to be driven by an electrostatic sandwich of C17 between carbonyl groups of the 5'-neighboring U and of G12 and involves cation-π interactions. Rate-pH profiles reveal that the major reactive channel under biological Mg²⁺ and pH involves a cationic C17 rather than a second metal ion. Substitution of a cationic base for a metal underscores the versatility of RNA.

Revisiting GNRA and UNCG folds: U-turns versus Z-turns in RNA hairpin loops

When thinking about RNA three-dimensional structures, coming across GNRA and UNCG tetraloops is perceived as a boon since their folds have been extensively described. Nevertheless, analyzing loop conformations within RNA and RNP structures led us to uncover several instances of GNRA and UNCG loops that do not fold as expected. We noticed that when a GNRA does not assume its "natural" fold, it adopts the one we typically associate with a UNCG sequence. The same folding interconversion may occur for loops with UNCG sequences, for instance within tRNA anticodon loops. Hence, we show that some structured tetranucleotide sequences starting with G or U can adopt either of these folds. The underlying structural basis that defines these two fold types is the mutually exclusive stacking of a backbone oxygen on either the first (in GNRA) or the last nucleobase (in UNCG), generating an oxygen-π contact. We thereby propose to refrain from using sequences to distinguish between loop conformations. Instead, we suggest using descriptors such as U-turn (for "GNRA-type" folds) and a newly described Z-turn (for "UNCG-type" folds). Because tetraloops adopt for the largest part only two (inter)convertible turns, we are better able to interpret from a structural perspective loop interchangeability occurring in ribosomes and viral RNA. In this respect, we propose a general view on the inclination for a given sequence to adopt (or not) a specific fold. We also suggest how long-noncoding RNAs may adopt discrete but transient structures, which are therefore hard to predict.

A Stably Protonated Adenine Nucleotide with a Highly Shifted pKa Value Stabilizes the Tertiary Structure of a GTP-Binding RNA Aptamer

RNA tertiary structure motifs are stabilized by a wide variety of hydrogen-bonding interactions. Protonated A and C nucleotides are normally not considered to be suitable building blocks for such motifs since their pK_a values are far from physiological pH. Here, we report the NMR solution structure of an in vitro selected GTP-binding RNA aptamer bound to GTP with an intricate tertiary structure. It contains a novel kind of base quartet stabilized by a protonated A residue. Owing to its unique structural environment in the base quartet, the pK_a value for the protonation of this A residue in the complex is shifted by more than 5 pH units compared to the pK_a for A nucleotides in single-stranded RNA. This is the largest pK_a shift for an A residue in structured nucleic acids reported so far, and similar in size to the largest pK_a shifts observed for amino acid side chains in proteins. Both RNA pre-folding and ligand binding contribute to the pK_a shift.

Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding

The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg²⁺ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl₂ were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.

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."