Diversity and evolution of nitric oxide reduction in bacteria and archaea
Edited by José J. Moura, Universidade Nova de Lisboa, Caparica, Portugal; received September 27, 2023; accepted April 24, 2024 by Editorial Board Member Marcetta Y. Darensbourg
Significance
With the advent of culture-independent techniques for studying environmental microbes, our knowledge of their diversity has exploded, uncovering unique organisms, pathways, and proteins carrying out important processes in the biosphere. Novel biochemical reactions are often proposed based on sequence data, but experimental validation is difficult and rare. In this work, we used environmental sequence data to find enzymes that produce the greenhouse gas N2O from NO and validated our hypothesis with experiments. These new enzymes likely contribute to global N2O fluxes and expand the breadth of nitrogen cycling. We also demonstrated that these enzymes evolved multiple times from oxygen reductases, indicating that the evolutionary histories of aerobic respiration and denitrification—and more broadly the oxygen and nitrogen cycles—are tightly connected.
Abstract
Nitrous oxide is a potent greenhouse gas whose production is catalyzed by nitric oxide reductase (NOR) members of the heme-copper oxidoreductase (HCO) enzyme superfamily. We identified several previously uncharacterized HCO families, four of which (eNOR, sNOR, gNOR, and nNOR) appear to perform NO reduction. These families have novel active-site structures and several have conserved proton channels, suggesting that they might be able to couple NO reduction to energy conservation. We isolated and biochemically characterized a member of the eNOR family from the bacterium Rhodothermus marinus and found that it performs NO reduction. These recently identified NORs exhibited broad phylogenetic and environmental distributions, greatly expanding the diversity of microbes in nature capable of NO reduction. Phylogenetic analyses further demonstrated that NORs evolved multiple times independently from oxygen reductases, supporting the view that complete denitrification evolved after aerobic respiration.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
The heme-copper oxidoreductase (HCO) superfamily is extremely diverse, with members playing crucial biogeochemical roles in both aerobic (oxygen reductases) and anaerobic [nitric oxide reductases (NORs)] respiration (1–3). While NO reduction can also be performed by fungal NORs (4) and flavodiiron proteins (5), in this paper we focus on NORs from the HCO superfamily. Fungal NO reduction is performed by cytochrome P450 (6), and flavodiiron proteins are primarily used for detoxification of NO. Respiratory denitrification in both Bacteria and Archaea involves NORs from the HCO superfamily. The HCO superfamily consists of three well-characterized oxygen reductase families (A, B, and C) and three NOR families (cNOR, qNOR, and qCuANOR) (1–3). The oxygen reductases catalyze the reduction of O2 to water (O2 + 4eout− + 4Hin+ + nHin+ → 2H2O + nHout+) and share a conserved reaction mechanism (3, 7), wherein three of the electrons required to reduce O2 are provided by the active-site metals, heme-Fe and CuB, while the fourth electron is derived from a unique redox-active cross-linked histidine-tyrosine cofactor (8) (Fig. 1). The free energy available from this reaction is converted into a transmembrane proton electrochemical gradient, allowing microbes to harness energy from aerobic respiration. The generation of electrochemical gradient occurs via two different mechanisms: charge separation across the membrane and proton pumping (9, 10). Both the protons used for chemistry (i.e., O2 reduction to water) and those separately pumped protons are taken up from the electrochemically negative side of the membrane (bacterial cytoplasm) by conserved proton-conducting channels that are composed of conserved polar residues and internal water molecules. The different oxygen reductase families exhibit differential proton pumping stoichiometries (n = 4 for the A-family, and n = 2 for the B and C-families) (10–12) and thus conserve energy differentially depending on their proton channels—though there is some dispute regarding the proton pumping stoichiometry of the C-family, with some studies reporting n = 4 (13). The oxygen reductases also vary in their secondary subunits that function as redox relays from electron donors in the electron transport chain (e.g., cytochrome c) to the protein complex active site, with the A and B-families utilizing a CuA-containing subunit (14–16) and the C-family containing one or more cytochrome c subunits (17) (Fig. 1).
Fig. 1.

NORs catalyze the reduction of NO to nitrous oxide (2NO + 2Hout+ + 2eout− + nHin+ → N2O + H2O + nHout+). NO reduction requires 2 molecules of NO to form nitrous oxide. With each N atom decreasing in oxidation state by 1, it is only a 2-electron reaction and does not require the cross-linked histidine-tyrosine cofactor for catalysis (18)— providing one metric for identifying putative NORs from environmental sequence data. There are currently three biochemically characterized NORs within the HCO superfamily, the cNOR, qNOR, and qCuANOR. The cNOR and qNOR families have a four amino acid coordinated FeB ion in their active sites, in contrast to the three amino acid coordinated CuB found in the HCO oxygen reductases (18–20). The cNOR and qNOR families are closely related to the C-family oxygen reductases (21). Like the C-family O2 reductase, cNOR has a secondary cytochrome c subunit, while qNOR appears to be the result of a gene fusion of the primary and secondary subunits forming a single polypeptide that lacks the heme c binding motif (22–24). cNOR does not conserve energy, with the enzyme taking both electrons and protons for NO reduction from the periplasmic side (25). Although qNOR is proposed to take up protons from the cytoplasm for NO reduction (20, 26), it does not have conserved residues that could form a proton channel from the cytoplasm, and it is not clear whether this enzyme conserves energy via either charge separation or proton pumping. The qCuANOR from Bacillus azotoformans (2, 27) is not closely related to cNOR and qNOR and is instead derived from within the B-family O2 reductases, leading it to be reclassified as the bNOR family (23). bNOR is fundamentally different from cNOR and qNOR, containing a CuA cofactor in the secondary subunit and a conserved proton channel for proton uptake from the cytoplasm. bNOR was shown to be electrogenic (27) and thus capable of generating more energy than previously characterized NORs. In earlier work, Hemp and Gennis demonstrated that the HCO superfamily was more diverse than previously thought, working with data from archaeal genomes (1). With recent work in a larger dataset including Bacteria and Archaea, we expanded that diversity to 12 families and demonstrated that quinol oxidation evolved within the HCO superfamily multiple times (23). In this work, we used phylogenomics of both isolates and environmental sequence data to study the diversity and evolution of multiple putative NOR families (eNOR, gNOR, nNOR, and sNOR) within the HCO superfamily and verified the biochemical NO reduction activity of eNOR from the bacterium, Rhodothermus marinus. We also identified a new family of putative N2O reductases. Our findings expanded the number of denitrification pathways in Bacteria and Archaea, increased the breadth of modern N2O production and further constrained the evolutionary history of one of the key protein scaffolds involved in aerobic and anaerobic respiration.
Results and Discussion
Expansion of the HCO Superfamily.
Phylogenomic analyses of genomic and metagenomic data identified at least six new families belonging to the HCO superfamily (Figs. 1 and 2) that are missing the active-site tyrosine, indicating that they do not catalyze O2 reduction. Analysis of structural models and sequences (SI Appendix) for each of these families showed no evidence for the sequence migration of a conserved tyrosine that could form an active-site cross-linked cofactor, as was observed in the C-family O2 reductases (8). Furthermore, their active sites exhibited structural features never before seen within the superfamily (Fig. 1). One of these families is closely related to qNOR and has been proposed to be a NO dismutase (NOD) contributing to intracellular O2 production in “Candidatus Methylomirabilis oxyfera” (28, 29). Another family is closely related to cNOR and might serve as a unique sulfide and acetylene-insensitive nitrous oxide reductase (N2O red) (30, 31). The remaining four families (eNOR, sNOR, nNOR, and gNOR) are closely related to the B-family of O2 reductases (Fig. 2) and encode homologs of the CuA-containing secondary subunits. This is consistent with the presence of CuA-containing subunits in the B-family of oxygen reductases (Fig. 1 and Dataset S1). Based on modeled active-site structures and genomic context, we proposed that these novel families perform NO reduction (Fig. 1).
Fig. 2.

Biochemical Characterization of eNOR.
To validate these predictions, we isolated and biochemically characterized a member of the eNOR family from R. marinus DSM 4252, a thermophilic member of the Bacteroidetes phylum. R. marinus was originally classified as a strict aerobe (34), but its genome encoded a periplasmic nitrate reductase (NapA), two nitrite reductases (NirS and NirK), and a N2O reductase (NosZ), suggesting that it may also be capable of denitrification (SI Appendix, Fig. S1). Denitrification was not observed under strictly anaerobic conditions, however, under microoxic conditions, we observed that isotopically labeled 15NO3− was converted to 30N2 (SI Appendix, Fig. S2), demonstrating that R. marinus DSM 4252 was capable of complete aerobic denitrification (NO3−→N2). Blockage of the N2O red (NosZ) with acetylene led to the accumulation of N2O (Fig. 3), implying that a NOR was also present in R. marinus DSM 4252. No known NORs (cNOR, qNOR, qCuANOR/bNOR, or flavodiiron proteins) were found in the genome. However, R. marinus DSM 4252 encoded a member of the eNOR family (SI Appendix, Fig. S1).
Fig. 3.

Isolation and biochemical characterization of the R. marinus DSM 4252 eNOR protein verified that it catalyzed NO reduction [at 25 °C, kcat = 0.68 ± 0.21 NO s−1 (n = 4)] (Fig. 3). This turnover number is lower than catalytic turnover rates reported for NORs purified from mesophilic bacteria such as Pseudmonas stutzeri [16 NO s−] (35) or Neisseria meningitidis [30 NO s−] (20) but is higher than activities reported for cNOR purified from other thermophilic microorganisms such as Thermus thermophilus [0.09 NO s−] (25). eNOR was unable to catalyze O2 reduction using a range of electron donors (SI Appendix, Fig. S3), showing that it functioned solely as a NOR. UV–Vis spectroscopy and heme characterization via mass spectrometry demonstrated that the R. marinus DSM 4252 eNOR contained a unique modified heme a that is used in both heme sites (Fig. 3 and SI Appendix, Figs. S3 and S4). Another member of the eNOR family was previously isolated from the aerobic denitrifier Magnetospirillum magnetotacticum MS-1 (36, 37); however, its function was never determined. The UV–Vis spectra of the M. magnetotacticum eNOR (36) were identical to the R. marinus eNOR, implying that the modified heme a is a general feature of the family. Mass spectroscopic analysis of the hemes extracted from eNOR revealed that this heme was As—a previously isolated heme a with a hydroxyethylgeranylgeranyl side chain first identified in the B-family oxygen reductase from Sulfolobus acidocaldarius (38). Many eNOR operons contain a CtaA homolog, an O2-dependent enzyme that converts heme o to heme a (39). This is consistent with the observation that eNOR required microoxic conditions to be expressed. Organisms performing denitrification with eNOR appear to be obligate aerobic denitrifiers, and future work will establish the extent of their role in environmental aerobic denitrification (40).
Unique Active-Site Features of recentiy identified NORs.
In addition to the experimental evidence that both eNOR and bNOR enzymes are NO reductases, there are several reasons to predict that the other recently identified families also perform NO reduction. The sNOR family has the same active-site structure as the bNOR family, strongly suggesting that it also performs NO reduction. However, the sNOR and bNOR families are not closely related and mark an example of convergent evolution of active-site structures within the HCO superfamily (Figs. 1 and 2). We identified another example of convergent evolution in the nNOR family. Members of this family have the same conserved active-site residues as the cNOR and qNOR families (Fig. 1) but are only very distantly related to them; nNOR is related to the B-family, whereas cNOR and qNOR are related to the C-family. Interestingly, the low-spin heme in nNOR is ligated by a histidine and methionine, which likely raises its redox potential by ~150 mV (41). This feature is similar to a modification found in some eNORs, wherein the low-spin heme is ligated by histidine and lysine. The gNOR is the first example of a HCO family member that has replaced one of the active-site histidines—residues completely conserved in all other families. The gNOR active site, with an aspartate in place of histidine, is likely capable of catalyzing NO reduction, since NO reduction was demonstrated in a bioinorganic mimic of the gNOR active-site (42). The gNOR has a secondary subunit with a cupredoxin fold that is missing the residues required to bind CuA, similar to the quinol-oxidizing oxygen reductase cytochrome bo3 from Escherichia coli. Conserved residues that could bind quinol have been identified in gNOR, so it may be a quinol-oxidizing NOR similar to qNOR (23).
The biochemically characterized eNOR and bNOR and proposed sNOR and gNOR families within the HCO superfamily have active sites that differ significantly from those found in the well-characterized cNOR and qNOR enzymes (Fig. 1). Importantly, while oxygen reduction chemistry is constrained to require a redox-active tyrosine cofactor, multiple HCO active-site structures appear to be compatible with NO reduction chemistry. Oxygen reductases from the A-family, B-family, and C-family (18, 43, 44) appear to catalyze NO reduction, albeit less efficiently and with a different mechanism than the NORs. This difference in biochemical constraints between NO and oxygen reduction chemistry suggests that the evolutionary transition of oxygen reduction to NO reduction is relatively simple within the HCO superfamily. Another useful chemical constraint that appears to differentiate the catalysis of O2 reduction and NO reduction is the active site metal: in the currently characterized HCOs, CuB is utilized for O2 reduction chemistry, whereas FeB is used for NO reduction chemistry. If this pattern is verified for the other predicted NOR families, it would indicate that the chemistry performed by HCOs is determined to a certain degree by the electronic properties of the active-site metal. It is important to note that the above biochemical constraints for NO and O2 reduction chemistry are only applicable within the HCO superfamily, since other enzymes such as the flavodiiron proteins (5) or cytochrome bd oxygen reductases (45) are capable of NO and O2 reduction with entirely different active site characteristics. Interestingly, tryptophan/tyrosine chains that are predicted to prevent oxidative damage in redox-active proteins (46) are conserved in both the A-family and B-family O2 reductases and several of the NO reductases that have evolved from the B-family (bNOR, eNOR, and gNOR) (Dataset S5). In these chains, radicals generated during substrate turnover move by hole hopping through a series of tryptophan and tyrosine residues to the surface of the protein where they are safely quenched by redox buffers within the cell (e.g., glutathione). Despite the difference in catalytic mechanisms between O2 reductases, 2 out of 3 residues implicated in oxidative protection are found in bNOR, eNOR, and gNOR; they are missing in sNOR and nNOR. Future investigation of the catalytic differences between these NORs will provide insight into the role played by these residues.
Bioenergetics of Denitrification Pathways with Recently Identified NORs.
Although both denitrification and aerobic respiration are highly exergonic processes, most of the enzymes in the denitrification pathway are not directly coupled to energy conservation in cells, making denitrification less energetically efficient than aerobic respiration (47). In the HCO oxygen reductases, conserved proton channels deliver protons from the cytoplasm to the active site for chemistry. These same channels are used to pump protons to the periplasmic side (9, 11, 12, 48). In contrast, previously characterized NORs do not appear to pump protons or conserve as much energy as the oxygen reductases. cNOR does not have conserved proton channels from the cytoplasm, which makes this enzyme incapable of conserving energy (25). The evidence regarding qNOR is currently ambiguous: although qNOR does not have conserved proton channels either, there is some evidence suggesting that it can conserve energy (20, 26).
We found that eNOR family has conserved hydrophilic residues, similar to the electrogenic bNOR, that closely resemble those found in the proton-conducting K-channel within the B-family of oxygen reductases (12, 27) (Dataset S1 and SI Appendix, Fig. S5). The sNOR family also has conserved residues in the K-channel region. However, this putative proton channel is slightly different from those found in the B1-subfamily of O2 reductases (that contains the T. thermophilus ba3) and the eNOR and bNOR families (Dataset S1). The conserved serine (S309, B1-subfamily T. thermophilus ba3 numbering) found in those families is missing in the sNOR, and instead, this enzyme has a conserved glutamate residue in a structurally different location from other HCO enzymes. Interestingly, the nNOR family, which has the same active site as cNOR and qNOR, also has a conserved proton channel (Dataset S1 and SI Appendix, Fig. S5). This implies that these recently identified NORs may be capable of energy conservation and the lack of a proton channel in the cNOR and qNOR may not be due to energetic constraints universal to NO reduction (49). The conserved proton channels in the eNOR, bNOR, sNOR, and nNOR families would allow them to conserve energy via charge separation and potentially by proton pumping. Detailed characterization of these new NOR families will be helpful for understanding the mechanism of proton pumping in the HCO superfamily—one of the longstanding questions in bioenergetics (50).
Environmental Distribution of NORs.
The recently identified HCO NOR families have broad phylogenetic and environmental distributions that substantially expand the scope of denitrification occurring in nature (Table 1 and Datasets S6 and S7). The eNOR, sNOR, gNOR, and nNOR families are all found in both Bacteria and Archaea, whereas the bNOR family was only found in the Bacillales order of Firmicutes (Dataset S6). Phylogenetic analysis of metagenomic data shows that the majority of eNOR, sNOR, gNOR, and nNOR enzymes appear in uncharacterized taxa, hinting at many more organisms capable of NO reduction than previously suspected. Furthermore, the new HCO NOR families were found in a wide variety of environments (Table 1 and Dataset S7). sNORs are broadly distributed in many environments, however, they are rarely found in Archaea. sNORs are found in most ammonia-oxidizing bacteria (AOB) sequenced to date, suggesting that the capability of NO reduction is an important difference in ammonia oxidation pathways between Bacteria and Archaea. Given the importance of AOB, it is likely that sNOR plays a role in this key biogeochemical process in nature (Dataset S7). The gNORs were predominantly found in microbes inhabiting sulfidic environments and as mentioned above may reflect an adaptation that allows for denitrification in the presence of free sulfide, which inhibits other NOR families (Dataset S7). Our analyses revealed that the eNOR family—the new HCO enzyme characterized in detail here—is extremely common in nature and has a broad distribution, similar to the cNOR and qNOR families (Datasets S6 and S7). eNORs were found in many strains of Candidatus Accumulibacter phosphatis, a critical microbe utilized in wastewater treatment plants for enhanced biological phosphorus removal. The eNOR is highly expressed in transcriptomic datasets from these facilities, demonstrating that Ca. Accumulibacter phosphatis is capable of complete denitrification in situ (51). eNOR has also been found in microbes capable of performing autotrophic nitrate reduction coupled to Fe(II) oxidation (NRFO). Gallionellaceae KS and related strains express an eNOR under denitrifying conditions, suggesting that an individual organism is capable of complete NRFO (52). eNOR is also common in hypersaline environments (Dataset S7), where it might play a role in the adaptation of denitrification to high salt conditions.
Table 1.
NCBI-Genomes | IMG-metagenomes | GTDB-genomes | |
---|---|---|---|
A-family | 20,290 | 102,368 | 45,135 |
B-family | 1,238 | 4,683 | 2,021 |
C-family | 13,976 | 23,015 | 14,981 |
qNOR | 4,388 | 7,680 | 3,458 |
cNOR | 2,801 | 4,824 | 2,594 |
eNOR | 68 | 2,709 | 547 |
sNOR | 95 | 872 | 344 |
bNOR | 51 | 12 | 200 |
nNOR | 6 | 289 | 32 |
gNOR | 10 | 913 | 156 |
NOD | 8 | 539 | 108 |
N2O red | 25 | 597 | 293 |
Distribution of NOR families in sequenced genomes vs. environmental datasets. The recently found NOR families account for approximately 2/3 of currently known diversity and 1/2 of the abundance of NORs in nature.
Many organisms encode NORs from multiple families (e.g., Candidatus Methylomirabilis oxyfera has qNOR, sNOR, and gNOR; B. azotoformans has qNOR, sNOR, and bNOR). The reasons underlying this apparent redundancy remain unclear, but it suggests that selection for different enzymatic properties (NO affinity, enzyme kinetics, energy conservation, or sensitivity to inhibitors) or the concentration of O2 may be important factors in determining their distribution and use, similar to what was observed for the HCO oxygen reductase families (10). Analysis of the presence of denitrification genes (nitrite reductases, NORs, and the NosZ-type N2O reds) within sequenced genomes revealed that many more organisms are capable of complete denitrification than previously realized (Dataset S8). Our current understanding of the diversity of organisms capable of performing denitrification in nature is far from complete but stands to grow with the recognition of these new families of NORs.
The Intertwined Evolutionary History of Aerobic Respiration and Denitrification.
Combining our biochemical results and insights with the phylogenetic relationships among different groups in the HCO superfamily, which contains both oxygen reductases and NORs, allowed us to better ordinate the evolutionary histories of aerobic respiration and denitrification. Previous work had demonstrated close evolutionary relationships between the A and B-family oxygen reductases (3), as well as close ancestry between the C-family oxygen reductases, cNOR, and qNOR (21). Yet, the question of which arose first—denitrification or aerobic respiration—has been harder to resolve.
Our analysis of the distribution of oxygen reductases and NORs across the wide diversity of microbial life revealed that oxygen reductases are far more widely distributed; over 30,238 of the 47,894 species in the genome taxonomy database (GTDB) encoded oxygen reductases, whereas NORs were only found in 6,626 species (Datasets S6 and S8). This distribution illustrated the massive impact that oxygen has had on the energetics of the biosphere. The A-family is by far the most widely distributed of the HCO enzymes. It is found in all three domains of life and in more phyla than any other enzyme of the HCO superfamily. This is consistent with the view that the A-family oxygen reductase holds greater antiquity (1, 10, 32, 53), wherein the B-family and C-family O2 reductases each evolved independently from within the A-family to facilitate specific metabolic and ecological challenges associated with exotic flavors of aerobic biology, like in hyperthermophiles (15) or in chemoautotrophic iron oxidation (54). An evolutionary transition from the A-family to B-family due to selection for higher oxygen affinity, which led to the loss of the D-proton channel to facilitate greater access of oxygen to the active site has been inferred from structural and phylogenetic data with a putative intermediate enzyme suggested in Nitrosopumilis maritimus (10). With that in mind, a clear transition from A to the C-family oxygen reductases has not been demonstrated.
The C-family branch of the HCO superfamily consists of the closely related C-family oxygen reductase, N2O red, cNOR, qNOR, and NOD (Fig. 2) (3, 21, 29, 33, 55). The transition from cNOR to qNOR as the result of a gene fusion of subunits I and II has been reasonably inferred based on the sequence similarity between the N-terminal domain of qNOR and subunit II of cNOR (24). That NOD was derived from qNOR is supported by the high level of sequence similarity between qNOR and NOD, as well as its branching topology within the qNOR clade (Fig. 2). These relative constraints support a simple interpretation of the evolutionary history within this branch—the C-family is the oldest, followed by the evolution of the N2O red and cNOR. The cNOR then is the ancestor of qNOR, followed by NOD. The sparse distribution of the C-family oxygen reductases, cNOR, and NOD in Archaea supports the hypothesis that these families evolved after the A-family oxygen reductases. While qNOR is widely distributed, it is rarely associated with energetically efficient denitrification and has been proposed to function in nature as a detoxification enzyme (21). Therefore, the presence of qNOR cannot be used as a robust constraint for the antiquity of denitrification. The wide distribution of the A-family, the indications of an evolutionary transition from the A- to B-family, and the relatively sparse distribution of the C-family branch members in Archaea, all suggest that the A-family likely hosts the root of the HCO superfamily (Fig. 2). In the future, the use of different comparative biological approaches—particularly those that might better capture the evolution of paralogs—to root the phylogenetic tree of HCOs could be used to test this idea.
What is clear from our new observations—regardless of the placement of the root of the HCO superfamily—is that NORs have evolved independently multiple times from the B-family and C-family oxygen reductases (Fig. 2). There are key underlying factors that enabled this, both chemical and environmental. It is biochemically straightforward to adapt an oxygen reductase (4 e− chemistry) for NO reduction (2 e− chemistry). B and C-family oxygen reductases can reduce NO at high concentrations in vitro (43, 56). It is thus unsurprising that small evolutionary modifications would lead to a cascade of enzyme descendants each capable of NO reduction at lower concentrations to enable more effective denitrification. It is also clear that in many environments, denitrification and aerobic respiration often co-occur (57, 58), and that many microorganisms display the respiratory flexibility to shift from aerobic to anaerobic respiration, especially at lower O2 concentrations (59). This respiratory flexibility is reflected in the fact that denitrification and oxygen respiration share much of the same bioenergetic logic, conserving energy via complex III or alternative complex III (60). Thus, the biochemical promiscuity of O2 reductases toward NO, the ecological proximity of NO3- and O2, and the close similarity between their respiratory pathways help explain why the evolutionary transition of O2 reduction to NO reduction is both favorable and readily achievable.
Nitrate is derived from biogeochemical processes involving oxygen (61) and consequently opportunities for denitrification prior to Earth’s great oxygenation event (GOE) were muted compared to those after the GOE with the rise of nitrate in seawater. Finally, both aerobic respiration and denitrification are constrained by the presence of copper. Copper is an essential bioinorganic component of the active site in the O2 reductases, nitrite reductase NirK, and for N2O red (NosZ) and is therefore essential for their biochemical activity. The environmental abundance of copper increased significantly after the GOE (62–64), suggesting that both of these metabolic pathways evolved and expanded after the accumulation of O2 in Earth’s surface environments.
Materials and Methods
This is an abbreviated version of the materials and methods used in this work. A detailed version is available as part of SI Appendix.
Purification of eNOR from R. marinus Grown under Denitrification Conditions.
R. marinus DSMZ 4252 was grown in DSM Medium 630 with 30 mM NO3− added and shaken at 75 rpm to induce denitrification. The microoxic conditions that result from slow shaking were essential for denitrification in R. marinus. We used labeled nitrate (15NO3−) to verify that R. marinus DSM 4252 was capable of complete denitrification (NO3− to N2). The experiments detailed below established that eNOR was expressed under these conditions and functions as a NOR.
Cultures of R. marinus were grown in 1L of this medium in 24 × 2L Erlenmeyer flasks for 36 h or to stationary phase, to generate sufficient biomass for protein purification. The cell pellet recovered from this culture was subject to lysis as described in SI Appendix, Materials and Methods and the membrane fraction was recovered by ultracentrifugation. eNOR was purified in a protocol similar to that described for caa3 from R. marinus (65). Purification of eNOR was improved when the membranes were first solubilized in 1% 3-[(3-Cholamidopropyl) dimethylammonio]-1-propanesulfonate (CHAPS), apparently recovering peripheral membrane proteins. The membranes not solubilized in this step were then pelleted with ultracentrifugation and solubilized in 1% N-dodecyl-β-D-maltoside (DDM). These solubilized membrane proteins were subject to a protein purification protocol detailed in SI Appendix, Materials and Methods. Purified protein was confirmed to be eNOR using electrophoretic analysis and mass spectrometric identification.
Biochemical Characterization of eNOR.
Purified eNOR was verified to perform NO reduction by measuring NO consumption using a Clark Electrode (World Precision Instruments) in a protocol previously described (25) and by measurement of the product, N2O using GC.
The heme cofactors of eNOR were first analyzed using a pyridine hemochrome assay (66) and then analyzed by Liquid Chromatography-Mass Spectrometry (LC-MS) after solvent extraction. Further details for both assessment of activity and cofactor identification are available in SI Appendix, Materials and Methods.
A detailed description of the phylogenomic analysis of NORs by taxonomy and environment is provided in SI Appendix, Materials and Methods.
Data, Materials, and Software Availability
All the protein accession numbers used for generation of trees in this study, as well as associated phylogenetic trees and multiple sequence alignments are included in the supporting information. The Hidden Markov Models (HMMs) used for identification of HCO sequences can be found at https://github.com/ranjani-m/HCO (67).
Acknowledgments
We thank the NIH for funding support (grant# U12AB123456 to PI: R.B.G.). This research was also supported by funding from the Agouron Institute (W.W.F. and J.H.) and by the Community Science Project 507064 (PI: R.H.) under the Joint Genome Institute (https://ror.org/04xm1d337), which is a Department of Energy (DOE) Office of Science User Facility. Resources were also used at Office of Biological and Environmental Research of the US Department of Energy Atmospheric System Research Program Interagency Agreement grant DE-AC02-05CH11231 (JGI). Resources were used at Office of Biological and Environmental Research of the United States Department of Energy Atmospheric System Research Program Interagency Agreement grant DE-AC05-76RL01830 Environmental Molecular Sciences Laboratory (EMSL). We thank Sylvia Choi for providing pure ba3 oxygen reductase from T. thermophilus to use as a control for oxygen reductase assays and for heme extraction, Paige Sheridan for providing purified qNOR from Persephonella marina, Lici Schurig-Briccio for guidance in performing NOR assays with the Clark Electrode, and Peter Yau at the University of Illinois’ Mass spectrometric facility for protein identification. We thank Alon Philosof and Connor Skennerton for valuable discussions on bioinformatics analysis. Finally, we would like to thank our reviewers for their valuable comments that led to the improvement of this manuscript.
Author contributions
R.M., W.W.F., R.B.G., and J.H. designed research; R.M., L.A.P., R.A.S., L.M.W., M.M.L., R.H., W.W.F., and J.H. performed research; R.M., R.A.S., L.M.W., W.W.F., R.B.G., and J.H. contributed new reagents/analytic tools; R.M., L.A.P., U.F.L., W.W.F., R.B.G., and J.H. analyzed data; and R.M., R.A.S., L.M.W., M.M.L., R.H., W.W.F., R.B.G., and J.H. wrote the paper.
Competing interests
The authors declare no competing interest.
Supporting Information
Appendix 01 (PDF)
- Download
- 19.32 MB
Dataset S01 (XLSX)
- Download
- 9.46 KB
Dataset S02 (PDF)
- Download
- 5.77 MB
Dataset S03 (RTF)
- Download
- 50.25 KB
Dataset S04 (XLSX)
- Download
- 17.18 KB
Dataset S05 (XLSX)
- Download
- 9.21 KB
Dataset S06 (XLSX)
- Download
- 5.77 MB
Dataset S07 (XLSX)
- Download
- 11.50 KB
Dataset S08 (XLSX)
- Download
- 829.27 KB
References
1
J. Hemp, R. B. Gennis, Diversity of the heme-copper superfamily in archaea: Insights from genomics and structural modeling. Results Probl. Cell Differ. 45, 1–31 (2008).
2
Suharti, M. J. F. Strampraad, I. Schröder, S. de Vries, A novel copper A containing menaquinol NO reductase from Bacillus azotoformans. Biochemistry 40, 2632–2639 (2001).
3
M. M. Pereira, F. L. Sousa, A. F. Veríssimo, M. Teixeira, Looking for the minimum common denominator in haem–copper oxygen reductases: Towards a unified catalytic mechanism. Biochim. Biophys. Acta 1777, 929–934 (2008).
4
K. Maeda et al., N2O production, a widespread trait in fungi. Sci. Rep. 5, 9697 (2015).
5
C. V. Romão, J. B. Vicente, P. T. Borges, C. Frazão, M. Teixeira, The dual function of flavodiiron proteins: Oxygen and/or nitric oxide reductases. J. Biol. Inorg. Chem. 21, 39–52 (2016).
6
N. Takaya, H. Shoun, Nitric oxide reduction, the last step in denitrification by Fusarium oxysporum, is obligatorily mediated by cytochrome P450nor. Mol. Gen. Genet. 263, 342–348 (2000).
7
J. Hemp, C. Christian, B. Barquera, R. B. Gennis, T. J. Martínez, Helix switching of a key active-site residue in the cytochrome cbb3 oxidases. Biochemistry 44, 10766–10775 (2005).
8
J. Hemp, D. E. Robinson, T. J. Martinez, N. L. Kelleher, R. B. Gennis, The evolutionary migration of a post-translationally modified active-site residue in the proton-pumping heme-copper oxygen reductases. Biochemistry 45, 15405–15410 (2006).
9
M. K. F. Wikstrom, Proton pump coupled to cytochrome c oxidase in mitochondria. Nature 266, 271–273 (1977).
10
H. Han et al., Adaptation of aerobic respiration to low O2 environments. Proc. Natl. Acad. Sci. U.S.A. 108, 14109–14114 (2011).
11
J. Hemp et al., Comparative genomics and site-directed mutagenesis support the existence of only one input channel for protons in the C-family (cbb3 oxidase) of heme-copper oxygen reductases. Biochemistry 46, 9963–9972 (2007).
12
H.-Y. Chang et al., Exploring the proton pump and exit pathway for pumped protons in cytochrome ba3 from Thermus thermophilus. Proc. Natl. Acad. Sci. U.S.A. 109, 5259–5264 (2012).
13
V. Rauhamäki, D. A. Bloch, M. Wikström, Mechanistic stoichiometry of proton translocation by cytochrome cbb3. Proc. Natl. Acad. Sci. U.S.A. 109, 7286–7291 (2012).
14
L. Qin, C. Hiser, A. Mulichak, R. M. Garavito, S. Ferguson-Miller, Identification of conserved lipid/detergent-binding sites in a high-resolution structure of the membrane protein cytochrome c oxidase. Proc. Natl. Acad. Sci. U.S.A. 103, 16117–16122 (2006).
15
T. Soulimane et al., Structure and mechanism of the aberrant ba3-cytochrome c oxidase from Thermus thermophilus. EMBO J. 19, 1766–1776 (2000).
16
T. Tiefenbrunn et al., High resolution structure of the ba3 cytochrome c oxidase from Thermus thermophilus in a lipidic environment. PLoS One 6, 1–12 (2011).
17
S. Buschmann et al., The structure of cbb3 cytochrome oxidase provides insights into proton pumping. Science 329, 327–330 (2010).
18
N. Lehnert et al., The biologically relevant coordination chemistry of iron and nitric oxide: Electronic structure and reactivity. Chem. Rev. 121, 14682–14905 (2021).
19
T. Hino et al., Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330, 1666–1670 (2010).
20
N. Gonska et al., Characterization of the quinol-dependent nitric oxide reductase from the pathogen Neisseria meningitidis, an electrogenic enzyme. Sci. Rep. 8, 3637 (2018).
21
W. G. Zumft, Nitric oxide reductases of prokaryotes with emphasis on the respiratory, heme–copper oxidase type. J. Inorg. Biochem. 99, 194–215 (2005).
22
Y. Matsumoto et al., Crystal structure of quinol-dependent nitric oxide reductase from Geobacillus stearothermophilus. Nat. Struct. Mol. Biol. 19, 238–245 (2012).
23
R. Murali, J. Hemp, R. B. Gennis, Evolution of quinol oxidation within the heme-copper oxidoreductase superfamily. Biochim. Biophys. Acta 1863, 148907 (2022).
24
J. Hendriks et al., Nitric oxide reductases in bacteria. Biochim. Biophys. Acta 1459, 266–273 (2000).
25
L. A. Schurig-Briccio et al., Characterization of the nitric oxide reductase from Thermus thermophilus. Proc. Natl. Acad. Sci. U.S.A. 110, 12613–12618 (2013).
26
C. C. Gopalasingam et al., Dimeric structures of quinol-dependent nitric oxide reductases (qNORs) revealed by cryo–electron microscopy. Sci. Adv. 5, eaax1803 (2019).
27
S. Al-Attar, S. de Vries, An electrogenic nitric oxide reductase. FEBS Lett. 589, 2050–2057 (2015).
28
K. F. Ettwig et al., Nitrite-driven anaerobic methane oxidation by oxygenic bacteria. Nature 464, 543–548 (2010).
29
K. F. Ettwig et al., Bacterial oxygen production in the dark. Front. Microbiol. 3, 273 (2012).
30
A. M. Jones, A. M. Adkins, R. Knowles, G. R. Rayat, Identification of a denitrifying gliding bacterium, isolated from soil and able to reduce nitrous oxide in the presence of sulfide and acetylene, as Flexibacter canadensis. Can. J. Microbiol. 36, 765–770 (1990).
31
A. M. Jones, T. C. Hollocher, R. Knowles, Nitrous oxide reductase of Flexibacter canadensis: A unique membrane-bound enzyme. FEMS Microbiol. Lett. 92, 205–209 (1992).
32
M. M. Pereira, M. Santana, M. Teixeira, A novel scenario for the evolution of haem–copper oxygen reductases. Biochim. Biophys. Acta 1505, 185–208 (2001).
33
A.-L. Ducluzeau et al., The evolution of respiratory O2/NO reductases: An out-of-the-phylogenetic-box perspective. J. R. Soc. Interface 11, 20140196 (2014).
34
G. A. Alfredsson, J. K. Kristjansson, S. Hjrleifsdottir, K. O. Stetter, Rhodothermus marinus, gen. nov., sp. nov., a Thermophilic, Halophilic Bacterium from submarine hot springs in Iceland. J. Gen. Microbiol. 134, 299–306 (1988).
35
B. Heiss, K. Frunzke, W. G. Zumft, Formation of the N-N bond from nitric oxide by a membrane-bound cytochrome bc complex of nitrate-respiring (denitrifying) Pseudomonas stutzeri. J. Bacteriol. 171, 3288–3297 (1989).
36
H. Tamegai, T. Yamanaka, Y. Fukumori, Purification and properties of a ‘cytochrome a1’-like hemoprotein from a magnetotactic bacterium, Aquaspirillum magnetotacticum. Biochim. Biophys. Acta 1158, 237–243 (1993).
37
Y. Tanimura, Y. Fukumori, Heme-copper oxidase family structure of Magnetospirillum magnetotacticum ‘cytochrome a1’-like hemoprotein without cytochrome c oxidase activity. J. Inorg. Biochem. 82, 73–78 (2000).
38
M. Lübben, K. Morand, Novel prenylated hemes as cofactors of cytochrome oxidases. Archaea have modified hemes A and O. J. Biol. Chem. 269, 21473–21479 (1994).
39
K. R. Brown, B. M. Allan, P. Do, E. L. Hegg, Identification of novel hemes generated by heme A synthase: Evidence for two successive monooxygenase reactions. Biochemistry 41, 10906–10913 (2002).
40
B. Ji et al., Aerobic denitrification: A review of important advances of the last 30 years. Biotechnol. Bioproc. E 20, 643–651 (2015).
41
C. Fufezan, J. Zhang, M. R. Gunner, Ligand preference and orientation in b- and c-type heme-binding proteins. Proteins 73, 690–704 (2008).
42
Y.-W. Lin et al., Introducing a 2-His-1-Glu nonheme iron center into myoglobin confers nitric oxide reductase activity. J. Am. Chem. Soc. 132, 9970–9972 (2010).
43
A. Giuffrè et al., The heme-copper oxidases of Thermus thermophilus catalyze the reduction of nitric oxide: Evolutionary implications. Proc. Natl. Acad. Sci. U.S.A. 96, 14718 (1999).
44
E. Forte et al., The cytochrome cbb3 from Pseudomonas stutzeri displays nitric oxide reductase activity. Eur. J. Biochem. 268, 6486–6491 (2001).
45
V. B. Borisov, R. B. Gennis, J. Hemp, M. I. Verkhovsky, The cytochrome bd respiratory oxygen reductases. Biochim. Biophys. Acta 1807, 1398–1413 (2011).
46
H. B. Gray, J. R. Winkler, Hole hopping through tyrosine/tryptophan chains protects proteins from oxidative damage. Proc. Natl. Acad. Sci. U.S.A. 112, 10920–10925 (2015).
47
J. Chen, M. Strous, Denitrification and aerobic respiration, hybrid electron transport chains and co-evolution. Biochim. Biophys. Acta 1827, 136–144 (2013).
48
H.-Y. Chang, J. Hemp, Y. Chen, J. A. Fee, R. B. Gennis, The cytochrome ba3 oxygen reductase from Thermus thermophilus uses a single input channel for proton delivery to the active site and for proton pumping. Proc. Natl. Acad. Sci. U.S.A. 106, 16169–16173 (2009).
49
M. R. A. Blomberg, P. E. M. Siegbahn, Why is the reduction of NO in cytochrome c dependent nitric oxide reductase (cNOR) not electrogenic? Biochim. Biophys. Acta 1827, 826–833 (2013).
50
M. Wikström, V. Sharma, Proton pumping by cytochrome c oxidase—A 40 year anniversary. Biochim. Biophys. Acta 1859, 692–698 (2018).
51
P. Y. Camejo, B. O. Oyserman, K. D. McMahon, D. R. Noguera, Integrated omic analyses provide evidence that a “Candidatus Accumulibacter phosphatis” strain performs denitrification under microaerobic conditions. mSystems 4, e00193-18 (2019).
52
Y.-M. Huang, D. Straub, N. Blackwell, A. Kappler, S. Kleindienst, Meta-omics reveal Gallionellaceae and Rhodanobacter species as interdependent key players for Fe(II) oxidation and nitrate reduction in the autotrophic enrichment culture KS. Appl. Environ. Microbiol. 87, e00496-21 (2021).
53
C. Brochier-Armanet, E. Talla, S. Gribaldo, The multiple evolutionary histories of dioxygen reductases: Implications for the origin and evolution of aerobic respiration. Mol. Biol. Evol. 26, 285–297 (2009).
54
C. Castelle et al., A new iron-oxidizing/O2-reducing supercomplex spanning both inner and outer membranes, isolated from the extreme acidophile acidithiobacillus ferrooxidans. J. Biol. Chem. 283, 25803–25811 (2008).
55
J. Castresana, M. Saraste, Evolution of energetic metabolism: The respiration-early hypothesis. Trends Biochem. Sci. 20, 443–448 (1995).
56
A. Loullis, E. Pinakoulaki, Probing the nitrite and nitric oxide reductase activity of cbb3 oxidase: Resonance Raman detection of a six-coordinate ferrous heme–nitrosyl species in the binuclear b3/CuB center. Chem. Commun. 51, 17398–17401 (2015).
57
F. J. Stewart, O. Ulloa, E. F. DeLong, Microbial metatranscriptomics in a permanent marine oxygen minimum zone. Environ. Microbiol. 14, 23–40 (2012).
58
T. Kalvelage et al., Aerobic microbial respiration in oceanic oxygen minimum zones. PLoS One 10, e0133526 (2015).
59
E. J. Zakem, A. Mahadevan, J. M. Lauderdale, M. J. Follows, Stable aerobic and anaerobic coexistence in anoxic marine zones. ISME J. 14, 288–301 (2020).
60
V. R. I. Kaila, M. Wikström, Architecture of bacterial respiratory chains. Nat. Rev. Microbiol. 19, 319–330 (2021).
61
P. G. Falkowski, T. Fenchel, E. F. Delong, The microbial engines that drive earth’s biogeochemical cycles. Science 320, 1034–1039 (2008).
62
J. J. R. Fraústo da Silva, R. J. P. Williams, The Biological Chemistry of the Elements (Clarendon-Press, 1991).
63
R. J. P. Williams, J. J. R. F. da Silva, The Natural Selection of the Chemical Elements (Clarendon-Press, 1996).
64
R. J. P. Williams, J. J. R. Fraústo da silva, Evolution was chemically constrained. J. Theor. Biol. 220, 323–343 (2003).
65
M. M. Pereira et al., The caa3 terminal oxidase of the thermohalophilic bacterium Rhodothermus marinus: A HiPIP:oxygen oxidoreductase lacking the key glutamate of the D-channel. Biochim. Biophys. Acta 1413, 1–13 (1999, 1413,).
66
E. A. Berry, B. L. Trumpower, Simultaneous determination of hemes a, b, and c from pyridine hemochrome spectra. Anal. Biochem. 161, 1–15 (1987).
67
R. Murali, Curated Hidden Markov Models for the heme-copper oxidoreductase superfamily. HCO HMMs. https://github.com/ranjani-m/HCO. Deposited 7 June 2022.
Information & Authors
Information
Published in
Classifications
Copyright
Copyright © 2024 the Author(s). Published by PNAS. This article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).
Data, Materials, and Software Availability
All the protein accession numbers used for generation of trees in this study, as well as associated phylogenetic trees and multiple sequence alignments are included in the supporting information. The Hidden Markov Models (HMMs) used for identification of HCO sequences can be found at https://github.com/ranjani-m/HCO (67).
Submission history
Received: September 27, 2023
Accepted: April 24, 2024
Published online: June 20, 2024
Published in issue: June 25, 2024
Keywords
Acknowledgments
We thank the NIH for funding support (grant# U12AB123456 to PI: R.B.G.). This research was also supported by funding from the Agouron Institute (W.W.F. and J.H.) and by the Community Science Project 507064 (PI: R.H.) under the Joint Genome Institute (https://ror.org/04xm1d337), which is a Department of Energy (DOE) Office of Science User Facility. Resources were also used at Office of Biological and Environmental Research of the US Department of Energy Atmospheric System Research Program Interagency Agreement grant DE-AC02-05CH11231 (JGI). Resources were used at Office of Biological and Environmental Research of the United States Department of Energy Atmospheric System Research Program Interagency Agreement grant DE-AC05-76RL01830 Environmental Molecular Sciences Laboratory (EMSL). We thank Sylvia Choi for providing pure ba3 oxygen reductase from T. thermophilus to use as a control for oxygen reductase assays and for heme extraction, Paige Sheridan for providing purified qNOR from Persephonella marina, Lici Schurig-Briccio for guidance in performing NOR assays with the Clark Electrode, and Peter Yau at the University of Illinois’ Mass spectrometric facility for protein identification. We thank Alon Philosof and Connor Skennerton for valuable discussions on bioinformatics analysis. Finally, we would like to thank our reviewers for their valuable comments that led to the improvement of this manuscript.
Author contributions
R.M., W.W.F., R.B.G., and J.H. designed research; R.M., L.A.P., R.A.S., L.M.W., M.M.L., R.H., W.W.F., and J.H. performed research; R.M., R.A.S., L.M.W., W.W.F., R.B.G., and J.H. contributed new reagents/analytic tools; R.M., L.A.P., U.F.L., W.W.F., R.B.G., and J.H. analyzed data; and R.M., R.A.S., L.M.W., M.M.L., R.H., W.W.F., R.B.G., and J.H. wrote the paper.
Competing interests
The authors declare no competing interest.
Notes
This article is a PNAS Direct Submission. J.J.M. is a guest editor invited by the Editorial Board.
Authors
Metrics & Citations
Metrics
Altmetrics
Citations
Cite this article
Diversity and evolution of nitric oxide reduction in bacteria and archaea, Proc. Natl. Acad. Sci. U.S.A.
121 (26) e2316422121,
https://doi.org/10.1073/pnas.2316422121
(2024).
Copied!
Copying failed.
Export the article citation data by selecting a format from the list below and clicking Export.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFLogin options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.