Facilitating RNA structure prediction with microarrays

Measuring intramolecular connectivity in long RNA molecules using two-dimensional DNA patch-probe arrays

We describe a simple method to infer intramolecular connections in a population of long RNA molecules in vitro. First we add DNA oligonucleotide "patches" that perturb the RNA connections, then we use a microarray containing a complete set of DNA oligonucleotide "probes" to record where perturbations occur. The pattern of perturbations reveals couplings between different regions of the RNA sequence, from which we infer connections as well as their prevalences in the population. We validate this patch-probe method using the 1,058-nucleotide RNA genome of satellite tobacco mosaic virus (STMV), which has previously been shown to have multiple long-range connections. Our results not only indicate long duplexes that agree with previous structures but also reveal the prevalence of competing connections. Together, these results suggest that globally-folded and locally-folded structures coexist in solution. We show that the prevalence of connections changes when pseudouridine, an important component of natural and synthetic RNA molecules, is substituted for uridine in STMV RNA.

Duplexed aptamers: history, design, theory, and application to biosensing

Nucleic acid aptamers are single stranded DNA or RNA sequences that specifically bind a cognate ligand. In addition to their widespread use as stand-alone affinity binding reagents in analytical chemistry, aptamers have been engineered into a variety of ligand-specific biosensors, termed aptasensors. One of the most common aptasensor formats is the duplexed aptamer (DA). As defined herein, DAs are aptasensors containing two nucleic acid elements coupled via Watson-Crick base pairing: (i) an aptamer sequence, which serves as a ligand-specific receptor, and (ii) an aptamer-complementary element (ACE), such as a short DNA oligonucleotide, which is designed to hybridize to the aptamer. The ACE competes with ligand binding, such that DAs generate a signal upon ligand-dependent ACE-aptamer dehybridization. DAs possess intrinsic advantages over other aptasensor designs. For example, DA biosensing designs generalize across DNA and RNA aptamers, DAs are compatible with many readout methods, and DAs are inherently tunable on the basis of nucleic acid hybridization. However, despite their utility and popularity, DAs have not been well defined in the literature, leading to confusion over the differences between DAs and other aptasensor formats. In this review, we introduce a framework for DAs based on ACEs, and use this framework to distinguish DAs from other aptasensor formats and to categorize cis- and trans-DA designs. We then explore the ligand binding dynamics and chemical properties that underpin DA systems, which fall under conformational selection and induced fit models, and which mirror classical SN1 and SN2 models of nucleophilic substitution reactions. We further review a variety of in vitro and in vivo applications of DAs in the chemical and biological sciences, including riboswitches and riboregulators. Finally, we present future directions of DAs as ligand-responsive nucleic acids. Owing to their tractability, versatility and ease of engineering, DA biosensors bear a great potential for the development of new applications and technologies in fields ranging from analytical chemistry and mechanistic modeling to medicine and synthetic biology.

Secondary structure of the segment 5 genomic RNA of influenza A virus and its application for designing antisense oligonucleotides

Influenza virus causes seasonal epidemics and dangerous pandemic outbreaks. It is a single stranded (-)RNA virus with a segmented genome. Eight segments of genomic viral RNA (vRNA) form the virion, which are then transcribed and replicated in host cells. The secondary structure of vRNA is an important regulator of virus biology and can be a target for finding new therapeutics. In this paper, the secondary structure of segment 5 vRNA is determined based on chemical mapping data, free energy minimization and structure-sequence conservation analysis for type A influenza. The revealed secondary structure has circular folding with a previously reported panhandle motif and distinct novel domains. Conservations of base pairs is 87% on average with many structural motifs that are highly conserved. Isoenergetic microarray mapping was used to additionally validate secondary structure and to discover regions that easy bind short oligonucleotides. Antisense oligonucleotides, which were designed based on modeled secondary structure and microarray mapping, inhibit influenza A virus proliferation in MDCK cells. The most potent oligonucleotides lowered virus titer by ~90%. These results define universal for type A structured regions that could be important for virus function, as well as new targets for antisense therapeutics.

Self-Folding of Naked Segment 8 Genomic RNA of Influenza A Virus

Influenza A is a negative sense RNA virus that kills hundreds of thousands of humans each year. Base pairing in RNA is very favorable, but possibilities for RNA secondary structure of the influenza genomic RNA have not been investigated. This work presents the first experimentally-derived exploration of potential secondary structure in an influenza A naked (protein-free) genomic segment. Favorable folding regions are revealed by in vitro chemical structure mapping, thermodynamics, bioinformatics, and binding to isoenergetic microarrays of an entire natural sequence of the 875 nt segment 8 vRNA and of a smaller fragment. Segment 8 has thermodynamically stable and evolutionarily conserved RNA structure and encodes essential viral proteins NEP and NS1. This suggests that vRNA self-folding may generate helixes and loops that are important at one or more stages of the influenza life cycle.

Experiment-Assisted Secondary Structure Prediction with RNAstructure

Experimental probing data can be used to improve the accuracy of RNA secondary structure prediction. The software package RNAstructure can take advantage of enzymatic cleavage data, FMN cleavage data, traditional chemical modification reactivity data, and SHAPE reactivity data for secondary structure modeling. This chapter provides protocols for using experimental probing data with RNAstructure to restrain or constrain RNA secondary structure prediction.

Efficient Design of Compact Unstructured RNA Libraries Covering All k-mers

Current microarray technologies to determine RNA structure or measure protein–RNA interactions rely on single-stranded, unstructured RNA probes on a chip covering together all k-mers. Since space on the array is limited, the problem is to efficiently design a compact library of unstructured ℓ-long RNA probes, where each k-mer is covered at least p times. Ray et al. designed such a library for specific values of k, ℓ, and p using ad-hoc rules. To our knowledge, there is no general method to date to solve this problem. Here, we address the problem of finding a minimum-size covering of all k-mers by ℓ-long sequences with the desired properties for any value of k, ℓ, and p. As we prove that the problem is NP-hard, we give two solutions: the first is a greedy algorithm with a logarithmic approximation ratio; the second, a heuristic greedy approach based on random walks in de Bruijn graphs. The heuristic algorithm works well in practice and produces a library of unstructured RNA probes that is only ∼1.1-times greater in size compared to the theoretical lower bound. We present results for typical values of k and probe lengths ℓ and show that our algorithm generates a library that is significantly smaller than the library of Ray et al.; moreover, we show that our algorithm outperforms naive methods. Our approach can be generalized and extended to generate RNA or DNA oligo libraries with other desired properties. The software is freely available online.

Structural determinants for alternative splicing regulation of the MAPT pre-mRNA

Alternative splicing at the MAPT gene exon 10 yields similar levels of the 3R and 4R tau protein isoforms. (1) The presence of mutations, particularly in exon 10 and intron 10-11, changes the quantity of tau isoforms. Domination each of the isoform yields tau protein aggregation and frontotemporal dementia and Parkinsonism linked to chromosome 17 (FTDP-17). Here, we report for the first time the secondary structure of the 194/195 nucleotide region for the wild type (WT) and 10 mutants of the MAPT gene pre-mRNA determined using both chemical and microarray mapping. Thermodynamic analyses indicate that single nucleotide mutations in the splicing regulatory element (SRE) that form a hairpin affect its stability by up to 4 and 7 kcal/mol. Moreover, binding the regulatory hairpin of small molecule ligands (neomycin, kanamycin, tobramycin and mitoxantrone) enhance its stability depending on the nature of the ligands and the RNA mutations. Experiments using the cos-7 cell line indicate that the presence of ligands and modified antisense oligonucleotides affect the quantity of 3R and 4R isoforms. This finding correlates with the thermodynamic stability of the regulatory hairpin. An alternative splicing regulation mechanism for exon 10 is postulated based on our experimental data and on published data.

A Conserved Secondary Structural Element in the Coding Region of the Influenza A Virus Nucleoprotein (NP) mRNA Is Important for the Regulation of Viral Proliferation

Influenza A virus is a threat to humans due to seasonal epidemics and infrequent, but dangerous, pandemics that lead to widespread infection and death. Eight segments of RNA constitute the genome of this virus and they encode greater than eight proteins via alternative splicing of coding (+)RNAs generated from the genomic (-)RNA template strand. RNA is essential in its life cycle. A bioinformatics analysis of segment 5, which encodes nucleoprotein, revealed a conserved structural motif in the (+)RNA. The secondary structure proposed by energy minimization and comparative analysis agrees with structure predicted based on experimental data using a 121 nucleotide in vitro RNA construct comprising an influenza A virus consensus sequence and also an entire segment 5 (+)RNA (strain A/VietNam/1203/2004 (H5N1)). The conserved motif consists of three hairpins with one being especially thermodynamically stable. The biological importance of this conserved secondary structure is supported in experiments using antisense oligonucleotides in cell line, which found that disruption of this motif led to inhibition of viral fitness. These results suggest that this conserved motif in the segment 5 (+)RNA might be a candidate for oligonucleotide-based antiviral therapy.

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

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

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.

Identifying and characterizing Hfq-RNA interactions

To regulate stress responses and virulence, bacteria use small regulatory RNAs (sRNAs). These RNAs can up or down regulate target mRNAs through base pairing by influencing ribosomal access and RNA decay. A large class of these sRNAs, called trans-encoded sRNAs, requires the RNA binding protein Hfq to facilitate base pairing between the regulatory RNA and its target mRNA. The resulting network of regulation is best characterized in Escherichia coli and Salmonella typhimurium, but the importance of Hfq dependent sRNA regulation is recognized in a diverse population of bacteria. In this review we present the approaches and methods used to discover Hfq binding RNAs, characterize their interactions and elucidate their functions.

Computational prediction of RNA secondary structure

The purpose of this section is to detail methods for the computational prediction of RNA secondary structure. This protocol is intended to provide an easy entry into the field of RNA structure prediction for those wishing to utilize it in their research and to suggest 'best practices' for going from sequence to secondary structure depending on the available data.

Mapping targetable sites on human telomerase RNA pseudoknot/template domain using 2'-OMe RNA-interacting polynucleotide (RIPtide) microarrays

Most cellular RNAs engage in intrastrand base-pairing that gives rise to complex three-dimensional folds. This self-pairing presents an impediment toward binding of the RNA by nucleic acid-based ligands. An important step in the discovery of RNA-targeting ligands is therefore to identify those regions in a folded RNA that are accessible toward the nucleic acid-based ligand. Because the folding of RNA targets can involve interactions between nonadjacent regions and employ both Watson-Crick and non-Watson-Crick base-pairing, screening of candidate binder ensembles is typically necessary. Microarray-based screening approaches have shown great promise in this regard and have suggested that achieving complete sequence coverage would be a valuable attribute of a next generation system. Here, we report a custom microarray displaying a library of RNA-interacting polynucleotides comprising all possible 2'-OMe RNA sequences from 4- to 8-nucleotides in length. We demonstrate the utility of this array in identifying RNA-interacting polynucleotides that bind tightly and specifically to the highly conserved, functionally essential template/pseudoknot domain of human telomerase RNA and that inhibit telomerase function in vitro.

Isoenergetic microarrays to study the structure and interactions of DsrA and OxyS RNAs in two- and three-component complexes

Information on the secondary structure and interactions of RNA is important to understand the biological function of RNA as well as in applying RNA as a tool for therapeutic purposes. Recently, the isoenergetic microarray mapping method was developed to improve the prediction of RNA secondary structure. Herein, for the first time, isoenergetic microarrays were used to study the binding of RNA to protein or other RNAs as well as the interactions of two different RNAs and protein in a three-component complex. The RNAs used as models were the regulatory DsrA and OxyS RNAs from Escherichia coli, the fragments of their target mRNAs (fhlA and rpoS), and their complexes with Hfq protein. The collected results showed the advantages and some limitations of microarray mapping.

Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain

Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).

MicroRNAs align with accessible sites in target mRNAs

The importance of microRNAs (miRs) in control of gene expression is now clearly recognized. While individual microRNAs are thought to target hundreds of disparate mRNAs via imperfect base pairing, little is known about the characteristics of miR target sites. Here we show that the miRs can be aligned with empirically identified accessible sites in a target RNA (Cytokeratin 19, KRT), and that some of the aligned miRs functionally down-regulate KRT expression post-transcriptionally. We employed an RNase-H-based random library selection protocol to identify accessible sites in KRT RNA. We then aligned the Sanger Institute database collection of human miRs to KRT mRNA, and also aligned them using the web-based MicroInspector program. Most miRs aligned with the accessible sites identified empirically; those not aligned with the empirically identified sites also functioned effectively in RNase-H-based assays. Similar results were obtained with a second target RNA (Mammoglobin). Transient transfection assays established that some of the miRs which aligned with KRT significantly down-regulated it at the protein level, with no effect on RNA level. The functionally effective miRs aligned within the coding region of KRT, whereas a number of miRs which aligned with the 3'-untranslated region did not produce down-regulation.

Rapid determination of RNA accessible sites by surface plasmon resonance detection of hybridization to DNA arrays

RNA accessible sites are the regions in an RNA molecule that are available for hybridization with cDNA or RNA molecules. The identification of these accessible sites is a critical first step in identifying antisense-mediated gene suppression sites, as well as in a variety of other RNA-based analysis methods. Here, we present a rapid, hybridization-based, label-free method of identifying RNA accessible sites with surface plasmon resonance imaging (SPRi) on in situ synthesized oligonucleotide arrays prepared on carbon-on-metal substrates. The accessible sites of three pre-miRNAs, miRNA precursors of approximately 75 nt in length, were determined by hybridizing the RNA molecules to RNA-specific tiling arrays. An array composed of all possible 6mer oligonucleotide sequences was also utilized in this work, offering a universal platform capable of studying RNA molecules in a high throughput manner.

Binding of short oligonucleotides to RNA: studies of the binding of common RNA structural motifs to isoenergetic microarrays

Binding of short oligonucleotides to RNA is important for many biological processes. On the basis of RNAi phenomena, antisense, and ribozyme approaches, it is useful in the inhibition of biological functions. To be considered as potential therapeutics, oligonucleotides must bind strongly and selectively to a complementary fragment of target RNA. Microarray technologies also involve the binding of oligonucleotide probes to DNA or RNA. Herein, the hybridization of common structural motifs of RNA, i.e., hairpins, internal loops, bulges, 3'- and 5'-dangling ends, and pseudoknots to isoenergetic microarray probes is presented. The analysis demonstrates that microarray probes bind to bulges, internal loops, and dangling ends as expected. Probes may also bind to terminal helixes, however, possibly due to the rearrangement of base pairs. These results suggest that isoenergetic microarray mapping can provide data to facilitate and improve RNA secondary structure prediction. However, optimal results require combination with chemical and/or enzymatic mapping.

Chemical synthesis of LNA-2-thiouridine and its influence on stability and selectivity of oligonucleotide binding to RNA

Hybridization to RNA is important for many applications, including antisense therapeutics, RNA interference, and microarray screening. Similar thermodynamic stabilities of A-U and G-U base pairs result in difficulties in selective binding to RNA. Moreover, A-U pairs are weaker than G-C pairs so that binding is sometimes weak when many A-U pairs are present. It is known, however, that replacement of uridine with 2-thiouridine significantly improves binding and selectivity. To test for additional improvement of binding and of the specificity for binding A over G, LNA-2-thiouridine was synthesized for the first time and incorporated into many LNA-2'-O-methyl-RNA/RNA duplexes. UV melting was used to measure the thermodynamic effect of replacing 2'-O-methyluridine with 2'-O-methyl-2-thiouridine or LNA-2-thiouridine. The 2-thiouridine usually enhances binding and selectivity. Selectivity is optimized when a single 2-thiouridine is placed at an internal position in a duplex.

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNA^Phe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules.

Contributions of stacking, preorganization, and hydrogen bonding to the thermodynamic stability of duplexes between RNA and 2'-O-methyl RNA with locked nucleic acids

Locked nucleic acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplexes. We report the thermodynamic stabilities of LNA-2'-O-methyl RNA/RNA duplexes designed to provide insight into the contributions of stacking and hydrogen bonding interactions to the enhanced stability. The results show that hydrogen bonding of LNA nucleotides is similar to that of 2'-O-methyl RNA nucleotides, whereas the 3'-stacking interactions are on average approximately 0.7 kcal/mol more favorable at 37 degrees C than for 2'-O-methyl or RNA nucleotides. Moreover, NMR spectra suggest helical preorganization of the single-stranded tetramer, C(L)A(M)A(L)U(M), probably due to restriction of some torsion angles. Thus, enhanced stacking interactions and helical preorganization of single-stranded oligonucleotides contribute to the extraordinary stabilization of duplexes by LNA nucleotides.

Secondary structures for 5' regions of R2 retrotransposon RNAs reveal a novel conserved pseudoknot and regions that evolve under different constraints

Sequences from the 5' region of R2 retrotransposons of four species of silk moth are reported. In Bombyx mori, this region of the R2 messenger RNA contains a binding site for R2 protein and mediates interactions critical to R2 element insertion into the host genome. A model of secondary structure for a segment of this RNA is proposed on the basis of binding to oligonucleotide microarrays, chemical mapping, and comparative sequence analysis. Five conserved secondary structures are identified, including a novel pseudoknot. There is an apparent transition from an entirely RNA structure coding function in most of the 5' segment to a protein coding function near the 3' end. This suggests that local regions evolved under separate functional constraints (structural, coding, or both).

LNA-modified primers drastically improve hybridization to target RNA and reverse transcription

Knowledge about the structure of RNA is crucial to understanding its biological activities. Very often, the presence of unusually thermodynamically stable structural fragments in RNAs, such as hairpins, makes it impossible to apply primer extension to visualize the results of chemical mapping experiments. However, replacement of DNA primers with LNA-modified primers overcomes this limitation. This approach was tested successfully on regulatory OxyS RNA and DsrA RNA from Escherichia coli.

Properties of pseudo-complementary DNA substituted with weakly pairing analogs of guanine or cytosine

A straightforward enzymatic protocol for converting regular DNA into pseudo-complementary DNA could improve the performance of oligonucleotide microarrays by generating readily hybridizable structure-free targets. Here we screened several highly destabilizing analogs of G and C for one that could be used with 2-aminoadenine (nA) and 2-thiothymine (sT) to generate structure-free DNA that is fully accessible to complementary probes. The analogs, which included bioactive bases such as 6-thioguanine (sG), 5-nitrocytosine (NitroC), 2-pyrimidinone (P; the free base of zebularine) and 6-methylfuranopyrimidinone (MefP), were prepared as dNTPs and evaluated as substrates for T7 and Phi29 DNA polymerases that lacked editor function. Pairing properties of the analogs were characterized by solution hybridization assays using modified oligonucleotides or primer extension products. P and MeP did not support robust primer extension whereas sG and NitroC did. In hybridization assays, however, sG lacked discrimination and NitroC paired too strongly to C. The dNTPs of two other base analogs, 7-nitro-7-deazahypoxanthine (NitrocH) and 2-thiocytosine (sC), exhibited the greatest promise. Either analog could be used with nA and sT to generate DNA that was nearly structure-free. Hybridization of probes to these modified DNAs will require the development of base analogs that pair strongly to NitrocH or sC.

NMR-assisted prediction of RNA secondary structure: identification of a probable pseudoknot in the coding region of an R2 retrotransposon

As the rate of functional RNA sequence discovery escalates, high-throughput techniques for reliable structural determination are becoming crucial for revealing the essential features of these RNAs in a timely fashion. Computational predictions of RNA secondary structure quickly generate reasonable models but suffer from several approximations, including overly simplified models and incomplete knowledge of significant interactions. Similar problems limit the accuracy of predictions for other self-folding polymers, including DNA and peptide nucleic acid (PNA). The work presented here demonstrates that incorporating unassigned data from simple nuclear magnetic resonance (NMR) experiments into a dynamic folding algorithm greatly reduces the potential folding space of a given RNA and therefore increases the confidence and accuracy of modeling. This procedure has been packaged into an NMR-assisted prediction of secondary structure (NAPSS) algorithm that can produce pseudoknotted as well as non-pseudoknotted secondary structures. The method reveals a probable pseudoknot in the part of the coding region of the R2 retrotransposon from Bombyx mori that orchestrates second-strand DNA cleavage during insertion into the genome.

Enzymatic synthesis of structure-free DNA with pseudo-complementary properties

Long single-stranded DNAs and RNAs possess considerable secondary structure under conditions that support stable hybrid formation with oligonucleotides. Consequently, different oligomeric probes can hybridize to the same target with efficiencies that vary by several orders of magnitude. The ability to enzymatically generate structure-free single-stranded copies of any nucleic acid without impairing Watson–Crick base pairing to short probes would eliminate this problem and significantly improve the performance of many oligonucleotide-based applications. Synthetic nucleic acids that exhibit these properties are defined as pseudo-complementary. Previously, we described a pseudo-complementary A-T couple consisting of 2-aminoadenine (nA) and 2-thiothymine (sT) bases. The nA-sT couple is a mismatch even though nA-T and A-sT are stable base pairs. Here we show that 7-alkyl-7-deazaguanine and N⁴-alkylcytosine (where alkyl = methyl or ethyl) can be used in conjunction with nA and sT to render DNA largely structure-free and pseudo-complementary. The deoxynucleoside triphosphates (dNTPs) of these bases are incorporated into DNA by selected mesophilic and thermophilic DNA polymerases and the resulting primer extension products hybridize with good specificity and stability to oligonucleotide probes composed of the standard bases. Further optimization and characterization of the synthesis and properties of pseudo-complementary DNA should lead to an ideal target for use with oligonucleotide probes that are <25 nt in length.

Study of the base discrimination ability of DNA and 2'-O-methylated RNA oligomers containing 2-thiouracil bases towards complementary RNA or DNA strands and their application to single base mismatch detection

Precise detection of target DNA and RNA sequences using chemically modified oligonucleotides is of crucial importance in gene analysis and gene silence. The hybridisation and base discrimination abilities of oligonucleotides containing 2'-O-methyl-2-thiouridine (s(2)Um) in homo- and hetero-duplexes composed of DNA and RNA strands have been studied in detail. When s(2)Um was incorporated into RNA or DNA strands, the hybridisation and base discrimination abilities of the modified RNA or DNA oligomers towards the complementary RNA strands were superior to those of the corresponding unmodified oligomers. On the other hand, their base discrimination abilities towards complementary DNA strands were almost the same as those of the unmodified ones. The base discrimination abilities of 2-thiouracil base-containing oligonucleotide probes on slide glass plates were also studied. These modified probes exhibited efficient detection of mismatched base pairing.

Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2′ oxygen and 4′ carbon of ribose) and 2,6-diaminopurine were incorporated into 2′-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5′ coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5′ RNA). This R2Bm 5′ RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5′ RNA.

Selective quenching of fluorescence from unbound oligonucleotides by gold nanoparticles as a probe of RNA structure

Binding of small oligonucleotides to the periphery of folded RNA can provide insight into the secondary structure of complex RNA in solution. To discriminate between bound and unbound fluorescein-labeled 2'-O-methyl RNA probes, we use ionically coated gold nanoparticles to selectively adsorb unbound probes and quench their fluorescence. The target is the 3' untranslated region of Bombyx mori R2 RNA. Fluorescence indicates that R2 sequences complementary to some of the probes are accessible for binding in the three-dimensional structure. Hybridization occurs under homogeneous conditions in the absence of the gold nanoparticles so that steric issues associated with chip-based assays are avoided. The assay is compatible with well plate formats, takes less than 5 min, and requires only 2 pmol or less of unlabeled target RNA per probe sequence tested.

Changes in the conformation of 5S rRNA cause alterations in principal functions of the ribosomal nanomachine

5S rRNA is an integral component of the large ribosomal subunit in virtually all living organisms. Polyamine binding to 5S rRNA was investigated by cross-linking of N¹-azidobenzamidino (ABA)-spermine to naked 5S rRNA or 50S ribosomal subunits and whole ribosomes from Escherichia coli cells. ABA-spermine cross-linking sites were kinetically measured and their positions in 5S rRNA were localized by primer extension analysis. Helices III and V, and loops A, C, D and E in naked 5S rRNA were found to be preferred polyamine binding sites. When 50S ribosomal subunits or poly(U)-programmed 70S ribosomes bearing tRNA^Phe at the E-site and AcPhe-tRNA at the P-site were targeted, the susceptibility of 5S rRNA to ABA-spermine was greatly reduced. Regardless of 5S rRNA assembly status, binding of spermine induced significant changes in the 5S rRNA conformation; loop A adopted an apparent ‘loosening’ of its structure, while loops C, D, E and helices III and V achieved a more compact folding. Poly(U)-programmed 70S ribosomes possessing 5S rRNA cross-linked with spermine were more efficient than control ribosomes in tRNA binding, peptidyl transferase activity and translocation. Our results support the notion that 5S rRNA serves as a signal transducer between regions of 23S rRNA responsible for principal ribosomal functions.

A chemical synthesis of LNA-2,6-diaminopurine riboside, and the influence of 2'-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes

Modified nucleotides are useful tools to study the structures, biological functions and chemical and thermodynamic stabilities of nucleic acids. Derivatives of 2,6-diaminopurine riboside (D) are one type of modified nucleotide. The presence of an additional amino group at position 2 relative to adenine results in formation of a third hydrogen bond when interacting with uridine. New method for chemical synthesis of protected 3'-O-phosphoramidite of LNA-2,6-diaminopurine riboside is described. The derivatives of 2'-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides were used to prepare complete 2'-O-methyl RNA and LNA-2'-O-methyl RNA chimeric oligonucleotides to pair with RNA oligonucleotides. Thermodynamic stabilities of these duplexes demonstrated that replacement of a single internal 2'-O-methyladenosine with 2'-O-methyl-2,6-diaminopurine riboside (D(M)) or LNA-2,6-diaminopurine riboside (D(L)) increases the thermodynamic stability (DeltaDeltaG degrees 37) on average by 0.9 and 2.3 kcal/mol, respectively. Moreover, the results fit a nearest neighbor model for predicting duplex stability at 37 degrees C. D-A and D-G but not D-C mismatches formed by D(M) or D(L) generally destabilize 2'-O-methyl RNA/RNA and LNA-2'-O-methyl RNA/RNA duplexes relative to the same type of mismatches formed by 2'-O-methyladenosine and LNA-adenosine, respectively. The enhanced thermodynamic stability of fully complementary duplexes and decreased thermodynamic stability of some mismatched duplexes are useful for many RNA studies, including those involving microarrays.

Prediction of RNA secondary structure by free energy minimization

RNA secondary structure is often predicted from sequence by free energy minimization. Over the past two years, advances have been made in the estimation of folding free energy change, the mapping of secondary structure and the implementation of computer programs for structure prediction. The trends in computer program development are: efficient use of experimental mapping of structures to constrain structure prediction; use of statistical mechanics to improve the fidelity of structure prediction; inclusion of pseudoknots in secondary structure prediction; and use of two or more homologous sequences to find a common structure.