Mechanism for N2O Generation in Bacterial Nitric Oxide Reductase: A Quantum Chemical Study

Bacterial nitric oxide reductase (NorBC) models employing click chemistry

Bacterial NO Reductase (NorBC or cNOR) is a membrane-bound enzyme found in denitrifying bacteria that catalyzes the two-electron reduction of NO to N₂O and water. The mechanism by which NorBC operates is highly debated, due to the fact that this enzyme is difficult to work with, and no intermediates of the NO reduction reaction could have been identified so far. The unique active site of NorBC consists of a heme b₃/non-heme Fe_B diiron center. Synthetic model complexes provide the opportunity to obtain insight into possible mechanistic alternatives for this enzyme. In this paper, we present three new synthetic model systems for NorBC, consisting of a tetraphenylporphyrin-derivative clicked to modified BMPA-based ligands (BMPA = bis(methylpyridyl)amine) that model the non-heme site in the enzyme. These complexes have been characterized by EPR, IR and UV-Vis spectroscopy. The reactivity with NO was then investigated, and it was found that the complex with the BMPA-carboxylate ligand as the non-heme component has a very low affinity for NO at the non-heme iron site. If the carboxylate functional group is replaced with a phenolate or pyridine group, reactivity is restored and formation of a diiron dinitrosyl complex was observed. Upon one-electron reduction of the nitrosylated complexes, following the semireduced pathway for NO reduction, formation of dinitrosyl iron complexes (DNICs) was observed in all three cases, but no N₂O could be detected.

Spectroscopic Characterization of a Diferric Mycobacterial Hemerythrin-Like Protein with Unprecedented Reactivity toward Nitric Oxide

Hemerythrin-like proteins (HLPs) are broadly distributed across taxonomic groups and appear to play highly diverse functional roles in prokaryotes. Mycobacterial HLPs contribute to the survival of these pathogenic bacteria in mammalian macrophages, but their modes of action remain unclear. A recent crystallographic characterization of Mycobacterium kansasii HLP (Mka-HLP) revealed the unexpected presence of a tyrosine sidechain (Tyr54) near the coordination sphere of one of the two iron centers. Here, we show that Tyr54 is a true ligand to the Fe2(III) ion which, in conjunction with the presence of a μ-oxo group bridging the two iron(III), brings unique reactivity toward nitric oxide (NO). Monitoring the titration of Mka-HLP with NO by Fourier-transform infrared and electron paramagnetic resonance spectroscopies shows that both diferric and diferrous forms of Mka-HLP accumulate an uncoupled high-spin and low-spin {FeNO}⁷ pair. We assign the reactivity of the diferric protein to an initial radical reaction between NO and the μ-oxo bridge to form nitrite and a mixed-valent diiron center that can react further with NO. Amperometric measurements of NO consumption by Mka-HLP confirm that this reactivity can proceed at low micromolar concentrations of NO, before additional NO consumption, supporting a NO scavenging role for mycobacterial HLPs.

NO Coupling at Copper to cis-Hyponitrite: N2O Formation via Protonation and H-Atom Transfer

Copper nitrite reductases (CuNIRs) convert NO2- to NO as well as NO to N2O under high NO flux at a mononuclear type 2 Cu center. While model complexes illustrate N-N coupling from NO that results in symmetric trans-hyponitrite [CuII]-ONNO-[CuII] complexes, we report NO assembly at a single Cu site in the presence of an external reductant Cp*2M (M = Co, Fe) to give the first copper cis-hyponitrites [Cp*2M]{[CuII](κ2-O2N2)[CuI]}. Importantly, the κ1-N-bound [CuI] fragment may be easily removed by the addition of mild Lewis bases such as CNAr or pyridine to form the spectroscopically similar anion {[CuII](κ2-O2N2)}-. The addition of electrophiles such as H+ to these anionic copper(II) cis-hyponitrites leads to N2O generation with the formation of the dicopper(II)-bis-μ-hydroxide [CuII]2(μ-OH)2. One-electron oxidation of the {[CuII](κ2-O2N2)}- core turns on H-atom transfer reactivity, enabling the oxidation of 9,10-dihydroanthracene to anthracene with concomitant formation of N2O and [CuII]2(μ-OH)2. These studies illustrate both the reductive coupling of NO at a single copper center and a way to harness the strong oxidizing power of nitric oxide via the neutral cis-hyponitrite [Cu](κ2-O2N2).

Reaction Mechanism of Human PAICS Elucidated by Quantum Chemical Calculations

Human PAICS is a bifunctional enzyme that is involved in the de novo purine biosynthesis, catalyzing the conversion of aminoimidazole ribonucleotide (AIR) into N-succinylcarboxamide-5-aminoimidazole ribonucleotide (SAICAR). It comprises two distinct active sites, AIR carboxylase (AIRc) where the AIR is initially converted to carboxyaminoimidazole ribonucleotide (CAIR) by reaction with CO₂ and SAICAR synthetase (SAICARs) in which CAIR then reacts with an aspartate to form SAICAR, in an ATP-dependent reaction. Human PAICS is a promising target for the treatment of various types of cancer, and it is therefore of high interest to develop a detailed understanding of its reaction mechanism. In the present work, density functional theory calculations are employed to investigate the PAICS reaction mechanism. Starting from the available crystal structures, two large models of the AIRc and SAICARs active sites are built and different mechanistic proposals for the carboxylation and phosphorylation–condensation mechanisms are examined. For the carboxylation reaction, it is demonstrated that it takes place in a two-step mechanism, involving a C–C bond formation followed by a deprotonation of the formed tetrahedral intermediate (known as isoCAIR) assisted by an active site histidine residue. For the phosphorylation–condensation reaction, it is shown that the phosphorylation of CAIR takes place before the condensation reaction with the aspartate. It is further demonstrated that the three active site magnesium ions are involved in binding the substrates and stabilizing the transition states and intermediates of the reaction. The calculated barriers are in good agreement with available experimental data.

The importance of exact exchange-A methodological investigation of NO reduction in heme-copper oxidases

Significant improvements of the density functional theory (DFT) methodology during the past few decades have made DFT calculations a powerful tool in studies of enzymatic reaction mechanisms. For metalloenzymes, however, there are still concerns about the reliability in the DFT-results. Therefore, a systematic study is performed where the fraction of exact exchange in a hybrid DFT functional is used as a parameter. By varying this parameter, a set of different but related functionals are obtained. The various functionals are applied to one of the reactions occurring in the enzyme family heme-copper oxidases, the reduction of nitric oxide (NO) to nitrous oxide (N₂O) and water. The results show that, even though certain parts of the calculated energetics exhibit large variations, the qualitative pictures of the reaction mechanisms are quite stable. Furthermore, it is found that the functional with 15% exact exchange (B3LYP^*) gives the best agreement with experimental data for the particular reactions studied. An important aspect of the procedure used is that the computational results are carefully combined with a few more general experimental data to obtain a complete description of the entire catalytic cycle of the reactions studied.

Activation of O2 and NO in heme-copper oxidases - mechanistic insights from computational modelling

Heme-copper oxidases are transmembrane enzymes involved in aerobic and anaerobic respiration. The largest subgroup contains the cytochrome c oxidases (CcO), which reduce molecular oxygen to water. A significant part of the free energy released in this exergonic process is conserved as an electrochemical gradient across the membrane, via two processes, electrogenic chemistry and proton pumping. A deviant subgroup is the cytochrome c dependent NO reductases (cNOR), which reduce nitric oxide to nitrous oxide and water. This is also an exergonic reaction, but in this case none of the released free energy is conserved. Computational studies applying hybrid density functional theory to cluster models of the bimetallic active sites in the heme-copper oxidases are reviewed. To obtain a reliable description of the reaction mechanisms, energy profiles of the entire catalytic cycles, including the reduction steps have to be constructed. This requires a careful combination of computational results with certain experimental data. Computational studies have elucidated mechanistic details of the chemical parts of the reactions, involving cleavage and formation of covalent bonds, which have not been obtainable from pure experimental investigations. Important insights regarding the mechanisms of energy conservation have also been gained. The computational studies show that the reduction potentials of the active site cofactors in the CcOs are large enough to afford electrogenic chemistry and proton pumping, i.e. efficient energy conservation. These results solve a conflict between different types of experimental data. A mechanism for the proton pumping, involving a specific and crucial role for the active site tyrosine, conserved in all CcOs, is suggested. For the cNORs, the calculations show that the low reduction potentials of the active site cofactors are optimized for fast elimination of the toxic NO molecules. At the same time, the low reduction potentials lead to endergonic reduction steps with high barriers. To prevent even higher barriers, which would lead to a too slow reaction, when the electrochemical gradient across the membrane is present, the chemistry must occur in a non-electrogenic manner. This explains why there is no energy conservation in cNOR.

Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling

Heme-copper oxidases (HCO), nitric oxide reductases (NOR), and sulfite reductases (SiR) catalyze the multi-electron and multi-proton reductions of O2, NO, and SO32-, respectively. Each of these reactions is important to drive cellular energy production through respiratory metabolism and HCO, NOR, and SiR evolved to contain heteronuclear active sites containing heme/copper, heme/nonheme iron, and heme-[4Fe-4S] centers, respectively. The complexity of the structures and reactions of these native enzymes, along with their large sizes and/or membrane associations, make it challenging to fully understand the crucial structural features responsible for the catalytic properties of these active sites. In this review, we summarize progress that has been made to better understand these heteronuclear metalloenzymes at the molecular level though study of the native enzymes along with insights gained from biomimetic models comprising either small molecules or proteins. Further understanding the reaction selectivity of these enzymes is discussed through comparisons of their similar heteronuclear active sites, and we offer outlook for further investigations.

Quantification of Ni-N-O Bond Angles and NO Activation by X-ray Emission Spectroscopy

A series of β-diketiminate Ni-NO complexes with a range of NO binding modes and oxidation states were studied by X-ray emission spectroscopy (XES). The results demonstrate that XES can directly probe and distinguish end-on vs side-on NO coordination modes as well as one-electron NO reduction. Density functional theory (DFT) calculations show that the transition from the NO 2s2s σ* orbital has higher intensity for end-on NO coordination than for side-on NO coordination, whereas the 2s2s σ orbital has lower intensity. XES calculations in which the Ni-N-O bond angle was fixed over the range from 80° to 176° suggest that differences in NO coordination angles of ∼10° could be experimentally distinguished. Calculations of Cu nitrite reductase (NiR) demonstrate the utility of XES for characterizing NO intermediates in metalloenzymes. This work shows the capability of XES to distinguish NO coordination modes and oxidation states at Ni and highlights applications in quantifying small molecule activation in enzymes.

Role of the Two Metals in the Active Sites of Heme Copper Oxidases-A Study of NO Reduction in cbb3 Cytochrome c Oxidase

The superfamily of heme copper oxidases reduces molecular oxygen or nitric oxide, and the active sites comprise a high-spin heme group (a₃ or b₃) and a non-heme metal (Cu_B or Fe_B). The cbb₃ C family of cytochrome c oxidases, with the high-spin heme b₃ and Cu_B in the active site, is a subfamily of the heme copper oxidases that can reduce both molecular oxygen, which is the main substrate, and nitric oxide. The mechanism for NO reduction in cbb₃ oxidase is studied here using hybrid density functional theory and compared to other cytochrome c oxidases (A and B families), with a high-spin heme a₃ and Cu_B in the active site, and to cytochrome c dependent NO reductase, with a high-spin heme b₃ and a non-heme Fe_B in the active site. It is found that the reaction mechanism and the detailed reaction energetics of the cbb₃ oxidases are not similar to those of cytochrome c dependent NO reductase, which has the same type of high-spin heme group but a different non-heme metal. This is in contrast to earlier expectations. Instead, the NO reduction mechanism in cbb₃ oxidases is very similar to that in the other cytochrome c oxidases, with the same non-heme metal, Cu_B, and is independent of the type of high-spin heme group. The conclusion is that the type of non-heme metal (Cu_B or Fe_B) in the active site of the heme copper oxidases is more important for the reaction mechanisms than the type of high-spin heme, at least for the NO reduction reaction. The reason is that the proton-coupled reduction potentials of the active site cofactors determine the energetics for the NO reduction reaction, and they depend to a larger extent on the non-heme metal. Observed differences in NO reduction reactivity among the various cytochrome c oxidases may be explained by differences outside the BNC, affecting the rate of proton transfer, rather than in the BNC itself.

Biological and Bioinspired Inorganic N-N Bond-Forming Reactions

The metallobiochemistry underlying the formation of the inorganic N-N-bond-containing molecules nitrous oxide (N2O), dinitrogen (N2), and hydrazine (N2H4) is essential to the lifestyles of diverse organisms. Similar reactions hold promise as means to use N-based fuels as alternative carbon-free energy sources. This review discusses research efforts to understand the mechanisms underlying biological N-N bond formation in primary metabolism and how the associated reactions are tied to energy transduction and organismal survival. These efforts comprise studies of both natural and engineered metalloenzymes as well as synthetic model complexes.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

Insights into the mechanism of nitric oxide reductase from a FeB -depleted variant

A key step of denitrification, the reduction of toxic nitric oxide to nitrous oxide, is catalysed by cytochrome c-dependent NO reductase (cNOR). cNOR contains four redox-active cofactors: three hemes and a nonheme iron (Fe_B ). Heme b₃ and Fe_B constitute the active site, but the specific mechanism of NO-binding events and reduction is under debate. Here, we used a recently constructed, fully folded and hemylated cNOR variant that lacks Fe_B to investigate the role of Fe_B during catalysis. We show that in the Fe_B -less cNOR, binding of both NO and O₂ to heme b₃ still occurs but further reduction is impaired, although to a lesser degree for O₂ than for NO. Implications for the catalytic mechanisms of cNOR are discussed.

Nitrosyl Linkage Isomers: NO Coupling to N2O at a Mononuclear Site

Linkage isomers of reduced metal-nitrosyl complexes serve as key species in nitric oxide (NO) reduction at monometallic sites to produce nitrous oxide (N₂O), a potent greenhouse gas. While factors leading to extremely rare side-on nitrosyls are unclear, we describe a pair of nickel-nitrosyl linkage isomers through controlled tuning of noncovalent interactions between the nitrosyl ligands and differently encapsulated potassium cations. Furthermore, these reduced metal-nitrosyl species with N-centered spin density undergo radical coupling with free NO and provide a N-N coupled cis-hyponitrite intermediate whose protonation triggers the release of N₂O. This report outlines a stepwise molecular mechanism of NO reduction to form N₂O at a mononuclear metal site that provides insight into the related biological reduction of NO to N₂O.

Nitric Oxide Reductase Activity in Heme-Nonheme Binuclear Engineered Myoglobins through a One-Electron Reduction Cycle

Fe_BMbs are structural and functional models of native bacterial nitric oxide reductases (NORs) generated through engineering of myoglobin. These biosynthetic models replicate the heme-nonheme diiron site of NORs and allow substitutions of metal centers and heme cofactors. Here, we provide evidence for multiple NOR turnover in monoformyl-heme-containing Fe_BMb1 proteins loaded with Fe^II, Co^II, or Zn^II metal ions at the Fe_B site (Fe^II/Co^II/Zn^II-Fe_BMb1(MF-heme)). FTIR detection of the ν(NNO) band of N₂O at 2231 cm^-1 provides a direct quantitative measurement of the product in solution. A maximum number of turnover is observed with Fe^II-Fe_BMb1(MF-heme), but the NOR activity is retained when the Fe_B site is loaded with Zn^II. These data support the viability of a one-electron semireduced pathway for the reduction of NO at binuclear centers in reducing conditions.

Mechanisms for enzymatic reduction of nitric oxide to nitrous oxide - A comparison between nitric oxide reductase and cytochrome c oxidase

Cytochrome c oxidases (CcO) reduce O₂ to H₂O in the respiratory chain of mitochondria and many aerobic bacteria. In addition, some species of CcO can also reduce NO to N₂O and water while others cannot. Here, the mechanism for NO-reduction in CcO is investigated using quantum mechanical calculations. Comparison is made to the corresponding reaction in a "true" cytochrome c-dependent NO reductase (cNOR). The calculations show that in cNOR, where the reduction potentials are low, the toxic NO molecules are rapidly reduced, while the higher reduction potentials in CcO lead to a slower or even impossible reaction, consistent with experimental observations. In both enzymes the reaction is initiated by addition of two NO molecules to the reduced active site, forming a hyponitrite intermediate. In cNOR, N₂O can then be formed using only the active-site electrons. In contrast, in CcO, one proton-coupled reduction step most likely has to occur before N₂O can be formed, and furthermore, proton transfer is most likely rate-limiting. This can explain why different CcO species with the same heme a₃-Cu active site differ with respect to NO reduction efficiency, since they have a varying number and/or properties of proton channels. Finally, the calculations also indicate that a conserved active site valine plays a role in reducing the rate of NO reduction in CcO.

Surface-Enhanced Infrared Absorption Spectroscopy of Bacterial Nitric Oxide Reductase under Electrochemical Control Using a Vibrational Probe of Carbon Monoxide

Nitric oxide reductases (NORs) reduce nitric oxide to nitrous oxide in the denitrification pathway of the global nitrogen cycle. NORs contain four iron cofactors and the NO reduction occurs at the heme b₃/nonheme Fe_B binuclear active site. The determination of reduction potentials of the iron cofactors will help us elucidate the enzymatic reaction mechanism. However, previous reports on these potentials remain controversial. Herein, we performed electrochemical and surface-enhanced infrared absorption (SEIRA) spectroscopic measurements of Pseudomonas aeruginosa NOR immobilized on gold electrodes. Cyclic voltammograms exhibited two reduction peaks at -0.11 and -0.44 V vs SHE, and a SEIRA spectrum using a vibrational probe of CO showed a characteristic band at 1972 cm^-1 at -0.4 V vs SHE, which was assigned to νCO of heme b₃-CO. Our results suggest that the reduction of heme b₃ initiates the enzymatic NO reduction.

A Theoretical Study of the Benzoylformate Decarboxylase Reaction Mechanism

Density functional theory calculations are used to investigate the detailed reaction mechanism of benzoylformate decarboxylase, a thiamin diphosphate (ThDP)-dependent enzyme that catalyzes the nonoxidative decarboxylation of benzoylformate yielding benzaldehyde and carbon dioxide. A large model of the active site is constructed on the basis of the X-ray structure, and it is used to characterize the involved intermediates and transition states and evaluate their energies. There is generally good agreement between the calculations and available experimental data. The roles of the various active site residues are discussed and the results are compared to mutagenesis experiments. Importantly, the calculations identify off-cycle intermediate species of the ThDP cofactor that can have implications on the kinetics of the reaction.

The insertion of the non-heme FeB cofactor into nitric oxide reductase from P. denitrificans depends on NorQ and NorD accessory proteins

Bacterial NO reductases (NOR) catalyze the reduction of NO into N₂O, either as a step in denitrification or as a detoxification mechanism. cNOR from Paracoccus (P.) denitrificans is expressed from the norCBQDEF operon, but only the NorB and NorC proteins are found in the purified NOR complex. Here, we established a new purification method for the P. denitrificans cNOR via a His-tag using heterologous expression in E. coli. The His-tagged enzyme is both structurally and functionally very similar to non-tagged cNOR. We were also able to express and purify cNOR from the structural genes norCB only, in absence of the accessory genes norQDEF. The produced protein is a stable NorCB complex containing all hemes and it can bind gaseous ligands (CO) to heme b₃, but it is catalytically inactive. We show that this deficient cNOR lacks the non-heme iron cofactor Fe_B. Mutational analysis of the nor gene cluster revealed that it is the norQ and norD genes that are essential to form functional cNOR. NorQ belongs to the family of MoxR P-loop AAA+ ATPases, which are in general considered to facilitate enzyme activation processes often involving metal insertion. Our data indicates that NorQ and NorD work together in order to facilitate non-heme Fe insertion. This is noteworthy since in many cases Fe cofactor binding occurs spontaneously. We further suggest a model for NorQ/D-facilitated metal insertion into cNOR.

Heme redox potentials hold the key to reactivity differences between nitric oxide reductase and heme-copper oxidase

Despite high structural homology between NO reductases (NORs) and heme-copper oxidases (HCOs), factors governing their reaction specificity remain to be understood. Using a myoglobin-based model of NOR (Fe_BMb) and tuning its heme redox potentials (E°') to cover the native NOR range, through manipulating hydrogen bonding to the proximal histidine ligand and replacing heme b with monoformyl (MF-) or diformyl (DF-) hemes, we herein demonstrate that the E°' holds the key to reactivity differences between NOR and HCO. Detailed electrochemical, kinetic, and vibrational spectroscopic studies, in tandem with density functional theory calculations, demonstrate a strong influence of heme E°' on NO reduction. Decreasing E°' from +148 to -130 mV significantly impacts electronic properties of the NOR mimics, resulting in 180- and 633-fold enhancements in NO association and heme-nitrosyl decay rates, respectively. Our results indicate that NORs exhibit finely tuned E°' that maximizes their enzymatic efficiency and helps achieve a balance between opposite factors: fast NO binding and decay of dinitrosyl species facilitated by low E°' and fast electron transfer facilitated by high E°'. Only when E°' is optimally tuned in Fe_BMb(MF-heme) for NO binding, heme-nitrosyl decay, and electron transfer does the protein achieve multiple (>35) turnovers, previously not achieved by synthetic or enzyme-based NOR models. This also explains a long-standing question in bioenergetics of selective cross-reactivity in HCOs. Only HCOs with heme E°' in a similar range as NORs (between -59 and 200 mV) exhibit NOR reactivity. Thus, our work demonstrates efficient tuning of E°' in various metalloproteins for their optimal functionality.

Characterization of the quinol-dependent nitric oxide reductase from the pathogen Neisseria meningitidis, an electrogenic enzyme

Bacterial nitric oxide reductases (NORs) catalyse the reduction of NO to N₂O and H₂O. NORs are found either in denitrification chains, or in pathogens where their primary role is detoxification of NO produced by the immune defense of the host. Although NORs belong to the heme-copper oxidase superfamily, comprising proton-pumping O₂-reducing enzymes, the best studied NORs, cNORs (cytochrome c-dependent), are non-electrogenic. Here, we focus on another type of NOR, qNOR (quinol-dependent). Recombinant qNOR from Neisseria meningitidis, a human pathogen, purified from Escherichia coli, showed high catalytic activity and spectroscopic properties largely similar to cNORs. However, in contrast to cNOR, liposome-reconstituted qNOR showed respiratory control ratios above two, indicating that NO reduction by qNOR was electrogenic. Further, we determined a 4.5 Å crystal structure of the N. meningitidis qNOR, allowing exploration of a potential proton transfer pathway from the cytoplasm by mutagenesis. Most mutations had little effect on the activity, however the E-498 variants were largely inactive, while the corresponding substitution in cNOR was previously shown not to induce significant effects. We thus suggest that, contrary to cNOR, the N. meningitidis qNOR uses cytoplasmic protons for NO reduction. Our results allow possible routes for protons to be discussed.

Mechanism of NO–CO reaction over highly dispersed cuprous oxide on γ-alumina catalyst using a metal–support interfacial site in the presence of oxygen: similarities to and differences from biological systems

The NO–CO reaction mechanism over the Cu/γ-Al₂O₃ catalyst was elucidated using DFT and a cluster model.

The mechanism for oxygen reduction in cytochrome c dependent nitric oxide reductase (cNOR) as obtained from a combination of theoretical and experimental results

Bacterial NO-reductases (NOR) belong to the heme-copper oxidase (HCuO) superfamily, in which most members are O₂-reducing, proton-pumping enzymes. This study is one in a series aiming to elucidate the reaction mechanisms of the HCuOs, including the mechanisms for cellular energy conservation. One approach towards this goal is to compare the mechanisms for the different types of HCuOs, cytochrome c oxidase (CcO) and NOR, reducing the two substrates O₂ and NO. Specifically in this study, we describe the mechanism for oxygen reduction in cytochrome c dependent NOR (cNOR). Hybrid density functional calculations were performed on large cluster models of the cNOR binuclear active site. Our results are used, together with published experimental information, to construct a free energy profile for the entire catalytic cycle. Although the overall reaction is quite exergonic, we show that during the reduction of molecular oxygen in cNOR, two of the reduction steps are endergonic with high barriers for proton uptake, which is in contrast to oxygen reduction in CcO, where all reduction steps are exergonic. This difference between the two enzymes is suggested to be important for their differing capabilities for energy conservation. An additional result from this study is that at least three of the four reduction steps are initiated by proton transfer to the active site, which is in contrast to CcO, where electrons always arrive before the protons to the active site. The roles of the non-heme metal ion and the redox-active tyrosine in the active site are also discussed.

Modulation of protein function in membrane mimetics: Characterization of P. denitrificans cNOR in nanodiscs or liposomes

For detailed functional characterization, membrane proteins are usually studied in detergent. However, it is becoming clear that detergent micelles are often poor mimics of the lipid environment in which these proteins function. In this work we compared the catalytic properties of the membrane-embedded cytochrome c-dependent nitric oxide reductase (cNOR) from Paracoccus (P.) denitrificans in detergent, lipid/protein nanodiscs, and proteoliposomes. We used two different lipid mixtures, an extract of soybean lipids and a defined mix of synthetic lipids mimicking the original P. denitrificans membrane. We show that the catalytic activity of detergent-solubilized cNOR increased threefold upon reconstitution from detergent into proteoliposomes with the P. denitrificans lipid mixture, and above two-fold when soybean lipids were used. In contrast, there was only a small activity increase in nanodiscs. We further show that binding of the gaseous ligands CO and O₂ are affected differently by reconstitution. In proteoliposomes the turnover rates are affected much more than in nanodiscs, but CO-binding is more significantly accelerated in liposomes with soybean lipids, while O₂-binding is faster with the P. denitrificans lipid mix. We also investigated proton-coupled electron transfer during the reaction between fully reduced cNOR and O₂, and found that the pK_a of the internal proton donor was increased in proteoliposomes but not in nanodiscs. Taking our results together, the liposome-reconstituted enzyme shows significant differences to detergent-solubilized protein. Nanodiscs show much more subtle effects, presumably because of their much lower lipid to protein ratio. Which of these two membrane-mimetic systems best mimics the native membrane is discussed.

Recent Trends in Quantum Chemical Modeling of Enzymatic Reactions

The quantum chemical cluster approach is a powerful method for investigating enzymatic reactions. Over the past two decades, a large number of highly diverse systems have been studied and a great wealth of mechanistic insight has been developed using this technique. This Perspective reviews the current status of the methodology. The latest technical developments are highlighted, and challenges are discussed. Some recent applications are presented to illustrate the capabilities and progress of this approach, and likely future directions are outlined.

Nitric oxide coupling to generate N2O promoted by a single-heme system as examined by density functional theory

Bacteria utilize a heme/non-heme enzyme system to detoxify nitric oxide (NO) to N₂O. In order to probe the capacity of a single-heme system to mediate this NO-to-N₂O transformation, various scenarios for addition of electrons, protons, and a second NO molecule to a heme nitrosyl to generate N₂O were explored by density functional theory calculations. We describe, utilizing this single-heme system, several stepwise intermediates along pathways that enable the critical N-N bond formation step yielding the desired Fe-N₂O product. We also report a hitherto unreported directional second protonation that results in either productive N₂O formation with loss of water, or formation of a non-productive hyponitrous acid adduct Fe{HONNOH}.

Investigating the Proton Donor in the NO Reductase from Paracoccus denitrificans

Variant nomenclature: the variants were made in the NorB subunit if not indicated by the superscript c, which are variants in the NorC subunit (e.g. E122A = exchange of Glu-122 in NorB for an Ala, E71cD; exchange of Glu-71 in NorC for an Asp). Bacterial NO reductases (NORs) are integral membrane proteins from the heme-copper oxidase superfamily. Most heme-copper oxidases are proton-pumping enzymes that reduce O2 as the last step in the respiratory chain. With electrons from cytochrome c, NO reductase (cNOR) from Paracoccus (P.) denitrificans reduces NO to N2O via the following reaction: 2NO+2e-+2H+→N2O+H2O. Although this reaction is as exergonic as O2-reduction, cNOR does not contribute to the electrochemical gradient over the membrane. This means that cNOR does not pump protons and that the protons needed for the reaction are taken from the periplasmic side of the membrane (since the electrons are donated from this side). We previously showed that the P. denitrificans cNOR uses a single defined proton pathway with residues Glu-58 and Lys-54 from the NorC subunit at the entrance. Here we further strengthened the evidence in support of this pathway. Our further aim was to define the continuation of the pathway and the immediate proton donor for the active site. To this end, we investigated the region around the calcium-binding site and both propionates of heme b3 by site directed mutagenesis. Changing single amino acids in these areas often had severe effects on cNOR function, with many variants having a perturbed active site, making detailed analysis of proton transfer properties difficult. Our data does however indicate that the calcium ligation sphere and the region around the heme b3 propionates are important for proton transfer and presumably contain the proton donor. The possible evolutionary link between the area for the immediate donor in cNOR and the proton loading site (PLS) for pumped protons in oxygen-reducing heme-copper oxidases is discussed.

Effect of Outer-Sphere Side Chain Substitutions on the Fate of the trans Iron-Nitrosyl Dimer in Heme/Nonheme Engineered Myoglobins (Fe(B)Mbs): Insights into the Mechanism of Denitrifying NO Reductases

Denitrifying NO reductases are transmembrane protein complexes that utilize a heme/nonheme diiron center at their active sites to reduce two NO molecules to the innocuous gas N2O. Fe(B)Mb proteins, with their nonheme iron sites engineered into the heme distal pocket of sperm whale myoglobin, are attractive models for studying the molecular details of the NO reduction reaction. Spectroscopic and structural studies of Fe(B)Mb constructs have confirmed that they reproduce the metal coordination spheres observed at the active site of the cytochrome c-dependent NO reductase from Pseudomonas aeruginosa. Exposure of Fe(B)Mb to excess NO, as examined by analytical and spectroscopic techniques, results primarily in the formation of a five-coordinate heme-nitrosyl complex without N2O production. However, substitution of the outer-sphere residue Ile107 with a glutamic acid (i.e., I107E) decreases the formation rate of the five-coordinate heme-nitrosyl complex and allows for the substoichiometric production of N2O. Here, we aim to better characterize the formation of the five-coordinate heme-nitrosyl complex and to explain why the level of N2O production increases with the I107E substitution. We follow the formation of the five-coordinate heme-nitrosyl inhibitory complex through the sequential exposure of Fe(B)Mb to different NO isotopomers using rapid-freeze-quench resonance Raman spectroscopy. The data show that the complex is formed by the displacement of the proximal histidine by a new NO molecule after the weakening of the Fe(II)-His bond in the intermediate six-coordinate low-spin (6cLS) heme-nitrosyl complex. These results lead us to explore diatomic migration within the scaffold of myoglobin and whether substitutions at residue 107 can be sufficient to control access to the proximal heme cavities. Results on a new Fe(B)Mb construct with an I107F substitution (Fe(B)Mb3) show an increased rate for the formation of the five-coordinate low-spin heme-nitrosyl complex without N2O production. Taken together, our results suggest that production of N2O from the [6cLS heme {FeNO}(7)/{Fe(B)NO}(7)] trans iron-nitrosyl dimer intermediate requires a proton transfer event facilitated by an outer-sphere residue such as E107 in Fe(B)Mb2 and E280 in P. aeruginosa cNOR.

Mechanism of Oxygen Reduction in Cytochrome c Oxidase and the Role of the Active Site Tyrosine

Cytochrome c oxidase, the terminal enzyme in the respiratory chain, reduces molecular oxygen to water and stores the released energy through electrogenic chemistry and proton pumping across the membrane. Apart from the heme-copper binuclear center, there is a conserved tyrosine residue in the active site (BNC). The tyrosine delivers both an electron and a proton during the O-O bond cleavage step, forming a tyrosyl radical. The catalytic cycle then occurs in four reduction steps, each taking up one proton for the chemistry (water formation) and one proton to be pumped. It is here suggested that in three of the reduction steps the chemical proton enters the center of the BNC, leaving the tyrosine unprotonated with radical character. The reproprotonation of the tyrosine occurs first in the final reduction step before binding the next oxygen molecule. It is also suggested that this reduction mechanism and the presence of the tyrosine are essential for the proton pumping. Density functional theory calculations on large cluster models of the active site show that only the intermediates with the proton in the center of the BNC and with an unprotonated tyrosyl radical have a high electron affinity of similar size as the electron donor, which is essential for the ability to take up two protons per electron and thus for the proton pumping. This type of reduction mechanism is also the only one that gives a free energy profile in accordance with experimental observations for the amount of proton pumping in the working enzyme.

A bridged di-iron porphyrin hyponitrite complex as a model for biological N2O production from hyponitrite

Heme-hyponitrites are intermediates that form at the bimetallic active sites of bacterial nitric oxide reductases. To probe a possible effect of the Fe-Fe distance on hyponitrite stability, we prepared a bridged bis-porphyrin Fe-hyponitrite compound, namely [(OEP-CH2)Fe]2(μ2,η(1),η(1)-ONNO). Its υNO of 992 cm(-1) (υ15NO of 976 cm(-1)) is close to the υNO of 983 cm(-1) reported previously by us for the crystallographically characterized [(OEP)Fe]2(μ2,η(1),η(1)-ONNO) compound. The bridged bis-porphyrin Fe-hyponitrite complex is unstable with respect to N2O production, supporting the role of the bis-Fe porphyrin system in hyponitrite conversion to N2O. The preparation and crystallographic determination of the bridging sulfato derivative is also reported.

Theoretical study of the reaction mechanism of phenolic acid decarboxylase

The cofactor-free phenolic acid decarboxylases (PADs) catalyze the non-oxidative decarboxylation of phenolic acids to their corresponding p-vinyl derivatives. Phenolic acids are toxic to some organisms, and a number of them have evolved the ability to transform these compounds, including PAD-catalyzed reactions. Since the vinyl derivative products can be used as polymer precursors and are also of interest in the food-processing industry, PADs might have potential applications as biocatalysts. We have investigated the detailed reaction mechanism of PAD from Bacillus subtilis using quantum chemical methodology. A number of different mechanistic scenarios have been considered and evaluated on the basis of their energy profiles. The calculations support a mechanism in which a quinone methide intermediate is formed by protonation of the substrate double bond, followed by C-C bond cleavage. A different substrate orientation in the active site is suggested compared to the literature proposal. This suggestion is analogous to other enzymes with p-hydroxylated aromatic compounds as substrates, such as hydroxycinnamoyl-CoA hydratase-lyase and vanillyl alcohol oxidase. Furthermore, on the basis of the calculations, a different active site residue compared to previous proposals is suggested to act as the general acid in the reaction. The mechanism put forward here is consistent with the available mutagenesis experiments and the calculated energy barrier is in agreement with measured rate constants. The detailed mechanistic understanding developed here might be extended to other members of the family of PAD-type enzymes. It could also be useful to rationalize the recently developed alternative promiscuous reactivities of these enzymes.

Recent advances in biosynthetic modeling of nitric oxide reductases and insights gained from nuclear resonance vibrational and other spectroscopic studies

This Forum Article focuses on recent advances in structural and spectroscopic studies of biosynthetic models of nitric oxide reductases (NORs). NORs are complex metalloenzymes found in the denitrification pathway of Earth’s nitrogen cycle where they catalyze the proton-dependent two-electron reduction of nitric oxide (NO) to nitrous oxide (N₂O). While much progress has been made in biochemical and biophysical studies of native NORs and their variants, a clear mechanistic understanding of this important metalloenzyme related to its function is still elusive. We report herein UV–vis and nuclear resonance vibrational spectroscopy (NRVS) studies of mononitrosylated intermediates of the NOR reaction of a biosynthetic model. The ability to selectively substitute metals at either heme or nonheme metal sites allows the introduction of independent ⁵⁷Fe probe atoms at either site, as well as allowing the preparation of analogues of stable reaction intermediates by replacing either metal with a redox inactive metal. Together with previous structural and spectroscopic results, we summarize insights gained from studying these biosynthetic models toward understanding structural features responsible for the NOR activity and its mechanism. The outlook on NOR modeling is also discussed, with an emphasis on the design of models capable of catalytic turnovers designed based on close mimics of the secondary coordination sphere of native NORs.

New insights into nitric oxide reductases (NORs) are obtained. Using nuclear resonance vibrational spectroscopy, we probe both iron atoms in mononitrosylated intermediates of the NOR reaction in a biosynthetic protein model that reveal new insights into the structural and electronic features responsible for the NOR activity and its likely mechanism.

Proton-Assisted Mechanism of NO Reduction on a Dinuclear Ruthenium Complex

Density-functional-theory (DFT) calculations are performed for the proposal of a plausible mechanism on the reduction of NO to N2O by a dinuclear ruthenium complex, reported by Arikawa and co-workers [J. Am. Chem. Soc. 2007, 129, 14160]. On the basis of the experimental fact that the reduction proceeds under strongly acidic conditions, the role of protons in the mechanistic pathways is investigated with model complexes, where one or two NO ligands are protonated. The reaction mechanism of the NO reduction is partitioned into three steps: reorientation of N2O2 (cis-NO dimer), O-N bond cleavage, and N2O elimination. A key finding is that the protonation of the NO ligand(s) significantly reduces the activation barrier in the rate-determining reorientation step. The activation energy of 43.1 kcal/mol calculated for the proton-free model is reduced to 30.2 and 17.6 kcal/mol for the mono- and diprotonated models, respectively. The protonation induces the electron transfer from the Ru(II)Ru(II) core to the O═N-N═O moiety to give a Ru(III)Ru(III) core and a hyponitrite (O-N═N-O)(2-) species. The formation of the hyponitrite species provides an alternative pathway for the N2O2 reorientation, resulting in the lower activation energies in the presence of proton(s). The protonation also has a marginal effect on the O-N bond cleavage and the N2O elimination steps. Our calculations reveal a remarkable role of protons in the NO reduction via N2O formation and provide new insights into the mechanism of NO reduction catalyzed by metalloenzymes such as nitric oxide reductase (NOR) that contains a diiron active site.

Protonation of the binuclear active site in cytochrome c oxidase decreases the reduction potential of CuB

One of the remaining mysteries regarding the respiratory enzyme cytochrome c oxidase is how proton pumping can occur in all reduction steps in spite of the low reduction potentials observed in equilibrium titration experiments for two of the active site cofactors, CuB(II) and Fea3(III). It has been speculated that, at least the copper cofactor can acquire two different states, one metastable activated state occurring during enzyme turnover, and one relaxed state with lower energy, reached only when the supply of electrons stops. The activated state should have a transiently increased CuB(II) reduction potential, allowing proton pumping. The relaxed state should have a lower reduction potential, as measured in the titration experiments. However, the structures of these two states are not known. Quantum mechanical calculations show that the proton coupled reduction potential for CuB is inherently high in the active site as it appears after reaction with oxygen, which explains the observed proton pumping. It is suggested here that, when the flow of electrons ceases, a relaxed resting state is formed by the uptake of one extra proton, on top of the charge compensating protons delivered in each reduction step. The extra proton in the active site decreases the proton coupled reduction potential for CuB by almost half a volt, leading to agreement with titration experiments. Furthermore, the structure for the resting state with an extra proton is found to have a hydroxo-bridge between CuB(II) and Fea3(III), yielding a magnetic coupling that can explain the experimentally observed EPR silence.

Bacterial nitric oxide reductase: a mechanism revisited by an ONIOM (DFT:MM) study

Bacterial nitric oxide reductase (cNOR) is an important binuclear iron enzyme responsible for the reduction of nitric oxide to nitrous oxide in the catalytic cycle of bacterial respiration. The reaction mechanism of cNOR as well as the key reactive intermediates of the reaction are still under debate. Here, we report a computational study based on ONIOM (DFT:MM) calculations aimed at investigating the reaction mechanism of cNOR. The results suggest that the reaction proceeds via the mono-nitrosyl mechanism which starts off by the binding of an NO molecule to the heme b3 center, N-N hyponitrite bond formation as a result of the reaction with a second NO molecule was found to proceed with an exothermic energy barrier to yield a hyponitrite adduct forming an open (incomplete) ring conformation with the non-heme FeB center (O-N-N-O-FeB). N-O bond cleavage to yield N2O was shown to be the rate-limiting step with an activation barrier of 22.6 kcal mol(-1). The dinitrosyl (trans) mechanism, previously proposed by several studies, was also examined and found unfavorable due to high activation barriers of the resulting intermediates.

How cytochrome c oxidase can pump four protons per oxygen molecule at high electrochemical gradient

Experiments have shown that the A-family cytochrome c oxidases pump four protons per oxygen molecule, also at a high electrochemical gradient. This has been considered a puzzle, since two of the reduction potentials involved, Cu(II) and Fe(III), were estimated from experiments to be too low to afford proton pumping at a high gradient. The present quantum mechanical study (using hybrid density functional theory) suggests a solution to this puzzle. First, the calculations show that the charge compensated Cu(II) potential for CuB is actually much higher than estimated from experiment, of the same order as the reduction potentials for the tyrosyl radical and the ferryl group, which are also involved in the catalytic cycle. The reason for the discrepancy between theory and experiment is the very large uncertainty in the experimental observations used to estimate the equilibrium potentials, mainly caused by the lack of methods for direct determination of reduced CuB. Second, the calculations show that a high energy metastable state, labeled EH, is involved during catalytic turnover. The EH state mixes the low reduction potential of Fe(III) in heme a3 with another, higher potential, here suggested to be that of the tyrosyl radical, resulting in enough exergonicity to allow proton pumping at a high gradient. In contrast, the corresponding metastable oxidized state, OH, is not significantly higher in energy than the resting state, O. Finally, to secure the involvement of the high energy EH state it is suggested that only one proton is taken up via the K-channel during catalytic turnover.