Disclosure of Key Stereoelectronic Factors for Efficient H2 Binding and Cleavage in the Active Site of [NiFe]-Hydrogenases

Expanding the scope of resonance Raman spectroscopy in hydrogenase research: New observable states and reporter vibrations

Oxygen-tolerant [NiFe] hydrogenases are valuable blueprints for the activation and evolution of molecular hydrogen under application-relevant conditions. Vibrational spectroscopic techniques play a key role in the investigation of these metalloenzymes. For instance, resonance Raman spectroscopy has been introduced as a site-selective approach for probing metal-ligand coordinates of the [NiFe] active site and FeS clusters. Despite its success, this approach is still challenged by a limited number of detectable active-site states - due to missing resonance enhancement or intrinsic light sensitivity - and difficulties in their assignment. Utilizing two oxygen-tolerant [NiFe] hydrogenases as model systems, we illustrate how these challenges can be met by extending excitation and detection wavelength regimes in resonance Raman spectroscopic studies. Specifically, we observe that this technique does not only probe low-frequency metal-ligand vibrations but also high-frequency intra-ligand modes of the diatomic CO/CN- ligands at the active site of [NiFe] hydrogenases. These reporter vibrations are routinely probed by infrared absorption spectroscopy, so that direct comparison of spectra from both techniques allows an unambiguous assignment of states detected by resonance Raman spectroscopy. Moreover, we find that a previously undetected state featuring a bridging hydroxo ligand between Ni and Fe can be probed using higher excitation wavelengths, as photoconversion occurring at lower wavelengths is avoided. In summary, this study expands the applicability of resonance Raman spectroscopy to hydrogenases and other complex metalloenzymes by introducing new strategies for probing and assigning redox-structural states of the active site.

Understanding the [NiFe] Hydrogenase Active Site Environment through Ultrafast Infrared and 2D-IR Spectroscopy of the Subsite Analogue K[CpFe(CO)(CN)2] in Polar and Protic Solvents

The [CpFe(CO)(CN)2]- unit is an excellent structural model for the Fe(CO)(CN)2 moiety of the active site found in [NiFe] hydrogenases. Ultrafast infrared (IR) pump-probe and 2D-IR spectroscopy have been used to study K[CpFe(CO)(CN)2] (M1) in a range of protic and polar solvents and as a dry film. Measurements of anharmonicity, intermode vibrational coupling strength, vibrational relaxation time, and solvation dynamics of the CO and CN stretching modes of M1 in H2O, D2O, methanol, dimethyl sulfoxide, and acetonitrile reveal that H-bonding to the CN ligands plays an important role in defining the spectroscopic characteristics and relaxation dynamics of the Fe(CO)(CN)2 unit. Comparisons of the spectroscopic and dynamic data obtained for M1 in solution and in a dry film with those obtained for the enzyme led to the conclusion that the protein backbone forms an important part of the bimetallic active site environment via secondary coordination sphere interactions.

Stepwise conversion of the Cys6[4Fe-3S] to a Cys4[4Fe-4S] cluster and its impact on the oxygen tolerance of [NiFe]-hydrogenase

The membrane-bound [NiFe]-hydrogenase of Cupriavidus necator is a rare example of a truly O2-tolerant hydrogenase. It catalyzes the oxidation of H2 into 2e- and 2H+ in the presence of high O2 concentrations. This characteristic trait is intimately linked to the unique Cys6[4Fe-3S] cluster located in the proximal position to the catalytic center and coordinated by six cysteine residues. Two of these cysteines play an essential role in redox-dependent cluster plasticity, which bestows the cofactor with the capacity to mediate two redox transitions at physiological potentials. Here, we investigated the individual roles of the two additional cysteines by replacing them individually as well as simultaneously with glycine. The crystal structures of the corresponding MBH variants revealed the presence of Cys5[4Fe-4S] or Cys4[4Fe-4S] clusters of different architecture. The protein X-ray crystallography results were correlated with accompanying biochemical, spectroscopic and electrochemical data. The exchanges resulted in a diminished O2 tolerance of all MBH variants, which was attributed to the fact that the modified proximal clusters mediated only one redox transition. The previously proposed O2 protection mechanism that detoxifies O2 to H2O using four protons and four electrons supplied by the cofactor infrastructure, is extended by our results, which suggest efficient shutdown of enzyme function by formation of a hydroxy ligand in the active site that protects the enzyme from O2 binding under electron-deficient conditions.

Reversible Glutamate Coordination to High-Valent Nickel Protects the Active Site of a [NiFe] Hydrogenase from Oxygen

NAD⁺-reducing [NiFe] hydrogenases are valuable biocatalysts for H₂-based energy conversion and the regeneration of nucleotide cofactors. While most hydrogenases are sensitive toward O₂ and elevated temperatures, the soluble NAD⁺-reducing [NiFe] hydrogenase from Hydrogenophilus thermoluteolus (HtSH) is O₂-tolerant and thermostable. Thus, it represents a promising candidate for biotechnological applications. Here, we have investigated the catalytic activity and active-site structure of native HtSH and variants in which a glutamate residue in the active-site cavity was replaced by glutamine, alanine, and aspartate. Our biochemical, spectroscopic, and theoretical studies reveal that at least two active-site states of oxidized HtSH feature an unusual architecture in which the glutamate acts as a terminal ligand of the active-site nickel. This observation demonstrates that crystallographically observed glutamate coordination represents a native feature of the enzyme. One of these states is diamagnetic and characterized by a very high stretching frequency of an iron-bound active-site CO ligand. Supported by density-functional-theory calculations, we identify this state as a high-valent species with a biologically unprecedented formal Ni(IV) ground state. Detailed insights into its structure and dynamics were obtained by ultrafast and two-dimensional infrared spectroscopy, demonstrating that it represents a conformationally strained state with unusual bond properties. Our data further show that this state is selectively and reversibly formed under oxic conditions, especially upon rapid exposure to high O₂ levels. We conclude that the kinetically controlled formation of this six-coordinate high-valent state represents a specific and precisely orchestrated stereoelectronic response toward O₂ that could protect the enzyme from oxidative damage.

Understanding 2D-IR Spectra of Hydrogenases: A Descriptive and Predictive Computational Study

[NiFe] hydrogenases are metalloenzymes that catalyze the reversible cleavage of dihydrogen (), a clean future fuel. Understanding the mechanism of these biocatalysts requires spectroscopic techniques that yield insights into the structure and dynamics of the [NiFe] active site. Due to the presence of CO and ligands at this cofactor, infrared (IR) spectroscopy represents an ideal technique for studying these aspects, but molecular information from linear IR absorption experiments is limited. More detailed insights can be obtained from ultrafast nonlinear IR techniques like IRpump−IRprobe and two-dimensional (2D-)IR spectroscopy. However, fully exploiting these advanced techniques requires an in-depth understanding of experimental observables and the encoded molecular information. To address this challenge, we present a descriptive and predictive computational approach for the simulation and analysis of static 2D-IR spectra of [NiFe] hydrogenases and similar organometallic systems. Accurate reproduction of experimental spectra from a first-coordination-sphere model suggests a decisive role of the [NiFe] core in shaping the enzymatic potential energy surface. We also reveal spectrally encoded molecular information that is not accessible by experiments, thereby helping to understand the catalytic role of the diatomic ligands, structural differences between [NiFe] intermediates, and possible energy transfer mechanisms. Our studies demonstrate the feasibility and benefits of computational spectroscopy in the 2D-IR investigation of hydrogenases, thereby further strengthening the potential of this nonlinear IR technique as a powerful research tool for the investigation of complex bioinorganic molecules.

Second and Outer Coordination Sphere Effects in Nitrogenase, Hydrogenase, Formate Dehydrogenase, and CO Dehydrogenase

Gases like H2, N2, CO2, and CO are increasingly recognized as critical feedstock in "green" energy conversion and as sources of nitrogen and carbon for the agricultural and chemical sectors. However, the industrial transformation of N2, CO2, and CO and the production of H2 require significant energy input, which renders processes like steam reforming and the Haber-Bosch reaction economically and environmentally unviable. Nature, on the other hand, performs similar tasks efficiently at ambient temperature and pressure, exploiting gas-processing metalloenzymes (GPMs) that bind low-valent metal cofactors based on iron, nickel, molybdenum, tungsten, and sulfur. Such systems are studied to understand the biocatalytic principles of gas conversion including N2 fixation by nitrogenase and H2 production by hydrogenase as well as CO2 and CO conversion by formate dehydrogenase, carbon monoxide dehydrogenase, and nitrogenase. In this review, we emphasize the importance of the cofactor/protein interface, discussing how second and outer coordination sphere effects determine, modulate, and optimize the catalytic activity of GPMs. These may comprise ionic interactions in the second coordination sphere that shape the electron density distribution across the cofactor, hydrogen bonding changes, and allosteric effects. In the outer coordination sphere, proton transfer and electron transfer are discussed, alongside the role of hydrophobic substrate channels and protein structural changes. Combining the information gained from structural biology, enzyme kinetics, and various spectroscopic techniques, we aim toward a comprehensive understanding of catalysis beyond the first coordination sphere.

Bioinorganic Chemistry on Electrodes: Methods to Functional Modeling

One of the major goals of bioinorganic chemistry has been to mimic the function of elegant metalloenzymes. Such functional modeling has been difficult to attain in solution, in particular, for reactions that require multiple protons and multiple electrons (nH⁺/ne^-). Using a combination of heterogeneous electrochemistry, electrode and molecule design one may control both electron transfer (ET) and proton transfer (PT) of these nH⁺/ne^- reactions. Such control can allow functional modeling of hydrogenases (H⁺ + e^- → 1/2 H₂), cytochrome c oxidase (O₂ + 4 e^- + 4 H⁺ → 2 H₂O), monooxygenases (RR'CH₂ + O₂ + 2 e^- + 2 H⁺ → RR'CHOH + H₂O) and dioxygenases (S + O₂ → SO₂; S = organic substrate) in aqueous medium and at room temperatures. In addition, these heterogeneous constructs allow probing unnatural bioinspired reactions and estimation of the inner- and outer-sphere reorganization energy of small molecules and proteins.

Hydrogen evolution, electron-transfer, and hydride-transfer reactions in a nickel-iron hydrogenase model complex: a theoretical study of the distinctive reactivities for the conformational isomers of nickel-iron hydride

Hydrogen fuel is a promising alternative to fossil fuel. Therefore, efficient hydrogen production is crucial to elucidate the distinctive reactivities of metal hydride species, the intermediates formed during hydrogen activation/evolution in the presence of organometallic catalysts. This study uses density functional theory (DFT) to investigate the isomerizations and reactivities of three nickel-iron (NiFe) hydride isomers synthesized by mimicking the active center of NiFe hydrogenase. Hydride transfer within these complexes, rather than a chemical reaction between the complexes, converts the three hydrides internally. Their reactivities, including their electron-transfer, hydride-transfer and proton-transfer reactions, are investigated. The bridging hydride complex exhibits a higher energy level for the highest occupied molecular orbital (HOMO) than the terminal hydride during the electron-transfer reaction. This energy level indicates that the bridging hydride is more easily oxidized and is more susceptible to electron transfer than the terminal hydride. Regarding the hydride-transfer reaction between the NiFe hydride complex and methylene blue, the terminal hydrides exhibit larger hydricity and lower reaction barriers than the bridging hydride complexes. The results of energy decomposition analysis indicate that the structural deformation energy of the terminal hydride in the transition state is smaller than that of the bridging hydride complex, which lowers the reaction barrier of hydride transfer in the terminal hydride. To produce hydrogen, the rate-determining step is represented by the protonation of the hydride, and the terminal hydrides are thermodynamically and kinetically superior to the bridging ones. The differences in the reactivities of the hydride isomers ensure the precise control of hydrogen, and the theoretical calculations can be applied to design catalysts for hydrogen activation/production.

EDA-NOCV Calculation for Efficient N2 Binding to the Reduced Ni3S8 Complex: Estimation of Ni-N2 Intrinsic Interaction Energies

The binding of the dinitrogen molecule to the metal center is the first and crucial step toward dinitrogen activation. Favorable interaction energies are desired by chemists and biochemists to study model complexes in the laboratory. An electrochemically reduced form of a previously isolated sulfur-bridged Ni₃S₈ complex is inferred to bind N₂ at multiple Ni centers, and this bonded N₂ undergoes reductive protonation to produce hydrazine (N₂H₄) as the product in the presence of a proton donor. Density functional theory (DFT) calculations and quantum theory of atoms in molecules (QTAIM) analysis have been carried out to shed light on the nature of N₂ binding to an anionic trinuclear Ni₃S₈ complex. Additionally, energy decomposition analysis with the combination of natural orbital for chemical valence (EDA-NOCV) analysis has been performed to estimate the pairwise interaction energies between the Ni center and the N₂ molecule under experimental conditions.

Computational Survey of Recent Experimental Developments in the Hydroxylation Mechanism of Kynurenine 3-Monooxygenase

Recently, two new mechanistic proposals for the kynurenine 3-monooxygenase (KMO) catalyzed hydroxylation reaction of l-Kynurenine (l-Kyn) have been proposed. According to the first proposal, instead of the distal oxygen, the proximal oxygen of the hydroperoxide intermediate of flavin adenine dinucleotide (FAD) is transferred to the substrate ring. The second study proposes that l-Kyn participates in its base form in the reaction. To address these proposals, the reaction was reconsidered with a 386 atom quantum cluster model that is based on a recent X-ray structure (PDB id: 6FOX). The computations were carried out at the UB3LYP/6-311+G(2d,2p)//UB3LYP/6-31G(d,p) level with solvation (polarizable continuum model) and dispersion (DFT-D3(BJ)) corrections. To supplement the results of the density functional theory (DFT) calculations, molecular dynamics (MD) simulations of the protein-substrate complex were employed. The comparison of a proximal oxygen transfer mechanism to the distal oxygen transfer mechanism revealed that the former requires too high of a barrier energy while the latter validated our previous results. According to the MD simulations, the hydroperoxy moiety does not favor an alignment that might promote the proximal oxygen transfer mechanism. In the second part of the study, hydroxylation reaction with the base form of l-Kyn was sought. Although DFT calculations confirmed a much more facile reaction with the base form of l-Kyn, a mechanism which would allow the deprotonation of the l-Kyn before the oxygen transfer could not be determined with the X-ray-based positions. A concerted mechanism with both the oxygen transfer and the deprotonation required a high barrier energy. A stepwise mechanism involving the deprotonation of l-Kyn was found, starting from an MD frame. The overall barrier of the oxygen transfer step of this model was found to be in the range of that of with neutral l-Kyn. MD simulations supported the idea of ineffectiveness of the nearby shell surrounding the utilized active site core on the deprotonation of l-Kyn.

H2 activation by hydrogenase-inspired NiFe catalyst using frustrated Lewis pair: effect of buffer and halide ion in the heterolytic H-H bond cleavage

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production. We investigate the reaction mechanism of H₂ activation in an aqueous solution by the recently developed NiFe complex (Ogo et al. Sci. Adv. 2020, 6, eaaz8181) using density functional theory (DFT) calculation. Our computational results showed that H₂ is activated using frustrated Lewis pair. That is, H₂ binds to the Fe site of the NiFe complex, acting as a Lewis acid, while the added buffer, acting as Lewis base, abstracts protons to form a hydride complex. Furthermore, the higher basicity in the proton abstraction reaction characterises reaction more exergonic and lowers the reaction barrier. In addition, in the proton abstraction by the water molecule, the reaction barrier was lowered when anion such as Cl⁻ is in the vicinity of the water. Understanding the chemical species that contribute to the catalytic process in cooperation with the metal catalyst at the atomic level should help to maximise the function of the catalyst.

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production.

Electrocatalytic Reduction of Nitrogen to Hydrazine Using a Trinuclear Nickel Complex

Activation and reduction of N₂ have been a major challenge to chemists and the focus since now has mostly been on the synthesis of NH₃. Alternatively, reduction of N₂ to hydrazine is desirable because hydrazine is an excellent energy vector that can release the stored energy very conveniently without the need for catalysts. To date, only one molecular catalyst has been reported to be able to reduce N₂ to hydrazine chemically. A trinuclear T-shaped nickel thiolate molecular complex has been designed to activate dinitrogen. The electrochemically generated all Ni(I) state of this molecule can reduce N₂ in the presence of PhOH as a proton donor. Hydrazine is detected as the only nitrogen-containing product of the reaction, along with gaseous H₂. The complex reported here is selective for the 4e^-/4H⁺ reduction of nitrogen to hydrazine with a minor overpotential of ∼300 mV.

[NiFe], [FeFe], and [Fe] hydrogenase models from isomers

The study of hydrogenase enzymes (H₂ases) is necessary because of their importance to a future hydrogen energy economy. These enzymes come in three distinct classes: [NiFe] H₂ases, which have a propensity toward H₂ oxidation; [FeFe] H₂ases, which have a propensity toward H₂ evolution; and [Fe] H₂ases, which catalyze H^- transfer. Modeling these enzymes has so far treated them as different species, which is understandable given the different cores and ligand sets of the natural molecules. Here, we demonstrate, using x-ray analysis and nuclear magnetic resonance, infrared, Mössbauer spectroscopies, and electrochemical measurement, that the catalytic properties of all three enzymes can be mimicked with only three isomers of the same NiFe complex.

The large subunit of the regulatory [NiFe]-hydrogenase from Ralstonia eutropha - a minimal hydrogenase?

Chemically synthesized compounds that are capable of facilitating the reversible splitting of dihydrogen into protons and electrons are rare in chemists' portfolio. The corresponding biocatalysts – hydrogenases – are, however, abundant in the microbial world. [NiFe]-hydrogenases represent a major subclass and display a bipartite architecture, composed of a large subunit, hosting the catalytic NiFe(CO)(CN)₂ cofactor, and a small subunit whose iron–sulfur clusters are responsible for electron transfer. To analyze in detail the catalytic competence of the large subunit without its smaller counterpart, we purified the large subunit HoxC of the regulatory [NiFe]-hydrogenase of the model H₂ oxidizer Ralstonia eutropha to homogeneity. Metal determination and infrared spectroscopy revealed a stoichiometric loading of the metal cofactor. This enabled for the first time the determination of the UV-visible extinction coefficient of the NiFe(CO)(CN)₂ cofactor. Moreover, the absence of disturbing iron–sulfur clusters allowed an unbiased look into the low-spin Fe²⁺ of the active site by Mössbauer spectroscopy. Isolated HoxC was active in catalytic hydrogen–deuterium exchange, demonstrating its capacity to activate H₂. Its catalytic activity was drastically lower than that of the bipartite holoenzyme. This was consistent with infrared and electron paramagnetic resonance spectroscopic observations, suggesting that the bridging position between the active site nickel and iron ions is predominantly occupied by water-derived ligands, even under reducing conditions. In fact, the presence of water-derived ligands bound to low-spin Ni²⁺ was reflected by the absorption bands occurring in the corresponding UV-vis spectra, as revealed by time-dependent density functional theory calculations conducted on appropriate in silico models. Thus, the isolated large subunits indeed represent simple [NiFe]-hydrogenase models, which could serve as blueprints for chemically synthesized mimics. Furthermore, our data point to a fundamental role of the small subunit in preventing water access to the catalytic center, which significantly increases the H₂ splitting capacity of the enzyme.

Spectroscopic investigation of an isolated [NiFe]-hydrogenase large subunit enables a unique view of the NiFe(CO)(CN)₂ cofactor.

Selective Oxidation of H2 and CO by NiIr Catalyst in Aqueous Solution: A DFT Mechanistic Study

One of the challenges in utilizing hydrogen gas (H₂) as a sustainable fossil fuel alternative is the inhibition of H₂ oxidation by carbon monoxide (CO), which is involved in the industrial production of H₂ sources. To solve this problem, a catalyst that selectively oxidizes either CO or H₂ or one that co-oxidizes H₂ and CO is needed. Recently, a NiIr catalyst [Ni^IICl(X)Ir^IIICl(η⁵-C₅Me₅)], (X = N,N'-dimethyl-3,7-diazanonane-1,9-dithiolate), which efficiently and selectively oxidizes either H₂ or CO depending on the pH, has been developed (Angew. Chem. Int. Ed. 2017, 56, 9723-9726). In the present work, density functional theory (DFT) calculations are employed to elucidate the pH-dependent reaction mechanisms of H₂ and CO oxidation catalyzed by this NiIr catalyst. During H₂ oxidation, our calculations suggest that dihydrogen binds to the Ir center and generates an Ir(III)-dihydrogen complex, followed by subsequent isomerization to an Ir(V)-dihydride species. Then, a proton is abstracted by a buffer base, CH₃COO^-, resulting in the formation of a hydride complex. The catalytic cycle completes with electron transfer from the hydride complex to a protonated 2,6-dichlorobenzeneindophenol (DCIP) and a proton transfer from the oxidized hydride complex to a buffer base. The CO oxidation mechanism involves three distinct steps, i.e., (1) formation of a metal carbonyl complex, (2) formation of a metallocarboxylic acid, and (3) conversion of the metallocarboxylic acid to a hydride complex. The formation of the metallocarboxylic acid involves nucleophilic attack of OH^- to the carbonyl-C followed by a large structural change with concomitant cleavage of the Ir-S bond and rotation of the COOH group along the NiIr axis. During the conversion of the metallocarboxylic acid to the hydride complex, intramolecular proton transfer followed by removal of CO₂ leads to the formation of the hydride complexes. In addition, the barrier heights for the binding of small molecules (H₂, OH^-, H₂O, and CO) to Ir were calculated, and the results indicated that dissociation from Ir is a faster process than the binding of H₂O and H₂. These calculations indicate that H₂ oxidation is inhibited by CO and OH^- and thus prefers acidic conditions. In contrast, the CO oxidation reactions occur more favorably under basic conditions, as the formation of the metallocarboxylic acid involves OH^- attack to a carbonyl-C and the binding of OH^- to Ni largely stabilizes the triplet spin state of the complex. Taken together, these calculations provide a rationale for the experimentally observed pH-dependent, selective oxidations of H₂ and CO.

Molecular Electrocatalysts for the Hydrogen Evolution Reaction: Input from Quantum Chemistry

In the pursuit of carbon-free fuels, hydrogen can be considered as an apt energy carrier. The design of molecular electrocatalysts for hydrogen production is important for the development of renewable energy sources that are abundant, inexpensive, and environmentally benign. Over the last 20 years, a large number of electrocatalysts have been developed, and considerable efforts have been directed toward the design of earth-abundant, first-row transition-metal complexes capable of promoting electrocatalytic hydrogen evolution reaction (HER). In this context, numerical approaches have emerged as powerful tools to study the catalytic performances of these complexes. This review covers some of the most significant theoretical mechanistic studies of biomimetic and bioinspired homogeneous HER catalysts. The approaches employed to study the free energy landscapes are discussed and methods used to obtain accurate estimates of relevant observables required to study the HER are presented. Furthermore, the structural and electronic parameters that govern the reactivity, and are necessary to achieve efficient hydrogen production, are discussed in view of future research directions.

X-ray Crystallography and Vibrational Spectroscopy Reveal the Key Determinants of Biocatalytic Dihydrogen Cycling by [NiFe] Hydrogenases

[NiFe] hydrogenases are complex model enzymes for the reversible cleavage of dihydrogen (H₂). However, structural determinants of efficient H₂ binding to their [NiFe] active site are not properly understood. Here, we present crystallographic and vibrational‐spectroscopic insights into the unexplored structure of the H₂‐binding [NiFe] intermediate. Using an F₄₂₀‐reducing [NiFe]‐hydrogenase from Methanosarcina barkeri as a model enzyme, we show that the protein backbone provides a strained chelating scaffold that tunes the [NiFe] active site for efficient H₂ binding and conversion. The protein matrix also directs H₂ diffusion to the [NiFe] site via two gas channels and allows the distribution of electrons between functional protomers through a subunit‐bridging FeS cluster. Our findings emphasize the relevance of an atypical Ni coordination, thereby providing a blueprint for the design of bio‐inspired H₂‐conversion catalysts.

Understanding the structure and dynamics of hydrogenases by ultrafast and two-dimensional infrared spectroscopy

Hydrogenases are valuable model enzymes for sustainable energy conversion approaches using H2, but rational utilization of these base-metal biocatalysts requires a detailed understanding of the structure and dynamics of their complex active sites. The intrinsic CO and CN- ligands of these metalloenzymes represent ideal chromophores for infrared (IR) spectroscopy, but structural and dynamic insight from conventional IR absorption experiments is limited. Here, we apply ultrafast and two-dimensional (2D) IR spectroscopic techniques, for the first time, to study hydrogenases in detail. Using an O2-tolerant [NiFe] hydrogenase as a model system, we demonstrate that IR pump-probe spectroscopy can explore catalytically relevant ligand bonding by accessing high-lying vibrational states. This ultrafast technique also shows that the protein matrix is influential in vibrational relaxation, which may be relevant for energy dissipation from the active site during fast reaction steps. Further insights into the relevance of the active site environment are provided by 2D-IR spectroscopy, which reveals equilibrium dynamics and structural constraints imposed on the H2-accepting intermediate of [NiFe] hydrogenases. Both techniques offer new strategies for uniquely identifying redox-structural states in complex catalytic mixtures via vibrational quantum beats and 2D-IR off-diagonal peaks. Together, these findings considerably expand the scope of IR spectroscopy in hydrogenase research, and new perspectives for the characterization of these enzymes and other (bio-)organometallic targets are presented.

Mechanism of Kynurenine 3-Monooxygenase-Catalyzed Hydroxylation Reaction: A Quantum Cluster Approach

The mechanism of the hydroxylation reaction between l-Kyn and model flavin adenine dinucleotide (FAD)-hydroperoxide was investigated via density functional theory (DFT) calculations in the absence and in the presence of the kynurenine 3-monooxygenase (KMO) enzyme by considering possible pathways that can lead to the product 3-hydroxykynurenine (3-HK). Crystal structure (pdb code: 5NAK )-based calculations involved a quantum cluster model in which the active site of the enzyme with the substrate l-Kyn was represented with 348 atoms. According to the deduced mechanism, KMO-catalyzed hydroxylation reaction takes place with four transformations. In the initial transition state, FAD delivers its peroxy hydroxyl to the l-Kyn ring, creating an sp³-hybridized carbon center. Then, the hydrogen on the hydroxyl moiety is immediately transferred back to the proximal oxygen that remained on FAD. These consequent transformations are in line with the somersault rearrangement previously described for similar enzymatic systems. The second step corresponds to a hydride shift from the sp³-hybridized carbon of the substrate ring to its adjacent carbon, producing the keto form of 3-HK. Then, keto-3-HK is transformed into its enol form (3-HK) with a water-assisted tautomerization. Lastly, FAD is oxidized with a water-assisted dehydration, which also involves 3-HK as a catalyst. In the proposed pathway, Asn54, Pro318, and a crystal water molecule were seen to play significant roles in the proton relays. The energies obtained via the cluster approach were calculated at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d,p) level with solvation (polarizable continuum model) and dispersion (DFT-D3(BJ)) corrections.

Reactivation of the Ready and Unready Oxidized States of [NiFe]-Hydrogenases: Mechanistic Insights from DFT Calculations

The apparently simple dihydrogen formation from protons and electrons (2H⁺ + 2e^- ⇄ H₂) is one of the most challenging reactions in nature. It is catalyzed by metalloenzymes of amazing complexity, called hydrogenases. A better understanding of the chemistry of these enzymes, especially that of the [NiFe]-hydrogenases subgroup, has important implications for production of H₂ as alternative sustainable fuel. In this work, reactivation mechanism of the oxidized and inactive Ni-B and Ni-A states of the [NiFe]-hydrogenases active site has been investigated using density functional theory. Results obtained from this study show that one-electron reduction and protonation of the active site promote the removal of the bridging hydroxide ligand contained in Ni-B and Ni-A. However, this process is sufficient to activate only the Ni-B state. H₂ binding to the active site is required to convert Ni-A to the active Ni-SI_a state. Here, we also propose a reasonable structure for the spectroscopically well-characterized Ni-SI_r and Ni-SU species, formed respectively from the one-electron reduction of Ni-B and Ni-A. Ni-SI_r, depending on the pH at which the reaction occurs, features a bridging hydroxide ligand or a water molecule terminally coordinated to the Ni atom, whereas in Ni-SU a water molecule is terminally coordinated to the Fe atom, and the Cys64 residue is oxidized to sulfenate. The sulfenate oxygen atom in the Ni-A state affects the stereoelectronic properties of the binuclear cluster by modifying the coordination geometry of Ni, and consequently, by switching the regiochemistry of H₂O and H₂ binding from the Ni to the Fe atom. This effect is predicted to be at the origin of the different reactivation kinetics of the oxidized and inactive Ni-B and Ni-A states.

Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods

[NiFe] hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons. This seemingly simple reaction has attracted much attention because of the prospective use of H₂ as a clean fuel. In this paper, we have studied the full reaction mechanism of this enzyme with various computational methods. Geometries were obtained with combined quantum mechanical and molecular mechanics (QM/MM) calculations. To get more accurate energies and obtain a detailed account of the surroundings, we performed big-QM calculations with 819 atoms in the QM region. Moreover, QM/MM thermodynamic cycle perturbation calculations were performed to obtain free energies. Finally, density matrix renormalisation group complete active space self-consistent field calculations were carried out to study the electronic structures of the various states in the reaction mechanism. Our calculations indicate that the Ni-L state is not involved in the reaction mechanism. Instead, the Ni-C state is reduced by one electron and then the bridging hydride ion is transferred to the sulfur atom of Cys546 as a proton and the two electrons transfer to the Ni ion. This step turned out to be rate-determining with an energy barrier of 58 kJ/mol, which is consistent with the experimental rate of 750 ± 90 s^-1 (corresponding to ∼52 kJ/mol). The cleavage of the H-H bond is facile with an energy barrier of 33 kJ/mol, according to our calculations. We also find that the reaction energies are sensitive to the size of the QM system, the basis set, and the density functional theory method, in agreement with previous studies.

Comprehensive reaction mechanisms at and near the Ni-Fe active sites of [NiFe] hydrogenases

[NiFe] hydrogenase (H2ase) catalyzes the oxidation of dihydrogen to two protons and two electrons and/or its reverse reaction. For this simple reaction, the enzyme has developed a sophisticated but intricate mechanism with heterolytic cleavage of dihydrogen (or a combination of a hydride and a proton), where its Ni-Fe active site exhibits various redox states. Recently, thermodynamic parameters of the acid-base equilibrium for activation-inactivation, a new intermediate in the catalytic reaction, and new crystal structures of [NiFe] H2ases have been reported, providing significant insights into the activation-inactivation and catalytic reaction mechanisms of [NiFe] H2ases. This Perspective provides an overview of the reaction mechanisms of [NiFe] H2ases based on these new findings.

QM/MM Investigation of the Role of a Second Coordination Shell Arginine in [NiFe]-Hydrogenases

[NiFe]-hydrogenases are highly efficient catalysts for the heterolytic splitting of molecular hydrogen (H₂). The heterobimetallic cysteine-coordinated active site of these enzymes is covered by a highly conserved arginine residue, whose role in the reaction is not fully resolved yet. The structural and catalytic role of this arginine is investigated here using QM/MM calculations with various exchange-correlation functionals. All of them give a very consistent picture of the thermodynamics of H₂ oxidation. The concept of the presence of a neutral arginine and its direct involvement as a Frustrated Lewis Pair (FLP) in the reaction is critically evaluated. The arginine, however, would exist in its standard protonation state and perform a critical role in positioning and slightly polarizing the substrate H₂. It is not directly involved in the heterolytic processing of H₂ but guides its approach and reduces its flexibility during binding. Upon substitution of the positively charged arginine by a charge-conserving lysine residue, the H₂ binding position remains unaffected. However, critical hydrogen bonding interactions with nearby aspartate residues are lost. In addition, the H₂ polarization is unfavorable and the reduced side-chain volume may negatively affect the kinetics of the catalytic process.

Hydrogen bonding effect between active site and protein environment on catalysis performance in H2-producing [NiFe] hydrogenases

The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol^-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H⁺/e^- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S^- in cysteine residues shows a larger influence on the energy profile than that via the CN^- ligands on Fe.

The oxygen reduction reaction on [NiFe] hydrogenases

Oxygen tolerance capacity is critical for hydrogen oxidation/evolution catalysts.

Jensen HJ, Hedegård ED. Exploration of H2 binding to the [NiFe]-hydrogenase active site with multiconfigurational density functional theory

The combination of density functional theory (DFT) with a multiconfigurational wave function is an efficient way to include dynamical correlation in calculations with multiconfiguration self-consistent field wave functions.

Theoretical investigation of aerobic and anaerobic oxidative inactivation of the [NiFe]-hydrogenase active site

The extraordinary capability of [NiFe]-hydrogenases to catalyse the reversible interconversion of protons and electrons into dihydrogen (H₂) has stimulated numerous experimental and theoretical studies addressing the direct utilization of these enzymes in H₂ production processes.

Catalytic Mechanism of Fungal Lytic Polysaccharide Monooxygenases Investigated by First-Principles Calculations

Lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes that facilitate the degradation of recalcitrant polysaccharides by the oxidative cleavage of glycosidic bonds. They are gaining rapidly increasing attention as key players in biomass conversion, especially for the production of second-generation biofuels. Elucidation of the detailed mechanism of the LPMO reaction is a major step toward the assessment and optimization of LPMO efficacy in industrial biotechnology, paving the way to utilization of sustainable fuel sources. Here, we used density functional theory calculations to study the reaction pathways suggested to date, exploiting a very large active-site model for a fungal AA9 LPMO and using a celloheptaose unit as a substrate mimic. We identify a copper oxyl intermediate as being responsible for H-atom abstraction from the substrate, followed by a rapid, water-assisted hydroxyl rebound, leading to substrate hydroxylation.

H2binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods

CCSD(T) and DMRG-CASPT2 calculations show that H₂prefers to bind to Ni rather than to Fe in [NiFe] hydrogenase.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

Structure and function of [NiFe] hydrogenases

Hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons via a heterolytic splitting mechanism. The active sites of [NiFe] hydrogenases comprise a dinuclear Ni-Fe center carrying CO and CN^- ligands. The catalytic activity of the standard (O₂-sensitive) [NiFe] hydrogenases vanishes under aerobic conditions. The O₂-tolerant [NiFe] hydrogenases can sustain H₂ oxidation activity under atmospheric conditions. These hydrogenases have very similar active site structures that change the ligand sphere during the activation/catalytic process. An important structural difference between these hydrogenases has been found for the proximal iron-sulphur cluster located in the vicinity of the active site. This unprecedented [4Fe-3S]-6Cys cluster can supply two electrons, which lead to rapid recovery of the O₂ inactivation, to the [NiFe] active site.

The working mechanism of the β-carbonic anhydrase degrading carbonyl sulphide (COSase): a theoretical study

In order to give insights into the working mechanism of the novel characterized enzyme carbonyl sulphide hydrolase (COSase), which efficiently converts COS into H2S and CO2, we have performed a detailed theoretical investigation using the framework of density functional theory (using B3LYP and M06 exchange-correlation functionals) by the cluster model approach. In the final part of the reaction the metal ion is unable to form a pentacoordinated species. The B3LYP-D3 and M06 potential energy surfaces have a very similar shape. The elucidation of the catalytic reduction of COS is important in view of its role in environmental chemistry.

A threonine stabilizes the NiC and NiR catalytic intermediates of [NiFe]-hydrogenase

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

Concepts in bio-molecular spectroscopy: vibrational case studies on metalloenzymes

Challenges and chances in bio-molecular spectroscopy are exemplified by vibrational case studies on metalloenzymes.

Stereoelectronic effects: a simple yet powerful tool to manipulate anion affinity

Stereoelectronic effects on anion binding were examined, IR spectroscopy was used to probe structures, and a well aligned non-interacting group can be more significant than a hydrogen bond donor.