DNA fragment conformations. 1H-NMR comparative studies of helix-coil transition and conformation of d(C-A-C-G-T-G) and d(G-T-G-C-A-C). Influence of helix formation on proton chemical shifts

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.

An oligodeoxyribonucleotide N3'--> P5' phosphoramidate duplex forms an A-type helix in solution

The solution conformations of the dinucleotide d(TT) and the modified duplex d(CGCGAATTCGCG)2 with N3'--> P5' phosphoramidate internucleoside linkages have been studied using circular dichroism (CD) and NMR spectroscopy. The CD spectra indicate that the duplex conformation is similar to that of isosequential phosphodiester RNA, a A-type helix, and is different from that of DNA, a B-type helix, NMR studies of model dimers d(TpT) and N3'--> P5' phosphoramidate d(TnpT) show that the sugar ring conformation changes from predominantly C2'-endo to C3'-endo when the 3'-phosphoester is replaced by a phosphoramidate group. Two-dimensional NMR (NOESY, DQF-COSY and TOCSY spectra) studies of the duplex provide additional details about the A-type duplex conformation of the oligonucleotide phosphoramidate and confirm that all furanose rings of 3'-aminonucleotides adopt predominantly N-type sugar puckering.

A 19F NMR study of the interaction of 3-fluoro-4-demethoxydaunomycin with the hexanucleotide d(TCCGGA)2

19F NMR spectroscopy has been used to study the binding of the fluorinated anthracycline 3-fluoro-4-demethoxydaunomycin with the hexanucleotide d(TCCGGA)2. In the spectrum of the 1:1 anthracycline-d(TCCGGA)2 complex four resonances of approximately equal intensity were observed. This indicated that 4-demethoxydaunomycin intercalated at all possible sites with similar affinity. This suggests that the specific high affinity binding sites that are observed in anthracycline-DNA footprinting experiments are strongly regulated by the local DNA conformation.

Transcriptional enhancer related DNA sequences: anomalous 1H NMR NOE crosspeaks

A dynamic heterogeneity which correlates with the function of the operator DNA in the lactose operon of E. coli. was previously observed (1) as a local minimum in the thymine imino proton T1 centered at a GTG/C-CAC sequence. Since this triplet occurs frequently in DNA regulatory regions, it was proposed that these sequences may be part of a structural element for specific protein interaction. We examine here three additional biologically significant 17 base pair duplexes containing GTG/CAC triplets: (1) a sequence from the mouse heavy chain immunoglobulin enhancer, (2) a sequence from the critical core of the Simian Virus 40 (SV40) enhancer, and (3) a sequence from pBR322 plasmid used as control for experiments with the SV40 DNA sequences. The 1H NMR resonance assignment for nearly all the nonexchangeable protons for both eukaryotic enhancer duplexes with the exception of the H5'/H5" protons was accomplished to use for structural analysis of these duplexes. The data presented show several NOE's associated with the GTG/CAC triplets which suggest structural variation from uniform B-DNA. In addition, anomalous broad crosspeaks for the fixed thymine methyl to its own H6 proton in combination with the imino proton kinetics associated with these triplets reinforces the original observation of a sequence dependent dynamic variation.

A 19F-NMR study of 3-fluoro-4-demethoxydaunomycin intercalation complexes with the hexanucleotide d(CTGCAG)2

Equilibrium systems containing intercalation complexes formed between the novel anthracycline drug, 3-fluoro-4-demethoxydaunomycin (3FD), and the hexanucleotide duplex d(CTGCAG)2 have been studied by 19F-NMR spectroscopy. Solutions containing a 1:1 molar ratio of 3FD/d(CTGCAG)2 gave four 19F signals which have been assigned to each of four possible intercalation isomers for the 1:1 3FD.d(CTGCAG)2 complex, which we denote by [d(CTGCAG)2][3FD]; these were where 3FD bound between the 5'-CT-3', 5'-TG-3', 5'-GC-3' or 5'-CA-3' base sequences, with the drug sugar moiety lying in the minor groove and pointed in the 3' direction in each case. Changes in temperature and NaCl concentration affecting the equilibrium distribution of these isomers were studied and indicated that no overriding binding site preference prevailed under standard biochemical conditions. Formation of some of the 2:1 3FD.d(CTGCAG)2 complex occurred when a solution of [d(CTGCAG)2][3FD] was exposed to excess 3FD; however, this complex was unstable to gel filtration and no co-operative binding of the second 3FD molecule was observed.

Conformations in crystals and solutions of d(CACGTG), d(CCGCGG) and d(GGCGCC) studied by vibrational spectroscopy

Crystals of self complementary DNA hexamers d(CACGTG), d(CCGCGG) and d(GGCGCC) were grown by vapour diffusion technique and studied by microRaman and microIR spectroscopies. The oligonucleotides were studied in parallel in solution by vibrational spectroscopy. A B- greater than Z transition was detected by Raman spectroscopy during the crystallization process for d(CACGTG). Vibrational spectroscopy shows that the d(GGCGCC) crystals adopt a B geometry. On the contrary the d(CCGCGG) sequence which is shown to be able to undergo in solution or in films quite easily the B- greater than Z transition, remains trapped in crystals in a geometry which may correspond to an intermediate conformation often proposed in models of the B- greater than Z transition. The crystals used in this study were characterized by X-ray diffraction. The unit cell and space group have been determined.

Proton and phosphorus NMR study of RNA self-complementary hexamers: influence of the 5-methylcytidine on the conformational transitions and the molecular motions

We present proton and phosphorus NMR data, which contribute to explain why the 5-methylcytidine (m5C) differently affects the conformational transitions of CGUAm5CG (Bloch et al., FEBS Letters 219, 464 (1987)) and CGm5CGCG (Ceolin et al., J. Am. Chem. Soc. 109, 2539 (1987)). The atypical intermediate form observed in the random coil to A-like duplex transition of CGm5CGCG is determined: this original duplex structure exhibits a frame-shift pairing of the two strands, with a partial conservation of the A-type structure. In the case of CGUAm5CG such a pairing process would lead to a complete mismatch pairing. This feature probably explains that the NMR data of CGUACG and CGUAm5CG are similar and consistent with a random coil to A-type helix transition. Nevertheless a significant difference between the correlation times observed for the molecular motions of the two duplexes is detected. Another example of unusual conformational transition of methylated RNA oligomers is given by CAm5CGUG.

Structure of d(CACGTG), a Z-DNA hexamer containing AT base pairs

The left-handed Z-DNA conformation has been observed in crystals made from the self-complementary DNA hexamer d(CACGTG). This is the first time that a non disordered Z form is found in the crystal structure of an alternating sequence containing AT base pairs without methylated or brominated cytosines. The structure has been determined and refined to an agreement factor R = 22.9% using 746 reflections in the resolution in the resolution shell 7 to 2.5 A. The overall shape of the molecule is very similar to the Z-structure of the related hexamer d(CG)3 confirming the rigidity of the Z form. No solvent molecules were detected in the minor groove of the helix near the A bases. The disruption of the spine of hydration in the AT step appears to be a general fact in the Z form in contrast with the B form. The biological relevance of the structure in relation to the CA genome repeats is discussed.

Resonance assignments of non-exchangeable protons in B type DNA oligomers, an overview

The chemical shifts of 1H resonances of non exchangeable protons (except H5', H5" and adenine H2) of over six hundred nucleotides have been collected. The influence which the base of the nucleotide itself as well as the bases on its 5' and 3' side exert on the chemical shifts of the various resonances has been investigated. Most of the resonances appear to be predominantly influenced by only one base. For H2', H2", H3', H4' and H6/H8 this is the base of the central nucleotide, for H5(C) and CH3(T) it is the one on the 5' side and for H1' it is the one on the 3' side. Chemical shift distribution profiles are presented which allow an estimation of the probability of finding a particular resonance at a particular position in the spectrum.

An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

By means of one- and two-dimensional NMR spectroscopy the solution structures of the partly self-complementary octamer d(m5C-G-m5C-G-A-G-m5C-G) were investigated. It is shown that this DNA fragment, under conditions of high DNA concentration (8 mM DNA) and/or high ionic strength prefers to adopt a duplex structure. At low DNA concentration (0.4 mM DNA), the duplex exists in a 1:1 slow equilibrium with a monomeric hairpin form. Addition of salt destabilizes the hairpin structure in favour of the dimer. At high temperatures the hairpin form, as well as the dimer structure, exist in a fast equilibrium with the random-coil form. For the hairpin/random-coil equilibrium a Tm of 329 K and a delta H degree of -121 kJ.mol-1 were deduced. These thermodynamic parameters are independent of the DNA concentration, as is expected for a monomeric structure. For the dimer to coil transition a Tm of 359 K (1 M DNA) and a delta H degree of -285 kJ.mol duplex-1 were derived. The thermodynamic data of the hairpin-coil transition mutually agree with those recently reported for the hairpin to random coil equilibrium of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) [Orbons, L. P. M., van der Marel, G. A., van Boom, J. H. & Altona, C. (1987) J. Biomol. Struct. Dyns. 4, 939-963]. It is demonstrated that the dimer structure exhibits B-DNA characteristics, as is witnessed by the NOESY experiments and the analysis of the proton-proton coupling data. It is shown that the base-pair formation of the G x A mismatches is anti-anti. A comparison of 1H and 31P chemical-shift data of the title compound with those of a well-characterized B-DNA structure reveals large differences in the dm5C(3)-dG(4)-dA(5) part of the mismatched dimer structure. These differences apparently indicate some major local structural changes due to the incorporation of the G x A mismatches. Under the most extreme conditions used (i.e. up to 3 M NaCl or 75% CH3OH in the presence of 10 mM MgCl2) no Z-DNA structure was observed. It is shown that the structural features of the hairpin form of the title compound mimic those of the hairpin structure of d(m5C-G-m5C-G-T-G-m5C-G). An energy-minimized model of the hairpin form is given.

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)

Z helix-coil transition of d(C-Br8G-C-G-C-Br8G) studied by CD, 1H-NMR and IR spectroscopies

The conformation of d(C-Br8G-C-G-C-Br8G) in aqueous solution was studied by CD and 1H-NMR spectroscopy and in condensed phase by IR spectroscopy. Whether in 0.1 M or 3 M NaCl solution or in film the only double helical structure adopted by brominated d(C-G)3 oligomer is the Z form. The IR spectrum of the film presents all the characteristic absorptions of the Z conformation and in particular is indicative of a syn conformation for the central guanosine as well as for the brominated one. Imino proton resonances of d(C-Br8G-C-G-C-Br8G) demonstrating the duplex formation were observed up to 60 degrees C. It is interesting to note that the significant highfield shifts of the dC H5" exocyclic sugar protons characteristic of the non exchangeable proton spectra of d(C-G)3 containing 5-methyl dC residues in the Z form were also detected in the proton spectrum of brominated oligomer. Whereas formation of the Z helix of methylated d(C-G)3 oligomers dependent on the salt concentration was found to occur via the preliminary formation of a B helix even in 4 M NaCl solution, the Z helix of d(C-Br8G-C-G-C-Br8G) is obtained directly from the coil form. However, IR data suggest that in the Z form of d(C-Br8G-C-G-C-Br8G), the overlapping of the base planes should be slightly different in comparison with the stacking observed in d(C-G)3 crystals. The kinetic data (activation energy and lifetime) of the Z helix-coil transition of brominated d(C-G)3 are compared to those of the B helix-coil transition observed for methylated d(C-G)3 in 0.1 M NaCl solution while the thermodynamic data of these two reactions (enthalpy and midpoint temperature) are slightly different.

Influence of uracil defect on DNA structure: 1H NMR investigation at 500 MHz

The local structure of two self complementary oligonucleotides d(GTAC-GTAC) and d(GTACGUAC) which differ only by the presence of uracil, not a normal component of DNA, have been investigated by 1H NMR at 500 MHz. The two octamers exhibit the same thermodynamical constants (t 1/2, delta H), their exchangeable protons broaden and disappear at the same temperature. The T-U substitution did not induce any significant changes on non exchangeable protons resonances from 2-D COSY and 2-D NOESY experiments. So the two octamers exhibit the same global structure. The only variation was detected by 1D NOE measurements: the base orientations around the N glycosidic bonds (chi angles) are different.

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.

Influence of dA X dT and d(2aminoA) X dT base pairs on the B in equilibrium Z transition of DNA fragments. 1H-NMR study of d(C-G-C-A-m5C-G-T-G-m5C-G), d(m5C-G-C-A-m5C-G-T-G-C-G) and d(C-2aminoA-C-G-T-G)

The Helical structures of d(C-G-C-A-m5C-G-T-G-m5C-G), d(m5C-G-C-A-m5C-G-T-G-C-G) and d(C-2aminoA-C-G-T-G) were studied in aqueous solution at various salt concentrations and temperatures by 1H-NMR spectroscopy. In 0.1 M NaCl solution only the B form was evidenced for these DNA fragments whereas in 4 M NaCl both B and Z forms, in slow exchange on the NMR time scale, were observed. Under these conditions the Z form accounted for less than 60% of the decamer conformation; conversely d(C-G)3 hexamers containing methylated cytidines were predominantly in the Z form (greater than 90%) [Tran-Dinh et al. (1984) Biochemistry 23, 1362; Cavaillès et al. (1984) J. Biomol. Struct. Dyn. 1, 1347-1371]. On the other hand, d(C-2aminoA-C-G-T-G) in which the d(2aminoA) X dT base pair forms three hydrogen bonds, was found to adopt the Z conformation in 4M NaCl solution which was not the case for d(C-A-C-G-T-G) (unpublished results). The present study shows that the B in equilibrium Z transition in solution is highly sequence-dependent and that correlation exists between the stability of the duplexes (essentially governed by the number of hydrogen bonds between complementary bases) and their ability to adopt the Z conformation.

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.

The B----Z transition in two synthetic oligonucleotides: d(C-2-amino-ACGTG) and d(m5CGCAm5CGTGCG) studied by IR, NMR and CD spectroscopies

The sequences CA'CGTG (where A' = 2-aminodeoxyadenosine) and m5CGCAm5CGTGCG are prepared and studied by IR, CD and 1H-NMR. Infrared spectra demonstrate the capacity of the modified hexamer and decamer to adopt a Z conformation. The influence of the NH2 substitution on the adenine or of the methylated terminal part of the decamer acting with the increase of the DNA concentration stabilizes the Z conformation at room temperature in low humidity films. Very weak proportion of Z conformation is detected in UV dilute solutions. In more concentrated NMR solutions, the Z proportion induced by high salt content is only 20-25%. The effects of the concentration and of the covalent modification of the bases are discussed.

DNA fragment conformations. A 1H-NMR conformational analysis of the d(G-G)-chelated platinum-oligonucleotide d(A-T-G-G)cisPt

The conformation of d(A-T-G-G) and d(A-T-G-G)cisPt has been investigated by 1H-NMR at 500 MHz and 90 MHz under various experimental conditions of temperature and concentration. Analysis of the coupling constants between the deoxyribose protons shows that all the sugar rings of d(A-T-G-G) adopt the S(C2'-endo) conformation most of the time. By contrast, in the platinated tetramer, d(A-T-G-G)cisPt, the N(C3'-endo) conformation is highly predominant for the internal dG residue while the S(C2'-endo) conformation is largely favoured for the other residues as in the case of the unplatinated compound. The relaxation time and nuclear Overhauser effect measurements indicate that the orientation of the two guanines of d(A-T-G-G)cisPt is anti in agreement with the previous results obtained for the dimers: r(G-G)cisPt, d(G-G)cisPt. On lowering the temperature from 80 degrees C to 20 degrees C, several proton resonances of d(A-T-G-G)cisPt exhibit large chemical shift and linewidth variations. The most spectacular temperature effect was observed for the internal dG(H1') and dT(H4') protons. All the delta = f(t) curves display a sigmoid form with the same mid-point temperature of 44 +/- 2 degrees C. This mid-point temperature together with the observed chemical shift and linewidth variations were found to be independent of the d(A-T-G-G)cisPt concentration. These results suggest that d(A-T-G-G)cisPt can adopt two different conformations depending on the temperature. The enthalpy for the transition between the high and low temperature conformations is about 84 kJ/mol.

Correlation of lac operator DNA imino proton exchange kinetics with its function

The kinetics for imino hydrogen exchange, at individual base pairs in the DNA sequence corresponding to the lactose operon operator of Escherichia coli, has been examined by NMR saturation recovery measurements as a function of temperature. Three 17-base-pair subsections of the lac operator DNA were chemically synthesized for these studies. The results support our previous observations in the 36-base-pair complete lac operator DNA fragment that has been used in our previous NMR studies. The results indicate faster opening kinetics at a GTG/CAC that is also the site of operator mutations leading to the highest level of constitutive beta-galactosidase synthesis. The GTG/CAC sequence occurs frequently and often symmetrically in prokaryotic and eukaryotic DNA sites where one anticipates specific protein interaction for gene regulation or recombination.

1H NMR and circular dichroism studies of the B and Z conformations of the self-complementary deoxyhexanucleotide d(m5C-G-C-G-m5-C-G): mechanism of the Z-B-coil transitions

The double-helical conformations of d(m5-C-G-C-G-m5-C-G) in aqueous solution were studied by circular dichroism and 1H NMR spectroscopy. In 0.1 M NaCl, only the B form is detected whereas the Z form is strongly predominant in 3 M NaCl. In the presence of 2 M NaCl, two resonance signals corresponding to the B and Z duplexes were observed for each proton below 50 degrees C, indicating a slow exchange between B and Z. However, the B-Z exchange becomes intermediate or fast in the 55-80 degrees C temperature interval. By contrast the exchange between B helix and single-stranded (or coil) forms is much faster for the same temperature conditions. The Z form is only detectable when the coil form is practically absent. With decreasing temperature the B form decreases in favor of the Z form. From proton line-width measurements under various experimental conditions, it was also shown that Z exchanges only with B, while the latter also exchanges with the single-stranded form (S): Z in equilibrium B in equilibrium S. The enthalpy value is about 8 +/- 1 kcal/mol for the B-Z transition and about 40 +/- 2 kcal/mol for the B-S dissociation (2 M NaCl solution). The activation energy is about 47 +/- 2 kcal/mol for the Z----B and 39 +/- 2 kcal/mol for the B----Z reaction. Very good agreement between the experimental results and computed data (based on the above kinetic reaction model) was found for the B, Z, and coil proportions. The B-Z transition of methylated d(C-G)n oligomers is only possible when the Watson-Crick hydrogen bonds between the CG base pairs are firmly maintained; otherwise, the transformation from B to Z would not occur, and B-S dissociation would take place instead.