Computational studies of [NiFe] and [FeFe] hydrogenases

Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H2 gas powered by a [NiFe] active site and transduces the free energy into proton pumping and Na+/H+ exchange across the membrane. Despite recent structural advances, the mechanistic principles of H2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H2 and how the catalysis couples to conformational dynamics in the membrane domain of Mbh. By combining large-scale quantum chemical density functional theory (DFT) and correlated ab initio wave function methods with atomistic molecular dynamics simulations, we show that the proton transfer reactions required for the catalysis are gated by electric field effects that direct the protons by water-mediated reactions from Glu21L toward the [NiFe] site, or alternatively along the nearby His75L pathway that also becomes energetically feasible in certain reaction steps. These local proton-coupled electron transfer (PCET) reactions induce conformational changes around the active site that provide a key coupling element via conserved loop structures to the ion transport activity. We find that H2 forms in a heterolytic proton reduction step, with spin crossovers tuning the energetics along key reaction steps. On a general level, our work showcases the role of electric fields in enzyme catalysis and how these effects are employed by the [NiFe] active site of Mbh to drive PCET reactions and ion transport.

Perspective: multi-configurational methods in bio-inorganic chemistry

Transition metal ions play crucial roles in the structure and function of numerous proteins, contributing to essential biological processes such as catalysis, electron transfer, and oxygen binding. However, accurately modeling the electronic structure and properties of metalloproteins poses significant challenges due to the complex nature of their electronic configurations and strong correlation effects. Multiconfigurational quantum chemistry methods are, in principle, the most appropriate tools for addressing these challenges, offering the capability to capture the inherent multi-reference character and strong electron correlation present in bio-inorganic systems. Yet their computational cost has long hindered wider adoption, making methods such as density functional theory (DFT) the method of choice. However, advancements over the past decade have substantially alleviated this limitation, rendering multiconfigurational quantum chemistry methods more accessible and applicable to a wider range of bio-inorganic systems. In this perspective, we discuss some of these developments and how they have already been used to answer some of the most important questions in bio-inorganic chemistry. We also comment on ongoing developments in the field and how the future of the field may evolve.

Molecular Catalysis of Energy Relevance in Metal-Organic Frameworks: From Higher Coordination Sphere to System Effects

The modularity and synthetic flexibility of metal-organic frameworks (MOFs) have provoked analogies with enzymes, and even the term MOFzymes has been coined. In this review, we focus on molecular catalysis of energy relevance in MOFs, more specifically water oxidation, oxygen and carbon dioxide reduction, as well as hydrogen evolution in context of the MOF-enzyme analogy. Similar to enzymes, catalyst encapsulation in MOFs leads to structural stabilization under turnover conditions, while catalyst motifs that are synthetically out of reach in a homogeneous solution phase may be attainable as secondary building units in MOFs. Exploring the unique synthetic possibilities in MOFs, specific groups in the second and third coordination sphere around the catalytic active site have been incorporated to facilitate catalysis. A key difference between enzymes and MOFs is the fact that active site concentrations in the latter are often considerably higher, leading to charge and mass transport limitations in MOFs that are more severe than those in enzymes. High catalyst concentrations also put a limit on the distance between catalysts, and thus the available space for higher coordination sphere engineering. As transport is important for MOF-borne catalysis, a system perspective is chosen to highlight concepts that address the issue. A detailed section on transport and light-driven reactivity sets the stage for a concise review of the currently available literature on utilizing principles from Nature and system design for the preparation of catalytic MOF-based materials.

Modeling of the Electrostatic Interaction and Catalytic Activity of [NiFe] Hydrogenases on a Planar Electrode

Hydrogenases are a group of enzymes that have caught the interest of researchers in renewable energies, due to their ability to catalyze the redox reaction of hydrogen. The exploitation of hydrogenases in electrochemical devices requires their immobilization on the surface of suitable electrodes, such as graphite. The orientation of the enzyme on the electrode is important to ensure a good flux of electrons to the catalytic center, through an array of iron-sulfur clusters. Here we present a computational approach to determine the possible orientations of a [NiFe] hydrogenase (PDB 1e3d) on a planar electrode, as a function of pH, salinity, and electrode potential. The calculations are based on the solution of the linearized Poisson-Boltzmann equation, using the PyGBe software. The results reveal that electrostatic interactions do not truly immobilize the enzyme on the surface of the electrode, but there is instead a dynamic equilibrium between different orientations. Nonetheless, after averaging over all thermally accessible orientations, we find significant differences related to the solution's salinity and pH, while the effect of the electrode potential is relatively weak. We also combine models for the protein adsoption-desorption equilibria and for the electron transfer between the proteins and the electrode to arrive at a prediction of the electrode's activity as a function of the enzyme concentration.

Catalytic Mechanism of Competing Proton Transfer Events from Water and Acetic Acid by [CoII(bpbH2)Cl2] for Water Splitting Processes

We performed first principles simulations to explore the water reduction process of the cobalt complex [Co^II(bpbH₂)Cl₂], where bpbH₂ = N,N'-bis(2'-pyridine carboxamide)-1,2-benzene. We considered the sequence steps of electron reduction followed by the proton addition process to observe the hydrogen evolution process. An experimental study of the catalyst showed that the increase in the acetic acid concentration triggers catalytic current and reduction of Co(II) to Co(I), and protonation occurred, yielding a Co(III)-H intermediate. Therefore, we used water and acetic acid as the proton sources. We compare the proton transfer kinetics from both the water and acetic acid. The reduction potentials and proton transfer kinetics from water or acetic acid to the reaction center were studied in a DMF solvent through the implicit solvent model. The first proton transfer from the acetic acid is more favorable, forming a Co^III-H complex and further reducing to Co^II-H. The second proton transfer from water to the Co^II-H moiety requires less free energy than acetic acid and is the rate-limiting step. The nature of the reduction process is also examined through the charge analysis, which reveals that the ligand becomes softer due to the C═O groups.

Activation of Small Molecules and Hydrogenation of CO2 Catalyzed by Frustrated Lewis Pairs

The chemistry of frustrated Lewis pair (FLP) is widely explored in the activation of small molecules, the hydrogenation of CO2, and unsaturated organic species. A survey of several experimental works on the activation of small molecules by FLPs and the related mechanistic insights into their reactivity from electronic structure theory calculation are provided in the present review, along with the catalytic hydrogenation of CO2. The mechanistic insight into H2 activation is thoroughly discussed, which may provide a guideline to design more efficient FLP for H2 activation. FLPs can activate other small molecules like, CO, NO, CO2, SO2, N2O, alkenes, alkynes, etc. by cooperative action of the Lewis centers of FLPs, as revealed by several computational analyses. The activation barrier of H2 and other small molecules by the FLP can be decreased by utilizing the aromaticity criterion in the FLP as demonstrated by the nucleus independent chemical shift (NICS) analysis. The term boron-ligand cooperation (BLC), which is analogous to the metal-ligand cooperation (MLC), is invoked to describe a distinct class of reactivity of some specific FLPs towards H2 activation.

Role of geochemical protoenzymes (geozymes) in primordial metabolism: specific abiotic hydride transfer by metals to the biological redox cofactor NAD. FEBS J 2021;289:3148-3162. [PMID: 34923745 PMCID: PMC9306933 DOI: 10.1111/febs.16329]

Hydrogen gas, H₂, is generated in serpentinizing hydrothermal systems, where it has supplied electrons and energy for microbial communities since there was liquid water on Earth. In modern metabolism, H₂ is converted by hydrogenases into organically bound hydrides (H^–), for example, the cofactor NADH. It transfers hydrides among molecules, serving as an activated and biologically harnessed form of H₂. In serpentinizing systems, minerals can also bind hydrides and could, in principle, have acted as inorganic hydride donors—possibly as a geochemical protoenzyme, a ‘geozyme’— at the origin of metabolism. To test this idea, we investigated the ability of H₂ to reduce NAD⁺ in the presence of iron (Fe), cobalt (Co) and nickel (Ni), metals that occur in serpentinizing systems. In the presence of H₂, all three metals specifically reduce NAD⁺ to the biologically relevant form, 1,4‐NADH, with up to 100% conversion rates within a few hours under alkaline aqueous conditions at 40 °C. Using Henry's law, the partial pressure of H₂ in our reactions corresponds to 3.6 mm, a concentration observed in many modern serpentinizing systems. While the reduction of NAD⁺ by Ni is strictly H₂‐dependent, experiments in heavy water (²H₂O) indicate that native Fe can reduce NAD⁺ both with and without H₂. The results establish a mechanistic connection between abiotic and biotic hydride donors, indicating that geochemically catalysed, H₂‐dependent NAD⁺ reduction could have preceded the hydrogenase‐dependent reaction in evolution.

Mechanism study of improving anaerobic co-digestion performance of waste activated sludge and food waste by Fe3O4

A large amount of waste activated sludge (WAS) and food waste (FW) are produced every year in China. Anaerobic co-digestion is considered to be an effective way to solve this problem. This study applied FW/WAS mixture as co-substrate to create different digestive environment, aiming to understand the mechanism of Fe₃O₄ particles in promoting AD performance. The results showed that the addition of Fe₃O₄ presented various performances when facing different digestive acidification stress brought by different mixing ratios of WAS and FW. Methanogenic pathways and microbial communities varied with substrates' properties. For group A (WAS mono-digestion), the acetoclastic methanogens dominated, 20 mg/g VS (according to the iron element) Fe₃O₄ could promote methane production, while 200 mg/g VS Fe₃O₄ would inhibit microbial activity. The promoted methane production by Fe₃O₄ was attributable to the promotion of sludge hydrolysis. For group B (WAS: FW = 1:0.5, based on VS addition, similarly hereinafter), Fe₃O₄ triggered direct interspecific electron transfer (DIET) between bacteria and methanogens. For group C (WAS: FW = 1:1), the hydrogenotrophic methanogens dominated, bacteria excreted more non-conductive polysaccharides in EPS to resist unfavorable environment, thereby it prevented their contact with Fe₃O₄ particles. So, it was difficult for Fe₃O₄ to trigger DIET and promote the digestive performance of batch experiments in such condition.

Functional Dynamics of an Ancient Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) is a redox-driven Na⁺/H⁺ transporter that employs the energy from hydrogen gas (H₂) production to catalyze proton pumping and Na⁺/H⁺ exchange across cytoplasmic membranes of archaea. Despite a recently resolved structure of this ancient energy-transducing enzyme [Yu et al. Cell2018, 173, 1636–1649], the molecular principles of its redox-driven ion-transport mechanism remain puzzling and of major interest for understanding bioenergetic principles of early cells. Here we use atomistic molecular dynamics (MD) simulations in combination with data clustering methods and quantum chemical calculations to probe principles underlying proton reduction as well as proton and sodium transport in Mbh from the hyperthermophilic archaeon Pyrococcus furiosus. We identify putative Na⁺ binding sites and proton pathways leading across the membrane and to the NiFe-active center as well as conformational changes that regulate ion uptake. We suggest that Na⁺ binding and protonation changes at a putative ion-binding site couple to proton transfer across the antiporter-like MbhH subunit by modulating the conformational state of a conserved ion pair at the subunit interface. Our findings illustrate conserved coupling principles within the complex I superfamily and provide functional insight into archaeal energy transduction mechanisms.

Photodynamics of Asymmetric Di-Iron-Cyano Hydrogenases Examined by Time-Resolved Mid-Infrared Spectroscopy

Two anionic asymmetric Fe-Fe hydrogenase model compounds containing a single cyano (CN) and five carboxyl (CO) ligands, [Et₄N][Fe₂(μ-S₂C₃H₆)(CO)₅(CN)₁] and [Et₄N][Fe₂(μ-S₂C₂H₄)(CO)₅(CN)₁], dissolved in room-temperature acetonitrile, are examined. The molecular asymmetry affects the redox potentials of the central iron atoms, thus changing the photophysics and possible catalytic properties of the compounds. Femtosecond ultraviolet excitation with mid-infrared probe spectroscopy of the model compounds was employed to better understand the ultrafast dynamics of the enzyme-active site. Continuous ultraviolet lamp excitation with Fourier transform infrared (FTIR) spectroscopy was also used to explore stable product formation on the second timescale. For both model compounds, two timescales are observed; a 20-30 ps decay and the formation of a long-lived photoproduct. The picosecond decay is assigned to vibrational cooling and rotational dynamics, while the residual spectra remain for up to 300 ps, suggesting the formation of new photoproducts. Static FTIR spectroscopy yielded a different stable photoproduct than that observed on the ultrafast timescale. Density functional theory calculations simulated photoproducts for CO-loss and CN-loss isomers, and the resulting photoproduct spectra suggest that the picosecond transients arise from a complex mixture of isomerization after CO-loss, while dimerization and formation of a CN-containing Fe-CO-Fe bridged species are also considered.

Switching Site Reactivity in Hydrogenase Model Systems by Introducing a Pendant Amine Ligand

Hydrogenases are versatile enzymatic catalysts with an unmet hydrogen evolution reactivity (HER) from synthetic bio-inspired systems. The binuclear active site only has one-site reactivity of the distal Fed atom. Here, binuclear complexes [Fe2(CO)5(μ-Mebdt)(P(4-C6H4OCH3)3)] 1 and [Fe2(CO)5(μ-Mebdt)(PPh2Py)] 2 are presented, which show electrocatalytic activity in the presence of weak acids as a proton source for the HER. Despite almost identical structural and spectroscopic properties (bond distances and angles from single-crystal X-ray; IR, UV/vis, and NMR), introduction of a nitrogen base atom in the phosphine ligand in 2 markedly changes site reactivity. The bridging benzenedithiolate ligand Mebdt interacts with the terminal ligand's phenyl aromatic rings and stabilizes the reduced states of the catalysts. Although 1 with monodentate phosphine terminal ligands only shows a distal iron atom HER activity by a sequence of electrochemical and protonation steps, the lone pair of pyridine nitrogen in 2 acts as the primary site of protonation. This swaps the iron atom catalytic activity toward the proximal iron for complex 2. Density-functional theory (DFT) calculations reveal the role of terminal phosphines ligands without/with pendant amines by directing the proton transfer steps. The reactivity of 1 is a thiol-based protonation of a dangling bond in 1- and distal iron hydride mechanism, which may follow either an ECEC or EECC sequence, depending on the choice of acid. The pendant amine in 2 enables a terminal ligand protonation and an ECEC reactivity. The introduction of a terminal nitrogen atom enables the control of site reactivity in a binuclear system.

QM/MM study of the binding of H2 to MoCu CO dehydrogenase: development and applications of improved H2 van der Waals parameters

The MoCu CO dehydrogenase enzyme not only transforms CO into CO₂ but it can also oxidise H₂. Even if its hydrogenase activity has been known for decades, a debate is ongoing on the most plausible mode for the binding of H₂ to the enzyme active site and the hydrogen oxidation mechanism. In the present work, we provide a new perspective on the MoCu-CODH hydrogenase activity by improving the in silico description of the enzyme. Energy refinement—by means of the BigQM approach—was performed on the intermediates involved in the dihydrogen oxidation catalysis reported in our previously published work (Rovaletti, et al. “Theoretical Insights into the Aerobic Hydrogenase Activity of Molybdenum–Copper CO Dehydrogenase.” Inorganics 7 (2019) 135). A suboptimal description of the H₂–HN(backbone) interaction was observed when the van der Waals parameters described in previous literature for H₂ were employed. Therefore, a new set of van der Waals parameters is developed here in order to better describe the hydrogen–backbone interaction. They give rise to improved binding modes of H₂ in the active site of MoCu CO dehydrogenase. Implications of the resulting outcomes for a better understanding of hydrogen oxidation catalysis mechanisms are proposed and discussed.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

A quantum chemical approach for the mechanisms of redox-active metalloenzymes

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes. In our group we have developed an approach that has been applied to a large number of such systems. Hybrid density functional theory (hybrid DFT) has from the start of these investigations been the leading electronic structure tool. An understanding of how the method works in practice has significantly improved the accuracy and applicability. During the past ten years, it has been found that the results for redox enzymes mainly depend on the chosen fraction of exact exchange in the functional, and that a choice of 15% has worked best. The idea has therefore been to vary that fraction over a reasonable range and study the relative energy dependence. For modeling the enzymes, a cluster approach has been developed. In the present review the development of the method we used is described from its start in work on photosystem II, fifteen years ago. Examples from a few recent applications are described, where the metals have been iron, nickel, copper, cobalt or manganese. The results are in excellent agreement with available experiments, and a large number of new predictions have been made.

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes.

Biochemical and artificial pathways for the reduction of carbon dioxide, nitrite and the competing proton reduction: effect of 2nd sphere interactions in catalysis

Reduction of oxides and oxoanions of carbon and nitrogen are of great contemporary importance as they are crucial for a sustainable environment.

Energetics for Proton Reduction in FeFe Hydrogenase

The energetics for proton reduction in FeFe-hydrogenase has been reinvestigated by theoretical modeling, in light of recent experiments. Two different mechanisms have been considered. In the first one, the bridging hydride position was blocked by the enzyme, which is the mechanism that has been supported by a recent spectroscopic study by Cramer et al. A major difficulty in the present study to agree with experimental energetics was to find the right position for the added proton in the first reduction step. It was eventually found that the best position was as a terminal hydride on the distal iron, which has not been suggested in any of the recent, experimentally based mechanisms. The lowest transition state was surprisingly found to be a bond formation between a proton on a cysteine and the terminal hydride. This type of TS is similar to the one for heterolytic H₂ cleavage in NiFe hydrogenase. The second mechanism investigated here is not supported by the present calculations or the recent experiments by Cramer et al., but was still studied as an interesting comparison. In that mechanism, the formation of the bridging hydride was allowed. The H–H formation barrier is only 3.6 kcal/mol higher than for the first mechanism, but there are severe problems concerning the motion of the protons.

Hydroxy-bridged resting states of a [NiFe]-hydrogenase unraveled by cryogenic vibrational spectroscopy and DFT computations

The catalytic mechanism of [NiFe]-hydrogenases is a subject of extensive research. Apart from at least four reaction intermediates of H₂/H⁺ cycling, there are also a number of resting states, which are formed under oxidizing conditions. Although not directly involved in the catalytic cycle, the knowledge of their molecular structures and reactivity is important, because these states usually accumulate in the course of hydrogenase purification and may also play a role in vivo during hydrogenase maturation. Here, we applied low-temperature infrared (cryo-IR) and nuclear resonance vibrational spectroscopy (NRVS) to the isolated catalytic subunit (HoxC) of the heterodimeric regulatory [NiFe]-hydrogenase (RH) from Ralstonia eutropha. Cryo-IR spectroscopy revealed that the HoxC protein can be enriched in almost pure resting redox states suitable for NRVS investigation. NRVS analysis of the hydrogenase catalytic center is usually hampered by strong spectral contributions of the FeS clusters of the small, electron-transferring subunit. Therefore, our approach to investigate the FeS cluster-free, ⁵⁷Fe-labeled HoxC provided an unprecedented insight into the [NiFe] site modes, revealing their contributions in a spectral range otherwise superimposed by FeS cluster-derived bands. Rationalized by density functional theory (DFT) calculations, our data provide structural descriptions of the previously uncharacterized hydroxy- and water-containing resting states. Our work highlights the relevance of cryogenic vibrational spectroscopy and DFT to elucidate the structure of barely defined redox states of the [NiFe]-hydrogenase active site.

Active site vibrations of a [NiFe]-hydrogenase catalytic subunit are selectively probed by IR and NRV spectroscopy in two Ni^IIFe^II and Ni^IIIFe^II resting states, contributing in combination with DFT modeling to rationalized structural candidates.

Rational Design of Fe2 (μ-PR2 )2 (L)6 Coordination Compounds Featuring Tailored Potential Inversion

It was recently discovered that some redox proteins can thermodynamically and spatially split two incoming electrons towards different pathways, resulting in the one-electron reduction of two different substrates, featuring reduction potential respectively higher and lower than the parent reductant. This energy conversion process, referred to as electron bifurcation, is relevant not only from a biochemical perspective, but also for the ground-breaking applications that electron-bifurcating molecular devices could have in the field of energy conversion. Natural electron-bifurcating systems contain a two-electron redox centre featuring potential inversion (PI), i. e. with second reduction easier than the first. With the aim of revealing key factors to tailor the span between first and second redox potentials, we performed a systematic density functional study of a 26-molecule set of models with the general formula Fe₂ (μ-PR₂ )₂ (L)₆ . It turned out that specific features such as i) a Fe-Fe antibonding character of the LUMO, ii) presence of electron-donor groups and iii) low steric congestion in the Fe's coordination sphere, are key ingredients for PI. In particular, the synergic effects of i)-iii) can lead to a span between first and second redox potentials larger than 700 mV. More generally, the "molecular recipes" herein described are expected to inspire the synthesis of Fe₂ P₂ systems with tailored PI, of primary relevance to the design of electron-bifurcating molecular devices.

A Dithiolato and Hydrido Bridged (CO/CN)Fe-Ni Complex with Unprotected CN: A Model for the [Ni-R] State of the [Ni-Fe] Hydrogenase Active Site

A dithiolate/hydride bridged Fe-Ni complex, [(CN)(CO)₂Fe^II(μ-pdt)(μ-H)Ni^II(CN)(PCy₃)]^- (2, pdt = propane-1,3-dithiolate) has been synthesized by the reaction of [(CN)₂(CO)₂Fe^II(pdt)]^2- with [Ni^II(Cl)(H)(PCy₃)₂] as a synthetic analogue of the Ni-R state of the active site of the [Ni-Fe] hydrogenase. X-ray crystallography of this model complex suggests that the hydride unsymmetrically binds to Ni and Fe similar to natural [Ni-Fe] hydrogenases.

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.

Hydrogen evolution reaction mediated by an all-sulfur trinuclear nickel complex

A trinuclear nickel complex with S-based ligands is reported as a bio-inspired model of the [NiFe] hydrogenases' active site. DFT calculations indicate that thiolate and thioether functions are involved as proton relays in the H₂ evolution mechanism.

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.

Conformational Effects of [Ni2 (μ-ArS)2 ] Cores on Their Electrocatalytic Activity

Two nickel complexes supported by tridentate NS₂ ligands, [Ni₂ (κ-N,S,S,S'-N^Ph {CH₂ (MeC₆ H₂ R')S}₂ )₂ ] (1; R'=3,5-(CF₃ )₂ C₆ H₃ ) and [Ni₂ (κ-N,S,S,S'-N^iBu {CH₂ C₆ H₄ S}₂ )₂ ] (2), were prepared as bioinspired models of the active site of [NiFe] hydrogenases. The solid-state structure of 1 reveals that the [Ni₂ (μ-ArS)₂ ] core is bent, with the planes of the nickel centers at a hinge angle of 81.3(5)°, whereas 2 shows a coplanar arrangement between both nickel(II) ions in the dimeric structure. Complex 1 electrocatalyzes proton reduction from CF₃ COOH at -1.93 (overpotential of 1.04 V, with i_cat /i_p ≈21.8) and -1.47 V (overpotential of 580 mV, with i_cat /i_p ≈5.9) versus the ferrocene/ferrocenium redox couple. The electrochemical behavior of 1 relative to that of 2 may be related to the bent [Ni₂ (μ-ArS)₂ ] core, which allows proximity of the two Ni⋅⋅⋅Ni centers at 2.730(8) Å; thus possibly favoring H⁺ reduction. In contrast, the planar [Ni₂ (μ-ArS)₂ ] core of 2 results in a Ni⋅⋅⋅Ni distance of 3.364(4) Å and is unstable in the presence of acid.

O2 sensitivity and H2 production activity of hydrogenases-A review

Hydrogenases are metalloproteins capable of catalyzing the interconversion between molecular hydrogen and protons and electrons. The iron-sulfur clusters within the enzyme enable rapid relay of electrons which are either consumed or generated at the active site. Their unparalleled catalytic efficiency has attracted attention, especially for potential use in H₂ production and/or fuel cell technologies. However, there are limitations to using hydrogenases, especially due to their high O₂ sensitivity. The subclass, called [FeFe] hydrogenases, are particularly more vulnerable to O₂ but proficient in H₂ production. In this review, we provide an overview of mechanistic and protein engineering studies focused on understanding and enhancing O₂ tolerance of the enzyme. The emphasis is on ongoing studies that attempt to overcome O₂ sensitivity of the enzyme while it catalyzes H₂ production in an aerobic environment. We also discuss pioneering attempts to utilize the enzyme in biological H₂ production and other industrial processes, as well as our own perspective on future applications.

Photodynamics of [FeFe]-Hydrogenase Model Compounds with Bidentate Heterocyclic Ligands

Two asymmetrically structured model compounds for the hydrogen-generating [Fe-Fe]-hydrogenase active site were investigated to determine the ultrafast photodynamics, structural intermediates, and photoproducts compared to more common symmetric di-iron species. The bidentate-ligand-containing compounds studied were Fe2(μ-S2C3H6)(CO)4(bipy), 1, and Fe2(μ-S2C3H6)(CO)4(phen), 2, in dilute room temperature acetonitrile solution and low-temperature 2Me-THF matrix isolation using static FTIR difference and time-resolved infrared spectroscopic methods (TRIR). Ultraviolet-visible spectra were also compared to time-dependent density functional theory (TD-DFT) to ascertain the orbital origins of long wavelength electronic absorption features. The spectroscopic evidence supports the conclusions that only a propyl-bridge flip occurs in low-temperature matrix, while early time CO ejection leads to the formation of solvated isomeric species on the 25 ps time scale in room temperature solution.

Can carbene decorated [FeFe]-hydrogenase model complexes catalytically produce dihydrogen? An insight from theory

Cyclic alkyl amino carbene (CAAC) anchored [FeFe]-hydrogenase model complex featuring rotated conformation at one of the iron centers are found to be promising candidate for effective production of dihydrogen. A stepwise comparison of the complete mechanism using the CAAC stabilized model complex [1]⁰ has been performed with that of an experimentally isolated one ([2]⁰). Interestingly, the reduction events involved in the catalytic cycles are found to be more favorable than those previously reported for a similar experimentally known system. Furthermore, the computed ΔpK_a values indicate that the distal iron center with a vacant coordination site is more basic compared to the amino nitrogen atom of the azadithiolate bridge. We also made an attempt to determine the oxidation states of the iron centers for the intermediates involved in the catalytic cycles on the basis of the computed Mössbauer isomer shift and Mulliken spin density values.

Theoretical Studies of Nickel-Dependent Enzymes

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.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Synthetic Designs and Structural Investigations of Biomimetic Ni-Fe Thiolates

Described are the syntheses of several Ni(μ-SR)₂Fe complexes, including hydride derivatives, in a search for improved models for the active site of [NiFe]-hydrogenases. The nickel(II) precursors include (i) nickel with tripodal ligands: Ni(PS₃)^- and Ni(NS₃)^- (PS₃^3- = tris(phenyl-2-thiolato)phosphine, NS₃^3- = tris(benzyl-2-thiolato)amine), (ii) traditional diphosphine-dithiolates, including chiral diphosphine R,R-DIPAMP, (iii) cationic Ni(phosphine-imine/amine) complexes, and (iv) organonickel precursors Ni( o-tolyl)Cl(tmeda) and Ni(C₆F₅)₂. The following new nickel precursor complexes were characterized: PPh₄[Ni(NS₃)] and the dimeric imino/amino-phosphine complexes [NiCl₂(PCH═N^An)]₂ and [NiCl₂(PCH₂NH^An)]₂ (P = Ph₂PC₆H₄-2-). The iron(II) reagents include [CpFe(CO)₂(thf)]BF₄, [Cp*Fe(CO)(MeCN)₂]BF₄, FeI₂(CO)₄, FeCl₂(diphos)(CO)₂, and Fe(pdt)(CO)₂(diphos) (diphos = chelating diphosphines). Reactions of the nickel and iron complexes gave the following new Ni-Fe compounds: Cp*Fe(CO)Ni(NS₃), [Cp(CO)Fe(μ-pdt)Ni(dppbz)]BF₄, [( R,R-DIPAMP)Ni(μ-pdt)(H)Fe(CO)₃]BAr^F₄, [(PCH═N^An)Ni(μ-pdt)(Cl)Fe(dppbz)(CO)]BF₄, [(PCH₂NH^An)Ni(μ-pdt)(Cl)Fe(dppbz)(CO)]BF₄, [(PCH═N^An)Ni(μ-pdt)(H)Fe(dppbz)(CO)]BF₄, [(dppv)(CO)Fe(μ-pdt)]₂Ni, {H[(dppv)(CO)Fe(μ-pdt)]₂Ni]}BF₄, and (C₆F₅)₂Ni(μ-pdt)Fe(CO)₂(dppv) (DIPAMP = (CH₂P(C₆H₄-2-OMe)₂)₂; BAr^F₄^- = [B(C₆H₃-3,5-(CF₃)₂]₄^-)) Within the context of Ni-(SR)₂-Fe complexes, these new complexes feature new microenvironments for the nickel center: tetrahedral Ni, chirality, imine, and amine coligands, and Ni-C bonds. In the case of {H[(dppv)(CO)Fe(μ-pdt)]₂Ni}⁺, four low-energy isomers are separated by ≤3 kcal/mol, one of which features a biomimetic HNi(SR)₄ site, as supported by density functional theory calculations.

A low-valent dinuclear ruthenium diazadiene complex catalyzes the oxidation of dihydrogen and reversible hydrogenation of quinones

The dinuclear ruthenium complex [Ru₂H(μ-H)(Me₂dad)(dbcot)₂] contains a 1,4-dimethyl-diazabuta-1,3-diene (Me₂dad) as a non-innocent bridging ligand between the metal centers to give a [Ru₂(Me₂dad)] core. In addition, each ruthenium is bound to one dibenzo[a,e]cyclooctatetraene (dbcot) ligand. This Ru dimer converts H₂ to protons and electrons. It also catalyzes reversibly under mild conditions the selective hydrogenation of vitamins K₂ and K₃ to their corresponding hydroquinone equivalents without affecting the C[double bond, length as m-dash]C double bonds. Mechanistic studies suggest that the [Ru₂(Me₂dad)] moiety, like hydrogenases, reacts with H₂ and releases electrons and protons stepwise.

Hydrogenases

Hydrogenases catalyze the simple yet important interconversion between H₂ and protons and electrons. Found throughout prokaryotes, lower eukaryotes, and archaea, hydrogenases are used for a variety of redox and signaling purposes and are found in many different forms. This diverse group of metalloenzymes is divided into [NiFe], [FeFe], and [Fe] variants, based on the transition metal contents of their active sites. A wide array of biochemical and spectroscopic methods has been used to elucidate hydrogenases, and this along with a general description of the main enzyme types and catalytic mechanisms is discussed in this chapter.

Pathways of selective catalytic CO2 two-step reduction on di-iron, di-cobalt and iron-cobalt disulfide carbonyls – an electronic structure study

The mixed iron–cobalt disulfide hexacarbonyl provides a selective route in the two step carbon dioxide reduction to formic acid.

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.

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.

Energy-Efficient Hydrogen Evolution by Fe-S Electrocatalysts: Mechanistic Investigations

The intrinsic catalytic property of a Fe-S complex toward H₂ evolution was investigated in a wide range of acids. The title complex exhibited catalytic events at -1.16 and -1.57 V (vs Fc⁺/Fc) in the presence of trifluoromethanesulfonic acid (HOTf) and trifluoroacetic acid (TFA), respectively. The processes corresponded to the single reduction of the Fe-hydride-S-proton and Fe-hydride species, respectively. When anilinium acid was used, the catalysis occurred at -1.16 V, identical with the working potential of the HOTf catalysis, although the employment of anilinium acid was only capable of achieving the Fe-hydride state on the basis of the spectral and calculated results. The thermodynamics and kinetics of individual steps of the catalysis were analyzed by density functional theory (DFT) calculations and electroanalytical simulations. The stepwise CCE or CE (C, chemical; E, electrochemical) mechanism was operative from the HOTf or TFA source, respectively. In contrast, the involvement of anilinium acid most likely initiated a proton-coupled electron transfer (PCET) pathway that avoided the disfavored intermediate after the initial protonation. Via the PCET pathway, the heterogeneous electron transfer rate was increased and the overpotential was decreased by 0.4 V in comparison with the stepwise pathways. The results showed that the PCET-involved catalysis exhibited substantial kinetic and thermodynamic advantages in comparison to the stepwise pathway; thus, an efficient catalytic system for proton reduction was established.

Ultrafast Photodynamics of Cyano-Functionalized [FeFe] Hydrogenase Model Compounds

[FeFe] hydrogenases are efficient enzymes that produce hydrogen gas under mild conditions. Synthetic model compounds containing all CO or mixed CO/PMe3 ligands were previously studied by us and others with ultrafast ultraviolet or visible pump-infrared probe spectroscopy in an effort to better understand the function and interactions of the active site with light. Studies of anionic species containing cyano groups, which more closely match the biological active site, have been elusive. In this work, two model compounds dissolved in room-temperature acetonitrile solution were examined: [Fe2(μ-S2C3H6)(CO)4(CN)2]2- (1) and [Fe2(μ-S2C2H4)(CO)4(CN)2]2- (2). These species exhibit long-lived transient signals consistent with loss of one CO ligand with potential isomerization of newly formed ground electronic state photoproducts, as previously observed with all-CO and CO/PMe3-containing models. We find no evidence for fast (ca. 150 ps) relaxation seen in the all-CO and CO/PMe3 compounds because of the absence of the metal-to-metal charge transfer band in the cyano-functionalized models. These results indicate that incorporation of cyano ligands may significantly alter the electronic properties and photoproducts produced immediately after photoexcitation, which may influence the catalytic activity of model compounds when attached to photosensitizers.