1
|
Sirohiwal A, Gamiz-Hernandez AP, Kaila VRI. Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase. J Am Chem Soc 2024; 146:18019-18031. [PMID: 38888987 PMCID: PMC11228991 DOI: 10.1021/jacs.4c04476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2024] [Revised: 05/22/2024] [Accepted: 05/23/2024] [Indexed: 06/20/2024]
Abstract
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.
Collapse
Affiliation(s)
- Abhishek Sirohiwal
- Department of Biochemistry
and Biophysics, Stockholm University, Stockholm 10691, Sweden
| | - Ana P. Gamiz-Hernandez
- Department of Biochemistry
and Biophysics, Stockholm University, Stockholm 10691, Sweden
| | - Ville R. I. Kaila
- Department of Biochemistry
and Biophysics, Stockholm University, Stockholm 10691, Sweden
| |
Collapse
|
2
|
Agbaglo DA, Summers TJ, Cheng Q, DeYonker NJ. The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study. Phys Chem Chem Phys 2024; 26:12467-12482. [PMID: 38618904 PMCID: PMC11090134 DOI: 10.1039/d3cp06100k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.
Collapse
Affiliation(s)
- Donatus A Agbaglo
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Thomas J Summers
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| | - Nathan J DeYonker
- Department of Chemistry, University of Memphis, Memphis, TN 38152, USA.
| |
Collapse
|
3
|
Isegawa M. Metal- and ligand-substitution-induced changes in the kinetics and thermodynamics of hydrogen activation and hydricity in a dinuclear metal complex. Dalton Trans 2024; 53:5966-5978. [PMID: 38462977 DOI: 10.1039/d4dt00361f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Catalytic function in organometallic complexes is achieved by carefully selecting their central metals and ligands. In this study, the effects of a metal and a ligand on the kinetics and thermodynamics of hydrogen activation, hydricity degree of the hydride complex, and susceptibility to electronic oxidation in bioinspired NiFe complexes, [NiIIX FeII(Cl)(CO)Y]+ ([NiFe(Cl)(CO)]+; X = N,N'-diethyl-3,7-diazanonane-1,9-dithiolato and Y = 1,2-bis(diphenylphosphino)ethane), were investigated. The density functional theory calculations revealed that the following order thermodynamically favored hydrogen activation: [NiFe(CO)]2+ > [NiRu(CO)]2+ > [NiFe(CNMe)]2+ ∼ [PdRu(CO)]2+ ∼ [PdFe(CO)]2+ ≫ [NiFe(NCS)]+. Moreover, the reverse order thermodynamically favored the hydricity degree.
Collapse
Affiliation(s)
- Miho Isegawa
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka, 819-0395, Japan.
| |
Collapse
|
4
|
Vysotskiy VP, Filippi C, Ryde U. 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. J Phys Chem A 2024; 128:1358-1374. [PMID: 38324717 PMCID: PMC10895656 DOI: 10.1021/acs.jpca.3c05808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 01/03/2024] [Accepted: 01/03/2024] [Indexed: 02/09/2024]
Abstract
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.
Collapse
Affiliation(s)
- Victor P. Vysotskiy
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University
of Twente, P.O. Box 217, Enschede 7500 AE, Netherlands
| | - Ulf Ryde
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| |
Collapse
|
5
|
Hagemann MM, Wieduwilt EK, Hedegård ED. Understanding the initial events of the oxidative damage and protection mechanisms of the AA9 lytic polysaccharide monooxygenase family. Chem Sci 2024; 15:2558-2570. [PMID: 38362420 PMCID: PMC10866358 DOI: 10.1039/d3sc05933b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Accepted: 12/31/2023] [Indexed: 02/17/2024] Open
Abstract
Lytic polysaccharide monooxygenase (LPMO) is a new class of oxidoreductases that boosts polysaccharide degradation employing a copper active site. This boost may facilitate the cost-efficient production of biofuels and high-value chemicals from polysaccharides such as lignocellulose. Unfortunately, self-oxidation of the active site inactivates LPMOs. Other oxidoreductases employ hole-hopping mechanisms as protection against oxidative damage, but little is generally known about the details of these mechanisms. Herein, we employ highly accurate theoretical models based on density functional theory (DFT) molecular mechanics (MM) hybrids to understand the initial steps in LPMOs' protective measures against self-oxidation; we identify several intermediates recently proposed from experiment, and quantify which are important for protective hole-hopping pathways. Investigations on two different LPMOs show consistently that a tyrosine residue close to copper is crucial for protection: this explains recent experiments, showing that LPMOs without this tyrosine are more susceptible to self-oxidation.
Collapse
Affiliation(s)
- Marlisa M Hagemann
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark Campusvej 55 5230 Odense Denmark
| | - Erna K Wieduwilt
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark Campusvej 55 5230 Odense Denmark
| | - Erik D Hedegård
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark Campusvej 55 5230 Odense Denmark
| |
Collapse
|
6
|
He N, Nakatani N, Hada M. How does multi-reference computation change the catalysis chemistry? DFT and CASPT2 studies of the Cu-catalysed coupling reactions between aryl iodides and β-diketones. Phys Chem Chem Phys 2023; 25:28871-28884. [PMID: 37853798 DOI: 10.1039/d3cp03418f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
The molecular mechanism of a Cu-catalysed coupling reaction was theoretically studied using density functional theory (DFT) and the complete active space self-consistent field method followed by the second-order perturbation theory (CASSCF/CASPT2) to investigate the effects of the strong electron correlation of the Cu centre on the reaction profile. Both DFT and CASSCF/CASPT2 calculations showed that the catalytic cycle proceeds via an oxidative addition (OA) reaction, followed by a reductive elimination (RE) reaction, where OA is the rate-determining step. Although the DFT-calculated activation energies of the OA and RE steps are highly dependent on the choice of functionals, the CASSCF/CASPT2 results are less affected by the choice of DFT-optimised geometries. Therefore, with a careful assessment based on the CASSCF/CASPT2 single-point energy evaluation, an optimal choice of the DFT geometry is of good qualitative use for energetics at the CASPT2 level of theory. Based on the changes in the electron populations of the 3d orbitals during the OA and RE steps, the characteristic features of the DFT-calculated electronic structure were qualitatively consistent with those calculated using the CASSCF method. Further electronic structure analysis by the natural orbital occupancy of the CASSCF wavefunction showed that the ground state is almost single-reference in this system and the strong electron correlation effect of the Cu centre can be dealt with using the MP2 or CCSD method, too. However, the slightly smaller occupation numbers of the 3dπ orbital in the course of reactions suggested that the electron correlation effect of the Cu(III) centre appears through the interaction between the 3dπ orbital and the C-I antibonding σ* orbital in the OA step, and between the 3dπ orbital and the Cu-C antibonding σ* orbital in the RE step.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry, Faculty of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo, 192-0397, Japan.
| | - Naoki Nakatani
- Department of Chemistry, Faculty of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo, 192-0397, Japan.
| | - Masahiko Hada
- Department of Chemistry, Faculty of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo, 192-0397, Japan.
| |
Collapse
|
7
|
Vysotskiy VP, Torbjörnsson M, Jiang H, Larsson ED, Cao L, Ryde U, Zhai H, Lee S, Chan GKL. Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art ab initio methods. J Chem Phys 2023; 159:044106. [PMID: 37486046 DOI: 10.1063/5.0152611] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 07/03/2023] [Indexed: 07/25/2023] Open
Abstract
We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.
Collapse
Affiliation(s)
- Victor P Vysotskiy
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Magne Torbjörnsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ernst D Larsson
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Lili Cao
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, SE-221 00 Lund, Sweden
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
8
|
Saitow M, Uemura K, Yanai T. A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals. J Chem Phys 2022; 157:084101. [DOI: 10.1063/5.0094777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field.
Collapse
Affiliation(s)
- Masaaki Saitow
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Kazuma Uemura
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
9
|
Affiliation(s)
- Milica Feldt
- Leibniz Institute for Catalysis: Leibniz-Institut fur Katalyse eV Theory & Catalysis Albert-Einstein-Str 29A 18059 Rostock GERMANY
| | - Quan Manh Phung
- Nagoya University: Nagoya Daigaku Department of Chemistry JAPAN
| |
Collapse
|
10
|
Vitillo JG, Cramer CJ, Gagliardi L. Multireference Methods are Realistic and Useful Tools for Modeling Catalysis. Isr J Chem 2022. [DOI: 10.1002/ijch.202100136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jenny G. Vitillo
- Department of Science and High Technology and INSTM Università degli Studi dell'Insubria Via Valleggio 9 I-22100 Como Italy
| | - Christopher J. Cramer
- Underwriters Laboratories Inc. 333 Pfingsten Road Northbrook Illinois 60602 United States
| | - Laura Gagliardi
- Department of Chemistry Pritzker School of Molecular Engineering James Franck Institute University of Chicago Chicago Illinois 60637 United States
| |
Collapse
|
11
|
Abstract
The level of carbon dioxide in the atmosphere has increased in a dangerous way during the past century. Methods to decrease this level are therefore of high interest at present. Inspiration to do so in an efficient way could come from biological systems. Molybdenum-containing formate dehydrogenase (Mo-FDH) is one of the most interesting enzymes in this respect. For example, the reduction potential required is not very low. The normal reaction catalyzed by Mo-FDH is actually the opposite one of oxidizing formate to CO2. However, recent electrochemical studies have shown that the reaction can be reversed by a moderate lowering of the reduction potential. The goal of the present study has been to study the full mechanism of Mo-FDH, particularly in the most interesting direction of reducing CO2, which has not been done before. The methods used are the same as those that have been shown to give excellent results for redox enzymes in all cases they have been tested. The results obtained for Mo-FDH are also in excellent agreement with the experimental results.
Collapse
Affiliation(s)
- Per E M Siegbahn
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden
| |
Collapse
|
12
|
Orio M, Pantazis DA. Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research. Chem Commun (Camb) 2021; 57:3952-3974. [DOI: 10.1039/d1cc00705j] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Maylis Orio
- Aix-Marseille Université
- CNRS
- iSm2
- Marseille
- France
| | - Dimitrios A. Pantazis
- Max-Planck-Institut für Kohlenforschung
- Kaiser-Wilhelm-Platz 1
- 45470 Mülheim an der Ruhr
- Germany
| |
Collapse
|
13
|
Sugisaki K, Toyota K, Sato K, Shiomi D, Takui T. A quantum algorithm for spin chemistry: a Bayesian exchange coupling parameter calculator with broken-symmetry wave functions. Chem Sci 2020; 12:2121-2132. [PMID: 34163976 PMCID: PMC8179312 DOI: 10.1039/d0sc04847j] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The Heisenberg exchange coupling parameter J (H = −2JSi · Sj) characterises the isotropic magnetic interaction between unpaired electrons, and it is one of the most important spin Hamiltonian parameters of multi-spin open shell systems. The J value is related to the energy difference between high-spin and low-spin states, and thus computing the energies of individual spin states are necessary to obtain the J values from quantum chemical calculations. Here, we propose a quantum algorithm, B̲ayesian ex̲change coupling parameter calculator with b̲roken-symmetry wave functions (BxB), which is capable of computing the J value directly, without calculating the energies of individual spin states. The BxB algorithm is composed of the quantum simulations of the time evolution of a broken-symmetry wave function under the Hamiltonian with an additional term jS2, the wave function overlap estimation with the SWAP test, and Bayesian optimisation of the parameter j. Numerical quantum circuit simulations for H2 under a covalent bond dissociation, C, O, Si, NH, OH+, CH2, NF, O2, and triple bond dissociated N2 molecule revealed that the BxB can compute the J value within 1 kcal mol−1 of errors with less computational costs than conventional quantum phase estimation-based approaches. A quantum algorithm “Bayesian exchange coupling parameter calculator with broken-symmetry wave function (BxB)” enables us to calculate Heisenberg exchange coupling parameter J without inspecting total energies of individual spin states, within 1 kcal mol−1 of energy tolerance.![]()
Collapse
Affiliation(s)
- Kenji Sugisaki
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan .,JST PRESTO 4-1-8 Honcho Kawaguchi Saitama 332-0012 Japan
| | - Kazuo Toyota
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan
| | - Kazunobu Sato
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan
| | - Daisuke Shiomi
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan
| | - Takeji Takui
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan .,Research Support Department, University Research Administrator Centre, University Administration Division, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan
| |
Collapse
|
14
|
Espinosa Ferao A, Rey Planells A, Streubel R. Between Oxirane and Phosphirane: The Spring‐loaded Oxaphosphirane Ring. Eur J Inorg Chem 2020. [DOI: 10.1002/ejic.202000881] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Arturo Espinosa Ferao
- Department of Organic Chemistry Faculty of Chemistry University of Murcia Campus de Espinardo 30100 Murcia
| | - Alicia Rey Planells
- Department of Organic Chemistry Faculty of Chemistry University of Murcia Campus de Espinardo 30100 Murcia
| | - Rainer Streubel
- Institute of Inorganic Chemistry Rheinische Friedrich-Wilhelms-Universität Bonn Gerhard-Domagk-Straße 1
| |
Collapse
|
15
|
Planells AR, Ferao AE. Accurate Ring Strain Energy Calculations on Saturated Three-Membered Heterocycles with One Group 13-16 Element. Inorg Chem 2020; 59:11503-11513. [PMID: 32799495 DOI: 10.1021/acs.inorgchem.0c01316] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate ring strain energy (RSE) data for parent (CH2)2X rings are reported, where X are group 13-16 elements (El) in their lowest oxidation state, from the second to the sixth row, with their covalence completed by bonds to H. They are obtained from appropriate homodesmotic and hyper-homodesmotic reactions at different levels up to the CCSD(T) level, thus providing a benchmark of high-quality reference RSE values, as well as acceptably accurate faster lower-level options. Derivatives of indium, thallium, and lead cannot be properly described by a three-member ring connectivity, because they display a unique donor-acceptor structure from an ethylene π(C═C) orbital to an empty p orbital of a metallylene subunit. The RSE of groups 13 and 14 heterocycles increases on descending in the group (except for Ga and Ge), while it decreases for groups 15 and 16. The latter is presumably due to a strain-releasing mechanism favored by the increase of p-character at the sp3-type atomic orbital used by El in the endocyclic El-C bonds, %p(El)El-C, originated by the tendency of the El lone pairs in groups 15-16 to increase their s-character. This strain-releasing mechanism does not exist in heavier tetrels, which keep almost unchanged the p-character in the endocyclic bonds at El, whereas for triels the p-character is still lower owing to their sp2-like hybridization. Remarkable linear correlations were found between the RSE and either the above-mentioned %p(El)El-C, the distal C-C bond distance or the relaxed force constants for the endocyclic bond angles.
Collapse
Affiliation(s)
- Alicia Rey Planells
- Depto. Quı́mica Orgánica, Facultad de Quı́mica, Campus de Espinardo, Universidad de Murcia, 30100 Murcia, Spain
| | - Arturo Espinosa Ferao
- Depto. Quı́mica Orgánica, Facultad de Quı́mica, Campus de Espinardo, Universidad de Murcia, 30100 Murcia, Spain
| |
Collapse
|
16
|
Ogo S, Kishima T, Yatabe T, Miyazawa K, Yamasaki R, Matsumoto T, Ando T, Kikkawa M, Isegawa M, Yoon KS, Hayami S. [NiFe], [FeFe], and [Fe] hydrogenase models from isomers. SCIENCE ADVANCES 2020; 6:eaaz8181. [PMID: 32577514 PMCID: PMC7286669 DOI: 10.1126/sciadv.aaz8181] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Accepted: 04/08/2020] [Indexed: 06/11/2023]
Abstract
The study of hydrogenase enzymes (H2ases) is necessary because of their importance to a future hydrogen energy economy. These enzymes come in three distinct classes: [NiFe] H2ases, which have a propensity toward H2 oxidation; [FeFe] H2ases, which have a propensity toward H2 evolution; and [Fe] H2ases, 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.
Collapse
Affiliation(s)
- Seiji Ogo
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Takahiro Kishima
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Takeshi Yatabe
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Keishi Miyazawa
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Ryunosuke Yamasaki
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Takahiro Matsumoto
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Tatsuya Ando
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Mitsuhiro Kikkawa
- Department of Chemistry and Biochemistry, Graduate School of Engineering, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Miho Isegawa
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Ki-Seok Yoon
- Center for Small Molecule Energy, Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- International Institute for Carbon-Neutral Energy Research (WPI-ICNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Shinya Hayami
- Graduate School of Science and Technology, Kumamoto University, 2-39-1 Kurokami, Chuo-ku, Kumamoto 860-8555, Japan
| |
Collapse
|
17
|
Affiliation(s)
- Per E. M. Siegbahn
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Rong-Zhen Liao
- Key Laboratory for Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Media, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 430074 Wuhan, People’s Republic of China
| |
Collapse
|
18
|
Farcaș AA, Bende A. Improving the Light-Induced Spin Transition Efficiency in Ni(II)-Based Macrocyclic-Ligand Complexes. Molecules 2019; 24:molecules24234249. [PMID: 31766599 PMCID: PMC6930591 DOI: 10.3390/molecules24234249] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 11/15/2019] [Accepted: 11/20/2019] [Indexed: 01/29/2023] Open
Abstract
The structural stability and photoabsorption properties of Ni(II)-based metal-organic complexes with octahedral coordination having different planar ligand ring structures were investigated employing density functional theory (DFT) and its time-dependent extension (TD-DFT) considering the M06 exchange-correlation functional and the Def2-TZVP basis set. The results showed that the molecular composition of different planar cyclic ligand structures had significant influences on the structural stability and photoabsorption properties of metal-organic complexes. Only those planar ligands that contained aromatic rings met the basic criteria (thermal stability, structural reversibility, and appropriate excitation frequency domain) for light-induced excited spin state trapping, but their spin transition efficiencies were very different. While, in all three aromatic cases, the singlet electronic excitations induced charge distribution that could help in the singlet-to-triplet spin transition, and triplet excitations, which could assist in the backward (triplet-to-singlet) spin transition, was found only for one complex.
Collapse
Affiliation(s)
- Alex-Adrian Farcaș
- National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania;
- Faculty of Physics, “Babeş-Bolyai” University, Mihail Kogalniceanu Street No. 1, Ro-400084 Cluj-Napoca, Romania
| | - Attila Bende
- National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania;
- Correspondence:
| |
Collapse
|
19
|
Matczak P, Mlostoń G, Hamera-Fałdyga R, Görls H, Weigand W. Structure of Diferrocenyl Thioketone: From Molecule to Crystal. MOLECULES (BASEL, SWITZERLAND) 2019; 24:molecules24213950. [PMID: 31683693 PMCID: PMC6864675 DOI: 10.3390/molecules24213950] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 10/28/2019] [Accepted: 10/29/2019] [Indexed: 01/07/2023]
Abstract
Ferrocenyl-functionalized thioketones have recently been recognized as useful building blocks for sulfur-containing compounds with potential applications in materials chemistry. This work is devoted to a single representative of such thioketones, namely diferrocenyl thioketone (Fc2CS), whose structure has been determined here for the first time. Both X-ray crystallography and a wide variety of quantum-chemical methods were used to explore the structure of Fc2CS. In addition to the X-ray structure determination, intermolecular interactions occurring in the crystal structure of Fc2CS were examined in detail by quantum-chemical methods. These methods were also an invaluable tool in studying the molecular structure of Fc2CS, from the gas phase to solutions and to its crystal. Intramolecular interactions governing the conformational behavior of an isolated Fc2CS molecule were deduced from quantum-chemical analyses carried out in orbital space and real space. Our experimental and theoretical results indicate that the main structural features of an isolated Fc2CS molecule in its lowest-energy geometry are retained both upon solvation and in the crystal. The tilt of ferrocenyl groups is only slightly affected by crystal packing forces that are dominated by dispersion. Nonetheless, a network of intermolecular interactions, such as H···H, C···H and S···H, was detected in the Fc2CS crystal but each of them is fairly weak.
Collapse
Affiliation(s)
- Piotr Matczak
- Department of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90236 Lodz, Poland.
| | - Grzegorz Mlostoń
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Róża Hamera-Fałdyga
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Helmar Görls
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| | - Wolfgang Weigand
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| |
Collapse
|
20
|
Gaggioli CA, Stoneburner SJ, Cramer CJ, Gagliardi L. Beyond Density Functional Theory: The Multiconfigurational Approach To Model Heterogeneous Catalysis. ACS Catal 2019. [DOI: 10.1021/acscatal.9b01775] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Carlo Alberto Gaggioli
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Samuel J. Stoneburner
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
21
|
Abstract
The advancements of quantum chemical methods and computer power allow detailed mechanistic investigations of metalloenzymes. In particular, both quantum chemical cluster and combined QM/MM approaches have been used, which have been proven to successfully complement experimental studies. This review starts with a brief introduction of nickel-dependent enzymes and then summarizes theoretical studies on the reaction mechanisms of these enzymes, including NiFe hydrogenase, methyl-coenzyme M reductase, nickel CO dehydrogenase, acetyl CoA synthase, acireductone dioxygenase, quercetin 2,4-dioxygenase, urease, lactate racemase, and superoxide dismutase.
Collapse
|
22
|
Qiu S, Li Q, Xu Y, Shen S, Sun C. Learning from nature: Understanding hydrogenase enzyme using computational approach. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1422] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Siyao Qiu
- Science & Technology Innovation Institute Dongguan University of Technology Dongguan China
| | - Qinye Li
- School of Chemical Engineering Monash University Clayton Victoria Australia
| | - Yongjun Xu
- Science & Technology Innovation Institute Dongguan University of Technology Dongguan China
| | - Shaohua Shen
- International Research Center for Renewable Energy, State Key Laboratory of Multiphase Flow in Power Engineering Xi'an Jiaotong University Shaanxi China
| | - Chenghua Sun
- Department of Chemistry and Biotechnology, and Center for Translational Atomaterials Swinburne University of Technology Hawthorn Victoria Australia
| |
Collapse
|
23
|
New insights into Fe–H$$_{2}$$ and Fe–H$$^{-}$$ bonding of a [NiFe] hydrogenase mimic: a local vibrational mode study. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2463-9] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
24
|
|
25
|
Vogiatzis KD, Polynski MV, Kirkland JK, Townsend J, Hashemi A, Liu C, Pidko EA. Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities. Chem Rev 2019; 119:2453-2523. [PMID: 30376310 PMCID: PMC6396130 DOI: 10.1021/acs.chemrev.8b00361] [Citation(s) in RCA: 214] [Impact Index Per Article: 42.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Indexed: 12/28/2022]
Abstract
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.
Collapse
Affiliation(s)
| | | | - Justin K. Kirkland
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Jacob Townsend
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ali Hashemi
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Chong Liu
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Evgeny A. Pidko
- TheoMAT
group, ITMO University, Lomonosova 9, St. Petersburg 191002, Russia
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| |
Collapse
|
26
|
Cao L, Ryde U. Extremely large differences in DFT energies for nitrogenase models. Phys Chem Chem Phys 2019; 21:2480-2488. [PMID: 30652711 DOI: 10.1039/c8cp06930a] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen avaiable for other organisms. It contains a complicated MoFe7S9C(homocitrate) cluster in its active site. Many computational studies with density-functional theory (DFT) of the nitrogenase enzyme have been presented, but they do not show any consensus - they do not even agree where the first four protons should be added, forming the central intermediate E4. We show that the prime reason for this is that different DFT methods give relative energies that differ by almost 600 kJ mol-1 for different protonation states. This is 4-30 times more than what is observed for other systems. The reason for this is that in some structures, the hydrogens bind to sulfide or carbide ions as protons, whereas in other structures they bind to the metals as hydride ions, changing the oxidation state of the metals, as well as the Fe-C, Fe-S and Fe-Fe distances. The energies correlate with the amount of Hartree-Fock exchange in the method, indicating a variation in the amount of static correlation in the structures. It is currently unclear which DFT method gives the best results for nitrogenase. We show that non-hybrid DFT functionals and TPSSh give the most accurate structures of the resting active site, whereas B3LYP and PBE0 give the best H2 dissociation energies. However, no DFT method indicates that a structure of E4 with two bridging hydride ions is lowest in energy, as spectroscopic experiments indicate.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Lund University, P. O. Box 124, SE-221 00 Lund, Sweden.
| | | |
Collapse
|
27
|
Qiu S, Azofra LM, Macfarlane DR, Sun C. Hydrogen Evolution in [NiFe] Hydrogenases: A Case of Heterolytic Approach between Proton and Hydride. Inorg Chem 2019; 58:2979-2986. [DOI: 10.1021/acs.inorgchem.8b02812] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Siyao Qiu
- Science & Technology Innovation Institute, Dongguan University of Technology, Dongguan 523808, China
- School of Chemistry, Faculty of Science, Monash University, Clayton, VIC 3800, Australia
| | - Luis Miguel Azofra
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Douglas R. Macfarlane
- School of Chemistry, Faculty of Science, Monash University, Clayton, VIC 3800, Australia
- ARC Centre of Excellence for Electromaterials Science (ACES), School of Chemistry, Faculty of Science, Monash University, Clayton, VIC 3800, Australia
| | - Chenghua Sun
- Science & Technology Innovation Institute, Dongguan University of Technology, Dongguan 523808, China
- Department of Chemistry and Biotechnology, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Hawthorn, VIC 3122, Australia
| |
Collapse
|
28
|
Dong G, Phung QM, Pierloot K, Ryde U. Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods. Inorg Chem 2018; 57:15289-15298. [PMID: 30500163 DOI: 10.1021/acs.inorgchem.8b02590] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
[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 H2 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.
Collapse
Affiliation(s)
- Geng Dong
- Department of Theoretical Chemistry, Chemical Centre , Lund University , P.O. Box 124, SE-221 00 Lund , Sweden
- Department of Biochemistry and Molecular Biology , Shantou University Medical College , Shantou 514041 , Guangdong , PR China
| | - Quan Manh Phung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Kristine Pierloot
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre , Lund University , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
29
|
Isegawa M, Sharma AK, Ogo S, Morokuma K. Electron and Hydride Transfer in a Redox-Active NiFe Hydride Complex: A DFT Study. ACS Catal 2018. [DOI: 10.1021/acscatal.8b02368] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Miho Isegawa
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Akhilesh K. Sharma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Seiji Ogo
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University, 744 Moto-oka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
30
|
Escorcia AM, Stein M. QM/MM Investigation of the Role of a Second Coordination Shell Arginine in [NiFe]-Hydrogenases. Front Chem 2018; 6:164. [PMID: 29868565 PMCID: PMC5962704 DOI: 10.3389/fchem.2018.00164] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 04/23/2018] [Indexed: 01/01/2023] Open
Abstract
[NiFe]-hydrogenases are highly efficient catalysts for the heterolytic splitting of molecular hydrogen (H2). 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 H2 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 H2. It is not directly involved in the heterolytic processing of H2 but guides its approach and reduces its flexibility during binding. Upon substitution of the positively charged arginine by a charge-conserving lysine residue, the H2 binding position remains unaffected. However, critical hydrogen bonding interactions with nearby aspartate residues are lost. In addition, the H2 polarization is unfavorable and the reduced side-chain volume may negatively affect the kinetics of the catalytic process.
Collapse
Affiliation(s)
- Andrés M Escorcia
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
| | - Matthias Stein
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
| |
Collapse
|
31
|
Isegawa M, Sharma AK, Ogo S, Morokuma K. DFT Study on Fe(IV)-Peroxo Formation and H Atom Transfer Triggered O2 Activation by NiFe Complex. Organometallics 2018. [DOI: 10.1021/acs.organomet.8b00098] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Miho Isegawa
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University 744 Moto-oka, Nishi-ku, Fukuoka 819-0385, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Akhilesh K. Sharma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Seiji Ogo
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER), Kyushu University 744 Moto-oka, Nishi-ku, Fukuoka 819-0385, Japan
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
32
|
Qiu S, Azofra LM, MacFarlane DR, Sun C. Hydrogen bonding effect between active site and protein environment on catalysis performance in H 2-producing [NiFe] hydrogenases. Phys Chem Chem Phys 2018; 20:6735-6743. [PMID: 29457815 DOI: 10.1039/c7cp07685a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Siyao Qiu
- School of Chemistry, Faculty of Science, Monash University, Clayton, VIC 3800, Australia.
| | | | | | | |
Collapse
|
33
|
Qiu S, Olsen S, MacFarlane DR, Sun C. The oxygen reduction reaction on [NiFe] hydrogenases. Phys Chem Chem Phys 2018; 20:23528-23534. [DOI: 10.1039/c8cp04160a] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Oxygen tolerance capacity is critical for hydrogen oxidation/evolution catalysts.
Collapse
Affiliation(s)
- Siyao Qiu
- Science & Technology Innovation Institute
- Dongguan University of Technology
- Dongguan
- China
- School of Chemistry
| | - Seth Olsen
- School of Chemistry
- Faculty of Science
- Monash University
- Clayton
- VIC 3800
| | | | - Chenghua Sun
- Science & Technology Innovation Institute
- Dongguan University of Technology
- Dongguan
- China
- Department of Chemistry and Biotechnology
| |
Collapse
|
34
|
Dong G, Ryde U, Aa. Jensen HJ, Hedegård ED. Exploration of H2 binding to the [NiFe]-hydrogenase active site with multiconfigurational density functional theory. Phys Chem Chem Phys 2018; 20:794-801. [DOI: 10.1039/c7cp06767d] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
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.
Collapse
Affiliation(s)
- Geng Dong
- Department of Theoretical Chemistry
- Lund University
- Chemical Centre
- SE-221 00 Lund
- Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry
- Lund University
- Chemical Centre
- SE-221 00 Lund
- Sweden
| | - Hans Jørgen Aa. Jensen
- Department of Physics, Chemistry and Pharmacy
- University of Southern Denmark
- DK-5230 Odense M
- Denmark
| | - Erik D. Hedegård
- Department of Theoretical Chemistry
- Lund University
- Chemical Centre
- SE-221 00 Lund
- Sweden
| |
Collapse
|
35
|
Breglia R, Greco C, Fantucci P, De Gioia L, Bruschi M. Theoretical investigation of aerobic and anaerobic oxidative inactivation of the [NiFe]-hydrogenase active site. Phys Chem Chem Phys 2018; 20:1693-1706. [DOI: 10.1039/c7cp06228a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The extraordinary capability of [NiFe]-hydrogenases to catalyse the reversible interconversion of protons and electrons into dihydrogen (H2) has stimulated numerous experimental and theoretical studies addressing the direct utilization of these enzymes in H2 production processes.
Collapse
Affiliation(s)
- Raffaella Breglia
- Department of Earth and Environmental Science
- University of Milano Bicocca
- 20126 Milan
- Italy
| | - Claudio Greco
- Department of Earth and Environmental Science
- University of Milano Bicocca
- 20126 Milan
- Italy
| | - Piercarlo Fantucci
- Department of Biotechnology and Biosciences
- University of Milano Bicocca
- 20126 Milan
- Italy
| | - Luca De Gioia
- Department of Biotechnology and Biosciences
- University of Milano Bicocca
- 20126 Milan
- Italy
| | - Maurizio Bruschi
- Department of Earth and Environmental Science
- University of Milano Bicocca
- 20126 Milan
- Italy
| |
Collapse
|
36
|
Alavi FS, Zahedi M, Safari N, Ryde U. QM/MM Study of the Conversion of Oxophlorin into Verdoheme by Heme Oxygenase. J Phys Chem B 2017; 121:11427-11436. [PMID: 29090581 DOI: 10.1021/acs.jpcb.7b08332] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Heme oxygenase is an enzyme that degrades heme, thereby recycling iron in most organisms, including humans. Pervious density functional theory (DFT) calculations have suggested that iron(III) hydroxyheme, an intermediate generated in the first step of heme degradation by heme oxygenase, is converted to iron(III) superoxo oxophlorin in the presence of dioxygen. In this article, we have studied the detailed mechanism of conversion of iron(III) superoxo oxophlorin to verdoheme by using combined quantum mechanics and molecular mechanics (QM/MM) calculations. The calculations employed the B3LYP method and the def2-QZVP basis set, considering dispersion effects with the DFT-D3 approach, obtaining accurate energies with large QM regions of almost 1000 atoms. The reaction was found to be exothermic by -35 kcal/mol, with a rate-determining barrier of 19 kcal/mol in the doublet state. The protein environment and especially water in the enzyme pocket significantly affects the reaction by decreasing the reaction activation energies and changing the structures by providing strategic hydrogen bonds.
Collapse
Affiliation(s)
- Fatemeh Sadat Alavi
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Mansour Zahedi
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Nasser Safari
- Department of Chemistry, Faculty of Sciences, Shahid Beheshti University , G.C., Evin, 19839-6313 Tehran, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University , Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
37
|
Quantum chemical approaches to [NiFe] hydrogenase. Essays Biochem 2017; 61:293-303. [PMID: 28487405 DOI: 10.1042/ebc20160079] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Revised: 02/22/2017] [Accepted: 03/01/2017] [Indexed: 11/17/2022]
Abstract
The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.
Collapse
|
38
|
Dong G, Phung QM, Hallaert SD, Pierloot K, Ryde U. H2binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods. Phys Chem Chem Phys 2017; 19:10590-10601. [DOI: 10.1039/c7cp01331k] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
CCSD(T) and DMRG-CASPT2 calculations show that H2prefers to bind to Ni rather than to Fe in [NiFe] hydrogenase.
Collapse
Affiliation(s)
- Geng Dong
- Department of Theoretical Chemistry
- Lund University
- SE-221 00 Lund
- Sweden
| | - Quan Manh Phung
- Department of Chemistry
- University of Leuven
- B-3001 Leuven
- Belgium
| | | | | | - Ulf Ryde
- Department of Theoretical Chemistry
- Lund University
- SE-221 00 Lund
- Sweden
| |
Collapse
|
39
|
Fouda A, Ryde U. Does the DFT Self-Interaction Error Affect Energies Calculated in Proteins with Large QM Systems? J Chem Theory Comput 2016; 12:5667-5679. [DOI: 10.1021/acs.jctc.6b00903] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Adam Fouda
- Department of Theoretical
Chemistry, Chemical Centre, Lund University, P. O. Box 124, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical
Chemistry, Chemical Centre, Lund University, P. O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
40
|
dePolo GE, Kaliakin DS, Varganov SA. Spin-Forbidden Transitions between Electronic States in the Active Site of Rubredoxin. J Phys Chem A 2016; 120:8691-8698. [DOI: 10.1021/acs.jpca.6b07717] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gwen E. dePolo
- Department of Chemistry, University of Nevada, Reno, 1664 North
Virginia Street, Reno, Nevada 89557-0216, United States
| | - Danil S. Kaliakin
- Department of Chemistry, University of Nevada, Reno, 1664 North
Virginia Street, Reno, Nevada 89557-0216, United States
| | - Sergey A. Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 North
Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
41
|
Phung QM, Wouters S, Pierloot K. Cumulant Approximated Second-Order Perturbation Theory Based on the Density Matrix Renormalization Group for Transition Metal Complexes: A Benchmark Study. J Chem Theory Comput 2016; 12:4352-61. [DOI: 10.1021/acs.jctc.6b00714] [Citation(s) in RCA: 78] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Quan Manh Phung
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Sebastian Wouters
- Center
for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Kristine Pierloot
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| |
Collapse
|
42
|
Schilter D, Camara JM, Huynh MT, Hammes-Schiffer S, Rauchfuss TB. Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides. Chem Rev 2016; 116:8693-749. [PMID: 27353631 PMCID: PMC5026416 DOI: 10.1021/acs.chemrev.6b00180] [Citation(s) in RCA: 397] [Impact Index Per Article: 49.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- David Schilter
- Center for Multidimensional Carbon Materials, Institute for Basic Science (IBS), Ulsan 44919, Republic of Korea
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Republic of Korea
| | - James M. Camara
- Department of Chemistry, Yeshiva University, 500 West 185th Street, New York, New York 10033, United States
| | - Mioy T. Huynh
- Department of Chemistry, University of Illinois at Urbana–Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana–Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Thomas B. Rauchfuss
- Department of Chemistry, University of Illinois at Urbana–Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
43
|
Qiu S, Azofra LM, MacFarlane DR, Sun C. Unraveling the Role of Ligands in the Hydrogen Evolution Mechanism Catalyzed by [NiFe] Hydrogenases. ACS Catal 2016. [DOI: 10.1021/acscatal.6b01359] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Siyao Qiu
- School
of Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
| | - Luis Miguel Azofra
- School
of Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
- ARC
Centre of Excellence for Electromaterials Science (ACES), School of
Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
| | - Douglas R. MacFarlane
- School
of Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
- ARC
Centre of Excellence for Electromaterials Science (ACES), School of
Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
| | - Chenghua Sun
- School
of Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
- ARC
Centre of Excellence for Electromaterials Science (ACES), School of
Chemistry, Faculty of Science, Monash University, Clayton, Victoria 3800, Australia
| |
Collapse
|
44
|
Mono- and binuclear non-heme iron chemistry from a theoretical perspective. J Biol Inorg Chem 2016; 21:619-44. [DOI: 10.1007/s00775-016-1357-8] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2016] [Accepted: 04/29/2016] [Indexed: 10/21/2022]
|
45
|
Protonation states of intermediates in the reaction mechanism of [NiFe] hydrogenase studied by computational methods. J Biol Inorg Chem 2016; 21:383-94. [DOI: 10.1007/s00775-016-1348-9] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2015] [Accepted: 02/23/2016] [Indexed: 10/22/2022]
|
46
|
|
47
|
Qiu S, Azofra LM, MacFarlane DR, Sun C. Why is a proton transformed into a hydride by [NiFe] hydrogenases? An intrinsic reactivity analysis based on conceptual DFT. Phys Chem Chem Phys 2016; 18:15369-74. [DOI: 10.1039/c6cp00948d] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The hydrogen evolution reaction (HER) catalysed by [NiFe] hydrogenases entails a series of chemical events involving great mechanistic interest.
Collapse
Affiliation(s)
- Siyao Qiu
- School of Chemistry
- Faculty of Science
- Monash University
- Clayton
- Australia
| | | | | | - Chenghua Sun
- School of Chemistry
- Faculty of Science
- Monash University
- Clayton
- Australia
| |
Collapse
|
48
|
Phung QM, Vancoillie S, Pierloot K. Theoretical Study of the Dissociation Energy of First-Row Metallocenium Ions. J Chem Theory Comput 2015; 10:3681-8. [PMID: 26588513 DOI: 10.1021/ct500376c] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The bond dissociation energy of a series of metallocenium ions, i.e., the energy difference of the reaction MCp2(+) → MCp(+) + Cp· (with M = Ti, V, Cr, Mn, Fe, Co, and Ni), was studied by means of multiconfigurational perturbation theory (CASPT2, RASPT2, NEVPT2) and restricted coupled cluster theory (CCSD(T)). From a comparison between the results obtained from these different methods, and a detailed analysis of their treatment of electron correlation effects, a set of MCp(+)-Cp binding energies are proposed with an accuracy of 5 kcal/mol. The computed results are in good agreement with the experimental data measured by threshold photoelectron photoion coincidence (TPEPICO) spectroscopy but disagree with the more recent threshold collision-induced dissociation (TCID) experiments.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Leuven, Flanders, Belgium
| | - Steven Vancoillie
- Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Leuven, Flanders, Belgium
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven , Celestijnenlaan 200F, B-3001 Leuven, Flanders, Belgium
| |
Collapse
|
49
|
Kaliakin DS, Zaari RR, Varganov SA. Effect of H2 Binding on the Nonadiabatic Transition Probability between Singlet and Triplet States of the [NiFe]-Hydrogenase Active Site. J Phys Chem A 2015; 119:1066-73. [DOI: 10.1021/jp510522z] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Danil S. Kaliakin
- Department
of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States
- Department
of Chemistry, Siberian Federal University, 79 Svobodnyi Prospect, Krasnoyarsk, Krasnoyarskiy kray 660041, Russia
| | - Ryan R. Zaari
- Department
of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States
| | - Sergey A. Varganov
- Department
of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
50
|
Gan L, Groy TL, Tarakeshwar P, Mazinani SKS, Shearer J, Mujica V, Jones AK. A Nickel Phosphine Complex as a Fast and Efficient Hydrogen Production Catalyst. J Am Chem Soc 2015; 137:1109-15. [DOI: 10.1021/ja509779q] [Citation(s) in RCA: 122] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Lu Gan
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| | - Thomas L. Groy
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| | - Pilarisetty Tarakeshwar
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| | - Shobeir K. S. Mazinani
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| | - Jason Shearer
- Department
of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Vladimiro Mujica
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| | - Anne K. Jones
- Department
of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287, United States
| |
Collapse
|