NMR relaxation in DNA. I. The contribution of torsional deformation modes of the elastic filament

31P NMR investigation of backbone dynamics in DNA binding sites

The backbone conformation of DNA plays an important role in the indirect readout mechanisms for protein--DNA recognition events. Thus, investigating the backbone dynamics of each step in DNA binding sequences provides useful information necessary for the characterization of these interactions. Here, we use 31P dynamic NMR to characterize the backbone conformation and dynamics in the Dickerson dodecamer, a sequence containing the EcoRI binding site, and confirm solid-state 2H NMR results showing that the C3pG4 and C9pG10 steps experience unique dynamics and that these dynamics are quenched upon cytosine methylation. In addition, we show that cytosine methylation affects the conformation and dynamics of neighboring nucleotide steps, but this effect is localized to only near neighbors and base-pairing partners. Last, we have been able to characterize the percent BII in each backbone step and illustrate that the C3pG4 and C9pG10 favor the noncanonical BII conformation, even at low temperatures. Our results demonstrate that 31P dynamic NMR provides a robust and efficient method for characterizing the backbone dynamics in DNA. This allows simple, rapid determination of sequence-dependent dynamical information, providing a useful method for studying trends in protein-DNA recognition events.

Sequence-dependent dynamics in duplex DNA

The submicrosecond bending dynamics of duplex DNA were measured at a single site, using a site-specific electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of the mean squared amplitude of bending relative to the end-to-end vector defined by the weakly bending rod model. The bending dynamics monitored at the single site varied when the length and position of a repeated AT sequence, distant from the spin probe, were changed. As the distance between the probe and the AT sequence was increased, the mean squared amplitude of bending seen by the probe due to that sequence decreased. A model for the sequence-dependent internal flexural motion of duplex DNA, which casts the mean squared bending amplitudes in terms of sequence-dependent bending parameters, has been developed. The best fit of the data to the model occurs when the (AT)(n) basepairs are assumed to be 20% more flexible than the average of the basepairs within the control sequence. These findings provide a quantitative basis for interpreting the kinetics of biological processes that depend on duplex DNA flexibility, such as protein recognition and chromatin packaging.

Manifestations of slow site exchange processes in solution NMR: a continuous Gaussian exchange model

The effects of site exchange due to slow conformational changes in rapidly rotating molecules in solution are examined in detail. Significant gaps in the currently available theory are filled. The effects of site exchange on the lineshape, decay of a simple spin-echo, decay of the even echoes in a Carr-Purcell-Meiboom-Gill (CPMG) pulse-sequence, and decay of the transverse magnetization in a resonant spin-locking field are investigated. Both trajectory and stochastic operator approaches are formulated and shown to be completely equivalent whenever the dynamics of population transfers among the inequivalent sites is governed by either a stationary or a nonstationary Markov process. A nonstationary Markov process may result from Brownian dynamics (a stationary Markov process) in a larger conformational space that contains the subspace of inequivalent sites. A continuous Gaussian exchange model is formulated in which a nucleus undergoes continuous one-dimensional motion in a harmonic potential well that is located in a linear chemical shift gradient. The effects of this Gaussian exchange model on the lineshape, simple spin-echo decay, and decay of the even echoes of a CPMG pulse train are treated rigorously via the trajectory approach. Compact analytical expressions are obtained for the relevant correlation functions in each case. The relevant decays are found to be exponential in the very short time and long time limits, which are not necessarily experimentally significant in any given case. In the fast exchange limit the relevant decays are exponential at all times, and explicit formulas are given for their decay rates. In the long time limit, all discrete multisite models with the same intrinsic Ro2 at every site are shown to be completely equivalent to a continuous Gaussian model with appropriate relaxation time and variance of the Larmor frequency. The effects of this Gaussian exchange model on the decay of the transverse magnetization in a resonant spin-locking field are treated heuristically by a trajectory approach. The intrinsic contribution (Ro1rho) of rapid rotations and dipole-dipole interactions to relax the transverse magnetizations of two nuclei of the same kind in the presence of a (nearly) resonant spin-locking field is also derived and found to be practically the same as the intrinsic contribution, Ro2, of those same rotations to the simple and CPMG spin-echo decay rates and linewidth. Literature data for the linewidth, decay rate of the CPMG even spin-echoes, and R(1rho) decay rate for the A9-H2 protons of adenines at the central TpA step in the sequence, 5'-GCAGGTTTAAACCTCG-3', are analyzed using the Gaussian exchange model to assess the time-scale and variance of the site exchange process as well as the intrinsic Ro2 rate. Although a single Gaussian exchange process with appropriate parameters can fit these three A9-H2 data rather well, this particular "solution" cannot be reconciled with NMR relaxation data on other protons in the same DNA molecule. Rather good agreement with all of the observations is obtained by using a model of two concurrent Gaussian exchange processes, whose relaxation times, tau = 7 and 460 micros, differ in time-scale by a factor of 65. The insensitivity of R1rho in the presence of a fast site exchange process to much slower concurrent site exchange processes is explicitly demonstrated. Protocols for detecting and characterizing a second slow site exchange process are suggested.

Large-amplitude picosecond anisotropy decay of the intrinsic fluorescence of double-stranded DNA

The conformational flexibility of the DNA double helix is of great interest because of its potential role in protein recognition, packaging into chromosomes, formation of photodefects, and interaction with drugs. Theory finds that DNA is very flexible; however, there is a scarcity of experimental results that examine intrinsic properties of the DNA bases for the inherent flexibility in solution. We have studied the dynamics of poly(dA).poly(dT) and (dA)20.(dT)20 in a 50 mM cacodylate, 0.1 M NaCl, pH 7 buffer by using the time-correlated picosecond fluorescence anisotropy of thymine selectively excited at 293 nm. For both nucleic acids, a large-amplitude biphasic decrease in the anisotropy is observed that has a very fast, large-amplitude component on the picosecond time scale and a slower, smaller-amplitude component on the nanosecond time scale. These modes are sensitive to sucrose concentration, and are greatly attenuated at 77% sucrose by volume. This observation suggests that motions of the bases make a significant contribution to the observed fluorescence depolarization (in the absence of sucrose). Measurements on the single-stranded systems poly(dT) and (dT)20 reveal a much smaller amplitude of the very fast depolarization mode. These observations are consistent with a mechanism that involves concerted motions in the interior of the double-stranded systems.

The amplitude of local angular motion of purines in DNA in solution

Nuclear magnetic resonance and optical experiments are combined to determine the rms amplitude of local angular motion of purines in DNA in solution. A 12 base-pair duplex DNA with the sequence d(CGCGAATTCGCG)2 is deuterated at the H8 positions of adenine and guanine by exchange with solvent at 55 degrees C. The deuterium nmr spectrum of this DNA is measured at 30 mg/mL at 30 degrees C in an 11.76 Tesla magnet (76.75 MHz). The time-resolved fluorescence polarization anisotropies (FPA) of this same sample and also a greatly diluted sample (0.215 mg/mL) were measured after addition of ethidium. FPA measurements of the dilute sample yield the hydrodynamic radius, RH = 9.94 +/- 0.2 A, while those at the nmr concentration are employed to characterize the collective motions in terms of either an enhanced viscosity or dimer formation. The rms amplitude of local angular motion was determined by analyzing the 2H-nmr spectrum, in particular the line width, using recently developed theory for the transverse relaxation rate (RQ2) together with essential information about the collective motions from these and other optical studies. When the principal-axis frame of the electric field gradient tensor is assumed to undergo overdamped libration around each of its three body-fixed axes in an isotropic deflection potential, then the rms amplitude of local angular motion around any single axis is found to lie in the range 10 degrees-11 degrees, provided the high DNA concentration acts to enhance the viscosity, and is about 9 degrees-11 degrees, if it acts to produce end-to-end dimers. The proton nmr relaxation data of Eimer et al. are reanalyzed and shown to yield an rms amplitude of angular motion of the cytosine H5-H6 internuclear vector of 9 degrees-10 degrees, depending upon its orientation with respect to the helix axis. In all of these analyses, full account is taken of the collective twisting and bending deformations, which have a small but significant effect on the results. It is shown that the rms amplitudes of local angular motion do not depend strongly on the model (potential), provided that isotropic rotation around the same number of axes is allowed and that one compares rms angles of the same dimensionality. The rms amplitudes of local angular motion in solution are comparable to those observed for the same sequence at low levels of hydration in the solid state.

Dynamics in synthetic oligonucleotides. A solid-state deuterium NMR study

Solid-state 2H NMR spectroscopy was used to investigate dynamics in the [methyl-2H]thymidine-labeled oligonucleotide [d(CGCGAAT*T*CGCG)]2. Quadrupole echo line shapes, spin-lattice relaxation, and quadrupolar echo decay rates were investigated as a function of hydration W (moles of water/moles of nucleotide) between 0 and approximately 30. The amplitude of the base motion, modeled as a fast four-site libration, or diffusion in a cone, increased slightly with higher levels of hydration. A slower component of motion about the helix axis appears at W approximately 10 and increases in rate and amplitude, leading to the intermediate rate line shape observed at W approximately 21.

Deuterium relaxation and internal motion in solid Li-DNA

As part of an effort to explore the nature of the internal motion in solid polynucleotides, the spectral densities of motion J1(omega 0) and J2(2 omega 0) have been measured for oriented, partially hydrated samples of calf thymus Li-DNA deuterated in the guanine and adenine 8-positions. Both spectral densities increase with increasing hydration level, J1 is found to be 2-5 times larger than J2, and their frequency dependence appears to be omega-1 and omega-3/2, respectively. The large values of the ratio J1/J2 rule out any in-plane torsional motion as the dominant relaxation mechanism in these samples, but a drop in this ratio at high hydration levels (G13 H2O/nucleotide) may indicate increasing contributions from such torsional motion. Although a satisfactory fit to a particular motional model has yet to be achieved, our findings show that the librational motion of the C8-D bond at or below a hydration level of 10 H2O/nucleotide is approximately uniaxial, with correlation times for the motion in the range 0.2-3.0 microseconds.

NMR distance measurements in DNA duplexes: sugars and bases have the same correlation times

To evaluate whether the sugar moieties of short DNA duplexes exhibit local motion of sufficient amplitude to affect interproton distance measurements, we have carried out a series of time-dependent NOESY experiments at increasingly shorter mixing times on dodecamer DNA duplexes. By use of the cytosine H5-H6 vector as a known distance in the bases and the geminal 2'H-2''H vector as a known distance in the sugars, the corresponding apparent cross-relaxation rates were sampled at various mixing times. While the ratio of the inverse sixth power of these two fixed distances is in the range 6-7, when the system is sampled at 100 ms the apparent initial rate of growth of the 2'H-2''H NOESY crosspeak is only 1.9-2.0 times faster than that of the H5-H6 crosspeak--in agreement with the results of Clore and Gronenborn [Clore, G. M., & Gronenborn, A. M. (1984) FEBS Lett. 172, 219; (1984) FEBS Lett. 175, 117] and of Gronenborn and Clore [Gronenborn, A. M., & Clore, G. M. (1985) Prog. NMR Spectrosc. 17, 1]. This observation was interpreted to indicate the existence of internal mobility with a 3-fold shorter correlation time for the sugar moieties in DNA and led to the use of this shorter correlation time to estimate sugar-sugar proton distances and many sugar-base proton distances in subsequent DNA structure determination. We have examined 2'H-2"H cross-relaxation and H5-H6 cross-relaxation at 100, 90, 60, 30, and 15 ms in dodecamer DNA duplexes.(ABSTRACT TRUNCATED AT 250 WORDS)

Interaction of chloroquine with linear and supercoiled DNAs. Effect on the torsional dynamics, rigidity, and twist energy parameter

The magnitude and uniformity of the torsion elastic constant (alpha) of linear pBR322 DNA and supercoiled pBR322 DNAs with high-twist (sigma = -0.083) and normal-twist (sigma = -0.48) are measured in 0.1 M NaCl as a function of added chloroquine/base-pair ratio (chl/bp) by studying the fluorescence polarization anisotrophy (FPA) of intercalated ethidium dye. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single-photon counting detection. A general theory is developed for the binding of ligands that unwind superhelical DNAs, and the simultaneous binding of two different intercalators is treated in detail. The equilibrium constant (K) for binding chloroquine to linear pBR322 DNA and the number (r) of bound chloroquines per base pair are determined from the relative amplitude ratio of the slow (normally intercalated) and fast (free) components in the decay of the (probe) ethidium fluorescence intensity as a function of chl/bp. For chloroquine binding to supercoiled pBR322 DNAs, the intrinsic binding constant is assumed to be the same as for the linear DNA, but the twist energy parameter ET (N times the free energy to change the linking number from 0 to 1 in units of kBT) is regarded as adjustable. Using the best-fit ET, the binding ratios r are calculated for each chl/bp ratio. Twist energy parameters are also determined for ethidium binding to these supercoiled DNAs by competitive dialysis. For chloroquine binding, we obtain ET = 360 and 460 respectively for the normal-twist and high-twist supercoiled DNAs. For ethidium binding the corresponding values are ET = 280 +/- 70 and 347 +/- 50. Like other dye-binding values, these are substantially lower than those obtained by ligation methods. In the absence of chloroquine, the torsion constants of all three DNAs are virtually identical, alpha = (5.0 +/- 0.4) x 10(-12) dyn.cm. For linear pBR322 DNA, the magnitude and uniformity of alpha remain unaltered by intercalated chloroquine up to r = 0.19. This finding argues that the FPA is not significantly relaxed by diffusion of any kinks or solitons. If alpha d denotes the torsion constant between a dye and a base pair and alpha 0 that between two base pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65-1.64, with a most probable value of 1.0.(ABSTRACT TRUNCATED AT 400 WORDS)

Three-dimensional structure of the wild-type lac Pribnow promoter DNA in solution. Two-dimensional nuclear magnetic resonance studies and distance geometry calculations

The solution structure of a 12 base-pair DNA duplex containing the wt-lac promoter Pribnow sequence TATGTT has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Proton assignments for the 24 sugar and base residues were obtained from two-dimensional correlated nuclear magnetic resonance and two-dimensional nuclear Overhauser effect spectra in both 2H2O and H2O, and by two-dimensional relayed coherence transfer nuclear magnetic resonance spectroscopy experiments. Time-dependent, two-dimensional nuclear Overhauser effect spectra were used to determine the initial cross-relaxation rates between 212 pairs of assigned protons, leading to 212 interproton distances in the double helix (8 to 9 per nucleotide). These distance constraints, and known bond lengths and angles, were entered into a distance matrix. After smoothing the bounds of the distance matrix, 12 trial matrices within the bounds constraints were independently generated and embedded in three-dimensional space using a distance geometry algorithm, to generate 12 trial structures. These trial structures were then refined until they no longer violated the distance matrix. The resulting structures are very similar at the local base-pair and nearest-neighbor base-pair level, but exhibit increasing variation at more distant and global levels. At the nearest-neighbor level, the A to T step and the G to T step within the Pribnow hexamer, as well as the G to T step preceding the hexamer, all exhibit very low screw pitch, i.e. 5(+/- 6) degrees. Conversely, the T to G step in the center of the promoter has a large screw pitch (47(+/- 2) degrees) and the T to G step at the 3' end of the promoter has a very large screw pitch (60(+/- 3) degrees). The limitations of nuclear magnetic resonance spectroscopy distance determination of structure are discussed in terms of resolution and spectral overlap of two-dimensional nuclear Overhauser effect crosspeaks. In the present duplex, the inability to measure several 1'-2' and 1'-2" distances resulted in underdetermination of the precise local sugar conformation for seven of the 24 residues, although the spatial position of all sugars was well defined.

A fluorescence photobleaching study of the microsecond reorientational motions of DNA

We have conducted a polarized fluorescence photobleaching recovery (FPR) study of the rotational dynamics of ethidium azide labeled DNA. Polarized photobleaching experiments provide data on microsecond and millisecond molecular reorientation that complement the information available from nanosecond fluorescence depolarization studies. In polarized FPR experiments an anisotropic angular concentration of fluorophore is created by bleaching dye molecules in a preferred orientation with a short, intense pulse of polarized light. The sample is then weakly illuminated, and the temporal variation in the emitted fluorescence is monitored. The fluorescence signal will systematically change as molecules undergo post-bleach reorientation and the angular distribution of dye tends toward isotropy. We have observed that the time dependence of our microsecond FPR curves is also determined in part by nonrotational phenomena. To isolate the reorientational recovery we conduct our FPR experiments in two modes (called parallel and perpendicular) that differ only in the polarization of the bleaching light. A quotient function, R(t), is constructed from the data obtained in these two modes; the variation with time of this new quantity is governed solely by processes that are sensitive to the polarization of the incident light (e.g., molecular rotation). It is found experimentally that R(t) remains constant, as expected, for rotationally restricted DNA systems despite a temporal recovery in the parallel and perpendicular FPR curves. We also follow the dynamics of solutions of phage lambda DNA as revealed in the temporal dependence of R(t). This DNA system rotationally relaxes after approximately 100 microseconds and the dye/DNA complex reorients substantially during the 10-microseconds bleach period. Our FPR data are interpreted in terms of dynamic models of DNA motion.

Local base dynamics and local structural features in RNA and DNA duplexes

Local base motion and local structural base information are derived with a simple motional model from site specifically spin-labeled polyribo- and polydeoxyribonucleotides. The model was developed earlier for some nucleic acids and has now been applied to analyze 22 different nucleic acid systems. We conclude that the base motion of the spin-labeled nucleotide in single-stranded RNA, DNA, or non-base-paired bases in duplexes is of the order of 1 ns and that its base mobility decreases by about a factor of 4 upon base pairing. Also, the tether motion of the probe is slower in an RNA than in a DNA duplex.

31P-NMR analysis of the B to Z transition in double-stranded (dC-dG)3 and (dC-dG)4 in high salt solution

In 4M NaCl solutions (dC-dG)n (n = 3,4; approximately 9 mM) exist as a mixture o +/- B and Z forms. The low and high field components of two 31P NMR resonances originating from internal phosphodiester groups are assigned to the GpC and CpG linkages, respectively. Low temperatures stabilize the Z-forms, which completely disappear above 50 degrees C (n = 3) and 65 degrees C (n = 4). delta H = -44 and -17 kJ/mol for B to Z transition in the hexamer and octamer duplexes, respectively. Temperature dependent changes (0-50 degrees C range) in the spin-lattice relaxation times at 145.7 MHz are distinctly different for the 31P nuclei o +/- GpC and CpG groups. The relaxation data can be explained by assuming that the GpC phosphodiester groups undergo more local internal motion than do the CpG groups.

NMR studies of conformational states and dynamics of DNA

The application of high resolution NMR techniques to the investigation of DNA double helices in solution is currently in a rapid state of change as a result of advances in three different fields. First, new methods (cloning, enzymatic degradation, sonication, and chemical synthesis) have been developed for producing large quantities of short DNA suitable for NMR studies. Second, there have been major advances in the field of NMR in terms of the introduction of new pulse techniques and improvements in instrumentation. Finally, as a result of recent X-ray diffraction studies on short DNA helices and the discovery of left-handed Z-DNA there is heightened interest in the study of DNA structures in solution and the effect of sequence on structure. In the present review, we discuss the way in which NMR techniques have been used to probe various aspects of the DNA properties, including base pairing structure, dynamics of breathing, effect of sequence on DNA structure, internal molecular motions, the effect of environment on the DNA, and the interaction of DNA with small ligands.

Understanding DNA conformational dynamics: answering questions and questioning answers

Phosphorus-31 and especially Carbon-13 NMR measurements have recently become primary input to the understanding of DNA solution dynamics. While the 31P measurements are inherently easier, the quality of 31P dynamics information is suspect and therefore 13C measurements are preferred. In fact, it is necessary to obtain several kinds of 13C data (T1s, NOE's, linewidths, integrated peak intensities) over a wide range of magnetic fields (13C NMR frequencies) in order to identify major features of DNA internal motions. Further information comes from variation of temperature and DNA fragment length and/or concentration. Most of our 13C measurements have been performed at 37.7-90.6 MHz on fully double stranded monomer size (147 base pair) DNA at concentrations in phosphate buffer of approximately less than 10 to approximately greater than 200 mg ml-1; temperatures studied range from 6 to 55 degrees C. Other measurements have been performed on monomer-size single-stranded DNA at 85 and 92 degrees. The large data set we have acquired appears to answer some important questions about the nature and extent of DNA overall and internal motional dynamics. However, the picture remains incomplete and a number of questions arise from these results: 1. Overall motion of the double stranded DNA fragments follows expected hydrodynamic behavior; 2. Restricted but rapid internal motion along the DNA structure is well represented by a spaghetti-like wobbling-in-a-cone model; 3. DNA-DNA Interactions and solvent ordering, present at relatively low DNA concentrations, partially quench the internal motion, consistent with hinge-model structural changes (and the spaghetti model above) but not as compatible with in-plane torsional motion models; 4. The deoxyribose C-2' sites undergo additional motion which is partially uncoupled from the internal wobbling motions: 5. At high DNA concentrations, a phase transition occurs, resulting in ordered structures which drastically affect DNA internal dynamics; 6. DNA interacting with ethidium does not greatly change its conformational mobility; 7. DNA interacting with Hg2+ ions shows less than anticipated change in internal DNA dynamics. The remaining challenge is to interpret our current results in terms of specific conformational processes and to understand why the conformational mobility of double stranded DNA is relatively unhindered by major structural perturbants such as intercalating ethidium and mercury ion.

New view of lipid bilayer dynamics from 2H and 13C NMR relaxation time measurements

Natural abundance 13C spin-lattice (T1) relaxation time measurements are reported for unilamellar vesicles of 1,2-dipalmitoylphosphatidylcholine (1,2-dipalmitoyl-sn-glycero-3-phosphocholine), in the liquid crystalline phase, at magnetic field strengths of 1.40, 1.87, 2.35, 4.23, 7.05, 8.45, and 11.7 tesla (resonance frequencies of 15.0, 20.0, 25.1, 45.3, 75.5, 90.5, and 126 MHz, respectively), and the results are compared to previous 2H T1 studies of multilamellar dispersions. For both the 13C and 2H T1 studies, a dramatic frequency dependence of the relaxation was observed. At superconducting magnetic field strengths (4.23-11.7 tesla), plots of the 13C T1(-1) relaxation rates as a function of acyl chain segment position clearly reveal the characteristic "plateau" signature of the liquid crystalline phase, as found previously from 2H NMR studies. The dependence of T1(1) on ordering, determined previously from 2H NMR, and the T1(-1) dependence on frequency, determined from both 13C and 2H NMR studies, suggest that a unified picture of the bilayer molecular dynamics can be provided by a simple relaxation law of the form T1(-1) approximately equal to A tau f + BS2C-H omega -1/2(0). In the above expression, A and B are constants, SC-H (= SC-D) is the bond segmental order parameter, and omega 0 is the nuclear Larmor frequency. The first (A) term includes contributions from fast, local segmental motions characterized by the effective correlation time tau f, whereas the second (B) term describes slower, collective fluctuations in the local ordering. The value of tau f approximately equal to 10(-11) sec, obtained by extrapolating T1(-1) to infinite frequency, suggests that the segmental microviscosity of the bilayer hydrocarbon region does not differ appreciably from that of the equivalent n-paraffinic liquids of similar chain length.