Abstract

Heme is an essential prosthetic group for many proteins involved in fundamental biological processes in all three domains of life. In Eukaryota and Bacteria heme is formed via a conserved and well-studied biosynthetic pathway. Surprisingly, in Archaea heme biosynthesis proceeds via an alternative route which is poorly understood. In order to formulate a working hypothesis for this novel pathway, we searched 59 completely sequenced archaeal genomes for the presence of gene clusters consisting of established heme biosynthetic genes and colocalized conserved candidate genes. Within the majority of archaeal genomes it was possible to identify such heme biosynthesis gene clusters. From this analysis we have been able to identify several novel heme biosynthesis genes that are restricted to archaea. Intriguingly, several of the encoded proteins display similarity to enzymes involved in heme d1 biosynthesis. To initiate an experimental verification of our proposals two Methanosarcina barkeri proteins predicted to catalyze the initial steps of archaeal heme biosynthesis were recombinantly produced, purified, and their predicted enzymatic functions verified.

1. Introduction

Heme, a modified tetrapyrrole, acts as an essential prosthetic group in many enzymes, sensory, and regulatory proteins. Hemes are also essential components of electron transport chains driving aerobic and anaerobic respiration and photosynthesis in almost all living organisms. Consequently, heme-containing proteins are found in all three domains of life, the Eukaryota, the Bacteria, and the Archaea. The biosynthesis of this important and ubiquitously distributed molecule has been intensively studied in eukaryotic and bacterial organisms, but little is known about heme biosynthesis in archaea. It is now well established for bacteria and eukarya that heme biosynthesis proceeds along a conserved pathway with highly related enzymes and identical biosynthetic intermediates (Figure 1(a)) [1]. Heme synthesis represents just one component of a larger, branched tetrapyrrole biosynthesis pathway, which is also responsible for the synthesis of chlorophylls, bacteriochlorophylls, cobalamin, siroheme, heme d1 and coenzyme F430 (Figure 1(b)) [2].

The common precursor for the formation of heme and all other tetrapyrroles is 5-aminolevulinic acid (ALA). Depending on the organism this molecule is either synthesized through the condensation of glycine and succinyl-CoA (Shemin pathway) by ALA synthase (HemAA) or in a two-step enzymatic process from glutamyl-tRNA via the intermediate glutamate-1-semialdehyde (GSA) by glutamyl-tRNA reductase (HemAB) and GSA-2,1-aminomutase (HemL) (C5-pathway) [3, 4]. Eight molecules of ALA are then converted into uroporphyrinogen III (UROGEN), the first cyclic tetrapyrrole of the pathway, in three consecutive enzymatic steps. First, two ALA molecules are condensed by porphobilinogen synthase (HemB) to the pyrrole derivative porphobilinogen (PBG) [5]. In the next step, four PBG molecules are oligomerized to the linear tetrapyrrole pre-uroporphyrinogen by PBG deaminase (HemC) and finally uroporphyrinogen III (UROGEN) is formed by cyclization of the pre-uroporphyrinogen by UROGEN synthase (HemD) [6]. The intermediate uroporphyrinogen III represents the last common precursor for all tetrapyrroles and is therefore an important branchpoint of the pathway. One of the diverting biosynthetic routes leads to the formation of hemes and (bacterio)chlorophylls via the intermediate coproporphyrinogen III (COPROGEN) and the other represents the first step of cobalamin, siroheme, heme d1, and coenzyme F430 biosyntheses via the common intermediate precorrin-2. COPROGEN and precorrin-2 are formed from UROGEN by the key branchpoint enzymes uroporphyrinogen III decarboxylase (HemE) and S-adenosyl-L-methionine-dependent uroporphyrinogen III methyltransferase (SUMT), respectively. Eukaryotic and bacterial heme biosynthesis further proceeds via the conversion of COPROGEN into protoporphyrinogen IX (PROTOGEN) by coproporphyrinogen III oxidase (HemF) or dehydrogenase (HemN) and the subsequent oxidation by protoporphyrinogen IX oxidase (HemY, HemG) to protoporphyrin IX (PROTO) [7]. Finally, the insertion of ferrous iron into PROTO by ferrochelatase (HemH) yields the end product heme [8]. All heme biosynthetic enzymes have been purified from many different eukaryotic and bacterial organisms and biochemically characterized [1]. The corresponding genes (hemA, L, B, C, D, E, F, N, Y, G, H) have all been cloned and sequenced [9]. In Figure 1(a) the commonly used gene designations for all bacterial heme biosynthesis genes are given together with their corresponding enzyme names.

In two independent bioinformatics studies in 2002 and 2008 the distribution of heme biosynthetic genes in prokaryotic organisms was investigated by analysis of the currently available sequenced microbial genomes [10, 11]. It was found that almost all bacteria which synthesize heme de novo possess the complete set of hem genes (i.e., hemAA or hemAB, hemL, hemB, hemC, hemD, hemE, hemF and/or hemN, hemY or hemG, hemH). In contrast, some bacteria (e.g., Clostridia and Desulfovibrio species) and almost all archaea were found to possess only the genes encoding the enzymes required for UROGEN formation (hemAB, L, B, C, D) and lacked the genes encoding the enzymes necessary for the conversion of UROGEN into heme [10, 11]. Such a finding can be explained by (i) the possibility that these organisms have no need for heme and require the initial genes for cobalamin, siroheme, coenzyme F430, or heme d1 formation, (ii) the possibility that they take up heme from the environment, or (iii) the existence of an alternative, yet unknown, heme biosynthesis pathway. For many free-living bacteria and archaea heme uptake is not very likely since heme is simply not available in their environment. However, a few examples of such cases exist in archaea and bacteria [12, 13]. Moreover, several pathogenic bacteria take up heme from their host and use it as an iron source. However, these bacteria often possess an intact heme biosynthetic apparatus [11].

It has been known for some time that Desulfovibrio species and many archaea contain cytochromes and other heme-containing proteins [1422], and therefore they must be able to synthesize their own heme. Indeed, for the sulfate-reducing bacterium Desulfovibrio vulgaris and the methanogenic archaeon Methanosarcina barkeri it was shown experimentally that an alternative heme biosynthesis pathway must exist. In these cases in vivo labeling studies demonstrated that their hemes contain methyl groups on rings A and B that are derived from methionine (via S-adenosyl-L-methionine) and not from ALA as is the case for hemes synthesized via the classical pathway [23, 24]. Further, in D. vulgaris sirohydrochlorin (the oxidized form of precorrin-2), 12,18-didecarboxysirohydrochlorin, coproporphyrin III, and PROTO were isolated as potential heme biosynthesis intermediates [25]. Thus, the alternative heme biosynthesis pathway seems to branch off the classical pathway at the stage of UROGEN. In the first step of the alternative route UROGEN is methylated at rings A and B by a SUMT-like enzyme to yield precorrin-2 (Figure 1(b)). This SUMT-dependent methylation of UROGEN is also required for the biosyntheses of cobalamin, siroheme, heme d1, and coenzyme F430 (Figure 1(b)). Recently, in D. vulgaris a bifunctional enzyme carrying both UROGEN synthase and SUMT activities and a precorrin-2 dehydrogenase (PC2-DH) catalyzing precorrin-2 oxidation to sirohydrochlorin were biochemically characterized [26]. Both enzymes are probably involved in the alternative heme biosynthesis pathway in this organism.

Besides the in vivo labeling study in M. barkeri the alternative heme biosynthesis pathway has not been investigated in archaea, so far. In the last few years the number of completely sequenced archaeal genomes has greatly increased and therefore we decided to start our investigation of archaeal heme biosynthesis with the search for potential heme biosynthesis gene clusters within these genomes. We found that many archaea indeed contain gene clusters consisting of the known early heme biosynthesis genes (hemAB, hemL, hemB, hemC, hemD) and of “nir-like” genes which encode proteins homologous to proteins involved in heme d1 biosynthesis in denitrifying bacteria. Moreover, very often the genes encoding a putative SUMT and a potential PC2-DH were found localized in these archaeal heme biosynthesis gene clusters. Here, the predicted SUMT and PC2-DH from M. barkeri were recombinantly produced, purified and shown in vitro to carry SUMT and PC2-DH activity, respectively.

2. Materials and Methods

2.1. Chemicals

All chemicals, reagents, and antibiotics were obtained from Sigma-Aldrich (Taufkirchen, Germany) or Merck (Darmstadt, Germany). DNA polymerase, restriction endonucleases, and PCR requisites were purchased from New England Biolabs (Frankfurt a.M., Germany). Oligonucleotide primers were obtained from metabion international AG (Martinsried, Germany). PCR purification and gel extraction Kits were purchased from Qiagen GmbH (Hilden, Germany). Ni Sepharose 6 Fast Flow was purchased from GE Healthcare (München, Germany). Uroporphyrin III was obtained from Frontier Scientific Europe (Carnforth, UK).

2.2. Construction of Vectors for Recombinant Protein Production

The gene mba_A1461 encoding a potential PC2-DH from M. barkeri was PCR amplified using the primers 01CysGN_Mba_BamHI_fw (GAA GGG ATC CGA TGA CCA AAA CCA ATA ATT TTC) and 02CysGN_Mba_NotI_rev (GAA CGC GGC CGC TTA ACG GTT GCT GTT CAC) containing BamHI and NotI restriction sites (underlined) and cloned into appropriately cut pET-Duet-1 (Novagen, Darmstadt, Germany) to generate pET-Duet_mba_A1461. The plasmid pMA_mba_A1461 (GeneART, Regensburg, Germany), which contains a synthetic copy of the mba_A1461 gene, codon-optimized for expression in Escherichia coli, was used as the DNA template for PCR.

The gene mba_A1791 encoding a putative SUMT from M. barkeri was PCR amplified using the primers MbarcobA-ATG (CAC ATA TGT CAG GAA ATT ACG GAA AAG) and MbarcobA-Stop (AGG ATC CAA AAC TAG TTA AAA GTC AAC TCC TGT CCG) containing NdeI and SpeI-BamHI restriction sites (underlined) from genomic M. barkeri DNA. The resulting PCR fragments and the vector pET14b (Novagen) were subsequently digested with NdeI and BamHI and ligated to generate pET14b_mba_A1791.

2.3. Bacterial Strains and Growth Conditions

E. coli DH10B was used as the host for cloning. For production of recombinant proteins the E. coli strains BL21 (DE3) and BL21 Star (DE3) pLysS were used, respectively. The expression vector pET14b_mba_A1791 was transformed into E. coli BL21 Star (DE3) pLysS. The vector pET-Duet_mba_A1461 was transformed into E. coli BL21 (DE3). For recombinant protein production the E. coli strains carrying the corresponding vectors were grown at 37°C in LB-medium containing appropriate antibiotics. Protein production was induced by adding 50 μM isopropyl isopropyl-β-D-thiogalactopyranosid (IPTG) to the cultures at an optical density at 578 nm of 0.6. The E. coli BL21 (DE3) strain containing pET-Duet_mba_A1461 was further cultivated at 37°C for 4 h. The E. coli BL21 Star (DE3) pLysS strain containing pET14b_mba_A1791 was further cultivated at 17°C for 18 h. The cells were harvested by centrifugation and stored at −20°C.

2.4. Purification of Enzymes and Tetrapyrroles

The affinity chromatographic purification of the recombinant enzymes was performed as described before [27] with minor changes. Briefly, for resuspension of the E. coli cells, harbouring the produced recombinant protein, buffer A (50 mM Tris/HCl (pH 7.5), 300 mM NaCl, 10% (w/v) glycerol) was used containing 1 mM phenylmethanesulfonyl fluoride. The cells were disrupted using a French press (1000 p.s.i.) and the soluble protein fraction was obtained by ultracentrifugation (60 min, 175000 × g, 4°C). The supernatant was applied to 1 mL of Ni Sepharose 6 Fast Flow (GE Healthcare). The flow-through containing the tetrapyrroles accumulated during in vivo protein production was applied to a 1 mL silica gel 100 C18-reversed phase column (Sigma-Aldrich) and the tetrapyrroles were extracted as described before [27]. The Ni-resin with bound proteins was washed extensively with buffer A. After a preelution step with buffer A containing 20 mM imidazole the recombinant protein was eluted with buffer A containing 300 mM imidazole. Immediately after elution a buffer exchange was performed in an anaerobic chamber (Coy Laboratories, Grass Lake, MI, USA) by passing the protein solution through a NAP-25 column (GE Healthcare) that had been equilibrated with degassed buffer A containing 5 mM dithiothreitol. The protein was stored at −20°C until required.

2.5. Determination of Protein Concentration

The Bradford Reagent (Sigma-Aldrich) was used to determine protein concentrations, according to the manufacturer’s instructions, using BSA as a standard.

2.6. Molecular Mass Determination

In order to determine the oligomeric state of proteins a gel permeation chromatography was performed using a Superdex 200 10/30 GL column with an ÄKTA Purifier system (GE Healthcare). The column was equilibrated with buffer A containing 5 mM dithiothreitol and calibrated using four standard proteins: cytochrome c, conalbumin, alcohol dehydrogenase, β-amylase (Kit Gel Filtration molecular weight markers (Sigma-Aldrich)). Protein samples (1 mg mL-1) were applied to the column and the elution of proteins was monitored by determination of the absorption of the eluate at 280 nm as described before [27].

2.7. In Vitro Enzyme Activity Assays

In vitro enzyme activities of the recombinantly produced and purified Mba_A1791 and Mba_A1461 proteins were measured using a coupled enzyme assay as described before [27]. The assay was performed in an anaerobic chamber (Coy Laboratories) under strictly anaerobic conditions (O2 = 0 ppm). The substrate uroporphyrinogen III was generated enzymatically from 1 mM ALA using purified HemB (0.14 μM) from Pseudomonas aeruginosa, HemC (0.15 μM), and HemD (0.17 μM) both from Bacillus megaterium in a final volume of 1 mL of degassed buffer B containing 50 mM Tris/HCl (pH 8.0), 100 mM KCl, 5 mM MgCl2, and 50 mM NaCl. In order to investigate the activity of Mba_A1791 the enzyme was added to a final concentration of 1.5 μM, and SAM as methyl donor was added to a final concentration of 200 μM. In order to determine the activity of Mba_1461, precorrin-2 was generated using P. aeruginosa NirE as SUMT [27] at a concentration of 1.5 μM. The Mba_1461 was added to a final concentration of 1.5 μM with 100 μM NAD+. The reaction mixtures were incubated overnight at 37°C in the dark. UV-visible spectra of the assay mixtures were recorded on a V-650 spectrophotometer (Jasco, Gross-Umstadt, Germany).

2.8. Bioinformatics Analysis

For the analysis and comparison of archaeal genomes the “Microbial Genome Database for Comparative Analysis” (http://mbgd.genome.ad.jp/) was used [2830]. This database contains a total of 68 completely sequenced archaeal genomes. Of these 68 genomes we initially chose one for each species analyzed, that is, different strains within one species were not included, which left us with 59 genomes. The archaeal species whose genomes were chosen are listed in Table 1. We also included the genomes of E. coli, P. aeruginosa, and D. vulgaris as positive and negative controls for our search. First, the database was searched for the known early heme biosynthesis genes of E. coli (hemAB, hemL, hemB, hemC, hemD) and then the database was used to find “orthologous clusters” in the archaeal genomes. With the “orthologous cluster” tool all homologous hem genes in the chosen genomes were displayed and a multiple genome map comparison could be viewed. For the clustering parameters we chose the default values of the database. Using the “multiple genome map comparison” tool we identified gene clusters with similar gene organization in the neighborhood of the known early heme biosynthesis genes in the archaeal genomes.

3. Results and Discussion

3.1. The Late Heme Biosynthesis Genes Are Missing in Archaeal Genomes

In order to identify potential heme biosynthesis gene clusters in the 59 archaeal genomes analyzed we first checked for the presence and genomic localization of the early heme biosynthesis genes hemAB, hemL, hemB, hemC, and hemD. Next, we inspected the genes located in their direct neighborhood using the MBGD database. Out of the 59 archaeal genomes included in this study we found 12 genomes which do not contain any obvious hem gene (Table 1). These organisms apparently do not synthesize tetrapyrroles de novo unless via a completely novel pathway. Alternatively, these members of the Archaea do not need heme and other tetrapyrroles, respectively, or they are able to take up these compounds from their environment as described previously [12]. For example, it was reported recently that many archaeal species possess genes that encode putative homologs of the prokaryotic BtuFCD system for cobalamin uptake [31]. Here, in the 12 genomes that lack the hem genes we also found btuFCD homologs with the exception of the Korarchaeum cryptofilum and the Nanoarchaeum equitans genomes (not shown).

However, in the majority (47) of the studied genomes we found all five hem genes (hemAB, hemL, hemB, hemC, hemD) whose encoded protein products are known to be responsible for formation of the tetrapyrrole precursor UROGEN. One exception is the genome of Aeropyrum pernix which is missing a recognizable hemD gene. As already observed before [10, 11] we failed to detect the genes hemE, hemF/N, hemG/Y, hemH encoding the known late heme biosynthesis enzymes catalyzing the conversion of UROGEN into heme. Exceptions from this rule came from the analysis of the Picrophilus torridus, Thermoplasma acidophilum, and Thermoplasma volcanium genomes. The genomes of these three species contain hemE and hemH genes encoding UROGEN decarboxylase and ferrochelatase, respectively, as highlighted in earlier studies [10, 11]. However, no genes encoding recognizable COPROGEN oxidases/dehydrogenases (hemF/N) or PROTOGEN oxidases (hemG/Y) were found. Thus, most archaea possess the genetic potential for synthesizing UROGEN from glutamyl-tRNA via the intermediates GSA, ALA, PBG, and pre-uroporphyrinogen and appear to have genes for heme-containing proteins. Therefore, to make heme they must transform the UROGEN by a novel pathway that differs from the known heme biosynthesis route. This is consistent with the observation that the methanogenic archaeon M. barkeri synthesizes its heme via the intermediate precorrin-2 [24].

3.2. Archaeal Hem Genes Are Clustered with SUMT and PC2-DH Genes

Upon closer bioinformatical inspection of the chromosomal organization of the detected archaeal hem genes we found that they are often located within gene clusters comprising two or more hem genes (Table 1 and Figure 2). Interestingly, within these hem gene clusters we also detected genes potentially encoding a SUMT and a precorrin-2 dehydrogenase (PC2-DH). SUMT proteins catalyze the S-adenosyl-L-methionine-dependent methylation of UROGEN on rings A and B at positions 2 and 7 to give precorrin-2 (Figure 1(b)). The PC2-DH proteins in turn oxidize the precorrin-2 to sirohydrochlorin in a NAD+-dependent reaction. Unfortunately, the nomenclature for these two genes in the MBGD database is quite inconsistent as the SUMT encoding gene is sometimes named cobA, cysG-1, cysG-2, cysG, uroM, or hemX, and the gene encoding PC2-DH is referred to as sirC, hemX, cysG, or cysG1. In the following we will refer to the genes encoding the methyltransferase and the dehydrogenase simply as the SUMT and PC2-DH genes, respectively. In archaea a SUMT is probably required for the synthesis of all tetrapyrroles including heme, cobalamin, siroheme, and coenzyme F430 [24, 32]. The PC2-DH has been shown conclusively to be involved in siroheme and anaerobic cobalamin formation in bacteria [33, 34] and will probably fulfill this function also in archaea. Its involvement in heme and coenzyme F430 biosyntheses has not yet been demonstrated and requires further experimental evidence.

Although both enzymes are required for the formation of all these different tetrapyrroles in archaea, it is interesting to note that almost all archaeal species possess only one SUMT and one PC2-DH gene. The only exception from this rule is Archaeoglobus fulgidus which possesses two SUMT genes. As outlined above, the SUMT and PC2-DH genes are often clustered on the genomes with the early hem genes. This clustering of genes encoding the enzymes responsible for the transformation of ALA into precorrin-2 or sirohydrochlorin provides the organisms with the possibility of coordinated gene expression and production of enzymes catalyzing consecutive biosynthetic steps. However, such heme biosynthesis gene clusters were not found in all archaeal genomes. For some of the investigated species, like Ignicoccus hospitalis and Caldivirga maquilingensis, the heme biosynthesis genes were found scattered randomly throughout the genome (Table 1).

3.3. Potential Involvement of nir-Like Genes in Archaeal Heme Biosynthesis

Interestingly, in 32 of the archaeal genomes that contain the early hem genes we also found so-called nir genes (nirD, nirH, nirJ) co-localized in large gene clusters with hemAB, hemL, hemB, hemC, hemD, SUMT, and PC2-DH genes (Table 1 and Figure 2). It was previously reported that D. vulgaris and some methanogenic archaea harbor these nir genes on their genomes. It was speculated that the nir genes might be involved in the alternative heme biosynthesis pathway in these organisms [26]. Here, we show not only that the methanogenic archaea contain nir genes, but also that the majority of archaea that synthesize heme de novo require these genes (Table 1). These nir genes encode proteins that are homologous to proteins involved in heme d1 biosynthesis in denitrifying bacteria such as P. aeruginosa. The dioxoisobacteriochlorin heme d1 serves as an essential prosthetic group in the cytochrome cd1 nitrite reductase which catalyzes the second step of denitrification [36]. However, based on amino acid sequence homology searches, only Pyrobaculum aerophilum, Pyrobaculum arsenaticum, and Pyrobaculum calidifontis possess a potential cytochrome cd1 nitrite reductase. All other archaeal genomes analyzed in this study do not. Consequently, the majority of archaeal nir genes are not involved in heme d1 biosynthesis. Rather, they are likely to be involved in heme biosynthesis. Therefore, we renamed these nir-like genes in the Archaea ahb(archaeal heme biosynthesis)-nir genes.

3.4. Structures of Potential Heme Biosynthesis Gene Clusters in Archaea

As mentioned above, the ahb-nir genes are often clustered with the hem, SUMT, and PC2-DH genes on the archaeal genomes. The most complete gene clusters, comprising ten out of the eleven potential heme biosynthesis genes, were found in the genomes of P. aerophilum and P. arsenaticum (Figure 2). In P. aerophilum these genes form one large, uninterrupted gene cluster. Another striking clustering of the potential heme biosynthesis genes was observed in the Methanosarcinales (Figure 2). For example, in the genomes of M. acetivorans and M. barkeri the genes hemAB, hemL, hemB, hemC, ahb-nirD, ahb-nirH, ahb-nirJ1, and PC2-DH are organized as one continuous gene cluster, while the genes hemD, ahb-nirJ2, and SUMT are localized together in a second gene cluster. In Halobacterium sp. NRC-1 three heme biosynthesis-related gene clusters were found. The first cluster comprises the genes hemL, hemB, hemC, hemD, and SUMT, the second consists of hemAB, PC2-DH, ahb-nirD, and ahb-nirH and the third contains ahb-nirJ1, and ahb-nirJ2 (Figure 2). In the other archaeal species that possess ahb-nir genes the clustering with the hem, SUMT, and PC2-DH genes is less distinct, but there is still often a colocalization of one or two ahb-nir genes with one or several hem genes (Table 1).

3.5. Proposed Function of the ahb-nir Genes during Heme Biosynthesis in the Archaea

As already mentioned above the ahb-nir genes encode proteins that are similar to proteins involved in heme d1 biosynthesis. We compared the amino acid sequences of the Ahb-Nir proteins from M. barkeri with the Nir proteins involved in heme d1 biosynthesis in P. aeruginosa. We found the following sequence identities: M. barkeri Ahb-NirD and P. aeruginosa NirD: 36.3%; Ahb-NirH and NirH: 40.1%; Ahb-NirJ1 and NirJ: 29.5%; Ahb-NirJ2 and NirJ: 38.8%. Thus, the Ahb-NirJ2 is more similar to the heme d1 biosynthesis protein NirJ than is the Ahb-NirJ1. Additionally, we found an amino acid sequence identity of 31.8% between the two Ahb-NirJ proteins.

Although the precise functions of the Nir proteins involved in heme d1 biosynthesis has not yet been established, several reasonable proposals were made [26, 3640]. First of all, it is known that heme d1 is biosynthesized from precorrin-2 [27, 35]. In order to obtain heme d1 from this precursor the following modifications have to take place: (a) decarboxylation of the acetate groups on rings C and D, (b) removal of the propionate side chains on rings A and B and replacement by oxo groups, (c) formation of an acrylate side chain on ring D, (d) oxidation of the tetrapyrrole macrocycle, and (e) iron insertion. The order of these reactions is not known. However, it was proposed that the oxidation reaction (b) might be catalyzed by the NirJ protein during heme d1 formation [37]. NirJ belongs to the so-called Radical SAM enzyme family whose members are known to catalyze chemically challenging reactions through radical-based mechanisms [41, 42]. It was also speculated that the NirD, NirL, NirG, and NirH proteins might be responsible for the decarboxylation reaction (a) [38].

In order to form heme from precorrin-2 the decarboxylation of the acetate groups on rings C and D as in reaction (a) is required and thus might be catalyzed by Ahb-NirD and Ahb-NirH. Further, the acetate side chains on rings A and B have to be removed, probably in a reaction which resembles the mechanism of reaction (b). Consequently, the Radical SAM enzyme Ahb-NirJ2 which shares 38.8% sequence identity with P. aeruginosa NirJ is a good candidate for catalyzing this reaction. This function was previously proposed for one of the NirJ-like proteins from D. vulgaris [26]. Another reaction which is required for heme formation from precorrin-2 that has, however, no equivalent in heme d1 biosynthesis is the oxidative decarboxylation of the propionate side chains on rings A and B to the corresponding vinyl groups. This reaction also takes place during the classical heme biosynthesis route in most bacteria and the Eukaryota. In bacteria it is catalyzed by either HemF or HemN (see Figure 1(a)). HemN also belongs to the Radical SAM enzyme family [43]. Thus, Ahb-NirJ1 (Radical SAM family member) might catalyze the formation of the required vinyl groups. In summary, we propose that the Ahb-Nir proteins catalyze some of the late reaction steps during archaeal heme biosynthesis from precorrin-2 (Figure 3).

3.6. Distribution of the ahb-nir Genes over the Archaeal Genomes

In accordance with the proposed function of the ahb-nir gene products during the late steps of archaeal heme biosynthesis we failed to detect any of the ahb-nir genes in those archaeal genomes without any hem genes (Table 1). However, the presence of the hem genes in an archaeal genome does not necessarily mean that the ahb-nir genes are also present. As mentioned above, out of 47 archaeal genomes containing all five early hem genes only 32 also contain the ahb-nir genes. The 15 archaeal species which possess the hem, but no ahb-nir genes probably synthesize their UROGEN solely as precursor for siroheme [44, 45], cobalamin [31] and, in the case of methanogens, for coenzyme F430 [46]. However, most likely they do not form heme. Accordingly, almost all of these 15 species also possess both a SUMT and a PC2-DH gene which are required for siroheme, cobalamin, and coenzyme F430 biosynthesis. One exception is P. torridus which does not contain a recognizable PC2-DH gene. Moreover, T. acidophilum and T. volcanium do not possess a SUMT gene. Interestingly, these three species are the only representatives of the Archaea for which hemE and hemH genes were found (see above). For P. torridus and T. acidophilum heme-containing proteins were biochemically characterized [18, 20]. However, considering the observation that their genomes lack recognizable ahb-nir genes and some of the late hem genes their route of heme biosynthesis remains currently unclear.

Within the group of the 32 archaeal species that contain ahb-nir genes several subgroups can be recognized. First of all, there are those species for which a complete set of ahb-nir genes (ahb-nirD, ahb-nirH, ahb-nirJ1, ahb-nirJ2) was found. Out of the 32 genomes containing ahb-nir genes 27 contain all four of them. Among the 5 genomes in which not all ahb-nir genes are present, those of Halorhabdus utahensis and Methanopyrus kandleri are missing ahb-nirD and ahb-nirH. The genomes of Sulfolobus acidocaldarius, Haloquadratum walsbyi, and Nirosopumilus maritimus do not contain ahb-nirJ1 and ahb-nirJ2. Thus, for these five species it is questionable whether they synthesize heme themselves, if at all required. For example, for M. kandleri and N. maritimus no indications were found in the literature or databases that they possess heme-containing proteins. Further, two subgroups of ahb-nir containing archaea can be distinguished depending on whether they possess two distinct ahb-nirD and ahb-nirH genes or whether they contain an ahb-nirDH gene fusion. In fact, almost all archaea possess the fused genes except for those methanogenic archaea which contain the ahb-nir genes (Table 1). However, in these methanogens the ahb-nirD and ahb-nirH genes are always located aside to each other on the genome with the only exception being Methanosaeta thermophila. Likewise, the two ahb-nirJ genes are also often (15 out of 29) co-localized on the genomes, either as direct neighbors or in close proximity to each other, indicating gene duplication as the origin of the two copies.

Our bioinformatics investigation of 59 archaeal genomes in combination with the experimental evidence that two methyl groups of archaeal heme are derived from S-adenosyl-L-methionine strongly suggests that heme biosynthesis in the Archaea follows a novel, yet mostly unknown route. It starts with the methylation of UROGEN to precorrin-2 catalyzed by SUMT, followed by the oxidation of precorrin-2 to sirohydrochlorin by PC2-DH and further transformations (decarboxylation of acetate groups, removal of acetate groups, oxidative decarboxylation of propionate to vinyl groups, and insertion of iron) of the macrocycle side chains which are most likely performed by the Ahb-Nir proteins (Figure 3). Clearly, these proposals need to be tested experimentally. Thus, we decided to first verify the predicted functions of the M. barkeri proteins Mba_1791 and Mba_1461 as SUMT and PC2-DH, respectively.

3.7. Production and Purification of Recombinant Mba_1791 and Mba_1461

The M. barkeri proteins Mba_1791 and Mba_1461 were recombinantly produced as N-terminal His-tagged fusion proteins in E. coli. In both cases the recombinant proteins were produced in a soluble form and in a high yield. We purified Mba_1791 and Mba_1461 to apparent homogeneity using a single affinity-chromatographic step on Ni Sepharose 6 Fast Flow (Figure 4(a)). The purified Mba_1791 exhibited a slight red-brown color. UV-visible absorption spectroscopy suggested the presence of a copurified tetrapyrrole, probably the reaction product of Mba_1791 (data not shown). For other SUMTs (e.g., P. aeruginosa NirE) the co-purification of their reaction product has been previously reported [27, 35, 47]. Therefore, the presence of a tetrapyrrole in the purified Mba_1791 was a first hint towards the function of this protein as a SUMT. In contrast, purified Mba_1461 appeared colorless.

The oligomeric state of Mba_1791 and Mba_1461 was determined by gel permeation chromatography. This experiment revealed a native relative molecular mass of 55,300 840 Da for Mba_1791 and 60,800 ± 7,300 Da for Mba_1461, respectively (Figure 4(a)). The calculated molecular masses based on the amino acid sequences of the proteins are 26,350 Da for Mba_1461 and 27,230 Da for Mba_1791. Thus, gel permeation chromatography suggests a dimeric structure for both proteins. Other SUMTs [48, 49] and PC2-DH [33, 50] are also thought to be dimeric proteins.

3.8. Mba_1791 Acts as a SUMT In Vivo

During production of Mba_1791 in E. coli a red compound accumulated and remained in the soluble protein fraction of the cell-free extract after disruption of the cells and ultracentrifugation. This compound was extracted using C18-reversed phase silica gel and analyzed by UV-visible absorption spectroscopy. The UV-visible absorption spectrum of the extracted compound exhibited an absorption maximum at 378 nm which strongly resembled the previously reported spectra of sirohydrochlorin (Figure 4(b)) [50]. Obviously, recombinantly produced M. barkeri Mba_1791 showed significant SUMT activity in the production host E. coli which led to the accumulation of sirohydrochlorin, the oxidized form of the SUMT reaction product precorrin-2. Such tetrapyrrole accumulation during recombinant SUMT production has been described before. Depending on the source of the enzyme the formation of either trimethylpyrrocorphin, which is a nonphysiological trimethylated reaction product, or sirohydrochlorin was reported [27, 35, 47, 5153]. Apparently, Mba_1791 belongs to the class of SUMT enzymes that accumulates sirohydrochlorin and does not catalyze the overmethylation of precorrin-2 to trimethylpyrrocorphin.

3.9. M. barkeri Mba_1791 Is a SUMT

In order to investigate the in vitro activity of Mba_1791 a coupled enzyme assay was performed. The substrate uroporphyrinogen III was produced enzymatically and incubated overnight with recombinant purified Mba_1791. The formation of reaction products was followed using UV-visible absorption spectroscopy (Figure 4(c)). The absorption spectrum of a control assay mixture containing only the uroporphyrinogen III producing enzymes showed no characteristic absorption features under anaerobic conditions. In contrast, the addition of purified Mba_1791 and SAM to the reaction mixture resulted in a yellow colored solution after overnight incubation. The corresponding absorption spectrum exhibited a broad absorption between 350–400 nm and 400–500 nm which is characteristic for precorrin-2 [50, 54]. Consequently, Mba_1791 is indeed a SUMT. Therefore, we will name the enzyme from now on Mba_SUMT. The specific activity of Mba_SUMT was determined with uroporphyrinogen III (produced by chemical reduction of uroporphyrin III) at a concentration of 17 μM, a SAM concentration of 200 μM and a Mba_SUMT concentration of 1.5 μM. Under these conditions we observed a specific activity of 616 nmol precorrin-2 × h−1 × mg−1 Mba_SUMT corresponding to a turnover of 17 h−1. This specific activity is in the same range as the activities observed for other SUMTs [27, 5557].

3.10. M. barkeri Mba_1461 is a PC2-Dehydrogenase

In order to verify the postulated PC2-DH activity for Mba_1461 in an in vitro assay the enzymatically produced uroporphyrinogen III was converted to the PC2-DH substrate precorrin-2 by addition of the SUMT NirE from P. aeruginosa. Additionally, purified M. barkeri Mba_1461 and NAD+ were added to the reaction mixture. After overnight incubation a UV-visible absorption spectrum of the now purple reaction mixture was measured. The spectrum (Figure 4(c)) corresponds to a typical absorption spectrum of sirohydrochlorin with an absorption maximum at 378 nm [50]. Thus, Mba_1461 exhibited strong PC2-DH activity in vitro and can be safely assigned as Mba_PC2-DH.

We also tested the activities of the M. barkeri enzymes Mba_SUMT and Mba_PC2-DH in a coupled assay. Both were added to the reaction mixture containing all enzymes for uroporphyrinogen III generation. After overnight incubation the UV-visible absorption spectrum of this assay mixture was similar to the spectrum of the reaction mixture containing P. aeruginosa NirE (SUMT) and Mba_PC2-DH showing again the formation of sirohydrochlorin (Figure 4(c)). When SAM and/or NAD+ were omitted from this activity assay no formation of sirohydrochlorin was observed (data not shown).

4. Conclusion

In this study we have identified gene clusters in many archaeal genomes that are likely required for the biosynthesis of heme via a novel pathway. These gene clusters consist of (i) the known hem genes (hemAB, hemL, hemB, hemC, hemD) necessary for the formation of the heme precursor UROGEN, (ii) the SUMT and PC2-DH genes required for the synthesis of the intermediates precorrin-2 and sirohydrochlorin, and (iii) the ahb-nir genes whose protein products are probably responsible for the conversion of sirohydrochlorin into heme (Figure 3). We propose that the detected ahb-nir genes are involved in archaeal heme biosynthesis and not in heme d1 biosynthesis since almost all archaea do not possess a cytochrome cd1 nitrite reductase. However, we failed to detect any obvious genes encoding potential ferrochelatases within the archaeal heme biosynthesis gene clusters. The Archaea do not usually possess a hemH gene encoding the bacterial-type ferrochelatase but do contain several copies of genes encoding putative cobalt- and/or magnesium chelatases, although they do not synthesize chlorophylls. These gene products might be involved in the archaeal heme biosynthesis. To confirm our bioinformatics findings and predictions further experimental verification will be required in order to determine the exact function of the Ahb-Nir proteins in the archaeal heme biosynthesis process.

Acknowledgments

The authors would like to thank Professor R. Thauer (MPI Marburg) for the gift of genomic DNA from M. barkeri. They also thank Dr. J. Moser (TU Braunschweig) for helpful discussions. This work was financially supported by grants from the Deutsche Forschungsgemeinschaft, the Fonds der Chemischen Industrie to G. Layer and the BBSRC to M. J. Warren.