1
|
Dance I. The activating capture of N 2 at the active site of Mo-nitrogenase. Dalton Trans 2024. [PMID: 39140218 DOI: 10.1039/d4dt01866d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2024]
Abstract
Dinitrogen is inherently inert. This report describes detailed density functional calculations (with a 485+ atom model) of mechanistic steps by which the enzyme nitrogenase activates unreactive N2 at the intact active site FeMo-co, to form a key intermediate with bound HNNH. This mechanism does not bind N2 first and then add H atoms, but rather captures N2 ('N2-ready') that diffuses in through the substrate channel and enters a strategic gallery of H atom donors in the reaction zone, between Fe2 and Fe6. This occurs at the E4 stage of the complete mechanism. Exploration of possible reactions of N2 in this space leads to the conclusion that the first reaction step is transfer of H on Fe7 to one end of N2-ready, soon followed by Fe-N bond formation, and then a second H transfer from bridging S2BH to the other N. Two H-N bonds and one or two N-Fe bonds are formed, in some cases with a single transition state. The variable positions and orientations of N2-ready lead to various reaction trajectories and products. The favourable products resulting from this capture, judged by the criteria of reaction energies, reaction barriers, and mechanistic competence for further hydrogenation reactions in the nitrogenase cycle, have Fe2-NH-NH bonding. The trajectory of one N2 capture reaction is described in detail, and calculations that separate the H atom component and the 'heavy atom' components of the classical activation energy are described, in the context of possible H atom tunneling in the activation of N2-ready. I present arguments for the activation of N2 by the pathway of concerted hydrogenation and binding of N2-ready, alternative to the commonly assumed pathway of binding N2 first, with subsequent hydrogenation. The active site of nitrogenase is well primed for the thermodynamic and kinetic advantages of N2 capture.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, NSW 2052, Australia.
| |
Collapse
|
2
|
Jiang H, Ryde U. Putative reaction mechanism of nitrogenase with a half-dissociated S2B ligand. Dalton Trans 2024; 53:11500-11513. [PMID: 38916132 DOI: 10.1039/d4dt00937a] [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: 06/26/2024]
Abstract
We have studied whether dissociation of the S2B sulfide ligand from one of its two coordinating Fe ions may affect the later parts of the reaction mechanism of nitrogenase. Such dissociation has been shown to be favourable for the E2-E4 states in the reaction mechanism, but previous studies have assumed that S2B either remains bridging or has fully dissociated from the active-site FeMo cluster. We employ combined quantum mechanical and molecular mechanical (QM/MM) calculations with two density-functional theory methods, r2SCAN and TPSSh. To make dissociation of S2B possible, we have added a proton to this group throughout the reaction. We study the reaction starting from the E4 state with N2H2 bound to the cluster. Our results indicate that half-dissociation of S2B is unfavourable in most steps of the reaction mechanism. We observe favourable half-dissociation of S2B only when NH or NH2 is bound to the cluster, bridging Fe2 and Fe6. However, the former state is most likely not involved in the reaction mechanism and the latter state is only an intermittent intermediate of the E7 state. Therefore, half-dissociation of S2B seems to play only a minor role in the later parts of the reaction mechanism of nitrogenase. Our suggested mechanism with a protonated S2B is alternating (the two N atoms of the substrate is protonated in an alternating manner) and the substrate prefers to bind to Fe2, in contrast to the preferred binding to Fe6 observed when S2B is unprotonated and bridging Fe2 and Fe6.
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
|
3
|
Hooper RX, Wertz AE, Shafaat HS, Holland PL. Evaluating Diazene to N 2 Interconversion at Iron-Sulfur Complexes. Chemistry 2024; 30:e202304072. [PMID: 38376370 PMCID: PMC11045311 DOI: 10.1002/chem.202304072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 02/19/2024] [Accepted: 02/20/2024] [Indexed: 02/21/2024]
Abstract
Biological N2 reduction occurs at sulfur-rich multiiron sites, and an interesting potential pathway is concerted double reduction/ protonation of bridging N2 through PCET. Here, we test the feasibility of using synthetic sulfur-supported diiron complexes to mimic this pathway. Oxidative proton transfer from μ-η1 : η1-diazene (HN=NH) is the microscopic reverse of the proposed N2 fixation pathway, revealing the energetics of the process. Previously, Sellmann assigned the purple metastable product from two-electron oxidation of [{Fe2+(PPr3)L1}2(μ-η1 : η1-N2H2)] (L1=tetradentate SSSS ligand) at -78 °C as [{Fe2+(PPr3)L1}2(μ-η1 : η1-N2)]2+, which would come from double PCET from diazene to sulfur atoms of the supporting ligands. Using resonance Raman, Mössbauer, NMR, and EPR spectroscopies in conjunction with DFT calculations, we show that the product is not an N2 complex. Instead, the data are most consistent with the spectroscopically observed species being the mononuclear iron(III) diazene complex [{Fe(PPr3)L1}(η2-N2H2)]+. Calculations indicate that the proposed double PCET has a barrier that is too high for proton transfer at the reaction temperature. Also, PCET from the bridging diazene is highly exergonic as a result of the high Fe3+/2+ redox potential, indicating that the reverse N2 protonation would be too endergonic to proceed. This system establishes the "ground rules" for designing reversible N2/N2H2 interconversion through PCET, such as tuning the redox potentials of the metal sites.
Collapse
Affiliation(s)
- Reagan X Hooper
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, CT-06511
| | - Ashlee E Wertz
- Department of Chemistry and Biochemistry, The Ohio State University, 100 W 18th Ave, Columbus, OH-43210
| | - Hannah S Shafaat
- Department of Chemistry and Biochemistry, The Ohio State University, 100 W 18th Ave, Columbus, OH-43210
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, CA-90095
| | - Patrick L Holland
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, CT-06511
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
Tzeliou CE, Mermigki MA, Tzeli D. Review on the QM/MM Methodologies and Their Application to Metalloproteins. Molecules 2022; 27:molecules27092660. [PMID: 35566011 PMCID: PMC9105939 DOI: 10.3390/molecules27092660] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 04/12/2022] [Accepted: 04/18/2022] [Indexed: 12/04/2022] Open
Abstract
The multiscaling quantum mechanics/molecular mechanics (QM/MM) approach was introduced in 1976, while the extensive acceptance of this methodology started in the 1990s. The combination of QM/MM approach with molecular dynamics (MD) simulation, otherwise known as the QM/MM/MD approach, is a powerful and promising tool for the investigation of chemical reactions’ mechanism of complex molecular systems, drug delivery, properties of molecular devices, organic electronics, etc. In the present review, the main methodologies in the multiscaling approaches, i.e., density functional theory (DFT), semiempirical methodologies (SE), MD simulations, MM, and their new advances are discussed in short. Then, a review on calculations and reactions on metalloproteins is presented, where particular attention is given to nitrogenase that catalyzes the conversion of atmospheric nitrogen molecules N₂ into NH₃ through the process known as nitrogen fixation and the FeMo-cofactor.
Collapse
Affiliation(s)
- Christina Eleftheria Tzeliou
- Laboratory of Physical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, 157 71 Athens, Greece; (C.E.T.); (M.A.M.)
| | - Markella Aliki Mermigki
- Laboratory of Physical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, 157 71 Athens, Greece; (C.E.T.); (M.A.M.)
| | - Demeter Tzeli
- Laboratory of Physical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, 157 71 Athens, Greece; (C.E.T.); (M.A.M.)
- Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, 48 Vassileos Constantinou Ave., 116 35 Athens, Greece
- Correspondence: ; Tel.: +30-210-727-4307
| |
Collapse
|
11
|
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
|
12
|
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
|
13
|
Liao RZ, Zhang JX, Lin Z, Siegbahn PE. Antiferromagnetically coupled [Fe8S9] cluster catalyzed acetylene reduction in a nitrogenase-like enzyme DCCPCh: Insights from QM/MM calculations. J Catal 2021. [DOI: 10.1016/j.jcat.2021.04.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
14
|
Palermo G, Spinello A, Saha A, Magistrato A. Frontiers of metal-coordinating drug design. Expert Opin Drug Discov 2020; 16:497-511. [PMID: 33874825 DOI: 10.1080/17460441.2021.1851188] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.
Collapse
Affiliation(s)
- Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, Riverside, United States
| | - Angelo Spinello
- National Research Council (CNR) of Italy, Institute of Material (IOM) @ International School for Advanced Studies (SISSA), Trieste, Italy
| | - Aakash Saha
- Department of Bioengineering, University of California Riverside, Riverside, United States
| | - Alessandra Magistrato
- National Research Council (CNR) of Italy, Institute of Material (IOM) @ International School for Advanced Studies (SISSA), Trieste, Italy
| |
Collapse
|
15
|
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]
|