skip to main content


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

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.

 
more » « less
Award ID(s):
1807809
NSF-PAR ID:
10469746
Author(s) / Creator(s):
; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Communications Chemistry
Volume:
6
Issue:
1
ISSN:
2399-3669
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Peptide nucleic acid (PNA) is a unique synthetic nucleic acid analog that has been adopted for use in many biological applications. These applications rely upon the robust Franklin–Watson–Crick base pairing provided by PNA, particularly at lower ionic strengths. However, our understanding of the relationship between the kinetics of PNA:DNA hybridization and ionic strength is incomplete. Here we measured the kinetics of association and dissociation of PNA with DNA across a range of ionic strengths and temperatures at single-molecule resolution using total internal reflection fluorescence imaging. Unlike DNA:DNA duplexes, PNA:DNA duplexes are more stable at lower ionic strength, and we demonstrate that this is due to a higher association rate. While the dissociation rate of PNA:DNA duplexes is largely insensitive to ionic strength, it is significantly lower than that of DNA:DNA duplexes having the same number and sequence of base pairing interactions. The temperature dependence of PNA:DNA kinetic rate constants indicate a significant enthalpy barrier to duplex dissociation, and to a lesser extent, duplex formation. This investigation into the kinetics of PNA:DNA hybridization provides a framework towards better understanding and design of PNA sequences for future applications. 
    more » « less
  2. Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism. 
    more » « less
  3. It has been shown previously in protonated, deprotonated and ionized guanine–cytosine base pairs that intra-base pair proton transfer from the N1–H at the Watson–Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure–kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e. , [9MG·9MG]˙ + , [9MOG·9MG]˙ + and [9MOG·9MOG]˙ + . Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion–molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1–H at the Watson–Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer. 
    more » « less
  4. Abstract

    Oligonucleotide hybridization is crucial in various biological, prebiotic and nanotechnological processes, including gene regulation, non-enzymatic primer extension and DNA nanodevice assembly. Although extensive research has focused on the thermodynamics and kinetics of nucleic acid hybridization, the behavior of complex mixtures and the outcome of competition for target binding remain less well understood. In this study, we investigate the impact of mismatches and bulges in a 12 bp DNA or RNA duplex on its association (kon) and dissociation (koff) kinetics. We find that such defects have relatively small effects on the association kinetics, while the dissociation kinetics vary in a position-dependent manner by up to 6 orders of magnitude. Building upon this observation, we explored a competition scenario involving multiple oligonucleotides, and observed a transient low specificity of probe hybridization to fully versus partially complementary targets in solution. We characterize these long-lived metastable states and their evolution toward equilibrium, and show that sufficiently long-lived mis-paired duplexes can serve as substrates for prebiotically relevant chemical copying reactions. Our results suggest that transient low accuracy states may spontaneously emerge within all complex nucleic acid systems comprising a large enough number of competing strands, with potential repercussions for gene regulation in the realm of modern biology and the prebiotic preservation of genetic information.

     
    more » « less
  5. Nucleic acid kinetic simulators aim to predict the kinetics of interacting nucleic acid strands. Many simulators model the kinetics of interacting nucleic acid strands as continuous-time Markov chains (CTMCs). States of the CTMCs represent a collection of secondary structures, and transitions between the states correspond to the forming or breaking of base pairs and are determined by a nucleic acid kinetic model. The number of states these CTMCs can form may be exponentially large in the length of the strands, making two important tasks challenging, namely, mean first passage time (MFPT) estimation and parameter estimation for kinetic models based on MFPTs. Gillespie’s stochastic simulation algorithm (SSA) is widely used to analyze nucleic acid folding kinetics, but could be computationally expensive for reactions whose CTMC has a large state space or for slow reactions. It could also be expensive for arbitrary parameter sets that occur in parameter estimation. Our work addresses these two challenging tasks, in the full state space of all non-pseudoknotted secondary structures of each reaction. In the first task, we show how to use a reduced variance stochastic simulation algorithm (RVSSA), which is adapted from SSA, to estimate the MFPT of a reaction’s CTMC. In the second task, we estimate model parameters based on MFPTs. To this end, first, we show how to use a generalized method of moments (GMM) approach, where we minimize a squared norm of moment functions that we formulate based on experimental and estimated MFPTs. Second, to speed up parameter estimation, we introduce a fixed path ensemble inference (FPEI) approach, that we adapt from RVSSA. We implement and evaluate RVSSA and FPEI using the Multistrand kinetic simulator. In our experiments on a dataset of DNA reactions, FPEI speeds up parameter estimation compared to inference using SSA, by more than a factor of three for slow reactions. Also, for reactions with large state spaces, it speeds up parameter estimation by more than a factor of two. 
    more » « less