A Method for the Determination of Furanose Ring Coordinates in its Pseudorotation Circuit for Different Amplitudes of Pucker

Abstract

Interconversion between energetically favored molecular conformations must proceed through less favored intermediate states. Thus, a knowledge of the nucleotide furanose ring conformations, other than the crystallographically well-determined ones, are of interest in investigating nucleotide conformational energies and dynamics. The sugar ring flexibility affects the conformation and dynamics of the monomer and determines the range of feasible nucleic acid secondary and tertiary structures. We have generated furanose geometries for varying amplitudes of pucker over its entire range of pseudorotation by making use of a ring closure procedure and the empirical dependence of endocyclic bond lengths and bond angles on sugar pucker. Atomic coordinates are tabulated here for the furanose ring at pseudorotation phase angle intervals of 9° for the average amplitude (τ_m) of pucker of 39° as well as for decreased (20° and 30°) and increased (44°) values of τ_m. However, the coordinates for any values of P and τ_m can be readily calculated.     

Stereochemical studies on nucleic acid analogues. I. Conformations of alpha-nucleosides and alpha-nucleotides: interconversion of sugar puckers via O4'-exo.

The preferred conformations of ribo and deoxyribo alpha-nucleosides and alpha-nucleotides, the stereoisomers of the naturally occurring beta-isomers, are worked out by minimizing the conformational energy as a function of all the major parameters including the sugar ring conformations along the pseudorotation path. The results of the studies bring out certain distinct conformational features that are at variance with their beta counterparts. The range of glycosyl conformations are found to be not only quite restricted here but favor predominantly the anti conformation. The syn glycosyl conformation for the entire region of P values are found to be energetically less favorable, with the barrier to anti in equilibrium with syn interconversion being higher especially in alpha-ribonucleosides. The energetically preferred range of pseudorotation phase angles (P) is also considerably restricted and P values corresponding to the C1'-exo range of sugars are highly unfavorable for alpha-nucleosides, in sharp contrast to the broad range of sugar ring conformations favored by beta-isomers. While both trans congruent to 180 degrees and skew congruent to 270 degrees conformations around the C3'-O3' (phi') bond are favored for alpha-3'-nucleotides with deoxyribose sugars, ribose sugars seem to favor only the skew values of phi'. Most interestingly and in sharp contrast to beta-stereoisomers, an energy barrier is encountered at P values corresponding to O4'-endo sugars. This suggests that the possible sugar pucker interconversion between C2'-endo/C3'-exo and C3'-endo/C2'-exo in alpha-anomers could take place only through the O4'-exo region. Likewise the possible path of anti in equilibrium with syn interconversion in alpha-nucleosides is not via high anti, in sharp contrast to beta-nucleosides. These observations should be borne in mind while assigning the sugar ring conformers in alpha-nucleosides and those containing them from nmr investigations. Comparison of the results with beta-anomers seem to suggest on the whole a lack of conformational variability or the restricted nature of alpha-stereoisomers. This could be one of the reasons for its nonselection in the naturally occurring nucleic acids.     

Tautomerism and conformation of the promutagenic analogue N6-methoxy-2'',3'',5''-tri-O-methyladenosine

N6-Methoxy-2',3',5'-tri-O-methyladenosine crystallizes in space group P2(1)2(1)2(1) with cell dimensions a = 4.693, b = 11.412, c = 31.741 A. Least-squares refinement of diffractometer data converged at R = 0.038. The location of a hydrogen atom at N1 and the observed bond lengths and bond angles indicate unequivocally the imino tautomer of the adenine moiety. The N6-methoxy group is oriented syn to N1 and the glycosidic torsion angle XCN is -3.6 degrees, i.e. in the anti range. The furanose ring has a C2'-exo/C3'- endo pucker (P = 0.9 degrees) and is unusually flattened (tau m = 30.0 degrees). The conformations of the O-methyl groups of the ribose ring are compared with those of monomethylated nucleosides, including the biologically important 2'-O-methyl nucleosides. Evidence is presented for the existence of C-H ... N intermolecular hydrogen bonds between adenine moieties. Bearing in mind that N6-methoxyadenosine is a promutagenic analogue, the results are compared with those for the corresponding promutagenic N4-methoxycytidine. They are also discussed in relation to the tautomerism, the conformation of the N6-methoxy group, and the associated base-pairing abilities in the absence and presence of polymerases.     

Crystal structure and conformation of 5-fluorouridine: conformational preferences for 5-fluorinated pyranosides

Crystals of 5-fluorouridine (5FUrd) have unit cell dimensions a = 7.716(1), b = 5.861(2), c = 13.041(1)A, alpha = gamma = 90 degrees, beta = 96.70 degrees (1), space group P2(1), Z = 2, rho obs = 1.56 gm/c.c and rho calc = 1574 gm/c.c The crystal structure was determined with diffractometric data and refined to a final reliability index of 0.042 for the observed 2205 reflections (I > or = 3sigma). The nucleoside has the anti conformation [chi = 53.1(4) degrees] with the furanose ring in the favorite C2'-endo conformation. The conformation across the sugar exocyclic bond is g+, with values of 49.1(4) and -69.3(4) degrees for phi(theta c) and phi (infinity) respectively. The pseudorotational amplitude tau(m) is 34.5 (2) with a phase angle of 171.6(4) degrees. The crystal structure is stabilized by a network of N-H...O and O-H...O involving the N3 of the uracil base and the sugar 03' and 02' as donors and the 02 and 04 of the uracil base and 03' oxygen as acceptors respectively. Fluorine is neither involved in the hydrogen bonding nor in the stacking interactions. Our studies of several 5-fluorinated nucleosides show the following preferred conformational features: 1) the most favored anti conformation for the nucleoside [chi varies from -20 to + 60 degrees] 2) an inverse correlation between the glycosyl bond distance and the chi angle 3) a wide variation of conformations of the sugar ranging froni C2'-endo through C3'-endo to C4'-exo 4) the preferred g+ across the exocyclic C4'-C5' bond and 5) no role for the fluorine atom in the hydrogen bonding or base stacking interactions.     

The conformationally constrained N-methanocarba-dT analogue adopts an unexpected C4'-exo sugar pucker in the structure of a DNA hairpin

Incorporation of a bicyclo[3.1.0]hexane scaffold into the nucleoside sugar was devised to lock the embedded cyclopentane ring in conformations that mimic the furanose North and South sugar puckers. To analyze the effects of North-methanocarba-2'-deoxythymidine (N-MCdT) on the B-form DNA, we crystallized d(CGCGAA[mcTmcT]CGCG) with two N-MCdTs. Instead of a duplex, the 12mer forms a tetraloop hairpin, whereby loop N-MCdTs adopt the C4'-exo pucker (NE; P = 50°). Thus, the bicyclic framework does not limit the pucker to the anticipated C2'-exo range (NNW; P = -18°).     

The crystal and molecular structure of 2'-deoxy-2'-fluoroinosine monohydrate.

The structure of the hydrate of 2'-deoxy-2'-fluoroinosine has been determined by single-crystal x-ray diffraction. The nucleoside crystallizes in space group P2(1)2(1)2(1) with unit cell dimensions, a = 33.291, b = 10. 871, c = 6.897A. There are two nucleosides and two water molecules in the asymmetric unit. The structure was solved by direct methods and refined to a residual R = 0.095. The two independent nucleosides in the asymmetric unit show different conformations about the glycosidic bond, while other structural details are similar. The base orientation to the sugar is syn in molecule A, whereas anti in molecule B. The exocyclic C(4')-C(5') bond conformation defined with respect to C(3')-C(4')-C(5')-O(5') is gauche+ in both molecules A and B. The sugar ring pucker defined by the pseudorotation phase angle P is a twisted conformation in both, C(3')-endo-C(4')-exo with P = 29 degrees in molecule A and C(4')-exo-C(3')-endo with P = 41 degrees in molecule B. It is shown by comparison with x-ray results of other 2'-fluoronucleosides and unmodified nucleosides including inosines that, in addition to a strong preference of the C(3')-endo type pucker, twisted conformations involving C(4')-exo puckering may be one of characteristic features of 2'-fluoronucleosides.     

Analysis of intrasugar interproton NOESY cross-peaks as an aid to determine sugar geometries in DNA fragments

A systematic analysis of the conformation of deoxyribofuranose rings in DNA fragments has been described using two-dimensional nuclear Overhauser effect spectroscopy (2D NOESY). The approach is based on the interpretation of the intrasugar proton-proton distances which can be estimated using a low-mixing-time pure-absorption mode w1-scaled NOESY spectrum. The experimental distances are compared with the theoretical values calculated as a function of pseudorotation phase angle (P) describing the sugar geometries. The approach can be used as a complementary aid to J couplings for establishing sugar conformations in individual nucleotide units of DNA fragments. Using this strategy on d-ACATCGATGT, we observed that individual nucleotides exhibit O4'-endo sugar pucker. The results rule out possibilities of the existence of a fast equilibrium (on the NMR time scale) between C2'-endo (or S-domain) and C3'-endo (or N-domain) sugar puckers.     

Crystal structure of formycin 5'-phosphate: an explanation for its tight binding to AMP nucleosidase

Formycin 5'-monophosphate (FMP) is a strong competitive inhibitor of AMP nucleosidase with Km/Kis from 1200 to 2600 depending on the source of the enzyme. The crystal structure of FMP has been determined in order to understand the basis for its high affinity for AMP nucleosidases and other biological properties. The key structural features of FMP are (1) the base is the N(7)-H tautomer, (2) the N(3) of the base forms an intramolecular hydrogen bond to the phosphate oxygen O(1), (3) the glycosyl torsion angle is syn with O(4')-C(1') relative to C(9)-C(4) being -6.43 degrees, and (4) the furanose ring pucker is C(3')-endo, with a pseudorotation angle of 20.3 degrees. The major difference between the AMP and FMP structures is that the glycosyl torsion angles differ by 190 degrees. The computed conformational energy necessary to distort AMP so that it has the same glycosyl torsion angle as FMP is 4.6 kcal/mol. This corresponds to a 2100-fold difference in binding energy, in good agreement with the observed interaction between AMP nucleosidase and FMP.     

Conformational Studies of Nucleic Acids I. A Rapid and Direct Method for Generating Furanose Coordinates from the Pseudorotation Angle

Abstract

The greatest difficulty in modeling a nucleic acid is generating the coordinates of its furanoses. This difficulty arises from constraints imposed by the closed ring geometries of these sugars. We have developed a new method for modeling these furanose rings. Using this method, the coordinates of a sugar can be obtained quickly and unambiguously for any point on the pseudorotational pathway from one parameter: the phase angle of pseudorotation P. The significant difference between this and previous sugar modeling schemes is that here the endocyclic bond lengths of the five-membered sugar ring are allowed to vary a small amount according to simple, explicit, and experimentally reasonable analytic functions of P. The coefficients of these functions follow from the empirical behavior of the endocyclic bond angles and from geometrical constraints due to ring closure. The ability to model the sugars directly from one parameter greatly facilitates carrying out the global conformational studies on nucleic acid constituents which will be attempted in subsequent papers of this series.     

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.     

Structure of triphosphoryl nucleotide bound at the active site of yeast hexokinase: 1H-nuclear magnetic resonance study

Conformation of a nonhydrolyzable adenosine triphosphate (ATP) analogue, adenylyl-(,-methylene)-diphosphonate (AMPPCP) bound at the active site of yeast hexokinase-PII was determined by proton two-dimensional transferred nuclear Overhauser effect spectroscopy (TRNOESY) and molecular dynamics simulations. The effect of the glucose-induced domain closure on the conformation of the nucleotide was evaluated by making measurements on two different complexes: PIIAMPPCPMg(II) and PIIGlcAMPPCPMg(II). TRNOE measurements were made at 500 MHz, 10°C, as a function of several mixing times varying in the range of 40 to 200 ms. Interproton distances derived from the analysis of NOE buildup curves were used as restraints in molecular dynamics simulations to determine the conformation of the enzyme bound nucleotide. The adenosine moiety was found to bind in high anti conformation with a glycosidic torsion angle = 48 ± 5 degrees in both complexes. However, significant differences in the conformations of the ribose and triphosphoryl chain of the nucleotide are observed between the two complexes. The phase angles of pseudorotation P in PIIAMPPCPMg(II) and PIIGlcAMPPCPMg(II) are 87 degrees and 77 degrees, describing a OE and OT₄ sugar pucker and the amplitudes of the sugar pucker () are 37 degrees and 61 degrees, respectively.     

Stereochemical effects of non-ionic methylphosphonates on nucleotide conformations.

The preferred conformations of deoxyribo and ribonucleoside 3'-methylphosphonates are analysed by minimizing the conformational energy as a function of all the major parameters including the sugar ring for both the S- and R-isomers. The results show that neither the substitution nor the nature of the diastereomer affects significantly the preferred conformations compared to the naturally occurring nucleoside 3'-phosphates. The preferred range of C3'-O3' bond torsions or the phase angles of pseudorotation (P) of the sugar are unaffected. The chiral substitution on the phosphate always adopts a conformation distal to the secondary C3' carbon atom in the minimum energy conformational state. Further, it introduces certain restrictions on the preferred range of P-O3' torsions depending on the methylphosphonate configuration. Methylphosphonate, especially the S-isomer, renders the normal gauche- range of P-O3' bond torsions responsible for the stacked helical duplexes to be energetically unfavourable besides introducing a high energy barrier between trans and gauche conformations. Therefore it is suggested that duplexes with S-methylphosphonate may favour extended phosphodiester conformations. These factors explain the observed lower melting temperature as well as the downfield shifts in the 31P signals in duplexes containing the S-isomer.     

Correlated variations of bond lengths in pseudorotating furanose rings

Correlated variations of bond lengths in pseudorotating furanose rings are investigated by a theoretical method. At first, matrix equations are proposed to determine the spatial coordinates of the ring atoms from the bond lengths, the bond angles, and the pseudorotation parameters. Secondly, a necessary functional form of the variations of the bond lengths of five-membered rings is derived from a consideration of symmetry. Finally, demonstrations are performed on a furanose ring whose bond angle variations have been precisely determined by experimental analyses. The resulting bond length variations are: delta Ri = beta icos(8/5 pi.(i-2)+2P) where delta Ri is the variation of the bond length between atoms i and i+1, P is the pseudorotation phase, and beta i is a negative constant about -0.01 A. These bond length variations are balanced on the apparent strains of the bond lengths and the bond angles.     

Why ribose was selected as the sugar component of nucleic acids

During evolution ribose was selected as the exclusive sugar component of nucleic acids. The selection is explained by using molecular models and by eliminating most of the other common sugars by looking at their chemical structure and envisioning how they would fit in a nucleic acid model. Comparisons of sugar pucker conformations and configurations of pentoses indicate that ribose was not randomly selected but the only choice, since beta-D-ribose fits best into the structure of physiological forms of nucleic acids. In other nucleotides containing arabinose, xylose, or lyxose, the C(2)'-OH and/or the C(3)'-OH are above the furanose ring, causing steric interference with the bulky base and the C(5)'-OH group.     

Conformational studies of nucleic acids. I. A rapid and direct method for generating furanose coordinates from the pseudorotation angle

The greatest difficulty in modeling a nucleic acid is generating the coordinates of its furanoses. This difficulty arises from constraints imposed by the closed ring geometries of these sugars. We have developed a new method for modeling these furanose rings. Using this method, the coordinates of a sugar can be obtained quickly and unambiguously for any point on the pseudorotational pathway from one parameter: the phase angle of pseudorotation P. The significant difference between this and previous sugar modeling schemes is that here the endocyclic bond lengths of the five-membered sugar ring are allowed to vary a small amount according to simple, explicit, and experimentally reasonable analytic functions of P. The coefficients of these functions follow from the empirical behavior of the endocyclic bond angles and from geometrical constraints due to ring closure. The ability to model the sugars directly from one parameter greatly facilitates carrying out the global conformational studies on nucleic acid constituents which will be attempted in subsequent papers of this series.     

Conformational analysis of the 2'-deoxyribofuranose ring from proton-proton coupling constants: analysis of a nucleoside-carcinogen adduct formed from 2-acetylaminofluorene utilizing a three-state model

The conformation of the sugar moiety of 8-(N-fluoren-2-ylamino)-2′-deoxyguanosine in solution has been examined as a function of temperature by ¹H-nmr spectroscopy. Analysis of coupling constants shows that lowering the temperature to ?50°C in methanol shifts the conformational equilibrium of the sugar ring resulting in a C2′-endo conformation at a mole fraction of 0.97. The computed phase angle of pseudorotation and amplitude of pucker are 154° and 36°, respectively, with very little discrepancy between the five calculated coupling constants and coupling constants extrapolated from the temperature profiles. A computer program has been written enabling a three-state best-fit analysis. The three-state analysis indicates an equilibrium between C2′-endo, C3′-endo, and 04′-endo conformations. In aqueous solution, the computed mole fraction of the 04′-endo form is 0.18 at 30°C. The conformation associated with the sugar ring and the C4′? C5′ bond is compared to that of 2′-deoxyguanosine.     

Deoxyribose conformation in [d(GTATATAC)]2: evaluation of sugar pucker by simulation of double-quantum-filtered COSY cross-peaks

Exchangeable and nonexchangeable proton and phosphorus resonances (11.75 T) of [d(GTATATAC)]2 in aqueous solution were assigned by using proton two-dimensional nuclear Overhauser effect (2D NOE) spectra, homonuclear proton double-quantum-filtered COSY (2QF-COSY) spectra, proton spin-lattice relaxation time measurements, and 31P1H heteronuclear shift correlation spectra. Due to the large line widths, it was not possible to directly extract vicinal proton coupling constant values from any spectrum including ECOSY or 2QF-COSY. However, comparison of quantitative 2QF-COSY spectral simulations with experimental spectra enabled elucidation of coupling constants. The scope and limitations of this approach were explored by computation and by use of experimental data. It was found that proton line widths exhibit some variability from one residue to the next as well as from one proton to the next within a residue and the exact line width is critical to accurate evaluation of coupling constants. Experimental 2QF-COSY spectra were not consistent with a rigid deoxyribose conformation for any of the nucleotide residues. A classical two-state model, with rapid jumps between C2'-endo (pseudorotation angle P = 162 degrees) and C3'-endo (P = 9 degrees) conformations, was able to account for the spectral characteristics of terminal residue sugars: 60% C2'-endo and 40% C3'-endo. However, the 2QF-COSY cross-peaks from the -TATATA- core could be simulated only if the classical two-state model was altered such that the dominant conformer had a pseudorotation angle at 144 degrees instead of 162 degrees. In this case, the major conformer amounted to 80-85%. Alternatively, the spectral data were consistent with a three-state model in which C2'-endo and C3'-endo conformations had the largest and smallest populations, respectively, but a third conformer corresponding to C1'-exo (P = 126 degrees) was present, consistent with recent molecular dynamics calculations. This alternative yielded populations of 50% (P = 162 degrees), 35% (P = 126 degrees), and 15% (P = 9 degrees) for the -TATATA- sugars. The spectral results indicate little variation of sugar pucker between T and A. Small differences in cross-peak component intensities and characteristic spectral distortions, however, do suggest some unquantified variation. 31P1H heteronuclear chemical shift correlation spectra manifested alternating chemical shifts and coupling constants suggestive of phosphodiester backbone conformational differences between TA and AT junctions.