Unveiling How Stereoelectronic Factors Affect Kinetics and Thermodynamics of Protonation Regiochemistry in [FeFe] Hydrogenase Synthetic Models: A DFT Investigation

Spectroscopic and Computational Evidence that [FeFe] Hydrogenases Operate Exclusively with CO-Bridged Intermediates

[FeFe] hydrogenases are extremely active H₂-converting enzymes. Their mechanism remains highly controversial, in particular, the nature of the one-electron and two-electron reduced intermediates called H_redH⁺ and H_sredH⁺. In one model, the H_redH⁺ and H_sredH⁺ states contain a semibridging CO, while in the other model, the bridging CO is replaced by a bridging hydride. Using low-temperature IR spectroscopy and nuclear resonance vibrational spectroscopy, together with density functional theory calculations, we show that the bridging CO is retained in the H_sredH⁺ and H_redH⁺ states in the [FeFe] hydrogenases from Chlamydomonas reinhardtii and Desulfovibrio desulfuricans, respectively. Furthermore, there is no evidence for a bridging hydride in either state. These results agree with a model of the catalytic cycle in which the H_redH⁺ and H_sredH⁺ states are integral, catalytically competent components. We conclude that proton-coupled electron transfer between the two subclusters is crucial to catalysis and allows these enzymes to operate in a highly efficient and reversible manner.

Differential Protonation at the Catalytic Six-Iron Cofactor of [FeFe]-Hydrogenases Revealed by 57Fe Nuclear Resonance X-ray Scattering and Quantum Mechanics/Molecular Mechanics Analyses

[FeFe]-hydrogenases are efficient biological hydrogen conversion catalysts and blueprints for technological fuel production. The relations between substrate interactions and electron/proton transfer events at their unique six-iron cofactor (H-cluster) need to be elucidated. The H-cluster comprises a four-iron cluster, [4Fe4S], linked to a diiron complex, [FeFe]. We combined ⁵⁷Fe-specific X-ray nuclear resonance scattering experiments (NFS, nuclear forward scattering; NRVS, nuclear resonance vibrational spectroscopy) with quantum-mechanics/molecular-mechanics computations to study the [FeFe]-hydrogenase HYDA1 from a green alga. Selective ⁵⁷Fe labeling at only [4Fe4S] or [FeFe], or at both subcomplexes was achieved by protein expression with a ⁵⁷Fe salt and in vitro maturation with a synthetic diiron site precursor containing ⁵⁷Fe. H-cluster states were populated under infrared spectroscopy control. NRVS spectral analyses facilitated assignment of the vibrational modes of the cofactor species. This approach revealed the H-cluster structure of the oxidized state (Hox) with a bridging carbon monoxide at [FeFe] and ligand rearrangement in the CO-inhibited state (Hox-CO). Protonation at a cysteine ligand of [4Fe4S] in the oxidized state occurring at low pH (HoxH) was indicated, in contrast to bridging hydride binding at [FeFe] in a one-electron reduced state (Hred). These findings are direct evidence for differential protonation either at the four-iron or diiron subcomplex of the H-cluster. NFS time-traces provided Mössbauer parameters such as the quadrupole splitting energy, which differ among cofactor states, thereby supporting selective protonation at either subcomplex. In combination with data for reduced states showing similar [4Fe4S] protonation as HoxH without (Hred') or with (Hhyd) a terminal hydride at [FeFe], our results imply that coordination geometry dynamics at the diiron site and proton-coupled electron transfer to either the four-iron or the diiron subcomplex discriminate catalytic and regulatory functions of [FeFe]-hydrogenases. We support a reaction cycle avoiding diiron site geometry changes during rapid H₂ turnover.

Models of the iron-only hydrogenase enzyme: structure, electrochemistry and catalytic activity of Fe2(CO)3(μ-dithiolate)(μ,κ1,κ2-triphos)

A series of Fe₂(triphos)(CO)₃(μ-dithiolate) complexes have been prepared and studied as models of the diiron centre in [FeFe]-hydrogenases.

The Molecular Proceedings of Biological Hydrogen Turnover

Over the past two decades, the bioinorganic chemistry of hydrogenases has attracted much interest from basic and applied research. Hydrogenases are highly efficient metalloenzymes that catalyze the reversible reduction of protons to molecular hydrogen (H₂) in all domains of life. Their iron- and nickel-based cofactors represent promising blueprints for the design of biomimetic, synthetic catalysts. In this Account, we address the molecular proceedings of hydrogen turnover with [FeFe]-hydrogenases. The active site cofactor of [FeFe]-hydrogenases ("H-cluster") comprises a unique diiron complex linked to a [4Fe-4S] cluster via a single cysteine. Since it was discovered that a synthetic analogue of the diiron site can be incorporated into apoprotein in vitro to yield active enzyme, significant progress has been made toward a comprehensive understanding of hydrogenase catalysis. The diiron site carries three to four carbon monoxide (CO) and two cyanide (CN^-) ligands that give rise to intense infrared (IR) absorption bands. These bands are sensitive reporters of the electron density across the H-cluster, which can be addressed by infrared spectroscopy to follow redox and protonation changes at the cofactor. Synthetic variation of the metal-bridging dithiolate ligand at the diiron site, as well as site-directed mutagenesis of amino acids, provides access to the proton pathways toward the cofactor. Quantum chemical calculations are employed to specifically assign IR bands to vibrational modes of the diatomic ligands and yield refined H-cluster geometries. Here, we provide an overview of recent research on [FeFe]-hydrogenases with emphasis on experimental and computational IR studies. We describe advances in attenuated total reflection Fourier transform infrared spectroscopy (ATR FTIR) and protein film electrochemistry, as well as density functional theory (DFT) calculations. Key cofactor species are discussed in terms of molecular geometry, redox state, and protonation. Isotope editing is introduced as a tool to evaluate the cofactor geometry beyond the limits of protein crystallography. In particular, the role of proton-coupled electron transfer (PCET) in the generation of catalytically relevant redox species is addressed. We propose that site-selective protonation of the H-cluster biases surplus electrons either to the [4Fe-4S] cluster or to the diiron site. Protonation of the [4Fe-4S] cluster prevents premature reduction at the diiron site and stabilizes a reactive, terminal hydride. The observed H-cluster species are assigned to rapid H₂ conversion or to reactions possibly involved in activity regulation and cellular H₂ sensing. In the catalytic cycle of [FeFe]-hydrogenases, an H-cluster geometry is preserved that features a bridging CO ligand. PCET levels the redox potential for two steps of sequential cofactor reduction. The concept of consecutive PCET at a geometrically inert cofactor with tight control of electron and proton localization may inspire the design of a novel generation of biomimetic catalysts for the production of H₂ as a fuel.

Proton-Coupled Reduction of the Catalytic [4Fe-4S] Cluster in [FeFe]-Hydrogenases

In nature, [FeFe]-hydrogenases catalyze the uptake and release of molecular hydrogen (H₂ ) at a unique iron-sulfur cofactor. The absence of an electrochemical overpotential in the H₂ release reaction makes [FeFe]-hydrogenases a prime example of efficient biocatalysis. However, the molecular details of hydrogen turnover are not yet fully understood. Herein, we characterize the initial one-electron reduction of [FeFe]-hydrogenases by infrared spectroscopy and electrochemistry and present evidence for proton-coupled electron transport during the formation of the reduced state Hred'. Charge compensation stabilizes the excess electron at the [4Fe-4S] cluster and maintains a conservative configuration of the diiron site. The role of Hred' in hydrogen turnover and possible implications on the catalytic mechanism are discussed. We propose that regulation of the electronic properties in the periphery of metal cofactors is key to orchestrating multielectron processes.

Mechanistic Insight into Electrocatalytic H2 Production by [Fe2(CN){μ-CN(Me)2}(μ-CO)(CO)(Cp)2]: Effects of Dithiolate Replacement in [FeFe] Hydrogenase Models

DFT has been used to investigate viable mechanisms of the hydrogen evolution reaction (HER) electrocatalyzed by [Fe₂(CN){μ-CN(Me)₂}(μ-CO)(CO)(Cp)₂] (1) in AcOH. Molecular details underlying the proposed ECEC electrochemical sequence have been studied, and the key functionalities of CN^- and amino-carbyne ligands have been elucidated. After the first reduction, CN^- works as a relay for the first proton from AcOH to the carbyne, with this ligand serving as the main electron acceptor for both reduction steps. After the second reduction, a second protonation occurs at CN^- that forms a Fe(CNH) moiety: i.e., the acidic source for the H₂ generation. The hydride (formally 2e/H⁺), necessary to the heterocoupling with H⁺ is thus provided by the μ-CN(Me)₂ ligand and not by Fe centers, as occurs in typical L₆Fe₂S₂ derivatives modeling the hydrogenase active site. It is remarkable, in this regard, that CN^- plays a role more subtle than that previously expected (increasing electron density at Fe atoms). In addition, the role of AcOH in shuttling protons from CN^- to CN(Me)₂ is highlighted. The incompetence for the HER of the related species [Fe₂{μ-CN(Me)₂}(μ-CO)(CO)₂(Cp)₂]⁺ (2⁺) has been investigated and attributed to the loss of proton responsiveness caused by CN^- replacement with CO. In the context of hydrogenase mimicry, an implication of this study is that the dithiolate strap, normally present in all synthetic models, can be removed from the Fe₂ core without loss of HER, but the redox and acid-base processes underlying turnover switch from a metal-based to a ligand-based chemistry. The versatile nature of the carbyne, once incorporated in the Fe₂ scaffold, could be exploited to develop more active and robust catalysts for the HER.

Interplay between Terminal and Bridging Diiron Hydrides in Neutral and Oxidized States

This study describes the structural, spectroscopic, and electrochemical properties of electronically unsymmetrical diiron hydrides. The terminal hydride Cp*Fe(pdt)Fe(dppe)(CO)H ([1(t-H)]⁰, Cp*^- = Me₅C₅^-, pdt^2- = CH₂(CH₂S^-)₂, dppe = Ph₂PC₂H₄PPh₂) was prepared by hydride reduction of [Cp*Fe(pdt)Fe(dppe)(CO)(NCMe)]⁺. As established by X-ray crystallography, [1(t-H)]⁰ features a terminal hydride ligand. Unlike previous examples of terminal diiron hydrides, [1(t-H)]⁰ does not isomerize to the bridging hydride [1(μ-H)]⁰. Oxidation of [1(t-H)]⁰ gives [1(t-H)]⁺, which was also characterized crystallographically as its BF₄^- salt. Density functional theory (DFT) calculations indicate that [1(t-H)]⁺ is best described as containing an Cp*Fe^III center. In solution, [1(t-H)]⁺ isomerizes to [1(μ-H)]⁺, as anticipated by DFT. Reduction of [1(μ-H)]⁺ by Cp₂Co afforded the diferrous bridging hydride [1(μ-H)]⁰. Electrochemical measurements and DFT calculations indicate that the couples [1(t-H)]^+/0 and [1(μ-H)]^+/0 differ by 210 mV. Qualitative measurements indicate that [1(t-H)]⁰ and [1(μ-H)]⁰ are close in free energy. Protonation of [1(t-H)]⁰ in MeCN solution affords H₂ even with weak acids via hydride transfer. In contrast, protonation of [1(μ-H)]⁰ yields 0.5 equiv of H₂ by a proposed protonation-induced electron transfer process. Isotopic labeling indicates that μ-H/D ligands are inert.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

Mechanisms and origins of the switchable regioselectivity of FeBr3-catalyzed [1,2]-aryl and [1,2]-alkyl shifts of α-aryl aldehydes

With the aid of DFT calculations, the FeBr3-catalyzed skeletal rearrangements of 2-cyclohexanal,2-p-C6H4OMe-propylaldehyde (1A) and 2-phenyl,2-p-C6H4OMe-propylaldehyde (1B) were investigated theoretically. As compared to mono-FeBr3 as a catalyst, the bis-FeBr3 serving as a catalyst is found to be not only enhancing the catalytic efficiency but also improving the product selectivity. For the reaction starting from 1A, the [1,2]-group shift (first step) is rate-determining, and why the Cy shift is the most favored is rationalized in comparison with the p-C6H4OMe and Me shifts. For the reaction starting from 1B, the [1,2]-H shift (second step) is rate-determining although the [1,2]-p-C6H4OMe shift is favored over the [1,2]-phenyl shift. In contrast to the experimental proposal, the newly established H2O/Br(-) joint-assisted H-shift mechanism explains the partial α-H source of the [1,2]-Cy shift product. In addition, we discussed the inherent mechanism that explains why both the [1,2]-p-C6H4OMe and [1,2]-p-C6H4CF3 shifts are more facile than the [1,2]-phenyl shift although the substituents -OMe and -CF3 have opposite electronic behaviors.

Hydrogen generation: aromatic dithiolate-bridged metal carbonyl complexes as hydrogenase catalytic site models

The design, syntheses and characteristics of metal carbonyl complexes with aromatic dithiolate linkers reported as bioinspired hydrogenase catalytic site models are described and reviewed. Among these the complexes capable of hydrogen generation have been discussed in detail. Comparisons have been made with carbonyl complexes having alkyl dithiolates as linkers between metal centers.

New reactions of terminal hydrides on a diiron dithiolate

Mechanisms for biological and bioinspired dihydrogen activation and production often invoke the intermediacy of diiron dithiolato dihydrides. The first example of such a Fe2(SR)2H2 species is provided by the complex [(term-H)(μ-H)Fe2(pdt)(CO)(dppv)2] ([H1H](0)). Spectroscopic and computational studies indicate that [H1H](0) contains both a bridging hydride and a terminal hydride, which, notably, occupies a basal site. The synthesis begins with [(μ-H)Fe2(pdt)(CO)2(dppv)2](+) ([H1(CO)](+)), which undergoes substitution to afford [(μ-H)Fe2(pdt)(CO)(NCMe)(dppv)2](+) ([H1(NCMe)](+)). Upon treatment of [H1(NCMe)](+) with borohydride salts, the MeCN ligand is displaced to afford [H1H](0). DNMR (EXSY, SST) experiments on this complex show that the terminal and bridging hydride ligands interchange intramolecularly at a rate of 1 s(-1) at -40 °C. The compound reacts with D2 to afford [D1D](0), but not mixed isotopomers such as [H1D](0). The dihydride undergoes oxidation with Fc(+) under CO to give [1(CO)](+) and H2. Protonation in MeCN solution gives [H1(NCMe)](+) and H2. Carbonylation converts [H1H](0) into [1(CO)](0).

[FeFe] hydrogenase active site modeling: a key intermediate bearing a thiolate proton and Fe hydride

The first di-protonated [FeFe] hydrogenase model relevant to key intermediates in catalytic hydrogen production is reported. The complex bearing the S-proton and Fe-hydride is structurally and spectroscopically characterized as well as studied by DFT calculations. The results show that the thiolate sulfur can accept protons during the catalytic routes.

Excited state properties of diiron dithiolate hydrides: implications in the unsensitized photocatalysis of H2 evolution

Density functional theory (DFT) and time-dependent DFT (TDDFT) have been used to investigate how visible light photons can excite an asymmetrically substituted diiron hydride, [Fe2(pdt)(μ-H)(CO)4dppv](+) (1(+), dppv = cis-1,2-C2H2(PPh2)2; pdt = 1,3-propanedithiolate), as well as the symmetric species [Fe2(pdt)(μ-H)(CO)4(PMe3)2](+) (2(+)), which are the first photocatalysts of proton reduction operating without employing sensitizers (Wang, W.; Rauchfuss, T. B.; Bertini, L.; Zampella, G.; J. Am. Chem. Soc., 2012, 134, 4525). Theoretical results illustrate that the peculiar reactivity associated to the excited states of 1(+) and 2(+) is compatible with three different scenarios: (i) it can arise from the movement of the hydride ligand from fully bridging to semibridging/terminal coordination, which is expected to be more reactive toward protons; (ii) reactivity could be related to cleavage of a Fe-S bond, which implies formation of a transient Fe penta-coordinate species that would trigger a facile turnstile hydride isomerization, if lifetime excitation is long enough; (iii) also in line with a Fe-S bond cleavage is the possibility that after excited state decay, a highly basic S center is protonated so that a species simultaneously containing S-H(δ+) and Fe-H(δ-) moieties is formed and, once reduced by a suitable electron donor, it can readily afford H2 plus an unprotonated form of the FeFe complex. This last possibility is consistent with (31)P NMR and IR solution data. All the three possibilities are compatible with the capability of 1(+) and 2(+) to perform photocatalysis of hydrogen evolving reaction (HER) without sensitizer. Moreover, even though it turned out difficult to discriminate among the three scenarios, especially because of the lack of experimental excitation lifetimes, it is worth underscoring that all of the three pathways represent a novelty regarding diiron carbonyl photoreactivity, which is usually associated with CO loss. Results provide also a rationale to the experimental observations which showed that the simultaneous presence of donor ligands (dppv in the case of 1(+)) and a H ligand in the coordination environment of diiron complexes is a key factor to prevent CO photodissociation and catalyze HER. Finally, the comparison of photoexcitation behavior of 1(+) and 2(+) allows a sort of generalization about the functioning of such hydride species.

Reactions of [FeFe]-hydrogenase models involving the formation of hydrides related to proton reduction and hydrogen oxidation

[FeFe]-hydrogenases are enzymes in nature that catalyze the reduction of protons and the oxidation of H2 at neutral pH with remarkably high activities and incredibly low overpotential. Structural and functional biomimicking of the active site of [FeFe]-hydrogenases can provide helpful hints for elucidating the mechanism of H2 evolution and uptake at the [FeFe]-hydrogenase active site and for designing bioinspired catalysts to replace the expensive noble metal catalysts for H2 generation and uptake. This perspective focuses on the recent progress in the formation and reactivity of iron hydrides closely related to the processes of proton reduction and hydrogen oxidation mediated by diiron dithiolate complexes. The second section surveys the bridging and terminal hydride species formed from various diiron complexes as well as the intramolecular proton transfer. The very recent progress in H2 activation by diiron dithiolate models are reviewed in the third section. In the concluding remarks and outlook, the differences in structure and catalytic mechanism between the synthetic models and the native [FeFe]-H2ase active site are compared and analyzed, which may cause the need for a significantly larger driving force and may lead to lower activities of synthetic models than the [FeFe]-H2ases for H2 generation and uptake.

Terminal vs bridging hydrides of diiron dithiolates: protonation of Fe2(dithiolate)(CO)2(PMe3)4

This investigation examines the protonation of diiron dithiolates, exploiting the new family of exceptionally electron-rich complexes Fe(2)(xdt)(CO)(2)(PMe(3))(4), where xdt is edt (ethanedithiolate, 1), pdt (propanedithiolate, 2), and adt (2-aza-1,3-propanedithiolate, 3), prepared by the photochemical substitution of the corresponding hexacarbonyls. Compounds 1-3 oxidize near -950 mV vs Fc(+/0). Crystallographic analyses confirm that 1 and 2 adopt C(2)-symmetric structures (Fe-Fe = 2.616 and 2.625 Å, respectively). Low-temperature protonation of 1 afforded exclusively [μ-H1](+), establishing the non-intermediacy of the terminal hydride ([t-H1](+)). At higher temperatures, protonation afforded mainly [t-H1](+). The temperature dependence of the ratio [t-H1](+)/[μ-H1](+) indicates that the barriers for the two protonation pathways differ by ∼4 kcal/mol. Low-temperature (31)P{(1)H} NMR measurements indicate that the protonation of 2 proceeds by an intermediate, proposed to be the S-protonated dithiolate [Fe(2)(Hpdt)(CO)(2)(PMe(3))(4)](+) ([S-H2](+)). This intermediate converts to [t-H2](+) and [μ-H2](+) by first-order and second-order processes, respectively. DFT calculations support transient protonation at sulfur and the proposal that the S-protonated species (e.g., [S-H2](+)) rearranges to the terminal hydride intramolecularly via a low-energy pathway. Protonation of 3 affords exclusively terminal hydrides, regardless of the acid or conditions, to give [t-H3](+), which isomerizes to [t-H3'](+), wherein all PMe(3) ligands are basal.

The reactions of pyridinyl thioesters with triiron dodecacarbonyl: their novel diiron carbonyl complexes and mechanistic investigations

Reaction of Fe(3)(CO)(12) with pyridinyl thioester ligand PyCH(2)SCOCH(3) (L(1), Py = pyridin-2-yl) produced complex, [Fe(2)(κ-COCH(3))(μ-SCH(2)Py)(CO)(5)] (1) (PyCH(2)S = pyridin-2-ylmethanethiolate). When complex 1 reacted with PPh(3), a monosubstituted complex, [Fe(2)(κ-COCH(3))(μ-SCH(2)Py)(CO)(4)PPh(3)] (2), was derived. Reaction of the same precursor with analogous thioester ligand PyCH(2)SCOPy (L(2)) generated three novel diiron complexes, [Fe(2)(κ-Py)(μ-SCH(2)Py)(CO)(5)] (3), [Fe(2)(κ-Py)'(μ-SCH(2)Py)(CO)(5)] (4), and [Fe(2)(κ-Py)(μ-SCH(2)Py)(CO)(6)] (5). Complexes 3 and 4 are structural isomers. Complex 5 could be converted into complex 4 but the conversion from complex 5 to the isomer 3 was not observed. All the five complexes were fully characterised using FTIR, NMR, and other techniques. Their structures were determined using X-ray single crystal diffraction analysis. The oxidative formation of complexes 1, 3, 4, and 5 involved C-S and/or C-C bonds cleavages. To probe possible mechanisms for these cleavages, DFT calculations were performed. From the calculations, viable reaction pathways leading to the formation of all the isolated products were delineated. The results of the theoretic calculations also allowed rationalisation of the experimental observations.

Electron delocalization from the fullerene attachment to the diiron core within the active-site mimics of [FeFe]hydrogenase

Attachment of the redox-active C(60)(H)PPh(2) group modulates the electronic structure of the Fe(2) core in [(μ-bdt)Fe(2)(CO)(5)(C(60)(H)PPh(2))]. The neutral complex is characterized by X-ray crystallography, IR, NMR spectroscopy, and cyclic voltammetry. When it is reduced by one electron, the spectroscopic and density functional theory results indicate that the Fe(2) core is partially spin-populated. In the doubly reduced species, extensive electron communication occurs between the reduced fullerene unit and the Fe(2) centers as displayed in the spin-density plot. The results suggest that the [4Fe4S] cluster within the H cluster provides an essential role in terms of the electronic factor.

On understanding proton transfer to the biocatalytic [Fe-Fe](H) sub-cluster in [Fe-Fe]H(2)ases: QM/MM MD simulations

Proton transfer to the [Fe-Fe](H) sub-cluster in the Desulfovibrio desulfuricans (DdH) and Clostridium pasteurianum (CpI) [Fe-Fe] hydrogenases was investigated by a combination of first principles and empirical molecular dynamics simulations. Pathways that can be inferred from the X-ray crystal structures of DdH and CpI, i.e., (Glu159→Ser198→Glu156→water460→Cys178→DTMA([Fe-Fe](H)) and (Glu282→Ser319→Glu279→water612→Cys299), respectively, were considered. Proton transfer from Cys178 to DTMA in the [Fe-Fe](H) sub-cluster in DdH was readily observed in our results, specifically when [Fe-Fe](H) was in the reduced state ([Fe(I)-Fe(I)]) or in the mixed valence state for the protonated distal iron Fe(d) ([Fe(I)-Fe(II)-H(-)](H)). A concerted mechanism is proposed, where proton transfer in DdH from Glu159 to Glu156 via Ser198 and Glu156 to Cys178 via water460 readily occurred, as well as from Glu282 to Glu279 via Ser319 and Glu279 to Cys299 via water612 in CpI. The theoretical prediction of the proton transfer characteristics is consistent with the assumed biocatalytic mechanism of the [Fe-Fe] hydrogenases in which the proton binds at Fe(d), providing confirmation that has not been explored so far. The computational results were qualitatively validated by the agreement with experimental hydrogen production activity data for mutated CpI enzymes, relative to the wild-type protein. Finally, the insight provided by the simulations, combined, in part, with experimental validation, are important for establishing an approach in future exploration of proton transfer to the active site in this class of enzymes, and possibly also for biomimetic analogs.