Computational Probing of Watson-Crick/Hoogsteen Breathing in a DNA Duplex Containing N1-Methylated Adenine

Isomerization Pathways of a Mismatched Base Pair of A:8OG in Free Duplex DNA

The A:8OG base pair (bp) is the outcome of DNA replication of the mismatched C:8OG bp. A high A:8OG bp population increases the C/G to A/T transversion mutation, which is responsible for various diseases. MutY is an important enzyme in the error-proof cycle and reverts A:8OG to C:8OG bp by cleaving adenine from the A:8OG bp. Several X-ray crystallography studies have determined the structure of MutY during the lesion scanning and lesion recognition stages. Interestingly, glycosidic bond (χ) angles of A:8OG bp in those two lesion recognition structures were found to differ, which implies that χ-torsion isomerization should occur during the lesion recognition process. In this study, as a first step to understanding this isomerization process, we characterized the intrinsic dynamic features of A:8OG in free DNAs by a free energy landscape simulation at the all-atom level. In this study, four isomerization states were assigned in the order of abundance: Aanti:8OGsyn > Aanti:8OGanti > Asyn:8OGanti ≈ Asyn:8OGsyn. Of these bp states, only 8OG in Asyn:8OGanti was located in the extrahelical space, whereas the purine bases (A and 8OG) in the other bp states remained inside the DNA helix. Also, free energy landscapes showed that the isomerization processes connecting these four bp states proceeded mostly in the intrahelical space via successive single glycosidic bond rotations of either A or 8OG.

Studies of DNA 'Breathing' by Polarization-Sweep Single-Molecule Fluorescence Microscopy of Exciton-Coupled (iCy3)2 Dimer-Labeled DNA Fork Constructs

Local fluctuations of the sugar-phosphate backbones and bases of DNA (often called DNA 'breathing') play a variety of critical roles in controlling the functional interactions of the DNA genome with the protein complexes that regulate it. Here, we present a single-molecule fluorescence method that we have used to measure and characterize such conformational fluctuations at and near biologically important positions in model DNA replication fork constructs labeled with exciton-coupled cyanine [(iCy3)2] dimer probes. Previous work has shown that the constructs that we tested here exhibit a broad range of spectral properties at the ensemble level, and these differences can be structurally and dynamically interpreted using our present methodology at the single-molecule level. The (iCy3)2 dimer has one symmetric (+) and one antisymmetric (-) exciton, with the respective transition dipole moments oriented perpendicular to one another. We excite single-molecule samples using a continuous-wave linearly polarized laser, with the polarization direction continuously rotated at the frequency of 1 MHz. The ensuing fluorescence signal is modulated as the laser polarization alternately excites the symmetric and antisymmetric excitons of the (iCy3)2 dimer probe. Phase-sensitive detection of the modulated signal provides information about the distribution of local conformations and the conformational interconversion dynamics of the (iCy3)2 probe. We find that at most construct positions that we examined, the (iCy3)2 dimer-labeled DNA fork constructs can adopt four topologically distinct conformational macrostates. These results suggest that in addition to observing DNA breathing at and near ss-dsDNA junctions, our new methodology should be useful to determine which of these pre-existing macrostates are recognized by, bind to, and are stabilized by various genome-regulatory proteins.

Conformational stability and order of Hoogsteen base pair induced by protein binding

Several experimental studies have shown that Hoogsteen (HG) base pair (bp) stabilizes in the presence of proteins. The molecular mechanism underlying this stabilization is not well known. This leads us to examine the stability of the HG bp in duplex DNA using all-atom molecular dynamics simulation in both the absence and presence of proteins. We use conformational thermodynamics to investigate the stability of a HG bp in duplex DNA at the molecular level. We compute the changes in the conformational free energy and entropy of DNA when DNA adopts a HG bp in its bp sequence rather than a Watson-Crick (WC) bp in both naked DNA and protein-bound DNA complex. We observe that the presence of proteins stabilizes and organizes the HG bp and the entire DNA duplex. Sugar-phosphate, sugar-base, and sugar-pucker torsion angles play key roles in stabilizing and ordering the HG bp in the protein-bound DNA complex.

Force-Field-Dependent DNA Breathing Dynamics: A Case Study of Hoogsteen Base Pairing in A6-DNA

The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.

Sequence dependence of transient Hoogsteen base pairing in DNA

Hoogsteen (HG) base pairing is characterized by a 180° rotation of the purine base with respect to the Watson-Crick-Franklin (WCF) motif. Recently, it has been found that both conformations coexist in a dynamical equilibrium and that several biological functions require HG pairs. This relevance has motivated experimental and computational investigations of the base-pairing transition. However, a systematic simulation of sequence variations has remained out of reach. Here, we employ advanced path-based methods to perform unprecedented free-energy calculations. Our methodology enables us to study the different mechanisms of purine rotation, either remaining inside or after flipping outside of the double helix. We study seven different sequences, which are neighbor variations of a well-studied A⋅T pair in A₆-DNA. We observe the known effect of A⋅T steps favoring HG stability, and find evidence of triple-hydrogen-bonded neighbors hindering the inside transition. More importantly, we identify a dominant factor: the direction of the A rotation, with the 6-ring pointing either towards the longer or shorter segment of the chain, respectively relating to a lower or higher barrier. This highlights the role of DNA’s relative flexibility as a modulator of the WCF/HG dynamic equilibrium. Additionally, we provide a robust methodology for future HG proclivity studies.

Recently, an alternative DNA base-pairing conformation, known as Hoogsteen (HG), has been found to coexist with the well-known Watson-Crick-Franklin (WCF) pairing. Several experimental and computational studies have focused on this heterogeneity, as it is involved in various recognition and replication processes. The WCF-to-HG transition mechanisms consist of a ±180° rotation of the purine base, occurring either inside of the double helix or while flipping temporarily outside of it. Even though molecular dynamics simulations can provide fine details about the transition pathways and their free-energy barriers, the computational cost has limited most studies to focus on only one particular chain (A₆-DNA). Here, we investigate the sequence dependence of the base-pairing transition, by systematically varying the direct neighbors of a transitioning A⋅T pair; probing inside and outside pathways in seven distinct systems. We discover that triple-hydrogen-bonded neighboring base-pairs hinder the inside rotation mechanism, due to the reduced flexibility needed for internal base rotation. Across all sequences, we confirm that outside transitions have a lower free-energy barrier. Most importantly, we observe that the direction of the A rotation, with the A 6-ring pointing either towards the long or short end of the modelled DNA chain, has a determinant effect on the height of the free-energy barrier. These results point to a critical role of DNA’s small- and medium-scale flexibility in modulating the proclivity of HG base pairs; providing a handle that might be employed by several biological mechanisms.

Free-Energy Landscape of a pH-Modulated G·C Base Pair Transition from Watson-Crick to Hoogsteen State in Duplex DNA

In double-helical DNAs, the most stable Watson-Crick (WC) base pair (bp) can be in thermal equilibrium with much less abundant Hoogsteen (HG) bp by the spontaneous rotation of the glycosidic angle in purine bases. Previous experimental studies showed that in the case of a G·C bp, the population of the transient HG is enhanced as a protonated form (HG⁺) through the protonation of the cytosine base under weakly acidic conditions. Hence, pH is a key factor that can modulate this transition event from the WC to HG⁺ bp. In this study, to computationally probe the overall free-energy landscapes of this pH-modulated G·C HG breathing, a comprehensive classical molecular dynamics (MD) simulation protocol is proposed using an enhanced sampling MD in conjunction with the standard thermodynamic integration method. From this MD protocol proposed, the free-energy surfaces of the G·C bp transition from the WC to HG bp were constructed successfully at any pH range, producing pH-dependent free-energy quantities in close agreement with previously reported experimental results. The simulation protocol is expected to provide valuable atomistic insight into the DNA bp transition events coupled with protonation or tautomeric shift in a target bp.

Effects of Adenine Methylation on the Structure and Thermodynamic Stability of a DNA Minidumbbell

DNA methylation is a prevalent regulatory modification in prokaryotes and eukaryotes. N1-methyladenine (m1A) and N6-methyladenine (m6A) have been found to be capable of altering DNA structures via disturbing Watson-Crick base pairing. However, little has been known about their influences on non-B DNA structures, which are associated with genetic instabilities. In this work, we investigated the effects of m1A and m6A on both the structure and thermodynamic stability of a newly reported DNA minidumbbell formed by two TTTA tetranucleotide repeats. As revealed by the results of nuclear magnetic resonance spectroscopic studies, both m1A and m6A favored the formation of a T·m1A and T·m6A Hoogsteen base pair, respectively. More intriguingly, the m1A and m6A modifications brought about stabilization and destabilization effects on the DNA minidumbbell, respectively. This work provides new biophysical insights into the effects of adenine methylation on the structure and thermodynamic stability of DNA.

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.