Insights into the Direct Oxidative Repair of Etheno Lesions: MD and QM/MM Study on the Substrate Scope of ALKBH2 and AlkB.

E. coli AlkB and human ALKBH2 belong to the AlkB family enzymes, which contain several α-ketoglutarate (α-KG)/Fe(II)-dependent dioxygenases that repair alkylated DNA. Specifically, the AlkB enzymes catalyze decarboxylation of α-KG to generate a high-valent Fe(IV)-oxo species that oxidizes alkyl groups on DNA adducts. AlkB and ALKBH2 have been reported to differentially repair select etheno adducts, with preferences for 1,N6-ethenoadenine (1,N6-εA) and 3,N4-ethenocytosine (3,N4-εC) over 1,N2-ethenoguanine (1,N2-εG). However, N2,3-ethenoguanine (N2,3-εG), the most common etheno adduct, is not repaired by the AlkB enzymes. Unfortunately, a structural understanding of the differential activity of E. coli AlkB and human ALKBH2 is lacking due to challenges acquiring atomistic details for a range of substrates using experiments. This study uses both molecular dynamics (MD) simulations and ONIOM(QM:MM) calculations to determine how the active site changes upon binding each etheno adduct and characterizes the corresponding catalytic impacts. Our data reveal that the preferred etheno substrates (1,N6-εA and 3,N4-εC) form favorable interactions with catalytic residues that situate the lesion near the Fe(IV)-oxo species and permit efficient oxidation. In contrast, although the damage remains correctly aligned with respect to the Fe(IV)-oxo moiety, repair of 1,N2-εG is mitigated by increased solvation of the active site and a larger distance between Fe(IV)-oxo and the aberrant carbons. Binding of non-substrate N2,3-εG in the active site disrupts key DNA-enzyme interactions, and positions the aberrant carbon atoms even further from the Fe(IV)-oxo species, leading to prohibitively high barriers for oxidative catalysis. Overall, our calculations provide the first structural insight required to rationalize the experimentally-reported substrate specificities of AlkB and ALKBH2 and thereby highlight the roles of several active site residues in the repair of etheno adducts that directly correlates with available experimental data. These proposed catalytic strategies can likely be generalized to other α-KG/Fe(II)-dependent dioxygenases that play similar critical biological roles, including epigenetic and post-translational regulation.

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.

DNA repair enzymes ALKBH2, ALKBH3, and AlkB oxidize 5-methylcytosine to 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine in vitro.

5-Methylcytosine (5mC) in DNA CpG islands is an important epigenetic biomarker for mammalian gene regulation. It is oxidized to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) by the ten-eleven translocation (TET) family enzymes, which are α-ketoglutarate (α-KG)/Fe(II)-dependent dioxygenases. In this work, we demonstrate that the epigenetic marker 5mC is modified to 5hmC, 5fC, and 5caC in vitro by another class of α-KG/Fe(II)-dependent proteins-the DNA repair enzymes in the AlkB family, which include ALKBH2, ALKBH3 in huamn and AlkB in Escherichia coli. Theoretical calculations indicate that these enzymes may bind 5mC in the syn-conformation, placing the methyl group comparable to 3-methylcytosine, the prototypic substrate of AlkB. This is the first demonstration of the AlkB proteins to oxidize a methyl group attached to carbon, instead of nitrogen, on a DNA base. These observations suggest a broader role in epigenetics for these DNA repair proteins.

Structural explanation for the tunable substrate specificity of an E. coli nucleoside hydrolase: insights from molecular dynamics simulations.

Parasitic protozoa rely on nucleoside hydrolases that play key roles in the purine salvage pathway by catalyzing the hydrolytic cleavage of the N-glycosidic bond that connects nucleobases to ribose sugars. Cytidine-uridine nucleoside hydrolase (CU-NH) is generally specific toward pyrimidine nucleosides; however, previous work has shown that replacing two active site residues with Tyr, specifically the Thr223Tyr and Gln227Tyr mutations, allows CU-NH to process inosine. The current study uses molecular dynamics (MD) simulations to gain atomic-level insight into the activity of wild-type and mutant E. coli CU-NH toward inosine. By examining systems that differ in the identity and protonation states of active site catalytic residues, key enzyme-substrate interactions that dictate the substrate specificity of CU-NH are identified. Regardless of the wild-type or mutant CU-NH considered, our calculations suggest that inosine binding is facilitated by interactions of the ribose moiety with active site residues and Ca2+, and π-interactions between two His residues (His82 and His239) and the nucleobase. However, the lack of observed activity toward inosine for wild-type CU-NH is explained by no residue being correctly aligned to stabilize the departing nucleobase. In contrast, a hydrogen-bonding network between hypoxanthine and a newly identified general acid (Asp15) is present when the two Tyr mutations are engineered into the active site. Investigation of the single CU-NH mutants reveals that this hydrogen-bonding network is only maintained when both Tyr mutations are present due to a π-interaction between the residues. These results rationalize previous experiments that show the single Tyr mutants are unable to efficiently hydrolyze inosine and explain how the Tyr residues work synergistically in the double mutant to stabilize the nucleobase leaving group during hydrolysis. Overall, our simulations provide a structural explanation for the substrate specificity of nucleoside hydrolases, which may be used to rationally develop new treatments for kinetoplastid diseases.

Quantum Chemical Studies of the Structure and Stability of N-Methylated DNA Nucleobase Dimers: Insights into the Mutagenic Base Pairing of Damaged DNA.

DNA is constantly under attack from exogenous and endogenous sources that modify the chemical structure of the nucleobases. A common type of nucleobase damage is N-methylation, which can result in mutagenesis. Nevertheless, these lesions are often repaired by the DNA repair enzyme AlkB, albeit at varying rates. Herein we use density functional theory (B3LYP-D3(BJ)/6-311++G(2df,2p)//B3LYP/6-31G(d,p)) to comprehensively examine the structural and energetic properties of base pairs between seven nucleobase lesions resulting from N-methylation on the Watson-Crick (WC) binding face and each canonical nucleobase. By characterizing 105 stable nucleobase dimers, we provide fundamental details regarding the preferred lesion base pairings. Specifically, we reveal that the flexibility of the methylamino group resulting from methylation of an exocyclic amino substituent allows the 2MeG, 4MeC, and 6MeA lesions to maintain a preference for canonical WC base pairing, which correlates with the experimentally reported lack of mutagenicity for these damage products. In contrast, calculated distortions in key structural parameters and altered binding energies for base pairs involving adducts formed upon methylation of a ring nitrogen (namely, 1MeG, 3MeT, 1MeA, and 3MeC) help rationalize the associated mutagenicity and repair efficiencies. Most importantly, our work provides molecular-level information about the interactions between N-methylated and canonical nucleobases that is critical for future large-scale modeling of damaged DNA and enzyme-DNA complexes that strive to further uncover the mutagenicity and repair propensities of these detrimental lesions.

QM/MM Study of the Reaction Catalyzed by Alkyladenine DNA Glycosylase: Examination of the Substrate Specificity of a DNA Repair Enzyme.

Human alkyladenine DNA glycosylase (AAG) functions as part of the base excision repair pathway to excise structurally diverse oxidized and alkylated DNA purines. Specifically, AAG uses a water molecule activated by a general base and a nonspecific active site lined with aromatic residues to cleave the N-glycosidic bond. Despite broad substrate specificity, AAG does not target the natural purines (adenine (A) and guanine (G)). Using the ONIOM(QM:MM) methodology, we provide fundamental atomic level details of AAG bound to DNA-containing a neutral substrate (hypoxanthine (Hx)), a nonsubstrate (G), or a cationic substrate (7-methylguanine (7MeG)) and probe changes in the reaction pathway that occur when AAG targets different nucleotides. We reveal that subtle differences in protein-DNA contacts upon binding different substrates within the flexible AAG active site can significantly affect the deglycosylation reaction. Notably, we predict that AAG excises Hx in a concerted mechanism that is facilitated through correct alignment of the (E125) general base due to hydrogen bonding with a neighboring aromatic amino acid (Y127). Hx departure is further stabilized by π-π interactions with aromatic amino acids and hydrogen bonds with active site water. Despite possessing a similar structure to Hx, G is not excised since the additional exocyclic amino group leads to misalignment of the general base due to disruption of the key E125-Y127 hydrogen bond, the catalytically unfavorable placement of water within the active site, and weakened π-contacts between aromatic amino acids and the nucleobase. In contrast, cationic 7MeG does not occupy the same position within the AAG active site as G due to steric clashes with the additional N7 methyl group, which results in the correct alignment of the general base and permits nucleobase excision as observed for neutral Hx. Overall, our structural data rationalizes the observed substrate specificity of AAG and contributes to our fundamental understanding of enzymes with flexible active sites and broad substrate specificities.

Hydrolytic Glycosidic Bond Cleavage in RNA Nucleosides: Effects of the 2'-Hydroxy Group and Acid-Base Catalysis.

Despite the inherent stability of glycosidic linkages in nucleic acids that connect the nucleobases to sugar-phosphate backbones, cleavage of these bonds is often essential for organism survival. The current study uses DFT (B3LYP) to provide a fundamental understanding of the hydrolytic deglycosylation of the natural RNA nucleosides (A, C, G, and U), offers a comparison to DNA hydrolysis, and examines the effects of acid, base, or simultaneous acid-base catalysis on RNA deglycosylation. By initially examining HCOO-···H2O mediated deglycosylation, the barriers for RNA hydrolysis were determined to be 30-38 kJ mol-1 higher than the corresponding DNA barriers, indicating that the 2'-OH group stabilizes the glycosidic bond. Although the presence of HCOO- as the base (i.e., to activate the water nucleophile) reduces the barrier for uncatalyzed RNA hydrolysis (i.e., unactivated H2O nucleophile) by ∼15-20 kJ mol-1, the extreme of base catalysis as modeled using a fully deprotonated water molecule (i.e., OH- nucleophile) decreases the uncatalyzed barriers by up to 65 kJ mol-1. Acid catalysis was subsequently examined by selectively protonating the hydrogen-bond acceptor sites of the RNA nucleobases, which results in an up to ∼80 kJ mol-1 barrier reduction relative to the corresponding uncatalyzed pathway. Interestingly, the nucleobase proton acceptor sites that result in the greatest barrier reductions match sites typically targeted in enzyme-catalyzed reactions. Nevertheless, simultaneous acid and base catalysis is the most beneficial way to enhance the reactivity of the glycosidic bonds in RNA, with the individual effects of each catalytic approach being weakened, additive, or synergistic depending on the strength of the base (i.e., degree of water nucleophile activation), the nucleobase, and the hydrogen-bonding acceptor site on the nucleobase. Together, the current contribution provides a greater understanding of the reactivity of the glycosidic bond in natural RNA nucleosides, and has fundamental implications for the function of RNA-targeting enzymes.

Evaluating the Substrate Selectivity of Alkyladenine DNA Glycosylase: The Synergistic Interplay of Active Site Flexibility and Water Reorganization.

Human alkyladenine DNA glycosylase (AAG) functions as part of the base excision repair (BER) pathway by cleaving the N-glycosidic bond that connects nucleobases to the sugar-phosphate backbone in DNA. AAG targets a range of structurally diverse purine lesions using nonspecific DNA-protein π-π interactions. Nevertheless, the enzyme discriminates against the natural purines and is inhibited by pyrimidine lesions. This study uses molecular dynamics simulations and seven different neutral or charged substrates, inhibitors, or canonical purines to probe how the bound nucleotide affects the conformation of the AAG active site, and the role of active site residues in dictating substrate selectivity. The neutral substrates form a common DNA-protein hydrogen bond, which results in a consistent active site conformation that maximizes π-π interactions between the aromatic residues and the nucleobase required for catalysis. Nevertheless, subtle differences in DNA-enzyme contacts for different neutral substrates explain observed differential catalytic efficiencies. In contrast, the exocyclic amino groups of the natural purines clash with active site residues, which leads to catalytically incompetent DNA-enzyme complexes due to significant reorganization of active site water. Specifically, water resides between the A nucleobase and the active site aromatic amino acids required for catalysis, while a shift in the position of the general base (E125) repositions (potentially nucleophilic) water away from G. Despite sharing common amino groups, the methyl substituents in cationic purine lesions (3MeA and 7MeG) exhibit repulsion with active site residues, which repositions the damaged bases in the active site in a manner that promotes their excision. Overall, we provide a structural explanation for the diverse yet discriminatory substrate selectivity of AAG and rationalize key kinetic data available for the enzyme. Specifically, our results highlight the complex interplay of many different DNA-protein interactions used by AAG to facilitate BER, as well as the crucial role of the general base and water (nucleophile) positioning. The insights gained from our work will aid the understanding of the function of other enzymes that use flexible active sites to exhibit diverse substrate specificity.

Glycosidic Bond Cleavage in DNA Nucleosides: Effect of Nucleobase Damage and Activation on the Mechanism and Barrier.

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.

Hydrolysis of the damaged deoxythymidine glycol nucleoside and comparison to canonical DNA.

Genomic integrity is continually under attack by both endogenous and exogenous sources. One of the most common forms of damage is oxidation of the thymine nucleobase to form (5R,6S)-dihydroxy-5,6-dihydro-thymine (thymine glycol or Tg), which stops DNA polymerases and is thus cytotoxic. Thymine glycol damage is repaired through a variety of mechanisms, including the multi-step base excision repair (BER) pathway. In the first BER step, the glycosidic bond of the dTg nucleotide is hydrolyzed by a DNA glycosylase. In order to understand the catalytic effect of the glycosylases, the corresponding uncatalyzed mechanisms and barriers are required, as well as an appreciation of the relative reactivity of the glycosidic bond with respect to the corresponding canonical nucleoside. To this end, the PCM-B3LYP/6-31+G(d) reaction potential energy surfaces (PES) for deoxythymidine (dT) and dTg hydrolysis are characterized in the present study using solvent-phase optimizations and a model containing three explicit water molecules. The surfaces are comparable to those generated using functionals that account for dispersion interactions (B3LYP-D3 and M06-2X). Mapping the PES as a function of the glycosidic bond length and nucleophile-sugar distance reveals a synchronous S(N)2 mechanism as the lowest energy pathway for damaged dTg hydrolysis, which contrasts the preferred dissociative S(N)1 mechanism isolated for the deglycosylation of natural dT. As proposed for other enzymes, the difference in excision pathway may at least in part help the enzyme selectively target the damaged base and discriminate against the natural counterpart. Interestingly, the barrier to dTg deglycosylation (ΔG(‡) = 138.0 kJ mol(-1)) is much higher than for dT deglycosylation (ΔG(‡) = 112.7 kJ mol(-1)), which supports the stability of this lesion and clarifies the catalytic feat presented to DNA repair enzymes that remove this detrimental damage from the genome. Although nucleotide excision repair (NER) typically targets bulky DNA lesions, the large calculated barrier for dTg deglycosylation rationalizes why the NER mechanism also excises this non-bulky lesion from cellular 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.