# A Novel Pathway for the Biosynthesis of Heme inArchaea: Genome-Based Bioinformatic Predictions and Experimental Evidence

Authors: Sonja Storbeck; Sarah Rolfes; Evelyne Raux-Deery; Martin J. Warren; Dieter Jahn; Gunhild Layer
Journal: Archaea
 (2010)

---

## Abstract

Heme is an essential prosthetic group for many proteins involved in fundamental biological processes in all three domains of life. InEukaryota 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.


---

## Body

## 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, theEukaryota, 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].Tetrapyrrole biosynthesis pathways. (a) Heme biosynthesis in most bacteria and theEukaryota. The first common precursor in the classical heme biosynthesis pathway is ALA of which eight molecules are converted into UROGEN in three consecutive enzymatic steps. UROGEN is then further converted into heme through successive modifications of the macrocycle side chains and finally iron insertion. The enzymes involved in the classical heme biosynthesis are  glutamyl-tRNA reductase (GluTR), glutamate-1-semialdehyde-2,1-aminomutase (GSAM), 5-aminolevulinic acid synthase (ALAS), porphobilinogen synthase (PBGS), porphobilinogen deaminase (PBGD), uroporphyrinogen III synthase (UROS), uroporphyrinogen III decarboxylase (UROD), oxygen-dependent coproporphyrinogen III oxidase (CPO), coproporphyrinogen III dehydrogenase (CPDH), oxygen-dependent and oxygen-independent protoporphyrinogen IX oxidase (PPO), and ferrochelatase (FC). The corresponding bacterial gene names are denoted in brackets below the enzyme names. (b) Overview of the different branches of the tetrapyrrole biosynthesis pathway. The last common precursor for the formation of all tetrapyrroles is UROGEN. Hemes and (bacterio)chlorophylls share PROTO as their last common intermediate. Siroheme, cobalamin, coenzyme F430, and heme d1 are all biosynthesized via precorrin-2. In the Archaea and some bacteria an alternative heme biosynthesis pathway exists in which the heme is biosynthesized from precorrin-2 via as yet unknown intermediates.
(a)(b)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 thatDesulfovibrio species and many archaea contain cytochromes and other heme-containing proteins [14–22], 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 thein 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 frommetabion 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 genemba_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 genemba_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 [28–30]. 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.Table 1
Heme biosynthesis genes in Archaea.
OrganismhemABahemLhemBhemChemDSUMTPC2-DHnirDnirHnirJ1nirJ2hemebhemed1B12csirohemedF4301Aeropyrum  pernix K1APE_2296e_2299.1_2300.1_2298.1—_0236.1_1491.1_1497.1_1497.1_1655_1652+———12Desulfurococcus  kamchatkensis 1221n—————————————23Ignicoccus  hospitalis KIN4/IIGNI_0512_0683_0324_0684_0090_0089_0483_0839_0839_0397_0630———34Staphylothermus  marinus F1——————————————45Hyperthermus  butylicus DSM 5456HBUT_0206_0837_1386_0836_0834_0835_0510_0036_0036_0035_0037———56Metallosphaera  sedula DSM 5348MSED_0214_0216_0215_0217_0218_0608_0213_0013_0013_0512_0511—++—67Sulfolobus  acidocaldarius DSM 639SACI_0777_0779_0778_0780_0781_0914_0776_0711_0711——+—++—78Sulfolobus  islandicus M.14.25M1425_1955_1953_1954_1952_1951_0226_1956_1894_1894_1048_1047—+—89Sulfolobus  solfataricus P2SSO018001820181018301842435—0245024516311632+—++—91840183910Sulfolobus  tokodaii 7ST02120215021402170218056302110293029301270126+—++—1011Thermofilum  pendens Hrk 5———————————————1112Caldivirga  maquilingensis IC-167CMAQ_1730_1413_1740_1121_1901_1313_1731_0150_0150_1900_1347——+—1213Pyrobaculum  aerophilum IM2PAE06010594058305800589059005850581058105960579++++—1314Pyrobaculum  arsenaticum DSM 13514PARS_2261_2254_2247_2245_2250_2251_2248_2246_2246_2255_2244++—+—1415Pyrobaculum  calidifontis JCM 11548PCAL_2034 _1481_1717_1709_1707_1712_1713_1710_1708_1708_1716_1706++++—1516Pyrobaculum  islandicum DSM 4184PISL_0096_0114_0042_0050_0116_0115_0119_0044_0044_0113_0051——+—1617Thermoproteus  neutrophilus V24StaTNEU_1917_1900_0957_0968_1898_1899_1896_0962_0962_1901_0969—+—1718Archaeoglobus  fulgidus DSM 4304AF19751241197412420116124315921593M1594M24131125+—++—181594M19Haloarcula  marismortui ATCC 43049RRNAC17082628261030863088308717091711171134891363+—++—1920Halobacterium  sp. NRC-1VNG1774G2326G2322G2330G2332G2331G1775C1776G1776G1185G1184GM+—+—2021Halomicrobium  mukohataei DSM 12286HMUK_1612_0914_0925_1105_1103_1104_1613_1614_1614_1679_1984+—+—2122Haloquadratum  walsbyi DSM 16790:HBSQ001HQ3336A3447A3443A3450A3452A3451A3335A3334A3334A——+—++—2223Halorhabdus  utahensis DSM 12940HUTA_2825_1970_1362_0927_0928_1761_1755——_0031_0734+—+—2324Halorubrum  lacusprofundi ATCC 49239HLAC_2132_2622_0015_2252_2254_2253_2131_2130_2130_1215_2081+—++—2425Natronomonas  pharaonis DSM 2160NP4502A1246A0920A1326A1330A1328A4500A4498A4498A1546A1542A+—++—2526Methanobrevibacter  smithii ATCC 35061MSM_0967_1233_1476_0881_1504_1550_0968—————++2627Methanosphaera  stadtmanae DSM 3091MSP_1408_1180_0416_1332_1191_1192_1407—————++2728Methanothermobacter  thermautotrophicus delta HeMTH10122287448741661671013—————+++2829Methanocaldococcus  fervens AG86MEFER_0496_1387_0258_0770_1494_0849_0978—————+2930Methanococcus  jannaschii DSM 2661MJ0143060306430569099409650140—————+++3031Methanocaldococcus  vulcanius M7METVU_1458_0867_0542_0749_0380_1448_0492—————+3132Methanococcus  aeolicus Nankai-3MAEO_0052_1208_0916_1324_0075_0861_0053—————++3233Methanococcus  maripaludis S2MMP0088022412580872039409660089—————++3334Methanococcus  vannielii SBMEVAN_1100_1237_0571_0048_1381_0289_1101—————++3435Methanocorpusculum  labreanum ZMLAB_0523_0525_0524_0526_0073_0526_0522—————++35_035936Methanoculleus  marisnigri JR1MEMAR_0980_0982_0981_0983_0539_0984_0979_0878_0877_0879_0876—++3637Methanospirillum  hungatei JF-1MHUN_2562_2560_2561_2559_2267_2558_2563—————++3738Methanoregula  boonei 6A8MBOO_1238_1235_1236_1234_0514_1233_1239_0957_0956_0958_0955—++3839Methanosphaerula  palustris E1-9cMPAL_1728_1726_1727_1725_2749_1724_1729_2627_2628_2626_2629—+3940Methanosaeta  thermophila PTMTHE_0049_1126_0050_1125_1124_1124—_1134_0047_1135_1123+—++4041Methanococcoides  burtonii DSM 6242MBUR_1229_1227_1228_1226_1106_1105_1230_1232_1231_1233_1236+—+++4142Methanosarcina  acetivorans C2AMA05770581057805823034303305760574057505733035+—++4243Methanosarcina  barkeri fusaroMBAR_A1462_A1464_A1463_A1465_A1792_A1791_A1461_A1459_A1460_A1458_A1793+—++4344Methanosarcina  mazei Goe1MM_1741_1743_1742_1744_0308_0307_1740_1738_1739_1737_0309+—++4445Methanopyrus  kandleri AV19MK0200MKT0801980746155015481495——08960980—+++4546Pyrococcus  abyssi GE5——————————————4647Pyrococcus  furiosus DSM 3638——————————————4748Pyrococcus  horikoshii OT3——————————————4849Thermococcus  gammatolerans EJ3—————————————4950Thermococcus  kodakarensis KOD1——————————————5051Thermococcus  onnurineus NA1—————————————5152Thermococcus  sibiricus MM 739—————————————5253Picrophilus  torridus DSM 9790PTO091802481311024902501435—————+—+—5354Thermoplasma  acidophilum DSM 1728TA05360571095505720573—0652————+—+—5455Thermoplasma  volcanium GSS1TVN05900635110006340633—0924—————+—5556uncultured methanogenic archaeon RC-IRCIX911913912914916915909—————+5657Korarchaeum  cryptofilum OPF8—————————————5758Nanoarchaeum  equitans Kin4-M——————————————5859Nitrosopumilus  maritimus SCM1NMAR_0510_0490_0509_0491_0493_0492_0511_0512_0512———+—5960Escherichia  coli K-12 MG1655B1210015403693805380433683368————+——+—6061Pseudomonas  aeruginosa PAO1PA466639775243526052590510f26110515f0512f—0511f++++—6162Desulfovibrio  vulgaris HildenboroughDVU14613168085618900734073414630854316708550857+—++—62ahemAB gene encoding glutamyl-tRNA reductase in the listed organisms. In the case of two or more copies for one gene, these copies are only listed if none of them is located in the gene cluster, otherwise only the copy within the cluster is listed. This is valid for all listed genes. Genes colocated within the same gene cluster in a certain species are highlighted in the same font (bold, italics, bold italics).bOrganisms possessing heme-containing proteins according to the literature [14–22] are marked with +. For organisms with empty fields the presence of heme-containing proteins was not clearly obvious from literature and bioinformatics data.cOrganisms possessing cobalamin biosynthesis genes according to [35] are marked with +. Organisms that do not synthesize cobalamin are marked with — [35]. For organisms with empty fields the ability to synthesize cobalamin was not obvious from the literature.dOrganisms possessing siroheme-containing sulfite or nitrite reductases based on sequence homology are marked with +. For organisms with empty fields the presence of siroheme-containing sulfite and nitrite reductases was not clearly obvious from literature and bioinformatics data.eNumbers represent the database gene number. The full database gene numbers including the strain abbreviation is given for the hemAB genes. For all other genes the database gene numbers are given without the strain abbreviation. For genes that are located as direct or near neighbors in an individual strain the numbers are written in the same font.fP. aeruginosa  nirE, nirD, nirH, and nirJ genes involved in heme d1 biosynthesis.


## 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 frommetabion 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 genemba_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 genemba_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 [28–30]. 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.Table 1
Heme biosynthesis genes in Archaea.
OrganismhemABahemLhemBhemChemDSUMTPC2-DHnirDnirHnirJ1nirJ2hemebhemed1B12csirohemedF4301Aeropyrum  pernix K1APE_2296e_2299.1_2300.1_2298.1—_0236.1_1491.1_1497.1_1497.1_1655_1652+———12Desulfurococcus  kamchatkensis 1221n—————————————23Ignicoccus  hospitalis KIN4/IIGNI_0512_0683_0324_0684_0090_0089_0483_0839_0839_0397_0630———34Staphylothermus  marinus F1——————————————45Hyperthermus  butylicus DSM 5456HBUT_0206_0837_1386_0836_0834_0835_0510_0036_0036_0035_0037———56Metallosphaera  sedula DSM 5348MSED_0214_0216_0215_0217_0218_0608_0213_0013_0013_0512_0511—++—67Sulfolobus  acidocaldarius DSM 639SACI_0777_0779_0778_0780_0781_0914_0776_0711_0711——+—++—78Sulfolobus  islandicus M.14.25M1425_1955_1953_1954_1952_1951_0226_1956_1894_1894_1048_1047—+—89Sulfolobus  solfataricus P2SSO018001820181018301842435—0245024516311632+—++—91840183910Sulfolobus  tokodaii 7ST02120215021402170218056302110293029301270126+—++—1011Thermofilum  pendens Hrk 5———————————————1112Caldivirga  maquilingensis IC-167CMAQ_1730_1413_1740_1121_1901_1313_1731_0150_0150_1900_1347——+—1213Pyrobaculum  aerophilum IM2PAE06010594058305800589059005850581058105960579++++—1314Pyrobaculum  arsenaticum DSM 13514PARS_2261_2254_2247_2245_2250_2251_2248_2246_2246_2255_2244++—+—1415Pyrobaculum  calidifontis JCM 11548PCAL_2034 _1481_1717_1709_1707_1712_1713_1710_1708_1708_1716_1706++++—1516Pyrobaculum  islandicum DSM 4184PISL_0096_0114_0042_0050_0116_0115_0119_0044_0044_0113_0051——+—1617Thermoproteus  neutrophilus V24StaTNEU_1917_1900_0957_0968_1898_1899_1896_0962_0962_1901_0969—+—1718Archaeoglobus  fulgidus DSM 4304AF19751241197412420116124315921593M1594M24131125+—++—181594M19Haloarcula  marismortui ATCC 43049RRNAC17082628261030863088308717091711171134891363+—++—1920Halobacterium  sp. NRC-1VNG1774G2326G2322G2330G2332G2331G1775C1776G1776G1185G1184GM+—+—2021Halomicrobium  mukohataei DSM 12286HMUK_1612_0914_0925_1105_1103_1104_1613_1614_1614_1679_1984+—+—2122Haloquadratum  walsbyi DSM 16790:HBSQ001HQ3336A3447A3443A3450A3452A3451A3335A3334A3334A——+—++—2223Halorhabdus  utahensis DSM 12940HUTA_2825_1970_1362_0927_0928_1761_1755——_0031_0734+—+—2324Halorubrum  lacusprofundi ATCC 49239HLAC_2132_2622_0015_2252_2254_2253_2131_2130_2130_1215_2081+—++—2425Natronomonas  pharaonis DSM 2160NP4502A1246A0920A1326A1330A1328A4500A4498A4498A1546A1542A+—++—2526Methanobrevibacter  smithii ATCC 35061MSM_0967_1233_1476_0881_1504_1550_0968—————++2627Methanosphaera  stadtmanae DSM 3091MSP_1408_1180_0416_1332_1191_1192_1407—————++2728Methanothermobacter  thermautotrophicus delta HeMTH10122287448741661671013—————+++2829Methanocaldococcus  fervens AG86MEFER_0496_1387_0258_0770_1494_0849_0978—————+2930Methanococcus  jannaschii DSM 2661MJ0143060306430569099409650140—————+++3031Methanocaldococcus  vulcanius M7METVU_1458_0867_0542_0749_0380_1448_0492—————+3132Methanococcus  aeolicus Nankai-3MAEO_0052_1208_0916_1324_0075_0861_0053—————++3233Methanococcus  maripaludis S2MMP0088022412580872039409660089—————++3334Methanococcus  vannielii SBMEVAN_1100_1237_0571_0048_1381_0289_1101—————++3435Methanocorpusculum  labreanum ZMLAB_0523_0525_0524_0526_0073_0526_0522—————++35_035936Methanoculleus  marisnigri JR1MEMAR_0980_0982_0981_0983_0539_0984_0979_0878_0877_0879_0876—++3637Methanospirillum  hungatei JF-1MHUN_2562_2560_2561_2559_2267_2558_2563—————++3738Methanoregula  boonei 6A8MBOO_1238_1235_1236_1234_0514_1233_1239_0957_0956_0958_0955—++3839Methanosphaerula  palustris E1-9cMPAL_1728_1726_1727_1725_2749_1724_1729_2627_2628_2626_2629—+3940Methanosaeta  thermophila PTMTHE_0049_1126_0050_1125_1124_1124—_1134_0047_1135_1123+—++4041Methanococcoides  burtonii DSM 6242MBUR_1229_1227_1228_1226_1106_1105_1230_1232_1231_1233_1236+—+++4142Methanosarcina  acetivorans C2AMA05770581057805823034303305760574057505733035+—++4243Methanosarcina  barkeri fusaroMBAR_A1462_A1464_A1463_A1465_A1792_A1791_A1461_A1459_A1460_A1458_A1793+—++4344Methanosarcina  mazei Goe1MM_1741_1743_1742_1744_0308_0307_1740_1738_1739_1737_0309+—++4445Methanopyrus  kandleri AV19MK0200MKT0801980746155015481495——08960980—+++4546Pyrococcus  abyssi GE5——————————————4647Pyrococcus  furiosus DSM 3638——————————————4748Pyrococcus  horikoshii OT3——————————————4849Thermococcus  gammatolerans EJ3—————————————4950Thermococcus  kodakarensis KOD1——————————————5051Thermococcus  onnurineus NA1—————————————5152Thermococcus  sibiricus MM 739—————————————5253Picrophilus  torridus DSM 9790PTO091802481311024902501435—————+—+—5354Thermoplasma  acidophilum DSM 1728TA05360571095505720573—0652————+—+—5455Thermoplasma  volcanium GSS1TVN05900635110006340633—0924—————+—5556uncultured methanogenic archaeon RC-IRCIX911913912914916915909—————+5657Korarchaeum  cryptofilum OPF8—————————————5758Nanoarchaeum  equitans Kin4-M——————————————5859Nitrosopumilus  maritimus SCM1NMAR_0510_0490_0509_0491_0493_0492_0511_0512_0512———+—5960Escherichia  coli K-12 MG1655B1210015403693805380433683368————+——+—6061Pseudomonas  aeruginosa PAO1PA466639775243526052590510f26110515f0512f—0511f++++—6162Desulfovibrio  vulgaris HildenboroughDVU14613168085618900734073414630854316708550857+—++—62ahemAB gene encoding glutamyl-tRNA reductase in the listed organisms. In the case of two or more copies for one gene, these copies are only listed if none of them is located in the gene cluster, otherwise only the copy within the cluster is listed. This is valid for all listed genes. Genes colocated within the same gene cluster in a certain species are highlighted in the same font (bold, italics, bold italics).bOrganisms possessing heme-containing proteins according to the literature [14–22] are marked with +. For organisms with empty fields the presence of heme-containing proteins was not clearly obvious from literature and bioinformatics data.cOrganisms possessing cobalamin biosynthesis genes according to [35] are marked with +. Organisms that do not synthesize cobalamin are marked with — [35]. For organisms with empty fields the ability to synthesize cobalamin was not obvious from the literature.dOrganisms possessing siroheme-containing sulfite or nitrite reductases based on sequence homology are marked with +. For organisms with empty fields the presence of siroheme-containing sulfite and nitrite reductases was not clearly obvious from literature and bioinformatics data.eNumbers represent the database gene number. The full database gene numbers including the strain abbreviation is given for the hemAB genes. For all other genes the database gene numbers are given without the strain abbreviation. For genes that are located as direct or near neighbors in an individual strain the numbers are written in the same font.fP. aeruginosa  nirE, nirD, nirH, and nirJ genes involved in heme d1 biosynthesis.


## 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 geneshemAB, 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 fivehem 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 volcaniumgenomes. 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 archaealhem 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.Figure 2
Putative heme biosynthesis gene clusters. The most complete gene clusters containing thehem, SUMT, PC2-DH, and ahb-nir genes were found in the genomes of different Pyrobaculum species, members of the Methanosarcinales and Halobacteria (see also Table 1). Genes located as direct or near neighbors are shown in the same color. The number of slashes indicates the number of genes encoding hypothetical proteins which interrupt the gene cluster. The database gene numbers are given in the gene arrow representation without the abbreviation for the strain. These abbreviations are PAE for P. aerophilum, PARS_ for P. arsenaticum, MA for M. acetivorans, MBAR_ for M. barkeri fusaro, and VNG for Halobacterium sp. NRC-1.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 oneSUMT 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 ofnir-Like Genes in Archaeal Heme Biosynthesis

Interestingly, in 32 of the archaeal genomes that contain the earlyhem 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, theahb-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 theahb-nir Genes during Heme Biosynthesis in the Archaea

As already mentioned above theahb-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 hemed1 biosynthesis has not yet been established, several reasonable proposals were made [26, 36–40]. 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 withP. 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).Figure 3
Proposal for the novel, alternative heme biosynthesis pathway in archaea. Archaeal heme biosynthesis starts with the SAM-dependent methylation of UROGEN to precorrin-2 by a SUMT and potentially proceedsvia oxidation of precorrin-2 to sirohydrochlorin by PC2-DH. The side chain modifications (highlighted by dashed circles) including acetate group decarboxylation on rings C and D, acetate group removal on rings A and B, and vinyl group formation on rings A and B are potentially catalyzed by the Ahb-Nir proteins.


### 3.6. Distribution of theahb-nir Genes over the Archaeal Genomes

In accordance with the proposed function of theahb-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 containahb-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 fromS-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

TheM. 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.Purification and characterization of Mba_1791 and Mba_1461. (a) SDS-PAGE analysis of purified Mba_1791 (lane 1) and Mba_1461 (lane 2). Gel permeation chromatography revealed native relative molecular masses of 55,300 ± 840 Da for Mba_1791 (dotted line) and 60,800 ± 7,300 Da for Mba_1461 (solid line), respectively. (b) UV-visible absorption spectrum of extracted tetrapyrroles which accumulated during production of recombinant Mba_1791 inE. coli. (c) UV-visible absorption spectra of enzyme assays after overnight incubation at 37°C in the anaerobic chamber. Uroporphyrinogen III was produced from ALA by the enzymes HemB, HemC, and HemD (dashed double dotted line). Addition of purified Mba_1791 and SAM to the assay mixture resulted in precorrin-2 formation (solid line). Addition of purified NirE, Mba_1461, and NAD+ to the assay resulted in formation of sirohydrochlorin (dotted line). In a coupled enzyme assay containing purified Mba_1791 and Mba_1461 the formation of sirohydrochlorin was also observed (dashed line). For exact details see Section 2.
(a)(b)(c)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 inE. 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, 51–53]. 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 thein 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, 55–57].


### 3.10.M. barkeri Mba_1461 is a PC2-Dehydrogenase

In order to verify the postulated PC2-DH activity for Mba_1461 in anin vitroassay 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 theM. 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).


## 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 geneshemAB, 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 fivehem 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 volcaniumgenomes. 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 archaealhem 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.Figure 2
Putative heme biosynthesis gene clusters. The most complete gene clusters containing thehem, SUMT, PC2-DH, and ahb-nir genes were found in the genomes of different Pyrobaculum species, members of the Methanosarcinales and Halobacteria (see also Table 1). Genes located as direct or near neighbors are shown in the same color. The number of slashes indicates the number of genes encoding hypothetical proteins which interrupt the gene cluster. The database gene numbers are given in the gene arrow representation without the abbreviation for the strain. These abbreviations are PAE for P. aerophilum, PARS_ for P. arsenaticum, MA for M. acetivorans, MBAR_ for M. barkeri fusaro, and VNG for Halobacterium sp. NRC-1.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 oneSUMT 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 ofnir-Like Genes in Archaeal Heme Biosynthesis

Interestingly, in 32 of the archaeal genomes that contain the earlyhem 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, theahb-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 theahb-nir Genes during Heme Biosynthesis in the Archaea

As already mentioned above theahb-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 hemed1 biosynthesis has not yet been established, several reasonable proposals were made [26, 36–40]. 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 withP. 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).Figure 3
Proposal for the novel, alternative heme biosynthesis pathway in archaea. Archaeal heme biosynthesis starts with the SAM-dependent methylation of UROGEN to precorrin-2 by a SUMT and potentially proceedsvia oxidation of precorrin-2 to sirohydrochlorin by PC2-DH. The side chain modifications (highlighted by dashed circles) including acetate group decarboxylation on rings C and D, acetate group removal on rings A and B, and vinyl group formation on rings A and B are potentially catalyzed by the Ahb-Nir proteins.


## 3.6. Distribution of theahb-nir Genes over the Archaeal Genomes

In accordance with the proposed function of theahb-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 containahb-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 fromS-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

TheM. 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.Purification and characterization of Mba_1791 and Mba_1461. (a) SDS-PAGE analysis of purified Mba_1791 (lane 1) and Mba_1461 (lane 2). Gel permeation chromatography revealed native relative molecular masses of 55,300 ± 840 Da for Mba_1791 (dotted line) and 60,800 ± 7,300 Da for Mba_1461 (solid line), respectively. (b) UV-visible absorption spectrum of extracted tetrapyrroles which accumulated during production of recombinant Mba_1791 inE. coli. (c) UV-visible absorption spectra of enzyme assays after overnight incubation at 37°C in the anaerobic chamber. Uroporphyrinogen III was produced from ALA by the enzymes HemB, HemC, and HemD (dashed double dotted line). Addition of purified Mba_1791 and SAM to the assay mixture resulted in precorrin-2 formation (solid line). Addition of purified NirE, Mba_1461, and NAD+ to the assay resulted in formation of sirohydrochlorin (dotted line). In a coupled enzyme assay containing purified Mba_1791 and Mba_1461 the formation of sirohydrochlorin was also observed (dashed line). For exact details see Section 2.
(a)(b)(c)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 inE. 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, 51–53]. 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 thein 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, 55–57].


## 3.10.M. barkeri Mba_1461 is a PC2-Dehydrogenase

In order to verify the postulated PC2-DH activity for Mba_1461 in anin vitroassay 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 theM. 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 hemevia 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.



---

*Source: 175050-2010-12-13.xml*