DNA hybridization kinetics: zippering, internal displacement and sequence dependence

DNA as a perfect quantum computer based on the quantum physics principles

DNA is a complex multi-resolution molecule whose theoretical study is a challenge. Its intrinsic multiscale nature requires chemistry and quantum physics to understand the structure and quantum informatics to explain its operation as a perfect quantum computer. Here, we present theoretical results of DNA that allow a better description of its structure and the operation process in the transmission, coding, and decoding of genetic information. Aromaticity is explained by the oscillatory resonant quantum state of correlated electron and hole pairs due to the quantized molecular vibrational energy acting as an attractive force. The correlated pairs form a supercurrent in the nitrogenous bases in a single band π -molecular orbital ( π -MO). The MO wave function ( Φ ) is assumed to be the linear combination of the n constituent atomic orbitals. The central Hydrogen bond between Adenine (A) and Thymine (T) or Guanine (G) and Cytosine (C) functions like an ideal Josephson Junction. The approach of a Josephson Effect between two superconductors is correctly described, as well as the condensation of the nitrogenous bases to obtain the two entangled quantum states that form the qubit. Combining the quantum state of the composite system with the classical information, RNA polymerase teleports one of the four Bell states. DNA is a perfect quantum computer.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Single-Molecule FRET at 10 MHz Count Rates

A bottleneck in many studies utilizing single-molecule Förster resonance energy transfer is the attainable photon count rate, as it determines the temporal resolution of the experiment. As many biologically relevant processes occur on time scales that are hardly accessible with currently achievable photon count rates, there has been considerable effort to find strategies to increase the stability and brightness of fluorescent dyes. Here, we use DNA nanoantennas to drastically increase the achievable photon count rates and observe fast biomolecular dynamics in the small volume between two plasmonic nanoparticles. As a proof of concept, we observe the coupled folding and binding of two intrinsically disordered proteins, which form transient encounter complexes with lifetimes on the order of 100 μs. To test the limits of our approach, we also investigated the hybridization of a short single-stranded DNA to its complementary counterpart, revealing a transition path time of 17 μs at photon count rates of around 10 MHz, which is an order-of-magnitude improvement compared to the state of the art. Concomitantly, the photostability was increased, enabling many seconds long megahertz fluorescence time traces. Due to the modular nature of the DNA origami method, this platform can be adapted to a broad range of biomolecules, providing a promising approach to study previously unobservable ultrafast biophysical processes.

Rational Design of DNA Hydrogels Based on Molecular Dynamics of Polymers

In recent years, DNA has emerged as a fascinating building material to engineer hydrogel due to its excellent programmability, which has gained considerable attention in biomedical applications. Understanding the structure-property relationship and underlying molecular determinants of DNA hydrogel is essential to precisely tailor its macroscopic properties at molecular level. In this review, the rational design principles of DNA molecular networks based on molecular dynamics of polymers on the temporal scale, which can be engineered via the backbone rigidity and crosslinking kinetics, are highlighted. By elucidating the underlying molecular mechanisms and theories, it is aimed to provide a comprehensive overview of how the tunable DNA backbone rigidity and the crosslinking kinetics lead to desirable macroscopic properties of DNA hydrogels, including mechanical properties, diffusive permeability, swelling behaviors, and dynamic features. Furthermore, it is also discussed how the tunable macroscopic properties make DNA hydrogels promising candidates for biomedical applications, such as cell culture, tissue engineering, bio-sensing, and drug delivery.

A rhythmically pulsing leaf-spring DNA-origami nanoengine that drives a passive follower

Molecular engineering seeks to create functional entities for modular use in the bottom-up design of nanoassemblies that can perform complex tasks. Such systems require fuel-consuming nanomotors that can actively drive downstream passive followers. Most artificial molecular motors are driven by Brownian motion, in which, with few exceptions, the generated forces are non-directed and insufficient for efficient transfer to passive second-level components. Consequently, efficient chemical-fuel-driven nanoscale driver-follower systems have not yet been realized. Here we present a DNA nanomachine (70 nm × 70 nm × 12 nm) driven by the chemical energy of DNA-templated RNA-transcription-consuming nucleoside triphosphates as fuel to generate a rhythmic pulsating motion of two rigid DNA-origami arms. Furthermore, we demonstrate actuation control and the simple coupling of the active nanomachine with a passive follower, to which it then transmits its motion, forming a true driver-follower pair.

Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.

Correlated Hybrid DNA Structures Explored by the oxDNA Model

Thermodynamically, perfect DNA hybridization can be formed between probes and their corresponding targets due to the favorable energy. However, this is not the case dynamically. Here, we use molecular dynamics (MD) simulations based on the oxDNA model to investigate the process of DNA microarray hybridization. In general, correlated hybrid DNA structures are formed, including one probe associated with several targets as well as one target hybrid with multiple probes leading to the target-mediated hybridization. The formation of these two types of correlated structures largely depends on the surface coverage of the DNA microarray. Moreover, DNA sequence, DNA length, and spacer length have an impact on the structural formation. Our findings shed light on the dynamics of DNA hybridization, which is important for the application of DNA microarray.

Peptide Nucleic Acids Containing Cationic/Amino-Alkyl Modified Bases Promote Enhanced Hybridization Kinetics and Thermodynamics with Single-Strand DNA

Peptide nucleic acids (PNAs) are antisense molecules with excellent polynucleotide hybridization properties; they are resistant to nuclease degradation but often have poor cell permeability leading to moderate cellular activity and limited clinical results. The addition of cationic substitutions (positive charges) to PNA molecules greatly increases cell permeability. In this report, we describe the synthesis and polynucleotide hybridization properties of a novel cationic/amino-alkyl nucleotide base-modified PNA (OPNA). This study was designed to quantitate the effect the cationic/amino-alkyl nucleotide base modification had on the kinetic and thermodynamic properties of OPNA-DNA hybridization using surface plasmon resonance and UV thermal melt studies. Kinetic studies reveal a favorable 10-30 fold increase in affinity for a single cationic modification on the base of an adenine, cytosine, or guanidine OPNA sequence compared to the nonmodified PNA strand. The increase in affinity is correlated directly with a favorable decrease in the dissociation rate constant and increase in the association rate constant. Introducing additional amino-alkyl base modifications further favors a decrease in the dissociation rate (3-10-fold per amino-alkyl). The thermodynamics driving the OPNA hybridization is promoted by an additional favorable -80 kJ/mol enthalpy of binding for a single amino-alkyl modification compared to the PNA strand. This increase in enthalpy is consistent with an ion-ion interaction with the DNA strand. These kinetic and thermodynamic hybridization studies reveal for the first time that this type of cationic/amino-alkyl base-modified PNA has favorable hybridization properties suitable for development as an antisense oligomer.

Unraveling the Kinetics of the 10-23 RNA-Cleaving DNAzyme

DNA-based enzymes, or DNAzymes, are single-stranded DNA sequences with the ability to catalyze various chemical reactions, including the cleavage of the bond between two RNA nucleotides. Lately, an increasing interest has been observed in these RNA-cleaving DNAzymes in the biosensing and therapeutic fields for signal generation and the modulation of gene expression, respectively. Additionally, multiple efforts have been made to study the effects of the reaction environment and the sequence of the catalytic core on the conversion of the substrate into product. However, most of these studies have only reported alterations of the general reaction course, but only a few have focused on how each individual reaction step is affected. In this work, we present for the first time a mathematical model that describes and predicts the reaction of the 10-23 RNA-cleaving DNAzyme. Furthermore, the model has been employed to study the effect of temperature, magnesium cations and shorter substrate-binding arms of the DNAzyme on the different kinetic rate constants, broadening the range of conditions in which the model can be exploited. In conclusion, this work depicts the prospects of such mathematical models to study and anticipate the course of a reaction given a particular environment.

Kinetics of radiation-induced DNA double-strand breaks through coarse-grained simulations

Double-strand breaks (DSBs), i.e., the covalent cut of the DNA backbone over both strands, are a detrimental outcome of cell irradiation, bearing chromosomal aberrations and leading to cell apoptosis. In the early stages of the evolution of a DSB, the disruption of the residual interactions between the DNA moieties drives the fracture of the helical layout; in spite of its biological significance, the details of this process are still largely uncertain. Here, we address the mechanical rupture of DNA by DSBs via coarse-grained molecular dynamics simulations: the setup involves a 3855-bp DNA filament and diverse DSB motifs, i.e., within a range of distances between strand breaks (or DSB distance). By employing a coarse-grained model of DNA, we access the molecular details and characteristic timescales of the rupturing process. A sequence-nonspecific, linear correlation is observed between the DSB distance and the internal energy contribution to the disruption of the residual (Watson-Crick and stacking) contacts between DNA moieties, which is seemingly driven by an abrupt, cooperative process. Moreover, we infer an exponential dependence of the characteristic rupture times on the DSB distances, which we associate to an Arrhenius-like law of thermally-activated processes. This work lays the foundations of a detailed, mechanistic assessment of DSBs in silico as a benchmark to both numerical simulations and data from single-molecule experiments.

Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.

Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

Low hybridization temperatures improve target capture success of invertebrate loci: a case study of leaf-footed bugs (Hemiptera: Coreoidea)

Target capture is widely used in phylogenomic, ecological and functional genomic studies. Bait sets that allow capture from a diversity of species can be advantageous, but high-sequence divergence from baits can limit yields. Currently, only four experimental comparisons of a critical target capture parameter, hybridization temperature, have been published. These have been in vertebrates, where bait divergences are typically low, and none include invertebrates where bait-target divergences may be higher. Most invertebrate capture studies use a fixed, high hybridization temperature to maximize the proportion of on-target data, but many report low locus recovery. Using leaf-footed bugs (Hemiptera: Coreoidea), we investigate the effect of hybridization temperature on capture success of ultraconserved elements targeted by (i) baits developed from divergent hemipteran genomes and (ii) baits developed from less divergent coreoid transcriptomes. Lower temperatures generally resulted in more contigs and improved recovery of targets despite a lower proportion of on-target reads, lower read depth and more putative paralogues. Hybridization temperatures had less of an effect when using transcriptome-derived baits, which is probably due to lower bait-target divergences and greater bait tiling density. Thus, accommodating low hybridization temperatures during target capture can provide a cost-effective, widely applicable solution to improve invertebrate locus recovery.

Purification of DNA oligonucleotides to improve hybridization chain reaction performance

Hybridization Chain Reaction (HCR) is a technique to generate a linear polymerization of oligonucleotide hairpins, used in multiple molecular biology methods. The HCR reaction is dependent on every hairpin being metastable in the absence of a triggering oligonucleotide and that every hairpin can continue the polymerization, which puts a strong demand on oligonucleotide quality. We show how further purification can greatly increase polymerization potential. It was found that a single extra PAGE-purification could greatly enhance hairpin polymerization both in solution and in situ. Purification using a ligation-based method further improved polymerization, yielding in situ immunoHCR stains at least 3.4-times stronger than a non-purified control. This demonstrates the importance of not only good sequence design of the oligonucleotide hairpins, but also the demand for high quality oligonucleotides to accomplish a potent and specific HCR.

Disruption of energetic and dynamic base pairing cooperativity in DNA duplexes by an abasic site

DNA duplex stability arises from cooperative interactions between multiple adjacent nucleotides that favor base pairing and stacking when formed as a continuous stretch rather than individually. Lesions and nucleobase modifications alter this stability in complex manners that remain challenging to understand despite their centrality to biology. Here, we investigate how an abasic site destabilizes small DNA duplexes and reshapes base pairing dynamics and hybridization pathways using temperature-jump infrared spectroscopy and coarse-grained molecular dynamics simulations. We show how an abasic site splits the cooperativity in a short duplex into two segments, which destabilizes small duplexes as a whole and enables metastable half-dissociated configurations. Dynamically, it introduces an additional barrier to hybridization by constraining the hybridization mechanism to a step-wise process of nucleating and zipping a stretch on one side of the abasic site and then the other.

Beneficial and detrimental effects of non-specific binding during DNA hybridization

DNA strands have to sample numerous states to find the alignment that maximizes Watson-Crick-Franklin base pairing. This process depends strongly on sequence, which affects the stability of the native duplex as well as the prevalence of non-native inter- and intramolecular helices. We present a theory that describes DNA hybridization as a three-stage process: diffusion, registry search, and zipping. We find that non-specific binding affects each of these stages in different ways. Mis-registered intermolecular binding in the registry search stage helps DNA strands sample different alignments and accelerates the hybridization rate. Non-native intramolecular structure affects all three stages by rendering portions of the molecule inert to intermolecular association, limiting mis-registered alignments to be sampled, and impeding the zipping process. Once in-register base pairs are formed, the stability of the native structure is important to hold the molecules together long enough for non-native contacts to break.

Mechanical DNA Origami to Investigate Biological Systems

The ability to self-assemble DNA nanodevices with programmed structural dynamics that can sense and respond to the local environment can enable transformative applications in fields including mechanobiology and nanomedicine. The responsive function of biomolecules is often driven by alterations in conformational distributions mediated by highly sensitive interactions with the local environment. In this review, the current state-of-the-art in constructing complex DNA geometries with dynamic and mechanical properties to enable a molecular scale force measurement is first summarized. Next, an overview of engineering modular DNA devices that interact with cell surfaces is highlighted detailing examples of mechanosensitive proteins and the force-induced dynamic molecular interaction on the downstream biochemical signaling. Finally, the challenges and an outlook on this promising class of DNA devices acting as nanomachines to operate at a low piconewton range suitable for a majority of biological effects or as hybrid materials to achieve higher tension exertion required for other biological investigations, are discussed.

Sequence-Dependent Melting and Refolding Dynamics of RNA UNCG Tetraloops Using Temperature-Jump/Drop Infrared Spectroscopy

Time-resolved temperature-jump/drop infrared (IR) spectroscopy has been used to measure the impact of stem base sequence on the melting and refolding dynamics of ribonucleic acid (RNA) tetraloops. A series of three 12-nucleotide RNA hairpin sequences were studied, each featuring a UACG tetraloop motif and a double-stranded stem containing four base pairs. In each case, the stem comprised three GC pairs plus a single AU base pair inserted at the closing point of the loop (RNA_loop), in the middle of the stem (RNA_mid), or at the stem terminus (RNA_end). Results from analogous DNA tetraloop (TACG) sequences were also obtained. Inclusion of AU or AT base pairs in the stem leads to faster melting of the stem-loop structure compared to a stem sequence featuring four GC base pairs while refolding times were found to be slower, consistent with a general reduction in stem-loop stability caused by the AU/AT pair. Independent measurement of the dynamic timescales for melting and refolding of ring vibrational modes of guanine (G_R) and adenine (A_R) provided position-specific insight into hairpin dynamics. The G_R-derived data showed that DNA sequences melted more quickly (0.5 ± 0.1 to 0.7 ± 0.1 μs at 70 °C) than analogous RNA sequences (4.3 ± 0.4 to 4.4 ± 0.3 μs at 70 °C). Position-sensitive data from the A_R modes suggests that DNA hairpins begin melting from the terminal end of the stem toward the loop while RNA sequences begin melting from the loop. Refolding timescales for both RNA and DNA hairpins were found to be similar (250 ± 50 μs at 70 °C) except for RNA_end and DNA_loop which refolded much more slowly (746 ± 36 and 430 ± 31 μs, respectively), showing that the refolding pathway is significantly impaired by the placement of AU/AT pairs at different points in the stem. We conclude that conformational changes of analogous pairs of RNA and DNA tetraloops proceed by different mechanisms.

Toehold-Mediated Strand Displacement in Random Sequence Pools

Toehold-mediated strand displacement (TMSD) has been used extensively for molecular sensing and computing in DNA-based molecular circuits. As these circuits grow in complexity, sequence similarity between components can lead to cross-talk, causing leak, altered kinetics, or even circuit failure. For small non-biological circuits, such unwanted interactions can be designed against. In environments containing a huge number of sequences, taking all possible interactions into account becomes infeasible. Therefore, a general understanding of the impact of sequence backgrounds on TMSD reactions is of great interest. Here, we investigate the impact of random DNA sequences on TMSD circuits. We begin by studying individual interfering strands and use the obtained data to build machine learning models that estimate kinetics. We then investigate the influence of pools of random strands and find that the kinetics are determined by only a small subpopulation of strongly interacting strands. Consequently, their behavior can be mimicked by a small collection of such strands. The equilibration of the circuit with the background sequences strongly influences this behavior, leading to up to 1 order of magnitude difference in reaction speed. Finally, we compare two established and one novel technique that speed up TMSD reactions in random sequence pools: a three-letter alphabet, protection of toeholds by intramolecular secondary structure, or by an additional blocking strand. While all of these techniques were useful, only the latter can be used without sequence constraints. We expect that our insights will be useful for the construction of TMSD circuits that are robust to molecular noise.

Measuring RNA UNCG Tetraloop Refolding Dynamics Using Temperature-Jump/Drop Infrared Spectroscopy

Determining the structural dynamics of RNA and DNA is essential to understanding their cellular function, but direct measurement of strand association or folding remains experimentally challenging. Here we illustrate a temperature-jump/drop method able to reveal refolding dynamics. Time-resolved temperature-jump/drop infrared spectroscopy is used to measure the melting and refolding dynamics of a 12-nucleotide RNA sequence comprising a UACG tetraloop and a four-base-pair double-stranded GC stem, comparing them to an equivalent DNA (TACG) sequence. Stem-loop melting occurred an order of magnitude more slowly in RNA than DNA (6.0 ± 0.1 μs versus 0.8 ± 0.1 μs at 70 °C). In contrast, the refolding dynamics of both sequences occurred on similar time scales (200 μs). While the melting and refolding dynamics of RNA and DNA hairpins both followed Arrhenius temperature dependences, refolding was characterized by an apparent negative activation energy, consistent with a mechanism involving multiple misfolded intermediates prior to zipping of the stem base pairs.

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.

The stability and number of nucleating interactions determine DNA hybridization rates in the absence of secondary structure

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson-Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.

Lattice model on the rate of DNA hybridization

We develop a lattice model on the rate of hybridization of the complementary single-stranded DNAs (c-ssDNAs). Upon translational diffusion mediated collisions, c-ssDNAs interpenetrate each other to form correct (cc), incorrect (icc), and trap correct contacts (tcc) inside the reaction volume. Correct contacts are those with exact registry matches, which leads to nucleation and zipping. Incorrect contacts are the mismatch contacts which are less stable compared to tcc, which can occur in the repetitive c-ssDNAs. Although tcc possess registry match within the repeating sequences, they are incorrect contacts in the view of the whole c-ssDNAs. The nucleation rate (k_{N}) is directly proportional to the collision rate and the average number of correct contacts (〈n_{cc}〉) formed when both c-ssDNAs interpenetrate each other. Detailed lattice model simulations suggest that 〈n_{cc}〉∝L/V where L is the length of c-ssDNAs and V is the reaction volume. Further numerical analysis revealed the scaling for the average radius of gyration of c-ssDNAs (R_{g}) with their length as R_{g}∝sqrt[L]. Since the reaction space will be approximately a sphere with radius equals to 2R_{g} and V∝L^{3/2}, one obtains k_{N}∝1/sqrt[L]. When c-ssDNAs are nonrepetitive, the overall renaturation rate becomes as k_{R}∝k_{N}L, and one finally obtains k_{R}∝sqrt[L] in line with the experimental observations. When c-ssDNAs are repetitive with a complexity of c, earlier models suggested the scaling k_{R}∝sqrt[L]/c, which breaks down at c=L. This clearly suggests the existence of at least two different pathways of renaturation in the case of repetitive c-ssDNAs, viz., via incorrect contacts and trap correct contacts. The trap correct contacts can lead to the formation of partial duplexes which can keep the complementary strands in the close proximity for a prolonged timescale. This is essential for the extended 1D slithering, inchworm movements, and internal displacement mechanisms which can accelerate the searching for the correct contacts. Clearly, the extent of slithering dynamics will be inversely proportional to the complexity. When the complexity is close to the length of c-ssDNAs, the pathway via incorrect contacts will dominate. When the complexity is much less than the length of c-ssDNA, pathway via trap correct contacts would be the dominating one.

Nucleic acid folding simulations using a physics-based atomistic free energy model

Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces-base stacking, counterion-induced backbone self-interactions, and base pairing-via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach.

DNA binding by polycomb-group proteins: searching for the link to CpG islands

Polycomb group proteins predominantly exist in polycomb repressive complexes (PRCs) that cooperate to maintain the repressed state of thousands of cell-type-specific genes. Targeting PRCs to the correct sites in chromatin is essential for their function. However, the mechanisms by which PRCs are recruited to their target genes in mammals are multifactorial and complex. Here we review DNA binding by polycomb group proteins. There is strong evidence that the DNA-binding subunits of PRCs and their DNA-binding activities are required for chromatin binding and CpG targeting in cells. In vitro, CpG-specific binding was observed for truncated proteins externally to the context of their PRCs. Yet, the mere DNA sequence cannot fully explain the subset of CpG islands that are targeted by PRCs in any given cell type. At this time we find very little structural and biophysical evidence to support a model where sequence-specific DNA-binding activity is required or sufficient for the targeting of CpG-dinucleotide sequences by polycomb group proteins while they are within the context of their respective PRCs, either PRC1 or PRC2. We discuss the current knowledge and open questions on how the DNA-binding activities of polycomb group proteins facilitate the targeting of PRCs to chromatin.

Electrochemical genosensor based on gold nanostars for the detection of Escherichia coli O157:H7 DNA

Escherichia coli O157:H7 (E. coli O157:H7) is an enterohemorrhagic E. coli (EHEC), which has been issued as a major threat to public health worldwide due to fatal contamination of water and food. Thus, its rapid and accurate detection has tremendous importance in environmental monitoring and human health. In this regard, we report a simple and sensitive electrochemical DNA biosensor by targeting Z3276 as a genetic marker in river water. The surface of the designed gold electrode was functionalized with gold nanostars and an aminated specific sensing probe of E. coli O157:H7 to fabricate the genosensor. Cyclic voltammetry (CV) and square wave voltammetry (SWV) techniques were applied for electrochemical characterization and detection. The synthesized gold nanostars were characterized using different characterization techniques. The fabricated DNA-based sensor exhibited a high selective ability for one, two, and three-base mismatched sequences. Regeneration, stability, selectivity, and kinetics of the bioassay were investigated. Under optimal conditions, the fabricated genosensor exhibited a linear response range of 10^-5 to 10^-17 μM in the standard sample and 7.3 to 1 × 10^-17 μM in water samples with a low limit of quantification of 0.01 zM in water samples. The detection strategy based on silver plated gold nanostars and DNA hybridization improved the sensitivity and specificity of the assay for E. coli O157:H7 detection in real water samples without filtration. The detection assay has the advantages of high selectivity, sensitivity, low amounts of reagents, short analysis time, commercialization, and potential application for the determination of other pathogenic bacteria.

Melting of confined DNA: static and dynamic properties

We study dsDNA (double strand DNA) melting in detail within varying strip-like confinement in a two-dimensional lattice model. The interplay between reduced configurational entropy and attractive base-pairing energy results in a non-monotonic melting profile of DNA. Structural transitions associated with confined DNA melting reveal a stretched or extended state for very strong confinement. By using the exact enumeration method, we investigate the emergence of a local denatured zone e.g. bubbles during DNA melting. The survival time of a single bubble within varying strip width is studied from the Fokker-Planck formalism by considering the bubble size as a reaction co-ordinate. We show that a simple lattice model can capture the sequence heterogeneity effect on DNA melting and bubble dynamics within the strip. Different time scales of bubble zipping for different DNA sequences are found, which may have potential applications in denaturation mapping.

Star-Polymer-DNA Gels Showing Highly Predictable and Tunable Mechanical Responses

Dynamically crosslinked gels are appealing materials for applications that require time-dependent mechanical responses. DNA duplexes are ideal crosslinkers for building such gels because of their excellent sequence addressability and flexible tunability in bond energy. However, the mechanical responses of most DNA gels are complicated and unpredictable. Here, a DNA gel with a highly homogeneous gel network and well predictable mechanical behaviors is demonstrated by using a pair of star-polymer-DNA precursors with presimulated DNA sequences showing the two-state transition. The melting curve analysis of the DNA gels reveals the good correspondence between the thermodynamic potentials of the DNA crosslinkers and the presimulated values by DNA calculators. Stress-relaxation tests and dissociation kinetics measurements show that the macroscopic relaxation time of the DNA gels is approximately equal to the lifetime of the DNA crosslinkers over 4 orders of magnitude from 0.1-2000 s. Furthermore, a series of durability tests find the DNA gels are hysteresis-less and self-healable after the applications of repeated temperature and mechanical stimuli. These results demonstrate the great potential of star-polymer-DNA precursors for building gels with predictable and tunable viscoelastic properties, suitable for applications such as stress-response extracellular matrices, injectable solids, and soft robotics.

OxDNA to Study Species Interactions

Molecular ecology uses molecular genetic data to answer traditional ecological questions in biogeography and biodiversity, among others. Several ecological principles, such as the niche hypothesis and the competitive exclusions, are based on the fact that species compete for resources. More in generally, it is now recognized that species interactions play a crucial role in determining the coexistence and abundance of species. However, experimentally controllable platforms, which allow us to study and measure competitions among species, are rare and difficult to implement. In this work, we suggest exploiting a Molecular Dynamics coarse-grained model to study interactions among single strands of DNA, representing individuals of different species, which compete for binding to other oligomers considered as resources. In particular, the well-established knowledge of DNA–DNA interactions at the nanoscale allows us to test the hypothesis that the maximum consecutive overlap between pairs of oligomers measure the species’ competitive advantages. However, we suggest that a more complex structure also plays a role in the ability of the species to successfully bind to the target resource oligomer. We complement the simulations with experiments on populations of DNA strands which qualitatively confirm our hypotheses. These tools constitute a promising starting point for further developments concerning the study of controlled, DNA-based, artificial ecosystems.

DNA-Based Biosensors for the Biochemical Analysis: A Review

In recent years, DNA-based biosensors have shown great potential as the candidate of the next generation biomedical detection device due to their robust chemical properties and customizable biosensing functions. Compared with the conventional biosensors, the DNA-based biosensors have advantages such as wider detection targets, more durable lifetime, and lower production cost. Additionally, the ingenious DNA structures can control the signal conduction near the biosensor surface, which could significantly improve the performance of biosensors. In order to show a big picture of the DNA biosensor's advantages, this article reviews the background knowledge and recent advances of DNA-based biosensors, including the functional DNA strands-based biosensors, DNA hybridization-based biosensors, and DNA templated biosensors. Then, the challenges and future directions of DNA-based biosensors are discussed and proposed.

The interplay of supercoiling and thymine dimers in DNA

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

Electrochemical Detection Platform Based on RGO Functionalized with Diazonium Salt for DNA Hybridization

In this paper, we propose an improved electrochemical platform based on graphene for the detection of DNA hybridization. Commercial screen-printed carbon electrodes (SPCEs) were used for this purpose due to their ease of functionalization and miniaturization opportunities. SPCEs were modified with reduced graphene oxide (RGO), offering a suitable surface for further functionalization. Therefore, aryl-carboxyl groups were integrated onto RGO-modified electrodes by electrochemical reduction of the corresponding diazonium salt to provide enough reaction sites for the covalent immobilization of amino-modified DNA probes. Our final goal was to determine the optimum conditions needed to fabricate a simple, label-free RGO-based electrochemical platform to detect the hybridization between two complementary single-stranded DNA molecules. Each modification step in the fabrication process was monitored by cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) using [Fe(CN)6]3-/4- as a redox reporter. Although, the diazonium electrografted layer displayed the expected blocking effect of the charge transfer, the next steps in the modification procedure resulted in enhanced electron transfer properties of the electrode interface. We suggest that the improvement in the charge transfer after the DNA hybridization process could be exploited as a prospective sensing feature. The morphological and structural characterization of the modified electrodes performed by scanning electron microscopy (SEM) and Raman spectroscopy, respectively, were used to validate different modification steps in the platform fabrication process.

Can DyeCycling break the photobleaching limit in single-molecule FRET?

Biomolecular systems, such as proteins, crucially rely on dynamic processes at the nanoscale. Detecting biomolecular nanodynamics is therefore key to obtaining a mechanistic understanding of the energies and molecular driving forces that control biomolecular systems. Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful technique to observe in real-time how a single biomolecule proceeds through its functional cycle involving a sequence of distinct structural states. Currently, this technique is fundamentally limited by irreversible photobleaching, causing the untimely end of the experiment and thus, a narrow temporal bandwidth of ≤ 3 orders of magnitude. Here, we introduce "DyeCycling", a measurement scheme with which we aim to break the photobleaching limit in smFRET. We introduce the concept of spontaneous dye replacement by simulations, and as an experimental proof-of-concept, we demonstrate the intermittent observation of a single biomolecule for one hour with a time resolution of milliseconds. Theoretically, DyeCycling can provide > 100-fold more information per single molecule than conventional smFRET. We discuss the experimental implementation of DyeCycling, its current and fundamental limitations, and specific biological use cases. Given its general simplicity and versatility, DyeCycling has the potential to revolutionize the field of time-resolved smFRET, where it may serve to unravel a wealth of biomolecular dynamics by bridging from milliseconds to the hour range.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.1007/s12274-022-4420-5 and is accessible for authorized users.

Kinetics of RNA and RNA:DNA Hybrid Strand Displacement

In nucleic acid nanotechnology, strand displacement is a widely used mechanism where one strand from a hybridized duplex is exchanged with an invading strand that binds to a toehold, a single-stranded region on the duplex. It is used to perform logic operations on a molecular level, initiate cascaded reactions, or even for in vivo diagnostics and treatments. While systematic experimental studies have been carried out to probe the kinetics of strand displacement in DNA with different toehold lengths, sequences, and mismatch positions, there has not been a comparable investigation of RNA or RNA-DNA hybrid systems. Here, we experimentally study how toehold length, toehold location (5' or 3' end of the strand), and mismatches influence the strand displacement kinetics. We observe reaction acceleration with increasing toehold length and placement of the toehold at the 5' end of the substrate. We find that mismatches closer to the interface of toehold and duplex slow down the reaction more than remote mismatches. A comparison of RNA and DNA displacement with hybrid displacement (RNA invading DNA or DNA invading RNA) is partly explainable by the thermodynamic stabilities of the respective toehold regions, but also suggests that the rearrangement from B-form to A-form helix in the case of RNA invading DNA might play a role in the kinetics.

Determining Sequence-Dependent DNA Oligonucleotide Hybridization and Dehybridization Mechanisms Using Coarse-Grained Molecular Simulation, Markov State Models, and Infrared Spectroscopy

A robust understanding of the sequence-dependent thermodynamics of DNA hybridization has enabled rapid advances in DNA nanotechnology. A fundamental understanding of the sequence-dependent kinetics and mechanisms of hybridization and dehybridization remains comparatively underdeveloped. In this work, we establish new understanding of the sequence-dependent hybridization/dehybridization kinetics and mechanism within a family of self-complementary pairs of 10-mer DNA oligomers by integrating coarse-grained molecular simulation, machine learning of the slow dynamical modes, data-driven inference of long-time kinetic models, and experimental temperature-jump infrared spectroscopy. For a repetitive ATATATATAT sequence, we resolve a rugged dynamical landscape comprising multiple metastable states, numerous competing hybridization/dehybridization pathways, and a spectrum of dynamical relaxations. Introduction of a G:C pair at the terminus (GATATATATC) or center (ATATGCATAT) of the sequence reduces the ruggedness of the dynamics landscape by eliminating a number of metastable states and reducing the number of competing dynamical pathways. Only by introducing a G:C pair midway between the terminus and the center to maximally disrupt the repetitive nature of the sequence (ATGATATCAT) do we recover a canonical "all-or-nothing" two-state model of hybridization/dehybridization with no intermediate metastable states. Our results establish new understanding of the dynamical richness of sequence-dependent kinetics and mechanisms of DNA hybridization/dehybridization by furnishing quantitative and predictive kinetic models of the dynamical transition network between metastable states, present a molecular basis with which to understand experimental temperature jump data, and furnish foundational design rules by which to rationally engineer the kinetics and pathways of DNA association and dissociation for DNA nanotechnology applications.

Measuring Excess Heat Capacities of Deoxyribonucleic Acid (DNA) Folding at the Single-Molecule Level

Measurements of the thermodynamic properties of biomolecular folding (ΔG°, ΔH°, ΔS°, etc.) provide a wealth of information on the folding process and have long played a central role in biophysical investigation. In particular, the excess heat capacity of folding (ΔC_P) is crucial, as typically measured in bulk ensemble studies by differential scanning calorimetry (DSC) and isothermal titration calorimetry (ITC). Here, we report the first measurements of ΔC_P at the single-molecule level using the single-molecule fluorescence resonance energy transfer (smFRET) as well as the very first measurements of the heat capacity change associated with achieving the transition state (ΔC^‡_P) for nucleic acid folding. The deoxyribonucleic acid (DNA) hairpin used in these studies exhibits an excess heat capacity for hybridization (ΔC_P = -340 ± 60 J/mol/K per base pair) consistent with the range of literature expectations (ΔC_P = -100 to -420 J/mol/K per base pair). Furthermore, the measured activation heat capacities (ΔC^‡_P) for such hairpin unfolding are consistent with a folding transition state containing few fully formed base pairs, in agreement with prevailing models of DNA hybridization.

Factors and methods to modulate DNA hybridization kinetics

DNA oligonucleotides are widely used in a diverse range of research fields from analytical chemistry, molecular biology, nanotechnology to drug delivery. In these applications, DNA hybridization is often the most important enabling reaction. Achieving control over hybridization kinetics and a high yield of hybridized products is needed to ensure high-quality and reproducible results. Since DNA strands are highly negatively charged and can also fold upon itself to form various intramolecular structures, DNA hybridization needs to overcome these barriers. Nucleation and diffusion are two main kinetic limiting steps although their relative importance differs in different conditions. The effects of length and sequence, temperature, pH, salt concentration, cationic polymers, organic solvents, freezing and crowding agents are summarized in the context of overcoming these barriers. This article will help researchers in the biotechnology-related fields to better understand and control DNA hybridization, as well as provide a landscape for future work in simulation and experiment to optimize DNA hybridization systems.

Folding Dynamics of DNA G-Quadruplexes Probed by Millisecond Temperature Jump Circular Dichroism

G-quadruplexes play important roles in cellular regulatory functions, but despite significant experimental and theoretical efforts, their folding mechanisms remain poorly understood. In this context, we developed a T-jump experiment to access the thermal denaturation and renaturation dynamics of short intramolecular G-quadruplexes in vitro, on the time scale of a few hundred milliseconds. With this new setup, we compared the thermal denaturation and renaturation kinetics of three antiparallel topologies made of the human telomeric sequences d[(5'-GGG(TTAGGG)₃-3']/Na⁺ and d[5'-AGGG(TTAGGG)₃-3']/Na⁺ and the thrombin-binding aptamer sequence d[5'-GGTTGGTGTGGTTGG-3']/K⁺, with those of the parallel topology made of the human CEB25 minisatellite d[5'-AAGGGTGGGTGTAAGTGTGGGTGGGT-3']/Na⁺. In all cases, exponential kinetics of the order of several hundred milliseconds were observed. Measurements performed for different initial temperatures revealed distinct denaturation and renaturation dynamics, ruling out a simple two-state mechanism. The parallel topology, in which all guanines adopt an anti conformation, displays much slower dynamics than antiparallel topologies associated with very low activation barriers. This behavior can be explained by the constrained conformational space due to the presence of the single-base propeller loops that likely hinders the movement of the coiled DNA strand and reduces the contribution of the entropy during the renaturation process at high temperatures.

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.

OxDNA.org: a public webserver for coarse-grained simulations of DNA and RNA nanostructures

OxDNA and oxRNA are popular coarse-grained models used by the DNA/RNA nanotechnology community to prototype, analyze and rationalize designed DNA and RNA nanostructures. Here, we present oxDNA.org, a graphical web interface for running, visualizing and analyzing oxDNA and oxRNA molecular dynamics simulations on a GPU-enabled high performance computing server. OxDNA.org automatically generates simulation files, including a multi-step relaxation protocol for structures exported in non-physical states from DNA/RNA design tools. Once the simulation is complete, oxDNA.org provides an interactive visualization and analysis interface using the browser-based visualizer oxView to facilitate the understanding of simulation results for a user's specific structure. This online tool significantly lowers the entry barrier of integrating simulations in the nanostructure design pipeline for users who are not experts in the technical aspects of molecular simulation. The webserver is freely available at oxdna.org.

Complex dependence of CRISPR-Cas9 binding strength on guide RNA spacer lengths

It is established that for CRISPR-Cas9 applications guide RNAs with 17-20 bp long spacer sequences are optimal for accurate target binding and cleavage. In this work we perform cell-free CRISPRa (CRISPR activation) and CRISPRi (CRISPR inhibition) experiments to demonstrate the existence of a complex dependence of CRISPR-Cas9 binding as a function of the spacer length and complementarity. Our results show that significantly truncated or mismatched spacer sequences can form stronger guide-target bonds than the conventional 17-20 bp long spacers. To explain this phenomenon, we take into consideration previous structural and single-molecule CRISPR-Cas9 experiments and develop a novel thermodynamic model of CRISPR-Cas9 target recognition.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Kinetic Microscale Thermophoresis for Simultaneous Measurement of Binding Affinity and Kinetics

Microscale thermophoresis (MST) is a versatile technique to measure binding affinities of binder-ligand systems, based on the directional movement of molecules in a temperature gradient. We extended MST to measure binding kinetics as well as binding affinity in a single experiment by increasing the thermal dissipation of the sample. The kinetic relaxation fingerprints were derived from the fluorescence changes during thermodynamic re-equilibration of the sample after local heating. Using this method, we measured DNA hybridization on-rates and off-rates in the range 10⁴ -10⁶  m^-1  s^-1 and 10^-4 -10^-1  s^-1 , respectively. We observed the expected exponential dependence of the DNA hybridization off-rates on salt concentration, strand length and inverse temperature. The measured on-rates showed a linear dependence on salt concentration and weak dependence on strand length and temperature. For biomolecular interactions with large enthalpic contributions, the kinetic MST technique offers a robust, cost-effective and immobilization-free determination of kinetic rates and binding affinity simultaneously, even in crowded solutions.

Nanomechanics of self-assembled DNA building blocks

DNA has become a powerful platform to design functional nanodevices. DNA nanodevices are often composed of self-assembled DNA building blocks that differ significantly from the structure of native DNA. In this study, we present Flow Force Microscopy as a massively parallel approach to study the nanomechanics of DNA self-assemblies on the single-molecular level. The high-throughput experiments performed in a simple microfluidic channel enable statistically meaningful studies with nanometer scale precision in a time frame of several minutes. A surprisingly high flexibility was observed for a typical construct used in DNA origami, reflected in a persistence length of 10.2 nm, a factor of five smaller than for native DNA. The enhanced flexibility is attributed to the discontinuous backbone of DNA self-assemblies that facilitate base pair opening by thermal fluctuations at the end of hybridized oligomers. We believe that the results will contribute to the fundamental understanding of DNA nanomechanics and help to improve the design of DNA nanodevices with applications in biological analysis and clinical research.

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.

DNA hybridisation kinetics using single-molecule fluorescence imaging

Deoxyribonucleic acid (DNA) hybridisation plays a key role in many biological processes and nucleic acid biotechnologies, yet surprisingly there are many aspects about the process which are still unknown. Prior to the invention of single-molecule microscopy, DNA hybridisation experiments were conducted at the ensemble level, and thus it was impossible to directly observe individual hybridisation events and understand fully the kinetics of DNA hybridisation. In this mini-review, recent single-molecule fluorescence-based studies of DNA hybridisation are discussed, particularly for short nucleic acids, to gain more insight into the kinetics of DNA hybridisation. As well as looking at single-molecule studies of intrinsic and extrinsic factors affecting DNA hybridisation kinetics, the influence of the methods used to detect hybridisation of single DNAs is considered. Understanding the kinetics of DNA hybridisation not only gives insight into an important biological process but also allows for further advancements in the growing field of nucleic acid biotechnology.

Programming folding cooperativity of the dimeric i-motif with DNA frameworks for sensing small pH variations

The response sensitivity of a molecular sensor is determined by the folding cooperativity of its responsive module. Using an H⁺-responsive dimeric DNA i-motif as a model, we demonstrate the enhancement of its folding cooperativity through preorganization by a DNA framework, and with it we fabricate robust intracellular pH sensors with high response sensitivity.