Solution structure of a DNA duplex containing 8-hydroxy-2'-deoxyguanosine opposite deoxyguanosine

MIM-ML: A Novel Quantum Chemical Fragment-Based Random Forest Model for Accurate Prediction of NMR Chemical Shifts of Nucleic Acids

We developed a random forest machine learning (ML) model for the prediction of 1H and 13C NMR chemical shifts of nucleic acids. Our ML model is trained entirely on reproducing computed chemical shifts obtained previously on 10 nucleic acids using a Molecules-in-Molecules (MIM) fragment-based density functional theory (DFT) protocol including microsolvation effects. Our ML model includes structural descriptors as well as electronic descriptors from an inexpensive low-level semiempirical calculation (GFN2-xTB) and trained on a relatively small number of DFT chemical shifts (2080 1H chemical shifts and 1780 13C chemical shifts on the 10 nucleic acids). The ML model is then used to make chemical shift predictions on 8 new nucleic acids ranging in size from 600 to 900 atoms and compared directly to experimental data. Though no experimental data was used in the training, the performance of our model is excellent (mean absolute deviation of 0.34 ppm for 1H chemical shifts and 2.52 ppm for 13C chemical shifts for the test set), despite having some nonstandard structures. A simple analysis suggests that both structural and electronic descriptors are critical for achieving reliable predictions. This is the first attempt to combine ML from fragment-based DFT calculations to predict experimental chemical shifts accurately, making the MIM-ML model a valuable tool for NMR predictions of nucleic acids.

Accurate and Cost-Effective NMR Chemical Shift Predictions for Nucleic Acids Using a Molecules-in-Molecules Fragmentation-Based Method

We have developed, implemented, and assessed an efficient protocol for the prediction of NMR chemical shifts of large nucleic acids using our molecules-in-molecules (MIM) fragment-based quantum chemical approach. To assess the performance of our approach, MIM-NMR calculations are calibrated on a test set of three nucleic acids, where the structure is derived from solution-phase NMR studies. For DNA systems with multiple conformers, the one-layer MIM method with trimer fragments (MIM1trimer) is benchmarked to get the lowest energy structure, with an average error of only 0.80 kcal/mol with respect to unfragmented full molecule calculations. The MIMI-NMRdimer calibration with respect to unfragmented full molecule calculations shows a mean absolute deviation (MAD) of 0.06 and 0.11 ppm, respectively, for 1H and 13C nuclei, but the performance with respect to experimental NMR chemical shifts is comparable to the more expensive MIM1-NMR and MIM2-NMR methods with trimer subsystems. To compare with the experimental chemical shifts, a standard protocol is derived using DNA systems with Protein Data Bank (PDB) IDs 1SY8, 1K2K, and 1KR8. The effect of structural minimizations is employed using a hybrid mechanics/semiempirical approach and used for computations in solution with implicit and explicit-implicit solvation models in our MIM1-NMRdimer methodology. To demonstrate the applicability of our protocol, we tested it on seven nucleic acids, including structures with nonstandard residues, heteroatom substitutions (F and B atoms), and side chain mutations with a size ranging from ∼300 to 1100 atoms. The major improvement for predicted MIM1-NMRdimer calculations is obtained from structural minimizations and implicit solvation effects. A significant improvement with the explicit-implicit solvation model is observed only for two smaller nucleic acid systems (1KR8 and 7NBK), where the expensive first solvation shell is replaced by the microsolvation model, in which a single water molecule is added for each solvent-exposed amino and imino protons, along with the implicit solvation. Overall, our target accuracy of ∼0.2-0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for large nucleic acids. The proposed MIM-NMR approach is accurate and cost-effective (linear scaling with system size), and it can aid in the structural assignments of a wide range of complex biomolecules.

Collision-induced dissociation of homodimeric and heterodimeric radical cations of 9-methylguanine and 9-methyl-8-oxoguanine: correlation between intra-base pair proton transfer originating from the N1-H at a Watson-Crick edge and non-statistical dissociation

It has been shown previously in protonated, deprotonated and ionized guanine-cytosine base pairs that intra-base pair proton transfer from the N1-H at the Watson-Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure-kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e., [9MG·9MG]˙⁺, [9MOG·9MG]˙⁺ and [9MOG·9MOG]˙⁺. Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion-molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1-H at the Watson-Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer.

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo-7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm ) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

Automated AFM analysis of DNA bending reveals initial lesion sensing strategies of DNA glycosylases

Base excision repair is the dominant DNA repair pathway of chemical modifications such as deamination, oxidation, or alkylation of DNA bases, which endanger genome integrity due to their high mutagenic potential. Detection and excision of these base lesions is achieved by DNA glycosylases. To investigate the remarkably high efficiency in target site search and recognition by these enzymes, we applied single molecule atomic force microscopy (AFM) imaging to a range of glycosylases with structurally different target lesions. Using a novel, automated, unbiased, high-throughput analysis approach, we were able to resolve subtly different conformational states of these glycosylases during DNA lesion search. Our results lend support to a model of enhanced lesion search efficiency through initial lesion detection based on altered mechanical properties at lesions. Furthermore, its enhanced sensitivity and easy applicability also to other systems recommend our novel analysis tool for investigations of diverse, fundamental biological interactions.

7,8-Dihydro-8-oxoguanosine Lesions Inhibit the Theophylline Aptamer or Change Its Selectivity

Aptamers are attractive constructs due to their high affinity/selectivity towards a target. Here 7,8-dihydro-8-oxoguanosine (8-oxoG) has been used, due in part to its unique H-bonding capabilities (Watson-Crick or Hoogsteen), to expand the "RNA alphabet". Its impact on the theophylline RNA aptamer was explored by modifying its binding pocket at positions G11, G25, or G26. Structural probing, with RNases A and T1 , showed that modification at G11 leads to a drastic structural change, whereas the G25-/G26-modified analogues exhibited cleavage patterns similar to that of the canonical construct. The recognition properties towards three xanthine derivatives were then explored through thermophoresis. Modifying the aptamer at position G11 led to binding inhibition. Modification at G25, however, changed the selectivity towards theobromine (Kd ≈160 μm), with a poor affinity for theophylline (Kd >1.5 mm) being observed. Overall, 8-oxoG can have an impact on the structures of aptamers in a position-dependent manner, leading to altered target selectivity.

Polymorphic G:G mismatches act as hotspots for inducing right-handed Z DNA by DNA intercalation

DNA mismatches are highly polymorphic and dynamic in nature, albeit poorly characterized structurally. We utilized the antitumour antibiotic Co^II(Chro)₂ (Chro = chromomycin A3) to stabilize the palindromic duplex d(TTGGCGAA) DNA with two G:G mismatches, allowing X-ray crystallography-based monitoring of mismatch polymorphism. For the first time, the unusual geometry of several G:G mismatches including syn–syn, water mediated anti–syn and syn–syn-like conformations can be simultaneously observed in the crystal structure. The G:G mismatch sites of the d(TTGGCGAA) duplex can also act as a hotspot for the formation of alternative DNA structures with a GC/GA-5′ intercalation site for binding by the GC-selective intercalator actinomycin D (ActiD). Direct intercalation of two ActiD molecules to G:G mismatch sites causes DNA rearrangements, resulting in backbone distortion to form right-handed Z-DNA structures with a single-step sharp kink. Our study provides insights on intercalators-mismatch DNA interactions and a rationale for mismatch interrogation and detection via DNA intercalation.

Reactivity and Specificity of RNase T1, RNase A, and RNase H toward Oligonucleotides of RNA Containing 8-Oxo-7,8-dihydroguanosine

Understanding how oxidatively damaged RNA interacts with ribonucleases is important because of its proposed role in the development and progression of disease. Thus, understanding structural aspects of RNA containing lesions generated under oxidative stress, as well as its interactions with other biopolymers, is fundamental. We explored the reactivity of RNase A, RNase T₁, and RNase H toward oligonucleotides of RNA containing 8-oxo-7,8-dihydroguanosine (8oxoG). This is the first example that addresses this relationship and will be useful for understanding (1) how these RNases can be used to characterize the structural impact that this lesion has on RNA and (2) how oxidatively modified RNA may be handled intracellularly. 8-OxoG was incorporated into 10-16-mers of RNA, and its reactivity with each ribonuclease was assessed via electrophoretic analyses, circular dichroism, and the use of other C8-purine-modified analogues (8-bromoguanosine, 8-methoxyguanosine, and 8-oxoadenosine). RNase T₁ does not recognize sites containing 8-oxoG, while RNase A recognizes and cleaves RNA at positions containing this lesion while differentiating if it is involved in H-bonding. The selectivity of RNase A followed the order C > 8-oxoG ≈ U. In addition, isothermal titration calorimetry showed that an 8-oxoG-C3'-methylphosphate derivative can inhibit RNase A activity. Cleavage patterns obtained from RNase H displayed changes in reactivity in a sequence- and concentration-dependent manner and displayed recognition at sites containing the modification in some cases. These data will aid in understanding how this modification affects reactivity with ribonucleases and will enable the characterization of global and local structural changes in oxidatively damaged RNA.

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

Molecular dynamics simulation of the opposite-base preference and interactions in the active site of formamidopyrimidine-DNA glycosylase

Background

Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available.

Results

Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition.

Conclusion

Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0075-y) contains supplementary material, which is available to authorized users.

8-Oxo-7,8-dihydroadenine and 8-Oxo-7,8-dihydroadenosine-Chemistry, Structure, and Function in RNA and Their Presence in Natural Products and Potential Drug Derivatives

A description and history of the role that 8-oxo-7,8-dihydroadenine (8-oxoAde) and 8-oxo-7,8-dihydroadenosine (8-oxoA) have in various fields has been compiled. This Review focusses on 1) the formation of this oxidatively generated modification in RNA, its interactions with other biopolymers, and its potential role in the development/progression of disease; 2) the independent synthesis and incorporation of this modified nucleoside into oligonucleotides of RNA to display the progress that has been made in establishing its behavior in biologically relevant systems; 3) reported synthetic routes, which date back to 1890, along with the progress that has been made in the total synthesis of the nucleobase, nucleoside, and their corresponding derivatives; and 4) the isolation, total synthesis, and biological activity of natural products containing these moieties as the backbone. The current state of research regarding this oxidatively generated lesion as well as its importance in the context of RNA, natural products, and potential as drug derivatives is illustrated using all available examples reported to date.

Oxidative DNA Damage and Repair in the Radioresistant Archaeon Thermococcus gammatolerans

The hyperthermophilic archaeon Thermococcus gammatolerans can resist huge doses of γ-irradiation, up to 5.0 kGy, without loss of viability. The potential to withstand such harsh conditions is probably due to complementary passive and active mechanisms, including repair of damaged chromosomes. In this work, we documented the formation and repair of oxidative DNA lesions in T. gammatolerans. The basal level of the oxidized nucleoside, 8-oxo-2'-deoxyguanosine (8-oxo-dGuo), was established at 9.2 (± 0.9) 8-oxo-dGuo per 10⁶ nucleosides, a higher level than those usually measured in eukaryotic cells or bacteria. A significant increase in oxidative damage, i.e., up to 24.2 (± 8.0) 8-oxo-dGuo/10⁶ nucleosides, was measured for T. gammatolerans exposed to a 5.0 kGy dose of γ-rays. Surprisingly, the yield of radiation-induced modifications was lower than those previously observed for human cells exposed to doses corresponding to a few grays. One hour after irradiation, 8-oxo-dGuo levels were significantly reduced, indicating an efficient repair. Two putative base excision repair (BER) enzymes, TGAM_1277 and TGAM_1653, were demonstrated both by proteomics and transcriptomics to be present in the cells without exposure to ionizing radiation. Their transcripts were moderately upregulated after gamma irradiation. After heterologous production and purification of these enzymes, biochemical assays based on electrophoresis and MALDI-TOF (matrix-assisted laser desorption ionization-time of flight) mass spectrometry indicated that both have a β-elimination cleavage activity. TGAM_1653 repairs 8-oxo-dGuo, whereas TGAM_1277 is also able to remove lesions affecting pyrimidines (1-[2-deoxy-β-d-erythro-pentofuranosyl]-5-hydroxyhydantoin (5-OH-dHyd) and 1-[2-deoxy-β-d-erythro-pentofuranosyl]-5-hydroxy-5-methylhydantoin (5-OH-5-Me-dHyd)). This work showed that in normal growth conditions or in the presence of a strong oxidative stress, T. gammatolerans has the potential to rapidly reduce the extent of DNA oxidation, with at least these two BER enzymes as bodyguards with distinct substrate ranges.

Crosslinking reactions of 4-amino-6-oxo-2-vinylpyrimidine with guanine derivatives and structural analysis of the adducts

DNA interstrand crosslinks (ICLs) are the primary mechanism for the cytotoxic activity of many clinical anticancer drugs, and numerous strategies for forming ICLs have been developed. One such method is using crosslink-forming oligonucleotides (CFOs). In this study, we designed a 4-amino-6-oxo-2-vinylpyrimidine (AOVP) derivative with an acyclic spacer to react selectively with guanine. The AOVP CFO exhibited selective crosslinking reactivity with guanine and thymine in DNA, and with guanine in RNA. These crosslinking reactions with guanine were accelerated in the presence of CoCl₂, NiCl₂, ZnCl₂ and MnCl₂. In addition, we demonstrated that the AOVP CFO was reactive toward 8-oxoguanine opposite AOVP in the duplex DNA. The structural analysis of each guanine and 8-oxoguanine adduct in the duplex DNA was investigated by high-resolution NMR. The results suggested that AOVP reacts at the N2 amine in guanine and at the N1 or N2 amines in 8-oxoguanine in the duplex DNA. This study demonstrated the first direct determination of the adduct structure in duplex DNA without enzyme digestion.

The structural impact of DNA mismatches

The structure and dynamics of all the transversion and transition mismatches in three different DNA environments have been characterized by molecular dynamics simulations and NMR spectroscopy. We found that the presence of mismatches produced significant local structural alterations, especially in the case of purine transversions. Mismatched pairs often show promiscuous hydrogen bonding patterns, which interchange among each other in the nanosecond time scale. This therefore defines flexible base pairs, where breathing is frequent, and where distortions in helical parameters are strong, resulting in significant alterations in groove dimension. Even if the DNA structure is plastic enough to absorb the structural impact of the mismatch, local structural changes can be propagated far from the mismatch site, following the expected through-backbone and a previously unknown through-space mechanism. The structural changes related to the presence of mismatches help to understand the different susceptibility of mismatches to the action of repairing proteins.

Chemistry and structural biology of DNA damage and biological consequences

The formation of adducts by the reaction of chemicals with DNA is a critical step for the initiation of carcinogenesis. The structural analysis of various DNA adducts reveals that conformational and chemical rearrangements and interconversions are a common theme. Conformational changes are modulated both by the nature of adduct and the base sequences neighboring the lesion sites. Equilibria between conformational states may modulate both DNA repair and error-prone replication past these adducts. Likewise, chemical rearrangements of initially formed DNA adducts are also modulated both by the nature of adducts and the base sequences neighboring the lesion sites. In this review, we focus on DNA damage caused by a number of environmental and endogenous agents, and biological consequences.

Synthesis of N2-alkyl-8-oxo-7,8-dihydro-2'-deoxyguanosine derivatives and effects of these modifications on RNA duplex stability

N(2)-alkyl analogues of 8-oxo-7,8-dihydro-2'-deoxyguanosine (OG) were synthesized (alkyl = propyl, benzyl) via reductive amination of the protected OG nucleoside and incorporated into various positions of an RNA strand. Thermal stability studies of duplexes containing A or C opposite a single modified base revealed only moderate destabilization. Both OG as well as its N(2)-alkyl analogues can pair opposite A or C with nearly equal stability, potentially offering a new means of modulating RNA-protein interactions in the minor vs major grooves.

Impact of conformational heterogeneity of OxoG lesions and their pairing partners on bypass fidelity by Y family polymerases

7,8-Dihydro-8-oxoguanine (oxoG), the predominant oxidative DNA damage lesion, is processed differently by high-fidelity and Y-family lesion bypass polymerases. Although high-fidelity polymerases extend predominantly from an A base opposite an oxoG, the Y-family polymerases Dpo4 and human Pol eta preferentially extend from the oxoG*C base pair. We have determined crystal structures of extension Dpo4 ternary complexes with oxoG opposite C, A, G, or T and the next nascent base pair. We demonstrate that neither template backbone nor the architecture of the active site is perturbed by the oxoG(anti)*C and oxoG*A pairs. However, the latter manifest conformational heterogeneity, adopting both oxoG(syn)*A(anti) and oxoG(anti)*A(syn) alignment. Hence, the observed reduced primer extension from the dynamically flexible 3'-terminal primer base A is explained. Because of homology between Dpo4 and Pol eta, such a dynamic screening mechanism might be utilized by Dpo4 and Pol eta to regulate error-free versus error-prone bypass of oxoG and other lesions.

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.

Effect of hyperoxidized guanine on DNA primer-template structures: spiroiminodihydantoin leads to strand slippage

Oxidation of guanine in DNA can lead to mutagenic lesions such as 7-hydro-8-oxoguanine (oG). Upon further oxidation, a more mutagenic lesion, spirominodihydantoin (Sp), can occur. In this study, nuclear magnetic resonance (NMR) investigations were performed to determine the structural features of DNA primer-template models with 5'-GG, 5'-G(oG), 5'-G(Sp) and 5'-T(Sp) templates, that mimic the situation in which the downstream G of the template has been oxidized to oG or hyperoxidized to Sp. Our results show that misalignment occurs only in the 5'-G(Sp) and 5'-T(Sp) templates, providing structural insights into the observed differences in mutagenicity of Sp and oG during DNA replication.

A distinct TthMutY bifunctional glycosylase that hydrolyzes not only adenine but also thymine opposite 8-oxoguanine in the hyperthermophilic bacterium, Thermus thermophilus

Oxidative damage represents a major threat to genomic stability because the major product of DNA oxidation, 8-oxoguanine (GO), frequently mispairs with adenine during replication. We were interested in finding out how hyperthermophilic bacteria under goes the process of excising mispaired adenine from A/GO to deal with genomic oxidative damage. Herein we report the properties of an Escherichia coli MutY (EcMutY) homolog, TthMutY, derived from a hyperthermophile Thermus thermophilus. TthMutY preferentially excises on A/GO and G/GO mispairs and has additional activities on T/GO and A/G mismatches. TthMutY has significant sequence homology to the A/G and T/G mismatch recognition motifs, respectively, of MutY and Mig.MthI. A substitution from Tyr112 to Ser or Ala (Y112S and Y112A) in the putative thymine-binding site of TthMutY showed significant decrease in DNA glycosylase activity. A mutant form of TthMutY, R134K, could form a Schiff base with DNA and fully retained its DNA glycosylase activity against A/GO and A/G mispair. Interestingly, although TthMutY cannot form a trapped complex with substrate in the presence of NaBH(4), it expressed AP lyase activity, suggesting Tyr112 in TthMutY may be the key residue for AP lyase activity. These results suggest that TthMutY may be an example of a novel class of bifunctional A/GO mismatch DNA glycosylase that can also remove thymine from T/GO mispair.

Pa-AGOG, the founding member of a new family of archaeal 8-oxoguanine DNA-glycosylases

Oxidative damage represents a major threat to genomic stability, as the major product of DNA oxidation, 8-oxoguanine (GO), frequently mispairs with adenine during replication. In order to prevent these mutagenic events, organisms have evolved GO-DNA glycosylases that remove this oxidized base from DNA. We were interested to find out how GO is processed in the hyperthermophilic archaeon Pyrobaculum aerophilum, which lives at temperatures around 100 degrees C. To this end, we searched its genome for open reading frames (ORFs) bearing the principal hallmark of GO-DNA glycosylases: a helix-hairpin-helix motif and a glycine/proline-rich sequence followed by an absolutely conserved aspartate (HhH-GPD motif). Interestingly, although the P.aerophilum genome encodes three such ORFs, none of these encodes the potent GO-processing activity detected in P.aerophilum extracts. Fractionation of the extracts, followed by analysis of the active fractions by denaturing polyacrylamide gel electrophoresis, showed that the GO-processing enzyme has a molecular size of approximately 30 kDa. Mass spectrometric analysis of proteins in this size range identified several peptides originating from P.aerophilum ORF PAE2237. We now show that PAE2237 encodes AGOG (Archaeal GO-Glycosylase), the founding member of a new family of DNA glycosylases, which can remove GO from single- and double-stranded substrates with great efficiency.

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

Synthesis, properties, and NMR studies of a C8-phenylguanine modified oligonucleotide that preferentially adopts the Z DNA conformation

Carcinogenic aryl hydrazines produce C8-arylated purine adducts. The effect of these adducts on DNA conformation and their role in hydrazine carcinogenesis are unknown. Here, we describe a new synthetic route to produce these adducts that is also compatible with the synthesis of the corresponding phosphoramidites needed for oligonucleotide synthesis. Two oligonucleotides were prepared, an unmodified oligonucleotide, d((5)(')CGCGCGCGCG(3)(')), and a C8-phenylguanine modified oligonucleotide, d((5)(')CGCGCGCGCG(3)(')) (G = 8-phenylguanine). These oligonucleotides were compared using thermal denaturation, circular dichroism, NMR, and molecular modeling. The phenyl modification destabilizes the B DNA form and stabilizes the Z DNA form such that the B:Z ratio is near one under physiological conditions. In light of recent studies that show a role for Z DNA in gene expression and cell transformation, Z DNA stabilization by C8-arylguanine formation from aryl hydrazines may be relevant to their role in carcinogenesis.