Hybridization kinetics of out-of-equilibrium mixtures of short RNA oligonucleotides

Hybridization and strand displacement kinetics determine the evolution of the base paired configurations of mixtures of oligonucleotides over time. Although much attention has been focused on the thermodynamics of DNA and RNA base pairing in the scientific literature, much less work has been done on the time dependence of interactions involving multiple strands, especially in RNA. Here we provide a study of oligoribonucleotide interaction kinetics and show that it is possible to calculate the association, dissociation and strand displacement rates displayed by short oligonucleotides (5nt–12nt) that exhibit no expected secondary structure as simple functions of oligonucleotide length, CG content, ΔG of hybridization and ΔG of toehold binding. We then show that the resultant calculated kinetic parameters are consistent with the experimentally observed time dependent changes in concentrations of the different species present in mixtures of multiple competing RNA strands. We show that by changing the mixture composition, it is possible to create and tune kinetic traps that extend by orders of magnitude the typical sub-second hybridization timescale of two complementary oligonucleotides. We suggest that the slow equilibration of complex oligonucleotide mixtures may have facilitated the nonenzymatic replication of RNA during the origin of life.     

DNA hairpins: fuel for autonomous DNA devices

We present a study of the hybridization of complementary DNA hairpin loops, with particular reference to their use as fuel for autonomous DNA devices. The rate of spontaneous hybridization between complementary hairpins can be reduced by increasing the neck length or decreasing the loop length. Hairpins with larger loops rapidly form long-lived kissed complexes. Hairpin loops may be opened by strand displacement using an opening strand that contains the same sequence as half of the neck and a "toehold" complementary to a single-stranded domain adjacent to the neck. We find loop opening via an external toehold to be 10-100 times faster than via an internal toehold. We measure rates of loop opening by opening strands that are at least 1000 times faster than the spontaneous interaction between hairpins. We discuss suitable choices for loop, neck, and toehold length for hairpin loops to be used as fuel for autonomous DNA devices.     

Kinetics of heterochiral strand displacement from PNA–DNA heteroduplexes

Dynamic DNA nanodevices represent powerful tools for the interrogation and manipulation of biological systems. Yet, implementation remains challenging due to nuclease degradation and other cellular factors. Use of l-DNA, the nuclease resistant enantiomer of native d-DNA, provides a promising solution. On this basis, we recently developed a strand displacement methodology, referred to as ‘heterochiral’ strand displacement, that enables robust l-DNA nanodevices to be sequence-specifically interfaced with endogenous d-nucleic acids. However, the underlying reaction – strand displacement from PNA–DNA heteroduplexes – remains poorly characterized, limiting design capabilities. Herein, we characterize the kinetics of strand displacement from PNA–DNA heteroduplexes and show that reaction rates can be predictably tuned based on several common design parameters, including toehold length and mismatches. Moreover, we investigate the impact of nucleic acid stereochemistry on reaction kinetics and thermodynamics, revealing important insights into the biophysical mechanisms of heterochiral strand displacement. Importantly, we show that strand displacement from PNA–DNA heteroduplexes is compatible with RNA inputs, the most common nucleic acid target for intracellular applications. Overall, this work greatly improves the understanding of heterochiral strand displacement reactions and will be useful in the rational design and optimization of l-DNA nanodevices that operate at the interface with biology.     

Exquisite allele discrimination by toehold hairpin primers

The ability to detect and monitor single nucleotide polymorphisms (SNPs) in biological samples is an enabling research and clinical tool. We have developed a surprising, inexpensive primer design method that provides exquisite discrimination between SNPs. The field of DNA computation is largely reliant on using so-called toeholds to initiate strand displacement reactions, leading to the execution of kinetically trapped circuits. We have now similarly found that the short toehold sequence to a target of interest can initiate both strand displacement within the hairpin and extension of the primer by a polymerase, both of which will further stabilize the primer:template complex. However, if the short toehold does not bind, neither of these events can readily occur and thus amplification should not occur. Toehold hairpin primers were used to detect drug resistance alleles in two genes, rpoB and katG, in the Mycobacterium tuberculosis genome, and ten alleles in the Escherichia coli genome. During real-time PCR, the primers discriminate between mismatched templates with Cq delays that are frequently so large that the presence or absence of mismatches is essentially a ‘yes/no’ answer.     

On the biophysics and kinetics of toehold-mediated DNA strand displacement

Dynamic DNA nanotechnology often uses toehold-mediated strand displacement for controlling reaction kinetics. Although the dependence of strand displacement kinetics on toehold length has been experimentally characterized and phenomenologically modeled, detailed biophysical understanding has remained elusive. Here, we study strand displacement at multiple levels of detail, using an intuitive model of a random walk on a 1D energy landscape, a secondary structure kinetics model with single base-pair steps and a coarse-grained molecular model that incorporates 3D geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Two factors explain the dependence of strand displacement kinetics on toehold length: (i) the physical process by which a single step of branch migration occurs is significantly slower than the fraying of a single base pair and (ii) initiating branch migration incurs a thermodynamic penalty, not captured by state-of-the-art nearest neighbor models of DNA, due to the additional overhang it engenders at the junction. Our findings are consistent with previously measured or inferred rates for hybridization, fraying and branch migration, and they provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.     

DNA Logic Gate Based on Metallo-Toehold Strand Displacement

DNA is increasingly being used as an ideal material for the construction of nanoscale structures, circuits, and machines. Toehold-mediated DNA strand displacement reactions play a very important role in these enzyme-free constructions. In this study, the concept of metallo-toehold was utilized to further develop a mechanism for strand displacement driven by Ag⁺ ions, in which the intercalation of cytosine–cytosine mismatched base pairs on the toeholds provides additional control by varying of the concentration of Ag⁺ ions. The characteristics of displacement reaction in response to different concentration of Ag⁺ ions are investigated by fluorescence spectral and non-denaturing polyacrylamide gel electrophoresis. The reaction can successfully occur when the concentration of Ag⁺ ions is suitabe; excess Ag⁺ ions block the reaction. Furthermore, the displacement reaction can be tuned and controlled most efficiently under the condition of two C:C mismatched base pairs placed on the six-nt toehold. Based on our research, a mechanism was developed to construct Boolean logic gate AND and OR by employing strand displacement reaction as a tool, Ag⁺ and Hg²⁺ as input.     

Sulfolobus chromatin proteins modulate strand displacement by DNA polymerase B1

Strand displacement by a DNA polymerase serves a key role in Okazaki fragment maturation, which involves displacement of the RNA primer of the preexisting Okazaki fragment into a flap structure, and subsequent flap removal and fragment ligation. We investigated the role of Sulfolobus chromatin proteins Sso7d and Cren7 in strand displacement by DNA polymerase B1 (PolB1) from the hyperthermophilic archaeon Sulfolobus solfataricus. PolB1 showed a robust strand displacement activity and was capable of synthesizing thousands of nucleotides on a DNA-primed 72-nt single-stranded circular DNA template. This activity was inhibited by both Sso7d and Cren7, which limited the flap length to 3–4 nt at saturating concentrations. However, neither protein inhibited RNA displacement on an RNA-primed single-stranded DNA minicircle by PolB1. Strand displacement remained sensitive to modulation by the chromatin proteins when PolB1 was in association with proliferating cell nuclear antigen. Inhibition of DNA instead of RNA strand displacement by the chromatin proteins is consistent with the finding that double-stranded DNA was more efficiently bound and stabilized than an RNA:DNA duplex by these proteins. Our results suggest that Sulfolobus chromatin proteins modulate strand displacement by PolB1, permitting efficient removal of the RNA primer while inhibiting excessive displacement of the newly synthesized DNA strand during Okazaki fragment maturation.     

DNA strand displacement, strand annealing and strand swapping by the Drosophila Bloom's syndrome helicase

Genetic analysis of the Drosophila Bloom's syndrome helicase homolog (mus309/DmBLM) indicates that DmBLM is required for the synthesis-dependent strand annealing (SDSA) pathway of homologous recombination. Here we report the first biochemical study of DmBLM. Recombinant, epitope-tagged DmBLM was expressed in Drosophila cell culture and highly purified protein was prepared from nuclear extracts. Purified DmBLM exists exclusively as a high molecular weight (~1.17 MDa) species, is a DNA-dependent ATPase, has 3′→5′ DNA helicase activity, prefers forked substrate DNAs and anneals complementary DNAs. High-affinity DNA binding is ATP-dependent and low-affinity ATP-independent interactions contribute to forked substrate DNA binding and drive strand annealing. DmBLM combines DNA strand displacement with DNA strand annealing to catalyze the displacement of one DNA strand while annealing a second complementary DNA strand.     

Characterization of the cold stress-induced cyanobacterial DEAD-box protein CrhC as an RNA helicase

We have shown previously that CrhC is a unique member of the DEAD-box family of RNA helicases whose expression occurs specifically under conditions of cold stress. Here we show that recombinant His-tagged CrhC, purified from Escherichia coli, is an ATP-independent RNA binding protein possessing RNA-dependent ATPase activity which is stimulated most efficiently by rRNA and polysome preparations. RNA strand displacement assays indicate that CrhC possesses RNA unwinding activity that is adenosine nucleotide specific. Unwinding of partially duplexed RNA proceeds in the 5′→3′ but not the 3′→5′ direction using standard assay conditions. Immunoprecipitation and far-western analysis indicate that CrhC is a component of a multisubunit complex, interacting specifically with a 37 kDa polypeptide. We propose that CrhC unwinds cold-stabilized secondary structure in the 5′-UTR of RNA during cold stress.     

Study of E. coli Hfq’s RNA annealing acceleration and duplex destabilization activities using substrates with different GC-contents

Folding of RNA molecules into their functional three-dimensional structures is often supported by RNA chaperones, some of which can catalyse the two elementary reactions helix disruption and helix formation. Hfq is one such RNA chaperone, but its strand displacement activity is controversial. Whereas some groups found Hfq to destabilize secondary structures, others did not observe such an activity with their RNA substrates. We studied Hfq’s activities using a set of short RNAs of different thermodynamic stabilities (GC-contents from 4.8% to 61.9%), but constant length. We show that Hfq’s strand displacement as well as its annealing activity are strongly dependent on the substrate’s GC-content. However, this is due to Hfq’s preferred binding of AU-rich sequences and not to the substrate’s thermodynamic stability. Importantly, Hfq catalyses both annealing and strand displacement with comparable rates for different substrates, hinting at RNA strand diffusion and annealing nucleation being rate-limiting for both reactions. Hfq’s strand displacement activity is a result of the thermodynamic destabilization of the RNA through preferred single-strand binding whereas annealing acceleration is independent from Hfq’s thermodynamic influence. Therefore, the two apparently disparate activities annealing acceleration and duplex destabilization are not in energetic conflict with each other.     

Choosing between DNA and RNA: the polymer specificity of RNA helicase NPH-II

NPH-II is a prototypical member of the DExH/D subgroup of superfamily II helicases. It exhibits robust RNA helicase activity, and a detailed kinetic framework for unwinding has been established. However, like most SF2 helicases, there is little known about its mode of substrate recognition and its ability to differentiate between RNA and DNA substrates. Here, we employ a series of chimeric RNA–DNA substrates to explore the molecular determinants for NPH-II specificity on RNA and to determine if there are conditions under which DNA is a substrate. We show that efficient RNA helicase activity depends exclusively on ribose moieties in the loading strand and in a specific section of the 3′-overhang. However, we also document the presence of trace activity on DNA polymers, showing that DNA can be unwound under extremely permissive conditions that favor electrostatic binding. Thus, while polymer-specific SF2 helicases control substrate recognition through specific interactions with the loading strand, alternative specificities can arise under appropriate reaction conditions.     

Regulation of yeast DNA polymerase δ-mediated strand displacement synthesis by 5′-flaps

The strand displacement activity of DNA polymerase δ is strongly stimulated by its interaction with proliferating cell nuclear antigen (PCNA). However, inactivation of the 3′–5′ exonuclease activity is sufficient to allow the polymerase to carry out strand displacement even in the absence of PCNA. We have examined in vitro the basic biochemical properties that allow Pol δ-exo⁻ to carry out strand displacement synthesis and discovered that it is regulated by the 5′-flaps in the DNA strand to be displaced. Under conditions where Pol δ carries out strand displacement synthesis, the presence of long 5′-flaps or addition in trans of ssDNA suppress this activity. This suggests the presence of a secondary DNA binding site on the enzyme that is responsible for modulation of strand displacement activity. The inhibitory effect of a long 5′-flap can be suppressed by its interaction with single-stranded DNA binding proteins. However, this relief of flap-inhibition does not simply originate from binding of Replication Protein A to the flap and sequestering it. Interaction of Pol δ with PCNA eliminates flap-mediated inhibition of strand displacement synthesis by masking the secondary DNA site on the polymerase. These data suggest that in addition to enhancing the processivity of the polymerase PCNA is an allosteric modulator of other Pol δ activities.     

Mobile D-loops are a preferred substrate for the Bloom's syndrome helicase

The Bloom's syndrome helicase, BLM, is a member of the highly conserved RecQ family, and possesses both DNA unwinding and DNA strand annealing activities. BLM also promotes branch migration of Holliday junctions. One role for BLM is to act in conjunction with topoisomerase IIIα to process homologous recombination (HR) intermediates containing a double Holliday junction by a process termed dissolution. However, several lines of evidence suggest that BLM may also act early in one or more of the recombination pathways to eliminate illegitimate or aberrantly paired DNA joint molecules. We have investigated whether BLM can disrupt DNA displacement loops (D-loops), which represent the initial strand invasion step of HR. We show that mobile D-loops created by the RecA recombinase are a highly preferred substrate for BLM with the invading strand being displaced from the duplex. We have identified structural features of the D-loop that determine the efficiency with which BLM promotes D-loop dissociation. We discuss these results in the context of models for the role of BLM as an ‘anti-recombinase’.     

Dengue virus 2 capsid protein chaperones the strand displacement of 5′-3′ cyclization sequences

By virtue of its chaperone activity, the capsid protein of dengue virus strain 2 (DENV2C) promotes nucleic acid structural rearrangements. However, the role of DENV2C during the interaction of RNA elements involved in stabilizing the 5′-3′ panhandle structure of DENV RNA is still unclear. Therefore, we determined how DENV2C affects structural functionality of the capsid-coding region hairpin element (cHP) during annealing and strand displacement of the 9-nt cyclization sequence (5CS) and its complementary 3CS. cHP has two distinct functions: a role in translation start codon selection and a role in RNA synthesis. Our results showed that cHP impedes annealing between 5CS and 3CS. Although DENV2C does not modulate structural functionality of cHP, it accelerates annealing and specifically promotes strand displacement of 3CS during 5′-3′ panhandle formation. Furthermore, DENV2C exerts its chaperone activity by favouring one of the active conformations of cHP. Based on our results, we propose mechanisms for annealing and strand displacement involving cHP. Thus, our results provide mechanistic insights into how DENV2C regulates RNA synthesis by modulating essential RNA elements in the capsid-coding region, that in turn allow for DENV replication.     

Accurate Strand-Specific Quantification of Viral RNA

The presence of full-length complements of viral genomic RNA is a hallmark of RNA virus replication within an infected cell. As such, methods for detecting and measuring specific strands of viral RNA in infected cells and tissues are important in the study of RNA viruses. Strand-specific quantitative real-time PCR (ssqPCR) assays are increasingly being used for this purpose, but the accuracy of these assays depends on the assumption that the amount of cDNA measured during the quantitative PCR (qPCR) step accurately reflects amounts of a specific viral RNA strand present in the RT reaction. To specifically test this assumption, we developed multiple ssqPCR assays for the positive-strand RNA virus o''nyong-nyong (ONNV) that were based upon the most prevalent ssqPCR assay design types in the literature. We then compared various parameters of the ONNV-specific assays. We found that an assay employing standard unmodified virus-specific primers failed to discern the difference between cDNAs generated from virus specific primers and those generated through false priming. Further, we were unable to accurately measure levels of ONNV (−) strand RNA with this assay when higher levels of cDNA generated from the (+) strand were present. Taken together, these results suggest that assays of this type do not accurately quantify levels of the anti-genomic strand present during RNA virus infectious cycles. However, an assay permitting the use of a tag-specific primer was able to distinguish cDNAs transcribed from ONNV (−) strand RNA from other cDNAs present, thus allowing accurate quantification of the anti-genomic strand. We also report the sensitivities of two different detection strategies and chemistries, SYBR® Green and DNA hydrolysis probes, used with our tagged ONNV-specific ssqPCR assays. Finally, we describe development, design and validation of ssqPCR assays for chikungunya virus (CHIKV), the recent cause of large outbreaks of disease in the Indian Ocean region.     

Unwinding of nucleic acids by HCV NS3 helicase is sensitive to the structure of the duplex

Hepatitis C virus (HCV) helicase, non-structural protein 3 (NS3), is proposed to aid in HCV genome replication and is considered a target for inhibition of HCV. In order to investigate the substrate requirements for nucleic acid unwinding by NS3, substrates were prepared by annealing a 30mer oligonucleotide to a 15mer. The resulting 15 bp duplex contained a single-stranded DNA overhang of 15 nt referred to as the bound strand. Other substrates were prepared in which the 15mer DNA was replaced by a strand of peptide nucleic acid (PNA). The PNA–DNA substrate was unwound by NS3, but the observed rate of strand separation was at least 25-fold slower than for the equivalent DNA–DNA substrate. Binding of NS3 to the PNA–DNA substrate was similar to the DNA–DNA substrate, due to the fact that NS3 initially binds to the single-stranded overhang, which was identical in each substrate. A PNA–RNA substrate was not unwound by NS3 under similar conditions. In contrast, morpholino–DNA and phosphorothioate–DNA substrates were utilized as efficiently by NS3 as DNA–DNA substrates. These results indicate that the PNA–DNA and PNA–RNA heteroduplexes adopt structures that are unfavorable for unwinding by NS3, suggesting that the unwinding activity of NS3 is sensitive to the structure of the duplex.