Next Article in Journal
Understanding the Role of the Glial Scar through the Depletion of Glial Cells after Spinal Cord Injury
Next Article in Special Issue
A Knockout of Poly(ADP-Ribose) Polymerase 1 in a Human Cell Line: An Influence on Base Excision Repair Reactions in Cellular Extracts
Previous Article in Journal
Developing Oncolytic Viruses for the Treatment of Cervical Cancer
Previous Article in Special Issue
The Activity of Natural Polymorphic Variants of Human DNA Polymerase β Having an Amino Acid Substitution in the Transferase Domain
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 12
Issue 14
10.3390/cells12141839
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2

by
Anastasiia T. Davletgildeeva
1,*,
Timofey E. Tyugashev
1,
Mingxing Zhao
2,
Nikita A. Kuznetsov
1,2,
Alexander A. Ishchenko
3,
Murat Saparbaev
3 and
Aleksandra A. Kuznetsova
1,*
1
Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
2
Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
3
Groupe Mechanisms of DNA Repair and Carcinogenesis, CNRS UMR9019, Gustave Roussy Cancer Campus, Université Paris-Saclay, CEDEX, F-94805 Villejuif, France
*
Authors to whom correspondence should be addressed.
Cells 2023, 12(14), 1839; https://doi.org/10.3390/cells12141839
Submission received: 28 March 2023 / Revised: 29 June 2023 / Accepted: 11 July 2023 / Published: 13 July 2023
(This article belongs to the Special Issue Epigenetics and DNA Repair: Regulatory Mechanisms and Therapeutic Opportunities)
Browse Figures
Versions Notes

Abstract

:
Human Fe(II)/α-ketoglutarate-dependent dioxygenase ABH2 plays a crucial role in the direct reversal repair of nonbulky alkyl lesions in DNA nucleobases, e.g., N1-methyladenine (m1A), N3-methylcytosine (m3C), and some etheno derivatives. Moreover, ABH2 is capable of a less efficient oxidation of an epigenetic DNA mark called 5-methylcytosine (m5C), which typically is a specific target of DNA dioxygenases from the TET family. In this study, to elucidate the mechanism of the substrate specificity of ABH2, we investigated the role of several active-site amino acid residues. Functional mapping of the lesion-binding pocket was performed through the analysis of the functions of Tyr122, Ile168, and Asp173 in the damaged base recognition mechanism. Interactions of wild-type ABH2, or its mutants Y122A, I168A, or D173A, with damaged DNA containing the methylated base m1A or m3C or the epigenetic marker m5C were analyzed by molecular dynamics simulations and kinetic assays. Comparative analysis of the enzymes revealed an effect of the substitutions on DNA binding and on catalytic activity. Obtained data clearly demonstrate the effect of the tested amino acid residues on the catalytic activity of the enzymes rather than the DNA-binding ability. Taken together, these data shed light on the molecular and kinetic consequences of the substitution of active-site residues for the mechanism of the substrate recognition.

1. Introduction

Exposure of the cell to alkylating agents originating from external and internal sources constantly damages DNA, thereby leading to the formation of various alkylated lesions that have a cytotoxic and/or mutagenic effect [1,2,3]. Some alkylating agents are also widely used in antitumor therapies and exert their action by generating cytotoxic lesions in tumor cells [4,5]. At the same time, targeted enzymatic DNA methylation plays an important role in the epigenetic regulation of gene expression [6,7]. There is no doubt in the extreme importance of thorough control over the alkylating damage to DNA bases. The participation of DNA methylation in the processes of DNA damage and in the regulation of DNA functions as an information carrier is well characterized; however, demethylation, which is also a very important process in the modulation of the DNA methylation status, remains to be elucidated.
In the past three decades, three main types of demethylation enzymes have been identified in the cells that protect the genome from any alkylated modifications. These types are methyltransferases, repair DNA glycosylases, and repair and epigenetic DNA dioxygenases [1,2,8,9]. Among these types of enzymes, repair DNA dioxygenases of the AlkB family are of particular interest [10,11] because these enzymes carry out a single enzymatic reaction to direct the catalytical repair of various alkyl DNA lesions, such as N1-methyladenine (m1A) and N3-methylcytosine (m3C), through an oxidative mechanism involving the Fe2+ ion, α-ketoglutarate (αKG), and oxygen [12,13,14,15,16,17]. After the catalytic oxidation of the methyl group, the decomposition of the hydroxylated intermediate proceeds with the formation of the undamaged base and formaldehyde. Nine orthologs of E. coli AlkB dioxygenase have been found in human cells [9]. Nonetheless, only two human enzymes homologous to AlkB have been identified, ABH2 and ABH3, which are capable of the oxidative demethylation of m1A and m3C [18,19,20]. It is noteworthy that enzymes AlkB and ABH3 have higher activity toward single-stranded DNA (ssDNA) than dsDNA, whereas ABH2 seems to prefer dsDNA substrates [19,21]. It is known that AlkB and ABH3, but not the ABH2 enzyme, can reverse m1A and m3C lesions in RNA substrates [19,22]. These data apparently indicate that AlkB and ABH3 predominantly repair DNA during replication, when segments of ssDNA are present, whereas the ABH2 enzyme is considered a DNA-repair enzyme acting primarily during the S phase through the interaction with the proliferating cell nuclear antigen (PCNA) [23,24].
To date, it has been thought that dioxygenase ABH2 is only a repair enzyme that protects genomic DNA from nonbulky alkyl lesions, such as m1A, m3C, 3-methylthymidine (m3T), 1,N6-ethenoadenine (εA), 3,N4-ethenocytosine (εC), and N2,3-ethenoguanine [25,26,27]. On the other hand, a recent study additionally revealed an ability of dioxygenase ABH2 to remove the epigenetic mark 5-methylcytosine (m5C) from ssDNA and dsDNA in vitro [28]. A major product of the m5C oxidation is 5-hydroxymethylcytosine (hm5C), but products of its further oxidation to 5-formylcytosine (f5C) and 5-carboxylcytosine (ca5C) form only in negligible amounts. Given that the demethylation of m5C is an important epigenetic process, the above findings suggest that a living organism can employ several redundant or complementary pathways to control an epigenetic modification. The possible biological significance of the discovered activity may be the implementation of a reserve pathway of epigenetic regulation in addition to the main mechanism, which is realized with the participation of TET family dioxygenases [29,30].
Nevertheless, the question of how the ABH2 enzyme discriminates alkylated DNA substrates from other substrates remains a subject of increased interest among researchers. As follows from X-ray data, an effective interaction of AlkB family dioxygenases with DNA and RNA strands is implemented by characteristic loop regions, so-called β-hairpins, whose size and amino acid composition differ substantially among homologs [31,32,33,34,35,36,37]. In particular, ABH2 contains two β-sheet hairpin motifs: the first one intercalates into the DNA duplex by means of residues Val101 and Phe102 and occupies the space left by the damaged base flipped out into the active site, whereas the other motif interacts with the intact DNA strand (Figure 1). To bind the opposite DNA strand, ABH2 uses a short positively charged loop containing Arg241, Lys242, and Lys243, as well as an additional flexible loop carrying the DNA-binding residues Arg198, Gly204, and Lys205. The Fe2+ ion is coordinated by the amino acid triad His171, Asp173, and His236.
Analysis of the available crystal structures of wild-type (WT) ABH2 with DNA containing m1A, m3C, or εA inside the active site (Figure 2) allowed us to compile a list of the amino acid residues (Table 1) that form the damaged-base-binding pocket, and thereby could be involved in the damage recognition and control of substrate specificity. It was found that the wall of the damaged-base pocket is formed by Val99, Arg110, Ser125, and Ile168. All types of flipped-out damaged bases are planarly placed between two aromatic residues, Phe124 and His171, resulting in stacking. Additionally, Tyr122, Asp173, and Glu175 come into direct contact with bases m1A and m3C, but not with εA. Residues Val101 and Phe102 of the DNA-intercalating loop interact with the 5′ and 3′ bases that flank the damaged base.
Several attempts have been made to conduct an extended analysis of the role of the key amino acid residues of the ABH2 enzyme [20,38]. Lee et al. report that substitution D173A inactivates the enzyme toward single-stranded substrates, whereas H236A leads only to some reduction in the activity toward both m1A and m3C bases. The phenylalanine “finger” Phe102 has been found to be crucial for ssDNA repair; however, the second intercalating residue Val101 shows only reduced demethylating activity toward ssDNA and has no effect on dsDNA. The substitution R110A completely eliminates the enzymatic activity, suggesting that Arg110 is involved in the base-flipping process. It has also been shown that Glu175 and Phe124 contribute to nucleotide-base-specific selection and stabilization in the active site for repair.
To elucidate the mechanism of the substrate specificity of the ABH2 enzyme, analyses of conformational changes occurring in the enzyme–substrate complex during the catalytic cycle have been performed by quantum mechanical/molecular mechanical (QM/MM) modeling [39,40] and by the stopped-flow fluorescence spectroscopy technique [41]. The results indicate that the conformational flexibility of protein domains and of the DNA substrate is important for their effective interaction.
Here, we report the mutational analysis of three of the eleven listed amino acid residues in Table 1, which form contacts with damaged bases inside the active site of ABH2. Tyr122, Ile168, and Asp173 amino acid residues were specifically chosen as the ones able to form direct contacts with the m1A or m3C lesions from three different spatial points. To elucidate the role of these residues in the mechanism of recognition by the ABH2 enzyme of the damaged methylated base m1A or m3C or the epigenetic mark m5C, we created single-point mutant forms containing the substitution of selected amino acids. In the current study, we decided to adhere to the ‘classic’ switching-off strategy to elucidate the impact of the individual amino acid residues into the catalytic process, thus making the substitution to alanine specifically. Combining computational molecular dynamics (MD) simulations with the experimental detection of DNA binding and enzyme activity allowed us to elucidate the role of the selected amino acid residues in the ABH2 action.

2. Results and Discussion

2.1. Choosing Potentially Important Amino Acid Residues Forming Specific Contacts with Damaged Bases

The mechanisms of target nucleotide recognition by site-specific enzymes that were elucidated in the last few years for a wide range of enzymes, such as DNA glycosylases [42,43,44,45,46,47] and AP-endonucleases [48,49,50,51,52], suggest that these enzymes utilize a combined approach, including both the indirect recognition of a lesion through the induction of the eversion of the target nucleotide in the active site by DNA bending and the unwinding and subsequent direct verification of the lesion in the active site via the formation of damage-specific contacts. It should be noted that crystal structures enable the analysis of only those contacts between amino acid residues and a damaged base that are already formed inside the binding pocket, whereas the trajectory of the base eversion from the DNA duplex in the active site may also reveal a very important contribution (if any) to the efficiency of damage recognition and substrate specificity.
In our analysis of the structures of the ABH2 enzyme in complex with DNA containing m3C (PDB ID: 3RZJ) or m1A (PDB ID: 3BTY) as damage, a number of amino acid residues were selected that were of potential interest for studying their role in the step-by-step damaged-base-recognition process. This analysis of the structures of the ABH2–DNA complexes (Figure 3) revealed that amino acid residues Tyr122, Ile168, and Asp173 engage in contacts with both m3C and m1A (damaged bases) placed inside the active site; these observations supported the possible importance of these residues for the damage-checking mechanism of substrate specificity. Therefore, these amino acid residues were chosen to verify their participation in DNA binding, catalytic complex formation, and catalysis.

2.2. MD Simulations of Mutant Forms of the Enzyme

The functions of the amino acid residues Tyr122, Ile168, and Asp173 in damage recognition were estimated by MD simulations. ABH2 cross-linked with DNA in the presence of the cofactor αKG and the Mn2+ ion (PDB ID: 3RZJ for m3C-containing DNA and PDB ID: 3BUC for m1A-containing DNA) were used as initial simulation structures [32,39]. The conformational equilibration of these initial structures helped to obtain models of WT ABH2 with DNA containing m3C or m1A (Figure 4a,b). It turned out that MD relaxation leads to the slight movement of damaged nucleotides within the active site in comparison with crystal structures. Such movement proceeds owing to the formation of set contacts between the damaged-base and active-site amino acid residues, thereby resulting in the optimized positioning of the damaged nucleotides. Then, complexes of ABH2 mutant forms Y122A (Figure 4c,d), I168A (Figure 4e,f), or D173A (Figure 4g,h) with DNA containing m1A or m3C were modeled to analyze conformational rearrangements induced by the selected amino acid substitution.
The obtained MD model of Y122A (Figure 4c,d) revealed that the substitution of Tyr122 causes a loss of the hydrogen bond with Arg254 and the disappearance of the salt bridge Arg254–Asp173, which is important for the coordination of the Fe2+ ion in the active site by Asp173. Such destabilization of the Fe2+ ion’s position disturbs the metal ion coordination sphere and most likely should affect the catalytic efficacy.
The substitution of I168A (Figure 4e,f) leads to the widening of the active-site pocket, thereby allowing an αKG molecule to move, rotate, and interact with the Fe2+ ion only via the carboxyl group. It should be pointed out that such orientations also occurred in our simulation trajectory for the WT enzyme, albeit for a much shorter period; therefore, the catalytic activity of ABH2 I168A most likely is low, owing to a catalytically “wrong” position of αKG.
In the case of D173A (Figure 4g,h), the stacking of the damaged bases with the Phe124 residue is disturbed owing to a different position of the target nucleotide. Moreover, the undereversion of both tested damaged bases leads to its greater mobility in the active-site pocket. Although the DNA conformation is imperfect for the catalytical state, the Fe2+ ion retains its catalytically competent position due to an additional water molecule in place of Asp173.
To clarify the specific features of the substrate-specificity mechanism that allows for the recognition of the epigenetic mark m5C in the active site of ABH2, as reported recently [28], MD models of WT ABH2 with DNA containing an m5C base were set up (Figure 5). These models enabled us to determine the importance of the contacts of the damaged base with active-site amino acids. It should be noted that the m5C base could be positioned in the active site in two alternative ways. In the first (“noncatalytic”) position, the m5C base is positioned similarly to the m3C base, and thus the methyl group is located far away from the catalytic amino acid triad in a position unsuitable for oxidation. In such an orientation, most of the stabilizing contacts with the damaged base are formed, which retain the base inside the active site. In the other (“catalytic”) position, the N-glycosidic bond of m5C is rotated to place the methyl group close to the catalytic amino acids. Nonetheless, in this position, there is a loss of many contacts that hold the everted base inside the active site. Taken together, these data suggest that an equilibrium between these positions will significantly affect the ability of ABH2 to catalyze the oxidation of the m5C base. The disturbance of the active-site architecture by the selected amino acid substitution could influence this equilibrium, owing to a loss of some specific contacts inside the active site.

2.3. Circular Dichroism (CD) Spectra

WT ABH2 and its mutants featuring the Y122A, I168A, or D173A substitutions were purified to experimentally evaluate the activity of these enzymes toward DNA containing m3C, m1A, or m5C. First of all, CD spectra of these recombinant proteins were recorded to compare the effects of the chosen substitutions on the structure of the enzyme (Figure 6). The results showed similar structural patterns and allowed us to conclude that the selected substitutions do not cause a global alteration of the secondary structure of the protein.

2.4. Equilibrium Binding of WT ABH2 or Its Mutant Forms to Methylated DNA

To assess the influence of the amino acid substitutions Y122A, I168A, and D173A on the ability of the enzyme to bind methylated DNA substrates, electrophoretic mobility shift assays (EMSAs) were performed in the presence of 5′-FAM-labeled oligonucleotides containing m1A or m3C as a lesion (Figure 7a,b). The enzymes were titrated with 17-nt dsDNA containing a methylated nucleotide. Equilibrium dissociation constants (Kd) for complexes of the enzyme with target dsDNA (Figure 7c,d) were calculated by means of Equation (1). Dissociation constants were in a narrow range of 2.1–2.7 µM for the binding to the m1A-containing DNA and 2.5–4.2 µM for the binding to the m3C-containing DNA (Figure 8). It should be mentioned that these values are in good agreement with data obtained previously in an equilibrium fluorescent titration assay [41].
The obtained Kd values indicate that the WT ABH2 enzyme—just as with its mutant forms—binds to dsDNA containing the methylated base m1A with an affinity similar to that for m3C-containing dsDNA, with a slight preference for m1A. It must be noted that, in the range of Kd values for m1A binding (D173A ≤ Y122A < WT ≤ I168A), the greatest value was observed in the case of I168A, and this value was only 1.3-fold greater than the smallest one (in case of D173A), meaning that all of the introduced amino acid substitutions have only a minor influence on the processes involving the binding of an m1A-containing substrate. For m3C, the impact of the substitution manifested itself as a 1.6-fold difference between the greatest value (for I168A) and the lowest one (for Y122A). It is also worth noting that in both cases—m1A- and m3C-containing DNA—the substitution of Ile168 weakens the binding affinity of the ABH2 enzyme.

2.5. Activity of WT ABH2 and Its Mutant Forms toward DNA Containing m1A or m3C

The efficiency of active demethylation of damaged dsDNA substrates by WT ABH2 or its mutant forms Y122A, I168A, or D173A was measured by a restriction-coupled assay followed by gel electrophoresis (Figure 9a,b). The reaction mixture contained equimolar amounts of the enzyme and of its substrate, which was preincubated with Fe2+ ions in the presence of αKG. After 30 min, the reaction was terminated, and the products were digested with a specific restriction endonuclease to cleave the product of the active demethylation. The shortened reaction product was then visualized in a gel, and thus the activity of each enzyme on both substrates, m1A and m3C DNA, was determined.
As shown in Figure 9c,d, WT ABH2 and all mutant forms were able to cleave the dsDNA substrate containing m1A or m3C. The highest level of product accumulation was observed during the active demethylation of the m1A- or m3C-containing dsDNA substrate by WT ABH2 (54.0% and 75.4%, respectively; Figure 9c,d). The incomplete substrate conversion could be associated with some degradation of the enzyme, substrate, or one of the cofactors during the long reaction process. Indeed, during the reaction, the concentration of Fe2+ ions or O2 can substantially decrease in the solution. Additionally, this effect could come out of the inhibition of the enzyme by the products of the reaction, such as succinate. However, this enzyme assay allows us to compare the activity of mutant forms. As such, the substitution of Ile168 has different consequences for the tested lesions: the level of the product accumulation was very similar to that of the WT enzyme for m1A (52.1% and 54.0%, respectively), whereas it was low in the case of m3C (50% and 75.4%, respectively). As determined by MD simulations (Figure 4), the I168A substitution causes the widening of the active-site pocket and increases the movements of the αKG molecule within the active site. It can be hypothesized that, in the case of the smaller m3C base, such increased mobility of αKG leads to the shorter average time of the catalytically optimal orientation of the substrate and cofactor. In the cases of the substitutions Y122A and D173A, the product accumulation for both m1A and m3C was significantly lower in comparison with the WT enzyme. Nonetheless, the decrease of activity toward m3C was again much more drastic, supporting the notion that a disturbance of the Fe2+ ion coordination environment plays a more important role for this damaged base. Accordingly, it can be concluded that the substitutions of Tyr122, Ile168, and Asp173 all weaken the preference for m3C over m1A, seen in WT ABH2.
To clarify the details of the impact of the amino acid substitutions on specific features of the reaction of the ABH2 enzyme with the methylated DNA substrates, the kinetics of product accumulation during the active demethylation of DNA substrates containing m1A or m3C by WT ABH2 or its mutant forms were investigated next. Time courses of the substrate demethylation are presented in Figure 10. A comparison of the time courses of the active demethylation of the m1A lesion between WT ABH2 and its mutant forms suggested that, despite similar levels of overall product accumulation, the repair of this lesion by the ABH2 I168A mutant form was slower than that by WT ABH2 (Figure 10a). The substitution of Tyr122 or Asp173 leads to the substantial deceleration of the enzyme, thus reducing the overall efficiency, and the effects are quite similar between the two substitutions, indicating that the effects on the catalytic activity of the ABH2 enzyme are nearly identical when these two amino acids are compared. Moreover, the much lower activity of Y122A and D173A mutant forms toward the m3C lesion was similar to that for the m1A lesion (Figure 10b). The decrease in the initial velocity and overall efficiency of the m3C demethylation by the I168A mutant form was well pronounced but not as dramatic as that seen with substitutions Tyr122 or Asp173.
It is interesting to note that the final amount of the overall product accumulation for the WT ABH2 interaction with m3C-containing DNA was higher than that for m1A-containing DNA, but the rate of the accumulation was slower. The introduction of one of three amino acid substitutions caused a greater decrease of enzymatic activity toward lesion m3C than toward m1A.
Based on the obtained time courses and the EMSA data, the observed rate constants kobs were determined with the help of Equation (2). The kobs values characterizing the catalytic repair of the m1A- and m3C-containing DNA substrates by WT ABH2 or its mutant forms Y122A, I168A, or D173A are presented as a diagram in Figure 11 and are listed in Table 2. Kinetic analysis revealed that, for the ABH2 I168A mutant form, the kobs values characterizing its catalytic activity toward m1A- and m3C-containing DNA substrates were 4- and 6-fold lower, respectively, as compared to the WT enzyme (Table 2). The kobs values of the m1A demethylation reactions depended on the enzyme and were in the range between 0.010 s−1 for ABH2 D173A and 0.071 s−1 for WT ABH2. As expected, the kobs values for D173A and Y122A (0.013 s−1) variants were very close. Just as in the case of m3C, the kobs values for D173A and Y122A were the lowest in the range, and their difference was within the margin of error (0.0045 ± 0.0009 and 0.005 ± 0.001 s−1, respectively).

3. Conclusions

Mutational analysis of Tyr122, Ile168, and Asp173 residues, which form contacts with damaged bases inside the active site of ABH2, was performed by combining of computational and experimental approaches. These amino acid residues form direct contacts with the m1A or m3C lesions and, thereby, allow us to estimate their possible role in the recognition of the epigenetic mark m5C. Molecular dynamics simulations revealed the impact of these amino acid residues in the catalytic process. In the case of Y122A, a loss of the hydrogen bond between Tyr122 and Arg254 was found, resulting in the disappearance of the salt bridge Arg254–Asp173, which is important for the coordination of the Fe2+ ion in the active site by Asp173. The substitution of I168A allowed an αKG molecule to move inside the active site, thereby providing a greater variety of αKG positions. D173A substitutions lead to an increase of the mobility of damaged bases in the active-site pocket, owing to the disturbance of the stacking of m1A or m3C with the Phe124 residue. The MD models of WT ABH2 with DNA containing an m5C base revealed that the m5C base could be positioned in the active site in alternative “noncatalytic” and “catalytic” ways. However, in the “catalytic” position, when the methyl group of m5C is close to the catalytic amino acids, there is a loss of many contacts that hold the everted base inside the active site.
Experimental EMSA analysis of the binding of WT ABH2 or its mutant forms to methylated DNA revealed only a minor influence of all the tested substitutions on the DNA-binding ability. Indeed, the obtained Kd values were varied between the enzymes within 1.3-fold and 1.6-fold for the DNA containing m1A or m3C, respectively. However, the analysis of the product accumulation data revealed that, in the case of mutant forms, the rate constant of the product accumulation kobs was significantly lower as compared to the WT enzyme. Therefore, the obtained data clearly demonstrate the effect of the tested amino acid residues on the catalytic activity of the enzymes rather than the DNA-binding ability. Taken together, these data shed light on the molecular and kinetic consequences of the substitution of active-site residues for the mechanism of the substrate recognition.

4. Materials and Methods

4.1. Oligodeoxyribonucleotides

The synthesis of the oligodeoxyribonucleotides (Table 3) was performed on an ASM-800DNA/RNA synthesizer (Biosset, Novosibirsk, Russia) by means of standard commercial phosphoramidites and CPG solid supports from Glen Research (Sterling, VA, USA). The oligonucleotides were deprotected according to the manufacturer’s protocols and were purified by high-performance liquid chromatography. Oligonucleotide homogeneity was checked by denaturing 20% PAGE. Concentrations of oligonucleotides were calculated from their absorbance at 260 nm (A260). Oligonucleotide duplexes were prepared via the annealing of oligonucleotide strands at a 1:1 molar ratio.

4.2. Site-Directed Mutagenesis

This procedure was used to mutate the ABH2 gene at sites of interest. Primers were synthetized as described above. Three mutant forms were designed by (1) the substitution of a tyrosine residue with alanine at position 122, (2) the substitution of an isoleucine residue with alanine at position 168, or (3) the substitution of an aspartate residue with alanine at position 173. Sequences of the forward primers in question were as follows (reverse primers were complementary to forward primers): 5′-GCTGGGCTGACCGCAACATTTTCAGG-3′ (Y122A), 5′-GATGGCTGTGACCACGCCGGGGAGCAC-3′ (I168A), and 5′-CGGGGAGCACCGAGCAGATGAAAGAGAACTG-3′ (D173A). After PCR, the methylated DNA was removed by MalI (SibEnzyme, Novosibirsk, Russia) digestion. These mixtures were transfected into the Escherichia coli strain ElectroMAX™ DH10B (Invitrogen, Waltham, MA, USA). Plasmid DNA was isolated from the resulting cell culture according to the protocol of the diaGene Isolation Kit (Novosibirsk, Russia). The concentration of the isolated plasmid DNA was measured on a NanoPhotometer (Implen, Munich, Germany), and its length was determined in an agarose gel. Sequences of the transformed clones were confirmed by direct sequencing.

4.3. Enzyme Purification

The WT ABH2 enzyme (ABH2_wt) and its mutant forms featuring substitutions Y122A, I168A, or D173A were isolated from E. coli ArcticExpress (DE3) cells (Invitrogen, Waltham, MA, USA) transformed with plasmid pET28c encoding an N-terminal His6 tag and the relevant gene. To purify the enzymes expressed as recombinant proteins, 1 L of culture (in YT broth) of E. coli cells carrying the desired vector construct was grown with 50 µg/mL kanamycin, 20 µg/mL gentamycin, and 10 µg/mL tetracycline at 37 °C until A600 reached 0.6–0.7; the expression of the enzyme was induced during 24 h at 16 °C with 0.2 mM isopropyl β-D-1-thiogalactopyranoside. The cells were harvested by centrifugation (5000× g, 20 min) and then resuspended in a buffer (20 mM HEPES-KOH (pH of 7.8), 40 mM NaCl) with the addition of a mixture of EDTA-free protease inhibitors (Inhibitor cocktail, Complete, Mannheim, Germany) followed by cell lysis by means of a French press. All the purification procedures were carried out at 4 °C. The homogenate was centrifuged at 40,000× g for 45 min. The supernatant was then passed through a 0.45 µm filter (Labfil, ALWSCI, Zhejiang, China) and the NaCl concentration in the supernatant was brought to 500 mM, whereupon the supernatant was mixed with Ni-Sepharose HP (Cytiva, GE Healthcare Life Sciences, Marlborough, MA, USA). The elution of the protein-containing fraction was carried out with a buffer composed of 500 mM NaCl, 20 mM HEPES-NaOH (pH of 7.8), and 600 mM imidazole. Then, the NaCl concentration in the eluent was brought to 150 mM by dilution with 20 mM HEPES-NaOH (pH of 7.8). The resulting solution was loaded on a 1 mL HiTrap-SP HPTM column (Cytiva GE Healthcare Life Sciences, United States). Bound protein was eluted with a linear 150→1200 mM gradient of NaCl. The homogeneity of the protein was verified by SDS-PAGE (Figure 12); the protein concentration was measured on the NanoPhotometer (Implen, Munich, Germany). Glycerol was added to the isolated fractions up to 50% (v/v). Enzyme stock solutions were stored at −70 °C.

4.4. CD Spectra

These spectra were recorded on a Jasco J-600 spectropolarimeter (Jacso, Tokyo, Japan) at 25 °C in quartz cells with a 1 cm light path length. The concentration of ABH2 in the device cell was 20 µM. The experiments were conducted in a buffer consisting of 50 mM Tris-HCl with a pH of 8.0. The spectra were acquired at a bandwidth of 1.0 nm and a resolution of 1.0 nm at a scan speed of 50 nm/min. The scans were accumulated and automatically averaged.

4.5. The EMSA

WT ABH2, or one of its mutant forms, was first diluted with EMSA buffer (50 mM Tris-HCl with a pH of 8.0, 50 mM KCl, 10 mM MgCl2, 0.5 mM αKG, and 2 mM sodium ascorbate). In parallel, 6-carboxyfluorescein (FAM)-5′-labeled oligodeoxyribonucleotides containing m1A or m3C as damage were heated with the respective complementary strand at 95 °C for 3 min, then cooled to 22 °C. The binding reactions were initiated by adding 5 μL of the DNA sample to 5 μL of the protein samples. Standard binding reactions contained 1.2 µM of FAM-labeled damaged DNA duplex with protein concentrations between 0.4 and 12.8 µM. The reaction mixtures were incubated for 10 min at 37 °C, then cooled for 5 min on ice. After the glycerol concentration in a sample was brought to 5%, the reaction mixtures were loaded on a nondenaturing 10% polyacrylamide gel (70:1 acrylamide:bis-acrylamide) in 0.5× TBE buffer (1× TBE buffer: 90 mM Tris, 90 mM boric acid, 2.4 mM EDTA), and electrophoresis was run for approximately 35 min at 200 V in a cool box (at ~4–8 °C). The gels were visualized using a Molecular Imager ®VersaDoc™ MP Imaging System (Bio-Rad, Hercules, CA, USA) and the bands were quantified by scanning densitometry in the Gel-Pro Analyzer software, v.4.0 (Media Cybernetics, Rockville, MD, USA). The fraction of bound DNA [F = DNAbound/(DNAbound + DNAunbound)] was plotted against the protein concentration, and the binding data were fitted to the Hill equation (Equation (1)) with the OriginPro 8 software (OriginLab, Northampton, MA, USA). Titration series (nine points) were implemented in triplicate.
F = F u + F b − F u 1 + K d E 0 h
where F is the fraction of bound DNA; Fu and Fb are the initial and final levels of F, respectively; h is the Hill coefficient; Kd is the effective dissociation constant of the complex.

4.6. Methylated DNA Repair Activity Assays

FAM-5′-labeled oligonucleotides were subjected to experiments on the separation of the cleavage products by PAGE. The repair activity assays were carried out in reaction buffer (50 mM Tris-HCl (pH of 8.0), 50 mM KCl, 10.0 mM MgCl2, 1.0 mM αKG, 2.0 mM sodium ascorbate, and 40 μM (NH4)2Fe(SO4)2·6H2O). Analyses of the active demethylation of the DNA substrates containing m1A or m3C as a lesion were performed at 37 °C in a 10 µL reaction mixture for 30 min. The substrate concentration was 1.0 µM and the concentration of the enzyme was 1.0 µM as well. The reaction was initiated by the addition of the enzyme. After 30 min, the reaction mixture was subjected to precipitation with a 2% solution of LiClO4 in acetone.
In the analysis of the kinetics of the active demethylation product accumulation, the reaction between m1A- or m3C-containing DNA substrates and the WT enzyme or mutant forms was stopped after a certain period ranging from 30 s to 1 h.
To visualize the product of DNA demethylation, each precipitated sample was treated with 5 U of restriction endonuclease BstMB I (SibEnzyme, Novosibirsk, Russia), which recognizes the GATC sequence and is sensitive to DNA methylation. The digestion of the ABH2 reaction product by BstMB I was carried out for 30 min at 37 °C. The reaction with BstMB I was quenched with 5 µL of a gel-loading dye containing 7 M urea and 50 mM EDTA; then, the samples were loaded onto 20% (w/v) polyacrylamide gel containing 7 M urea. PAGE was performed at 55 °C and a voltage of 200–300 V. The gels were visualized using the Molecular Imager VersaDoc™ MP Imaging System (Bio-Rad, Hercules, CA, USA), and the bands were quantified by scanning densitometry in Gel-Pro Analyzer v.4.0 (Media Cybernetics, Rockville, MD, USA).
The data were analyzed via the following equation, under the assumption that an observed catalytic rate constant is a ratio of the initial velocity, V0, to the equilibrium concentration of the precatalytic complex [E·S]:
k o b s = 2 V 0 E 0 + S 0 + K d + E 0 + S 0 + K d 2 − 4 E 0 S 0 − 1
where V0 is initial velocity estimated as the initial slope of the kinetic curve obtained under steady-state reaction conditions, E0 and S0 are the total concentrations of the enzyme and DNA, and Kd is an equilibrium dissociation constant calculated in the EMSA assay by means of Equation (1).

4.7. MD Simulations

Initial simulation structures were derived from available crystal structures of ABH2 cross-linked with DNA containing modified nucleotides everted into the enzyme active site in the presence of cofactor αKG and the Mn2+ ion: PDB ID: 3RZJ for m3C- and m5C-containing models and PDB ID: 3BUC and 3RZK for m1A models [32,39]. Poorly resolved amino acid residues were filled in using Modeller [53], cross-linking modifications were removed, and DNA duplex sequences were adjusted to match the experimental ones.
Amino acid sidechain protonation states were assigned on the H++ server [54]. Simulation setup and simulations were performed using the GROMACS MD package [55]. The starting structures were solvated and neutralized in a dodecahedral PBC box by means of TIP3P model water and 50 mM KCl JC ions [56,57]. An AMBER 14SB forcefield with OL15 corrections was used to describe the protein and the substrate DNA [58,59,60,61]. The AMBER forcefield parameters for αKG, m3C, and m1A were generated by Antechamber, with RESP charges calculated by the REDD server [62,63,64,65,66]. Conversion to GROMACS-compatible formats was performed with the ACPYPE tool [67]. The active-site Fe2+ ion was modeled as a nonbonded cationic dummy model [68].
The cutoff of nonbonded interactions was set to 1.0 nm, and long-range electrostatic interactions were processed by the PME method [69,70]. Covalent bonds involving hydrogen atoms were constrained using the LINCS solver [71]. Steepest-descent energy minimization was followed by 1 ns NVT and NPT equilibrations with restrained positions of all solute heavy atoms, with the help of the Bussi thermostat and Parrinello–Rahman barostat [72,73]. Postequilibration unrestrained MD simulations were run for 100 ns in triplicate, with one set each for m3C and m1A, and were extended up to 500 ns. Trajectory processing was performed using the integrated GROMACS toolset. Images were generated in the open-source version of the PyMOL viewer.

Author Contributions

Conceptualization, N.A.K.; methodology, A.T.D. and T.E.T.; validation, A.T.D., T.E.T., A.A.K. and N.A.K.; formal analysis, A.T.D. and T.E.T.; investigation, A.T.D. and T.E.T.; resources, A.A.I., M.S., A.A.K. and N.A.K.; data curation, A.T.D., T.E.T., M.Z., A.A.I., M.S., A.A.K. and N.A.K.; writing—original draft preparation, A.T.D., T.E.T. and N.A.K.; writing—review and editing, A.T.D., T.E.T. and N.A.K.; supervision, A.A.K. and N.A.K.; project administration, A.A.K.; funding acquisition, A.A.I., M.S., A.A.K. and N.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the Russian Federal Ministry of Science and Higher Education (project no. 121031300041-4 for N.A.K.), by Electricité de France (grant nos. RB 2020-02 and RB 2021-05 to M.S.), by the French National Research Agency (grant no. ANR-18-CE44-0008 to A.A.I.), and Fondation ARC (grant no. PJA-2021060003796 to A.A.I.). The part of the work involving the analysis of the ABH2 conformational dynamics and activity was specifically funded by the Russian Science Foundation grant no. 21-14-00018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Experimental data are available upon request to N.A.K. Tel. +7-383-363-5174, E-mail: [email protected].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sedgwick, B. Repairing DNA-Methylation Damage. Nat. Rev. Mol. Cell Biol. 2004, 5, 148–157. [Google Scholar] [CrossRef]
  2. Mishina, Y.; Duguid, E.M.; He, C. Direct Reversal of DNA Alkylation Damage. Chem. Rev. 2006, 106, 215–232. [Google Scholar] [CrossRef]
  3. Hecht, S.S. DNA Adduct Formation from Tobacco-Specific N-Nitrosamines. Mutat. Res. 1999, 424, 127–142. [Google Scholar] [CrossRef] [PubMed]
  4. Rajski, S.R.; Williams, R.M. DNA Cross-Linking Agents as Antitumor Drugs. Chem. Rev. 1998, 98, 2723–2796. [Google Scholar] [CrossRef]
  5. Hurley, L.H. DNA and Its Associated Processes as Targets for Cancer Therapy. Nat. Rev. Cancer 2002, 2, 188–200. [Google Scholar] [CrossRef] [PubMed]
  6. Walport, L.J.; Hopkinson, R.J.; Schofield, C.J. Mechanisms of Human Histone and Nucleic Acid Demethylases. Curr. Opin. Chem. Biol. 2012, 16, 525–534. [Google Scholar] [CrossRef] [PubMed]
  7. Cantoni, G.L. The Nature of the Active Methyl Donor Formed Enzymatically from L-Methionine and Adenosinetriphosphate. J. Am. Chem. Soc. 1952, 74, 2942–2943. [Google Scholar] [CrossRef]
  8. Lindahl, T.; Sedgwick, B.; Sekiguchi, M.; Nakabeppu, Y. Regulation and Expression of the Adaptive Response to Alkylating Agents. Annu. Rev. Biochem. 1988, 57, 133–157. [Google Scholar] [CrossRef] [PubMed]
  9. Kuznetsov, N.A.; Kanazhevskaya, L.Y.; Fedorova, O.S. DNA Demethylation in the Processes of Repair and Epigenetic Regulation Performed by 2-Ketoglutarate-Dependent DNA Dioxygenases. Int. J. Mol. Sci. 2021, 22, 10540. [Google Scholar] [CrossRef]
  10. Trewick, S.C.; Henshaw, T.F.; Hausinger, R.P.; Lindahl, T.; Sedgwick, B. Oxidative Demethylation by Escherichia coli AlkB Directly Reverts DNA Base Damage. Nature 2002, 419, 174–178. [Google Scholar] [CrossRef]
  11. Falnes, P.Ø.; Johansen, R.F.; Seeberg, E. AlkB-Mediated Oxidative Demethylation Reverses DNA Damage in Escherichia coli. Nature 2002, 419, 178–182. [Google Scholar] [CrossRef] [PubMed]
  12. Teo, I.; Sedgwick, B.; Kilpatrick, M.W.; McCarthy, T.V.; Lindahl, T. The Intracellular Signal for Induction of Resistance to Alkylating Agents in E. coli. Cell 1986, 45, 315–324. [Google Scholar] [CrossRef] [PubMed]
  13. Aravind, L.; Koonin, E.v. The DNA-Repair Protein AlkB, EGL-9, and Leprecan Define New Families of 2-Oxoglutarate- and Iron-Dependent Dioxygenases. Genome Biol. 2001, 2, RESEARCH0007. [Google Scholar] [CrossRef]
  14. Delaney, J.C.; Smeester, L.; Wong, C.; Frick, L.E.; Taghizadeh, K.; Wishnok, J.S.; Drennan, C.L.; Samson, L.D.; Essigmann, J.M. AlkB Reverses Etheno DNA Lesions Caused by Lipid Oxidation in Vitro and in Vivo. Nat. Struct. Mol. Biol. 2005, 12, 855–860. [Google Scholar] [CrossRef]
  15. Delaney, J.C.; Essigmann, J.M. Mutagenesis, Genotoxicity, and Repair of 1-Methyladenine, 3-Alkylcytosines, 1-Methylguanine, and 3-Methylthymine in AlkB Escherichia Coli. Proc. Natl. Acad. Sci. USA 2004, 101, 14051–14056. [Google Scholar] [CrossRef] [PubMed]
  16. Frick, L.E.; Delaney, J.C.; Wong, C.; Drennan, C.L.; Essigmann, J.M. Alleviation of 1,N6-Ethanoadenine Genotoxicity by the Escherichia Coli Adaptive Response Protein AlkB. Proc. Natl. Acad. Sci. USA 2007, 104, 755–760. [Google Scholar] [CrossRef]
  17. Yi, C.; Yang, C.-G.; He, C. A Non-Heme Iron-Mediated Chemical Demethylation in DNA and RNA. Acc. Chem. Res. 2009, 42, 519–529. [Google Scholar] [CrossRef] [Green Version]
  18. Duncan, T.; Trewick, S.C.; Koivisto, P.; Bates, P.A.; Lindahl, T.; Sedgwick, B. Reversal of DNA Alkylation Damage by Two Human Dioxygenases. Proc. Natl. Acad. Sci. USA 2002, 99, 16660–16665. [Google Scholar] [CrossRef]
  19. Aas, P.A.; Otterlei, M.; Falnes, P.O.; Vågbø, C.B.; Skorpen, F.; Akbari, M.; Sundheim, O.; Bjørås, M.; Slupphaug, G.; Seeberg, E.; et al. Human and Bacterial Oxidative Demethylases Repair Alkylation Damage in Both RNA and DNA. Nature 2003, 421, 859–863. [Google Scholar] [CrossRef]
  20. Lee, D.-H.; Jin, S.-G.; Cai, S.; Chen, Y.; Pfeifer, G.P.; O’Connor, T.R. Repair of Methylation Damage in DNA and RNA by Mammalian AlkB Homologues. J. Biol. Chem. 2005, 280, 39448–39459. [Google Scholar] [CrossRef] [Green Version]
  21. Dinglay, S.; Trewick, S.C.; Lindahl, T.; Sedgwick, B. Defective Processing of Methylated Single-Stranded DNA by E. coli AlkB Mutants. Genes Dev. 2000, 14, 2097–2105. [Google Scholar] [CrossRef] [PubMed]
  22. Ougland, R.; Zhang, C.-M.; Liiv, A.; Johansen, R.F.; Seeberg, E.; Hou, Y.-M.; Remme, J.; Falnes, P.Ø. AlkB Restores the Biological Function of MRNA and TRNA Inactivated by Chemical Methylation. Mol. Cell 2004, 16, 107–116. [Google Scholar] [CrossRef] [PubMed]
  23. Fu, D.; Samson, L.D.; Hübscher, U.; van Loon, B. The Interaction between ALKBH2 DNA Repair Enzyme and PCNA Is Direct, Mediated by the Hydrophobic Pocket of PCNA and Perturbed in Naturally-Occurring ALKBH2 Variants. DNA Repair Amst. 2015, 35, 13–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gilljam, K.M.; Feyzi, E.; Aas, P.A.; Sousa, M.M.L.; Müller, R.; Vågbø, C.B.; Catterall, T.C.; Liabakk, N.B.; Slupphaug, G.; Drabløs, F.; et al. Identification of a Novel, Widespread, and Functionally Important PCNA-Binding Motif. J. Cell Biol. 2009, 186, 645–654. [Google Scholar] [CrossRef] [Green Version]
  25. Falnes, P.O. Repair of 3-Methylthymine and 1-Methylguanine Lesions by Bacterial and Human AlkB Proteins. Nucleic Acids Res. 2004, 32, 6260–6267. [Google Scholar] [CrossRef] [Green Version]
  26. Ringvoll, J.; Moen, M.N.; Nordstrand, L.M.; Meira, L.B.; Pang, B.; Bekkelund, A.; Dedon, P.C.; Bjelland, S.; Samson, L.D.; Falnes, P.Ø.; et al. AlkB Homologue 2-Mediated Repair of Ethenoadenine Lesions in Mammalian DNA. Cancer Res. 2008, 68, 4142–4149. [Google Scholar] [CrossRef] [Green Version]
  27. Ringvoll, J.; Nordstrand, L.M.; Vagbo, C.B.; Talstad, V.; Reite, K.; Aas, P.A.; Lauritzen, K.H.; Liabakk, N.B.; Bjork, A.; Doughty, R.W.; et al. Repair Deficient Mice Reveal MABH2 as the Primary Oxidative Demethylase for Repairing 1meA and 3meC Lesions in DNA. Embo J. 2006, 25, 2189–2198. [Google Scholar] [CrossRef] [Green Version]
  28. Bian, K.; Lenz, S.A.P.; Tang, Q.; Chen, F.; Qi, R.; Jost, M.; Drennan, C.L.; Essigmann, J.M.; Wetmore, S.D.; Li, D. DNA Repair Enzymes ALKBH2, ALKBH3, and AlkB Oxidize 5-Methylcytosine to 5-Hydroxymethylcytosine, 5-Formylcytosine and 5-Carboxylcytosine in Vitro. Nucleic Acids Res. 2019, 47, 5522–5529. [Google Scholar] [CrossRef] [Green Version]
  29. Seethy, A.; Pethusamy, K.; Chattopadhyay, I.; Sah, R.; Chopra, A.; Dhar, R.; Karmakar, S. TETology: Epigenetic Mastermind in Action. Appl. Biochem. Biotechnol. 2021, 193, 1701–1726. [Google Scholar] [CrossRef]
  30. Ito, S.; Shen, L.; Dai, Q.; Wu, S.C.; Collins, L.B.; Swenberg, J.A.; He, C.; Zhang, Y. Tet Proteins Can Convert 5-Methylcytosine to 5-Formylcytosine and 5-Carboxylcytosine. Science 2011, 333, 1300–1303. [Google Scholar] [CrossRef] [Green Version]
  31. Yu, B.; Edstrom, W.C.; Benach, J.; Hamuro, Y.; Weber, P.C.; Gibney, B.R.; Hunt, J.F. Crystal Structures of Catalytic Complexes of the Oxidative DNA/RNA Repair Enzyme AlkB. Nature 2006, 439, 879–884. [Google Scholar] [CrossRef]
  32. Yang, C.-G.; Yi, C.; Duguid, E.M.; Sullivan, C.T.; Jian, X.; Rice, P.A.; He, C. Crystal Structures of DNA/RNA Repair Enzymes AlkB and ABH2 Bound to DsDNA. Nature 2008, 452, 961–965. [Google Scholar] [CrossRef] [Green Version]
  33. Sundheim, O.; Talstad, V.A.; Vågbø, C.B.; Slupphaug, G.; Krokan, H.E. AlkB Demethylases Flip out in Different Ways. DNA Repair (Amst.) 2008, 7, 1916–1923. [Google Scholar] [CrossRef]
  34. Holland, P.J.; Hollis, T. Structural and Mutational Analysis of Escherichia Coli AlkB Provides Insight into Substrate Specificity and DNA Damage Searching. PLoS ONE 2010, 5, e8680. [Google Scholar] [CrossRef] [Green Version]
  35. Sundheim, O.; Vågbø, C.B.; Bjørås, M.; Sousa, M.M.L.; Talstad, V.; Aas, P.A.; Drabløs, F.; Krokan, H.E.; Tainer, J.A.; Slupphaug, G. Human ABH3 Structure and Key Residues for Oxidative Demethylation to Reverse DNA/RNA Damage. EMBO J. 2006, 25, 3389–3397. [Google Scholar] [CrossRef] [Green Version]
  36. Giri, N.C.; Sun, H.; Chen, H.; Costa, M.; Maroney, M.J. X-Ray Absorption Spectroscopy Structural Investigation of Early Intermediates in the Mechanism of DNA Repair by Human ABH2. Biochemistry 2011, 50, 5067–5076. [Google Scholar] [CrossRef] [Green Version]
  37. Monsen, V.T.; Sundheim, O.; Aas, P.A.; Westbye, M.P.; Sousa, M.M.L.; Slupphaug, G.; Krokan, H.E. Divergent SS-Hairpins Determine Double-Strand versus Single-Strand Substrate Recognition of Human AlkB-Homologues 2 and 3. Nucleic Acids Res. 2010, 38, 6447–6455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Chen, B.; Liu, H.; Sun, X.; Yang, C.-G. Mechanistic Insight into the Recognition of Single-Stranded and Double-Stranded DNA Substrates by ABH2 and ABH3. Mol. Biosyst. 2010, 6, 2143–2149. [Google Scholar] [CrossRef] [PubMed]
  39. Yi, C.; Chen, B.; Qi, B.; Zhang, W.; Jia, G.; Zhang, L.; Li, C.J.; Dinner, A.R.; Yang, C.-G.; He, C. Duplex Interrogation by a Direct DNA Repair Protein in Search of Base Damage. Nat. Struct. Mol. Biol. 2012, 19, 671–676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Waheed, S.O.; Ramanan, R.; Chaturvedi, S.S.; Lehnert, N.; Schofield, C.J.; Christov, C.Z.; Karabencheva-Christova, T.G. Role of Structural Dynamics in Selectivity and Mechanism of Non-Heme Fe(II) and 2-Oxoglutarate-Dependent Oxygenases Involved in DNA Repair. ACS Cent. Sci. 2020, 6, 795–814. [Google Scholar] [CrossRef]
  41. Kanazhevskaya, L.Y.; Smyshliaev, D.A.; Timofeyeva, N.A.; Ishchenko, A.A.; Saparbaev, M.; Kuznetsov, N.A.; Fedorova, O.S. Conformational Dynamics of Human ALKBH2 Dioxygenase in the Course of DNA Repair as Revealed by Stopped-Flow Fluorescence Spectroscopy. Molecules 2022, 27, 4960. [Google Scholar] [CrossRef]
  42. Kuznetsova, A.A.; Kuznetsov, N.A.; Ishchenko, A.A.; Saparbaev, M.K.; Fedorova, O.S. Step-by-Step Mechanism of DNA Damage Recognition by Human 8-Oxoguanine DNA Glycosylase. Biochim. Biophys. Acta 2014, 1840, 387–395. [Google Scholar] [CrossRef]
  43. Yakovlev, D.A.; Kuznetsova, A.A.; Fedorova, O.S.; Kuznetsov, N.A. Search for Modified DNA Sites with the Human Methyl-CpG-Binding Enzyme MBD4. Acta Naturae 2017, 9, 88–98. [Google Scholar] [CrossRef]
  44. Kuznetsova, A.A.; Fedorova, O.S.; Kuznetsov, N.A. Lesion Recognition and Cleavage of Damage-Containing Quadruplexes and Bulged Structures by DNA Glycosylases. Front. Cell Dev. Biol. 2020, 1462. [Google Scholar] [CrossRef] [PubMed]
  45. Kuznetsov, N.A.; Kuznetsova, A.A.; Vorobjev, Y.N.; Krasnoperov, L.N.; Fedorova, O.S. Thermodynamics of the DNA Damage Repair Steps of Human 8-Oxoguanine DNA Glycosylase. PLoS ONE 2014, 9, e98495. [Google Scholar] [CrossRef] [Green Version]
  46. Tyugashev, T.E.; Vorobjev, Y.N.; Kuznetsova, A.A.; Lukina, M.V.; Kuznetsov, N.A.; Fedorova, O.S. Roles of Active-Site Amino Acid Residues in Specific Recognition of DNA Lesions by Human 8-Oxoguanine-DNA Glycosylase (OGG1). J. Phys. Chem. B 2019, 123, 4878–4887. [Google Scholar] [CrossRef]
  47. Kuznetsov, N.A.; Kiryutin, A.S.; Kuznetsova, A.A.; Panov, M.S.; Barsukova, M.O.; Yurkovskaya, A.V.; Fedorova, O.S. The Formation of Catalytically Competent Enzyme-Substrate Complex Is Not a Bottleneck in Lesion Excision by Human Alkyladenine DNA Glycosylase. J. Biomol. Struct. Dyn. 2017, 35, 950–967. [Google Scholar] [CrossRef]
  48. Bulygin, A.A.; Kuznetsova, A.A.; Vorobjev, Y.N.; Fedorova, O.S.; Kuznetsov, N.A. The Role of Active-Site Plasticity in Damaged-Nucleotide Recognition by Human Apurinic/Apyrimidinic Endonuclease APE1. Molecules 2020, 25, 3940. [Google Scholar] [CrossRef] [PubMed]
  49. Kuznetsova, A.A.; Senchurova, S.I.; Ishchenko, A.A.; Saparbaev, M.; Fedorova, O.S.; Kuznetsov, N.A. Common Kinetic Mechanism of Abasic Site Recognition by Structurally Different Apurinic/Apyrimidinic Endonucleases. Int. J. Mol. Sci. 2021, 22, 8874. [Google Scholar] [CrossRef] [PubMed]
  50. Kuznetsova, A.A.; Matveeva, A.G.; Milov, A.D.; Vorobjev, Y.N.; Dzuba, S.A.; Fedorova, O.S.; Kuznetsov, N.A. Substrate Specificity of Human Apurinic/Apyrimidinic Endonuclease APE1 in the Nucleotide Incision Repair Pathway. Nucleic Acids Res. 2018, 46, 11454–11465. [Google Scholar] [CrossRef] [Green Version]
  51. Senchurova, S.I.; Syryamina, V.N.; Kuznetsova, A.A.; Novopashina, D.S.; Ishchenko, A.A.; Saparbaev, M.; Dzuba, S.A.; Fedorova, O.S.; Kuznetsov, N.A. The Mechanism of Damage Recognition by Apurinic/Apyrimidinic Endonuclease Nfo from Escherichia coli. BBA General. Subj. 2022, 1866, 130216. [Google Scholar] [CrossRef]
  52. Kuznetsov, N.A.; Fedorova, O.S. Thermodynamic Analysis of Fast Stages of Specific Lesion Recognition by DNA Repair Enzymes. Biochemistry Mosc. 2016, 81, 1136–1152. [Google Scholar] [CrossRef]
  53. Šali, a MODELLER: A Program for Protein Structure Modeling Release 9.12, R9480; Rockefeller University: New York, NY, USA, 2013.
  54. Anandakrishnan, R.; Aguilar, B.; Onufriev, A.V. H++ 3.0: Automating PK Prediction and the Preparation of Biomolecular Structures for Atomistic Molecular Modeling and Simulations. Nucleic Acids Res. 2012, 40, W537–W541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindah, E. Gromacs: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef] [Green Version]
  56. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926. [Google Scholar] [CrossRef]
  57. Joung, I.S.; Cheatham, T.E. Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations. J. Phys. Chem. B 2008, 112, 9020–9041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Cornell, W.D.; Cieplak, P.; Bayly, C.I.; Gould, I.R.; Merz, K.M.; Ferguson, D.M.; Spellmeyer, D.C.; Fox, T.; Caldwell, J.W.; Kollman, P.A. A Second Generation All Atom Force Field for the Simulation of Proteins, Nucleic Acids and Organic Molecules. J. Am. Chem. Soc. 1995, 117, 5179–5197. [Google Scholar] [CrossRef] [Green Version]
  59. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. Ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from Ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  60. Zgarbová, M.; Otyepka, M.; Šponer, J.; Mládek, A.; Banáš, P.; Cheatham, T.E.; Jurečka, P. Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles. J. Chem. Theory Comput. 2011, 7, 2886–2902. [Google Scholar] [CrossRef]
  61. Zgarbová, M.; Šponer, J.; Otyepka, M.; Cheatham, T.E.; Galindo-Murillo, R.; Jurečka, P. Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA. J. Chem. Theory Comput. 2015, 11, 5723–5736. [Google Scholar] [CrossRef]
  62. Bayly, C.I.; Cieplak, P.; Cornell, W.D.; Kollman, P.A. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges: The RESP Model. J. Phys. Chem. 1993, 97, 10269–10280. [Google Scholar] [CrossRef]
  63. Wang, J.M.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General Amber Force Field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  64. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic Atom Type and Bond Type Perception in Molecular Mechanical Calculations. J. Mol. Graph. Model. 2006, 25, 247–260. [Google Scholar] [CrossRef] [PubMed]
  65. Meagher, K.L.; Redman, L.T.; Carlson, H.A. Development of Polyphosphate Parameters for Use with the AMBER Force Field. J. Comput. Chem. 2003, 24, 1016–1025. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Vanquelef, E.; Simon, S.; Marquant, G.; Garcia, E.; Klimerak, G.; Delepine, J.C.; Cieplak, P.; Dupradeau, F.Y. RED Server: A Web Service for Deriving RESP and ESP Charges and Building Force Field Libraries for New Molecules and Molecular Fragments. Nucleic Acids Res. 2011, 39, W511–W517. [Google Scholar] [CrossRef] [Green Version]
  67. Sousa da Silva, A.W.; Vranken, W.F. ACPYPE—AnteChamber PYthon Parser InterfacE. BMC Res. Notes 2012, 5, 367. [Google Scholar] [CrossRef] [Green Version]
  68. Jiang, Y.; Zhang, H.; Tan, T. Rational Design of Methodology-Independent Metal Parameters Using a Nonbonded Dummy Model. J. Chem. Theory Comput. 2016, 12, 3250–3260. [Google Scholar] [CrossRef] [Green Version]
  69. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577. [Google Scholar] [CrossRef] [Green Version]
  70. Wennberg, C.L.; Murtola, T.; Hess, B.; Lindahl, E. Lennard-Jones Lattice Summation in Bilayer Simulations Has Critical Effects on Surface Tension and Lipid Properties. J. Chem. Theory Comput. 2013, 9, 3527–3537. [Google Scholar] [CrossRef]
  71. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  72. Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1998, 52, 7182. [Google Scholar] [CrossRef]
Cells 12 01839 g001 550
Figure 1. An overview of the structure of a complex of WT ABH2 with DNA containing m3C (Protein Data Bank (PDB) ID 3RZJ). Functionally important loops and elements are shown: an intercalating loop containing Val101 and Phe102 (red); a positively charged loop containing Arg241, Lys242, and Lys243 (orange, pointed with an orange arrow); a flexible unresolved loop carrying DNA-binding residues Arg198, Gly204, and Lys205 (gray, pointed with a gray arrow); a metal ion (dark blue); the flipped damaged base inside the active site (blue); αKG (yellow).
Figure 1. An overview of the structure of a complex of WT ABH2 with DNA containing m3C (Protein Data Bank (PDB) ID 3RZJ). Functionally important loops and elements are shown: an intercalating loop containing Val101 and Phe102 (red); a positively charged loop containing Arg241, Lys242, and Lys243 (orange, pointed with an orange arrow); a flexible unresolved loop carrying DNA-binding residues Arg198, Gly204, and Lys205 (gray, pointed with a gray arrow); a metal ion (dark blue); the flipped damaged base inside the active site (blue); αKG (yellow).
Cells 12 01839 g001
Cells 12 01839 g002 550
Figure 2. Active-site architecture in the case of m1A ((a), PDB ID: 3BTY), m3C ((b), PDB ID: 3RZJ), or εA ((c), PDB ID: 3RZK). Amino acid residues coming into contact with the damaged base are shown.
Figure 2. Active-site architecture in the case of m1A ((a), PDB ID: 3BTY), m3C ((b), PDB ID: 3RZJ), or εA ((c), PDB ID: 3RZK). Amino acid residues coming into contact with the damaged base are shown.
Cells 12 01839 g002
Cells 12 01839 g003 550
Figure 3. Residues targeted by site-directed mutagenesis in the active site of ABH2 that is in complex with m1A-containing DNA ((a), PDB ID: 3BTY) or m3C-containing DNA ((b), PDB ID: 3RZJ).
Figure 3. Residues targeted by site-directed mutagenesis in the active site of ABH2 that is in complex with m1A-containing DNA ((a), PDB ID: 3BTY) or m3C-containing DNA ((b), PDB ID: 3RZJ).
Cells 12 01839 g003
Cells 12 01839 g004a 550Cells 12 01839 g004b 550
Figure 4. Close-up view of m3C (a,c,e,g) and m1A (b,d,f,h) damage-recognition workspace in the active site of WT ABH2 (a,b) or its mutant forms Y122A (c,d), I168A (e,f), or D173A (g,h), as determined by MD simulations.
Figure 4. Close-up view of m3C (a,c,e,g) and m1A (b,d,f,h) damage-recognition workspace in the active site of WT ABH2 (a,b) or its mutant forms Y122A (c,d), I168A (e,f), or D173A (g,h), as determined by MD simulations.
Cells 12 01839 g004aCells 12 01839 g004b
Cells 12 01839 g005 550
Figure 5. Close-up view of possible positions (“noncatalytic” (a) and “catalytic” (b)) of an m5C base in the active site of WT ABH2.
Figure 5. Close-up view of possible positions (“noncatalytic” (a) and “catalytic” (b)) of an m5C base in the active site of WT ABH2.
Cells 12 01839 g005
Cells 12 01839 g006 550
Figure 6. CD spectra of WT ABH2 (green) and of mutant forms featuring the Y122A (pink), I168A (orange), or D173A (purple) substitutions.
Figure 6. CD spectra of WT ABH2 (green) and of mutant forms featuring the Y122A (pink), I168A (orange), or D173A (purple) substitutions.
Cells 12 01839 g006
Cells 12 01839 g007 550
Figure 7. The EMSA of the WT ABH2 enzyme and its mutant forms Y122A, I168A, and D173A. (a) EMSA of WT ABH2 or of Y122A, I168A, or D173A with 1.2 µM FAM-labeled m1A-containing DNA. (b) EMSA of WT ABH2 or Y122A, I168A, or D173A with 1.2 µM FAM-labeled m3C-containing DNA. (c) Fractions of bound FAM-labeled m1A-containing DNA plotted against ABH2 concentrations, with the curve fitting of the obtained data. (d) The fractions of bound FAM-labeled m3C-containing DNA plotted against ABH2 concentrations, with the curve fitting of the obtained data. In each panel, the concentration of DNA is 1.2 µM, and concentrations of the enzyme are 0, 0.4, 0.8, 1.2, 1.6, 3.2, 6.4, 9.6, and 12.8 µM.
Figure 7. The EMSA of the WT ABH2 enzyme and its mutant forms Y122A, I168A, and D173A. (a) EMSA of WT ABH2 or of Y122A, I168A, or D173A with 1.2 µM FAM-labeled m1A-containing DNA. (b) EMSA of WT ABH2 or Y122A, I168A, or D173A with 1.2 µM FAM-labeled m3C-containing DNA. (c) Fractions of bound FAM-labeled m1A-containing DNA plotted against ABH2 concentrations, with the curve fitting of the obtained data. (d) The fractions of bound FAM-labeled m3C-containing DNA plotted against ABH2 concentrations, with the curve fitting of the obtained data. In each panel, the concentration of DNA is 1.2 µM, and concentrations of the enzyme are 0, 0.4, 0.8, 1.2, 1.6, 3.2, 6.4, 9.6, and 12.8 µM.
Cells 12 01839 g007
Cells 12 01839 g008 550
Figure 8. Analysis of the binding affinity of WT ABH2 or its mutant forms Y122A, I168A, and D173A in an EMSA. Dissociation constants (Kd) were determined by the fitting of the experimental data points to the Hill equation (Equation (1)). Error bars represent standard deviation of three technical replicates.
Figure 8. Analysis of the binding affinity of WT ABH2 or its mutant forms Y122A, I168A, and D173A in an EMSA. Dissociation constants (Kd) were determined by the fitting of the experimental data points to the Hill equation (Equation (1)). Error bars represent standard deviation of three technical replicates.
Cells 12 01839 g008
Cells 12 01839 g009a 550Cells 12 01839 g009b 550
Figure 9. PAGE analysis of the DNA demethylation by ABH2 enzymes during an interaction with the dsDNA substrate containing m1A (a) or m3C (b). A comparison of efficacy of demethylation of damaged DNA substrates by ABH2 enzymes during interaction with m1A-containing (c) or m3C-containing DNA (d). [Enzyme] = [DNA] = 1.0 µM. T = 37 °C, reaction time = 30 min. S is a substrate, P is a product of the cleavage of the demethylated DNA chain by the restriction enzyme. Levels of accumulation of the product are presented as the average of three experimental values ± SD.
Figure 9. PAGE analysis of the DNA demethylation by ABH2 enzymes during an interaction with the dsDNA substrate containing m1A (a) or m3C (b). A comparison of efficacy of demethylation of damaged DNA substrates by ABH2 enzymes during interaction with m1A-containing (c) or m3C-containing DNA (d). [Enzyme] = [DNA] = 1.0 µM. T = 37 °C, reaction time = 30 min. S is a substrate, P is a product of the cleavage of the demethylated DNA chain by the restriction enzyme. Levels of accumulation of the product are presented as the average of three experimental values ± SD.
Cells 12 01839 g009aCells 12 01839 g009b
Cells 12 01839 g010 550
Figure 10. Active demethylation of m1A-containing (a) or m3C-containing (b) DNA substrates by WT ABH2, Y122A, I168A, or D173A. [Enzyme] = [DNA] = 1.0 µM. T = 37 °C. Levels of accumulation of the product at each time point are presented as the average of three experimental values ± SD.
Figure 10. Active demethylation of m1A-containing (a) or m3C-containing (b) DNA substrates by WT ABH2, Y122A, I168A, or D173A. [Enzyme] = [DNA] = 1.0 µM. T = 37 °C. Levels of accumulation of the product at each time point are presented as the average of three experimental values ± SD.
Cells 12 01839 g010
Cells 12 01839 g011 550
Figure 11. Observed rate constants kobs characterizing the catalytic activity of WT ABH2 and its mutant forms Y122A, I168A, and D173A toward m1A- and m3C-containing DNA substrates. The kobs values were determined by means of Equation (2). Error bars represent the SD determined by mathematical data processing.
Figure 11. Observed rate constants kobs characterizing the catalytic activity of WT ABH2 and its mutant forms Y122A, I168A, and D173A toward m1A- and m3C-containing DNA substrates. The kobs values were determined by means of Equation (2). Error bars represent the SD determined by mathematical data processing.
Cells 12 01839 g011
Cells 12 01839 g012 550
Figure 12. Analysis of protein purity of the WT, Y122A, I168A, and D173A ABH2 enzymes by SDS-PAGE. The proteins were stained with Coomassie Blue R-250 (Sangon Biotech, Shanghai, China).
Figure 12. Analysis of protein purity of the WT, Y122A, I168A, and D173A ABH2 enzymes by SDS-PAGE. The proteins were stained with Coomassie Blue R-250 (Sangon Biotech, Shanghai, China).
Cells 12 01839 g012
Table
Table 1. The list of amino acid residues forming contacts with three types of damaged bases inside the active site of WT ABH2.
Table 1. The list of amino acid residues forming contacts with three types of damaged bases inside the active site of WT ABH2.
m1Am3CεA
Val99Wall of damaged-base pocketWall of damaged-base pocketWall of damaged-base pocket
Val101DNA-intercalating residueDNA-intercalating residueDNA-intercalating residue
Phe102DNA-intercalating residueDNA-intercalating residueDNA-intercalating residue
Arg110Coordination of damaged base (4.38 Å)Coordination of damaged base (3.9 Å)Coordination of damaged base (4.52 Å)
Tyr122Hydrogen bond with damaged baseHydrogen bond with damaged baseNo direct interaction with damaged base
Phe124Stacking with damaged baseStacking with damaged baseStacking with damaged base
Ser125Wall of damaged-base pocketWall of damaged-base pocketCoordination of damaged base (4.25 Å)
Ile168Wall of damaged-base pocketWall of damaged-base pocketWall of damaged-base pocket
His171Stacking with damaged baseStacking with damaged baseStacking with damaged base
Asp173Hydrogen bond with damaged baseHydrogen bond with damaged baseNo direct interaction with damaged base
Glu175Hydrogen bond with damaged baseHydrogen bond with damaged baseNo direct interaction with damaged base
Table
Table 2. A summary of kinetic parameters of DNA demethylation by ABH2 enzymes, as determined under steady-state conditions.
Table 2. A summary of kinetic parameters of DNA demethylation by ABH2 enzymes, as determined under steady-state conditions.
Kd, µMkobs, s−1
WT2.45 ± 0.070.071 ± 0.007m1A
Y122A2.11 ± 0.080.013 ± 0.006
I168A2.7 ± 0.40.02 ± 0.01
D173A2.10 ± 0.090.010 ± 0.003
WT3.5 ± 0.10.045 ± 0.007m3C
Y122A2.49 ± 0.060.005 ± 0.001
I168A4.1 ± 0.10.008 ± 0.002
D173A3.09 ± 0.070.0045 ± 0.0009
Table
Table 3. DNA substrates used in this study.
Table 3. DNA substrates used in this study.
ShorthandSequence
m1A5′-FAM-AGTTCAATG-m1A-TCTTCAT-3′
3′-TCAAGTTACTAGAAGTA-5′
m3C5′-FAM-AGTTCAATGAT-m3C-TTCAT-3′
3′-TCAAGTTACTAGAAGTA-5′
GATC5′-FAM-AGTTCAATGATCTTCAT-3′
3′-TCAAGTTACTAGAAGTA-5′
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Davletgildeeva, A.T.; Tyugashev, T.E.; Zhao, M.; Kuznetsov, N.A.; Ishchenko, A.A.; Saparbaev, M.; Kuznetsova, A.A. Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2. Cells 2023, 12, 1839. https://doi.org/10.3390/cells12141839

AMA Style

Davletgildeeva AT, Tyugashev TE, Zhao M, Kuznetsov NA, Ishchenko AA, Saparbaev M, Kuznetsova AA. Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2. Cells. 2023; 12(14):1839. https://doi.org/10.3390/cells12141839

Chicago/Turabian Style

Davletgildeeva, Anastasiia T., Timofey E. Tyugashev, Mingxing Zhao, Nikita A. Kuznetsov, Alexander A. Ishchenko, Murat Saparbaev, and Aleksandra A. Kuznetsova. 2023. "Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2" Cells 12, no. 14: 1839. https://doi.org/10.3390/cells12141839

APA Style

Davletgildeeva, A. T., Tyugashev, T. E., Zhao, M., Kuznetsov, N. A., Ishchenko, A. A., Saparbaev, M., & Kuznetsova, A. A. (2023). Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2. Cells, 12(14), 1839. https://doi.org/10.3390/cells12141839

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop