Abasic frameshift in DNA. Solution conformation determined by proton NMR and molecular mechanics calculations

We have determined the three-dimensional structure of a non-self-complementary oligodeoxynucleotide duplex that contains a model abasic site. The duplex contains six GC base pairs plus the abasic site at the center of one strand and corresponds to an abasic frameshift. Two-dimensional NMR studies on the nonexchangeable protons show that the guanine bases on either side of the abasic site are stacked over each other and that the abasic site is rotated out of the helix. Close proton-proton interactions are observed between the H4' proton of the abasic site and sugar protons of the guanosine in the 5' direction, which allows the position of the free sugar to be well-defined. NOE buildup curves from NOESY spectra recorded at very short mixing times were used to calculate a set of interproton distances. This data set was incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. Two conformations that differ in the sugar conformation of the guanosine next to the abasic site in the 3' direction were necessary to fit all the NMR data. One of these two conformations could only be stabilized by addition of counterions at specific sites.     

Extrahelical adenosine stacks into right-handed DNA: solution conformation of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

This paper reports on features of the three-dimensional structure of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) self-complementary duplex (designated adenosine 13-mer), which contains symmetrical extrahelical adenosines in the interior of the helix. The majority of the protons have been assigned from two-dimensional nuclear Overhauser effect (NOESY) spectra of the adenosine 13-mer in H2O and D2O solution. The measurement of NOESY cross-peak volume integrals as a function of mixing time has yielded a set of 96 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds, which have served as input parameters for a distance geometry analysis of one symmetric half of the adenosine 13-mer duplex. We demonstrate that the extrahelical adenosine stacks into the duplex for all refined structures without disruption of base pairing on either side of the modification site. The distance geometry refinement yields two classes of conformations consistent with distance measurements but which differ in orientation of the stacked extrahelical adenosine at the modification site.     

Conformation of a bulge-containing oligomer from a hot-spot sequence by NMR and energy minimization

Two-dimensional nmr data on a bulge-containing oligodeoxyribonucleotide, 5'dGATGGGCAG.dCTGACCCATC, and a regular oligomer of similar sequence, 5'dGATGGCAG.dCTGCCATC, are presented. The nonexchangeable protons are assigned from sequential nuclear Overhauser effect spectroscopy (NOESY) connectivities. The two-dimensional NOE (NOESY) and correlated (COSY) spectra of the bulge-containing oligomer are compared to those of the perfect 8-mer. Experimental proton-proton distances are determined from NOESY spectra acquired with mixing times of 100, 150, and 200 ms, using comparable distances in the B-DNA region of the molecule as a calibration. With this approach, measured distances do not depend systematically on mixing time. Energy minimization techniques are used to calculate a three-dimensional structure for the bulge-containing oligomer in agreement with the nmr data. The helix is of the B family, with the extra adenine stacked into the helix, and the helix axis is bent by 20 degrees.     

Sequence-specific solution structure of d-GGTACGCGTACC

Complete resonance assignments of nonexchangeable base protons and sugar protons in d-GGTACGCGTACC at 500 MHz have been obtained by two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY). The characteristic phase-sensitive multiplet patterns of the ntrasugar cross peaks in the omega 1-scaled COSY spectrum have been used to estimate several scalar coupling constants (J). These coupling constants combined with the intranucleotide COSY cross peak intensities have been used to identify the sugar pucker of individual nucleotide units. In most cases, the deoxyribose rings adopt a conformation close to O4'-endo. Spin-diffusion has been monitored from the buildup of the normalized volumes of NOE cross peaks in NOESY spectra as a function of mixing time. A set of 55 intranucleotide and internucleotide interproton distances have been estimated from the low mixing time NOESY spectrum (tau m = 75 ms). The estimated intranucleotide proton-proton distances have been used to determine the individual glycosidic dihedral angles of the nucleotide units which lie in the anti domain. It is observed that the molecule adopts an overall conformation close to that of the B-form although there are differences in the intricate details.     

Wobble dG X dT pairing in right-handed DNA: solution conformation of the d(C-G-T-G-A-A-T-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

We report below on features of the three-dimensional structure of the d(C-G-T-G-A-A-T-T-C-G-C-G) self-complementary duplex (designated 12-mer GT) containing symmetrical G X T mismatches in the interior of the helix. The majority of the base and sugar protons in the 12-mer GT duplex were assigned by two-dimensional nuclear Overhauser effect (NOESY) spectra in H2O and D2O solution. A set of 92 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds for one symmetrical half of the 12-mer GT duplex were estimated from NOESY data sets recorded as a function of mixing time. These experimental distances combined with nucleotide bond length parameters were embedded into Cartesian space; several trial structures were refined to minimize bond geometry and van der Waals and chirality error. Confidence in this approach is based on the similarity of the refined structures for the solution conformation of the 12-mer GT duplex. The G and T bases pair through two imino-carbonyl hydrogen bonds, and stacking is maintained between the G X T wobble pair and adjacent Watson-Crick G X C pairs. The experimental distance information is restricted to base and sugar protons, and hence structural features such as base pair overlap, glycosidic torsion angles, and sugar pucker are well-defined by this combination of NMR and distance geometry methods. By contrast, we are unable to define the torsion angles about the bonds C3'-O3'-P-O5'-C5'-C4' in the backbone of the nucleic acid.     

NMR structural refinement of an extrahelical adenosine tridecamer d(CGCAGAATTCGCG)2 via a hybrid relaxation matrix procedure

Until very recently interproton distances from NOESY experiments have been derived solely from the two-spin approximation method. Unfortunately, even at short mixing times, there is a significant error in many of these distances. A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations avoids the approximation of the two-spin method. We have calculated the structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, by an iterative refinement approach using a hybrid relaxation matrix method combined with restrained molecular dynamics calculations. Distances from the 2D NOESY spectra have been calculated from the relaxation rate matrix which has been evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix derived distances have then been used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure is then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. Although the crystal structure of the tridecamer clearly shows the extrahelical adenosine looped out way from the duplex, the NOESY distance restrained hybrid matrix/molecular dynamics structural refinement establishes that the extrahelical adenosine stacks into the duplex.     

NMR studies of abasic sites in DNA duplexes: deoxyadenosine stacks into the helix opposite acyclic lesions

Proton and phosphorus NMR studies are reported for two complementary nonanucleotide duplexes containing acyclic abasic sites. The first duplex, d(C-A-T-G-A-G-T-A-C).d(G-T-A-C-P-C-A-T-G), contains an acyclic propanyl moiety, P, located opposite a deoxyadenosine at the center of the helix (designated APP 9-mer duplex). The second duplex, d(C-A-T-G-A-G-T-A-C).d(G-T-A-C-E-C-A-T-G), contains a similarly located acyclic ethanyl moiety, E (designated APE 9-mer duplex). The ethanyl moiety is one carbon shorter than the natural carbon-phosphodiester backbone of a single nucleotide unit of DNA. The majority of the exchangeable and nonexchangeable base and sugar protons in both the APP 9-mer and APE 9-mer duplexes, including those at the abasic site, have been assigned by recording and analyzing two-dimensional phase-sensitive NOESY data sets in H2O and D2O solution between -5 and 5 degrees C. These spectroscopic observations establish that A5 inserts into the helix opposite the abasic site (P14 and E14) and stacks between the flanking G4.C15 and G6.C13 Watson-Crick base pairs in both the APP 9-mer and APE 9-mer duplexes. The helix is right-handed at and adjacent to the abasic site, and all glycosidic torsion angles are anti in both 9-mer duplexes. Proton NMR parameters for the APP 9-mer and APE 9-mer duplexes are similar to those reported previously for the APF 9-mer duplex (F = furan) in which a cyclic analogue of deoxyribose was embedded in an otherwise identical DNA sequence [Kalnik, M. W., Chang, C. N., Grollman, A. P., & Patel, D. J. (1988) Biochemistry 27, 924-931]. These proton NMR experiments demonstrate that the structures at abasic sites are very similar whether the five-membered ring is open or closed or whether the phosphodiester backbone is shortened by one carbon atom. Phosphorus spectra of the APP 9-mer and APE 9-mer duplexes (5 degrees C) indicate that the backbone conformation is similarly perturbed at three phosphodiester backbone torsion angles. These same torsion angles are also distorted in the APF 9-mer but assume a different conformation than those in the APP 9-mer and APE 9-mer duplexes.     

NMR studies of abasic sites in DNA duplexes: deoxyadenosine stacks into the helix opposite the cyclic analogue of 2-deoxyribose

Proton and phosphorus NMR studies are reported for the complementary d(C-A-T-G-A-G-T-A-C).d(G-T-A-C-F-C-A-T-G) nonanucleotide duplex (designated APF 9-mer duplex) which contains a stable abasic site analogue, F, in the center of the helix. This oligodeoxynucleotide contains a modified tetrahydrofuran moiety, isosteric with 2-deoxyribofuranose, which serves as a structural analogue of a natural apurinic/apyrimidinic site [Takeshita, M., Chang, C.N., Johnson, F., Will, S., & Grollman, A.P. (1987) J. Biol. Chem. 262, 10171-10179]. Exchangeable and nonexchangeable base and sugar protons, including those located at the abasic site, have been assigned in the complementary APF 9-mer duplex by recording and analyzing two-dimensional phase-sensitive NOESY data sets in H2O and D2O solution at low temperature (0 degrees C). These studies indicate that A5 inserts into the helix opposite the abasic site F14 and stacks with flanking G4.C15 and G6.C13 Watson-Crick base pairs. Base-sugar proton NOE connectivities were measured through G4-A5-G6 on the unmodified strand and between the base protons of C15 and the sugar protons of the 5'-flanking residue F14 on the modified strand. These studies establish that all glycosidic torsion angles are anti and that the helix is right-handed at and adjacent to the abasic site in the APF 9-mer duplex. Two of the 16 phosphodiester groups exhibit phosphorus resonances outside the normal spectral dispersion indicative of altered torsion angles at two of the phosphate groups in the backbone of the APF 9-mer duplex.     

Experimental and theoretical studies of the conformational perturbations induced by an abasic site

The three-dimensional structure of the natural undecamer duplex d(CGCACACACGC). d(GCGTGTGTGCG) has been determined by the combined use of NMR spectroscopy and restrained molecular dynamics (rMD) and also by molecular mechanics calculations using the JUMNA program without experimental distance constraints. Both procedures have also been used to model the abasic structure d(CGCACOCACGC).d(GCGTGTGTGCG), where 'O' indicates a modified abasic site: 3-hydroxy-2-(hydroxymethyl) tetrahydrofuran. For the natural duplex, 134 interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peaks intensities, using MARDIGRAS software. These distances along with 100 torsion angles for sugar ring and additional data derived from canonical A and B-DNA, have been used for structures refinement by restrained molecular dynamics. Comparison of the natural oligomer with the abasic structure obtained earlier by NMR/rMD (Y. Coppel, N. Berthet, C. Coulombeau, Ce. Coulombeau, J. Garcia and J. Lhomme, Biochemistry 36, 4817-4830, 1997) confirms that the creation of an abasic site, in this sequence context, leads to marked helix kinking. It is also shown that the JUMNA procedure is capable of reproducing the overall structural features of the natural and damaged DNA conformations without the use of experimental constraints.     

Structural model for an oligonucleotide containing a bulged guanosine by NMR and energy minimization

We present three-dimensional structural models for a DNA oligomer containing a bulged guanosine based on proton NMR data and energy minimization computations. The nonexchangeable proton resonances of the duplex 5'd(GATGGGCAG).d(CTGCGCCATC) are assigned by nuclear Overhauser effect spectroscopy (NOESY) and correlated spectroscopy connectivities, and the NMR spectrum is compared with that of a regular 8-mer of similar sequence, 5'd(GATGGCAG).d(CTGCCATC). Experimental proton-proton distances are obtained from NOESY spectra acquired with mixing times of 100, 150, and 200 ms. A refined three-dimensional structure for the bulge-containing duplex is calculated from regular B DNA starting coordinates by using the AMBER molecular mechanics program [Weiner, S. J., Kollman, P. A., Case, D. A., Singh, U. C., Ghio, C., Alagona, G., Profeta, S., & Weiner, P. (1984) J. Am. Chem. Soc. 106, 765-784]. We compare structures obtained by building the helix in three and four base pair increments with structures obtained by direct minimization of the entire nine base sequence, with and without experimental distance constraints. The general features of all the calculated structures are very similar. The helix is of the B family, with the extra guanine stacked into the helix, and the helix axis is bent by 18-23 degrees, in agreement with gel mobility data for bulge-containing sequences [Rice, J. A. (1987) Ph.D. Thesis, Yale University].     

NMR solution structures of bistranded abasic site lesions in DNA

Ionizing radiation produces clustered lesions in DNA. Since the orientation of bistranded lesions affects their recognition by DNA repair enzymes, clustered damages are more difficult to process and thus more toxic than single oxidative lesions. In order to understand the structural determinants that lead to differential recognition, we used NMR spectroscopy and restrained molecular dynamics to solve the structure of two DNA duplexes, each containing two stable abasic site analogues positioned on opposite strands of the duplex and staggered in the 3' (-1 duplex, (AP) 2-1 duplex) or 5' (+1 duplex, (AP) 2+1 duplex) direction. Cross-peak connectivities observed in the nonexchangeable NOESY spectra indicate compression of the helix at the lesion site of the duplexes, resulting in the formation of two abasic bulges. The exchangeable proton spectra show the AP site partner nucleotides forming interstrand hydrogen bonds that are characteristic of a Watson-Crick G.C base pairs, confirming the extra helical nature of the AP residues. Restrained molecular dynamics simulations generate a set of converging structures in full agreement with the spectroscopic data. In the (AP) 2-1 duplex, the extra helical abasic site residues reside in the minor groove of the helix, while they appear in the major groove in the (AP) 2+1 duplex. These structural differences are consistent with the differential recognition of bistranded abasic site lesions by human AP endonuclease.     

Solution structure of the EcoRI DNA sequence: refinement of NMR-derived distance geometry structures by NOESY spectrum back-calculations

The solution structure of the self-complementary DNA duplex [d(CGCGAATTCGCG)]2, which contains the EcoRI restriction site sequence GAATTC at the center, has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Time-dependent nuclear Overhauser effect spectra were used to obtain the initial cross-relaxation rates between 155 pairs of protons. These initial cross-relaxation rates were converted into interproton distances and entered into a distance (bounds) matrix. A distance geometry algorithm (DSPACE) was used to create embedded starting structures and to refine these structures until they showed good agreement with the distance matrix; symmetry constraints were included in the refinement procedure, making the two strands in the refined distance geometry structures virtually identical and significantly improving the agreement with the distance matrix. The NOESY spectrum for one of these distance geometry structures was then calculated from the explicit coordinates by numerically integrating all the z-magnetization transfer pathways among neighboring protons within a specified radius. Distances in this distance geometry structure that did not agree with the experimental NOESY time course were then adjusted accordingly. This process was iterated until a good agreement between calculated and experimental NOESY spectra was reached. The final structure, which generates good agreement with the experimental NOESY spectrum, displays kinks at the C3-G4 base step and at the A6-T7 base step that appear to be similar to those reported for the EcoRI restriction site DNA bound to its endonuclease. The solution structure is not the same as the crystal structure of this DNA duplex.     

High-resolution NMR study of a synthetic DNA-RNA hybrid dodecamer containing the consensus pribnow promoter sequence: d(CGTTATAATGCG).r(CGCAUUAUAACG)

The nonsymmetrical double-helical hybrid dodecamer d(CGTTATAATGCG).r(CGCAUUAUAACG) was synthesized with solid-phase phosphoramidite methods and studied by high-resolution 2D NMR. The imino protons were assigned by one-dimensional nuclear Overhauser methods. All the base protons and H1', H2', H2", H3', and H4' sugar protons of the DNA strand and the base protons, H1', H2', and most of the H3'-H4' protons of the RNA strand were assigned by 2D NMR techniques. The well-resolved spectra allowed a qualitative analysis of relative proton-proton distances in both strands of the dodecamer. The chemical shifts of the hybrid duplex were compared to those of the pure DNA double helix with the same sequence (Wemmer et al., 1984). The intrastrand and cross-strand NOEs from adenine H2 to H1' resonances of neighboring base pairs exhibited characteristic patterns that were very useful for checking the spectral assignments, and their highly nonsymmetric nature reveals that the conformations of the two strands are quite different. Detailed analysis of the NOESY and COSY spectra, as well as the chemical shift data, indicate that the RNA strand assumes a normal A-type conformation (C3'-endo) whereas the DNA strand is in the general S domain but not exactly in the normal C2'-endo conformation. The overall structure of this RNA-DNA duplex is different from that reported for hybrid duplexes in solution by other groups (Reid et al., 1983a; Gupta et al., 1985) and is closer to the C3'-endo-C2'-endo hybrid found in poly(dA).poly(dT) and poly(rU).poly(dA) in the fiber state (Arnott et al., 1983, 1986).     

2'-deoxyribonolactone lesion in DNA: refined solution structure determined by nuclear magnetic resonance and molecular modeling

The solution conformation of the DNA duplex d(C1G2C3A4C5L6C7A8C9G10C11).d(G12C13G14T15G16T17G18T19G20C21G22 ) containing the 2'-deoxyribonolactone lesion (L6) in the middle of the sequence has been investigated by NMR spectroscopy and restrained molecular dynamics calculations. Interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peak intensities. These distances, along with torsion angles for sugar rings and additional data derived from canonical A- and B-DNA, have been used for structure refinement by restrained molecular dynamics (rMD). Six rMD simulations have been carried out starting from both regular A- and B-DNA forms. The pairwise rms deviations calculated for each refined structure are <1 A, indicating convergence to essentially the same geometry. The accuracy of the rMD structures has been assessed by complete relaxation matrix back-calculation. The average sixth-root residual index (Rx = 0.052 +/- 0.003) indicated that a good fit between experimental and calculated NOESY spectra has been achieved. Detailed analysis revealed a right-handed DNA conformation for the duplex in which both the T17 nucleotide opposite the abasic site and the lactone ring are located inside the helix. No kinking is observed for this molecule, even at the abasic site step. This structure is compared to that of the oligonucleotide with the identical sequence containing the stable tetrahydrofuran abasic site analogue that we reported previously [Coppel, Y., Berthet, N., Coulombeau, C., Coulombeau, Ce., Garcia, J., and Lhomme, J. (1997) Biochemistry 36, 4817-4830].     

A Renaissance Man: In Memoriam of Jon Widom (1955–2011)

Abstract

Assignment of the ¹H and ³¹P resonances of a decamer DNA duplex, d(CGCTTAAGCG)₂ was determined by two-dimensional COSY, NOESY and ¹H- ³¹P Pure Absorption phase Constant time (PAC) heteronuclear correlation spectroscopy. The solution structure of the decamer was calculated by an iterative hybrid relaxation matrix method combined with NOESY-distance restrained molecular dynamics. The distances from the 2D NOESY spectra were calculated from the relaxation rate matrix which were evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix-derived distances were then used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure was then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. J_H3′-P coupling constants for each of the phosphates of the decamer were obtained from ¹H-³¹P J-resolved selective proton flip 2D spectra. By using a modified Karplus relationship the C4′-C3′-03′-P torsional angles (?) were obtained. Comparison of the ³¹P chemical shifts and J_H3′-P coupling constants of this sequence has allowed a greater insight into the various factors responsible for ³¹P chemical shift variations in oligonucleotides. It also provides an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These correlations are consistent with the hypothesis that changes in local helical structure perturb the deoxyribose phosphate backbone. The variation of the ³¹P chemical shift, and the degree of this variation from one base step to the next is proposed as a potential probe of local helical conformation within the DNA double helix. The pattern of calculated ? and ζ torsional angles from the restrained molecular dynamics refinement agrees quite well with the measured J_H3′-P coupling constants. Thus, the local helical parameters determine the length of the phosphodiester backbone which in turn constrains the phosphate in various allowed conformations.     

Solution conformation of d-CTCGAGCTCGAG by two-dimensional NMR: conformational heterogeneity at XhoI cleavage site

Nonexchangeable proton resonances in the 500-MHz NMR spectrum of d-CTCGAGCTCGAG have been assigned by using two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY). 1H-1H coupling constants (J) in the deoxyribose rings have been measured by analyzing intensity and multiplet patterns in the phase-sensitive omega 1-scaled COSY spectra. A modification of the J-resolved technique, called amplitude-modulated J-resolved spectroscopy, has been described and used to increase the accuracy of J measurements. Absorption mode omega 1-scaled NOESY spectra at mixing times in the range 50-200 ms have been analyzed to monitor spin diffusion. A 50-ms spectrum has been used to estimate several interproton distances. The coupling constant and distance data have been used to arrive at sequence-specific sugar geometries and glycosidic torsion angles. The backbone structure has been refined by model building using the FRODO program, employing the sugar geometries and glycosidic torsion angles discussed above. The molecule shows interesting sequence-dependent variations in the structure. The cleavage site of the restriction enzyme XhoI exhibits unique differences in the sugar geometry and backbone torsion angles.     

The conformation of the d(ACCCGGGT) duplex in aqueous solution

The nonexchangeable base and sugar protons of the octanucleotide d(ACCCGGGT)2 have been assigned using two dimensional homonuclear Hartmann-Hahn relayed spectroscopy (HOHAHA), double quantum filtered homonuclear correlation spectroscopy (DQFCOSY) and nuclear Overhauser spectroscopy (NOESY) in D2O at 12 degrees C. The observed NOE's between the base protons and their own H2' protons and between the base protons and the H2' protons of the 5' adjacent nucleotide and the observed coupling constants between the deoxyribose 1' and 2',2' protons indicate that this duplex assumes a right-handed B-type helix conformation in solution.     

Conformational analysis of d(C3G3), a B-family duplex in solution

NMR and circular dichroism studies of the duplex formed by the self-complementary DNA hexanucleotide d(C3G3) indicate that it is a B-type structure but differs from standard B-form. An analysis of NMR coupling constants within the deoxyribose moieties yields a 70% or greater contribution from pseudorotation phase angles corresponding to the C3'-exo conformation, a conformation similar to the C2'-endo conformation associated with B-form DNA. Intranucleotide interproton distances are consistent with a B-form structure, but some internucleotide distances are intermediate between A- and B-form structures. Circular dichroism spectra have B-form characteristics but also include an unusual negative band at 282 nm. The solution spectroscopic results are in contrast with X-ray crystallographic studies which find A-form structures for similar sequences.     

Quantification of DNA structure from NMR data: conformation of d-ACATCGATGT

Resonance assignments of nonexchangeable base and sugar protons have been obtained in double-helical d-ACATCGATGT by using two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY). The exchangeable imino protons have been assigned on the basis of their chemical shifts. The characteristic phase-sensitive multiplet patterns of the intrasugar cross-peaks in the omega 1-scaled COSY spectrum have been used to estimate several scalar coupling constants (J). The information on the J values combined with the intranucleotide COSY cross-peak intensities has been used to identify sugar puckers of individual nucleotide units. In most cases, the deoxyribofuranose rings are found to adopt a conformation close to O4'-endo. Spin diffusion has been monitored from the buildup of the normalized volumes of NOE cross-peaks in NOESY spectra as a function of mixing time. A set of 52 intranucleotide and internucleotide proton-proton distances have been estimated by using low mixing time NOESY spectra (tau m = 40 and 80 ms). The estimated intrasugar proton-proton distances rule out possibilities of existence of a fast equilibrium between C2'-endo and C3'-endo conformations. Intranucleotide proton-proton distances combined with the knowledge of sugar puckers have been used to fix the glycosidic bond torsion angle (chi). For this purpose, simulated distance contours depicting the dependence of intranucleotide proton-proton distances on pseudorotational phase angle (P) and glycosidic bond torsion angle (chi) have been used. Further, the proton homonuclear (J, delta) spectrum has been used to monitor the 31P-1H heteronuclear couplings, which are preserved in the omega 2 projection.(ABSTRACT TRUNCATED AT 250 WORDS)     

Solution structure of a trinucleotide A-T-A bulge loop within a DNA duplex.

We have synthesized an oligodeoxynucleotide duplex, d(G-C-A-T-C-G-A-T-A-G-C-T-A-C-G).d(C-G-T-A-G-C-C-G-A-T-C-G), with a three-base bulge loop (A-T-A) at a central site in the first strand. Nuclear Overhauser experiments (NOESY) in H2O indicate that the GC base pairs flanking the bulge loop are intact between 0 and 25 degrees C. Nuclear Overhauser effects in both H2O and D2O indicate that all bases within the bulge loop are stacked into the helix. These unpaired bases retain an anti conformation about their glycosidic bonds as they stack within the duplex. The absence of normal sequential connectivities between the two cytosine residues flanking the bulge site on the opposite strand indicates a disruption in the geometry of this base step upon insertion of the bulged bases into the helix. This conformational perturbation is more akin to a shearing apart of the bases, which laterally separates the two halves of the molecule, rather than the "wedge" model often invoked for single-base bulges. Using molecular dynamics calculations, with both NOE-derived proton-proton distances and relaxation matrix-calculated NOESY cross peak volumes as restraints, we have determined the solution structure of an A-T-A bulge loop within a DNA duplex. The bulged bases are stacked among themselves and with the guanine bases on either side of the loop. All three of the bulged bases are displaced by 2-3 A into the major groove, increasing the solvent accessibility of these residues. The ATA-bulge duplex is significantly kinked at the site of the lesion, in agreement with previously reported electron microscopy and gel retardation studies on bulge-containing duplexes [Hsieh, C.-H., & Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A 86, 4833-4837; Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840]. Bending occurs in a direction away from the bulge-containing strand, and we find a significant twist difference of 84 degrees between the two base pairs flanking the bulge loop site. This value represents 58% of the twist difference for base pairs four steps apart in B-DNA. These results suggest a structural mechanism for the bending of DNA induced by unpaired bases, as well as accounting for the effect bulge loops may have on the secondary and tertiary structures of nucleic acids.