How [FeFe]-Hydrogenase Facilitates Bidirectional Proton Transfer

Hydrogenases are metalloenzymes that catalyze the conversion of protons and molecular hydrogen, H2. [FeFe]-hydrogenases show particularly high rates of hydrogen turnover and have inspired numerous compounds for biomimetic H2 production. Two decades of research on the active site cofactor of [FeFe]-hydrogenases have put forward multiple models of the catalytic proceedings. In comparison, our understanding of proton transfer is poor. Previously, residues were identified forming a hydrogen-bonding network between active site cofactor and bulk solvent; however, the exact mechanism of catalytic proton transfer remained inconclusive. Here, we employ in situ infrared difference spectroscopy on the [FeFe]-hydrogenase from Chlamydomonas reinhardtii evaluating dynamic changes in the hydrogen-bonding network upon photoreduction. While proton transfer appears to be impaired in the oxidized state (Hox), the presented data support continuous proton transfer in the reduced state (Hred). Our analysis allows for a direct, molecular unique assignment to individual amino acid residues. We found that transient protonation changes of glutamic acid residue E141 and, most notably, arginine R148 facilitate bidirectional proton transfer in [FeFe]-hydrogenases.


■ INTRODUCTION
Hydrogenases are gas-processing iron−sulfur enzymes that catalyze the reversible reduction of protons to molecular hydrogen in all kingdoms of life. 1,2 Most hydrogenases are biased toward H 2 oxidation, for example, in the context of energy metabolism and H 2 sensing. 3−5 The [FeFe]-hydrogenases from bacteria and algae, in contrast, are truly bidirectional and catalyze H 2 oxidation and H 2 evolution with similar efficiency. 6−8 Combining high turnover frequencies (10 000 H 2 s −1 ) and a catalytic midpoint potential close to the H + /H 2 redox couple, 9−11 the active site cofactor of [FeFe]hydrogenases ("H-cluster") inspired the design of numerous biomimetic complexes for H 2 production. 12−14 The H-cluster comprises a conventional [4Fe-4S] center linked to a bimetallic iron−sulfur complex (Figure 1a). 15−17 The diiron site carries two terminal carbonyl and cyanide ligands (CO, CN − ) as well as a single carbonyl ligand in Fe− Fe bridging position (μCO). 18−20 An aminodithiolate (ADT) group connects the proximal and distal iron ion (Fe p and Fe d , relative to the [4Fe-4S] center) 21 and functions as proton relay between active site cofactor and protein environment. 22 While prokaryotic [FeFe]-hydrogenases like CPI from Clostridium pasteurianum and DDH from Desulfovibrio desulf uricans hydrogenase) bind additional iron−sulfur clusters, the enzyme from Chlamydomonas reinhardtii (HYDA1) exclusively carries the H-cluster. 6 During hydrogen turnover, the H-cluster adopts different redox and protonation states. The oxidized resting state (Hox) 23−25 can be distinguished from intermediates with a reduced [4Fe-4S] center (Hred′, Hhyd) 26−32 or a reduced diiron site (Hred, Hsred). 33−35 These states are formed upon concerted proton and electron transfer. The active-ready geometry of Hox is characterized by a square-pyramidal configuration of both metal ions, a μCO ligand, and an open coordination site at Fe d . 15−17 While this geometry is conserved in Hred′ and Hhyd, 27,30 the structural changes upon reduction of the diiron site are under debate. 36 The H-cluster may undergo rigorous ligand rearrangement forming a μH geometry, which would exclude both Hred and Hsred from catalytic turnover. 35 Alternatively, diiron site geometries with a bridging 37,38 or "semi-bridging" CO ligand 20 have been suggested. On the basis of the assumption of a protonated ADT group in these species, the nomenclature HredH + and HsredH + can be found in literature. 34 [FeFe]-hydrogenases exchange protons with the bulk solvent via a trajectory of conserved, polar amino acid residues that connect active site cofactor and protein surface ( Figure  1a). 39−41 Making use of protein crystallography and infrared spectroscopy, previously we were able to identify the residues that render catalytic proton transfer possible. 42 In the [FeFe]hydrogenase from C. reinhardtii, this includes R148, E144, S189, E141, C169 (numbers 1−5 in Figure 1a), and water cluster W1. A second, mostly aqueous proton transfer pathway facilitates protonation of the [4Fe-4S] center. 27 These protons are not consumed in the H 2 release reaction but stabilize the active-ready geometry and compensate for the drop in redox potential after a first reduction step. 26,32 The dynamics of catalytic proton transfer have been addressed by molecular dynamics simulations before; 43−45 however, no experimental data exist on the changes in the hydrogen-bonding network, for example, when switching from H 2 evolution to H 2 oxidation. Such data are key to understanding the bidirectional catalysis of [FeFe]-hydrogenases.
The high extinction coefficient of the active site CO and CN − ligands facilitated numerous investigations of the Hcluster by Fourier-transform infrared (FTIR) spectroscopy. 25−34 In chromophoric proteins, FTIR difference spectroscopy is routinely used for the analysis of light-induced proton transfer steps. 46−48 However, compared to the H-cluster ligands, the involved amino acid side chains exhibit a drastically lower extinction coefficient and overlap with the intense absorption bands of liquid water and the protein backbone.
[FeFe]-hydrogenases lack a natural chromophore to selectively trigger catalytic activity and protonation changes by light, which renders an analysis of proton transfer challenging. Redox dyes provides an opportunity to characterize hydrogenbonding networks in visibly transparent enzymes. 49 Here, we explore the dynamics of proton transfer in HYDA1 using in situ attenuated total reflection (ATR) FTIR difference spectroscopy. We report steady-state photoreduction of the Hcluster in the presence of 5′-carboxy eosin Y (5CE) 50 that resulted in a fast (t 1/2 ≈ 20 s), near complete (>90%), and highly selective redox conversion (<5% other species) of the one-electron reduced state Hred over the oxidized resting state, Hox (Figure 1b). These parameters were not achieved by in situ gas treatments 27,30 or electrochemical titrations. 32,35 The proton uptake associated with formation of Hred induced spectral differences in the IR regime from 1750 to 1650 cm −1 . Exploiting in situ H/D exchange and site-directed mutagenesis, these differences are assigned to the CO stretching vibrations of carboxylic acid side chains (COOH) and the coupled vibrational mode of an protonated arginine side chain (C(NH 2 ) 3 + ). Infrared spectroscopy provides evidence for changes in hydrogen bonding involving glutamic acid E141 and serine S189 close to the active site as well as glutamic acid E144 and arginine R148 near the protein surface. This work presents the first direct, experimental characterization of the hydrogen-bonding network that facilitates catalytic proton transfer in [FeFe]-hydrogenases.

■ EXPERIMENTAL SECTION
The [FeFe]-hydrogenase from C. reinhardtii HYDA1 was expressed and synthesized in Escherichia coli, purified by strep-tactin affinity chromatography, and activated in vitro with synthetic ADT-containing diiron complex under anaerobic conditions. 51,52 After removal of excess complex, the protein concentration was adjusted to ∼3 mM (∼150 g/L). 5′-Carboxy eosin Y (5CE) 50 and ethylenediaminetetraacetic acid (EDTA) were prepared in aqueous stock solutions of 6 and 90 mM, respectively. One part of each component was mixed to yield a HYDA1/5CE/EDTA ratio of 1:2:30.
All spectroscopic experiments were performed under anaerobic conditions, at room temperature, ambient pressure, and on hydrated protein films of physiological pH values. First, 1−2 μL of the reaction mix was pipetted onto the silicon crystal of the ATR unit (DuraDisc SamplIR-2, Smiths Detection) in the FTIR spectrometer (Tensor 27, Bruker). Spectra from 3900−1300 cm −1 were recorded with a narrowband mercury cadmium telluride (MCT) detector with a spectral resolution of 2 cm −1 and 25 interferometer scans each (80 kHz). The solution was protected from stray light, dried under N 2 , and rehydrated via the gas phase with 10 mM 2-(N-morpholino)ethanesulfonic acid (MES) buffer (pH 6). Traces of reduced and CO-inhibited species were lost in favor of Hox upon auto-oxidation. 27 For the pH jump experiments, guanidine-HCl was solved in H 2 O or D 2 O. Ten microliters of a 1000 g L −1 guanidine-HCl solution (pH ≈∼ 8) was pipetted onto the ATR crystal and continuously measured by FTIR spectroscopy. Addition of 1 μL of NaOH solution (∼100 g L −1 , in H 2 O or D 2 O) increased the bulk pH above the pK a of guanidine-HCl (∼13.5) deprotonating the guanidinium ion. Steady-state illumination of 5′-carboxy eosin Y (5CE, photosensitizer) at 505 nm populated the excited triplet state, 5CE*. The latter was quenched by the oxidized enzyme to enrich Hred over Hox. Increasing the basicity upon reduction induces a protonation of the H-cluster. The associated changes in the hydrogen-bonding network of the catalytic proton transfer pathway are followed by in situ ATR FTIR difference spectroscopy. The sacrificial electron donor EDTA rereduced 5CE + to 5CE.

■ RESULTS
The experiment was initiated upon steady-state illumination of the film at 505 nm and followed by FTIR spectroscopy with a time resolution of 5 s (up to 60 s, see Figure S1). After ∼20 s, half of the Hox population converted into Hred. Importantly, no other reduced species than Hred were observed, although H 2 released upon reduction of HYDA1 50 may rereact with the enzyme forming Hred′, Hsred, and Hhyd. Therefore, the continuous exchange of gas in our setup was found to be of particular importance, as it precludes reoxidation of H 2 and a buildup of multiple reduced species. 22 No photoreduction was observed under off-resonant conditions (590 nm), whereas illumination at 455 nm induced notable H-cluster corruption ( Figure S2), as noted earlier. 53 Figure 2a shows an overlay of absorbance spectra in the range of 3900−1300 cm −1 in the dark (black) and light (magenta). The absorbance ratio of ∼1.3 for amide I (1635 cm −1 ) to amide II (1545 cm −1 ) suggests a well-hydrated protein film. 27 From 2700 to 1800 cm −1 neither liquid water (H 2 O) nor protein solution show strong IR intensities, which allows analyzing the CO/CN − bands of the H-cluster in absolute spectra. The inset highlights the IR signature of the H-cluster from 2150 to 1750 cm −1 . The difference spectrum in Figure 2b emphasizes how the hydrogenase adopted Hox in the dark while Hred clearly dominated upon illumination. The respective IR band patterns have been identified earlier. 34,35 Fitting these patterns to absolute spectra before and after illumination indicated a redox conversion larger than 90%.
In contrast to the CO/CN − bands, the intense absorbance of liquid water (HOH bending) and protein backbone (amide I, amide II) overlaps with signals in the COOH regime from 1750 to 1650 cm −1 and precludes any meaningful analysis in absolute spectra. Figure 2c shows a "dark−light" difference spectrum computed from the single channel spectra that generated the absorbance spectra in Figure 2a. The cofactor bands clearly dominate the spectrum. Efficient photoreduction prevents an accumulation of unspecific changes in the film (i.e., hydration level, protein concentration) and allows analyzing the full spectrum. This includes the OH, SH, and COOH regime as well as frequencies less than 1600 cm −1 comprising marker bands of the photosensitizer 5CE ( Figure S2).
All difference bands in the COOH regime are specific for functional HYDA1. The inset in Figure 2c shows that no such changes were observed when HYDA1 apoprotein was probed (apo-HYDA1 lacks the diiron site and is catalytically unreactive 51,52 ). Moreover, difference spectra of HYDA1 recorded upon exposure to CO (Hox-CO over Hox) or in the presence of zinc porphyrin as an alternative redox dye 54 (Hred′ over Hox) confirmed that all bands in the COOH regime are specific for the formation of Hred ( Figure S3).
It is important to point out that our results do not suggest protonation or hydrogen-bonding differences involving OH or SH groups ( Figure S4). The former would give rise to sharp absorbance bands around 3650 cm −1 indicative of "dangling", weakly hydrogen-bonded water. 55,56 The SH group absorbs around 2550 cm −1 and is very sensitive to changes in hydrophilicity. 56−58 Typically, this frequency regime is addressed to analyze hydrogen-bonding changes involving the side chain of a cysteine, for example, C169. Comparing our data with recent work by Hirota et al. on [NiFe]-hydrogenases 58 shows that the signal-to-noise ratio of the Hred − Hox difference spectrum would allow identifying potential changes in the OH and SH regime ( Figure S4). However, the lack thereof suggests an invariable hydrogen-bonding network between the water cluster W1, C169, and the ADT headgroup of the H-cluster (Figure 1a).
Band Fitting and Tentative Assignments. The Hred − Hox spectrum in the COOH regime was described by a fit routine including a minimum of nine Gaussians with a fixed half-max width of 6−8 cm −1 and third-order polynomial baseline correction (Figure 3a). The temporal evolution of these bands is in excellent agreement with those of the H- (c) The full difference spectrum allows analyzing the OH, SH, and COOH regime as well as frequencies less than 1500 cm −1 comprising vibrational marker bands of photosensitizer 5CE. (inset) Magnification of the COOH regime. The band changes are specific for functional HYDA1 (black) and not observed in HYDA1 apoprotein (red). * 2337 cm −1 , assigned to CO 2 .

Journal of the American Chemical Society
Article cluster ( Figure S1). Vibrations at frequencies greater than 1700 cm −1 are typically assigned to the CO stretches of the COOH side chain from aspartic acid or glutamic acid residues (E, D). 57,59 The CO stretching frequency is inversely proportional to the hydrogen-bonding strength and can vary between none to multiple hydrogen-bonding partners from 1750 to 1700 cm −1 . Besides clearly discriminable bands at 1721, 1715, and 1700 cm −1 our fit routine suggested additional contributions centered at 1710, 1696, 1690, and 1681 ( Figure  3a). The negative band at 1681 cm −1 may be attributed to the asymmetric C(NH 2 ) 3 + vibration of the protonated arginine side chain, 60−62 while the broader features at ∼1670 and 1655 cm −1 potentially arise from changes in amide I absorbance. 63 The latter likely reflects minor changes in secondary structure induced upon reduction of the H-cluster.
H/D Exchange. To achieve an experimental band assignment in the COOH regime, we performed photoreduction on hydrated and deuterated hydrogenases films. Bands indicative of hydrogen bonding or protonation changes involving the carboxylic side chains are supposed to shift to lower frequencies in deuterated sample. 63 Absorbance spectra of the HYDA1/5CE/EDTA reaction mixture show a complete exchange of solvent in the presence of either H 2 O or D 2 O, and Hred − Hox difference spectra prove that deuteration did not affect the H-clusters' CO/CN − band position ( Figure S5).
However, in the COOH regime, the spectra show significant changes (Figure 3b). The prominent H/D shift of 1715 and 1700 cm −1 to 1709 and 1694 cm −1 immediately supports an assignment of this motif to a titratable group, for example, an aspartic or glutamic acid side chain. Bands at 1696 and 1690 cm −1 were affected by the H/D shift as well. While a dissection of components is not immediately possible here, the mean frequency downshift by 12 ± 2 cm −1 suggests strongly hydrogen-bonded aspartic or glutamic acid side chains. 57,59 The positive band at 1721 cm −1 is insensitive to H/D exchange.
To achieve an unambiguous experimental band assignment, we analyzed three amino acid variants of the proton transfer pathway, namely, R148A, E144D, and S189A. The enrichment of Hred over Hox depends on functional proton transfer. 33−35 In particular, amino acid residues C169 and E141 close to the H-cluster were found to be susceptible to variations of the hydrogen-bonding network, slowly accumulating the hydride state Hhyd over Hox rather than Hred ( Figure S6). This impedes a direct comparison, and only a limited number of variants allowed screening the hydrogen-bonding changes associated with catalytic proton transfer. An invariable hydrogen-bonding network between H-cluster and E141 is in striking agreement with the aforementioned lack of hydrogenbonding changes around W1 and C169 ( Figure S4). Relative to the H-cluster, we will refer to C169, W1, and E141 as "inner core" of the proton transfer pathway.
By contrast, Figure 4 shows Hred − Hox difference spectra of HYDA1 variants that constitute the "outer core" of the proton transfer pathway, namely, R148, E144, and S189. Sitedirected mutagenesis at these positions included an exchange against alanine (A). No accumulation of Hred was observed upon photoreduction of E144A ( Figure S6) so that the conservative variation of glutamic to aspartic acid was analyzed instead. Variants R148A, E144D, and S189A adopted Hred upon illumination but showed only 15−25% of the native conversion efficiency ( Figure S7). For comparison with native HYDA1, difference spectra were normalized to the amplitude of the band pair at 1715 and 1700 cm −1 that was found to be prominently conserved in all spectra. The resulting scaling factors were in good agreement with the amplitudes observed for the CO difference bands of the conversion of Hred over Hox. Figure S6 depicts the spectral transitions over time for each variant including an evaluation of signal-to-noise in the COOH regime.
Amino Acid Variant R148A. The H 2 evolution activity of ∼50% for HYDA1 variant R148A indicates that glutamic acid E144 can partially replace R148 as proton loading site. 42 Figure  4a shows an overlay of Hred − Hox difference spectra in the COOH regime for R148A and native HYDA1 (left panel). The right panel depicts an overlay of the respective CPI crystal structures. For the sake of convenience, we will use HYDA1 numbering here. Site-directed mutagenesis resulted in spectra with missing features at 1721, 1696, or 1681 cm −1 (marked "X"), while the shift from 1715 to 1700 cm −1 and a negative band at 1690 cm −1 was conserved.
Above, we tentatively assigned the negative feature at 1681 cm −1 to the asymmetric C(NH 2 ) 3 + vibration of an arginine residue. 60−62 The evident lack of this band in amino acid variant R148A supports this assignment. Moreover, the band at 1681 cm −1 shifted to 1607 cm −1 in deuterated sample, which is in excellent agreement to guanidine hydrochloride H/D reference spectra in Figure 5. Therefore, we conclude

Journal of the American Chemical Society
Article deprotonation of R148 + upon formation of Hred in native HYDA1. Poisson−Boltzmann calculations predicted a pK a of ∼3.5 for E144 in native HYDA1 (Table S1). Although our experiments were conducted at pH 6, we suggest hydrogen bonding of E144 to R148 (∼2.8 Å) and S189 (∼3.1 Å) stabilizing the carboxylic over the carboxylate form of E144. However, the carboxylate form likely prevails in the absence of the arginine side chain. The lack of spectral features at 1721 and 1696 cm −1 in the R148A difference spectrum (X) therefore facilitates the assignment to hydrogen-bonding changes involving E144 in native HYDA1.
Glutamic Acid Variant E144D. As observed for the arginine variant, the interaction between R148 and E144 (E144D) is not strictly essential for catalytic activity. Amino acid variant E144D is reported with ∼50% H 2 evolution activity. 42 Figure 4b shows an overlay of Hred − Hox difference spectra in the COOH regime for R144D and native HYDA1 (left panel). The right panel depicts an overlay of the respective CPI crystal structures. 42 Site-directed mutagenesis resulted in spectra with a pronounced band upshift from 1721 to 1742 cm −1 and an intensity inversion of the band at ∼1690 cm −1 (negative in native HYDA1, slightly shifted and positive in E144D). The band pair at 1715 and 1700 cm −1 and negative bands at 1694 and 1681 cm −1 are conserved among variant and native HYDA1.
Shortening of the alkyl side chain at position 144 causes a different hydrogen-bonding situation (Figure 4b, right panel). Instead of forming a hydrogen bond with R148 (∼5.0 Å) rotation of the aspartic acid side chain forces D144 into a weak complex with S189 in Hox (3.0 and 3.7 Å) that reflects in a pK a increase of nearly three units compared to native HYDA1 (Table S1). The upshift of the E144 band from 1721 to 1742 cm −1 suggests significantly weaker hydrogen bonding in reduced enzyme. Furthermore, variant E144D allows differentiating the 1696/1690 cm −1 peak doublet. The latter band appears positive in the spectrum, and thus only the 1696 cm −1 band is assigned to E144.
Serine Variant S189A. Figure 4c shows an overlay of Hred − Hox difference spectra in the COOH regime for S189A and native HYDA1 (left panel). The right panel depicts an overlay of the respective CPI crystal structures. 42 Despite the relatively low H 2 evolution activity of ∼10%, site-directed mutagenesis resulted in spectra indicative of only minor differences to native HYDA1. The crystal structure of the S189A variant revealed an additional water molecule (W*) between E144 and E141 (Figure 4c, right panel). This arrangement was proposed to compensate the lack of the serine side chain. 42 E144 is in fair hydrogen-bonding distance to W* (∼2.4 Å) and R148 (∼2.9 Å), which largely restores the spectral phenotype of native HYDA1. However, the distance of 5.8 Å between E141 and W* is clearly out of range for hydrogen bonding or proton transfer.
Glutamic Acid E141. The band pair at 1715 and 1700 cm −1 is prominently conserved in all protein samples that accumulate Hred over Hox. The H/D specific band shift hints at a carboxylic group and suggests efficient proton exchange;

Journal of the American Chemical Society
Article however, the motif could not be assigned to E144. Residing in a hydrophobic pocket at the interface of inner and outer core of the proton transfer pathway, glutamic acid E141 has been calculated to adopt the carboxylic acid form for pH values less than 8. 42 Accordingly, changes associated with E141 will be visible in the COOH regime. 57 Any variation of E141 abolished catalytic activity and the formation of Hred ( Figure  S5) hinting at the central role of E141 in proton transfer. On the basis of this line of evidence, we assign the band pair at 1715 and 1700 cm −1 to E141. Figure S7 provides a conclusive overview on the observed frequencies and experiment band assignment.
The CO stretching frequencies of glutamic acid E141 indicate strong hydrogen-bonding contacts, irrespective of redox state. 57,59 To this end, the crystal structure of oxidized enzyme supports a trans complex between E141 and W1 (distances 2.4 and 3.4 Å). 64−66 The 15 cm −1 frequency downshift upon reduction may reflect a release of the E141/ W1 complex in favor of hydrogen bonding with S189 (or W* in serine variant S189A). This demands a certain level of structural flexibility, as the distance between E141 and S189 amounts to ∼3.8 Å in oxidized enzyme. Molecular dynamics simulations showed that E141 and S189 change between hydrogen-bonding donor and acceptor when switching from proton uptake to proton release. 43,45 Apparently, smaller structural changes at the interface of inner and outer core are well within the thermodynamic range of functional [FeFe]hydrogenases. The large distance between E141 and W* (5.8 Å) reduces the probability of proton transfer in serine variant S189A to ∼10% H 2 release activity. 42 The experimental assignment of our data is summed up in Figure 6 and Table 1. Figure 7a depicts the progression of amino acid residues involved in catalytic proton transfer as identified in the crystal structure of oxidized [FeFe]-hydrogenase (Hox). 42 This arrangement favors proton uptake and H 2 evolution. Arginine R148 + donates a hydrogen bond to glutamic acid E144 (3.1 Å), the latter forming a hydrogen bond with serine S189 (2.8 Å). 43−45 On the basis of pK a calculations, we previously favored an ionic bond between R148 + and the carboxylate of E144; 42 however, such stabilization fails to explain the FTIR band changes. Trapped between R148 and S189, E144 likely persists in protonated, carboxylic acid form, even at pH values well above the predicted pK a of 3.5 (Table S1).

■ DISCUSSION
Serine S189 is located at the interface of inner and outer core of the proton transfer pathway, as defined above. The distance of ∼3.8 Å between S189 and E141 does reflect discontinued hydrogen bonding, and the probability of proton transfer appears insufficient to justify turnover frequencies greater than 10 000 H 2 s −1 . 9−11 The crystal structure supports

Journal of the American Chemical Society
Article E141 and water molecule W1 forming a trans complex that represents the most stable configuration of COOH groups in aqueous solution. 64−66 This arrangement would interrupt the catalytic hydrogen-bonding network as indicated by the yellow boxes in Figure 7. The E141/W1 complex is weakened by the relatively long donor distance (α = 3.4 Å), which may alleviate changing from W1 to S189 as hydrogen-bond acceptor, for example, upon reduction of the H-cluster. Our IR data are compatible with subtle structural changes upon formation of Hred. Figure 7b illustrates how the catalytic proton transfer pathway between active site cofactor and solvent is significantly more continuous in the reduced enzyme. The downshift of the E141 band may reflect dissolution of the E141/W1 complex and hydrogen bonding to S189 with α > β (in the oxidized crystal, β is 0.4 Å larger than α; see Figure 7a). Serine S189 donates a single bond to E144; however, no second hydrogen bond is formed by E144 due to deprotonation of R148 + . This is in agreement with the pronounced upshift of the E144 band. We assume that a reprotonation of R148 is precluded by the unfavorable local electrostatics betweenE144 and R148 + . Our data support a model in which the E144/R148 interface is defined by the OH donor group of the E144 side chain, held in place by a hydrogen bond from S189 (Figure 7b). Formation of R148 + appears to be hindered in the presence of a potential hydrogenbonding donor. The stability of deprotonated arginine residues in proteins has been questioned 67 and, to our knowledge, has not been observed before. Yet, our spectroscopic investigation on hydrated and deuterated [FeFe]-hydrogenase clearly allows concluding deprotonation of R148 + upon reduction of the enzyme. This conserved arginine functions as proton donor to the active site cofactor, fine-tuning proton transfer efficiency and catalytic bias.
The mechanism of discontinuous proton transfer conceptualized above likely includes a transient step we can speculate about now. Reduction of the H-cluster by one electron leads to an increase in basicity and the formation of Hred upon protonation via the catalytic proton transfer pathway. On the basis of the lack of difference signals in the SH and OH regime, we consider a rigid donor/acceptor conformation between the H-cluster, C169, and W1. Protonation of the ADT headgroup in Hred was discussed 33−35 and computed earlier, 43−45 but the present data do not support steady-state protonation changes at the ADT headgroup, cysteine thiolate, or water cluster. Transiently, however, formation of Hred may trigger deprotonation of E141 ( Figure S8). In a second step, the high basicity of E141 − would induce steady-state deprotonation of R148 + , proton transfer via E144 and S189, and reprotonation of E141, now hydrogen-bonded to S189 instead of W1 (Figure 7b).

■ CONCLUSIONS
In this work, we demonstrate how in situ infrared spectroscopy was applied to analyze the hydrogen-bonding network of the catalytic proton transfer pathway in [FeFe]-hydrogenases. Discontinuous proton transfer was triggered by the enhanced basicity of the active site cofactor (H-cluster) upon photoreduction of a highly active iron−sulfur enzyme lacking a natural chromophore. Infrared spectroscopy provides a direct read-out for changes in hydrogen bonding perfectly complementary to X-ray crystallography. Thereby, the first experimental description of the dynamic hydrogen-bonding changes in the catalytic proton transfer pathway of [FeFe]-hydrogenases was accomplished.
Discontinuous proton transfer is common in nature, for example, in retinal proteins, 68,69 photosystem II, 70 cytochrome c oxidase, 71 and other systems, 48,72 but it has not yet been considered in hydrogenases. Our model rationalizes how ions transcend the gap between inner and outer core of the proton transfer pathway. Furthermore, it provides a reasonable explanation for the catalytic bidirectionality of [FeFe]hydrogenases. 6−8 Glutamic acid E141 switches as hydrogenbonding donor between W1 (proton uptake) and S189 (proton release). Deviations in distance less than 0.5 Å suggest a flexible hydrogen-bonding network that facilitates both H 2 evolution (proton uptake) and H 2 oxidation (proton release).

* S Supporting Information
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10  All distances refer to pdb coordinates 4XDC (HYDA1 numbering). Our data suggest hydrogen-bonding contacts between S189, E144, and R148 + of the catalytic proton transport pathway. Furthermore, invariable hydrogen bonding between ADT, C169, and W1 was observed. The W1/E141 complex establishes the difference between inner and outer core of the catalytic proton transport pathway. (b) In reduced [FeFe]-hydrogenase (Hred), the W1/141 complex is terminated in favor of hydrogen bonding between E141, S189, and E144, which facilitates continuous proton transfer (see main text for details). Our data indicate deprotonation of R148 + upon formation of Hred.

Journal of the American Chemical Society
Article