Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations

NMR structure of the full-length linear dimer of stem-loop-1 RNA in the HIV-1 dimer initiation site

The packaging signal of HIV-1 RNA contains a stem-loop structure, SL1, which serves as the dimerization initiation site for two identical copies of the genome and is important for packaging of the RNA genome into the budding virion and for overall infectivity. SL1 spontaneously dimerizes via a palindromic hexanucleotide sequence in its apical loop, forming a metastable kissing dimer form. Incubation with nucleocapsid protein causes this form to refold to a thermodynamically stable mature linear dimer. Here, we present an NMR structure of the latter form of the full-length SL1 sequence of the Lai HIV-1 isolate. The structure was refined using nuclear Overhauser effect and residual dipolar coupling data. The structure presents a symmetric homodimer of two RNA strands of 35 nucleotides each; it includes five stems separated by four internal loops. The central palindromic stem is surrounded by two symmetric adenine-rich 1-2 internal loops, A-bulges. All three adenines in each A-bulge are stacked inside the helix, consistent with the solution structures of shorter SL1 constructs determined previously. The outer 4-base pair stems and, proximal to them, purine-rich 1-3 internal loops, or G-bulges, are the least stable parts of the molecule. The G-bulges display high conformational variability in the refined ensemble of structures, despite the availability of many structural restraints for this region. Nevertheless, most conformations share a similar structural motif: a guanine and an adenine from opposite strands form a GA mismatch stacked on the top of the neighboring stem. The two remaining guanines are exposed, one in the minor groove and another in the major groove side of the helix, consistent with secondary structure probing data for SL1. These guanines may be recognized by the nucleocapsid protein, which binds tightly to the G-bulge in vitro.

Changes in DNA bending induced by restricting nucleotide ring pucker studied by weak alignment NMR spectroscopy

Changes in bending of the DNA helix axis caused by the introduction of conformationally locked nucleotide analogs into the center region of the palindromic Dickerson dodecamer, d(CGCGAATTCGCG)(2), have been studied by NMR measurement of residual one-bond (13)C-(1)H dipolar couplings. Thymidine analogs, in which the deoxyribose was substituted by bicyclo[3.1.0]hexane, were incorporated in the T7, T8, and T7T8 positions. These nucleotide analogs restrict the ring pucker to the C2'-exo or "north" conformation, instead of C2'-endo or "south," which dominates in regular B-form DNA. For all three oligomers, bending toward the major groove is found relative to the native molecule. The effects are additive with bending of 5 +/- 1 degrees per locked nucleotide. Measurement of the change in bending is more accurate than measurement of the bending angle itself and requires far fewer experimental data.

Structure of DNA Hexamer Sequence d-CGATCG by Two-dimensional Nuclear Magnetic Resonance Spectroscopy and Restrained Molecular Dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.

Force field dependence of NMR-Based, restrained molecular dynamics DNA structure calculations including an analysis of the influence of residual dipolar coupling restraints

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.

Incorporating residual dipolar couplings into the NMR solution structure determination of nucleic acids

NMR solution structures of nucleic acids are generally less well defined than similar-sized proteins. Most NMR structures of nucleic acids are defined only by short-range interactions, such as intrabase-pair or sequential nuclear Overhauser effects (NOEs), and J-coupling constants, and there are no long-range structural data on the tertiary structure. Residual dipolar couplings represent an extremely valuable source of distance and angle information for macromolecules but they average to zero in isotropic solutions. With the recent advent of general methods for partial alignment of macromolecules in solution, residual dipolar couplings are rapidly becoming indispensable constraints for solution NMR structural studies. These residual dipolar couplings give long-range global structural information and thus complement the strictly local structural data obtained from standard NOE and torsion angle constraints. Such global structural data are especially important in nucleic acids due to the more elongated, less-globular structure of many DNAs and RNAs. Here we review recent progress in application of residual dipolar couplings to structural studies of nucleic acids. We also present results showing how refinement procedures affect the final solution structures of nucleic acids.

Relationship of DNA structure to internal dynamics: correlation of helical parameters from NOE-based NMR solution structures of d(GCGTACGC)(2) and d(CGCTAGCG)(2) with (13)C order parameters implies conformational coupling in dinucleotide units

The coupling between the conformational properties of double-stranded DNA and its internal dynamics has been examined. The solution structures of the isomeric DNA oligomers d(GCGTACGC)(2) (UM) and d(CGCTAGCG)(2) (CTSYM) were determined with (1)H NMR spectroscopy by utilizing distance restraints from total relaxation matrix analysis of NOESY cross-peak intensities in restrained molecular dynamics calculations. The root-mean-square deviation of the coordinates for the ensemble of structures was 0.13 A for UM and 0.49 A for CTSYM, with crystallographic equivalent R(c)=0.41 and 0.39 and sixth-root residual R(x)=0.11 and 0.10 for UM and CTSYM, respectively. Both UM and CTSYM are B-form with straight helical axes and show sequence-dependent variations in conformation. The internal dynamics of UM and CTSYM were previously determined by analysis of (13)C relaxation parameters in the context of the Lipari & Szabo model-free formalism. Helical parameters for the two DNA oligomers were examined for linear correlations with the order parameters (S(2)) of groups of (13)C spins in base-pairs and dinucleotide units of UM and CTSYM. Correlations were found for six interstrand base-pair parameters tip, y-displacement, inclination, buckle and stretch with various combinations of S(2) for atoms in Watson-Crick base-pairs and for two inter-base-pair parameters, rise and roll with various combinations of S(2) for atoms in dinucleotides. The correlations for the interstrand base-pair helical parameters indicate that the conformations of the deoxyribose residues of each strand are dynamically coupled. Also, the inter-base-pair separation has a profound effect on the local internal motions available to the DNA, supporting the idea that rise is a principal degree of freedom for DNA conformational variability. The correlations indicate collective atomic motions of spins that may represent specific motional modes in DNA, and that base sequence has a predictable effect on the relative order of groups of spins both in the bases and in the deoxyribose ring of the DNA backbone. These observations suggest that an important functional outcome of DNA base sequence is the modulation of both the conformation and dynamic behavior of the DNA backbone.

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.

DNA bending and sequence-dependent backbone conformation NMR and computer experiments

Although DNA bending plays a crucial role in several biological processes, very little is known experimentally about the relationship between sugar phosphate conformation and sequence directed bending. In this paper, we determine the coupling constants for a nonself-complementary 11-mer A-tract DNA duplex from 2D NMR experiments and along each chain of the duplex, we report the sugar pucker, torsional preferences and conformational averaging about the C3'-O3', C4'-C5' and C5'-O5' bonds for each nucleotide. The A-tract exists as an equilibrium blend of canonical B-form and noncanonical B-form in which the exocyclic C4'-C5' bond is in trans conformation as in the original Watson-Crick model [Crick, F.H.C. & Watson, J.D. (1954) Proc. Roy. Soc. (London), A223, 80-96]. The trans conformation at the C4'-C5' can increase the interphosphate distance and lead to local unwinding of the duplex and rolling of the base pair into the major groove. This will create a kink or hinge. At the 3'-end of the A-tract in the purine-thymine step, the duplex is compressed by the presence of a junction between A and B forms of DNA exclusively in one strand, with consequent reduction of the phosphate-phosphate distance. The coupling constant data seriously disagree with the A-tract DNA bending model of Crothers [Koo, H.-S., Wu, H.-M. & Crothers, D.M. (1986) Nature 320, 501-506], but is in agreement with the finding of Leroy et al. [Leroy, J.-L., Charretier, E., Kochoyan, M. & Gueron, M. (1988) Biochemistry 27, 8894-8898] that the structure that drives bending in the A-tract is locally different from B-DNA. Structural distortions are extremely localized with little or no propagation. It is likely that transcription factor proteins recognize these preexisting deformations in the free DNA itself and mold it into the matrix of the protein.

GC rich DNA oligonucleotides with narrow minor groove width

Investigation of the width of the minor groove using 500 MHz NMR spectroscopy in three closely related 11-mer B-DNA duplexes shows that the minor groove is narrow in a GC rich oligonucleotide, and that a narrow minor groove is not something endemic to DNAs with persistent repetitions of adenine nucleotides (A-tract DNA). The width of the groove is dictated by local sequence contexts and independent of neighboring A-tract DNA.

Structure of an anti-HIV-1 hammerhead ribozyme complex with a 17-mer DNA substrate analog of HIV-1 gag RNA and a mechanism for the cleavage reaction: 750 MHz NMR and computer experiments

The structure of an anti-HIV-1 ribozyme-DNA abortive substrate complex was investigated by 750 MHz NMR and computer modeling experiments. The ribozyme was a chimeric molecule with 30 residues-18 DNA nucleotides, and 12 RNA residues in the conserved core. The DNA substrate analog had 17 residues. The chimeric ribozyme and the DNA substrate formed a shortened ribozyme-abortive substrate complex of 47 nucleotides with two DNA stems (stems I and III) and a loop consisting of the conserved core residues. Circular dichroism spectra showed that the DNA stems assume A-family conformation at the NMR concentration and a temperature of 15 degrees C, contrary to the conventional wisdom that DNA duplexes in aqueous solution populate entirely in the B-form. It is proposed that the A-family RNA residues at the core expand the A-family initiated at the core into the DNA stems because of the large free energy requirement for the formation of A/B junctions. Assignments of the base H8/H6 protons and H1' of the 47 residues were made by a NOESY walk. In addition to the methyl groups of all T's, the imino resonances of stems I and III and AH2's were assigned from appropriate NOESY walks. The extracted NMR data along with available crystallographic data, were used to derive a structural model of the complex. Stems I and III of the final model displayed a remarkable similarity to the A form of DNA; in stem III, a GC base pair was found to be moving into the floor of the minor groove defined by flanking AT pairs; data suggest the formation of a buckled rhombic structure with the adjacent pair; in addition, the base pair at the interface of stem III and the loop region displayed deformed geometry. The loop with the catalytic core, and the immediate region of the stems displayed conformational multiplicity within the NMR time scale. A catalytic mechanism for ribozyme action based on the derived structure, and consistent with biochemical data in the literature, is proposed. The complex between the anti HIV-1 gag ribozyme and its abortive DNA substrate manifests in the detection of a continuous track of A.T base pairs; this suggests that the interaction between the ribozyme and its DNA substrate is stronger than the one observed in the case of the free ribozyme where the bases in stem I and stem III regions interact strongly with the ribozyme core region (Sarma, R. H., et al. FEBS Letters 375, 317-23, 1995). The complex formation provides certain guidelines in the design of suitable therapeutic ribozymes. If the residues in the ribozyme stem regions interact with the conserved core, it may either prevent or interfere with the formation of a catalytically active tertiary structure.

Sensitivity of NMR internucleotide distances to B-DNA conformation: underlying mechanics

Nuclear magnetic resonance (NMR) spectroscopy, combining correlated spectroscopy (COSY) coupling constant measurements with nuclear Overhauser effect spectroscopy (NOESY) interatomic distances, should make it possible to determine an averaged solution structure for DNA oligomers. However, even if such data could be obtained with high accuracy, it is not clear which structural parameters of DNA would be determined. Here, the relationships between measurable internucleotide distances and helical parameters are systematically studied through molecular modelling. Investigations are carried out using four representative sequences, (ACGT)n, (TCGA)n, (AGCT)n and (TGCA)n, composed of repeated tetranucleotides belonging to oligomers previously studied by NMR. Correlations between interatomic distances become evident and strong connections between distances and inter-base helical parameters are observed. Results imply that twist, roll, shift and slide values can be accurately determined from NMR data. Sequence independent mechanical coupling which link backbone and sugar conformations to helical twist are also described.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

Sequence dependent effects of CpG cytosine methylation. A joint 1H-NMR and 31P-NMR study

The impact of cytosine methylation in the central CpG step of two closely related octanucleotide duplexes d(CATCGATG)2 and d(CTTCGAAG)2 was examined by 1H-NMR and 31P-NMR experiments, and a quantitative structural analysis was performed using the NOE-derived distances, the sugar puckers and the epsilon torsion angles. The two starting oligonucleotides displayed a B-DNA conformation with, however, significant local structure differences. Although the methylated oligonucleotides retained their B-DNA conformation, different structural and thermal stability effects were observed. The magnitude of the methylation effects was to depend on the initial conformation of the CpG site, which is governed by the nature of the dinucleotide AT or TT located on the CpG flanks. As an example of sequence dependence, the methylation of CpG entailed larger conformational variation in d(CATCGATG)2 than in d(CTTCGAAG)2. In this study, the 1H and 31P chemical-shift parameters averred as extremely sensitive probes for detecting subtle conformational changes. Finally, our comparative results may aid our understanding of the structural and related biological effects produced by cytosine methylation in DNA.

Secondary structure in solution of two anti-HIV-1 hammerhead ribozymes as investigated by two-dimensional 1H 500 MHz NMR spectroscopy in water

Two hammerhead chimeric RNA/DNA ribozymes (HRz) were synthesized in pure form. Both were 30 nucleotides long, and the sequences were such that they could be targeted to cleave the HIV-1 gag RNA. Named HRz-W and HRz-M, the former had its invariable core region conserved, the latter had a uridine in the invariable region replaced by a guanine. Their secodary structures were determined by 2D NOESY 1H 500 MHz NMR spectroscopy in 90% water and 10% D2(0), following the imino protons. The data show that both HRz-M and HRz-W form identical secondary structures with stem regions consisting of continuous stacks of AT and GT pairs. An energy minimized computer model of this stem region is provided. The results suggest that the loss of catalytic activity that is known to result when an invariant core residue is replaced is not related to the secondary structure of the ribozymes in the absence of substrate.

Probability assessment of conformational ensembles: sugar repuckering in a DNA duplex in solution

Conformational flexibility of molecules in solution implies that different conformers contribute to the NMR signal. This may lead to internal inconsistencies in the 2D NOE-derived interproton distance restraints and to conflict with scalar coupling-based torsion angle restraints. Such inconsistencies have been revealed and analyzed for the DNA octamer GTATAATG.CATATTAC, containing the Pribnow box consensus sequence. A number of subsets of distance restraints were constructed and used in the restrained Monte Carlo refinement of different double-helical conformers. The probabilities of conformers were then calculated by a quadratic programming algorithm, minimizing a relaxation rate-base residual index. The calculated distribution of conformers agrees with the experimental NOE data as an ensemble better than any single structure. A comparison with the results of this procedure, which we term PARSE (Probability Assessment via Relaxation rates of a Structural Ensemble), to an alternative method to generate solution ensembles showed, however, that the detailed multi-conformational description of solution DNA structure remains ambiguous at this stage. Nevertheless, some ensemble properties can be deduced with confidence, the most prominent being a distribution of sugar puckers with minor populations in the N-region and major populations in the S-region. Importantly, such a distribution is in accord with the analysis of independent experimental data--deoxyribose proton-proton scalar coupling constants.

Solution structure of the octamer motif in immunoglobulin genes via restrained molecular dynamics calculations

The solution structure of the DNA decamer d(CATTTGCATC)-d(GATGCAAATG), comprising the octamer motif of immunoglobulin genes, is determined by restrained molecular dynamics (rMD) simulations. The restraint data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were previously obtained by a complete relaxation matrix analysis of the two-dimensional nuclear Overhauser enhancement (2D NOE) intensities and by the quantitative simulation of cross-peaks in double-quantum-filtered correlated (2QF-COSY) spectra. The influence of torsion angles and the number of experimental distance restraints on the structural refinement has been systematically examined. Omitting part of the experimental NOE-derived distances results in reduced restraint violations and lower R factors but impairs structural convergence in the rMD refinement. Eight separate restrained molecular dynamics simulations were carried out for 20 ps each, starting from either energy-minimized A- or B-DNA. Mutual atomic root-mean-square (rms) differences among the refined structures are well below 1 A and comparable to the rms fluctuations of the atoms about their average position, indicating convergence to essentially identical structures. The average refined structure was subjected to an additional 100 ps of rMD simulations and analyzed in terms of average torsion angles and helical parameters. The B-type duplex exhibits clear sequence-dependent variations in its geometry with a narrow minor groove at the T3.A3 tract and a large positive roll at the subsequent TG.CA step. This is accompanied by a noticeable bend of the global helix axis into the major groove. There is also evidence of significant flexibility of the sugar-phosphate backbone with rapid interconversion among different conformers.

Base dependence of B-DNA sugar conformation in solution and in the solid state

Analysis of 1H-NOESY solution data for eight short DNA duplexes has revealed pronounced differences between the sugar conformations of purine and pyrimidine nucleotides. It was found that the H1'-H4' interproton distance is less than ca. 3.0 A in pyrimidine sugars, while in purine sugars it is more than ca. 3.0A. This difference has been analyzed by comparison with the sugar conformations of highly resolved B-DNA crystal structures and model sugar conformations. The conclusion can be drawn that the deoxyribose conformation is of the general C2'-endo type but pyrimidine sugars are characterized by smaller phase angles of pseudorotation P (90 degrees < P < 150 degrees), while purine sugars have larger P values that are greater than ca. 140 degrees (140 degrees < P < 180 degrees). There is no such clear base dependence of sugar conformation in highly resolved B-DNA crystal structures; however the similar trend can be seen as in the solution studies. Based on B-type DNA crystal structures, J-coupling constants have been calculated, and the applicability of experimental coupling measurements to the determination of sugar conformation is discussed.

Investigation of solution structure of d(GAATTTAAATTC)2 by 1H NMR, molecular dynamics, mechanics, refinement by back-calculation of the NOESY spectrum and analysis of this structure using X-ray data

1H NMR spectroscopy, restrained molecular mechanics and dynamics and refinements after back-calculation of the NOESY spectrum have been performed to study the structure of the d(GAATTTAAATTC)2 duplex and to determine whether it is bent or not. It is found that the duplex adopts a B-type conformation; all sugar conformations belong to the C2' endo region and purines have a larger pseudorotation angle as compared to pyrimidines. The cross-strand AH2(n)-AH1' (m + 1) distance (where (n) and (m) are complementary residues), crucial for an anomalous A/T tract structure, is large on the TA step and gradually decreases at the 3' and 5' ends of the TTTAAA tract and follows the rules proposed previously [Chuprina, V.P., Lipanov, A.A., Fedoroff, O.Yu, Kim, S.G., Kintanar, A., and Reid, B.R., Proc. Natl. Acad. Sci. U.S.A. 88, 9087 (1991)]. The changes in this distance correlate with those in the T1 value for AH2 protons which we measured for several oligonucleotide sequences. A total number of about 250 interproton distance constraints were determined from NOESY spectra and were used for structure determination by molecular mechanics, dynamics and refinement by back-calculations. It is shown that these data are not enough to determine whether the duplex is bent or not. The whole family of B-type conformations including bent and straight structures fit well with the available NMR data. In principle, additional non-NMR data could be used in order to reduce the number of the allowable structures. The refinement of the structures with additional different non-NMR constraints (used as a driving force) on P-P or H1'-H1' minor groove width distances in the TA region shows a very good correlation between these distances and the angle of bending of the dodecamer. The more the minor groove width increases in the TA region the more the duplex is bent at the major groove of this region. On the other hand, there is also a very good correlation between P-P, H1'-H1' and AH2-H1' cross-strand distances as follows from analysis of X-ray B-type structures. These two correlations, together with the increased AH2-H1' cross-strand NMR distance in the TA region of the dodecamer indicate that the duplex should be characterized by a wider minor groove in the TA region and be bent in the major groove in this region.