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.     

Conformational studies of nucleic acids: III. Empirical multiple correlation functions for nucleic acid torsion angles

There are seven significantly variable torsion angles in each monomer unit of a polynucleotide. Because of this, it is computationally infeasible to consider the energetics of all conformations available to a nucleic acid without the use of simplifications. In this paper, we develop functions suggested by and regression fit to crystallographic data which allow three of these torsion angles, alpha (O3'-P-O5'-C5'), delta (C5'-C4'-C3'-O3') and epsilon (C4'-C3'-O3'-P), to be calculated as dependent variables of those remaining. Using these functions, the seven independent torsions are reduced to four, a reduction in complexity sufficient to allow an examination of the global conformational energetics of a nucleic acid for the remaining independent torsion angles. These functions are the first to quantitatively relate a dependent nucleic acid torsion angle to several different independent angles. In all three cases the data are fit reasonably well, and in one case, alpha, the fit is exceptionally good, lending support for the suitability of the functions in conformational searches. In addition, an examination of the most significant terms in each of the correlation functions allows insight into the physical basis for the correlations.     

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.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

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.     

DFT studies of the disaccharide, alpha-maltose: relaxed isopotential maps

The disaccharide, alpha-maltose, forms the molecular basis for the analysis of the structure of starch, and determining the conformational energy landscape as the molecule oscillates around the glycosidic bonds is of importance. Thus, it is of interest to determine, using density functionals and a medium size basis set, a relaxed isopotential contour map plotted as a function of the phi(H) and psi(H) dihedral angles. The technical aspects include the method of choosing the starting conformations, the choice of scanning step size, the method of constraining the specific dihedral angles, and the fitting of data to obtain well defined contour maps. Maps were calculated at the B3LYP/6-31+G( *) level of theory in 5 degrees intervals around the (phi(H),psi(H))=(0 degrees ,0 degrees ) position, out to approximately +/-30 degrees or greater, for gg-gg'-c, gg-gg'-r, gt-gt'-c, gt-gt'-r, tg-tg'-c, and tg-tg'-r conformers, as well as one-split gg(c)-gg'(r) conformer. The results show that the preferred conformation of alpha-maltose in vacuo depends strongly upon the hydroxyl group orientations ('c'/'r'), but the energy landscape moving away from the minimum-energy position is generally shallow and transitions between conformational positions can occur without the addition of significant energy. Mapped deviations of selected parameters such as the dipole moment; the C1-O1-C4', H1-C1-O1, and H4'-C4'-O1 bond angles; and deviations in hydroxymethyl rotamers, O5-C5-C6-O6, O5'-C5'-C6'-O6', C5-C6-O6-H, and C5'-C6'-O6'-H', are presented. These allow visualization of the structural and energetic changes that occur upon rotation about the glycosidic bonds. Interactions across the bridge are visualized by deviations in H(O2)...O3', H(O3')...O2, and H1...H4' distances and the H(O2)-O2-C2-C1 and H'(O3')-O3'-C3'-C4' hydroxyl dihedral angles.     

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 3. The magnitude of torsional angle epsilon in d(TpA) from CCOP and HCOP NMR coupling constants.

Carbon-13 NMR spectra of the deoxyribonucleotide d(TpA), 3',5'-cyclic AMP and 3',5'-cyclic dAMP were measured. It is shown that the different substitution of C2' in deoxyribonucleotides versus ribonucleotides does not affect the vicinal C2'-C3'-O3'-P coupling to a measurable extent. Therefore, the same set of Karplus parameters may be used for the C2'-C3-O3'-P couplings in ribonucleotides and in deoxyribonucleotides. Vicinal carbon-phosphorus and proton-phosphorus coupling constants are used to calculate the magnitude of the torsion angle epsilon (C4'-C3'-O3'-P), which amounts to 195(0) in the trans conformer and to 261(0) in the gauche(-) conformer.     

Vibrational analysis of nucleic acids. V. Force field and conformation-dependent modes of the phosphodiester backbone modeled by diethyl phosphate.

A generalized valence force field is derived for the diethyl phosphate anion [(CH3CH2O)2PO2-] and its deuterium [(CH3CD2O)2PO2-, (CD3CH2O)2PO2- and (CD3CD2O)2PO2-] and carbon-13 [(CH3 13CH2O)2PO2-] derivatives in the stable trans-gauche-gauche-trans conformation. Normal coordinate analysis of the trans-gauche-gauche-trans conformer, which serves as a structural analog of the nucleic acid phosphodiester group, is based on comprehensive infrared and Raman spectroscopic data and vibrational assignments obtained for the diethyl phosphate anion. The generalized valence force field is in good agreement with the scaled ab initio force field of diethyl phosphate and represents significant improvement over earlier modeling of the phosphodiester moiety with dimethyl phosphate. The conformational dependence of skeletal C-C-O-P(O2-)-O-C-C stretching vibrations is also explored. Starting with the trans-gauche-gauche-trans conformation, the frequency dependence of skeletal stretching modes has been obtained by stepwise rotation of the torsion angles of the P-O and C-O bonds corresponding to nucleic acid torsions alpha (P-O5'), beta (O5'-C5'), epsilon (C3'-O3'), and zeta (O3'-P). Both symmetric and antisymmetric phosphoester stretching modes are highly sensitive to P-O and C-O torsions, whereas symmetric and antisymmetric phosphodioxy (PO2-) stretching modes are less sensitive. The present results provide an improved structural basis for understanding previously developed empirical correlations between vibrational marker bands and nucleic acid backbone conformation.     

Flexibility of 3',5' deoxyribonucleoside diphosphates

In 3',5' deoxyribonucleoside diphosphates, in addition to the nature of the base and the sugar puckering, there are six single bond rotations. However, from the analysis of crystal structure data on the constituents of nucleic acids, only three rotational angles, that are about glycosyl bond, about C4'-C5' and about C3'-O3' bonds, are flexible. For a given sugar puckering and a base, potential energy calculations using non-bonded, electrostatic and torsional functions were carried out by varying the three torsion angles. The energies are represented as isopotential energy surfaces. Since the availability of the real-time color graphics, it is possible to analyse these isopotential energy surfaces. The calculations were carried out for C3' exo and C3' endo puckerings for deoxyribose and also for four bases. These calculations throw more light not only on the allowed regions for the three rotational angles but also on the relationships among them. The dependence of base and the puckering of the sugar on these rotational angles and thereby the flexibility of the 3',5' deoxyribonucleoside diphosphates is discussed. From our calculations, it is now possible to follow minimum energy path for interconversion among various conformers.     

Effect of distortions in the deoxyribose phosphate backbone conformation of duplex oligodeoxyribonucleotide dodecamers containing GT, GG, GA, AC, and GU base-pair mismatches on 31P NMR spectra

We have previously suggested that variations in the 31P chemical shifts of individual phosphates in duplex oligonucleotides are attributable to torsional angle changes in the deoxyribose phosphate backbone. This hypothesis is not directly supported by analysis of the 1H/31P two-dimensional J-resolved spectra of a number of mismatch dodecamer oligonucleotide duplexes including the following sequences: d-(CGTGAATTCGCG), d(CGUGAATTCGCG), d(CGGGAATTCGCG), d(CGAGAATTCGCG), and d(CGCGAATTCACG). The 31P NMR signals of the dodecamer mismatch duplexes were assigned by 2D 1H/31P pure absorption phase constant time (PAC) heteronuclear correlation spectra. From the assigned H3' and H4' signals, the 31P signals of the base-pair mismatch dodecamers were identified. JH3'-P coupling constants for each of the phosphates of the dodecamers were obtained from 1H/31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. JH3'-P coupling constants were measured for many of the oligonucleotides as a function of temperature. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. This correlation can be further extended to the C3'-O3'-P-O5' torsional angle (zeta) by using a linear relationship between epsilon and zeta obtained from crystal structure studies. The 31P chemical shifts follow the general observation that the more internally the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. In addition, 31P chemical shifts show sequence- and site-specific variations. Analysis of the backbone torsional angle variations from the coupling constant analysis has provided additional information regarding the origin of these variations in 31P chemical shifts.     

Visualisation of a 2'-5' parallel stranded double helix at atomic resolution: crystal structure of cytidylyl-2',5'-adenosine

X-ray crystallographic studies on 3'-5' oligomers have provided a great deal of information on the stereochemistry and conformational flexibility of nucleic acids and polynucleotides. In contrast, there is very little information available on 2'-5' polynucleotides. We have now obtained the crystal structure of Cytidylyl-2',5'-Adenosine (C2'p5'A) at atomic resolution to establish the conformational differences between these two classes of polymers. The dinucleoside phosphate crystallises in the monoclinic space group C2, with a = 33.912(4)A, b = 16.824(4)A, c = 12.898(2)A and beta = 112.35(1) with two molecules in the asymmetric unit. Spectacularly, the two independent C2'p5'A molecules in the asymmetric unit form right handed miniature parallel stranded double helices with their respective crystallographic two fold (b axis) symmetry mates. Remarkably, the two mini duplexes are almost indistinguishable. The cytosines and adenines form self-pairs with three and two hydrogen bonds respectively. The conformation of the C and A residues about the glycosyl bond is anti same as in the 3'-5' analog but contrasts the anti and syn geometry of C and A residues in A2'p5'C. The furanose ring conformation is C3' endo, C2' endo mixed puckering as in the C3'p5'A-proflavine complex. A comparison of the backbone torsion angles with other 2'-5' dinucleoside structures reveals that the major deviations occur in the torsion angles about the C3'-C2' and C4'-C3' bonds. A right-handed 2'-5' parallel stranded double helix having eight base pairs per turn and 45 degrees turn angle between them has been constructed using this dinucleoside phosphate as repeat unit. A discussion on 2'-5' parallel stranded double helix and its relevance to biological systems is presented.     

Structural studies on LNA quadruplexes

LNAs (locked nucleic acids) are new DNA analogues with higher binding affinities toward nucleic acids than the canonical counterparts mainly due to the characteristic conformational restriction arising from the 2'-O, 4'-C methylene bridge. In light of the promising therapeutic applications and considering the advantageous characteristics of LNAs, such as their high water solubility, easy handling, and synthetic accessibility through the conventional phosphoramidite chemistry, we undertook a study concerning the capability of these nucleic acid analogues to form quadruplex structures. Particularly, we have been investigating the LNA/DNA chimeras corresponding to the well-known DNA sequences 5-GGTTGGTGTGGTTGG-3', capable of forming an unimolecular quadruplex. This article deals with the study of the sequence 5'-ggTTggTGTggTTgg-3' (upper and lower case letters represent DNA and LNA residues, respectively), which, according to CD spectroscopy, is able to fold into a quadruplex structure.     

Stacking-unstacking of the dinucleoside monophosphate guanylyl-3'',5''-uridine studied with molecular dynamics.

Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.     

Polypurine/polypyrimidine hairpins form a triple helix structure at low pH.

1D and 2D NMR investigations of the 15 residue deoxynucleotide sequence d(TCTCTC-TTT-GAGAGA) show that above pH = 6.5 the molecule adopts a B-form hairpin conformation. As the pH is lowered below 6.5 molecules progressively associate in pairs to form a partially triple helical, partially single stranded structure in which the bases of the oligopyrimidine d(TC)3 tract from one molecule form Hoogsteen pairs with the d(G-A)3 tract of the other. Imino protons of protonated cytosines can be observed at very low field (approximately 15 ppm). The enthalpy of triplex formation was estimated by NMR techniques to be -16 kcal mol-1. Intense H6 to H3' cross peaks from residues in all three strands suggest the presence of N-type sugars at some but not at all possible sites. Surprisingly strong cross peaks between H5' or H5" and non-exchangeable base protons are also observed. These suggest that certain of the O5'-C5'-C4'-C3' phosphate backbone torsion angles (gamma) are unusual.     

Molecular conformation of thyroid hormones: structure and binding interactions of thyroxine

In an accurately determined X-ray diffraction study of the thyroid hormone thyroxine (T4), the two independent conformations in the crystal lattice show significant differences in the outer phenyl ring geometry when compared with that of 3,5,3'-triiodothyronine (T3). The major differences between the T4 and T3 structures are a shortened C4'-O4' bond, contraction of the C3'-C4'-C5' angle and an increase in the C3' and C5' angles of T4. These changes can be correlated with the difference in acidity of the 4'-OH of T4 and T3 and help to explain binding affinity differences among thyroactive compounds. The hydrogen bond directionality observed in T4 and other thyroid structures offers an insight into the molecular details of the hormone-receptor site. The conformation of one T4 molecule is cisoid, that of the other transoid, the first such instance of different overall conformations to be found in the same crystal lattice. One T4 molecule has the side chain nearly coplanar with the inner ring, an unusual conformation among thyroid structures.     

The effect of (2''-5'') and (3''-5'') phosphodiester linkages on conformational and stacking properties of cytidylyl-cytidine in aqueous solution.

Conformational properties of (2'-5') and (3'-5') CpC have been determined by proton magnetic resonance spectroscopy at 220 MHz. The ribose ring structures are predominantly 3E with the exception of the ring from the 2'-phosphate fragment of C(2'-5')pC which exhibits an 2E pucker. Bases are oriented anti with respect to the ribose and the conformations about C4'-C5', C5'-O5', C3'-O3' (C2'-O2') are gg, g'g', and g+ in equilibrium g-, respectively. The dimers exist as mixtures of stacked (g+g+ and g-g- about the P-O(C) bonds) and unstacked species at 20 degrees C. Stacking is estimated to be 35% in both dimers.     

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.     

Solid-state NMR determination of sugar ring pucker in (13)C-labeled 2'-deoxynucleosides

The H3'-C3'-C4'-H4' torsional angles of two microcrystalline 2'-deoxynucleosides, thymidine and 2'-deoxycytidine.HCl, doubly (13)C-labeled at the C3' and C4' positions of the sugar ring, have been measured by solid-state magic-angle-spinning nuclear magnetic resonance (NMR). A double-quantum heteronuclear local field experiment with frequency-switched Lee-Goldberg homonuclear decoupling was used. The H3'-C3'-C4'-H4' torsional angles were obtained by comparing the experimental curves with numerical simulations, including the two (13)C nuclei, the directly bonded (1)H nuclei, and five remote protons. The H3'-C3'-C4'-H4' angles were converted into sugar pucker angles and compared with crystallographic data. The delta torsional angles determined by solid-state NMR and x-ray crystallography agree within experimental error. Evidence is also obtained that the proton positions may be unreliable in the x-ray structures. This work confirms that double-quantum solid-state NMR is a feasible tool for studying sugar pucker conformations in macromolecular complexes that are unsuitable for solution NMR or crystallography.