Implications of Alternative Substrate Binding Modes for Catalysis by Uracil-DNA Glycosylase: An Apparent Discrepancy Resolved
- 28 October 2006
- journal article
- research article
- Published by American Chemical Society (ACS) in Biochemistry
- Vol. 45 (46) , 13687-13696
- https://doi.org/10.1021/bi061061y
Abstract
A theoretical study showed that the base excision repair enzyme uracil-DNA glycosylase (UDG) exploits electrostatic interactions with backbone phosphate groups in the substrate for catalysis. Although experiments performed to test the calculated results confirmed the predicted importance of the -2, -1, and +1 phosphate groups, there was an apparent disagreement with regard to the +2 phosphate group. The calculations indicated that it made an important contribution, while experimentally, the effect of its deletion or neutralization was small. The +2 phosphate group interacts directly with an active site histidine (H148 in humans) in the crystal structure of UDG in complex with double-stranded (ds) DNA. We previously calculated that H148 has a strong anticatalytic effect due to its protonation, and here we use alchemical free energy simulations to estimate its site-specific pKa. We find that it is positively charged over the entire experimental pH range (4-10), so its deprotonation cannot compensate for deletion or neutralization of the +2 phosphate group. The free energy simulations are facilitated by an efficient charge-scaling procedure that allows quantitative correction for the implicit treatment of solvent far from the active site; improvements are made to that method to account carefully for differences in the truncation of electrostatic interactions in the contributing molecular-mechanical and continuum-electrostatic approaches. Additional simulations are used to demonstrate that the +2 phosphate group is fully solvent exposed in complexes with single-stranded DNA substrates like those used in the experiments. In contrast, it is well-structured and buried in the dsDNA complex used in the original simulations. Differences in solvent shielding thus account for the apparent lack of an effect observed experimentally upon neutralization or deletion of this group.Keywords
This publication has 16 references indexed in Scilit:
- Dynamic coupling between coordinates in a model for biomolecular isomerizationThe Journal of Chemical Physics, 2006
- Comment on "Uracil DNA Glycosylase Activity Is Dispensable for Immunoglobulin Class Switch"Science, 2004
- Response to Comment on "Uracil DNA Glycosylase Activity Is Dispensable for Immunoglobulin Class Switch"Science, 2004
- Influence of the treatment of electrostatic interactions on the results of free energy calculations of dipolar systemsThe Journal of Chemical Physics, 2000
- All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target dataJournal of Computational Chemistry, 2000
- Continuum Treatment of Long-Range Interactions in Free Energy Calculations. Application to Protein−Ligand Binding.The Journal of Physical Chemistry B, 1997
- Free energy calculations: Applications to chemical and biochemical phenomenaChemical Reviews, 1993
- THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The methodJournal of Computational Chemistry, 1992
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983