1
|
Dance I. What triggers the coupling of proton transfer and electron transfer at the active site of nitrogenase? Dalton Trans 2024; 53:7996-8004. [PMID: 38651170 DOI: 10.1039/d4dt00474d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
In converting N2 to NH3 the enzyme nitrogenase utilises 8 electrons and 8 protons in the complete catalytic cycle. The source of the electrons is an Fe4S4 reductase protein (Fe-protein) which temporarily docks with the MoFe-protein that contains the catalytic active cofactor, FeMo-co, and an electron transfer cluster called the P cluster. The overall mechanism involves 8 repetitions of a cycle in which reduced Fe-protein docks with the MoFe-protein, one electron transfers to the P-cluster, and then to FeMo-co, followed by dissociation of the two proteins and re-reduction of the Fe-protein. Protons are supplied serially to FeMo-co by a Grotthuss proton translocation mechanism from the protein surface along a conserved chain of water molecules (a proton wire) that terminates near S atoms of the FeMo-co cluster [CFe7S9Mo(homocitrate)] where the multiple steps of the chemical conversions are effected. It is assumed that the chemical mechanisms use proton-coupled electron-transfer (PCET) and that H atoms (e- + H+) are involved in each of the hydrogenation steps. However there is neither evidence for, or mechanism proposed, for this coupling. Here I report calculations of cluster charge distribution upon electron addition, revealing that the added negative charge is on the S atoms of FeMo-co, which thereby become more basic, and able to trigger proton transfer from H3O+ waiting at the near end of the proton wire. This mechanism is supported by calculations of the dynamics of the proton transfer step, in which the barrier is reduced by ca. 3.5 kcal mol-1 and the product stabilised by ca. 7 kcal mol-1 upon electron addition. H tunneling is probable in this step. In nitrogenase it is electron transfer that triggers proton transfer.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, NSW 2052, Australia.
| |
Collapse
|
2
|
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
|
3
|
Jiang H, Ryde U. H 2 formation from the E 2-E 4 states of nitrogenase. Phys Chem Chem Phys 2024; 26:1364-1375. [PMID: 38108422 DOI: 10.1039/d3cp05181a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 × 102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).
Collapse
Affiliation(s)
- Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| |
Collapse
|
4
|
Warmack RA, Rees DC. Nitrogenase beyond the Resting State: A Structural Perspective. Molecules 2023; 28:7952. [PMID: 38138444 PMCID: PMC10745740 DOI: 10.3390/molecules28247952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 11/30/2023] [Accepted: 12/01/2023] [Indexed: 12/24/2023] Open
Abstract
Nitrogenases have the remarkable ability to catalyze the reduction of dinitrogen to ammonia under physiological conditions. How does this happen? The current view of the nitrogenase mechanism focuses on the role of hydrides, the binding of dinitrogen in a reductive elimination process coupled to loss of dihydrogen, and the binding of substrates to a binuclear site on the active site cofactor. This review focuses on recent experimental characterizations of turnover relevant forms of the enzyme determined by cryo-electron microscopy and other approaches, and comparison of these forms to the resting state enzyme and the broader family of iron sulfur clusters. Emerging themes include the following: (i) The obligatory coupling of protein and electron transfers does not occur in synthetic and small-molecule iron-sulfur clusters. The coupling of these processes in nitrogenase suggests that they may involve unique features of the cofactor, such as hydride formation on the trigonal prismatic arrangement of irons, protonation of belt sulfurs, and/or protonation of the interstitial carbon. (ii) Both the active site cofactor and protein are dynamic under turnover conditions; the changes are such that more highly reduced forms may differ in key ways from the resting-state structure. Homocitrate appears to play a key role in coupling cofactor and protein dynamics. (iii) Structural asymmetries are observed in nitrogenase under turnover-relevant conditions by cryo-electron microscopy, although the mechanistic relevance of these states (such as half-of-sites reactivity) remains to be established.
Collapse
Affiliation(s)
- Rebeccah A. Warmack
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
- Howard Hughes Medical Institute, California Institute of Technology, Pasadena, CA 91125, USA
| | - Douglas C. Rees
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
- Howard Hughes Medical Institute, California Institute of Technology, Pasadena, CA 91125, USA
| |
Collapse
|
5
|
Zhai H, Lee S, Cui ZH, Cao L, Ryde U, Chan GKL. Multireference Protonation Energetics of a Dimeric Model of Nitrogenase Iron-Sulfur Clusters. J Phys Chem A 2023; 127:9974-9984. [PMID: 37967028 PMCID: PMC10694817 DOI: 10.1021/acs.jpca.3c06142] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 10/31/2023] [Accepted: 11/01/2023] [Indexed: 11/17/2023]
Abstract
Characterizing the electronic structure of the iron-sulfur clusters in nitrogenase is necessary to understand their role in the nitrogen fixation process. One challenging task is to determine the protonation state of the intermediates in the nitrogen fixing cycle. Here, we use a dimeric iron-sulfur model to study relative energies of protonation at C, S, or Fe. Using a composite method based on coupled cluster and density matrix renormalization group energetics, we converge the relative energies of four protonated configurations with respect to basis set and correlation level. We find that accurate relative energies require large basis sets as well as a proper treatment of multireference and relativistic effects. We have also tested ten density functional approximations for these systems. Most of them give large errors in their relative energies. The best performing functional in this system is B3LYP, which gives mean absolute and maximum deviations of only 10 and 13 kJ/mol with respect to our correlated wave function estimates, respectively, comparable to the uncertainty in our correlated estimates. Our work provides benchmark results for the calibration of new approximate electronic structure methods and density functionals for these problems.
Collapse
Affiliation(s)
- Huanchen Zhai
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Zhi-Hao Cui
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Lili Cao
- Department
of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department
of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
6
|
Barchenko M, O’Malley PJ, de Visser SP. Mechanism of Nitrogen Reduction to Ammonia in a Diiron Model of Nitrogenase. Inorg Chem 2023; 62:14715-14726. [PMID: 37650683 PMCID: PMC10498488 DOI: 10.1021/acs.inorgchem.3c02089] [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: 06/23/2023] [Indexed: 09/01/2023]
Abstract
Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P2P'PhFeH)2(μ-N2)] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.
Collapse
Affiliation(s)
- Maxim Barchenko
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester M1 7DN, U.K.
- Department
of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, U.K.
| | - Patrick J. O’Malley
- Department
of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, U.K.
| | - Sam P. de Visser
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester M1 7DN, U.K.
- Department
of Chemical Engineering, The University
of Manchester, Oxford
Road, Manchester M13 9PL, U.K.
| |
Collapse
|
7
|
Pang Y, Bjornsson R. The E3 state of FeMoco: one hydride, two hydrides or dihydrogen? Phys Chem Chem Phys 2023; 25:21020-21036. [PMID: 37522223 DOI: 10.1039/d3cp01106b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/01/2023]
Abstract
Hydrides are present in the reduced states of the iron-molybdenum cofactor (FeMoco) of Mo nitrogenase and are believed to play a key mechanistic role in the dinitrogen reduction reaction catalyzed by the enzyme. Two hydrides are present in the E4 state according to 1H ENDOR and there is likely a single hydride in the E2 redox state. The 2-hydride E4 state has been experimentally observed to bind N2 and it has been speculated that E3 may bind N2 as well. However, the E3 state has not been directly observed and very little is known about its molecular and electronic structure or reactivity. In recent computational studies, we have explored the energy surfaces of the E2 and E4 by QM/MM modelling, and found that the most stable hydride isomers contain bridging or partially bridging hydrides with an open protonated belt sulfide-bridge. In this work we systematically explore the energy surface of the E3 redox state, comparing single hydride and two-hydride isomers with varying coordination and bridging vs. terminal sulfhydryl groups. We also include a model featuring a triply protonated carbide. The results are only mildly dependent on the QM-region size. The three most stable E3 isomers at the r2SCAN level of theory have in common: an open belt sulfide-bridge (terminal sulfhydryl group on Fe6) and either 2 bridging hydrides (between Fe2 and Fe6), 1 bridging-1-terminal hydride (around Fe2 and Fe6) or a dihydrogen ligand bound at the Fe2 site. Analyzing the functional dependency of the results, we find that functionals previously found to predict accurate structures of spin-coupled Fe/Mo dimers and FeMoco (TPSSh, B97-D3, r2SCAN, and B3LYP*) are in generally good agreement about the stability of these 3 E3 isomers. However, B3LYP*, similar to its parent B3LYP method, predicts a triply protonated carbide isomer as the most stable isomer, an unlikely scenario in view of the lack of experimental evidence for carbide protonation occurring in reduced FeMoco states. Distinguishing further between the 3 hydride isomers is difficult and this flexible coordination nature of hydrides suggests that multiple hydride isomers could be present during experimental conditions. N2 binding was explored and resulted in geometries with 2 bridging hydrides and N2 bound to either Fe2 or Fe6 with a local low-spin state on the Fe. N2 binding is predicted to be mildly endothermic, similar to the E2 state, and it seems unlikely that the E3 state is capable of binding N2.
Collapse
Affiliation(s)
- Yunjie Pang
- College of Chemistry, Beijing Normal University, 100875, Beijing, China
- Max-Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Ragnar Bjornsson
- Max-Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
- Univ. Grenoble Alpes, CNRS, CEA, IRIG, Laboratoire de Chimie et Biologie des Métaux, 17 Rue des Martyrs, F-38054 Grenoble, Cedex, France.
| |
Collapse
|
8
|
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: 0] [Impact Index Per Article: 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
|
9
|
Wang KY, Zhang J, Hsu YC, Lin H, Han Z, Pang J, Yang Z, Liang RR, Shi W, Zhou HC. Bioinspired Framework Catalysts: From Enzyme Immobilization to Biomimetic Catalysis. Chem Rev 2023; 123:5347-5420. [PMID: 37043332 PMCID: PMC10853941 DOI: 10.1021/acs.chemrev.2c00879] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2022] [Indexed: 04/13/2023]
Abstract
Enzymatic catalysis has fueled considerable interest from chemists due to its high efficiency and selectivity. However, the structural complexity and vulnerability hamper the application potentials of enzymes. Driven by the practical demand for chemical conversion, there is a long-sought quest for bioinspired catalysts reproducing and even surpassing the functions of natural enzymes. As nanoporous materials with high surface areas and crystallinity, metal-organic frameworks (MOFs) represent an exquisite case of how natural enzymes and their active sites are integrated into porous solids, affording bioinspired heterogeneous catalysts with superior stability and customizable structures. In this review, we comprehensively summarize the advances of bioinspired MOFs for catalysis, discuss the design principle of various MOF-based catalysts, such as MOF-enzyme composites and MOFs embedded with active sites, and explore the utility of these catalysts in different reactions. The advantages of MOFs as enzyme mimetics are also highlighted, including confinement, templating effects, and functionality, in comparison with homogeneous supramolecular catalysts. A perspective is provided to discuss potential solutions addressing current challenges in MOF catalysis.
Collapse
Affiliation(s)
- Kun-Yu Wang
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- Department
of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry
(MOE) and Renewable Energy Conversion and Storage Center (RECAST),
College of Chemistry, Nankai University, Tianjin 300071, China
| | - Jiaqi Zhang
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- Department
of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry
(MOE) and Renewable Energy Conversion and Storage Center (RECAST),
College of Chemistry, Nankai University, Tianjin 300071, China
| | - Yu-Chuan Hsu
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| | - Hengyu Lin
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| | - Zongsu Han
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- Department
of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry
(MOE) and Renewable Energy Conversion and Storage Center (RECAST),
College of Chemistry, Nankai University, Tianjin 300071, China
| | - Jiandong Pang
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- School
of Materials Science and Engineering, Tianjin Key Laboratory of Metal
and Molecule-Based Material Chemistry, Nankai
University, Tianjin 300350, China
| | - Zhentao Yang
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- Department
of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry
(MOE) and Renewable Energy Conversion and Storage Center (RECAST),
College of Chemistry, Nankai University, Tianjin 300071, China
| | - Rong-Ran Liang
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| | - Wei Shi
- Department
of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry
(MOE) and Renewable Energy Conversion and Storage Center (RECAST),
College of Chemistry, Nankai University, Tianjin 300071, China
| | - Hong-Cai Zhou
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| |
Collapse
|
10
|
Dance I. The binding of reducible N 2 in the reaction domain of nitrogenase. Dalton Trans 2023; 52:2013-2026. [PMID: 36691966 DOI: 10.1039/d2dt03599e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The binding of N2 to FeMo-co, the catalytic site of the enzyme nitrogenase, is central to the conversion to NH3, but also has a separate role in promoting the N2-dependent HD reaction (D2 + 2H+ + 2e- → 2HD). The protein surrounding FeMo-co contains a clear channel for ingress of N2, directly towards the exo-coordination position of Fe2, a position which is outside the catalytic reaction domain. This led to the hypothesis [I. Dance, Dalton Trans., 2022, 51, 12717] of 'promotional' N2 bound at exo-Fe2, and a second 'reducible' N2 bound in the reaction domain, specifically the endo-coordination position of Fe2 or Fe6. The range of possibilities for the binding of reducible N2 in the presence of bound promotional N2 is described here, using density functional simulations with a 486 atom model of the active site and surrounding protein. The pathway for ingress of the second N2 through protein, past the first N2 at exo-Fe2, and tumbling into the binding domain between Fe2 and Fe6, is described. The calculations explore 24 structures involving 6 different forms of hydrogenated FeMo-co, including structures with S2BH unhooked from Fe2 but tethered to Fe6. The calculations use the most probable electronic states. End-on (η1) binding of N2 at the endo position of either Fe2 or Fe6 is almost invariably exothermic, with binding potential energies ranging up to -18 kcal mol-1. Many structures have binding energies in the range -6 to -14 kcal mol-1. The relevant entropic penalty for N2 binding from a diffusible position within the protein is estimated to be 4 kcal mol-1, and so the binding free energies for reducible N2 are suitably negative. N2 binding at endo-Fe2 is stronger than at endo-Fe6 in three of the six structure categories. In many cases the reaction domain containing reducible N2 is expanded. These results inform computational simulation of the subsequent steps in which surrounding H atoms transfer to reducible N2.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, Australia.
| |
Collapse
|
11
|
Near ambient N2 fixation on solid electrodes versus enzymes and homogeneous catalysts. Nat Rev Chem 2023; 7:184-201. [PMID: 37117902 DOI: 10.1038/s41570-023-00462-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/31/2022] [Indexed: 02/04/2023]
Abstract
The Mo/Fe nitrogenase enzyme is unique in its ability to efficiently reduce dinitrogen to ammonia at atmospheric pressures and room temperature. Should an artificial electrolytic device achieve the same feat, it would revolutionize fertilizer production and even provide an energy-dense, truly carbon-free fuel. This Review provides a coherent comparison of recent progress made in dinitrogen fixation on solid electrodes, homogeneous catalysts and nitrogenases. Specific emphasis is placed on systems for which there is unequivocal evidence that dinitrogen reduction has taken place. By establishing the cross-cutting themes and synergies between these systems, we identify viable avenues for future research.
Collapse
|
12
|
Dance I. The HD Reaction of Nitrogenase: a Detailed Mechanism. Chemistry 2023; 29:e202202502. [PMID: 36274057 PMCID: PMC10099629 DOI: 10.1002/chem.202202502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Indexed: 11/06/2022]
Abstract
Nitrogenase is the enzyme that converts N2 to NH3 under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe7 S9 CMo(homocitrate)] is unknown. An obligatory co-product is H2 , while exogenous H2 is a competitive inhibitor. Isotopic substitution using exogenous D2 revealed the N2 -dependent reaction D2 +2H+ +2e- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D2 binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D2 , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N2 -dependence results from promotional N2 bound at exo-Fe2; (c) different N2 binding Km for the HD reaction and the NH3 formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H2 +D2 results from the synchronicity of the two transfers of H to D2 ; (f) inhibition of HD production at high pN2 is by competitive binding of N2 at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW, Sydney, Australia
| |
Collapse
|
13
|
Threatt SD, Rees DC. Biological nitrogen fixation in theory, practice, and reality: a perspective on the molybdenum nitrogenase system. FEBS Lett 2023; 597:45-58. [PMID: 36344435 PMCID: PMC10100503 DOI: 10.1002/1873-3468.14534] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Revised: 10/30/2022] [Accepted: 10/31/2022] [Indexed: 11/09/2022]
Abstract
Nitrogenase is the sole enzyme responsible for the ATP-dependent conversion of atmospheric dinitrogen into the bioavailable form of ammonia (NH3 ), making this protein essential for the maintenance of the nitrogen cycle and thus life itself. Despite the widespread use of the Haber-Bosch process to industrially produce NH3 , biological nitrogen fixation still accounts for half of the bioavailable nitrogen on Earth. An important feature of nitrogenase is that it operates under physiological conditions, where the equilibrium strongly favours ammonia production. This biological, multielectron reduction is a complex catalytic reaction that has perplexed scientists for decades. In this review, we explore the current understanding of the molybdenum nitrogenase system based on experimental and computational research, as well as the limitations of the crystallographic, spectroscopic, and computational techniques employed. Finally, essential outstanding questions regarding the nitrogenase system will be highlighted alongside suggestions for future experimental and computational work to elucidate this essential yet elusive process.
Collapse
Affiliation(s)
- Stephanie D Threatt
- Division of Chemistry and Chemical Engineering, Howard Hughes Medical Institute, California Institute of Technology, Pasadena, CA, USA
| | - Douglas C Rees
- Division of Chemistry and Chemical Engineering, Howard Hughes Medical Institute, California Institute of Technology, Pasadena, CA, USA
| |
Collapse
|
14
|
Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase. MOLECULES (BASEL, SWITZERLAND) 2022; 28:molecules28010065. [PMID: 36615260 PMCID: PMC9822455 DOI: 10.3390/molecules28010065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/18/2022] [Revised: 12/14/2022] [Accepted: 12/19/2022] [Indexed: 12/24/2022]
Abstract
We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron-sulfur clusters with 1-4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400-500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0-E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4-E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.
Collapse
|
15
|
Jiang H, Svensson OKG, Ryde U. QM/MM Study of Partial Dissociation of S2B for the E 2 Intermediate of Nitrogenase. Inorg Chem 2022; 61:18067-18076. [DOI: 10.1021/acs.inorgchem.2c02488] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hao Jiang
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| | - Oskar K. G. Svensson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| |
Collapse
|
16
|
Dance I. Understanding the tethered unhooking and rehooking of S2B in the reaction domain of FeMo-co, the active site of nitrogenase. Dalton Trans 2022; 51:15538-15554. [PMID: 36168836 DOI: 10.1039/d2dt02571j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The active site of the nitrogen fixing enzyme nitrogenase is an Fe7MoS9C cluster, and investigations of the enigmatic chemical mechanism of the enzyme have focussed on a pair of Fe atoms, Fe2 and Fe6, and the S2B atom that bridges them. There are three proposals for the status of the Fe2-S2B-Fe6 bridge during the catalytic cycle: one that it remains intact, another that it is completely labile and absent during catalysis, and a third that S2B is hemilabile, unhooking one of its bonds to Fe2 or Fe6. This report examines the tethered unhooking of S2B and factors that affect it, using DFT calculations of 50 geometric/electronic possibilities with a 485 atom model including all relevant parts of surrounding protein. The outcomes are: (a) unhooking the S2B-Fe2 bond is feasible and favourable, but alternative unhooking of the S2B-Fe6 bond is unlikely for steric reasons, (b) energy differences between hooked and unhooked isomers are generally <10 kcal mol-1, usually with unhooked more stable, (c) ligation at the exo-Fe6 position inhibits unhooking, (d) unhooking of hydrogenated S2B is more favourable than that of bare S2B, (e) hydrogen bonding from the NεH function of His195 to S2B occurs in hooked and unhooked forms, and possibly stabilises unhooking, (f) unhooking is reversible with kinetic barriers ranging 10-13 kcal mol-1. The conclusion is that energetically accessible reversible unhooking of S2B or S2BH, as an intrinsic property of FeMo-co, needs to be considered in the formulation of mechanisms for the reactions of nitrogenase.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, Australia.
| |
Collapse
|
17
|
Jiang H, Svensson OKG, Cao L, Ryde U. Proton Transfer Pathways in Nitrogenase with and without Dissociated S2B. Angew Chem Int Ed Engl 2022; 61:e202208544. [PMID: 35920055 PMCID: PMC9804283 DOI: 10.1002/anie.202208544] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Indexed: 01/05/2023]
Abstract
Nitrogenase is the only enzyme that can convert N2 to NH3 . Crystallographic structures have indicated that one of the sulfide ligands of the active-site FeMo cluster, S2B, can be replaced by an inhibitor, like CO and OH- , and it has been suggested that it may be displaced also during the normal reaction. We have investigated possible proton transfer pathways within the FeMo cluster during the conversion of N2 H2 to two molecules of NH3 , assuming that the protons enter the cluster at the S3B, S4B or S5A sulfide ions and are then transferred to the substrate. We use combined quantum mechanical and molecular mechanical (QM/MM) calculations with the TPSS and B3LYP functionals. The calculations indicate that the barriers for these reactions are reasonable if the S2B ligand remains bound to the cluster, but they become prohibitively high if S2B has dissociated. This suggests that it is unlikely that S2B reversibly dissociates during the normal reaction cycle.
Collapse
Affiliation(s)
- Hao Jiang
- Theoretical ChemistryLund UniversityChemical CentreP. O. Box 12422100LundSweden
| | | | - Lili Cao
- Theoretical ChemistryLund UniversityChemical CentreP. O. Box 12422100LundSweden
| | - Ulf Ryde
- Theoretical ChemistryLund UniversityChemical CentreP. O. Box 12422100LundSweden
| |
Collapse
|
18
|
Dance I. Calculating the chemical mechanism of nitrogenase: new working hypotheses. Dalton Trans 2022; 51:12717-12728. [PMID: 35946501 DOI: 10.1039/d2dt01920e] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The enzyme nitrogenase converts N2 to NH3 with stoichiometry N2 + 8H+ + 8e- → 2NH3 + H2. The mechanism is chemically complex with multiple steps that must be consistent with much accumulated experimental information, including exchange of H2 and N2 and the N2-dependent hydrogenation of D2 to HD. Previous investigations have developed a collection of working hypotheses that guide ongoing density functional investigations of mechanistic steps and sequences. These include (i) hypotheses about the serial provision of protons and their conversion to H atoms bonded to S and Fe atoms of the FeMo-co catalytic site, (ii) the migration of H atoms over the surface of FeMo-co, (iii) the roles of His195, (iv) identification of three protein channels, one for the ingress of N2, a separate pathway for the passage of exogenous H2 (D2) and product H2 (HD), and a hydrophilic pathway for egress of product NH3. Two additional working hypotheses are described in this paper. N2 passing along the N2 channel approaches and binds end-on to the exo coordination position of Fe2, with favourable energetics when FeMo-co is pre-hydrogenated. This exo-Fe2-N2 is apparently not reduced but has a promotional role by expanding the reaction zone. A second N2 can enter via the N2 ingress channel and bind at the endo-Fe6 position, where it is surrounded by H atom donors suitable for the N2 → NH3 conversion. It is proposed that this endo-Fe6 position is also the binding site for H2 (generated or exogenous), accounting for the competitive inhibition of N2 reduction by H2. The HD reaction occurs at the endo-Fe6 site, promoted by N2 at the exo-Fe2 site. The second hypothesis concerns the most stable electronic states of FeMo-co with ligands bound at Fe2 and Fe6, and provides a protocol for management of electronic states in mechanism calculations.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, NSW 2051, Australia.
| |
Collapse
|
19
|
Jiang H, Svensson OKG, Cao L, Ryde U. Proton Transfer Pathways in Nitrogenase with and without Dissociated S2B. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202208544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Hao Jiang
- Lund University: Lunds Universitet Theoretical Chemistry P. O. Box 124 22100 Lund SWEDEN
| | - Oskar K. G. Svensson
- Lund University: Lunds Universitet Theoretical Chemistry P. O. Box 124 22100 Lund SWEDEN
| | - Lili Cao
- Lund University: Lunds Universitet Theoretical Chemistry P. O. Box 124 Lund SWEDEN
| | - Ulf Ryde
- Lund university Theoretical Chemistry P. O. Box 124 S-221 00 Lund SWEDEN
| |
Collapse
|
20
|
Stripp ST, Duffus BR, Fourmond V, Léger C, Leimkühler S, Hirota S, Hu Y, Jasniewski A, Ogata H, Ribbe MW. Second and Outer Coordination Sphere Effects in Nitrogenase, Hydrogenase, Formate Dehydrogenase, and CO Dehydrogenase. Chem Rev 2022; 122:11900-11973. [PMID: 35849738 PMCID: PMC9549741 DOI: 10.1021/acs.chemrev.1c00914] [Citation(s) in RCA: 49] [Impact Index Per Article: 24.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Gases like H2, N2, CO2, and CO are increasingly recognized as critical feedstock in "green" energy conversion and as sources of nitrogen and carbon for the agricultural and chemical sectors. However, the industrial transformation of N2, CO2, and CO and the production of H2 require significant energy input, which renders processes like steam reforming and the Haber-Bosch reaction economically and environmentally unviable. Nature, on the other hand, performs similar tasks efficiently at ambient temperature and pressure, exploiting gas-processing metalloenzymes (GPMs) that bind low-valent metal cofactors based on iron, nickel, molybdenum, tungsten, and sulfur. Such systems are studied to understand the biocatalytic principles of gas conversion including N2 fixation by nitrogenase and H2 production by hydrogenase as well as CO2 and CO conversion by formate dehydrogenase, carbon monoxide dehydrogenase, and nitrogenase. In this review, we emphasize the importance of the cofactor/protein interface, discussing how second and outer coordination sphere effects determine, modulate, and optimize the catalytic activity of GPMs. These may comprise ionic interactions in the second coordination sphere that shape the electron density distribution across the cofactor, hydrogen bonding changes, and allosteric effects. In the outer coordination sphere, proton transfer and electron transfer are discussed, alongside the role of hydrophobic substrate channels and protein structural changes. Combining the information gained from structural biology, enzyme kinetics, and various spectroscopic techniques, we aim toward a comprehensive understanding of catalysis beyond the first coordination sphere.
Collapse
Affiliation(s)
- Sven T Stripp
- Freie Universität Berlin, Experimental Molecular Biophysics, Berlin 14195, Germany
| | | | - Vincent Fourmond
- Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Institut Microbiologie, Bioénergies et Biotechnologie, CNRS, Aix Marseille Université, Marseille 13402, France
| | - Christophe Léger
- Laboratoire de Bioénergétique et Ingénierie des Protéines, Institut de Microbiologie de la Méditerranée, Institut Microbiologie, Bioénergies et Biotechnologie, CNRS, Aix Marseille Université, Marseille 13402, France
| | - Silke Leimkühler
- University of Potsdam, Molecular Enzymology, Potsdam 14476, Germany
| | - Shun Hirota
- Nara Institute of Science and Technology, Division of Materials Science, Graduate School of Science and Technology, Nara 630-0192, Japan
| | - Yilin Hu
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States
| | - Andrew Jasniewski
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States
| | - Hideaki Ogata
- Nara Institute of Science and Technology, Division of Materials Science, Graduate School of Science and Technology, Nara 630-0192, Japan.,Hokkaido University, Institute of Low Temperature Science, Sapporo 060-0819, Japan.,Graduate School of Science, University of Hyogo, Hyogo 678-1297, Japan
| | - Markus W Ribbe
- Department of Molecular Biology & Biochemistry, University of California, Irvine, California 92697-3900, United States.,Department of Chemistry, University of California, Irvine, California 92697-2025, United States
| |
Collapse
|
21
|
Jiang H, Ryde U. Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand. Chemistry 2022; 28:e202103933. [PMID: 35006641 PMCID: PMC9305431 DOI: 10.1002/chem.202103933] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2021] [Indexed: 12/16/2022]
Abstract
We have used combined quantum mechanical and molecular mechanical (QM/MM) calculations to study the reaction mechanism of nitrogenase, assuming that none of the sulfide ligands dissociates. To avoid the problem that there is no consensus regarding the structure and protonation of the E4 state, we start from a state where N2 is bound to the cluster and is protonated to N2H2, after dissociation of H2. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH2) and H2NNH2 as intermediates and the two NH3 products dissociate at the E7 and E8 levels. For all intermediates, coordination to Fe6 is preferred, but for the E4 and E8 intermediates, binding to Fe2 is competitive. For the E4, E5 and E7 intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH2, H2NNH2 and NH3 intermediates. This may explain why homocitrate is a mandatory component of nitrogenase. All steps in the suggested reaction mechanism are thermodynamically favourable compared to protonation of the nearby His‐195 group and in all cases, protonation of the NE2 atom of the latter group is preferred.
Collapse
Affiliation(s)
- Hao Jiang
- Department of Theoretical Chemistry, Lund University Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| |
Collapse
|
22
|
Dance I. Structures and reaction dynamics of N 2 and H 2 binding at FeMo-co, the active site of nitrogenase. Dalton Trans 2021; 50:18212-18237. [PMID: 34860237 DOI: 10.1039/d1dt03548g] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The chemical reactions occurring at the Fe7MoS9C(homocitrate) cluster, FeMo-co, the active site of the enzyme nitrogenase (N2 → NH3), are enigmatic. Experimental information collected over a long period reveals aspects of the roles of N2 and H2, each with more than one type of reactivity. This paper reports investigations of the binding of H2 and N2 at intact FeMo-co, using density functional simulations of a large 486 atom relevant portion of the protein, resulting in 27 new structures containing H2 and/or N2 bound at the exo and endo coordination sites of the participating Fe atoms, Fe2 and Fe6. Binding energies and transition states for association/dissociation are determined, and trajectories for the approach, binding and separation of H2/N2 are described, including diffusion of these small molecules through proximal protein. Influences of surrounding amino acids are identified. FeMo-co deforms geometrically when binding H2 or N2, and a procedure for calculating the energy cost involved, the adaptation energy, is introduced here. Adaptation energies, which range from 7 to 36 kcal mol-1 for the reported structures, are influenced by the protonation state of the His195 side chain. Seven N2 structures and three H2 structures have negative binding free energies, which include the estimated entropy penalties for binding of N2, H2 from proximal protein. These favoured structures have N2 bound end-on at exo-Fe2, exo-Fe6 and endo-Fe2 positions of FeMo-co, and H2 bound at the endo-Fe2 position. Various postulated structures with N2 bridging Fe2 and Fe6 revert to end-on-N2 at endo positions. The structures are also assessed via the calculated potential energy barriers for association and dissociation. Barriers to the binding of H2 range from 1 to 20 kcal mol-1 and barriers to dissociation of H2 range from 3 to 18 kcal mol-1. Barriers to the binding of N2, in either side-on or end-on mode, range from 2 to 18 kcal mol-1, while dissociation of bound N2 encounters barriers of 3 to 8 kcal mol-1 for side-on bonding and 7 to 18 kcal mol-1 for end-on bonding. These results allow formulation of mechanisms for the H2/N2 exchange reaction, and three feasible mechanisms for associative exchange and three for dissociative exchange are identified. Consistent electronic structures and potential energy surfaces are maintained throughout. Changes in the spin populations of Fe2 and Fe6 connected with cluster deformation and with metal-ligand bond formation are identified.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, NSW 2051, Australia.
| |
Collapse
|
23
|
Wang YT, Yang CH, Huang KS, Shaw JF. Chlorophyllides: Preparation, Purification, and Application. Biomolecules 2021; 11:biom11081115. [PMID: 34439782 PMCID: PMC8392590 DOI: 10.3390/biom11081115] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 07/20/2021] [Accepted: 07/20/2021] [Indexed: 12/18/2022] Open
Abstract
Chlorophyllides can be found in photosynthetic organisms. Generally, chlorophyllides have a-, b-, c-, d-, and f-type derivatives, and all chlorophyllides have a tetrapyrrole structure with a Mg ion at the center and a fifth isocyclic pentanone. Chlorophyllide a can be synthesized from protochlorophyllide a, divinyl chlorophyllide a, or chlorophyll. In addition, chlorophyllide a can be transformed into chlorophyllide b, chlorophyllide d, or chlorophyllide f. Chlorophyllide c can be synthesized from protochlorophyllide a or divinyl protochlorophyllide a. Chlorophyllides have been extensively used in food, medicine, and pharmaceutical applications. Furthermore, chlorophyllides exhibit many biological activities, such as anti-growth, antimicrobial, antiviral, antipathogenic, and antiproliferative activity. The photosensitivity of chlorophyllides that is applied in mercury electrodes and sensors were discussed. This article is the first detailed review dedicated specifically to chlorophyllides. Thus, this review aims to describe the definition of chlorophyllides, biosynthetic routes of chlorophyllides, purification of chlorophyllides, and applications of chlorophyllides.
Collapse
Affiliation(s)
- Yi-Ting Wang
- Department of Biological Science and Technology, I-Shou University, Kaohsiung 82445, Taiwan; (Y.-T.W.); (C.-H.Y.)
| | - Chih-Hui Yang
- Department of Biological Science and Technology, I-Shou University, Kaohsiung 82445, Taiwan; (Y.-T.W.); (C.-H.Y.)
- Pharmacy Department of E-Da Hospital, Kaohsiung 82445, Taiwan
- Taiwan Instrument Research Institute, National Applied Research Laboratories, Taipei 106214, Taiwan
| | - Keng-Shiang Huang
- The School of Chinese Medicine for Post-Baccalaureate, I-Shou University, Kaohsiung 82445, Taiwan
- Correspondence: (K.-S.H.); (J.-F.S.); Tel.: +886-7-6151100 (ext. 7063) (K.-S.H.); +886-7-6151100 (ext. 7310) (J.-F.S.); Fax: +886-7-6151959 (J.-F.S.)
| | - Jei-Fu Shaw
- Department of Biological Science and Technology, I-Shou University, Kaohsiung 82445, Taiwan; (Y.-T.W.); (C.-H.Y.)
- Correspondence: (K.-S.H.); (J.-F.S.); Tel.: +886-7-6151100 (ext. 7063) (K.-S.H.); +886-7-6151100 (ext. 7310) (J.-F.S.); Fax: +886-7-6151959 (J.-F.S.)
| |
Collapse
|
24
|
Wu L, Qin L, Nie Y, Xu Y, Zhao YL. Computer-aided understanding and engineering of enzymatic selectivity. Biotechnol Adv 2021; 54:107793. [PMID: 34217814 DOI: 10.1016/j.biotechadv.2021.107793] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Revised: 04/26/2021] [Accepted: 06/28/2021] [Indexed: 12/26/2022]
Abstract
Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.
Collapse
Affiliation(s)
- Lunjie Wu
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China
| | - Lei Qin
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China
| | - Yao Nie
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China; Suqian Industrial Technology Research Institute of Jiangnan University, Suqian 223814, China.
| | - Yan Xu
- School of Biotechnology and Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi 214122, China; State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi 214122, China.
| | - Yi-Lei Zhao
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, MOE-LSB & MOE-LSC, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
25
|
Martín-Fernández C, Harvey JN. On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT. J Phys Chem A 2021; 125:4639-4652. [PMID: 34018759 DOI: 10.1021/acs.jpca.1c01290] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.
Collapse
Affiliation(s)
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan, 200F 3001 Leuven, Belgium
| |
Collapse
|
26
|
Bergmann J, Oksanen E, Ryde U. Quantum-refinement studies of the bidentate ligand of V‑nitrogenase and the protonation state of CO-inhibited Mo‑nitrogenase. J Inorg Biochem 2021; 219:111426. [PMID: 33756394 DOI: 10.1016/j.jinorgbio.2021.111426] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 02/19/2021] [Accepted: 03/10/2021] [Indexed: 10/21/2022]
Abstract
Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to plants (although the enzyme itself is strictly microbial). It has been studied extensively with both experimental and computational methods, but many details of the reaction mechanism are still unclear. X-ray crystallography is the main source of structural information for biomacromolecules, but it has problems to discern hydrogen atoms or to distinguish between elements with the same number of electrons. These problems can sometimes be alleviated by introducing quantum chemical calculations in the refinement, providing information about the ideal structure (in the same way as the empirical restraints used in standard crystallographic refinement) and comparing different interpretations of the structure with normal crystallographic and quantum mechanical quality measures. We have performed such quantum-refinement calculations to address two important issues for nitrogenase. First, we show that the bidentate ligand of the active-site FeV cluster in V‑nitrogenase is carbonate, rather than bicarbonate or nitrate. Second, we study the CO-inhibited structure of Mo‑nitrogenase. CO binds to a reduced and protonated state of the enzyme by replacing one of the sulfide ions (S2B) in the active-site FeMo cluster. We examined if it is possible to deduce from the crystal structure the location of the protons. Our results indicates that the crystal structure is best modelled as fully deprotonated.
Collapse
Affiliation(s)
- Justin Bergmann
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Esko Oksanen
- European Spallation Source ESS ERIC, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden.
| |
Collapse
|
27
|
Masero F, Perrin MA, Dey S, Mougel V. Dinitrogen Fixation: Rationalizing Strategies Utilizing Molecular Complexes. Chemistry 2021; 27:3892-3928. [PMID: 32914919 PMCID: PMC7986120 DOI: 10.1002/chem.202003134] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Indexed: 02/06/2023]
Abstract
Dinitrogen (N2 ) is the most abundant gas in Earth's atmosphere, but its inertness hinders its use as a nitrogen source in the biosphere and in industry. Efficient catalysts are hence required to ov. ercome the high kinetic barriers associated to N2 transformation. In that respect, molecular complexes have demonstrated strong potential to mediate N2 functionalization reactions under mild conditions while providing a straightforward understanding of the reaction mechanisms. This Review emphasizes the strategies for N2 reduction and functionalization using molecular transition metal and actinide complexes according to their proposed reaction mechanisms, distinguishing complexes inducing cleavage of the N≡N bond before (dissociative mechanism) or concomitantly with functionalization (associative mechanism). We present here the main examples of stoichiometric and catalytic N2 functionalization reactions following these strategies.
Collapse
Affiliation(s)
- Fabio Masero
- Department of Chemistry and Applied BiosciencesLaboratory of Inorganic ChemistryETH ZürichVladimir Prelog Weg 1–58093ZürichSwitzerland
| | - Marie A. Perrin
- Department of Chemistry and Applied BiosciencesLaboratory of Inorganic ChemistryETH ZürichVladimir Prelog Weg 1–58093ZürichSwitzerland
| | - Subal Dey
- Department of Chemistry and Applied BiosciencesLaboratory of Inorganic ChemistryETH ZürichVladimir Prelog Weg 1–58093ZürichSwitzerland
| | - Victor Mougel
- Department of Chemistry and Applied BiosciencesLaboratory of Inorganic ChemistryETH ZürichVladimir Prelog Weg 1–58093ZürichSwitzerland
| |
Collapse
|
28
|
Cao L, Ryde U. Putative reaction mechanism of nitrogenase after dissociation of a sulfide ligand. J Catal 2020. [DOI: 10.1016/j.jcat.2020.08.028] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
29
|
Luxem KE, Leavitt WD, Zhang X. Large Hydrogen Isotope Fractionation Distinguishes Nitrogenase-Derived Methane from Other Methane Sources. Appl Environ Microbiol 2020; 86:e00849-20. [PMID: 32709722 PMCID: PMC7499036 DOI: 10.1128/aem.00849-20] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 07/16/2020] [Indexed: 02/01/2023] Open
Abstract
Biological nitrogen fixation is catalyzed by the enzyme nitrogenase. Two forms of this metalloenzyme, the vanadium (V)- and iron (Fe)-only nitrogenases, were recently found to reduce small amounts of carbon dioxide (CO2) into the potent greenhouse gas methane (CH4). Here, we report carbon (13C/12C) and hydrogen (2H/1H) stable isotopic compositions and fractionations of methane generated by V- and Fe-only nitrogenases in the metabolically versatile nitrogen fixer Rhodopseudomonas palustris The stable carbon isotope fractionation imparted by both forms of alternative nitrogenase are within the range observed for hydrogenotrophic methanogenesis (13αCO2/CH4 = 1.051 ± 0.002 for V-nitrogenase and 1.055 ± 0.001 for Fe-only nitrogenase; values are means ± standard errors). In contrast, the hydrogen isotope fractionations (2αH2O/CH4 = 2.071 ± 0.014 for V-nitrogenase and 2.078 ± 0.018 for Fe-only nitrogenase) are the largest of any known biogenic or geogenic pathway. The large 2αH2O/CH4 shows that the reaction pathway nitrogenases use to form methane strongly discriminates against 2H, and that 2αH2O/CH4 distinguishes nitrogenase-derived methane from all other known biotic and abiotic sources. These findings on nitrogenase-derived methane will help constrain carbon and nitrogen flows in microbial communities and the role of the alternative nitrogenases in global biogeochemical cycles.IMPORTANCE All forms of life require nitrogen for growth. Many different kinds of microbes living in diverse environments make inert nitrogen gas from the atmosphere bioavailable using a special enzyme, nitrogenase. Nitrogenase has a wide substrate range, and, in addition to producing bioavailable nitrogen, some forms of nitrogenase also produce small amounts of the greenhouse gas methane. This is different from other microbes that produce methane to generate energy. Until now, there was no good way to determine when microbes with nitrogenases are making methane in nature. Here, we present an isotopic fingerprint that allows scientists to distinguish methane from microbes making it for energy versus those making it as a by-product of nitrogen acquisition. With this new fingerprint, it will be possible to improve our understanding of the relationship between methane production and nitrogen acquisition in nature.
Collapse
Affiliation(s)
- Katja E Luxem
- Department of Geosciences, Princeton University, Princeton, New Jersey, USA
- Princeton Environmental Institute, Princeton University, Princeton, New Jersey, USA
| | - William D Leavitt
- Department of Earth Sciences, Dartmouth College, Hanover, New Hampshire, USA
- Department of Chemistry, Dartmouth College, Hanover, New Hampshire, USA
- Department of Biological Sciences, Dartmouth College, Hanover, New Hampshire, USA
| | - Xinning Zhang
- Department of Geosciences, Princeton University, Princeton, New Jersey, USA
- Princeton Environmental Institute, Princeton University, Princeton, New Jersey, USA
| |
Collapse
|
30
|
Cao L, Caldararu O, Ryde U. Does the crystal structure of vanadium nitrogenase contain a reaction intermediate? Evidence from quantum refinement. J Biol Inorg Chem 2020; 25:847-861. [PMID: 32856107 PMCID: PMC7511287 DOI: 10.1007/s00775-020-01813-z] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Accepted: 08/14/2020] [Indexed: 12/16/2022]
Abstract
Recently, a crystal structure of V-nitrogenase was presented, showing that one of the µ2 sulphide ions in the active site (S2B) is replaced by a lighter atom, suggested to be NH or NH2, i.e. representing a reaction intermediate. Moreover, a sulphur atom is found 7 Å from the S2B site, suggested to represent a storage site for this ion when it is displaced. We have re-evaluated this structure with quantum refinement, i.e. standard crystallographic refinement in which the empirical restraints (employed to ensure that the final structure makes chemical sense) are replaced by more accurate quantum-mechanical calculations. This allows us to test various interpretations of the structure, employing quantum-mechanical calculations to predict the ideal structure and to use crystallographic measures like the real-space Z-score and electron-density difference maps to decide which structure fits the crystallographic raw data best. We show that the structure contains an OH--bound state, rather than an N2-derived reaction intermediate. Moreover, the structure shows dual conformations in the active site with ~ 14% undissociated S2B ligand, but the storage site seems to be fully occupied, weakening the suggestion that it represents a storage site for the dissociated ligand.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Octav Caldararu
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
31
|
Hu Y, Ribbe MW. Special Issue on Nitrogenases and Homologous Systems. Chembiochem 2020; 21:1668-1670. [PMID: 32426925 DOI: 10.1002/cbic.202000279] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Indexed: 12/16/2022]
Abstract
The nitrogenase superfamily comprises homologous enzyme systems that carry out fundamentally important processes, including the reduction of N2 and CO, and the biosynthesis of bacteriochlorophyll and coenzyme F430. This special issue provides a cross-disciplinary overview of the ongoing research in this highly diverse and unique research area of metalloprotein biochemistry.
Collapse
Affiliation(s)
- Yilin Hu
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, 92697-3900, USA
| | - Markus W Ribbe
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, 92697-3900, USA.,Department Chemistry, University of California, Irvine, Irvine, CA, 92697-2025, USA
| |
Collapse
|
32
|
Jin WT, Yang M, Zhu SS, Zhou ZH. Bond-valence analyses of the crystal structures of FeMo/V cofactors in FeMo/V proteins. ACTA CRYSTALLOGRAPHICA SECTION D-STRUCTURAL BIOLOGY 2020; 76:428-437. [PMID: 32355039 DOI: 10.1107/s2059798320003952] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Accepted: 03/18/2020] [Indexed: 12/15/2022]
Abstract
The bond-valence method has been used for valence calculations of FeMo/V cofactors in FeMo/V proteins using 51 crystallographic data sets of FeMo/V proteins from the Protein Data Bank. The calculations show molybdenum(III) to be present in MoFe7S9C(Cys)(HHis)[R-(H)homocit] (where H4homocit is homocitric acid, HCys is cysteine and HHis is histidine) in FeMo cofactors, while vanadium(III) with a more reduced iron complement is obtained for FeV cofactors. Using an error analysis of the calculated valences, it was found that in FeMo cofactors Fe1, Fe6 and Fe7 can be unambiguously assigned as iron(III), while Fe2, Fe3, Fe4 and Fe5 show different degrees of mixed valences for the individual Fe atoms. For the FeV cofactors in PDB entry 5n6y, Fe4, Fe5 and Fe6 correspond to iron(II), iron(II) and iron(III), respectively, while Fe1, Fe2, Fe3 and Fe7 exhibit strongly mixed valences. Special situations such as CO-bound and selenium-substituted FeMo cofactors and O(N)H-bridged FeV cofactors are also discussed and suggest rearrangement of the electron configuration on the substitution of the bridging S atoms.
Collapse
Affiliation(s)
- Wan Ting Jin
- State Key Laboratory for Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Min Yang
- State Key Laboratory for Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shuang Shuang Zhu
- State Key Laboratory for Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Zhao Hui Zhou
- State Key Laboratory for Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|