skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Free Energy Landscape and Conformational Kinetics of Hoogsteen Base Pairing in DNA vs. RNA
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
Award ID(s):
2028443
PAR ID:
10227752
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Biophysical journal
Volume:
119
Issue:
8
ISSN:
0006-3495
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; (, Frontiers in Molecular Biosciences)
    Numerous biological processes and mechanisms depend on details of base pairing and hydrogen bonding in DNA. Hydrogen bonds are challenging to quantify by X-ray crystallography and cryo-EM due to difficulty of visualizing hydrogen atom locations but can be probed with site specificity by NMR spectroscopy in solution and the solid state with the latter particularly suited to large, slowly tumbling DNA complexes. Recently, we showed that low-temperature dynamic nuclear polarization (DNP) enhanced solid-state NMR is a valuable tool for distinguishing Hoogsteen base pairs (bps) from canonical Watson-Crick bps in various DNA systems under native-like conditions. Here, using a model 12-mer DNA duplex containing two central adenine-thymine (A-T) bps in either Watson-Crick or Hoogsteen confirmation, we demonstrate DNP solid-state NMR measurements of thymine N3-H3 bond lengths, which are sensitive to details of N-H···N hydrogen bonding and permit hydrogen bonds for the two bp conformers to be systematically compared within the same DNA sequence context. For this DNA duplex, effectively identical TN3-H3 bond lengths of 1.055 ± 0.011 Å and 1.060 ± 0.011 Å were found for Watson-Crick A-T and Hoogsteen A (syn)-T base pairs, respectively, relative to a reference amide bond length of 1.015 ± 0.010 Å determined for N-acetyl-valine under comparable experimental conditions. Considering that prior quantum chemical calculations which account for zero-point motions predict a somewhat longer effective peptide N-H bond length of 1.041 Å, in agreement with solution and solid-state NMR studies of peptides and proteins at ambient temperature, to facilitate direct comparisons with these earlier studies TN3-H3 bond lengths for the DNA samples can be readily scaled appropriately to yield 1.083 Å and 1.087 Å for Watson-Crick A-T and Hoogsteen A (syn)-T bps, respectively, relative to the 1.041 Å reference peptide N-H bond length. Remarkably, in the context of the model DNA duplex, these results indicate that there are no significant differences in N-H···N A-T hydrogen bonds between Watson-Crick and Hoogsteen bp conformers. More generally, high precision measurements of N-H bond lengths by low-temperature DNP solid-state NMR based methods are expected to facilitate detailed comparative analysis of hydrogen bonding for a range of DNA complexes and base pairing environments. 
    more » « less
  2. ; ; ; (, Physical Chemistry Chemical Physics)
    null (Ed.)
    We investigated the collision-induced dissociation (CID) reactions of a protonated Hoogsteen 9-methylguanine–1-methylcytosine base pair (HG-[9MG·1MC + H] + ), which aims to address the mystery of the literature reported “anomaly” in product ion distributions and compare the kinetics of a Hoogsteen base pair with its Watson-Crick isomer WC-[9MG·1MC + H] + (reported recently by Sun et al. ; Phys. Chem. Chem. Phys. , 2020, 22 , 24986). Product ion cross sections and branching ratios were measured as a function of center-of-mass collision energy using guided-ion beam tandem mass spectrometry, from which base-pair dissociation energies were determined. Product structures and energetics were assessed using various theories, of which the composite DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97XD/6-311++G(d,p) was adopted as the best-performing method for constructing a reaction potential energy surface. The statistical Rice–Ramsperger–Kassel–Marcus theory was found to provide a useful framework for rationalizing the dominating abundance of [1MC + H] + over [9MG + H] + in the fragment ions of HG-[9MG·1MC + H] + . The kinetics analysis proved the necessity for incorporating into kinetics modeling not only the static properties of reaction minima and transition states but more importantly, the kinetics of individual base-pair conformers that have formed in collisional activation. The analysis also pinpointed the origin of the statistical kinetics of HG-[9MG·1MC + H] + vs. the non-statistical behavior of WC-[9MG·1MC + H] + in terms of their distinctively different intra-base-pair hydrogen-bonds and consequently the absence of proton transfer between the N1 position of 9MG and the N3′ of 1MC in the Hoogsteen base pair. Finally, the Hoogsteen base pair was examined in the presence of a water ligand, i.e. , HG-[9MG·1MC + H] + ·H 2 O. Besides the same type of base-pair dissociation as detected in dry HG-[9MG·1MC + H] + , secondary methanol elimination was observed via the S N 2 reaction of water with nucleobase methyl groups. 
    more » « less
  3. ; ; ; ; ; ; ; (, Proceedings of the National Academy of Sciences)
    The majority of base pairs in double-stranded DNA exist in the canonical Watson-Crick geometry. However, they can also adopt alternate Hoogsteen conformations in various complexes of DNA with proteins and small molecules, which are key for biological function and mechanism. While detection of Hoogsteen base pairs in large DNA complexes and assemblies poses considerable challenges for traditional structural biology techniques, we show here that multidimensional dynamic nuclear polarization–enhanced solid-state NMR can serve as a unique spectroscopic tool for observing and distinguishing Watson-Crick and Hoogsteen base pairs in a broad range of DNA systems based on characteristic NMR chemical shifts and internuclear dipolar couplings. We illustrate this approach using a model 12-mer DNA duplex, free and in complex with the antibiotic echinomycin, which features two central adenine-thymine base pairs with Watson-Crick and Hoogsteen geometry, respectively, and subsequently extend it to the ∼200 kDa Widom 601 DNA nucleosome core particle. 
    more » « less
  4. ; ; ; ; ; (, RSC Chemical Biology)
    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
  5. ; ; (, Communications Chemistry)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email