Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Fluorescent Nanoparticles Synthesized from DNA, RNA, and Nucleotides

Ubiquitous on Earth, DNA and other nucleic acids are being increasingly considered as promising biomass resources. Due to their unique chemical structure, which is different from that of more common carbohydrate biomass polymers, materials based on nucleic acids may exhibit new, attractive characteristics. In this study, fluorescent nanoparticles (biodots) were prepared by a hydrothermal (HT) method from various nucleic acids (DNA, RNA, nucleotides, and nucleosides) to establish the relationship between the structure of precursors and fluorescent properties of biodots and to optimize conditions for preparation of the most fluorescent product. HT treatment of nucleic acids results in decomposition of sugar moieties and depurination/depyrimidation of nucleobases, while their consequent condensation and polymerization gives fluorescent nanoparticles. Fluorescent properties of DNA and RNA biodots are drastically different from biodots synthesized from individual nucleotides. In particular, biodots synthesized from purine-containing nucleotides or nucleosides show up to 50-fold higher fluorescence compared to analogous pyrimidine-derived biodots. The polymeric nature of a precursor disfavors formation of a bright fluorescent product. The reported effect of the structure of the nucleic acid precursor on the fluorescence properties of biodots should help designing and synthesizing brighter fluorescent nanomaterials with broader specification for bioimaging, sensing, and other applications.

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

The importance of hydration and DNA conformation in interpreting infrared spectra of cells and tissues

Since Watson and Crick's historical papers on the structure and function of DNA based on Rosalind Franklin's and Maurice Wilkin's X-ray diffraction patterns tremendous scientific curiosity has been aroused by the unique and dynamic structure of the molecule of life. A-DNA and B-DNA represent different conformations of the DNA molecule, which is stabilised by hydrogen interactions between base pairs, stacking interactions between neighboring bases and long-range intra- and inter-backbone forces. This review highlights the contribution Fourier transform infrared (FTIR) spectroscopy has made to the understanding of DNA conformation in relation to hydration and its potential role in clinical diagnostics. The review will first begin by elucidating the main forms of DNA conformation found in nature and the general structures of the A, B and Z forms. This is followed by a detailed critique on infrared spectroscopy applied to DNA conformation highlighting pivotal studies on isolated DNA, polynucleotides, nucleoprotein and nucleohistone complexes. A discussion on the potential of diagnosing cancer using FTIR spectroscopy based on the detection of DNA bands in cells and tissues will ensue, highlighting the recent studies investigating the conformation of DNA in hydrated and dehydrated cells. The method of hydration as a way to facilitate DNA conformational band assignment will be discussed and the conformational change to the A-form upon dehydration will be used to explain the reason for the apparent lack of FTIR DNA signals observed in fixed or air-dried cells and tissues. The advantages of investigating B-DNA in the hydrated state, as opposed to A-DNA in the dehydrated state, are exemplified in a series of studies that show: (1) improved quantification of DNA in cells; (2) improved discrimination and reproducibility of FTIR spectra recorded of cells progressing through the cell cycle; (3) insights into the biological significance of A-DNA as evidenced by an interesting study on bacteria, which can survive desiccation and at the same time undergo the B-A-B transition. Finally, the importance of preserving the B-DNA conformation for the diagnosis of cancer is put forward as way to improve the sensitivity of this powerful technique.

8-Oxoguanine Affects DNA Backbone Conformation in the EcoRI Recognition Site and Inhibits Its Cleavage by the Enzyme

8-oxoguanine is one of the most abundant and impactful oxidative DNA lesions. However, the reasons underlying its effects, especially those not directly explained by the altered base pairing ability, are poorly understood. We report the effect of the lesion on the action of EcoRI, a widely used restriction endonuclease. Introduction of 8-oxoguanine inside, or adjacent to, the GAATTC recognition site embedded within the Drew—Dickerson dodecamer sequence notably reduced the EcoRI activity. Solution NMR revealed that 8-oxoguanine in the DNA duplex causes substantial alterations in the sugar—phosphate backbone conformation, inducing a BI→BII transition. Moreover, molecular dynamics of the complex suggested that 8-oxoguanine, although does not disrupt the sequence-specific contacts formed by the enzyme with DNA, shifts the distribution of BI/BII backbone conformers. Based on our data, we propose that the disruption of enzymatic cleavage can be linked with the altered backbone conformation and dynamics in the free oxidized DNA substrate and, possibly, at the protein—DNA interface.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Structural basis of the RNase H1 activity on stereo regular borano phosphonate DNA/RNA hybrids

Numerous DNA chemistries for improving oligodeoxynucleotide (ODN)-based RNA targeting have been explored. The majority of the modifications render the ODN/RNA target insensitive to RNase H1. Borano phosphonate ODN's are among the few modifications that are tolerated by RNase H1. To understand the effect of the stereochemistry of the BH(3) modification on the nucleic acid structure and RNase H1 enzyme activity, we have investigated two DNA/RNA hybrids containing either a R(P) or S(P) BH(3) modification by nuclear magnetic resonance (NMR) spectroscopy. T(M) studies show that the stabilities of R(P) and S(P) modified DNA/RNA hybrids are essentially identical (313.8 K) and similar to that of an unmodified control (312.9 K). The similarity is also reflected in the imino proton spectra. To characterize such similar structures, we used a large number of NMR restraints (including dipolar couplings and backbone torsion angles) to determine structural features that were important for RNase H1 activity. The final NMR structures exhibit excellent agreement with the data (total R(x) values of <6%) with helical properties between those of an A and B helix. Subtle backbone variations are observed in the DNA near the modification, while the RNA strands are relatively unperturbed. In the case of the S(P) modification, for which more perturbations are recorded, a slightly narrower minor groove is also obtained. Unique NOE base contacts localize the S(P) BH(3) group in the major groove while the R(P) BH(3) group points away from the DNA. However, this creates a potential clash of the R(P) BH(3) groups with important RNase H1 residues in a complex, while the S(P) BH(3) groups could be tolerated. We therefore predict that on the basis of our NMR structures a fully R(P) BH(3) DNA/RNA hybrid would not be a substrate for RNase H1.

Structural analysis of the binding of the diquaternary pyridophenazine derivative dqdppn to B-DNA oligonucleotides

The interaction of the ethylene-bipyridyldiylium-naphthaphenazine dication, dqdppn, with several hexa- and octanucleotide duplexes has been studied using CD and NMR. Taken together, these studies reveal that with the hexanucleotide, dqdppn intercalates into the terminal base pair, and causes a large twisting of the terminal base pair. In contrast, with all three octanucleotides, dqdppn intercalates more centrally within the sequence. The NMR-derived structures of two of the binding complexes demonstrate that dqdppn intercalates from the major groove in an unusual 'side-on' geometry, rather than threading through the helix. An analysis of these results indicates that the preferred binding site is not sequence-specific, but primarily at the most conformationally flexible DNA step.

Nucleic acid folding determined by mesoscale modeling and NMR spectroscopy: solution structure of d(GCGAAAGC)

Determination of DNA solution structure is a difficult task even with the high-sensitivity method used here based on simulated annealing with 35 restraints/residue (Cryoprobe 750 MHz NMR). The conformations of both the phosphodiester linkages and the dinucleotide segment encompassing the sharp turn in single-stranded DNA are often underdetermined. To obtain higher quality structures of a DNA GNRA loop, 5'-d(GCGAAAGC)-3', we have used a mesoscopic molecular modeling approach, called Biopolymer Chain Elasticity (BCE), to provide reference conformations. By construction, these models are the least deformed hairpin loop conformation derived from canonical B-DNA at the nucleotide level. We have further explored this molecular conformation at the torsion angle level with AMBER molecular mechanics using different possible (epsilon,zeta) constraints to interpret the 31P NMR data. This combined approach yields a more accurate molecular conformation, compatible with all the NMR data, than each method taken separately, NMR/DYANA or BCE/AMBER. In agreement with the principle of minimal deformation of the backbone, the hairpin motif is stabilized by maximal base-stacking interactions on both the 5'- and 3'-sides and by a sheared G.A mismatch base pair between the first and last loop nucleotides. The sharp turn is located between the third and fourth loop nucleotides, and only two torsion angles beta6 and gamma6 deviate strongly with respect to canonical B-DNA structure. Two other torsion angle pairs epsilon3,zeta3 and epsilon5,zeta5 exhibit the newly recognized stable conformation BIIzeta+ (-70 degrees, 140 degrees). This combined approach has proven to be useful for the interpretation of an unusual 31P chemical shift in the 5'-d(GCGAAAGC)-3' hairpin.

NMR methods for studying quadruplex nucleic acids

Solution NMR spectroscopy has traditionally played a central role in examining quadruplex structure, dynamics, and interactions. Here, an overview is given of the methods currently applied to structural, dynamics, thermodynamics, and kinetics studies of nucleic acid quadruplexes and associated cations.

NMR studies of the conformational effect of single and double 3'-S-phosphorothiolate substitutions within deoxythymidine trinucleotides

NMR spectroscopy has been used to investigate the conformational effects of single and two consecutive 3'-S-phosphorothiolate modifications within a deoxythymidine trinucleotide. The presence of a single 3'-phosphorothioate modification shifts the conformation of the sugar ring it is attached to, from a mainly south to north pucker; this effect is also transmitted to the 3'-neighbour deoxyribose. This transmission is thought to be caused by favourable stacking of the heterocyclic bases. Similar observations have been made previously by this group. When two adjacent modifications are present, the conformations of the attached deoxyribose rings are again shifted almost completely to the north, however, there is no transmission to the 3' deoxyribose ring. Base proton chemical shift analysis and molecular modelling have been used to aid elucidation of the origin of this feature. The observation for the dimodified sequence is consistent with our previously reported results for a related system in which spaced modifications are more thermodynamically stable than consecutive ones.

Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.

Mispairing of a site specific major groove (2S,3S)-N6-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl DNA Adduct of butadiene diol epoxide with deoxyguanosine: formation of a dA(anti).dG(anti) pairing interaction

The (2S,3S)-N6-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl (BDT) adduct arising from alkylation of adenine N6 by butadiene diol epoxide (BDE) was placed opposite a mismatched deoxyguanosine nucleotide in the complementary strand of the oligodeoxynucleotide 5'-d(CGGACXAGAAG)-3'.5'-d(CTTCTGGTCCG)-3'. This oligodeoxynucleotide contains codon 61 (underlined) of the human N-ras protooncogene. The BDT adduct was at the second position of codon 61, and this was named the ras61 S,S-BDT-(61,2) A.G adduct. NMR spectroscopy revealed the presence of two conformations of the adducted mismatched duplex. In the major conformation, the mismatched base pair X6.G17 was oriented in a "face-to-face" orientation, in which both the modified nucleotide X6 and its complement G17 were intrahelical and in the anti conformation about the glycosyl bond. Hydrogen bonding was suggested between X6 N1 and G17 N1H and between X6 N6H and G17 O6. The presence of the BDT moiety allowed formation of a stable A.G mismatch pair. The identity of the minor conformation could not be determined. If not repaired, the resulting mismatch pair would generate A-->C mutations, which have been associated with this adenine N6 BDT adduct [Carmical, J. R., Nechev, L. N., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2000) Env. Mol. Mutagen. 35, 48-56].

Structure of a site specific major groove (2S,3S)-N6-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl DNA adduct of butadiene diol epoxide

The solution structure of the (2S,3S)-N(6)-(2,3,4-trihydroxybutyl)-2'-deoxyadenosyl adduct arising from the alkylation of adenine N(6) at position X(6) in d(CGGACXAGAAG).d(CTTCTTGTCCG), by butadiene diol epoxide, was determined. This oligodeoxynucleotide contains codon 61 (underlined) of the human N-ras protooncogene. This oligodeoxynucleotide, containing the adenine N(6) adduct butadiene triol (BDT) adduct at the second position of codon 61, was named the ras61 S,S-BDT-(61,2) adduct. NMR spectroscopy revealed modest structural perturbations localized to the site of adduction at X(6).T(17), and its nearest-neighbor base pairs C(5).G(18) and A(7).T(16). All sequential NOE connectivities arising from DNA protons were observed. Torsion angle analysis from COSY data suggested that the deoxyribose sugar at X(6) remained in the C2'-endo conformation. Molecular dynamics calculations using a simulated annealing protocol restrained by a total of 442 NOE-derived distances and J coupling-derived torsion angles refined structures in which the BDT moiety oriented in the major groove. Relaxation matrix analysis suggested hydrogen bonding between the hydroxyl group located at the beta-carbon of the BDT moiety and the T(17) O(4) of the modified base pair X(6).T(17). The minimal perturbation of DNA induced by this major groove adduct correlated with its facile bypass by three Escherichia coli DNA polymerases in vitro and its weak mutagenicity [Carmical, J. R., Nechev, L. V., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2000) Environ. Mol. Mutagen. 35, 48-56]. Overall, the structure of this adduct is consistent with an emerging pattern in which major groove adenine N(6) alkylation products of styrene and butadiene oxides that do not strongly perturb DNA structure are not strongly mutagenic.

Random coil phosphorus chemical shift of deoxyribonucleic acids

Random coil phosphorus chemical shift has been studied using 16 17-nucleotide DNA sequences. Due to the presence of residual base stacking in these sequences, the temperature and sequence effects were investigated at 50 and 55 degrees C. The phosphorus chemical shifts of random coil DNA sequences have been found to be independent of temperature. Sequence effect analysis shows that the phosphorus chemical shift of a nucleotide in a random coil DNA sequence depends on both its 5'- and 3'-nearest neighbors. A trimer model has been used to establish the random coil 31P chemical shift prediction protocol which shows an accuracy of 0.02 ppm.

Discrimination between BI and BII conformational substates of B-DNA based on sugar-base interproton distances

Molecular dynamics (MD) simulations of four water-solvated DNA duplexes were used to generate a database of approximately 27000 dinucleotide conformations. Analyzing this database, we investigated the relationship between so-called BI-BII transitions and short-range interproton distances. Four H-H distances were found particularly sensitive to BI-BII transitions: internucleotide H1'(n)-H68(n+1), H2'(n)-H68(n+1), and H2"(n)-H68(n+1), and intranucleotide H2"(n)-H68(n). Determination of these distances using classical NOESY spectroscopy can thus provide valuable indications on the existence of BII substates, complementing the existing method based on (31)P chemical shifts and (31)P-(1)H spin-spin coupling constants.

Sequence context effect on the structure of nitrous acid induced DNA interstrand cross-links

In the preceding paper in this journal, we described the solution structure of the nitrous acid cross-linked dodecamer duplex [d(GCATCCGGATGC)]2 (the cross-linked guanines are underlined). The structure revealed that the cross-linked guanines form a nearly planar covalently linked 'G:G base pair', with the complementary partner cytidines flipped out of the helix. Here we explore the flanking sequence context effect on the structure of nitrous acid cross-links in [d(CG)]2 and the factors allowing the extrahelical cytidines to adopt such fixed positions in the minor groove. We have used NMR spectroscopy to determine the solution structure of a second cross-linked dodecamer duplex, [d(CGCTACGTAGCG)]2, which shows that the identity of the flanking base pairs significantly alters the stacking patterns and phosphate backbone conformations. The cross-linked guanines are now stacked well on adenines preceding the extrahelical cytidines, illustrating the importance of purine- purine base stacking. Observation of an imino proton resonance at 15.6 p.p.m. provides evidence for hydrogen bonding between the two cross-linked guanines. Preliminary structural studies on the cross-linked duplex [d(CGCGACGTCGCG)]2 show that the extrahelical cytidines are very mobile in this sequence context. We suggest that favorable van der Waals interactions between the cytidine and the adenine 2 bp away from the cross-link localize the cytidines in the previous cross-linked structures.

Solution structure of a nitrous acid induced DNA interstrand cross-link

Nitrous acid is a mutagenic agent. It can induce interstrand cross-links in duplex DNA, preferentially at d(CpG) steps: two guanines on opposite strands are linked via a single shared exocyclic imino group. Recent synthetic advances have led to the production of large quantities of such structurally homogenous cross-linked duplex DNA. Here we present the high resolution solution structure of the cross-linked dodecamer [d(GCATCCGGATGC)]2 (the cross-linked guanines are underlined), determined by 2D NMR spectroscopy, distance geometry, restrained molecular dynamics and iterative NOE refinement. The cross-linked guanines form a nearly planar covalently linked 'G:G base pair' with only minor propeller twisting, while the cytidine bases of their normal base pairing partners have been flipped out of the helix and adopt well defined extrahelical positions in the minor groove. On the 5'-side of the cross-link, the minor groove is widened to accommodate these extrahelical bases, and the major groove becomes quite narrow at the cross-link. The cross-linked 'G:G base pair' is well stacked on the spatially adjacent C:G base pairs, particularly on the 3'-side guanines. In addition to providing the first structure of a nitrous acid cross-link in DNA, these studies could be of major importance to the understanding of the mechanisms of nitrous acid cross-linking and mutagenicity, as well as the mechanisms responsible for its repair in intracellular environments. It is also the shortest DNA cross-link structure to be described.

Solution structure of a DNA duplex containing 8-hydroxy-2'-deoxyguanosine opposite deoxyguanosine

Deoxyguanosine residues are hydroxylated by reactive oxygen species at the C-8 position to form 8-hydroxy-2'-deoxyguanosine (8-OG), one of the most important mutagenic lesions in DNA. Though the spontaneous G:C to C:G transversions are rare events, the pathways leading to this mutation are not established. An 8-OG:G mispair, if not corrected by DNA repair enzymes, could lead to G:C to C:G transversions. NMR spectroscopy and restrained molecular dynamics calculations are used to refine the solution structure of the base mismatch formed by the 8-OG:G pair on a self complementary DNA dodecamer duplex d(CGCGAATT(8-O)GGCG)(2). The results reveal that the 8-OG base is inserted into the helix and forms Hoogsteen base-pairing with the G on the opposite strand. The 8-OG:G base-pairs are seen to be stabilized by two hydrogen bonding interactions, one between the H7 of the 8-OG and the O6 of the G, and a three-center hydrogen bonding between the O8 of the 8-OG and the imino and amino protons of the G. The 8-OG:G base-pairs are very well stacked between the Watson-Crick base-paired flanking bases. Both strands of the DNA duplex adopt right-handed conformations. All of the unmodified bases, including the G at the lesion site, adopt anti glycosidic torsion angles and form Watson-Crick base-pairs. At the lesion site, the 8-OG residues adopt syn conformations. The structural studies demonstrate that 8-OG(syn):G(anti) forms a stable pair in the interior of the duplex, providing a basis for the in vivo incorporation of G opposite 8-OG. Calculated helical parameters and backbone torsional angles, and the observed 31P chemical shifts, indicate that the structure of the duplex is perturbed near lesion sites, with the local unwinding of the double helix. The melting temperature of the 8-OG:G containing duplex is only 2.6 deg. C less than the t(m) of the unmodified duplex.

Synthesis, characterization and solution structure of tethered oligonucleotides containing an internal 3'-phosphoglycolate, 5'-phosphate gapped lesion

Bleomycins (BLMs) are antitumor antibiotics that in the presence of iron and oxygen mediate DNA damage by 4'-hydrogen atom abstraction of pyrimidines 3' to guanines. The resulting 4'-deoxyribose radicals can be trapped by O2 and ultimately result in the formation of base-propenal and gapped DNA with 3'-phosphoglycolate (3'-PG) and 5'-phosphate (5'-P) ends. The role of this lesion in triggering double-strand cleavage of duplex DNA by a single BLM molecule and the mechanism by which this lesion is repaired in vivo remain unsolved problems. The structure of these lesions is an essential step in addressing both of these problems. Duplex DNAs (13mers containing tethered hexaethylene glycol linkers) with GTAC and GGCC cleavage sites have been synthesized in which gaps containing 3'-PG and 5'-P ends at the sites of BLM cleavage have been inserted. The former sequence represents a hot spot for double-strand cleavage, while the latter is a hot spot for single-strand cleavage. Analytical methods to characterize the lesioned products have been developed. These oligonucleotides have been examined using 2D NMR methods and molecular modeling. The studies reveal that the lesioned DNAs are B-form and the 3'-PG and 5'-P are extrahelical. The base opposite the gap and the base pairs adjacent to the gap remain well stacked in the DNA duplex. Titrations of the lesioned GGCC oligomer with HOO-CoBLM leads to a mixture of complexes, in contrast to results of a similar titration with the lesioned GTAC oligomer.

The role of the phosphorus BI-BII transition in protein-DNA recognition: the NF-kappaB complex

We examined, by 1H and 31P NMR, the solution structure of a 16 bp non-palindromic DNA fragment (16M2) containing the HIV-1 NF-kappaB-binding site, in which the sequences flanking the kappaB site had been mutated. 31P NMR was particularly useful for obtaining structural information on the phosphodiester backbone conformation. Structural features were then compared with those of the two previously studied DNA fragments corresponding, respectively, to the native kappaB fragment (16N) and a fragment in which mutations have been introduced at the 5' end of the kappaB site (16M1). For the mutated 16M2 duplex, NMR data showed that the BI-BII equilibrium, previously reported for the native fragment (16N) at the kappaB flanking steps, was lost. The role of the BI-BII equilibrium in NF-kappaB recognition by DNA was then investigated by electrophoretic mobility shift assay. We found that the isolated kappaB site has the potential to bind efficiently due to the BI-BII equilibrium of the kappaB flanking sequences.

NMR analysis of helix I from the 5S RNA of Escherichia coli

The structure of helix I of the 5S rRNA from Escherichia coli has been determined using a nucleolytic digest fragment of the intact molecule. The fragment analyzed, which corresponds to bases (-1)-11 and 108-120 of intact 5S rRNA, contains a G-U pair and has unpaired bases at its termini. Its proton resonances were assigned by two-dimensional NMR methods, and both NOE distance and coupling constant information have been used to calculate structural models for it using the full relaxation matrix algorithm of the molecular dynamics program XPLOR. Helix I has A-type helical geometry, as expected. Its most striking departure from regular helical geometry occurs at its G-U, which stacks on the base pair to the 5' side of its G but not on the base pair to its 3' side. This stacking pattern maximizes interstrand guanine-guanine interactions and explains why the G-U in question fails to give imino proton NOE's to the base pair to 5' side of its G. These results are consistent with the crystal structures that have been obtained for wobble base pairs in tRNAPhe [Mizuno, H., & Sundaralingam, M. (1978) Nucleic Acids Res. 5, 4451-4461] and A-form DNA [Rabbinovich, D., Haran, T., Eisenstein, M., & Shakked, Z. (1988) J. Mol. Biol. 200, 151-161]. The conformations of the terminal residues of helix I, which corresponds to bases (-1)-11 and 108-120 of native 5S RNA, are less well-determined, and their sugar puckers are intermediate between C2' and C3'-endo, on average.

Quadruple intercalated G-6 stack: a possible motif in the fold-back structure of the Drosophila centromeric dodeca-satellite?

The purine-rich strand d(GTACGGGACCGA)(n) of the Drosophila centromeric dodeca-satellite sequence is highly conserved and was found to form stable fold-back structures in which the homopurine 5'-GGGA-3' sequence was determined to play a crucial role. Here, we report the stable formation of the d(GGGA)(2) motif in the stem of a DNA hairpin closed by a single-residue d(ACC) loop. Similar to the zipper-like d(GGA)(2) motif observed in the human centromeric (TGGAA)(n) sequence, the central four guanosine bases in the d(GGGA)(2) motif do not pair, but interdigitate to form an elongated zipper-like quadruple-intercalated G-6 stack bracketed by sheared G.A base-pairs. Comparison between the current d(GGGA)(2) structure and the published crystal d(GAAA)(2) structure implies that the alignment of the unpaired purine bases plays an important role in determining the minor groove width of the purine-rich d(GPuPuA)(2) motif. Similarity between the zipper-like motifs possibly present in the Drosophila centromeric dodeca-satellite sequence and in the human centromeric (TGGAA)(n) sequence led us to propose that these special zipper-like motifs may constitute common cores in organizing eukaryotic centromeres.

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.

Solution conformation of an RNA--DNA hybrid duplex containing a pyrimidine RNA strand and a purine DNA strand

RNA--DNA hybrid duplexes are involved in transcription, replication and reverse transcription of nucleic acids. Information on such duplexes may shed some light on the mechanism of these processes. For this purpose, the influence of base composition on the structure of a polypyrimidine--polypurine RNA--DNA duplex r(cucuccuucucuu). d(GAGAGGAAGAGAA) has been studied using 1H, 31P and 13C NMR experiments, molecular modeling (JUMNA program) and NOE back-calculation methods. The resulting structure of the 13-mer hybrid duplex shows that the RNA strand is in the expected A-type conformation while the DNA strand is in a very flexible conformation. In the DNA strand, the desoxyribose sugars retain the C2'-endo B-type conformation. The duplex helical parameters (such as inclination, twist and displacement of the bases) are close to the A-type conformation. No bending was observed for the global axis curvature. The major groove width is close to the B-form value and the minor groove width is intermediate between standard values for A and B-forms. These results are in favour of the independence of minor groove size (where RNase H interacts) and the base composition of the hybrid duplexes.

BII nucleotides in the B and C forms of natural-sequence polymeric DNA: A new model for the C form of DNA

A combination of solid-state (31)P and (13)C NMR, X-ray diffraction, and model building is used to show that the B and C forms of fibrous macromolecular DNA consist of two distinct nucleotide conformations, which correspond closely to the BI and BII nucleotide conformations known from oligonucleotide crystals. The proportion of the BII conformation is higher in the C form than in the B form. We show structural models for a 10(1) double helix involving BI nucleotides and a 9(1) double helix involving BII nucleotides. The 10(1) BI model is similar to a previous model of B-form DNA, while the 9(1) BII model is novel. The BII model has a very deep and narrow minor groove, a shallow and wide major groove, and highly inclined bases. This work shows that the B to C transition in fibers corresponds to BI to BII conformational changes of the individual nucleotides.

Solution structure of the orphan nuclear receptor rev-erb beta response element by 1H, 31P NMR and molecular simulation*

Rev-erb beta is an orphan receptor that binds as a homodimer or as a monomer to DNA. The solution structure of the non-palindromic 15 bp DNA duplex d(TAGAATGTAGGTCAG), the response element of Rev-erb beta for monomeric binding, was determined by 1H and 31P NMR, energy minimization with NMR-derived restraints for distances and NOE back-calculation methods. The refined final structures have the typical overall features of B-type DNA. However, titration of this 15 bp duplex with ReDBD, the DNA binding domain of Rev-erb beta, showed large shifts of imino protons and 31P signals, suggesting major conformational changes.

Cross-strand purine-pyrimidine stack and sheared purine.pyrimidine pairing in the human HIV-1 reverse transcriptase inhibitors

Cross-strand homo purine-purine (G-G or A-A) stacks and sheared purine.purine pairing have been found to be important motifs in nucleic acid duplex structures. We now report novel cross-strand purine-pyrimidine (A-C) and hetero purine-purine (G-A) stacks that are established from a sheared purine.pyrimidine (A.C) pair adjacent to a sheared G.A pair in the 5'-AA/GC-3' sequence. This "internal loop" sequence is conserved in two families of single-stranded DNA inhibitors of the reverse transcriptase of type 1 human immunodeficiency virus. The distorted backbone of these inhibitors, resulting from the unique helical twists and kinks in the 5'-AA/GC-3' sequence, may be responsible for the increased affinities of these single-stranded DNA inhibitors as compared with other regular B-form duplex substrates. Two simple rules have been generalized to account for all reported cross-strand stacks.

Structure determination and analysis of helix parameters in the DNA decamer d(CATGGCCATG)2 comparison of results from NMR and crystallography

The solution structure of the DNA decamer (CATGGCCATG)2 has been determined by NMR spectroscopy and restrained molecular dynamic and distance geometry calculations. The restrainted data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were obtained by nuclear Overhauser enhancement intensities and quantitative simulation of cross-peaks from double quantum filtered correlation spectroscopy. The backbone torsion angles were constrained using experimental data from NOE cross-peaks, 1H-1H and 1H-31P-coupling constants. The NMR structure and the crystal structure of the DNA decamer deviates from the structure of the canonical form of B-DNA in a number of observable characteristics. Particularly, both structures display a specific pattern of stacking interaction in the central GGC base triplet. Furthermore, a specific local conformation of the TG/CA base-pair step is present in NMR and crystal structure, highlighting the unusually high flexibility of this DNA duplex part. The solution structure of the TG/CA base-pair step obtained by our high resolution NMR study is characterized by a positive roll angle, whereas in crystal this base-pair step tends to adopt remarkably high twist angles.

31P NMR study of the interactions between oligodeoxynucleotides containing (6-4) photoproduct and Fab fragments of monoclonal antibodies specific for (6-4) photoproduct

A 31P nuclear magnetic resonance (NMR) study of the interactions between oligonucleotides containing the (6-4) photoproduct and the Fab fragments of monoclonal antibodies (64M3 and 64M5) recognizing the (6-4) photoproduct is reported. The 31P chemical shift data indicate that backbone conformation of (64) adduct is affected by the presence of flanking oligodeoxynucleotides, and (6-4) adducts with different backbone conformations are accommodated in the antigen binding sites of these antibodies. It was also revealed that epitopes for these antibodies consist of not only the (6-4) adduct but the flanking di- or tri-deoxynucleotides on both the 5' and 3' sides as well.

Centromeric pyrimidine strands fold into an intercalated motif by forming a double hairpin with a novel T:G:G:T tetrad: solution structure of the d(TCCCGTTTCCA) dimer

The solution structures of the oligodeoxynucleotides d(CCCGTTTCC) and d(TCCCGTTTCCA) have been determined by two-dimensional NMR spectroscopy. These oligomers are part of a DNA box in human centromeric alpha satellite targeted by the centromere protein B (CENP-B). Both CENP-B and its recognition box in alphoid DNA are conserved in mammals, suggesting an important biological role. At acidic pH, d(CCCGTTTCC), d(TCCCGTTTCCA) and the full d(TCCCGTTTCCAACGAAG) CENP-B box strand all fold and dimerize in solution forming a stable bimolecular structure containing two GTTT hairpin loops that interact through a novel T : G : G : T tetrad. The stem region of the dimer is a four-stranded intercalated motif in which the hairpin monomers are parallel and held together by C : C+ hydrogen-bonding and intercalation. The loops are at the same end of the dimer and lie across the narrow grooves of the tetraplex. They are remarkably structured and stabilized by base-base cross-stacking, sugar-base stacking, and parallel G:G and antiparallel G:T pairing. In the d(TCCCGTTTCCA)2 structure, the intercalated motif is continued at the other end of the dimer with unpaired but stacked adenine and thymine bases. The possible biological implications of these structures are discussed.

Sheared purine x purine pairing in biology

The Watson-Crick G x C and A x T base-paired DNA duplex has been the single most important milestone in modem molecular biology. However, it is possible that other types of stable DNA structures besides the double helix might exist, since only about 5% of the human chromosome is transcribed and expressed. Stable, four-stranded G-tetraplex DNA structures occur in the extensive tandem repeated sequences at the telomeres of chromosome. Formation of stable triplexes of the Py x Pu x Py or Pu x Pu x Py type have been implicated at the control regions of certain human genes. We review and discuss the various types of DNA duplex structures containing stable sheared base-pairs and compare their structural characteristics with that of B-DNA. Pu x Pu structural motifs are found in the highly conserved sequences at the replication origins of several single-stranded DNA viruses and in the peri-centromeric regions of human chromosomes, and may be involved in important biological functions, such as viral DNA replication and centromere formation.

Sequence-specific local structural variations in solution structures of d(CGXX'CG)2 and d(CAXX'TG)2 self-complementary deoxyribonucleic acids

In this study, the solution structures of eight self-complementary deoxyribonucleic acid molecules, d(CGXX'CG)2 and d(CAXX'TG)2 (where X = C, G, T or A and X' is complementary to X), have been determined using NMR interproton distances, endocyclic sugar torsion angles, backbone torsion angles and hydrogen bond constraints. Based on the structural data obtained in solution, a novel sequence-specific local structure function, sigma LS, composed of the sum of the contributions from the helix twist omega, base roll rho, base-pair slide delta delta and propeller twist omega, is introduced to describe their sequence-specific local structures. Sigma LS is found to produce an acceptable correlation (r = 0.96) with the relative local stability (delta G0 10) of a base-pair-step. This result demonstrates that in addition to the inter-strand purine-purine clashes, the base morphology of nearest-neighbor base-pairs is also important in defining the local geometry of base-pairs. Thus, in analyzing the base-pair structural parameter blocks of trimers are used as the basic unit whereas for the base-pair-step structural parameters, the basic unit is composed of blocks of tetramers. The omega LS-delta G(zero) 10 correlation is the first experimental evidence demonstrating the relationship between the relative local stabilities and the sequence-specific local structures of DNA duplexes. It also forms the basis for using the trimer-tetramer model for the prediction of sequence-specific local structures of deoxyribonucleic acid molecules.

Hairpin loops consisting of single adenine residues closed by sheared A.A and G.G pairs formed by the DNA triplets AAA and GAG: solution structure of the d(GTACAAAGTAC) hairpin

The DNA undecamers GTACAAAGTAC (AAA 11-mer) and GTACGAGGTAC (GAG 11-mer) have been studied in solution by high-resolution NMR spectroscopy. Both duplexes form stable hairpins containing single deoxyadenosine loops and stems containing five base-pairs that are closed at the loop end by sheared AxA and GxC pairs, respectively. These molecules thus contain new AAA and GAG loop turn motifs. All protons, including the chiral H5'/H5" protons of the loop residues, were assigned using NOESY, DQF-COSY and heteronuclear 1H-31P COSY experiments. The backbone torsion angles were constrained using experimental data from NOE crosspeaks, three-bond 1H-1H coupling constants and four-bond 1H-31P coupling constants and four-bond 1H-31P coupling constants. The AAA and GAG 11-mers form similar structures in solution. The detailed structure of the AAA 11-mer was determined by the combined use of NMR, distance geometry and energy minimization methods. This structure exhibits good stacking of the loop adenosine base on the closing 5Ax7A sheared pair, with the 6A base stacking on the 5A base and the 6A deoxyribose stacking with the 7A base. All sugars in the AAA 11-mer hairpin adopt the typical DNA C2'-endo conformation and a sharp backbone turn occurs between residues 6A and 7A. This loop turn is brought about mainly by a change in the backbone phosphate torsion angles from zeta(g-) alpha(g-) to zeta(g+) alphat(g+) at the turn. The gamma torsion angle of residue 7A in the closing sheared pair also changes from gauche+ to trans. In Pu1NPu2 loop turns of the GCA, AAA and GAG types, the chemical shift of the H4' proton of the loop deoxyribose depends on the nature of Pu2; this reflects the stacking of the loop sugar on the Pu2 base and the different ring current effects of A or G in this position.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.

Conformational states governing the rates of spontaneous transition mutations

We report the results of a theoretical study, combining the results of sequence analysis and integral equation structural methods for nucleic acids in aqueous solutions, on the effects of nearest neighbors on the (T.G) mispair in solution, for 12 nearest neighbor contexts. Attempts have been made to classify the structural and energetic effects of the 5' and 3' neighbors with respect to the observed spontaneous mutation rates in vertebrates. It is found that 5' nearest neighbor is probably the most critical structural factor in facilitating or discouraging mutations. Local conformational states correlate with discrimination of bases to be excised in mispairs. Our study confirms the role of the flexibility of the DNA molecule in governing the rates of spontaneous mutations.

Sequence dependent effects of CpG cytosine methylation. A joint 1H-NMR and 31P-NMR study

The impact of cytosine methylation in the central CpG step of two closely related octanucleotide duplexes d(CATCGATG)2 and d(CTTCGAAG)2 was examined by 1H-NMR and 31P-NMR experiments, and a quantitative structural analysis was performed using the NOE-derived distances, the sugar puckers and the epsilon torsion angles. The two starting oligonucleotides displayed a B-DNA conformation with, however, significant local structure differences. Although the methylated oligonucleotides retained their B-DNA conformation, different structural and thermal stability effects were observed. The magnitude of the methylation effects was to depend on the initial conformation of the CpG site, which is governed by the nature of the dinucleotide AT or TT located on the CpG flanks. As an example of sequence dependence, the methylation of CpG entailed larger conformational variation in d(CATCGATG)2 than in d(CTTCGAAG)2. In this study, the 1H and 31P chemical-shift parameters averred as extremely sensitive probes for detecting subtle conformational changes. Finally, our comparative results may aid our understanding of the structural and related biological effects produced by cytosine methylation in DNA.

Modifying the helical structure of DNA by design: recruitment of an architecture-specific protein to an enforced DNA bend

BACKGROUND

Proteins can force DNA to adopt distorted helical structures that are rarely if ever observed in naked DNA. The ability to synthesize DNA that contains defined helical aberrations would offer a new avenue for exploring the structural and energetic plasticity of DNA. Here we report a strategy for the enforcement of non-canonical helical structures through disulfide cross-linking; this approach is exemplified by the design and synthesis of an oligonucleotide containing a pronounced bend.

RESULTS

A localized bend was site-specifically introduced into DNA by the formation of a disulfide cross-link between the 5' adenines of a 5'-AATT-3' region in complementary strands of DNA. The DNA bend was characterized by high-resolution NMR structure determination of a cross-linked dodecamer and electrophoretic mobility assays on phased multimers, which together indicate that the cross-linked tetranucleotide induces a helical bend of approximately 30 degrees and a modest degree of unwinding. The enforced bend was found to stimulate dramatically the binding of an architecture-specific protein, HMG-D, to the DNA. DNase I foot-printing analysis revealed that the protein is recruited to the section of DNA that is bent.

CONCLUSIONS

The present study reports a novel approach for the investigation of non-canonical DNA structures and their recognition by architecture-specific proteins. The mode of DNA bending induced by disulfide cross-linking resembles that observed in structures of protein-DNA complexes. The results reveal common elements in the DNA-binding mode employed by sequence-specific and architecture-specific HMG proteins.

Molecular dynamics simulations of chlorambucil/DNA adducts. A structural basis for the 5'-GNC interstrand DNA crosslink formed by nitrogen mustards

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG3G2CGC).-d(GCG1CCCG) duplex, in which the N7 atoms of G1, G2 and G3 are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5'-GNC site (a 1,3 crosslink, G1-G3) rather than at a 5'-GC sites (a 1,2 crosslink, G1-G2). In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G1-G3 diadduct formation. It proved possible to construct a G1-G3 diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four BII conformations on the 3' side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G1. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G3. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.