Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

Estimation of strength in different extra Watson-Crick hydrogen bonds in DNA double helices through quantum chemical studies

It was shown earlier, from database analysis, model building studies, and molecular dynamics simulations that formation of cross-strand bifurcated or Extra Watson-Crick hydrogen (EWC) bonds between successive base pairs may lead to extra rigidity to DNA double helices of certain sequences. The strengths of these hydrogen bonds are debatable, however, as they do not have standard linear geometry criterion. We have therefore carried out detailed ab initio quantum chemical studies using RHF/6-31G(2d,2p) and B3LYP/6-31G(2p,2d) basis sets to determine strengths of several bent hydrogen bonds with different donor and acceptors. Interaction energy calculations, corrected for the basis set superposition errors, suggest that N-H...O type bent EWC hydrogen bonds are possible along same strands or across the strands between successive base pairs, leading to significant stability (ca. 4-9 kcal/mol). The N-H...N and C-H...O type interactions, however, are not so stabilizing. Hence, consideration of EWC N-H...O H-bonds can lead to a better understanding of DNA sequence directed structural features.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

NMR studies of the stable mismatch purine-thymine in the self-complementary d(CGPuAATTTCG) duplex in solution

One- and two-dimensional nuclear Overhauser effect experiments demonstrate that a single hydrogen bond between a T imino proton and purine N3 is sufficient to hold the base pair dPu.dT in d(CGPuAATTTCG) by a Watson-Crick fashion rather than a Hoogsteen type. In addition, the dPu.dT base pair is well stacked with neighboring base pairs. The spin-lattice relaxation measurements at 30 and 35 degrees C of two decamers, d(CGPuAATTTCG) and d(CGAAATTTCG), reveal that the elimination of two single hydrogen bonds of dA.dT base pairs (due to the substitution of adenine for purine) in the sequence results in an increase in the overall imino proton exchange rate from 7 to 36 s-1 at the site of mismatch.     

Kinetics for exchange of imino protons in the d(C-G-C-G-A-A-T-T-C-G-C-G) double helix and in two similar helices that contain a G . T base pair, d(C-G-T-G-A-A-T-T-C-G-C-G), and an extra adenine, d(C-G-C-A-G-A-A-T-T-C-G-C-G)

The relaxation lifetimes of imino protons from individual base pairs were measured in (I) a perfect helix, d(C-G-C-G-A-A-T-T-C-G-C-G), (II) this helix with a G . C base pair replaced with a G . T base pair, d(C-G-T-G-A-A-T-T-C-G-C-G), and (III) the perfect helix with an extra adenine base in a mismatch, d(C-G-C-A-G-A-A-T-T-C-G-C-G). The lifetimes were measured by saturation recovery proton nuclear magnetic resonance experiments performed on the imino protons of these duplexes. The measured lifetimes of the imino protons were shown to correspond to chemical exchange lifetimes at higher temperatures and spin-lattice relaxation times at lower temperatures. Comparison of the lifetimes in these duplexes showed that the destabilizing effect of the G . T base pair in II affected the opening rate of only the nearest-neighbor base pairs. For helix III, the extra adenine affected the opening rates of all the base pairs in the helix and thus was a larger perturbation for opening of the base pairs than the G . T base pair. The temperature dependence of the exchange rates of the imino proton in the perfect helix gives values of 14-15 kcal/mol for activation energies of A . T imino protons. These relaxation rates were shown to correspond to exchange involving individual base pair opening in this helix, which means that one base-paired imino proton can exchange independent of the others. For the other two helices that contain perturbations, much larger activation energies for exchange of the imino protons were found, indicating that a cooperative transition involving exchange of at least several base pairs was the exchange mechanism of the imino protons. The effects of a perturbation in a helix on the exchange rates and the mechanisms for exchange of imino protons from oligonucleotide helices are discussed.     

The contribution of pseudouridine to stabilities and structure of RNAs

Thermodynamic data are reported revealing that pseudouridine (Ψ) can stabilize RNA duplexes when replacing U and forming Ψ-A, Ψ-G, Ψ-U and Ψ-C pairs. Stabilization is dependent on type of base pair, position of Ψ within the RNA duplex, and type and orientation of adjacent Watson–Crick pairs. NMR spectra demonstrate that for internal Ψ-A, Ψ-G and Ψ-U pairs, the N3 imino proton is hydrogen bonded to the opposite strand nucleotide and the N1 imino proton may also be hydrogen bonded. CD spectra show that general A-helix structure is preserved, but there is some shifting of peaks and changing of intensities. Ψ has two hydrogen donors (N1 and N3 imino protons) and two hydrogen bond acceptors because the glycosidic bond is C-C rather than C-N as in uridine. This greater structural potential may allow Ψ to behave as a kind of structurally driven universal base because it can enhance stability relative to U when paired with A, G, U or C inside a double helix. These structural and thermodynamic properties may contribute to the biological functions of Ψ.     

Can G-C Hoogsteen-wobble pairs contribute to the stability of d(G. C-C) triplexes?

Quantum mechanics, molecular dynamics and statistical mechanics methods are used to analyze the importance of neutral Hoogsteen-wobble G.C pairing in the stabilization of triple helices based on the poly-(G.C-C) trio at neutral pH and low ionic strength. In spite of the existence of a single hydrogen bond, the Hoogsteen-wobble G.C pair is found to be quite stable both in gas phase and solvated DNA. Molecular dynamics simulations of different triplexes based on the d(G.C-C) trio leads to stable structures if the neutral d(G.C-C) steps stabilized by Hoogsteen-wobble pairs are mixed with d(G.C-C+) steps. Finally, high level ab initio calculations and thermodynamic integration techniques are used to determine the relative stability of G.C wobble and G.C imino pairings. It is found that triplexes containing the imino pairing are slightly more stable structures than those with the wobble one, due mainly to a better stacking.     

Studies on the cross stand hydrogen bonds in DNA double helices.

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross-strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may in fact be considered as yet another general structural feature of DNA helices.     

Proton nuclear magnetic resonance investigation of the conformation and dynamics in the synthetic deoxyribonucleic acid decamers d(ATATCGATAT) and d(ATATGCATAT)

A variety of one-dimensional proton NMR methods have been used to investigate the properties of two synthetic DNA decamers, d(ATATCGATAT) and d(ATATGCATAT). These results, in conjunction with the results of two-dimensional NMR experiments, permit complete assignment of the base proton resonances. Low-field resonances were assigned by sequential "melting" of the A . T base pairs and by comparison of the spectra of the two decamers. Below 20 degree C spin-lattice relaxation is dominated by through-space dipolar interactions. A substantial isotope effect on the G imino proton relaxation is observed in 75% D2O, confirming the importance of the exchangeable amino protons in the relaxation process. A somewhat smaller isotope effect is observed on the T imino proton relaxation. At elevated temperatures spin-lattice relaxation of the imino protons is due to proton exchange with solvent. Apparent activation energies for exchange vary from 36 kcal/base pair for base pairs (3,8) to 64 kcal/mol for the most interior base pairs (5,6), indicating that disruption of part, or all, of the double helix contributes significantly to the exchange of the imino protons in these decamers. By contrast, single base pair opening events are the major low-temperature pathways for exchange from A X T and G X C base pairs in the more stable higher molecular weight DNA examined in other studies. The temperature dependence of the chemical shifts and line widths of certain aromatic resonances indicates that the interconversion between the helix and coil states is not in fast exchange below the melting temperature, Tm. Within experimental error, no differential melting of base pairs was found in either molecule, and both exhibited melting points Tm = 50-52 degrees C. Spin-spin and spin-lattice relaxation rates of the nonexchangeable protons (TH6, AH8, and AH2) are consistent with values calculated by using an isotropic rotor model with a rotational correlation time of 6 ns and interproton distances appropriate for B-family DNA. The faster decay of AH8 compared with GH8 is attributed to an interaction between the thymine methyl protons and the AH8 protons in adjacent adenines (5'ApT3'). The base protons (AH8, GH8, and TH6) appear to be located close (1.9-2.3 A) to sugar H2',2" protons.(ABSTRACT TRUNCATED AT 400 WORDS)     

Structural Studies of DNA Fragments: The G·T Wobble Base Pair in A,B and Z DNA; The G·A Base Pair in B-DNA

Abstract

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G·T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with 02 of T and 06 of G with N3 of T. The X-ray analyses establish that the G·T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G·A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Melting of Cross-Linked DNA V. Cross-Linking Effect Caused by Local Stabilization of the Double Helix

Abstract

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this “cross-linking effect” (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the “cross-linking effect” (18–20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20–40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Melting of cross-linked DNA v. cross-linking effect caused by local stabilization of the double helix

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this "cross-linking effect" (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the "cross-linking effect" (18 - 20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20- 40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Structural and dynamic properties of a fluorouracil-adenine base pair in DNA studied by proton NMR

An oligodeoxynucleotide duplex containing the chemotherapeutic agent 5-fluorouracil (FU) has been constructed by solid phase phosphotriester synthesis and has been studied in solution by proton NMR. In this study, we provide the first structural characterization of a DNA complex containing a FU.A base pair. It has been determined that the 7-mer duplex containing a central FU.A base pair adopts a normal right-handed configuration and the A residue in the FU.A pair is oriented in the normal anticonfiguration giving a Watson-Crick base pair. The significant difference between T.A and FU.A base pairs is dynamic, not structural: the FU.A base pair opens faster than normal base pairs in the oligonucleotide studied. We provide evidence that the FU.A base pair has a significantly enhanced opening rate resulting form decreased stacking of the 5-fluorouracil residue and not from the enhanced acidity of the 5'-fluorouracil imino proton.     

High base pair opening rates in tracts of GC base pairs

Sequence-dependent structural features of the DNA double helix have a strong influence on the base pair opening dynamics. Here we report a detailed study of the kinetics of base pair breathing in tracts of GC base pairs in DNA duplexes derived from 1H NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia. In the limit of infinite exchange catalyst concentration, the exchange times of the guanine imino protons of the GC tracts extrapolate to much shorter base pair lifetimes than commonly observed for isolated GC base pairs. The base pair lifetimes in the GC tracts are below 5 ms for almost all of the base pairs. The unusually rapid base pair opening dynamics of GC tracts are in striking contrast to the behavior of AT tracts, where very long base pair lifetimes are observed. The implication of these findings for the structural principles governing spontaneous helix opening as well as the DNA-binding specificity of the cytosine-5-methyltransferases, where flipping of the cytosine base has been observed, are discussed.     

Structural energetics and base-pair opening dynamics in sarcin-ricin domain RNA

The sarcin-ricin domain is a universal element of the RNA from the large ribosomal subunit. The domain is part of the binding site for elongation factors and is specifically cleaved by the toxins alpha-sarcin and ricin. In this work, we have mapped the energetics and dynamics of individual structural motifs in a 29-mer RNA oligomer containing the sarcin-ricin domain. The stability of individual base pairs in the structure was characterized from measurements of the exchange rates of imino protons using nuclear magnetic resonance spectroscopy at 10 degrees C. The measurements also provided the rates of opening and closing for selected base pairs. The results reveal that the structural stabilization free energies in the sarcin-ricin domain are broadly distributed between 2.9 and 10.6 kcal/mol. One of the least stable sites in the structure is the noncanonical G-A base pair located next to the phosphodiester bond that is cleaved by alpha-sarcin. The low stability of this base pair supports the proposal that cleavage by alpha-sarcin occurs by a base flipping mechanism. The opening dynamics of other base pairs is affected by elements of the structure such as the bulged-G motif and its cross-strand stacking. Participation in these motifs increases the lifetimes of the bases in an open, solvent-accessible conformation.     

Chemical consequences of incorporation of 5-fluorouracil into DNA as studied by NMR

The cytotoxic analogue of thymine, 5-fluorouracil (Uf), is known to be incorporated into DNA in biological systems. This abnormal base has been synthetically incorporated into short DNA oligomers. The ionization of the N-3 proton of this base within DNA oligomers was measured by observation of the 19F chemical shift at varying pH values. The pKa values for the Uf ring of dTpdUfpdT and dApdUfpdA were determined to be 7.84 and 7.9, respectively. The self-complementary 12-mers d(G-C-G-C-A-A-T-Uf-G-C-G-C) and d(C-G-A-T-Uf-A-T-A-A-T-C-G) were synthesized, and 1H NMR was used to compare the helix dynamics and stability of the interstrand imino proton hydrogen bonds with those of the 12-mers d(G-C-G-C-A-A-T-T-G-C-G-C) and d(C-G-A-T-T-A-T-A-A-T-C-G). The N-3 hydrogen bond of the A-Uf base pair was less stable than the corresponding hydrogen bond in A-T base pairs in the same helix, and the A-Uf base pair was less stable than the A-T base pair in the analogous position of the control helix. The observed temperature-dependent dynamics and NMR melting temperatures of the control and dUf-containing oligomers were similar.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Stability of XGCGCp, GCGCYp, and XGCGCYp helixes: an empirical estimate of the energetics of hydrogen bonds in nucleic acids

The stabilizing effects of dangling ends and terminal base pairs on the core helix GCGC are reported. Enthalpy and entropy changes of helix formation were measured spectrophotometrically for AGCGCU, UGCGCA, GGCGCCp, CGCGCGp, and the corresponding pentamers XGCGCp and GCGCYp containing the GCGC core plus a dangling end. Each 5' dangling end increases helix stability at 37 degrees C roughly 0.2 kcal/mol and each 3' end from 0.8 to 1.7 kcal/mol. The free energy increments for dangling ends on GCGC are similar to the corresponding increments reported for the GGCC core [Freier, S. M., Alkema, D., Sinclair, A., Neilson, T., & Turner, D. H. (1985) Biochemistry 24, 4533-4539], indicating a nearest-neighbor model is adequate for prediction of stabilization due to dangling ends. Nearest-neighbor parameters for prediction of the free energy effects of adding dangling ends and terminal base pairs next to G.C pairs are presented. Comparison of these free energy changes is used to partition the free energy of base pair formation into contributions of "stacking" and "pairing". If pairing contributions are due to hydrogen bonding, the results suggest stacking and hydrogen bonding make roughly comparable favorable contributions to the stability of a terminal base pair. The free energy increment associated with forming a hydrogen bond is estimated to be -1 kcal/mol of hydrogen bond.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.