Combining structural and bioinformatics methods for the analysis of functionally important residues in DNA glycosylases

Molecular dynamics simulation of the opposite-base preference and interactions in the active site of formamidopyrimidine-DNA glycosylase

Background

Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available.

Results

Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition.

Conclusion

Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0075-y) contains supplementary material, which is available to authorized users.

A discontinuous DNA glycosylase domain in a family of enzymes that excise 5-methylcytosine

DNA cytosine methylation (5-meC) is a widespread epigenetic mark associated to gene silencing. In plants, DEMETER-LIKE (DML) proteins typified by Arabidopsis REPRESSOR OF SILENCING 1 (ROS1) initiate active DNA demethylation by catalyzing 5-meC excision. DML proteins belong to the HhH-GPD superfamily, the largest and most functionally diverse group of DNA glycosylases, but the molecular properties that underlie their capacity to specifically recognize and excise 5-meC are largely unknown. We have found that sequence similarity to HhH-GPD enzymes in DML proteins is actually distributed over two non-contiguous segments connected by a predicted disordered region. We used homology-based modeling to locate candidate residues important for ROS1 function in both segments, and tested our predictions by site-specific mutagenesis. We found that amino acids T606 and D611 are essential for ROS1 DNA glycosylase activity, whereas mutations in either of two aromatic residues (F589 and Y1028) reverse the characteristic ROS1 preference for 5-meC over T. We also found evidence suggesting that ROS1 uses Q607 to flip out 5-meC, while the contiguous N608 residue contributes to sequence-context specificity. In addition to providing novel insights into the molecular basis of 5-meC excision, our results reveal that ROS1 and its DML homologs possess a discontinuous catalytic domain that is unprecedented among known DNA glycosylases.

A "moving metal mechanism" for substrate cleavage by the DNA repair endonuclease APE-1

Apurinic/apyrimidinic endonuclease (APE-1) is essential for base excision repair (BER) of damaged DNA. Here molecular dynamics (MD) simulations of APE1 complexed with cleaved and uncleaved damaged DNA were used to determine the role and position of the metal ion(s) in the active site before and after DNA cleavage. The simulations started from an energy minimized wild-type structure of the metal-free APE1/damaged-DNA complex (1DE8). A grid search with one Mg2+ ion located two low energy clusters of Mg2+ consistent with the experimentally determined metal ion positions. At the start of the longer MD simulations, Mg2+ ions were placed at different positions as seen in the crystal structures and the movement of the ion was followed over the course of the trajectory. Our analysis suggests a "moving metal mechanism" in which one Mg2+ ion moves from the B- (more buried) to the A-site during substrate cleavage. The anticipated inversion of the phosphate oxygens occurs during the in-line cleavage reaction. Experimental results, which show competition between Ca2+ and Mg2+ for catalyzing the reaction, and high concentrations of Mg2+ are inhibitory, indicate that both sites cannot be simultaneously occupied for maximal activity.

Catalytic mechanism of Escherichia coli endonuclease VIII: roles of the intercalation loop and the zinc finger

Endonuclease VIII (Nei) excises oxidatively damaged pyrimidines from DNA and shares structural and functional homology with formamidopyrimidine-DNA glycosylase. Although the structure of Escherichia coli Nei is solved [Zharkov et al. (2002) EMBO J. 21, 789-800], the functions of many of its amino acid residues involved in catalysis and substrate specificity are not known. We constructed a series of Nei mutants that interfere with eversion of the damaged base from the helix (QLY69-71AAA, DeltaQLY69-71) or perturb the conserved zinc finger (R171A, Q261A). Steady-state kinetics were measured with these mutant enzymes using substrates containing 5,6-dihydrouracil, two enantiomers of thymine glycol, 8-oxo-7,8-dihydroguanine, and an abasic site positioned opposite each of the four canonical DNA bases. To some extent, all Nei mutants were deficient in processing damaged DNA, with mutations in the zinc finger generally having a more profound effect. Wild-type Nei showed prominent opposite-base specificity (G > C approximately = T > A) when the lesion was 5,6-dihydrouracil or cis-(5S,6R)-thymine glycol but not for other lesions tested. Mutations in the Q69-Y71 loop eliminated this effect. Only wild-type Nei and Nei-Q261A mutants could be reductively cross-linked to damaged base-containing DNA. Experiments involving trapping with NaBH4 and the kinetics of DNA cleavage catalyzed by Nei-Q261A suggested that this mutant was deficient in regenerating free enzyme from the Nei-DNA covalent complex formed during the reaction. We conclude that the opposite-base specificity of Nei is primarily governed by residues in the Q69-Y71 loop and that both this loop and the zinc finger contribute significantly to the substrate specificity of Nei.

The DNA trackwalkers: principles of lesion search and recognition by DNA glycosylases

DNA glycosylases, the pivotal enzymes in base excision repair, are faced with the difficult task of recognizing their substrates in a large excess of unmodified DNA. We present here a kinetic analysis of DNA glycosylase substrate specificity, based on the probability of error. This novel approach to this subject explains many features of DNA surveillance and catalysis of lesion excision by DNA glycosylases. This approach also is applicable to the general issue of substrate specificity. We discuss determinants of substrate specificity in damaged DNA and in the enzyme, as well as methods by which these determinants can be identified.

Recognition of cyclonucleoside lesions by the Lactococcus lactis FPG protein

Several purine and pyrimidine cyclonucleosides were found to be not recognized by several Escherichia coli and yeast DNA N-glycosylases. Interestingly, a non covalent complex was observed between the Lactoccocus lactis formamidopyrimidine-DNA glycosylases (Fpg-Ll) and the cyclonucleosides. This may provide new information on the mechanism involved in the activity of the latter enzyme.

Substrate discrimination by formamidopyrimidine-DNA glycosylase: a mutational analysis

Formamidopyrimidine-DNA glycosylase (Fpg) is a primary participant in the repair of 8-oxoguanine, an abundant oxidative DNA lesion. Although the structure of Fpg has been established, amino acid residues that define damage recognition have not been identified. We have combined molecular dynamics and bioinformatics approaches to address this issue. Site-specific mutagenesis coupled with enzyme kinetics was used to test our predictions. On the basis of molecular dynamics simulations, Lys-217 was predicted to interact with the O8 of extrahelical 8-oxoguanine accommodated in the binding pocket. Consistent with our computational studies, mutation of Lys-217 selectively reduced the ability of Fpg to excise 8-oxoguanine from DNA. Dihydrouracil, also a substrate for Fpg, served as a nonspecific control. Other residues involved in damage recognition (His-89, Arg-108, and Arg-109) were identified by combined conservation/structure analysis. Arg-108, which forms two hydrogen bonds with cytosine in Fpg-DNA, is a major determinant of opposite-base specificity. Mutation of this residue reduced excision of 8-oxoguanine from thermally unstable mispairs with guanine or thymine, while excision from the stable cytosine and adenine base pairs was less affected. Mutation of His-89 selectively diminished the rate of excision of 8-oxoguanine, whereas mutation of Arg-109 nearly abolished binding of Fpg to damaged DNA. Taken together, these results suggest that His-89 and Arg-109 form part of a reading head, a structural feature used by the enzyme to scan DNA for damage. His-89 and Lys-217 help determine the specificity of Fpg in recognizing the oxidatively damaged base, while Arg-108 provides specificity for bases positioned opposite the lesion.

Structural characterization of the Fpg family of DNA glycosylases

Until recently, the Fpg family was the only major group of DNA glycosylases for which no structural data existed. Prototypical members of this family, found in eukaryotes as well as prokaryotes, have now been crystallized as free proteins and as complexes with DNA. In this review, we analyze the available structural information for formamidopyrimidine-DNA glycosylase (Fpg) and endonuclease VIII (Nei). Special emphasis is placed on mechanisms by which these enzymes recognize and selectively excise cognate lesions from oxidatively damaged DNA. The problem of lesion recognition is considered in two parts: how the enzyme efficiently locates a single lesion embedded in a vast excess of DNA; and how the lesion is accommodated in a pocket near the active site of the enzyme. Although all crystal structures reported to date for the Fpg family lack the damaged base, functionally important residues that participate in DNA binding and enzyme catalysis have been clearly identified and other residues, responsible for substrate specificity, have been inferred.

Structure of formamidopyrimidine-DNA glycosylase covalently complexed to DNA

Formamidopyrimidine-DNA glycosylase (Fpg) is a DNA repair enzyme that excises oxidized purines from damaged DNA. The Schiff base intermediate formed during this reaction between Escherichia coli Fpg and DNA was trapped by reduction with sodium borohydride, and the structure of the resulting covalently cross-linked complex was determined at a 2.1-A resolution. Fpg is a bilobal protein with a wide, positively charged DNA-binding groove. It possesses a conserved zinc finger and a helix-two turn-helix motif that participate in DNA binding. The absolutely conserved residues Lys-56, His-70, Asn-168, and Arg-258 form hydrogen bonds to the phosphodiester backbone of DNA, which is sharply kinked at the lesion site. Residues Met-73, Arg-109, and Phe-110 are inserted into the DNA helix, filling the void created by nucleotide eversion. A deep hydrophobic pocket in the active site is positioned to accommodate an everted base. Structural analysis of the Fpg-DNA complex reveals essential features of damage recognition and the catalytic mechanism of Fpg.