New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Quantitative and Systematic NMR Measurements of Sequence-Dependent A-T Hoogsteen Dynamics in the DNA Double Helix

The dynamic properties of DNA depend on the sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging, especially when they involve lowly populated and short-lived conformational states. Using 1H CEST supplemented by targeted 13C R1ρ NMR experiments, we quantitatively measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in 13 trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. Our results uncover a unique conformational basis for sequence-specificity in the DNA double helix and establish the utility of NMR to quantitatively and comprehensively measure sequence-dependent DNA dynamics.

Binding of Homeodomain Proteins to DNA with Hoogsteen Base Pair

In DNA double helices, Hoogsteen (HG) base pairing is an alternative mode of Watson-Crick (WC) base pairing. HG bp has a different hydrogen bonding pattern than WC bp. We investigate here the binding energy of homeodomain proteins with a HG-DNA duplex, where DNA adopts a HG bp in its sequence. We observe that the presence of the HG bp increases the binding energy of both the specific and nonspecific homeodomain proteins compared to WC-DNA bps. The neutral mutation in the N-terminal basic arm of the nonspecific protein significantly changes the binding energy between nonspecific protein and HG-DNA only, while the acidic mutation significantly changes the binding energy of both the specific and nonspecific proteins with HG-DNA. The significant variation in the binding energy of the homeodomains within distinct DNA-protein complexes can be ascribed to the differences in the number of intermolecular contacts between DNA bases and protein residues. Our conformational thermodynamics calculations based on the fluctuation of microscopic conformational variables at the interface show that with increasing conformational stability and order at the interface, the binding of the homeodomain protein gets stronger.

Nanopore Discriminates Watson-Crick and Hoogsteen Hydrogen Bonds in Multiple DNA Contexts

Watson-Crick and Hoogsteen hydrogen bonds aid the formation of highly ordered structures in genomic DNA that dynamically govern genetic modes such as gene regulation and replication. Hence, measuring and distinguishing these two types of hydrogen bonds in different DNA contexts are essential for understanding DNA architectures. However, due to their transient nature and minimal structure differences at the sub-nanometer scale, differentiating Watson-Crick hydrogen bonds from Hoogsteen hydrogen bonds is difficult. Relying on nanopore technology, we successfully discriminated the two types of hydrogen bonds in multiple DNA contexts in the presence of epigenetic modification, changes in DNA structures, and proton strength in the environment. Our results indicate that Watson-Crick and Hoogsteen hydrogen bonds show different susceptibilities to changes in physicochemical characteristics that matter in stabilizing DNA hydrogen bonds. This work provides insight into the features of hydrogen bonds at the nanoscale and may benefit profiling complex DNA architectures by measuring subtle structural changes.

CRISPR-Cas12a bends DNA to destabilize base pairs during target interrogation

RNA-guided endonucleases are involved in processes ranging from adaptive immunity to site-specific transposition and have revolutionized genome editing. CRISPR-Cas9, -Cas12 and related proteins use guide RNAs to recognize ∼20-nucleotide target sites within genomic DNA by mechanisms that are not yet fully understood. We used structural and biochemical methods to assess early steps in DNA recognition by Cas12a protein-guide RNA complexes. We show here that Cas12a initiates DNA target recognition by bending DNA to induce transient nucleotide flipping that exposes nucleobases for DNA-RNA hybridization. Cryo-EM structural analysis of a trapped Cas12a-RNA-DNA surveillance complex and fluorescence-based conformational probing show that Cas12a-induced DNA helix destabilization enables target discovery and engagement. This mechanism of initial DNA interrogation resembles that of CRISPR-Cas9 despite distinct evolutionary origins and different RNA-DNA hybridization directionality of these enzyme families. Our findings support a model in which RNA-mediated DNA interference begins with local helix distortion by transient CRISPR-Cas protein binding.

Selective Recognition of the Dimeric NG16 Parallel G-Quadruplex Structure Using Synthetic Turn-On Red Fluorescent Protein Chromophore

Red fluorescent protein (RFP)-based fluorescent probes that can selectively interact with specific nucleic acids are of great importance for therapeutic and bioimaging applications. Herein, we have reported the synthesis of RFP chromophores for selective recognition of G-quadruplex nucleic acids in vitro and ex vivo. We identified DFHBI-DM as a fluorescent turn-on probe that binds to the dimeric NG16 parallel quadruplex with superior selectivity and sensitivity over various parallel, antiparallel, and hybrid topologies. The binding of DFHBI-DM to NG16 exhibited excellent photophysical properties, including high binding affinity, large Stokes shift, high photostability, and quantum yield. The MD simulation study supports the 1:1 binding stoichiometry. It confirms the planar conformation of DFHBI-DM, which makes strong binding interactions with a flat quartet of NG16 compared to other antiparallel and hybrid topologies. The cell imaging and MTT assays revealed that DFHBI-DM is a biocompatible and efficient fluorescent probe for intracellular imaging of NG16. Overall, these results demonstrated that DFHBI-DM could be an effective fluorescent G4-stabilizing agent for the dimeric NG16 parallel quadruplex, and it could be a promising candidate for further exploration of bioimaging and therapeutic applications.

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

Structural basis of water-mediated cis Watson-Crick/Hoogsteen base-pair formation in non-CpG methylation

Non-CpG methylation is associated with several cellular processes, especially neuronal development and cancer, while its effect on DNA structure remains unclear. We have determined the crystal structures of DNA duplexes containing -CGCCG- regions as CCG repeat motifs that comprise a non-CpG site with or without cytosine methylation. Crystal structure analyses have revealed that the mC:G base-pair can simultaneously form two alternative conformations arising from non-CpG methylation, including a unique water-mediated cis Watson-Crick/Hoogsteen, (w)cWH, and Watson-Crick (WC) geometries, with partial occupancies of 0.1 and 0.9, respectively. NMR studies showed that an alternative conformation of methylated mC:G base-pair at non-CpG step exhibits characteristics of cWH with a syn-guanosine conformation in solution. DNA duplexes complexed with the DNA binding drug echinomycin result in increased occupancy of the (w)cWH geometry in the methylated base-pair (from 0.1 to 0.3). Our structural results demonstrated that cytosine methylation at a non-CpG step leads to an anti→syntransition of its complementary guanosine residue toward the (w)cWH geometry as a partial population of WC, in both drug-bound and naked mC:G base pairs. This particular geometry is specific to non-CpG methylated dinucleotide sites in B-form DNA. Overall, the current study provides new insights into DNA conformation during epigenetic regulation.

CRISPR-Cas12a bends DNA to destabilize base pairs during target interrogation

RNA-guided endonucleases are involved in processes ranging from adaptive immunity to site-specific transposition and have revolutionized genome editing. CRISPR-Cas9, -Cas12 and related proteins use guide RNAs to recognize ~20-nucleotide target sites within genomic DNA by mechanisms that are not yet fully understood. We used structural and biochemical methods to assess early steps in DNA recognition by Cas12a protein-guide RNA complexes. We show here that Cas12a initiates DNA target recognition by bending DNA to induce transient nucleotide flipping that exposes nucleobases for DNA-RNA hybridization. Cryo-EM structural analysis of a trapped Cas12a-RNA-DNA surveillance complex and fluorescence-based conformational probing show that Cas12a-induced DNA helix destabilization enables target discovery and engagement. This mechanism of initial DNA interrogation resembles that of CRISPR-Cas9 despite distinct evolutionary origins and different RNA-DNA hybridization directionality of these enzyme families. Our findings support a model in which RNA-mediated DNA engineering begins with local helix distortion by transient CRISPR-Cas protein binding.

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

High-Resolution NMR Structures of Intrastrand Hairpins Formed by CTG Trinucleotide Repeats

The CAG and CTG trinucleotide repeat expansions cause more than 10 human neurodegenerative diseases. Intrastrand hairpins formed by trinucleotide repeats contribute to repeat expansions, establishing them as potential drug targets. High-resolution structural determination of CAG and CTG hairpins poses as a long-standing goal to aid drug development, yet it has not been realized due to the intrinsic conformational flexibility of repetitive sequences. We herein investigate the solution structures of CTG hairpins using nuclear magnetic resonance (NMR) spectroscopy and found that four CTG repeats with a clamping G-C base pair was able to form a stable hairpin structure. We determine the first solution NMR structure of dG(CTG)4C hairpin and decipher a type I folding geometry of the TGCT tetraloop, wherein the two thymine residues form a T·T loop-closing base pair and the first three loop residues continuously stack. We further reveal that the CTG hairpin can be bound and stabilized by a small-molecule ligand, and the binding interferes with replication of a DNA template containing CTG repeats. Our determined high-resolution structures lay an important foundation for studying molecular interactions between native CTG hairpins and ligands, and benefit drug development for trinucleotide repeat expansion diseases.

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.

Hybrid nucleobases as new and efficient unnatural genetic letters

To expand the existing genetic letters beyond the natural four nucleotides, such as G, C, A, and T, it is necessary to design robust nucleotides that can not only produce stable and unperturbed DNA but also function naturally in living cells. Although hydrophobic bases, such as d5SICS (2,6-dimethyl-2H-isoquiniline-1-thione) and dNaM (2-methoxy-3-methylnaphthalene) were shown to be replicated in bacterial cells, the d5SICS:dNaM base-pair was found to perturb the structure of the duplex DNA. Therefore, it is necessary to design nucleobases that can form base pairs like the natural G:C and A:T pairs. Here, a reliable dispersion-corrected density functional theory has been used to design several nucleobases that can produce three-hydrogen-bonded base pairs like the G:C pair. In doing so, the Watson-Crick faces of d5SICS and dNaM were modified by replacing the hydrophobic groups with hydrogen bond donors and acceptors. As dNaM contains an unnatural C-glycosidic bond (C-dNaM), it was also modified to contain the natural N-glycosidic bond (N-dNaM). This technique produced 91 new bases (N-d5SICS-X (X = 1-33), C-dNaM-X (X = 1-35), and N-dNaM-X (X = 1-23), where X is the different types of modifications applied to d5SICS and dNaM) and 259 base-pairs. Among these base pairs, 76 base pairs are found to be more stable than the G:C pair. Interestingly, the N-d5SICS-32:C-dNaM-32 and N-d5SICS-32:N-dNaM-20 pairs are found to be the most stable with binding energies of about -28.0 kcal/mol. The base-pair patterns of these pairs are also analogous to that of the G:C pair. Hence, it is proposed that N-d5SICS-32, C-dNaM-32, and N-dNaM-20 would act as efficient new genetic letters to produce stable and unperturbed artificial DNA.Communicated by Ramaswamy H. Sarma.

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.

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.

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.

Bacterial Biofilms Utilize an Underlying Extracellular DNA Matrix Structure That Can Be Targeted for Biofilm Resolution

Bacterial biofilms contribute significantly to the antibiotic resistance, pathogenesis, chronicity and recurrence of bacterial infections. Critical to the stability and survival of extant biofilms is the extracellular DNA (eDNA)-dependent matrix which shields the resident bacteria from hostile environments, allows a sessile metabolic state, but also encourages productive interactions with biofilm-inclusive bacteria. Given the importance of the eDNA, approaches to this area of research have been to target not just the eDNA, but also the additional constituent structural components which appear to be widespread. Chief among these is a ubiquitous two-member family of bacterial nucleoid associated proteins (the DNABII proteins) responsible for providing structural integrity to the eDNA and thereby the biofilm. Moreover, this resultant novel eDNA-rich secondary structure can also be targeted for disruption. Here, we provide an overview of both what is known about the eDNA-dependent matrix, as well as the resultant means that have resulted in biofilm resolution. Results obtained to date have been highly supportive of continued development of DNABII-targeted approaches, which is encouraging given the great global need for improved methods to medically manage, or ideally prevent biofilm-dependent infections, which remains a highly prevalent burden worldwide.

Nucleic Acids and Their Analogues for Biomedical Applications

Nucleic acids are emerging as powerful and functional biomaterials due to their molecular recognition ability, programmability, and ease of synthesis and chemical modification. Various types of nucleic acids have been used as gene regulation tools or therapeutic agents for the treatment of human diseases with genetic disorders. Nucleic acids can also be used to develop sensing platforms for detecting ions, small molecules, proteins, and cells. Their performance can be improved through integration with other organic or inorganic nanomaterials. To further enhance their biological properties, various chemically modified nucleic acid analogues can be generated by modifying their phosphodiester backbone, sugar moiety, nucleobase, or combined sites. Alternatively, using nucleic acids as building blocks for self-assembly of highly ordered nanostructures would enhance their biological stability and cellular uptake efficiency. In this review, we will focus on the development and biomedical applications of structural and functional natural nucleic acids, as well as the chemically modified nucleic acid analogues over the past ten years. The recent progress in the development of functional nanomaterials based on self-assembled DNA-based platforms for gene regulation, biosensing, drug delivery, and therapy will also be presented. We will then summarize with a discussion on the advanced development of nucleic acid research, highlight some of the challenges faced and propose suggestions for further improvement.

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.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Discovering Collective Variables of Molecular Transitions via Genetic Algorithms and Neural Networks

With the continual improvement of computing hardware and algorithms, simulations have become a powerful tool for understanding all sorts of (bio)molecular processes. To handle the large simulation data sets and to accelerate slow, activated transitions, a condensed set of descriptors, or collective variables (CVs), is needed to discern the relevant dynamics that describes the molecular process of interest. However, proposing an adequate set of CVs that can capture the intrinsic reaction coordinate of the molecular transition is often extremely difficult. Here, we present a framework to find an optimal set of CVs from a pool of candidates using a combination of artificial neural networks and genetic algorithms. The approach effectively replaces the encoder of an autoencoder network with genes to represent the latent space, i.e., the CVs. Given a selection of CVs as input, the network is trained to recover the atom coordinates underlying the CV values at points along the transition. The network performance is used as an estimator of the fitness of the input CVs. Two genetic algorithms optimize the CV selection and the neural network architecture. The successful retrieval of optimal CVs by this framework is illustrated at the hand of two case studies: the well-known conformational change in the alanine dipeptide molecule and the more intricate transition of a base pair in B-DNA from the classic Watson-Crick pairing to the alternative Hoogsteen pairing. Key advantages of our framework include the following: optimal interpretable CVs, avoiding costly calculation of committor or time-correlation functions, and automatic hyperparameter optimization. In addition, we show that applying a time-delay between the network input and output allows for enhanced selection of slow variables. Moreover, the network can also be used to generate molecular configurations of unexplored microstates, for example, for augmentation of the simulation data.

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.

Dinucleotides as simple models of the base stacking-unstacking component of DNA 'breathing' mechanisms

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.

Slow motions in A·T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R_1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using ¹³C/¹⁵ N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

Rapid and accurate determination of atomistic RNA dynamic ensemble models using NMR and structure prediction

Biomolecules form dynamic ensembles of many inter-converting conformations which are key for understanding how they fold and function. However, determining ensembles is challenging because the information required to specify atomic structures for thousands of conformations far exceeds that of experimental measurements. We addressed this data gap and dramatically simplified and accelerated RNA ensemble determination by using structure prediction tools that leverage the growing database of RNA structures to generate a conformation library. Refinement of this library with NMR residual dipolar couplings provided an atomistic ensemble model for HIV-1 TAR, and the model accuracy was independently supported by comparisons to quantum-mechanical calculations of NMR chemical shifts, comparison to a crystal structure of a substate, and through designed ensemble redistribution via atomic mutagenesis. Applications to TAR bulge variants and more complex tertiary RNAs support the generality of this approach and the potential to make the determination of atomic-resolution RNA ensembles routine.

Determining dynamic ensembles of biomolecules is still challenging. Here the authors present an approach for rapid RNA ensemble determination that combines RNA structure prediction tools and NMR residual dipolar coupling data and use it to determine atomistic ensemble models for a variety of RNAs.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

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.

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC⁺ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.

DNA-specific selectivity in pairing of model nucleobases in the solid state

Solid-state may serve as the reaction medium for selective recognition between model nucleobases.

Atomistic insight into the kinetic pathways for Watson-Crick to Hoogsteen transitions in DNA

DNA predominantly contains Watson-Crick (WC) base pairs, but a non-negligible fraction of base pairs are in the Hoogsteen (HG) hydrogen bonding motif at any time. In HG, the purine is rotated ∼180° relative to the WC motif. The transitions between WC and HG may play a role in recognition and replication, but are difficult to investigate experimentally because they occur quickly, but only rarely. To gain insight into the mechanisms for this process, we performed transition path sampling simulations on a model nucleotide sequence in which an AT pair changes from WC to HG. This transition can occur in two ways, both starting with loss of hydrogen bonds in the base pair, followed by rotation around the glycosidic bond. In one route the adenine base converts from WC to HG geometry while remaining entirely within the double helix. The other route involves the adenine leaving the confines of the double helix and interacting with water. Our results indicate that this outside route is more probable. We used transition interface sampling to compute rate constants and relative free energies for the transitions between WC and HG. Our results agree with experiments, and provide highly detailed insights into the mechanisms of this important process.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements

NMR relaxation dispersion studies have shown that Watson-Crick G-C and A-T base pairs in duplex DNA exist in dynamic equilibrium with their Hoogsteen counterparts. Hoogsteen base pairs form through concurrent rotation of the purine base about the glycosidic bond from an anti to a syn conformation and constriction of the C1'-C1' distance across the base pair by ∼2 Å to allow Hoogsteen type hydrogen bonding. Owing to their unique structure, Hoogsteen base pairs can play important roles in DNA recognition, the accommodation, recognition, and repair of DNA damage, and in DNA replication. NMR relaxation dispersion experiments targeting imino nitrogen and protonated base and sugar carbons have provided insights into many structural features of transient Hoogsteen base pairs, including one of two predicted hydrogen bonds involving (G)N7···H-N3(C)⁺ and (A)N7···H-N3(T). Here, through measurement of cytosine amino (N4) R_1ρ relaxation dispersion, we provide direct evidence for the second (G)O6···H₂-N4(C)⁺ hydrogen bond in G(syn)-C⁺ transient Hoogsteen base pairs. The utility of cytosine N4 R_1ρ relaxation dispersion as a new sensitive probe of transient Hoogsteen base pairs, and cytosine dynamics in general, is further demonstrated by measuring G(syn)-C⁺ Hoogsteen exchange near neutral pH and in the context of the naturally occurring DNA modification 5-methyl cytosine (m⁵C), in DNA samples prepared using chemical synthesis and a ¹⁵N labeled m⁵C phosphoramidite.

Infrared Spectroscopic Observation of a G-C+ Hoogsteen Base Pair in the DNA:TATA-Box Binding Protein Complex Under Solution Conditions

Hoogsteen DNA base pairs (bps) are an alternative base pairing to canonical Watson-Crick bps and are thought to play important biochemical roles. Hoogsteen bps have been reported in a handful of X-ray structures of protein-DNA complexes. However, there are several examples of Hoogsteen bps in crystal structures that form Watson-Crick bps when examined under solution conditions. Furthermore, Hoogsteen bps can sometimes be difficult to resolve in DNA:protein complexes by X-ray crystallography due to ambiguous electron density and by solution-state NMR spectroscopy due to size limitations. Here, using infrared spectroscopy, we report the first direct solution-state observation of a Hoogsteen (G-C+ ) bp in a DNA:protein complex under solution conditions with specific application to DNA-bound TATA-box binding protein. These results support a previous assignment of a G-C+ Hoogsteen bp in the complex, and indicate that Hoogsteen bps do indeed exist under solution conditions in DNA:protein complexes.

Bending of DNA duplexes with mutation motifs

Mutations can be induced by environmental factors but also arise spontaneously during DNA replication or due to deamination of methylated cytosines at CpG dinucleotides. Sites where mutations occur with higher frequency than would be expected by chance are termed hotspots while sites that contain mutations rarely are termed coldspots. Mutations are permanently scanned and repaired by repair systems. Among them, the mismatch repair targets base pair mismatches, which are discriminated from canonical base pairs by probing altered elasticity of DNA. Using biased molecular dynamics simulations, we investigated the elasticity of coldspots and hotspots motifs detected in human genes associated with inherited disorders, and also of motifs with Czech population hotspots and de novo mutations. Main attention was paid to mutations leading to G/T and A+/C pairs. We observed that hotspots without CpG/CpHpG sequences are less flexible than coldspots, which indicates that flexible sequences are more effectively repaired. In contrary, hotspots with CpG/CpHpG sequences exhibited increased flexibility as coldspots. Their mutability is more likely related to spontaneous deamination of methylated cytosines leading to C > T mutations, which are primarily targeted by base excision repair. We corroborated conclusions based on computer simulations by measuring melting curves of hotspots and coldspots containing G/T mismatch.

Why are Hoogsteen base pairs energetically disfavored in A-RNA compared to B-DNA?

A(syn)-U/T and G(syn)-C⁺ Hoogsteen (HG) base pairs (bps) are energetically more disfavored relative to Watson–Crick (WC) bps in A-RNA as compared to B-DNA by >1 kcal/mol for reasons that are not fully understood. Here, we used NMR spectroscopy, optical melting experiments, molecular dynamics simulations and modified nucleotides to identify factors that contribute to this destabilization of HG bps in A-RNA. Removing the 2′-hydroxyl at single purine nucleotides in A-RNA duplexes did not stabilize HG bps relative to WC. In contrast, loosening the A-form geometry using a bulge in A-RNA reduced the energy cost of forming HG bps at the flanking sites to B-DNA levels. A structural and thermodynamic analysis of purine-purine HG mismatches reveals that compared to B-DNA, the A-form geometry disfavors syn purines by 1.5–4 kcal/mol due to sugar-backbone rearrangements needed to sterically accommodate the syn base. Based on MD simulations, an additional penalty of 3–4 kcal/mol applies for purine-pyrimidine HG bps due to the higher energetic cost associated with moving the bases to form hydrogen bonds in A-RNA versus B-DNA. These results provide insights into a fundamental difference between A-RNA and B-DNA duplexes with important implications for how they respond to damage and post-transcriptional modifications.

Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA

A( syn)-T and G( syn)-C⁺ Hoogsteen base pairs in protein-bound DNA duplexes can be difficult to resolve by X-ray crystallography due to ambiguous electron density and by nuclear magnetic resonance (NMR) spectroscopy due to poor chemical shift dispersion and size limitations with solution-state NMR spectroscopy. Here we describe an NMR strategy for characterizing Hoogsteen base pairs in protein-DNA complexes, which relies on site-specifically incorporating ¹³C- and ¹⁵N-labeled nucleotides into DNA duplexes for unambiguous resonance assignment and to improve spectral resolution. The approach was used to resolve the conformation of an A-T base pair in a crystal structure of an ∼43 kDa complex between a 34 bp duplex DNA and the integration host factor (IHF) protein. In the crystal structure (Protein Data Bank entry 1IHF ), this base pair adopts an unusual Hoogsteen conformation with a distorted sugar backbone that is accommodated by a nearby nick used to aid in crystallization. The NMR chemical shifts and interproton nuclear Overhauser effects indicate that this base pair predominantly adopts a Watson-Crick conformation in the intact DNA-IHF complex under solution conditions. Consistent with these NMR findings, substitution of 7-deazaadenine at this base pair resulted in only a small (∼2-fold) decrease in the IHF-DNA binding affinity. The NMR strategy provides a new approach for resolving crystallographic ambiguity and more generally for studying the structure and dynamics of protein-DNA complexes in solution.

The Adaptive Path Collective Variable: A Versatile Biasing Approach to Compute the Average Transition Path and Free Energy of Molecular Transitions

In the past decade, great progress has been made in the development of enhanced sampling methods, aimed at overcoming the time-scale limitations of molecular dynamics (MD) simulations. Many sampling schemes rely on adding an external bias to favor the sampling of transitions and to estimate the underlying free energy landscape. Nevertheless, sampling molecular processes described by many order parameters, or collective variables (CVs), such as complex biomolecular transitions, remains often very challenging. The computational cost has a prohibitive scaling with the dimensionality of the CV-space. Inspiration can be taken from methods that focus on localizing transition pathways: the CV-space can be projected onto a path-CV that connects two stable states, and a bias can be exerted onto a one-dimensional parameter that captures the progress of the transition along the path-CV. In principle, such a sampling scheme can handle an arbitrarily large number of CVs. A standard enhanced sampling technique combined with an adaptive path-CV can then locate the mean transition pathway and obtain the free energy profile along the path. In this chapter, we discuss the adaptive path-CV formalism and its numerical implementation. We apply the path-CV with several enhanced sampling methods-steered MD, metadynamics, and umbrella sampling-to a biologically relevant process: the Watson-Crick to Hoogsteen base-pairing transition in double-stranded DNA. A practical guide is provided on how to recognize and circumvent possible pitfalls during the calculation of a free energy landscape that contains multiple pathways. Examples are presented on how to perform enhanced sampling simulations using PLUMED, a versatile plugin that can work with many popular MD engines.

Dual coding potential of a 2',5'-branched ribonucleotide in DNA

Branchpoints in RNA templates are highly mutagenic, but it is not known yet whether this also applies to branchpoints in DNA templates. Here, we report how nucleic acid polymerases replicate a 2',5'-branched DNA (bDNA) molecule. We constructed long-chained bDNA templates containing a branch guanosine and T7 promoters at both arms by splinted ligation. Quantitative real-time PCR analysis was used to investigate whether a branchpoint blocks DNA synthesis from the two arms in the same manner. We find that the blocking effect of a branchpoint is arm-specific. DNA synthesis from the 2'-arm is more than 20,000-fold decreased, whereas from the 3'-arm only 15-fold. Our sequence analysis of full-length nucleic acid generated by Taq DNA polymerase, Moloney murine leukemia virus reverse transcriptase, and T7 RNA polymerase from the 2'-arm of bDNA shows that the branched guanine has a dual coding potential and can base-pair with cytosine and guanine. We find that branchpoint templating is influenced by the type of the surrounding nucleic acid and is probably modulated by polymerase and RNase H active sites. We show that the branchpoint bypass by the polymerases from the 3'-arm of bDNA is predominantly error-free, indicating that bDNA is not as highly mutagenic as 2',5'-branched RNA.

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

DNA breathing is a local conformational fluctuation spontaneously occurring in double-stranded DNAs. In particular, the possibility of individual base pairs (bps) in duplex DNA to flip between alternate bp modes, i.e., Watson-Crick (WC)-like and Hoogsteen (HG)-like, at relevant time scales has impacted DNA research fields for many years. In this study, to computationally probe effects of chemical modification on the DNA breathing, we present a free energy landscape of spontaneous thermal transitions between WC and HG bps in a free DNA duplex containing N1-methylated adenine (m1A). For the current free energy computation, a variant of well-tempered metadynamics simulation was extensively performed for a total of 40 μs to produce free energy surfaces. The free energy profile indicated that, upon the chemical modification of adenine, the HG bp (m1A·T) was located in the most favorable conformation (96.7%); however, the canonical WC bp (m1A·T) was distorted into two WC-like bps of WC* (2.8%) and WC** (0.5%). The conformational exchange between these two minor WC-like bps occurs with the first hundred nanoseconds. The transition between WC-like and HG bp features multiple transition pathways displaying various extents of base flipping in combination with glycosidic rotation. Analysis of the simulated ensemble showed that the m1A-induced changes of the backbone and sugar pucker were in a reasonable agreement with previous results inferred from NMR experiments. Also, this study revealed that the formation of the stable HG bp upon the mutation alters the characteristics of dynamic fluctuations of the neighboring WC residues of m1A. We expect this simulation approach to be a robust computational scheme to complement and guide future high-resolution experiments on many outstanding issues of duplex DNA breathing.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Dynamics of physiologically relevant noncanonical DNA structures: an overview from experimental and theoretical studies

DNA is a complex molecule with phenomenal inherent plasticity and the ability to form different hydrogen bonding patterns of varying stabilities. These properties enable DNA to attain a variety of structural and conformational polymorphic forms. Structurally, DNA can exist in single-stranded form or as higher-order structures, which include the canonical double helix as well as the noncanonical duplex, triplex and quadruplex species. Each of these structural forms in turn encompasses an ensemble of dynamically heterogeneous conformers depending on the sequence composition and environmental context. In vivo, the widely populated canonical B-DNA attains these noncanonical polymorphs during important cellular processes. While several investigations have focused on the structure of these noncanonical DNA, studying their dynamics has remained nontrivial. Here, we outline findings from some recent advanced experimental and molecular simulation techniques that have significantly contributed toward understanding the complex dynamics of physiologically relevant noncanonical forms of DNA.

New Insights into the Role of DNA Shape on Its Recognition by p53 Proteins

The tumor suppressor p53 acts as a transcription factor recognizing diverse DNA response elements (REs). Previous structural studies of p53-DNA complexes revealed non-canonical Hoogsteen geometry of A/T base pairs at conserved CATG motifs leading to changes in DNA shape and its interface with p53. To study the effects of DNA shape on binding characteristics, we designed REs with modified base pairs "locked" into either Hoogsteen or Watson-Crick form. Here we present crystal structures of these complexes and their thermodynamic and kinetic parameters, demonstrating that complexes with Hoogsteen base pairs are stabilized relative to those with all-Watson-Crick base pairs. CATG motifs are abundant in p53REs such as GADD45 and p53R2 related to cell-cycle arrest and DNA repair. The high-resolution structures of these complexes validate their propensity to adopt the unique Hoogsteen-induced structure, thus providing insights into the functional role of DNA shape and broadening the mechanisms that contribute to DNA recognition by proteins.

Non-dissociative structural transitions of the Watson-Crick and reverse Watson-Crick А·Т DNA base pairs into the Hoogsteen and reverse Hoogsteen forms

In this study it was theoretically shown that discovered by us recently (Brovarets' et al., Frontiers in Chemistry, 2018, 6:8; doi: 10.3389/fchem.2018.00008) high-energetical, significantly non-planar (symmetry C₁), short-lived wobbled conformers of the classical Watson-Crick А·Т(WC), reverse Watson-Crick А·Т(rWC), Hoogsteen А·Т(Н) and reverse Hoogsteen А·Т(rН) DNA base pairs are the intermediates of their pairwise А∙Т(WC)/А∙Т(rWC) ↔ А∙Т(H)/А∙Т(rH) conformational transformations. These transitions do not require for their realization the energy-consumable anisotropic rotation of the amino group of A around the exocyclic C6-N6 bond. They are controlled by the non-planar transition states with quasi-orthogonal geometry (symmetry C₁) joined by the single intermolecular (Т)N3H···N6(А) H-bond (~4 kcal∙mol^-1). The Gibbs free energies of activation for these non-dissociative, dipole-active conformational transitions consist 7.33 and 7.81 kcal∙mol^-1, accordingly. Quantum-mechanical (QM) calculations in combination with Bader's quantum theory of "Atoms in Molecules" (QTAIM) have been performed at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory in the continuum with ε = 4 under normal conditions.