How accurately can oligonucleotide structures be determined from the hybrid relaxation rate matrix/NOESY distance restrained molecular dynamics approach?

The accuracy and precision of structures derived from a combined hybrid relaxation rate matrix/NOESY distance restrained molecular dynamics methodology were examined with simulations that included typical experimental errors. NOESY data were simulated for a DNA dodecamer duplex, d-(CGCGAATTCGCG)2, with added volume error of approximately 20% and low-level thermal noise. Distances derived from a hybrid relaxation matrix analysis of the NOE data were used as constraints in molecular dynamics driven structural refinements of several initial model geometries. The final structures were compared against results obtained from the traditional isolated two-spin approximation treatment of these NOESY volumes and also against refined structures that employed error-free data. Results show that the structures derived from the relaxation rate matrix analysis of the NOESY data are more accurate than those derived from a simple two-spin approximation analysis and it is possible to achieve refinement to the level of simulated experimental error. Results may be significantly improved with the use of either more accurately measured NOESY volumes or additional matrix-derived constraints. Many of the helical parameters and backbone torsional angles may be accurately reproduced by the hybrid matrix methodology.     

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of a covalent CPI-CDPI2-oligodeoxyribonucleotide decamer complex

CPI-CDPI2 is a synthetic analogue of CC-1065, which is a naturally occurring antitumor antibiotic. Assignment of the 1H NMR spectra of a CPI-CDPI2-oligodeoxyribonucleotide decamer, d-(CGCTTAAGCG)2, complex has been made by two-dimensional 1H/1H spectroscopy. The solution structure of the complex was calculated by an iterative hybrid relaxation matrix method combined with NOESY distance restrained molecular dynamics. Refinement proceeded in two steps in which the decamer was initially refined alone and then CPI-CDPI2 was added to the structure to allow initial estimates of drug-DNA contacts. A hybrid matrix/MD refinement was used to better take into account problems associated with spin diffusion. Thus the distances from the 2D NOESY spectra were calculated from the relaxation rate matrix which were evaluated from a hybrid NOESY volume matrix comprising elements from the experimental spectrum and those calculated from an initial structure. The hybrid matrix derived distances were then used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure was then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. The efficacy of CC-1065 has been attributed to its minor groove binding and alkylation to the N3 position of adenosine. CPI-CDPI2 appears to bind to the decamer in a similar manner. The effect of CPI-CDPI2 on the decamer's 1H and 31P spectrum was consistent with a minor groove binding motif with the drug alkylating at A17 with the CDPI rings oriented toward the 5'-end of the alkylated strand. In addition, the NMR data support one major adduct but also indicate the presence of a minor adduct. The latter could represent a drug alkylation of the DNA at a secondary site (or alternative orientation of the rings).     

NMR structural refinement of a tandem G.A mismatched decamer d(CCAAGATTGG)2 via the hybrid matrix procedure

A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations is used to evaluate the NMR solution structure of a tandemly positioned G.A double mismatch decamer oligodeoxyribonucleotide duplex, d(CCAAGATTGG)2. An iterative refinement method using a hybrid relaxation matrix combined with restrained molecular dynamics calculations is shown to provide structures having good agreement with the experimentally derived structures. Distances incorporated into the MD simulations have been calculated from the relaxation rate matrix evaluated from a hybrid NOESY volume matrix whose elements are obtained from the merging of experimental and calculated NOESY intensities. Starting from both A- and B-DNA and mismatch syn and anti models, it is possible to calculate structures that are in good atomic RMS agreement with each other (less than 1.6 A RMS) but differ from the reported crystal structure (greater than 3.6 A). Importantly, the hybrid matrix derived structures are in excellent agreement with the experimental solution conformation as determined by comparison of the 200-ms simulated and experimental NOESY spectra, while the crystallographic data provide spectra that are grossly different.     

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex

Assignment of the 1H and 31P NMR spectra of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. The downfield 31P resonance previously noted by Patel et al. (1982) has been assigned by both 17O labeling of the phosphate as well as a pure absorption phase constant-time heteronuclear 31P-1H correlated spectrum and has been associated with the phosphate on the 3' side of the extrahelical adenosine. JH3'-P coupling constants for each of the phosphates of the tridecamer were obtained from the 1H-31P J-resolved selective proton-flip 2D spectrum. By use of a modified Karplus relationship the C4-C3'-O3-P torsional angles (epsilon) were obtained. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. The 31P chemical shifts and epsilon torsional angles follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. Because the extrahelical adenosine significantly distorts the deoxyribose phosphate backbone conformation even several bases distant from the extrahelical adenosine, 31P chemical shifts show complex site- and sequence-specific variations. Modeling and NOESY distance-restrained energy minimization and restrained molecular dynamics suggest that the extrahelical adenosine stacks into the duplex. However, a minor conformation is also observed in the 1H NMR, which could be associated with a structure in which the extrahelical adenosine loops out into solution.     

A Renaissance Man: In Memoriam of Jon Widom (1955–2011)

Abstract

Assignment of the ¹H and ³¹P resonances of a decamer DNA duplex, d(CGCTTAAGCG)₂ was determined by two-dimensional COSY, NOESY and ¹H- ³¹P Pure Absorption phase Constant time (PAC) heteronuclear correlation spectroscopy. The solution structure of the decamer was calculated by an iterative hybrid relaxation matrix method combined with NOESY-distance restrained molecular dynamics. The distances from the 2D NOESY spectra were calculated from the relaxation rate matrix which were evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix-derived distances were then used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure was then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. J_H3′-P coupling constants for each of the phosphates of the decamer were obtained from ¹H-³¹P J-resolved selective proton flip 2D spectra. By using a modified Karplus relationship the C4′-C3′-03′-P torsional angles (?) were obtained. Comparison of the ³¹P chemical shifts and J_H3′-P coupling constants of this sequence has allowed a greater insight into the various factors responsible for ³¹P chemical shift variations in oligonucleotides. It also provides an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These correlations are consistent with the hypothesis that changes in local helical structure perturb the deoxyribose phosphate backbone. The variation of the ³¹P chemical shift, and the degree of this variation from one base step to the next is proposed as a potential probe of local helical conformation within the DNA double helix. The pattern of calculated ? and ζ torsional angles from the restrained molecular dynamics refinement agrees quite well with the measured J_H3′-P coupling constants. Thus, the local helical parameters determine the length of the phosphodiester backbone which in turn constrains the phosphate in various allowed conformations.     

Extrahelical adenosine stacks into right-handed DNA: solution conformation of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

This paper reports on features of the three-dimensional structure of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) self-complementary duplex (designated adenosine 13-mer), which contains symmetrical extrahelical adenosines in the interior of the helix. The majority of the protons have been assigned from two-dimensional nuclear Overhauser effect (NOESY) spectra of the adenosine 13-mer in H2O and D2O solution. The measurement of NOESY cross-peak volume integrals as a function of mixing time has yielded a set of 96 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds, which have served as input parameters for a distance geometry analysis of one symmetric half of the adenosine 13-mer duplex. We demonstrate that the extrahelical adenosine stacks into the duplex for all refined structures without disruption of base pairing on either side of the modification site. The distance geometry refinement yields two classes of conformations consistent with distance measurements but which differ in orientation of the stacked extrahelical adenosine at the modification site.     

2'-deoxyribonolactone lesion in DNA: refined solution structure determined by nuclear magnetic resonance and molecular modeling

The solution conformation of the DNA duplex d(C1G2C3A4C5L6C7A8C9G10C11).d(G12C13G14T15G16T17G18T19G20C21G22 ) containing the 2'-deoxyribonolactone lesion (L6) in the middle of the sequence has been investigated by NMR spectroscopy and restrained molecular dynamics calculations. Interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peak intensities. These distances, along with torsion angles for sugar rings and additional data derived from canonical A- and B-DNA, have been used for structure refinement by restrained molecular dynamics (rMD). Six rMD simulations have been carried out starting from both regular A- and B-DNA forms. The pairwise rms deviations calculated for each refined structure are <1 A, indicating convergence to essentially the same geometry. The accuracy of the rMD structures has been assessed by complete relaxation matrix back-calculation. The average sixth-root residual index (Rx = 0.052 +/- 0.003) indicated that a good fit between experimental and calculated NOESY spectra has been achieved. Detailed analysis revealed a right-handed DNA conformation for the duplex in which both the T17 nucleotide opposite the abasic site and the lactone ring are located inside the helix. No kinking is observed for this molecule, even at the abasic site step. This structure is compared to that of the oligonucleotide with the identical sequence containing the stable tetrahydrofuran abasic site analogue that we reported previously [Coppel, Y., Berthet, N., Coulombeau, C., Coulombeau, Ce., Garcia, J., and Lhomme, J. (1997) Biochemistry 36, 4817-4830].     

Experimental and theoretical studies of the conformational perturbations induced by an abasic site

The three-dimensional structure of the natural undecamer duplex d(CGCACACACGC). d(GCGTGTGTGCG) has been determined by the combined use of NMR spectroscopy and restrained molecular dynamics (rMD) and also by molecular mechanics calculations using the JUMNA program without experimental distance constraints. Both procedures have also been used to model the abasic structure d(CGCACOCACGC).d(GCGTGTGTGCG), where 'O' indicates a modified abasic site: 3-hydroxy-2-(hydroxymethyl) tetrahydrofuran. For the natural duplex, 134 interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peaks intensities, using MARDIGRAS software. These distances along with 100 torsion angles for sugar ring and additional data derived from canonical A and B-DNA, have been used for structures refinement by restrained molecular dynamics. Comparison of the natural oligomer with the abasic structure obtained earlier by NMR/rMD (Y. Coppel, N. Berthet, C. Coulombeau, Ce. Coulombeau, J. Garcia and J. Lhomme, Biochemistry 36, 4817-4830, 1997) confirms that the creation of an abasic site, in this sequence context, leads to marked helix kinking. It is also shown that the JUMNA procedure is capable of reproducing the overall structural features of the natural and damaged DNA conformations without the use of experimental constraints.     

An abasic site in DNA. Solution conformation determined by proton NMR and molecular mechanics calculations

We have determined the three-dimensional structure of a non-selfcomplementary nonanucleotide duplex which contains an abasic (apyrimidinic) site in the centre, i.e. a deoxyribose residue opposite an adenosine. The majority of the base and sugar proton resonances were assigned by NOESY, COSY and 2DQF spectra in D2O and H2O. We have measured the initial slope of buildup of NOEs in NOESY spectra at very short mixing times (25 to 50 ms), and from these were able to establish interproton distances for the central part of the duplex. We propose a different strategy for proton-proton distance determinations which takes into account the observed variations in correlation times for particular proton-proton vectors. A set of 31 measured interproton distances was incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. Two structures were obtained which retain all aspects of a classical B DNA in which the unpaired adenine and the abasic deoxyribose lie inside the helix. We observe that the non-hydrogen bonded adenine is held well in the helix, the Tm of this base being the same as that of the A.T base pairs in the same duplex.     

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

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

Solution structures of human transforming growth factor alpha derived from 1H NMR data

The 600-MHz 1H NMR spectrum of the des-Val-Val mutant of human transforming growth factor alpha (TGF-alpha) was reassigned at pH = 6.3. The conformation space of des-Val-Val TGF-alpha was explored by distance geometry embedding followed by restrained molecular dynamics refinement using NOE distance constraints and some torsion angle constraints derived from J-couplings. Over 80 long-range NOE constraints were found by completely assigning all resolved cross-peaks in the NOESY spectra. Low NOE constraint violations were observed in structures obtained with the following three different refinement procedures: interactive annealing in DSPACE, AMBER 3.0 restrained molecular dynamics, and dynamic simulated annealing in XPLOR. The segment from Phe15 to Asp47 was found to be conformationally well-defined. Back-calculations of NOESY spectra were used to evaluate the quality of the structures. Our calculated structures resemble the ribbon diagram presentations that were recently reported by other groups. Several side-chain conformations appear to be well-defined as does the relative orientation of the C loop to the N-terminal half of the protein.     

Solution structure of a highly stable DNA duplex conjugated to a minor groove binder.

The tripeptide 1,2-dihydro-(3 H )-pyrrolo[3,2- e ]indole-7-carboxylate (CDPI3) binds to the minor groove of DNA with high affinity. When this minor groove binder is conjugated to the 5'-end of short oligonucleotides the conjugates form unusually stable hybrids with complementary DNA and thus may have useful diagnostic and/or therapeutic applications. In order to gain an understanding of the structural interactions between the CDPI3minor groove binding moiety and the DNA, we have determined and compared the solution structure of a duplex consisting of oligodeoxyribonucleotide 5'-TGATTATCTG-3' conjugated at the 5'-end to CDPI3 and its complementary strand to an unmodified control duplex of the same sequence using nuclear magnetic resonance techniques. Thermal denaturation studies indicated that the hybrid of this conjugate with its complementary strand had a melting temperature that was 30 degrees C higher compared with the unmodified control duplex. Following restrained molecular dynamics and relaxation matrix refinement, the solution structure of the CDPI3-conjugated DNA duplex demonstrated that the overall shape of the duplex was that of a straight B-type helix and that the CDPI3moiety was bound snugly in the minor groove, where it was stabilized by extensive van der Waal's interactions.     

Refined solution structure and ligand-binding properties of PDC-109 domain b. A collagen-binding type II domain.

We have determined, via 1H-n.m.r., the solution conformation of the collagen-binding b-domain of the bovine seminal fluid protein PDC-109 (PDC-109/b). The structure determination is based on 341 interproton distance estimates and 42 dihedral angle estimates: a set of 24 initial structures were computed; 12 using the variable target function program DIANA, and 12 using the metric matrix program DISGEO. These structures were optimized by restrained energy minimization and dynamic simulated annealing using the CHARMM and X-PLOR programs. The average pairwise root-mean-square difference (r.m.s.d) between the optimized DIANA (DISGEO) structures is 0.71 A (0.82 A) for the backbone atoms, and 1.73 A (2.03 A) for all atoms. Both sets of structures exhibit the same global fold, secondary structure and placement of most non-polar side-chains. Two central antiparallel beta-sheets, which lie roughly perpendicular to each other, and two irregular loops support a large, partially exposed, hydrophobic surface that defines a putative binding site. A test of a hybrid relaxation matrix-based distance refinement protocol (MIDGE program) was performed using a normalized 250 millisecond NOESY spectrum. The resulting distances were input to the molecular mechanics/dynamics procedures mentioned above in order to optimize the DIANA structures. Our results indicate that relaxation matrix refinement of distances is most useful when used conservatively for identifying underestimated distance constraints. 1H-n.m.r. monitored ligand titration experiments revealed definite, albeit weak, binding interactions for phenethylamine and leucine analogs (Ka less than or equal to 25 M-1). Residues perturbed by ligand binding include Tyr7, Trp26, Tyr33, Asp34 and Trp39. These results suggest that PDC-109/b may recognize specific leucine and/or isoleucine-containing sequences within collagen.     

Relaxation matrix refinement of the solution structure of squash trypsin inhibitor

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.     

Dynamic NMR structures of [Rp]- and [Sp]-phosphorothioated DNA-RNA hybrids: is flexibility required for RNase H recognition?

Chemically modified DNA oligonucleotides have been crucial to the development of antisense therapeutics. High-resolution structural studies of pharmaceutically relevant derivatives have been limited to only a few molecules. We have used NMR to elucidate the structure in solution of two DNA-RNA hybrids with the sequence d(CCTATAATCC).r(GGAUUAUAGG). The two hybrids contain an unmodified RNA target strand, whereas the DNA strand contains one of two different stereoregular sugar-phosphate backbone linkages at each nucleotide: 1), [Rp]-phosphorothioate or 2), [Sp]-phosphorothioate. Homonuclear two-dimensional spectroscopy afforded nearly complete nonlabile proton assignments. Distance bounds, calculated from the nuclear Overhauser effect (NOE) crosspeak intensities via a complete relaxation matrix approach with the program MARDIGRAS, were used to restrain the structure of the two hybrids during simulations of molecular dynamics. Analysis of restrained molecular dynamics trajectories suggests that both hybrids are flexible, requiring the use of molecular dynamics with time-averaged restraints (MDtar) to generate ensembles of structures capable of satisfying the NMR data. In particular, the deoxyribose sugars of the DNA strand show strong evidence of repuckering. Furthermore, deoxyribose sugar repuckering is accompanied by increased flexibility of overall helical geometry. These observations, together with the analysis of the crystal structure of a hybrid duplex in complex with ribonuclease H (RNase H), suggested that this flexibility may be required for recognition by RNase H.     

Direct NOE refinement of biomolecular structures using 2D NMR data

Summary A set of computer programs called DINOSAUR has been developed, which allows the refinement of biomolecular structures directly from 2D NOE intensities. The NOE restraining potential implemented emphasises the weak intensities corresponding to larger distances which are more likely to determine the three-dimensional structure. An approximation based on a two-spin approximation is proposed for the gradient of the NOE intensities instead of the exact solution which is extremely time-consuming. The DINOSAUR routines have been implemented in various refinement programs (Distance bound Driven Dynamics, Molecular Dynamics and Energy Minimisation) and tested on an eight-residue model peptide.     

Solution structure of the EcoRI DNA octamer containing 5-fluorouracil via restrained molecular dynamics using distance and torsion angle constraints extracted from NMR spectral simulations.

The self-complementary DNA octamer [d(GGAATUFCC)]2, containing the EcoRI recognition sequence with one of the thymines replaced by 5-fluorouracil (UF), was synthesized. Proton homonuclear two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered correlation (2QF-COSY) spectra, as well as one-dimensional spectra at different temperatures, were recorded for the octamer. Consequently, all proton resonances were assigned. The thermally induced transition from the duplex to single strands has been followed, demonstrating the stability of the duplex containing 5-fluorouracil. Simulations of the 2QF-COSY cross-peaks by means of the programs SPHINX and LINSHA were compared with experimental data, establishing scalar coupling constants for the sugar ring protons and hence sugar pucker parameters. The deoxyribose rings exhibit a dynamic equilibrium of N- and S-type conformers with 75-95% populations of the latter. Two programs used for complete relaxation matrix analysis 2D NOE spectra, CORMA and MARDIGRAS, were modified to account for the influence of the fluorines on dipolar interactions in the proton system. Quantitative assessment of the 2D NOE cross-peak intensities for different mixing times, in conjunction with the program MARDIGRAS, gave a set of interproton distances for each mixing time. The largest and smallest values of each of the interproton distances were chosen as the upper and lower bounds for each distance constraint. The distance bounds define the size of a flat-well potential function term, incorporated into the AMBER force field, which was employed for restrained molecular dynamics calculations. Torsion angle constraints in the form of a flat-well potential were also constructed from the analysis of the sugar pucker data. Several restrained molecular dynamics runs of 35 ps were performed, utilizing 284 experimental distance and torsion angle constraints and two different starting structures, energy-minimized A- and B-DNA. Convergence to similar structures with a root-mean-square deviation of 1.2 A was achieved for the central hexamer of the octamer, starting from A- and B-DNA. The average structure from six different molecular dynamics runs was subjected to final restrained energy minimization. The resulting final structure was in good agreement with the structures derived from different molecular dynamics runs and showed a substantial improvement of the 2D NOE sixth-root residual index in comparison with classical and energy-minimized B-DNA. A detailed analysis of the conformation of the final structure and comparison with structures of similar sequences, obtained by different methods, were performed.     

Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement

Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.     

Solution structure of a nitrous acid induced DNA interstrand cross-link

Nitrous acid is a mutagenic agent. It can induce interstrand cross-links in duplex DNA, preferentially at d(CpG) steps: two guanines on opposite strands are linked via a single shared exocyclic imino group. Recent synthetic advances have led to the production of large quantities of such structurally homogenous cross-linked duplex DNA. Here we present the high resolution solution structure of the cross-linked dodecamer [d(GCATCCGGATGC)]2 (the cross-linked guanines are underlined), determined by 2D NMR spectroscopy, distance geometry, restrained molecular dynamics and iterative NOE refinement. The cross-linked guanines form a nearly planar covalently linked 'G:G base pair' with only minor propeller twisting, while the cytidine bases of their normal base pairing partners have been flipped out of the helix and adopt well defined extrahelical positions in the minor groove. On the 5'-side of the cross-link, the minor groove is widened to accommodate these extrahelical bases, and the major groove becomes quite narrow at the cross-link. The cross-linked 'G:G base pair' is well stacked on the spatially adjacent C:G base pairs, particularly on the 3'-side guanines. In addition to providing the first structure of a nitrous acid cross-link in DNA, these studies could be of major importance to the understanding of the mechanisms of nitrous acid cross-linking and mutagenicity, as well as the mechanisms responsible for its repair in intracellular environments. It is also the shortest DNA cross-link structure to be described.     

Solution structure of duplex DNA containing the mutagenic lesion 1,N(2)-etheno-2'-deoxyguanine

We have used high-resolution NMR spectroscopy and molecular dynamics simulations to determine the solution structure of DNA containing the genotoxic lesion 1, N (2)-etheno-2'-deoxyguanosine (epsilonG), paired to dC. The NMR data suggest the presence of a major, minimally perturbed structure at neutral pH. NOESY spectra indicate the presence of a right-handed helix with all nucleotides in anti, 2'-deoxyribose conformations within the C2'-endo/C1'-exo range and proper Watson-Crick base pair alignments outside the lesion site. The epsilonG residue remains deeply embedded inside the helix and stacks between the flanking base pairs. The lesion partner dC is extrahelical and is located in the minor groove of the duplex, where it is highly exposed to solvent. Upon acidification of the sample, a second conformation at the lesion site of the duplex emerges, with protonation of the lesion partner dC and possible formation of a Hoogsteen base pair. Restrained molecular dynamics simulations of the neutral-pH structure generated a set of three-dimensional models that show epsilonG inside the helix, where the lesion is stabilized by stacking interactions with flanking bases but without participating in hydrogen bonding. The lesion counterbase dC is displaced in the minor groove of the duplex where it can form a hydrogen bond with the sugar O4' atom of a residue 2 bp away.