Geometry and Electronic Structure of the P-Cluster in Nitrogenase Studied by Combined Quantum Mechanical and Molecular Mechanical Calculations and Quantum Refinement

Hydrogenase and Nitrogenase: Key Catalysts in Biohydrogen Production

Hydrogen with high energy content is considered to be a promising alternative clean energy source. Biohydrogen production through microbes provides a renewable and immense hydrogen supply by utilizing raw materials such as inexhaustible natural sunlight, water, and even organic waste, which is supposed to solve the two problems of "energy supply and environment protection" at the same time. Hydrogenases and nitrogenases are two classes of key enzymes involved in biohydrogen production and can be applied under different biological conditions. Both the research on enzymatic catalytic mechanisms and the innovations of enzymatic techniques are important and necessary for the application of biohydrogen production. In this review, we introduce the enzymatic structures related to biohydrogen production, summarize recent enzymatic and genetic engineering works to enhance hydrogen production, and describe the chemical efforts of novel synthetic artificial enzymes inspired by the two biocatalysts. Continual studies on the two types of enzymes in the future will further improve the efficiency of biohydrogen production and contribute to the economic feasibility of biohydrogen as an energy source.

Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase

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 P¹⁺/P^N redox potential is nearly the same independently of whether P¹⁺ 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 E₀ state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E₀-E₄ 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 (E₄-E₈) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E₀/E₁ couple.

Role of Serine Coordination in the Structural and Functional Protection of the Nitrogenase P-Cluster

Nitrogenase catalyzes the multielectron reduction of dinitrogen to ammonia. Electron transfer in the catalytic protein (MoFeP) proceeds through a unique [8Fe-7S] cluster (P-cluster) to the active site (FeMoco). In the reduced, all-ferrous (P^N) state, the P-cluster is coordinated by six cysteine residues. Upon two-electron oxidation to the P²⁺ state, the P-cluster undergoes conformational changes in which a highly conserved oxygen-based residue (a Ser or a Tyr) and a backbone amide additionally ligate the cluster. Previous studies of Azotobacter vinelandii (Av) MoFeP revealed that when the oxygen-based residue, βSer188, was mutated to a noncoordinating residue, Ala, the P-cluster became redox-labile and reversibly lost two of its eight Fe centers. Surprisingly, the Av strain with a MoFeP variant that lacked the serine ligand (Av βSer188Ala MoFeP) displayed the same diazotrophic growth and in vitro enzyme turnover rates as wild-type Av MoFeP, calling into question the necessity of this conserved ligand for nitrogenase function. Based on these observations, we hypothesized that βSer188 plays a role in protecting the P-cluster under nonideal conditions. Here, we investigated the protective role of βSer188 both in vivo and in vitro by characterizing the ability of Av βSer188Ala cells to grow under suboptimal conditions (high oxidative stress or Fe limitation) and by determining the tendency of βSer188Ala MoFeP to be mismetallated in vitro. Our results demonstrate that βSer188 (1) increases Av cell survival upon exposure to oxidative stress in the form of hydrogen peroxide, (2) is necessary for efficient Av diazotrophic growth under Fe-limiting conditions, and (3) may protect the P-cluster from metal exchange in vitro. Taken together, our findings suggest a structural adaptation of nitrogenase to protect the P-cluster via Ser ligation, which is a previously unidentified functional role of the Ser residue in redox proteins and adds to the expanding functional roles of non-Cys ligands to FeS clusters.

Statistical analysis of PN clusters in Mo/VFe protein crystals using a bond valence method toward their electronic structures

Iron valences of 129 P-clusters from FeMo/V proteins were analyzed using a bond valence method, supposing the existence of Fe3+ in a generally considered all-ferrous PN cluster in solution with excess reducing agent.

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.

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

Nitrogenase is the only enzyme that can cleave the triple bond in N₂, 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.

Quantum refinement with multiple conformations: application to the P-cluster in nitrogenase

X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms of electron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe8S7Cys6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced PN state (Fe82+, 20%) and the doubly oxidized P2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P1+. Here, it is shown that it is rather a 50/50% mixture of the P1+ and P2+ states and that many of the Fe-Fe and Fe-S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.

Does the crystal structure of vanadium nitrogenase contain a reaction intermediate? Evidence from quantum refinement

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.

Electron Transfer in Nitrogenase

Nitrogenase is the only enzyme capable of reducing N2 to NH3. This challenging reaction requires the coordinated transfer of multiple electrons from the reductase, Fe-protein, to the catalytic component, MoFe-protein, in an ATP-dependent fashion. In the last two decades, there have been significant advances in our understanding of how nitrogenase orchestrates electron transfer (ET) from the Fe-protein to the catalytic site of MoFe-protein and how energy from ATP hydrolysis transduces the ET processes. In this review, we summarize these advances, with focus on the structural and thermodynamic redox properties of nitrogenase component proteins and their complexes, as well as on new insights regarding the mechanism of ET reactions during catalysis and how they are coupled to ATP hydrolysis. We also discuss recently developed chemical, photochemical, and electrochemical methods for uncoupling substrate reduction from ATP hydrolysis, which may provide new avenues for studying the catalytic mechanism of nitrogenase.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.