Characterizing the protonation state of cytosine in transient G·C Hoogsteen base pairs in duplex DNA

The bisulfite reaction with cytosine and genomic DNA structure

The bisulfite reaction with native DNA has been extensively employed in the detection of non-B DNA structures that can form spontaneously in DNA. These sequences are dynamic in that they can adopt both normal Watson-Crick paired B-DNA or unusual structures like the Triplex, G-Quadruplex, i-motif and Cruciform or Hairpin. Considerable evidence now suggests that these dynamic sequences play roles in both epigenetics and mutagenesis. The bisulfite reaction with native DNA offers a key approach to their detection. In this application whole cells, isolated nuclei or isolated DNA are treated with bisulfite under non-denaturing conditions in order to detect bisulfite accessible regions DNA that are associated with these structures. Here I review the stereochemistry of the bisulfite reaction, the electronic structure of its DNA cytosine substrates and its application in the detection of unusual structures in native DNA.

Harnessing Cytosine for Tunable Nanoparticle Self-Assembly Behavior Using Orthogonal Stimuli

Nucleobases control the assembly of DNA, RNA, etc. due to hydrogen bond complementarity. By combining these unique molecules with state-of-the-art synthetic polymers, it is possible to form nanoparticles whose self-assembly behavior could be altered under orthogonal stimuli (pH and temperature). Herein, we report the synthesis of cytosine-containing nanoparticles via aqueous reversible addition-fragmentation chain transfer polymerization-induced self-assembly. A poly(N-acryloylmorpholine) macromolecular chain transfer agent (mCTA) was chain-extended with cytosine acrylamide, and a morphological phase diagram was constructed. By exploiting the ability of cytosine to form dimers via hydrogen bonding, the self-assembly behavior of cytosine-containing polymers was altered when performed under acidic conditions. Under these conditions, stable nanoparticles could be formed at longer polymer chain lengths. Furthermore, the resulting nanoparticles displayed different morphologies compared to those at pH 7. Additionally, particle stability post-assembly could be controlled by varying pH and temperature. Finally, small-angle X-ray scattering was performed to probe their dynamic behavior under thermal cycling.

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.

Enhancement of intrinsic guanine fluorescence by protonation in DNA of various structures

Understanding the diversity of DNA structure and functions in biology requires tools to study this biomolecule selectively and thoroughly. Fluorescence methods are powerful technique for non-invasive research. Due to the low quantum yield, the intrinsic fluorescence of nucleotides has not been considered for use in the detection and differentiation of nucleic acid bases. Here, we have studied the influence of protonation of nucleotides on their fluorescence properties. We show that protonation of ATP and GTP leads to enhanced intrinsic fluorescence. Fluorescence enhancement at acidic pH has been observed for double-stranded DNA and single-stranded oligonucleotides. The formation of G4 secondary structures apparently protected certain nucleotides from protonation, resulting in less pronounced fluorescence enhancement. Furthermore, acid-induced depurination under protonation was less noticeable in G4 structures than in double-stranded and single-stranded DNA. We show that changes in the intrinsic fluorescence of guanine can be used as a sensitive sensor for changes in the structure of the DNA and for the protonation of specific nucleotides.

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.

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

The photophysics of protonated cytidine and hemiprotonated cytidine base pair: A computational study

We here study the effect that a lowering of the pH has on the excited state processes of cytidine and a cytidine/cytidine pair in solution, by integrating time-dependent density functional theory and CASSCF/CASPT2 calculations, and including solvent by a mixed discrete/continuum model. Our calculations reproduce the effect of protonation at N3 on the steady-state infrared and absorption spectra of a protonated cytidine (CH+ ), and predict that an easily accessible non-radiative deactivation route exists for the spectroscopic state, explaining its sub-ps lifetime. Indeed, an extremely small energy barrier separates the minimum of the lowest energy bright state from a crossing region with the ground electronic state, reached by out-of-plane motion of the hydrogen substituents of the CC double bond, the so-called ethylenic conical intersection typical of cytidine and other pyrimidine bases. This deactivation route is operative for the two bases forming an hemiprotonated cytidine base pair, [CH·C]+ , the building blocks of I-motif secondary structures, whereas interbase processes play a minor role. N3 protonation disfavors instead the nπ* transitions, associated with the long-living components of cytidine photoactivated dynamics.

Ag+-stabilized DNA triplex coupled with catalytic hairpin assembly and CRISPR/Cas12a amplifications for sensitive metallothionein assay

Metallothionein (MT) is a protein biomarker secreted by liver in response to the treatment for heavy metal toxicity and oncological diseases. On the basis of a new Ag+-stabilized DNA triplex probe (Ag+-SDTP), we establish a fluorescent biosensing system for high sensitivity detection of MT by combining catalytic hairpin assembly (CHA) and the CRISPR/Cas12a signal enhancements. The MT analyte complexes with Ag+ in Ag+-SDTP to disrupt the triplex structure and to release the ssDNA strands, which trigger subsequent CHA formation of many protospacer adjacent motif (PAM)-containing dsDNAs from two hairpins. Cas12a/crRNA further recognizes these PAM sequences to activate its trans-catalytic activity to cyclically cleave the fluorescently quenched ssDNA reporters to recovery drastically amplified fluorescence for detecting MT down to 0.34 nM within the dynamic range of 1∼800 nM. Moreover, the sensing method is able to selectively discriminate MT from other non-specific molecules and can realize low level detection of MT in diluted human serums, manifesting its potentiality for monitoring the disease-specific MT biomarker at trace levels.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

Measuring Local Electrostatic Potentials Around Nucleic Acids by Paramagnetic NMR Spectroscopy

Electrostatic potentials around macromolecules in the presence of mobile charges are difficult to assess especially for highly charged systems. Here, we report measurements of local electrostatic potentials around DNA by paramagnetic NMR. Through quantitative analysis of NMR paramagnetic relaxation enhancement arising from positively charged or neutral paramagnetic cosolutes, we were able to determine local electrostatic potentials around 1H nuclei at >100 sites in major and minor grooves of 13C,15N-labeled 15-bp DNA at 100 mM NaCl. Our experimental electrostatic potential data directly confirmed the Coulombic end effects of DNA. The effective near-surface electrostatic potentials from the NMR data were in good agreement with the theoretical predictions with the Poisson-Boltzmann equation. This NMR method allows for unprecedented experimental investigations into the electrostatic properties of nucleic acids.

Bifurcation of signalling in human innate immune pathways to NF-kB and IRF family activation

The human innate immune response can be activated through a wide range of stimuli. This multi-faceted system can be triggered by a range of immunostimulants including pathogen-associated molecular patterns (PAMPs) and damage-associated molecular patterns (DAMPs). These stimuli drive intracellular signalling pathways that branch off downstream to activate several distinct transcription factors. The two most impactful of which in innate immune outcomes are the NF-κB and the IRF family members. Both transcription factor families play defining roles in driving inflammation as well as the antiviral response. Pathways leading to their simultaneous activation share common upstream components but eventually distinct regulators which directly facilitate their activation. This review will discuss the current state of knowledge about what is known about how these pathways bifurcate to activate NF-κB and IRF family members.

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.

Slow Rotational Dynamics of Cytosine NH2 Groups in Double-Stranded DNA

Aromatic NH₂ groups are essential as hydrogen-bond donors in secondary structures of DNA and RNA. Although rapid rotations of NH₂ groups of adenine and guanine bases were previously characterized, there has been a lack of quantitative information about slow rotations of cytosine NH₂ groups in Watson-Crick base pairs. In this study, using an NMR method we had recently developed, we determined the kinetic rate constants and energy barriers for cytosine NH₂ rotations in a 15-base-pair DNA duplex. Our data show that the rotational dynamics of cytosine NH₂ groups depend on local environments. Qualitative correlation between the ranges of ¹⁵N chemical shifts and rotational time scales for various NH₂ groups of nucleic acids and proteins illuminates a relationship between the partial double-bond character of the C-N bond and the time scale for NH₂ rotations.

Fluorescence Sensing of the Panhandle Structure of the Influenza A Virus RNA Promoter by Thiazole Orange Base Surrogate-Carrying Peptide Nucleic Acid Conjugated with Small Molecule

We have developed a new class of triplex-forming peptide nucleic acid (PNA)-based fluorogenic probes for sensing of the panhandle structure of the influenza A virus (IAV) RNA promoter region. Here, a small molecule (DPQ) capable of selectively binding to the internal loop structure was conjugated with triplex-forming forced intercalation of the thiazole orange (tFIT) probe with natural PNA nucleobases. The resulting conjugate, tFIT-DPQ, showed a significant light-up response (83-fold) upon strong (K_d = 107 nM) and structure-selective binding to the IAV RNA promoter region under physiological conditions (pH 7.0, 100 mM NaCl). We demonstrated the conjugation of these two units through the suitable spacer was key to show useful binding and fluorogenic signaling functions. tFIT-DPQ facilitated the sensitive and selective detection of IAV RNA based on its binding to the promoter region. Furthermore, we found that tFIT-DPQ could work as a sensitive indicator for screening of test compounds targeting the IAV RNA promoter region in the fluorescence indicator displacement assay.

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.

Changes of C≡C Triple Bond Vibration that Disclosed Non-Canonical Cytosine Protonation in i-Motif-Forming Oligodeoxynucleotides

Non-canonical protonation at cytosine (C) in DNA is related to a formation of second order DNA structures such as i-motif, which has a role in gene regulation. Although the detailed structural information is indispensable for comprehension of their functions in cells, the protonation status of C in complicated environments is still elusive. To provide a reporter system of non-canonical protonation, we focused on the molecular vibration that could be monitored using the Raman spectroscopy. We prepared a cytosine derivative (^PC) with an acetylene unit as a Raman tag, and found that the Raman signal of acetylene in ^PC in oligodeoxynucleotides (ODNs) changed due to protonation at the cytosine ring which shortened an acetylene bond. The signal change in i-motif-forming ODNs was also observed in crowded environments with polyethylene glycol, evidencing protonation in i-motif DNA in complicated environments. This system would be one of tracking tools for protonation in DNA structures.

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 13 C or 15 N isotopically enriched DNA samples for NMR relaxation dispersion (RD) experiments. Here, using SELective Optimized Proton Experiments (SELOPE) 1 H CEST experiments employing high-power radiofrequency fields (B 1  >  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 1 H CEST were in very good agreement with counterparts measured using off-resonance 13 C /  15 N spin relaxation in the rotating frame (R 1 ρ ). It is shown that 1 H-1 H NOE effects which typically introduce artifacts in 1 H-based measurements of chemical exchange can be effectively suppressed by selective excitation, provided that the relaxation delay is short (≤  100 ms). The 1 H CEST experiment can be performed with ∼  10× higher throughput and ∼  100× lower cost relative to 13 C /  15 N R 1 ρ and enabled Hoogsteen chemical exchange measurements undetectable by R 1 ρ . The results reveal an increased propensity to form Hoogsteen bp's near terminal ends and a diminished propensity within A-tract motifs. The 1 H CEST experiment provides a basis for rapidly screening Hoogsteen breathing in duplex DNA, enabling identification of unusual motifs for more in-depth characterization.

Spectroscopic, thermodynamic and kinetic analysis of selective triplex formation by peptide nucleic acid with double-stranded RNA over its DNA counterpart

Unlike conventional triplex-forming oligonucleotide (TFO), triplex-forming peptide nucleic acid (PNA) can tightly bind with double-stranded RNA (dsRNA) than double-stranded DNA (dsDNA). Here, we performed spectroscopic, thermodynamic and kinetic experiments for triplex formation by PNA to examine different binding behaviors between PNA - dsRNA and PNA - dsDNA triplexes. We found 9-mer PNA (cytosine content of 66%) formed the thermally stable triplex with dsRNA compared to dsDNA over a wide range of pH (5.5-8.0), salt concentration (50-500 mM NaCl). Both the calorimetric binding constant and the association rate constant for dsRNA were larger than those for dsDNA, indicating the favorable association process for the PNA - dsRNA triplex formation. Comparison with the DNA/RNA heteroduplexes revealed that the DNA strand was detrimental to the triplex stability for PNA, a contrasting result for conventional TFO. The keys underlying the difference in the triplex formation of PNA with different duplexes appear to be the conformational adoptability and the geometric compatibility of PNA to fit the deep, narrow major groove of dsRNA and the helical rigidity difference of the duplexes. Our results emphasize the importance of both the sugar puckering of the duplex and the appropriate conformational flexibility of PNA for the triplex formation.

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.

The Molecular Tête-à-Tête between G-Quadruplexes and the i-motif in the Human Genome

G-Quadruplex (GQ) and i-motif structures are the paradigmatic examples of nonclassical tetrastranded nucleic acids having multifarious biological functions and widespread applications in therapeutics and material science. Recently, tetraplexes emerged as promising anticancer targets due to their structural robustness, gene-regulatory roles, and predominant distribution at specific loci of oncogenes. However, it is arguable whether the i-motif evolves in the complementary single-stranded region after GQ formation in its opposite strand and vice versa. In this review, we address the prerequisites and significance of the simultaneous and/or mutually exclusive formation of GQ and i-motif structures at complementary and sequential positions in duplexes in the cellular milieu. We discussed how their dynamic interplay Sets up cellular homeostasis and exacerbates carcinogenesis. The review gives insights into the spatiotemporal formation of GQ and i-motifs that could be harnessed to design different types of reporter systems and diagnostic platforms for potential bioanalytical and therapeutic intervention.

Structural dynamics of double-stranded DNA with epigenome modification

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson-Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson-Crick base pair to allow base-flipping in living cells.

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.

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.

Microhydration of protonated biomolecular building blocks: protonated pyrimidine

Protonation and hydration of biomolecules govern their structure, conformation, and function. Herein, we explore the microhydration structure in mass-selected protonated pyrimidine-water clusters (H⁺Pym-W_n, n = 1-4) by a combination of infrared photodissociation spectroscopy (IRPD) between 2450 and 3900 cm^-1 and density functional theory (DFT) calculations at the dispersion-corrected B3LYP-D3/aug-cc-pVTZ level. We further present the IR spectrum of H⁺Pym-N₂ to evaluate the effect of solvent polarity on the intrinsic molecular parameters of H⁺Pym. Our combined spectroscopic and computational approach unequivocally shows that protonation of Pym occurs at one of the two equivalent basic ring N atoms and that the ligands in H⁺Pym-L (L = N₂ or W) preferentially form linear H-bonds to the resulting acidic NH group. Successive addition of water ligands results in the formation of a H-bonded solvent network which increasingly weakens the NH group. Despite substantial activation of the N-H bond upon microhydration, no intracluster proton transfer occurs up to n = 4 because of the balance of relative proton affinities of Pym and W_n and the involved solvation energies. Comparison to neutral Pym-W_n clusters reveals the drastic effects of protonation on microhydration with respect to both structure and interaction strength.

5-Carboxylcytosine and Cytosine Protonation Distinctly Alter the Stability and Dehybridization Dynamics of the DNA Duplex

Applications associated with nucleobase protonation events are grounded in their fundamental impact on DNA thermodynamics, structure, and hybridization dynamics. Of the canonical nucleobases, N3 protonation of cytosine (C) is the most widely utilized in both biology and nanotechnology. Naturally occurring C derivatives that shift the N3 pK_a introduce an additional level of tunability. The epigenetic nucleobase 5-carboxylcytosine (caC) presents a particularly interesting example since this derivative forms Watson-Crick base pairs of similar stability and displays pH-dependent behavior over the same range as the canonical nucleobase. However, the titratable group in caC corresponds to the exocyclic carboxyl group rather than N3, and the implications of these divergent protonation events toward DNA hybridization thermodynamics, kinetics, and base pairing dynamics remain poorly understood. Here, we study the pH dependence of these physical properties using model oligonucleotides containing C and caC with FTIR and temperature-jump IR spectroscopy. We demonstrate that N3 protonation of C completely disrupts duplex stability, leading to large shifts in the duplex/single-strand equilibrium, a reduction in the cooperativity of melting, and an acceleration in the rate of duplex dissociation. In contrast, while increasing 5-carboxyl protonation in caC-containing duplexes induces an increase in base pair fluctuations, the DNA duplex can tolerate substantial protonation without significant perturbation to the duplex/single-strand equilibrium. However, 5-carboxyl protonation has a large impact on hybridization kinetics by reducing the transition state free energy. Our thermodynamic and kinetic analysis provides new insight on the impact of two divergent protonation mechanisms in naturally occurring nucleobases on the biophysical properties of DNA.

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.

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 micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

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.

Protonation of Nucleobases in Single- and Double-Stranded DNA

Single-stranded model oligodeoxyribonucleotides, each containing a single protonatable base-cytosine, adenine, guanine, or 5-methylcytosine-centrally located in a background of non-protonatable thymine residues, were acid-titrated in aqueous solution, with UV monitoring. The basicity of the central base was shown to depend on the type of the central base and its nearest neighbours and to rise with increasing oligonucleotide length and decreasing ionic strength of the solution. More complex model oligonucleotides, each containing a centrally located 5-methylcytosine base, were comparatively evaluated in single-stranded and double-stranded form, by UV spectroscopy and high-field NMR. The N³ protonation of the 5-methylcytosine moiety in the double-stranded case occurred at much lower pH, at which the duplex was already experiencing general dissociation, than in the single-stranded case. The central guanine:5-methylcytosine base pair remained intact up to this point, possibly due to an unusual alternative protonation on O² of the 5-methylcytosine moiety, already taking place at neutral or weakly basic pH, as indicated by UV spectroscopy, thus suggesting that 5-methylcytosine sites in double-stranded DNA might be protonated to a significant extent under physiological conditions.

Apurinic/apyrimidinic (AP) endonuclease 1 processing of AP sites with 5' mismatches

Despite the DNA duplex being central to biological functions, many intricacies of this molecule, including the dynamic nature of mismatched base pairing, are still unknown. The unique conformations adopted by DNA mismatches can provide insight into the forces at play between nucleotides. Moreover, DNA-binding proteins apply their own individualized steric and electrochemical influences on the nucleotides that they interact with, further altering base-pairing conformations. Here, seven X-ray crystallographic structures of the human nuclease apurinic/apyrimidinic (AP) endonuclease 1 (APE1) in complex with its substrate target flanked by a 5' mismatch are reported. The structures reveal how APE1 influences the conformations of a variety of different mismatched base pairs. Purine-purine mismatches containing a guanine are stabilized by a rotation of the guanine residue about the N-glycosidic bond to utilize the Hoogsteen edge for hydrogen bonding. Interestingly, no rotation of adenine, the other purine, is observed. Mismatches involving both purine and pyrimidine bases adopt wobble conformations to accommodate the mismatch. Pyrimidine-pyrimidine mismatches also wobble; however, the smaller profile of a pyrimidine base results in a gap between the Watson-Crick faces that is reduced by a C1'-C1' compression. These results advance our understanding of mismatched base pairing and the influence of a bound protein.

(C2G4)n repeat expansion sequences from the C9orf72 gene form an unusual DNA higher-order structure in the pH range of 5-6

Massive expansion of a DNA hexanucleotide sequence repeat (C₂G₄) within the human C9orf72 gene has been linked to a number of neurodegenerative diseases. In sodium or potassium salt solutions, single-stranded d(C₂G₄)_n DNAs fold to form G-quadruplexes. We have found that in magnesium or lithium salt solutions, especially under slightly acidic conditions, d(C₂G₄)_n oligonucleotides fold to form a distinctive higher order structure whose most striking feature is an “inverted” circular dichroism spectrum, which is distinguishable from the spectrum of the left handed DNA double-helix, Z-DNA. On the basis of CD spectroscopy, gel mobility as well as chemical protection analysis, we propose that this structure, which we call “iCD-DNA”, may be a left-handed Hoogsteen base-paired duplex, an unorthodox G-quadruplex/i-motif composite, or a non-canonical, “braided” DNA triplex. Given that iCD-DNA forms under slightly acidic solution conditions, we do not know at this point in time whether or not it forms within living cells.

Modulation of Hoogsteen dynamics on DNA recognition

In naked duplex DNA, G–C and A–T Watson-Crick base pairs exist in dynamic equilibrium with their Hoogsteen counterparts. Here, we used nuclear magnetic resonance (NMR) relaxation dispersion and molecular dynamics (MD) simulations to examine how Watson-Crick/Hoogsteen dynamics are modulated upon recognition of duplex DNA by the bisintercalator echinomycin and monointercalator actinomycin D. In both cases, DNA recognition results in the quenching of Hoogsteen dynamics at base pairs involved in intermolecular base-specific hydrogen bonds. In the case of echinomycin, the Hoogsteen population increased 10-fold for base pairs flanking the chromophore most likely due to intermolecular stacking interactions, whereas actinomycin D minimally affected Hoogsteen dynamics at other sites. Modulation of Hoogsteen dynamics at binding interfaces may be a general phenomenon with important implications for DNA–ligand and DNA–protein recognition.

DNA is found in a dynamic equilibrium between standard Watson-Crick (WC) base pairs and non-standard Hoogsteen (HG) base pairs. Here the authors describe the influence of echinomycin and actinomycin D ligands binding on the HG-WC base pair dynamics in DNA.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Detection of protonated non-Watson-Crick base pairs using electrospray ionization mass spectrometry

Many studies have shown that protonated nucleic acid base pairs are involved in a wide variety of nucleic acid structures. However, little information is available on relative stability of hemiprotonated self- and non-self-dimers at monomer level. We used electrospray ionization mass spectrometry (ESI-MS) to evaluate the relative stability under various concentrations of hydrogen ion. These enable conjecture of the formation of protonated non-Watson-Crick base pairs based on DNA and RNA base sequence. In the present study, we observed that ESI-MS peaks corresponded to respective self-dimers for all examined nucleosides except for adenosine. Peak heights depended on the concentration of hydrogen ion. The ESI-MS peak heights of the hemiprotonated cytidine dimers and the hemiprotonated thymidine dimer sharply increased with increased concentration of hydrogen ion, suggesting direct participation of hydrogen ion in dimer formations. In ESI-MS measurements of the solutions containing adenosine, cytidine, thymidine and guanosine, we observed protonated cytidine-guanosine dimer (CH+-G) and protonated cytidine-thymidine dimer (CH+-T) in addition to hemiprotonated cytidine-cytidine dimer (CH+-C) with following relative peak height, (CH+-C) > (CH+-G) ≈ (CH+-T) > (CH+-A). Additionally, in the ESI-MS measurements of solutions containing adenosine, thymidine and guanosine, we observed a considerable amount of protonated adenosine-guanosine (AH+-G) and protonated adenosine-thymidine (AH+-T).

Surprising Conformers of the Biologically Important A·T DNA Base Pairs: QM/QTAIM Proofs

For the first time novel high-energy conformers-A·T(wWC) (5.36), A·T(wrWC) (5.97), A·T(wH) (5.78), and A·T(wrH) (ΔG = 5.82 kcal·mol-1) (See Graphical Abstract) were revealed for each of the four biologically important A·T DNA base pairs - Watson-Crick A·T(WC), reverse Watson-Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of quantum-mechanical theory in the continuum with ε = 4 under normal conditions. Each of these conformers possesses substantially non-planar wobble (w) structure and is stabilized by the participation of the two anti-parallel N6H/N6H'…O4/O2 and N3H…N6 H-bonds, involving the pyramidalized amino group of the A DNA base as an acceptor and a donor of the H-bonding. The transition states - TSA·T(WC)↔A·T(wWC), TSA·T(rWC)↔A·T(wrWC), TSA·T(H)↔A·T(wH), and TSA·T(rH)↔A·T(wrH), controlling the dipole-active transformations of the conformers from the main plane-symmetric state into the high-energy, significantly non-planar state and vice versa, were localized. They also possess wobble structures similarly to the high-energy conformers and are stabilized by the participation of the N6H/N6H'…O4/O2 and N3H…N6 H-bonds. Discovered conformers of the A·T DNA base pairs are dynamically stable short-lived structures [lifetime τ = (1.4-3.9) ps]. Their possible biological significance and future perspectives have been briefly discussed.

Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA

The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.

Mechanism Underlying the Nucleobase-Distinguishing Ability of Benzopyridopyrimidine (BPP)

Benzopyridopyrimidine (BPP) is a fluorescent nucleobase analogue capable of forming base pairs with adenine (A) and guanine (G) at different sites. When incorporated into oligodeoxynucleotides, it is capable of differentiating between the two purine nucleobases by virtue of the fact that its fluorescence is largely quenched when it is base-paired to guanine, whereas base-pairing to adenine causes only a slight reduction of the fluorescence quantum yield. In the present article, the photophysics of BPP is investigated through computer simulations. BPP is found to be a good charge acceptor, as demonstrated by its positive and appreciably large electron affinity. The selective quenching process is attributed to charge transfer (CT) from the purine nucleobase, which is predicted to be efficient in the BPP-G base pair, but essentially inoperative in the BPP-A base pair. The CT process owes its high selectivity to a combination of two factors: the ionization potential of guanine is lower than that of adenine, and less obviously, the site occupied by guanine enables a greater stabilization of the CT state through electrostatic interactions than the one occupied by adenine. The case of BPP illustrates that molecular recognition via hydrogen bonding can enhance the selectivity of photoinduced CT processes.

An Extraordinarily Stable DNA Minidumbbell

The minidumbbell (MDB) is a new type of native DNA structure. At neutral pH, two TTTA or CCTG repeats can fold into the highly compact MDB with a melting temperature of ∼22 °C. Owing to the relatively low thermodynamic stability, MDBs have been proposed to be the structural intermediates that lead to efficient DNA repair escape and thus repeat expansions. In this study, we reveal that two CCTG repeats can also form an extraordinarily stable MDB with a melting temperature of ∼46 °C at pH 5.0. This unusual stability predominantly results from the formation of a three hydrogen bond C⁺·C mispair between the two minor groove cytosine residues. Due to the drastic stability change, the CCTG MDB, when combined with its complementary sequence, shows instant and complete structural conversions when the pH switches between 5.0 and 7.0, making the system serve as a simple and efficient pH-controlled molecular switch.

Robust IR-based detection of stable and fractionally populated G-C+ and A-T Hoogsteen base pairs in duplex DNA

Noncanonical G-C+ and A-T Hoogsteen base pairs can form in duplex DNA and play roles in recognition, damage repair, and replication. Identifying Hoogsteen base pairs in DNA duplexes remains challenging due to difficulties in resolving syn versus antipurine bases with X-ray crystallography; and size limitations and line broadening can make them difficult to characterize by NMR spectroscopy. Here, we show how infrared (IR) spectroscopy can identify G-C+ and A-T Hoogsteen base pairs in duplex DNA across a range of different structural contexts. The utility of IR-based detection of Hoogsteen base pairs is demonstrated by characterizing the first example of adjacent A-T and G-C+ Hoogsteen base pairs in a DNA duplex where severe broadening complicates detection with NMR.

Cooperative Interactions in the Hammerhead Ribozyme Drive pKa Shifting of G12 and Its Stacked Base C17

General acid-base catalysis is a key mechanistic strategy in protein and RNA enzymes. Ribozymes use hydrated metal ions, nucleobases, and organic cofactors to carry this out. In most small ribozymes, a guanosine is positioned to participate in proton transfer with the nucleophilic 2'-OH. The unshifted pK_a values for nucleobases and solvated metal ions are far from neutrality, however, and thus nonideal for general acid-base catalysis. Herein, evidence is provided for cooperative interaction in the hammerhead ribozyme among the guanine that interacts with the nucleophilic 2'-OH, G12, the -1 nucleobase C17, and Mg²⁺ ions. We introduce global fitting for analyzing ribozyme rate-pH data parametric in Mg²⁺ concentration and benchmark this method on data from the hepatitis delta virus ribozyme. We then apply global fitting to new rate-pH data for the hammerhead ribozyme using a minimal three-dimensional, four-channel cooperative model. The value for the pK_a of G12 that we obtain is channel-dependent and varies from 8.1 to 9.9, shifting closest toward neutrality in the presence of two cationic species: C17H⁺ and a Mg²⁺ ion. The value for the pK_a of the -1 nucleotide, C17, is increased a remarkable 3.5-5 pK_a units toward neutrality. Shifting of the pK_a of C17 appears to be driven by an electrostatic sandwich of C17 between carbonyl groups of the 5'-neighboring U and of G12 and involves cation-π interactions. Rate-pH profiles reveal that the major reactive channel under biological Mg²⁺ and pH involves a cationic C17 rather than a second metal ion. Substitution of a cationic base for a metal underscores the versatility of RNA.

Constant-pH Hybrid Nonequilibrium Molecular Dynamics-Monte Carlo Simulation Method

A computational method is developed to carry out explicit solvent simulations of complex molecular systems under conditions of constant pH. In constant-pH simulations, preidentified ionizable sites are allowed to spontaneously protonate and deprotonate as a function of time in response to the environment and the imposed pH. The method, based on a hybrid scheme originally proposed by H. A. Stern (J. Chem. Phys. 2007, 126, 164112), consists of carrying out short nonequilibrium molecular dynamics (neMD) switching trajectories to generate physically plausible configurations with changed protonation states that are subsequently accepted or rejected according to a Metropolis Monte Carlo (MC) criterion. To ensure microscopic detailed balance arising from such nonequilibrium switches, the atomic momenta are altered according to the symmetric two-ends momentum reversal prescription. To achieve higher efficiency, the original neMD-MC scheme is separated into two steps, reducing the need for generating a large number of unproductive and costly nonequilibrium trajectories. In the first step, the protonation state of a site is randomly attributed via a Metropolis MC process on the basis of an intrinsic pKa; an attempted nonequilibrium switch is generated only if this change in protonation state is accepted. This hybrid two-step inherent pKa neMD-MC simulation method is tested with single amino acids in solution (Asp, Glu, and His) and then applied to turkey ovomucoid third domain and hen egg-white lysozyme. Because of the simple linear increase in the computational cost relative to the number of titratable sites, the present method is naturally able to treat extremely large systems.