Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods

Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH)4H]- Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.

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

Gases like H₂, N₂, CO₂, 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 N₂, CO₂, and CO and the production of H₂ 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 N₂ fixation by nitrogenase and H₂ production by hydrogenase as well as CO₂ 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.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

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.

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.

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.