The solution structure of the circular trinucleotide cr(GpGpGp) determined by NMR and molecular mechanics calculation.

The 3'-5' circular trinucleotide cr(GpGpGp) was studied by means of 1D and 2D high resolution NMR techniques and molecular mechanics calculations. Analysis of the J-couplings, obtained from the 1H and 13C-NMR spectra, allowed the determination of the conformation of the sugar rings and of the 'circular' phosphate backbone. In the course of the investigations it was found that the Karplus-equation most recently parametrized for the CCOP J-coupling constants could not account for the measured J(C4'P) of 11.1 Hz and a new parametrization for both HCOP and CCOP coupling constants is therefore presented. Subsequent analysis of the coupling constants yielded 'fixed' values for the torsion angles beta and delta (with beta = 178 degrees and delta = 139 degrees). The value of the latter angle corresponds to an S-type sugar conformation. The torsion angles gamma and epsilon are involved in a rapid equilibrium in which they are converted between the gauche(+) and trans and between the trans and gauche(-) domain respectively. We show that the occurrence of epsilon in the gauche(-) domain necessitates S-type sugar conformations. Given the aforementioned values for beta, gamma, delta and epsilon the ring closure constraints for the ring, formed by the phosphate backbone can only be fulfilled if alpha and zeta adopt some special values. After energy minimization with the CHARMm force field only two combinations of alpha and zeta result in energetically favourable structures, i.e. the combination alpha (t)/zeta(g-) in case gamma is in a gauche(+) and epsilon is in a trans conformation, and the combination alpha (t)/zeta (g+) for the combination gamma (t)/epsilon (g-). The results are discussed in relation to earlier findings obtained for cd(ApAp) and cr(GpGp), the latter molecule being a regulator of the synthesis of cellulose in Acetobacter xylinum.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

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 crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) shows significant differences to that found in solution.

The crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) has been solved and refined at 2.5 A resolution. The refinement procedure converged at R = 0.181 for all reflections in the range 20.0-2.5 A. In the crystal, the RNA/DNA hybrid duplex has an A' conformation with all but one of the nucleotide sugar moieties adopting a C3'- endo (N) conformation. Both strands in the double helix adopt a global conformation close to the A-form and the width of the minor groove is typical of that found in the crystal structures of other A-form duplexes. However, differences are observed between the RNA and DNA strands that make up the hybrid at the local level. In the central portion of the duplex, the RNA strand has backbone alpha, beta and gamma torsion angles that alternate between the normal gauche -/ trans / gauche + conformation and an unusual trans / trans / trans conformation. Coupled with this so-called 'alpha/gamma flipping' of the backbone torsion angles, the distance between adjacent phosphorous atoms on the RNA strand systematically varies. Neither of these phenomena are observed on the DNA strand. The structure of the RNA/DNA hybrid presented here differs significantly from that found in solution for this and other sequences. Possible reasons for these differences and their implications for the current model of RNase H activity are discussed.     

Conformational analysis of the 3'-5'-cyclic dinucleotide d less than pApA greater than by means of molecular mechanics

The conformational properties of the cyclic dinucleotide d less than pApA greater than were studied by means of molecular mechanics calculations in which a multiconformation analysis was combined with minimum energy calculations. In this approach models of possible conformers are built by varying the torsion angles of the molecule systematically. These models are then subjected to energy minimization; in the present investigation use was made of the AMBER Force field. It followed that the lowest energy conformer has a pseudo-two-fold axis of symmetry. In this conformer the deoxyribose sugars adopt a N-type conformation. The conformation of the sugar-phosphate backbone is determined by the following torsion angles: alpha +, beta t, gamma +, epsilon t and zeta +. The conformation of this ringsystem corresponds to the structure derived earlier by means of NMR spectroscopy and X-ray diffraction. The observation of a preference for N-type sugar conformations in this molecule can be explained by the steric hindrance induced between opposite H3' atoms when one sugar is switched from N- to S-type puckers. The sugars can in principle switch from N- to S-type conformations, but this requires at least the transition of gamma + to gamma -. In this process the molecule obtains an extended shape in which the bases switch from a pseudo-axial to a pseudo-equatorial position. The calculations demonstrate that, apart from the results obtained for the lowest energy conformation, the 180 degrees change in the propagation direction of the phosphate backbone can be achieved by several different combinations of the backbone torsion angles. It appeared that in the low energy conformers five higher order correlations are found. The combination of torsion angles which are involved in changes in the propagation direction of the sugar-phosphate backbone in DNA-hairpin loops and in tRNA, are found in the dataset obtained for cyclic d less than pApA greater than. It turns out, that in the available examples, 180 degrees changes in the backbone direction are localized between two adjacent nucleotides.     

Boltzmann-type distribution of side-chain conformation in proteins

We analyze packing imperfections in globular proteins as reflected in deviations of torsion angles from the equilibrium values for the isolated side chains. The distribution of conformations of methionine and lysine residues in a database of high-resolution structures is compared with energies of model compounds calculated with high-level quantum-mechanics. The distribution of the C-C and C-S torsion angles (chi(3)) correlates well with the Boltzmann factor of the torsion energy, exp(-betaE) of the model compounds C(2)H(5)-C(2)H(5) and C(2)H(5)-S-CH(3). An exponential relation was again found between the relative occurrence of g+, g- and t conformations for C(alpha)-C(beta) bonds in long side chains and the energy differences of rotamers of alpha-amino n-butyric acid, when dependence on backbone conformation was taken into account. The distribution of all 27 rotamers of methionine was correlated with the energy differences between the model's rotamers, corrected for clashes with nearby residues, the correlation being good for a set with backbone in the beta-conformation, but less clear for backbone alpha-conformation. In all correlations, the value of the coefficient beta corresponds to a temperature of circa 300 K. These results can be interpreted with a model that considers the structure of a folded protein as resulting from packing imperfectly complementary parts, with a requirement of an overall low energy. Compromises are required to optimize the fit of nonbonded contacts with surrounding groups, and side chains assume conformations away from the energy minimum. An exponential distribution is a most probable distribution, and this can be established easily under conditions other than thermal equilibrium.     

Solution conformation of a deoxynucleotide containing tandem G.A mismatched base pairs and 3'-overhanging ends in d(GTGAACTT)2.

We have used 31P and 1H NMR spectroscopy and circular dichroism to define the solution conformation of d(GTGAACTT)2 which contains tandem G.A mismatched base pairs and 3'-overhanging TT ends. Measurements of coupling constants and NOE intensities show that the sugar puckers of the nucleotides are predominantly in the south domain (i.e., near C2'-endo) and that the glycosidic torsion angles are anti. The sequential NOE intensities indicate the presence of a right-handed helix. Analysis of the 31P and 1H NMR spectra of the duplex shows that the tandem mismatch forms a block in which there are unusual backbone torsion angles (i.e., in the BII state), within an otherwise B-like structure. The chemical shift of the N1H of the mismatched guanosine and NOEs between the mismatched base pairs and their nearest neighbors are inconsistent with the imino pairing present in single A.G mismatches or in the X-ray structure of a tandem mismatch [Privé, G. G., et al. (1987) Science 238, 498-503] but the data are consistent with the amino pairing found by Li et al. (1991) [Li, Y., et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The strong base-base stacking both within the tandem G.A block and between the G.A mismatches and their other nearest neighbors offsets the intrinsic destabilizing effects of the mismatch. Further, the 3'-TT overhangs stack onto the ends of the helix and stabilize the duplex against fraying, which accounts for the observed increase in the melting temperature compared with the flush-ended duplex.     

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.     

Dinucleoside monophosphate having a high anti conformation: the crystal structure of 8,2'-S-cycloadenylyl-(3'-5')-8,2'-S-cycloadenosine hydrochloride

The crystal and molecular structure of 8,2'-S-cycloadenylyl-(3'-5')-8,2'-S-cycloadenosine (AspAs) hydrochloride has been determined by X-ray method. The conformation of two independent AspAs molecules found in an asymmetric unit are almost identical to each other. The torsion angles concerning the sugar-phosphate backbone are different from those in crystalline dinucleoside monophosphates so far determined by X-rays. Both AspAs molecules are in the sharp bend conformations, i.e. each rotation around P-O bond (omega', omega) is (g-, t) rather than the preferred (g-, g-) or (g+, g+) conformation. There is no intramolecular base stacking or base-pairing but the intermolecular base stacking was found.     

The solution structure of the intramolecular photoproduct of d(TpA) derived with the use of NMR and a combination of distance geometry and molecular dynamics.

One and two dimensional NMR techniques have been used together with molecular modelling to obtain the solution structure for the photoproduct d(TpA)*. The NMR data confirm that the cyclobutane linkage is formed between the bonds thymine C6-C5 and adenine C5-C6. The 2D NOE data are used as constraints in a distance geometry calculation. The structures obtained show a trans-syn cyclobutane linkage and the glycosidic angles are SYN and ANTI for thymidine and deoxyadenosine, respectively. The coupling constant data are used to check the backbone torsion angles of the obtained structures. Typical torsion angles are a gamma+ and beta t for the deoxyadenosine residue. A free molecular dynamics simulation of a trans-syn d(TpA) photoproduct confirmed all these structural characteristics.     

Structural characterization of d(CAACCCGTTG) and d(CAACGGGTTG) mini-hairpin loops by heteronuclear NMR: the effects of purines versus pyrimidines in DNA hairpins.

The DNA decamers, d(CAACCCGTTG) and d(CAACGGGTTG) were studied in solution by proton and heteronuclear NMR. Under appropriate conditions of pH, temperature, salt concentration and DNA concentration, both decamers form hairpin conformations with similar stabilities [Avizonis and Kearns (1995) Biopolymers, 35, 187-200]. Both decamers adopt mini-hairpin loops, where the first and last four nucleotides are involved in Watson-Crick hydrogen bonding and the central two nucleotides, CC or GG respectively, form the loop. Through the use of proton-proton, proton-phosphorus and natural abundance proton-carbon NMR experiments, backbone torsion angles (beta, gamma and epsilon), sugar puckers and interproton distances were measured. The nucleotides forming the loops of these decamers were found to stack upon one another in an L1 type of loop conformation. Both show gamma tr and unusual beta torsion angles in the loop-closing nucleotide G7, as expected for mini-hairpin loop formation. Our results indicate that the beta and epsilon torsion angles of the fifth and sixth nucleotides that form the loop and the loop-closing nucleotide G7 are not in the standard trans conformation as found in B-DNA. Although the loop structures calculated from NMR-derived constraints are not well defined, the stacking of the bases in the two different hairpins is different. This difference in the base stacking of the loop may provide an explanation as to why the cytosine-containing hairpin is thermodynamically more stable than the guanine-containing hairpin.     

Conformation of the glycotripeptide repeating unit of antifreeze glycoprotein of polar fish as determined from the fully assigned proton n.m.r. spectrum

With its simple glycotripeptide repeating structure the antifreeze glycoprotein of polar fish may be an especially simple conformational mode for mucin glycoproteins with similar but more complex structures. The fully assigned proton n.m.r. spectrum confirms the anomeric configurations of the hexapyranosidic sugars of the side chains and the coupling constants of the alpha GalNAc and the beta Gal residues show both to be in the expected 4C1 chair conformation. The assignment of a single resonance for each proton of the (Ala-Thr-Ala)n repeat unit coupled with the observation of long range nuclear Overhauser effects (n.O.e.) implies a three-fold repeating conformation. The resonances of the two alanines are distinct and can be assigned to their correct positions in the peptide sequence by n.O.e. observed at the amide proton resonances on saturation of the alpha proton signals. The amide proton coupling constants of all three peptide residues are similar and imply a limited range of peptide backbone torsion angles, phi CN. The large n.O.e. which has been observed between the amide proton and the alpha proton of the residue preceding it in the sequence implies large positive values for the peptide dihedral angle, psi CC. Limits are placed on possible values of side chain dihedral angles by the observation of the coupling constant between the alpha and beta protons of the threonyl residue. The observation of n.O.e. between the anomeric proton of GalNAc and the threonyl side chain protons gives information on the conformation of the alpha glycosidic linkage between the disaccharide and the peptide. n.O.e. observed between the protons of the beta glycosidic linkage indicates the conformation of the disaccharide and the large amide proton coupling constant of the GalNAc residue shows a trans proton relationship. The spectroscopically derived data have been combined with conformational energy calculations to give a conformational model for antifreeze glycoprotein in which the hydrophobic surfaces of the disaccharide side chains are wrapped closely against a three-fold left handed helical peptide backbone. The hydrophilic sides of the disaccharides are aligned so that they may bind to the ice crystal face, which is perpendicular to the fast growth axis inhibiting normal crystal growth.     

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 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.     

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 and molecular modeling studies on two glycopeptides from the carbohydrate-protein linkage region of connective tissue proteoglycans.

Complete 1H and 13C NMR assignments are reported for two glycopeptides representing the carbohydrate-protein linkage region of connective tissue proteoglycans. These glycopeptides are the octasaccharide hexapeptide, Ser(GlcpAbeta(1-->3) Galpbeta(1-->3)Galpbeta(1-->4)Xylpbeta)-Gly-Ser-Gly-Se r (GlcpAbeta(1-->3)Galpbeta(1-->3)Galpbeta(1-->4)Xylp beta)-Gly (1), and the tetrasaccharide dipeptide, Ser(GlcpAbeta(1-->3)Galpbeta(1-->3)Galpbeta(1-->4)X ylpbeta)-Gly (2). The vicinal coupling constant data show that the monosaccharide residues adopt4 C 1 chair conformations. Distance geometry/simulated annealing calculations using 2D NOESY derived distance constraints yielded a single family of structures for the tetrasaccharide moiety, with well defined interglycosidic linkage conformations. The straight phi torsion angles of the glycosidic C1'-O1 bonds showed a strict preference for the -sc range whereas the psi torsion angles (O1-Cn) exhibited dependence upon the interglycosidic linkage position (-ac for beta(1-->3) linkage, +ac for beta(1-->4) linkage). The predominant conformation about the glycopeptide bond is straight phi = -sc and psi = +ac. The presence of strong daN (i, i+1) NOE contacts, and the general absence of dNN (i, i+1) contacts (except for a weak Ser-5/Gly-6 dNN contact) and the dbN (i, i+1) contacts (except for Ser-1/Gly-2) in the ROESY spectrum, suggest that the backbone for 1 is predominantly in an extended conformation. A comparison of the ROESY data for 1 with those obtained from the unglycosylated hexapeptide (3) of the same sequence suggests that glycosylation has only a marginal influence on the backbone conformation of the hexapeptide.     

C-C-A-G-G-C-m5C-T-G-G. Helical fine structure, hydration, and comparison with C-C-A-G-G-C-C-T-G-G.

The self-complementary DNA duplex C-C-A-G-G-C-m5C-T-G-G has been refined against 1.75-A x-ray diffraction data to an R value of 17.4%. In the crystal of space group P6, 10-base pair DNA fragments with characteristic sequence-related fine structure stack end to end to form long antiparallel B-type double helices. As shown by a structure analysis at lower resolution (Heinemann, U., and Alings, C. (1991) EMBO J. 10, 35-43), the overall geometry of C-C-A-G-G-C-m5C-T-G-G is similar to that of the unmethylated analog C-C-A-G-G-C-C-T-G-G despite a different crystal environment. The present high resolution structure analysis permits a detailed comparison of the two duplexes and their hydration spheres. Helical parameters are significantly correlated between both molecules, with the exception of the base pair propeller. Sugar pucker and backbone torsion angles alpha, gamma, delta, and chi show similar mean values, but their individual values deviate significantly between duplexes. In contrast, torsion angles beta, epsilon, and zeta change along the strands of both duplexes in much the same way. The effect of single-site methylation on DNA conformation appears to be small and limited to the base pairs directly involved. Methylation tends to push base pairs toward the minor groove of the helix. A regular minor groove hydration pattern involves dual hydrogen bonding of water molecules to O-4' and base atoms of C-C-A-G-G-C-m5C-T-G-G.     

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.     

Conformational analysis of a dinucleotide photodimer with the aid of the genetic algorithm.

The solution structure of the photodimer cis,syn-dUp[]dT is derived with the aid of the genetic algorithm. The conformational space available for the molecule is sampled efficiently using the computer program DENISE and tested against a set of constraints available from nmr experiments. The dominant conformation in solution found with this approach can be described by the following combinations of sugar-phosphate backbone torsion angles: epsilon(t), zeta(t), alpha(+), beta(-ac), and gamma(t). The conformation of the sugars and glycosidic torsion angles are S type and syn, respectively. The cyclobutane ring and pyrimidines are puckered. In addition, other conformations that exist in equilibrium with the first are found. It is concluded that the cyclobutane-pyrimidine system is rigid, whereas the sugar-phosphate backbone is flexible. The solution structures are compared with the crystal structure of the strongly related cyano-ethyl ester of cis,syn-dTp[]dT.     

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

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