Figures
Abstract
Human 8-oxoguanine DNA glycosylase (hOGG1) is a key enzyme responsible for initiating the base excision repair of 7,8-dihydro-8-oxoguanosine (oxoG). In this study a thermodynamic analysis of the interaction of hOGG1 with specific and non-specific DNA-substrates is performed based on stopped-flow kinetic data. The standard Gibbs energies, enthalpies and entropies of specific stages of the repair process were determined via kinetic measurements over a temperature range using the van’t Hoff approach. The three steps which are accompanied with changes in the DNA conformations were detected via 2-aminopurine fluorescence in the process of binding and recognition of damaged oxoG base by hOGG1. The thermodynamic analysis has demonstrated that the initial step of the DNA substrates binding is mainly governed by energy due to favorable interactions in the process of formation of the recognition contacts, which results in negative enthalpy change, as well as due to partial desolvation of the surface between the DNA and enzyme, which results in positive entropy change. Discrimination of non-specific G base versus specific oxoG base is occurring in the second step of the oxoG-substrate binding. This step requires energy consumption which is compensated by the positive entropy contribution. The third binding step is the final adjustment of the enzyme/substrate complex to achieve the catalytically competent state which is characterized by large endothermicity compensated by a significant increase of entropy originated from the dehydration of the DNA grooves.
Citation: Kuznetsov NA, Kuznetsova AA, Vorobjev YN, Krasnoperov LN, Fedorova OS (2014) Thermodynamics of the DNA Damage Repair Steps of Human 8-Oxoguanine DNA Glycosylase. PLoS ONE 9(6): e98495. https://doi.org/10.1371/journal.pone.0098495
Editor: Fenfei Leng, Florida International University, United States of America
Received: February 3, 2014; Accepted: May 3, 2014; Published: June 9, 2014
Copyright: © 2014 Kuznetsov et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Program of the Russian Academy of Sciences “Molecular and Cell Biology” (6.11); President Grants (nos. SS-1205.2014.4 and SP-4012.2013.4); the Grant from Russian Scientific Foundation (14-14-00063); O.S.F., A.A.K., N.A.K. and Y.N.V. are supported by Russian Foundation for Basic Research (nos. 13-04-00013, 14-04-31174, 14-04-00531 and 12-04-00135, respectively). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
DNA is continuously damaged by reactive oxygen species (ROS) generated by UV light, ionizing radiation and during metabolism [1]–[5]. Among the various products of oxidative stress, 7,8-dihydro-8-oxoguanosine (oxoG) is the most commonly found. It is a pre-mutagenic DNA lesion since oxoG is able to mispair with adenine, thus generating G/C to T/A transversion mutations [6]. As is the case of most oxidized bases, oxoG is primarily removed in the base excision repair (BER) pathway. This pathway is initiated by the recognition of the modified bases by specific DNA glycosylases. In human cells, 8-oxoguanine DNA glycosylase hOGG1 accomplishes the excision of oxoG residues. HOGG1 belongs to the HhH-GPD superfamily of DNA glycosylases containing conserved structural helix-hairpin-helix and GPD motifs [7]–[10]. HOGG1 is a bifunctional enzyme, possessing DNA glycosylase activity (hydrolysis of the N-glycosidic bond of the damaged nucleotide resulting in formation of the abasic product, Fig. 1) and AP lyase activity (elimination of the 3′-phosphate, often referred to as β-elimination resulting in formation of the nicked product, Fig. 1). In contrast to several other bifunctional DNA glycosylases, the AP lyase activity of hOGG1 is much weaker than its glycosylase activity [11]–[14].
Step 1: oxoG-base removal and formation of abasic product, step 2: β-elimination of 3′-phosphate resulting in formation of the nicked product.
The 3D structures of free hOGG1 and its covalent complexes with DNA have been determined by X-ray crystallography [9], [10], [15]–[19]. The structural data indicate that the conformations of both enzyme and DNA undergo significant changes upon binding. The enzyme forms non-specific electrostatic and hydrophobic contacts with the sugar-phosphate backbone. In the specific lesion recognition complex the DNA is kinked by 70°. The oxoG base is fully flipped out from the DNA helix and deeply inserted into the active site of the enzyme (Fig. 2A). Amino acids Arg-154 and Arg-204 inserted in the DNA helix participate in bidentate hydrogen-bonding interactions with the cytosine base, located opposite to the damaged base (Fig. 2B). Additionally, amino acids Asn-149 and Tyr-203 are inserted into the abasic void released due to flipping out of the oxoG residue.
(A) The lesion recognition complex of hOGG1 bound to DNA with the oxoG base inserted in the active site (PDB ID: 1YQR, [17]), (B) specific contacts of enzyme with the cytosine which is located opposite to oxoG base (PDB ID: 1EBM, [9]), (C) the complex of hOGG1 bound to the DNA with the G base inserted in the exo-site (PDB ID: 1YQK, [17]), (D) the complex of hOGG1 bound to the DNA, with the G base inserted in the active site (PDB ID: 3IH7, [19]). The overlays of the trapped catalytic complex hOGG1/oxoG-substrate (red, 1HU0, [16]) with either (E) the complex of hOGG1 bound to DNA with the G base inserted in the exo-site (blue, 1YQK, [17]) or (F) the complex of hOGG1 bound to DNA, with the G base inserted in the active site (green, 3IH7, [19]).
According to the structure of the complex of hOGG1 with a non-specific DNA, the non-damaged G base is held in the semi-everted transition state (Fig. 2C) [17]. In this state the base is located in the exo-site at the distance of ca. 5 Å from the catalytic enzyme pocket. Nevertheless, in a recent paper from Verdine’s laboratory, a different structure of the complex hOGG1 with non-damaged DNA is proposed (Fig. 2D) [19]. It was shown that although a non-damaged G base can be placed into the active site of hOGG1 approximately in the same position as an oxoG base in the lesion recognition complex (Fig. 2E and 2F), the catalytic step does not occur with DNAs containing non-damaged G base [16]. Therefore the question arises - what is the mechanism of recognition of the damaged nucleotide?
To address this issue, in the previous studies we investigated the conformational dynamics of hOGG1 and DNA substrates during the catalytic cycles using fluorescence detection [14], [20]. The enzyme molecule contains 10 fluorescent tryptophan amino acids, located in several parts of the spatial structure. The fluorescence changes of Trp residues in hOGG1 as well as 2-aminopurine (aPu) residues in the DNA substrates have been detected during catalytic cycle at single-turnover conditions. It was shown that the interaction proceeds through multiple conformational transitions, attributed to sequential lesion recognition process, catalytic steps and release of the product from the DNA complex with enzyme [14], [20]. A kinetic scheme (Fig. 3) containing three initial reversible recognition steps, two irreversible chemical steps and final reversible step was proposed to explain the observations. The structural nature of these steps was suggested based on the analysis of several mutant forms of hOGG1 as discussed in detail in [21]. These studies resulted in atomic-level structural and kinetic information to understand the sequence of hOGG1-DNA interactions during the recognition of oxoG base, and provided insights into the roles of individual amino acids in this process. The analysis of the fluorescence kinetic data lead to the conclusion that the first step of the oxoG-DNA repair is a non-specific binding process which consists of C base pulling by Arg-154/Arg-204, pushing of Tyr-203 into DNA helix, electrostatic interaction between Lys-249 and the damaged nucleotide, oxoG flipping out in the exo-site, and the interaction of His-270 with the everted base [21]. During the second step the oxoG base is everted into the active site and is held there due to the π–π stacking with the aromatic ring of Phe-319. After that, His-270 interacts with phosphate group of the everted nucleotide and forms hydrogen bonds between Arg-154, Arg-204, Asn-149 and the C base. The third DNA damage recognition step is most likely associated with the full insertion of Tyr-203 in the duplex and oxoG interaction with Gly-42.
Here E is hOGG1; OG is the oxoG-substrate; (E•OG)n are different enzyme-substrate recognition complexes; E•AP is the complex of E with the abasic site formed in the course of N-glycosylase reaction; E•P is the enzyme-product complex formed in the AP-lyase reaction; P is the final free reaction product; ki and k–i are the rate constants of the forward and backward individual processes.
Understanding of the mechanism of the DNA damaged base recognition by hOGG1 may be improved by elucidation of the thermodynamics of the binding and cleavage steps via determination of their Gibbs free energy (ΔG°), enthalpy (ΔH°) and entropy (ΔS°). Correlation of the thermodynamic data with the proposed sequential mechanism of oxoG-substrate binding and cleavage might result in the quantitative description of the driving forces that govern the structural adjustment of the enzyme and DNA. In the present work, the thermodynamic parameters of the specific steps occurring in the time range 1 ms–1000 s during oxoG base recognition and removal by hOGG1 were obtained via the temperature-dependent fluorescence kinetic data using the van’t Hoff and Eyring equations [22]. The rate constants of the forward and reverse reactions as well as the resulting equilibrium constants of the individual steps in a multi-step reaction mechanism were determined based on the temporal 2-aminopurine (aPu) fluorescence traces recorded at different temperatures. This fluorescent analog of purine bases is sensitive to the DNA conformation and was widely used for the studies of structural changes in DNA (see [23], [24] for details).
Materials and Methods
Protein Expression and Purification
To purify hOGG1 expressed as a recombinant protein, 1 liter of Escherichia coli strain Rosetta II(DE3) (Invitrogen, France) carrying the pET28c(+) hOGG1 construct was grown in LB broth with 25 µg/ml kanamycin at 25°C until A600 = 0.6–0.7 and induced overnight with 0.2 mM isopropyl-β-D-thiogalacto-pyranoside. The method of the isolation of the wild-type hOGG1 was described previously [14], [21]. The protein concentration was measured according to the Bradford method [25]; the stock solution was stored at −20°C.
Oligodeoxynucleotides
The sequences of ODNs used in this work are listed in Table 1. The ODNs were synthesized by standard phosphoramidite methods using an ASM-700 synthesizer (BIOSSET Ltd., Novosibirsk, Russia) in the Laboratory of Bionanotechnology of ICBFM from phosphoramidites purchased from Glen Research (Sterling, VA). Synthetic oligonucleotides were unloaded from the solid support with ammonium hydroxide according to the manufacturer protocols. Deprotected oligonucleotides were purified by HPLC. The purity of ODNs exceeded 98% as estimated by electrophoresis in 20% denaturing PAGE after staining with the Stains-All dye (Sigma-Aldrich). Concentrations of oligonucleotides were determined from their absorbances at 260 nm. The ODN duplexes were prepared by annealing of modified and complementary strands at the 1∶1 molar ratio in the reaction buffer (50 mM Tris-HCl (pH 7.5), 50 mM KCl, 1 mM EDTA, 1 mM dithiothreitol, 9% glycerol).
Stopped-flow Fluorescence Measurements
Stopped-flow measurements with fluorescence detection were carried out essentially as described previously [14], [20], [21]. A SX.18MV stopped-flow spectrometer (Applied Photophysics Ltd, UK) equipped with a 150 W Xe arc lamp and 2 mm path length optical cell was used. The dead time of the instrument is 1.4 ms. The fluorescence of aPu was excited at λex = 310 nm and monitored at λem>370 nm as transmitted by a Corion filter LG370. The concentration of substrates containing aPu in the reaction cell was 1 µM, and the concentration of hOGG1 protein was varied from 0.5 to 2.0 µM. Typically, each trace shown in the Figures is the average of four or more fluorescence traces recorded in individual experiments. All experiments were carried out in a buffer containing 50 mM Tris-HCl, pH 7.5, 50 mM KCl, 1 mM EDTA, 1 mM DDT and 9% glycerol (v/v) at different temperatures over the range of 10–30°C.
Product Analysis
To analyze the products formed by hOGG1, the oxoG-substrate was 5′-32P-labeled using T4 polynucleotide kinase and [γ-32P] ATP. The reaction mixtures contained reaction buffer, 1.0 µM 32P-labeled oxoG-substrate and 2.0 µM enzyme. The reaction was initiated by adding the enzyme and allowed to proceed at different temperatures (10–30°C). Aliquots (2 µl) were withdrawn as required, mixed with 3 µl of gel-loading dye containing 7 M urea, and analyzed by 20% denaturing PAGE. The gels were exposed to Agfa CP-BU X-ray film (Agfa-Geavert), and the autoradiograms were scanned and quantified using Gel-Pro Analyzer, version 4.0.
Data Processing
The processing of the fluorescence intensity traces was performed by numerical fitting of a solution of a system of ordinary differential equations corresponding to the reaction mechanism using DynaFit software (BioKin, Pullman, WA) [26] as described previously in [27]–[30].
The approach is based on the fluorescence intensity variation in the course of the reaction due to the sequential formation and further transformation of DNA-enzyme complex and its conformers. The software performs numerical integration of a system of ODE with subsequent non-linear least-squares regression analysis. The response factors of the intermediates (that are essentially the products of the molar extinction coefficients and the fluorescence quantum yields) are treated as fitting parameters in the data processing. In the evaluated mechanisms, except for the first bimolecular step, all other reactions are first order reactions. In the data processing, the kinetic information is obtained from the ‘temporal behavior’ of the fluorescence intensity, not from the ‘amplitudes’ of the specific signal contributions. Only the rate constants obtained in the non-linear fits were used to derive the thermochemical parameters. The ‘response factors’ for different conformers resulting from the fits were not used in the determination of the equilibrium constants, but rather provided additional information on the fluorescence intensity variations in different complexes and conformers.
Processing of individual kinetic curves does not unambiguously provide the kinetic parameters; therefore, global fits of sets of kinetic curves obtained at different concentrations of the reactants at each temperature was used. In the fits all relevant rate constants for the forward and reverse reactions, as well as the specific molar responses for all intermediate complexes were optimized.
Using the measured rate constants the equilibrium constants Ki (ki/k–i, i is the step number) were determined for G- and oxoG-substrate. The standard thermodynamic functions of the i-th equilibrium step were determined using the van’t Hoff equation, which represents the relation between the true thermodynamic equilibrium constant (Ki) and the Standard Gibbs Energy, the Standard Enthalpy, and the Standard Entropy of the i-th reaction step [22]:(1)
The typical dependence ln(Ki) versus 1/T was linear, as expected for the relatively narrow temperature range of the study. The Gibbs free energies ΔG°i at 25°C were calculated from ln(Ki) = −ΔG°i/RT. The validity of this approach was discussed elsewhere (e.g., Ref. [31]).
Analysis of the temperature dependence of the rate constant of a chemical reaction ki permits determination of the standard enthalpy of activation (ΔH°,‡) and the standard entropy of activation (ΔS°,‡) based on the Transition State Theory (Eyring equation [22]). For unimolecular reactions, such as the catalytic step 5 in the reaction mechanism,(2)where kB and h are Boltzmann and Planck’s constants, respectively, R is the gas constant, T is absolute temperature in Kelvins, ki is the rate constant of the chemical step i.
Molecular Modeling
The X-ray data obtained for complexes of hOGG1 with G and oxoG containing DNA ([17], PDB ID 1YQK and 1YQR, respectively), are taken as models of atomic structures of protein-DNA complexes (E•OG)1 and (E•OG)3 in Figure 3.
Inspection of these structures revealed, that structure 1YQK represents an example of the initial protein-DNA complex (E•OG)1, with bended DNA interacting with the protein and the G base placed in the exo-site. Structure 1YQR represents the specific protein-DNA complex (E•OG)3 with oxoG base everted from the double helix into the enzyme’s active site, with filled void in the double helix by Arg-154, Arg-204, Tyr-203 and Asn-149 and the enzyme forming a tight complex with the DNA-substrate in the catalytically active conformation. The contact areas ΔMS of the protein-DNA interface were calculated as the difference of the molecular surface (MS) confining the solvent-excluded volume [32] of isolated protein and the DNA and the protein-DNA complexes, ΔMS = MS (protein-DNA complex) – MS (protein) – MS (DNA), where the buried MS has negative sign. The MS surface was calculated by the SIMS method [33]. The free energy of the cavity formation in water in the process of desolvation has entropic nature and it is approximated by the linear equation ΔGcav = γΔMS, where the parameter γ is in the range of 70–117 cal/mol/Å2 [32]. For the evaluation of the entropy loss, γ = 80 cal/mol/Å2 was accepted (which approximates the entropy of desolvation of non-polar species in water [34]).
Results
Interactions of hOGG1 with oxoG-substrate
The thermodynamics of the oxoG recognition by hOGG1 provided very important information for identification of the key steps in the DNA lesion search process. To quantitatively characterize the thermodynamics of the damage recognition process the changes in the fluorescence of an aPu label on 5′-side from oxoG base were used to reveal the conformational changes of DNA in the close proximity to the damaged base. The fluorescence traces measured for various concentrations of enzyme at each temperature (Fig. 4A and 4B) were fitted by the kinetic mechanism shown in Figure 3, as in the previous studies [14], [20], [21]. The rate constants obtained from the fits are summarized in Table 2. The thermodynamic parameters (Table 3) were then calculated using the van’t Hoff analysis (Fig. 5) as described in Materials and Methods.
Changes in aPu fluorescence intensity during interaction of hOGG1 with oxoG-substrate at different concentrations of enzyme at 25°C (A) and at different temperatures (B). Solid lines represent the fitted curves. The kinetics of the accumulation of the abasic (C) and nicked (D) products, formed in the N-glycosylase and AP-lyase reactions, respectively, as detected in the PAGE experiments. The concentrations of hOGG1 and DNA for (B), (C) and (D) panels were 2.0 µM and 1.0 µM, respectively.
The effect of temperature on the rate of the chemical reactions was determined in the direct measurements of the rate of the products formation (abasic product in N-glycosylase reaction and DNA chain nick in AP-lyase reaction) using PAGE analysis of 32P-labeled oxoG-substrate. The results are shown in Figures 4C and 4D. The characteristic times of the abasic intermediate and nicked DNA product accumulation coincide with the characteristic times of aPu fluorescence increase in the 10–1000 s interval, indicating that the fluorescence changes in this time interval characterize the chemical steps of the enzymatic reaction.
Interactions of hOGG1 with Non-specific DNA
The aPu fluorescence increases substantially during the process of hOGG1 binding of the non-specific G-ligand; the process is essentially completed in about 0.2 s (Fig. 6). The observed increase in the fluorescence intensity of aPu indicates local melting of the DNA chain interior of the non-specific enzyme-DNA complex. As it is known from the structural data obtained for covalently linked complex hOGG1 with a non-specific DNA, the enzyme is capable to flip out the undamaged G base not only in the exo-site [17] but also in the active site [19]. According to the kinetic data obtained for a set of hOGG1 mutant forms [21], the non-specific binding includes pulling of the opposite C base with Arg-154 and Arg-204, pushing by the Tyr-203 residue into the duplex between cytosine and the next base, the H-bond formation between the His-270 residue and the phosphate group of the damaged nucleotide, and stacking of the G base with the Phe-319 residue. All these factors lead to disruption of the stacking between aPu and G bases due to the flipping out of the latter and the DNA bending, which are expected to enhance aPu fluorescence. The aPu fluorescence traces obtained for G-ligand do not exhibit decreases in the fluorescence intensity as would be indicative of the void-filling processes. This means that the full insertion of the G base in the active site of hOGG1 as well as interaction of the opposite C with Arg-154, Arg-204, Asn-149 and Tyr-203 and its insertion into the DNA void, as suggested by the structural studies [19], is kinetically unfeasible.
Changes in aPu fluorescence intensity during interaction of hOGG1 with G-ligand at different concentrations of enzyme at 25°C (A) and at different temperatures (B). Solid lines represent the fitted curves. The concentrations of hOGG1 and G-ligand were 2.0 µM and 1.0 µM, respectively. (C) Van’t Hoff analysis of the temperature dependence of Ki for G-ligand.
The measured fluorescence traces at various concentrations of enzyme at each temperature (Fig. 6) were fitted by a kinetic mechanism containing two binding steps (Fig. 7). The forward and the reverse rate constants as well as the dissociation constant obtained at different temperatures are listed in Table 4. The thermodynamic parameters calculated as described in Materials and Methods are listed in Table 5.
Discussion
Comparative analysis of the thermodynamic data presented in Tables 3 and 5 together with the structural information for free hOGG1 and various hOGG1/DNA complexes [9], [10], [15]–[19] provides important information on the energy consumption in the stages of the DNA binding, specific sites recognition and catalysis.
The process of G-ligand binding proceeds in two steps, which could belong to the formation of non-specific electrostatic and hydrophobic contacts with the sugar-phosphate backbone, DNA bending, Arg-154, Arg-204, Asn-149 and Tyr-203 amino acids insertion into the DNA and the G base flipping out into the exo- or active site of the hOGG1, as follows from the X-ray studies [17], [19].
The formation of the 1st complex is favorable in terms of the Gibbs energy change, due to both negative enthalpy and positive entropy changes (Table 5). The second detectable binding step for G-ligand is characterized by the increase of aPu fluorescence and completed within ca. 0.2 s. As follows from Figure 6, no decrease in the aPu fluorescence was detected, indicating that the full insertion of mentioned above amino acids into the DNA void does not occur. The formation of a tighter 2nd enzyme-ligand complex is characterized by the decrease of entropy compensated by a favorable change of enthalpy, indicating formation of new contacts, which stabilize this complex and make it more rigid (Table 5). Overall, binding of hOGG1 with G-ligand is thermodynamically favorable process with ΔG° = −6.9 kcal/mol.
Therefore, hOGG1 binding to G-ligand as monitored by direct fluorescence comprises melting of the duplex and eversion of the G base, with the G base left mainly outside the active center. This follows from the absence of kinetic evidence of full insertion of the enzyme amino acids into the DNA duplex, which is expected in the case when the G base is placed in the active site similar to that as for the oxoG base. Therefore, it is reasonable to conclude that in the case of true equilibrium and kinetically restricted state for the hOGG1/G-DNA complex the location of G base in the spatial structure could be different.
The fluorescence traces obtained for oxoG-substrate (Fig. 4A and 4B) show that at 0.2 s the aPu fluorescence starts to decrease, which characterizes the insertion of the amino acids. It means that the insertion of amino acids occurs immediately after the eversion of oxoG base from DNA-helix. Therefore, the second step is very important to discriminate G bases vs. oxoG bases. According to [21], specific interactions between Arg-154, Arg-204, Asn-149 and DNA chain occur during this step. Additionally, the oxoG is placed in the active site where it forms stacking structure with Phe-319.
For oxoG-substrate, the formation of the second (E•OG)2 and the third (E•OG)3 complexes in Figure 3 is characterized by the positive enthalpy and entropy changes (Table 3). A several processes during the protein/DNA complexation can lead to a significant unfavorable contribution to the enthalpy. The desolvation of polar groups at protein/DNA interface is expected to be energetically unfavorable [35]. Moreover, the protein bindings in the minor groove exhibit positive enthalpy and positive entropy changes that compensates the process endothermicity. Therefore, these processes are completely entropy driven [36]. The unfavorable enthalpy changes for the minor groove binding proteins result from the displacement of water molecules that are in the minor groove. Because these water molecules are highly ordered, their release into the bulk solution leads to a significant entropy growth. In the case of hOGG1, the aryl ring of Tyr-203 is wedged into the DNA duplex from the side of the minor groove [9]. Therefore these results are in accord with the conclusion that the enzyme interacts with the minor groove of oxoG-substrate [9].
Another potential source of unfavorable enthalpy ΔH° is the stress of the DNA structure induced by the ligand attachment. The comparison of the enthalpy and entropy changes for a wide range of DNA-binding proteins [35] shows that complexes in which the DNA is bent or severely distorted exhibit an unfavorable enthalpy changes compensated by favorable entropy changes. According to the structural data, in the catalytically active complex hOGG1 bends the DNA (ca. 70°) at the lesion site. Therefore, DNA bending can be also associated with these steps of oxoG-substrate binding by hOGG1.
The X-ray structures of the complexes of hOGG1 with G (1YQK) and oxoG (1YQR) containing DNA [17] might be considered as the atomic models for the hOGG1-DNA complexes corresponding to the initial complex (E•OG)1 and the catalytically competent complex (E•OG)3, respectively. Therefore, it was possible to estimate the change of the area of the protein-DNA interface during the formation of (E•OG)3 from (E•OG)1. The calculations show that the mutual adjustment of the enzyme and DNA in the complex (E•OG)3 is accompanied by the increase of the contact surface interface by 672 Å2 from 764 to 1436 Å2 in complexes (E•OG)1 and (E•OG)3, respectively. Burring of 672 Å2 of the molecular surface of the tight protein-DNA complex results in the entropy increase of about 180 cal mol−1 K−1 [34]. This entropy increase is partially compensated by the decrease of entropy due to the formation of the new tight protein-DNA contacts in the course of oxoG-substrate binding steps (Table 3).
The transition state for the catalytic step of the base excision (N-glycosylase reaction, step 4) as well as the 3′-phosphodiester bond cleavage (β-elimination reaction, step 5) has unfavorable enthalpy ΔH°,‡ and favorable entropy ΔS°,‡ changes (Table 3). The results of the thermodynamic analysis of the catalytic steps obtained using aPu fluorescence data were qualitatively confirmed by the direct PAGE analysis of the products accumulation at different temperatures.
As it was recently demonstrated for another oxoG excision DNA glycosylase Fpg from E. coli, belonging to the Fpg/Nei family of DNA glycosylases, the determination of the thermodynamic parameters of the substrate binding and lesion removal steps provided important information about oxoG recognition process [27], [29], [30]. The pre-steady-state kinetic data were obtained using pyrrolocytosine fluorescence in DNA [27]. These data allowed determination of the Gibbs energy, enthalpy and entropy for the five rapidly proceeding steps of sequential damage recognition mechanism and to develop a detailed model of the Fpg enzymatic process. The first step had similar thermodynamic parameters both for undamaged G-ligand and oxoG-specific substrate, as it was obtained for hOGG1 catalysis in the present study. This step is highly exothermic. In this step the favorable enthalpy change is accompanied by an increase in entropy most probably due to the DNA melting in the point of contact and the conformational changes of the DNA binding site in the Fpg molecule with the emphasis on the contribution of the Phe-110 wedge movement. The next 2nd, 3rd and 4th steps are characterized by the negative or slightly positive ΔG° changes and reflect the formation of a kink in the DNA chain; eversion of oxoG base from the double helix into the enzyme’s active site; and filling the resulting void in the double helix by Arg-108 and Met-73. The 5th reversible recognition step was characterized by a highly unfavorable enthalpy that was compensated by favorable entropic contribution. The large positive entropy of the protein binding originated from the dehydration of the DNA grooves and isomerization of the enzyme to form a tight complex with the DNA-substrate and to produce catalytically active conformation.
Conclusions
Summarizing, in the present study a thermodynamic analysis of the reaction of human 8-oxoguanine DNA glycosylase with DNA duplexes containing undamaged and damaged nucleotides was performed. It was shown that the first binding step only weakly depends on the nature of the damaged nucleotide, supporting the hypothesis that the main contribution to the thermodynamics of this step is due to the non-specific interactions between the enzyme and the DNA backbone. The results of this study clearly demonstrate that the G base does not evert into the active site, however, do not exclude a possibility for the evertion of the G base in the exo-site suggested earlier by Verdine et al. [17]. The second step is crucial for the discrimination between oxoG and G bases. In the case of oxoG-substrate the second and the third steps have unfavorable enthalpy, which is compensated by favorable desolvation entropy. During these steps displacement of water molecules from the minor groove as well as DNA bending occurs. The last binding step also includes the final adjustment of the enzyme/substrate complex to the catalytically competent state. This step might be attributed to the full insertion of Tyr-203 in the DNA helix accompanied with the formation contacts between oxoG and Gly-42.
The thermodynamic analyses of the catalytic processes performed by functionally homologous but structurally different DNA-glycosylases, human OGG1 (this work) and E. coli Fpg [27], belonging to two distinct structural families, HhH-GPD and Fpg/Nei, respectively, demonstrate that they employ the common energetic features in the main steps of 8-oxoguanine lesion recognition. The first step of the DNA binding process is the formation of a non-specific complex, resulting in the double helix melting and insertion of wedging aminoacids (Tyr-203 for hOGG1 and Phe-110 for Fpg), characterized by favorable changes of both enthalpy and entropy. The following steps of lesion recognition are favorable in enthalpy but may be unfavorable in entropy due to the formation of structurally rigid complexes. The last step of the lesion recognition where enzyme-substrate complexes are rearranged to catalytically competent conformations is characterized by large endothermicity compensated by a significant increase of entropy originated from the dehydration of DNA grooves.
Acknowledgments
We thank Dr. Alexander A. Ishchenko and Dr. Murat K. Saparbaev, Groupe «Réparation de l’ADN», Université Paris-Sud XI, UMR8200 CNRS, Institute Gustave Roussy, Villejuif Cedex F-94805, France, for providing the pET28c-hOgg1 plasmid and assistance with protein isolation.
Author Contributions
Conceived and designed the experiments: NAK. Performed the experiments: NAK AAK YNV. Analyzed the data: NAK YNV LNK OSF. Wrote the paper: NAK LNK OSF.
References
- 1. Evans MD, Dizdaroglu M, Cooke MS (2004) Oxidative DNA damage and disease: induction, repair and significance. Mutat Res 567: 1–61.
- 2. Bjelland S, Seeberg E (2003) Mutagenicity, toxicity and repair of DNA base damage induced by oxidation. Mutat Res 531: 37–80.
- 3. Wallace SS (2002) Biological consequences of free radical-damaged DNA bases. Free Radic Biol Med 33: 1–14.
- 4. Rumyantseva GV, Weiner LM, Frolova EI, Fedorova OS (1989) Hydroxyl radical generation and DNA strand scission mediated by natural anticancer and synthetic quinones. FEBS Lett 242: 397–400.
- 5. Kuznetsova AA, Knorre DG, Fedorova OS (2009) Oxidation of DNA and its components with reactive oxygen species. Russian Chemical Reviews 78: 659–678.
- 6. Cheng KC, Cahill DS, Kasai H, Nishimura S, Loeb LA (1992) 8-Hydroxyguanine, an abundant form of oxidative DNA damage, causes G–T and A–C substitutions. J Biol Chem 267: 166–172.
- 7. Thayer MM, Ahern H, Xing D, Cunningham RP, Tainer JA (1995) Novel DNA binding motifs in the DNA repair enzyme endonuclease III crystal structure. EMBO J 14: 4108–4120.
- 8. Nash HM, Bruner SD, Scharer OD, Kawate T, Addona TA, et al. (1996) Cloning of a yeast 8-oxoguanine DNA glycosylase reveals the existence of a base-excision DNA-repair protein superfamily. Curr Biol 6: 968–980.
- 9. Bruner SD, Norman DP, Verdine GL (2000) Structural basis for recognition and repair of the endogenous mutagen 8-oxoguanine in DNA. Nature 403: 859–866.
- 10. Bjoras M, Seeberg E, Luna L, Pearl LH, Barrett TE (2002) Reciprocal “flipping” underlies substrate recognition and catalytic activation by the human 8-oxo-guanine DNA glycosylase. J Mol Biol 317: 171–177.
- 11. Bjoras M, Luna L, Johnsen B, Hoff E, Haug T, et al. (1997) Opposite base-dependent reactions of a human base excision repair enzyme on DNA containing 7,8-dihydro-8-oxoguanine and abasic sites. EMBO J 16: 6314–6322.
- 12. Zharkov DO, Rosenquist TA, Gerchman SE, Grollman AP (2000) Substrate specificity and reaction mechanism of murine 8-oxoguanine-DNA glycosylase. J Biol Chem 275: 28607–28617.
- 13. Radicella JP, Dherin C, Desmaze C, Fox MS, Boiteux S (1997) Cloning and characterization of hOGG1, a human homolog of the OGG1 gene of Saccharomyces cerevisiae. Proc Natl Acad Sci USA 94: 8010–8015.
- 14. Kuznetsov NA, Koval VV, Zharkov DO, Nevinsky GA, Douglas KT, et al. (2005) Kinetics of substrate recognition and cleavage by human 8-oxoguanine-DNA glycosylase. Nucleic Acids Res 33: 3919–3931.
- 15. Norman DP, Chung SJ, Verdine GL (2003) Structural and biochemical exploration of a critical amino acid in human 8-oxoguanine glycosylase. Biochemistry 42: 1564–1572.
- 16. Fromme JC, Bruner SD, Yang W, Karplus M, Verdine GL (2003) Product-assisted catalysis in base-excision DNA repair. Nat Struct Biol 10: 204–211.
- 17. Banerjee A, Yang W, Karplus M, Verdine GL (2005) Structure of a repair enzyme interrogating undamaged DNA elucidates recognition of damaged DNA. Nature 434: 612–618.
- 18. Radom CT, Banerjee A, Verdine GL (2007) Structural characterization of human 8-oxoguanine DNA glycosylase variants bearing active site mutations. J Biol Chem 282: 9182–9194.
- 19. Crenshaw CM, Nam K, Oo K, Kutchukian PS, Bowman BR, et al. (2012) Enforced presentation of an extrahelical guanine to the lesion recognition pocket of human 8-oxoguanine glycosylase, hOGG1. J Biol Chem 287: 24916–24928.
- 20. Kuznetsov NA, Koval VV, Nevinsky GA, Douglas KT, Zharkov DO, et al. (2007) Kinetic conformational analysis of human 8-oxoguanine-DNA glycosylase. J Biol Chem 282: 1029–1038.
- 21. Kuznetsova AA, Kuznetsov NA, Ishchenko AA, Saparbaev MK, Fedorova OS (2014) Step-by-Step Mechanism of DNA Damage Recognition by Human 8-Oxoguanine DNA Glycosylase. Biochim Biophys Acta 1840: 387–395.
- 22.
Atkins P, Paula J (2006) Atkins’ Physical Chemistry: Oxford university press.
- 23. Dunlap CA, Tsai MD (2002) Use of 2-aminopurine and tryptophan fluorescence as probes in kinetic analyses of DNA polymerase β. Biochemistry 41: 11226–11235.
- 24. Hardman SJ, Thompson KC (2006) Influence of base stacking and hydrogen bonding on the fluorescence of 2-aminopurine and pyrrolocytosine in nucleic acids. Biochemistry 45: 9145–9155.
- 25. Bradford MM (1976) A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72: 248–254.
- 26. Kuzmic P (1996) Program DYNAFIT for the analysis of enzyme kinetic data: application to HIV proteinase. Anal Biochem 237: 260–273.
- 27. Kuznetsov NA, Vorobjev YN, Krasnoperov LN, Fedorova OS (2012) Thermodynamics of the multi-stage DNA lesion recognition and repair by formamidopyrimidine-DNA glycosylase using pyrrolocytosine fluorescence–stopped-flow pre-steady-state kinetics. Nucleic Acids Res 40: 7384–7392.
- 28. Kuznetsov NA, Koval VV, Zharkov DO, Fedorova OS (2012) Conformational dynamics of the interaction of Escherichia coli endonuclease VIII with DNA substrates. DNA Repair (Amst) 11: 884–891.
- 29. Kuznetsov NA, Zharkov DO, Koval VV, Buckle M, Fedorova OS (2009) Reversible Chemical Step and Rate-Limiting Enzyme Regeneration in the Reaction Catalyzed by Formamidopyrimidine-DNA Glycosylase. Biochemistry 48: 11335–11343.
- 30. Koval VV, Kuznetsov NA, Ishchenko AA, Saparbaev MK, Fedorova OS (2010) Real-time studies of conformational dynamics of the repair enzyme E. coli formamidopyrimidine-DNA glycosylase and its DNA complexes during catalytic cycle. Mutat Res 685: 3–10.
- 31. Ragone R, Colonna G, Ambrosone L (1995) Reliability of the van’t Hoff Plots. J Phys Chem 99: 13050.
- 32. Vorobjev YN (2011) Advances in implicit models of water solvent to compute conformational free energy and molecular dynamics of proteins at constant pH. Adv Protein Chem Struct Biol 85: 281–322.
- 33. Vorobjev YN, Hermans J (1997) SIMS: computation of a smooth invariant molecular surface. Biophys J 73: 722–732.
- 34. Rose GD, Wolfenden R (1993) Hydrogen bonding, hydrophobicity, packing, and protein folding. Annu Rev Biophys Biomol Struct 22: 381–415.
- 35. Jen-Jacobson L, Engler LE, Jacobson LA (2000) Structural and thermodynamic strategies for site-specific DNA binding proteins. Structure 8: 1015–1023.
- 36. Privalov PL, Dragan AI, Crane-Robinson C (2011) Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components. Nucleic Acids Res 39: 2483–2491.