A Solid-State Deuterium NMR Study of the Localized Dynamics at the C9pG10 Step in the DNA Dodecamer [d(CGCGAATTCGCG)]2

Solid-state NMR studies of nucleic acid components

Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.

Investigation of the effect of a variety of pulse errors on spin I=1 quadrupolar alignment echo spectroscopy

We report on an analysis of a well known three-pulse sequence for generating and detecting spin I=1 quadrupolar order when various pulse errors are taken into account. In the situation of a single quadrupolar frequency, such as the case found in a single crystal, we studied the potential leakage of single and/or double quantum coherence when a pulse flip error, finite pulse width effect, RF transient or a resonance offset is present. Our analysis demonstrates that the four-step phase cycling scheme studied is robust in suppressing unwanted double and single quantum coherence as well as Zeeman order that arise from the experimental artifacts, allowing for an unbiased measurement of the quadrupolar alignment relaxation time, T(1Q). This work also reports on distortions in quadrupolar alignment echo spectra in the presence of experimental artifacts in the situation of a powdered sample, by simulation. Using our simulation tool, it is demonstrated that the spectral distortions associated with the pulse artifacts may be minimized, to some extent, by optimally choosing the time between the first two pulses. We highlight experimental results acquired on perdeuterated hexamethylbenzene and polyethylene that demonstrate the efficacy of the phase cycling scheme for suppressing unwanted quantum coherence when measuring T(1Q). It is suggested that one employ two separate pulse sequences when measuring T(1Q) to properly analyze the short time behavior of quadrupolar alignment relaxation data.

Solid state 2H NMR analysis of furanose ring dynamics in DNA containing uracil

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions, such as uracil, in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Increased dynamic properties in lesion-containing DNAs have been suggested to assist recognition and specificity. Deuterium solid-state nuclear magnetic resonance (SSNMR) has been used to directly observe local dynamics of the furanose ring within a uracil:adenine (U:A) base pair and compared to a normal thymine:adenine (T:A) base pair. Quadrupole echo lineshapes, <T(1Z)>, and <T(2e)> relaxation data were collected, and computer modeling was performed. The results indicate that the relaxation times are identical within the experimental error, the solid lineshapes are essentially indistinguishable above the noise level, and our lineshapes are best fit with a model that does not have significant local motions. Therefore, U:A base pair furanose rings appear to have essentially identical dynamic properties as a normal T:A base pair, and the local dynamics of the furanose ring are unlikely to be the sole arbiter for uracil recognition and specificity in U:A base pairs.

An FTIR investigation of flanking sequence effects on the structure and flexibility of DNA binding sites

Fourier transform infrared (FTIR) spectroscopy and a library of FTIR marker bands have been used to examine the structure and relative flexibilities conferred by different flanking sequences on the EcoRI binding site. This approach allowed us to examine unique peaks and subtle changes in the spectra of d(AAAGAATTCTTT)(2), d(TTCGAATTCGAA)(2), and d(CGCGAATTCGCG)(2) and thereby identify local changes in base pairing, base stacking, backbone conformation, glycosidic bond rotation, and sugar puckering in the studied sequences. The changes in flanking sequences induce differences in the sugar puckers, glycosidic bond rotation, and backbone conformations. Varying levels of local flexibility are observed within the sequences in agreement with previous biological activity assays. The results also provide supporting evidence for the presence of a splay in the G(4)-C(9) base pair of the EcoRI binding site and a potential pocket of flexibility at the G(4) cleavage site that have been proposed in the literature. In sum, we have demonstrated that FTIR is a powerful methodology for studying the effect of flanking sequences on DNA structure and flexibility, for it can provide information about the local structure of the nucleic acid and the overall relative flexibilities conferred by different flanking sequences.

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.

Conformational dynamics of an intact virus: order parameters for the coat protein of Pf1 bacteriophage

This study has examined the atomic-level dynamics of the protein in the capsid of filamentous phage Pf1. This capsid consists of approximately 7,300 small subunits of only 46 aa in a helical array around a highly extended, circular single-stranded DNA molecule of 7,349 nt. Measurements were made of site-specific, solid-state NMR order parameters, S, the values which are dimensionless quantities between 0 (mobile) and 1 (static) that characterize the amplitudes of molecular bond angular motions that are faster than microseconds. It was found that the protein subunit backbone is very static, and of particular interest, it appears to be static at residues glycine 15 and glutamine 16 where it had been previously thought to be mobile. In contrast to the backbone, several side chains display large-amplitude angular motions. Side chains on the virion exterior that interact with solvent are highly mobile, but surprisingly, the side chains of residues arginine 44 and lysine 45 near the DNA deep in the interior of the virion are also highly mobile. The large-amplitude dynamic motion of these positively charged side chains in their interactions with the DNA were not previously expected. The results reveal a highly dynamic aspect of a DNA-protein interface within a virus.

13C relaxation studies of the DNA target sequence for hhai methyltransferase reveal unique motional properties

The goal of this work was to examine if sequence-dependent conformational flexibility in DNA plays a role in base extrusion, a common conformational change induced by many DNA-modifying enzymes. We studied the dynamics of the double-stranded DNA target of the HhaI methyltransferase by recording an extensive set of (13)C NMR relaxation parameters. We observe that the cytidine furanose rings experience fast (picosecond to nanosecond) motions that are not present in other nucleotides; the methylation site experiences particularly high mobility. We also observe that the bases of guanosine and cytidine residues within the HhaI recognition sequence GCGC experience motions on a much slower (1-100 micros) time scale. We compare these observations with previous solution and solid-state NMR studies of the EcoRI nuclease target sequence, and solid-state NMR studies of a similar HhaI target construct. While an increased mobility of cytidine furanose rings compared to those of other nucleotides is observed for both sequences, the slower motions are only observed in the HhaI target DNA. We propose that this inherent flexibility lowers the energetic barriers that must occur when the DNA binds to the HhaI methyltransferase and for extrusion of the cytidine prior to its methylation.

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.

Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site

The dynamics of the furanose rings in the GCGC moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied by using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs selectively deuterated on the furanose rings of nucleotides within the 5'-GCGC-3' moiety indicated that all of these positions are structurally flexible. The furanose ring within the deoxycytidine that is the methylation target displays the largest-amplitude structural changes according to the observed deuterium NMR line shapes, whereas the furanose rings of nucleotides more remote from the methylation site have less-mobile furanose rings (i.e., with puckering amplitudes < 0.3 A). Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through the furanose ring. These NMR data indicate that the 5'-GCGC-3' is dynamic, with the largest-amplitude motions occurring nearest the methylation site. The inherent flexibility of this moiety in DNA makes the molecule more amenable to the large-amplitude structural rearrangements that must occur when the DNA binds to the HhaI methyltransferase.

Simulations of molecular dynamics in solid-state NMR spectra of spin-1 nuclei including effects of CSA- and EFG-terms up to second order

By numerical simulations MAS and QCPMG methods for acquiring spectra of spin-1 nuclei were compared in order to determine the most sensitive experiment for analysis of molecular dynamics. To comply with the large quadrupolar constants for 14N and the CSA reported for 6Li both of these interactions are included up to second order. For 2H and 6Li both QCPMG and single-pulse MAS experiments were suitable for dynamics studies whereas the single-pulse MAS experiment were the method of choice for investigation of 14N dynamics for C(Q)'s larger than 750kHz at 14.1T. This property prohibits excitation of the 14N lineshape using either single hard or softer composite rf-pulses. Focusing on 14N it was demonstrated that the centerband lineshape is sensitive toward both off-MAS and CSA effects. In addition, excitation by real-time pulses showed that proper lineshapes corresponding to a site with a C(Q) of 3MHz may be excited by a very short pulse.

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.

Bending and flexibility of methylated and unmethylated EcoRI DNA

We used cyclization kinetics experiments and Monte Carlo simulations to determine a structural model for a DNA decamer containing the EcoRI restriction site. Our findings agree well with recent crystal and NMR structures of the EcoRI dodecamer, where an overall bend of seven degrees is distributed symmetrically over the molecule. Monte Carlo simulations indicate that the sequence has a higher flexibility, assumed to be isotropic, compared to that of a "generic" DNA sequence. This model was used as a starting point for the investigation of the effect of cytosine methylation on DNA bending and flexibility. While methylation did not affect bend magnitude or direction, it resulted in a reduction in bending flexibility and under-winding of the methylated nucleotides. We demonstrate that our approach can augment the understanding of DNA structure and dynamics by adding information about the global structure and flexibility of the sequence. We also show that cyclization kinetics can be used to study the properties of modified nucleotides.

Modeling furanose ring dynamics in DNA

Determination of the conformational flexibility of the furanose ring is of vital importance in understanding the structure of DNA. In this work we have applied a model of furanose ring motion to the analysis of deuterium line shape data obtained from sugar rings in solid hydrated DNA. The model describes the angular trajectories of the atoms in the furanose ring in terms of pseudorotation puckering amplitude (q) and the pseudorotation puckering phase phi. Fixing q, the motion is thus treated as Brownian diffusion through an angular-dependent potential U(phi). We have simulated numerous line shapes varying the adjustable parameters, including the diffusion coefficient D, pseudorotation puckering amplitude q, and the form of the potential U(phi). We have used several forms of the potential, including equal double-well potentials, unequal double-well potentials, and a potential truncated to "second order" in the Fourier series. To date, we have obtained best simulations for both equilibrium and nonequilibrium (partially relaxed) solid-state deuterium NMR line shapes for the sample [2' '-2H]-2'-deoxycytidine at the position C3 (underlined) in the DNA sequence [d(CGCGAATTCGCG)]2, using a double-well potential with an equal barrier height of U(0) = 5.5k(B)T ( approximately 3.3 kcal/mol), a puckering amplitude of q = 0.4 A, and a diffusion coefficient characterizing the underlying stochastic jump rate D = 9.9 x 10(8) Hz. Then the rate of flux for the C-D bond over the barrier, i.e., the escape velocity or the overall rate of puckering between modes, was found to be 0.7 x 10(7) Hz.

Dynamic impact of methylation at the M. Hhai target site: a solid-state deuterium NMR study

Base methylation plays an important role in numerous biological functions of DNA, from inhibition of cleavage by endonucleases to inhibition of transcription factor binding. Studies of nucleic acid structure have shown little differences in unmethylated DNAs and the identical sequence containing methylated analogues. We have investigated changes in the local dynamics of DNA upon substitution of a methylated cytosine analogue for cytosine using solid-state deuterium NMR. In particular, we have observed changes in the local dynamics at the target site of the M. HhaI restriction system. These studies observe changes in the amplitudes of the local backbone dynamics at the actual target site of the HhaI methyltransferase. This conclusion is another indication that the significant result of base methylation is to perturb the local dynamics, and therefore the local conformational flexibility, of the DNA helix, inhibiting or restricting the protein's ability to manipulate the DNA helix in order to perform its chemical alterations.

The dynamic impact of CpG methylation in DNA

Solid-state deuterium NMR is used to investigate perturbations of the local, internal dynamics in the EcoRI restriction binding site, -GAATTC- induced by cytidine methylation. Methylation of the cytidine base in this sequence is known to suppress hydrolysis by the EcoRI restriction enzyme. Previous solid-state deuterium NMR studies have detected large amplitude motions of the phosphate-sugar backbone at the AT-CG junction of the unmethylated DNA sequence. This study shows that methylation of the cytidine base in a CpG dinucleotide reduces the amplitudes of motions of the phosphate-sugar backbone. These observations suggest a direct link between suppression of the amplitudes of localized, internal motions of the sugar-phosphate backbone of the DNA and inhibition of restriction enzyme cleavage.

Helix morphology changes in B-DNA induced by spontaneous B(I)<==>B(II) substrate interconversion

Investigations of spontaneous, i.e. not forced, B-DNA's B(I)<==>B(II) substate transitions are carried out on the d(CGCGAATTCGCG)2 EcoRI dodecamer sequence using Molecular Dynamics Simulations. Analysis of the resulting transition processes with respect to the backbone angles reveals concerted changes not only for backbone angles epsilon, zeta, and beta, but also for the 5'-delta and 5'-chi angles. For alpha and delta inside the interconverting base step, a change is seen in short lived B(II) conformers. With respect to base morphology distinct changes are observed for buckle, propeller twist, shift, roll and twist, as well as x-displacement and tip. The base mainly involved in the changes is identified as the base preceding the interconverting phosphate. Altogether single B(I)<==>B(II) interconversions result only in local distortions represented by the larger spread of most parameters. Comparison of the atomic positional fluctuations derived from the simulation with those obtained from the static X-ray structure results in striking similarities.