Hoogsteen Research Articles

KCNQ1OT1 promotes genome-wide transposon repression by guiding RNA-DNA triplexes and HP1 binding.

Transposon (de)repression and heterochromatin reorganization are dynamically regulated during cell fate determination and are hallmarks of cellular senescence. However, whether they are sequence specifically regulated remains unknown. Here we uncover that the KCNQ1OT1 lncRNA, by sequence-specific Hoogsteen base pairing with double-stranded genomic DNA via its repeat-rich region and binding to the heterochromatin protein HP1α, guides, induces and maintains epigenetic silencing at specific repetitive DNA elements. Repressing KCNQ1OT1 or deleting its repeat-rich region reduces DNA methylation and H3K9me3 on KCNQ1OT1-targeted transposons. Engineering a fusion KCNQ1OT1 with an ectopically targeting guiding triplex sequence induces de novo DNA methylation at the target site. Phenotypically, repressing KCNQ1OT1 induces senescence-associated heterochromatin foci, transposon activation and retrotransposition as well as cellular senescence, demonstrating an essential role of KCNQ1OT1 to safeguard against genome instability and senescence.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions.

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Alchemical Free-Energy Calculations of Watson-Crick and Hoogsteen Base Pairing Interconversion in DNA.

Hoogsteen (HG) base pairs have a transient nature and can be structurally similar to Watson-Crick (WC) base pairs, making their occurrence and thermodynamic stability difficult to determine experimentally. Herein, we employed the restrain-free-energy perturbation-release (R-FEP-R) method to calculate the relative free energy of the WC and HG base pairing modes in isolated and bound DNA systems and predict the glycosyl torsion conformational preference of purine bases. Notably, this method does not require prior knowledge of the transition pathway between the two end states. Remarkably, relatively fast convergence was reached, with results in excellent agreement with experimental data for all the examined DNA systems. The R-REP-R method successfully determined the stability of HG base pairing and more generally, the conformational preference of purine bases, in these systems. Therefore, this computational approach can help to understand the dynamic equilibrium between the WC and HG base pairing modes in DNA.

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.

Insight into achirality and chirality effects in interactions of an racemic ruthenium(II) polypyridyl complex and its Δ- and Λ-enantiomers with an RNA triplex

RNA triplexes have a variety of potential applications in molecular biology, diagnostics and therapeutics, while low stabilization of the third strand hinders their practical utilities under physiological conditions. In this regard, achieving the third-strand stabilization by binding small molecules is a promising strategy. Chirality is one of the basic properties of nature. To clarify achirality and chirality effects on the binding and stabilizing effects of RNA triplexes by small molecules, we report for the first time the RNA interactions of an racemic ruthenium(II) polypyridyl complex [Ru(bpy)2(11-CN-dppz)]2+ (rac-Ru1) and its two enantiomers Δ/Λ-[Ru(bpy)2(11-CN-dppz)]2+ (Δ/Λ-Ru1) with an RNA triplex poly(U-A*U) (where “-” represents Watson-Crick base pairing, and “*” denotes Hoogsteen base pairing, respectively) in this work. Research shows that although rac-Ru1 and its two enantiomers Δ/Λ-Ru1 bind to the RNA triplex through the same mode of intercalation, the binding affinity for enantiomer Δ-Ru1 is much higher than that for rac-Ru1 and enantiomer Λ-Ru1. However, compared to enantiomer Λ-Ru1, the binding affinity for rac-Ru1 does not show much of an advantage, which is slightly greater than that for the former. Thermal denaturation measurements reveal both rac-Ru1 and Δ-Ru1 to have a preference for stabilizing the third strand rather than the template duplex of the RNA triplex, while Λ-Ru1 stabilizes the RNA triplex without significant selectivity. Besides, the third-strand stabilizing effects by rac-Ru1 and Δ-Ru1 are not markedly different from each other, but more marked than that by Λ-Ru1. This work shows that the binding properties of the racemic Ru(II) polypyridyl complex with the RNA triplex are not simply an average of its two enantiomers, indicating potentially complicated binding events.

Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson–Crick base pair at the two drug-binding interfaces can result in a single-stranded ‘chair-shaped’ DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a ‘curved-shaped’ duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs’ potential synergetic use.

Probing Watson-Crick and Hoogsteen base pairing in duplex DNA using dynamic nuclear polarization solid-state NMR spectroscopy

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.

Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA.

Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.

Theoretical insights into water network of B-DNA duplex with Watson-Crick and Hoogsteen base pairing geometries

Solvation is essential for the proper structure and functioning of biomolecules. It has a major role in establishing the B-form as the dominant structure of DNA and strongly affects its biological functions. In our genome, two modes of base pairing (Watson-Crick (WC) and Hoogsteen (HG)) exist in equilibrium and are related to the functional properties of the DNA including ligand recognition and gene regulation. Here, we use molecular dynamics simulation along with topological analysis of the H-bond network to investigate the hydration patterns of DNA structure with WC and HG base pairings. We showed that DNA most significantly disturbs the H-bonding properties of the first hydration shell and slightly the second shell compared to bulk water and that the water network is considerably different in the minor and major grooves. We found that HG base pairing leads to increased solvation of the major groove and a less ordered water structure in the minor groove accompanied with increasing hydrophobicity of the minor groove which could also contribute to the site recognition of specific DNA binding proteins. Our newly developed method to selectively study the solvation of the minor and major grooves opens the door to study in finer details the binding between various form of DNA and different ligands such drugs, RNA and proteins.

Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting

Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C+ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.

1-Deazaguanosine-Modified RNA: The Missing Piece for Functional RNA Atomic Mutagenesis.

Atomic mutagenesis is the key to advance our understanding of RNA recognition and RNA catalysis. To this end, deazanucleosides are utilized to evaluate the participation of specific atoms in these processes. One of the remaining challenges is access to RNA-containing 1-deazaguanosine (c1G). Here, we present the synthesis of this nucleoside and its phosphoramidite, allowing first time access to c1G-modified RNA. Thermodynamic analyses revealed the base pairing parameters for c1G-modified RNA. Furthermore, by NMR spectroscopy, a c1G-triggered switch of Watson-Crick into Hoogsteen pairing in HIV-2 TAR RNA was identified. Additionally, using X-ray structure analysis, a guanine–phosphate backbone interaction affecting RNA fold stability was characterized, and finally, the critical impact of an active-site guanine in twister ribozyme on the phosphodiester cleavage was revealed. Taken together, our study lays the synthetic basis for c1G-modified RNA and demonstrates the power of the completed deazanucleoside toolbox for RNA atomic mutagenesis needed to achieve in-depth understanding of RNA recognition and catalysis.

Understanding the role of non‐Watson‐Crick base pairs in DNA–protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

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 A6-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.

Role of Hoogsteen interaction in the stability of different phases of triplex DNA.

A simple coarse-grained model of DNA which includes both Watson-Crick and Hoogsteen base pairing has been used to study the melting and unzipping of triplex DNA. Using Langevin dynamics simulations, we reproduce the qualitative features of one-step and two-step thermal melting of triplex as seen in experiments. The thermal melting phase diagram shows the existence of a stable interchain three-strand complex (bubble-bound state). Our studies based on the mechanical unzipping of a triplex revealed that it is mechanically more stable compared to an isolated duplex-DNA.

Curiosity: a tool for exploring density features in crystallographic and cryoEM maps

A key component of macromolecular structure determination from crystallographic or cryoEM density maps is the discovery and repair of modeling errors. In the case of chemically and physically unreasonable arrangements of atoms, we are informed by the clashscore, Ramachandran plots, bond distance/angle/torsion/chirality/planarity outliers, and so forth. Some tools for identifying chemically unlikely but biologically important structural motifs such as Hoogsteen base pairing have also been developed.

The "guiding" principles of noncoding RNA function.

The human genome is pervasively transcribed and yet only a small fraction of these RNAs (less than 2%) are known to code for proteins. The vast majority of the RNAs are classified as noncoding RNAs (ncRNAs) and are further subgrouped as small (shorter than 200 bases) and long noncoding RNAs. The ncRNAs have been identified in all three domains of life and regulate diverse cellular processes through transcriptional and posttranscriptional gene regulation. Most of these RNAs work in conjunction with proteins forming a wide array of base pairing interactions. The determinants of these base pairing interactions are now becoming more evident and show striking similarities among the diverse group of ncRNAs. Here we present a mechanistic overview of pairing between RNA-RNA or RNA-DNA that dictates the function of ncRNAs; we provide examples to illustrate that ncRNAs work through shared evolutionary mechanisms that encompasses a guide-target interaction, involving not only classical Watson-Crick but also noncanonical Wobble and Hoogsteen base pairing. We also highlight the similarities in target selection, proofreading, and the ruler mechanism of ncRNA-protein complexes that confers target specificity and target site selection. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA-Based Catalysis > RNA-Mediated Cleavage RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA.

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.

RNA:DNA triple helices: from peculiar structures to pervasive chromatin regulators.

The genomes of complex eukaryotes largely contain non-protein-coding DNA, which is pervasively transcribed into a plethora of non-coding RNAs (ncRNAs). The functional importance of many of these ncRNAs has been investigated in the last two decades, revealing their crucial and multifaceted roles in chromatin regulation. A common mode of action of ncRNAs is the recruitment of chromatin modifiers to specific regions in the genome. Whereas many ncRNA-protein interactions have been characterised in detail, binding of ncRNAs to their DNA target sites is much less understood. Recently developed RNA-centric methods have mapped the genome-wide distribution of ncRNAs, however, how ncRNAs achieve locus-specificity remains mainly unresolved. In terms of direct RNA-DNA interactions, two kinds of triple-stranded structures can be formed: R-loops consisting of an RNA:DNA hybrid and a looped out DNA strand, and RNA:DNA triple helices (triplexes), in which the RNA binds to the major groove of the DNA double helix by sequence-specific Hoogsteen base pairing. In this essay, we will review the current knowledge about RNA:DNA triplexes, summarising triplex formation rules, detection methods, and ncRNAs reported to engage in triplexes. While the functional characterisation of RNA:DNA triplexes is still anecdotal, recent advances in high-throughput and computational analyses indicate their widespread distribution in the genome. Thus, we are witnessing a paradigm shift in the appreciation of RNA:DNA triplexes, away from exotic structures towards a prominent mode of ncRNA-chromatin interactions.

Revealing the Hydrogen Bonding Interaction of DNA with Unnatural Bases via Plasmonic Antenna Enhanced Infrared Spectroscopy.

The incorporation of unnatural bases in DNA programming can break through the limits of Watson-Crick and Hoogsteen base pairing to expand the diversity of DNA structures. Thus, understanding the interaction between DNA and unnatural bases is of great importance in DNA nanotechnology. Here, we propose an approach of plasmonic antenna enhanced infrared spectroscopy to study the hydrogen bonding interaction between poly(thymine) DNA (poly T DNA) and melamine. The formation of multiple hydrogen bonds between melamine and thymine of poly T DNA is revealed by the appearance of a new infrared (IR) feature of the NH2 deformation vibration at 1680 cm-1. The binding rate constant (kb) and the dissociation rate constant (kd) of the affinity reaction reach 39.70 M-1·s-1 and 4.49 × 10-5 s-1, respectively. This work offers a valuable IR technique to study DNA nanostructures at the molecular level, providing unique physicochemical views of the interaction mechanism between DNA and unnatural bases in DNA programming.

Rapid assessment of Watson-Crick to Hoogsteen exchange in unlabeled DNA duplexes using high-power SELOPE imino 1H CEST.

In duplex DNA, Watson-Crick A-T and G-C base pairs (bp's) exist in dynamic equilibrium with an alternative Hoogsteen conformation, which is low in abundance and short-lived. Measuring how the Hoogsteen dynamics varies across different DNA sequences, structural contexts and physiological conditions is key for identifying potential Hoogsteen hot spots and for understanding the potential roles of Hoogsteen base pairs in DNA recognition and repair. However, such studies are hampered by the need to prepare C or N isotopically enriched DNA samples for NMR relaxation dispersion (RD) experiments. Here, using SELective Optimized Proton Experiments (SELOPE) H CEST experiments employing high-power radiofrequency fields ( 250 Hz) targeting imino protons, we demonstrate accurate and robust characterization of Watson-Crick to Hoogsteen exchange, without the need for isotopic enrichment of the DNA. For 13 residues in three DNA duplexes under different temperature and pH conditions, the exchange parameters deduced from high-power imino H CEST were in very good agreement with counterparts measured using off-resonance C N spin relaxation in the rotating frame (). It is shown that H-H NOE effects which typically introduce artifacts in H-based measurements of chemical exchange can be effectively suppressed by selective excitation, provided that the relaxation delay is short ( 100 ms). The H CEST experiment can be performed with 10 higher throughput and 100 lower cost relative to C N and enabled Hoogsteen chemical exchange measurements undetectable by . The results reveal an increased propensity to form Hoogsteen bp's near terminal ends and a diminished propensity within A-tract motifs. The H CEST experiment provides a basis for rapidly screening Hoogsteen breathing in duplex DNA, enabling identification of unusual motifs for more in-depth characterization.