Investigation of the solution structure of a DNA octamer [d(GGAATTCC)]2 using two-dimensional nuclear Overhauser enhancement spectroscopy

Proton two-dimensional nuclear Overhauser enhancement (2D NOE) spectra in the pure absorption phase were obtained at 500 MHz for [d(GGAATTCC)]2 in aqueous solution at a series of mixing times. The experimental data were analyzed by comparison with theoretical spectra calculated using the complete 70 X 70 relaxation matrix including all proton dipole-dipole interactions and spin diffusion [Keepers, J. W. & James, T. L. (1984) J. Magn. Reson. 57, 404-426]. The theoretical spectra at each mixing time were calculated using two structures: a standard B-form DNA structure and an energy-minimized structure based on the similarity of the six internal residues of the title octamer with those of the dodecamer [d(CGCGAATTCGCG)]2, for which the crystal structure has been determined. Neither the standard B-form nor the energy-minimized structure will yield theoretical 2D NOE spectra which accurately reproduce all peak intensities in the experimental spectra. However, many features of the experimental spectra can be represented by both the B-form and the energy-minimized structure. Sequence-dependent structural characteristics are manifest in the 2D NOE spectra, in particular at the purine-pyrimidine junction as noted previously in the crystal structure. On the whole, the energy-minimized structure appears to yield theoretical 2D NOE spectra which mimic many, if not all, aspects of the experimental spectra. All 2D NOE data were consistent with nanosecond correction times as implied by proton spin-lattice relaxation time measurements. But better fits of some of the 2D NOE data using small variations in an effective isotropic correlation time suggest that there may be some local variations in mobility within the octamer duplex structure in solution.     

Solution structure of [d(GGTATACC)]2: wrinkled D structure of the TATA moiety

Phase-sensitive two-dimensional nuclear Overhauser effect spectra of [d(GGTATACC)]2 in aqueous deuterium oxide solution at four mixing times were quantified to give all nonoverlapping cross-peak intensities. A structural model for [d(GGTATACC)]2 was built in which the GG- and -CC moieties were in the B-DNA form, while the middle -TATA- moiety was in the wrinkled-D form (BDB model). This model was subjected to energy refinement by molecular mechanics calculations with the program AMBER. Counterions (Na+) were added to neutralize the charges, and water molecules were placed bridging across the minor groove. A complete relaxation matrix analysis was used to calculate two-dimensional nuclear Overhauser effect spectra of [d(GGTATACC)]2 from the above models (before and after energy refinement) and from four other [d(GGTATACC)]2 structural models: regular A, crystalline A, regular B, and energy-minimized B. Among them, the energy-minimized BDB model yielded a set of theoretical spectra that gave the best fit to the experimental spectra. It was also the energetically most stable. Therefore, it is a good representation of the ensemble- and time-averaged structure of the octamer in solution. This model has backbone torsion angles similar to those of B-form DNA in the GG- and -CC moieties and torsion angles similar to those of wrinkled D form DNA in the -TATA- moiety. The base stacking and base pairing are not interrupted at the junctions between the two structural moieties. Its minor groove is narrower than that of B DNA, and the solvent-accessible surface of the minor groove forms a closed hydration tunnel in the middle -TATA- segment.     

Two-dimensional nuclear Overhauser enhancement investigation of the solution structure and dynamics of the DNA octamer [d(GGTATACC)]2

The resonances of nearly all 70 of the non-exchangeable protons of the duplex [d(GGTATACC)]2 in aqueous solution are assigned by proton two-dimensional nuclear Overhauser enhancement (2D NOE) spectra obtained in pure absorption phase at 500 MHz. Experimental and theoretical 2D NOE spectra are compared at each mixing time (100, 175, 250 and 400 ms) using two B-DNA structures: a standard B-form and an energy-minimized form. The GG and CC ends of the octamer duplex are well represented by the regular B-DNA structure. But large discrepancies from these models are observed for the 'TATA' box. All 2D NOE data are consistent with nanosecond correlation times, as indicated by non-selective proton spin-lattice relaxation times, but small variations in the correlation time are observed, suggesting that there are some local differences in mobility within the octamer duplex structure in solution.     

Solution structure of the EcoRI DNA octamer containing 5-fluorouracil via restrained molecular dynamics using distance and torsion angle constraints extracted from NMR spectral simulations.

The self-complementary DNA octamer [d(GGAATUFCC)]2, containing the EcoRI recognition sequence with one of the thymines replaced by 5-fluorouracil (UF), was synthesized. Proton homonuclear two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered correlation (2QF-COSY) spectra, as well as one-dimensional spectra at different temperatures, were recorded for the octamer. Consequently, all proton resonances were assigned. The thermally induced transition from the duplex to single strands has been followed, demonstrating the stability of the duplex containing 5-fluorouracil. Simulations of the 2QF-COSY cross-peaks by means of the programs SPHINX and LINSHA were compared with experimental data, establishing scalar coupling constants for the sugar ring protons and hence sugar pucker parameters. The deoxyribose rings exhibit a dynamic equilibrium of N- and S-type conformers with 75-95% populations of the latter. Two programs used for complete relaxation matrix analysis 2D NOE spectra, CORMA and MARDIGRAS, were modified to account for the influence of the fluorines on dipolar interactions in the proton system. Quantitative assessment of the 2D NOE cross-peak intensities for different mixing times, in conjunction with the program MARDIGRAS, gave a set of interproton distances for each mixing time. The largest and smallest values of each of the interproton distances were chosen as the upper and lower bounds for each distance constraint. The distance bounds define the size of a flat-well potential function term, incorporated into the AMBER force field, which was employed for restrained molecular dynamics calculations. Torsion angle constraints in the form of a flat-well potential were also constructed from the analysis of the sugar pucker data. Several restrained molecular dynamics runs of 35 ps were performed, utilizing 284 experimental distance and torsion angle constraints and two different starting structures, energy-minimized A- and B-DNA. Convergence to similar structures with a root-mean-square deviation of 1.2 A was achieved for the central hexamer of the octamer, starting from A- and B-DNA. The average structure from six different molecular dynamics runs was subjected to final restrained energy minimization. The resulting final structure was in good agreement with the structures derived from different molecular dynamics runs and showed a substantial improvement of the 2D NOE sixth-root residual index in comparison with classical and energy-minimized B-DNA. A detailed analysis of the conformation of the final structure and comparison with structures of similar sequences, obtained by different methods, were performed.     

Solution structure of [d(A-T)5]2 via complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect spectra and molecular mechanics calculations: evidence for a hydration tunnel

Pure absorption phase proton two-dimensional nuclear Overhauser effect (2D NOE) spectra at 500 MHz have been obtained for [d(5'ATATATATAT3')]2 in deuterium oxide solution at several mixing times. The 100 nonexchangeable proton resonances have been assigned. The experimental 2D NOE spectra were compared with theoretical spectra calculated by using the complete relaxation matrix analysis method [Keepers, J. W., & James, T. L. (1984) J. Magn. Reson. 57, 404-426] and x-ray diffraction determined molecular coordinates of A, B, alternating B, left-handed B, C, D, and wrinkled D forms of DNA and of energy-minimized structures calculated from the most promising X-ray crystal structures by using the molecular mechanics program AMBER, in which all hydrogens, counterions, and hydration water molecules were included. The analysis of all features of the 2D NOE spectra played an important role in extracting the promising structures, and it was concluded that the wrinkled D form yields the best fit for the 2D NOE data of the A-T decamer. The molecular mechanics calculation indicated that this model structure, whose minor groove is comparatively deep and narrow, may be energetically more stable than the B form for alternating d(A-T) DNA. Interesting features of the structure include possible intra- and interchain sugar-phosphate attractions and a hydration tunnel inside the minor groove capable of accommodating three types of water molecules that aid in helix stabilization via hydrogen bonding. Counterions (sodium) serve to reduce interchain phosphate-phosphate repulsive effects.     

Conformational studies of d-(AAAAATTTTT)2 using constraints from nuclear overhauser effects and from quantitative analysis of the cross-peak fine structures in two-dimensional 1H nuclear magnetic resonance spectra

The conformation at the dA-dT junction in d-(AAAAATTTTT)2 was investigated by using a variety of phase-sensitive two-dimensional nuclear magnetic resonance experiments at 500 MHz for detailed studies of the deoxyribose ring puckers. Conformational constraints were collected from two-dimensional nuclear Overhauser enhancement spectra recorded with short mixing times and from quantitative simulations of the cross-peaks in two-dimensional correlated spectra. Overall, the decamer duplex adopts a conformation of the B-DNA type, and for dA4 and dA5 the pseudorotation phase angle P is in the standard range 150-180 degrees. The deoxyribose puckers for the other nucleotides deviate significantly from the standard B-DNA structure. Spectrum simulations assuming either static deviations from standard B-DNA or a simple two-state dynamic equilibrium between the C2'-endo and C3'-endo forms of the deoxyribose were used to analyze the experimental data. It was thus found that the ring pucker for dT6 deviates from the regular C2'-endo form of B-DNA by a static distortion, with the pseudorotation phase angle P in the range 100-130 degrees, and a similar value of P is indicated for dT7. For the peripheral base pairs dynamic distortions of the C2'-endo form of the deoxyribose were found. In agreement with recent papers on related duplexes containing (dA)n tracts, we observed prominent nuclear Overhauser effects between adenine-2H and deoxyribose-1'H, which could be largely due to pronounced propeller twisting as observed in the crystal structures of (dA)n-containing compounds.     

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.     

Solution structure of [d(GCGTATACGC)]2.

The solution structure of the alternating pyrimidine-purine DNA duplex [d(GCGTATACGC)]2 has been determined using two-dimensional nuclear magnetic resonance techniques and distance geometry methods. Backbone distance constraints derived from experimental nuclear Overhauser enhancement and J-coupling torsion angle constraints were required to adequately define the conformation of the inter-residue backbone linkages and to avoid underwinding of the duplex. The distance geometry structures were further refined by back-calculation of the two-dimensional nuclear Overhauser enhancement spectra to correct spin-diffusion distance errors. Fifteen final structures for [d(GCGTATACGC)]2 were generated from the refined experimental distance bounds. These structures all exhibit fully wound B-form geometry with small penalty values (< 1.5 A) against the distance bounds and small pair-wise root-mean-square deviation values (typically 0.6 A to 1.5 A). The final structures exhibit positive base-pair inclination with respect to the helix axis, a marked alternation in rise and twist, and are shorter and wider than classical fiber B-form DNA. The purines were found to adopt a sugar pucker close to the C-2'-endo conformation while pyrimidine sugars exhibited significantly lower pseudorotation phase angles in the C-1'-exo to C-2'-endo range. The minor groove cross-strand steric clashes at pyrimidine-purine steps that would exist in pure B-DNA are attenuated by an increased rise at these steps (and an increased roll angle at TpA steps). Concomitantly the backbone torsion angles of the pyrimidine moieties have larger gamma values, larger epsilon values, and smaller zeta values than the purines. The structures generated by distance geometry methods were also compared with those obtained from restrained molecular dynamics with empirical force-field potentials. The results indicate that the nuclear magnetic resonance/distance geometry approach alone is capable of elucidating most of the salient structural features of double-stranded helical nucleic acids in solution without resorting to empirical energy potentials and without using any structural assumptions from crystallographic data.     

Sequential resonance assignment of the 500-MHz nmr spectrum of d-(CG)6 and its structure under low-salt conditions

Two-dimensional (2D) nmr methods (correlated spectroscopy, nuclear Overhauser enhancement spectroscopy, and relayed correlated spectroscopy) have been used to obtain resonance assignment of the nonexchangeable base and sugar protons of a double-helical DNA segment, d-(CG)₆ in D₂O solutions under conditions of low ionic strength. Detailed information about the glycosidic torsion angle, sugar geometry, stacking patterns of the bases, and the overall solution structure of the dodecanucleotide has been obtained from the relative intensities of cross-peaks in the 2D spectra. The molecule shows general features of B-DNA under the experimental conditions employed. However, in spite of the repeating base sequence, there are subtle and detectable variations in the structure along the double helix. The terminal residues show considerable conformational flexibility.     

Solution conformation of purine-pyrimidine DNA octamers using nuclear magnetic resonance, restrained molecular dynamics and NOE-based refinement

The solution structures of two alternating purine-pyrimidine octamers, [d(G-T-A-C-G-T-A-C)]2 and the reverse sequence [d(C-A-T-G-C-A-T-G)]2, are investigated by using nuclear magnetic resonance spectroscopy and restrained molecular dynamics calculations. Chemical shift assignments are obtained for non-exchangeable protons by a combination of two-dimensional correlation and nuclear Overhauser enhancement (NOE) spectroscopy experiments. Distances between protons are estimated by extrapolating distances derived from time-dependent NOE measurements to zero mixing time. Approximate dihedral angles are determined within the deoxyribose ring from coupling constants observed in one and two-dimensional spectra. Sets of distance and dihedral determinations for each of the duplexes form the bases for structure determination. Molecular dynamics is then used to generate structures that satisfy the experimental restraints incorporated as effective potentials into the total energy. Separate runs start from classical A and B-form DNA and converge to essentially identical structures. To circumvent the problems of spin diffusion and differential motion associated with distance measurements within molecules, models are improved by NOE-based refinement in which observed NOE intensities are compared to those calculated using a full matrix analysis procedure. The refined structures generally have the global features of B-type DNA. Some, but not all, variations in dihedral angles and in the spatial relationships of adjacent base-pairs are observed to be in synchrony with the alternating purine-pyrimidine sequence.     

1H nuclear magnetic resonance assignments for d-(GCATTAATGC)2 using experimental refinements of established procedures

Sequence-specific 1H nuclear magnetic resonance assignments are presented for d-(GCATTAATGC)2. Using omega 1-scaled double quantum-filtered correlated spectroscopy, two-quantum spectroscopy, relayed coherence transfer spectroscopy and detailed analysis of the fine structure in these phase-sensitive spectra, the spin system of the bases and deoxyribose rings were identified entirely via scalar proton-proton couplings. The sequential connectivities were established with two-dimensional nuclear Overhauser enhancement spectra recorded with a short mixing time of 60 milliseconds. These spectra contain only a small number of cross-peaks, corresponding to the shortest proton-proton distances prevailing in the DNA. They are thus easy to interpret, and therefore the presently proposed modifications of the established assignment procedures should enable studies of larger DNA duplexes with intrinsically more complex nuclear magnetic resonance spectra, and they also provided an improved basis for conformational studies of DNA fragments.     

Relaxation matrix analysis of two-dimensional nuclear Overhauser effect spectra

The large number of interproton distances extracted from two-dimensional nuclear Overhauser effect spectra has enabled determination of biomolecular structures in solution. The accuracy of those distances is increased substantially and the number of distances increased significantly by analysis of the experimental peak intensities using a complete relaxation matrix approach. More distances and more accurate distances both lead to a higher resolution structure. A complete relaxation matrix analysis also enables simulation of peak intensities for any postulated structure; comparison of these intensities with experimental intensities can provide a guide for structure refinement as well as a measure of the quality of the structure derived.     

Solution conformation of an oligonucleotide containing a G.G mismatch determined by nuclear magnetic resonance and molecular mechanics.

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three dimensional solution structure of the non-selfcomplementary oligonucleotide, d(GAGGAGGCACG). d(CGTGCGTCCTC) in which the central base pair is G.G. This is the first structural determination of a G.G mismatch in a oligonucleotide. Two dimensional nuclear magnetic resonance spectra show that the bases of the mismatched pair are stacked into the helix and that the helix adopts a classical B-DNA form. Spectra of the exchangeable protons show that the two guanosines are base paired via their imino protons. For the non-exchangeable protons and for some of the exchangeable protons nuclear Overhauser enhancement build up curves at short mixing times have been measured. These give 84 proton-proton distances which are sensitive to the helix conformation. One of the guanosines adopts a normal anti conformation while the other is syn or close to syn. All non-terminal sugars are C2' endo. These data sets were incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. The G.G mismatch shows a symmetrical base pairing structure. Although the mismatch is very bulky many of its features are close to that of normal B-DNA. The mismatch induces a small lateral shift in the helix axis and the sum of the helical twist above and below the mismatch is close to that of B-DNA.     

Sequential resonance assignments in protein 1H nuclear magnetic resonance spectra: Computation of sterically allowed proton-proton distances and statistical analysis of proton-proton distances in single crystal protein conformations

Two different, theoretical studies of intramolecular proton-proton distances in polypeptide chains are described. Firstly, the distances between amide, C^α and C^β protons of neighbouring residues in the amino acid sequence, which correspond to the sterically allowed values for the dihedral angles φ_i, ψ_i and χ_i¹, were computed. Secondly, the frequency with which short distances occur between amide, C^α and C^β protons of neighbouring and distant residues in the amino acid sequence were statistically evaluated in a representative sample of globular protein crystal structures. Both approaches imply that semi-quantitative measurements of short, non-bonding proton-proton distances, e.g. by nuclear Overhauser experiments, should present a reliable and generally applicable method for sequential, individual resonance assignments in protein ¹H nuclear magnetic resonance spectra. Similar calculations imply that corresponding distance measurements can be used for resonance assignments in the side-chains of the aromatic amino acid residues, asparagine and glutamine, where the complete spin systems cannot usually be identified from through-bond spin-spin coupling connectivities.     

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.     

Solution structure of [d(GTATATAC)]2 via restrained molecular dynamics simulations with nuclear magnetic resonance constraints derived from relaxation matrix analysis of two-dimensional nuclear Overhauser effect experiments

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.     

Determination of the conformation of d(GGAAATTTCC)2 in solution by use of 1H NMR and restrained molecular dynamics

The conformation of the putative bent DNA d(GGAAATTTCC)2 in solution was studied by use of 1H NMR and restrained molecular dynamics. Most of the resonances were assigned sequentially. A total of 182 interproton distance restraints were determined from two-dimensional nuclear Overhauser effect spectra with short mixing times. Torsion angle restraints for each sugar moiety were determined by qualitative analysis of a two-dimensional correlated spectrum. Restrained molecular dynamics was carried out with the interproton distances and torsion angles incorporated into the total energy function of the system in the form of effective potential terms. As initial conformations for restrained molecular dynamics, classical A-DNA and B-DNA were adopted. The root mean square deviation (rmsd) between these two conformations is 5.5 A. The conformations obtained by use of restrained molecular dynamics are very similar to each other, the rmsd being 0.8 A. On the other hand, the conformations obtained by use of molecular dynamics without experimental restraints or restrained energy minimization depended heavily on the initial conformations, and convergence to a similar conformation was not attained. The conformation obtained by use of restrained molecular dynamics exhibits a few remarkable features. The second G residue takes on the BII conformation [Fratini, A. V., Kopka, M. L., Drew, H. R., & Dickerson, R. E. (1982) J. Biol. Chem. 257, 14686-14707] rather than the standard BI conformation. There is discontinuity of the sugar puckering between the eighth T and ninth C. The minor groove of the oligo(dA) tract is rather compressed. As a result, d(GGAAATTTCC)2 is bent.     

Structural analysis of the decanucleotide duplex d(GGTAATTACC)2 using 1H NMR spectroscopy.

The conformation of the decanucleotide duplex d(GGTAATTACC)2 has been investigated in solution by one- and two-dimensional proton NMR spectroscopy. Intra- and inter-nucleotide two-dimensional nuclear Overhauser enhancement data, recorded at mixing times between 15 and 250 ms, reveal a right-handed B-DNA structure. The data also show that the A-T basepairs of the TAATTA tract are highly propeller twisted and the minor groove is particularly narrow.     

The structure of the double-stranded RNA pentamer 5'(CACAG) . 5'(CUGUG) determined by nuclear Overhauser enhancement measurements: interproton distance determination and structure refinement on the basis of X-ray coordinates

The structure in solution of the duplex RNA pentamer 5'(CACAG) . 5'(CUGUG), comprising the stem of the T psi C loop of yeast tRNAPhe, has been investigated by means of one- and two-dimensional nuclear Overhauser enhancement measurements. All non-exchangeable base and sugar proton resonances with the exception of the H5'/H5" sugar resonances are assigned in a sequential manner. From the relative intensities of the cross-peaks obtained in the pure-phase absorption two-dimensional nuclear Overhauser enhancement spectra at several mixing times, it is deduced that the RNA pentamer adopts an A-type conformation in solution. Cross-relaxation rates and interproton distances are determined from the time dependence of the nuclear Overhauser effects, principally by one-dimensional measurements. The structure of the RNA pentamer is then refined by restrained least-squares minimization on the basis of both distance and planarity restraints using fibre diffraction data as an initial model. The refined structure of the RNA pentamer is of the A type but exhibits local structural variations in glycosidic bond and backbone torsion angles as well as in propeller twist, base roll and base tilt angles.     

Structure refinement of oligonucleotides by molecular dynamics with nuclear Overhauser effect interproton distance restraints: application to 5' d(C-G-T-A-C-G)2

The solution structure of the self-complementary DNA hexamer 5' d(C-G-T-A-C-G)2 is refined by restrained molecular dynamics in which 192 interproton distances, determined from pre-steady-state nuclear Overhauser enhancement measurements, are incorporated into the total energy of the system in the form of effective potentials. First the method is tested by applying an idealized set of distance restraints taken from classical B-DNA to a simulation starting off from A-DNA and vice versa. It is shown that in both cases the expected transition between A- and B-DNA occurs. Second, a set of restrained molecular dynamics calculations is carried out starting from both A- and B-DNA with the experimental interproton distances for 5' d(C-G-T-A-C-G)2 as restraints. Convergence to the same B-type structure is achieved with the interproton distances equal to the measured values within experimental error. The root-mean-square atomic difference between the two average restrained dynamics structures (less than 1 A) is approximately the same as the root-mean-square fluctuations of the atoms.