Stereospecificity of Oligonucleotide Interactions Revisited: No Evidence for Heterochiral Hybridization and Ribozyme/DNAzyme Activity

A major challenge for the application of RNA- or DNA-oligonucleotides in biotechnology and molecular medicine is their susceptibility to abundant nucleases. One intriguing possibility to tackle this problem is the use of mirror-image (l-)oligonucleotides. For aptamers, this concept has successfully been applied to even develop therapeutic agents, so-called Spiegelmers. However, for technologies depending on RNA/RNA or RNA/DNA hybridization, like antisense or RNA interference, it has not been possible to use mirror-image oligonucleotides because Watson-Crick base pairing of complementary strands is (thought to be) stereospecific. Many scientists consider this a general principle if not a dogma. A recent publication proposing heterochiral Watson-Crick base pairing and sequence-specific hydrolysis of natural RNA by mirror-image ribozymes or DNAzymes (and vice versa) prompted us to systematically revisit the stereospecificity of oligonucleotides hybridization and catalytic activity. Using hyperchromicity measurements we demonstrate that hybridization only occurs among homochiral anti-parallel complementary oligonucleotide strands. As expected, achiral PNA hybridizes to RNA and DNA irrespective of their chirality. In functional assays we could not confirm an alleged heterochiral hydrolytic activity of ribozymes or DNAzymes. Our results confirm a strict stereospecificity of oligonucleotide hybridization and clearly argue against the possibility to use mirror-image oligonucleotides for gene silencing or antisense applications.     

Oligonucleotides form a duplex with non-helical properties on a positively charged surface

The double helix is known to form as a result of hybridization of complementary nucleic acid strands in aqueous solution. In the helix the negatively charged phosphate groups of each nucleic acid strand are distributed helically on the outside of the duplex and are available for interaction with cationic groups. Cation-coated glass surfaces are now widely used in biotechnology, especially for covalent attachment of cDNAs and oligonucleotides as surface-bound probes on microarrays. These cationic surfaces can bind the nucleic acid backbone electrostatically through the phosphate moiety. Here we describe a simple method to fabricate DNA microarrays based upon adsorptive rather than covalent attachment of oligonucleotides to a positively charged surface. We show that such adsorbed oligonucleotide probes form a densely packed monolayer, which retains capacity for base pair-specific hybridization with a solution state DNA target strand to form the duplex. However, both strand dissociation kinetics and the rate of DNase digestion suggest, on symmetry grounds, that the target DNA binds to such adsorbed oligonucleotides to form a highly asymmetrical and unwound duplex. Thus, it is suggested that, at least on a charged surface, a non-helical DNA duplex can be the preferred structural isomer under standard biochemical conditions.     

Energetic signatures of single base bulges: thermodynamic consequences and biological implications

DNA bulges are biologically consequential defects that can arise from template-primer misalignments during replication and pose challenges to the cellular DNA repair machinery. Calorimetric and spectroscopic characterizations of defect-containing duplexes reveal systematic patterns of sequence-context dependent bulge-induced destabilizations. These distinguishing energetic signatures are manifest in three coupled characteristics, namely: the magnitude of the bulge-induced duplex destabilization (ΔΔG_Bulge); the thermodynamic origins of ΔΔG_Bulge (i.e. enthalpic versus entropic); and, the cooperativity of the duplex melting transition (i.e. two-state versus non-two state). We find moderately destabilized duplexes undergo two-state dissociation and exhibit ΔΔG_Bulge values consistent with localized, nearest neighbor perturbations arising from unfavorable entropic contributions. Conversely, strongly destabilized duplexes melt in a non-two-state manner and exhibit ΔΔG_Bulge values consistent with perturbations exceeding nearest-neighbor expectations that are enthalpic in origin. Significantly, our data reveal an intriguing correlation in which the energetic impact of a single bulge base centered in one strand portends the impact of the corresponding complementary bulge base embedded in the opposite strand. We discuss potential correlations between these bulge-specific differential energetic profiles and their overall biological implications in terms of DNA recognition, repair and replication.     

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2′-O-methyl (2′-OMe) oligonucleotides GCGAAUUCGC, (UA)₆ and (CG)₃ using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Δn_W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Δn_W was 3–4 for RNA duplexes and 2–3 for DNA and 2′-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2′-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2′-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.     

Thermodynamic calculations and statistical correlations for oligo-probes design

Optimization of probe design for array-based experiments requires improved predictability of oligonucleotide hybridization behavior. Currently, designing oligonucleotides capable of interacting efficiently and specifically with the relevant target is not a routine procedure. Multiple examples demonstrate that oligonucleotides targeting different regions of the same RNA differ in their hybridization ability. The present work shows how thermodynamic evaluations of oligo-target duplex or oligo self-structure stabilities can facilitate probe design. Statistical analysis of large sets of hybridization data reveals that thermodynamic evaluation of oligonucleotide properties can be used to avoid poor RNA binders. Thermodynamic criteria for the selection of 20 and 21mers, which, with high probability, interact efficiently and specifically with their targets, are suggested. The design of longer oligonucleotides can also be facilitated by the same calculations of ΔG°_T values for oligo-target duplex or oligo self-structure stabilities and similar selection schemes.     

Hybridization of 2'-ribose modified mixed-sequence oligonucleotides: thermodynamic and kinetic studies

In this study, we characterize the thermodynamics of hybridization, binding kinetics and conformations of four ribose-modified (2′-fluoro, 2′-O-propyl, 2′-O-methoxyethyl and 2′-O-aminopropyl) decameric mixed-sequence oligonucleotides. Hybridization to the complementary non-modified DNA or RNA decamer was probed by fluorescence and circular-dichroism spectroscopy and compared to the same duplex formed between two non-modified strands. The thermal melting points of DNA–DNA duplexes were increased by 1.8, 2.2, 0.3 and 1.3°C for each propyl, methoxyethyl, aminopropyl and fluoro modification, respectively. In the case of DNA–RNA duplexes, the melting points were increased by 3.1, 4.1 and 1.0°C for each propyl, methoxyethyl and aminopropyl modification, respectively. The high stability of the duplexes formed with propyl-, methoxyethyl- and fluoro-modified oligonucleotides correlated with high preorganization in these single-strands. Despite higher thermodynamic duplex stability, hybridization kinetics to complementary DNA or RNA was slower for propyl- and methoxyethyl-modified oligonucleotides than for the non-modified control. In contrast, the positively-charged aminopropyl-modified oligonucleotide showed rapid binding to the complementary DNA or RNA.     

Design considerations for array CGH to oligonucleotide arrays.

BACKGROUND: Representational oligonucleotide microarray analysis has been developed for detection of single nucleotide polymorphisms and/or for genome copy number changes. In this process, the intensity of hybridization to oligonucleotides arrays is increased by hybridizing a polymerase chain reaction (PCR)-amplified representation of reduced genomic complexity. However, hybridization to some oligonucleotides is not sufficiently high to allow precise analysis of that portion of the genome. METHODS: In an effort to identify aspects of oligonucleotide hybridization affecting signal intensity, we explored the importance of the PCR product strand to which each oligonucleotide is homologous and the sequence of the array oligonucleotides. We accomplished this by hybridizing multiple PCR-amplified products to oligonucleotide arrays carrying two sense and two antisense 50-mer oligonucleotides for each PCR amplicon. RESULTS: In some cases, hybridization intensity depended more strongly on the PCR amplicon strand (i.e., sense vs. antisense) than on the detection oligonucleotide sequence. In other cases, the oligonucleotide sequence seemed to dominate. CONCLUSION: Oligonucleotide arrays for analysis of DNA copy number or for single nucleotide polymorphism content should be designed to carry probes to sense and antisense strands of each PCR amplicon to ensure sufficient hybridization and signal intensity.     

First passage time study of DNA strand displacement

DNA strand displacement, in which a single-stranded nucleic acid invades a DNA duplex, is pervasive in genomic processes and DNA engineering applications. The kinetics of strand displacement have been studied in bulk; however, the kinetics of the underlying strand exchange were obfuscated by a slow bimolecular association step. Here, we use a novel single-molecule fluorescence resonance energy transfer approach termed the “fission” assay to obtain the full distribution of first passage times of unimolecular strand displacement. At a frame time of 4.4 ms, the first passage time distribution for a 14-nucleotide displacement domain exhibited a nearly monotonic decay with little delay. Among the eight different sequences we tested, the mean displacement time was on average 35 ms and varied by up to a factor of 13. The measured displacement kinetics also varied between complementary invaders and between RNA and DNA invaders of the same base sequence, except for T → U substitution. However, displacement times were largely insensitive to the monovalent salt concentration in the range of 0.25–1 M. Using a one-dimensional random walk model, we infer that the single-step displacement time is in the range of ～30–300 μs, depending on the base identity. The framework presented here is broadly applicable to the kinetic analysis of multistep processes investigated at the single-molecule level.     

Thermodynamic criteria for high hit rate antisense oligonucleotide design

Antisense oligonucleotides are used for therapeutic applications and in functional genomic studies. In practice, however, many of the oligonucleotides complementary to an mRNA have little or no antisense activity. Theoretical strategies to improve the ‘hit rate’ in antisense screens will reduce the cost of discovery and may lead to identification of antisense oligonucleotides with increased potency. Statistical analysis performed on data collected from more than 1000 experiments with phosphorothioate-modified oligonucleotides revealed that the oligo-probes, which form stable duplexes with RNA (ΔG^o₃₇ ≤ –30 kcal/mol) and have small self-interaction potential, are more frequently efficient than molecules that form less stable oligonucleotide–RNA hybrids or more stable self-structures. To achieve optimal statistical preference, the values for self-interaction should be (ΔG^o₃₇) ≥ –8 kcal/mol for inter-oligonucleotide pairing and (ΔG^o₃₇) ≥ –1.1 kcal/mol for intra-molecular pairing. Selection of oligonucleotides with these thermodynamic values in the analyzed experiments would have increased the ‘hit rate’ by as much as 6-fold.     

Analysis of Damage Tolerance Pathways in Saccharomyces cerevisiae: a Requirement for Rev3 DNA Polymerase in Translesion Synthesis

The replication of double-stranded plasmids containing a single N-2-acetylaminofluorene (AAF) adduct located in a short, heteroduplex sequence was analyzed in Saccharomyces cerevisiae. The strains used were proficient or deficient for the activity of DNA polymerase ζ (REV3 and rev3Δ, respectively) in a mismatch and nucleotide excision repair-defective background (msh2Δ rad10Δ). The plasmid design enabled the determination of the frequency with which translesion synthesis (TLS) and mechanisms avoiding the adduct by using the undamaged, complementary strand (damage avoidance mechanisms) are invoked to complete replication. To this end, a hybridization technique was implemented to probe plasmid DNA isolated from individual yeast transformants by using short, ³²P-end-labeled oligonucleotides specific to each strand of the heteroduplex. In both the REV3 and rev3Δ strains, the two strands of an unmodified heteroduplex plasmid were replicated in ~80% of the transformants, with the remaining 20% having possibly undergone prereplicative MSH2-independent mismatch repair. However, in the presence of the AAF adduct, TLS occurred in only 8% of the REV3 transformants, among which 97% was mostly error free and only 3% resulted in a mutation. All TLS observed in the REV3 strain was abolished in the rev3Δ mutant, providing for the first time in vivo biochemical evidence of a requirement for the Rev3 protein in TLS.     

Formation of an intramolecular triple-stranded DNA structure monitored by fluorescence of 2-aminopurine or 6-methylisoxanthopterin

The parallel (recombination) ‘R-triplex’ can accommodate any nucleotide sequence with the two identical DNA strands in parallel orientation. We have studied oligonucleotides able to fold back into such a recombination-like structure. We show that the fluorescent base analogs 2-aminopurine (2AP) and 6-methylisoxanthopterin (6MI) can be used as structural probes for monitoring the integrity of the triple-stranded conformation and for deriving the thermodynamic characteristics of these structures. A single adenine or guanine base in the third strand of the triplex-forming and the control oligonucleotides, as well as in the double-stranded (ds) and single-stranded (ss) reference molecules, was substituted with 2AP or 6MI. The 2AP*(T·A) and 6MI*(C·G) triplets were monitored by their fluorescence emission and the thermal denaturation curves were analyzed with a quasi-two-state model. The fluorescence of 2AP introduced into an oligonucleotide sequence unable to form a triplex served as a negative control. We observed a remarkable similarity between the thermodynamic parameters derived from melting of the secondary structures monitored through absorption of all bases at 260 nm or from fluorescence of the single base analog. The similarity suggests that fluorescence of the 2AP and 6MI base analogs may be used to monitor the structural disposition of the third strand. We consider the data in the light of alternative ‘branch migration’ and ‘strand exchange’ structures and discuss why these are less likely than the R-type triplex.     

Pestivirus NS3 (p80) protein possesses RNA helicase activity.

The pestivirus bovine viral diarrhea virus (BVDV) p80 protein (referred to here as the NS3 protein) contains amino acid sequence motifs predictive of three enzymatic activities: serine proteinase, nucleoside triphosphatase, and RNA helicase. We have previously demonstrated that the former two enzymatic activities are associated with this protein. Here, we show that a purified recombinant BVDV NS3 protein derived from baculovirus-infected insect cells possesses RNA helicase activity. BVDV NS3 RNA helicase activity was specifically inhibited by monoclonal antibodies to the p80 protein. The activity was dependent on the presence of nucleoside triphosphate and divalent cation, with a preference for ATP and Mn2+. Hydrolysis of the nucleoside triphosphate was necessary for strand displacement. The helicase activity required substrates with an un-base-paired region on the template strand 3' of the duplex region. As few as three un-base-paired nucleotides were sufficient for efficient oligonucleotide displacement. However, the enzyme did not act on substrates having a single-stranded region only to the 5' end of the duplex or on substrates lacking single-stranded regions altogether (blunt-ended duplex substrates), suggesting that the directionality of the BVDV RNA helicase was 3' to 5' with respect to the template strand. The BVDV helicase activity was able to displace both RNA and DNA oligonucleotides from RNA template strands but was unable to release oligonucleotides from DNA templates. The possible role of this activity in pestivirus replication is discussed.     

Thermodynamic and Kinetic Characterization of Antisense Oligodeoxynucleotide Binding to a Structured mRNA

Antisense oligonucleotides act as exogenous inhibitors of gene expression by binding to a complementary sequence on the target mRNA, preventing translation into protein. Antisense technology is being applied successfully as a research tool and as a molecular therapeutic. However, a quantitative understanding of binding energetics between short oligonucleotides and longer mRNA targets is lacking, and selecting a high-affinity antisense oligonucleotide sequence from the many possibilities complementary to a particular RNA is a critical step in designing an effective antisense inhibitor. Here, we report measurements of the thermodynamics and kinetics of hybridization for a number of oligodeoxynucleotides (ODNs) complementary to the rabbit β-globin (RBG) mRNA using a binding assay that facilitates rapid separation of bound from free species in solution. A wide range of equilibrium dissociation constants were observed, and association rate constants within the measurable range correlated strongly with binding affinity. In addition, a significant correlation was observed of measured binding affinities with binding affinity values predicted using a thermodynamic model involving DNA and RNA unfolding, ODN hybridization, and RNA restructuring to a final free energy minimum. In contrast to the behavior observed for hybridization of short strands, the association rate constant increased with temperature, suggesting that the kinetics of association are related to disrupting the native structure of the target RNA. The rate of cleavage of the RBG mRNA in the presence of ribonuclease H and ODNs of varying association kinetics displayed apparent first-order kinetics, with the rate constant exhibiting binding-limited behavior at low association rates and reaction-limited behavior at higher rates. Implications for the rational design of effective antisense reagents are discussed.     

Measurement of steady-state kinetic parameters for DNA unwinding by the bacteriophage T4 Dda helicase: use of peptide nucleic acids to trap single-stranded DNA products of helicase reactions

Measurement of steady-state rates of unwinding of double-stranded oligonucleotides by helicases is hampered due to rapid reannealing of the single-stranded DNA products. Including an oligonucleotide in the reaction mixture which can hybridize with one of the single strands can prevent reannealing. However, helicases bind to single-stranded DNA, therefore the additional oligonucleotide can sequester the enzyme, leading to slower observed rates for unwinding. To circumvent this problem, the oligonucleotide that serves as a trap was replaced with a strand of peptide nucleic acid (PNA). Fluorescence polarization was used to determine that a 15mer PNA strand does not bind to the bacteriophage T4 Dda helicase. Steady-state kinetic parameters of unwinding catalyzed by Dda were determined by using PNA as a trapping strand. The substrate consisted of a partial duplex with 15 nt of single-stranded DNA and 15 bp. In the presence of 250 nM substrate and 1 nM Dda, the rate of unwinding in the presence of the DNA trapping strand was 0.30 nM s^–1 whereas the rate was 1.34 nM s^–1 in the presence of the PNA trapping strand. PNA prevents reannealing of single-stranded DNA products, but does not sequester the helicase. This assay will prove useful in defining the complete kinetic mechanism for unwinding of oligonucleotide substrates by this helicase.     

The synaptic complex of RecA protein participates in hybridization and inverse strand exchange reactions

RecA protein catalyzes strand exchange between homologous single-stranded and double-stranded DNAs. In the presence of ATPgammaS, the post-strand exchange synaptic complex is a stable end product that can be studied. Here we ask whether such complexes can hybridize to or exchange with DNA, 2'-OMe RNA, PNA, or LNA oligonucleotides. Using a gel mobility shift assay, we show that the displaced strand of a 45 bp synaptic complex can hybridize to complementary oligonucleotides with different backbones to form a four-stranded (double D-loop) joint that survives removal of the RecA protein. This hybridization reaction, which confirms the single-stranded character of the displaced strand in a synaptic complex, might initiate recombination-dependent DNA replication if it occurs in vivo. We also show that either strand of the heteroduplex in a 30 bp synaptic complex can be replaced with a homologous DNA oligonucleotide in a strand exchange reaction that is mediated by the RecA filament. Consistent with the important role that deoxyribose plays in strand exchange, oligonucleotides with non-DNA backbones did not participate in this reaction. The hybridization and strand exchange reactions reported here demonstrate that short synaptic complexes are dynamic structures even in the presence of ATPgammaS.     

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.     

A General Framework for Thermodynamically Consistent Parameterization and Efficient Sampling of Enzymatic Reactions

Kinetic models provide the means to understand and predict the dynamic behaviour of enzymes upon different perturbations. Despite their obvious advantages, classical parameterizations require large amounts of data to fit their parameters. Particularly, enzymes displaying complex reaction and regulatory (allosteric) mechanisms require a great number of parameters and are therefore often represented by approximate formulae, thereby facilitating the fitting but ignoring many real kinetic behaviours. Here, we show that full exploration of the plausible kinetic space for any enzyme can be achieved using sampling strategies provided a thermodynamically feasible parameterization is used. To this end, we developed a General Reaction Assembly and Sampling Platform (GRASP) capable of consistently parameterizing and sampling accurate kinetic models using minimal reference data. The former integrates the generalized MWC model and the elementary reaction formalism. By formulating the appropriate thermodynamic constraints, our framework enables parameterization of any oligomeric enzyme kinetics without sacrificing complexity or using simplifying assumptions. This thermodynamically safe parameterization relies on the definition of a reference state upon which feasible parameter sets can be efficiently sampled. Uniform sampling of the kinetics space enabled dissecting enzyme catalysis and revealing the impact of thermodynamics on reaction kinetics. Our analysis distinguished three reaction elasticity regions for common biochemical reactions: a steep linear region (0> ΔG_r >-2 kJ/mol), a transition region (-2> ΔG_r >-20 kJ/mol) and a constant elasticity region (ΔG_r <-20 kJ/mol). We also applied this framework to model more complex kinetic behaviours such as the monomeric cooperativity of the mammalian glucokinase and the ultrasensitive response of the phosphoenolpyruvate carboxylase of Escherichia coli. In both cases, our approach described appropriately not only the kinetic behaviour of these enzymes, but it also provided insights about the particular features underpinning the observed kinetics. Overall, this framework will enable systematic parameterization and sampling of enzymatic reactions.     

Formation of the double helix: a mutational study

To investigate the mechanisms by which oligonucleotides hybridize to target molecules, the binding of two oligodeoxynucleotide probes to RNA targets was measured over a broad range of temperatures. Mutations were then scanned across each DNA/RNA hybrid to map, at single base resolution, sequences important for hybridization. Despite being unrelated in sequence, each hybrid formed by a similar mechanism. In the absence of secondary structure, two stretches of bases, termed nucleation regions, cooperated with one another by a looping mechanism to nucleate hybridization. Mutations inside each nucleation region strongly decreased hybridization rates, even at temperatures well below the melting temperature (T_m) of the hybridized duplex. Surprisingly, nucleation regions were detected in a RNA target but not a corresponding DNA target. When either nucleation region was sequestered in secondary structure, the hybridization rate fell and the mechanism of hybridization changed. Single-stranded bases within the nucleation region of the probe and target first collided to form a double helix. If sufficiently G + C rich, the double helix then propagated throughout the oligonucleotide by a strand invasion process. On the basis of these results, general mechanisms for the hybridization of oligonucleotides to complementary and mutant targets are proposed.     

DNA assembly using bis-peptide nucleic acids (bisPNAs)

DNA nanostructures are ordered oligonucleotide arrangements that have applications for DNA computers, crystallography, diagnostics and material sciences. Peptide nucleic acid (PNA) is a DNA/RNA mimic that offers many advantages for hybridization, but its potential for application in the field of DNA nanotechnology has yet to be thoroughly examined. We report the synthesis and characterization of tethered PNA molecules (bisPNAs) designed to assemble two individual DNA molecules through Watson–Crick base pairing. The spacer regions linking the PNAs were varied in length and contained amino acids with different electrostatic properties. We observed that bisPNAs effectively assembled oligonucleotides that were either the exact length of the PNA or that contained overhanging regions that projected outwards. In contrast, DNA assembly was much less efficient if the oligonucleotides contained overhanging regions that projected inwards. Surprisingly, the length of the spacer region between the PNA sequences did not greatly affect the efficiency of DNA assembly. Reasons for inefficient assembly of inward projecting DNA oligonucleotides include non-sequence-specific intramolecular interactions between the overhanging region of the bisPNA and steric conflicts that complicate simultaneous binding of two inward projecting strands. These results suggest that bisPNA molecules can be used for self-assembling DNA nanostructures provided that the arrangement of the hybridizing DNA oligonucleotides does not interfere with simultaneous hybridization to the bisPNA molecule.