Four hundred twenty-eight high-resolution DNA–protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π–contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π–π stacking, π–π T-shaped or sugar–π contacts. Three hundred forty-four nucleobase–amino acid π–π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA–protein sugar–π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C–H···π and/or lone–pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π–π and sugar–π interactions display a range of relative monomer orientations and therefore interaction energies (up to –50 (–70) kJ mol−1 for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA–protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

The β/δ-lyase activity of bifunctional glycosylases on damaged nucleotides in DNA involves the formation of a covalent linkage between the protein (lysine or N-terminal proline) and DNA (C1' of the damaged nucleotide). In the present study, the conformational and mechanistic flexibility of the cross-link is examined. Repair of 8-oxoguanine damage by hOgg1 is considered as a representative system, and the glycosylase through β-lyase steps are investigated using density functional theory. (PCM/SMD)-M06-2X/6-311+G(2df,2p)//PCM-B3LYP/6-31G(d) energetics were determined for eight unique mechanisms differing in the conformation of the imine linkage (E/Z), the proton (pro-S/R) abstracted during elimination, and whether the ring-opening step is base catalyzed. This initial study used a model system limited to the damaged nucleoside 3'-monophosphate and a model nucleophile to investigate this series of complex reaction steps. The great flexibility exhibited by the linkage and clustered β-elimination energetics indicate sterics will play a large role in predicting the preferred lyase mechanism for a given enzyme. The stationary points identified herein can be overlaid into a protein structure to assist in generating initial guesses for large model systems. By comparing the characterized geometries and enzyme active sites, methods for catalysis of the various chemical steps can be identified, and these possibilities are discussed in detail for hOgg1. Interestingly, the most stable structure on the potential energy surface occurs before elimination of the 3'-phosphate. Hydrolysis of the protein-DNA cross-link at this point would yield an AP-site, which provides support for the recently observed monofunctional activity of hOgg1.

Density functional theory (B3LYP) was employed to examine the hydrolysis of the canonical 2'-deoxynucleotides in varied environments (gas phase or water) using different computational models for the sugar residue (methyl or phosphate group at C5') and nucleophile (water activated through full or partial proton abstraction). Regardless of the degree of nucleophile activation, our results show that key geometrical parameters along the reaction pathway are notably altered upon direct inclusion of solvent effects in the optimization routine, which leads to significant changes in the reaction energetics and better agreement with experiment. Therefore, despite the wide use of gas-phase calculations in the literature, small model computational work, as well as large-scale enzyme models, that strive to understand nucleotide deglycosylation must adequately describe the environment. Alternatively, although inclusion of the phosphate group at C5' also affects the geometries of important stationary points, the effects cancel to yield unchanged deglycosylation barriers, and therefore smaller computational models can be used to estimate the energy associated with nucleotide deglycosylation, with the 5' phosphate group included if full (geometric) details of the reaction are desired. Hydrogen-bonding interactions with the nucleobase can significantly reduce the barrier to deglycosylation, which supports suggestions that discrete hydrogen-bonding interactions with active-site amino acid residues can play a significant role in enzyme-catalyzed nucleobase excision. Taken together with previous studies, the present work provides vital clues about the components that must be included in future studies of the deglycosylation of isolated noncanonical nucleotides, as well as the corresponding enzyme-catalyzed reactions.

8-Oxoguanine (OG) is one of the most frequently occurring forms of DNA damage and is particularly deleterious since it forms a stable Hoogsteen base pair with adenine (A). The repair of an OG:A mispair is initiated by adenine-DNA glycosylase (MutY), which hydrolyzes the sugar-nucleobase bond of the adenine residue before the lesion is processed by other proteins. MutY has been proposed to use a two-part chemical step involving protonation of the adenine nucleobase, followed by SN1 hydrolysis of the glycosidic bond. However, differences between a recent (fluorine recognition complex, denoted as the FLRC) crystal structure and the structure on which most mechanistic conclusions have been based to date (namely, the lesion recognition complex or LRC) raise questions regarding the mechanism used by MutY and the discrete role of various active-site residues. The present work uses both molecular dynamics (MD) and quantum mechanical (ONIOM) models to compare the active-site conformational dynamics in the two crystal structures, which suggests that only the understudied FLRC leads to a catalytically competent reactant. Indeed, all previous computational studies on MutY have been initiated from the LRC structure. Subsequently, for the first time, various mechanisms are examined with detailed ONIOM(M06-2X:PM6) reaction potential energy surfaces (PES) based on the FLRC structure, which significantly extends the mechanistic picture. Specifically, our work reveals that the reaction proceeds through a different route than the commonly accepted mechanism and the catalytic function of various active-site residues (Geobacillus stearothermophilus numbering). Specifically, contrary to proposals based on the LRC, E43 is determined to solely be involved in the initial adenine protonation step and not the deglycosylation reaction as the general base. Additionally, a novel catalytic role is proposed for Y126, whereby this residue plays a significant role in stabilizing the highly charged active site, primarily through interactions with E43. More importantly, D144 is found to explicitly catalyze the nucleobase dissociation step through partial nucleophilic attack. Although this is a more direct role than previously proposed for any other DNA glycosylase, comparison to previous work on other glycosylases justifies the larger contribution in the case of MutY and allows us to propose a unified role for the conserved Asp/Glu in the DNA glycosylases, as well as other enzymes that catalyze nucleotide deglycosylation reactions.

Although DNA damage can have a variety of deleterious effects on cells (e.g., senescence, death, and rapid growth), the base excision repair (BER) pathway combats the effects by removing several types of damaged DNA. Since the first BER step involves cleavage of the bond between the damaged nucleobase and the DNA sugar-phosphate backbone, we have used density functional theory to compare the intrinsic stability of the glycosidic bond in a number of common DNA oxidation, deamination, and alkylation products to the corresponding natural nucleosides. Our calculations predict that the dissociative (SN1) and associative (SN2) pathways are nearly isoenergetic, with the dissociative pathway only slightly favored on the Gibbs reaction surface for all canonical and damaged nucleosides, which suggests that DNA damage does not affect the inherently most favorable deglycosylation pathway. More importantly, with the exception of thymine glycol, all DNA lesions exhibit reduced glycosidic bond stability relative to the undamaged nucleosides. Furthermore, the trend in the magnitude of the deglycosylation barrier reduction directly correlates with the relative nucleobase acidity (at N9 for purines or N1 for pyrimidines), which thereby provides a computationally efficient, qualitative measure of the glycosidic bond stability in DNA damage. The effect of nucleobase activation (protonation) at different sites predicts that the positions leading to the largest reductions in the deglycosylation barrier are typically used by DNA glycosylases to facilitate base excision. Finally, deaza purine derivatives are found to have greater glycosidic bond stability than the canonical counterparts, which suggests that alterations to excision rates measured using these derivatives to probe DNA glycosylase function must be interpreted in reference to the inherent differences in the nucleoside reactivity. Combined with previous studies of the deglycosylation of DNA nucleosides, the current study provides a greater fundamental understanding about the reactivity of the glycosidic bond in damaged DNA, which has direct implications to the function of critical DNA repair enzymes.