The influence of terminal 3', 5' phosphates on conformations of dApdA

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 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) 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(CXnG) 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(6) 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.

Conformational characteristics of mixed sugar puckered deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp from energy minimization studies

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.

Base sequence effects in double helical DNA. I. Potential energy estimates of local base morphology

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.

Visualization of an AAF induced frameshift mutation: molecular views of base displacement in B-DNA from minimized potential energy calculations

Energy minimized structures of base displacement in an AAF modified B-DNA dodecamer are presented. A rational search strategy, beginning with a global search of the conformation space of the modified deoxydinucleoside monophosphate, together with model building by computer graphics, has been employed. A number of different minimum energy conformations have been located which reveal base displaced structures. These show fluorene interstrand stacking, fluorene inter- and intrastrand stacking, and non-stacked fluorene situated in the denatured bulge. The local helix axis is bent to various extents in the different forms, and one or two base pairs are fully denatured. One structure of special interest offers a molecular view that suggests how AAF can induce the -2 deletion mutation observed in AAF modified E. coli.

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 studies of nucleic acids. II. The conformational energetics of commonly occurring nucleosides

We have examined the conformational energetics of the eight most commonly occurring nucleosides--A, U, G, C, dA, dT, dG, dC--as monitored by a semi-empirical energy force field. These are the first reported calculations to completely explore the entire conformational spaces available to all eight major nucleosides using experimentally consistent furanose geometries and an appropriate force field. Central to our approach is the ability to model an experimentally reasonable furanose for each nucleoside directly from only one parameter, the phase angle of pseudorotation P, as described in the previous paper (D.A. Pearlman, and S.-H. Kim, preceeding paper in this issue). This allows us to specify the conformation of a nucleoside by three variables: torsion angle gamma (O5'-C5'-C4'-C3'); torsion angle chi (O4'-C1'-N9/N1-C4/C2); and P. In our study each of these parameters was allowed to vary independently and in small increments over the range 0-360 degrees. The empirically observed preferences for C3'-endo and C2'-endo sugar conformations, for anti and syn values of chi and for staggered (g+, t, g-) values of gamma can be explained on the basis of the energy maps so obtained. Finer details, such as the different conformational preferences of ribonucleosides and deoxyribonucleosides and of purines and pyrimidines, can also be extracted from these maps and are consistent with experiment. The calculations support previous descriptions of pseudorotation as hindered. Statistical Boltzmann population factors for different conformational ranges in gamma, chi, and P, as predicted by the calculations, are consistent with factors obtained from crystallographic data. The excellent results here provide additional support for the suitability of the new sugar modeling technique used.

Backbone conformation in nucleic acids: an analysis of local helicity through heminucleotide scheme and a proposal for a unified conformational plot

A relationship has been established to express the local helicity of a polynucleotide backbone directly in terms of the virtual bonds spanning the conformationally equivalent heminucleotide repeats, with a view to provide a better understanding of the cumulative effects of all the chemical bond rotational variations on local helicity. Using this, an analysis made with a few oligodeoxynucleotide crystal structures clearly brings forth that it is the concerted movements manifested in the near neighbour correlations between the pair of chemical bonds C4'-C5' and P-O5' and C4'-C3' and P-O3' of the 5' and 3' heminucleotides respectively that are primarily responsible for the observed non-uniform helical twists both in A and B type helical backbones. That these need not be restricted to oligodeoxynucleotides but may be a feature of oligoribonucleotides backbone also is shown from an analysis of helical segments of yeast tRNA(Phe). A proposal of a unified or a grand two dimensional conformational plot which would help visualise succinctly the overall effect of the variations in all the repeating six chemical bonds of a polynucleotide backbone is made. Apart from considerable simplification, the plot affords identification on it regions characteristic of helical, and loop and bend conformations of nucleic acid backbone chain.

Conformational characteristics of the dinucleoside triphosphate pCpGp from energy-minimization studies

The influence of the 3'- and 5'- terminal phosphates on the conformational characteristics of the dinucleoside monophosphate CpG is described in this paper. The computed potential energy of the system is minimized with respect to the relevant 10 dihedral angles permitting the two sugar rings to adopt the alternative puckering states, 2E and 3E. Of the 84 conformations considered, 22 become energetically accessible. The familiar A-, B-, Z- and Watson-Crick-type backbone states of DNA subunits become low-energy forms for this RNA unit pCpGp also. The Watson-Crick-type backbone is invariably preferred in all the four sugar pucker sequences, indicating its importance in the dynamics of sugar pucker fluctuations and in the DNA-RNA association. The interphosphate geometries and the possible hydrogen-bonded states are discussed in relation to the varied folded/extended polynucleotide structures.

Conformation of the deoxydinucleoside monophosphate dCpdG modified at carbon 8 of guanine with 2-(acetylamino)fluorene

Minimized conformational potential energy calculations were performed for dCpdG modified with the carcinogen 2-(acetylamino)fluorene (AAF). The major adduct, linked via a covalent bond between guanine C-8 and N-2 of AAF, was investigated. The 12 variable torsion angles and both deoxyribose puckers were independent flexible parameters in the energy minimizations. Three categories of low-energy conformers were calculated in which the guanine was syn and nearly perpendicular to the plane of the fluorene: (1) forms in which fluorene is stacked with cytidine (included among these is the global minimum energy conformation); (2) conformers which preserve guanine-cytidine stacking while placing the fluorene in a base-pair obstructing position; (3) conformers which maintain guanine-cytidine stacking and place the fluorene at the helix exterior, without interfering with base pairing. The Z form is important in this group. In addition, a low-energy conformation with guanine anti, but still nearly perpendicular to fluorene, was computed. Molecular models were constructed for the most important conformations incorporated into larger polymers. These indicated that the fluorene-cytidine stacked forms induce a severe kink in the B helix. Conformers with guanine-cytidine stacking and AAF in a base-pair obstructing position place the AAF at the B-type helix interior with little distortion in the helix direction. Conformers with the guanine-cytidine stack in which AAF does not affect base pairing place the fluorene at the Z or alternate helix exterior. It is suggested that base sequence, extent of modification, and external conditions such as salt concentration determine which of a number of possible conformational effects is actually induced by AAF. The variety of observed experimental results with AAF-modified DNA may reflect there various conformational possibilities.

Computational studies of polynucleotide flexibility

Details of polynucleotide flexibility may be probed through a combination of semiempirical potential energy calculations and statistical mechanical analyses. The pseudorotational motions of the furanose and the long-range correlated rotations of the chain backbone are described briefly here.

Conformational characteristics of dimeric subunits of RNA from energy minimization studies. Mixed sugar-puckered ApG, ApU, CpG, and CpU

Following the procedure described in the preceding article, the low energy conformations located for the four dimeric subunits of RNA, ApG, ApU, CpG, and CpU are presented. The A-RNA type and Watson-Crick type helical conformations and a number of different kinds of loop promoting ones were identified as low energy in all the units. The 3E-3E and 3E-2E pucker sequences are found to be more or less equally preferred; the 2E-2E sequence is occasionally preferred, while the 2E-3E is highly prohibited in all the units. A conformation similar to the one observed in the drug-dinucleoside monophosphate complex crystals becomes a low energy case only for the CpG unit. The low energy conformations obtained for the four model units were used to assess the stability of the conformational states of the dinucleotide segments in the four crystal models of the tRNAPhe molecule. Information on the occurrence of the less preferred sugar-pucker sequences in the various loop regions in the tRNAPhe molecule has been obtained. A detailed comparison of the conformational characteristics of DNA and RNA subunits at the dimeric level is presented on the basis of the results.

Conformational characteristics of the dimeric subunits of DNA from energy minimization studies. Mixed sugar-puckered dApdA, dApdT, dTpdA, and dTpdT

An extensive investigation on the conformational characteristics of four deoxydinucleoside monophosphates, namely, dApdA, dApdT, dTpdA, and dTpdT was carried through calculation of the classical potential energy of the systems. The four major types of sugar-pucker sequences, namely, 3E-3E, 3E-2E, 2E-3E, 2E-2E, were included in the study. For each of the units, energies were computed for 96 starting conformations that resulted from the consideration of all possible low energy regions for the relevant seven dihedral angles and the four sugar-pucker sequences, and minimized by permitting all the seven dihedral angles to vary simultaneously. The number and the order of preference of low energy conformations obtained were found to be characteristic of the base sequence of the unit considered. The conformational states close to the A-DNA, B-DNA, C-DNA, and Watson-Crick DNA structures are noted to be preferred for all the units except dTpdT. The 3E-2E sugar-pucker sequence is the most favored and the 2E-3E sequence is the least favored state in terms of the associated number of local minima. For each unit, there exists a set of specific conformational states with more or less equal stabilities but different sugar-pucker sequences. The mixed sugar-pucker states 2E-3E and 3E-2E, when incorporated, in the conventional A-DNA and B-DNA conformational states, respectively, have energies that allow them to act as intermediates in the B form in equilibrium with A form transitions. Such transitions are most likely to occur at sites with a Thymine-Adenine base sequence. Available experimental results were interpreted in terms of their stabilities.

Yeast tRNAPhe conformation wheels: a novel probe of the monoclinic and orthorhombic models

A series of conformation wheels is constructed from the recently refined X-ray crystallographic data of monoclinic and orthorhombic yeast tRNAPhe. These circular plots relate the primary chemical structure (i.e., base sequence) directly to the secondary and tertiary structure of the molecule. The circular sequence of backbone torsion angles displays a unique pattern that is useful both in distinguishing the ordered and disordered regions of the molecule and in comparing the three sets of experimental data. Composite conformation wheels describe the fluctuations in the "fixed" parameters (phi', phi, chi) and independent conformation wheels reveal the changes in the "variable" parameters (omega', omega, psi, psi') of the three different yeast tRNAPhe models. Additional plots of base-stacking parameters help to visualize the intimate interrelationship between chemical sequence and three-dimensional folding of yeast tRNAPhe. The composite data illustrate several conformational schemes that position the bases of adjacent nucleosides in a parallel stacked array and reveal an even larger number of conformations that introduce bends or turns in the polynucleotide chain.