1
|
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
|
2
|
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
|
3
|
Lukoyanov DA, Yang ZY, Pérez-González A, Raugei S, Dean DR, Seefeldt LC, Hoffman BM. 13C ENDOR Characterization of the Central Carbon within the Nitrogenase Catalytic Cofactor Indicates That the CFe 6 Core Is a Stabilizing "Heart of Steel". J Am Chem Soc 2022; 144:18315-18328. [PMID: 36166637 DOI: 10.1021/jacs.2c06149] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Substrates and inhibitors of Mo-dependent nitrogenase bind and react at Fe ions of the active-site FeMo-cofactor [7Fe-9S-C-Mo-homocitrate] contained within the MoFe protein α-subunit. The cofactor contains a CFe6 core, a carbon centered within a trigonal prism of six Fe, whose role in catalysis is unknown. Targeted 13C labeling of the carbon enables electron-nuclear double resonance (ENDOR) spectroscopy to sensitively monitor the electronic properties of the Fe-C bonds and the spin-coupling scheme adopted by the FeMo-cofactor metal ions. This report compares 13CFe6 ENDOR measurements for (i) the wild-type protein resting state (E0; α-Val70) to those of (ii) α-Ile70, (iii) α-Ala70-substituted proteins; (iv) crystallographically characterized CO-inhibited "hi-CO" state; (v) E4(4H) Janus intermediate, activated for N2 binding/reduction by accumulation of 4[e-/H+]; (vi) E4(2H)* state containing a doubly reduced FeMo-cofactor without Fe-bound substrates; and (vii) propargyl alcohol reduction intermediate having allyl alcohol bound as a ferracycle to FeMo-cofactor Fe6. All states examined, both S = 1/2 and 3/2 exhibited near-zero 13C isotropic hyperfine coupling constants, Ca = [-1.3 ↔ +2.7] MHz. Density functional theory computations and natural bond orbital analysis of the Fe-C bonds show that this occurs because a (3 spin-up/3 spin-down) spin-exchange configuration of CFe6 Fe-ion spins produces cancellation of large spin-transfers to carbon in each Fe-C bond. Previous X-ray diffraction and DFT both indicate that trigonal-prismatic geometry around carbon is maintained with high precision in all these states. The persistent structure and Fe-C bonding of the CFe6 core indicate that it does not provide a functionally dynamic (hemilabile) "beating heart"─instead it acts as "a heart of steel", stabilizing the structure of the FeMo-cofactor-active site during nitrogenase catalysis.
Collapse
Affiliation(s)
- Dmitriy A Lukoyanov
- Departments of Chemistry and Molecular Biosciences, Northwestern University, Evanston, Illinois60208, United States
| | - Zhi-Yong Yang
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah84322, United States
| | - Ana Pérez-González
- Biochemistry Department, Virginia Tech, Blacksburg, Virginia24061, United States
| | - Simone Raugei
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Dennis R Dean
- Biochemistry Department, Virginia Tech, Blacksburg, Virginia24061, United States
| | - Lance C Seefeldt
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah84322, United States
| | - Brian M Hoffman
- Departments of Chemistry and Molecular Biosciences, Northwestern University, Evanston, Illinois60208, United States
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
Benediktsson B, Bjornsson R. Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco. J Chem Theory Comput 2022; 18:1437-1457. [PMID: 35167749 PMCID: PMC8908755 DOI: 10.1021/acs.jctc.1c00753] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
![]()
The open-shell electronic
structure of iron–sulfur clusters
presents considerable challenges to quantum chemistry, with the complex
iron–molybdenum cofactor (FeMoco) of nitrogenase representing
perhaps the ultimate challenge for either wavefunction or density
functional theory. While broken-symmetry density functional theory
has seen some success in describing the electronic structure of such
cofactors, there is a large exchange–correlation functional
dependence in calculations that is not fully understood. In this work,
we present a geometric benchmarking test set, FeMoD11, of synthetic
spin-coupled Fe–Fe and Mo–Fe dimers, with relevance
to the molecular and electronic structure of the Mo-nitrogenase FeMo
cofactor. The reference data consists of high-resolution crystal structures
of metal dimer compounds in different oxidation states. Multiple density
functionals are tested on their ability to reproduce the local geometry,
specifically the Fe–Fe/Mo–Fe distance, for both antiferromagnetically
coupled and ferromagnetically coupled dimers via the broken-symmetry
approach. The metal–metal distance is revealed not only to
be highly sensitive to the amount of exact exchange in the functional
but also to the specific exchange and correlation functionals. For
the antiferromagnetically coupled dimers, the calculated metal–metal
distance correlates well with the covalency of the bridging metal–ligand
bonds, as revealed via the corresponding orbital analysis, Hirshfeld
S/Fe charges, and Fe–S Mayer bond order. Superexchange via
bridging ligands is expected to be the dominant interaction in these
dimers, and our results suggest that functionals that predict accurate
Fe–Fe and Mo–Fe distances describe the overall metal–ligand
covalency more accurately and in turn the superexchange of these systems.
The best performing density functionals of the 16 tested for the FeMoD11
test set are revealed to be either the nonhybrid functionals r2SCAN and B97-D3 or hybrid functionals with 10–15% exact
exchange: TPSSh and B3LYP*. These same four functionals are furthermore
found to reproduce the high-resolution X-ray structure of FeMoco well
according to quantum mechanics/molecular mechanics (QM/MM) calculations.
Almost all nonhybrid functionals systematically underestimate Fe–Fe
and Mo–Fe distances (with r2SCAN and B97-D3 being
the sole exceptions), while hybrid functionals with >15% exact
exchange
(including range-separated hybrid functionals) overestimate them.
The results overall suggest r2SCAN, B97-D3, TPSSh, and
B3LYP* as accurate density functionals for describing the electronic
structure of iron–sulfur clusters in general, including the
complex FeMoco cluster of nitrogenase.
Collapse
Affiliation(s)
- Bardi Benediktsson
- Science Institute, University of Iceland, Dunhagi 3, 107 Reykjavik, Iceland
| | - Ragnar Bjornsson
- Science Institute, University of Iceland, Dunhagi 3, 107 Reykjavik, Iceland.,Max-Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
8
|
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
|