Mechanistic and conformational flexibility of the covalent linkage formed during β-lyase activity on an AP-site: application to hOgg1

Electrostatic fingerprints of catalytically active amino acids in enzymes

The computed electrostatic and proton transfer properties are studied for 20 enzymes that represent all six major enzyme commission classes and a variety of different folds. The properties of aspartate, glutamate, and lysine residues that have been previously experimentally determined to be catalytically active are reported. The catalytic aspartate and glutamate residues studied here are strongly coupled to at least one other aspartate or glutamate residue and often to multiple other carboxylate residues with intrinsic pK_a differences less than 1 pH unit. Sometimes these catalytic acidic residues are also coupled to a histidine residue, such that the intrinsic pK_a of the acidic residue is higher than that of the histidine. All catalytic lysine residues studied here are strongly coupled to tyrosine or cysteine residues, wherein the intrinsic pK_a of the anion-forming residue is higher than that of the lysine. Some catalytic lysines are also coupled to other lysines with intrinsic pK_a differences within 1 pH unit. Some evidence of the possible types of interactions that facilitate nucleophilicity is discussed. The interactions reported here provide important clues about how side chain functional groups that are weak Brønsted acids or bases for the free amino acid in solution can achieve catalytic potency and become strong acids, bases or nucleophiles in the enzymatic environment.

Molecular dynamics approach to identification of new OGG1 cancer-associated somatic variants with impaired activity

DNA of living cells is always exposed to damaging factors. To counteract the consequences of DNA lesions, cells have evolved several DNA repair systems, among which base excision repair is one of the most important systems. Many currently used antitumor drugs act by damaging DNA, and DNA repair often interferes with chemotherapy and radiotherapy in cancer cells. Tumors are usually extremely genetically heterogeneous, often bearing mutations in DNA repair genes. Thus, knowledge of the functionality of cancer-related variants of proteins involved in DNA damage response and repair is of great interest for personalization of cancer therapy. Although computational methods to predict the variant functionality have attracted much attention, at present, they are mostly based on sequence conservation and make little use of modern capabilities in computational analysis of 3D protein structures. We have used molecular dynamics (MD) to model the structures of 20 clinically observed variants of a DNA repair enzyme, 8-oxoguanine DNA glycosylase. In parallel, we have experimentally characterized the activity, thermostability, and DNA binding in a subset of these mutant proteins. Among the analyzed variants of 8-oxoguanine DNA glycosylase, three (I145M, G202C, and V267M) were significantly functionally impaired and were successfully predicted by MD. Alone or in combination with sequence-based methods, MD may be an important functional prediction tool for cancer-related protein variants of unknown significance.

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.

Structural Insight into the Discrimination between 8-Oxoguanine Glycosidic Conformers by DNA Repair Enzymes: A Molecular Dynamics Study of Human Oxoguanine Glycosylase 1 and Formamidopyrimidine-DNA Glycosylase

hOgg1 and FPG are the primary DNA repair enzymes responsible for removing the major guanine (G) oxidative product, namely, 7,8-dihydro-8-oxoguanine (OG), in humans and bacteria, respectively. While natural G adopts the anti conformation and forms a Watson-Crick pair with cytosine (C), OG can also adopt the syn conformation and form a Hoogsteen pair with adenine (A). hOgg1 removes OG paired with C but is inactive toward the OG:A pair. In contrast, FPG removes OG from OG:C pairs and also exhibits appreciable (although diminished) activity toward OG:A pairs. As a first step toward understanding this difference in activity, we have employed molecular dynamics simulations to examine how the anti and syn conformers of OG are accommodated in the hOgg1 and FPG active sites. When anti-OG is bound, hOgg1 active site residues are properly aligned to initiate catalytic base departure, while geometrical parameters required for the catalytic reaction are not conserved for syn-OG. On the other hand, the FPG catalytic residues are suitably aligned for both OG conformers, with anti-OG being more favorably bound. Thus, our data suggests that the differential ability of hOgg1 and FPG to accommodate the anti- and syn-OG glycosidic conformations is an important factor that contributes to the relative experimental excision rates. Nevertheless, the positions of the nucleophiles with respect to the lesion in the active sites suggest that the reactant complex is poised to initiate catalysis through a similar mechanism for both repair enzymes and supports a recently proposed mechanism in which sugar-ring opening precedes nucleoside deglycosylation.

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^-···H₂O 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 H₂O 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.

Quantum mechanical study of the β- and δ-lyase reactions during the base excision repair process: application to FPG

The β- and δ-elimination reactions catalyzed by FPG during the base excision repair of 8-oxoguanine are intrinsically different.

Stimulation of N--glycoside transfer in deoxythymidine glycol: mechanism of the initial step in base excision repair

Thymine glycol (Tg), a toxic oxidative DNA lesion, is preferentially removed by endonuclease III (Endo III). To investigate the glycosylase activity of Endo III, the N--glycoside transfer mechanism in deoxythymidine glycol (dTg) is examined in this theoretical study based on the BHandHLYP/6-311++G(d,p) level of theory. Two controversial mechanisms were characterized, i.e., the displacement and endocyclic mechanisms. For each mechanism, three types of reaction models were established, including the direct reaction, local microhydration and protonated models. The calculated results indicate that (i) all three reaction models favor the displacement mechanism more than the endocyclic mechanism; (ii) the local microhydration model allows for discrete proton transfer and contributes to the reduction of activation energies, nevertheless, large activation energies are still involved; (iii) the O4'-protonated endocyclic model can efficiently promote the nucleophilic attack of lysine residue and an amino acid residue other than the nucleophilic lysine should be responsible for the opening of the sugar ring; (iv) the O2-protonated displacement model facilitates the leaving group (Tg) stabilization and therefore is the preferred mechanism for the N--glycoside transfer of dTg, whose activation energy of 17.7 kcal mol⁻¹ is in good agreement with the experimental estimate of 19.0 kcal mol⁻¹. As a result, the protonation of nucleobase plays a significant role in predicting the preferred glycosylase mechanism. Our findings can propose appropriate mechanisms for future large-scale enzymatic modeling of Endo III and provide more fundamental information about the important residues that may be included in the enzyme-catalyzed reactions.

The activation of N-glycosidic bond cleavage performed by base-excision repair enzyme hOGG1; theoretical study of the role of Lys 249 residue in activation of G, OxoG and FapyG

NLMOs of lone-pair electrons at N9 nitrogen and Fukui indexesf²of N9.

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.