Spatial configurations of polynucleotide chains. I. Steric interactions in polyribonucleotides: a virtual bond model

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.

Theoretical studies of the structure and energies of base-paired nucleotides and the dissociation kinetics of a proflavine-dinucleotide complex

We presented calculations of base-paired dinucleoside phosphates and hexanucleoside pentaphosphates of varying compositions. Complete energy minimizations were performed for (a) the ten base-pair combinations of dinucleoside phosphates, starting from a B-DNA conformation, (b) six hexanucleoside pentaphosphates--base-paired CGCGCG, GCGCGC, G6-C6, TATATA, ATATAT, and A6-T6--starting with a B-DNA geometry, and (c) the four hexanucleoside pentaphosphates that have alternating pyrimidine-purine sequences, starting with a Z-DNA geometry. In addition, we studied the proflavine-base-paired CpG complex, using both complete energy minimization and energetic constraints to force the drug to dissociate from the dinucleoside phosphate. In many of these calculations, we examined the dependence of the calculated energies and structures on the potential function, focusing mainly on the effect of nonbonded potentials, the effective dielectric constant, and the role of counterions. These calculations allow us to explain why pur-(3',5')-pyr sequence isomers are more stable than pyr-(3'-5')-pur isomers. Both base-base and base-backbone energies are important in this differentiation, with the former being mainly van der Waals attraction and the latter mainly electrostatic energies. The calculations also allow us to understand the differences in double helical stabilities found by Wells et al. These differences, caused by electrostatic interactions between those bases not Watson-Crick hydrogen bonded, allow us to explain the following experimental data: poly(dG-dC) melts 12 degrees C higher than poly dG-poly dC, poly(dA-dT) melts 6 degrees C lower than poly dA-poly dT, and poly(dA-dG)-poly(dC-dT) melts 6 degrees C lower than poly(dA-dC)-poly(dT-dG). These results have interesting implications for drug binding: they imply that simple intercalators, such as ethidium, will exhibit a greater affinity for hetero- than for homopolymers and that this preference will be greater in the AT polymers than it is in the GC polymers. Our calculations allow us to explain the fact that Z-DNA is more stable than B-DNA under high salt conditions and to suggest some sequence dependence for the Z to B transition. We found that the activation energy for proflavine dissociating from dCpG is almost equal to the dissociation energy.

A model for the single stranded random coil form of polydeoxyadenylic acid from minimum energy conformations of the dimeric subunit

The minimum energy conformations of dApdA have been examined for their suitability as buildings blocks of the single stranded coil form of polynucleotides. Calculations of the characteristic ratio C difference = less than ro greater than 2/n liter2 were made for a polymer generated from all the low energy conformers, as well as for selected combinations. A polymer composed of a conformer with omega', omega = t*,g+,(skewed) psi = t, C-(2)-endo type pucker, in combination with the 'B' form, has a C difference equal to that observed in coils of apurinic acid (6) when the fraction of 'B' form conformers is approximately 25% and approximately 91%. The t*,g+ conformer is the second lowest energy form in the C-(2)-endo puckering domain, following the 'B' form.

RNA structure

This review is concerned primarily with the physical structure and changes in the structure of RNA molecules. It will be evident that we have not attempted comprehensive coverage of what amounts to a vast literature. We have tried to stay away from particular areas that have been recently reviewed elsewhere. Citations to and information from them are included, however, so that access to the literature is available. Much of what we treat in depth deals with the crystal structures and solution behaviour of model RNA compounds, including synthetic polymers and molecular fragments such as dinucleoside phosphates. Sequence data on natural RNA are cited, but not in detail. Similarly, apart from tRNA, natural RNAs the structural determinations of which are presently not so far advanced, are not dwelt upon. We have tried to present in detail the available structural data with scaled drawings that permit facile comparisons of molecular geometries.

Conformational energy calculations for dinucleotide molecules. A study of the component mononucleotide adenosine 3'-monophosphate

Semi-empirical conformational energy calculations were performed for the mononucleotides 5'-AMP, NMN+ and 3'-AMP. Only intramolecular forces are considered. Essentially all conformational states were explored to investigate the population distribution likely to be found in a non-crystalline environment. The calculations suggest that 5'-AMP and 3'-AMP are relatively flexible and a mixture of conformational states is expected. In contrast, the results for NMN+ suggest that a strong electrostatic interaction between the positively charged nicotinamide nitrogen atom and negatively charged phosphate oxygen is possible, stabilizing a few specific states. This interaction will be most significant in a solvent-free situation or an apolar environment.