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.     

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.     

Structural characterization of a new subclass of panicum mosaic virus-like 3′ cap-independent translation enhancer

Canonical eukaryotic mRNA translation requires 5′cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5′cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3′UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson–Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson–Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.     

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.     

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.     

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.     

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.     

DNA ligases ensure fidelity by interrogating minor groove contacts

DNA ligases, found in both prokaryotes and eukaryotes, covalently link the 3′-hydroxyl and 5′-phosphate ends of duplex DNA segments. This reaction represents a completion step for DNA replication, repair and recombination. It is well established that ligases are sensitive to mispairs present on the 3′ side of the ligase junction, but tolerant of mispairs on the 5′ side. While such discrimination would increase the overall accuracy of DNA replication and repair, the mechanisms by which this fidelity is accomplished are as yet unknown. In this paper, we present the results of experiments with Tth ligase from Thermus thermophilus HB8 and a series of nucleoside analogs in which the mechanism of discrimination has been probed. Using a series of purine analogs substituted in the 2 and 6 positions, we establish that the apparent base pair geometry is much more important than relative base pair stability and that major groove contacts are of little importance. This result is further confirmed using 5-fluorouracil (FU) mispaired with guanine. At neutral pH, the FU:G mispair on the 3′ side of a ligase junction is predominantly in a neutral wobble configuration and is poorly ligated. Increasing the solution pH increases the proportion of an ionized base pair approximating Watson–Crick geometry, substantially increasing the relative ligation efficiency. These results suggest that the ligase could distinguish Watson–Crick from mispaired geometry by probing the hydrogen bond acceptors present in the minor groove as has been proposed for DNA polymerases. The significance of minor groove hydrogen bonding interactions is confirmed with both Tth and T4 DNA ligases upon examination of base pairs containing the pyrimidine shape analog, difluorotoluene (DFT). Although DFT paired with adenine approximates Watson–Crick geometry, a minor groove hydrogen bond acceptor is lost. Consistent with this hypothesis, we observe that DFT-containing base pairs inhibit ligation when on the 3′ side of the ligase junction. The NAD⁺-dependent ligase, Tth, is more sensitive to the DFT analog on the unligated strand whereas the ATP-dependent T4 ligase is more sensitive to substitutions in the template strand. Electrophoretic gel mobility-shift assays demonstrate that the Tth ligase binds poorly to oligonucleotide substrates containing analogs with altered minor groove contacts.     

Riboswitch structure: an internal residue mimicking the purine ligand

The adenine and guanine riboswitches regulate gene expression in response to their purine ligand. X-ray structures of the aptamer moiety of these riboswitches are characterized by a compact fold in which the ligand forms a Watson–Crick base pair with residue 65. Phylogenetic analyses revealed a strict restriction at position 39 of the aptamer that prevents the G39–C65 and A39–U65 combinations, and mutational studies indicate that aptamers with these sequence combinations are impaired for ligand binding. In order to investigate the rationale for sequence conservation at residue 39, structural characterization of the U65C mutant from Bacillus subtilis pbuE adenine riboswitch aptamer was undertaken. NMR spectroscopy and X-ray crystallography studies demonstrate that the U65C mutant adopts a compact ligand-free structure, in which G39 occupies the ligand-binding site of purine riboswitch aptamers. These studies present a remarkable example of a mutant RNA aptamer that adopts a native-like fold by means of ligand mimicking and explain why this mutant is impaired for ligand binding. Furthermore, this work provides a specific insight into how the natural sequence has evolved through selection of nucleotide identities that contribute to formation of the ligand-bound state, but ensures that the ligand-free state remains in an active conformation.     

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.     

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.     

Structure of the cytosine–cytosine mismatch in the thymidylate synthase mRNA binding site and analysis of its interaction with the aminoglycoside paromomycin

The structure of a cytosine–cytosine (CC) mismatch-containing RNA molecule derived from a hairpin structure in the thymidylate synthase mRNA that binds the aminoglycoside paromomycin with high affinity was determined using nuclear magnetic resonance (NMR) spectroscopy. The cytosines in the mismatch form a noncanonical base pair where both cytosines are uncharged and stack within the stem of the RNA structure. Binding to paromomycin was analyzed using isothermal titration calorimetry (ITC) to demonstrate the necessity of the CC mismatch and to determine the affinity dissociation constant of this RNA to paromomycin to be 0.5 ± 0.3 μM. The CC mismatch, and the neighboring GC base pairs experienced the highest degree of chemical shift changes in their H6 and H5 resonances indicating that paromomycin binds in the major groove at the CC mismatch site. In comparing the structure of CC mismatch RNA with a fully Watson–Crick GC base paired stem, the CC mismatch is shown to confer a widening of the major groove. This widening, combined with the dynamic nature of the CC mismatch, enables binding of paromomycin to this RNA molecule.     

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.     

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.     

Probing the structural hierarchy and energy landscape of an RNA T-loop hairpin

The T-loop motif is an important recurrent RNA structural building block consisting of a U-turn sub-motif and a UA trans Watson–Crick/Hoogsteen base pair. In the presence of a hairpin stem, the UA non-canonical base pair becomes part of the UA-handle motif. To probe the hierarchical organization and energy landscape of the T-loop, we performed replica exchange molecular dynamics (REMD) simulations of the T-loop in isolation and as part of a hairpin. Our simulations reveal that the isolated T-loop adopts coil conformers stabilized by base stacking. The T-loop hairpin shows a highly rugged energy landscape featuring multiple local minima with a transition state for folding consisting of partially zipped states. The U-turn displays a high conformational flexibility both when the T-loop is in isolation and as part of a hairpin. On the other hand, the stability of the UA non-canonical base pair is enhanced in the presence of the UA-handle. This motif is apparently a key component for stabilizing the T-loop, while the U-turn is mostly involved in long-range interaction. Our results suggest that the stability and folding of small RNA motifs are highly dependent on local context.     

The solution structure of an oligonucleotide duplex containing a 2'-deoxyadenosine-3-(2-hydroxyethyl)- 2'-deoxyuridine base pair determined by NMR and molecular dynamics studies

Determination of the solution structure of the duplex d(GCAAGTC(HE)AAAACG)·d(CGTTTTAGACTTGC) containing a 3-(2-hydroxyethyl)-2′-deoxyuridine·deoxyadenine (HE·A) base pair is reported. The three-dimensional solution structure, determined starting from 512 models via restrained molecular mechanics using inter-proton distances and torsion angles, converged to two final families of structures. For both families the HE and the opposite A residues are intrahelical and in the anti conformation. The hydroxyethyl chain lies close to the helix axis and for one family the hydroxyl group is above the HE·A plane and in the other case it is below. These two models were used to start molecular dynamic calculations with explicit solvent to explore the hydrogen bonding possibilities of the HE·A base pair. The dynamics calculations converge finally to one model structure in which two hydrogen bonds are formed. The first is formed all the time and is between HEO4 and the amino group of A, and the second, an intermittent one, is between the hydroxyl group and the N1 of A. When this second hydrogen bond is not formed a weak interaction CH···N is possible between HEC7H2 and N1A21. All the best structures show an increase in the C1′–C1′ distance relative to a Watson–Crick base pair.