Prediction of DNA structure from sequence: a build-up technique

A build-up technique has been devised that permits prediction of DNA structure from sequence. No experimental information is employed other than the force field parameters. This strategy for dealing with the multiple minimum problem requires a supercomputer to make the necessary global searches. The number of energy minimization trials that were made for each of the 16 deoxydinucleoside monophosphate conformational building blocks of DNA was 1944. As a test case, the minimum energy conformations of d(GpC) and d(CpG) to 5.5 kcal/mole were then combined to generate energy-minimized structures for d(CpGpC). The number of trials that were made for d(CpGpC) was 3752. Minima for this single-stranded trimer to 15 kcal/mole were then employed to search for minimum energy conformations of the duplex d(CpGpC).d(GpCpG). The number of starting conformations that were utilized at this stage was 1514. The lowest energy duplex had a Z-II-DNA conformation, followed by a B-DNA form at 1.2 kcal/mole. The A- and Z-I-forms as well as many novel Watson-Crick base-paired structures were found at higher energy. Finally, energy-minimized structures of d(CG)6.d(CG)6 in Z-II and B-DNA conformations were computed using torsion angles from the analogous duplex trimer minima.     

Conformations of dinucleoside monophosphates in relation to duplex DNA structures

Mononucleotide conformations are important in understanding the structural aspects of nucleic acids and polynucleotides. In order to study the influence of stacking interactions between adjacent bases in a polynucleotide on the preferred conformations of mononucleotides, conformational energy calculations have been carried out on dinucleoside monophosphate fragments. Four base sequences—d(ApT), d(TpA), d(CpG), and d(GpC)— have been analyzed in the framework of helical structures. Flexibility of the furanose ring has been incorporated in the investigations. Energetically favored conformers of the four compounds correspond to a variety of left- and right-handed uniform helical structures, similar to those of the commonly observed polymorphous forms. Implications of these investigations on the further understanding of double-helical polynucleotide conformations are briefly discussed.     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

Minimum energy conformations of proline-containing helices.

Proline occurs frequently in transmembrane alpha-helices of transport and receptor proteins even though statistical surveys demonstrate the overwhelming preference of this residue for a non-alpha-helical, hydrophilic environment. As a result, membrane-buried proline has been proposed to be functionally important, with function arising from structural discontinuity or destabilization of the helix. Destabilization may occur by Pro-mediated conformational transitions between discrete states, and may be manifested in membrane protein systems through reversible processes such as channel opening and closing or signal transduction. In this study, computer modeling of a model transmembrane alpha-helix, (Ala)8-Leu-Pro-Phe-(Ala)8, in a medium of low polarity (dielectric = 2), is used to examine the occurrence and energetic accessibility of Pro-mediated conformational interconversions. Leu psi and chi 1, Pro psi, and Phe phi and chi 1 torsion angles were assigned random values so that a data base of 200 conformations for each of the cis and trans states was generated. The conformations were minimized and low-energy structures organized into families. This analysis demonstrated that the most populated lowest energy family is the Trans-I conformation, corresponding to proline in a kinked alpha-helix. Two additional trans structures, Trans-II and Trans-III, as well as a cis conformation, Cis-I, are also energetically competitive. Interconversions between the trans states could thus be mediated by changes at a single torsion angle, accompanied by minor local hydrogen-bonding rearrangements. This work substantiates that membrane-buried proline can provide the basis for conformational transitions between discrete alpha-helix-based structures in a nonpolar environment.     

Determining local conformational variations in DNA. Nuclear magnetic resonance structures of the DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC) generated using back-calculation of the nuclear Overhauser effect spectra, a distance geometry algorithm and constrained molecular dynamics

Two-dimensional nuclear magnetic resonance (n.m.r.) spectroscopy and a variety of computational techniques have been used to generate three-dimensional structures of the two DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC). The central six base-pairs in these two decamers contain all ten dinucleotide pairs in DNA and thus, represent a model system for investigating how the local structure of DNA varies with base sequence. Resonance assignments were made for the non-exchangeable base protons and most of the C-1'-C-4' sugar protons in both decamers. Three-dimensional structures were generated using a distance geometry algorithm and these initial structures were refined by optimizing the fit of back-calculated spectra against the experimental two-dimensional nuclear Overhauser effect (NOE) spectra. This back-calculation procedure consists of calculating NOE cross relaxation rates for a given structure by solution of the Bloch equations, and directly accounts for spin diffusion effects. Use of this refinement procedure eliminates some assumptions that have been invoked when generating structures of DNA oligomers from n.m.r. data. Constrained energy minimization and constrained quenched molecular dynamics calculation were also performed on both decamers to help generate energetically favorable structures consistent with the experimental data. Analysis of the local conformational parameters of helical twist, helical rise, propeller twist, displacement and the alpha, beta, gamma, epison and zeta backbone torsion angles in these structures shows that these parameters span a large range of values relative to the X-ray data of nucleic acids. However, the glycosidic and pseudorotation angles are quite well defined in these structures. The implications that these results have for determination of local structural variations of DNA in solution, such as those predicted by Callidine's rules, are discussed. Our results differ significantly from some previous studies on determining local conformations of nucleic acids and comparisons with these studies are made.     

Raman studies of nucleic acids. XII. Conformations of oligonucleotides and deuterated polynucleotides

Laser Raman spectra of the trinucleoside diphoshate ApApA and dinucleoside phosphates ApU, UpA, GpC, CpG, and GpU are reported and discussed. Assignments of conformationally sensitive frequencies are-facilitated by comparison with spectra reported here of poly(rA), poly(rC), and poly(rU) in deuterium oxide solutions. The significant spectral differences between ApU and UpA, and between GpC and CpG, reveal that the sequence isomers have nonidentical conformations in aqueous solution. In UpA at low temperature the bases are stacked and the backbone conformation is similar to that found in ordered polynucleotide structures and RNA. In ApU no base stacking can be detected and the backbone conformation differs from that found in UpA, both in the orientation of phosphodiester linkages and in the internal conformation of ribose. At the conditions employed neither ApU nor UpA exhibits base pairing in aqueous solutions. In both GpC and CpG the bases are stacked and the phosphodiester conformations are similar to those encountered for UpA and RNA. However, major differences between spectra of GpC and CpG indicate that the geometries of stacking and ribosyl conformations are different. In GpC the Raman data favor the formation of hydrogen bonded dimers containing GC pairs. Protonation of C in GpC is sufficient to eliminate the ordered conformation detected by Raman spectroscopy. Despite the ordered backbone conformation evident in GpU, this dinucleoside apparently contains neither stacked nor hydrogen bonded bases at the conditions employed here. The Raman data also confirm the stacking interactions in ApApA, poly(rA), and poly(rC) but suggest that the backbone conformation in poly(rC) differs qualitatively from that found in most ordered polynucleotide structures and is thermally more stable. The present results demonstrate the sensitivity of the Raman technique to sequence-related structural differences in oligonucleotides and provide additional spectra–structure correlations for future conformational studies of RNA by laser Raman spectroscopy.     

Role of the closing base pair for d(GCA) hairpin stability: free energy analysis and folding simulations

Hairpin loops belong to the most important structural motifs in folded nucleic acids. The d(GNA) sequence in DNA can form very stable trinucleotide hairpin loops depending, however, strongly on the closing base pair. Replica-exchange molecular dynamics (REMD) were employed to study hairpin folding of two DNA sequences, d(gcGCAgc) and d(cgGCAcg), with the same central loop motif but different closing base pairs starting from single-stranded structures. In both cases, conformations of the most populated conformational cluster at the lowest temperature showed close agreement with available experimental structures. For the loop sequence with the less stable G:C closing base pair, an alternative loop topology accumulated as second most populated conformational state indicating a possible loop structural heterogeneity. Comparative-free energy simulations on induced loop unfolding indicated higher stability of the loop with a C:G closing base pair by ~3 kcal mol(-1) (compared to a G:C closing base pair) in very good agreement with experiment. The comparative energetic analysis of sampled unfolded, intermediate and folded conformational states identified electrostatic and packing interactions as the main contributions to the closing base pair dependence of the d(GCA) loop stability.     

Structure and dynamics of phosphate linkages and sugars in an abasic hexaloop RNA hairpin

Hairpins containing hexaloops are well represented among the diverse conformations adopted by the RNA molecules. To investigate the intrinsic properties of a backbone submitted to a hexaloop fold, we present here a molecular dynamics study of an abasic hexaloop closed by an A-form 6 basepair stem. The analysis of the 23 ns trajectory made in explicit solvent shows that both the sugars and the torsion angles in the loop undergo numerous conformational transitions. The south sugars, although not in a majority, are the major actors of the loop stretching. The five torsion angles, epsilon, zeta, alpha, beta, and gamma, are unequally variable, and only zeta and alpha exhibit trimodal distributions. The analysis of the phosphate linkages in terms of epsilonzeta'-alpha'-beta'-gamma-combinations allows us to define five conformational families, each one composed of one major substate in equilibrium with several less populated ones. The transitions between the substates within a family follow specific pathways involving the angles epsilon, zeta, and alpha. Thus, this work reveals that the backbone conformational space is both reduced and ordered even in a hexaloop devoid of bases.     

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5′-GGGCGNNAGCCU-3′ sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction. Received: 20 January 1999 / Revised version: 4 June 1999 / Accepted: 10 June 1999     

Wrinkled DNA.

The B form of poly d(GC):poly d(GC) in orthorhombic microcrystallites in oriented fibers has a secondary structure in which a dinucleotide is the repeated motif rather than a mononucleotide as in standard, smooth B DNA. One set of nucleotides (probably GpC) has the same conformations as the smooth form but the alternate (CpG) nucleotides have a different conformation at C3'-O3'. This leads to a distinctive change in the orientation of the phosphate groups. Similar perturbations can be detected in other poly d(PuPy):poly d(PuPy) DNAs such as poly d(IC):poly d(IC) and poly d(AT):poly d(AT) in their D forms which have tetragonal crystal environments. This suggests that such perturbations are intrinsic to all stretches of duplex DNA where purines and pyrimidines alternate and may play a role in the detection and exploitation of such sequences by regulatory proteins.     

Alpha/gamma transitions in the B-DNA backbone

In the crystal structures of protein complexes with B-DNA, α and γ DNA backbone torsion angles often exhibit non-canonical values. It is not known if these alternative backbone conformations are easily accessible in solution and can contribute to the specific recognition of DNA by proteins. We have analysed the coupled transition of the α and γ torsion angles within the central GpC step of a B-DNA dodecamer by computer simulations. Five stable or metastable non-canonical α/γ sub-states are found. The most favourable pathway from the canonical α/γ structure to any unusual form involves a counter-rotation of α and γ, via the trans conformation. However, the corresponding free energy indicates that spontaneous flipping of the torsions is improbable in free B-DNA. This is supported by an analysis of the available high resolution crystallographic structures showing that unusual α/γ states are only encountered in B-DNA complexed to proteins. An analysis of the structural consequences of α/γ transitions shows that the non-canonical backbone geometry influences essentially the roll and twist values and reduces the equilibrium dispersion of structural parameters. Our results support the hypothesis that unusual α/γ backbones arise during protein–DNA complexation, assisting the fine structural adjustments between the two partners and playing a role in the overall complexation free energy.     

Efficient Search on Energy Minima for Structure Prediction of Nucleic Acid Motifs

Abstract

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during con-formational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

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.     

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.     

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.     

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.     

The post-SCF quantum chemistry characteristics of inter- and intra-strand stacking interactions in d(CpG) and d(GpC) steps found in B-DNA, A-DNA and Z-DNA crystals

The energies of intra- and inter-strand stacking interactions in model d(GpC) and d(CpG) two-base-pair steps were estimated by MP2/aug-cc-pVDZ single point calculations corrected for basis superposition errors. The stacked two-nucleobase pairs were constructed using experimental values of base pair and base step parameters taken from Nucleic Acid Database (). Three distinct polymorphic forms were analysed, namely A-, B- and Z-DNA. The applied methodology enables statistical analysis of structural and energetic diversities. The structural relationships between polymorphic forms are quite complex and depend on the sequence of pairs. The variability of parameters such as shift and tilt is almost the same irrespective of the polymorphic form and sequence of steps analysed. In contrast, shift and twist distributions easily discriminate all three polymorphic forms of DNA. Interestingly, despite significant structural diversities, the energies of the most frequent energy ranges are comparable irrespective of the polymorphic form and base sequence. There was observed compensation of inter- and intra-strand interactions, especially for d(GpC) and d(CpG) steps found in A- and B-DNA. Thus, among many other roles, these pairs act as a kind of energetic buffer, balancing the double helix. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conformational Studies of Nucleic Acids: IV. The Conformational Energetics of Oligonucleotides: d(ApApApA) and ApApApA

Abstract

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: δ (C5′-C4′-C3′-03′)χ (04′-C1′-N9- C4), where the energetic basis for an approximately linear δ-χ correlation can be seen; ζ (C3′- 03′-P-05′)-δ, in which the experimentally observed linear correlation between ζ and δ in DNA (220° < ζ <280°) is clearly predicted; ζ-ε (C4′-C3′-03′-P), which shows that e increases with decreasing ζ <260°; α (03′-P-05′-C5′)-γ (05′-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 Ô is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in a and y. 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.