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

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein–nucleic acid complexes.

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker. J-couplings were calculated in 16 distinct nucleic acid conformations, two principal double-helical DNAs, B- and A-, the main RNA form, A-RNA, as well as in 13 other RNA conformations. High-level quantum mechanics calculations used a baseless dinucleoside phosphate as a molecular model, and the effect of solvent was included. The predicted J-couplings correlate reliably with available experimental data from the literature.

Thermodynamics of melting of the circular dumbbell d<pCGC-TT-GCG-TT>

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson-Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5'-d<pCGC-TT-GCG-TT>-3' (I) and reference compounds 5'-d<pGGC-TT-GCC-TT>-3' (II) and 5'-d(GCG-TC-CGC)-3' (III) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin loop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4. The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson-Crick C-G base pairs. The high-temperature form, L2L4, contains a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop, with the cytosine base placed in a syn orientation. The opposite 5'-CTTG-3' loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5'-GTTC-3' loop as found in I; a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point Tm of II was estimated to lie above 365 K. The Tm value of the duplex stem and the 5'-CTTG-3' loop of the L2L4 form of I is 352 +/- 2 K. The delta H0 is calculated as -89 +/- 10 kJ/mol. The Tm value determined for the individual residues of the 5'-GTTC-3' loop lies 4 degrees-11 degrees lower. The enthalpy delta H0 of melting the thymine residues in the 5'-GTTC-3' loop is calculated to be -61 +/- 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for delta H0 of +30 kJ/mol in going from a two-membered to a four-membered loop, in agreement with the enthalpy difference of +28 kJ/mol between the two loops at the Tm midpoint. For hairpin III the upper limit for delta H0 in going from a two-membered to a four-membered loop amounts to +/- 21 kJ/mol.(ABSTRACT TRUNCATED AT 400 WORDS)

Comparative NMR analysis of the decadeoxynucleotide d-(GCATTAATGC)2 and an analogue containing 2-aminoadenine

The dodecadeoxynucleotide duplex d-(GCATTAATGC)2 has been prepared with all adenine bases replaced by 2-NH2-adenine. This modified duplex has been characterized by nuclear magnetic resonance (NMR) spectroscopy. Complete sequence-specific 1H resonance assignments have been obtained by using a variety of 2D NMR methods. Multiple quantum-filtered and multiple quantum experiments have been used to completely assign all sugar ring protons, including 5'H and 5'H resonances. The assignments form the basis for a detailed comparative analysis of the 1H NMR parameters of the modified and parent duplex. The structural features of both decamer duplexes in solution are characteristic of the B-DNA family. The spin-spin coupling constants in the sugar rings and the relative spatial proximities of protons in the bases and sugars (as determined from the comparison of corresponding nuclear Overhauser effects) are virtually identical in the parent and modified duplexes. Thus, substitution by this adenine analogue in oligonucleotides appears not to disturb the global or local conformation of the DNA duplex.

1H-, 13C-, 31P-NMR studies and conformational analysis of NADP+, NADPH coenzymes and of dimers from electrochemical reduction of NADP+

All H,H, H,P and several C,P coupling constants, including those between C-4' and the vicinal phosphorus atom, have been determined for NADP+, NADPH coenzymes and for a 4,4-dimer obtained from one-electron electrochemical reduction of NADP+. From these data the preferred conformation of the ribose, that of the 1,4-dihydronicotinamide rings, and the conformation about bonds C(4')-C(5') and C(5')-O(5') were deduced. The preferred form of the 1,4- and 1,6-dihydropyridine rings and the conformation about the ring-ring junction were also obtained for all the other 4,4- and 4,6-dimers formed in the same reduction. All the dimers show a puckered structure, i.e., a boat form for the 1,4- and a twist-boat for the 1,6-dihydronicotinamide ring; both protons at the ring-ring junctions are equatorial and have preferred gauche orientation. On the contrary, the reduced coenzyme NADPH displays a planar or highly flexible conformation, rapidly flipping between two limiting boat structures. The conformation of the ribose rings, already suggested for the NADP coenzymes to be an equilibrium mixture of C(2')-endo (S-type) and C(3')-endo (N-type) puckering modes, has been reexamined by using the Altona procedure and the relative proportion of the two modes has been obtained. The S and N families of conformers have almost equal population for the adenine-ribose, whereas for the nicotinamide-ribose rings the S-type reaches the 90%. The rotation about the ester bond C(5')-O(5') and about C(4')-C(5'), defined by torsion angles beta and gamma respectively, displays a constant high preference for the trans conformer beta t (75-80%), whereas the rotamers gamma are spread out in a range of different populations. The values are distributed between the gauche gamma + (48-69%) and the trans gamma t forms (28-73%). The gamma + conformer reaches a 90% value in the case of NADP+ and NMN+. The conformations of the mononucleotides 5'-AMP, NMN+ and NMNH were also calculated from the experimental coupling constant values of the literature.

Conformational consequences of the incorporation of arabinofuranosylcytidine in DNA. An NMR study of the DNA fragments d(CGCTAGCG) and d(CGaCTAGCG) in solution

The self-complementary octamers d(CGCTAGCG) and d(CGaCTAGCG) (aC, arabinofuranosylcytidine) were studied by means of NMR spectroscopy. It is shown that d(CGaCTAGCG), under suitable conditions of oligonucleotide concentration, ionic strength and temperature, exclusively adopts a hairpin structure. However, under the same experimental conditions (5 mM DNA, no added salt, 295 K) d(CGCTAGCG) mainly adopts a B-DNA-type duplex. At lower temperatures (less than or equal to 290 K) the hairpin form of d(CGaCTAGCG) occurs in slow exchange with an intact B-DNA-type duplex. When the DNA concentration of d(CGCTAGCG) is dramatically reduced (less than or equal to 0.5 mM) the hairpin form becomes highly favoured at the expense of the dimer. Moreover, proton-chemical-shift considerations indicate that the structural features of the hairpin structure of d(CGCTAGCG) mimic, in part, those of the modified octamer d(CGaCTAGCG), i.e. a loop comprising only the two central residues with the thymine located into the minor groove (Pieters, J. M. L., de Vroom, E., van der Marel, G. A., van Boom, J. H., Koning, T. M. G., Kaptein, R. and Altona, C. unpublished results). Thermodynamic analysis of d(CGCTAGCG) yields an average Tmd value of 342 K (1 M DNA) and a delta Hod value of -266 kJ/mol for the dimer/coil transition and an average Tmh value of 321 K and delta Hoh - 102 kJ/mol for the hairpin/coil equilibrium. For the duplex/coil equilibrium of d(CGaCTAGCG) an average Tmd value of 336 K (1 M DNA) and delta Hod value of -253 kJ/mol are deduced. The hairpin/coil transition of d(CGaCTAGCG) is characterized by a delta Hoh value of -104 kJ/mol and an average Tmh value of 331 K. It is concluded that incorporation of an arabinofuranosylcytidine in the octamer d(CGaCTAGCG) results in stabilization of the hairpin form, whereas the dimer is destablized by two aC.dG base pairs.

Conformational and thermodynamic consequences of the introduction of a nick in duplexed DNA fragments: an NMR study augmented by biochemical experiments

NMR studies were carried out on various equimolar mixtures consisting of a combination of oligomers: d(ACGGCT) (I). d(pACGGCT) (Ia), d(TGCAGT) (II), d(AGCCGTACTGCA) (III), d(TGCAGTACGGCT) (IV). It is shown that I + II + III (MI) and Ia + II + III (M2) form stable duplexes with nicks in the centre of the respective double helices. A close analysis of the NOESY experiments of M1 and M2 revealed that these fragments form B-DNA type duplex structures. A comparison of the chemical-shift data of the nicked duplexes with those of the intact duplex of III + IV (M3) demonstrated that only small local distortions occur when a nick is introduced. The chemical-shift profiles of M1 and M3 were used to obtain the thermodynamic data for the duplex/coil transitions. The profiles of M1 were analysed by means of a new thermodynamic model (TRIDUP). From the calculated thermodynamic data of M1 and M3 it is concluded that the melting behaviour of M1 occurs cooperatively. A ligation experiment demonstrated that the relatively small substrate (M2) was almost completely joined after an overnight incubation at 14 degrees C.

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.

Alterations in the d(CpGpT) structure in solution as a result of [PtCl(diethylenetriamine)]+ binding

The trinucleotide d(CpGpT) reacts with [PtCl(dien)]Cl (dien = diethylenetriamine) to yield as a single adduct Pt(dien)[d(CpGpT)-N7(2)]. The structure of this adduct in solution has been analysed with the aid of NMR spectroscopy and compared with that of the unmodified trinucleotide. A change in the population of the S conformer of the guanosine deoxyribose ring and a syn preference of the guanine residue are the most important changes occurring upon platination. As a result the dC-dG stack disappears, whereas the dG-dT stack is hardly affected. The CD spectra of both platinated and free d(CpGpT) confirm the different nature of the two molecules.

Bulge-out structures in the single-stranded trimer AUA and in the duplex (CUGGUGCGG).(CCGCCCAG). A model-building and NMR study

Model-building studies were carried out on the trimer AUA. Bulge-out structures which allow incorporation into a continuous RNA helix were generated and energy-minimized. All geometrical features obtained by previous NMR studies on purine-pyrimidine-purine sequences are accounted for in these models. One of the models was used to fit into a double helical fragment. Only minor changes were necessary to construct a central bulge-out in an otherwise intact duplex. NMR and model-building studies were performed on the duplex (CUGGUGCGG).(CCGCCCAG) which contains an unpaired uridine residue. NOE data, chemical-shift profiles and imino-proton resonances provided evidence that the extra U is bulged out of the duplex. The relatively small dispersion in 31P chemical shifts (approximately equal to 0.7 ppm) indicate the absence of t/g or g/t combinations for the phosphodiester angles zeta/alpha. An energy-minimized model of the duplex, which fits the present collection of data, is presented.

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.

Synthesis and proton-NMR studies of oligonucleotides containing an apurinic (AP) site

In order to elucidate the conformational properties of base-deleted oligodeoxyribonucleotides, the molecules d-CpS(pCpG)n (n = 1,2; S = sugar) were synthesized by the phosphotriester method and characterized by 1H-NMR spectroscopy. Complete assignment of all non-exchangeable proton resonances of both compounds was obtained by 1D- and 2D-NMR techniques. In combination with computer simulation, these spectra yielded proton-proton and proton-phosphorus coupling constants of high accuracy. These data provide valuable information about the sugar and the backbone conformation. It appears that d-Cp1Sp2Cp3G4 does not form a duplex under any of the conditions studied. On the contrary, the base-deleted hexamer d-Cp1Sp2Cp3Gp4Cp5G6 occurs as a right-handed' staggered' DNA duplex at 280 K: the core of this duplex is formed by the residues C(3)-G(6); two 'dangling' residues C(1) and S(2) are located at the two 5'-ends of the duplex. The assignment of the corresponding imino proton resonances for [d-CpS(pCpG)2]2 was based on their thermal behavior: the line broadening of these resonances was studied as a function of temperature. The chemical shift and the number of imino proton resonances accord well with the number and type of Watson-Crick base pairs which can be formed in the staggered duplex described above. Thermodynamic parameters of duplex formation were obtained from an analysis of the chemical shift versus temperature profiles of aromatic base and H-1' protons. It is suggested that the cytosine ring of C(1) stacks, at least part of the time, with the guanine ring on the nucleotide residue, G(6), situated in the complementary strand. The binding of Lys-Trp-Lys to [d-CpS(pCpG)2]2 as well as to [d-CpGpCpG]1 was investigated. It is concluded that the indole ring of the tryptophan residue probably stacks on top of the 3'-terminal guanine base of both duplexes, but not on the nucleic acid bases next to the apurinic (AP) site.

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.

Influence of N6-methylation of residue A(5) on the conformational behaviour of d(C-C-G-A-A-T-T-C-G-G) in solution studied by 1H-NMR spectroscopy. 1. The duplex form

One- and two-dimensional NMR studies at 300 MHz and 500 MHz were carried out on the two oligonucleotides d(C-C-G-A-A-T-T-C-G-G) and d(C-C-G-A-m6A-T-T-C-G-G) in aqueous solution. NMR spectra were observed at 10 mM sample concentration over the temperature range 273-368 K. Assignments are given of the base, H1', H2', H2", H3' and of some H4' resonances, based upon a combination of two-dimensional correlation spectra (COSY) and two-dimensional nuclear Overhauser effect spectra (NOESY); imino-proton resonances were assigned with the aid of a two-dimensional NOE experiment. Chemical shift vs temperature profiles were constructed in order to gain insight into the influence of N6-methylation of residue A(5) on the temperature-dependent conformational behaviour of the decamer and to determine thermodynamic parameters for the duplex-to-coil transition. The NOESY spectra, the imino-proton spectra and the shift profiles of the two compounds, under conditions where each forms a B-DNA-type duplex, are very similar. This is taken to indicate that the influence of N6-methylation of residue A(5) on the local structure of the duplex must be small. However, the temperature dependence of the (non-)exchangeable proton resonances of the two compounds reveals that methylation slows down the duplex-single-strand exchange. Furthermore, a thermodynamic analysis of the two compounds indicates that N6-methylation slightly decreases the stability of the duplex relative to the monomeric forms (Tm is reduced from 332 K down to 325 K at 10 mM sample concentration). Proton-proton couplings were obtained by means of one-dimensional and two-dimensional NMR experiments and were used in a conformational analysis of the sugar ring of each residue of the two compounds in the duplex form. The analysis indicated that all sugar rings display conformational flexibility in the intact duplex: population S-type sugar conformation ranges from 70% to 100%. A more refined analysis of the sugar rings of the parent compound revealed a sequence-dependent variation of the sugar geometry. This variation does not follow well the trend predicted by the Calladine/Dickerson sigma 3-sum rule [Dickerson, R. E. (1983) J. Mol. Biol. 166, 419-441; Calladine, C. R. (1982) J. Mol. Biol. 161, 343-352]; moreover the actual variations appear to be smaller in solution than those expected on the basis of known X-ray structures.

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.

Conformational analysis of the B and Z forms of the d(m5C-G)3 and d(br5C-G)3 hexamers in solution. A 300-MHz and 500-MHz NMR study

The B and the Z forms of the DNA hexamers d(m5C-G)3 and d(br5C-G)3 were investigated by means of NMR spectroscopy. It is demonstrated that the low-salt form of d(m5C-G)3 is a B DNA structure. The form, which becomes increasingly predominant when increasing amounts of MgCl2 and/or methanol are added to the solution, has Z DNA characteristics. It is shown that the major geometrical features of the Z form of d(m5C-G)3 in the crystal structure are maintained in solution, with the dC residues S sugar conformation, gamma + and the base in the anti orientation and the dG residues N (except the 3'-terminal residue), gamma t and syn. Neither the Z form of the methylated nor that of the brominated compound resembles the Z' form, in which the deoxy guanosine sugar rings adopt a C1'-exo conformation. Substitution of m5C by br5C causes no perceptible conformational changes in either the B or in the Z forms.

Conformational analysis of the octamer d(G-G-C-C-G-G-C-C) in aqueous solution. A one-dimensional and two-dimensional proton NMR study at 300 MHz and 500 MHz

Proton NMR studies in H2O and 2H2O were carried out on the self-complementary DNA octamer d(G-G-C-C-G-G-C-C) and compared with similar studies on the dimethylated analogue d(G-G-m5C-m5C-G-G-C-C) [Sanderson, M. R., Mellema, J.-R., van der Marel, G. A., Wille, G., van Boom, J. H. & Altona, C. (1983) Nucleic Acids Res. 11, 3333-3346]. Base, H1', H2' and H2" resonances were assigned by means of two-dimensional nuclear Overhauser spectroscopy (NOESY) and correlated spectroscopy (COSY) experiments. Chemical shift vs temperature profiles were used to analyze the temperature-dependent conformational behaviour and to extract thermodynamic parameters for the duplex-to-coil transition. Analysis of proton-proton couplings were used to discriminate between J1'2' and J1'2" and between the H2' and H2" resonances as well as to obtain conformational parameters of the sugar rings. From the NOESY and COSY experiments, the temperature profiles and the analysis of the coupling data it is concluded that: (a) the compound adopts a B-DNA-type helix in solution; (b) the sugar rings in the intact duplex display limited conformational freedom; (c) methylation of the cytosine bases at the C5 position has no measurable effect on the conformational behaviour of the octamer, nor on the conformation of the sugar rings; however, methylation does increase the stability of the duplex in aqueous solution under conditions of low salt concentration.

Influence of 2-aminoadenosine, A', on the conformational behaviour of d(T-A'-T-A'). A one-dimensional proton NMR study at 300 MHz and 500 MHz

Proton NMR studies at 300 MHz and 500 MHz in aqueous solution were carried out on the modified self-complementary tetranucleoside trisphosphate d(T-A'-T-A'), in which d(A') represents 2-aminodeoxyadenosine. The NMR spectra were observed at two sample concentrations over the temperature range 2-70 degrees C. Assignments, based on homonuclear decoupling and nuclear Overhauser enhancement (NOE) experiments, are given. The concentration dependence of the chemical shift vs temperature profiles was used to extract information concerning duplex formation. The proton-proton and proton-phosphorus coupling constants were obtained at four temperatures and yielded accurate conformational data on the sugar ring and on the back-bone angles beta, gamma and delta. From the observed line broadening, shift profiles, NOEs, and from the observation of imino proton resonances, it is concluded that the compound exists as a miniduplex at low temperature. Comparison of these observations with similar observations on the parent compound d(T-A-T-A) indicate that substitution of dA by dA' increases the tendency towards duplex formation. At low temperature the compound adopts a stacked B-DNA type structure: destacking occurs on raising the temperature. The sequence-dependent sugar ring geometry [Mellema, J.-R., Pieters, J. M. L., van der Marel, G. A., van Boom, J. H., Haasnoot, C. A. G. & Altona, C. (1984) Eur. J. Biochem. 143, 285-301] is present in this molecule. The conformational parameters of d(T-A'-T-A') and d(T-A-T-A) are quite similar, thus substitution of dA by dA' has no measurable influence on the geometry of the sugar ring or on the backbone angles beta, gamma and delta.

Conformational characteristics of the hexanucleoside pentaphosphate AUAUAU: a 2D NMR study at 500 MHz

All 36 ribose proton resonances and most of the base proton resonances of the hexanucleoside pentaphosphate AUAUAU have been assigned unequivocally using 2D J-resolved spectroscopy, spin echo correlated spectroscopy (SECSY) and 2D NOE spectroscopy (NOESY). The NMR parameters of AUAUAU are compared with those of smaller fragments that contain methylated adenine bases: m62AU, m62AUm62A, m62AUm62AU and m62AUm62AUm62A. Previous studies on this series of compounds have shown that in all these cases purine-pyrimidine-purine sequences prefer to adopt a mixture of states which have as common feature that the interior pyrimidine residues are bulged out, whereas the purine residues stack upon each other. Chemical shift data, proton-proton coupling constants, as well as the observation of imino-proton resonances for AUAUAU show unambiguously that upon lowering the temperature the high-temperature "bulged out" situation reverts to a normal A-RNA-like double helix.

A 1H and 31P NMR study of cis-Pt(NH3)2[d(CpGpG)-N7(2),N7(3)]. The influence of a 5'-terminal cytosine, on the structure of the cis-Pt(NH3)2[d(GpG)-N7,N7] intrastrand cross-link

Proton NMR studies at 300 MHz and 500 MHz have been carried out on the trinucleoside bisphosphate d(CpGpG) and on cis-Pt(NH3)2[d(CpGpG)-N7(2),N7(3)] [abbreviated as d(CpGpGp) . cisPt]. For the Pt adduct, 13C and 31P NMR was also used for characterizing the oligonucleotide. d(CpGpG) appears to revert to a B-DNA-type single helix at lower temperatures. The relatively small concentration dependence of the proton chemical shifts, in comparison with shifts due to intramolecular stacking effects, indicates that the compound is essentially single-stranded. In d(CpGpGp) . cisPt, the first nucleoside, C(1), stacks well on top of the second, G(2), despite the N conformation of the G(2) sugar ring. The platinated GpG part in this trimer adopts largely the same structure as in cis-Pt(NH3)2[d(GpGpG)-N7(1),N7(2)] [den Hartog, J. H. J., et al. (1982) Nucleic Acids Res. 10, 4715-4730]. Main differences however, are changes in H8 chemical shifts and a 0.6-ppm downfield shift of the third nucleotide phosphorus, P(3), in d(CpGpGp) . cisPt with respect to P(2) in d(GpG) . cisPt. The latter shift change is likely to be induced by a structural alteration, caused by stacking of C(1) on top of G(2). Also, the large chemical shift differences between the two H8 protons in d(NpGpG) . cisPt fragments is discussed; the deviation from a mirror symmetry of the two guanine bases seems to be the main origin of this effect. The chemical shift changes, observed in the proton and phosphorus NMR chemical shift temperature and chemical shift pH profiles have been explained in terms of stack-destack equilibria changes.

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 4. The torsion angle distribution about the C3'-O3' bond in DNA constituents

Carbon-13 and proton NMR spectra of a series of oligodeoxynucleotides (d(CT), d(CC), d(TA), d(AT), d(CG), d(GC), d(AG), d(AAA), d(TATA) and d(GGTAAT] were measured at various temperatures. The three coupling constants that are related to the magnitude of backbone angle epsilon (J(C4'-P), J(C2'-P) and J(H3'-P] are analyzed in terms of a three-state equilibrium about this bond. Two epsilon (trans) angles occur, which differ in magnitude depending on the conformation (N or S) of the adjoining deoxyribose ring. The S-type deoxyribose ring is associated with a smaller epsilon (trans) angle: epsilon (t,S) = 192 degrees. The N-type deoxyribose ring is associated with a larger epsilon (trans) angle epsilon (t,N) = 212 degrees. The third rotamer participating in the conformational equilibrium, is a gauche(-) (epsilon (-] conformer and occurs exclusively in combination with the S-type sugar ring (epsilon (-,S) = 266 degrees). Within the limits of experimental error, the magnitude of these three angles appears to be independent of the particular base sequence, except in the case of d(CG) where a slightly larger epsilon (t,S) angle (197 degrees) is indicated. A simple equation is proposed which may be used to calculate the population of epsilon (t,S) conformer in cases where only J(H3'-P) is known.