Sequence dependence of hydrogen exchange kinetics in DNA duplexes at the individual base pair level in solution

A novel transition pathway of ligand-induced topological conversion from hybrid forms to parallel forms of human telomeric G-quadruplexes

The folding topology of DNA G-quadruplexes (G4s) depends not only on their nucleotide sequences but also on environmental factors and/or ligand binding. Here, a G4 ligand, 3,6-bis(1-methyl-4-vinylpyridium iodide)-9-(1-(1-methyl-piperidinium iodide)-3,6,9-trioxaundecane) carbazole (BMVC-8C3O), can induce topological conversion of non-parallel to parallel forms in human telomeric DNA G4s. Nuclear magnetic resonance (NMR) spectroscopy with hydrogen-deuterium exchange (HDX) reveals the presence of persistent imino proton signals corresponding to the central G-quartet during topological conversion of Tel23 and Tel25 G4s from hybrid to parallel forms, implying that the transition pathway mainly involves local rearrangements. In contrast, rapid HDX was observed during the transition of 22-CTA G4 from an anti-parallel form to a parallel form, resulting in complete disappearance of all the imino proton signals, suggesting the involvement of substantial unfolding events associated with the topological transition. Site-specific imino proton NMR assignments of Tel23 G4 enable determination of the interconversion rates of individual guanine bases and detection of the presence of intermediate states. Since the rate of ligand binding is much higher than the rate of ligand-induced topological conversion, a three-state kinetic model was evoked to establish the associated energy diagram for the topological conversion of Tel23 G4 induced by BMVC-8C3O.

Structural basis of sodium-potassium exchange of a human telomeric DNA quadruplex without topological conversion

Understanding the mechanism of Na(+)/K(+)-dependent spectral conversion of human telomeric G-quadruplex (G4) sequences has been limited not only because of the structural polymorphism but also the lack of sufficient structural information at different stages along the conversion process for one given oligonucleotide. In this work, we have determined the topology of the Na(+) form of Tel23 G4, which is the same hybrid form as the K(+) form of Tel23 G4 despite the distinct spectral patterns in their respective nuclear magnetic resonance (NMR) and circular dichroism spectra. The spectral difference, particularly the well-resolved imino proton NMR signals, allows us to monitor the structural conversion from Na(+) form to K(+) form during Na(+)/K(+) exchange. Time-resolved NMR experiments of hydrogen-deuterium exchange and hybridization clearly exclude involvement of the global unfolding for the fast Na(+)/K(+) spectral conversion. In addition, the K(+) titration monitored by NMR reveals that the Na(+)/K(+) exchange in Tel23 G4 is a two-step process. The addition of K(+) significantly stabilizes the unfolding kinetics of Tel23 G4. These results offer a possible explanation of rapid spectral conversion of Na(+)/K(+) exchange and insight into the mechanism of Na(+)/K(+) structural conversion in human telomeric G4s.

Molecular electrostatic potentials of DNA base-base pairing and mispairing

An understanding of why adenine (A) pairs with thymine (T) and cytosine (C) with guanine (G) in DNA is very useful in the design of sensors and other related devices. We report the use of dissociation energies, geometries and molecular electrostatic potentials (MEPs) to justify the canonical (AT and CG) Watson-Crick pairs. We also analyze all mismatches in both configurations-cis and trans-with respect to their glycoside bonds. As expected, we found that the most stable pair configuration corresponds to CG, providing an energy criterion for that preferred configuration. The reason why A gets together with T is much more difficult to explain as the energy of this pair is smaller than the energy of some other mismatched pairs. We tested MEPs to see if they could shed light on this problem. Interestingly, MEPs yield a unique pattern (shape) for the two canonical cases but different shapes for the mismatches. A tunnel of positive potential surrounded by a negative one is found interconnecting the three H-bonds of CG and the two of AT. This MEP tunnel, assisted partially by energetics and geometrical criteria, unambiguously determine a distinctive feature of the affinity between A and T as well as that between G and C.

Comparison of the thermodynamics and base-pair dynamics of a full LNA:DNA duplex and of the isosequential DNA:DNA duplex

Locked nucleic acids (LNA), conformationally restricted nucleotide analogues, are known to enhance pairing stability and selectivity toward complementary strands. With the aim to contribute to a better understanding of the origin of these effects, the structure, thermal stability, hybridization thermodynamics, and base-pair dynamics of a full-LNA:DNA heteroduplex and of its isosequential DNA:DNA homoduplex were monitored and compared. CD measurements highlight differences in the duplex structures: the homoduplex and heteroduplex present B-type and A-type helical conformations, respectively. The pairing of the hybrid duplex is characterized, at all temperatures monitored (between 15 and 37 degrees C), by a larger stability constant but a less favorable enthalpic term. A major contribution to this thermodynamic profile emanates from the presence of a hairpin structure in the LNA single strand which contributes favorably to the entropy of interaction but leads to an enthalpy penalty upon duplex formation. The base-pair opening dynamics of both systems was monitored by NMR spectroscopy via imino protons exchange measurements. The measurements highlight that hybrid G-C base-pairs present a longer base-pair lifetime and higher stability than natural G-C base-pairs, but that an LNA substitution in an A-T base-pair does not have a favorable effect on the stability. The thermodynamic and dynamic data confirm a more favorable stacking of the bases in the hybrid duplex. This study emphasizes the complementarities between dynamic and thermodynamical studies for the elucidation of the relevant factors in binding events.

Profiling the thermodynamic softness of adenoviral promoters

We showed previously that anharmonic DNA dynamical features correlate with transcriptional activity in selected viral promoters, and hypothesized that areas of DNA softness may represent loci of functional significance. The nine known promoters from human adenovirus type 5 were analyzed for inherent DNA softness using the Peyrard-Bishop-Dauxois model and a statistical mechanics approach, using a transfer integral operator. We found a loosely defined pattern of softness peaks distributed both upstream and downstream of the transcriptional start sites, and that early transcriptional regions tended to be softer than late promoter regions. When reported transcription factor binding sites were superimposed on our calculated softness profiles, we observed a close correspondence in many cases, which suggests that DNA duplex breathing dynamics may play a role in protein recognition of specific nucleotide sequences and protein-DNA binding. These results suggest that genetic information is stored not only in explicit codon sequences, but also may be encoded into local dynamic and structural features, and that it may be possible to access this obscured information using DNA dynamics calculations.

NMR analysis of duplex d(CGCGATCGCG)2 modified by Λ- and Δ-[Ru(bpy)2(m-GHK)]Cl2 and DNA photocleavage study

The interaction of the diastereomeric complexes Lambda-[Ru(bpy)2(m-GHK)]Cl2 and Delta-[Ru(bpy)2(m-GHK)]Cl2 (bpy is 2,2'-bipyridine, GHK is glycine-L-histidine-L-lysine) with the deoxynucleotide duplex d(5'-CGCGATCGCG)2 was studied by means of 1H NMR spectroscopy. At a Delta-isomer to DNA ratio of 1:1, significant shifts for the metal complex are observed, whereas there is negligible effect on the oligonucleotide protons and only one intermolecular nuclear Overhauser effect (NOE) is present at the 2D nuclear Overhauser enhancement spectroscopy spectrum. The 1Eta NMR spectrum at ratio 2:1 is characterized by a slight shift for the Delta-isomer's bpy aromatic protons as well as significant shifts for the decanucleotide G4 H1' and Eta2'', A5 H2, G10 H1', T6 NH and G2 NH protons. Furthermore, at ratio 2:1, 11 intermolecular NOEs are observed. The majority of the NOEs involve the sugar Eta2' and Eta2'' protons sited in the major groove of the decanucleotide. Increasing the Delta-isomer to d(CGCGATCGCG)2 ratio to 5:1 results in noteworthy spectral changes. The Delta-isomer's proton shifts are reduced, whereas significant shifts are observed for the decanucleotide protons, especially the sugar protons, as well as for the exchangeable protons. Interaction is characterized by the presence of only one intermolecular NOE. Furthermore, there is significant broadening of the imino proton signals as the ratio of the Delta-isomer to DNuAlpha increases, which is attributed to the opening of the two strands of the duplex. The Lambda-isomer, on the other hand, approaches the minor groove of the oligonucleotide and interacts only weakly, possibly by electrostatic interactions. Photocleavage studies were also conducted with the plasmid pUC19 and a 158-bp restriction fragment, showing that both diastereomers cleave DNA with similar efficiency, attacking mainly the guanines of the sequence probably by generating active oxygen species.

Trans-hydrogen bond deuterium isotope effects of A:T base pairs in DNA

The chemical shifts of (13)C2 of adenosine residues of DNA were observed to experience a through-space or trans-hydrogen bond isotope effect as a result of deuterium substitution at the imino hydrogen site of base-paired thymidine residues. NMR measurements of several self-complementary DNA duplexes at natural abundance (13)C in 50% H(2)O, 50% D(2)O solvent mixtures yielded an average trans-hydrogen bond isotope effect, (2h)Delta(13)C2, of -47 ppb. The data suggest that stronger hydrogen bonds have more negative (2h)Delta(13)C2 values, which means that A:T N1.H3 hydrogen bonds increase the anharmonicity of the effective vibrational potential of H3. However, (2h)Delta(13)C2 values do not correlate with intra-residue (2)Delta(13)C4 values of thymidine observed here and earlier (Vakonakis et al., 2003), which suggests that (2h)Delta(13)C2 is not determined entirely by hydrogen bond strength. Instead, the variations observed in (2h)Delta(13)C2 values suggest that they may also be sensitive to base pair geometry.

DNA dynamically directs its own transcription initiation

It has long been known that double-stranded DNA is subject to temporary, localized openings of its two strands. Particular regions along a DNA polymer are destabilized structurally by available thermal energy in the system. The localized sequence of DNA determines the physical properties of a stretch of DNA, and that in turn determines the opening profile of that DNA fragment. We show that the Peyrard-Bishop nonlinear dynamical model of DNA, which has been used to simulate denaturation of short DNA fragments, gives an accurate representation of the instability profile of a defined sequence of DNA, as verified using S1 nuclease cleavage assays. By comparing results for a non-promoter DNA fragment, the adenovirus major late promoter, the adeno-associated viral P5 promoter and a known P5 mutant promoter that is inactive for transcription, we show that the predicted openings correlate almost exactly with the promoter transcriptional start sites and major regulatory sites. Physicists have speculated that localized melting of DNA might play a role in gene transcription and other processes. Our data link sequence-dependent opening behavior in DNA to transcriptional activity for the first time.

Deuterium isotope effects and fractionation factors of hydrogen-bonded A:T base pairs of DNA

Deuterium isotope effects and fractionation factors of N1.H3-N3 hydrogen bonded Watson-Crick A:T base pairs of two DNA dodecamers are presented here. Specifically, two-bond deuterium isotope effects on the chemical shifts of (13)C2 and (13)C4, (2)delta(13)C2 and (2)delta(13)C4, and equilibrium deuterium/protium fractionation factors of H3, Phi, were measured and seen to correlate with the chemical shift of the corresponding imino proton, delta(H3). Downfield-shifted imino protons associated with larger values of (2)delta(13)C2 and (2)delta(13)C4 and smaller Phi values, which together suggested that the effective H3-N3 vibrational potentials were more anharmonic in the stronger hydrogen bonds of these DNA molecules. We anticipate that (2)delta(13)C2, (2)delta(13)C4 and Phi values can be useful gauges of hydrogen bond strength of A:T base pairs.

Correlation between normal modes in the 20-200 cm-1 frequency range and localized torsion motions related to certain collective motions in proteins

In certain biologically relevant collective motions, such as protein domain motions and sub-domain motions, large amplitude movements are localized in one or a few flexible regions consisting of a small number of residues. This paper explores the possible use of normal mode analysis in probing localized vibrational torsion motions in these flexible regions that may be related to certain collective motions. The normal modes of 10 structures of five proteins in different conformation (TRP repressor, calmodulin, calbindin D(9k), HIV-1 protease and troponin C), known to have shear or hinge domain or sub-domain motion, respectively, are analyzed. Our study identifies, for each structure, unique normal modes in the 20-200 cm-1 frequency range, whose corresponding motions are primarily concentrated in the region where large amplitude torsion movements of a known domain or sub-domain motion occur. This suggests possible correlation between normal modes at 20-200 cm-1 frequency range and initial fluctuational motions leading to localized collective motions in proteins, and thus the potential application of normal mode analysis in facilitating the study of biologically important localized motions in biomolecules.

NMR investigations of duplex stability of phosphorothioate and phosphorodithioate DNA analogues modified in both strands

Duplex formation from the self-complementary 12mer d(CGCGAATTCGCG) (Dickerson dodecamer) in which all phosphodiester linkages were replaced by phosphorothioate or phosphorodithioate linkages was studied using variable-temperature 1H and 31P NMR spectroscopy. Melting temperatures of the dodecamer, measured spectrophotometrically, showed significant decrease upon sulfur substitution (Tm 49 degrees C for the phosphorothioate and 21 degrees C for the phosphorodithioate, compared with 68 degrees C for the unmodified oligomer, in 1 M salt). Hyperchromicity observed upon melting of the dithioate was surprisingly low. NOESY spectra of the monothioate showed a cross-peak pattern characteristic for a right-handed duplex. Imino proton resonances of the duplex, shown by the mono- and the dithioate, were similar to those of the parent compound. In spite of monophasic melting curves, temperature dependence of the imino proton resonances and phosphorus resonances of the phosphorodithioate indicated heterogeneity with respect to base-pairing, compatible with the presence of a hairpin loop. Relaxation times (T1) of the imino protons in the phosphorothioate, determined by the saturation recovery method, were considerably shorter than in the unmodified oligomer. Base-pair lifetimes in the unmodified Dickerson dodecamer, determined by catalyst-dependent changes in relaxation rates of imino protons, were in the range of 2-30 ms at 20 degrees C. Strongly reduced base-pair lifetimes were found in the phosphorothioate analogue.

Direct estimation of base-pair exchange kinetics in oligo-DNA by a combination of NOESY and ROESY experiments

A new method for the determination of the kinetics of exchange of the imino protons of DNA duplex is reported using a combination NOESY and ROESY experiments at short mixing times (< or = 20 ms). These results have been compared with the commonly used longitudinal relaxation approach through the T1 measurement. To calculate kex and pi ex by ROESY-NOESY experiment, the volume of the cross-peaks between imino protons and water in the NOESY and ROESY spectra have been measured separately from the magnetization term. This work shows that the present approach for the measurement of the kinetics of slow exchanging imino protons of DNA duplex is comparable to the saturation recovery experiment in which the exchange rate can be accelerated by the addition of a base catalyst. The present ROESY-NOESY approach has been found to be particularly useful and reasonably accurate for the measurement of exchange kinetics of both the fast- and slow-exchanging imino protons in DNA duplex both under non-physiological and physiological condition where the saturation recovery method can not be used.

100ps molecular dynamic simulation of d(TATCACC)2

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.

Interaction of berenil with the EcoRI dodecamer d(CGCGAATTCGCG)2 in solution studied by NMR

The conformation of the EcoRI dodecamer d(CGCGAATTCGCG)2 has been examined in solution by 1H and 31P NMR. Spin-spin coupling constants and nuclear Overhauser (NOE) enhancement spectroscopy show that all deoxyriboses lie in the south domain, with a small admixture of the north conformation (0-20%). The time dependence of the nuclear Overhauser enhancements also reveals a relatively uniform conformation at the glycosidic bonds (average angle, chi = -114 degrees). The average helical twist is 36.5 degrees (9.8 base pairs per turn). Tilt angles are small (in the range 0 to -10 degrees), and roll angles are poorly determined. Unlike single-crystal X-ray studies of the same sequence, there is no evidence for asymmetry in the structure. Both the NOE intensities and 31P relaxation data imply conformational anomalies at the C3-G4/C9-G10 and the A5-A6/T7-T8 steps. Berenil binds in 1:1 stoichiometry to the dodecamer with high affinity (Kd = 1 microM at 298 K) and causes substantial changes in chemical shifts of the sugar protons of nucleotides Ado 5-Cyt 9 and of the H2 resonances of the two Ado residues. No significant asymmetry appears to be induced in the DNA conformation on binding, and there is no evidence for intercalation, although the binding site is not centrosymmetric. NOEs are observed between the aromatic protons of berenil and the H1' of both Thy 7 and Thy 8, as well as to Ado 5 and Ado 6 H2. These results firmly establish that berenil binds via the minor groove and closely approaches the nucleotides Ado 6, Thy 7, and Thy 8. On the basis of quantitative NOE spectroscopy and measurements of spin-spin coupling constants, changes in the conformations of the nucleotides are found to be small. Using the observed NOEs between the ligand and the DNA together with the derived glycosidic torsion angles, we have built models that satisfy all of the available solution data. The berenil molecule binds at the 5'-AAT (identical to 5'-ATT on the complementary strand) site such that (i) favorable hydrogen bonds are formed between the charged amidinium groups and the N3 atoms of Ado 6 and Ado 18 and (ii) the ligand is closely isohelical with the floor of the minor groove.

Solution structure of the EcoRI DNA sequence: refinement of NMR-derived distance geometry structures by NOESY spectrum back-calculations

The solution structure of the self-complementary DNA duplex [d(CGCGAATTCGCG)]2, which contains the EcoRI restriction site sequence GAATTC at the center, has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Time-dependent nuclear Overhauser effect spectra were used to obtain the initial cross-relaxation rates between 155 pairs of protons. These initial cross-relaxation rates were converted into interproton distances and entered into a distance (bounds) matrix. A distance geometry algorithm (DSPACE) was used to create embedded starting structures and to refine these structures until they showed good agreement with the distance matrix; symmetry constraints were included in the refinement procedure, making the two strands in the refined distance geometry structures virtually identical and significantly improving the agreement with the distance matrix. The NOESY spectrum for one of these distance geometry structures was then calculated from the explicit coordinates by numerically integrating all the z-magnetization transfer pathways among neighboring protons within a specified radius. Distances in this distance geometry structure that did not agree with the experimental NOESY time course were then adjusted accordingly. This process was iterated until a good agreement between calculated and experimental NOESY spectra was reached. The final structure, which generates good agreement with the experimental NOESY spectrum, displays kinks at the C3-G4 base step and at the A6-T7 base step that appear to be similar to those reported for the EcoRI restriction site DNA bound to its endonuclease. The solution structure is not the same as the crystal structure of this DNA duplex.

N.m.r. assignments and temperature-dependent conformational transitions of a mutant trp operator-promoter in solution

A total of 145 protons in the mutant trp operator-promoter sequence CGTACTGATTAATCAGTACG were assigned by one-dimensional and two-dimensional n.m.r. methods. Except at the sites of mutation (underlined), the chemical shifts and other n.m.r. parameters are very similar to those observed in the symmetrized wild-type sequence [Lefèvre, Lane & Jardetzky (1987) Biochemistry 26, 5076-5090]. Spin-spin-relaxation rate constants of the resolved base protons and intra- and inter-nucleotide nuclear-Overhauser-enhancement intensities argue for a sequence-dependent structure similar to that of the wild-type, except at and close to the sites of the mutation. The overall tumbling time as a function of temperature was determined from cross-relaxation rate constants for the H-6-H-5 vectors of the four cytosine residues. The values are consistent with the oligonucleotide maintaining a double-helical conformation over the entire temperature range 5-45 degrees C, and that internal motions of the bases are of small amplitude on the subnanosecond time scale. The temperature-dependence of chemical shifts, spin-spin-relaxation rate constants and cross-relaxation rate constants show the occurrence of two conformational transitions localized to the TTAA sequence in the centre of the molecule. The thermodynamics of the transition at the lower temperature (tm = 16 degrees C) were analysed according to a two-state process. The mid-point temperature is about 6 degrees C higher than in the wild-type sequence. The conformational transition does not lead to rupture of the Watson-Crick hydrogen bonds, but probably involves changes in the propellor twists of T.A-9 and T.A-10. The second transition occurs at about 40 degrees C, but cannot be fully characterized. This conformational variability seems to be a property of the sequence TTAA, and may have functional significance in bacterial promoters.

Phosphate catalyzed exchange rates of NH-N protons of DNA

We have studied the catalysis of the exchange of the hydrogen-bonded NH-N protons of the short DNA helix (d-CCAAGCTTGG)2 by phosphate addition. The NH exchange rates were monitored by the line widths of the corresponding NH resonances in the 1H nmr spectra. The exchange catalyst phosphate is most effective on the exchange rate of the terminal CG1 base pairs. However, all internal base pairs are also moderately affected by phosphate which suggests an exchange mechanism governed by a fast equilibrium between opened and closed states of the duplex. Within the limits of error the same effectiveness of phosphate on the exchange rate of all internal NH-N protons has been observed. With the exception of the terminal base pairs, no sequence and/or position specificity of the exchange rates of the NH-N protons of the base pairs has been found.

Sequence-specific assignments and their use in NMR studies of DNA structure

There has been a surge of recent interest, reflected by a sharp increase in the number of publications, in the area of high-resolution nuclear magnetic resonance (NMR) studies of DNA. The goal of many of these studies is to monitor the structure of biologically important DNA sequences directly in solution; the impetus for such studies was the realization, from early single-crystal X-ray structures, that nearest-neighbor context effects are a major determinant of local structure in short double-helical DNAs (Dickerson & Drew, 1981; Dickerson, 1983). Thus, instead of the previously assumed regular averaged structure of the double helix derived from fibre diffraction analysis, the more interesting concept emerged that specific sequence-dependent distortions from ‘classical’ DNA structure might be responsible for the recognition of such sequences by a variety of ligands such as repressors, polymerases, drugs, etc.

DNA and RNA: NMR studies of conformations and dynamics in solution

The early NMR research on nucleic acids was of a qualitative nature and was restricted to partial characterization of short oligonucleotides in aqueous solution. Major advances in magnet design, spectrometer electronics, pulse techniques, data analysis and computational capabilities coupled with the availability of pure and abundant supply of long oligonucleotides have extended these studies towards the determination of the 3-D structure of nucleic acids in solution.

Dissociation kinetics of 19 base paired oligonucleotide-DNA duplexes containing different single mismatched base pairs

The dissociation kinetics of 19 base paired oligonucleotide-DNA duplex containing a various single mismatched base pair are studied on dried agarose gels. The kinetics of the dissociation are first order under our experimental conditions. The incorporation of a single mismatched base pair destabilizes the DNA duplexes to some extent, the amount depending on the nature of the mismatched base pair. G-T and G-A mismatches slightly destabilize a duplex, while A-A, T-T, C-T and C-A mismatches significantly destabilize it. The activation energy for the overall dissociation processes for these oligonucleotide-DNA duplexes containing 19 base pairs is 52 +/- 2 Kcal mol-1 as determined from the slope of Arrhenius plot.

Influence of the polyamines spermine and spermidine on yeast tRNAPhe as revealed from its imino proton NMR spectrum

A comparison of imino proton NMR spectra of yeast tRNAPhe recorded at various solution conditions indicates, that polyamines have a limited effect on the structure of this tRNA molecule. Polyamines are found to catalyse the solvent exchange of several imino protons in yeast tRNAPhe not only of non hydrogen bonded imino protons, but also of imino protons of the GU and of some AU and tertiary base pairs. It is concluded that at low levels of catalysing components the exchange rates of the latter protons are not determined by the base pair lifetime. In the presence of high levels of spermidine the solvent exchange rates of imino protons of several base pairs in the molecule were assessed as a function of the temperature. Apparent activation energies derived from these rates were found to be less than 80 kJ/mol, which is indicative for (transient) independent opening of the corresponding base pairs. In the acceptor helix the GU base pair acts as a dynamic dislocation. The AU base pairs at one side of the GU base pair exhibit faster transient opening than the GC base pairs on the other side of this wobble pair. The base pairs m2GC10 and GC11 from the D stem and GC28 from the anticodon stem show relatively slow opening up to high temperatures. Model studies suggest that 1-methyladenosine, an element of tRNA itself, catalyses imino proton solvent exchange in a way similar to polyamines.

An NMR study of A-T base pair opening rates in oligonucleotides. Influence of sequence and of adenine methylation

We report relaxation time measurements by semi-selective and totally selective NMR techniques on the thymidine imino protons of d(GGATATCC) and d(GGm6ATATCC). For these oligonucleotides helix fraying, rather than single base pair opening, is the major exchange mechanism even 25 degrees C below the Tm. We have therefore applied a new saturation transfer technique to measure exchange rates at temperatures where fraying has a very small or negligible contribution. Measurements of exchange rates as a function of temperature give significantly different activation energies for base pairs 3 and 4 in d(GGATATCC). Adenine methylation results in a slowing down of the opening rate for the m6A-T base pair but surprisingly has an even greater effect upon the adjacent non-methylated A-T base pair.

Conformation, dynamics, and structural transitions of the TATA box region of self-complementary d[(C-G)n-T-A-T-A-(C-G)n] duplexes in solution

Structural and kinetic features of the TATA box located in the center of the alternating self-complementary d(C-G-C-G-T-A-T-A-C-G-C-G) duplex (TATA 12-mer) and d(C-G-C-G-C-G-T-A-T-A-C-G-C-G-C-G) duplex (TATA 16-mer) have been probed by high-resolution proton and phosphorus NMR spectroscopy in aqueous solution. The imino exchangeable Watson-Crick protons and the nonexchangeable base protons in the TATA box of the TATA 12-mer and TATA 16-mer duplexes have been assigned from intra and inter base pair nuclear Overhauser effect (NOE) measurements. Imino proton line-width and hydrogen exchange saturation recovery measurements demonstrate that the dA X dT base pairs in the TATA box located in the center of the TATA 12-mer and TATA 16-mer duplexes are kinetically more labile than flanking dG X dC base pairs. The proton and phosphorus NMR parameters of the TATA 12-mer monitor a cooperative premelting transition in the TATA box prior to the onset of the melting transition to unstacked strands. Phosphorus NMR studies have been unable to detect any indication of a right-handed B DNA to a left-handed Z DNA transition for the TATA 12-mer duplex in saturated NaCl solution. By contrast, we do detect the onset of the B to Z transition for the TATA 16-mer in saturated NaCl solution. Proton and phosphorus NMR studies demonstrate formation of a loop conformation with chain reversal at the TATA segment for the TATA 12-mer and TATA 16-mer duplexes on lowering the DNA and counterion concentration. The imino protons (10-11 ppm) and phosphorus resonances (3.5-4.0 ppm; 4.5-5.0 ppm) of the loop segment fall in spectral windows well resolved from the corresponding markers in fully paired segments so tha it should be possible to identify loops in longer DNA helixes. The equilibrium between the loop and fully paired duplex conformations of the TATA 12-mer and TATA 16-mer is shifted toward the latter on addition of moderate salt.

Conformation and dynamics of the Pribnow box region of the self-complementary d(C-G-A-T-T-A-T-A-A-T-C-G) duplex in solution

Nuclear magnetic resonance (NMR) has been used to monitor the conformation and dynamics of the d(C1-G2-A3-T4-T5-A6-T6-A5-A4-T3-C2-G1) self-complementary dodecanucleotide duplex (henceforth called Pribnow 12-mer), which contains a TATAAT Pribnow box and a central core of eight dA X dT base pairs. The exchangeable imino and nonexchangeable base protons have been assigned from one-dimensional intra and inter base pair nuclear Overhauser effect (NOE) measurements. Premelting conformational changes are observed at all the dA X dT base pairs in the central octanucleotide core in the Pribnow 12-mer duplex with the duplex to strand transition occurring at 55 degrees C in 0.1 M phosphate solution. The magnitude of the NOE measurements between minor groove H-2 protons of adjacent adenosines demonstrates that the base pairs are propeller twisted with the same handedness as observed in the crystalline state. The thymidine imino proton hydrogen exchange at the dA X dT base pairs has been measured from saturation recovery measurements as a function of temperature. The exchange rates and activation barriers show small variations among the four different dA X dT base pairs in the Pribnow 12-mer duplex.(ABSTRACT TRUNCATED AT 250 WORDS)

Deoxyguanosine-deoxyadenosine pairing in the d(C-G-A-G-A-A-T-T-C-G-C-G) duplex: conformation and dynamics at and adjacent to the dG X dA mismatch site

Nuclear magnetic resonance (NMR) has been used to monitor the conformation and dynamics of the d-(C1-G2-A3-G4-A5-A6-T6-T5-C4-G3-C2-G1) self-complementary dodecanucleotide (henceforth called 12-mer GA) that contains a dG X dA purine-purine mismatch at position 3 in the sequence. These results are compared with the corresponding d(C-G-C-G-A-A-T-T-C-G-C-G) dodecamer duplex (henceforth called 12-mer) containing standard Watson-Crick base pairs at position 3 [Patel, D.J., Kozlowski, S.A., Marky, L.A., Broka, C., Rice, J.A., Itakura, K., & Breslauer, K.J. (1982) Biochemistry 21, 428-436]. The dG X dA interaction at position 3 was monitored at the guanosine exchangeable H-1 and nonexchangeable H-8 protons and the nonexchangeable adenosine H-2 proton. We demonstrate base-pair formation between anti orientations of the guanosine and adenosine rings on the basis of nuclear Overhauser effects (NOE) observed between the H-2 proton of adenosine 3 and the imino protons of guanosine 3 (intra base pair) and guanosines 2 and 4 (inter base pair). The dG(anti) X dA(anti) pairing should result in hydrogen-bond formation between the guanosine imino H-1 and carbonyl O-6 groups and the adenosine N-1 and NH2-6 groups, respectively. The base pairing on either side of the dG X dA pair remains intact at low temperature, but these dG X dC pairs at positions 2 and 4 are kinetically destabilized in the 12-mer GA compared to the 12-mer duplex. We have estimated the hydrogen exchange kinetics at positions 4-6 from saturation-recovery measurements on the imino protons of the 12-mer GA duplex between 5 and 40 degrees C. The measured activation energies for imino proton exchange in the 12-mer GA are larger by a factor of approximately 2 compared to the corresponding values in the 12-mer duplex. This implies that hydrogen exchange in the 12-mer GA duplex results from a cooperative transition involving exchange of several base pairs as was previously reported for the 12-mer containing a G X T wobble pair at position 3 [Pardi, A., Morden, K.M., Patel, D.J., & Tinoco, I., Jr. (1982) Biochemistry 21, 6567-6574]. We have assigned the nonexchangeable base protons by intra and inter base pair NOE experiments and monitored these assigned markers through the 12-mer GA duplex to strand transition.(ABSTRACT TRUNCATED AT 400 WORDS)

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.

Possible molecular detent in the DNA structure at regulatory sequences

A common feature that appears in a number of DNA sites where proteins interact is the sequence GTG/CAC. In the lac operator this sequence leads to a region with a higher imino proton exchange rate well below the optical melting temperature. It is suggested that this reflects a structural feature recognized by proteins that bind specific sites on the DNA molecule.