Dynamics of the [4Fe-4S] cluster in Pyrococcus furiosus D14C ferredoxin via nuclear resonance vibrational and resonance Raman spectroscopies, force field simulations, and density functional theory calculations

Nuclear Resonance Vibrational Spectroscopy: A Modern Tool to Pinpoint Site-Specific Cooperative Processes

Nuclear resonant vibrational spectroscopy (NRVS) is a synchrotron radiation (SR)-based nuclear inelastic scattering spectroscopy that measures the phonons (i.e., vibrational modes) associated with the nuclear transition. It has distinct advantages over traditional vibration spectroscopy and has wide applications in physics, chemistry, bioinorganic chemistry, materials sciences, and geology, as well as many other research areas. In this article, we present a scientific and figurative description of this yet modern tool for the potential users in various research fields in the future. In addition to short discussions on its development history, principles, and other theoretical issues, the focus of this article is on the experimental aspects, such as the instruments, the practical measurement issues, the data process, and a few examples of its applications. The article concludes with introduction to non-57Fe NRVS and an outlook on the impact from the future upgrade of SR rings.

Development of parameters compatible with the CHARMM36 force field for [Fe4S4]2+ clusters and molecular dynamics simulations of adenosine-5'-phosphosulfate reductase in GROMACS 2019

DFT calculations were used to obtain parameters compatible with the CHARMM36 force field for iron-sulfur clusters (Fe-S) of the type [Fe₄S₄]²⁺ that are coordinated to dissimilatory adenosine-5'-phosphosulfate reductase (APSrAB). Classical molecular dynamics (MD) simulations were performed on two APSrAB systems to validate the parameters and verify the stability of the studied systems. The time analysis of the parameters inserted into the force field was in reasonable agreement with the experimental X-ray diffraction data. The analysis of the time evolution of the studied systems indicated that these systems and, in particular, the clusters in their respective cavities had a good stability and were in agreement with what was observed in previous works. The parameters obtained provide the basis for the study of APSrAB as well as other systems that contain [Fe₄S₄]²⁺ through the CHARMM36 force field.

Insights from 125Te and 57Fe nuclear resonance vibrational spectroscopy: a [4Fe-4Te] cluster from two points of view

Can sulfur-to-tellurium exchange serve as a method to understand iron–sulfur clusters of enzymatic systems?

Iron–sulfur clusters are common building blocks for electron transport and active sites of metalloproteins. Their comprehensive investigation is crucial for understanding these enzymes, which play important roles in modern biomimetic catalysis and biotechnology applications. We address this issue by utilizing (Et₄N)₃[Fe₄Te₄(SPh)₄], a tellurium modified version of a conventional reduced [4Fe–4S]⁺ cluster, and performed both ⁵⁷Fe- and ¹²⁵Te-NRVS to reveal its characteristic vibrational features. Our analysis exposed major differences in the resulting ⁵⁷Fe spectrum profile as compared to that of the respective [4Fe–4S] cluster, and between the ⁵⁷Fe and ¹²⁵Te profiles. DFT calculations are applied to rationalize structural, electronic, vibrational, and redox-dependent properties of the [4Fe–4Te]⁺ core. We herein highlight the potential of sulfur/tellurium exchange as a method to isolate the iron-only motion in enzymatic systems.

Terminal Hydride Species in [FeFe]-Hydrogenases Are Vibrationally Coupled to the Active Site Environment

A combination of nuclear resonance vibrational spectroscopy (NRVS), FTIR spectroscopy, and DFT calculations was used to observe and characterize Fe-H/D bending modes in CrHydA1 [FeFe]-hydrogenase Cys-to-Ser variant C169S. Mutagenesis of cysteine to serine at position 169 changes the functional group adjacent to the H-cluster from a -SH to -OH, thus altering the proton transfer pathway. The catalytic activity of C169S is significantly reduced compared to that of native CrHydA1, presumably owing to less efficient proton transfer to the H-cluster. This mutation enabled effective capture of a hydride/deuteride intermediate and facilitated direct detection of the Fe-H/D normal modes. We observed a significant shift to higher frequency in an Fe-H bending mode of the C169S variant, as compared to previous findings with reconstituted native and oxadithiolate (ODT)-substituted CrHydA1. On the basis of DFT calculations, we propose that this shift is caused by the stronger interaction of the -OH group of C169S with the bridgehead -NH- moiety of the active site, as compared to that of the -SH group of C169 in the native enzyme.

Sterically Stabilized Terminal Hydride of a Diiron Dithiolate

The kinetically robust hydride [t-HFe₂(Me₂pdt)(CO)₂(dppv)₂]⁺ ([t-H1]⁺) (Me₂pdt^2- = Me₂C(CH₂S^-)₂; dppv = cis-1,2-C₂H₂(PPh₂)₂) and related derivatives were prepared with ⁵⁷Fe enrichment for characterization by NMR, FT-IR, and NRVS. The experimental results were rationalized using DFT molecular modeling and spectral simulations. The spectroscopic analysis was aimed at supporting assignments of Fe-H vibrational spectra as they relate to recent measurements on [FeFe]-hydrogenase enzymes. The combination of bulky Me₂pdt^2- and dppv ligands stabilizes the terminal hydride with respect to its isomerization to the 5-16 kcal/mol more stable bridging hydride ([μ-H1]⁺) with t_1/2(313.3 K) = 19.3 min. In agreement with the nOe experiments, the calculations predict that one methyl group in [t-H1]⁺ interacts with the hydride with a computed CH···HFe distance of 1.7 Å. Although [t-H⁵⁷1]⁺ exhibits multiple NRVS features in the 720-800 cm^-1 region containing the bending Fe-H modes, the deuterated [t-D⁵⁷1]⁺ sample exhibits a unique Fe-D/CO band at ∼600 cm^-1. In contrast, the NRVS spectra for [μ-H⁵⁷1]⁺ exhibit weaker bands near 670-700 cm^-1 produced by the Fe-H-Fe wagging modes coupled to Me₂pdt^2- and dppv motions.

NRVS for Fe in Biology: Experiment and Basic Interpretation

For over 20 years, nuclear resonance vibrational spectroscopy (NRVS) has been used to study vibrational dynamics of iron-containing materials. With the only selection rule being iron motion, ⁵⁷Fe NRVS has become an excellent tool to study iron-containing enzymes. Over the past decade, considerable progress has been made in the study of complex metalloenzymes using NRVS. Iron cofactors in heme-containing globins; [2Fe2S], [3Fe4S], [4Fe4S] proteins; the [NiFe] and [FeFe] hydrogenases; and nitrogenases have been explored in a fashion not possible through traditional vibrational spectroscopy. In this chapter, we discuss the basics of NRVS, a strategy to perform NRVS, and a discussion of the application of NRVS on rubredoxin and [FeFe] hydrogenase.

Reaction Coordinate Leading to H2 Production in [FeFe]-Hydrogenase Identified by Nuclear Resonance Vibrational Spectroscopy and Density Functional Theory

[FeFe]-hydrogenases are metalloenzymes that reversibly reduce protons to molecular hydrogen at exceptionally high rates. We have characterized the catalytically competent hydride state (Hhyd) in the [FeFe]-hydrogenases from both Chlamydomonas reinhardtii and Desulfovibrio desulfuricans using 57Fe nuclear resonance vibrational spectroscopy (NRVS) and density functional theory (DFT). H/D exchange identified two Fe-H bending modes originating from the binuclear iron cofactor. DFT calculations show that these spectral features result from an iron-bound terminal hydride, and the Fe-H vibrational frequencies being highly dependent on interactions between the amine base of the catalytic cofactor with both hydride and the conserved cysteine terminating the proton transfer chain to the active site. The results indicate that Hhyd is the catalytic state one step prior to H2 formation. The observed vibrational spectrum, therefore, provides mechanistic insight into the reaction coordinate for H2 bond formation by [FeFe]-hydrogenases.

Enzymatic and spectroscopic properties of a thermostable [NiFe]‑hydrogenase performing H2-driven NAD+-reduction in the presence of O2

Biocatalysts that mediate the H₂-dependent reduction of NAD⁺ to NADH are attractive from both a fundamental and applied perspective. Here we present the first biochemical and spectroscopic characterization of an NAD⁺-reducing [NiFe]‑hydrogenase that sustains catalytic activity at high temperatures and in the presence of O₂, which usually acts as an inhibitor. We isolated and sequenced the four structural genes, hoxFUYH, encoding the soluble NAD⁺-reducing [NiFe]‑hydrogenase (SH) from the thermophilic betaproteobacterium, Hydrogenophilus thermoluteolus TH-1^T (Ht). The HtSH was recombinantly overproduced in a hydrogenase-free mutant of the well-studied, H₂-oxidizing betaproteobacterium Ralstonia eutropha H16 (Re). The enzyme was purified and characterized with various biochemical and spectroscopic techniques. Highest H₂-mediated NAD⁺ reduction activity was observed at 80°C and pH6.5, and catalytic activity was found to be sustained at low O₂ concentrations. Infrared spectroscopic analyses revealed a spectral pattern for as-isolated HtSH that is remarkably different from those of the closely related ReSH and other [NiFe]‑hydrogenases. This indicates an unusual configuration of the oxidized catalytic center in HtSH. Complementary electron paramagnetic resonance spectroscopic analyses revealed spectral signatures similar to related NAD⁺-reducing [NiFe]‑hydrogenases. This study lays the groundwork for structural and functional analyses of the HtSH as well as application of this enzyme for H₂-driven cofactor recycling under oxic conditions at elevated temperatures.

Electron-transfer chain in respiratory complex I

Complex I is a part of the respiration energy chain converting the redox energy into the cross-membrane proton gradient. The electron-transfer chain of iron-sulfur cofactors within the water-soluble peripheral part of the complex is responsible for the delivery of electrons to the proton pumping subunit. The protein is porous to water penetration and the hydration level of the cofactors changes when the electron is transferred along the chain. High reaction barriers and trapping of the electrons at the iron-sulfur cofactors are prevented by the combination of intense electrostatic noise produced by the protein-water interface with the high density of quantum states in the iron-sulfur clusters caused by spin interactions between paramagnetic iron atoms. The combination of these factors substantially lowers the activation barrier for electron transfer compared to the prediction of the Marcus theory, bringing the rate to the experimentally established range. The unique role of iron-sulfur clusters as electron-transfer cofactors is in merging protein-water fluctuations with quantum-state multiplicity to allow low activation barriers and robust operation. Water plays a vital role in electron transport energetics by electrowetting the cofactors in the chain upon arrival of the electron. A general property of a protein is to violate the fluctuation-dissipation relation through nonergodic sampling of its landscape. High functional efficiency of redox enzymes is a direct consequence of nonergodicity.

Spectroscopic and Computational Investigations of Ligand Binding to IspH: Discovery of Non-diphosphate Inhibitors

Isoprenoid biosynthesis is an important area for anti-infective drug development. One isoprenoid target is (E)-1-hydroxy-2-methyl-but-2-enyl 4-diphosphate (HMBPP) reductase (IspH), which forms isopentenyl diphosphate and dimethylallyl diphosphate from HMBPP in a 2H+ /2e- reduction. IspH contains a 4 Fe-4 S cluster, and in this work, we first investigated how small molecules bound to the cluster by using HYSCORE and NRVS spectroscopies. The results of these, as well as other structural and spectroscopic investigations, led to the conclusion that, in most cases, ligands bound to IspH 4 Fe-4 S clusters by η1 coordination, forming tetrahedral geometries at the unique fourth Fe, ligand side chains preventing further ligand (e.g., H2 O, O2 ) binding. Based on these ideas, we used in silico methods to find drug-like inhibitors that might occupy the HMBPP substrate binding pocket and bind to Fe, leading to the discovery of a barbituric acid analogue with a Ki value of ≈500 nm against Pseudomonas aeruginosa IspH.

Nitrosylation of Nitric-Oxide-Sensing Regulatory Proteins Containing [4Fe-4S] Clusters Gives Rise to Multiple Iron-Nitrosyl Complexes

The reaction of protein-bound iron-sulfur (Fe-S) clusters with nitric oxide (NO) plays key roles in NO-mediated toxicity and signaling. Elucidation of the mechanism of the reaction of NO with DNA regulatory proteins that contain Fe-S clusters has been hampered by a lack of information about the nature of the iron-nitrosyl products formed. Herein, we report nuclear resonance vibrational spectroscopy (NRVS) and density functional theory (DFT) calculations that identify NO reaction products in WhiD and NsrR, regulatory proteins that use a [4Fe-4S] cluster to sense NO. This work reveals that nitrosylation yields multiple products structurally related to Roussin's Red Ester (RRE, [Fe2 (NO)4 (Cys)2 ]) and Roussin's Black Salt (RBS, [Fe4 (NO)7 S3 ]. In the latter case, the absence of 32 S/34 S shifts in the Fe-S region of the NRVS spectra suggest that a new species, Roussin's Black Ester (RBE), may be formed, in which one or more of the sulfide ligands is replaced by Cys thiolates.

Characterization of the [3Fe-4S](0/1+) cluster from the D14C variant of Pyrococcus furiosus ferredoxin via combined NRVS and DFT analyses

The D14C variant of Pyrococcus furiosus ferredoxin provides an extraordinary framework to investigate a [3Fe-4S] cluster at two oxidation levels and compare the results to its physiologic [4Fe-4S] counterpart in the very same protein. Our spectroscopic and computational study reveals vibrational property changes related to the electronic and structural aspects of both Fe-S clusters.

Hydride bridge in [NiFe]-hydrogenase observed by nuclear resonance vibrational spectroscopy

The metabolism of many anaerobes relies on [NiFe]-hydrogenases, whose characterization when bound to substrates has proven non-trivial. Presented here is direct evidence for a hydride bridge in the active site of the ⁵⁷Fe-labelled fully reduced Ni-R form of Desulfovibrio vulgaris Miyazaki F [NiFe]-hydrogenase. A unique ‘wagging' mode involving H⁻ motion perpendicular to the Ni(μ-H)⁵⁷Fe plane was studied using ⁵⁷Fe-specific nuclear resonance vibrational spectroscopy and density functional theory (DFT) calculations. On Ni(μ-D)⁵⁷Fe deuteride substitution, this wagging causes a characteristic perturbation of Fe–CO/CN bands. Spectra have been interpreted by comparison with Ni(μ-H/D)⁵⁷Fe enzyme mimics [(dppe)Ni(μ-pdt)(μ-H/D)⁵⁷Fe(CO)₃]⁺ and DFT calculations, which collectively indicate a low-spin Ni(II)(μ-H)Fe(II) core for Ni-R, with H⁻ binding Ni more tightly than Fe. The present methodology is also relevant to characterizing Fe–H moieties in other important natural and synthetic catalysts.

Understanding the catalytic mechanism of redox-active hydrogenases is a key to efficient hydrogen production and consumption. Here, the authors use nuclear resonance vibrational spectroscopy to study [NiFe]-hydrogenase, and observe a bridging hydride structure in an EPR silent intermediate.

Low frequency dynamics of the nitrogenase MoFe protein via femtosecond pump probe spectroscopy - Observation of a candidate promoting vibration

We have used femtosecond pump-probe spectroscopy (FPPS) to study the FeMo-cofactor within the nitrogenase (N2ase) MoFe protein from Azotobacter vinelandii. A sub-20-fs visible laser pulse was used to pump the sample to an excited electronic state, and a second sub-10-fs pulse was used to probe changes in transmission as a function of probe wavelength and delay time. The excited protein relaxes to the ground state with a ~1.2ps time constant. With the short laser pulse we coherently excited the vibrational modes associated with the FeMo-cofactor active site, which are then observed in the time domain. Superimposed on the relaxation dynamics, we distinguished a variety of oscillation frequencies with the strongest band peaks at ~84, 116, 189, and 226cm(-1). Comparison with data from nuclear resonance vibrational spectroscopy (NRVS) shows that the latter pair of signals comes predominantly from the FeMo-cofactor. The frequencies obtained from the FPPS experiment were interpreted with normal mode calculations using both an empirical force field (EFF) and density functional theory (DFT). The FPPS data were also compared with the first reported resonance Raman (RR) spectrum of the N2ase MoFe protein. This approach allows us to outline and assign vibrational modes having relevance to the catalytic activity of N2ase. In particular, the 226cm(-1) band is assigned as a potential 'promoting vibration' in the H-atom transfer (or proton-coupled electron transfer) processes that are an essential feature of N2ase catalysis. The results demonstrate that high-quality room-temperature solution data can be obtained on the MoFe protein by the FPPS technique and that these data provide added insight to the motions and possible operation of this protein and its catalytic prosthetic group.

Isoprenoid Biosynthesis in Pathogenic Bacteria: Nuclear Resonance Vibrational Spectroscopy Provides Insight into the Unusual [4Fe-4S] Cluster of the E. coli LytB/IspH Protein

The LytB/IspH protein catalyzes the last step of the methylerythritol phosphate (MEP) pathway which is used for the biosynthesis of essential terpenoids in most pathogenic bacteria. Therefore, the MEP pathway is a target for the development of new antimicrobial agents as it is essential for microorganisms, yet absent in humans. Substrate-free LytB has a special [4Fe-4S](2+) cluster with a yet unsolved structure. This motivated us to use synchrotron-based nuclear resonance vibrational spectroscopy (NRVS) in combination with quantum chemical-molecular mechanical (QM/MM) calculations to gain more insight into the structure of substrate-free LytB. The apical iron atom of the [4Fe-4S](2+) is clearly linked to three water molecules. We additionally present NRVS data of LytB bound to its natural substrate, (E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate (HMBPP) and to the inhibitors (E)-4-amino-3-methylbut-2-en-1-yl diphosphate and (E)-4-mercapto-3-methylbut-2-en-1-yl diphosphate.

Concepts in bio-molecular spectroscopy: vibrational case studies on metalloenzymes

Challenges and chances in bio-molecular spectroscopy are exemplified by vibrational case studies on metalloenzymes.

Nuclear resonance vibrational spectroscopy reveals the FeS cluster composition and active site vibrational properties of an O2-tolerant NAD+-reducing [NiFe] hydrogenase

Nuclear resonance vibrational spectroscopy is used to characterize all Fe-containing cofactors in a complex multicofactor enzyme.

Hydrogenases are complex metalloenzymes that catalyze the reversible splitting of molecular hydrogen into protons and electrons essentially without overpotential. The NAD⁺-reducing soluble hydrogenase (SH) from Ralstonia eutropha is capable of H₂ conversion even in the presence of usually toxic dioxygen. The molecular details of the underlying reactions are largely unknown, mainly because of limited knowledge of the structure and function of the various metal cofactors present in the enzyme. Here, all iron-containing cofactors of the SH were investigated by ⁵⁷Fe specific nuclear resonance vibrational spectroscopy (NRVS). Our data provide experimental evidence for one [2Fe2S] center and four [4Fe4S] clusters, which is consistent with the amino acid sequence composition. Only the [2Fe2S] cluster and one of the four [4Fe4S] clusters were reduced upon incubation of the SH with NADH. This finding explains the discrepancy between the large number of FeS clusters and the small amount of FeS cluster-related signals as detected by electron paramagnetic resonance spectroscopic analysis of several NAD⁺-reducing hydrogenases. For the first time, Fe–CO and Fe–CN modes derived from the [NiFe] active site could be distinguished by NRVS through selective ¹³C labeling of the CO ligand. This strategy also revealed the molecular coordinates that dominate the individual Fe–CO modes. The present approach explores the complex vibrational signature of the Fe–S clusters and the hydrogenase active site, thereby showing that NRVS represents a powerful tool for the elucidation of complex biocatalysts containing multiple cofactors.

Energy calibration issues in nuclear resonant vibrational spectroscopy: observing small spectral shifts and making fast calibrations

The conventional energy calibration for nuclear resonant vibrational spectroscopy (NRVS) is usually long. Meanwhile, taking NRVS samples out of the cryostat increases the chance of sample damage, which makes it impossible to carry out an energy calibration during one NRVS measurement. In this study, by manipulating the 14.4 keV beam through the main measurement chamber without moving out the NRVS sample, two alternative calibration procedures have been proposed and established: (i) an in situ calibration procedure, which measures the main NRVS sample at stage A and the calibration sample at stage B simultaneously, and calibrates the energies for observing extremely small spectral shifts; for example, the 0.3 meV energy shift between the 100%-(57)Fe-enriched [Fe4S4Cl4](=) and 10%-(57)Fe and 90%-(54)Fe labeled [Fe4S4Cl4](=) has been well resolved; (ii) a quick-switching energy calibration procedure, which reduces each calibration time from 3-4 h to about 30 min. Although the quick-switching calibration is not in situ, it is suitable for normal NRVS measurements.

Characterization of [4Fe-4S] cluster vibrations and structure in nitrogenase Fe protein at three oxidation levels via combined NRVS, EXAFS, and DFT analyses

Azotobacter vinelandii nitrogenase Fe protein (Av2) provides a rare opportunity to investigate a [4Fe-4S] cluster at three oxidation levels in the same protein environment. Here, we report the structural and vibrational changes of this cluster upon reduction using a combination of NRVS and EXAFS spectroscopies and DFT calculations. Key to this work is the synergy between these three techniques as each generates highly complementary information and their analytical methodologies are interdependent. Importantly, the spectroscopic samples contained no glassing agents. NRVS and DFT reveal a systematic 10-30 cm(-1) decrease in Fe-S stretching frequencies with each added electron. The "oxidized" [4Fe-4S](2+) state spectrum is consistent with and extends previous resonance Raman spectra. For the "reduced" [4Fe-4S](1+) state in Fe protein, and for any "all-ferrous" [4Fe-4S](0) cluster, these NRVS spectra are the first available vibrational data. NRVS simulations also allow estimation of the vibrational disorder for Fe-S and Fe-Fe distances, constraining the EXAFS analysis and allowing structural disorder to be estimated. For oxidized Av2, EXAFS and DFT indicate nearly equal Fe-Fe distances, while addition of one electron decreases the cluster symmetry. However, addition of the second electron to form the all-ferrous state induces significant structural change. EXAFS data recorded to k = 21 Å(-1) indicates a 1:1 ratio of Fe-Fe interactions at 2.56 Å and 2.75 Å, a result consistent with DFT. Broken symmetry (BS) DFT rationalizes the interplay between redox state and the Fe-S and Fe-Fe distances as predominantly spin-dependent behavior inherent to the [4Fe-4S] cluster and perturbed by the Av2 protein environment.

Nuclear resonance vibrational spectroscopy and electron paramagnetic resonance spectroscopy of 57Fe-enriched [FeFe] hydrogenase indicate stepwise assembly of the H-cluster

The [FeFe] hydrogenase from Clostridium pasteurianum (CpI) harbors four Fe-S clusters that facilitate the transfer of an electron to the H-cluster, a ligand-coordinated six-iron prosthetic group that catalyzes the redox interconversion of protons and H(2). Here, we have used (57)Fe nuclear resonance vibrational spectroscopy (NRVS) to study the iron centers in CpI, and we compare our data to that for a [4Fe-4S] ferredoxin as well as a model complex resembling the [2Fe](H) catalytic domain of the H-cluster. To enrich the hydrogenase with (57)Fe nuclei, we used cell-free methods to post-translationally mature the enzyme. Specifically, inactive CpI apoprotein with (56)Fe-labeled Fe-S clusters was activated in vitro using (57)Fe-enriched maturation proteins. This approach enabled us to selectively label the [2Fe](H) subcluster with (57)Fe, which NRVS confirms by detecting (57)Fe-CO and (57)Fe-CN normal modes from the H-cluster nonprotein ligands. The NRVS and iron quantification results also suggest that the hydrogenase contains a second (57)Fe-S cluster. Electron paramagnetic resonance (EPR) spectroscopy indicates that this (57)Fe-enriched metal center is not the [4Fe-4S](H) subcluster of the H-cluster. This finding demonstrates that the CpI hydrogenase retained an (56)Fe-enriched [4Fe-4S](H) cluster during in vitro maturation, providing unambiguous evidence of stepwise assembly of the H-cluster. In addition, this work represents the first NRVS characterization of [FeFe] hydrogenases.

Observation of the Fe-CN and Fe-CO vibrations in the active site of [NiFe] hydrogenase by nuclear resonance vibrational spectroscopy

Nuclear inelastic scattering of (57)Fe labeled [NiFe] hydrogenase is shown to give information on different states of the enzyme. It was thus possible to detect and assign Fe-CO and Fe-CN bending and stretching vibrations of the active site outside the spectral range of the Fe-S cluster normal modes.

Effects of quantum coherence in metalloprotein electron transfer

Many intramolecular electron transfer (ET) reactions in biology are mediated by metal centers in proteins. This process is commonly described by a model of diffusive hopping according to the semiclassical theories of Marcus and Hopfield. However, recent studies have raised the possibility that nontrivial quantum mechanical effects play a functioning role in certain biomolecular processes. Here, we investigate the potential effects of quantum coherence in biological ET by extending the semiclassical model to allow for the possibility of quantum coherent phenomena using a quantum master equation based on the Holstein Hamiltonian. We test the model on the structurally defined chain of seven iron-sulfur clusters in nicotinamide adenine dinucleotide plus hydrogen:ubiquinone oxidoreductase (complex I), a crucial respiratory enzyme and one of the longest chains of metal centers in biology. Using experimental parameters where possible, we find that, in limited circumstances, a small quantum mechanical contribution can provide a marked increase in the ET rate above the semiclassical diffusive-hopping rate. Under typical biological conditions, our model reduces to well-known diffusive behavior.

Force-Field Development and Molecular Dynamics of [NiFe] Hydrogenase

Classical molecular force-field parameters describing the structure and motion of metal clusters in [NiFe] hydrogenase enzymes can be used to compare the dynamics and thermodynamics of [NiFe] under different oxidation, protonation, and ligation circumstances. Using density functional theory (DFT) calculations of small model clusters representative of the active site and the proximal, medial, and distal Fe/S metal centers and their attached protein side chains, we have calculated classical force-field parameters for [NiFe] in reduced and oxidized states, including internal coordinates, force constants, and atom-centered charges. Derived force constants revealed that cysteinate ligands bound to the metal ions are more flexible in the Ni-B active site, which has a bridging hydroxide ligand, than in the Ni-C active site, which has a bridging hydride. Ten nanosecond all-atom, explicit-solvent MD simulations of [NiFe] hydrogenase in oxidized and reduced catalytic states established the stability of the derived force-field parameters in terms of Cα and metal cluster fluctuations. Average active site structures from the protein MD simulations are consistent with [NiFe] structures from the Protein Data Bank, suggesting that the derived force-field parameters are transferrable to other hydrogenases beyond the structure used for testing. A comparison of experimental H2-production rates demonstrated a relationship between cysteinate side chain rotation and activity, justifying the use of a fully dynamic model of [NiFe] metal cluster motion.

Mössbauer spectroscopy on respiratory complex I: the iron-sulfur cluster ensemble in the NADH-reduced enzyme is partially oxidized

In mitochondria, complex I (NADH:quinone oxidoreductase) couples electron transfer to proton translocation across an energy-transducing membrane. It contains a flavin mononucleotide to oxidize NADH, and an unusually long series of iron–sulfur (FeS) clusters that transfer the electrons to quinone. Understanding electron transfer in complex I requires spectroscopic and structural data to be combined to reveal the properties of individual clusters and of the ensemble. EPR studies on complex I from Bos taurus have established that five clusters (positions 1, 2, 3, 5, and 7 along the seven-cluster chain extending from the flavin) are (at least partially) reduced by NADH. The other three clusters, positions 4 and 6 plus a cluster on the other side of the flavin, are not observed in EPR spectra from the NADH-reduced enzyme: they may remain oxidized, have unusual or coupled spin states, or their EPR signals may be too fast relaxing. Here, we use Mössbauer spectroscopy on ⁵⁷Fe-labeled complex I from the mitochondria of Yarrowia lipolytica to show that the cluster ensemble is only partially reduced in the NADH-reduced enzyme. The three EPR-silent clusters are oxidized, and only the terminal 4Fe cluster (position 7) is fully reduced. Together with the EPR analyses, our results reveal an alternating profile of higher and lower potential clusters between the two active sites in complex I; they are not consistent with the consensus picture of a set of isopotential clusters. The implications for intramolecular electron transfer along the extended chain of cofactors in complex I are discussed.