DNA duplex dynamics: NMR relaxation studies of a decamer with uniformly 13C-labeled purine nucleotides

Base Dynamics in the HhaI Protein Binding Site

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

Structure and Dynamics of the Tetra-A Loop and (A-A)-U Sequence Motif within the Coliphage GA Replicase RNA Operator

The three-dimensional structure of a RNA hairpin containing the RNA operator binding site for bacteriophage GA coat protein is presented. The phage GA operator contains the asymmetric (A-A)-U sequence motif and is capped by a four-adenine (tetra-A) loop. The uridine of the (A-A)-U motif preferentially pairs with the 5'-proximal cross-strand adenine, and the 3'-proximal adenine stacks into the helix. The tetra-A loop is well-ordered with adenine residues 2-4 forming a 3' stack. This loop conformation stands in contrast to the structure of the 5'-AUUA loop of the related phage MS2 operator in which residues 1 and 2 form a 5' stack. The context dependence of the (A-A)-U sequence motif conformation was examined using structures of 76 unique occurrences from the Protein Data Bank. The motif almost always has one adenine bulged and the other adenine adopting an A-U base pair. In the case in which the (A-A)-U motif is flanked by only one Watson-Crick base pair, the adenine adjacent to the flanking base pair tends to bulge; 80% of motifs with a 3' flanking pair have a 3' bulged adenine, and 84% of motifs with a 5' flanking pair have a 5' bulged adenine. The frequencies of 3'- and 5'-proximal adenines bulging are 33 and 67%, respectively, when the (A-A)-U motif is flanked by base pairs on both sides. Although a 3' flanking cytidine correlates (88%) with bulging of the 5'-proximal adenine, no strict dependence on flanking nucleotide identity was identified for the 5' side.

Recent Methods for Purification and Structure Determination of Oligonucleotides

Aptamers are single-stranded DNA or RNA oligonucleotides that can interact with target molecules through specific three-dimensional structures. The excellent features, such as high specificity and affinity for target proteins, small size, chemical stability, low immunogenicity, facile chemical synthesis, versatility in structural design and engineering, and accessible for site-specific modifications with functional moieties, make aptamers attractive molecules in the fields of clinical diagnostics and biopharmaceutical therapeutics. However, difficulties in purification and structural identification of aptamers remain a major impediment to their broad clinical application. In this mini-review, we present the recently attractive developments regarding the purification and identification of aptamers. We also discuss the advantages, limitations, and prospects for the major methods applied in purifying and identifying aptamers, which could facilitate the application of aptamers.

Stable isotope labeling methods for DNA

NMR is a powerful method for studying proteins and nucleic acids in solution. The study of nucleic acids by NMR is far more challenging than for proteins, which is mainly due to the limited number of building blocks and unfavorable spectral properties. For NMR studies of DNA molecules, (site specific) isotope enrichment is required to facilitate specific NMR experiments and applications. Here, we provide a comprehensive review of isotope-labeling strategies for obtaining stable isotope labeled DNA as well as specifically stable isotope labeled building blocks required for enzymatic DNA synthesis.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Enzymatic preparation of multimilligram amounts of pure single-stranded DNA samples for material and analytical sciences

We present a method for high-yield production of multimilligram amounts of pure single-stranded DNA employing rolling circle amplification (RCA) and processing by restriction enzymes. Pure and homogeneous samples are produced with minimal handling time, reagents, and waste products. The RCA method is more than twice as efficient in dNTP incorporation than conventional polymerase chain reaction in producing end product. The validity and utility of the method are demonstrated in the production of a uniformly (13)C/(15)N-labeled 38-nt cocaine aptamer DNA used in nanosensing devices.

Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC)

BACKGROUND

The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA.

METHODS

Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms.

RESULTS

These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44μs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5μs timescale.

CONCLUSIONS

We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the μs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine.

GENERAL SIGNIFICANCE

With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.

Role of microscopic flexibility in tightly curved DNA

The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari–Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Conformational analysis by quantitative NOE measurements of the β-proton pairs across individual disulfide bonds in proteins

NOEs between the β-protons of cysteine residues across disulfide bonds in proteins provide direct information on the connectivities and conformations of these important cross-links, which are otherwise difficult to investigate. With conventional [U-(13)C, (15)N]-proteins, however, fast spin diffusion processes mediated by strong dipolar interactions between geminal β-protons prohibit the quantitative measurements and thus the analyses of long-range NOEs across disulfide bonds. We describe a robust approach for alleviating such difficulties, by using proteins selectively labeled with an equimolar mixture of (2R, 3S)-[β-(13)C; α,β-(2)H(2)] Cys and (2R, 3R)-[β-(13)C; α,β-(2)H(2)] Cys, but otherwise fully deuterated. Since either one of the prochiral methylene protons, namely β2 (proS) or β3 (proR), is always replaced with a deuteron and no other protons remain in proteins prepared by this labeling scheme, all four of the expected NOEs for the β-protons across disulfide bonds could be measured without any spin diffusion interference, even with long mixing times. Therefore, the NOEs for the β2 and β3 pairs across each of the disulfide bonds could be observed at high sensitivity, even though they are 25% of the theoretical maximum for each pair. With the NOE information, the disulfide bond connectivities can be unambiguously established for proteins with multiple disulfide bonds. In addition, the conformations around disulfide bonds, namely χ(2) and χ(3), can be determined based on the precise proton distances of the four β-proton pairs, by quantitative measurements of the NOEs across the disulfide bonds. The feasibility of this method is demonstrated for bovine pancreatic trypsin inhibitor, which has three disulfide bonds.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.

RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis

We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.

Preparation, resonance assignment, and preliminary dynamics characterization of residue specific 13C/15N-labeled elongated DNA for the study of sequence-directed dynamics by NMR

DNA is a highly flexible molecule that undergoes functionally important structural transitions in response to external cellular stimuli. Atomic level spin relaxation NMR studies of DNA dynamics have been limited to short duplexes in which sensitivity to biologically relevant fluctuations occurring at nanosecond timescales is often inadequate. Here, we introduce a method for preparing residue-specific (13)C/(15)N-labeled elongated DNA along with a strategy for establishing resonance assignments and apply the approach to probe fast inter-helical bending motions induced by an adenine tract. Preliminary results suggest the presence of elevated A-tract independent end-fraying internal motions occurring at nanosecond timescales, which evade detection in short DNA constructs and that penetrate deep (7 bp) within the DNA helix and gradually fade away towards the helix interior.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

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.

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.

NMR and MD studies of the temperature-dependent dynamics of RNA YNMG-tetraloops

In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.

Chemical shift tensors of protonated base carbons in helical RNA and DNA from NMR relaxation and liquid crystal measurements

Knowledge of (13)C chemical shift anisotropy (CSA) tensors in nucleotide bases is important for interpretation of NMR relaxation data in terms of local dynamic properties of nucleic acids and for analysis of residual chemical shift anisotropy (RCSA) resulting from weak alignment. CSA tensors for protonated nucleic acid base carbons have been derived from measurements on a uniformly (13)C-enriched helical A-form RNA segment and a helical B-form DNA dodecamer at natural (13)C abundance. The magnitudes of the derived CSA principal values are tightly restricted by the magnetic field dependencies of the (13)C transverse relaxation rates, whereas the tensor orientation and asymmetry follow from quantitative measurements of interference between (13)C-{(1)H} dipolar and (13)C CSA relaxation mechanisms. Changes in the chemical shift between the isotropic and aligned states, Deltadelta, complement these measurements and permit cross-validation. The CSA tensors are determined from the experimental Deltadelta values and relaxation rates, under the assumption that the CSA tensor of any specific carbon in a given type of base is independent of the base position in either the RNA or DNA helix. However, the experimental data indicate that for pyrimidine C(6) carbons in A-form RNA the CSA magnitude is considerably larger than in B-form DNA. This result is supported by quantum chemical calculations and is attributed in part to the close proximity between intranucleotide C(6)H and O(5)' atoms in RNA. The magnitudes of the measured CSA tensors, on average, agree better with previous solid-state NMR results obtained on powdered nucleosides than with prior results from quantum chemical calculations on isolated bases, which depend rather strongly on the level of theory at which the calculations are carried out. In contrast, previously computed orientations of the chemical shift tensors agree well with the present experimental results and exhibit less dependence on the level of theory at which the computations are performed.

Site-specific dynamics of strands in ss- and dsDNA as revealed by time-domain fluorescence of 2-aminopurine

It is well recognized that structure and dynamics of DNA strands guide proteins toward their cognate sites in DNA. While the dynamics is controlled primarily by the nucleotide sequence, the context of a particular sequence in relation to an open end could also play a significant role. In this work we have used the fluorescent analogue of adenine, 2-aminopurine (2-AP), to extract information on site-specific dynamics of DNA strands associated with 30-70 nucleotides length. Measurement of fluorescence lifetime and anisotropy decay kinetics in various types of DNA strands in which 2-AP was located in specific positions revealed novel insights into the dynamics of strands. We find that in single-stranded (ss) DNA, the extent of motional dynamics of the bases falls off sharply from the very end toward the middle of the strand. In contrast, the flexibility of the backbone decreases more gradually in the same direction. In double-stranded (ds) DNA, the level of base-pair fraying increases toward the ends in a graded manner. Surprisingly, the same is countered by the presence of ss-overhangs emanating from dsDNA ends. Moreover, the extent of concerted motion of bases in duplex DNA increased from the end to the middle of the duplex, a result which is both striking and counterintuitive. Most surprisingly, the two complementary strands of a duplex that were unequal in length exhibited differential dynamics: the longer one with overhangs showed a distinctly higher level of flexibility than the recessed shorter strand in the same duplex. All these results, taken together, provoke newer insights in our understanding of how different bases in DNA strands are endowed with specific dynamic properties as a function of their positions. These properties are likely to be used in facilitating specific recognitions of DNA bases by proteins during various DNA-protein interaction systems.

A physical picture of atomic motions within the Dickerson DNA dodecamer in solution derived from joint ensemble refinement against NMR and large-angle X-ray scattering data

The structure and dynamics of the Dickerson DNA dodecamer [5'd(CGCGAATTCGCG)2] in solution have been investigated by joint simulated annealing refinement against NMR and large-angle X-ray scattering data (extending from 0.25 to 3 A-1). The NMR data comprise an extensive set of hetero- and homonuclear residual dipolar coupling and 31P chemical shift anisotropy restraints in two alignment media, supplemented by NOE and 3J coupling data. The NMR and X-ray scattering data cannot be fully ascribed to a single structure representation, indicating the presence of anisotropic motions that impact the experimental observables in different ways. Refinement with ensemble sizes (Ne) of >or=2 to represent the atomic motions reconciles all the experimental data within measurement error. Cross validation against both the dipolar coupling and X-ray scattering data suggests that the optimal ensemble size required to account for the current data is 4. The resulting ensembles permit one to obtain a detailed view of the conformational space sampled by the dodecamer in solution and permit one to analyze fluctuations in helicoidal parameters, sugar puckers, and BI-BII backbone transitions and to obtain quantitative metrics of atomic motion such as generalized order parameters and thermal B factors. The calculated order parameters are in good agreement with experimental order parameters obtained from 13C relaxation measurements. Although DNA behaves as a relatively rigid rod with a persistence length of approximately 150 bp, dynamic conformational heterogeneity at the base pair level is functionally important since it readily permits optimization of intermolecular protein-DNA interactions.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

Ultrafast dynamics in DNA: "fraying" at the end of the helix

The dynamics of the electric fields in the interior of DNA are measured by using oligonucleotides in which a native base pair is replaced by a dye molecule (coumarin 102) whose emission spectrum is sensitive to the local electric field. Time-resolved measurements of the emission spectrum have been extended to a six decade time range (40 fs to 40 ns) by combining results from time-correlated photon counting, fluorescence up-conversion, and transient absorption. Recent results showed that when the reporter is placed in the center of the oligonucleotide, the dynamics are very broadly distributed over this entire time range and do not show specific time constants associated with individual processes (Andreatta, D.; et al. J. Am. Chem. Soc. 2005, 127, 7270). This paper examines an oligonucleotide with the reporter near its end. The broadly distributed relaxation seen before remains with little attenuation. In addition, a new relaxation with a well-defined relaxation time of 5 ps appears. This process is assigned to the rapid component of "fraying" at the end of the helix.

Insight into the CSA tensors of nucleobase carbons in RNA polynucleotides from solution measurements of residual CSA: towards new long-range orientational constraints

Using residual chemical shift anisotropies (RCSAs) measured in a weakly aligned stem-loop RNA, we examined the carbon chemical shift anisotropy (CSA) tensors of nucleobase adenine C2, pyrimidine C5 and C6, and purine C8. The differences between the measured RCSAs and the values back-calculated using three nucleobase carbon CSA sets [D. Stueber, D.M. Grant, 13C and 15N chemical shift tensors in adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine, J. Am. Chem. Soc. 124 (2002) 10539-10551; D. Sitkoff, D.A. Case, Theories of chemical shift anisotropies in proteins and nucleic acids, Prog. NMR Spectrosc. 32 (1998) 165-190; R. Fiala, J. Czernek, V. Sklenar, Transverse relaxation optimized triple-resonance NMR experiments for nucleic acids, J. Biomol. NMR 16 (2000) 291-302] reported previously for mononucleotides (1.4 Hz) is significantly smaller than the predicted RCSA range (-10-10 Hz) but remains larger than the RCSA measurement uncertainty (0.8 Hz). Fitting of the traceless principal CSA values to the measured RCSAs using a grid search procedure yields a cytosine C5 CSA magnitude (CSAa=(3/2.(delta11(2)+delta22(2)+delta33(2)))1/2=173+/-21 ppm), which is significantly higher than the reported mononucleotide values (131-138 ppm) and a guanine C8 CSAa (148+/-13 ppm) that is in very good agreement with the mononucleotide value reported by solid-state NMR [134 ppm, D. Stueber, D.M. Grant, 13C and (15)N chemical shift tensors in adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine, J. Am. Chem. Soc. 124 (2002) 10539-10551]. Owing to a unique sensitivity to directions normal to the base plane, the RCSAs can be translated into useful long-range orientational constraints for RNA structure determination even after allowing for substantial uncertainty in the nucleobase carbon CSA tensors.

Carbon-13 chemical shift anisotropy in DNA bases from field dependence of solution NMR relaxation rates

Knowledge of (13)C chemical shift anisotropy (CSA) in nucleotide bases is important for the interpretation of solution-state NMR relaxation data in terms of local dynamic properties of DNA and RNA. Accurate knowledge of the CSA becomes particularly important at high magnetic fields, prerequisite for adequate spectral resolution in larger oligonucleotides. Measurement of (13)C relaxation rates of protonated carbons in the bases of the so-called Dickerson dodecamer, d(CGCGAATTCGCG)(2), at 500 and 800 MHz (1)H frequency, together with the previously characterized structure and diffusion tensor yields CSA values for C5 in C, C6 in C and T, C8 in A and G, and C2 in A that are closest to values previously reported on the basis of solid-state FIREMAT NMR measurements, and mostly larger than values obtained by in vacuo DFT calculations. Owing to the noncollinearity of dipolar and CSA interactions, interpretation of the NMR relaxation rates is particularly sensitive to anisotropy of rotational diffusion, and use of isotropic diffusion models can result in considerable errors.

The structure of a mixed LNA/DNA:RNA duplex is driven by conformational coupling between LNA and deoxyribose residues as determined from 13C relaxation measurements

A study of the internal dynamics of an LNA/DNA:RNA duplex has been performed to further characterize the conformational changes associated with the incorporation of locked nucleic acid (LNA) nucleotides in a DNA:RNA duplex. In general, it was demonstrated that the LNA/DNA:RNA duplex has a very high degree of order compared to dsDNA and dsRNA duplexes. The order parameters of the aromatic carbon atoms in the LNA/DNA strand are uniformly high, whereas a sharp drop in the degree of order was seen in the RNA strand in the beginning of the AUAU stretch in the middle of the strand. This can be related to a return to normal dsRNA dynamics for the central A:U base pair. The high order of the heteroduplex is consistent with preorganization of the chimera strand for an A-form duplex conformation. These results partly explain the dramatic increase in T(m) of the chimeric heteroduplex over dsDNA and DNA:RNA hybrids of the same sequence.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

Partial 13C isotopic enrichment of nucleoside monophosphates: useful reporters for NMR structural studies

Analysis of the ¹³C isotopic labeling patterns of nucleoside monophosphates (NMPs) extracted from Escherichia coli grown in a mixture of C-1 and C-2 glucose is presented. By comparing our results to previous observations on amino acids grown in similar media, we have been able to rationalize the labeling pattern based on the well-known biochemistry of nucleotide biosynthesis. Except for a few notable absences of label (C4 in purines and C3′ in ribose) and one highly enriched site (C1′ in ribose), most carbons are randomly enriched at a low level (an average of 13%). These sparsely labeled NMPs give less complex NMR spectra than their fully isotopically labeled analogs due to the elimination of most ¹³C–¹³C scalar couplings. The spectral simplicity is particularly advantageous when working in ordered systems, as illustrated with guanosine diphosphate (GDP) bound to ADP ribosylation factor 1 (ARF1) aligned in a liquid crystalline medium. In this system, the absence of scalar couplings and additional long-range dipolar couplings significantly enhances signal to noise and resolution.

Residue specific ribose and nucleobase dynamics of the cUUCGg RNA tetraloop motif by MNMR 13C relaxation

The dynamics of the nucleobase and the ribose moieties in a 14-nt RNA cUUCGg hairpin-loop uniformly labeled with 13C and 15N were studied by 13C spin relaxation experiments. R1, R1rho and the 13C-[1H] steady-state NOE of C6 and C1' in pyrimidine and C8 and C1' in purine residues were obtained at 298 K. The relaxation data were analyzed by the model-free formalism to yield dynamic information on timescales of pico-, nano- and milli-seconds. An axially symmetric diffusion tensor with an overall rotational correlation time tau(c) of 2.31 +/- 0.13 ns and an axial ratio of 1.35 +/- 0.02 were determined. Both findings are in agreement with hydrodynamic calculations. For the nucleobase carbons, the validity of different reported 13C chemical shift anisotropy values (Stueber, D. and Grant, D. M., 2002 J. Am. Chem. Soc. 124, 10539-10551; Fiala et al., 2000 J. Biomol. NMR 16, 291-302; Sitkoff, D. and Case, D. A., 1998 Prog. NMR Spectroscopy 32, 165-190) is discussed. The resulting dynamics are in agreement with the structural features of the cUUCGg motif in that all residues are mostly rigid (0.82 < S2 < 0.96) in both the nucleobase and the ribose moiety except for the nucleobase of U7, which is protruding into solution (S2 = 0.76). In general, ribose mobility follows nucleobase dynamics, but is less pronounced. Nucleobase dynamics resulting from the analysis of 13C relaxation rates were found to be in agreement with 15N relaxation data derived dynamic information (Akke et al., 1997 RNA 3, 702-709).

Solid-phase synthesis of selectively labeled DNA: applications for multidimensional nuclear magnetic resonance spectroscopy

The solid-phase chemical synthesis method has a strong advantage over the enzymatic method for preparing selectively labeled DNA oligomers. Atom-specific and fully labeled 2'-deoxynucleosides are economically prepared with routinely available isotope precursors using this synthetic route. Special DNA oligomers prepared by advanced labeling techniques are needed for advanced NMR applications, and chemical synthesis is the method of choice to respond to such demands. As a summary of this chapter, two tables are given. Table I lists the labeled nucleosides reported to be available by chemical syntheses. Table II lists the NMR studies using labeled DNA oligomers that were prepared by chemical syntheses.

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.

Solution structure and dynamics of yeast elongin C in complex with a von Hippel-Lindau peptide

Elongin is a transcription elongation factor that stimulates the rate of elongation by suppressing transient pausing by RNA polymerase II at many sites along the DNA. It is heterotrimeric in mammals, consisting of elongins A, B and C subunits, and bears overall similarity to a class of E3 ubiquitin ligases known as SCF (Skp1-Cdc53 (cullin)-F-box) complexes. A subcomplex of elongins B and C is a target for negative regulation by the von Hippel-Lindau (VHL) tumor-suppressor protein. Elongin C from Saccharomyces cerevisiae, Elc1, exhibits high sequence similarity to mammalian elongin C. Using NMR spectroscopy we have determined the three-dimensional structure of Elc1 in complex with a human VHL peptide, VHL(157-171), representing the major Elc1 binding site. The bound VHL peptide is entirely helical. Elc1 utilizes two C-terminal helices and an intervening loop to form a binding groove that fits VHL(157-171). Chemical shift perturbation and dynamics analyses reveal that a global conformational change accompanies Elc1/VHL(157-171) complex formation. Moreover, the disappearance of conformational exchange phenomena on the microsecond to millisecond time scale within Elc1 upon VHL peptide binding suggests a role for slow internal motions in ligand recognition.

NMR evidence for mechanical coupling of phosphate B(I)-B(II) transitions with deoxyribose conformational exchange in DNA

The conformational exchange of the phosphate and deoxyribose groups of the DNA oligomers d(GCGTACGC)(2) and d(CGCTAGCG)(2) have been investigated using a combination of homonuclear and heteronuclear NMR techniques. Two-state exchange between phosphate B(I) and B(II) conformations and deoxyribose N and S conformations was expressed as percent population of the major conformer, %B(I) or %S. Sequence context-dependent variations in %B(I) and %S were observed. The positions of the phosphate and deoxyribose equilibria provide a quantitative measure of the ps to ns timescale dynamic exchange processes in the DNA backbone. Linear correlations between %B(I), %S, and previously calculated model free (13)C order parameters (S(2)) were observed. The %B(I) of the phosphates were found to be correlated to the S(2) of the flanking C3' and C4' atoms. The %B(I) was also found to be correlated with the %S and C1' S(2) of the deoxyribose ring 5' of the phosphates. The %B(I) of opposing phosphates is correlated, while the %B(I) of sequential phosphates is anti-correlated. These correlations suggest that conformational exchange processes in DNA are coupled to each other and are modulated by DNA base sequence, which may have important implications for DNA-protein interactions.

Slow motion in the CAA*TTG sequence of a DNA decamer duplex studied by NMR

In DNA duplexes, pyrimidine-purine steps are believed to be flexible or conformationally unstable. Indeed, several DNA crystal structures exhibit a multitude of conformations for CpA*TpG steps. The question arises of whether this structural flexibility is accompanied by dynamical flexibility, i.e., a question pertaining to the energy barrier between conformations. Except for TpA steps, slow motions on the microsecond-to-millisecond time scale have not been detected in duplexes until now. In the present study, such slow motion was investigated by 1H, 13C, and 15N NMR relaxation measurements on a DNA decamer d(CATTTGCATC)*d(GATGCAAATG). The DNA decamer was enriched with 15% 13C and 98% 15N isotopes for each adenosine and guanosine residue. Three lines of evidence support the notion of slow motion in the CAA*TTG moiety. Analysis of (15)N relaxation showed that the order parameter, S2, of guanosine imino NH groups was about 0.8, similar to that of CH groups for this oligomer. The strong temperature dependence of guanosine NH S2 in the CAA*TTG sequence indicated the presence of a large-amplitude motion. Signals of adenosine H8 protons in the CAA*TTG sequence were broadened in 2D 1H NOESY spectra, which also suggested the existence of slow motion. As well as being smaller than for other adenine residues, the 1H T2 values exhibited a magnetic field strength dependence for all adenosine H8 signals in the ATTTG*CAAAT region, suggesting slow motions more pronounced at the first adenosine in the CAA*TTG sequence but extending over the CAAAT*ATTTG region. This phenomenon was further examined by the pulse field strength dependence of the 1H, 13C, and 15N T1rho values. 1H and 13C T1rho values showed a pulse field strength dependence, but 15N T1rho did not. Assuming a two-site exchange process, an exchange time constant of 20-300 micros was estimated for the first adenosine in the CAA sequence. The exact nature of this motion remains unknown.

A nicked duplex decamer DNA with a PEG(6) tether

A dumbbell double-stranded DNA decamer tethered with a hexaethylene glycol linker moiety (DDSDPEG), with a nick in the centre of one strand, has been synthesised. The standard NMR methods, E.COSY, TOCSY, NOESY and HMQC, were used to measure (1)H, (31)P and T:(1) spectral parameters. Molecular modelling using rMD-simulated annealing was used to compute the structure. Scalar couplings and dipolar contacts show that the molecule adopts a right-handed B-DNA helix in 38 mM phosphate buffer at pH 7. Its high melting temperature confirms the good base stacking and stability of the duplex. This is partly attributed to the presence of the PEG(6) linker at both ends of the duplex that restricts the dynamics of the stem pentamers and thus stabilises the oligonucleotide. The inspection of the global parameters shows that the linker does not distort the B-DNA geometry. The computed structure suggests that the presence of the nick is not disturbing the overall tertiary structure, base pair geometry or duplex base pairing to a substantial extent. The nick has, however, a noticeable impact on the local geometry at the nick site, indicated clearly by NMR analysis and reflected in the conformational parameters of the computed structure. The (1)H spectra also show much sharper resonances in the presence of K(+) indicating that conformational heterogeneity of DDSDPEG is reduced in the presence of potassium as compared to sodium or caesium ions. At the same time the (1)H resonances have longer T:(1) times. This parameter is suggested as a sensitive gauge of stabilisation.

Active site dynamics in the lead-dependent ribozyme

Conformational dynamics are an important property of ribozymes and other RNA molecules but there is currently only limited information on the relationship between dynamics and RNA function. A recent structural study of the lead-dependent ribozyme, known as the leadzyme, showed significant dynamics at the active site and indicated that a structural rearrangement is required for the reaction to proceed from the ground to the transition state. In this work, microsecond-to-millisecond dynamics of the leadzyme are probed by analysis of the power dependence of (13)C NMR relaxation times in the rotating frame (T(1)(rho)). These results revealed a wide range of conformational dynamics for various residues in the leadzyme. For residue A25 in the active site, the power dependence of T(1)(rho) yielded an exchange lifetime similar to that previously measured by line-shape analysis, and provides an important calibration of this T(1)(rho) methodology for probing the dynamics of macromolecules. Strong evidence was also found for a previously suggested dynamic network of hydrogen bonds stabilizing the GAAA tetraloop motif. Within the active site of the leadzyme, internal motions are observed on a wide variety of time scales, suggesting a complex landscape of accessible states, and potential correlations between observed motions and catalytic function are discussed. These results demonstrate that the power dependence of (13)C T(1)(rho) relaxation times provides a valuable method for probing dynamics in nucleic acids.