Probing Site-Specific Energetics in Proteins and Nucleic Acids by Hydrogen Exchange and Nuclear Magnetic Resonance Spectroscopy

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

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.

Systems NMR: single-sample quantification of RNA, proteins and metabolites for biomolecular network analysis

Cellular behavior is controlled by the interplay of diverse biomolecules. Most experimental methods, however, can only monitor a single molecule class or reaction type at a time. We developed an in vitro nuclear magnetic resonance spectroscopy (NMR) approach, which permitted dynamic quantification of an entire 'heterotypic' network-simultaneously monitoring three distinct molecule classes (metabolites, proteins and RNA) and all elementary reaction types (bimolecular interactions, catalysis, unimolecular changes). Focusing on an eight-reaction co-transcriptional RNA folding network, in a single sample we recorded over 35 time points with over 170 observables each, and accurately determined five core reaction constants in multiplex. This reconstruction revealed unexpected cross-talk between the different reactions. We further observed dynamic phase-separation in a system of five distinct RNA-binding domains in the course of the RNA transcription reaction. Our Systems NMR approach provides a deeper understanding of biological network dynamics by combining the dynamic resolution of biochemical assays and the multiplexing ability of 'omics'.

NMR characterization of RNA small molecule interactions

Exciting discoveries of naturally occurring ligand-sensing and disease-linked noncoding RNAs have promoted significant interests in understanding RNA-small molecule interactions. NMR spectroscopy is a powerful tool for characterizing intermolecular interactions. In this review, we describe protocols and approaches for applying NMR spectroscopy to investigate interactions between RNA and small molecules. We review protocols for RNA sample preparation, methods for identifying RNA-binding small molecules, approaches for mapping RNA-small molecule interactions, determining complex structures, and characterizing binding kinetics. We hope this review will provide a guideline to streamline NMR applications in studying RNA-small molecule interactions, facilitating both basic mechanistic understandings of RNA functions and translational efforts in developing RNA-targeted therapeutics.

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.

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.

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.

Characterizing slow chemical exchange in nucleic acids by carbon CEST and low spin-lock field R(1ρ) NMR spectroscopy

Quantitative characterization of dynamic exchange between various conformational states provides essential insights into the molecular basis of many regulatory RNA functions. Here, we present an application of nucleic-acid-optimized carbon chemical exchange saturation transfer (CEST) and low spin-lock field R(1ρ) relaxation dispersion (RD) NMR experiments in characterizing slow chemical exchange in nucleic acids that is otherwise difficult if not impossible to be quantified by the ZZ-exchange NMR experiment. We demonstrated the application on a 47-nucleotide fluoride riboswitch in the ligand-free state, for which CEST and R(1ρ) RD profiles of base and sugar carbons revealed slow exchange dynamics involving a sparsely populated (p ~ 10%) and shortly lived (τ ~ 10 ms) NMR "invisible" state. The utility of CEST and low spin-lock field R(1ρ) RD experiments in studying slow exchange was further validated in characterizing an exchange as slow as ~60 s(-1).

Opening dynamics of 8-oxoguanine in DNA

8-oxoguanine is a major lesion of genomic DNA that results from oxidation of guanine by reactive oxygen species. The repair of this lesion is initiated by 8-oxoguanine glycosylases, which excise the damaged base by "flipping" it outside the DNA double helix. The molecular mechanisms involved in the specific recognition of the damaged base by the enzyme are not yet fully understood. Several models have proposed that, in DNA, the base pair between 8-oxoguanine and cytosine may possess altered dynamic properties that could help the enzyme locate the lesion and could favor the selective extra-helical flipping of the damaged base. To test this proposal, we have characterized the spontaneous opening of the base pair between 8-oxoguanine and cytosine in a DNA double helix using NMR spectroscopy and proton exchange. The results show that the rate of spontaneous opening of 8-oxoguanine and the lifetime of the base in the extra-helical state are the same as those of a canonical guanine-cytosine base pair, in the same base sequence context. This finding suggests that the opening dynamics of 8-oxoguanine, when paired with cytosine in DNA, does not play a significant role in the recognition of the lesion by glycosylases.

Individual basepair stability of DNA and RNA studied by NMR-detected solvent exchange

In this study, we have optimized NMR methodology to determine the thermodynamic parameters of basepair opening in DNA and RNA duplexes by characterizing the temperature dependence of imino proton exchange rates of individual basepairs. Contributions of the nuclear Overhauser effect to exchange rates measured with inversion recovery experiments are quantified, and the influence of intrinsic and external catalysis exchange mechanisms on the imino proton exchange rates is analyzed. Basepairs in DNA and RNA have an approximately equal stability, and the enthalpy and entropy values of their basepair dissociation are correlated linearly. Furthermore, the compensation temperature, T(c), which is derived from the slope of the correlation, coincides with the melting temperature, and duplex unfolding occurs at that temperature where all basepairs are equally thermodynamically stable. The impact of protium-deuterium exchange of the imino hydrogen on the free energy of RNA basepair opening is investigated, and it is found that two A·U basepairs show distinct fractionation factors.

Probing protein interactions with hydrogen/deuterium exchange and mass spectrometry-a review

Assessing the functional outcome of protein interactions in structural terms is a goal of structural biology, however most techniques have a limited capacity for making structure-function determinations with both high resolution and high throughput. Mass spectrometry can be applied as a reader of protein chemistries in order to fill this void, and enable methodologies whereby protein structure-function determinations may be made on a proteome-wide level. Protein hydrogen/deuterium exchange (H/DX) offers a chemical labeling strategy suitable for tracking changes in "dynamic topography" and thus represents a powerful means of monitoring protein structure-function relationships. This review presents the exchange method in the context of interaction analysis. Applications involving interface detection, quantitation of binding, and conformational responses to ligation are discussed, and commentary on recent analytical developments is provided.

Characterizing RNA dynamics at atomic resolution using solution-state NMR spectroscopy

Many recently discovered noncoding RNAs do not fold into a single native conformation but sample many different conformations along their free-energy landscape to carry out their biological function. Here we review solution-state NMR techniques that measure the structural, kinetic and thermodynamic characteristics of RNA motions spanning picosecond to second timescales at atomic resolution, allowing unprecedented insights into the RNA dynamic structure landscape. From these studies a basic description of the RNA dynamic structure landscape is emerging, bringing new insights into how RNA structures change to carry out their function as well as applications in RNA-targeted drug discovery and RNA bioengineering.

Enhanced base-pair opening in the adenine tract of a RNA double helix

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the kinetics and energetics of base-pair opening in two nucleic acid double helices. One is the RNA duplex 5'-r(GCGAUAAAAAGGCC)-3'/5'-r(GGCCUUUUUAUCGC)-3', which contains a central tract of five AU base pairs. The other is the homologous DNA duplex with a central tract of five AT base pairs. The rates and the equilibrium constants of the opening reaction of each base pair are measured from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results reveal that the tract of AU base pairs in the RNA duplex differs from the homologous tract of AT base pairs in DNA in several ways. The rates of opening of AU base pairs in RNA are high and increase progressively along the tract, reaching their largest values at the 3'-end of the tract. In contrast, the opening rates of AT base pairs in DNA are much lower than those of AU base pairs. Within the tract, the largest opening rate is observed for the AT base pair at the 5'-end of the tract. These differences in opening kinetics are paralleled by differences in the stabilities of individual base pairs. All AU base pairs in the RNA are less stable than the AT base pairs in the DNA. The presence of the tract enhances these differences by increasing the stability of AT base pairs in DNA while decreasing the stability of AU base pairs in RNA. Due to these divergent trends, along the tracts, the AU base pairs become progressively less stable than AT base pairs. These findings demonstrate that tracts of AU base pairs in RNA have specific dynamic and energetic signatures that distinguish them from similar tracts of AT base pairs in DNA.

Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5'-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine-Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14-C25 base pair. The mismatch base pair in the wt RNA thermometer (A8-G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.

Structural energetics of the adenine tract from an intrinsic transcription terminator

Intrinsic transcription termination sites generally contain a tract of adenines in the DNA template that yields a tract of uracils at the 3' end of the nascent RNA. To understand how this base sequence contributes to termination of transcription, we have investigated two nucleic acid structures. The first is the RNA-DNA hybrid that contains the uracil tract 5'-rUUUUUAU-3' from the tR2 intrinsic terminator of bacteriophage lambda. The second is the homologous DNA-DNA duplex that contains the adenine tract 5'-dATAAAAA-3'. This duplex is present at the tR2 site when the DNA is not transcribed. The opening and the stability of each rU-dA/dT-dA base pair in the two structures are characterized by imino proton exchange and nuclear magnetic resonance spectroscopy. The results reveal concerted opening of the central rU-dA base pairs in the RNA-DNA hybrid. Furthermore, the stability profile of the adenine tract in the RNA-DNA hybrid is very different from that of the tract in the template DNA-DNA duplex. In the RNA-DNA hybrid, the stabilities of rU-dA base pairs range from 4.3 to 6.5 kcal/mol (at 10 degrees C). The sites of lowest stability are identified at the central positions of the tract. In the template DNA-DNA duplex, the dT-dA base pairs are more stable than the corresponding rU-dA base pairs in the hybrid by 0.9 to 4.6 kcal/mol and, in contrast to the RNA-DNA hybrid, the central base pairs have the highest stability. These results suggest that the central rU-dA/dT-dA base pairs in the adenine tract make the largest energetic contributions to transcription termination by promoting both the dissociation of the RNA transcript and the closing of the transcription bubble. The results also suggest that the high stability of dT-dA base pairs in the DNA provides a signal for the pausing of RNA polymerase at the termination site.

Dynamics and stability of individual base pairs in two homologous RNA-DNA hybrids

Nuclear magnetic resonance spectroscopy and proton exchange have been used to characterize two RNA-DNA hybrids from the tR2 intrinsic transcription terminator site of phage lambda. The hybrids have the same base sequence [5'-GGCGCAGGCC(T/U)(T/U)CC-3'/5'-GGAAGGCC(T/U)GCGCC-3'] but differ from each other by an interchange of DNA and RNA strands. The opening of single base pairs in the two hybrids is characterized by measuring the rates of exchange of imino protons with solvent protons as a function of the concentration of a proton acceptor (ammonia base) at 10 degrees C. The free energy change in the opening reaction provides a measure of the stability of the base pair, while the rates of opening and closing define the base pair dynamics. The results demonstrate that, within the same base sequence context, dA-rU base pairs are less stable than dT-rA base pairs. The differences in stability are enhanced when two dA-rU base pairs are located next to each other in the hybrid structure. For the G-C base pairs, the rates of opening and closing and the stability are affected by the base sequence context and by the nature of the sugar moiety attached to the guanine. The dominant feature of the base sequence is the proximity of the dA-rU base pair, which destabilizes the G-C base pair when the guanine is located on the DNA strand. Two G-C base pairs (namely, those in the fourth and 10th positions) exhibit large differences in their opening and closing rates between the two hybrids, while maintaining the same stability. These results provide the first demonstration that, for RNA-DNA hybrid structures with the same base sequence, the opening dynamics and the stability of individual base pairs are strongly influenced by the chemical nature of each strand.

Bending modes of DNA directly addressed by cryo-electron microscopy of DNA minicircles

We use cryo-electron microscopy (cryo-EM) to study the 3D shapes of 94-bp-long DNA minicircles and address the question of whether cyclization of such short DNA molecules necessitates the formation of sharp, localized kinks in DNA or whether the necessary bending can be redistributed and accomplished within the limits of the elastic, standard model of DNA flexibility. By comparing the shapes of covalently closed, nicked and gapped DNA minicircles, we conclude that 94-bp-long covalently closed and nicked DNA minicircles do not show sharp kinks while gapped DNA molecules, containing very flexible single-stranded regions, do show sharp kinks. We corroborate the results of cryo-EM studies by using Bal31 nuclease to probe for the existence of kinks in 94-bp-long minicircles.

Influence of magnesium ions on spontaneous opening of DNA base pairs

A large amount of experimental evidence is available for the effects of magnesium ions on the structure and the stability of the DNA double helix. Less is known, however, on how these ions affect the dynamics of the molecule and the stability of each individual base pair. The present work addresses these questions by a study of the DNA duplex [dCGCAGATCTGCG]2, and its interactions with magnesium ions using nuclear magnetic resonance (NMR) spectroscopy and proton exchange. Two-dimensional NMR experiments indicate that binding of magnesium to this DNA duplex does not affect its structure. However, even in the absence of structural changes, magnesium ions specifically affect the exchange properties of imino protons in the four GC/CG base pairs that are located in the interior of the double helix. These specific changes do not result from alterations in the rates of spontaneous opening of these base pairs. Instead, the changes most likely reflect an enhancement in the energetic propensity for spontaneous opening of the GC/CG base pairs that is induced by the binding of magnesium ions.

Exchange rate constants of invisible protons in proteins determined by NMR spectroscopy

Although labile protons that are exchanging rapidly with those of the solvent cannot be observed directly, their exchange rate constants can be determined by indirect detection of scalar-coupled neighboring nuclei. We have used heteronuclear NMR spectroscopy to measure the exchange rate constants of labile protons in the side chains of lysine and arginine residues in ubiquitin enriched in carbon-13 and nitrogen-15 at neutral pH. Exchange rate constants as fast as 40x10(3) s(-1) were thus measured. These results demonstrate that NMR spectroscopy is a powerful tool for the characterization of lysine NH3(+) and arginine NH groups in proteins at physiologically relevant pH values.

Kinking the double helix by bending deformation

DNA bending and torsional deformations that often occur during its functioning inside the cell can cause local disruptions of the regular helical structure. The disruptions created by negative torsional stress have been studied in detail, but those caused by bending stress have only been analyzed theoretically. By probing the structure of very small DNA circles, we determined that bending stress disrupts the regular helical structure when the radius of DNA curvature is smaller than 3.5 nm. First, we developed an efficient method to obtain covalently closed DNA minicircles. To detect structural disruptions in the minicircles we treated them by single-strand-specific endonucleases. The data showed that the regular DNA structure is disrupted by bending deformation in the 64-65-bp minicircles, but not in the 85-86-bp minicircles. Our results suggest that strong DNA bending initiates kink formation while preserving base pairing.

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.

Examination of the long-range effects of aminofluorene-induced conformational heterogeneity and its relevance to the mechanism of translesional DNA synthesis

Adduct-induced conformational heterogeneity complicates the understanding of how DNA adducts exert mutation. A case in point is the N-deacetylated AF lesion [N-(2'-deoxyguanosin-8-yl)-2-aminofluorene], the major adduct derived from the strong liver carcinogen N-acetyl-2-aminofluorene. Three conformational families have been previously characterized and are dependent on the positioning of the aminofluorene rings: B is in the "B-DNA" major groove, S is "stacked" into the helix with base-displacement, and W is "wedged" into the minor groove. Here, we conducted (19)F NMR, CD, T(m), and modeling experiments at various primer positions with respect to a template modified by a fluorine tagged AF-adduct (FAF). In the first set, the FAF-G was paired with C and in the second set it was paired with A. The FAF-G:C oligonucleotides were found to preferentially adopt the B or S-conformers while the FAF-G:A mismatch ones preferred the B and W-conformers. The conformational preferences of both series were dependent on temperature and complementary strand length; the largest differences in conformation were displayed at lower temperatures. The CD and T(m) results are in general agreement with the NMR data. Molecular modeling indicated that the aminofluorene moiety in the minor groove of the W-conformer would impose a steric clash with the tight-packing amino acid residues on the DNA binding area of the Bacillus fragment (BF), a replicative DNA polymerase. In the case of the B-type conformer, the carcinogenic moiety resides in the solvent-exposed major groove throughout the replication/translocation process. The present dynamic NMR results, combined with previous primer extension kinetic data by Miller & Grollman, support a model in which adduct-induced conformational heterogeneities at positions remote from the replication fork affect polymerase function through a long-range DNA-protein interaction.

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.

A nuclear magnetic resonance investigation of the energetics of basepair opening pathways in DNA

The opening of basepairs plays a key role in DNA replication and transcription, and in the action of DNA repair and modification enzymes. In this article, we have used proton exchange to define the energetics of the pathways for basepair opening in two DNA 17-mer duplexes. The rates of exchange of imino protons with solvent protons were measured by NMR spectroscopy for each DNA duplex, as a function of the concentration of exchange catalyst and of temperature. The measurements provided the rates and the equilibrium constants of the opening reactions for individual basepairs at different temperatures. These temperature dependences were used to calculate the enthalpies and the free energies of the barrier to opening and of the open state for each basepair. The results reveal the existence of three distinct patterns of enthalpy changes in the opening reactions. The patterns differ from each other in the location of the kinetic opening barrier relative to the open state. Neighboring bases, which are one or more positions removed from the opening basepair, influence the enthalpic pattern of the opening pathway. The free energies of the opening barriers are found to be linearly related to the free energies of the open state. This correlation is analyzed in terms of rate-equilibrium free energy relationships previously observed in other systems, and suggests that the transition state in the opening reaction is closer to the native closed state of the basepair than to its open state.

NMR study of the folding-unfolding mechanism for the thrombin-binding DNA aptamer d(GGTTGGTGTGGTTGG)

Hydrogen exchange rates of the imino protons of the thrombin-binding 15 mer DNA aptamer d(G(1)G(2)T(3)T(4)G(5)G(6)T(7)G(8)T(9)G(10)G(11)T(12)T(13)G(14)G(15)) in the presence of Sr(2+) were measured. In the temperature range 15-35 degrees C, the exchange rates of the eight iminos in the quadruplex core were not uniform, with the G(2), G(11) and G(15) iminos exchanging faster, the G(1), G(5), G(10) and G(14) iminos exchanging slower, and the G(6) imino exchanging at a medium rate. In the quadruplex G(1), G(5), G(10) and G(14) adopted syn glycosidic conformation, while G(2), G(6), G(11) and G(15) adopted anti-conformation. It was found that the four slowly exchanging iminos, which were all the syn-iminos, happened to be located in the TT loops that were not easy to open to the solvent. The anti-iminos exchanged faster, but the G(6) imino exchanged slower than other anti-iminos, because its hydrogen bond with the G(10)O6 was stabilized by the TGT loop. The fact that the G(6) imino exchanged at a faster rate than those syn-iminos in the TT loops suggested that the TGT loop was less stable than the TT loops. Unfolding mechanism for the quadruplex was thus proposed: The quadruplex first uncoupled the three base pairs: G(1)-G(15), G(2)-G(14) and G(5)-G(11), which were not protected by any loops. Then it opened the TGT loop. Finally, it opened the TT loops and the sequence became an unstructured random coil that exchanged with the quadruplex conformation. The conformational exchange between the quadruplex and random coil had been detected.

Base Pair Opening in Three DNA-unwinding Elements

DNA-unwinding elements are specific base sequences that are located in the origin of DNA replication where they provide the start point for strand separation and unwinding of the DNA double helix. In the present work we have obtained the first characterization of the opening of individual base pairs in DNA-unwinding elements. The three DNA molecules investigated reproduce the 13-mer DNA-unwinding elements present in the Escherichia coli chromosome. The base sequences of the three 13-mers are conserved in the origins of replication of enteric bacterial chromosomes. The exchange of imino protons with solvent protons was measured for each DNA as a function of the concentration of exchange catalyst using nuclear magnetic resonance spectroscopy. The exchange rates provided the rates and the equilibrium constants for opening of individual base pairs in each DNA at 20 degrees C. The results reveal that the kinetics and energetics of the opening reactions for AT/TA base pairs are different in the three DNA-unwinding elements due to long range effects of the base sequence. These differences encompass the AT/TA base pairs that are conserved in various bacterial genomes. Furthermore, a qualitative correlation is observed between the kinetics and energetics of opening of AT/TA base pairs and the location of the corresponding DNA-unwinding element in the origin of DNA replication.

Sequence-dependence of the energetics of opening of at basepairs in DNA

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the energetics of opening of AT/TA basepairs in the DNA dodecamer 5'-d(GCTATAAAAGGG)-3'/5'-d(CCCTTTTATAGC)-3'. The dodecamer contains the TATA box of the adenovirus major late promoter. The equilibrium constants for opening of each basepair are measured from the dependence of the exchange rates of imino protons on ammonia concentration. The enthalpy, entropy, and free energy changes in the opening reaction of each basepair are determined from the temperature dependence of the exchange rates. The results reveal that the opening enthalpy changes encompass a wide range of values, namely, from 17 to 29 kcal/mol. The largest values are observed for the AT basepairs in 7th and 8th positions. These values, and the exchange rates of the corresponding imino protons, suggest that these two basepairs open in a single concerted reaction. The enthalpy changes for opening of the central six basepairs are correlated to the opening entropy changes. This enthalpy-entropy compensation minimizes the variations in the opening free energies among these central basepairs. Deviations from the enthalpy-entropy compensation pattern are observed for basepairs located close to the ends of the duplex structure, suggesting a different mode of opening for these basepairs.