Studies of base pair kinetics by NMR measurement of proton exchange

Optimising in-cell NMR acquisition for nucleic acids

Understanding the structure and function of nucleic acids in their native environment is crucial to structural biology and one focus of in-cell NMR spectroscopy. Many challenges hamper in-cell NMR in human cell lines, e.g. sample decay through cell death and RNA degradation. The resulting low signal intensities and broad line widths limit the use of more complex NMR experiments, reducing the possible structural and dynamic information that can be extracted. Here, we optimize the detection of imino proton signals, indicators of base-pairing and therefore secondary structure, of a double-stranded DNA oligonucleotide in HeLa cells, using selective excitation. We demonstrate the reproducible quantification of in-cell selective longitudinal relaxation times (selT1), which are reduced compared to the in vitro environment, as a result of interactions with the complex cellular environment. By measuring the intracellular selT1, we optimize the existing proton pulse sequences, and shorten measurement time whilst enhancing the signal gained per unit of time. This exemplifies an advantage of selective excitation over conventional methods like jump-return water suppression for in-cell NMR. Furthermore, important experimental controls are discussed, including intracellular quantification, supernatant control measurements, as well as the processing of lowly concentrated in-cell NMR samples. We expect that robust and fast in-cell NMR experiments of nucleic acids will facilitate the study of structure and dynamics and reveal their functional correlation.

Dynamics of 8-Oxoguanine in DNA: Decisive Effects of Base Pairing and Nucleotide Context

8-Oxo-7,8-dihydroguanine (oxoG), an abundant DNA lesion, can mispair with adenine and induce mutations. To prevent this, cells possess DNA repair glycosylases that excise either oxoG from oxoG:C pairs (bacterial Fpg, human OGG1) or A from oxoG:A mispairs (bacterial MutY, human MUTYH). Early lesion recognition steps remain murky and may include enforced base pair opening or capture of a spontaneously opened pair. We adapted the CLEANEX-PM NMR protocol to detect DNA imino proton exchange and analyzed the dynamics of oxoG:C, oxoG:A, and their undamaged counterparts in nucleotide contexts with different stacking energy. Even in a poorly stacking context, the oxoG:C pair did not open easier than G:C, arguing against extrahelical base capture by Fpg/OGG1. On the contrary, oxoG opposite A significantly populated the extrahelical state, which may assist recognition by MutY/MUTYH.

Shedding light on the base-pair opening dynamics of nucleic acids in living human cells

Base-pair opening is a fundamental property of nucleic acids that plays important roles in biological functions. However, studying the base-pair opening dynamics inside living cells has remained challenging. Here, to determine the base-pair opening kinetics inside living human cells, the exchange rate constant ([Formula: see text]) of the imino proton with the proton of solvent water involved in hairpin and G-quadruplex (GQ) structures is determined by the in-cell NMR technique. It is deduced on determination of [Formula: see text] values that at least some G-C base pairs of the hairpin structure and all G-G base-pairs of the GQ structure open more frequently in living human cells than in vitro. It is suggested that interactions with endogenous proteins could be responsible for the increase in frequency of base-pair opening. Our studies demonstrate a difference in dynamics of nucleic acids between in-cell and in vitro conditions.

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.

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.

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.

Impact of 3-deazapurine nucleobases on RNA properties

Deazapurine nucleosides such as 3-deazaadenosine (c3A) are crucial for atomic mutagenesis studies of functional RNAs. They were the key for our current mechanistic understanding of ribosomal peptide bond formation and of phosphodiester cleavage in recently discovered small ribozymes, such as twister and pistol RNAs. Here, we present a comprehensive study on the impact of c3A and the thus far underinvestigated 3-deazaguanosine (c3G) on RNA properties. We found that these nucleosides can decrease thermodynamic stability of base pairing to a significant extent. The effects are much more pronounced for 3-deazapurine nucleosides compared to their constitutional isomers of 7-deazapurine nucleosides (c7G, c7A). We furthermore investigated base pair opening dynamics by solution NMR spectroscopy and revealed significantly enhanced imino proton exchange rates. Additionally, we solved the X-ray structure of a c3A-modified RNA and visualized the hydration pattern of the minor groove. Importantly, the characteristic water molecule that is hydrogen-bonded to the purine N3 atom and always observed in a natural double helix is lacking in the 3-deazapurine-modified counterpart. Both, the findings by NMR and X-ray crystallographic methods hence provide a rationale for the reduced pairing strength. Taken together, our comparative study is a first major step towards a comprehensive understanding of this important class of nucleoside modifications.

Magnetization Transfer to Enhance NOE Cross-Peaks among Labile Protons: Applications to Imino-Imino Sequential Walks in SARS-CoV-2-Derived RNAs

2D NOESY plays a central role in structural NMR spectroscopy. We have recently discussed methods that rely on solvent‐driven exchanges to enhance NOE correlations between exchangeable and non‐exchangeable protons in nucleic acids. Such methods, however, fail when trying to establish connectivities within pools of labile protons. This study introduces an alternative that also enhances NOEs between such labile sites, based on encoding a priori selected peaks by selective saturations. The resulting selective magnetization transfer (SMT) experiment proves particularly useful for enhancing the imino–imino cross‐peaks in RNAs, which is a first step in the NMR resolution of these structures. The origins of these enhancements are discussed, and their potential is demonstrated on RNA fragments derived from the genome of SARS‐CoV‐2, recorded with better sensitivity and an order of magnitude faster than conventional 2D counterparts.

Magnetization Transfer to Enhance NOE Cross-Peaks among Labile Protons: Applications to Imino-Imino Sequential Walks in SARS-CoV-2-Derived RNAs

2D NOESY plays a central role in structural NMR spectroscopy. We have recently discussed methods that rely on solvent-driven exchanges to enhance NOE correlations between exchangeable and non-exchangeable protons in nucleic acids. Such methods, however, fail when trying to establish connectivities within pools of labile protons. This study introduces an alternative that also enhances NOEs between such labile sites, based on encoding a priori selected peaks by selective saturations. The resulting selective magnetization transfer (SMT) experiment proves particularly useful for enhancing the imino-imino cross-peaks in RNAs, which is a first step in the NMR resolution of these structures. The origins of these enhancements are discussed, and their potential is demonstrated on RNA fragments derived from the genome of SARS-CoV-2, recorded with better sensitivity and an order of magnitude faster than conventional 2D counterparts.

Duplexes Formed by G4C2 Repeats Contain Alternate Slow- and Fast-Flipping G·G Base Pairs

Aberrant expansion of the hexanucleotide GGGGCC (or G4C2) repeat in the human C9ORF72 gene is the most common genetic factor found behind amyotrophic lateral sclerosis and frontotemporal dementia. The hypothesized pathways, through which the repeat expansions contribute to the pathology, involve one or more secondary structural forms of the DNA and/or RNA sequences, such as G-quadruplexes, duplexes, and hairpins. Here, we study the structures of DNA and RNA duplexes formed by G4C2 repeats, which contain G(syn)·G(anti) base pairs flanked by either G·C or C·G base pairs. We show that duplexes formed by G4C2 repeats contain alternately two types of G·G pair contexts exhibiting different syn-anti base flipping dynamics (∼100 ms vs ∼2 ms for DNA and ∼50 ms vs ∼20 ms for RNA at 10 °C, respectively) depending on the flanking bases, with the slow-flipping G·G pairs being flanked by a guanine at the 5'-end and the fast-flipping G·G pairs being flanked by a cytosine at the 5'-end. Our findings on the structures and dynamics of G·G base pairs in DNA and RNA duplexes formed by G4C2 repeats provide a foundation for further studies of the functions and targeting of such biologically relevant motifs.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

A domain-level DNA strand displacement reaction enumerator allowing arbitrary non-pseudoknotted secondary structures

Information technologies enable programmers and engineers to design and synthesize systems of startling complexity that nonetheless behave as intended. This mastery of complexity is made possible by a hierarchy of formal abstractions that span from high-level programming languages down to low-level implementation specifications, with rigorous connections between the levels. DNA nanotechnology presents us with a new molecular information technology whose potential has not yet been fully unlocked in this way. Developing an effective hierarchy of abstractions may be critical for increasing the complexity of programmable DNA systems. Here, we build on prior practice to provide a new formalization of 'domain-level' representations of DNA strand displacement systems that has a natural connection to nucleic acid biophysics while still being suitable for formal analysis. Enumeration of unimolecular and bimolecular reactions provides a semantics for programmable molecular interactions, with kinetics given by an approximate biophysical model. Reaction condensation provides a tractable simplification of the detailed reactions that respects overall kinetic properties. The applicability and accuracy of the model is evaluated across a wide range of engineered DNA strand displacement systems. Thus, our work can serve as an interface between lower-level DNA models that operate at the nucleotide sequence level, and high-level chemical reaction network models that operate at the level of interactions between abstract species.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Kinetic Mechanism of RNA Helix-Terminal Basepairing-A Kinetic Minima Network Analysis

RNA functions are often kinetically controlled. The folding kinetics of RNAs involves global structural changes and local nucleotide movement, such as base flipping. The most elementary step in RNA folding is the closing and opening of a basepair. By integrating molecular dynamics simulation, master equation, and kinetic Monte Carlo simulation, we investigate the kinetics mechanism of RNA helix-terminal basepairing. The study reveals a six-state folding scheme with three dominant folding pathways of tens, hundreds, and thousands of nanoseconds of folding timescales, respectively. The overall kinetics is rate limited by the detrapping of a misfolded state with the overall folding time of 10^-5 s. Moreover, the analysis examines the different roles of the various driving forces, such as the basepairing and stacking interactions and the ion binding/dissociation effects on structural changes. The results may provide useful insights for developing a basepair opening/closing rate model and further kinetics models of large RNAs.

Base-pair Opening Dynamics of Nucleic Acids in Relation to Their Biological Function

Base-pair opening is a conformational transition that is required for proper biological function of nucleic acids. Hydrogen exchange, observed by NMR spectroscopic experiments, is a widely used method to study the thermodynamics and kinetics of base-pair opening in nucleic acids. The hydrogen exchange data of imino protons are analyzed based on a two-state (open/closed) model for the base-pair, where hydrogen exchange only occurs from the open state. In this review, we discuss examples of how hydrogen exchange data provide insight into several interesting biological processes involving functional interactions of nucleic acids: i) selective recognition of DNA by proteins; ii) regulation of RNA cleavage by site-specific mutations; iii) intermolecular interaction of proteins with their target DNA or RNA; iv) formation of PNA:DNA hybrid duplexes.

•

This review systematically summarizes hydrogen exchange theory for base-paired imino protons of nucleic acids.

•

Base-pair opening kinetics explain how the DNA can be selectively recognized by its target proteins.

•

Base-pair opening kinetics explain the mechanisms by which site-specific mutations regulate RNA cleavage.

•

Hydrogen exchange studies can elucidate the intermolecular interaction of proteins with their target DNA or RNA.

Substituting Inosine for Guanosine in DNA: Structural and Dynamic Consequences

Inosine differs from the guanosine nucleoside only by the absence of the N2 amino group. Both nucleosides also have similar electrostatic potentials. Therefore, substituting I for G has been used to probe various properties of nucleic acids and to facilitate the interpretation of binding studies. In particular, the absence of the amino group permits the assessment of its importance in the binding of ligands to the minor groove of duplex DNA. It has been known for some time that an I-C base pair is of lower stability than a regular G-C base pair, which needs to be considered when making DNA constructs containing inosine. However, it is generally assumed that both base pairs are structurally highly similar. To test this assumption in an identical sequence environment, we have determined the fine structure of two hairpin DNA substrates that differ only in the substitution of an I-C base pair for a G-C base pair. The structures have been solved using nuclear magnetic resonance (NMR) restraints in conjunction with Mardigras and molecular dynamics. The structural data are complemented with thermodynamic and dynamic data to get a comprehensive evaluation of the consequences of G-C vs I-C base pair substitutions. Our data show a strong similarity in the structures of the hairpins, but a significant difference in the melting temperatures, T m. This difference is also reflected in the drastically decreased base pair lifetime of 7.4 milliseconds compared to the G-C base pair lifetime of 155 milliseconds. The substitution of I-C for G-C is to probe for specific effect due to the amino group is satisfactory, as long as the lowered thermal stability and the drastically increased local dynamics are considered.

Energetics of base flipping at a DNA mismatch site confined at the latch constriction of α-hemolysin

Unique, two-state modulating current signatures are observed when a cytosine-cytosine mismatch pair is confined at the 2.4 nm latch constriction of the α-hemolysin (αHL) nanopore. We have previously speculated that the modulation is due to base flipping at the mismatch site. Base flipping is a biologically significant mechanism in which a single base is rotated out of the DNA helical stack by 180°. It is the mechanism by which enzymes are able to access bases for repair operations without disturbing the global structure of the helix. Here, temperature dependent ion channel recordings of individual double-stranded DNA duplexes inside αHL are used to derive thermodynamic (ΔH, ΔS) and kinetic (E_A) parameters for base flipping of a cytosine at an unstable cytosine-cytosine mismatch site. The measured activation energy for flipping a cytosine located at the latch of αHL out of the helix (18 ± 1 kcal mol^-1) is comparable to that previously reported for base flipping at mismatch sites from NMR measurements and potential mean force calculations. We propose that the αHL nanopore is a useful tool for measuring conformational changes in dsDNA at the single molecule level.

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.

Enhanced basepair dynamics pre-disposes protein-assisted flips of key bases in DNA strand separation during transcription initiation

Localized separation of strands of duplex DNA is a necessary step in many DNA-dependent processes, including transcription and replication.

Dynamic and Energetic Signatures of Adenine Tracts in a rA-dT RNA-DNA Hybrid and in Homologous RNA-DNA, RNA-RNA, and DNA-DNA Double Helices

Nuclear magnetic resonance spectroscopy and proton exchange are being used to characterize the opening reactions of individual base pairs in the RNA-DNA hybrid 5'-rGCGAUAAAAAGGCC-3'/5'-dGGCCTTTTTATCGC-3'. The hybrid contains a central tract of five rA-dT base pairs. The rates and the equilibrium constant of the opening reaction for each base pair are determined from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results are compared to those previously obtained by our laboratory for three homologous duplexes of the same base sequence (except for the appropriate T/U substitution), containing tracts of dA-rU, rA-rU, or dA-dT base pairs. The rA-dT tract is distinguished by an enhanced propensity of the base pairs to exist in the extrahelical state. The opening rates of rA-dT base pairs also exhibit a strong dependence on the location of the base pair in the structure; namely, as one advances into the tract, the opening rates of rA-dT base pairs gradually decrease. The local stability of each rA-dT base pair within the tract is the same as that of the corresponding rA-rU base pair in the homologous RNA-only duplex but differs from the stabilities of dA-dT and dA-rU base pairs in the other two duplexes (namely, dA-dT > rA-dT > dA-rU). These results demonstrate that, in nucleic acid double helices with the same base sequence, the opening dynamics and the energetics of individual base pairs are strongly influenced by the nature of the strand and by the structural context of the base pair.

5-Formylcytosine does not change the global structure of DNA

The mechanism by which the recently identified DNA modification 5-formylcytosine (^fC) is recognized by epigenetic writer and reader proteins is not known. Recently, an unusual DNA structure, F-DNA, has been proposed as the basis for enzyme recognition of clusters of ^fC. We used NMR and X-ray crystallography to compare several modified DNA duplexes with unmodified analogs and found that in the crystal state the duplexes all belong to the A family, whereas in solution they are all members of the B family. We found that, contrary to previous findings, ^fC does not significantly affect the structure of DNA, although there are modest local differences at the modification sites. Hence, global conformation changes are unlikely to account for the recognition of this modified base, and our structural data favor a mechanism that operates at base-pair resolution for the recognition of ^fC by epigenome-modifying enzymes.

Comparison of the free and ligand-bound imino hydrogen exchange rates for the cocaine-binding aptamer

Using NMR magnetization transfer experiments, the hydrogen exchange rate constants (k _ex ) of the DNA imino protons in the cocaine-binding aptamer have been determined for the free, cocaine-bound, and quinine-bound states. The secondary structure of the cocaine-binding aptamer is composed of three stems built around a three-way junction. In the free aptamer the slowest exchanging imino protons are located in the middle of the stems. The highest k _ex values were found for a nucleotide in the GAA loop of stem 3 and for nucleotides at the end of the stems that form the three-way junction structure and in the tandem GA mismatch. Upon ligand binding, the k _ex values of nucleotides at the ligand binding site are reduced, indicating that these base pairs become more stable or less solvent accessible in the bound state. The imino proton k _ex values of nucleotides located away from the binding site are only minimally affected by ligand binding.

Sequence dependency of canonical base pair opening in the DNA double helix

The flipping-out of a DNA base from the double helical structure is a key step of many cellular processes, such as DNA replication, modification and repair. Base pair opening is the first step of base flipping and the exact mechanism is still not well understood. We investigate sequence effects on base pair opening using extensive classical molecular dynamics simulations targeting the opening of 11 different canonical base pairs in two DNA sequences. Two popular biomolecular force fields are applied. To enhance sampling and calculate free energies, we bias the simulation along a simple distance coordinate using a newly developed adaptive sampling algorithm. The simulation is guided back and forth along the coordinate, allowing for multiple opening pathways. We compare the calculated free energies with those from an NMR study and check assumptions of the model used for interpreting the NMR data. Our results further show that the neighboring sequence is an important factor for the opening free energy, but also indicates that other sequence effects may play a role. All base pairs are observed to have a propensity for opening toward the major groove. The preferred opening base is cytosine for GC base pairs, while for AT there is sequence dependent competition between the two bases. For AT opening, we identify two non-canonical base pair interactions contributing to a local minimum in the free energy profile. For both AT and CG we observe long-lived interactions with water and with sodium ions at specific sites on the open base pair.

The DNA double helix, a molecule that stores biological information, has become an iconic image of biomedical research. In order to use or repair the information it carries, the bases that are stacked in the helix need to be chemically exposed. This can happen either by separating the two strands in the helix or by flipping out individual bases. Here, we focus on the latter process. Usually proteins are involved in interactions with bases, but it is still unclear if bases are pulled out actively by proteins or if they act on spontaneously flipped bases. Although experiments can detect base pair opening, it is difficult to detect which base moves in which direction. Here, we present results from molecular dynamics simulations using a recently developed sampling method which improves the statistics in the simulations by enhancing the probability of the base pair opening event. We observe differences in probability, modes and mechanism of opening that depend not only on the types of the bases in the pair, but also strongly on their neighbors. This provides essential information for understanding how DNA functions.

Purity assessment of ginsenoside Rg1 using quantitative 1H nuclear magnetic resonance

Ginseng herbs comprise a group of the most popular herbs, including Panax ginseng, P. notoginseng and P. quinquefolius (Family Araliaceae), which are used as traditional Chinese medicine (TCM) and are some of the best-selling natural products in the world. The accurate quantification of ginsenoside Rg1 is one of the major aspects of its quality control. However, the purity of the commercial Rg1 chemical reference substance (CRS) is often measured with high-performance chromatography coupled with an ultraviolet detector (HPLC-UV), which is a selective detector with unequal responses to different compounds; thus, this detector introduces probable error to purity assessments. In the present study, quantitative nuclear magnetic resonance (qNMR), due to its absolute quantification ability, was applied to accurately assess the purity of Rg1 CRS. Phenylmethyl phthalate was used as the internal standard (IS) to calibrate the purity of Rg1 CRS. The proton signal of Rg1 CRS in methanol-d4 at 4.37ppm was selected to avoid interfering signals, enabling accurate quantitative analysis. The relaxation delay, number of scans, and NMR windowing were optimized for data acquisition. For post-processing, the Lorentz/Gauss deconvolution method was employed to increase the signal accuracy by separating the impurities and noise in the integrated region of the quantitative proton. The method validation showed that the developed method has acceptable sensitivity, linearity, precision, and accuracy. The purity of the commercial Rg1 CRS examined with the method developed in this research was 90.34±0.21%, which was obviously lower than that reported by the manufacturer (>98.0%, HPLC-UV). The cross-method validation shows that the commonly used HPLC-UV, HPLC-ELSD (evaporative light scattering detector) and even LC-MS (mass spectrometry) methods provide significantly higher purity values of Rg1 CRS compared with the qNMR method, and the accuracy of these LC-based methods largely depend on the amount of the sample that was loaded and the properties of the impurities.

Glycan OH Exchange Rate Determination in Aqueous Solution: Seeking Evidence for Transient Hydrogen Bonds

Hydrogen bonds (Hbonds) are important stabilizing forces in biomolecules. However, for glycans in aqueous solution, direct NMR detection of Hbonds is elusive because of their transient nature. Here, we present Isotope-based Natural-abundance TOtal correlation eXchange SpectroscopY (INTOXSY), a new ¹H-¹³C heteronuclear single quantum coherence-total correlation spectroscopy based method, to extract OH groups' exchange rate constants (k_ex) for molecules in natural ¹³C abundance and show that OH Hbonds can be inferred from "slower" H/D k_ex. We evaluate k_ex measured with INTOXSY in light of those extracted with line-shape analysis. Subsequently, we use a set of common glycans to establish a k_ex reference basis set and to infer the existence of transient Hbonds involving OH donor groups. Then, we report k_ex values for a series of mono- and disaccharides, as well as for oligosaccharides sialyl Lewis X and β-cyclodextrin, and compare the results with those from the reference set to extract Hbond information. Finally, we utilize NMR experimental data in conjunction with molecular dynamics simulations to establish donor and acceptor Hbond pairs. Our exchange rate measurements indicate that OH/OD exchange rates, k_HD, values <10 s^-1 are consistent with transient Hbond OH groups and potential acceptor groups can be uncovered through MD simulations.

An intermolecular G-quadruplex as the basis for GTP recognition in the class V-GTP aptamer

Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V-GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.

Base-pair opening dynamics of primary miR156a using NMR elucidates structural determinants important for its processing level and leaf number phenotype in Arabidopsis

MicroRNAs originate from primary transcripts containing hairpin structures. The levels of mature miR156 influence the leaf number prior to flowering in the life cycle of plants. To understand the molecular mechanism of biogenesis of primary miR156a (pri-miR156a) to mature miR156, a base-pair opening dynamics study was performed using model RNAs mimicking the cleavage site of wild type and B5 bulge-stabilizing mutant pri-miR156a constructs. We also determined the mature miR156 levels and measured leaf numbers at flowering of plants overexpressing the wild type and mutant constructs. Our results suggest that the stabilities and/or opening dynamics of the C15·G98 and U16·A97 base-pairs at the cleavage site are essential for formation of the active conformation and for efficient processing of pri-miR156a, and that mutations of the B5 bulge can modulate mature miR156 levels as well as miR156-driven leaf number phenotypes via changes in the base-pair stability of the cleavage site.

Using NMR and molecular dynamics to link structure and dynamics effects of the universal base 8-aza, 7-deaza, N8 linked adenosine analog

A truly universal nucleobase enables a host of novel applications such as simplified templates for PCR primers, randomized sequencing and DNA based devices. A universal base must pair indiscriminately to each of the canonical bases with little or preferably no destabilization of the overall duplex. In reality, many candidates either destabilize the duplex or do not base pair indiscriminatingly. The novel base 8-aza-7-deazaadenine (pyrazolo[3,4-d]pyrimidin- 4-amine) N⁸-(2'deoxyribonucleoside), a deoxyadenosine analog (UB), pairs with each of the natural DNA bases with little sequence preference. We have utilized NMR complemented with molecular dynamic calculations to characterize the structure and dynamics of a UB incorporated into a DNA duplex. The UB participates in base stacking with little to no perturbation of the local structure yet forms an unusual base pair that samples multiple conformations. These local dynamics result in the complete disappearance of a single UB proton resonance under native conditions. Accommodation of the UB is additionally stabilized via heightened backbone conformational sampling. NMR combined with various computational techniques has allowed for a comprehensive characterization of both structural and dynamic effects of the UB in a DNA duplex and underlines that the UB as a strong candidate for universal base applications.

Structural Insights into Conformation Differences between DNA/TNA and RNA/TNA Chimeric Duplexes

Threose nucleic acid (TNA) is an artificial genetic polymer capable of heredity and evolution, and is studied in the context of RNA chemical etiology. It has a four-carbon threose backbone in place of the five-carbon ribose of natural nucleic acids, yet forms stable antiparallel complementary Watson-Crick homoduplexes and heteroduplexes with DNA and RNA. TNA base-pairs more favorably with RNA than with DNA but the reason is unknown. Here, we employed NMR, ITC, UV, and CD to probe the structural and dynamic properties of heteroduplexes of RNA/TNA and DNA/TNA. The results indicate that TNA templates the structure of heteroduplexes, thereby forcing an A-like helical geometry. NMR measurement of kinetic and thermodynamic parameters for individual base pair opening events reveal unexpected asymmetric "breathing" fluctuations of the DNA/TNA helix. The results suggest that DNA is unable to fully adapt to the conformational constraints of the rigid TNA backbone and that nucleic acid breathing dynamics are determined from both backbone and base contributions.

Observation of coherent delocalized phonon-like modes in DNA under physiological conditions

Underdamped terahertz-frequency delocalized phonon-like modes have long been suggested to play a role in the biological function of DNA. Such phonon modes involve the collective motion of many atoms and are prerequisite to understanding the molecular nature of macroscopic conformational changes and related biochemical phenomena. Initial predictions were based on simple theoretical models of DNA. However, such models do not take into account strong interactions with the surrounding water, which is likely to cause phonon modes to be heavily damped and localized. Here we apply state-of-the-art femtosecond optical Kerr effect spectroscopy, which is currently the only technique capable of taking low-frequency (GHz to THz) vibrational spectra in solution. We are able to demonstrate that phonon modes involving the hydrogen bond network between the strands exist in DNA at physiologically relevant conditions. In addition, the dynamics of the solvating water molecules is slowed down by about a factor of 20 compared with the bulk.

Fifty years of DNA "breathing": Reflections on old and new approaches

The coding sequences for genes, and much other regulatory information involved in genome expression, are located 'inside' the DNA duplex. Thus the "macromolecular machines" that read-out this information from the base sequence of the DNA must somehow access the DNA "interior." Double-stranded (ds) DNA is a highly structured and cooperatively stabilized system at physiological temperatures, but is also only marginally stable and undergoes a cooperative "melting phase transition" at temperatures not far above physiological. Furthermore, due to its length and heterogeneous sequence, with AT-rich segments being less stable than GC-rich segments, the DNA genome 'melts' in a multistate fashion. Therefore the DNA genome must also manifest thermally driven structural ("breathing") fluctuations at physiological temperatures that should reflect the heterogeneity of the dsDNA stability near the melting temperature. Thus many of the breathing fluctuations of dsDNA are likely also to be sequence dependent, and could well contain information that should be "readable" and useable by regulatory proteins and protein complexes in site-specific binding reactions involving dsDNA "opening." Our laboratory has been involved in studying the breathing fluctuations of duplex DNA for about 50 years. In this "Reflections" article we present a relatively chronological overview of these studies, starting with the use of simple chemical probes (such as hydrogen exchange, formaldehyde, and simple DNA "melting" proteins) to examine the local stability of the dsDNA structure, and culminating in sophisticated spectroscopic approaches that can be used to monitor the breathing-dependent interactions of regulatory complexes with their duplex DNA targets in "real time."

A historical account of Hoogsteen base-pairs in duplex DNA

In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.

The formation of catalytically competent enzyme-substrate complex is not a bottleneck in lesion excision by human alkyladenine DNA glycosylase

Human alkyladenine DNA glycosylase (AAG) protects DNA from alkylated and deaminated purine lesions. AAG flips out the damaged nucleotide from the double helix of DNA and catalyzes the hydrolysis of the N-glycosidic bond to release the damaged base. To understand better, how the step of nucleotide eversion influences the overall catalytic process, we performed a pre-steady-state kinetic analysis of AAG interaction with specific DNA-substrates, 13-base pair duplexes containing in the 7th position 1-N6-ethenoadenine (εA), hypoxanthine (Hx), and the stable product analogue tetrahydrofuran (F). The combination of the fluorescence of tryptophan, 2-aminopurine, and 1-N6-ethenoadenine was used to record conformational changes of the enzyme and DNA during the processes of DNA lesion recognition, damaged base eversion, excision of the N-glycosidic bond, and product release. The thermal stability of the duplexes characterized by the temperature of melting, T_m, and the rates of spontaneous opening of individual nucleotide base pairs were determined by NMR spectroscopy. The data show that the relative thermal stability of duplexes containing a particular base pair in position 7, (T_m(F/T) < T_m(εA/T) < T_m(Hx/T) < T_m(A/T)) correlates with the rate of reversible spontaneous opening of the base pair. However, in contrast to that, the catalytic lesion excision rate is two orders of magnitude higher for Hx-containing substrates than for substrates containing εA, proving that catalytic activity is not correlated with the stability of the damaged base pair. Our study reveals that the formation of the catalytically competent enzyme-substrate complex is not the bottleneck controlling the catalytic activity of AAG.

Solution structure of the Z-DNA binding domain of PKR-like protein kinase from Carassius auratus and quantitative analyses of the intermediate complex during B-Z transition

Z-DNA binding proteins (ZBPs) play important roles in RNA editing, innate immune response and viral infection. Structural and biophysical studies show that ZBPs initially form an intermediate complex with B-DNA for B-Z conversion. However, a comprehensive understanding of the mechanism of Z-DNA binding and B-Z transition is still lacking, due to the absence of structural information on the intermediate complex. Here, we report the solution structure of the Zα domain of the ZBP-containing protein kinase from Carassius auratus(caZαPKZ). We quantitatively determined the binding affinity of caZαPKZ for both B-DNA and Z-DNA and characterized its B-Z transition activity, which is modulated by varying the salt concentration. Our results suggest that the intermediate complex formed by caZαPKZ and B-DNA can be used as molecular ruler, to measure the degree to which DNA transitions to the Z isoform.

Direct ¹³C-detected NMR experiments for mapping and characterization of hydrogen bonds in RNA

In RNA secondary structure determination, it is essential to determine whether a nucleotide is base-paired and not. Base-pairing of nucleotides is mediated by hydrogen bonds. The NMR characterization of hydrogen bonds relies on experiments correlating the NMR resonances of exchangeable protons and can be best performed for structured parts of the RNA, where labile hydrogen atoms are protected from solvent exchange. Functionally important regions in RNA, however, frequently reveal increased dynamic disorder which often leads to NMR signals of exchangeable protons that are broadened beyond (1)H detection. Here, we develop (13)C direct detected experiments to observe all nucleotides in RNA irrespective of whether they are involved in hydrogen bonds or not. Exploiting the self-decoupling of scalar couplings due to the exchange process, the hydrogen bonding behavior of the hydrogen bond donor of each individual nucleotide can be determined. Furthermore, the adaption of HNN-COSY experiments for (13)C direct detection allows correlations of donor-acceptor pairs and the localization of hydrogen-bond acceptor nucleotides. The proposed (13)C direct detected experiments therefore provide information about molecular sites not amenable by conventional proton-detected methods. Such information makes the RNA secondary structure determination by NMR more accurate and helps to validate secondary structure predictions based on bioinformatics.

Base Flipping within the α-Hemolysin Latch Allows Single-Molecule Identification of Mismatches in DNA

A method for identifying and differentiating DNA duplexes containing the mismatched base pairs CC and CA at single molecule resolution with the protein pore α-hemolysin (αHL) is presented. Unique modulating current signatures are observed for duplexes containing the CC and CA mismatches when the mismatch site in the duplex is situated in proximity to the latch constriction of αHL during DNA residence inside the pore. The frequency and current amplitude of the modulation states are dependent on the mismatch type (CC or CA) permitting easy discrimination of these mismatches from one another, and from a fully complementary duplex that exhibits no modulation. We attribute the modulating current signatures to base flipping and subsequent interaction with positively charged lysine residues at the latch constriction of αHL. Our hypothesis is supported by the extended residence times of DNA duplexes within the pore when a mismatch is in proximity to the latch constriction, and by the loss of the two-state current signature in low pH buffers (<6.3), where the protonation of one of the cytosine bases increases the stability of the intrahelical state.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

Influence of 5-N-carboxamide modifications on the thermodynamic stability of oligonucleotides

We have recently shown that the incorporation of modified nucleotides such as 5-N-carboxamide-deoxyuridines into random nucleic acid libraries improves success rates in SELEX experiments and facilitates the identification of ligands with slow off-rates. Here we report the impact of these modifications on the thermodynamic stability of both duplexes and intramolecular ‘single-stranded’ structures. Within duplexes, large, hydrophobic naphthyl groups were destabilizing relative to the all natural DNA duplex, while the hydrophilic groups exhibited somewhat improved duplex stability. All of the significant changes in stability were driven by opposing contributions from the enthalpic and entropic terms. In contrast, both benzyl and naphthyl modifications stabilized intramolecular single-stranded structures relative to their natural DNA analogs, consistent with the notion that intramolecular folding allows formation of novel, stabilizing hydrophobic interactions. Imino proton NMR data provided evidence that elements of the folded structure form at temperatures well below the T_m, with a melting transition that is distinctly less cooperative when compared to duplex DNA. Although there are no data to suggest that the unmodified DNA sequences fold into structures similar to their modified analogs, this still represents clear evidence that these modifications impart thermodynamic stability to the folded structure not achievable with unmodified DNA.

Centromeric Alpha-Satellite DNA Adopts Dimeric i-Motif Structures Capped by AT Hoogsteen Base Pairs

Human centromeric alpha-satellite DNA is composed of tandem arrays of two types of 171 bp monomers; type A and type B. The differences between these types are concentrated in a 17 bp region of the monomer called the A/B box. Here, we have determined the solution structure of the C-rich strand of the two main variants of the human alpha-satellite A box. We show that, under acidic conditions, the C-rich strands of two A boxes self-recognize and form a head-to-tail dimeric i-motif stabilized by four intercalated hemi-protonated C:C(+) base pairs. Interestingly, the stack of C:C(+) base pairs is capped by T:T and Hoogsteen A:T base pairs. The two main variants of the A box adopt a similar three-dimensional structure, although the residues involved in the formation of the i-motif core are different in each case. Together with previous studies showing that the B box (known as the CENP-B box) also forms dimeric i-motif structures, our finding of this non-canonical structure in the A box shows that centromeric alpha satellites in all human chromosomes are able to form i-motifs, which consequently raises the possibility that these structures may play a role in the structural organization of the centromere.

Uncovering the thermodynamics of monomer binding for RNA replication

The nonenzymatic replication of primordial RNA is thought to have been a critical step in the origin of life. However, despite decades of effort, the poor rate and fidelity of model template copying reactions have thus far prevented an experimental demonstration of nonenzymatic RNA replication. The overall rate and fidelity of template copying depend, in part, on the affinity of free ribonucleotides to the RNA primer-template complex. We have now used (1)H NMR spectroscopy to directly measure the thermodynamic association constants, Kas, of the standard ribonucleotide monophosphates (rNMPs) to native RNA primer-template complexes. The binding affinities of rNMPs to duplexes with a complementary single-nucleotide overhang follow the order C > G > A > U. Notably, these monomers bind more strongly to RNA primer-template complexes than to the analogous DNA complexes. The relative binding affinities of the rNMPs for complementary RNA primer-template complexes are in good quantitative agreement with the predictions of a nearest-neighbor analysis. With respect to G:U wobble base-pairing, we find that the binding of rGMP to a primer-template complex with a 5'-U overhang is approximately 10-fold weaker than to the complementary 5'-C overhang. We also find that the binding of rGMP is only about 2-fold weaker than the binding of rAMP to 5'-U, consistent with the poor fidelity observed in the nonenzymatic copying of U residues in RNA templates. The accurate Ka measurements for ribonucleotides obtained in this study will be useful for designing higher fidelity, more effective RNA replication systems.

Mechanistic insights into temperature-dependent regulation of the simple cyanobacterial hsp17 RNA thermometer at base-pair resolution

The cyanobacterial hsp17 ribonucleicacid thermometer (RNAT) is one of the smallest naturally occurring RNAT. It forms a single hairpin with an internal 1×3-bulge separating the start codon in stem I from the ribosome binding site (RBS) in stem II. We investigated the temperature-dependent regulation of hsp17 by mapping individual base-pair stabilities from solvent exchange nuclear magnetic resonance (NMR) spectroscopy. The wild-type RNAT was found to be stabilized by two critical CG base pairs (C14-G27 and C13-G28). Replacing the internal 1×3 bulge by a stable CG base pair in hsp17^rep significantly increased the global stability and unfolding cooperativity as evidenced by circular dichroism spectroscopy. From the NMR analysis, remote stabilization and non-nearest neighbour effects exist at the base-pair level, in particular for nucleotide G28 (five nucleotides apart from the side of mutation). Individual base-pair stabilities are coupled to the stability of the entire thermometer within both the natural and the stabilized RNATs by enthalpy–entropy compensation presumably mediated by the hydration shell. At the melting point the Gibbs energies of the individual nucleobases are equalized suggesting a consecutive zipper-type unfolding mechanism of the RBS leading to a dimmer-like function of hsp17 and switch-like regulation behaviour of hsp17^rep. The data show how minor changes in the nucleotide sequence not only offset the melting temperature but also alter the mode of temperature sensing. The cyanobacterial thermosensor demonstrates the remarkable adjustment of natural RNATs to execute precise temperature control.

Differential stabilities and sequence-dependent base pair opening dynamics of Watson-Crick base pairs with 5-hydroxymethylcytosine, 5-formylcytosine, or 5-carboxylcytosine

5-Hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) form during active demethylation of 5-methylcytosine (5mC) and are implicated in epigenetic regulation of the genome. They are differentially processed by thymine DNA glycosylase (TDG), an enzyme involved in active demethylation of 5mC. Three modified Dickerson-Drew dodecamer (DDD) sequences, amenable to crystallographic and spectroscopic analyses and containing the 5'-CG-3' sequence associated with genomic cytosine methylation, containing 5hmC, 5fC, or 5caC placed site-specifically into the 5'-T(8)X(9)G(10)-3' sequence of the DDD, were compared. The presence of 5caC at the X(9) base increased the stability of the DDD, whereas 5hmC or 5fC did not. Both 5hmC and 5fC increased imino proton exchange rates and calculated rate constants for base pair opening at the neighboring base pair A(5):T(8), whereas 5caC did not. At the oxidized base pair G(4):X(9), 5fC exhibited an increase in the imino proton exchange rate and the calculated kop. In all cases, minimal effects to imino proton exchange rates occurred at the neighboring base pair C(3):G(10). No evidence was observed for imino tautomerization, accompanied by wobble base pairing, for 5hmC, 5fC, or 5caC when positioned at base pair G(4):X(9); each favored Watson-Crick base pairing. However, both 5fC and 5caC exhibited intranucleobase hydrogen bonding between their formyl or carboxyl oxygens, respectively, and the adjacent cytosine N(4) exocyclic amines. The lesion-specific differences observed in the DDD may be implicated in recognition of 5hmC, 5fC, or 5caC in DNA by TDG. However, they do not correlate with differential excision of 5hmC, 5fC, or 5caC by TDG, which may be mediated by differences in transition states of the enzyme-bound complexes.

NMR analysis of base-pair opening kinetics in DNA

Base pairing in nucleic acids plays a crucial role in their structure and function. Differences in the base-pair opening and closing kinetics of individual double-stranded DNA sequences or between chemically modified base pairs provide insight into the recognition of these base pairs by DNA processing enzymes. This unit describes how to quantify the kinetics for localized base pairs by observing changes in the imino proton signals by nuclear magnetic resonance spectroscopy. The determination of all relevant parameters using state-of-the art techniques and NMR instrumentation, including cryoprobes, is discussed.

Biosynthetic modularity rules in the bisintercalator family of antitumor compounds

Diverse actinomycetes produce a family of structurally and biosynthetically related non-ribosomal peptide compounds which belong to the chromodepsipeptide family. These compounds act as bisintercalators into the DNA helix. They give rise to antitumor, antiparasitic, antibacterial and antiviral bioactivities. These compounds show a high degree of conserved modularity (chromophores, number and type of amino acids). This modularity and their high sequence similarities at the genetic level imply a common biosynthetic origin for these pathways. Here, we describe insights about rules governing this modular biosynthesis, taking advantage of the fact that nowadays five of these gene clusters have been made public (thiocoraline, triostin, SW-163 and echinomycin/quinomycin). This modularity has potential application for designing and producing novel genetic engineered derivatives, as well as for developing new chemical synthesis strategies. These would facilitate their clinical development.

Mechanistic aspects of hydration of guanine radical cations in DNA

The mechanistic aspects of hydration of guanine radical cations, G(•+) in double- and single-stranded oligonucleotides were investigated by direct time-resolved spectroscopic monitoring methods. The G(•+) radical one-electron oxidation products were generated by SO4(•-) radical anions derived from the photolysis of S2O8(2-) anions by 308 nm laser pulses. In neutral aqueous solutions (pH 7.0), after the complete decay of SO4(•-) radicals (∼5 μs after the actinic laser flash) the transient absorbance of neutral guanine radicals, G(-H)(•) with maximum at 312 nm, is dominant. The kinetics of decay of G(-H)(•) radicals depend strongly on the DNA secondary structure. In double-stranded DNA, the G(-H)(•) decay is biphasic with one component decaying with a lifetime of ∼2.2 ms and the other with a lifetime of ∼0.18 s. By contrast, in single-stranded DNA the G(-H)(•) radicals decay monophasically with a ∼ 0.28 s lifetime. The ms decay component in double-stranded DNA is correlated with the enhancement of 8-oxo-7,8-dihydroguanine (8-oxoG) yields which are ∼7 greater than in single-stranded DNA. In double-stranded DNA, it is proposed that the G(-H)(•) radicals retain radical cation character by sharing the N1-proton with the N3-site of C in the [G(•+):C] base pair. This [G(-H)(•):H(+)C ⇆ G(•+):C] equilibrium allows for the hydration of G(•+) followed by formation of 8-oxoG. By contrast, in single-stranded DNA, deprotonation of G(•+) and the irreversible escape of the proton into the aqueous phase competes more effectively with the hydration mechanism, thus diminishing the yield of 8-oxoG, as observed experimentally.