Sequence-dependent dynamics in duplex DNA

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

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.     

Nuclear magnetic resonance structural studies and molecular modeling of duplex DNA containing normal and 4'-oxidized abasic sites

A 4'-oxidized abasic site (X) has been synthesized in a defined duplex DNA sequence, 5'-d(CCAAAGXACCGGG)-3'/3'-d(GGTTTCATGGCCC)-5' (1). Its structure has been determined by two-dimensional NMR methods, molecular modeling, and molecular dynamics simulations. 1 is globally B-form with the base (A) opposite X intrahelical and well-stacked. Only the alpha anomer of X is observed, and the abasic site deoxyribose is largely intrahelical. These results are compared with a normal abasic site (Y) in the same sequence context (2). Y is composed of a 60:40 mixture of alpha and beta anomers (2alpha and 2beta). In both 2alpha and 2beta, the base (A) opposite Y is intrahelical and well-stacked and the abasic site deoxyribose is predominantly extrahelical, consistent with the reported structures of the normal abasic site in a similar sequence context [Hoehn, S. T., Turner, C. J., and Stubbe, J. (2001) Nucleic Acids Res. 29, 3413-3423]. Molecular dynamics simulations reveal that the normal abasic site appears to be conformationally more flexible than the 4'-oxidized abasic site. The importance of the structure and flexibility of the abasic site in the recognition by the DNA repair enzyme Ape1 is discussed.     

Solution structure of duplex DNA containing an extrahelical abasic site analog determined by NMR spectroscopy and molecular dynamics.

Translesional DNA synthesis past abasic sites proceeds with the preferential incorporation of dAMP opposite the lesion and, depending on the sequence context, one or two base deletions. High-resolution NMR spectroscopy and molecular dynamics simulations were used to determine the three-dimensional structure of a DNA heteroduplex containing a synthetic abasic site (tetrahydrofuran) residue positioned in a sequence that promotes one base deletions. Analysis of NMR spectra indicates that the stem region of the duplex adopts a right-handed helical structure and the glycosidic torsion angle is in anti orientation for all residues. NOE interactions establish Watson-Crick alignments for all canonical base pairs of the duplex. Measurement of distance interactions at the lesion site shows the abasic residue excluded from the helix. Restrained molecular dynamics simulations generated three-dimensional models in excellent agreement with the spectroscopic data. These structures show a regular duplex region and a slight bend at the lesion site. The tetrahydrofuran residue extrudes from the helix and is highly flexible. The model reported here, in conjunction with a previous study performed on abasic sites, explains the structural bias of one-base deletion mutations.     

Structure of the DNA interstrand cross-link of 4,5',8-trimethylpsoralen.

4,5',8-Trimethylpsoralen (TMP) cross-links a 5' TpA or a 5' ApT site by photoreacting with one thymine moiety in each DNA strand. We are interested in whether psoralen interstrand cross-links all share one structure or whether there are significant differences. In this paper, we employed a rapid method for probing the structure of the cross-link by making a series of TMP cross-linked duplexes containing specific base-pair mismatches. The relative stability provided by a base pair can be correlated with neighboring base pairs by comparing the extents of gel retardation when base-pair mismatches happen in each position. From our studies, we infer that with respect to the furan-side strand, the 5'T.A base pair of the two T.A base pairs in the TpA site is not hydrogen bonded. Immediately on each side of the cross-linked TpA site is a highly stabilized base pair. Next, a region of decreased stability occurs in each arm of a cross-linked duplex and these base pairs of least stability are located farther away from the cross-linked thymines as the lengths of the arms of the cross-linked helix increase. Finally, even in 7 M urea at 49 degrees C the cross-linked helix is hydrogen bonded at both ends of a duplex of 22 base pairs. We propose that the structures of interstrand cross-links in DNA vary appreciably with the DNA sequence, the length of the DNA duplex, and the structures of the DNA cross-linking agents.     

Cleavage of oligoribonucleotides by a ribozyme derived from the hepatitis delta virus RNA sequence.

A self-cleaving RNA sequence from hepatitis delta virus was modified to produce a ribozyme capable of catalyzing the cleavage of RNA in an intermolecular (trans) reaction. The delta-derived ribozyme cleaved substrate RNA at a specific site, and the sequence specificity could be altered with mutations in the region of the ribozyme proposed to base pair with the substrate. A substrate target size of approximately 8 nucleotides in length was identified. Octanucleotides containing a single ribonucleotide immediately 5' to the cleavage site were substrates for cleavage, and cleavage activity was significantly reduced only with a guanine base at that position. A deoxyribose 5' to the cleavage site blocked the reaction. These data are consistent with a proposed secondary structure for the self-cleaving form of the hepatitis delta virus ribozyme in which a duplex forms with sequences 3' to the cleavage site, and they support a proposed mechanism in which cleavage involves attack on the phosphorus at the cleavage site by the adjacent 2'-hydroxyl group.     

Nitrosamine-induced cancer: O4-alkylthymine produces sites of DNA hyperflexibility.

The carcinogenic properties of N-nitroso compounds are associated with their ability to alkylate DNA, in particular to form O6-alkylguanine and O4-alkylthymine. DNA duplexes containing either O6-alkylguanine or O4-alkylthymine were synthesized, and each duplex was ligated to form a set of DNAs of increasing length with the alkylated base out of phase (16 base-pairs apart) or in phase (21 base-pairs apart) with the helical repeat of the DNA. The DNA contained the sequence 5' CAA 3', which is the 61st codon of the K-ras gene, because this codon is a preferred site of mutation for a number of carcinogens including the methylating carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1- butanone (NNK). O4-Methylthymine or O4-ethylthymine replaced thymine in either of the two A.T base-pairs of this codon (normally CAA), and O6-methylguanine replaced the guanine in the G.C pair. All the sequences containing O4-alkylthymine exhibited anomalous, slow, gel migration and ligated to form circles of unusually small diameter. In general, the effect was seen when the alkylated base-pair was out of phase with the helical repeat as well as when it was in phase, suggesting that the alkylated base-pair confers flexibility which is largely isotropic, i.e., has no preferred direction, rather than anisotropic flexibility or bending. However, at pH 8.3 the 21-base-pair set containing O4-alkylT.A had significantly greater anomalous migration than the 16-base-pair set, suggesting that the flexibility produced by this base-pair has a significant anisotropic component and thus resembles true bending.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure and dynamics of a DNA duplex containing single alpha-anomeric deoxyadenosine residue.

Structure and dynamics of an undecamer DNA duplex containing a single alpha-anomeric deoxyadenosine residue placed in opposition to a thymidine unit have been studied using simulation of molecular dynamics in aqueous solution. Despite several noticeable deviations from the B-DNA duplex structure caused by the anomerisation, such as: West type puckering of the alpha-anomeric sugar, disrupted base stacking pattern and unstable duplex bending, the formation of a non-classical alpha-dA-T pair was observed. A novel way of visual presentation of trajectory data allowing high throughput screening of the conformational parameters is presented.     

Structure of the 1,4-bis(2'-deoxyadenosin-N6-yl)-2R,3R-butanediol cross-link arising from alkylation of the human N-ras codon 61 by butadiene diepoxide

The solution structure of the 1,4-bis(2'-deoxyadenosin-N(6)-yl)-2R,3R-butanediol cross-link arising from N(6)-dA alkylation of nearest-neighbor adenines by butadiene diepoxide (BDO(2)) was determined in the oligodeoxynucleotide 5'-d(CGGACXYGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained codon 61 (underlined) of the human N-ras protooncogene. The cross-link was accommodated in the major groove of duplex DNA. At the 5'-side of the cross-link there was a break in Watson-Crick base pairing at base pair X(6).T(17), whereas at the 3'-side of the cross-link at base pair Y(7).T(16), base pairing was intact. Molecular dynamics calculations carried out using a simulated annealing protocol, and restrained by a combination of 338 interproton distance restraints obtained from (1)H NOESY data and 151 torsion angle restraints obtained from (1)H and (31)P COSY data, yielded ensembles of structures with good convergence. Helicoidal analysis indicated an increase in base pair opening at base pair X(6).T(17), accompanied by a shift in the phosphodiester backbone torsion angle beta P5'-O5'-C5'-C4' at nucleotide X(6). The rMD calculations predicted that the DNA helix was not significantly bent by the presence of the four-carbon cross-link. This was corroborated by gel mobility assays of multimers containing nonhydroxylated four-carbon N(6),N(6)-dA cross-links, which did not predict DNA bending. The rMD calculations suggested the presence of hydrogen bonding between the hydroxyl group located on the beta-carbon of the four-carbon cross-link and T(17) O(4), which perhaps stabilized the base pair opening at X(6).T(17) and protected the T(17) imino proton from solvent exchange. The opening of base pair X(6).T(17) altered base stacking patterns at the cross-link site and induced slight unwinding of the DNA duplex. The structural data are interpreted in terms of biochemical data suggesting that this cross-link is bypassed by a variety of DNA polymerases, yet is significantly mutagenic [Kanuri, M., Nechev, L. V., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580].     

1H NMR and optical spectroscopic investigation of the sequence-dependent dimerization of a symmetrical cyanine dye in the DNA minor groove

A symmetrical cyanine dye was previously shown to bind as a cofacial dimer to alternating A-T sequences of duplex DNA. Indirect evidence suggested that dimerization of the dye occurred in the minor groove. 1H NMR experiments reported here verify this model based on broadening and shifting of signals due to protons on carbon 2 of adenine and imino protons at the central five A-T pairs of the 11 base pair duplex: 5'-GCGTATATGCG-3'/3'-CGCATATACGC-5'. This binding mode is similar to that of distamycin A, even though the dye lacks the hydrogen-bonding groups used by distamycin for sequence-specific recognition. Surprisingly, the third base pair (G-C) was also implicated in the binding site. UV-vis experiments were used to compare the extent of dimerization of the dye for 11 different sequence variants. These experiments verified the importance of a G-C pair at the third position: replacing this pair with A-T suppressed dimerization. These results indicate that the dye binding site spans six base pairs: 5'-GTATAT-3'. The initial G-C pair seems to be important for widening the minor groove rather than for making important contacts with the dye molecules since inverting its orientation to C-G or replacing it with I-C still led to favorable dimerization of the dye.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)₂ shows unique stacking patterns of a stable B↔A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)₂ duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B↔A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.     

Using NMR and molecular dynamics to link structure and dynamics effects of the universal base 8-aza, 7-deaza,N8 linked adenosine analog

A truly universal nucleobase enables a host of novel applications such as simplified templates for PCR primers, randomized sequencing and DNA based devices. A universal base must pair indiscriminately to each of the canonical bases with little or preferably no destabilization of the overall duplex. In reality, many candidates either destabilize the duplex or do not base pair indiscriminatingly. The novel base 8-aza-7-deazaadenine (pyrazolo[3,4-d]pyrimidin- 4-amine) N⁸-(2′deoxyribonucleoside), a deoxyadenosine analog (UB), pairs with each of the natural DNA bases with little sequence preference. We have utilized NMR complemented with molecular dynamic calculations to characterize the structure and dynamics of a UB incorporated into a DNA duplex. The UB participates in base stacking with little to no perturbation of the local structure yet forms an unusual base pair that samples multiple conformations. These local dynamics result in the complete disappearance of a single UB proton resonance under native conditions. Accommodation of the UB is additionally stabilized via heightened backbone conformational sampling. NMR combined with various computational techniques has allowed for a comprehensive characterization of both structural and dynamic effects of the UB in a DNA duplex and underlines that the UB as a strong candidate for universal base applications.     

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.