The 31P-NMR spectrum of the dodecamer d(GACGATATCGTC)

Solution structure of the octamer motif in immunoglobulin genes via restrained molecular dynamics calculations

The solution structure of the DNA decamer d(CATTTGCATC)-d(GATGCAAATG), comprising the octamer motif of immunoglobulin genes, is determined by restrained molecular dynamics (rMD) simulations. The restraint data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were previously obtained by a complete relaxation matrix analysis of the two-dimensional nuclear Overhauser enhancement (2D NOE) intensities and by the quantitative simulation of cross-peaks in double-quantum-filtered correlated (2QF-COSY) spectra. The influence of torsion angles and the number of experimental distance restraints on the structural refinement has been systematically examined. Omitting part of the experimental NOE-derived distances results in reduced restraint violations and lower R factors but impairs structural convergence in the rMD refinement. Eight separate restrained molecular dynamics simulations were carried out for 20 ps each, starting from either energy-minimized A- or B-DNA. Mutual atomic root-mean-square (rms) differences among the refined structures are well below 1 A and comparable to the rms fluctuations of the atoms about their average position, indicating convergence to essentially identical structures. The average refined structure was subjected to an additional 100 ps of rMD simulations and analyzed in terms of average torsion angles and helical parameters. The B-type duplex exhibits clear sequence-dependent variations in its geometry with a narrow minor groove at the T3.A3 tract and a large positive roll at the subsequent TG.CA step. This is accompanied by a noticeable bend of the global helix axis into the major groove. There is also evidence of significant flexibility of the sugar-phosphate backbone with rapid interconversion among different conformers.

1H- and 31P-NMR studies of ditercalinium binding to a d(GCGC)2 and d(CCTATAGG)2 minihelices: a sequence specificity study

The structures of the complexes formed in aqueous solution between ditercalinium, a bis-intercalating drug, and both the self-complementary tetranucleotide d(GCGC)2 and octanucleotide d(CCTATAGG)2, have been investigated by 400-MHz 1H-nmr and 162-MHz 31P-nmr. All the nonexchangeable protons, as well as the exchangeable imino protons and the phosphorus signals, have been assigned. Both oligonucleotides have been shown to adopt a right-handed B-DNA type structure. The addition of ditercalinium to the oligonucleotides lead to the formation of complexes in slow exchange at the nmr time scale with the free helices. At all drug-to-helix ratios studied, the ditercalinium was found in the bound form, whereas free and complexed oligonucleotides were in slow exchange, allowing resonance assignments through two-dimensional chemical exchange experiments. for d(GCGC)2 the strong upfield shifts induced on all aromatic protons of both the bases and the drug by complexation with ditercalinium suggest an interaction by intercalation of the two rings. However, the loss of twofold symmetry upon binding, as well as the chemical shift variation of the drug proton signals of one of the chromophores with temperature and concentration, favor a model in which the drug-nucleotide complexes have one ring of the drug intercalated and the other stacked on top of the external base pair. The intermolecular contacts between drug protons and nucleotide protons give a defined geometry for complexation that is consistent with the proposed model. In contrast, with d(CCTATAGG)2 several drug-nucleotide complexes were formed and a large increase in line broadening was observed at high drug-to-DNA ratios, precluding a detailed analysis of these complexes. However, the large upfield shift in the imino proton resonances together with the shielding of the ditercalinium ring protons favor a model with bis-intercalation of ditercalinium. This model is supported by the downfield shift of at least 4 out of 14 phosphorus signals. The results are compared with those obtained on ditercalinium binding to the homologous sequences d(CGCG)2 and d(TTCGCGAA)2, and discussed in terms of sequence specificity.

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of an oligodeoxyribonucleotide duplex refined via a hybrid relaxation matrix procedure

Assignment of the 1H and 31P resonances of a decamer DNA duplex, d(CGCTTAAGCG)2 was determined by two-dimensional COSY, NOESY and 1H-31P Pure Absorption phase Constant time (PAC) heteronuclear correlation spectroscopy. The solution structure of the decamer was calculated by an iterative hybrid relaxation matrix method combined with NOESY-distance restrained molecular dynamics. The distances from the 2D NOESY spectra were calculated from the relaxation rate matrix which were evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix-derived distances were then used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure was then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. JH3'-P coupling constants for each of the phosphates of the decamer were obtained from 1H-31P J-resolved selective proton flip 2D spectra. By using a modified Karplus relationship the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. Comparison of the 31P chemical shifts and JH3'-P coupling constants of this sequence has allowed a greater insight into the various factors responsible for 31P chemical shift variations in oligonucleotides. It also provides an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These correlations are consistent with the hypothesis that changes in local helical structure perturb the deoxyribose phosphate backbone. The variation of the 31P chemical shift, and the degree of this variation from one base step to the next is proposed as a potential probe of local helical conformation within the DNA double helix. The pattern of calculated epsilon and zeta torsional angles from the restrained molecular dynamics refinement agrees quite well with the measured JH3'-P coupling constants. Thus, the local helical parameters determine the length of the phosphodiester backbone which in turn constrains the phosphate in various allowed conformations.

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of a mismatched GA decamer oligodeoxyribonucleotide duplex

Assignment of the 1H and 31P NMR spectra of a tandem G.A mismatched base pair decamer oligodeoxyribonucleotide duplex, d(CCAAGATTGG)2, has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. Unusual downfield 31P resonances have been assigned by a pure absorption phase constant-time heteronuclear 31P-1H correlated spectrum to be associated with the phosphates on the 5'- and 3'-sides of the mismatched guanosine residue. JH3'-P coupling constants for each of the phosphates of the decamer were obtained from the 1H-31P J-resolved selective proton-flip 2D spectrum. The two most downfield-shifted 31P resonances each appear to consist of two overlapping signals that can be resolved into two distinct doublets with different coupling constants in the J-resolved spectrum. This as well as the temperature dependence of the 31P spectra demonstrates that two distinct conformations exist at lower temperatures. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. A linear correlation between 31P chemical shifts and the measured coupling constants is quite good (only when the larger set of coupling constants of the two most downfield 31P signals is included). The 31P chemical shifts as well as the measured coupling constants tend to follow the positional variation seen in other duplexes of interior phosphates resonating more upfield than terminal residues and of interior phosphates exhibiting smaller coupling constants; however, this pattern is disrupted at the site of the mismatch. Modeling and initial NOESY distance restrained molecular mechanics energy minimization and restrained molecular dynamics support previous observations that the mismatched guanine and adenine bases are both in anti conformations. Most significantly, the epsilon backbone torsional angle variaions calculated from the NOESY distance restrained structures are in agreement with both the crystal structure values and the measured JH3'-P coupling constants.

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

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

Two-dimensional 1H and 31P NMR spectra of a decamer oligodeoxyribonucleotide duplex and a quinoxaline ((MeCys3, MeCys7)(TANDEM) drug duplex complex

Assignment of the 1H and 31P NMR spectra of a decamer oligodeoxyribonucleotide duplex, d(CCCGATCGGG), and its quinoxaline ((MeCys3, MeCys7]TANDEM) drug duplex complex has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. The 31P chemical shifts of this 10 base pair oligonucleotide follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. While the 31P chemical shifts show sequence-specific variations, they also do not generally follow the Calladine "rules" previously demonstrated. 31P NMR also provides a convenient monitor of the phosphate ester backbone conformational changes upon binding of the drug to the duplex. Although the quinoxaline drug, [MeCys3, MeCys7]TANDEM, is generally expected to bind to duplex DNA by bis-intercalation, only small 31P chemical shift changes are observed upon binding the drug to duplex d(CCCGATCGGG). Additionally, only small perturbations in the 1H NMR and UV spectra are observed upon binding the drug to the decamer, although association of the drug stabilizes the duplex form relative to the other states. These results are consistent with a non-intercalative mode of association of the drug. Modeling and molecular mechanics energy minimization demonstrate that a novel structure in which the two quinoxaline rings of the drug binds in the minor groove of the duplex is possible.

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex

Assignment of the 1H and 31P NMR spectra of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. The downfield 31P resonance previously noted by Patel et al. (1982) has been assigned by both 17O labeling of the phosphate as well as a pure absorption phase constant-time heteronuclear 31P-1H correlated spectrum and has been associated with the phosphate on the 3' side of the extrahelical adenosine. JH3'-P coupling constants for each of the phosphates of the tridecamer were obtained from the 1H-31P J-resolved selective proton-flip 2D spectrum. By use of a modified Karplus relationship the C4-C3'-O3-P torsional angles (epsilon) were obtained. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. The 31P chemical shifts and epsilon torsional angles follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. Because the extrahelical adenosine significantly distorts the deoxyribose phosphate backbone conformation even several bases distant from the extrahelical adenosine, 31P chemical shifts show complex site- and sequence-specific variations. Modeling and NOESY distance-restrained energy minimization and restrained molecular dynamics suggest that the extrahelical adenosine stacks into the duplex. However, a minor conformation is also observed in the 1H NMR, which could be associated with a structure in which the extrahelical adenosine loops out into solution.

Sequence-dependent variations in the 31P NMR spectra and backbone torsional angles of wild-type and mutant Lac operator fragments

Assignment of the 31P resonances of a series of six sequenced-related tetradecamer DNA duplexes, d(TGTGAGCGCTCACA)2, d(TATGAGCGCTCATA)2, d(TCTGAGCGCTCAGA)2, d(TGTGTGCGCACACA)2, d(TGTGACGCGTCACA)2 and d(CACAGTATACTGTG)2, related to the lac operator DNA sequence was determined either by site-specific 17O labeling of the phosphoryl groups or by two-dimensional 1H-31P pure absorption phase constant time (PAC) heteronuclear correlation spectroscopy. J(H3'-P) coupling constants for each of the phosphates of the tetradecamers were obtained from 1H-31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. Comparison of the 31P chemical shifts and J(H3'-P) coupling constants of these sequences has allowed greater insight into those various factors responsible for 31P chemical shift variations in oligonucleotides and provided an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These sequence-specific variations in the conformation of the DNA sugar phosphate backbone of various lac operator DNA sequences can possibly explain the sequence-specific recognition of DNA by DNA binding proteins, as mediated through direct contacts between the phosphates and the protein.

N.m.r. and c.d. studies of the DNA fragments d(TATATATA) and d(TATATA) in solution

DNA fragments d(TATATATA) and d(TATATA) were studied in low-salt aqueous solutions and found to coexist in more than one conformer. 1H-n.m.r. demonstrates that single-stranded and double-stranded states are involved in the conformational coexistence. Circular dichroism spectroscopy indicates a global B-DNA stacking of bases in the fragments. 31P-n.m.r. resonances of the TpA and ApT phosphodiester bonds are substantially separated in the spectra of both d(TATATATA) and d(TATATA) duplexes to suggest an alternating architecture of their backbones. In fact, the oligonucleotide duplexes are much more alternating than the corresponding polynucleotide under the same solution conditions. The alternating character of the d(TATATATA) double helix is further enhanced in molar caesium fluoride solutions. The oligonucleotide isomerization into X-DNA is, however, accompanied by gel formation, which makes high resolution n.m.r. measurements impossible.

NMR studies of tris-intercalation: solution structure and interaction of d(CTTCGCGCGAAG) with an acridine trimer

Tris-intercalation of an acridine trimer into the self-complementary dodecanucleotide d(CTTCGCGCGAAG) has been studied, in solution, by means of 1H and 31P nuclear magnetic resonance. In a first step all the non-exchangeable protons (except H5', H5"), the imino protons and seven of the eleven phosphorus have been assigned. The dodecanucleotide is shown to adopt a double helical B-type structure. Most of the sugar puckers are in the O1'endo range, those of the internal guanosines being closer to C2'endo. Deviations from the canonical B structure are observed in the base stacking and the phosphodiester torsional angles at the 3T4C5G stretch. The addition of an acridine trimer to the base-paired dodecanucleotide leads to the conclusion that the trimer, which is in slow exchange at the NMR time scale, tris-intercalates into the three C(3'-5')G sites of the central core, according to the excluded site model. This is evidenced by the large (1.4 ppm) upfield shift experienced by the imino protons of the three internal guanines and the shielding undergone by the acridine ring protons. Tris-intercalation is also supported by the downfield shift experienced by 6 out of the 22 phosphorus. Two of them are shifted by nearly 2 ppm, a shift range reported for oligonucleotides complexed to actinomycin D; this suggests that the structure of the backbone of the dodecanucleotide is altered.

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.

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)

Oligodeoxyribonucleoside methylphosphonates. NMR and UV spectroscopic studies of Rp-Rp and Sp-Sp methylphosphonate (Me) modified duplexes of (d[GGAATTCC])2

1H NMR chemical shift assignments for the title compounds were made for most of the 1H signals using two-dimensional nuclear Overhauser effect (2D-NOE) data, which were also used to establish the absolute configuration at the modified phosphorus. The chemical shifts were similar to those reported [Broido, M.S., et al. (1985) Eur. J. Biochem. 150, 117-128] for the unmodified, parent, B-type duplex [d(GGAATTCC)]2. Differences in chemical shifts were mostly localized to the nucleotides on the 5'- and 3'-sides of the modified phosphorus. The Rp-Rp isomers exhibited UV-derived Tm values similar to that of the parent duplex. On the other hand, the Sp-Sp isomers generally exhibited lower Tm values which correlated with P-CH3--H3' (n-1 nucleotide) cross peak intensities and 31P spectral parameters. The combined data argue for increased steric interactions with the Sp-P-Me methyl group as the modification position is moved toward the center of the oligomer. All of the Tm results can be explained in terms of three factors which result from replacement of a phosphate by a methylphosphonate group: reduction of oligomer charge; electronic and other substituent effects; steric interactions.