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 more » 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. « less
Award ID(s):
Publication Date:
Journal Name:
Biophysical journal
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 pinpointedmore »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.« less
  2. 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 methylationmore »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.« less
  3. 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.
  4. 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.
  5. 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 andmore »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.« less