Quantum theoretical study of cleavage of the glycosidic bond of 2'-deoxyadenosine: base excision-repair mechanism of DNA by MutY

Preorganized Internal Electric Field Promotes a Double-Displacement Mechanism for the Adenine Excision Reaction by Adenine DNA Glycosylase

Adenine DNA glycosylase (MutY) is a monofunctional glycosylase, removing adenines (A) misinserted opposite 8-oxo-7,8-dihydroguanine (OG), a common product of oxidative damage to DNA. Through multiscale calculations, we decipher a detailed adenine excision mechanism of MutY that is consistent with all available experimental data, involving an initial protonation step and two nucleophilic displacement steps. During the first displacement step, N-glycosidic bond cleavage is accompanied by the attack of the carboxylate group of residue Asp144 at the anomeric carbon (C1'), forming a covalent glycosyl-enzyme intermediate to stabilize the fleeting oxocarbenium ion. After departure of the excised base, water nucleophiles can be recruited to displace Asp144, completing the catalytic cycle with retention of stereochemistry at the C1' position. The two displacement reactions are found to mostly involve the movement of the oxocarbenium ion, occurring with large charge reorganization and thus sensitive to the internal electric field (IEF) exerted by the polar protein environment. Intriguingly, we find that the negatively charged carboxylate group is a good nucleophile for the oxocarbenium ion, yet an unactivated water molecule is not, and that the electric field catalysis strategy is used by the enzyme to enable its unique double-displacement reaction mechanism. A strong IEF, pointing toward 5' direction of the substrate sugar ring, greatly facilitates the second displacement reaction at the expense of elevating the barrier of the first one, thereby allowing both reactions to occur. These findings not only increase our understanding of the strategies used by DNA glycosylases to repair DNA lesions, but also have important implications for how internal/external electric field can be applied to modulate chemical reactions.

Distinctive Formation of a DNA-Protein Cross-Link during the Repair of DNA Oxidative Damage: Insights into Human Disease from MD Simulations and QM/MM Calculations

Reactive oxygen species damage DNA and result in health issues. The major damage product, 8-oxo-7,8-dihydroguanine (8oG), is repaired by human adenine DNA glycosylase homologue (MUTYH). Although MUTYH misfunction is associated with a genetic disorder called MUTYH-associated polyposis (MAP) and MUTYH is a potential target for cancer drugs, the catalytic mechanism required to develop disease treatments is debated in the literature. This study uses molecular dynamics simulations and quantum mechanics/molecular mechanics techniques initiated from DNA-protein complexes that represent different stages of the repair pathway to map the catalytic mechanism of the wild-type MUTYH bacterial homologue (MutY). This multipronged computational approach characterizes a DNA-protein cross-linking mechanism that is consistent with all previous experimental data and is a distinct pathway across the broad class of monofunctional glycosylase repair enzymes. In addition to clarifying how the cross-link is formed, accommodated by the enzyme, and hydrolyzed for product release, our calculations rationalize why cross-link formation is favored over immediate glycosidic bond hydrolysis, the accepted mechanism for all other monofunctional DNA glycosylases to date. Calculations on the Y126F mutant MutY highlight critical roles for active site residues throughout the reaction, while investigation of the N146S mutant rationalizes the connection between the analogous N224S MUTYH mutation and MAP. In addition to furthering our knowledge of the chemistry associated with a devastating disorder, the structural information gained about the distinctive MutY mechanism compared to other repair enzymes represents an important step for the development of specific and potent small-molecule inhibitors as cancer therapeutics.

Development of a QM/MM(ABEEM) method combined with a polarizable force field to investigate the excision reaction mechanism of damaged thymine

This paper focuses on the development of a quantum mechanics/molecular mechanics method using the ABEEM polarizable force field (QM/MM(ABEEM) method) to investigate the excision reaction mechanism of damaged thymine. This method does not simply combine the QM method with the polarizable force field. A valence electronegativity piecewise function with the distance between atoms as a variable is introduced to describe the atomic partial charges, and changes greatly during the reaction process. At the same time, the charge transfer effect is treated using the condition of local charge conservation. Compared with the traditional QM/MM method, the QM/MM(ABEEM) method can more accurately simulate the polarization effect and charge transfer effect in the reaction process. Focusing on the controversial problems of the excision of damaged bases, six reaction pathways were designed for monofunctional and difunctional deglycosylation of neutral bases and protonated bases. The results show that the QM/MM(ABEEM) method accurately simulates the polarization effect, charge transfer effect, activation energy and other properties of the reaction process. The process in which the active residue Asp activates the nucleophile H₂O to attack the protonated base is the preferred path. The average activation energy and free activation energy of the protonated base are 7.00-14.00 kcal mol^-1 lower than that of the neutral base. The study in this paper is helpful to understand the mechanism of repair enzymes in repairing bases.

Development of QM/MM (ABEEM polarizable force field) method to simulate the excision reaction mechanism of damaged cytosine

DNA damages are regarded as having harmful effects on cell. The base excision repair mechanism combats these effects by removing damaged bases. The deglycosylation mechanism of excising damaged bases by DNA glycosylase and the state of the leaving base have been controversial. The enzymatic reaction of DNA glycosylase to remove the damaged bases involves not only the formation and breaking of chemical bonds, but also complex polarization effect and charge transfer, which cannot be accurately simulated by the QM/MM method combined with the fixed charge force field. This work has developed the ABEEM fluctuating polarizable force field combining with the QM method, that is (QM/MM[ABEEM]), to accurately simulate the proton transfer, charge transfer and the charge distribution. The piecewise function is used as the valence-state electronegativity in the QM/MM (ABEEM) to realize the accurate fitting of the charge distribution in reaction. And the charge transfer is accurately simulated by the local charge conservation conditions. Four deglycosylation mechanisms including the monofunctional and difunctional mechanisms of four neutral and protonated cytosine derivatives are explored. It is confirmed that the monofunctional mechanism of Asp-activated nucleophile water is a better deglycosylation mechanism and the base is protonated before the reaction occurs. Protonization of the base reduced the activation energy by 10.00-17.00 kcal/mol. Asp provides the necessary charge for the reaction, and DNA glycosylase preferentially cleaves ɛC. This work provides a theoretical basis for the research of excising damaged bases by DNA glycosylase.

QM/MM and MM MD Simulations on the Pyrimidine-Specific Nucleoside Hydrolase: A Comprehensive Understanding of Enzymatic Hydrolysis of Uridine

The pyrimidine-specific nucleoside hydrolase Yeik (CU-NH) from Escherichia coli cleaves the N-glycosidic bond of uridine and cytidine with a 10²-10⁴-fold faster rate than that of purine nucleoside substrates, such as inosine. Such a remarkable substrate specificity and the plausible hydrolytic mechanisms of uridine have been explored by using QM/MM and MM MD simulations. The present calculations show that the relatively stronger hydrogen-bond interactions between uridine and the active-site residues Gln227 and Tyr231 in CU-NH play an important role in enhancing the substrate binding and thus promoting the N-glycosidic bond cleavage, in comparison with inosine. The estimated energy barrier of 30 kcal/mol for the hydrolysis of inosine is much higher than 22 kcal/mol for uridine. Extensive MM MD simulations on the transportation of substrates to the active site of CU-NH indicate that the uridine binding is exothermic by ∼23 kcal/mol, more remarkable than inosine (∼12 kcal/mol). All of these arise from the noncovalent interactions between the substrate and the active site featured in CU-NH, which account for the substrate specificity. Quite differing from other nucleoside hydrolases, here the enzymatic N-glycosidic bond cleavage of uridine is less influenced by its protonation.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated adenine nucleosides: N3 protonation induces base rotation and enhances N-glycosidic bond stability

N3 protonation induces base rotation and stabilizes the syn orientation of the adenine nucleobase of [dAdo+H]⁺ and [Ado+H]⁺via formation of a strong intramolecular N3H⁺⋯O5′ hydrogen-bonding interaction, which in turn influences the mechanisms and energetics for N-glycosidic bond cleavage.

Monitoring potential molecular interactions of adenine with other amino acids using Raman spectroscopy and DFT modeling

We report on the modes of inter-molecular interaction between adenine (Ade) and the amino acids: glycine (Gly), lysine (Lys) and arginine (Arg) using Raman spectroscopy of binary mixtures of adenine and each of the three amino acids at varying molar ratios in the spectral region 1550-550 cm(-1). We focused our attention on certain specific changes in the Raman bands of adenine arising due to its interaction with the amino acids. While the changes are less apparent in the Ade/Gly system, in the Ade/Lys or Ade/Arg systems, significant changes are observed, particularly in the Ade Raman bands that involve the amino group moiety and the N7 and N1 atoms of the purine ring. The ν(N1-C6), ν(N1-C2), δ(C8-H) and δ(N7-C8-N9) vibrations at 1486, 1332, 1253 and 948 cm(-1) show spectral changes on varying the Ade to amino acid molar ratio, the extent of variation being different for the three amino acids. This observation suggests a specific interaction mode between Ade and Lys or Arg, which is due to the hydrogen bonding. The measured spectral changes provide a clear indication that the interaction of Ade depends strongly on the structures of the amino acids, especially their side chains. Density functional theory (DFT) calculations were carried out to elucidate the most probable interaction modes of Ade with the different amino acids.

Standard role for a conserved aspartate or more direct involvement in deglycosylation? An ONIOM and MD investigation of adenine-DNA glycosylase

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.

Theoretical study on the degradation of ADP-ribose polymer catalyzed by poly(ADP-ribose) glycohydrolase

Poly(ADP-ribose) glycohydrolase (PARG) is the only enzyme responsible for the degradation of ADP-ribose polymers. Very recently, the first crystal structure of PARG was reported (Dea Slade, et al., Nature 477 (2011) 616), and a possible SN1-type-like mechanism was proposed. In this work, we present a computational study on the hydrolysis of glycosidic ribose-ribose bond catalyzed by PARG using hybrid density functional theory (DFT) methods. Based on the crystal structure of PARG, three models of the active site were constructed. The calculation results suggest that the degradation of poly(ADP-ribose) follows an SN2 mechanism, and the oxocarbenium expected by Dea Slade is a possible transition state but not an intermediate. The calculated reaction pathway agrees with the proposed mechanism. According to the computational models with different sizes, the roles of key residues are elucidated. Our results may provide useful information for the subsequent experimental and theoretical studies on the structure and functional relationships of PARG.

MUTYH DNA glycosylase: the rationale for removing undamaged bases from the DNA

Maintenance of genetic stability is crucial for all organisms in order to avoid the onset of deleterious diseases such as cancer. One of the many proveniences of DNA base damage in mammalian cells is oxidative stress, arising from a variety of endogenous and exogenous sources, generating highly mutagenic oxidative DNA lesions. One of the best characterized oxidative DNA lesion is 7,8-dihydro-8-oxoguanine (8-oxo-G), which can give rise to base substitution mutations (also known as point mutations). This mutagenicity is due to the miscoding potential of 8-oxo-G that instructs most DNA polymerases (pols) to preferentially insert an Adenine (A) opposite 8-oxo-G instead of the appropriate Cytosine (C). If left unrepaired, such A:8-oxo-G mispairs can give rise to CG→AT transversion mutations. A:8-oxo-G mispairs are proficiently recognized by the MutY glycosylase homologue (MUTYH). MUTYH can remove the mispaired A from an A:8-oxo-G, giving way to the canonical base-excision repair (BER) that ultimately restores undamaged Guanine (G). The importance of this MUTYH-initiated pathway is illustrated by the fact that biallelic mutations in the MUTYH gene are associated with a hereditary colorectal cancer syndrome termed MUTYH-associated polyposis (MAP). In this review, we will focus on MUTYH, from its discovery to the most recent data regarding its cellular roles and interaction partners. We discuss the involvement of the MUTYH protein in the A:8-oxo-G BER pathway acting together with pol λ, the pol that can faithfully incorporate C opposite 8-oxo-G and thus bypass this lesion in a correct manner. We also outline the current knowledge about the regulation of MUTYH itself and the A:8-oxo-G repair pathway by posttranslational modifications (PTM). Finally, to achieve a clearer overview of the literature, we will briefly touch on the rather confusing MUTYH nomenclature. In short, MUTYH is a unique DNA glycosylase that catalyzes the excision of an undamaged base from DNA.

Glycosidic bond cleavage in deoxynucleotides: effects of solvent and the DNA phosphate backbone in the computational model

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.