Modelling Toehold-Mediated RNA Strand Displacement

We study the thermodynamics and kinetics of an RNA toehold-mediated strand displacement reaction with a recently developed coarse-grained model of RNA. Strand displacement, during which a single strand displaces a different strand previously bound to a complementary substrate strand, is an essential mechanism in active nucleic acid nanotechnology and has also been hypothesized to occur in vivo. We study the rate of displacement reactions as a function of the length of the toehold and temperature and make two experimentally testable predictions: that the displacement is faster if the toehold is placed at the 5′ end of the substrate; and that the displacement slows down with increasing temperature for longer toeholds.     

Reassociation rate limited displacement of DNA strands by branch migration.

Large branched DNA structures are constructed by two-step reassociation of separated complementary strands from restriction fragments of different lengths. The displacement of DNA strands initially annealed to longer complementary DNA sequences, a process mediated by branch migration, is very rapid and has thus far been followed only under conditions which are second order, DNA reassociation rate limiting. The average lifetime of branched DNA leading to displacement of 1.6 Kb strands is estimated to be less than 10 seconds under conditions of DNA reassociation, Tm-25 degrees C. Several DNA-binding drugs, including intercalating dyes, have been tested to determine their influence, if any, on the kinetics of DNA strand displacements by branch migration. Only actinomycin D was found to have significant effect under the conditions we have described. The kinetics of the strand displacement in the presence of low concentrations of actinomycin D remain second order and slower rate of strand displacement must be attributed to decreased rate of reassociation of DNA strands to form the branched intermediates. Consideration is given to the potential manipulation of DNA structures at site-directed branches and the limitations due to rapid strand displacements. The feasibility of constructing sufficiently large branched DNA regions to approach first order, branch migration rate limiting kinetics is also discussed.     

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.     

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.     

Replication Protein A Stimulates the Werner Syndrome Protein Branch Migration Activity

Loss of the RecQ DNA helicase WRN protein causes Werner syndrome, in which patients exhibit features of premature aging and increased cancer. WRN deficiency induces cellular defects in DNA replication, mitotic homologous recombination (HR), and telomere stability. In addition to DNA unwinding activity, WRN also possesses exonuclease, strand annealing, and branch migration activities. The single strand binding proteins replication protein A (RPA) and telomere-specific POT1 specifically stimulate WRN DNA unwinding activity. To determine whether RPA and POT1 also modulate WRN branch migration activity, we examined biologically relevant mobile D-loops that mimic structures in HR strand invasion and at telomere ends. Both RPA and POT1 block WRN exonuclease digestion of the invading strand by loading on the strand. However, only RPA robustly stimulates WRN branch migration activity and increases the percentage of D-loops that are disrupted. Our results are consistent with cellular data that support RPA enhancement of branch migration during HR repair and, conversely, POT1 limitation of inappropriate recombination and branch migration at telomeric ends. This is, to our knowledge, the first evidence that RPA can stimulate branch migration activity.     

In vitro strand exchange promoted by the herpes simplex virus type-1 single strand DNA-binding protein (ICP8) and DNA helicase-primase

The genome of herpes simplex virus type-1 undergoes a high frequency of homologous recombination in the absence of a virus-encoded RecA-type protein. We hypothesized that viral homologous recombination is mediated by the combined action of the viral single strand DNA-binding protein (ICP8) and helicase-primase. Our results show that ICP8 catalyzes the formation of recombination intermediates (joint molecules) between circular single-stranded acceptor and linear duplex donor DNA. Joint molecules formed by invasion of a 3'-terminal strand displaces the non-complementary 5'-terminal strand, thereby creating a loading site for the helicase-primase. Helicase-primase acts on these joint molecules to promote ATP-dependent branch migration. Finally, we have reconstituted strand exchange by the synchronous action of ICP8 and helicase-primase. Based on these data, we present a recombination mechanism for a eukaryotic DNA virus in which a single strand DNA-binding protein and helicase cooperate to promote homologous pairing and branch migration.     

Effects of the transciption inhibitory protein, TF1, on phage SP01 promoter complex formation and stability.

The uptake of a homologous single-stranded fragment by superhelical DNA produces a complex that contains a stable displacement loop. When the circular DNA was relaxed by the random action of pancreatic DNAase, complexes dissociated by a process which requires that the single-stranded arm of the D-loop be intact. We attribute the dissociation to branch migration, the exchange of like strands at a branch point. The kinetics of dissociation were biphasic. A fraction of the nicked complexes dissociated in a few seconds, the rest dissociated much more slowly. The fraction of molecules that dissociated slowly was directly related to the length of the third strand, and inversely related to temperature. Salt also inhibited dissociation. Under physiological conditions, 37 °C and 0.15 m-NaCl, more than half of complexes containing a third strand of 1000-nucleotide residues survived for at least one minute. These observations provide a guide to handling certain natural or synthetic branched derivatives of DNA. Analyzing our data by the method of Thompson et al. (1976), we have estimated that the time for the exchange of one nucleotide for another at a single-stranded branch is 12 microseconds; but the calculated value depends strongly upon the assumption that single-strand branch migration occurs by a random walk.     

Theoretical analysis of DNA branch migration in the presence of a slow reversible initiation step

Branched DNA structures include several DNA regions connected by three- or four-way DNA junctions. Branched DNAs can be intermediates in DNA replication and recombination in living organisms and in sequence-specific DNA targeting in vitro. Branched DNA structures are usually metastable and irreversibly dissociate to non-branched products via a DNA strand exchange process commonly known as DNA branch migration. The key parameter in the DNA dissociation process is its characteristic time, which depends on the length of the dissociating DNA structure. Here, we predict that the presence of a slow reversible initiation step, which precedes DNA branch migration, can alter, to almost linear dependence, the "classic" quadratic dependence of the dissociation time on the length of the dissociating DNA structure. This prediction can be applied to dissociation of Y-like DNA structures and double D-loop DNA hybrids, which are DNA structures similar to replication bubbles. In addition, the slow initiation step can increase the effect of DNA sequence heterologies within the structure on its kinetic stability. Applications of our analysis for genetic manipulations with branched DNA structures are discussed.     

The structure of replicating adenovirus 2 DNA molecules

Adenovirus 2 (Ad2)-infected KB cells were exposed to a 2.5 min pulse of ³H-thymidine at 19 hr after infection. The labeled DNA molecules were separated from cell DNA and mature Ad2 DNA by sucrose gradient sedimentation and CsCI equilibrium centrifugation under conditions designed to minimize branch migration and hybridization of single strands. Electron microscopy-of fractions containing radioactivity revealed two basic types of putative replicating molecules: Ad2 length duplex DNA molecules with one or more single-stranded branches (type I) and Ad2 length linear DNA molecules with a single-stranded region extending a variable distance from one end (type II). Length measurements, partial denaturation studies and 3′ terminal labeling experiments were consistent with the following model for Ad2 DNA replication. Initiation of DNA synthesis occurs at or near an end of the Ad2 duplex. Following initiation, a daughter strand is synthesized in the 5′ to 3′ direction, displacing the parental strand with the same polarity. This results in the formation of a branched replicating molecule (type I). Initiations at the right and left molecular ends are approximately equal in frequency, and multiple initiations on the same replicating molecule are common. At any given displacement fork in a type I molecule, only one of the two parental strands is replicated. Two nonexclusive mechanisms are proposed to account for the replication of the other parental strand. In some cases, before completion of a round of displacement synthesis initiated at one end of the Ad2 duplex, a second initiation will occur at the opposite end. In these doubly initiated molecules, both parental strands serve as templates for displacement synthesis. Two type II molecules are generated when the oppositely moving displacement forks meet. Alternatively, displacement synthesis may proceed to the end of the Ad2 duplex, resulting in the formation of a daughter duplex and a parental single strand. Replication of the displaced parental strand is then initiated at or near its 3′ terminus, producing a type II molecule. Daughter strand synthesis proceeds in the 5′ to 3′ direction in type II molecules generated by either mechanism, and completion of synthesis results in the formation of a daughter duplex.     

Characterization of strand displacement synthesis catalyzed by bacteriophage T7 DNA polymerase

The DNA polymerase induced after infection of Escherichia coli by bacteriophage T7 can exist in two forms. One distinguishing property of Form I, the elimination of nicks in double-stranded DNA templates, strongly suggests that this form of the polymerase catalyzes limited DNA synthesis at nicks, resulting in displacement of the downstream strand. In this paper, we document this reaction by a detailed characterization of the DNA product. DNA synthesis on circular, duplex DNA templates containing a single site-specific nick results in circular molecules bearing duplex branches. Analysis of newly synthesized DNA excised from the product shows that the majority of the branches are less than 500 base pairs in length and that they arise from a limited number of sites. The branches have fully base-paired termini but are attached by two noncomplementary DNA strands that have a combined length of less than 30 nucleotides. The product molecules are topologically constrained as a result of the duplex branch. DNA sequence analysis has provided an unequivocal structure of one such product molecule. We conclude that strand displacement synthesis catalyzed by Form I of T7 DNA polymerase is terminated by a template-switching reaction. We propose two distinct models for template-switching that we call primer relocation and rotational strand exchange. Strand displacement synthesis catalyzed by Form I of T7 DNA polymerase effectively converts T7 DNA circles that are held together by hydrogen bonds in their 160-nucleotide-long terminal redundancy to T7-length linear molecules. We suggest that strand displacement synthesis catalyzed by T7 DNA polymerase is essential in vivo to the processing of a T7 DNA concatemer to mature T7 genomes.     

Preparation and melting of single strand circular DNA loops.

A method for preparation of single strand DNA circles of almost arbitrary sequence is described. By ligating two sticky ended hairpins together a linear duplex is formed, closed at both ends by single stranded loops. The melting characteristics of such loops are investigated using optical absorbance and NMR. It is shown by comparison with the corresponding linear sequence (closed circle minus the end loops) that the effects of end fraying and the strand concentration dependence of the melting temperature are eliminated in the circular form. Over the concentration range examined (0.5 to 2.0 micromolar strands), the circular DNA has a monophasic melting curve, while the linear duplex is biphasic, probably due to hairpin formation. Since effects of duplex to single strands dissociation do not contribute to melting of the circular molecules (dumbells), these DNAs present a realistic experimental model for examining local thermal stability in DNA.     

Branch migration of three-strand recombination intermediates by RecG, a possible pathway for securing exchanges initiated by 3''-tailed duplex DNA.

RecG protein is required for normal levels of recombination and DNA repair in Escherichia coli. This 76 kDa polypeptide is a junction-specific DNA helicase that acts post-synaptically to drive branch migration of Holliday junction intermediates made by RecA during the strand exchange stage of recombination. To gain further insight into the role of RecG, we studied its activity on three-strand intermediates formed by RecA between circular single-stranded and linear duplex DNAs. Once RecA is removed, RecG drives branch migration of these intermediates by a junction-targeted activity that depends on hydrolysis of ATP. RuvAB has a similar activity. However, when RecG is added to a RecA strand exchange reaction it severely reduces the accumulation of joint molecule intermediates by driving branch migration of junctions in the reverse direction to that catalysed by RecA strand exchange. In comparison, RuvAB has little effect on the reaction. We discuss how reverse branch migration by RecG, which acts counter of the 5'-->3' polarity of RecA binding and strand exchange, could serve to promote or abort the early stages of recombination, depending on the orientation of the single DNA strand initiating the exchange relative to the adjacent duplex region.     

A euryarchaeal histone modulates strand displacement synthesis by replicative DNA polymerases

Euryarchaeota and Crenarchaeota, the two main lineages of the domain Archaea, encode different chromatin proteins and differ in the use of replicative DNA polymerases. Crenarchaea possess a single family B DNA polymerase (PolB), which is capable of strand displacement modulated by the chromatin proteins Cren7 and Sul7d. Euryarchaea have two distinct replicative DNA polymerases, PolB and PolD, a family D DNA polymerase. Here we characterized the strand displacement activities of PolB and PolD from the hyperthermophilic euryarchaeon Pyrococcus furiosus and investigated the influence of HPfA1, a homolog of eukaryotic histones from P. furiosus, on these activities. We showed that both PolB and PolD were efficient in strand displacement. HPfA1 inhibited DNA strand displacement by both DNA polymerases but exhibited little effect on the displacement of a RNA strand annealed to single-stranded template DNA. This is consistent with the finding that HPfA1 bound more tightly to double-stranded DNA than to a RNA:DNA hybrid. Our results suggest that, although crenarchaea and euryarchaea differ in chromosomal packaging, they share similar mechanisms in modulating strand displacement by DNA polymerases during lagging strand DNA synthesis.     

Branch migration of Holliday junctions: identification of RecG protein as a junction specific DNA helicase.

The product of the recG gene of Escherichia coli is needed for normal recombination and DNA repair in E. coli and has been shown to help process Holliday junction intermediates to mature products by catalysing branch migration. The 76 kDa RecG protein contains sequence motifs conserved in the DExH family of helicases, suggesting that it promotes branch migration by unwinding DNA. We show that RecG does not unwind blunt ended duplex DNA or forked duplexes with short unpaired single-strand ends. It also fails to unwind a partial duplex (52 bp) classical helicase substrate containing a short oligonucleotide annealed to circular single-stranded DNA. However, unwinding activity is detected when the duplex region is reduced to 26 bp or less, although this requires high levels of protein. The unwinding proceeds with a clear 3' to 5' polarity with respect to the single strand bound by RecG. Substantially higher levels of unwinding are observed with substrates containing a three-way duplex branch. This is attributed to RecG's particular affinity for junction DNA which we demonstrate would be heightened by single-stranded DNA binding protein in vivo. Reaction requirements for unwinding are the same as for branch migration of Holliday junctions, with a strict dependence on hydrolysis of ATP. These results define RecG as a new class of helicase that has evolved to catalyse the branch migration of Holliday junctions.     

Rad51 uses one mechanism to drive DNA strand exchange in both directions

The Rad51 protein of Saccharomyces cerevisiae, like its bacterial counterpart RecA, promotes strand exchange between circular single-stranded DNA (ssDNA) and linear double-stranded DNA (dsDNA) in vitro. However, the two proteins differ in the requirement for initiating joint molecules and in the polarity of branch migration. Whereas RecA initiates joint molecules from any type of ends on the dsDNA and branch migration proceeds exclusively in the 5'- to 3'-direction with respect to the single strand DNA substrate, initiation mediated by Rad51 requires a complementary 3' or 5' overhanging end of the linear dsDNA and branch migration proceeds in either direction. Here we report that the rates of Rad51-mediated branch migration in either the 5'- to 3'- or 3'- to 5'-directions are affected to the same extent by temperature and MgCl(2). Furthermore, branch migration in both directions is equally impeded by insertions of non-homologous sequences in the dsDNA, inserts of 6 base pairs or more being completely inhibitory. We have also found that the preference of strand exchange in the 5'- to 3'-direction does not change if RPA is replaced by Escherichia coli SSB or T4 gene 32 proteins, suggesting that the preference for the direction of strand exchange is intrinsic to Rad51. Based on these results, we conclude that Rad51-promoted branch migration in either direction occurs fundamentally by the same mechanism, quite probably by stabilizing successively formed heteroduplex base pair.     

Inhibition of recA-mediated strand exchange by adducts of azacytosine-containing DNA and the EcoRII methylase

Wild type Escherichia coli cells containing elevated levels of DNA (cytosine-5)methyltransferases have increased sensitivity to the toxic effects of 5-azacytidine. The methyltransferases form tight binding complexes with azacytosine in DNA which could interfere with the recA recBCD repair pathway which is largely responsible for cell survival after treatment with the drug. We therefore determined if these complexes interfered with recA-mediated strand exchange in vitro. 32P-Labeled DNA fragments containing a single EcoRII site, with cytosine in the (-) strand replaced by 5-azacytosine, were prepared. We investigated the effect of the EcoRII methyltransferase on recA-mediated strand exchange with homologous M13 DNA by electrophoresis in agarose gels. In the absence of the methylase the rate and extent of strand exchange of azacytosine-containing DNA is the same as control DNA. In the presence of the methyltransferase strand exchange is inhibited, but some incorporation of duplexes into recA-single-stranded DNA (ssDNA) complexes still occurs. The formation of these complexes is dependent on the length of the fragment 3' to the methylase binding site on the strand complementary to the ssDNA. The greater the length the greater the number of complexes that form. S-Adenosyl-L-methionine, which enhances binding of the methyltransferase to azacytosine-containing DNA, causes an increase in the inhibition of strand exchange and an increase in the number of inactive complexes formed. The complexes can be dissociated with guanidinium chloride which denatures the methyltransferase and leads to release of the (+) strand. The (-) strand remains associated with the ssDNA. This result implies that a plectonemic joint is formed between recA-ssDNA complexes and azacytosine-containing DNA-methyltransferase complexes. However, branch migration in these complexes is inhibited. Denaturation of the methyltransferase allows branch migration to proceed to completion, releasing the (+) strand.     

A PAM-free CRISPR/Cas12a ultra-specific activation mode based on toehold-mediated strand displacement and branch migration

CRISPR (clustered regularly interspaced short palindromic repeats) technology has achieved great breakthroughs in terms of convenience and sensitivity; it is becoming the most promising molecular tool. However, only two CRISPR activation modes (single and double stranded) are available, and they have specificity and universality bottlenecks that limit the application of CRISPR technology in high-precision molecular recognition. Herein, we proposed a novel CRISPR/Cas12a unrestricted activation mode to greatly improve its performance. The new mode totally eliminates the need for a protospacer adjacent motif and accurately activates Cas12a through toehold-mediated strand displacement and branch migration, which is highly universal and ultra-specific. With this mode, we discriminated all mismatch types and detected the EGFR T790M and L858R mutations in very low abundance. Taken together, our activation mode is deeply incorporated with DNA nanotechnology and extensively broadens the application boundaries of CRISPR technology in biomedical and molecular reaction networks.     

Modeling translocation dynamics of strand displacement DNA synthesis by DNA polymerase I

A model is presented for the translocation dynamics of the strand displacement DNA synthesis by DNA polymerases such as polymerase I family. (i) The model gives an explanation to the experimental results which showed that the rate of strand displacement DNA synthesis is nearly consistent with that of single stranded primer extension synthesis, although the two are expected to have substantial differences in their energetics. (ii) During strand displacement DNA synthesis, the pausing at the specific sequence is considered to be due to an affinity of the fingers subdomain for the specific sequence of dsDNA downstream of the single strand. The theoretical results on the sequence-dependent pausing dynamics such as the mean pausing lifetimes and the distribution of the pausing lifetime are consistent with the experimental data. Moreover, predicted results are presented for the binding affinity of the fingers subdomain for the specific sequence of dsDNA and the dependence of the mean sequence-dependent pausing lifetime on the external force acting on the polymerase.