Ab-inito Quantum Mechanical Calculations of NMR Chemical Shifts in Nucleic Acids Constituents II Conformational Dependence of the lH and 13C Chemical Shifts in the Ribose

Abstract

The magnetic shielding constant of the different ¹³C and ¹³H nuclei of a deoxyribose are calculated for the C2′ endo and C3′ endo puckerings of the furanose ring as a function of the conformation about the C4′C5′ bond. For the carbons the calculated variations are of several ppm, the C3′ endo puckering corresponding in most cases to a larger shielding than the C2′ endo one. For the protons the calculated variations of chemical shifts are all smaller than 1.3 ppm, that is of the order of magnitude of the variation of the geometrical shielding produced on these protons by the other units of a DNA double helix, with a change of the overall structure of the helix. The computations carried out on the deoxyribose ?3′ and 5′ phosphates for several conformations of the phosphate group tend to show that the changes of conformation of the charged group of atoms produce chemical shift variations smaller than the two conformational parameters of the deoxyribose itself. The calculations carried out for a ribose do give the general features of the differences between the carbon and proton spectra of deoxynucleosides and nucleosides.

The comparison of the measured and calculated phosphorylation shifts tend to show that the counterion contributes significantly, for some nuclei of the deoxyribose, to the shifts measured. The calculated magnitude of this polarization effect on carbon shifts suggests a tentative qualitative interpretation of carbon spectra of the ribose part of DNA double helices.     

Secondary-structure dependent chemical shifts in proteins

Chemical shift data have been collected on eight proteins that have the same conformation in solution as in their crystal structures. Ring-current shifts have been calculated and subtracted from the exerimentally measured shifts, to leave shifts that depend only on local conformation. Overall, the shifts show an approximately normal distribution with no appreciable skewness, thus confirming that ring-current shifts have the overall effect of skewing the distribution to high field. In helices, NH and C(alpha)H have a high significant tendency to resonate to high field, whereas they resonate to low field in beta-sheets. Side-chain protons resonate slightly to high field in beta-sheets. Chemical shift distributions are narrowest for side-chain protons, and widest for amide protons. When only slowly exchanging amide protons are considered, the high field shift for amide protons in helices is more pronounced, but there is only a small difference in sheets. C(alpha)H signals at the N-terminal end of helices tend to resonate to higher field than those at the C-terminal end, whereas for NH signals it is the C-terminal end that resonates to higher field. There is no significant effect of position within the helix on side-chain signals, implying that the helix dipole has little effect on shifts within the helix.     

Nuclear magnetic spectra of self complementary decanucleotides in solution; base sequence effect on the chemical shifts of nonexchangeable protons

The aim of this study was to attempt to determine the extent to which the chemical shifts of the nonexchangeable base protons of a DNA helix depend upon the base sequence. We measured the proton NMR spectra of twelve decadeoxynucleotides in order to carry out a "statistical" treatment. In the helices, the chemical shifts were found to be determined within +/- 0.04 ppm, largely by the nearest neighbor residues on the 5'-side, and to a smaller extent by the residue on the 3'-side. The theoretical chemical shift calculations reproduced very well the polymerization shifts measured for H2 protons of adenosines if the electrostatic field effect was taken into account. A fair agreement was also obtained for H8 protons of the adenosine and guanosine residues. However, theory underestimates the polarization effects of the base protons of cytidine. This discrepancy suggests that the conformation of this residue is different in the mononucleotides relative to double helices.     

Parameters for the calculation of proton NMR chemical shifts of oligoribonucleotides

A set of empirical parameters which allows the prediction of the proton NMR chemical shifts at 70 C of non-exchangeable heterobase and anomeric protons in oligoribonucleotides has been constructed. The set is based on the highly flexible nature of oligoribonucleotide single strands and the wide range of conformational states which can be populated at relatively high temperatures (70 C or greater). A pairwise subtractive procedure, using 129 ribonucleotide oligomers (all 16 dimers, all 64 trimers, 37 tetramers, and 12 pentamers), shows that significant contributions to the observed chemical shift of protons in a given nucleoside residue are made by first, second, and third neighbors on the 3' and the 5' sides. The majority of the neighbors cause shielding effects with the exception of some first neighbors on the 5' side of a given residue. The magnitude of the shielding effects is greatest for the purine heterobases and follows the order A greater than G greater than C greater than U, with first neighbors on the 3'side showing more pronounced effects than second neighbors and these in turn showing larger effects than third neighbors. Second neighbors on the 5' side showed consistently greater shieldings than first neighbors, a result attributed to the deshielding effects of the first 5' neighbor phosphate group. The parameter Tables are applied to the prediction of proton chemical shifts in one heptamer, four hexamers, and two pentamers and give average absolute differences between predicted and observed shifts less than 0.030 ppm. The parameter approach represents an excellent method of generating initial assignments of proton chemical shifts for any single strand oligoribonucleotide.     

The conformation of the d(ACCCGGGT) duplex in aqueous solution

The nonexchangeable base and sugar protons of the octanucleotide d(ACCCGGGT)2 have been assigned using two dimensional homonuclear Hartmann-Hahn relayed spectroscopy (HOHAHA), double quantum filtered homonuclear correlation spectroscopy (DQFCOSY) and nuclear Overhauser spectroscopy (NOESY) in D2O at 12 degrees C. The observed NOE's between the base protons and their own H2' protons and between the base protons and the H2' protons of the 5' adjacent nucleotide and the observed coupling constants between the deoxyribose 1' and 2',2' protons indicate that this duplex assumes a right-handed B-type helix conformation in solution.     

Determination of the DNA sugar pucker using 13C NMR spectroscopy

Solid-state 13C NMR spectroscopy of a series of crystalline nucleosides and nucleotides allows direct measurement of the effect of the deoxyribose ring conformation on the carbon chemical shift. It is found that 3'-endo conformers have 3' and 5' chemical shifts significantly (5-10 ppm) upfield of comparable 3'-exo and 2'-endo conformers. The latter two conformers may be distinguished by smaller but still significant differences in the carbon chemical shifts at the C-2' and C-4' positions. High-resolution solid-state NMR of three modifications of fibrous calf thymus DNA shows that these trends are maintained in high-molecular-weight DNA and confirms that the major ring pucker in A-DNA is 3'-endo, while both B-DNA and C-DNA are largely 2'-endo. The data show that 13C NMR spectroscopy is a straightforward and useful probe of DNA ring pucker in both solution and the solid state.     

Ab-initio quantum mechanical calculations of NMR chemical shifts in nucleic acids constituents. III. Chemical shift variations due to base stacking

Ab inito computations of the different contributions to chemical shift variations due to intra and interstrand stacking are reported for the GC, CG, AT and TA sequences of a B DNA helix. The results obtained for the non hydrogen atoms of the GC stacks show that the chemical shift variations are mainly due to the polarization contribution, the term which decreases slowly with the intermolecular distance. Because of the weaker polarity of adenine and thymine the geometric and polarization contributions are of closer absolute magnitude for the non hydrogen atoms of the intrastrand stacks but the polarization term is the determining contribution in the corresponding interstrand stacks. For the protons which undergo smaller shifts due to the polarization (or electric field effects) the role of the geometric contribution is more important and is even the leading one for the hydrogens of cytosine and thymine in the case of intrastrand stacking. The charge transfer plus exchange term has a non negligeable value for a limited number of cases corresponding to the shortest intermolecular interatomic distances. These results are discussed in relation with the qualitative differences observed between the proton and carbon spectra of dinucleotides and B-DNA duplexes.     

Two dimensional NMR studies on the solution structure of d-CTCGAGCTCGAG

Two dimensional (2D) FT-NMR investigations have been carried out on the self-complementary dodecanucleotide d-CTCGAGCTCGAG, which has cleavage sites for the restriction enzyme Xho I (between C and T). The central TCG portion is also known to show a preference for DNAase activity. Complete resonance assignments have been obtained for the non-exchangeable sugar and base protons of the oligonucleotide. Information regarding sugar geometries, glycosidic torsion angles and other structural parameters has been obtained from the relative intensities of the cross peaks in the COSY and NOESY spectra. The results indicate that deoxyribose rings of C1 and C7 adopt a conformation different from the remaining sugars in the double helical oligonucleotide. The central TCG portion also exhibits variations in the backbone structure. The base stacking in the double helix shows interesting sequence dependent effects suggesting that the sequence effects are not localised to nearest neighbours but extended over longer stretches.     

The 1H-NMR assignments of the aromatic resonances in complexes of Lactobacillus casei dihydrofolate reductase and the origins of their chemical shifts

All the aromatic proton resonances in the 500-MHz NMR spectra of Lactobacillus casei dihydrofolate reductase have been assigned for several of its complexes with inhibitors. For the complexes with methotrexate and trimethoprim this was achieved by using a combination of NMR techniques in conjunction with a selectively deuterated protein designed to simplify the spectra such that nuclear Overhauser effect (NOE) connections could be detected with greater ease and certainty. By correlating these NOE data with crystal structure data on related complexes it was possible to assign all the aromatic resonances and to extend these assignments to spectra of other complexes of dihydrofolate reductase. The conformation-dependent chemical shifts observed for many of the resonances could be explained qualitatively, but not quantitatively, in terms of ring-current shifts. The discrepancies between calculated ring-current shifts and the observed conformation-dependent shifts could not in general be accounted for satisfactorily in terms of carbonyl-group anisotropic shielding contributions calculated using presently available models. In the case of the H delta 1, delta 2 protons of Phe30 some of the discrepancy probably results from a difference in the conformation of the Phe ring between the solution and crystal states.     

Flexibility of the furanose ring in nucleic acids: a compilation of crystal data

A compilation of crystal structure data on deoxyribo- and ribonucleosides and their higher derivatives is presented. The aim of this paper is to highlight the flexibility of deoxyribose and ribose rings. So far, the conformational parameters of nucleic acids constituents of ribose and deoxyribose have not been analysed separately. This paper aims to correlate the conformational parameters with the nature and puckering of the sugar. Deoxyribose puckering occurs in the C2′ endo region while ribose puckering is observed both in the C3′ endo and C2′ endo regions. A few endocyclic and exocyclic bond angles depend on the puckering and the nature of the sugar. The majority of structures have an anti conformation about the glycosyl bond. There appears to be a puckering dependence on the torsion angle about the C4′C5′ bonds. Such stereochemical information is useful in model building studies of polynucleotides and nucleic acids.     

Phosphorus-31 nuclear magnetic resonance of double- and triple-helical nucleic acids. Phosphorus-31 chemical shifts as a probe of phosphorus-oxygen ester bond torsional angles

The temperature dependence to the 31P NMR spectra of poly[d(GC)] . poly [d(GC)],d(GC)4, phenylalanine tRNA (yeast) and mixtures of poly(A) + oligo(U) is presented. The 31P NMR spectra of mixtures of complementary RNA and of the poly d(GC) self-complementary DNA provide torsional information on the phosphate ester conformation in the double, triple, and "Z" helix. The increasing downfield shift with temperature of the single-strand nucleic acids provides a measure of the change in the phosphate ester conformation in the single helix to coil conversion. A separate upfield peak (20-60% of the total phosphates) is observed at lower temperatures in the oligo(U) . poly(A) mixtures which is assigned to the double helix/triple helix. Proton NMR and UV spectra confirm the presence of the multistrand forms. The 31P chemical shift for the double helix/triple helix is 0.2-0.5 ppm upfield from the chemical shift for the single helix which in turn is 1.0 ppm upfield from the chemical shift for the random coil conformation.     

Conformational Analysis of a Hybrid DNA-RNA Double Helical Oligonucleotide in Aqueous Solution: d(CG)r(CG)d(CG) Studied by 1D- and 2D-1H NMR Spectroscopy

Abstract

The double helical structure of the self-complementary DNA-RNA-DNA hybrid d(CG)r(CG) d(CG) was studied in solution by 500 MHz ¹H-NMR spectroscopy. The non-exchangeable base protons and the (deoxy)ribose H1′, H2′ and H2″ protons were unambiguously assigned using 2D-J-correlated (COSY) and 2D-NOE (NOESY) spectroscopy techniques. A general strategy for the sequential assignment of ¹H-NMR spectra of (double) helical DNA and RNA fragments by means of 2D-NMR methods is presented.

Conformational analysis of the sugar rings of d(CG)r(CG)d(CG) at 300 K shows that the central ribonucleotide part of the helix adopts an A-type double helical conformation. The 5′- and 3′-terminal deoxyribose base pairs, however, take up the normal DNA-type conformation. The A-to-B transition in this molecule involves only one (deoxyribose) base pair. It is shown that this A-to-B conformational transition can only be accomodated by two specific sugar pucker combinations for the junction base pair, i.e. N·S (C3′-endo-C2′-endo, 60%, where the pucker given first is that assigned to the junction nucleotide residue of the strand running 5′ → 3′ from A-RNA to B-DNA) and S·S (C2′-endo-C2′-endo, 40%).     

H8 chemical shifts in oligonucleotide cross-linked at a GpG sequence by cis-Pt(NH3)2(2+): a clue to the adduct structure.

The origin of the anomalous H8 chemical shifts observed in 1H-NMR spectra of oligonucleotides cross-linked at a GpG sequence with cis-[Pt(NH3)2]2+ has been investigated and clarified. The main contributions that distinguish the H8 resonances of the two platinum-ligating guanines from other GH8 signals and from each other are: (a) the inductive effect of platinum binding which we have recently quantified as a downfield shift of 0.48 +/- 0.07 ppm (M. H. Fouchet, D. Lemaire, J. Kozelka and J.-C. Chottard, unpublished results); (b) the ring-current effect of one GpG guanine on the H8 resonance of the other guanine, which is negative (shielding) for the 5'-H8 and positive (deshielding) for the 3'-H8 in single-stranded adducts, but has the opposite sign in double-stranded adducts; (c) a deshielding polarization effect of the phosphate 5' to the GpG unit. The different signs of the ring-current effects in single-stranded and double-stranded oligonucleotides originate from the orientation of the guanines in the cis-[Pt(NH3)2(Gua)2]2+ moiety (Gua, guanine), which is left-handed helicoidal in single strands and right-handed helicoidal in double strands. In the platinated dinucleotides (cis-[Pt(NH3)2(GpG)]+, cis-[Pt(NH3)2(d(GpG))]+ and cis-[Pt(NH3)2(d(GpG))]), the guanines assume either the left-handed or the right-handed arrangement, depending on the sugar moiety (ribose or deoxyribose), protonation state at N1 and, in the solid state, on crystal forces. This work shows that chemical shifts contain valuable structural information which is complementary to that extracted from correlated spectroscopy and nuclear Overhauser spectroscopy data.     

Prediction of 1H NMR chemical shifts of DNA oligomers

A set of parameters, devised for the prediction of 1H NMR chemical shifts of heterobase and anomeric protons in the high temperature (greater than 70 degrees C) spectra of RNA oligomers has been found to be applicable to the corresponding DNA oligomers. Fifteen examples of DNA oligomers that have had high temperature spectra recorded and assigned show a mean absolute difference between predicted and assigned shifts of 0.045 ppm. The parameters for uridine H-5 are applied to the calculation of thymidine methyl group shifts and give excellent agreement with experimental assigned shifts. The RNA parameter set is a practical means of assigning heterobase and anomeric protons in DNA oligomers. A programme using the RNA parameter set has been written which enables the sequence of short DNA oligomers to be predicted from their 1H NMR spectra.     

The effect of the trp repressor from Escherichia coli on the structure and dynamics of dA20dT20

We have determined the effect of the tryptophan (trp) repressor from Escherichia coli on the structure and dynamics of dA20dT20. The structure was determined using time-dependent nuclear Overhauser effects and spin-lattice relaxation times. The deoxyribose conformation is near C3' endo for the thymine residues, and a mixture of about 30% C3' endo and 70% C2' endo for the adenine residues. The glycosidic torsion angles are -50 degrees for T and -60 degrees for A. The roll is 20 degrees and the propellor twist is about 29 degrees. The conformation is consistent with recent calculations (Rao, K. and Kollman, P.A. (1985) J. Am. Chem. Soc. 107, 1507-1511). The rate constant for exchange of the imino protons is similar to that usually found for AT base-pairs, with an activation energy of 20 +/- 2 kcal/mol, and an activation entropy of 17 +/- 7 cal/mol per K. The repressor greatly retards the exchange of imino protons, and the activation energy increases to 38 kcal/mol. There are small changes in the structure of the DNA on forming the complex, with the adenine and thymidine residues becoming more similar in conformation.     

Ab initio quantum mechanical calculations of the variation of the 1H and 13C nuclear magnetic shielding constants in proline as a function of the angle psi

The variation of the nuclear magnetic shielding constant of the different protons and carbons of trans HCO-L-Pro-NH2 with the value of the angle psi is calculated by a non-empirical method for three conformations of the proline ring. The results concerning the CH protons show that the chemical shift of the alpha, beta and gamma endo hydrogens can vary by more than 1 ppm when psi goes from -30 degrees to 180 degrees. The theoretical variation of the chemical shift difference between alpha and gamma or beta and gamma carbons is found to be sensitive to the puckering of the proline ring. For the second of these differences the theoretical results are in agreement with Siemion's relation only for a limited range of molecular conformations. Additional calculations show that the variations of the proton shifts with the value of psi are due to the magnetic anisotropy of the proline carbonyl group and to the polarization of the CH bonds by the multipolar charge distribution carried by this carbonyl. The results are discussed in relation to experiment and the possibility of using 1H and 13C chemical shifts for the determination of the value of the torsion angle about the C alpha C' bond.     

Prediction of fluorine chemical shifts in proteins

Molecular dynamics calculations have been used in an effort to estimate the change in fluorine nmr shielding when a fluorine nucleus enters the tertiary structure of a protein. Considerations of the possible interactions that can define the shift parameter change suggest that van der Waals interactions are the leading determinant of fluorine shifts in proteins, although aromatic ring currents, other magnetic anisotropies, and electrostatic field effects could result in shift distinctions of 1 ppm or smaller. Results of our studies of a fluorine-containing analogue of the ribonuclease A S-protein/S-peptide complex indicate that static structures such as those implied by crystallographic data lead to overestimates of the magnitude of the van der Waals shielding term; molecular dynamics simulations provide indications of the effects of conformational averaging in defining this term. The treatment used predicts the correct direction of the shift change when the fluorine enters this protein environment from aqueous solution and, with an experimentally supported choice of adjustable parameters, gives agreement with the magnitude of the shift.     

Conformation and interaction of short nucleic acid double-stranded helices. I. Proton magnetic resonance studies on the nonexchangeable protons of ribosyl ApApGpCpUpU.

1H nuclear magnetic resonance (NMR) spectra of a self-complementary ribosyl hexanucleotide, A2GCU2, are investigated as a function of temperature and ionic strength in D2O. Seventeen nonexchangeable base and ribose-H1' resonances are resolved, and unequivocally assigned by a systematic comparison with the spectra of a series of oligonucleotide fragments of the A2GCU2 sequence varying in chain length from 2 to 5. Changes in the chemical shifts of the 17 protons from the hexamer as well as the six H1'-H2' coupling constants are followed throughout a thermally induced helix-coil transition. These sigma vs. T and J vs. T (degrees C) profiles indicate that the transition is not totally cooperative and that substantial populations of partially bonded structures must exist at intermediate temperatures, with the central G-C region being most stable. Transitions in chemical shift for protons in the same base pair exhibit considerable differences in their Tm values as the data reflect both thermodynamic and local magnetic field effects in the structural transition, which are not readily separable. However, an average of the Tm values agrees well with the value predicted from studies of the thermally induced transition made by optical methods. The values of J1'-2' for all six residues become very small (less than 1.5 Hz) at low temperatures indicating that C3'-endo is the most heavily populated furanose conformation in the helix. The sigma values of protons in the duplex were compared with those calculated from the ring current magnetic anisotropies of nearest and next-nearest neighboring bases using the geometrical parameters of the A'-RNA and B-DNA models. The sigma values of the base protons in the duplex calculated assuming the A'-RNA geometry agree (+/- approximately 0.1 ppm) with the observed values much more accurately than those calculated on the basis of B-DNA geometry. The measured sigma values of the H1' are not accurately predicted from either model. The synthesis of 35 mg of A2GCU2 using primer-dependent polynucleotide phosphorylase is described in detail with extensive discussion in the microfilm edition.     

Solution conformation of an oligonucleotide containing a G.G mismatch determined by nuclear magnetic resonance and molecular mechanics.

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three dimensional solution structure of the non-selfcomplementary oligonucleotide, d(GAGGAGGCACG). d(CGTGCGTCCTC) in which the central base pair is G.G. This is the first structural determination of a G.G mismatch in a oligonucleotide. Two dimensional nuclear magnetic resonance spectra show that the bases of the mismatched pair are stacked into the helix and that the helix adopts a classical B-DNA form. Spectra of the exchangeable protons show that the two guanosines are base paired via their imino protons. For the non-exchangeable protons and for some of the exchangeable protons nuclear Overhauser enhancement build up curves at short mixing times have been measured. These give 84 proton-proton distances which are sensitive to the helix conformation. One of the guanosines adopts a normal anti conformation while the other is syn or close to syn. All non-terminal sugars are C2' endo. These data sets were incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. The G.G mismatch shows a symmetrical base pairing structure. Although the mismatch is very bulky many of its features are close to that of normal B-DNA. The mismatch induces a small lateral shift in the helix axis and the sum of the helical twist above and below the mismatch is close to that of B-DNA.     

High-resolution NMR study of a synthetic DNA-RNA hybrid dodecamer containing the consensus pribnow promoter sequence: d(CGTTATAATGCG).r(CGCAUUAUAACG)

The nonsymmetrical double-helical hybrid dodecamer d(CGTTATAATGCG).r(CGCAUUAUAACG) was synthesized with solid-phase phosphoramidite methods and studied by high-resolution 2D NMR. The imino protons were assigned by one-dimensional nuclear Overhauser methods. All the base protons and H1', H2', H2", H3', and H4' sugar protons of the DNA strand and the base protons, H1', H2', and most of the H3'-H4' protons of the RNA strand were assigned by 2D NMR techniques. The well-resolved spectra allowed a qualitative analysis of relative proton-proton distances in both strands of the dodecamer. The chemical shifts of the hybrid duplex were compared to those of the pure DNA double helix with the same sequence (Wemmer et al., 1984). The intrastrand and cross-strand NOEs from adenine H2 to H1' resonances of neighboring base pairs exhibited characteristic patterns that were very useful for checking the spectral assignments, and their highly nonsymmetric nature reveals that the conformations of the two strands are quite different. Detailed analysis of the NOESY and COSY spectra, as well as the chemical shift data, indicate that the RNA strand assumes a normal A-type conformation (C3'-endo) whereas the DNA strand is in the general S domain but not exactly in the normal C2'-endo conformation. The overall structure of this RNA-DNA duplex is different from that reported for hybrid duplexes in solution by other groups (Reid et al., 1983a; Gupta et al., 1985) and is closer to the C3'-endo-C2'-endo hybrid found in poly(dA).poly(dT) and poly(rU).poly(dA) in the fiber state (Arnott et al., 1983, 1986).