Molecular dynamics simulation of 7, 8-dihydro-8-oxoguanine DNA

Hidden in plain sight: subtle effects of the 8-oxoguanine lesion on the structure, dynamics, and thermodynamics of a 15-base pair oligodeoxynucleotide duplex

The base lesion 8-oxoguanine is formed readily by oxidation of DNA, potentially leading to G → T transversion mutations. Despite the apparent similarity of 8-oxoguanine-cytosine base pairs to normal guanine-cytosine base pairs, cellular base excision repair systems effectively recognize the lesion base. Here we apply several techniques to examine a single 8-oxoguanine lesion at the center of a nonpalindromic 15-mer duplex oligonucleotide in an effort to determine what, if anything, distinguishes an 8-oxoguanine-cytosine (8oxoG-C) base pair from a normal base pair. The lesion duplex is globally almost indistinguishable from the unmodified parent duplex using circular dichroism spectroscopy and ultraviolet melting thermodynamics. The DNA mismatch-detecting photocleavage agent Rh(bpy)(2)chrysi(3+) cleaves only weakly and nonspecifically, revealing that the 8oxoG-C pair is locally stable at the level of the individual base pairs. Nuclear magnetic resonance spectra are also consistent with a well-conserved B-form duplex structure. In the two-dimensional nuclear Overhauser effect spectra, base-sugar and imino-imino cross-peaks are strikingly similar between parent and lesion duplexes. Changes in chemical shift due to the 8oxoG lesion are localized to its complementary cytosine and to the 2-3 bp immediately flanking the lesion on the lesion strand. Residues further removed from the lesion are shown to be unperturbed by its presence. Notably, imino exchange experiments indicate that the 8-oxoguanine-cytosine pair is strong and stable, with an apparent equilibrium constant for opening equal to that of other internal guanine-cytosine base pairs, on the order of 10(-6). This collection of experiments shows that the 8-oxoguanine-cytosine base pair is incredibly stable and similar to the native pair.

Structural and energetic consequences of oxidation of d(ApGpGpGpTpT) telomere repeat unit in complex with TRF1 protein

The configuration hyperspace of canonical and oxidized 14-mers of B-DNA comprising telomere repeat units d(ApGpGpGpTpT) was sampled over 40 ns via molecular dynamic (MD) simulations. The energetic and structural consequences of TRF1 binding to telomere B-DNA were compared with non-complexed systems. Energetic properties of analyzed pairs, di- and tri-nucleotide steps occurring in central telomere repeat unit were estimated by means of advanced quantum chemistry computations including not only BSSE corrections, electron correlation contributions but also non-negligible many-body terms. These data along with bases pair and base step parameters distributions allow for quantization of consequences of oxidation and/or TRF1 binding to telomere repeat units. Occurrence of 8-oxoguanine in central telomeric triad (CTT) is the source of high stiffness if compared to non-modified oligomer. The origin of this property comes from significantly alteration of intermolecular interactions introduced by 8-oxoguanine. The increased stability observed for base-base interactions are accumulated and characterizes also di- and tri-nucleotides. The observed changes in the intermolecular interactions originate from structural alterations imposed by TRF1 binding to canonical and oxidized telomere B-DNA. First and most direct consequence of TRF1 binding to oxidized telomere repeat unit is alteration of shift-slide correlations if compared to canonical system. This in turn leads to large differences in purine-purine overlapping in oxidized structures. Thus, oxidized telomere B-DNA double strands are sensitive to interactions with protein ligands and numerous structural and energetic changes are imposed on base pairs forming CTT.

Impact of 7,8-dihydro-8-oxoguanine on methylation of the CpG site by Dnmt3a

7,8-Dihydro-8-oxoguanine (8-oxoG) is a ubiquitous oxidative DNA lesion resulting from injury to DNA via reactive oxygen species. 8-oxoG lesions may play a role in the formation of aberrant DNA methylation patterns during carcinogenesis. In this study, we assessed the effects of 8-oxoG on methylation and complex formation of nine 30-mer oligodeoxynucleotide duplexes by the catalytic domain of murine Dnmt3a DNA methyltransferase (Dnmt3a-CD). The effects of 8-oxoG on the methylation rate of hemimethylated duplexes varied from a 25-fold decrease to a 1.8-fold increase, depending on the position of the lesion relative to the Dnmt3a-CD recognition site (CpG) and target cytosine (C). The most significant effect was observed when 8-oxoG replaced guanine within the recognition site immediately downstream of the target cytosine. Fluorescence polarization experiments with fluorescein-labeled duplexes revealed that two molecules of Dnmt3a-CD bind per duplex, generating sigmoid binding curves. Duplexes exhibiting the highest apparent binding cooperativity formed the least stable 1:2 complexes with Dnmt3a-CD and were methylated at the lowest rate. Kinetic analyses disclosed the formation of very stable nonproductive enzyme-substrate complexes with hemimethylated duplexes that act as suicide substrates of Dnmt3a-CD. The presence of 8-oxoG within the CpG site downstream of the target cytosine markedly diminished productive versus nonproductive binding. We propose that 8-oxoG located adjacent to the target cytosine interferes with methylation by weakening the affinity of DNA for Dnmt3a-CD, thereby favoring a nonproductive binding mode.

How does modification of adenine by hydroxyl radical influence the stability and the nature of stacking interactions in adenine-cytosine complex?

This study reports on ab initio calculations of adenine-cytosine complexes in two different context alignments appearing in B-DNA. The influence of adenine modification by hydroxyl radical on the stability of the complexes is also discussed. The analysis was performed on over 40 crystallographic structures for each of the sequence contexts. In most cases, modification of adenine by hydroxyl radical leads to less negative intermolecular interaction energies. The issue of the influence of alteration of structural base step parameters on the stability of modified and unmodified adenine-cytosine complexes is also addressed. Analysis of the dependence of intermolecular interaction energy on base step parameters reveals that for twist and shift modification of adenine by hydroxyl radical leads to quite different interaction energy profiles in comparison with unmodified complexes. In order to elucidate the physical origins of this phenomenon, i.e. to analyze how the modification of adenine by hydroxyl radical is reflected in the change of intermolecular interaction energy components, a variational-perturbational decomposition scheme was applied at the MP2/aug-cc-pVDZ level of theory.

Theoretical analysis of the effects of guanine oxidative damage on the properties of B-DNA telomere fragments

The molecular dynamics as well as ab initio MP2/6-31G(d = 0.25) single point calculations were performed for native and oxidized B-DNA telomeric fragments. The structural, dynamic, energetic and electrostatic properties along with frontier orbitals distribution were described for the central triad consisting of three guanine molecules in its canonical or oxidized forms. Although the average structural parameters characterizing all of the studied telomere fragments are close each to other, the significant consequence on angular and displacement flexibilities are observed. Namely, the increase of mutual displacement of two successive base pairs along either axis and increase of the rotation of two bases of opposite strand are main dynamic consequences of presence of 8-oxo-guanine in the central triad of telomeric B-DNA. Besides, the significant increase of stacking energies in case of 8-oxo-guanine were found. Furthermore, the guanine pattern visible from the major groove may be described as donor-acceptor-acceptor formed by H8-N7-O6 atoms, respectively. To the contrary the presence of 8-oxo-guanine changes the electrostatic properties of the major groove into acceptor-donor-acceptor coming from O8-H7-O6 atoms. This results in significant alteration of ESP characteristics. Finally, the molecular orbital properties are also significantly affected by oxidation of telomeric B-DNA fragments. All these factors contribute to decrease of binding of telomere proteins.

Structural and dynamic effects of single 7-hydro-8-oxoguanine bases located in a frameshift target DNA sequence

DNA 7-hydro-8-oxoguanine (8-oxoG) is implicated in frameshift formation in an G(6) sequence of the HPRT gene in mismatch repair (MMR) defective cells. Using oligonucleotides based on this frameshift hotspot, we investigated how a single 8-oxoG modified the structural and dynamic properties of the G(6) tract. A 30 ns molecular dynamics (MD) simulation indicated compression of the minor groove in the immediate vicinity of the lesion. Fluorescence polarization anisotropy (FPA) and MD demonstrated that 8-oxoG increases DNA torsional rigidity and also constrains the movement of the single-stranded region at the single/double stranded DNA junction of model DNA replication template/primer. These constraints influenced the efficiency of primer extension by Klenow (exo(-)) DNA polymerase.

Dynamic behavior of DNA base pairs containing 8-oxoguanine

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.

8-Oxoguanine enhances bending of DNA that favors binding to glycosylases

Molecular dynamics (MD) simulations were carried out on the DNA oligonucleotide GGGAACAACTAG:CTAGTTGTTCCC in its native form and with guanine in the central G(19):C(6) base pair replaced by 8-oxoguanine (8oxoG). A box of explicit water molecules was used for solvation, and Na(+) counterions were added to neutralize the system. The direction and magnitude of global bending were assessed by a technique used previously to analyze simulations of DNA containing a thymine dimer. The presence of 8oxoG did not greatly affect the magnitude of DNA bending; however, bending into the major groove was significantly more probable when 8oxoG replaced G(19). Crystal structures of glycosylases bound to damaged-DNA substrates consistently show a sharp bend into the major groove at the damage site. We conclude that changes in bending dynamics that assist the formation of this kink are a part of the mechanism by which glycosylases of the base excision repair pathway recognize the presence of 8oxoG in DNA.