Determination of nucleic acid backbone conformation by 1H NMR.

In DNA or RNA duplexes, the six-bond C3'-O3'-P-O5'-C5'-C4'-C3' backbone linkage connecting adjacent residues contains six torsion angles (epsilon, zeta, alpha, beta, gamma, delta) but only four protons. This seriously limits the ability to define the backbone conformation by NMR using purely 1H-1H distance geometry (DG) methods. The problem is further compounded by the inability to assign two of the four backbone protons, namely the poorly resolved H5' and H5' protons, and invariably leads to DG structures with poorly defined backbone conformations. We have developed and tested a reliable method to constrain the beta, gamma, and epsilon (and indirectly alpha and zeta) backbone torsion angles by lower-bound NOE distances to unassigned H5'/H5' resonances combined with either 1H line widths or the conservative use of sigma J measurements; the method relies only on 1H 2-D NMR data, does not involve any structural assumptions, and leads to much improved backbone convergence among DG structures. The C4'-C5' torsion angle gamma is constrained by lower-bound NOE distances from H2' and from H6/H8 to any H5'/H5', as well as by sigma JH4, coupling measurements in the 3.9-4.4 ppm region; delta is constrained by H1'-H4' NOE distances and by H3'-H4' and H3'-H2' J couplings in COSY data; epsilon is partially constrained by H3' line width and/or further constrained by subtracting the minimum possible sigma JH3'-H from the observed sigma JH3' (COSY) to arrive at the maximum possible JH3'-P, which is then converted to H3'-P distance bounds. The angle beta is partially constrained via H5'-P and H5'-P distance bounds consistent with the maximum H5'-P and H5'-P J couplings derived from the observed H5' and H5' line widths, while alpha and zeta are indirectly constrained by lower distance bounds on the observed (n)H1' to (n + 1)H5'/H5' NOEs combined with the prior partial constraints on beta, gamma, delta, and epsilon. The combined effects of these additional constraints in determining distance geometry structures have been demonstrated using a 12-base duplex, [d(GCCGTTAACGGC)]2. Coordinate RMSDs per atom between structures refined with these constraints from random-embedded DG structures, from ideal A-DNA, and from B-DNA starting structures were less than 0.4 A for the central 8 base pairs indicating good convergence. All backbone angles for the central 8 base pairs are very well constrained with less than 10 degrees variation in any of the 48 torsion angles.     

The occurence of the syn-C3' endo conformation and the distorted backbone conformations for C4'-C5' and P-O5' in oligo and polynucleotides.

The nucleoside constituents of nucleic acids prefer the anti conformation (1). When the sugar pucker is taken into account the nucleosides prefer the C2'endo-anti conformation. Of the nearly 300 nucleosides known, about 250 are in the anti conformation and 50 are in the syn-conformation, i.e., anti to syn conformation is 5:1. The nucleotide building blocks of nucleic acids show the same trend as nucleosides. Both the deoxy-guanosine and riboguanosine residues in nucleosides and nucleotides prefer the syn-C2'endo conformation with an intra-molecular hydrogen bond (for nucleosides) between the O5'-H and the N3 of the base and, a few syn-C3'endo conformations are also observed. Evidence is presented for the occurrence of the C3'endo-syn conformation for guanines in mis-paired double helical right-handed structures with the distorted sugar phosphate C4'-C5' and P-O5' bonds respectively, from g+ (gg) and g- to trans. Evidence is also provided for guanosine nucleotides in left-handed double-helical (Z-DNA) oligo and polynucleotides which has the same syn-C3'endo conformation and the distorted backbone sugar-phosphate bonds (C4'-C5' and P-O5') as in the earlier right-handed case.     

Molecular structures of two new anti-HIV nucleoside analogs: 9-(2,3-dideoxy-2-fluoro-beta-D-threo-pentofuranosyl)adenine and 9-(2,3-dideoxy-2-fluoro-beta-D-threo-pentofuranosyl)hypoxanthine.

The x-ray crystal structures of two new anti-HIV compounds, 9-(2,3-dideoxy-2-fluoro-beta-D-threo-pentofuranosyl)adenine (2'-F-dd-araA) and 9-(2,3-dideoxy-2-fluoro-beta-D-threo- pentofuranosyl)hypoxanthine (2'-F-dd-aral), have been determined at two temperatures. Both crystals are in the space group P2(1)2(1)2(1), and their structures were solved by direct methods. Least-squares refinement produced final R-factors of 0.027 for the 2'-F-dd-araA structure and of 0.044 for the 2'-F-dd-aral structure, respectively. The latter structure contains a two-fold disordered conformation of the sugar moiety. All three conformers (one for 2'-F-dd-araA and two for 2'-F-dd-aral) adopt an anti chi CN glycosyl torsion angle. The sugar in the 2'-F-dd-araA structure has a C2'-endo pucker conformation, whereas the sugar in the 2'-F-dd-aral structure has a mixture of C2'-endo and C3'-endo pucker conformations. When the sugar adopts the C2'-endo conformation, the torsion angle about the C4'-C5' bond is in a transgauche+ conformation. In contrast, when the sugar adopts the C3'-endo conformation, the torsion angle about the C4'-C5' bond is in a gauche(+)-gauche- conformation. The C2'-F bond distance is 1.406(3) A, similar to that found in other aliphatic C-F bonds. The results suggest that the 2'-fluoro-2',3'-dideoxyarabinosyl nucleosides do not have a strong preference for either C2'-endo or C3'-endo sugar pucker.     

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.     

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.     

The structure of 3'-O-anthraniloyladenosine, an analogue of the 3'-end of aminoacyl-tRNA.

3'-O-Anthraniloyladenosine, an analogue of the 3'- terminal aminoacyladenosine residue in aminoacyl-tRNAs, was prepared by chemical synthesis, and its crystal structure was determined. The sugar pucker of 3'-O-anthraniloyladenosine is 2'-endo resulting in a 3'-axial position of the anthraniloyl residue. The nucleoside is insynconformation, which is stabilized by alternating stacking of adenine and benzoyl residues of the neighboring molecules in the crystal lattice. The conformation of the 5'-hydroxymethylene in 3'-O- anthraniloyladenosine is gauche-gauche. There are two intramolecular and two intermolecular hydrogen bonds and several H-bridges with surrounding water molecules. The predominant structure of 3'-O-anthraniloyladenosine in solution, as determined by NMR spectroscopy, is 2'-endo,gauche-gauche and anti for the sugar ring pucker, the torsion angle around the C4'-C5'bond and the torsion angle around the C1'-N9 bond, respectively. The 2'-endo conformation of the ribose in 2'(3')-O-aminoacyladenosine, which places the adenine and aminoacyl residues in equatorial and axial positions, respectively, could serve as a structural element that is recognized by enzymes that interact with aminoacyl-tRNA or by ribosomes to differentiate between aminoacylated and non-aminoacylated tRNA.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Flexibility of 3',5' deoxyribonucleoside diphosphates

In 3',5' deoxyribonucleoside diphosphates, in addition to the nature of the base and the sugar puckering, there are six single bond rotations. However, from the analysis of crystal structure data on the constituents of nucleic acids, only three rotational angles, that are about glycosyl bond, about C4'-C5' and about C3'-O3' bonds, are flexible. For a given sugar puckering and a base, potential energy calculations using non-bonded, electrostatic and torsional functions were carried out by varying the three torsion angles. The energies are represented as isopotential energy surfaces. Since the availability of the real-time color graphics, it is possible to analyse these isopotential energy surfaces. The calculations were carried out for C3' exo and C3' endo puckerings for deoxyribose and also for four bases. These calculations throw more light not only on the allowed regions for the three rotational angles but also on the relationships among them. The dependence of base and the puckering of the sugar on these rotational angles and thereby the flexibility of the 3',5' deoxyribonucleoside diphosphates is discussed. From our calculations, it is now possible to follow minimum energy path for interconversion among various conformers.     

DFT studies of the disaccharide, alpha-maltose: relaxed isopotential maps

The disaccharide, alpha-maltose, forms the molecular basis for the analysis of the structure of starch, and determining the conformational energy landscape as the molecule oscillates around the glycosidic bonds is of importance. Thus, it is of interest to determine, using density functionals and a medium size basis set, a relaxed isopotential contour map plotted as a function of the phi(H) and psi(H) dihedral angles. The technical aspects include the method of choosing the starting conformations, the choice of scanning step size, the method of constraining the specific dihedral angles, and the fitting of data to obtain well defined contour maps. Maps were calculated at the B3LYP/6-31+G( *) level of theory in 5 degrees intervals around the (phi(H),psi(H))=(0 degrees ,0 degrees ) position, out to approximately +/-30 degrees or greater, for gg-gg'-c, gg-gg'-r, gt-gt'-c, gt-gt'-r, tg-tg'-c, and tg-tg'-r conformers, as well as one-split gg(c)-gg'(r) conformer. The results show that the preferred conformation of alpha-maltose in vacuo depends strongly upon the hydroxyl group orientations ('c'/'r'), but the energy landscape moving away from the minimum-energy position is generally shallow and transitions between conformational positions can occur without the addition of significant energy. Mapped deviations of selected parameters such as the dipole moment; the C1-O1-C4', H1-C1-O1, and H4'-C4'-O1 bond angles; and deviations in hydroxymethyl rotamers, O5-C5-C6-O6, O5'-C5'-C6'-O6', C5-C6-O6-H, and C5'-C6'-O6'-H', are presented. These allow visualization of the structural and energetic changes that occur upon rotation about the glycosidic bonds. Interactions across the bridge are visualized by deviations in H(O2)...O3', H(O3')...O2, and H1...H4' distances and the H(O2)-O2-C2-C1 and H'(O3')-O3'-C3'-C4' hydroxyl dihedral angles.     

Molecular mechanics and dynamics of an abasic frameshift in DNA and comparison to NMR data

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide: [sequence; see text] where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2' endo or C3' endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4. In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium. Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.     

Determination of the DNA sugar pucker using 13C NMR spectroscopy

Solid-state 13C NMR spectroscopy of a series of crystalline nucleosides and nucleotides allows direct measurement of the effect of the deoxyribose ring conformation on the carbon chemical shift. It is found that 3'-endo conformers have 3' and 5' chemical shifts significantly (5-10 ppm) upfield of comparable 3'-exo and 2'-endo conformers. The latter two conformers may be distinguished by smaller but still significant differences in the carbon chemical shifts at the C-2' and C-4' positions. High-resolution solid-state NMR of three modifications of fibrous calf thymus DNA shows that these trends are maintained in high-molecular-weight DNA and confirms that the major ring pucker in A-DNA is 3'-endo, while both B-DNA and C-DNA are largely 2'-endo. The data show that 13C NMR spectroscopy is a straightforward and useful probe of DNA ring pucker in both solution and the solid state.     

High-resolution pyrimidine- and ribose-specific 4D HCCH-COSY spectra of RNA using the filter diagonalization method

The NMR spectra of nucleic acids suffer from severe peak overlap, which complicates resonance assignments. 4D NMR experiments can overcome much of the degeneracy in 2D and 3D spectra; however, the linear increase in acquisition time with each new dimension makes it impractical to acquire high-resolution 4D spectra using standard Fourier transform (FT) techniques. The filter diagonalization method (FDM) is a numerically efficient algorithm that fits the entire multi-dimensional time-domain data to a set of multi-dimensional oscillators. Selective 4D constant-time HCCH-COSY experiments that correlate the H5-C5-C6-H6 base spin systems of pyrimidines or the H1'-C1'-C2'-H2' spin systems of ribose sugars were acquired on the (13)C-labeled iron responsive element (IRE) RNA. FDM-processing of these 4D experiments recorded with only 8 complex points in the indirect dimensions showed superior spectral resolution than FT-processed spectra. Practical aspects of obtaining optimal FDM-processed spectra are discussed. The results here demonstrate that FDM-processing can be used to obtain high-resolution 4D spectra on a medium sized RNA in a fraction of the acquisition time normally required for high-resolution, high-dimensional spectra.     

Effect of distortions in the deoxyribose phosphate backbone conformation of duplex oligodeoxyribonucleotide dodecamers containing GT, GG, GA, AC, and GU base-pair mismatches on 31P NMR spectra

We have previously suggested that variations in the 31P chemical shifts of individual phosphates in duplex oligonucleotides are attributable to torsional angle changes in the deoxyribose phosphate backbone. This hypothesis is not directly supported by analysis of the 1H/31P two-dimensional J-resolved spectra of a number of mismatch dodecamer oligonucleotide duplexes including the following sequences: d-(CGTGAATTCGCG), d(CGUGAATTCGCG), d(CGGGAATTCGCG), d(CGAGAATTCGCG), and d(CGCGAATTCACG). The 31P NMR signals of the dodecamer mismatch duplexes were assigned by 2D 1H/31P pure absorption phase constant time (PAC) heteronuclear correlation spectra. From the assigned H3' and H4' signals, the 31P signals of the base-pair mismatch dodecamers were identified. JH3'-P coupling constants for each of the phosphates of the dodecamers were obtained from 1H/31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. JH3'-P coupling constants were measured for many of the oligonucleotides as a function of temperature. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. This correlation can be further extended to the C3'-O3'-P-O5' torsional angle (zeta) by using a linear relationship between epsilon and zeta obtained from crystal structure studies. The 31P chemical shifts follow the general observation that the more internally the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. In addition, 31P chemical shifts show sequence- and site-specific variations. Analysis of the backbone torsional angle variations from the coupling constant analysis has provided additional information regarding the origin of these variations in 31P chemical shifts.     

Conformational microheterogeneity in a DNA double helix: structure of restriction endonuclease Bam H1 recognition site

Structural studies using 500 MHz 1H NMR spectroscopy on Bam H1 recognition site d(GGATCC)2 in solution at 19 degrees is reported. The resonances from the sugar ring and base protons have been assigned from the 2D-COSY and NOESY spectra. Analyses of the NOESY cross-peaks between the base protons H8/H6 and sugar protons H2'/H2", H3' reveal that the nucleotide units G2, A3 and C6 adopt (C3'-endo, chi = 200 degrees-220 degrees) conformation while G1, T4 and C5 exhibit (C2'-endo, chi = 240 degrees-260 degrees) conformation. NMR data clearly suggest that the two strands of d(GGATCC)2 are conformationally equivalent and there is a structural two-fold between the two A-T pairs. The above information and the NOESY data are used to generate a structural model of d(GGATCC)2. The important features are: (i) G1-G2 stack, the site of cleavage, shows an alternation in sugar pucker i.e. C2'-endo, C3'-endo as in a B-A junction, (ii) G2-A3 stack adopts a mini A-DNA, both the sugars being C3'-endo, (iii) A3-T4 stack, the site of two-fold, displays an A-B junction with alternation in sugar pucker as C3'-endo, C2'-endo, (iv) T4-C5 stack adopts a mini B-DNA both the sugars being C2'-endo and (v) C5-C6 stack exhibits a B-A junction with C2'-endo, C3'-endo sugar puckers. Thus, our studies demonstrate that conformational microheterogeneity with a structural two fold, is present in the Bam H1 recognition site.     

Polypurine/polypyrimidine hairpins form a triple helix structure at low pH.

1D and 2D NMR investigations of the 15 residue deoxynucleotide sequence d(TCTCTC-TTT-GAGAGA) show that above pH = 6.5 the molecule adopts a B-form hairpin conformation. As the pH is lowered below 6.5 molecules progressively associate in pairs to form a partially triple helical, partially single stranded structure in which the bases of the oligopyrimidine d(TC)3 tract from one molecule form Hoogsteen pairs with the d(G-A)3 tract of the other. Imino protons of protonated cytosines can be observed at very low field (approximately 15 ppm). The enthalpy of triplex formation was estimated by NMR techniques to be -16 kcal mol-1. Intense H6 to H3' cross peaks from residues in all three strands suggest the presence of N-type sugars at some but not at all possible sites. Surprisingly strong cross peaks between H5' or H5" and non-exchangeable base protons are also observed. These suggest that certain of the O5'-C5'-C4'-C3' phosphate backbone torsion angles (gamma) are unusual.     

Simulations of the B-DNA molecular dynamics of d(CGCGAATTCGCG)2 and d(GCGCGCGCGC)2: an analysis of the role of initial geometry and a comparison of united and all-atom models

Molecular dynamics simulations on the sequence d(CGCGAATTCGCG)2 have been carried out using both united atom and all-atom representations, and starting the simulations both from a regular repeating B-DNA structure and from the x-ray single crystal B-DNA structure. An all-atom B-DNA simulation on the sequence d(GCGCGCGCGC)2 has also been carried out, in order to compare it with a previous united atom simulation. The helix repeats, H-bonding, sugar pucker profiles, and average torsional angles are all in the range observed in crystallographic and nmr studies for B-DNA helices. In some of the sequences, there is a significant bend in the DNA helices. The individual helix repeats, with focus on 3'CpG5' and 3'GpC5' units, show the opposite helix repeat to that suggested by Calladine's rules.     

Mechanisms of chain folding in nucleic acids. The (omega, omega) plot and its correlation to the nucleotide geometry in yeast tRNAPhe1.

The (omega', omega) polot depicting the internucleotide P-O bond rotation angles in yeast phenylalanyl transfer RNA has established the interdependence of the phosphodiesters and the nucleotide geometries in the folding of the polynucleotide backbone. The plot distinguishes the regions characteristic of secondary helical structures and tertiary structural loops and bends. The folding of the polynucleotide chain is accomplished either solely by rotations around the P-O bonds or in concert with rotations around the nucleotide C4'-C5' bond with or without changes in the sugar ring pucker. In spite of differences in nucleotide sequence and intraloop tertiary interactions in the anticodon and pseudouridine loops, a characteristic repeating structural unit is found for the sugar-phosphate backbone of the tetranucleotide segment around the sharp turns.     

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.