Proton Transfer Pathways between Active Sites and Proximal Clusters in the Membrane-Bound [NiFe] Hydrogenase

Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H2 gas powered by a [NiFe] active site and transduces the free energy into proton pumping and Na+/H+ exchange across the membrane. Despite recent structural advances, the mechanistic principles of H2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H2 and how the catalysis couples to conformational dynamics in the membrane domain of Mbh. By combining large-scale quantum chemical density functional theory (DFT) and correlated ab initio wave function methods with atomistic molecular dynamics simulations, we show that the proton transfer reactions required for the catalysis are gated by electric field effects that direct the protons by water-mediated reactions from Glu21L toward the [NiFe] site, or alternatively along the nearby His75L pathway that also becomes energetically feasible in certain reaction steps. These local proton-coupled electron transfer (PCET) reactions induce conformational changes around the active site that provide a key coupling element via conserved loop structures to the ion transport activity. We find that H2 forms in a heterolytic proton reduction step, with spin crossovers tuning the energetics along key reaction steps. On a general level, our work showcases the role of electric fields in enzyme catalysis and how these effects are employed by the [NiFe] active site of Mbh to drive PCET reactions and ion transport.

Conformational and mechanical stability of the isolated large subunit of membrane-bound [NiFe]-hydrogenase from Cupriavidus necator

Comprising at least a bipartite architecture, the large subunit of [NiFe]-hydrogenase harbors the catalytic nickel-iron site while the small subunit houses an array of electron-transferring Fe-S clusters. Recently, some [NiFe]-hydrogenase large subunits have been isolated showing an intact and redox active catalytic cofactor. In this computational study we have investigated one of these metalloproteins, namely the large subunit HoxG of the membrane-bound hydrogenase from Cupriavidus necator (CnMBH), targeting its conformational and mechanical stability using molecular modelling and long all-atom Gaussian accelerated molecular dynamics (GaMD). Our simulations predict that isolated HoxG is stable in aqueous solution and preserves a large portion of its mechanical properties, but loses rigidity in regions around the active site, in contrast to the MBH heterodimer. Inspired by biochemical data showing dimerization of the HoxG protein and IR measurements revealing an increased stability of the [NiFe] cofactor in protein preparations with higher dimer content, corresponding simulations of homodimeric forms were also undertaken. While the monomeric subunit contains several flexible regions, our data predicts a regained rigidity in homodimer models. Furthermore, we computed the electrostatic properties of models obtained by enhanced sampling with GaMD, which displays a significant amount of positive charge at the protein surface, especially in solvent-exposed former dimer interfaces. These data offer novel insights on the way the [NiFe] core is protected from de-assembly and provide hints for enzyme anchoring to surfaces, which is essential information for further investigations on these minimal enzymes.

Proton Transfer Mechanisms in Bimetallic Hydrogenases

Hydrogenases are metalloenzymes that catalyze proton reduction and H₂ oxidation with outstanding efficiency. They are model systems for bioinorganic chemistry, including low-valent transition metals, hydride chemistry, and proton-coupled electron transfer. In this Account, we describe how photochemistry and infrared difference spectroscopy can be used to identify the dynamic hydrogen-bonding changes that facilitate proton transfer in [NiFe]- and [FeFe]-hydrogenase.[NiFe]-hydrogenase binds a heterobimetallic nickel/iron site embedded in the protein by four cysteine ligands. [FeFe]-hydrogenase carries a homobimetallic iron/iron site attached to the protein by only a single cysteine. Carbon monoxide and cyanide ligands in the active site facilitate detailed investigations of hydrogenase catalysis by infrared spectroscopy because of their strong signals and redox-dependent frequency shifts. We found that specific redox-state transitions in [NiFe]- and [FeFe]-hydrogenase can be triggered by visible light to record extremely sensitive "light-minus-dark" infrared difference spectra monitoring key amino acid residues. As these transitions are coupled to protonation changes, our data allowed investigation of dynamic hydrogen-bonding changes that go well beyond the resolution of protein crystallography.In [NiFe]-hydrogenase, photolysis of the bridging hydride ligand in the Ni-C state was followed by infrared difference spectroscopy. Our data clearly indicate the formation of a protonated cysteine residue as well as hydrogen-bonding changes involving a glutamic acid residue and a "dangling water" molecule. These findings are in excellent agreement with crystallographic analyses of [NiFe]-hydrogenase. In [FeFe]-hydrogenase, an external redox dye was used to accumulate the Hred state. Infrared difference spectra indicate hydrogen-bonding changes involving two glutamic acid residues and a conserved arginine residue. While crystallographic analyses of [FeFe]-hydrogenase in the oxidized state failed to explain the rapid proton transfer because of a breach in the succession of residues, our findings facilitated a precise molecular model of discontinued proton transfer.Comparing both systems, our data emphasize the role of the outer coordination sphere in bimetallic hydrogenases: we suggest that protonation of a nickel-ligating cysteine in [NiFe]-hydrogenase causes the notable preference toward H₂ oxidation. On the contrary, proton transfer in [FeFe]-hydrogenase involves an adjacent cysteine as a relay group, promoting both H₂ oxidation and proton reduction. These observations may guide the design of organometallic compounds that mimic the catalytic properties of hydrogenases.

Hydroxy-bridged resting states of a [NiFe]-hydrogenase unraveled by cryogenic vibrational spectroscopy and DFT computations

The catalytic mechanism of [NiFe]-hydrogenases is a subject of extensive research. Apart from at least four reaction intermediates of H₂/H⁺ cycling, there are also a number of resting states, which are formed under oxidizing conditions. Although not directly involved in the catalytic cycle, the knowledge of their molecular structures and reactivity is important, because these states usually accumulate in the course of hydrogenase purification and may also play a role in vivo during hydrogenase maturation. Here, we applied low-temperature infrared (cryo-IR) and nuclear resonance vibrational spectroscopy (NRVS) to the isolated catalytic subunit (HoxC) of the heterodimeric regulatory [NiFe]-hydrogenase (RH) from Ralstonia eutropha. Cryo-IR spectroscopy revealed that the HoxC protein can be enriched in almost pure resting redox states suitable for NRVS investigation. NRVS analysis of the hydrogenase catalytic center is usually hampered by strong spectral contributions of the FeS clusters of the small, electron-transferring subunit. Therefore, our approach to investigate the FeS cluster-free, ⁵⁷Fe-labeled HoxC provided an unprecedented insight into the [NiFe] site modes, revealing their contributions in a spectral range otherwise superimposed by FeS cluster-derived bands. Rationalized by density functional theory (DFT) calculations, our data provide structural descriptions of the previously uncharacterized hydroxy- and water-containing resting states. Our work highlights the relevance of cryogenic vibrational spectroscopy and DFT to elucidate the structure of barely defined redox states of the [NiFe]-hydrogenase active site.

Active site vibrations of a [NiFe]-hydrogenase catalytic subunit are selectively probed by IR and NRV spectroscopy in two Ni^IIFe^II and Ni^IIIFe^II resting states, contributing in combination with DFT modeling to rationalized structural candidates.

Probing the Structure of [NiFeSe] Hydrogenase with QM/MM Computations

The geometry and vibrational behavior of selenocysteine [NiFeSe] hydrogenase isolated from Desulfovibrio vulgaris Hildenborough have been investigated using a hybrid quantum mechanical (QM)/ molecular mechanical (MM) approach. Structural models have been built based on the three conformers identified in the recent crystal structure resolved at 1.3 Å from X-ray crystallography. In the models, a diamagnetic Ni2+ atom was modeled in combination with both Fe2+ and Fe3+ to investigate the effect of iron oxidation on geometry and vibrational frequency of the nonproteic ligands, CO and CN-, coordinated to the Fe atom. Overall, the QM/MM optimized geometries are in good agreement with the experimentally resolved geometries. Analysis of computed vibrational frequencies, in comparison with experimental Fourier-transform infrared (FTIR) frequencies, suggests that a mixture of conformers as well as Fe2+ and Fe3+ oxidation states may be responsible for the acquired vibrational spectra.

Cysteine SH and Glutamate COOH Contributions to [NiFe] Hydrogenase Proton Transfer Revealed by Highly Sensitive FTIR Spectroscopy

A [NiFe] hydrogenase (H₂ ase) is a proton-coupled electron transfer enzyme that catalyses reversible H₂ oxidation; however, its fundamental proton transfer pathway remains unknown. Herein, we observed the protonation of Cys546-SH and Glu34-COOH near the Ni-Fe site with high-sensitivity infrared difference spectra by utilizing Ni-C-to-Ni-L and Ni-C-to-Ni-SI_a photoconversions. Protonated Cys546-SH in the Ni-L state was verified by the observed SH stretching frequency (2505 cm^-1 ), whereas Cys546 was deprotonated in the Ni-C and Ni-SI_a states. Glu34-COOH was double H-bonded in the Ni-L state, as determined by the COOH stretching frequency (1700 cm^-1 ), and single H-bonded in the Ni-C and Ni-SI_a states. Additionally, a stretching mode of an ordered water molecule was observed in the Ni-L and Ni-C states. These results elucidate the organized proton transfer pathway during the catalytic reaction of a [NiFe] H₂ ase, which is regulated by the H-bond network of Cys546, Glu34, and an ordered water molecule.