OligoDesign: Optimal design of LNA (locked nucleic acid) oligonucleotide capture probes for gene expression profiling

We report the development of new software, OligoDesign, which provides optimal design of LNA (locked nucleic acid) substituted oligonucleotides for functional genomics applications. LNAs constitute a novel class of bicyclic RNA analogs having an exceptionally high affinity and specificity toward their complementary DNA and RNA target molecules. The OligoDesign software features recognition and filtering of the target sequence by genome-wide BLAST analysis in order to minimize cross-hybridization with non-target sequences. Furthermore it includes routines for prediction of melting temperature, self-annealing and secondary structure for LNA substituted oligonucleotides, as well as secondary structure prediction of the target nucleotide sequence. Individual scores for all these properties are calculated for each possible LNA oligonucleotide in the query gene and the OligoDesign program ranks the LNA capture probes according to a combined fuzzy logic score and finally returns the top scoring probes to the user in the output. We have successfully used the OligoDesign tool to design a Caenorhabditis elegans LNA oligonucleotide microarray, which allows monitoring of the expression of a set of 120 potential marker genes for a variety of stress and toxicological processes and toxicologically relevant pathways. The OligoDesign program is freely accessible at http://lnatools.com/.     

Characterization of an RNA bulge structure by Fourier transform infrared spectroscopy

There may be several advantages associated with an antisense oligonucleotide that induces a bulged structure into its RNA target molecule. Many structures of RNA bulges are elucidated from single-stranded RNA models. However, a two-component system is the minimum requirement for a realistic antisense model. We have used Fourier transform infrared spectroscopy to investigate a single-stranded RNA oligonucleotide with known NMR solution structure, constructed to model a five nucleotide bulge, and its two-component oligonucleotide counterpart. The infrared spectra show A-helical base-paired stems and non-base-paired loops in both systems. The nucleosides are mainly in an anti-conformation. Both N-type and S-type of sugar puckers can be inferred from the infrared region sensitive to sugar conformations. The S-type of sugar pucker is likely to be associated with the nucleotides in the bulge. The FTIR results display an overall structural similarity between the two model systems.     

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.     

A shell program for the design of PCR primers using genetics computer group (GCG) software (7.1) on VAX/VMS systems.

A number of software analysis packages for the design of PCR primers are available for PCs; however, software for users that depend on VAX/VMS operating systems is not available. By treating oligonucleotides as RNA molecules, I have designed an alternative means toward studying oligonucleotide interactions using software that is currently available from The Genetics Computer Group (GCG, Madison, WI). The oligonucleotide interactions with self and non-self are analyzed by the GCG FOLD program, a program which finds a secondary structure of minimum free energy for an RNA molecule. This approach allows the identification of self-priming primer pairs, and the interaction energies provide a guideline for the prediction of optimal PCR primers.     

Some Properties of Site-Specific Nickase BspD6I and the Possibility of Its Use in Hybridization Analysis of DNA

A new method for hybridization analysis of nucleic acids is proposed on the basis of the ability of site-specific nickases to cleave only one DNA strand. The method is based on the use of a labeled oligonucleotide with the recognition site of the nickase hybridized with the target (DNA or RNA) at an optimal temperature of the enzyme (55°C). The two shorter oligonucleotides formed after the cleavage with the nickase do not complex with the target. Thus, a multiple cleavage of the labeled oligonucleotide takes place on one target molecule. The cleavage of the nucleotide is recorded either by polyacrylamide gel electrophoresis (when a radioactive labeled oligonucleotide is used) or by fluorescence measurements (if the oligonucleotide has the structure of a molecular beacon). The new method was tested on nickase BspD6I and a radioactive oligonucleotide complementary to the polylinker region of the viral DNA strand in bacteriophage M13mp19. Unfortunately, nickase BspD6I does not cleave DNA in the RNA–DNA duplexes and therefore cannot be used for detection of RNA targets.     

The effect of structure in a long target RNA on ribozyme cleavage efficiency.

Inhibition of gene expression by catalytic RNA (ribozymes) requires that ribozymes efficiently cleave specific sites within large target RNAs. However, the cleavage of long target RNAs by ribozymes is much less efficient than cleavage of short oligonucleotide substrates because of higher order structure in the long target RNA. To further study the effects of long target RNA structure on ribozyme cleavage efficiency, we determined the accessibility of seven hammerhead ribozyme cleavage sites in a target RNA that contained human immunodeficiency virus type 1 (HIV-1) vif - vpr . The base pairing-availability of individual nucleotides at each cleavage site was then assessed by chemical modification mapping. The ability of hammerhead ribozymes to cleave the long target RNA was most strongly correlated with the availability of nucleotides near the cleavage site for base pairing with the ribozyme. Moreover, the accessibility of the seven hammerhead ribozyme cleavage sites in the long target RNA varied by up to 400-fold but was directly determined by the availability of cleavage sites for base pairing with the ribozyme. It is therefore unlikely that steric interference affected hammerhead ribozyme cleavage. Chemical modification mapping of cleavage site structure may therefore provide a means to identify efficient hammerhead ribozyme cleavage sites in long target RNAs.     

STUDIES IN OLIGONUCLEOTIDE-BASED ARTIFICIAL NUCLEASE SYSTEMS. INTRAMOLECULAR COPPER (II) COMPLEX FORMATION IN AN OLIGONUCLEOTIDE BIS-PHENANTHROLINE CONJUGATE

We have recently developed oligonucleotide based artificial nuclease (OBAN) systems based on 2′-O-methyloligoribonucleotides carrying a 2,9-dimethylphenanthroline · Zn(II) complex. These hybridize to an RNA molecule with bulge formation in the central region of the target and cleave the RNA target in a catalytic manner. When studying an 11-mer 2′-O-methyloligoribonucleotide carrying two 2,9-dimethylphenanthroline moieties, located 5 base pairs apart from each other, we found that this forms a cyclic structure in the presence of Cu²⁺ ions. This is due to intramolecular Cu(2,9-dimethylphenanthroline)₂ complex formation, i.e., with the two ligands conjugated to the same oligonucleotide.     

Specificity of site directed psoralen addition to RNA.

We describe the attachment of a psoralen derivative (site specific psoralen, SSP) to the 5' end of a DNA oligonucleotide and the hybridization and the photoreaction of this reagent with a complementary target site on an RNA molecule. SSP was coupled to a variety of DNA oligonucleotides to investigate the structural requirements for addition to the RNA. Efficient SSP photoadducts were made on specific uridines by designing an intercalation site at an unpaired nucleotide in the RNA strand within the heteroduplex region. The optimal location for this site was five nucleotides from the oligonucleotide 5' end and just 5' to the target uridine residue. Because the attachment of the SSP to the oligonucleotide is through a disulfide bond, the DNA oligonucleotide can be removed with reduction to leave SSP attached to the RNA strand. The SSP adduct made in this way will be useful for subsequent biochemical and biophysical experiments.     

OligoMatcher

OligoMatcher is a web-based tool for analysis and selection of unique oligonucleotide sequences for gene silencing by antisense oligonucleotides (ASOs) or small interfering RNA (siRNA). A specific BLAST server was built for analysing sequences of ASOs that target pre-mRNA in the cell nucleus. Tissue- and cell-specific expression data of potential cross-reactive genes are integrated in the OligoMatcher program, which allows biologists to select unique oligonucleotide sequences for their target genes in specific experimental systems. AVAILABILITY: The OligoMatcher web server is available at http://shelob.cs.iupui.edu:18081/oligomatch.php. The source code is freely available for non-profit use on request to the authors. CONTACT: Mathew Palakal (mpalakal@cs.iupui.edu) or Shuyu Li (li_shuyu_dan@lilly.com).     

Mapping of accessible sites for oligonucleotide hybridization on hepatitis delta virus ribozymes

Semi-random libraries of DNA 6mers and RNase H digestion were applied to search for sites accessible to hybridization on the genomic and antigenomic HDV ribozymes and their 3′ truncated derivatives. An approach was proposed to correlate the cleavage sites and most likely sequences of oligomers, members of the oligonucleotide libraries, which were engaged in the formation of RNA–DNA hybrids. The predicted positions of oligomers hybridizing to the genomic ribozyme were compared with the fold of polynucleotide chain in the ribozyme crystal structure. The data exemplified the crucial role of target RNA structural features in the binding of antisense oligonucleotides. It turned out that cleavages were induced if the bound oligomer could adapt an ordered helical conformation even when it required partial penetration of an adjacent double-stranded region. The major features of RNA structure disfavoring hybridization and/or RNase H hydrolysis were sharp turns of the polynucleotide chain and breaks in stacking interactions of bases. Based on the predicted positions of oligomers hybridizing to the antigenomic ribozyme we chose and synthesized four antisense DNA 6mers which were shown to direct hydrolysis in the desired, earlier predicted regions of the molecule.     

Construction of new ribozymes requiring short regulator oligonucleotides as a cofactor

A hairpin loop and an oligonucleotide bound to the loop form one-half of the pseudoknot structure. We have designed an allosteric hammerhead ribozyme, which is activated by the introduction of this motif by using a short complementary oligonucleotide as a cofactor. Stem II of the hammerhead ribozyme was substituted with a non-self-complementary loop sequence (loop II) to abolish the cleavage activity. The new ribozyme had almost no cleavage activity of the target RNA. However, it exhibited the cleavage activity in the presence of a cofactor oligoribonucleotide, which is complementary to loop II of the ribozyme. The activity is assumed to be derived from the formation of a pseudo-stem structure between the cofactor oligonucleotide and loop II. The structure including the loop may be similar to the pseudo-half-knot structure. The activation efficiencies of the cofactor oligonucleotides were decreased as the lengths of the oligonucleotides increased, and the ribozyme with a longer loop II was more active than that with a short loop II. Oligoribonucleotides with 3'-dangling purine bases served as efficient cofactors of the ribozyme, and a 2'-O-methyloligonucleotide enhanced the cleavage activity of the ribozyme most efficiently, by as much as about 750-fold as compared with that in the absence of the oligonucleotide. Cofactor oligonucleotides with a cytidine base at the 3'-end also activated a ribozyme with the G10.1.G11.1 mutation, which eliminates the cleavage activity in the wild-type. The binding sites of the oligonucleotide were identified by photo-crosslinking experiments and were found to be the predicted sites in the loop. This is the first report of a design aimed at positively controlling the activity of ribozymes by employing a structural motif. This method can be applied to control the activities of other functional RNAs with hairpin loops.     

Predicting Helical Topologies in RNA Junctions as Tree Graphs

RNA molecules are important cellular components involved in many fundamental biological processes. Understanding the mechanisms behind their functions requires knowledge of their tertiary structures. Though computational RNA folding approaches exist, they often require manual manipulation and expert intuition; predicting global long-range tertiary contacts remains challenging. Here we develop a computational approach and associated program module (RNAJAG) to predict helical arrangements/topologies in RNA junctions. Our method has two components: junction topology prediction and graph modeling. First, junction topologies are determined by a data mining approach from a given secondary structure of the target RNAs; second, the predicted topology is used to construct a tree graph consistent with geometric preferences analyzed from solved RNAs. The predicted graphs, which model the helical arrangements of RNA junctions for a large set of 200 junctions using a cross validation procedure, yield fairly good representations compared to the helical configurations in native RNAs, and can be further used to develop all-atom models as we show for two examples. Because junctions are among the most complex structural elements in RNA, this work advances folding structure prediction methods of large RNAs. The RNAJAG module is available to academic users upon request.     

Efficient reduction of target RNAs by small interfering RNA and RNase H-dependent antisense agents. A comparative analysis

RNA interference can be considered as an antisense mechanism of action that utilizes a double-stranded RNase to promote hydrolysis of the target RNA. We have performed a comparative study of optimized antisense oligonucleotides designed to work by an RNA interference mechanism to oligonucleotides designed to work by an RNase H-dependent mechanism in human cells. The potency, maximal effectiveness, duration of action, and sequence specificity of optimized RNase H-dependent oligonucleotides and small interfering RNA (siRNA) oligonucleotide duplexes were evaluated and found to be comparable. Effects of base mismatches on activity were determined to be position-dependent for both siRNA oligonucleotides and RNase H-dependent oligonucleotides. In addition, we determined that the activity of both siRNA oligonucleotides and RNase H-dependent oligonucleotides is affected by the secondary structure of the target mRNA. To determine whether positions on target RNA identified as being susceptible for RNase H-mediated degradation would be coincident with siRNA target sites, we evaluated the effectiveness of siRNAs designed to bind the same position on the target mRNA as RNase H-dependent oligonucleotides. Examination of 80 siRNA oligonucleotide duplexes designed to bind to RNA from four distinct human genes revealed that, in general, activity correlated with the activity to RNase H-dependent oligonucleotides designed to the same site, although some exceptions were noted. The one major difference between the two strategies is that RNase H-dependent oligonucleotides were determined to be active when directed against targets in the pre-mRNA, whereas siRNAs were not. These results demonstrate that siRNA oligonucleotide- and RNase H-dependent antisense strategies are both valid strategies for evaluating function of genes in cell-based assays.     

No evidence that mRNAs have lower folding free energies than random sequences with the same dinucleotide distribution

This work investigates whether mRNA has a lower estimated folding free energy than random sequences. The free energy estimates are calculated by the mfold program for prediction of RNA secondary structures. For a set of 46 mRNAs it is shown that the predicted free energy is not significantly different from random sequences with the same dinucleotide distribution. For random sequences with the same mononucleotide distribution it has previously been shown that the native mRNA sequences have a lower predicted free energy, which indicates a more stable structure than random sequences. However, dinucleotide content is important when assessing the significance of predicted free energy as the physical stability of RNA secondary structure is known to depend on dinucleotide base stacking energies. Even known RNA secondary structures, like tRNAs, can be shown to have predicted free energies indistinguishable from randomized sequences. This suggests that the predicted free energy is not always a good determinant for RNA folding.     

ECHO-LNA conjugates: hybridization-sensitive fluorescence and its application to fluorescent detection of various RNA strands

Hybridization-sensitive fluorescent DNA probes containing the nucleotide units of locked nucleic acid (LNA) have been developed. Exciton-controlled hybridization-sensitive fluorescent oligonucleotide (ECHO) probes that incorporated LNA nucleotides achieved high thermostability of the hybrid with target RNA strands. The appropriately designed ECHO-LNA chimeric probes exhibited an effective on-off switching property of fluorescence depending on hybridization with RNA and facilitated fluorescent detection of the TAR RNA strand forming a hairpin structure and distinction of one base difference in PLAC4 RNA sequence.