Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.     

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.     

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.     

RNA aptamers that specifically bind to a 16S ribosomal RNA decoding region construct

RNA–RNA recognition is a critical process in controlling many key biological events, such as translation and ribozyme functions. The recognition process governing RNA–RNA interactions can involve complementary Watson–Crick (WC) base pair binding, or can involve binding through tertiary structural interaction. Hence, it is of interest to determine which of the RNA–RNA binding events might emerge through an in vitro selection process. The A-site of the 16S rRNA decoding region was chosen as the target, both because it possesses several different RNA structural motifs, and because it is the rRNA site where codon/anticodon recognition occurs requiring recognition of both mRNA and tRNA. It is shown here that a single family of RNA molecules can be readily selected from two different sizes of RNA library. The tightest binding aptamer to the A-site 16S rRNA construct, 109.2-3, has its consensus sequences confined to a stem–loop region, which contains three nucleotides complementary to three of the four nucleotides in the stem–loop region of the A-site 16S rRNA. Point mutations on each of the three nucleotides on the stem–loop of the aptamer abolish its binding capacity. These studies suggest that the RNA aptamer 109.2-3 interacts with the simple 27 nt A-site decoding region of 16S rRNA through their respective stem–loops. The most probable mode of interaction is through complementary WC base pairing, commonly referred to as a loop–loop ‘kissing’ motif. High affinity binding to the other structural motifs in the decoding region were not observed.     

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.     

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.     

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.     

Analysis of intermolecular base pair formation of prohead RNA of the phage phi29 DNA packaging motor using NMR spectroscopy

The bacteriophage ø29 DNA packaging motor that assembles on the precursor capsid (prohead) contains an essential 174-nt structural RNA (pRNA) that forms multimers. To determine the structural features of the CE- and D-loops believed to be involved in multimerization of pRNA, 35- and 19-nt RNA molecules containing the CE-loop or the D-loop, respectively, were produced and shown to form a heterodimer in a Mg²⁺-dependent manner, similar to that with full-length pRNA. It has been hypothesized that four intermolecular base pairs are formed between pRNA molecules. Our NMR study of the heterodimer, for the first time, proved directly the existence of two intermolecular Watson–Crick G–C base pairs. The two potential intermolecular A–U base pairs were not observed. In addition, flexibility of the D-loop was found to be important since a Watson–Crick base pair introduced at the base of the D-loop disrupted the formation of the intermolecular G–C hydrogen bonds, and therefore affected heterodimerization. Introduction of this mutation into the biologically active 120-nt pRNA (U80C mutant) resulted in no detectable dimerization at ambient temperature as shown by native gel and sedimentation velocity analyses. Interestingly, this pRNA bound to prohead and packaged DNA as well as the wild-type 120-nt pRNA.     

Solution structure and dynamics of DNA duplexes containing the universal base analogues 5-nitroindole and 5-nitroindole 3-carboxamide

Universal bases hybridize with all other natural DNA or RNA bases, and have applications in PCR and sequencing. We have analysed by nuclear magnetic resonance spectroscopy the structure and dynamics of three DNA oligonucleotides containing the universal base analogues 5-nitroindole and 5-nitroindole-3-carboxamide. In all systems studied, both the 5-nitroindole nucleotide and the opposing nucleotide adopt a standard anti conformation and are fully stacked within the DNA duplex. The 5-nitroindole bases do not base pair with the nucleotide opposite them, but intercalate between this base and an adjacent Watson–Crick pair. In spite of their smooth accommodation within the DNA double-helix, the 5-nitroindole-containing duplexes exist as a dynamic mixture of two different stacking configurations exchanging fast on the chemical shift timescale. These configurations depend on the relative intercalating positions of the universal base and the opposing base, and their exchange implies nucleotide opening motions on the millisecond time range. The structure of these nitroindole-containing duplexes explains the mechanism by which these artificial moieties behave as universal bases.     

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.     

A bridged nucleic acid, 2′,4′-BNACOC: synthesis of fully modified oligonucleotides bearing thymine, 5-methylcytosine, adenine and guanine 2′,4′-BNACOC monomers and RNA-selective nucleic-acid recognition

Recently, we synthesized pyrimidine derivatives of the 2′-O,4′-C-methylenoxymethylene-bridged nucleic-acid (2′,4′-BNA^COC) monomer, the sugar conformation of which is restricted in N-type conformation by a seven-membered bridged structure. Oligonucleotides (BNA^COC) containing this monomer show high affinity with complementary single-stranded RNA and significant resistance to nuclease degradation. Here, BNA^COC consisting of 2′,4′-BNA^COC monomers bearing all four bases, namely thymine, 5-methylcytosine, adenine and guanine was efficiently synthesized and properties of duplexes containing the 2′,4′-BNA^COC monomers were investigated by UV melting experiments and circular dichroism (CD) spectroscopy. The UV melting curve analyses showed that the BNA^COC/BNA^COC duplex possessed excellent thermal stability and that the BNA^COC increased thermal stability with a complementary RNA strand. On the other hand, BNA^COC/DNA heteroduplexes showed almost the same thermal stability as RNA/DNA heteroduplexes. Furthermore, mismatched sequence studies showed that BNA^COC generally improved the sequence selectivity with Watson–Crick base-pairing compared to the corresponding natural DNA and RNA. A CD spectroscopic analysis indicated that the BNA^COC formed duplexes with complementary DNA and RNA in a manner similar to natural RNA.     

Closing loop base pairs in RNA loop–loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop–loop RNA duplexes neutralized by Mg²⁺ counterions in aqueous phases. Mismatches GA, CU and Watson–Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop–loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.     

Frequency and isostericity of RNA base pairs

Most of the hairpin, internal and junction loops that appear single-stranded in standard RNA secondary structures form recurrent 3D motifs, where non-Watson–Crick base pairs play a central role. Non-Watson–Crick base pairs also play crucial roles in tertiary contacts in structured RNA molecules. We previously classified RNA base pairs geometrically so as to group together those base pairs that are structurally similar (isosteric) and therefore able to substitute for each other by mutation without disrupting the 3D structure. Here, we introduce a quantitative measure of base pair isostericity, the IsoDiscrepancy Index (IDI), to more accurately determine which base pair substitutions can potentially occur in conserved motifs. We extract and classify base pairs from a reduced-redundancy set of RNA 3D structures from the Protein Data Bank (PDB) and calculate centroids (exemplars) for each base combination and geometric base pair type (family). We use the exemplars and IDI values to update our online Basepair Catalog and the Isostericity Matrices (IM) for each base pair family. From the database of base pairs observed in 3D structures we derive base pair occurrence frequencies for each of the 12 geometric base pair families. In order to improve the statistics from the 3D structures, we also derive base pair occurrence frequencies from rRNA sequence alignments.     

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.     

Integrity of duplex structures without hydrogen bonding: DNA with pyrene paired at abasic sites

DNA polymerases specifically insert the hydrophobic pyrene deoxynucleotide (P) opposite tetrahydrofuran (F), an stable abasic site analog, and DNA duplexes containing this non-hydrogen-bonded pair possess a high degree of thermodynamic stability. These observations support the hypothesis that steric complementarity and stacking interactions may be sufficient for maintaining stability of DNA structure and specificity of DNA replication, even in the absence of hydrogen bonds across the base pair. Here we report the NMR characterization and structure determination of two DNA molecules containing pyrene residues. The first is a 13mer duplex with a pyrene·tetrahydrofuran pair (P·F pair) at the ninth position and the second mimics a replication intermediate right after incorporation of a pyrene nucleoside opposite an abasic site. Our data indicate that both molecules adopt right-handed helical conformations with Watson– Crick alignments for all canonical base pairs. The pyrene ring stays inside the helix close to its baseless partner in both molecules. The single-stranded region of the replication intermediate folds back over the opposing strand, sheltering the hydrophobic pyrene moiety from water exposure. The results support the idea that the stability and replication of a P·F pair is due to its ability to mimic Watson–Crick structure.     

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.     

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.     

Minimising the secondary structure of DNA targets by incorporation of a modified deoxynucleoside: implications for nucleic acid analysis by hybridisation

Some regions of nucleic acid targets are not accessible to heteroduplex formation with complementary oligonucleotide probes because they are involved in secondary structure through intramolecular Watson–Crick pairing. The secondary conformation of the target may be destabilised to assist its interaction with oligonucleotide probes. To achieve this, we modified a DNA target, which has self-complementary sequence able to form a hairpin loop, by replacing dC with N4-ethyldeoxycytidine (d^4EtC), which hybridises specifically with natural dG to give a G:^4EtC base pair with reduced stability compared to the natural G:C base pair. Substitution by d^4EtC greatly reduced formation of the target secondary structure. The lower level of secondary structure allowed hybridisation with complementary probes made with natural bases. We confirmed that hybridisation could be further enhanced by modifying the probes with intercalating groups which stabilise the duplex.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.