QM/MM Study of the Uracil DNA Glycosylase Reaction Mechanism: A Competition between Asp145 and His148

Uridine Embedded within DNA is Repaired by Uracil DNA Glycosylase via a Mechanism Distinct from That of Ribonuclease H2

Uridine (rU) and 2'-deoxyuridine (dU) are common DNA lesions. dU is repaired through a base excision repair (BER) pathway initiated by uracil DNA glycosylase (UDG), while rU is typically removed from DNA via ribonucleotide excision repair, mediated by RNase H2. In this study, we report that rU is also repaired through the UDG-mediated BER pathway. We found that UDG catalyzes the removal of uracil from rU embedded in DNA, but exhibits no activity toward rU in RNA. Biochemical and crystallographic analyses revealed that the 2'-OH group of rU is effectively accommodated by UDG and directly participates in catalyzing the hydrolysis of the N-glycosidic bond. The abasic site product generated upon removal of uracil from rU by UDG is further processed by downstream BER enzymes to restore undamaged DNA. Our findings suggest that UDG-initiated BER constitutes a previously unrecognized pathway for the repair of rU-specific ribonucleotides. Additionally, we developed a method for selectively quantifying rU content in DNA. Using this method, we determined that rU repair by UDG is not a major pathway in human cells. This discovery expands our understanding of the diverse biological functions of UDG, and inspires further investigation to determine the role of its rU-removal in cells.

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.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

Structural Insights into the Mechanism of Base Excision by MBD4

DNA glycosylases remove damaged or modified nucleobases by cleaving the N-glycosyl bond and the correct nucleotide is restored through subsequent base excision repair. In addition to excising threatening lesions, DNA glycosylases contribute to epigenetic regulation by mediating DNA demethylation and perform other important functions. However, the catalytic mechanism remains poorly defined for many glycosylases, including MBD4 (methyl-CpG binding domain IV), a member of the helix-hairpin-helix (HhH) superfamily. MBD4 excises thymine from G·T mispairs, suppressing mutations caused by deamination of 5-methylcytosine, and it removes uracil and modified uracils (e.g., 5-hydroxymethyluracil) mispaired with guanine. To investigate the mechanism of MBD4 we solved high-resolution structures of enzyme-DNA complexes at three stages of catalysis. Using a non-cleavable substrate analog, 2'-deoxy-pseudouridine, we determined the first structure of an enzyme-substrate complex for wild-type MBD4, which confirms interactions that mediate lesion recognition and suggests that a catalytic Asp, highly conserved in HhH enzymes, binds the putative nucleophilic water molecule and stabilizes the transition state. Observation that mutating the Asp (to Gly) reduces activity by 2700-fold indicates an important role in catalysis, but probably not one as the nucleophile in a double-displacement reaction, as previously suggested. Consistent with direct-displacement hydrolysis, a structure of the enzyme-product complex indicates a reaction leading to inversion of configuration. A structure with DNA containing 1-azadeoxyribose models a potential oxacarbenium-ion intermediate and suggests the Asp could facilitate migration of the electrophile towards the nucleophilic water. Finally, the structures provide detailed snapshots of the HhH motif, informing how these ubiquitous metal-binding elements mediate DNA binding.

Multiscale Modeling of Wobble to Watson-Crick-Like Guanine-Uracil Tautomerization Pathways in RNA

Energetically unfavorable Watson–Crick (WC)-like tautomeric forms of nucleobases are known to introduce spontaneous mutations, and contribute to replication, transcription, and translation errors. Recent NMR relaxation dispersion techniques were able to show that wobble (w) G•U mispair exists in equilibrium with the short-lived, low-population WC-like enolic tautomers. Presently, we have investigated the wG•U → WC-like enolic reaction pathway using various theoretical methods: quantum mechanics (QM), molecular dynamics (MD), and combined quantum mechanics/molecular mechanics (QM/MM). The previous studies on QM gas phase calculations were inconsistent with experimental data. We have also explored the environmental effects on the reaction energies by adding explicit water. While the QM-profile clearly becomes endoergic in the presence of water, the QM/MM-profile remains consistently endoergic in the presence and absence of water. Hence, by including microsolvation and QM/MM calculations, the experimental data can be explained. For the G•Uenol→ Genol•U pathway, the latter appears to be energetically more favorable throughout all computational models. This study can be considered as a benchmark of various computational models of wG•U to WC-like tautomerization pathways with and without the environmental effects, and may contribute on further studies of other mispairs as well.

DFT Study on the Deglycosylation of Methylated, Oxidized, and Canonical Pyrimidine Nucleosides in Water: Implications for Epigenetic Regulation and DNA Repair

Density functional theory (B3LYP) was used to characterize the kinetics and thermodynamics of the (nonenzymatic) deglycosylation in water for a variety of 2'-deoxycytidine (dC) and 2'-deoxyuridine (dU) nucleoside derivatives that differ in methylation and subsequent oxidation of the C5 substituent. A range of computational models are considered that combine implicit and explicit solvation of the nucleophile and nucleobase. Regardless of the model implemented, our calculations reveal that the glycosidic bond in dC is inherently more stable than that in dU. Furthermore, C5 methylation of either pyrimidine and subsequent oxidation of the methyl group yield overall small changes to the Gibbs reaction energy profiles and thereby preserve lower deglycosylation barriers for the dC compared to those for the dU nucleoside derivatives. However, hydrolytic deglycosylation becomes significantly more energetically favorable when 5-methyl-dC (5m-dC) undergoes two or three rounds of oxidation, with the Gibbs energy barrier decreasing and the reaction becoming more exergonic by up to 40 kJ/mol. In fact, two or three oxidation reactions from 5m-dC result in a deglycosylation barrier similar to that for dU, as well as those for the associated C5-methylated (2'-deoxythymidine) and oxidized (5-hydroxymethyl-dU) derivatives. These predicted trends in the inherent deglycosylation energetics in water directly correlate with the previously reported activity of thymine DNA glycosylase (TDG), which cleaves the glycosidic bond in select dC nucleosides as part of epigenetic regulation and in dU variants as part of DNA repair. Thus, our data suggests that fundamental differences in the intrinsic reactivity of the pyrimidine nucleosides help regulate the function of human enzymes that maintain cellular integrity.