Influence of thermodynamically unfavorable secondary structures on DNA hybridization kinetics

Rapid in situ RNA imaging based on Cas12a thrusting strand displacement reaction

RNA In situ imaging through DNA self-assembly is advantaged in illustrating its structures and functions with high-resolution, while the limited reaction efficiency and time-consuming operation hinder its clinical application. Here, we first proposed a new strand displacement reaction (SDR) model (Cas12a thrusting SDR, CtSDR), in which Cas12a could overcome the inherent reaction limitation and dramatically enhance efficiency through energy replenishment and by-product consumption. The target-initiated CtSDR amplification was established for RNA analysis, with order of magnitude lower limit of detection (LOD) than the Cas13a system. The CtSDR-based RNA in situ imaging strategy was developed to monitor intra-cellular microRNA expression change and delineate the landscape of oncogenic RNA in 66 clinic tissue samples, possessing a clear advantage over classic in situ hybridization (ISH) in terms of operation time (1 h versus 14 h) while showing comparable sensitivity and specificity. This work presents a promising approach to developing advanced molecular diagnostic tools.

Generation of DNA oligomers with similar chemical kinetics via in-silico optimization

Networks of interacting DNA oligomers are useful for applications such as biomarker detection, targeted drug delivery, information storage, and photonic information processing. However, differences in the chemical kinetics of hybridization reactions, referred to as kinetic dispersion, can be problematic for some applications. Here, it is found that limiting unnecessary stretches of Watson-Crick base pairing, referred to as unnecessary duplexes, can yield exceptionally low kinetic dispersions. Hybridization kinetics can be affected by unnecessary intra-oligomer duplexes containing only 2 base-pairs, and such duplexes explain up to 94% of previously reported kinetic dispersion. As a general design rule, it is recommended that unnecessary intra-oligomer duplexes larger than 2 base-pairs and unnecessary inter-oligomer duplexes larger than 7 base-pairs be avoided. Unnecessary duplexes typically scale exponentially with network size, and nearly all networks contain unnecessary duplexes substantial enough to affect hybridization kinetics. A new method for generating networks which utilizes in-silico optimization to mitigate unnecessary duplexes is proposed and demonstrated to reduce in-vitro kinetic dispersions as much as 96%. The limitations of the new design rule and generation method are evaluated in-silico by creating new oligomers for several designs, including three previously programmed reactions and one previously engineered structure.

Enzyme-Free Nucleic Acid Circuits for Fold-Change Detection

Fold-change detection is widespread in sensory systems of various organisms. Dynamic DNA nanotechnology provides an important toolbox for reproducing structures and responses of cellular circuits. In this work, we construct an enzyme-free nucleic acid circuit based on the incoherent feed-forward loop using toehold-mediated DNA strand displacement reactions and explore its dynamic behaviors. The mathematical model based on ordinary differential equations is used to evaluate the parameter regime required for fold-change detection. After selecting appropriate parameters, the constructed synthetic circuit exhibits approximate fold-change detection for multiple rounds of inputs with different initial concentrations. This work is anticipated to shed new light on the design of DNA dynamic circuits in the enzyme-free environment.

FRET Imaging of Nonuniformly Distributed DNA SAMs on Gold Reveals the Role Played by the Donor/Acceptor Ratio and the Local Environment in Measuring the Rate of Hybridization

Mixed DNA SAMs labeled with a fluorophore (either AlexaFluor488 or AlexaFluor647) were prepared on a single crystal gold bead electrode using potential-assisted thiol exchange and studied using Förster resonance energy transfer (FRET). A measure of the local environment of the DNA SAM (e.g., crowding) was possible using FRET imaging on these surfaces since electrodes prepared this way have a range of surface densities (Γ_DNA). The FRET signal was strongly dependent on Γ_DNA and on the ratio of AlexaFluor488 to AlexaFluor647 used to make the DNA SAM, which were consistent with a model of FRET in 2D systems. FRET was shown to provide a direct measure of the local DNA SAM arrangement on each crystallographic region of interest providing a direct assessment of the probe environment and its influence on the rate of hybridization. The kinetics of duplex formation for these DNA SAMs was also studied using FRET imaging over a range of coverages and DNA SAM compositions. Hybridization of the surface-bound DNA increased the average distance between the fluorophore label and the gold electrode surface and decreased the distance between the donor (D) and acceptor (A), both of which result in an increase in FRET intensity. This increase in FRET was modeled using a second order Langmuir adsorption rate equation, reflecting the fact that both D and A labeled DNA are required to become hybridized to observe a FRET signal. The self-consistent analysis of the hybridization rates on low and high coverage regions on the same electrode showed that the low coverage regions achieved full hybridization 5× faster than the higher coverage regions, approaching rates typically found in solution. The relative increase in FRET intensity from each region of interest was controlled by manipulating the donor to acceptor composition of the DNA SAM without changing the rate of hybridization. The FRET response can be optimized by controlling the coverage and the composition of the DNA SAM sensor surface and could be further improved with the use of a FRET pair with a larger (e.g., > 5 nm) Förster radius.

Predicting DNA kinetics with a truncated continuous-time Markov chain method

Predicting the kinetics of reactions involving nucleic acid strands is a fundamental task in biology and biotechnology. Reaction kinetics can be modeled as an elementary step continuous-time Markov chain, where states correspond to secondary structures and transitions correspond to base pair formation and breakage. Since the number of states in the Markov chain could be large, rates are determined by estimating the mean first passage time from sampled trajectories. As a result, the cost of kinetic predictions becomes prohibitively expensive for rare events with extremely long trajectories. Also problematic are scenarios where multiple predictions are needed for the same reaction, e.g., under different environmental conditions, or when calibrating model parameters, because a new set of trajectories is needed multiple times. We propose a new method, called pathway elaboration, to handle these scenarios. Pathway elaboration builds a truncated continuous-time Markov chain through both biased and unbiased sampling. The resulting Markov chain has moderate state space size, so matrix methods can efficiently compute reaction rates, even for rare events. Also the transition rates of the truncated Markov chain can easily be adapted when model or environmental parameters are perturbed, making model calibration feasible. We illustrate the utility of pathway elaboration on toehold-mediated strand displacement reactions, show that it well-approximates trajectory-based predictions of unbiased elementary step models on a wide range of reaction types for which such predictions are feasible, and demonstrate that it performs better than alternative truncation-based approaches that are applicable for mean first passage time estimation. Finally, in a small study, we use pathway elaboration to optimize the Metropolis kinetic model of Multistrand, an elementary step simulator, showing that the optimized parameters greatly improve reaction rate predictions. Our framework and dataset are available at https://github.com/DNA-and-Natural-Algorithms-Group/PathwayElaboration.

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.

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.

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.

Nearest-neighbour transition-state analysis for nucleic acid kinetics

We used stopped-flow to monitor hypochromicity for 43 oligonucleotide duplexes to study nucleic acid kinetics and extract transition-state parameters for association and dissociation. Reactions were performed in 1.0 M NaCl (for literature comparisons) and 2.2 mM MgCl2 (PCR conditions). Dissociation kinetics depended on sequence, increased exponentially with temperature, and transition-state parameters inversely correlated to thermodynamic parameters (r = -0.99). Association had no consistent enthalpic component, varied little with temperature or sequence, and poorly correlated to thermodynamic parameters (r = 0.28). Average association rates decreased 78% in MgCl2 compared to NaCl while dissociation was relatively insensitive to ionic conditions. A nearest-neighbour kinetic model for dissociation predicted rate constants within 3-fold of literature values (n = 11). However, a nearest-neighbour model for association appeared overparameterized and inadequate for predictions. Kinetic predictions were used to simulate published high-speed (<1 min) melting analysis and extreme (<2 min) PCR experiments. Melting simulations predicted apparent melting temperatures increase on average 2.4°C when temperature ramp rates increased from 0.1 to 32°C/s, compared to 2.8°C reported in the literature. PCR simulations revealed that denaturation kinetics are dependent on the thermocycling profile. Simulations overestimated annealing efficiencies at shorter annealing times and suggested that polymerase interactions contribute to primer-template complex stability at extension temperatures.

An alternative approach to nucleic acid memory

DNA is a compelling alternative to non-volatile information storage technologies due to its information density, stability, and energy efficiency. Previous studies have used artificially synthesized DNA to store data and automated next-generation sequencing to read it back. Here, we report digital Nucleic Acid Memory (dNAM) for applications that require a limited amount of data to have high information density, redundancy, and copy number. In dNAM, data is encoded by selecting combinations of single-stranded DNA with (1) or without (0) docking-site domains. When self-assembled with scaffold DNA, staple strands form DNA origami breadboards. Information encoded into the breadboards is read by monitoring the binding of fluorescent imager probes using DNA-PAINT super-resolution microscopy. To enhance data retention, a multi-layer error correction scheme that combines fountain and bi-level parity codes is used. As a prototype, fifteen origami encoded with 'Data is in our DNA!\n' are analyzed. Each origami encodes unique data-droplet, index, orientation, and error-correction information. The error-correction algorithms fully recover the message when individual docking sites, or entire origami, are missing. Unlike other approaches to DNA-based data storage, reading dNAM does not require sequencing. As such, it offers an additional path to explore the advantages and disadvantages of DNA as an emerging memory material.

How Surfaces Affect Hybridization Kinetics

Hybridization between nucleic acid strands immobilized on a solid support with partners in solution is widely practiced in bioanalytical technologies and materials science. An important fundamental aspect of understanding these reactions is the role played by immobilization in the dynamics of duplex formation and disassembly. This report reviews and analyzes literature kinetic data to identify commonly observed trends and to correlate them with probable molecular mechanisms. The analysis reveals that while under certain conditions impacts from immobilization are minimal so that surface and solution hybridization kinetics are comparable, it is more typical to observe pronounced offsets between the two scenarios. In the forward (hybridization) direction, rates at the surface commonly decrease by one to two decades relative to solution, while in the reverse direction rates of strand separation at the surface can exceed those in solution by tens of decades. By recasting the deviations in terms of activation barriers, a consensus of how immobilization impacts nucleation, zipping, and strand separation can be conceived within the classical mechanism in which duplex formation is rate limited by preassembly of a nucleus a few base pairs in length, while dehybridization requires the cumulative breakup of base pairs along the length of a duplex. Evidence is considered for how excess interactions encountered on solid supports impact these processes.

Molybdenum disulfide@5-carboxyfluorescein-probe biosensor for unamplified specific fragment detection in long nucleic acids based on magnetic composite probe-actuated deblocking of secondary structure

Secondary structures in long circulating tumor nucleic acids have potential obstacles for specific location point hybridized detection of gene fragments. The exploration of biosensing strategies requires selectively changing the nucleic acids conformation and inducing signal switching. Herein, a self-assembled magnetic composite probe (MCP) was fabricated by the hybridization reaction of Linker DNA and a "Y"-junction-DNA nanostructure on the surface of magnetic beads, contributing to the capture, secondary structure unlocking, and enrichment of the PML/RARα DNA "L" subtype. Then, by integrating the MCP-actuated reactor, a one-step "off-on" signal switching MoS2@FAM-probe biosensing method was developed for the efficient detection of the PML/RARα DNA "L" subtype. The proposed biosensor was capable of detecting 100 bases PML/RARα DNA "L" subtype with a wide linear range of 1 pM to 200 nM and a limit of detection down to 0.223 pM without signal amplification. In addition, the biosensing method was successfully applied for the detection of target in serum samples. It is worth pointing out that this simple biosensing strategy could enable long nucleic acids fragments with secondary structures from ctDNA and ctRNA to be quantitatively assayed based on direct hybridization.

Tethered multifluorophore motion reveals equilibrium transition kinetics of single DNA double helices

Understanding cellular functions and dysfunctions often begins with quantifying the interactions between the binding partners involved in the processes. Learning about the kinetics of the interactions is of particular importance to understand the dynamics of cellular processes. We created a tethered multifluorophore motion assay using DNA origami that enables over 1-hour-long recordings of the statistical binding and unbinding of single pairs of biomolecules directly in equilibrium. The experimental concept is simple and the data interpretation is very direct, which makes the system easy to use for a wide variety of researchers. Due to the modularity and addressability of the DNA origami-based assay, our system may be readily adapted to study various other molecular interactions.

We describe a tethered multifluorophore motion assay based on DNA origami for revealing bimolecular reaction kinetics on the single-molecule level. Molecular binding partners may be placed at user-defined positions and in user-defined stoichiometry; and binding states are read out by tracking the motion of quickly diffusing fluorescent reporter units. Multiple dyes per reporter unit enable singe-particle observation for more than 1 hour. We applied the system to study in equilibrium reversible hybridization and dissociation of complementary DNA single strands as a function of tether length, cation concentration, and sequence. We observed up to hundreds of hybridization and dissociation events per single reactant pair and could produce cumulative statistics with tens of thousands of binding and unbinding events. Because the binding partners per particle do not exchange, we could also detect subtle heterogeneity from molecule to molecule, which enabled separating data reflecting the actual target strand pair binding kinetics from falsifying influences stemming from chemically truncated oligonucleotides. Our data reflected that mainly DNA strand hybridization, but not strand dissociation, is affected by cation concentration, in agreement with previous results from different assays. We studied 8-bp-long DNA duplexes with virtually identical thermodynamic stability, but different sequences, and observed strongly differing hybridization kinetics. Complementary full-atom molecular-dynamics simulations indicated two opposing sequence-dependent phenomena: helical templating in purine-rich single strands and secondary structures. These two effects can increase or decrease, respectively, the fraction of strand collisions leading to successful nucleation events for duplex formation.