Conformational analysis of an RNA triplet in aqueous solution: m6 2 A-U-m6 2 A studied by two-dimensional nuclear magnetic resonance at 500 MHz

Influence of uracil on the conformational behaviour of RNA oligonucleotides in solution

NMR studies were carried out on some alternating pyrimidine-purine sequences: the single-stranded tetramers CACA and UGUG and the self-complementary octamer CACAUGUG. Assignments, based upon COSY, homonuclear Hartmann-Hahn, and NOESY experiments, are given for the resonances of all base protons and of several sugar protons. Chemical shift vs temperature profiles were used to obtain thermodynamic parameters for the single-stranded stack in equilibrium with random coil and the duplex in equilibrium with random coil equilibria. The populations of N-type conformer of the ribose rings were estimated from the observed J1'2'. Comparisons with another alternating pyrimidine-purine sequence Um2(6)AUm2(6)A and with the deoxyribose counterparts d(CACA), d(TGTG) and d(CACATGTG) are given. Previous 1H-NMR investigations of Um2(6)AUm2(6)A revealed that the population of bulge-out structure diminishes compared to m2(6)AUm2(6)A due to the U(1)-m2(6)A(2) stacking interaction. In CACA a strong stacking proclivity (Tm = 310 K) together with a clear preference for N-type ribose is observed. However, the stacking interactions in UGUG are relatively less stable (Tm = 288 K) and a bias towards S-type sugar is present. Besides a small amount of stack, a significant contribution of bulge out structure is proposed for UGUG. We conclude that the nature of the pyrimidine base mainly determines the formation of bulge-out structures. The poor stacking properties of uracil now appear to be mainly responsible for this phenomenon. Comparison with the deoxyribose counterparts shows a reasonable agreement between the Tm values of CACA and d(CACA), whereas the Tm of UGUG (288 K) is much lower than the Tm of d(TGTG) (315 K). It is suggested that the absence of bulge-out structures in DNA purine-pyrimidine-purine sequences is related to the relatively strong stacking proclivity of dT residues compared to that of U residues. The Tm values (average 341 K) for the duplex in equilibrium with random coil transition obtained for each residue of CACAUGUG appear very similar. All ribose rings, except the G(8), adopt a pure N conformer in the duplex. This is taken to mean that the differences in conformational behaviour of the constituent tetramers disappear upon duplex formation.

Conformational analysis of the tetranucleotides m6(2)A-m6(2)A-U-m6(2)A(m6(2)A = N6-dimethyladenosine) and U-m6(2)A-U-m6(2)A and of the hybrid dA-r(U-A). A one- and two-dimensional NMR study

A 1H-NMR investigation was carried out on the tetranucleotides U-m6(2)A-U-m6(2)A and m6(2)A-m6(2)A-U-m6(2)A (m6(2) = N6-dimethyladenosine) as well as on the hybrid trinucleotide dA-r(U-A). An extensive comparison with m6(2)A-U-m6(2)A and other relevant compounds is made. Previous proton NMR studies on trinucleotides have shown that purine-pyrimidine-purine sequences prefer to adopt a mixture of states which have as a common feature that the interior pyrimidine residue bulges out, whereas the flanking purine residues stack upon each other. A stacking interaction on the 3' side of the bulge is known to have no measurable effect on the bulge population. Chemical-shift data, ribose ring conformational analysis and information from NOE experiments now show unambiguously that the moderate U(1)-m6(2)A(2) stack in U-m6(2)A-U-m6(2)A diminishes the population of bulged-out structures in favour of a regular stack. This tendency towards conformational transmission in the downstream 5'----3' direction is fully confirmed by the fact that the strong m6(2)A(1)-m6(2)A(2) stack in the tetranucleotide m6(2)A-m6(2)A-U-m6(2)A virtually precludes the formation of bulged-out structures. The conformational characteristics of dA-r(U-A) appear comparable with those of m6(2)A-U-m6(2)A, which indicates that the presence of a 2'-hydroxyl group in the first purine residue is not a necessary prerequisite for the formation of a bulge.

Conformational analysis of the dinucleotides 5'-methylphospho-N6-dimethyladenylyl-uridine (mpm62A-U) and 5'methylphospho-uridylyl-N6-dimethyladenosine (mpU-m62A) and of the trinucleotide U-m62A-U. A nuclear magnetic resonance and circular dichroic study

NMR and CD studies were carried out on the dinucleotides 5'-methylphospho-N6-dimethyladenylyl-uridine (mpm62-U) and 5'-methylphospho-uridylyl-N6-dimethyladenosine (mpU-m62A) and on the trinucleotide U-m62A-U. A detailed comparison is given of the conformational features of mpm62A-U and mpU-m62A with the corresponding 5'-nonphosphorylated dinucleotides m62A-U and U-m62A, respectively. The behaviour of the trinucleotide U-m62A-U is compared with the properties of the constituent dinucleotides U-m62A and mpm62A-U. Chemical-shift and CD data were used to determine the amount of stacking interactions. For each compound NMR spectra were recorded at two or three sample concentrations in order to separate intermolecular and intramolecular base-base interactions. The coupling constants of the ribose ring are interpreted in terms of the N/S equilibrium, and population distributions along the backbone angles beta, gamma and epsilon are presented. The combined data indicate a strong similarity between mpm62A-U and m62A-U both in degree and in mode of stacking. In contrast, the existence of different types of stacking interactions in mpU-m62A and U-m62A is suggested in order to explain the NMR and CD data. It is concluded that dinucleoside bisphosphates serve better as a model for the behaviour of trinucleotides than dinucleoside monophosphates. The trinucleotide U-m62A-U adopts a regular single-stranded stacked RNA structure with preference for N-type ribose and gamma+ and beta t backbone torsion angles. The difference in behaviour between the U-m62A- part of U-m62A-U and the dimer U-m62A is seen as a typical example of conformational transmission.

Conformational analysis of the deoxyribofuranose ring in DNA by means of sums of proton-proton coupling constants: a graphical method

A graphical method is presented for the conformational analysis of the sugar ring in DNA fragments by means of proton-proton couplings. The coupling data required for this analysis consist of sums of couplings, which are referred to as sigma 1' (= J1'2' + J1'2''), sigma 2' (= J1'2' + J2'3' + J2'2''), sigma 2'' (= J1'2'' + J2''3' + J2'2'') and sigma 3' (= J2'3' + J2''3' + J3'4'). These sums of couplings correspond to the distance between the outer peaks of the H1', H2', H2'' and H3' [31P] resonances, respectively, (except for sigma 2' and sigma 2'' in the case of a small chemical shift difference between the H2' and H2'' resonances) and can often be obtained from 1H-NMR spectra via first-order measurement, obviating the necessity of a computer-assisted simulation of the fine structure of these resonances. Two different types of graphs for the interpretation of the coupling data are discussed: the first type of graph serves to probe as to whether or not the sugar ring occurs as a single conformer, and if so to analyze the coupling data in terms of the geometry of this sugar ring. In cases where the sugar ring does not occur as a single conformer, but as a blend of N- and S-type sugar puckers, the second type of graph is used to analyze the coupling data in terms of the geometry and population of the most abundant form. It is shown that the latter type of analysis can be carried out on the basis of experimental values for merely sigma 1',sigma 2' and sigma 2'', without any assumptions or restrictions concerning a relation between the geometry of the N- and S-type conformer. In addition, the question is discussed as to how insight can be gained into the conformational purity of the sugar ring from the observed fine structure of the H1' resonance. Finally, a comparison is made between experimental coupling data reported for single-stranded and duplex DNA fragments and covalent RNA-DNA hybrids on the one hand and the predicted couplings and sums of couplings presented in this paper on the other hand.

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.

13C-NMR of ribosyl ApApA, ApApG and ApUpG

The chemical shifts as well as the 13C-31P coupling constants of the carbon-13 nuclei in single-stranded ApApA, ApApG, and ApUpG are sensitive to sequence and temperature. ApApA and ApApG have similar properties with large shielding (up to 1.7 ppm) of many of the base carbons upon decreasing the temperature from 70 degrees C to 11 degrees C; the base carbons have smaller shielding changes in ApUpG. Large shielding and deshielding effects are observed for the 1', 3', 4' and 5'-carbons over this temperature range. Analysis of the 13C-31P couplings measured at the 4' ribose carbons show that the population of the anti rotamer about O5'-C5' varies from 98 to 75%, and is higher in ApApA and ApApG than in ApUpG. The CCOP coupling data at 2' and 4' is consistent with a blend of the -antiperiplanar/-synclinal nonclassical rotamers about the C3'-O3' bond, varying from 89/11% in ApApG to 55/45% in ApUpG. The coupling and chemical shift data support the thesis that ApUpG is stacked much less than the other two molecules. The stacked forms of all three trinucleotides is most easily interpreted by a standard A-RNA model. It is not necessary to invoke the "bulged base" hypothesis [Lee, C.-H. and Tinoco, Jr., I. (1981) Biophysical Chemistry 1, 283-294; Lankhorst, P.P., Wille, G., van Boom., J.H., Altona, C., and Haasnoot, C.A.G. (1983) Nucleic Acids Research 11, 2839-2856] to explain the contrast in 13C spectroscopic properties of ApUpG in comparison to ApApG and ApApA.

Sequence-dependent structural variation in single-helical DNA. Proton NMR studies of d(T-A-T-A) and d(A-T-A-T) in aqueous solution

The two deoxyribotetranucleoside triphosphates d(T-A-T-A) and d(A-T-A-T) were investigated in aqueous solution by one- and two-dimensional proton NMR at 300 and 500 MHz. It is demonstrated that both compounds occur predominantly in the single-helical form. Accurate coupling constants are obtained by computer simulation of several 500-MHz spectra. The data are interpreted in terms of N and S pseudorotational ranges. The geometry of the major S-type conformers displays a clear sequence dependence, as expressed by variation of the endocyclic backbone angle delta (C5'-C4'-C3'-O3'). A simple sum rule is proposed to predict delta variation in single-helical DNA fragments. Comparisons are made with other sequence-dependent geometries as observed in a double-helical B-DNA fragment in the crystalline state. Furthermore, one- and two-dimensional nuclear Overhauser effect (NOE) spectroscopy was carried out on d(T-A-T-A). An inventory is made of the observed intra- and inter-residue NOEs. The NOE data confirm the presence of a highly stacked single-helical conformation of d(T-A-T-A) in solution. No indications are found for the formation of a bulge-out structure as observed for analogous alternating purine-pyrimidine oligoribonucleotides.

Conformational studies of trinucleoside bisphosphates. 2. Potential energy calculations

Semi-empirical potential energy calculations were carried out for the conformations of ApApA, a trinucleoside bisphosphate (trimer), as proposed in the preceding paper in this journal. Energy minimals were obtained for these conformers, i. e. conformers I-I, I-II, III-I' and I'-I' (the first Roman numeral represents the conformation of the ApAp- moiety and the second one the -pApA moiety of ApApA), of the nearest-neighbor interactions, and conformers B1, B2 and B3 of bulged structures. The results show that the detailed conformations of the component dinuclotide (dimer) moieties in I-I, I-II, III-I' and I'-I' are close (+/- 10 degrees) to the corresponding stable conformations of dinucleotides in solution proposed by Lee and coworkers. The bulged conformations B1, B2 and B3 have open forms for each dinucleotide moiety. Their preferred local conformations are anti, gamma +, beta +, epsilon + (except for Ap- in B3, epsilon-), and 3' endo (except for Ap- in B3, 2' endo) for the glycosidic bond, C4'-C5', C5'-O5', C3'-O3' and the ribose moiety, respectively. The conformation of the phosphodiester linkage in -pApA moiety is zeta t2, alpha t3 for all the three bulged conformers, while that in ApAp- is zeta- 1, alpha t2 for B1, zeta+ 1 alpha- 2 for B2, and zeta+ 1, alpha+ 2 for B3. The calculated composite chemical shifts of the base protons from ring-current effects agree with the corresponding observed trimerization shifts. Conformers I-I and I'-I' are helix-like structures and correspond to the right-handed and the 'vertical-stacked' helices, respectively. The structure of I-II is similar to that observed in the ApApA single crystals. The bulged conformations may serve as models for frameshift mutations.

Conformational studies of 13 trinucleoside bisphosphates by 360-MHz 1H-NMR spectroscopy. 1. Ribose protons

The ribose protons of 13 trinucleoside bisphosphates (trimers) were studied, using 360-MHz proton nuclear magnetic resonance spectroscopy. Complete assignments and analyses of the NMR signals of these protons were carried out by the methods of homonuclear decoupling and computer line-shape simulations. It was shown that the trinucleotides preferred the anti, 3' endo, gamma +, beta t and epsilon t/epsilon- conformations for the glycosidic torsions, the ribose rings, the C4'-C5' bonds, the C5'-O5' bonds, and the C3'-O3' bonds, respectively. It was also found that the trimers, especially those which had noticeable population of 'bulged' structures, did not necessarily have a higher population of these preferred local conformations than their component dimers. The overall conformations of the trinucleotides are classified into two categories. The conformations in the first category involve the nearest-neighbor interactions. Each dinucleotide moiety can assume one of the four stable conformations (I, I', II and III) or the open forms of dinucleoside monophosphates. However, due to steric hindrance, there are only four cases in which both dinucleotide moieties can assume one of the four stable conformations at the same time. These four combinations of conformations are I-I, I'-I', I-II and III-I', where the first Roman numeral represents the conformation of the NpN'p-moiety and the second one, that of the -pN'pN'' moiety of the trimers. Among them, I-I and I'-I' are helical structures, capable of forming a double helix. The second category contains conformations with bulged structures which have the two dinucleotide moieties in open forms (i.e. no nearest-neighbor interactions) and the bases of the two terminal residues stacking on each other while the middle residue is bulged out. These bulged conformations may serve as structural models for frame-shift mutations.

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

The double helical structure of the self-complementary DNA-RNA-DNA hybrid d(CG)r(CG) d(CG) was studied in solution by 500 MHz1H-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 1H-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 accommodated 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%).

cis-Platinum induced distortions in DNA. Conformational analysis of d(GpCpG) and cis-pt(NH3)2[d(GpCpG)], studied by 500-MHz NMR

Proton NMR studies at 500 MHz in aqueous solution were carried out on the G-G chelated deoxytrinucleosidediphosphate platinum complex cis-Pt(NH3)2[d(GpCpG], on the uncoordinated trinucleotide d(GpCpG) and on the constituent monomers cis-Pt(NH3)2[d(Gp)]2, cis-Pt(NH3)2[d(pG)]2, d(Gp), d(pCp) and d(pG). Complete NMR spectral assignments are given and chemical shifts and coupling constants are analysed to obtain an impression of the detailed structure of d(GpCpG) and the distortion of the structure due to chelation with [cis-Pt(NH3)2]2+. Platination of the guanosine monophosphates affects the sugar conformational equilibrium to favour the N conformation of the deoxyribose ring. This feature is also apparent in ribose mononucleotides and is possibly caused by an increased anomeric effect. In cis-Pt(NH3)2[d(pG)]2 the phase angle of pseudorotation of the S-type sugar ring is 20 degrees higher than in 'free' d(pG) which might be an indication for an ionic interaction between the positive platinum and the negatively charged phosphate. It appears that d(GpCpG) reverts from a predominantly random coil to a normal right-handed B-DNA-like single-helical structure at lower temperatures, whereas the conformational features of cis-Pt(NH3)2[d(GpCpG)] are largely temperature-independent. In the latter compound much conformational freedom along the backbone angles is seen. The cytosine protons and deoxyribose protons exhibit almost no shielding effect as should normally be exerted by the guanine bases in stacking positions. This is interpreted in terms of a 'turning away' of the cytosine residue from both chelating guanines. Conformational features of cis-Pt(NH3)2[d(GpCpG)[ are compared with the 'bulge-out' of the ribose-trinucleotide m6(2)ApUpm6(2)A.

Conformational analysis of the single-stranded ribonucleic acid A-A-C-C. A one-dimensional and two-dimensional proton NMR study at 500 MHz

Proton NMR studies at 300 MHz and 500 MHz are reported on the ribotetranucleotide A-A-C-C. The complete 1H-NMR spectral assignment at 20 degrees C is given. Two-dimensional NMR was used to elucidate spin multiplets in 'crowded' regions. Nuclear Overhauser enhancement (NOE) experiments made an unambiguous spectral assignment possible and yielded information on interproton distances. The N/S equilibrium of the riboses and the rotamer populations around some backbone torsion angles are presented. A large preference for N-type ribose and gamma+ and beta t backbone torsions is observed, in particular in the central A-C unit of A-A-C-C. Information on distances between protons of different nucleotide units, obtained from NOE experiments, constitutes a probe of base-base stacking. It is concluded that A-A-C-C offers a good model for RNA single-strand conformation.

Thermodynamics of stacking and of self-association of the dinucleoside monophosphate m2(6)A-U from proton NMR chemical shifts: differential concentration temperature profile method

Chemical shifts of base and sugar protons of the modified ribodinucleoside monophosphate N6-dimethyladenylyl(3'-5')uridine (m2(6)A-U) were measured at 100, 360 and 400 MHz in aqueous solution. Seven different samples were used with concentrations ranging from 0.28 mM to 32.7 mM. The temperature was varied from -5 degrees C to 105 degrees C. An internal temperature calibration was used. The effects of intermolecular self-association and of intramolecular stacking on the chemical shifts were quantitatively separated by means of a new approach: differential concentration/temperature profiles (DCTP). Several computational models were tested and the analysis allowed deeper insight into the behaviour of m2(6)A-U at the molecular level. The simple two-state approach for both self-association and stacking already afforded a significant improvement over models in which the association is entirely neglected. A computer least-squares analysis of the chemical shift behaviour of each individual proton yielded thermodynamic parameters for self-association and stacking. However, the two-state model did not suffice to reproduce accurately all of the observations. A satisfactory fit required two additional assumptions: (a) the aromatic protons experience different association shifts in stacked and in unstacked molecules: (b) a temperature-dependent conformational equilibrium exists between sets of unstacked microstates. The stacked state is taken to represent a single conformational species. The implementation of this extended model in the least-squares optimization allowed the reproduction of over one thousand chemical shift observations within experimental error. Thermodynamic equilibrium parameters deduced for intramolecular stacking are: delta H degrees x = -28.8 kJ mol-1, delta S degrees x = -93 J mol-1 K-1. These numbers agree well with those obtained earlier by us from circular dichroism spectra. The equilibrium enthalpy and entropy values deduced for the association process are: delta H degrees A = -35 kJ mol-1 and delta S degrees A = -95 J mol-1 K-1.