Next Article in Journal
The Effect of Storage-Induced Changes in Ammonia Borane on Hydrogen Release during Its Low-Temperature Thermolysis
Next Article in Special Issue
Preliminary Characterization of a Ni2+-Activated and Mycothiol-Dependent Glyoxalase I Enzyme from Streptomyces coelicolor
Previous Article in Journal
Effect of the (Nd,Dy)-Double Doping on the Structural Properties of Ceria
Previous Article in Special Issue
pH Dependent Reversible Formation of a Binuclear Ni2 Metal-Center within a Peptide Scaffold
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Inorganics
Volume 7
Issue 8
10.3390/inorganics7080095
inorganics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Theoretical Studies of Nickel-Dependent Enzymes

by
Per E. M. Siegbahn
1,*,
Shi-Lu Chen
2 and
Rong-Zhen Liao
3
1
Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
2
School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China
3
School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Inorganics 2019, 7(8), 95; https://doi.org/10.3390/inorganics7080095
Submission received: 10 June 2019 / Revised: 19 July 2019 / Accepted: 23 July 2019 / Published: 29 July 2019
(This article belongs to the Special Issue Bioinorganic Chemistry of Nickel)
Browse Figures
Review Reports Versions Notes

Abstract

:
The advancements of quantum chemical methods and computer power allow detailed mechanistic investigations of metalloenzymes. In particular, both quantum chemical cluster and combined QM/MM approaches have been used, which have been proven to successfully complement experimental studies. This review starts with a brief introduction of nickel-dependent enzymes and then summarizes theoretical studies on the reaction mechanisms of these enzymes, including NiFe hydrogenase, methyl-coenzyme M reductase, nickel CO dehydrogenase, acetyl CoA synthase, acireductone dioxygenase, quercetin 2,4-dioxygenase, urease, lactate racemase, and superoxide dismutase.

Graphical Abstract

1. Introduction

Nature harnesses transition metals to catalyze many different types of biological reactions that are crucial to life. Although about one-third of the enzymes are metalloenzymes, the use of nickel as a cofactor is not very common [1,2,3,4,5,6,7,8,9]. It was in 1975 that Zerner discovered the first nickel-dependent enzyme, namely urease [10]. Until now, ten different types of nickel enzymes (Table 1) have been reported with diverse biological functions. Among the known nickel enzymes, four of them are involved in the processing of anaerobic gases, namely, H2, CO, CO2, and CH4, which have been proposed to be crucial to the origin of life [7]. Ni–Fe hydrogenase mediates the reversible two-electron reduction of protons to H2 [11]. Methyl-coenzyme M reductase catalyzes the reversible transformation of methyl-coenzyme M plus coenzyme B into CH4 and heterodisulfide CoM–S–S–CoB [12]. CO dehydrogenase catalyzes the reversible oxidation of CO to CO2 using water as the oxygen source [13]. Acetyl-CoA synthase interacts tightly with CO dehydrogenase and uses the CO generated by CO dehydrogenase, CoA-SH, and a methyl group from a corrinoid/FeS protein to synthesize acetyl-CoA [14]. Acireductone dioxygenase [15], and quercetin 2,4-dioxygenase [16] are both nickel-dependent dioxygenases, and produce CO as one of the products. Ni-dependent superoxide dismutase targets the toxic reactive oxygen species superoxide and catalyzes its conversion to O2 and H2O2 [17]. All these enzymes belong to the oxidoreductase class, while the remaining three enzymes are involved in the hydrolysis of urea (urease) [10], the isomerization of methylglyoxal to lactate (glyoxylase I) [18], and the racemization of lactate (lactate racemase, or LarA) [19].
Understanding the reaction mechanisms of nickel-dependent enzymes is of fundamental and practical importance. A particularly interesting question is why these enzymes select nickel as a catalytic center. One important advantage of using nickel has been suggested to be related to its flexible coordination geometry [8]. For redox reactions, the ligand environment has been demonstrated to play a crucial role in tuning the redox potentials of nickel, ranging from +0.89 V in superoxide dismutase to −0.60 V in Methyl-coenzyme M reductase and CO dehydrogenase [6]. Sulfur-donor ligands, such as sulfides, thiolates of cysteine residues, are commonly used to fulfill this purpose. In addition, the Lewis acidic character of nickel has also been suggested to be important for urease, glyoxalase I, acireductone dioxygenase, and quercetin 2,4-dioxygenase [8]. Importantly, these enzymes preferentially adopt O/N-donor ligands. Another important question is how the enzyme catalyzes the reaction. To understand the mechanistic details, especially the structures of transition states and intermediates for all elementary steps, quantum chemical calculations have proven to be useful to complement experimental work.
Two popular but different approaches have been developed in the modeling of metalloenzymes. The first one is termed the quantum chemical cluster approach, developed by Siegbahn and Himo et al. [20,21,22,23,24,25,26]. More than 30 years of experience has demonstrated that a relatively small model of the active site, from about 30 atoms initially to more than 200 atoms nowadays, is capable of describing all important mechanistic features of the catalysis. The alternative approach is called the combined quantum mechanics/molecular mechanics (QM/MM) method, which was first proposed by Warshel and Levitt in 1976 [27]. It has been shown by Liao and Thiel that with proper selection and an increase of the number of atoms in the QM region, both approaches gave similar results and the same conclusion [28,29].
In the present review, we summarize the mechanistic studies on nickel enzymes by quantum chemical calculations.

2. NiFe Hydrogenase

In nature, hydrogenases are the enzymes designed to form dihydrogen from protons and electrons in a reversible process:
2H+ + 2e− ⇌ H2
The most common forms of hydrogenases are NiFe- and FeFe-hydrogenases. The NiFe-enzymes are primarily used for hydrogen oxidation, and the FeFe-enzymes for proton reduction. NiFe-hydrogenase is probably the Ni-enzyme that has been most intensively studied by theoretical methods. The active site [30] is shown in Figure 1. Apart from the four cysteines, there are three additional ligands on iron, two cyanides, and one carbon monoxide, which are very unusual in biological systems, since they are potentially poisonous. A complicated machinery involving many enzymes is, therefore, used to construct this complex.
The theoretical work prior to 2007 has been reviewed [31,32]. A variety of different models had been used, ranging from a minimal one, with only the direct ligands of the NiFe complex, to larger ones, with some charged second shell residues. All studies agreed on the general mechanism described in Figure 2. It contains three intermediate states observed experimentally, Nia–C*, Nia–S and Nia–SR. Nia–C* is the electron paramagnetic resonance (EPR)-active resting state, which has a bridging hydride between the metals. The EPR-silent Nia–SR state is reached by reduction of Nia–C*, while Nia–S, which is also EPR-silent, is obtained by oxidation. To complete the cycle, Nia–SR can be reached by adding dihydrogen to Nia–S. Several other states have also been observed experimentally under varying conditions. The TS for H–H cleavage was generally described as heterolytic, leading to a bridging hydride and a protonation of a cysteine, usually Cys543. This mechanism could be regarded as being in consensus among the theoretical groups [31,32], and had strong support from experiments. For the energy profile, it could be concluded that entropy and dispersion played important roles [31], even though the dispersion effect could not be routinely included in the calculations, as it can be now.
In the years following the early period before 2007, some new possibilities and aspects were investigated. This development has also been briefly reviewed [26,33]. First, from experiments, it was suggested that a variant of the mechanism in Figure 2 might be operative, leading to a mechanism in Figure 3 [34]. The main difference to the mechanism of Figure 2 described above was that Nia–S does not participate in the active cycle. Instead, catalysis starts only after an initial reaction between Nia–S and dihydrogen, where Nia–C* is created. Nia–C* is then only in equilibrium with Nia–SR during turnover. Some support for the alternative mechanism was given in one study [35]. An interesting state in the context of that mechanism was a Ni(I)-state (Nia–R*), appearing at the end of the single step heterolytic mechanism. A TS for oxidative addition was found leading to two hydrides, one bridging and one bound to nickel (Nia–X*).
Another development was that the interest switched over to the oxidized states. The reason for this interest is the possibility to generate fuel by combining water oxidation with H2 formation. A major problem in this context is that most hydrogenases are inactivated with O2 present. Two main oxidized states of the NiFe cofactor have been identified. In Ni–B, there is a bridging hydroxide between Ni(III) and Fe(II). It is inactive for H2 cleavage, but is relatively easy to reactivate. The Ni–A state is more difficult to reactivate. At the highest resolution (1.1 Å) it was initially concluded that Ni–A is a peroxide with one oxygen bridging between the metals and the other oxygen on nickel [30,36]. The early modeling studies were reviewed in 2007 [31]. While there was no problem obtaining the Ni–B structure with a bridging hydroxide [37,38], the modeling of Ni–A was more problematic. A structure similar to the one found experimentally could rather easily be obtained, but the energy was too high. A detailed analysis of the density of the experimental structure showed that it was likely to be a mixture of states, and the structure was probably also over reduced [39]. Calculations on the mechanism for the reactivation of the initially suggested structure for Ni–A were made [37]. It was suggested that the full reactivation of the enzyme requires two additional electrons and protons and this is the reason that it was slow. A reactivation mechanism for an oxidized cysteine ligand was also investigated. In 2015, a re-analysis of the experimental structure for Ni–A was performed, and it was concluded that Ni–A instead has a bridging hydroxide and an oxidized sulfenated form of Cys75 [40]. DFT calculations were then done for the new structure using larger models than before [41,42], and a revised mechanism for the reactivation was made.
Several oxygen insensitive NiFe hydrogenases have been found. In these enzymes, the Ni–A state is activated by two electrons supplied from the proximal FeS cluster [43]. A few DFT studies have been performed to determine the details of the mechanism of how the electrons are released from the FeS cluster [44,45,46]. The oxygen tolerant hydrogenases have been found to have specifically designed proximal Fe–S clusters, where a sulfide in a corner is missing and replaced by two cysteines. This allows the Fe–S cluster to deliver two electrons to the NiFe-cluster. Together with two electrons from the active site, the accumulation of the Ni–A state is diminished. DFT calculations were used to determine electronic and structural changes of the Fe–S cluster.
After these initial modeling studies, the interest in the mechanism of NiFe-hydrogenase has decreased. The studies mentioned above gave mechanisms that are quite similar with minor differences, and in good agreement with available experiments [31]. However, in recent years there was a set of studies by Ryde et al. in 2014–2018, where many of the results obtained previously were seriously criticized [47,48,49,50,51]. In particular, the accuracy of the Becke, 3-parameter, Lee–Yang–Parr exchange-correlation functional (B3LYP), used in practically all previous studies, was claimed to be insufficient. Comparisons between B3LYP and more advanced methods of ab initio type seemed to give results, showing that there were large errors using B3LYP. For the bare cofactor without hydrogen, the singlet–triplet splitting (S–T) error was tolerable with a usual error for B3LYP of a few kcal/mol compared to the advanced methods [47]. However, when H2 was added, the S–T splitting showed large errors for B3LYP [48], with a value of 13.2 kcal/mol compared to high level studies that gave 2–3 kcal/mol [49], but there are severe problems with that comparison. The most important one is that the structures were not optimized, but taken directly from a QM/MM study. This meant that the structures, as used in the comparison, are about 40 kcal/mol higher in energy than the optimized ones. The result most emphasized in the new DFT studies was that Nia–R* should not be an intermediate in the reaction sequence [51], as it was in some previous studies [35]. However, there is an interesting experimental result in this context using electro-chemical methods, where the conclusion was that Nia–R* was indeed observed as an intermediate [52]. Another more important aspect of the Ryde et al. studies is that they showed that it is actually possible to do calculations with the most accurate methods in quantum chemistry, such as the coupled cluster single-double and perturbative triple method (CCSD(T)) and the density matrix renormalization group method (DMRG). They have not yet been used for the reaction mechanism but it will be interesting to follow that development in the future.

3. Methyl-Coenzyme M Reductase

Methane, acting as the second most important anthropogenic greenhouse gas after CO2 [53,54], has a significant influence on the global carbon cycling and the climate [53]. With a methane emission annually of around 500–600 Tg, about 69% is biologically released by microbial metabolisms, i.e., methanogenic archaea [55]. Meanwhile, a large amount of methane is oxidatively consumed by aerobic methanotrophic bacteria and anaerobic methane oxidizing archaea (ANME) [56,57]. The anaerobic methane oxidation is usually coupled with the reduction of sulfate (sulfate-reducing δ-proteobacteria) [58,59] iron [60,61], manganese [61], nitrate [62], soluble metal complexes [63], etc. [64], which may critically affect the Proterozoic and Neoproterozoic global climate systems [65] and the breeding of animals [66]. The final step of microbial methane formation and the first step of anaerobic methane oxidation are both performed by the Ni-dependent methyl-coenzyme M reductase (MCR) [56,67], which is capable of reversibly catalyzing the conversion of methyl-coenzyme M (CH3–SCoM) and coenzyme B (CoB–SH) to a heterodisulfide (CoM–S–S–CoB) and methane using a Ni(I) porphyrinoid cofactor (F430) (Figure 4) [6,9]. It is further found that MCR-like enzymes may constitute a family of alkyl-coenzyme M reductases, which are essential for the microbial anaerobic oxidation of ethane [68], propane [56,68], and butane [69] following the same reversible reaction pattern as MCR, i.e., alkane + CoM–S–S–CoB ⇌ alkyl–SCoM + CoB–SH. Thus, the MCR mechanism may serve as a template for the family of alkyl-coenzyme M reductases [68], rendering it significant for the understanding of the alkane biosynthesis and the development of biomimetic catalysts for the alkane conversion.
Mainly four mechanisms for MCR have been proposed and examined by experiments or theoretical calculations (Figure 5). Accompanying the crystallization of MCR [12], mechanism I (Figure 5a) was hypothesized to include an organometallic Ni(III)-methyl intermediate (referred to as MCRMe) resulting from the Ni-activated C–S bond dissociation of CH3–SCoM and the proton transfer from CoB–SH to the sulfur of CH3–SCoM. In mechanism I, the first step of C–S bond dissociation is most likely not rate-limiting, since the following electron transfer was suggested to lead to an unstable CoM–S·H radical. Methane is formed in the next step of S–S bond formation. This mechanism received the strongest support from the direct synthesis of MCRMe from the active MCR (MCRred1) with methyl bromide [70] or iodide [71], although its formation from the native substrate (CH3–SCoM) has never been found. However, a rough estimation by means of the B3LYP bond dissociation energies (BDEs) predicted a very large endothermicity for the MCRMe formation in the native MCR reaction [72]. A more elaborate B3LYP investigation based on a chemical model built from the crystal structure (Figure 6) showed that the MCRMe formation is dependent on the acidity of the methyl-X leaving group, i.e., the C–X bond dissociation energy [73]. The MCRMe formation is facile with a barrier of a few kcal/mol when X is a halide (I−, Br−, or Cl−), while it has a very high barrier when starting from the native substrate with a stronger C–S bond, which was predicted to have a large endothermicity of 23.5 kcal/mol [73]. This endothermicity was then updated to 21.8 kcal/mol with dispersion and entropy effects added [74], which is sufficient to rule out mechanism I, considering the possible additional barriers in the subsequent steps (proton transfer and electron transfer).
A more detailed DFT cluster modeling of mechanism II (Figure 5b) was then performed [72,76]. Mechanism II proceeds through a C–S bond dissociation forming Ni(II)–SCoM and a methyl radical, the immediate quenching of the methyl radical by CoB–SH leading to methane and a CoB–S· radical. This is followed by a rebound mechanism between Ni(II)–SCoM and CoB–S· to form a CoM–S–S–CoB and a regeneration of Ni(I) F430. This mechanism has the following characteristics that are compatible with most of the mechanistic experiments: (i) The first step of C–S bond dissociation is rate-limiting with a barrier of 19.5 kcal/mol, which is supported by a substantial carbon kinetic isotope effects (12CH3–SCoM/13CH3–SCoM) of 1.04 ± 0.01 [77]. (ii) The methyl radical is captured rapidly by CoB–SH, so that methane formation happens almost simultaneously with the C–S bond dissociation and the whole MCR reaction can be regarded as two chemical steps via two transition states. (iii) An inversion of configuration at the methyl carbon takes place via a trigonal planar configuration during the methane formation; step (iii) combined with step (ii) is consistent with the finding using an isotopically chiral form of ethyl-coenzyme M [78] and a large α-secondary kinetic H/D isotope effect (kH/kD = 1.19 ± 0.01), which indicates a geometric change of the methyl group from tetrahedral to trigonal planar at the highest transition state [77]. However, the energetics obtained in 2002 and 2003 [72,76] are not accurate enough to give a quantitatively correct description of the MCR mechanism, mainly because of the use of methanol as the model of Tyr333 and Tyr367, the lack of dispersion effects, and the omission of entropy effects in the second step of the S–S bond formation [74]. This leads to an excessively high barrier of ~26 kcal/mol for the reverse reaction (i.e., the anaerobic oxidation of methane) and an excessively large barrier difference of 10.0 kcal/mol between the two transition states of C–S bond dissociation and S–S bond formation [72,76]. This process faced additional challenges posed by two new isotope exchange experiments. These new experiments proved the feasibility of the reverse oxidation of methane [79] and showed a small barrier difference between the two chemical steps [80]. With dispersion and more accurate entropy effects included, a new DFT investigation with methylphenols as tyrosine models gave more reasonable barriers of 15.1 and 21.1 kcal/mol for the forward and reverse reactions, respectively, and a smaller barrier difference of 3.6 kcal/mol for the two transition states [74], which demonstrated that mechanism II is still acceptable and in agreement with mechanistic experiments at that time. Interestingly, an ingenious spectroscopic experiment using a shorter CoB–SH with one methylene less in order to elongate the distance of CoB–SH to the reacting core, lowering the MCR reaction rate, leads to a direct observation of the Ni(II)–SCoM intermediate [81]. This provided the clearest experimental evidence for mechanism II and made mechanism II a consensus MCR mechanism [82].
An oxidative addition mechanism (i.e., Mechanism III in Figure 5c) proposed by DFT calculations [83] and experiments [80,84] involves a side-on C–S (for CH3–SCoM) [84] or C–H (for the reverse methane oxidation) [80] coordination to the Ni ion, forming a ternary intermediate with both methyl and CoM–S−/hydride coordinating to the metal. However, from a DFT study based on the same cluster model that was used in the previous calculations of mechanism II, the formation of a ternary intermediate was predicted to have a very large endothermicity of 44.5 kcal/mol starting from CH3–SCoM, and also a very high barrier starting from methane [75]. This can safely rule out mechanism III. Another proposal (mechanism IV in Figure 5d) suggested a protonation at the sulfur position of CH3–SCoM, which should facilitate the C–S bond activation and the subsequent methyl–Ni(III)F430 formation. Nevertheless, the DFT calculations showed that the substrate sulfur has a much smaller proton affinity than the F430 cofactor, in particular the carbonyl oxygen in F430, which means that a proton entering the MCR active site should prefer binding at F430 instead of CH3–SCoM, and that an additional energy of ~23 kcal/mol is required to move the proton from F430 to the substrate sulfur before the C–S bond activation [75]. This makes mechanism IV infeasible.
In summary, mechanism II (Figure 5b) is finally the most acceptable one for the MCR enzyme after two decades of computational and experimental studies. It fits reasonably well with several structural characteristics of the MCR active site. For example, a deeply-buried active site and the strict arrangement of F430, CH3–SCoM, and CoB–SH from bottom to top [12,85], ensuring the quick quenching of the methyl radical. In fact, the absence of CoB–SH cannot initiate the C–S bond cleavage and a CoB–SH analogue with one methylene less makes the MCR reaction much slower [86], which may be explained by the fast rebound between Ni(II)–SCoM and the unstable methyl radical. Even if the methyl radical escapes from CoB–SH, it cannot easily cause irreparable damage to the MCR active site, since the latter is hydrophobic and four active-site residues (His257, Arg271, Gln400, and Cys452) have been uncommonly methylated [12,85]. With this, MCR presents a perfect example of a mechanistic proposal by quantum chemical modeling finally leading to the convergence between theory and experiment. The understanding of the MCR mechanism could lead to new ways of producing methane/alkanes as fuels [87], of utilizing methane/alkanes as chemical feedstocks [82,87], and of breaking the C–S bond [88].

4. Nickel CO Dehydrogenase and Acetyl CoA Synthase

Carbon dioxide fixation is an extremely important process in nature, with implications on present day discussions on the greenhouse effect. The dominant enzymes under anaerobic conditions are Ni–CODH and Acetyl CoA synthase, which both have an active site containing Ni. Structures of the enzymes were first determined in 2002–2003 [14,89,90]. In the latest structure [13], obtained at a high resolution of 1.1 Å, an intermediate with bound CO2 was found. The Ni–CODH and Acetyl CoA synthase enzymes are contained in one complex, and are connected by a 140 Å long channel in which CO is transported [91,92]. In Ni–CODH, CO can be reversibly transformed to CO2. The reduction potential used for oxidation of CO is −0.3 V, while for reduction of CO2 it has to be decreased to −0.6 V. Experimentally, the rate for reduction is 10 s−1 at pH 7, and at pH 8 the oxidation rate is as high as 3.9 × 104.
The active site of Ni–CODH has an unusual NiFe3S4 complex, termed the C-cluster, connected to another iron termed Feu. The experimental X-ray structure with a bound CO2 [13], termed Cred1, has carbon bound to nickel and one of the oxygens bound to Feu, as shown in Figure 7. There are two important hydrogen bonds to CO2, with one being from the positively charged Lys563, and the other one from a His93–Asp219 couple, which is most probably protonated. The spin state is a doublet [93] indicating an oxidation state of Ni(II)Fe3(III)Fe(II) for the NiFe4 cluster. Theoretical studies have supported this assignment [94,95,96]. After adding two electrons, one C–O bond is cleaved and the state termed Cred2 is reached, which could be either Ni(I)Fe4(II) or Ni(0)Fe3(III)Fe(II). In the X-ray structure of Cred2, the CO, resulting from the C–O cleavage of CO2, has disappeared from the cluster leaving an empty site on nickel [13].
The first computational study of the mechanism of Ni–CODH was made by Amara et al. [96]. To explain the X-ray structures, they suggested an unusual mechanism, where in Cred2 a hydride occupies the empty site on nickel. The hydride would then have a role as an electron reservoir in Cred2. This would mean that nickel always stays in the Ni(II) state. A bound hydride would not be seen in the X-ray structure and might explain the observation of small amounts of H2 in the experiments. In the CO2 reduction direction from Cred2, CO2 was suggested to insert into the Ni–hydride bond. The protein matrix should prevent formate formation. Mössbauer parameter calculations, and an orbital and a charge distribution analysis, were also performed to support their mechanism.
Quite recently, the energetics for the full catalytic cycle of Ni–CODH were computed for the first time (Figure 8) [97]. Energy diagrams for both the reduction of CO2, using a redox potential of −0.6 V, and oxidation of CO, using a redox potential of −0.3 V, were calculated. In line with experiments, the oxidation is faster with a barrier difference of 4.8 kcal/mol, compared to experiments of 6.2 kcal/mol [92]. The rate-limiting step in both directions is the cleavage (formation) of the C–O bond of CO2. The strong hydrogen bond from Lys563 and the proton transfer from the His93–Asp219 couple are very important. To obtain accurate energetics, it was found that one of the sulfides in the NiFeS cluster should be protonated. A water molecule in the active site was also found to be reasonably important.
An important difference between the computed and experimental structures found was that CO was bound to nickel in the reduced structure, in contrast to the experimental structure. To explain this result, it was suggested that the X-ray structure was over-reduced. The calculations for the over-reduced structure showed that CO is initially exchanged with a hydride. The hydride and the water molecule then form a hydrogen molecule, suggested to explain the presence of small amounts of H2.
As mentioned above, acetyl CoA synthase (ACS) is located in the same complex as Ni–CODH, at the other end of the long channel through which CO is transported. The active catalyst in ACS, termed the A-cluster, is a nickel dimer connected by two sulfide bridges to an Fe4S4 cluster (Figure 9). The nickel closest to the FeS cluster is termed Nip and the other one Niu. Nip is, apart from the sulfide bridges to the FeS cluster, connected by two bridging cysteines to Niu. Niu is furthermore ligated to two backbone nitrogens of two connected cysteines. Nip obtains a square planar coordination with a Ni(II) oxidation state before the substrates enter. Nip is five-coordinated and has been suggested to be the active metal in the dimer. There are three substrates, CO, CH3, and coenzyme A (CoA). CH3 is donated to the cluster from a cobalamine. The order of the entering substrates is debated. There are reasons to believe that the order might matter, since the protein has at least two different conformations, of which only one—the open conformation—allows CO to be transported all the way to ACS. In the closed conformation the transport of CO is blocked.
There has been one theoretical study of the mechanism of ACS [98]. A cluster QM approach was used. A problem in the modeling is that no X-ray structures with bound substrates are available. Possible intermediates in the catalytic cycle were studied, including their oxidation states. Ni(0) had earlier been suggested as a possible state of Nip, but the study shows that this is unlikely. The zero-valent state is stabilized only if the FeS cluster is reduced to a state with S = 1/2 ([Fe4S4]1+). As CO binds, Nip(0) is still favored. If the FeS cluster is oxidized, Nip(I) is instead favored. It is suggested that CO and CH3 both bind to Nip in a mononuclear mechanism on Nip. There are major remaining issues in the mechanism, mostly because of the limited structural information. In particular, the transfer of the methyl from cobalamine to the nickel dimer would be very interesting to study in the future, but also many questions remain regarding the final stages of CoA synthase.
A theoretical study of the Ni-dimer in ACS with emphasis on pKa values and redox potentials has also been performed [99]. Based on these computed values, the initial part of a mechanism for ACS was suggested. The low reduction potential computed indicates that the reduced form is protonated. Before or during methylation of the A-cluster, the proton could be detached by a base. A proton coupled electron transfer (PCET) mechanism should be operational in the methylation step. Furthermore, a short early note on the ACS mechanism by another group has been published [100]. Their calculations suggest that Nip is better described by the oxidation state Ni(I) than Ni(0).

5. Acireductone Dioxygenase

Acireductone dioxygenases (ARDs) are a family of metalloenzymes that are capable of using different transition metal ions (Ni2+ and Fe2+) to catalyze the dioxygenation of acireductone [101,102]. Fe2+-ARD and Ni2+-ARD share exactly the same protein sequence. However, Fe2+-ARD catalyzes the formation of methylthioketobutyrate and formate [101], whereas Ni2+-ARD oxidizes acireductone into three products, namely formate, 3-methylthiopropionic acid, along with carbon monoxide [102]. It has been previously suggested that the substrate binding mode for Fe2+-ARD is a five-membered ring, while there is a six-membered ring for Ni2+-ARD (Scheme 1) [103], and this should be the main reason for the formation of different products.
Sparta et al. performed discrete molecular dynamics simulations (QM/DMD) to investigate the hypothesis of the different substrate binding modes in Fe2+-ARD and Ni2+-ARD [104]. The QM/DMD simulations suggested that both Fe2+-ARD and Ni2+-ARD exhibit a six-membered ring coordination, with the O1 and O3 atoms binding to the central metal ion (Scheme 2, React). A second-shell residue Arg154 forms a hydrogen bond to the O2 atom of the substrate, which stabilizes this coordination mode. One snapshot was then selected for DFT calculations in order to understand the mechanisms of Fe2+-ARD and Ni2+-ARD. The energy barrier for the first dioxygen attack was not considered. The starting complex is a dioxygen adduct, with a peroxy bridging C1 and C3 of the substrate (reaction, Scheme 2). For Ni2+-ARD, the reaction takes place in the triplet state. After a concerted cleavage of the O–O, C1–C2, and C2–C3 bonds, carbon monoxide and two carboxylates are generated (production, Scheme 2). The energy barrier is only 3.7 kcal/mol.
Recently, based on the structure of human acireductone dioxygenase, Borowski et al. have performed QM/MM calculations to revisit the reaction mechanisms of ARDs [105]. For the Ni2+-ARD, the calculations suggested that the mechanism involves four major steps (Scheme 3). A proton is first transferred from the substrate to His88 when the substrate enters into the active site, and His88 dissociates from the nickel ion. Dioxygen binds to the metal coupled with one electron transfer from the substrate to the dioxygen moiety to generate a Ni2+-superoxide intermediate (complex 1), which is quite similar to quercetin 2,3-dioxygenases [106,107]. The ground state of complex 1 is a quintet state (high-spin nickel(II) ions, one electron on the dioxygen, one electron on the substrate), while the following reaction is in the triplet state (two unpaired electrons on the nickel (II) ion). Hence, the energy of the triplet complex 1 was taken as the reference. After the distal oxygen atom attacks the C2 atom, a peroxo intermediate (complex 2) is formed. This step was found to be barrier-less. The next step is the migration of the C2-bound oxygen atom to the C3 atom (complex 3), which is ready for the following formation of the O–O bridge between C1 and C3 (complex 4). The final concerted cleavages of the O–O, C1–C2, and C2–C3 bonds were found to be the rate-limiting steps (TS45) for the substrate oxidation, with a barrier of only 6.8 kcal/mol. It is more likely that the substrate binding or product release is the rate-limiting step for the whole reaction.
Alternative pathways have also been considered from complex 2. The C2-bound oxygen atom may attack the C1 atom to form complex 6. After the formation of a dioxethane intermediate (complex 7) and the cleavage of the O–O and C1–C2 bonds (TS78), methylthioketobutyrate and formate are generated (complex 8). However, the barrier of TS78 was found to be 11.2 kcal/mol, higher than that of TS45. Other pathways that pass through a Baeyer-Villiger type rearrangement (2 → 6 → 10 → 11) and the homolytic O–O cleavage (2 → 6 → 9) were also ruled out due to their higher energy barriers compared with the most favorable mechanism.

6. Quercetin 2,4-Dioxygenase

Quercetin 2,4-dioxygenase from Streptomyces sp. strain FLA is a nickel-dependent enzyme catalyzing the oxidative ring-cleavage of quercetin to produce CO and 2-protocatechuoylphloroglucinol carboxylic acid using O2 as the oxidant [16,108]. The mechanism for the dioxygen activation and substrate oxidation was first investigated by Liu and co-workers [109]. From the 19 ns molecular dynamics simulations, one snapshot was chosen for the following QM/MM calculations. The first-shell ligand Glu76 set to be deprotonated as a proton is transferred from the neutral quercetin substrate. For the reactant complex with the O2 binding to the Ni ion in an end-on fashion, their calculations showed that the quintet is the ground state, the triplet lies 2.3 kcal/mol higher in energy, while the singlet was found to be much higher in energy. It is likely that the closed-shell singlet was obtained, while the open-shell singlet was not considered. In the triplet state, one electron is transferred from the substrate to the O2 moiety, leading to the formation of a Ni2+-superoxide-substrate radical species. Importantly, the reaction takes place in the triplet state involving four major steps (Scheme 4). First, the terminal oxygen of the superoxide attacks C2 of the substrate to form the first C–O bond, which has a barrier of only 5.0 kcal/mol. This is followed by a conformational change by a rotation around the newly-formed C2–O bond, and this step turned out to be rate-limiting, with a barrier of 19.9 kcal/mol. Subsequently, the terminal oxygen of the peroxide attacks C4 of the substrate, resulting in a peroxide-bridging intermediate. Finally, simultaneous cleavage of the O–O bond and two C–C bonds produces a CO molecule and the 2-protocatechuoylphloroglucinol carboxylic acid. When the central nickel ion is replaced by an iron ion, the reaction proceeds via a similar pathway in the quintet state. However, the last step becomes rate-limiting, associated with a very high barrier of 29.9 kcal/mol.
Liao and co-workers independently investigated the reaction mechanism of this enzyme almost at the same time [107]. In their QM/MM calculations, the model started directly from the X-ray conformation. Importantly, the first-shell ligand Glu74 was considered to be both neutral and deprotonated. The aim of that study was to elucidate the mechanism and also to rationalize the chemoselectivity, in which only 2,4-dioxygenolytic cleavage takes place, but not 2,3-dioxygenolytic cleavage. They found a similar mechanism as Liu et al. [109] when Glu74 is protonated. However, the last step was found to be rate-limiting with a barrier of 24.8 kcal/mol. Unexpectedly, the 2,3-dioxygenolytic pathway was found to be more favorable, with a barrier of 21.8 kcal/mol. Therefore, the model with a protonated Glu74 residue could not reproduce the observed chemoselectivity. Instead, a model with a deprotonated Glu74 was demonstrated to reproduce the chemoselectivity (Scheme 5). In this case, the total barrier for the 2,4-dioxygenolytic pathway decreases to 17.4 kcal/mol, while it increases to 30.6 kcal/mol for the 2,3-dioxygenolytic pathway. The calculated barrier of 17.4 kcal/mol is in good agreement with the experimental kinetic data [108], which gave a barrier of about 15 kcal/mol.

7. Urease

Urease is a nickel-containing enzyme that catalyzes the hydrolysis of urea to yield ammonia and carbon dioxide in the last step of nitrogen mineralization [110]. The uncatalyzed reaction has been suggested by Merz and Estiu [111] to proceed via an elimination pathway. For the enzyme catalyzed reaction, there have been three different mechanistic proposals. First, Karplus et al. suggested a mechanism in which a Nickel-coordinated terminal hydroxide performs the nucleophilic attack, while an adjacent protonated His320 residue delivers a proton to the leaving amino group [110,112]. Second, Benini et al. proposed a mechanism involving a nucleophilic attack of a bridging hydroxide on the carbonyl carbon of urea, followed by the protonation of the leaving amino group by this hydroxide [113]. Third, Lippard et al. suggested an elimination mechanism, similarly to the uncatalyzed reaction [114]. However, this elimination pathway could not explain the experimental observation of carbamate as the first product found by Zerner et al. [115]. A number of computational studies have, therefore, been conducted to elucidate the reaction mechanism of this enzyme [116,117,118,119,120,121,122,123].
The two different hydrolytic mechanisms were investigated by Merz et al. using DFT calculations (Scheme 6 and Scheme 7) [120]. For the bridging hydroxide attack mechanism, the calculations suggested that the reaction starts from a NiII2 complex with a bridging hydroxide (Scheme 6). In the reactant complex S3, the urea substrate binds in a bi-dentate fashion, where the carbonyl oxygen coordinates to Ni1, while one of the amino groups coordinate to Ni2. Both NiII ions are in the triplet state, and the complex is, thus, a quintet for the ferromagnetically-coupled state. The bridging hydroxide performs a nucleophilic attack on the carbonyl carbon of the urea substrate via TS1A, leading to the formation of a tetrahedral intermediate (I1A) with the oxyanion stabilized by one of the two nickel ions. This pathway is consistent with the proposal by Ciurli et al., on the basis of docking studies [120]. Then, a proton is delivered from the bridging hydroxide to the amino group of urea (TS2A), facilitated by Asp360. This leads to the production of an ammonia molecule and a carbamate. The second step was found to be rate-limiting, with a barrier of 19.7 kcal/mol. Protonation of the leaving ammonia group by a protonated His320 residue has also been taken into account, and this step is close to isothermic. It should be noted that the role of this histidine residue has very recently been confirmed by modulating its protonation state and conformation of the mobile flap [124]. After the release of the ammonia molecule, a water molecule binds to Ni2, and this water molecule then performs a nucleophilic attack on the carbon atom of the carbamate anion (TS4A), concertedly with a proton transfer from the water molecule to the nearby Asp360 residue. Finally, the protonated Asp360 residue delivers a proton to the bridging oxyanion (TS6A), coupled with the C–O bond cleavage. The di-nickel complex (I7A) with a bridging hydroxide is regenerated. It should be pointed out that this kind of mechanism agrees very well with the recent crystal structure of the complex between urea and fluoride-inhibited urease, in which an unreactive fluoride replaces the bridging hydroxide [125].
Nucleophilic attack by the nickel bound water molecule was also taken into consideration by Merz et al. [120], and the starting complex was a monodentate urea-bound adduct with the carbonyl oxygen coordinated to Ni1 (S1B, Scheme 7). The calculations suggested that this mechanism involves only two steps (Scheme 7): the Ni2-bound water molecule performs a nucleophilic attack on the urea carbonyl carbon, which is coupled with a proton transfer from the water molecule to the leaving ammonia, assisted by the bridging hydroxide (TS1B). The next step is the C–N bond cleavage with the release of an ammonia molecule. The second step is rate-limiting, with a barrier of 24.6 kcal/mol, which is higher than the one of the bridging attack by a hydroxide.
They have also performed DFT calculations to analyze the elimination mechanism of urease (Scheme 8) [122]. They started from the enzyme without the urea substrate, and the substrate binding was calculated to be endothermic by 5.0 kcal/mol. Then, deprotonation of the amino group of urea by the bridging hydroxide should take place (TS3uw) and a new bridging water molecule is generated (I3uw). In the next step, a proton is transferred from the bridging water to the other amino group of urea (TS4uw), with a barrier of 30.5 kcal/mol. Subsequently, a water molecule enters into the active site and forms hydrogen bonds with the two Ni-bound water molecules, which facilitates the release of NH4+ and OCN−. This mechanism has a substantially higher barrier than those for the two other mechanisms.

8. Lactate Racemase

Lactate racemase (LarA) has been identified to be the enzyme that catalyzes the final step in the racemization between d- and l-lactic acids [126], using a previously unknown Ni cofactor (Ni-PTTMN), pyridinium-3-thioamide-5-thiocarboxylic acid mononucleotide Ni pincer complex [19,127] (Figure 10). Ni-PTTMN is formed by the coordination between a Ni(II) ion and a prosthetic group derived from nicotinic acid (Figure 10). The latter has a thioamide (modified from the nicotinic moiety) covalently bound to Lys184 and a thiocarboxylate, and acts as a tridentate pincer ligand with the two sulfurs and one pyridine carbon ligating to the nickel, where a Ni–C bond is formed. The His200 works as the fourth nickel ligand. The Ni-PTTMN cofactor in LarA serves as a rare example of a pincer complex found in enzymes and its Ni–C bond is the first case of a metal-carbon bond originating from nicotinic acid in biology [128]. Considering the common and widespread NAD+-type cofactors (nicotinamide adenine dinucleotide and its phosphorylated form, NADP+) derived from nicotinic acid, it is surprising that such a complicated Ni-PTTMN cofactor has been biologically evolved with extra consumption of resources. A question is whether it has evolutional advantages compared with NAD+ [128].
A proton-coupled hydride transfer (PCHT) mechanism using a Ni(II) ion (Figure 10) has been verified by two DFT cluster modeling investigations by Chung et al. [130] and Chen et al. [129]. In this mechanism, a hydride is transferred to the Ni-PTTMN pyridine carbon from the substrate α-carbon, coupled with a proton transfer from the substrate hydroxyl to a histidine, followed by a transfer back to the α-carbon of the resultant pyruvate from the opposite side. This is achieved via a rotation of the pyruvate acetyl moiety. Chung et al. used an active-site model of 139 atoms to give a barrier of ~18 kcal/mol for the hydride transfer at the B3LYP-D3 level [130]. Based on small models (mainly involving cofactors and a proton-acceptor imidazole) and low-polarity solvation effects with a dielectric constant (ε) of 4, the further comparison of Ni-PTTMN with NAD+-like cofactors showed that the Ni-PTTMN cofactor has a higher racemization barrier than NAD+ [130], which makes it interesting to explore the LarA activity in the full active-site environment. In the subsequent DFT calculations by Chen et al. using a 200-atom model with Tyr294 and full side chains of the residues included, a lower rate-limiting barrier of 12 kcal/mol for the hydride transfer was predicted. The catalytic acceleration effect is achieved via the stabilization of the transition states by Tyr294 and Lys298 [129]. Further calculations were performed with various modified Ni-PTTMN cofactors in the LarA active site. They showed that the barrier of 12 kcal/mol for the native LarA with a Ni-PTTMN cofactor is much lower than those (at least by 19 kcal/mol) for the LarA enzymes with NAD+-like cofactors [129], indicating an enhanced racemization activity of Ni-PTTMN. To further analyze the acceleration effects of Ni-PTTMN, the hydride affinities of the Ni-PTTMN and NAD+-like cofactors were calculated [129]. It was revealed that compared with NAD+-like cofactors, Ni-PTTMN has a hydride affinity less sensitive to the environmental polarity, thus keeping a stronger hydride-addition reactivity in moderately and highly polar surroundings (ε ≥ 8), while a weaker one in a low-polarity environment (ε < 8) [129]. This not only explains why a higher barrier for Ni-PTTMN was obtained in the calculations by Chung et al. [130], but also indicates the evolutionary advantage of Ni-PTTMN using an architecture with a Ni pincer, which fits perfectly in the moderately polar active site of LarA, in which a dielectric constant of around 14–18 was estimated by a barrier analysis [130]. Inspired by the Ni-PTTMN-dependent LarA reaction, a series of Ni-PTTMN-like metal pincer complexes (metal = Ni and Pd) were proposed by Yang et al. using DFT calculations to be potential catalysts for lactate racemization with a barrier as low as ~26 kcal/mol, adopting the same PCHT mechanism as LarA [131], which is a step closer to utilizing the novel Ni-PTTMN cofactor and its related chemistry in man-made catalysis.
Another mechanism starting from a Ni(III) state was proposed in a QM/MM study by Shaik et al. [132]. It includes the electron transfer from the lactate carboxylate to Ni(III), the C–C bond dissociation leading to a CO2 radical anion and an acetaldehyde, the rotation of acetaldehyde, and then a C–C rebound step and the electron transfer back to the carboxylate. In this mechanism, the Ni-PTTMN cofactor works as a reversible electrode to accept and donate back an electron and does not result in covalent bonding to the lactate moiety. However, it is not clear if a Ni(II) Ni-PTTMN cofactor can be oxidized to the Ni(III) state in the LarA active site [133].

9. Superoxide Dismutase

Superoxide dismutases (SOD) catalyze the disproportionation of the toxic superoxide radicals into the less toxic hydrogen peroxide and molecular oxygen,
2O2−· + 2H+ → O2 + H2O2
They have an important protective role, particularly in aerobic metabolism and photosynthesis. There are at least three different types of SOD’s, and in one of them there is an active nickel complex.
The metal complex in Ni-SOD has a mononuclear nickel, coordinated to two cysteines and two backbone nitrogens in a square planar configuration (Figure 11). This is very similar to the coordination of Niu in acetyl CoA synthase (see above). However, there is one additional, axial, histidine ligand, which stabilizes a Ni(III) oxidation state. The disproportionation reaction occurs in two consecutive half-reactions, each one with one superoxide radical substrate. In the oxidative phase, the first superoxide radical binds and is transformed to dioxygen,
Ni(III)-L + O2−· + H+ → Ni(II)-LH + O2
In the reductive phase, the second superoxide substrate binds and is transformed to hydrogen peroxide,
Ni(II)-LH + O2−· + H+ → Ni(III)-L + H2O2
Shortly after the first high-resolution X-ray structure appeared in 2004 [17], there were two independent quantum chemical studies of the mechanism [134,135]. Essentially two different mechanisms had been suggested based on experiments. The most popular one was based on the findings on the structure formed by X-ray reduction [17,136]. It was found that upon reduction, the axial histidine lost its coordination to nickel. Therefore, it was suggested that the His became protonated as the proton mediator in the mechanism. The alternative is that one of the cysteines becomes protonated instead. The calculations used essentially minimal cluster models, including only the directly coordinating ligands and the backbone connecting the two cysteines.
The mechanism obtained in the first study [134] is shown in Figure 11. It was suggested that the superoxide substrate enters as an O2H· radical and binds in the empty axial site, where the cost for the protonation was included in the energetics. Anti-ferromagnetic coupling between nickel and the radical was found to be best. After substrate binding, the proton moves over to Cys2 and O2 leaves. Nickel has now been reduced to Ni(II). The second substrate radical then enters and binds at the same axial position. The proton on Cys2 is now well suited to move over to the O2H− ligand to form H2O2, and the disproportionation is completed. The protonation of His1 was found to be much less favorable. The barrier for the first half-reaction was found to be 9.7 kcal/mol and for the second one 11.5 kcal/mol, in good agreement with expectations based on experiments.
In the second study, a different mechanism was found for the first half-reaction, which was the only one studied [135]. A few second shell residues were included in the model apart from the directly binding ones. It was suggested that Cys6 started out protonated and donated its proton to the first superoxide. H2O2 can then be formed by a donation of a proton from the medium via His1 and Tyr9. A proton donating role of His1 was, therefore, suggested in line with earlier experimental suggestions, based on the X-ray reduced structure [17,136]. Transition states were not determined due to the large size of the model.
Two additional studies can be mentioned. Several years later in 2011, a study was made [137] that reached similar conclusions as the earlier ones. An energy diagram was calculated showing a rate-limiting step with a barrier for one of them of around 19 kcal/mol. In another study on a mimic of SOD, multireference effects were studied with the use of the complete active space self-consistent field method (CASSCF) [138]. It was found that the near-degeneracy effects were small, which was considered surprising at the time. Multi-reference effects were expected, since a radical became bound to a transition metal.

10. Conclusions

In this review, we have presented the progress of theoretical studies on the reaction mechanisms of nickel-dependent enzymes. Both quantum chemical cluster and QM/MM approaches have been successfully used to answer mechanistic questions in these enzymes. Detailed information can be obtained on the structures, energies, electron densities, and spin densities of all relevant intermediates and transition states.
Nickel is redox active and the common oxidation states are NiI, NiII, and NiIII in Ni enzymes. Due to its flexible coordination geometry, its redox potential can be tuned by different ligand coordinations with a span from +0.89 V to −0.60 V. Therefore, it can be used for challenging redox reactions, such as for superoxide dismutase, selective reduction of CO2, and methane formation. In addition, the nickel ion can also function as a Lewis acid to stabilize anionic intermediates and transition states during the reactions.

Author Contributions

Writing—original draft preparation, P.E.M.S., S.-L.C., and R.-Z.L.; writing—review and editing, P.E.M.S., S.-L.C., and R.-Z.L.

Funding

This research was funded by the Swedish Research Council, the Knut and Alice Wallenberg Foundation, the National Natural Science Foundation of China, grant numbers 21673019 and 21873031.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Walsh, C.T.; Orme-Johnson, W.H. Nickel Enzymes. Biochemistry 1987, 26, 4901–4906. [Google Scholar] [CrossRef] [PubMed]
  2. Ermler, U.; Grabarse, W.; Shima, S.; Goubeaud, M.; Thauer, R.K. Active sites of transition-metal enzymes with a focus on nickel. Curr. Opin. Struct. Biol. 1998, 8, 749–758. [Google Scholar] [CrossRef]
  3. Maroney, M.J. Structure/function relationships in nickel metallobiochemistry. Curr. Opin. Chem. Biol. 1999, 3, 188–199. [Google Scholar] [CrossRef]
  4. Mulrooney, S.B.; Hausinger, R.P. Nickel uptake and utilization by microorganisms. FEMS Microbiol. Rev. 2003, 27, 239–261. [Google Scholar] [CrossRef]
  5. Li, Y.; Zamble, D.B. Nickel Homeostasis and Nickel Regulation: An Overview. Chem. Rev. 2009, 109, 4617–4643. [Google Scholar] [CrossRef]
  6. Ragsdale, S.W. Nickel-based Enzyme Systems. J. Biol. Chem. 2009, 284, 18571–18575. [Google Scholar] [CrossRef] [Green Version]
  7. Fontecilla-Camps, J.C.; Amara, P.; Cavazza, C.; Nicolet, Y.; Volbeda, A. Structure–function relationships of anaerobic gas-processing metalloenzymes. Nature 2009, 460, 814–822. [Google Scholar] [CrossRef] [PubMed]
  8. Maroney, M.J.; Ciurli, S. Nonredox Nickel Enzymes. Chem. Rev. 2014, 114, 4206–4228. [Google Scholar] [CrossRef] [PubMed]
  9. Boer, J.L.; Mulrooney, S.B.; Hausinger, R.P. Nickel-dependent metalloenzymes. Arch. Biochem. Biophys. 2014, 544, 142–152. [Google Scholar] [CrossRef] [PubMed]
  10. Dixon, N.E.; Gazzola, C.; Blakeley, R.L.; Zerner, B. Jack Bean Urease (EC 3.5.1.5). A Metalloenzyme. A Simple Biological Role for Nickel? J. Am. Chem. Soc. 1975, 97, 4131–4133. [Google Scholar] [CrossRef]
  11. Volbeda, A.; Charon, M.H.; Piras, C.; Hatchikian, E.C.; Frey, M.; Fontecilla-Camps, J.C. Crystal structure of the nickel–iron hydrogenase from Desulfovibrio gigas. Nature 1995, 373, 580–587. [Google Scholar] [CrossRef] [PubMed]
  12. Ermler, U.; Grabarse, W.; Shima, S.; Goubeaud, M.; Thauer, R.K. Crystal Structure of Methyl-Coenzyme M Reductase: The Key Enzyme of Biological Methane Formation. Science 1997, 278, 1457–1462. [Google Scholar] [CrossRef] [PubMed]
  13. Jeoung, J.H.; Dobbek, H. Carbon dioxide activation at the Ni, Fe-cluster of anaerobic carbon monoxide dehydrogenase. Science 2007, 318, 1461–1464. [Google Scholar] [CrossRef] [PubMed]
  14. Doukov, T.; Iverson, T.M.; Seravalli, J.; Ragsdale, S.W.; Drennan, C.L. A Ni–Fe–Cu Center in a Bifunctional Carbon Monoxide Dehydrogenase/Acetyl-CoA Synthase. Science 2002, 298, 567–572. [Google Scholar] [CrossRef] [PubMed]
  15. Pochapsky, T.C.; Pochapsky, S.S.; Ju, T.; Hoefler, C.; Liang, J. A refined model for the structure of acireductone dioxygenase from Klebsiella ATCC 8724 incorporating residual dipolar couplings. J. Biomol. NMR. 2006, 34, 117–127. [Google Scholar] [CrossRef] [PubMed]
  16. Jeoung, J.H.; Nianios, D.; Fetzner, S.; Dobbek, H. Quercetin 2,4-Dioxygenase Activates Dioxygen in a Side-On O2–Ni Complex. Angew. Chem. Int. Ed. 2016, 55, 3281–3284. [Google Scholar] [CrossRef]
  17. Wuerges, J.; Lee, J.W.; Yim, Y.I.; Yim, H.S.; Kang, S.O.; Carugo, K.D. Crystal structure of nickel-containing superoxide dismutase reveals another type of active site. Proc. Natl Acad. Sci. USA 2004, 101, 8569–8574. [Google Scholar] [CrossRef] [Green Version]
  18. He, M.M.; Clugston, S.L.; Honek, J.F.; Matthews, B.W. Determination of the Structure of Escherichia coli Glyoxalase I Suggests a Structural Basis for Differential Metal Activation. Biochemistry 2000, 39, 8719–8727. [Google Scholar] [CrossRef]
  19. Desguin, B.; Zhang, T.; Soumillion, P.; Hols, P.; Hu, J.; Hausinger, R.P. A tethered niacin-derived pincer complex with a nickel–carbon bond in lactate racemase. Science 2015, 349, 66–69. [Google Scholar] [CrossRef]
  20. Siegbahn, P.E.M.; Blomberg, M.R. Transition-metal systems in biochemistry studied by high-accuracy quantumchemical methods. Chem. Rev. 2000, 100, 421–437. [Google Scholar] [CrossRef]
  21. Himo, F.; Siegbahn, P.E.M. Quantum chemical studies of radical containing enzymes. Chem. Rev. 2003, 103, 2421–2456. [Google Scholar] [CrossRef]
  22. Siegbahn, P.E.M.; Borowski, T. Modeling enzymatic reactions involving transition metals. Acc. Chem. Res. 2006, 39, 729–738. [Google Scholar] [CrossRef]
  23. Siegbahn, P.E.M.; Himo, F. Recent developments of the quantum chemical cluster approach for modeling enzyme reactions. J. Biol. Inorg. Chem. 2009, 14, 643–651. [Google Scholar] [CrossRef]
  24. Siegbahn, P.E.M.; Himo, F. The quantum chemical cluster approach for modeling enzyme reactions. WIREs Comput. Mol. Sci. 2011, 1, 323–336. [Google Scholar] [CrossRef]
  25. Himo, F. Recent trends in quantum chemical modeling of enzymatic reactions. J. Am. Chem. Soc. 2017, 139, 6780–6786. [Google Scholar] [CrossRef]
  26. Blomberg, M.R.; Borowski, T.; Himo, F.; Liao, R.Z.; Siegbahn, P.E.M. Quantum chemical studies of mechanisms for metalloenzymes. Chem. Rev. 2014, 114, 3601–3658. [Google Scholar] [CrossRef]
  27. Warshel, A.; Levitt, M. Theoretical studies of enzymatic reactions dielectric, electrostatic, and steric stabilization of carbonium-ion in reaction of lysozyme. J. Mol. Biol. 1976, 103, 227–249. [Google Scholar] [CrossRef]
  28. Liao, R.Z.; Thiel, W. Comparison of QM-only and QM/MM models for the mechanism of tungsten-dependent acetylene hydratase. J. Chem. Theory Comput. 2012, 8, 3793–3803. [Google Scholar] [CrossRef]
  29. Liao, R.Z.; Thiel, W. Convergence in the QM-only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase. J. Comput. Chem. 2013, 34, 2389–2397. [Google Scholar] [CrossRef] [Green Version]
  30. Volbeda, A.; Martin, L.; Cavazza, C.; Matho, M.; Faber, B.; Roseboom, W.; Albracht, S.; Garcin, E.; Rousset, M.; Fontecilla-Camps, J.C. Structural differences between the ready and unready oxidized states of [NiFe] hydrogenases. J. Biol. Inorg. Chem. 2005, 10, 239–249. [Google Scholar] [CrossRef]
  31. Siegbahn, P.E.M.; Tye, J.W.; Hall, M.B. Computational studies of [NiFe] and [FeFe] hydrogenases. Chem. Rev. 2007, 107, 4414–4435. [Google Scholar] [CrossRef]
  32. Bruschi, M.; Zampella, G.; Greco, Y.; Bertini, L.; Fantucci, P.; Gioia, D.L. Computational Inorganic and Bioinorganic Chemistry; Solomon, E.I., King, B., Scott, R., Eds.; John and Wiley and Sons: Chichester, UK, 2009. [Google Scholar]
  33. Siegbahn, P.E.M.; Blomberg, M.R.A. Quantum Chemical Studies of Proton-Coupled Electron Transfer in Metalloenzymes. Chem. Rev. 2010, 110, 7040–7061. [Google Scholar] [CrossRef]
  34. Kurkin, S.; George, S.J.; Thorneley, R.N.F.; Albracht, S.P.J. Hydrogen-induced activation of the [NiFe]-hydrogenase from Allochromatium vinosum as studied by stopped-flow infrared spectroscopy. Biochemistry 2004, 43, 6820–6831. [Google Scholar] [CrossRef]
  35. Nilsson-Lill, S.O.; Siegbahn, P.E.M. An autocatalytic Mechanism for NiFe-Hydrogenase: Reduction to Ni(I) Followed by Oxidative Addition. Biochemistry 2009, 48, 1056–1066. [Google Scholar] [CrossRef]
  36. Ogata, H.; Hirota, S.; Nakahara, A.; Komori, H.; Shibata, N.; Kato, T.; Kano, K.; Higuchi, Y. Activation process of [NiFe] hydrogenase elucidated by high-resolution X-ray analyses: Conversion of the ready to the unready state. Structure 2005, 13, 1635–1642. [Google Scholar] [CrossRef]
  37. Siegbahn, P.E.M. Hybrid density functional study of the oxidized states of NiFe-hydrogenase. C. R. Chim. 2007, 10, 766–774. [Google Scholar] [CrossRef]
  38. Pardo, A.; Lacey, A.L.D.; Fernandez, V.M.; Fan, Y.B.; Hall, M.B. Characterization of the active site of catalytically inactive forms of [NiFe] hydrogenases by density functional theory. J. Biol. Inorg. Chem. 2007, 12, 751–760. [Google Scholar] [CrossRef]
  39. Söderhjelm, P.; Ryde, U. Combined computational and crystallographic study of the oxidised states of [NiFe] hydrogenase. Theochem 2006, 770, 108–116. [Google Scholar] [CrossRef]
  40. Volbeda, A.; Martin, L.; Barbier, E.; Gutiérrez-Sanz, O.; De Lacey, A.L.; Liebgott, P.P.; Dementin, S.; Rousset, M.; Fontecilla-Camps, J.C. Crystallographic studies of [NiFe]-hydrogenase mutants: Towards consensus structures for the elusive unready oxidized states. J. Biol. Inorg. Chem. 2015, 20, 11–22. [Google Scholar] [CrossRef]
  41. Breglia, R.; Greco, C.; Fantucci, P.; De Gioia, L.; Bruschi, M. Theoretical investigation of aerobic and anaerobic oxidative inactivation of the [NiFe]-hydrogenase active site. Phys. Chem. Chem. Phys. 2018, 20, 1693–1706. [Google Scholar] [CrossRef]
  42. Breglia, R.; Greco, C.; Fantucci, P.; De Gioia, L. Reactivation of the Ready and Unready Oxidized States of [NiFe]-Hydrogenases: Mechanistic Insights from DFT Calculations. Inorg. Chem. 2019, 58, 279–293. [Google Scholar] [CrossRef]
  43. Goris, T.; Wait, A.F.; Saggu, M.; Fritsch, J.; Heidary, N.; Stein, M.; Zebger, I.; Lendzian, F.; Armstrong, F.A.; Friedrich, B.; et al. A unique iron–sulfur cluster is crucial for oxygen tolerance of a [NiFe]-hydrogenase. Nat. Chem. Biol. 2011, 7, 310–318. [Google Scholar] [CrossRef]
  44. Mouesca, J.M.; Fontecilla-Camps, J.C.; Amara, P. The Structural Plasticity of the Proximal [4Fe3S] Cluster is Responsible for the O2 Tolerance of Membrane-Bound [NiFe] Hydrogenases. Angew. Chem. Int. Ed. 2013, 52, 2002–2006. [Google Scholar] [CrossRef]
  45. Pandelia, M.E.; Bykov, D.; Izsak, R.; Infossi, P.; Giudici-Orticoni, M.T.; Bill, E.; Neese, F.; Lubitz, W. Electronic structure of the unique [4Fe–3S] clusterin O2-tolerant hydrogenases characterized by 57Fe Mössbauer and EPR spectroscopy. Proc. Natl Acad. Sci. USA 2013, 110, 483–488. [Google Scholar] [CrossRef]
  46. Pelmenschikov, V.; Kaupp, M. Redox-Dependent Structural Transformations of the [4Fe–3S] Proximal Cluster in O2-Tolerant Membrane-Bound [NiFe]-Hydrogenase: A DFT Study. J. Am. Chem. Soc. 2013, 135, 11809–11823. [Google Scholar] [CrossRef]
  47. Delcey, M.G.; Pierloot, K.; Quan, M.P.; Vancoille, S.; Lind, R.; Ryde, U. Accurate calculations of geometries and singlet–triplet energy differences for active-site models of [NiFe] hydrogenase. Phys. Chem. Chem. Phys. 2014, 16, 7927–7938. [Google Scholar] [CrossRef]
  48. Dong, G.; Quan, M.P.; Hallaert, S.D.; Pierloot, K.; Ryde, U. H2 binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods. Phys. Chem. Chem. Phys. 2017, 19, 10590–10601. [Google Scholar] [CrossRef]
  49. Dong, G.; Ryde, U. Protonation states of intermediates in the reaction mechanism of [NiFe] hydrogenase studied by computational methods. J. Biol. Inorg. Chem. 2016, 21, 383–394. [Google Scholar] [CrossRef] [Green Version]
  50. Dong, G.; Ryde, U.; Jensen, H.J.; Hedegård, E.D. Exploration of H2 binding to the [NiFe]-hydrogenase active site with multiconfigurational density functional methods. Phys. Chem. Chem. Phys. 2018, 20, 794–801. [Google Scholar] [CrossRef]
  51. Dong, G.; Quan, M.P.; Pierloot, K.; Ryde, U. Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods. Inorg. Chem. 2018, 57, 15289–15298. [Google Scholar] [CrossRef]
  52. Murphy, B.J.; Hidalgo, R.; Roessler, M.M.; Evans, R.M.; Ash, P.A.; Myers, W.K.; Vincent, K.A.; Armstrong, F.A. Discovery of Dark pH-Dependent H+ Migration in a [NiFe]-Hydrogenase and Its Mechanistic Relevance: Mobilizing the Hydrido Ligand of the Ni–C Intermediate. J. Am. Chem. Soc. 2015, 137, 8484–8489. [Google Scholar] [CrossRef]
  53. Murrell, J.C.; Jetten, M.S.M. The microbial methane cycle. Environ. Microbiol. Rep. 2009, 1, 279–284. [Google Scholar] [CrossRef]
  54. Hogan, K.B.; Hoffman, J.S.; Thompson, A.M. Methane on the greenhouse agenda. Nature 1991, 354, 181–182. [Google Scholar] [CrossRef]
  55. Conrad, R. The global methane cycle: Recent advances in understanding the microbial processes involved. Environ. Microbiol. Rep. 2009, 1, 285–292. [Google Scholar] [CrossRef]
  56. Evans, P.N.; Boyd, J.A.; Leu, A.O.; Woodcroft, B.J.; Parks, D.H.; Hugenholtz, P.; Tyson, G.W. Methane on the greenhouse agenda. Nat. Rev. Microbiol. 2019, 17, 219–232. [Google Scholar] [CrossRef]
  57. Thauer, R.K.; Kaster, A.K.; Seedorf, H.; Buckel, W.; Hedderich, R. Methanogenic archaea: Ecologically relevant differences in energy conservation. Nat. Rev. Microbiol. 2008, 6, 579–591. [Google Scholar] [CrossRef]
  58. Nauhaus, K.; Albrecht, M.; Elvert, M.; Boetius, A.; Widdel, F. In vitro cell growth of marine archaeal-bacterial consortia during anaerobic oxidation of methane with sulfate. Environ. Microbiol. 2007, 9, 187–196. [Google Scholar] [CrossRef]
  59. Nauhaus, K.; Boetius, A.; Krüger, M.; Widdel, F. In vitro demonstration of anaerobic oxidation of methane coupled to sulphate reduction in sediment from a marine gas hydrate area. Environ. Microbiol. 2002, 4, 296–305. [Google Scholar] [CrossRef]
  60. Ettwig, K.F.; Zhu, B.; Speth, D.; Keltjens, J.T.; Jetten, M.S.M.; Kartal, B. Archaea catalyze iron-dependent anaerobic oxidation of methane. Proc. Natl Acad. Sci. USA 2016, 113, 12792–12796. [Google Scholar] [CrossRef] [Green Version]
  61. Beal, E.J.; House, C.H.; Orphan, V.J. Manganese-and iron-dependent marine methane oxidation. Science 2009, 325, 184–187. [Google Scholar] [CrossRef]
  62. Haroon, M.F.; Hu, S.; Shi, Y.; Imelfort, M.; Keller, J.; Hugenholtz, P.; Yuan, Z.; Tyson, G.W. Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature 2013, 500, 567–570. [Google Scholar] [CrossRef]
  63. Scheller, S.; Yu, H.; Chadwick, G.L.; McGlynn, S.E.; Orphan, V.J. Artificial electron acceptors decouple archaeal methane oxidation from sulfate reduction. Science 2016, 351, 703–707. [Google Scholar] [CrossRef] [Green Version]
  64. Rotaru, A.E.; Thamdrup, B. A new diet for methane oxidizers. Science 2016, 351, 658–659. [Google Scholar] [CrossRef]
  65. Shen, B.; Dong, L.; Xiao, S.H.; Lang, X.G.; Huang, K.J.; Peng, Y.B.; Zhou, C.M.; Ke, S.; Liu, P.J. Molar tooth carbonates and benthic methane fluxes in Proterozoic oceans. Nat. Commun. 2016, 7, 10317. [Google Scholar] [CrossRef]
  66. Planavsky, N.J.; Reinhard, C.T.; Wang, X.; Thomson, D.; McGoldrick, P.; Rainbird, R.H.; Johnson, T.; Fischer, W.W.; Lyons, T.W. Low Mid-Proterozoic atmospheric oxygen levels and the delayed rise of animals. Science 2014, 346, 635–638. [Google Scholar] [CrossRef] [Green Version]
  67. Thauer, R.K. The Wolfe cycle comes full circle. Proc. Natl Acad. Sci. USA 2012, 109, 15084–15085. [Google Scholar] [CrossRef] [Green Version]
  68. Chen, S.C.; Musat, N.; Lechtenfeld, O.J.; Paschke, H.; Schmidt, M.; Said, N.; Popp, D.; Calabrese, F.; Stryhanyuk, H.; Jaekel, U.; et al. Anaerobic oxidation of ethane by archaea from a marine hydrocarbon seep. Nature 2019, 568, 108–111. [Google Scholar] [CrossRef]
  69. Pérez, R.L.; Wegener, G.; Knittel, K.; Widdel, F.; Harding, K.J.; Krukenberg, V.; Meier, D.V.; Richter, M.; Tegetmeyer, H.E.; Riedel, D.; et al. Thermophilic archaea activate butane via alkyl-coenzyme M formation. Nature 2016, 539, 396–401. [Google Scholar] [CrossRef]
  70. Yang, N.; Reiher, M.; Wang, M.; Harmer, J.; Duin, E.C. Formation of a nickel-methyl species in methyl-coenzyme M reductase, an enzyme catalyzing methane formation. J. Am. Chem. Soc. 2007, 129, 11028–11029. [Google Scholar] [CrossRef]
  71. Dey, M.; Telser, J.; Kunz, R.C.; Lees, N.S.; Ragsdale, S.W.; Hoffman, B.M. Biochemical and spectroscopic studies of the electronic structure and reactivity of a methyl-Ni species formed on methyl-coenzyme m reductase. J. Am. Chem. Soc. 2007, 129, 11030–11032. [Google Scholar] [CrossRef]
  72. Pelmenschikov, V.; Blomberg, M.R.A.; Siegbahn, P.E.M.; Crabtree, R.H. A mechanism from quantum chemical studies for methane formation in methanogenesis. J. Am. Chem. Soc. 2002, 124, 4039–4049. [Google Scholar] [CrossRef]
  73. Chen, S.L.; Pelmenschikov, V.; Blomberg, M.R.A.; Siegbahn, P.E.M. Is There a Ni-Methyl Intermediate in the Mechanism of Methyl-Coenzyme M Reductase? J. Am. Chem. Soc. 2009, 131, 9912–9913. [Google Scholar] [CrossRef]
  74. Chen, S.L.; Blomberg, M.R.A.; Siegbahn, P.E.M. How is methane formed and oxidized reversibly when catalyzed by Ni-containing methyl-coenzyme M reductase? Chem. Eur. J. 2012, 18, 6309–6315. [Google Scholar] [CrossRef]
  75. Chen, S.L.; Blomberg, M.R.A.; Siegbahn, P.E.M. An investigation of possible competing mechanisms for Ni-containing methyl-coenzyme M reductase. Phys. Chem. Chem. Phys. 2014, 16, 14029–14035. [Google Scholar] [CrossRef]
  76. Pelmenschikov, V.; Siegbahn, P.E.M. Catalysis by methyl-coenzyme M reductase: A theoretical study for heterodisulfide product formation. J. Biol. Inorg. Chem. 2003, 8, 653. [Google Scholar] [CrossRef]
  77. Scheller, S.; Goenrich, M.; Thauer, R.K.; Jaun, B. Methyl-coenzyme M reductase from methanogenic archaea: Isotope effects on label exchange and ethane formation with the homologous substrate ethyl-coenzyme M. J. Am. Chem. Soc. 2013, 135, 14975–14984. [Google Scholar] [CrossRef]
  78. Ahn, Y.; Krzycki, J.A.; Floss, H.G. Steric course of the reduction of ethyl coenzyme M to ethane catalyzed by methyl coenzyme m reductase from methanosarcina barkeri. J. Am. Chem. Soc. 1991, 113, 4700–4701. [Google Scholar] [CrossRef]
  79. Scheller, S.; Goenrich, M.; Boecher, R.; Thauer, R.K.; Jaun, B. The key nickel enzyme of methanogenesis catalyses the anaerobic oxidation of methane. Nature 2010, 465, 606–609. [Google Scholar] [CrossRef]
  80. Scheller, S.; Goenrich, M.; Mayr, S.; Thauer, R.K.; Jaun, B. Intermediates in the Catalytic Cycle of Methyl Coenzyme? M Reductase: Isotope Exchange is Consistent with Formation of a σ-Alkane-Nickel Complex. Angew. Chem. Int. Ed. 2010, 49, 8112–8115. [Google Scholar] [CrossRef]
  81. Wongnate, T.; Sliwa, D.; Ginovska, B.; Smith, D.; Wolf, M.W.; Lehnert, N.; Raugei, S.; Ragsdale, S.W. The radical mechanism of biological methane synthesis by methyl-coenzyme M reductase. Science 2016, 352, 953–958. [Google Scholar] [CrossRef] [Green Version]
  82. Lawton, T.J.; Rosenzweig, A.C. Methane-make it or break it. Science 2016, 352, 892–893. [Google Scholar] [CrossRef]
  83. Duin, E.C.; Mckee, M.L. A new mechanism for methane production from methyl-coenzyme M reductase as derived from density functional calculations. J. Phys. Chem. B 2008, 112, 2466–2482. [Google Scholar] [CrossRef]
  84. Harmer, J.; Finazzo, C.; Piskorski, R.; Ebner, S.; Duin, E.C.; Goenrich, M.; Thauer, R.K.; Reiher, M.; Schweiger, A.; Hinderberger, D.; et al. A nickel hydride complex in the active site of methyl-coenzyme m reductase: Implications for the catalytic cycle. J. Am. Chem. Soc. 2008, 130, 10907–10920. [Google Scholar] [CrossRef]
  85. Grabarse, W.; Mahlert, F.; Duin, E.C.; Goubeaud, M.; Shima, S.; Thauer, R.K.; Lamzin, V.; Ermler, U. On the mechanism of biological methane formation: Structural evidence for conformational changes in methyl-coenzyme M reductase upon substrate binding. J. Mol. Biol. 2001, 309, 315–330. [Google Scholar] [CrossRef]
  86. Horng, Y.C.; Becker, D.F.; Ragsdale, S.W. Mechanistic studies of methane biogenesis by methyl-coenzyme M reductase: Evidence that coenzyme B participates in cleaving the C–S bond of methyl-coenzyme M. Biochemistry 2001, 40, 12875–12885. [Google Scholar] [CrossRef]
  87. Borman, S. Carbon dioxide hydrogenated to methanol on large scale. Chem. Eng. News 2016, 94, 7. [Google Scholar]
  88. Shima, S. The Biological Methane-Forming Reaction: Mechanism Confirmed Through Spectroscopic Characterization of a Key Intermediate. Angew. Chem. Int. Ed. 2016, 55, 13648–13649. [Google Scholar] [CrossRef]
  89. Dobbek, H.; Svetlitchnyi, V.; Gremer, L.; Huber, R.; Meyer, O. Crystal Structure of a Carbon Monoxide Dehydrogenase Reveals a [Ni–4Fe–5S] Cluster. Science 2001, 293, 1281–1285. [Google Scholar] [CrossRef]
  90. Darnault, C.; Volbeda, A.; Kim, E.J.; Legrand, P.; Vernede, X.; Lindahl, P.A.; Fontecilla-Camps, J.C. Ni–Zn–[Fe4–S4] and Ni–Ni–[Fe4–S4] clusters in closed and open α subunits of acetyl-CoA synthase/carbon monoxide dehydrogenase. Nat. Struct. Biol. 2003, 10, 271–279. [Google Scholar] [CrossRef]
  91. Appel, A.M.; Bercaw, J.E.; Bocarsly, A.B.; Dobbek, H.; DuBois, D.L.; Dupuis, M.; Ferry, J.G.; Fujita, E.; Hille, R.; Kenis, P.J.A.; et al. Frontiers, Opportunities, and Challenges in Biochemical and Chemical Catalysis of CO2 Fixation. Chem. Rev. 2013, 113, 6621–6658. [Google Scholar] [CrossRef]
  92. Can, M.; Armstrong, F.A.; Ragsdale, S.W. Structure, Function, and Mechanism of the Nickel Metalloenzymes, CO Dehydrogenase, and Acetyl-CoA Synthase. Chem. Rev. 2014, 114, 4149–4174. [Google Scholar] [CrossRef]
  93. Lindahl, P.A.; Munck, E.; Ragsdale, S.W. CO Dehydrogenase from Clostridium thermoaceticum: EPR and electrochemical studies in CO2 and argon atmospheres. J. Biol. Chem. 1990, 265, 3873–3879. [Google Scholar]
  94. Cao, Z.-X.; Mo, Y.-R. Computational Characterization of the Elusive C-Cluster of Carbon Monoxide Dehydrogenase. J. Theor. Comput. Chem. 2008, 7, 473–484. [Google Scholar] [CrossRef]
  95. Xie, H.-J.; Cao, Z.-X. Computational Characterization of Reactive Intermediates of Carbon Monoxide Dehydrogenase. Chin. J. Struct. Chem. 2009, 28, 1525–1532. [Google Scholar]
  96. Amara, P.; Mouesca, J.-M.; Volbeda, A.; Fontecilla-Camps, J.C. Carbon Monoxide Dehydrogenase Reaction Mechanism: A Likely Case of Abnormal CO2 Insertion to a Ni–H− Bond. Inorg. Chem. 2011, 50, 1868–1878. [Google Scholar] [CrossRef]
  97. Liao, R.-Z.; Siegbahn, P.E.M. Energetics for the mechanism of nickel-containing CO-dehydrogenase. Inorg. Chem. 2019. [Google Scholar] [CrossRef]
  98. Amara, P.; Volbeda, A.; Fontecilla-Camps, J.C.; Field, M.J. A quntum chemical study of the reaction mechanism of Acetyl-Coenzyme A synthase. J. Am. Chem. Soc. 2005, 127, 2776–2784. [Google Scholar] [CrossRef]
  99. Chmielowska, A.; Lodowski, P.; Jaworska, M. Redox Potentials and Protonation of the A-Cluster from Acetyl-CoA Synthase. A Density Functional Theory Study. J. Phys. Chem. 2013, 117, 12484–12496. [Google Scholar] [CrossRef]
  100. Schenker, R.P.; Brunold, T.C. Computational studies on the A cluster of acetyl-coenzyme A synthase: Geometric and electronic properties of the NiFeC species and mechanistic implications. J. Am. Chem. Soc. 2003, 125, 13962–13963. [Google Scholar] [CrossRef]
  101. Dai, Y.; Wensink, P.C.; Abeles, R.H.J. One protein, two enzymes. Biol. Chem. 1999, 274, 1193–1195. [Google Scholar] [CrossRef]
  102. Wray, J.W.; Abeles, R.H. The Methionine Salvage Pathway in Klebsiella pneumoniae and Rat Liver identification and characterization of two novel dioxygenases. J. Biol. Chem. 1995, 270, 3147–3153. [Google Scholar] [CrossRef]
  103. Dai, Y.; Pochapsky, T.C.; Abeles, R.H. Mechanistic studies of two dioxygenases in the methionine salvage pathway of Klebsiella pneumonia. Biochemistry 2001, 40, 6379–6387. [Google Scholar] [CrossRef]
  104. Sparta, M.; Valdez, C.E.; Alexandrova, A.N. Metal-dependent activity of Fe and Ni acireductone dioxygenases: How two electrons reroute the catalytic pathway. J. Mol. Biol. 2013, 425, 3007–3018. [Google Scholar] [CrossRef]
  105. Miłaczewska, A.; Kot, E.; Amaya, J.A.; Makris, T.M.; Zając, M.; Korecki, J.; Chumakov, A.; Trzewik, B.; Kędracka-Krok, S.; Minor, W.; et al. On the Structure and Reaction Mechanism of Human Acireductone Dioxygenase. Chem. Eur. J. 2018, 24, 5225–5237. [Google Scholar] [CrossRef]
  106. Siegbahn, P.E.M. Hybrid DFT study of the mechanism of quercetin 2,3-dioxygenase. Inorg. Chem. 2004, 43, 5944–5953. [Google Scholar] [CrossRef]
  107. Wang, W.J.; Wei, W.J.; Liao, R.Z. Deciphering the chemoselectivity of nickel-dependent quercetin 2, 4-dioxygenase. Phys. Chem. Chem. Phys. 2018, 20, 15784–15794. [Google Scholar] [CrossRef]
  108. Merkens, H.; Kappl, R.; Jakob, R.P.; Schmid, F.X.; Fetzner, S. Quercetinase QueD of Streptomyces sp. FLA, a Monocupin Dioxygenase with a Preference for Nickel and Cobalt. Biochemistry 2008, 47, 12185–12196. [Google Scholar] [CrossRef]
  109. Li, H.; Wang, X.; Tian, G.; Liu, Y. Insights into the dioxygen activation and catalytic mechanism of the nickel-containing quercetinase. Catal. Sci. Technol. 2018, 8, 2340–2351. [Google Scholar] [CrossRef]
  110. Karplus, P.A.; Pearson, M.A.; Hausinger, R.P. 70 years of crystalline urease: What have we learned? Acc. Chem. Res. 1997, 30, 330–337. [Google Scholar] [CrossRef]
  111. Estiu, G.; Merz, K.M. The hydrolysis of urea and the proficiency of urease. J. Am. Chem. Soc. 2004, 126, 6932–6944. [Google Scholar] [CrossRef]
  112. Pearson, M.A.; Park, I.S.; Schaller, R.A.; Michel, L.O.; Karplus, P.A.; Hausinger, R.P. Kinetic and Structural Characterization of Urease Active Site Variants. Biochemistry 2000, 39, 8575–8584. [Google Scholar] [CrossRef]
  113. Benini, S.; Rypniewski, W.R.; Wilson, K.S.; Miletti, S.; Ciurli, S.; Mangani, S. A new proposal for urease mechanism based on the crystal structures of the native and inhibited enzyme from Bacillus pasteurii: Why urea hydrolysis costs two nickels. Structure 1999, 7, 205–216. [Google Scholar] [CrossRef]
  114. Barrios, A.M.; Lippard, S.J. Interaction of urea with a hydroxide-bridged dinuclear nickel center: An alternative model for the mechanism of urease. J. Am. Chem. Soc. 2000, 122, 9172–9177. [Google Scholar] [CrossRef]
  115. Blakeley, R.L.; Hinds, J.A.; Kunze, H.E.; Webb, E.C.; Zerner, B. Jack bean urease (EC 3.5.1.5) demonstration of a carbamoyl-transfer reaction and inhibition by hydroxamic acids. Biochemistry 1969, 8, 1991–2000. [Google Scholar] [CrossRef]
  116. Smyj, R.P. A conformational analysis study of a nickel(II) enzyme: Urease. J. Mol. Struct. THEOCHEM 1997, 391, 207–223. [Google Scholar] [CrossRef]
  117. Zimmer, M. Molecular mechanics evaluation of the proposed mechanisms for the degradation of urea by urease. J. Biomol. Struct. Dyn. 2000, 17, 787–797. [Google Scholar] [CrossRef]
  118. Zimmer, M. Are classical molecular mechanics calculations still useful in bioinorganic simulations? Coord. Chem. Rev. 2009, 253, 817–826. [Google Scholar] [CrossRef]
  119. Musiani, F.; Arnofi, E.; Casadio, R.; Ciurli, S. Structure-based computational study of the catalytic and inhibition mechanisms of urease. J. Biol. Inorg. Chem. 2001, 6, 300–314. [Google Scholar] [CrossRef]
  120. Suárez, D.; Díaz, N.; Merz, K.M. Ureases: Quantum chemical calculations on cluster models. J. Am. Chem. Soc. 2003, 125, 15324–15337. [Google Scholar] [CrossRef]
  121. Estiu, G.; Merz, K.M. Catalyzed decomposition of urea. Molecular dynamics simulations of the binding of urea to urease. Biochemistry 2006, 45, 4429–4443. [Google Scholar] [CrossRef]
  122. Estiu, G.; Merz, K.M. Competitive hydrolytic and elimination mechanisms in the urease catalyzed decomposition of urea. J. Phys. Chem. B 2007, 111, 10263–10274. [Google Scholar] [CrossRef]
  123. Carlsson, H.; Nordlander, E. Computational Modeling of the Mechanism of Urease. Bioinorg. Chem. Appl. 2010, 2010, 8. [Google Scholar] [CrossRef]
  124. Mazzei, L.; Cianci, M.; Benini, S.; Ciurli, S. The impact of pH on catalytically critical protein conformational changes; the case of the urease, a nickel enzyme. Chem. Eur. J. 2019. [Google Scholar] [CrossRef]
  125. Mazzei, L.; Cianci, M.; Benini, S.; Ciurli, S. The structure of the elusive urease-urea complex unveils the mechanism of a paradigmatic nickel-dependent enzyme. Angew. Chem. Int. Ed. 2019, 131, 7493–7497. [Google Scholar] [CrossRef]
  126. Desguin, B.; Goffin, P.; Viaene, E.; Kleerebezem, M.; Diaconescu, V.M.; Maroney, M.J.; Declercq, J.P.; Soumillion, P.; Hols, P. Lactate racemase is a nickel-dependent enzyme activated by a widespread maturation system. Nat. Commun. 2014, 5, 3615. [Google Scholar] [CrossRef]
  127. Xu, T.; Bauer, G.; Hu, X. A novel nickel pincer complex in the active site of lactate racemase. ChemBioChem 2016, 17, 31–32. [Google Scholar] [CrossRef]
  128. Zamble, D. It costs more than a nickel. Science 2015, 349, 35–36. [Google Scholar] [CrossRef]
  129. Yu, M.J.; Chen, S.L. From NAD+ to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase. Chem. Eur. J. 2017, 23, 7545–7557. [Google Scholar] [CrossRef]
  130. Zhang, X.; Chung, L.W. Alternative Mechanistic Strategy for Enzyme Catalysis in a Ni-Dependent Lactate Racemase (LarA): Intermediate Destabilization by the Cofactor. Chem. Eur. J. 2017, 23, 3623–3630. [Google Scholar] [CrossRef]
  131. Qiu, B.; Yang, X.Z. A bio-inspired design and computational prediction of scorpion-like SCS nickel pincer complexes for lactate racemization. Chem. Commun. 2017, 53, 11410–11413. [Google Scholar] [CrossRef]
  132. Wang, B.; Shaik, S. The Nickel-Pincer Complex in Lactate Racemase Is an Electron Relay and Sink that acts through Proton-Coupled Electron Transfer. Angew. Chem. Int. Ed. 2017, 56, 10098–10102. [Google Scholar] [CrossRef]
  133. Zhang, T.H.; Zhang, X.Y.; Chung, L.W. Computational Insights into the Reaction Mechanisms of Nickel-Catalyzed Hydrofunctionalizations and Nickel-Dependent Enzymes. Asian J. Organ. Chem. 2018, 7, 522–536. [Google Scholar] [CrossRef]
  134. Pelmenschikov, V.; Siegbahn, P.E.M. Nickel superoxide dismutase reaction mechanism studied by hybrid density functional methods. J. Am. Chem. Soc. 2006, 128, 7466–7475. [Google Scholar] [CrossRef]
  135. Prabhakar, R.; Morokuma, K.; Musaev, D.G. A DFT study of the mechanism of Ni superoxide dismutase (NiSOD): Role of the active site cysteine-6 residue in the oxidative half-reaction. J. Comput. Chem. 2006, 27, 1438–1445. [Google Scholar] [CrossRef]
  136. Barondeau, D.P.; Kassmann, C.J.; Bruns, C.K.; Tainer, J.A.; Getzoff, E.D. Nickel superoxide dismutase structure and mechanism. Biochemistry 2004, 43, 8038–8047. [Google Scholar] [CrossRef]
  137. Wang, Q.; Chen, D.; Liu, X.; Zhang, L.F. Theoretical mechanisms of the superoxide radical anion catalyzed by the nickel superoxide dismutase. Comput. Theor. Chem. 2011, 966, 357–363. [Google Scholar] [CrossRef]
  138. Shearer, J.; Schmitt, J.C.; Clewett, H.S. Adiabaticity of the Proton-Coupled Electron-Transfer Step in the Reduction of Superoxide Effected by Nickel-Containing Superoxide Dismutase Metallopeptide-Based Mimics. J. Phys. Chem. B, 2015, 119, 5453–5461. [Google Scholar] [CrossRef]
Inorganics 07 00095 g001 550
Figure 1. The X-ray structure of the active site of NiFe hydrogenase (Protein Data Bank (PDB) code: 1YRQ) [30].
Figure 1. The X-ray structure of the active site of NiFe hydrogenase (Protein Data Bank (PDB) code: 1YRQ) [30].
Inorganics 07 00095 g001
Inorganics 07 00095 g002 550
Figure 2. Schematic drawing of the catalytic cycles for NiFe-hydrogenases, suggested by experiments. Adapted from [31], Copyright © 2007, American Chemical Society.
Figure 2. Schematic drawing of the catalytic cycles for NiFe-hydrogenases, suggested by experiments. Adapted from [31], Copyright © 2007, American Chemical Society.
Inorganics 07 00095 g002
Inorganics 07 00095 g003 550
Figure 3. A variant of the mechanism for NiFe-hydrogenases. Adapted from [31], Copyright © 2007, American Chemical Society.
Figure 3. A variant of the mechanism for NiFe-hydrogenases. Adapted from [31], Copyright © 2007, American Chemical Society.
Inorganics 07 00095 g003
Inorganics 07 00095 g004 550
Figure 4. The reaction catalyzed by methyl-coenzyme M reductase (MCR) using a F430 cofactor.
Figure 4. The reaction catalyzed by methyl-coenzyme M reductase (MCR) using a F430 cofactor.
Inorganics 07 00095 g004
Inorganics 07 00095 g005 550
Figure 5. Four proposed mechanisms for the MCR reaction. (a) Mechanism I, including a methyl-Ni(III)F430 intermediate. (b) Mechanism II, involving a methyl radical and a CoM–S–Ni(II)F430 intermediate. (c) Mechanism III ,including a ternary intermediate with side-on methyl and CoM–S− coordinating to Ni. (d) Mechanism IV, involving a proton binding at the sulfur of CH3–SCoM. The energies for minima and transition states are given in brackets and parentheses, respectively. The energies for mechanisms I and II are adapted from [74], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, while the ones for mechanisms III and IV are adapted from [75], with permission from the PCCP Owner Societies.
Figure 5. Four proposed mechanisms for the MCR reaction. (a) Mechanism I, including a methyl-Ni(III)F430 intermediate. (b) Mechanism II, involving a methyl radical and a CoM–S–Ni(II)F430 intermediate. (c) Mechanism III ,including a ternary intermediate with side-on methyl and CoM–S− coordinating to Ni. (d) Mechanism IV, involving a proton binding at the sulfur of CH3–SCoM. The energies for minima and transition states are given in brackets and parentheses, respectively. The energies for mechanisms I and II are adapted from [74], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, while the ones for mechanisms III and IV are adapted from [75], with permission from the PCCP Owner Societies.
Inorganics 07 00095 g005
Inorganics 07 00095 g006 550
Figure 6. Model used in the calculations of the methyl-coenzyme M reductase (MCR). Adapted from [74], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Figure 6. Model used in the calculations of the methyl-coenzyme M reductase (MCR). Adapted from [74], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Inorganics 07 00095 g006
Inorganics 07 00095 g007 550
Figure 7. The X-ray structure of the active site of Ni–CODH with bound CO2 [13].
Figure 7. The X-ray structure of the active site of Ni–CODH with bound CO2 [13].
Inorganics 07 00095 g007
Inorganics 07 00095 g008 550
Figure 8. Mechanism of nickel CO dehydrogenase suggested on the basis of density functional calculations. Adapted from [97], Copyright © 2019, American Chemical Society.
Figure 8. Mechanism of nickel CO dehydrogenase suggested on the basis of density functional calculations. Adapted from [97], Copyright © 2019, American Chemical Society.
Inorganics 07 00095 g008
Inorganics 07 00095 g009 550
Figure 9. Mechanism of Acetyl CoA suggested on the basis of quantum mechanics (QM) cluster calculations. Adapted from [98], Copyright © 2005, American Chemical Society.
Figure 9. Mechanism of Acetyl CoA suggested on the basis of quantum mechanics (QM) cluster calculations. Adapted from [98], Copyright © 2005, American Chemical Society.
Inorganics 07 00095 g009
Inorganics 07 00095 sch001 550
Scheme 1. The previously suggested mechanisms for ARDs. Adapted from [103], Copyright © 2001, American Chemical Society.
Scheme 1. The previously suggested mechanisms for ARDs. Adapted from [103], Copyright © 2001, American Chemical Society.
Inorganics 07 00095 sch001
Inorganics 07 00095 sch002 550
Scheme 2. Reaction mechanism for Ni-acireductone dioxygenase supported by the results of a DFT study. Adapted from [104], Copyright © 2013 Elsevier Ltd.
Scheme 2. Reaction mechanism for Ni-acireductone dioxygenase supported by the results of a DFT study. Adapted from [104], Copyright © 2013 Elsevier Ltd.
Inorganics 07 00095 sch002
Inorganics 07 00095 sch003 550
Scheme 3. Reaction mechanisms for Ni-ARD supported by the results of a quantum mechanics/molecular mechanics (QM/MM) study. Adapted from [105], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Scheme 3. Reaction mechanisms for Ni-ARD supported by the results of a quantum mechanics/molecular mechanics (QM/MM) study. Adapted from [105], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Inorganics 07 00095 sch003
Inorganics 07 00095 sch004 550
Scheme 4. Reaction mechanism for Ni-QueD suggested by Liu et al. Adapted from [109], with permission from The Royal Society of Chemistry.
Scheme 4. Reaction mechanism for Ni-QueD suggested by Liu et al. Adapted from [109], with permission from The Royal Society of Chemistry.
Inorganics 07 00095 sch004
Inorganics 07 00095 sch005 550
Scheme 5. Reaction mechanism and chemoselectivity for Ni-QueD suggested by Liao et al. Adapted from [107], with permission from the PCCP Owner Societies.
Scheme 5. Reaction mechanism and chemoselectivity for Ni-QueD suggested by Liao et al. Adapted from [107], with permission from the PCCP Owner Societies.
Inorganics 07 00095 sch005
Inorganics 07 00095 sch006 550
Scheme 6. The bridging hydroxide attack mechanism suggested by Merz et al. on the basis of DFT calculations. Adapted from [120], Copyright © 2003, American Chemical Society. Energies are given in kcal/mol relative to S3.
Scheme 6. The bridging hydroxide attack mechanism suggested by Merz et al. on the basis of DFT calculations. Adapted from [120], Copyright © 2003, American Chemical Society. Energies are given in kcal/mol relative to S3.
Inorganics 07 00095 sch006
Inorganics 07 00095 sch007 550
Scheme 7. The nickel bound water nucleophilic attack mechanism considered by Merz et al. on the basis of DFT calculations [120], Copyright © 2003, American Chemical Society. Energies are given in kcal/mol relative to S3 in Scheme 6.
Scheme 7. The nickel bound water nucleophilic attack mechanism considered by Merz et al. on the basis of DFT calculations [120], Copyright © 2003, American Chemical Society. Energies are given in kcal/mol relative to S3 in Scheme 6.
Inorganics 07 00095 sch007
Inorganics 07 00095 sch008 550
Scheme 8. The elimination mechanism of urea considered by Merz et al. on the basis of DFT calculations. Adapted from [122], Copyright © 2007, American Chemical Society.
Scheme 8. The elimination mechanism of urea considered by Merz et al. on the basis of DFT calculations. Adapted from [122], Copyright © 2007, American Chemical Society.
Inorganics 07 00095 sch008
Inorganics 07 00095 g010 550
Figure 10. A modified proton-coupled hydride transfer mechanism for lactate racemase (LarA) proposed by Chen et al. Adapted from [129], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim. Tyr294 and Lys298 play crucial roles in orienting substrates and stabilizing the negative charge developing at the substrate hydroxyl oxygen in the transition states, thus lowering reaction barriers significantly.
Figure 10. A modified proton-coupled hydride transfer mechanism for lactate racemase (LarA) proposed by Chen et al. Adapted from [129], copyright © WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim. Tyr294 and Lys298 play crucial roles in orienting substrates and stabilizing the negative charge developing at the substrate hydroxyl oxygen in the transition states, thus lowering reaction barriers significantly.
Inorganics 07 00095 g010
Inorganics 07 00095 g011 550
Figure 11. Reaction mechanism for superoxide dismutase (Ni-SOD) suggested by Pelmenschikov et al. Adapted from [134], Copyright © 2006, American Chemical Society.
Figure 11. Reaction mechanism for superoxide dismutase (Ni-SOD) suggested by Pelmenschikov et al. Adapted from [134], Copyright © 2006, American Chemical Society.
Inorganics 07 00095 g011
Table
Table 1. Reactions catalyzed by Ni-dependent enzymes.
Table 1. Reactions catalyzed by Ni-dependent enzymes.
Enzyme NameReaction
Hydrogenase Inorganics 07 00095 i001
Methyl-coenzyme M reductase Inorganics 07 00095 i002
CO dehydrogenase Inorganics 07 00095 i003
Acetyl-CoA synthase Inorganics 07 00095 i004
Acireductone dioxygenase Inorganics 07 00095 i005
Quercetin 2,4-dioxygenase Inorganics 07 00095 i006
Superoxide dismutase Inorganics 07 00095 i007
Urease Inorganics 07 00095 i008
Glyoxylase I Inorganics 07 00095 i009
Lactate racemase Inorganics 07 00095 i010

Share and Cite

MDPI and ACS Style

Siegbahn, P.E.M.; Chen, S.-L.; Liao, R.-Z. Theoretical Studies of Nickel-Dependent Enzymes. Inorganics 2019, 7, 95. https://doi.org/10.3390/inorganics7080095

AMA Style

Siegbahn PEM, Chen S-L, Liao R-Z. Theoretical Studies of Nickel-Dependent Enzymes. Inorganics. 2019; 7(8):95. https://doi.org/10.3390/inorganics7080095

Chicago/Turabian Style

Siegbahn, Per E. M., Shi-Lu Chen, and Rong-Zhen Liao. 2019. "Theoretical Studies of Nickel-Dependent Enzymes" Inorganics 7, no. 8: 95. https://doi.org/10.3390/inorganics7080095

APA Style

Siegbahn, P. E. M., Chen, S. -L., & Liao, R. -Z. (2019). Theoretical Studies of Nickel-Dependent Enzymes. Inorganics, 7(8), 95. https://doi.org/10.3390/inorganics7080095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop