Quantum-refinement studies of the bidentate ligand of V‑nitrogenase and the protonation state of CO-inhibited Mo‑nitrogenase

Quantum refinement in real and reciprocal space using the Phenix and ORCA software

X-ray and neutron crystallography, as well as cryogenic electron microscopy (cryo-EM), are the most common methods to obtain atomic structures of biological macromolecules. A feature they all have in common is that, at typical resolutions, the experimental data need to be supplemented by empirical restraints, ensuring that the final structure is chemically reasonable. The restraints are accurate for amino acids and nucleic acids, but often less accurate for substrates, inhibitors, small-molecule ligands and metal sites, for which experimental data are scarce or empirical potentials are harder to formulate. This can be solved using quantum mechanical calculations for a small but interesting part of the structure. Such an approach, called quantum refinement, has been shown to improve structures locally, allow the determination of the protonation and oxidation states of ligands and metals, and discriminate between different interpretations of the structure. Here, we present a new implementation of quantum refinement interfacing the widely used structure-refinement software Phenix and the freely available quantum mechanical software ORCA. Through application to manganese superoxide dismutase and V- and Fe-nitrogenase, we show that the approach works effectively for X-ray and neutron crystal structures, that old results can be reproduced and structural discrimination can be performed. We discuss how the weight factor between the experimental data and the empirical restraints should be selected and how quantum mechanical quality measures such as strain energies should be calculated. We also present an application of quantum refinement to cryo-EM data for particulate methane monooxygenase and show that this may be the method of choice for metal sites in such structures because no accurate empirical restraints are currently available for metals.

Reaction Mechanism for CO Reduction by Mo-Nitrogenase Studied by QM/MM

We have studied the conversion of two molecules of carbon monoxide to ethylene catalyzed by nitrogenase. We start from a recent crystal structure showing the binding of two carbon monoxide molecules to nitrogenase and employ the combined quantum mechanics and molecular mechanics approach. Our results indicate that the reaction is possible only if S2B dissociates as H2S (i.e., the charge of the FeMo cluster remains the same as in the E0 state, indicating that the Fe ions are formally reduced two steps when CO binds). Eight electrons and protons are needed for the reaction, and our mechanism suggests that the first four bind alternatively to the two carbon atoms. The C-C bond formation takes place already after the first protonation (in the E3 state). The next two protons bind to the same O atom, which then dissociates as water. In the same state (E8), the second C-O bond is cleaved, forming the ethylene product. The last two electrons and protons are used to form a water molecule that can be exchanged by S2B or by two CO molecules to start a new reaction cycle.

Analyses of the electronic structures of FeFe-cofactors compared with those of FeMo- and FeV-cofactors and their P-clusters

The electronic structures of FeFe-cofactors (FeFe-cos) in resting and turnover states, together with their PN clusters from iron-only nitrogenases, have been calculated using the bond valence method, and their crystallographic data were reported recently and deposited in the Protein Data Bank (PDB codes: 8BOQ and 8OIE). The calculated results have also been compared with those of their homologous Mo- and V-nitrogenases. For FeFe-cos in the resting state, Fe1/2/4/5/6/7/8 atoms are prone to Fe3+, while the Fe3 atom shows different degrees of mixed valences. The results support that the Fe8 atom at the terminal positions of FeFe-cos possesses the same oxidation states as the Mo3+/V3+ atoms of FeMo-/FeV-cos. In the turnover state, the overall oxidation state of FeFe-co is slightly reduced than those in the resting species, and its electronic configuration is rearranged after the substitution of S2B with OH, compatible with those found in CO-bound FeV-co. Moreover, the calculations give the formal oxidation states of 6Fe2+-2Fe3+ for the electronic structures of PN clusters in Fe-nitrogenases. By the comparison of Mo-, V- and Fe-nitrogenases, the overall oxidation levels of 7Fe atoms (Fe1-Fe7) for both FeFe- and FeMo-cos in resting states are found to be higher than that of FeV-co. For the PN clusters in MoFe-, VFe- and FeFe-proteins, they all exhibit a strong reductive character.

The energetics of N2 reduction by vanadium containing nitrogenase

The main class of nitrogenases has a molybdenum in its cofactor. A mechanism for Mo-nitrogenase has recently been described. In the present study, another class of nitrogenases has been studied, the one with a vanadium instead of a molybdenum in its cofactor. It is generally believed that these classes use the same general mechanism to activate nitrogen. The same methodology has been used here as the one used for Mo-nitrogenase. N2 activation is known to occur after four reductions in the catalytic cycle, in the E4 state. The main features of the mechanism for Mo-nitrogenase found in the previous study are an activation process in four steps prior to catalysis, the release of a sulfide during the activation steps and the formation of H2 from two hydrides in E4, just before N2 is activated. The same features have been found here for V-nitrogenase. A difference is that five steps are needed in the activation process, which explains why the ground state of V-nitrogenase is a triplet (even number) and the one for Mo-nitrogenase is a quartet (odd number). The reason an additional step is needed for V-nitrogenase is that V3+ can be reduced to V2+, in contrast to the case for Mo3+ in Mo-nitrogenase. The fact that V3+ is Jahn-Teller active has important consequences. N2H2 is formed in E4 with reasonably small barriers.

Fe/Zr-MOFs constructed by a sunlight-responsive ligand for efficient photocatalytic nitrogen fixation under ambient condition

Photocatalytic nitrogen fixation opens new opportunities for sustainable and healthier futures, and developing effective and inexpensive photocatalysts is the key. We use the ligand 3,3',5,5'-azomellitic acid (H₄abtc) to connect with Fe clusters and Zr clusters to form stable metal-organic frameworks (MOFs) Fe-abtc and Zr-abtc, both of which are responsive to visible lights for nitrogen fixation. It is worth noting that the presence of NN in the ligand makes it respond to visible lights. The tetracarboxyl group is connected to the metal cluster to form a stable structure. The field-only surface integral method verified that the ligands were successfully applied into the synthesized MOF particles, which expanded the photoresponse range and enhanced the photonic interactions of the synthesized photocatalysts compared with pure MOF particles. The best photocatalytic nitrogen fixation performance of Fe-abtc and Zr-abtc is 49.8 μmol·g_(cat.)^-1·h^-1 and 35.7 μmol·g_(cat.)^-1·h^-1, respectively, the apparent quantum efficiency (AQY) of the sample Fe-abtc is 0.56 %, and the reliability of the source of N element is proved by the isotope ¹⁵N₂. This work provides a new idea for the design of cheap and effective MOFs for photocatalytic nitrogen fixation.

Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

Combining crystallography with quantum mechanics

In standard crystallographic refinement of biomacromolecules, the crystallographic raw data are supplemented by empirical restraints that ensure that the structure makes chemical sense. These restraints are typically accurate for amino acids and nucleic acids, but less so for cofactors, substrates, inhibitors, ligands and metal sites. In quantum refinement, this potential is replaced by more accurate quantum mechanical (QM) calculations. Several implementations have been presented, differing in the level of QM and whether it is used for the entire structure or only for a site of particular interest. It has been shown that the method can improve and correct errors in crystal structures and that it can be used to determine protonation and tautomeric states of various ligands and to decide what is really seen in the structure by refining different interpretations and using standard crystallographic and QM quality measures to decide which fits the structure best.