skip to main content


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
NSF-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. 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
  2. null (Ed.)
    Abstract The N4-methylation of cytidine (m4C and m42C) in RNA plays important roles in both bacterial and eukaryotic cells. In this work, we synthesized a series of m4C and m42C modified RNA oligonucleotides, conducted their base pairing and bioactivity studies, and solved three new crystal structures of the RNA duplexes containing these two modifications. Our thermostability and X-ray crystallography studies, together with the molecular dynamic simulation studies, demonstrated that m4C retains a regular C:G base pairing pattern in RNA duplex and has a relatively small effect on its base pairing stability and specificity. By contrast, the m42C modification disrupts the C:G pair and significantly decreases the duplex stability through a conformational shift of native Watson-Crick pair to a wobble-like pattern with the formation of two hydrogen bonds. This double-methylated m42C also results in the loss of base pairing discrimination between C:G and other mismatched pairs like C:A, C:T and C:C. The biochemical investigation of these two modified residues in the reverse transcription model shows that both mono- or di-methylated cytosine bases could specify the C:T pair and induce the G to T mutation using HIV-1 RT. In the presence of other reverse transcriptases with higher fidelity like AMV-RT, the methylation could either retain the normal nucleotide incorporation or completely inhibit the DNA synthesis. These results indicate the methylation at N4-position of cytidine is a molecular mechanism to fine tune base pairing specificity and affect the coding efficiency and fidelity during gene replication. 
    more » « less
  3. null (Ed.)
    Abstract Regulatory protein access to the DNA duplex ‘interior’ depends on local DNA ‘breathing’ fluctuations, and the most fundamental of these are thermally-driven base stacking-unstacking interactions. The smallest DNA unit that can undergo such transitions is the dinucleotide, whose structural and dynamic properties are dominated by stacking, while the ion condensation, cooperative stacking and inter-base hydrogen-bonding present in duplex DNA are not involved. We use dApdA to study stacking-unstacking at the dinucleotide level because the fluctuations observed are likely to resemble those of larger DNA molecules, but in the absence of constraints introduced by cooperativity are likely to be more pronounced, and thus more accessible to measurement. We study these fluctuations with a combination of Molecular Dynamics simulations on the microsecond timescale and Markov State Model analyses, and validate our results by calculations of circular dichroism (CD) spectra, with results that agree well with the experimental spectra. Our analyses show that the CD spectrum of dApdA is defined by two distinct chiral conformations that correspond, respectively, to a Watson–Crick form and a hybrid form with one base in a Hoogsteen configuration. We find also that ionic structure and water orientation around dApdA play important roles in controlling its breathing fluctuations. 
    more » « less
  4. 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. 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