Discrimination among individual Watson-Crick base pairs at the termini of single DNA hairpin molecules

Nanoscale α-hemolysin pores can be used to analyze individual DNA or RNA molecules. Serial examination of hundreds to thousands of molecules per minute is possible using ionic current impedance as the measured property. In a recent report, we showed that a nanopore device coupled with machine learning algorithms could automatically discriminate among the four combinations of Watson–Crick base pairs and their orientations at the ends of individual DNA hairpin molecules. Here we use kinetic analysis to demonstrate that ionic current signatures caused by these hairpin molecules depend on the number of hydrogen bonds within the terminal base pair, stacking between the terminal base pair and its nearest neighbor, and 5′ versus 3′ orientation of the terminal bases independent of their nearest neighbors. This report constitutes evidence that single Watson–Crick base pairs can be identified within individual unmodified DNA hairpin molecules based on their dynamic behavior in a nanoscale pore.     

The interaction networks of structured RNAs

All pairwise interactions occurring between bases which could be detected in three-dimensional structures of crystallized RNA molecules are annotated on new planar diagrams. The diagrams attempt to map the underlying complex networks of base–base interactions and, especially, they aim at conveying key relationships between helical domains: co-axial stacking, bending and all Watson–Crick as well as non-Watson–Crick base pairs. Although such wiring diagrams cannot replace full stereographic images for correct spatial understanding and representation, they reveal structural similarities as well as the conserved patterns and distances between motifs which are present within the interaction networks of folded RNAs of similar or unrelated functions. Finally, the diagrams could help devising methods for meaningfully transforming RNA structures into graphs amenable to network analysis.     

The UA_handle: a versatile submotif in stable RNA architectures

Stable RNAs are modular and hierarchical 3D architectures taking advantage of recurrent structural motifs to form extensive non-covalent tertiary interactions. Sequence and atomic structure analysis has revealed a novel submotif involving a minimal set of five nucleotides, termed the UA_handle motif (5′XU/AN_nX3′). It consists of a U:A Watson–Crick: Hoogsteen trans base pair stacked over a classic Watson–Crick base pair, and a bulge of one or more nucleotides that can act as a handle for making different types of long-range interactions. This motif is one of the most versatile building blocks identified in stable RNAs. It enters into the composition of numerous recurrent motifs of greater structural complexity such as the T-loop, the 11-nt receptor, the UAA/GAN and the G-ribo motifs. Several structural principles pertaining to RNA motifs are derived from our analysis. A limited set of basic submotifs can account for the formation of most structural motifs uncovered in ribosomal and stable RNAs. Structural motifs can act as structural scaffoldings and be functionally and topologically equivalent despite sequence and structural differences. The sequence network resulting from the structural relationships shared by these RNA motifs can be used as a proto-language for assisting prediction and rational design of RNA tertiary structures.     

Requirement for canonical base pairing in the short pseudoknot structure of genomic hepatitis delta virus ribozyme

The tertiary structure of the 3′-cleaved product of the genomic hepatitis delta virus (HDV) ribozyme was solved by X-ray crystallographic analysis. In this structure, three single-stranded regions (SSrA, -B and -C) interact intricately with one another via hydrogen bonds between nucleotide bases, phosphate oxygens and 2′-OHs to form a nested double pseudoknot structure. Among these interactions, two Watson–Crick (W–C) base pairs, 726G–710C and 727G–709C, that form between SSrA and SSrC (P1.1) seem to be especially important for compact folding. To characterize the importance of these base pairs, ribozymes were subjected to in vitro selection from a pool of RNA molecules randomly substituted at positions 709, 710, 726 and 727. The results establish the importance of the two W–C base pairs for activity, although some mutants are active with one G–C base pair. In addition, the kinetic parameters were analyzed in all 16 combinations with two canonical base pairs. Comparison of variant ribozymes with the wild-type ribozyme reveals that the difference in reaction rates for these variants (ΔΔG^‡) is not simply accounted for by the differences in the stability of P1.1 (ΔΔG⁰₃₇). The role played by Mg²⁺ ions in formation of the P1.1 structure is also discussed.     

Studies of translational misreading in vivo show that the ribosome very efficiently discriminates against most potential errors

Protein synthesis must rapidly and repeatedly discriminate between a single correct and many incorrect aminoacyl-tRNAs. We have attempted to measure the frequencies of all possible missense errors by tRNA, tRNA and tRNA. The most frequent errors involve three types of mismatched nucleotide pairs, U•U, U•C, or U•G, all of which can form a noncanonical base pair with geometry similar to that of the canonical U•A or C•G Watson–Crick pairs. Our system is sensitive enough to measure errors at other potential mismatches that occur at frequencies as low as 1 in 500,000 codons. The ribosome appears to discriminate this efficiently against any pair with non-Watson–Crick geometry. This extreme accuracy may be necessary to allow discrimination against the errors involving near Watson–Crick pairing.     

Genetic and structural analysis of base substitutions in the central pseudoknot of Thermus thermophilus 16S ribosomal RNA

Characterization of base substitutions in rRNAs has provided important insights into the mechanism of protein synthesis. Knowledge of the structural effects of such alterations is limited, and could be greatly expanded with the development of a genetic system based on an organism amenable to both genetics and structural biology. Here, we describe the genetic analysis of base substitutions in 16S ribosomal RNA of the extreme thermophile Thermus thermophilus, and an analysis of the conformational effects of these substitutions by structure probing with base-specific modifying agents. Gene replacement methods were used to construct a derivative of strain HB8 carrying a single 16S rRNA gene, allowing the isolation of spontaneous streptomycin-resistant mutants and subsequent genetic mapping of mutations by recombination. The residues altered to give streptomycin resistance reside within the central pseudoknot structure of 16S rRNA comprised of helices 1 and 27, and participate in the U13–U20–A915 base triple, the G21–A914 type II sheared G–A base pair, or the G885–C912 Watson–Crick base pair closing helix 27. Substitutions at any of the three residues engaged in the base triple were found to confer resistance. Results from structure probing of the pseudoknot are consistent with perturbation of RNA conformation by these substitutions, potentially explaining their streptomycin-resistance phenotypes.     

NMR structures of double loops of an RNA aptamer against mammalian initiation factor 4A

A high affinity RNA aptamer (APT58, 58 nt long) against mammalian initiation factor 4A (eIF4A) requires nearly its entire nucleotide sequence for efficient binding. Since splitting either APT58 or eIF4A into two domains diminishes the affinity for each other, it is suggested that multiple interactions or a global interaction between the two molecules accounts for the high affinity. To understand the structural basis of APT58's global recognition of eIF4A, we determined the solution structure of two essential nucleotide loops (AUCGCA and ACAUAGA) within the aptamer using NMR spectroscopy. The AUCGCA loop is stabilized by a U-turn motif and contains a non-canonical A:A base pair (the single hydrogen bond mismatch: Hoogsteen/Sugar-edge). On the other hand, the ACAUAGA loop is stabilized by an AUA tri-nucleotide loop motif and contains the other type of A:A base pair (single hydrogen bond mismatch: Watson–Crick/Watson–Crick). Considering the known structural and functional properties of APT58, we propose that the AUCGCA loop is directly involved in the interaction with eIF4A, while the flexibility of the ACAUAGA loop is important to support this interaction. The Watson–Crick edges of C7 and C9 in the AUCGCA loop may directly interact with eIF4A.     

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem–loop I (SLI) substrate and stem–loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem–loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8–3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7–8 kcal/mol than predicted for a comparable duplex containing three Watson–Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6–7 Watson–Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2–3 Watson–Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

A novel structural rearrangement of hepatitis delta virus antigenomic ribozyme

A bioinformatic covariation analysis of a collection of 119 novel variants of the antigenomic, self-cleaving hepatitis delta virus (HDV) RNA motif supported the formation of all of the Watson–Crick base pairs (bp) of the catalytic centre except the C19–G81 pair located at the bottom of the P2 stem. In fact, a novel Watson–Crick bp between C19 and G80 is suggested by the data. Both chemical and enzymatic probing demonstrated that initially the C19–G81 pair is formed in the ribozyme (Rz), but upon substrate (S) binding and the formation of the P1.1 pseudoknot C19 switches its base-pairing partner from G81 to G80. As a result of this finding, the secondary structure of this ribozyme has been redrawn. The formation of the C19–G80 bp results in a J4/2 junction composed of four nucleotides, similar to that seen in the genomic counterpart, thereby increasing the similarities between these two catalytic RNAs. Additional mutagenesis, cleavage activity and probing experiments yield an original characterization of the structural features involving the residues of the J4/2 junction.     

Crystal structures of complexes between aminoglycosides and decoding A site oligonucleotides: role of the number of rings and positive charges in the specific binding leading to miscoding

The crystal structures of six complexes between aminoglycoside antibiotics (neamine, gentamicin C1A, kanamycin A, ribostamycin, lividomycin A and neomycin B) and oligonucleotides containing the decoding A site of bacterial ribosomes are reported at resolutions between 2.2 and 3.0 Å. Although the number of contacts between the RNA and the aminoglycosides varies between 20 and 31, up to eight direct hydrogen bonds between rings I and II of the neamine moiety are conserved in the observed complexes. The puckered sugar ring I is inserted into the A site helix by stacking against G1491 and forms a pseudo base pair with two H-bonds to the Watson–Crick sites of the universally conserved A1408. This central interaction helps to maintain A1492 and A1493 in a bulged-out conformation. All these structures of the minimal A site RNA complexed to various aminoglycosides display crystal packings with intermolecular contacts between the bulging A1492 and A1493 and the shallow/minor groove of Watson–Crick pairs in a neighbouring helix. In one crystal, one empty A site is observed. In two crystals, two aminoglycosides are bound to the same A site with one bound specifically and the other bound in various ways in the deep/major groove at the edge of the A sites.     

Recognition of T*G mismatched base pairs in DNA by stacked imidazole-containing polyamides: surface plasmon resonance and circular dichroism studies

An imidazole-containing polyamide trimer, f-ImImIm, where f is a formamido group, was recently found using NMR methods to recognize T·G mismatched base pairs. In order to characterize in detail the T·G recognition affinity and specificity of imidazole-containing polyamides, f-ImIm, f-ImImIm and f-PyImIm were synthesized. The kinetics and thermodynamics for the polyamides binding to Watson–Crick and mismatched (containing one or two T·G, A·G or G·G mismatched base pairs) hairpin oligonucleotides were determined by surface plasmon resonance and circular dichroism (CD) methods. f-ImImIm binds significantly more strongly to the T·G mismatch-containing oligonucleotides than to the sequences with other mismatched or with Watson–Crick base pairs. Compared with the Watson–Crick CCGG sequence, f-ImImIm associates more slowly with DNAs containing T·G mismatches in place of one or two C·G base pairs and, more importantly, the dissociation rate from the T·G oligonucleotides is very slow (small k_d). These results clearly demonstrate the binding selectivity and enhanced affinity of side-by-side imidazole/imidazole pairings for T·G mismatches and show that the affinity and specificity increase arise from much lower k_d values with the T·G mismatched duplexes. CD titration studies of f-ImImIm complexes with T·G mismatched sequences produce strong induced bands at ~330 nm with clear isodichroic points, in support of a single minor groove complex. CD DNA bands suggest that the complexes remain in the B conformation.     

Structural effect of the anticancer agent 6-thioguanine on duplex DNA

The incorporation of 6-thioguanine (S6G) into DNA is an essential step in the cytotoxic activity of thiopurines. However, the structural effects of this substitution on duplex DNA have not been fully characterized. Here, we present the solution structures of DNA duplexes containing S6G opposite thymine (S6G·T) and opposite cytosine (S6G·C), solved by high-resolution NMR spectroscopy and restrained molecular dynamics. The data indicate that both duplexes adopt right-handed helical conformations with all Watson–Crick hydrogen bonding in place. The S6G·T structures exhibit a wobble-type base pairing at the lesion site, with thymine shifted toward the major groove and S6G displaced toward the minor groove. Aside from the lesion site, the helices, including the flanking base pairs, are not highly perturbed by the presence of the lesion. Surprisingly, thermal dependence experiments suggest greater stability in the S6G-T mismatch than the S6G-C base pair.     

Molecular dynamics simulations and coupled nucleotide substitution experiments indicate the nature of A·A base pairing and a putative structure of the coralyne-induced homo-adenine duplex

Coralyne is an alkaloid drug that binds homo-adenine DNA (and RNA) oligonucleotides more tightly than it does Watson–Crick DNA. Hud’s laboratory has shown that poly(dA) in the presence of coralyne forms an anti-parallel duplex, however attempts to determine the structure by NMR spectroscopy and X-ray crystallography have been unsuccessful. Assuming adenine–adenine hydrogen bonding between the two poly(dA) strands, we constructed 40 hypothetical homo-(dA) anti-parallel duplexes and docked coralyne into the six most favorable duplex structures. The two most stable structures had trans glycosidic bonds, but distinct pairing geometries, i.e. either Watson–Crick Hoogsteen (transWH) or Watson–Crick Watson–Crick (transWW) with stability of transWH > transWW. To narrow down the possibilities, 7-deaza adenine base substitutions (dA→7) were engineered into homo-(dA) sequences. These substitutions significantly reduced the thermal stability of the coralyne-induced homo-(dA) structure. These experiments strongly suggest the involvement of N7 in the coralyne-induced A·A base pairs. Moreover, due to the differential effect on melting as a function of the location of the dA→7 mutations, these results are consistent with the N1–N7 base pairing of the transWH pairs. Together, the simulation and base substitution experiments predict that the coralyne-induced homo-(dA) duplex structure adopts the transWH geometry.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

The contribution of pseudouridine to stabilities and structure of RNAs

Thermodynamic data are reported revealing that pseudouridine (Ψ) can stabilize RNA duplexes when replacing U and forming Ψ-A, Ψ-G, Ψ-U and Ψ-C pairs. Stabilization is dependent on type of base pair, position of Ψ within the RNA duplex, and type and orientation of adjacent Watson–Crick pairs. NMR spectra demonstrate that for internal Ψ-A, Ψ-G and Ψ-U pairs, the N3 imino proton is hydrogen bonded to the opposite strand nucleotide and the N1 imino proton may also be hydrogen bonded. CD spectra show that general A-helix structure is preserved, but there is some shifting of peaks and changing of intensities. Ψ has two hydrogen donors (N1 and N3 imino protons) and two hydrogen bond acceptors because the glycosidic bond is C-C rather than C-N as in uridine. This greater structural potential may allow Ψ to behave as a kind of structurally driven universal base because it can enhance stability relative to U when paired with A, G, U or C inside a double helix. These structural and thermodynamic properties may contribute to the biological functions of Ψ.     

Isostericity and tautomerism of base pairs in nucleic acids

The natural bases of nucleic acids form a great variety of base pairs with at least two hydrogen bonds between them. They are classified in twelve main families, with the Watson–Crick family being one of them. In a given family, some of the base pairs are isosteric between them, meaning that the positions and the distances between the C1′ carbon atoms are very similar. The isostericity of Watson–Crick pairs between the complementary bases forms the basis of RNA helices and of the resulting RNA secondary structure. Several defined suites of non-Watson–Crick base pairs assemble into RNA modules that form recurrent, rather regular, building blocks of the tertiary architecture of folded RNAs. RNA modules are intrinsic to RNA architecture are therefore disconnected from a biological function specifically attached to a RNA sequence. RNA modules occur in all kingdoms of life and in structured RNAs with diverse functions. Because of chemical and geometrical constraints, isostericity between non-Watson–Crick pairs is restricted and this leads to higher sequence conservation in RNA modules with, consequently, greater difficulties in extracting 3D information from sequence analysis. Nucleic acid helices have to be recognised in several biological processes like replication or translational decoding. In polymerases and the ribosomal decoding site, the recognition occurs on the minor groove sides of the helical fragments. With the use of alternative conformations, protonated or tautomeric forms of the bases, some base pairs with Watson–Crick-like geometries can form and be stabilized. Several of these pairs with Watson–Crick-like geometries extend the concept of isostericity beyond the number of isosteric pairs formed between complementary bases. These observations set therefore limits and constraints to geometric selection in molecular recognition of complementary Watson–Crick pairs for fidelity in replication and translation processes.     

Self-association of short DNA loops through minor groove C:G:G:C tetrads

In addition to the better known guanine-quadruplex, four-stranded nucleic acid structures can be formed by tetrads resulting from the association of Watson–Crick base pairs. When such association occurs through the minor groove side of the base pairs, the resulting structure presents distinctive features, clearly different from quadruplex structures containing planar G-tetrads. Although we have found this unusual DNA motif in a number of cyclic oligonucleotides, this is the first time that this DNA motif is found in linear oligonucleotides in solution, demonstrating that cyclization is not required to stabilize minor groove tetrads in solution. In this article, we have determined the solution structure of two linear octamers of sequence d(TGCTTCGT) and d(TCGTTGCT), and their cyclic analogue d<pCGCTCCGT>, utilizing 2D NMR spectroscopy and restrained molecular dynamics. These three molecules self-associate forming symmetric dimers stabilized by a novel kind of minor groove C:G:G:C tetrad, in which the pattern of hydrogen bonds differs from previously reported ones. We hypothesize that these quadruplex structures can be formed by many different DNA sequences, but its observation in linear oligonucleotides is usually hampered by competing Watson–Crick duplexes.     

Nearest neighbor parameters for Watson-Crick complementary heteroduplexes formed between 2'-O-methyl RNA and RNA oligonucleotides

Results from optical melting studies of Watson–Crick complementary heteroduplexes formed between 2′-O-methyl RNA and RNA oligonucleotides are used to determine nearest neighbor thermodynamic parameters for predicting the stabilities of such duplexes. The results are consistent with the physical model assumed by the individual nearest neighbor-hydrogen bonding model, which contains terms for helix initiation, base pair stacking and base pair composition. The sequence dependence is similar to that for Watson–Crick complementary RNA/RNA duplexes, which suggests that the sequence dependence may also be similar to that for other backbones that favor A-form RNA conformations.