Enzymatic Deamination of the Epigenetic Base N-6-Methyladenine

Cog1001

A sequence similarity network for cog1001 is presented in Figure 1 at an E-value cutoff of 10⁻⁸⁰.^3,4 Adenine deaminase (ADE) is the predominant enzyme in cog1001 and is represented in Figure 1 by subgroups 1 and 4. Subgroups 3 and 5 consist of enzymes that are currently annotated as isoaspartyl dipeptidases and enamidases, respectively.^5,6 The functional annotations for most members of subgroup 2 in Figure 1 are currently listed as unknown in the NCBI database but some of them have been labeled as adenine deaminases without apparent experimental validation. Proteins in subgroup 2 have less than 35% sequence identity to the prototypical adenine deaminases found in subgroups 1 and 4.

There are sequences available for ~32 proteins from subgroup 2. These enzymes are found in various species of Bacillus. At least a dozen of these microorganisms also have an authentic adenine deaminase from subgroup 1 in their genome. Most of the genes in subgroup 2 from cog1001 lie adjacent to the pur operon for the biosynthesis of purine nucleotides. Moreover, in Bacillus halodurans, the genes for ADE from subgroup 1 (Bh0640) and Bh0637 from subgroup 2 are nearly adjacent to one another. The co-localization of these two genes suggests that these enzymes may use the same or similar substrates and that the products of their reactions are similar. Based on these observations, we have thus predicted that Bh0637 is likely to have an adenine-like substrate and that the reaction product resembles hypoxanthine.

There are a limited number of purine moieties in Nature that closely resemble adenine. A prime substrate candidate for enzymes in subgroup 2 is N-6-methyladenine because this compound is found in the genomes of all bacteria.⁷ This compound is often considered as the 6^th DNA base with 5-methylcytosine being the 5^th base. The occurrence of N-6-methyladenine in most bacterial genomes varies from 1.4 – 2.0% of the total adenine content.⁸ This nucleotide modification functions to protect host DNA against restriction endonucleases and is critical for cell viability.⁸ More recently, N-6-methyladenine has been implicated in DNA repair and replication, regulation of gene expression, control of transposition, and host-pathogen interactions.^7,8 In bacteria, DNA methylases and cell-cycle regulated methyltransferases are responsible for the methylation of adenine in DNA.^7,8 Previous studies with cell free extracts from various sources have implicated a cryptic enzyme activity for the conversion of N-6-methyladenine to hypoxanthine.^9,10 However, no enzyme is currently known to metabolize N-6-methyladenine. Other modified adenine moieties that could serve as potential substrates include 7-methyladenine and 2-hydroxyadenine (isoguanine).^11,12

Enzyme Purification

Bh0637 from B. halodurans C-125 was selected from enzymes clustered in subgroup 2 for substrate interrogation. The gene for Bh0637 was cloned into a pET30a(+) plasmid (Novagen) with NdeI and HindIII restriction sites, expressed in Escherichia coli and purified to homogeneity. Bsu06560 from Bacillus subtilis (gi|16077724) was also purified to homogeneity. These two proteins share 71% sequence identity. The gene for Bsu06560 was codon-optimized for expression in E. coli, synthesized, and then cloned into an expression vector that yields a protein with a non-cleavable C-terminal histidine tag. The sequence and expression/purification details for Bsu06560 are available from PepcDB (http://pepcdb.pdb.org/) under target ID “NYSGXRC-9208a”, clone 9208a4BCt13p1. The proteins were expressed with an iron chelator, 2,2′-dipyridyl, to prevent oxidation of active site sidechains. ICP-MS of the purified enzymes established the presence of ~2 equivalents of manganese per protein monomer. Apo-Bh0637 was prepared by dialysis against 10 mM o-phenanthroline at pH 7.5. Zn²⁺ and Fe²⁺ were added to reconstitute the enzymatic activity of Bh0637. Iron was added anaerobically to prevent air-oxidation prior to reconstitution.

Enzymatic Assays

Bh0637 and Bsu06560 were screened against a small highly focused library of approximately 20 adenine-like compounds (Supporting Information). All of the compounds in the compound library were monitored for changes in the UV spectrum. The compounds (100-150 μM) were dissolved in 20 mM HEPES, pH 7.5, and incubated with 1.0 μM Bh0637 or Bsu06560 for 2 hours. The only compounds showing any significant absorbance changes were adenine, N-6-methyladenine, 6-methoxypurine and 6-methylmercaptopurine.

The kinetic constants for the deamination of adenine were measured using a coupled assay for the formation of ammonia.¹³ The deamination of N-6-methyladenine was followed directly at 270 nm (Δε = 15,000 M⁻¹ cm⁻¹).⁸ The cleavage of methanol from 6-methoxypurine was monitored at 270 nm (Δε = 3760 M⁻¹ cm⁻¹) and the loss of methylthiol from 6-methylmercaptopurine was determined at 290 nm (Δε = 20,200 M⁻¹ cm⁻¹).^14,15 The formation of hypoxanthine was confirmed using UV absorbance spectroscopy. The production of methylamine from N-6-methyladenine was demonstrated by mass spectrometry and the release of methanol from 6-methoxypurine was verified with alcohol dehydrogenase. DTNB (5,5′-dithiobis(2-nitro-benzoic acid)) was utilized to monitor the formation of methylthiol from 6-methylmercaptopurine. The kinetic constants for the various substrates with the two enzymes are presented in Table 1. Since similar kinetic constants are observed for the Zn-, Mn-, and Fe-reconstituted forms of Bh0637 and Bsu06560, it is difficult to predict the metal used by either enzyme in vivo. The overall reaction for the deamination of N-6-methyladenine is represented in Scheme 1.

Structural Model for Bh0637

A comparative model of the three-dimensional structure for Bh0637 was constructed. The amino acid sequence of Bh0637 was first submitted to the PSI-BLAST server at NCBI to identify suitable template structures.¹⁶ The crystal structure of adenine deaminase (Atu4426, PDB code: 3NQB) was the only template that produced a significant alignment with an E-value (6 × 10⁻²⁰) better than the default threshold of 6 × 10⁻³ (Supporting Information). The sequence alignment between Bh0637 and Atu4426 was computed by profile-profile alignment as implemented in the “profile.scan” routine of MODELLER-9v8.¹⁷ A total of 500 comparative models were generated with the standard “automodel” class in MODELLER, and the model with the best DOPE score¹⁸ was selected for refinement of backbone and side-chain conformations of residues near the binuclear metal center (Ser-464 to Gly-469, Ser-552 to Pro-557) using the “loopmodel” class in MODELLER. Sidechains of residues in these two loops and Phe-92, Arg-253, Thr-467, Phe-550 and His-555 were subsquently refined using the “sidechain prediction” protocol in PLOP.¹⁹ This refinement protocol resulted in eight representative models of Bh0637.

Computational Docking to Bh0637

The high-energy intermediate (HEI) library of KEGG molecules was docked to all eight comparative models of Bh0637 by DOCK 3.5.54.^20-26 The model yielding the highest enrichment for N-6-methyladenine was selected to represent the most likely binding mode of this substrate in Bh0637. A comparison of the binding site in the ligand-free template structure of Atu4426 (PDB: 3NQB) and the binding site in the comparative model of Bh0637 is presented in Figure 2. In the Bh0637 model, Glu-233 and the six residues that coordinate the binuclear metal center (His-88, His-90, Glu-180, His-211, His-232 and Asp-284) overlap with the aligned residues from the structural template. Several residues that probably interact with N-6-methyladenine are not conserved between the target and the modeling template. Ser-97, His-122, Glu-123, Asp-290 and Gly-563 in Atu4426 are replaced by Phe-92, Leu-117, Met-118, Gly-285 and Phe-550 in Bh0637, respectively. The substrate binding site of Bh0637 is therefore more hydrophobic than that of Atu4426.

A model for the binding of the intermediate that would be formed by the attack of hydroxide on N-6-methyladenine in the active site of Bh0637 is presented in Figure 2. The tetrahedral intermediate was formed by attack of the bridging hydroxide on the C6 carbon from the re-face of N-6-methyladenine.²⁶ In this model, the N6 amino group hydrogen bonds with the side chain of Asp-284. The protonated N1 nitrogen hydrogen bonds with the sidechain of Glu-233. In addition, the N7 nitrogen coordinates the α-metal ion and the N9 nitrogen hydrogen bonds with the sidechain of Thr-467. The N3 nitrogen is adjacent to the Nε nitrogen from His-555 (3.9 Å), which was treated as protonated during the virtual screening calculations. The methyl group of the docked substrate is not predicted to be in close contact with any of the neighboring residues. However, if Gly-285 is replaced by a larger residue, such as the aspartate found at the equivalent position in Atu4426, the methyl group of the substrate would fall within 4 Å of the sidechain of the larger residue.

In addition to N-6-methyladenine, the top ranked compounds that could undergo nucleophilic attack by the hydroxide on the re-face were inspected visually. The docked pose of adenine is nearly identical to that of N-6-methyladenine but the relative rank of adenine (ranked 282) is worse than that of N-6-methyladenine (ranked 176). The prioritization of N-6-methyladenine over adenine appears to reflect the more favorable van der Waals interactions between the N-6-methyl group and hydrophobic residues in the substrate binding site, such as Phe-92 and Phe-550. Five other compounds that scored well in the computational docking were subsequently tested as substrates but none were found to be deaminated at measureable rates. These compounds included 5-amino-4-imidazolecarboxyamide, toxopyrimidine, 5-acetylamino-6-amino-3-methyluracil, 2-hydroxy-4-amino-6-(cyclopropylamino)-1,3,5-triazine, and deisopropylhydroxyatrazine (Supporting Information).

Physiological Relevance

For Bh0637 and Bsu06560 N-6-methyladenine is a 100-fold better substrate than the next best candidate compound based on the relative values of k_cat/K_m. N-6-methyladenine is a known metabolite and is widely prevalent in all bacteria. Bh0637 and Bsu06560 therefore can support the recycling of N-6-methyladenine to form hypoxanthine, which is presumably taken up by the purine salvage pathway to form guanine nucleotides from inosine.²⁷ These two enzymes can also use adenine to form hypoxanthine but not nearly at the same rate as N-6-methyladenine. We conclude that N-6-methyadenine is the likely physiological substrate for members of subgroup 2 in cog1001 and thus these enzymes should now be annotated as 6-methyladenine deaminases (6-MAD). Consistent with this view, the prototypical adenine deaminases in B. halodurans (Bh0640) and B. subtilus (Bsu14520) cannot use N-6-methyladenine as substrates.²⁷ The kinetic constants for the deamination of adenine by Bh0640 are 120 ± 4 s⁻¹, 0.30 ± 0.04 mM, and (4 ± 0.3) × 10⁵ M⁻¹ s⁻¹ for k_cat, K_m and k_cat/K_m respectively. The kinetic constants for the deamination of adenine by Bsu14520 are 140 ± 6 s⁻¹, 0.38 ± 0.04 mM, and (3.8 ± 0.3) × 10⁵ M⁻¹ s⁻¹ for k_cat, K_m and k_cat/K_m respectively.

Strategy for Functional Annotation

It has been estimated that nearly one third of the currently sequenced genes have an unknown, uncertain, or incorrect functional annotation.²⁸ We have demonstrated that Bh0437 and Bsu06560 catalyze the deamination of N-6-methyladenine. Conservation of residues suggest that the rest of the enzymes depicted in subgroup 2 of Figure 1 will catalyze the same enzymatic reaction (see Supporting Information for a list of annotated N-6-methyladenine deaminases in cog1001). The key observations that led to the successful experimental annotation of these proteins included: recognition of adenine deaminase as the closest functional homologue, localization of the gene for this enzyme near the operon for purine synthesis, and identification of residues within the substrate binding site that could accommodate larger substituents attached to the N6 nitrogen. These insights enabled experimental identification of the enzyme substrate using a relatively modest compound library. The construction of a comparative model and the docking of N-6-methyl adenine in the active site of Bh0637 elucidated the structural rationale for the preferential deamination of this compound relative to the poorer substrate, adenine.

Supplementary Material