Jun 1, 1974 - Department of Biochemistry, College of Agricultural and Life Sciences, University of Wisconsin, Madison, Wisc. 53706. Communicated by ...
Proc. Nat. Acad. Sci. USA Vol. 71, No. 9, pp. 3325-3328, September 1974
Backbone Conformations in Secondary and Tertiary Structural Units of Nucleic Acids. Constraint in the Phosphodiester Conformation* (energy calculation/polynucleotide backbone/constrained phosphodiester conformation/ dinucleoside mono-, di-, and triphosphates)
N. YATHINDRA AND M. SUNDARALINGAM Department of Biochemistry, College of Agricultural and Life Sciences, University of Wisconsin, Madison, Wisc. 53706
Communicated by Nelson J. Leonard, June 1, 1974 ABSTRACT The possible backbone phosphodiester conformations in a dinucleoside monophosphate and a dinucleoside triphosphate have been investigated by semiempirical energy calculations. Conformational energies have been computed as a function of the rotations w' and cu about the internucleotide P.0(3') and P-0(5') linkages, with the nucleotide residues themselves assumed to be in one of the preferred [C(3')-endoj conformations. The terminal phosphates in a dinueleoside triphosphate greatly limit the possible conformations for the backbone (in a polynucleotide) compared to a dinucleoside monophosphate. There appear to be two major types of conformations that are favored for the backbone. The phosphodiester conformation (,',,) - (2900,2900) characteristic of helical structures is one of them, indicating that the polynucleotide backbone shows an inherent tendency for the helical conformation. The other favored conformation is centered at (c', ) (1900,3000) and results in an extended backbone structure with unstacked bases. A third possible conformation centered at (ca', w) - (2000, 600) and the (1900, 3000) conformation appear to be important for the folding of a polynucleotide chain. The conformation (V',c) (800,800), observed in a dinucleoside monophosphate and believed to be a candidate for producing an abrupt turn in a polynucleotide chain, is found to be stereochemically unfavorable in a dinucleoside triphosphate and a polynucleotide. -
and a dinucleoside monophosphate has provided considerable insight into the folding of a polynucleotide chain in nucleic acids. METHODS Conformational energies for a dinucleoside monophosphate and- a dinucleoside triphosphate have been computed as a
function of the rotations a' and X about the P-0(3') and P-0(5') bonds (Fig. 1) by using partitioned potential functions that take into account the van der Waals, electrostatic, and torsional potentials. The van der Waals energies have been evaluated by use of the 6-12 Lennard-Jones function with the parameters reported earlier (15, 16). Electrostatic interactions have been estimated in the monopole approximation, with the charges given by Renugopalakrishnan et al.
Considerable attention is currently being focused on the backbone conformations in the tertiary folds and loops of nucleic acids such as transfer RNA. It has been found that the nucleotide residues in a polynucleotide exhibit essentially two favored conformations differing in the sugar pucker, and the major flexibility in a polynucleotide backbone resides primarily in the internucleotide P-O links (1-3). Results from both x-ray and potential energy calculations have indicated that there can be several possible conformations about the P-O bonds for a dinucleoside monophosphate (1, 3-14). In continuation of our studies on the conformations of nucleic acids and their structural units, we recently reported (15) results of our semiempirical energy calculations on mononucleotides. Here we present results of similar calculations on a dinucleoside monophosphate and a dinucleoside triphosphate. It is shown that in contrast to a dinucleoside monophosphate, the presence of the terminal 3'- and 5'-phosphate groups greatly restricts the number of possible phosphodiester conformations in a dinucleoside triphosphate. This distinction in the conformational properties of a dinucleoside triphosphate FIG. 1. Section of a polynucleotide chain showing the backbone atoms. The notations and.definitions used for the torsion angles are the same as given in ref. 1.
This work was presented at the ACS Meeting, Chicago, Aug. 26-31, 1973 and the ACA Meeting, Berkeley, Mar. 24-28, 1974. *
3325
3326
Proc. Nat. Acad. Sci. USA 71
Biochemistry: Yathindra and Sundaralingam
(1974)
TABLE 1. Conformational angles for the nucleotides*
C(3')-endo 600 1800 2100
Sugar pucker * 4.'
* These angles show a range of 15-30'.
sugars in addition to varying the nucleotide torsion angles over the observed range (1). Possible backbone conformations for a dinucleoside monophosphate The coiformational energy map obtained for the rotations
about the adjacent P-O(3') (w') and P-0(5') (X) bonds in a dinucleoside monophosphate is shown in Fig. 2. The global (600,1700) corresponding to the (g+,t) minimum at (W',W) (800,800) = conformation and the minimum at (w',) (g+,g+) have nearly the same energies. The latter conformation has been observed in the x-ray structures of the dinucleoside monophosphate, uridylyl (3',5')-adenosine, UpA (4,5), and in the A2pA3 fragment of the trinucleoside diphosphate, adenylyl-(3',5')-adenylyl-(3',5')-adenosine A'pA2pA3 (9). The low energy domain at (WIco,) (2900,2900) = (g-, g-) is characteristic of helical structures and has been observed in the dinucleoside monophosphates, adenylyl-(3',5')uridine ApU (8), guanylyl-(3',5')-cytidine GpC (7), and in the A'pA*2 fragment of the trinucleoside diphosphate, A pA2pA8 (9). The only other conformation (t,g-) that has been observed corresponds to the domain (w',w) = (1900, 3000) and is exhibited by the second molecule of UpA (4, 5) and thymidyl-(5',3')-thymidylate-5' (pTpT) (10). It is worthwhile to -
-
180 -_. W
FiG. 2. Equienergy contours for a dinucleoside monophosphate drawn as a function of the rotations w' and w about the P-0(3') and P-0(5') bonds, respectively. The conformations along the lines 600, 1800, and 300° are referred to as gauche (+) (g+), trans (t), and gauche (-), (gj).
(17). A 3-fold symmetric potential with a barrier height of 1 kcal/mole has been assumed for P-O bond torsions. In this paper the calculations are restricted to nucleotides with C(3')-endo sugars, and the torsion angles used for the nucleotide units are given in Table 1. These angles are similar to the values observed for nucleotide residues in di- and oligonucleoside structures (4, 5, 7-10). Calculations are also being carried out for dinucleoside phosphates comprising C(2')-endo sugars and alternate C(3')-endo and C(2')-endo
-
0 -
FIG. 3. Equienergy contours for a dinucleoside triphosphate drawn as a function of the rotations co' and co about the P-0(3') and P-O(5') bonds, respectively.
(g 99)
FIG. 4. Perspective view of the conformation (w',cw) = (800,800) (g+,g+), which produces the sharpest turn in the backbone of a dinucleoside monophosphate, but appears to be unfavorable for a polynucleotide (see text), as illustrated with the molecule UpA.
Proc. Nat. Acad. Sci. USA 71
(1974)
Constraint in Phosphodiester Conformations
b
a
g9-w9-)
v
3327
C
(twg
(tg)
FIG. 5. Perspective views of the three favored backbone conformations for a dinucleoside triphosphate and a polynucleotide. (a) (w',c) (2900,2900) (gPgj) is the most favored backbone conformation for a helix and it can also be a part of loop structures; (b) (c',o) (1900,3000) (t,g-) generates extended helix-like conformation for the backbone and appears to be important in single-stranded regions of nucleic acids; (c) (C',w) (200°,60°) (tg+) appears to be the conformation favored for the folding of a polynucleotide chain in hairpin loops. Note that the terminal phosphates are not shown in the above figures. -
-
mention that in the latter (pTpT), the only known structure of a deoxydinucleotide, the sugar rings exhibit the C(2')-endo pucker, while in the former they are C(3')-endo. In addition to the above four conformations there are a few other probable conformations for a dinucleoside monophosphate with (w',co) values around (3000,1800), (1900,600), and (800,2400). It may be noted that there is one other minimum (1800,1800) (t,t). This correin the energy map at (co',w) sponds to the fully extended conformation for the phosphodiester group and is likely to be less favored due to eclipsed interaction of the lone pair electrons on the adjacent ester oxygens (1, 3). In the semiempirical calculations, the effects due to the lone pair electrons and their conjugation with the phosphorous d orbitals are not accounted for, and consideration of these effects may reduce the importance of the (t,t) conformation. It is interesting that ab initio calculations on dimethyl phosphate (although not an exact model of a dinucleoside monophosphate) show that the fully extended conformation (t,t) turns out to be a high energy conformation
(18, 19). Possible backbone conformations for a dinucleoside triphosphate and a polynucleotide. Role of the phosphates
The above studies on dinucleoside monophosphates (see also refs. 1, 3, 6, and 11-14) show that the rotations about the ester P4) bonds are quite flexible and suggest that there can be several backbone conformations for a polynucleotide. This situation would have made the conformational analysis of polynucleotides quite difficult. But in a polynuoleotide, unlike the situation in a dinucleoside monophosphate, the 3'- and 5'-hydroxyl groups are phosphorylated. It was anticipated that the terminal phosphates of a dinucleoside triphosphate would render some of the folded conformations for the back-
bone sterically unfavorable. Consequently, calculations were carried out for a dinucleoside triphosphate similar to those performed for the dinucleoside monophosphate. The (',co) map thus obtained is shown in Fig. 3. It is seen that the number of possible conformations for the dinucleoside triphosphate is strikingly limited, in sharp contrast to the dinucleoside monophosphate (Fig. 2). The minima corresponding to the folded backbone conformations (co',w) (800,800) (g+, g+) and (600,1700) (g+,t) are conspicuously absent in Fig. 3 due to steric and electrostatic interactions introduced by the terminal phosphates. These interactions also render the region near (w',w) (800,2400) a higher energy conformation. Thus, the conformations for the range O'00 to 800 are expected to be energetically unfavorable in dinucleoside triphosphates and polynucleotides. As mentioned above, the fully extended conformation (1800,1800) (t,t) is also not likely to be favored (1,3). There appear to be two predominant conformations for the dinucleoside triphosphate pXpXp. The conformation centered at (w',w) (2900,2900) (g-,g-) (Fig. 5a) corresponding to the global minimum is characteristic of all known double helical nucleic acids and polynucleotides (1, 20) and appears to be the most favored conformation for the helical backbone. It is interesting that despite the omission of both the base stacking and base pairing energies, the helical conformation turns out to be a favored conformation. This result substantiated an earlier observation (1) that the polynucleotide backbone shows an inherent tendency for the helical conformation and that base stacking and base pairing interactions provide additional stabilization. The other favored conformation is centered at (co',co) (1900,300°) (t,g-) and corresponds to an unstacked conformation with an extended backbone (Fig. 5b). The above results indicate that in solution, polynucleo-
3328
Biochemistry: Yathindra and Sundaralingam
tides are predominantly composed of a mixture of the stacked (ggj) and unstacked (t,g-) conformations. The relative populations of these two conformations will be expected to depend on the nature of the bases and their sequence. The conformation centered at (co',c) (200°,600) (tg+) (Fig. 5c), which is considerably less favored compared to the above two, turns out to be a possible candidate for effecting a hairpin turn in a polynucleotide chain. In addition to the above three conformations, the (2900,1900) (gt) and (1000,2700) (g+,g) conformations may also be important in regulating the backbone conformations of single-stranded regions of nucleic acids. The (g+,g+) conformation
Although the conformation around (w',w) = (800,80°) (g+, g+) is energetically unfavorable within a polynucleotide chain, it can be entertained by the terminal dinucleoside diphosphate residues in a nucleic acid when stabilized by an intramolecular hydrogen bond involving the free hydroxyl group 0(3')-H or 0(5')-H and the corresponding penultimate phosphate at the 3'- or 5-end, respectively. In fact, this conformation has been observed in the trinucleoside diphosphate A'pA2pA' (9), which exhibits an intramolecular hydrogen bond between the hydroxyl group 0(3')-H and the distal phosphate group. By analogy, a similar intramolecular hydrogen bond can stabilize the (g+,g+) conformation in a dinucleoside diphosphate. Therefore, the (g+,g+) backbone conformation seems to be restricted to the dinucleoside monophosphates, dinucleoside diphosphates, and the terminal dinucleoside diphosphate moieties of nucleic acids (Fig 4). CONCLUSIONS
The main result that has emerged from this work is that the rotations about the internucleotide P-0 bonds in a dinucleoside triphosphate and a polynucleotide are considerably more restricted than in a dinucleoside monophosphate. It must be stressed that the P-0 bonds still show greater rotational flexibility (besides the sugar pucker) than the other backbone bonds C(4')-C(5'), C(5')-O(59, and C(3')-(3'). It is found that certain phosphodiester conformations that are likely in dinucleoside monophosphates turn out to be energetically unfavorable in a polynucleotide, suggesting that extrapolation of the results derived from dinucleoside monophosphates for the analysis of polynucleotide conformation (5()0,800) can be misleading. The conformation (w',w) (g+,g+), which produces the sharpest turn in a dinucleoside monophosphate, is found to be stereochemically unfavorable for a dinucleoside triphosphate and a polynucleotide. It appears that the phosphodiester conformations (t,g+) and (grt) are the most probable candidates for the folding of a polynucleotide chain. The distinction between the conforma-
Proc. Nat. Acad. Sci. USA 71
(1974)
tional properties of a dinucleoside triphosphate and a dinucleoside monophosphate is reminiscent of the conformational rigidity found in a nucleotide in contrast to a nucleoside (3). The presence of the phosphate group in a nucleotide strikingly restricts the conformations about the backbone C(4')-C(5'), C(5')-O(59, and C(3')-O(3') bonds and the side-chain glycosyl C(19)-N bond (15). Therefore, the phosphate groups play a major role in limiting not only the conformations of the nucleotide building blocks themselves, but also their relative orientations about the internucleotide P-O linkages in a polynucleotide. The observations that the terminal phosphates in a dinucleoside triphosphate greatly limit the possible conformations for the phosphodiester group, taken in conjunction with the rigid nucleotide concept, further simplifies the conformational analysis of polynucleotides. We gratefully thank the National Institutes of Health of the United States Public Health Service for support of this work through Grants GM-17378 and GM-18455. 1. Sundaralingam, M. (1969) Biopolymers 7, 821-869. 2. Rao, S. T. & Sundaralingam, M. (1970) J. Amer. Chem. Soc. 92, 4963-4970. 3. Sundaralingam, M. (1973) in Conformations of Biological Molecules and Polymers, The Jerusalem Symposia on Quantum Chemistry and Biochemistry, Israel Academy of Sciences and Humanities, Jerusalem, 5, 417-453. 4. Rubin, J., Brennan, T. & Sundaralingam, M. (1972) Biochemistry 11, 3112-3127. 5. Sussman, J. L., Seeman, N. C., Kim, S. H. & Berman, H. M. (1972) J. Mol. Biol. 66, 403-421. 6. Kim, S. H., Berman, H. M., Seeman, N. C. & Newton, M. D. (1973) Acta CrystaUogr. B29, 703-710 7. Day, R. O., Seeman, N. C., Rosenberg, J. M. & Rich, A. (1973) J. Mot. Biol. 70, 849-53. 8. Rosenberg, J. M., Seeman, N. C., Kim, J. J. P., Suddath, F. L., Nicholas, H. B. & Rich, A. (1973) Nature 243, 150-154. 9. Suck, D., Manor, P. C., Germain, G., Schwalbe, C. H., Weimann, G. & Saenger, W. (1973) Nature New Biol. 246, 161-165. 10. Camerman, N., Fawcett, J. K. & Camerman, A. (1973) Science 182, 1142-1143. 11. Saran, A. & Govil, G. (1971) J. Theor. Biol. 33, 407-418. 12. Pullman, B., Perahia, D. & Saran, A. (1972) Biochim. Biophys. Acta 269, 1-14. 13. Olson, W. K. & Flory, P. J. (1972) Biopolymers 11, 25-56. 14. Stellaman, S., D., Hingerty, B., Broyde, S. B., Subramanian, E., Sato, T. & Langridge, R. (1973) Biopolymers 12, 2731-2750. 15. Yathindra, N. & Sundaralingam, M. (1973) Biopolymers 12, 297-314; 2075-2082; 2261-2277. 16. Yathindra, N. & Sundaralingam, M. (1973) Biochim. Biophys. Acta 308, 17-21. 17. Renugopalakrishnan, V., Lakshminarayanan A. V. & Sasisekharan, V. (1971) Biopolymers 10, 1159-1167. 18. Newton, M. (1973) J. Amer. Chem. Soc. 95, 256-258. 19. Perahia, D., Pullman, B. & Saran, A. (1973) C.R.H. Acad. Sci. 277, 2257-2260. 20. Arnott, S. (1970) Progr. Biophys. Mol. Biol. 21, 267-319.