Proton Transport Pathways in [NiFe]-Hydrogenase

A Chimeric NiFe Hydrogenase Heterodimer to Assess the Role of the Electron Transfer Chain in Tuning the Enzyme's Catalytic Bias and Oxygen Tolerance

The observation that some homologous enzymes have the same active site but very different catalytic properties demonstrates the importance of long-range effects in enzyme catalysis, but these effects are often difficult to rationalize. The NiFe hydrogenases 1 and 2 (Hyd 1 and Hyd 2) from E. coli both consist of a large catalytic subunit that embeds the same dinuclear active site and a small electron-transfer subunit with a chain of three FeS clusters. Hyd 1 is mostly active in H2 oxidation and resistant to inhibitors, whereas Hyd 2 also catalyzes H2 production and is strongly inhibited by O2 and CO. Based on structural and site-directed mutagenesis data, it is currently believed that the catalytic bias and tolerance to O2 of Hyd 1 are defined by the distal and proximal FeS clusters, respectively. To test these hypotheses, we produced and characterized a hybrid enzyme made of the catalytic subunit of Hyd 1 and the electron transfer subunit of Hyd 2. We conclude that catalytic bias and sensitivity to CO are set by the catalytic subunit rather than by the electron transfer chain. We confirm the importance of the proximal cluster in making the enzyme Hyd 1 resist long-term exposure to O2, but we show that other structural determinants, in both subunits, contribute to O2 tolerance. A similar strategy based on the design of chimeric heterodimers could be used in the future to elucidate various structure-function relationships in hydrogenases and other multimeric metalloenzymes and to engineer useful hydrogenases that combine the desirable properties of distinct, homologous enzymes.

Structure of the membrane-bound formate hydrogenlyase complex from Escherichia coli

The prototypical hydrogen-producing enzyme, the membrane-bound formate hydrogenlyase (FHL) complex from Escherichia coli, links formate oxidation at a molybdopterin-containing formate dehydrogenase to proton reduction at a [NiFe] hydrogenase. It is of intense interest due to its ability to efficiently produce H₂ during fermentation, its reversibility, allowing H₂-dependent CO₂ reduction, and its evolutionary link to respiratory complex I. FHL has been studied for over a century, but its atomic structure remains unknown. Here we report cryo-EM structures of FHL in its aerobically and anaerobically isolated forms at resolutions reaching 2.6 Å. This includes well-resolved density for conserved loops linking the soluble and membrane arms believed to be essential in coupling enzymatic turnover to ion translocation across the membrane in the complex I superfamily. We evaluate possible structural determinants of the bias toward hydrogen production over its oxidation and describe an unpredicted metal-binding site near the interface of FdhF and HycF subunits that may play a role in redox-dependent regulation of FdhF interaction with the complex.

New cryo-EM structures of the formate hydrogenlyase complex from the model bacterium E. coli clarify how electrons and protons move through the complex and are combined to make H2 gas. The complex shows important similarities and differences to related bioenergetic complexes across the tree of life.

Functional Dynamics of an Ancient Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) is a redox-driven Na⁺/H⁺ transporter that employs the energy from hydrogen gas (H₂) production to catalyze proton pumping and Na⁺/H⁺ exchange across cytoplasmic membranes of archaea. Despite a recently resolved structure of this ancient energy-transducing enzyme [Yu et al. Cell2018, 173, 1636–1649], the molecular principles of its redox-driven ion-transport mechanism remain puzzling and of major interest for understanding bioenergetic principles of early cells. Here we use atomistic molecular dynamics (MD) simulations in combination with data clustering methods and quantum chemical calculations to probe principles underlying proton reduction as well as proton and sodium transport in Mbh from the hyperthermophilic archaeon Pyrococcus furiosus. We identify putative Na⁺ binding sites and proton pathways leading across the membrane and to the NiFe-active center as well as conformational changes that regulate ion uptake. We suggest that Na⁺ binding and protonation changes at a putative ion-binding site couple to proton transfer across the antiporter-like MbhH subunit by modulating the conformational state of a conserved ion pair at the subunit interface. Our findings illustrate conserved coupling principles within the complex I superfamily and provide functional insight into archaeal energy transduction mechanisms.

Unifying Activity, Structure, and Spectroscopy of [NiFe] Hydrogenases: Combining Techniques To Clarify Mechanistic Understanding

Achieving a unified understanding of the mechanism of a multicenter redox enzyme such as [NiFe] hydrogenase is complicated by difficulties in reconciling information obtained by using different techniques and on samples in different physical forms. Measurements of the activity of the enzyme, and of factors which perturb activity, are generally carried out using biochemical assays in solution or with electrode-immobilized enzymes using protein film electrochemistry (PFE). Conversely, spectroscopy aimed at reporting on features of the metalloclusters in the enzyme, such as electron paramagnetic resonance (EPR) or X-ray absorption spectroscopy (XAS), is often conducted on frozen samples and is thus difficult to relate to catalytically relevant states as information about turnover and activity has been lost. To complicate matters further, most of our knowledge of the atomic-level structure of metalloenzymes comes from X-ray diffraction studies in the solid, crystalline state, which are again difficult to link to turnover conditions. Taking [NiFe] hydrogenases as our case study, we show here how it is possible to apply infrared (IR) spectroscopic sampling approaches to unite direct spectroscopic study with catalytic turnover. Using a method we have named protein film IR electrochemistry (PFIRE), we reveal the steady-state distribution of intermediates during catalysis and identify catalytic "bottlenecks" introduced by site-directed mutagenesis. We also show that it is possible to study dynamic transitions between active site states of enzymes in single crystals, uniting solid state and solution spectroscopic information. In all of these cases, the spectroscopic data complement and enhance interpretation of purely activity-based measurements by providing direct chemical insight that is otherwise hidden. The [NiFe] hydrogenases possess a bimetallic [NiFe] active site, coordinated by CO and CN^- ligands, linked to the protein via bridging and terminal cysteine sulfur ligands, as well as an electron relay chain of iron sulfur clusters. Infrared spectroscopy is ideal for probing hydrogenases because the CO and CN^- ligands are strong IR absorbers, but the suite of IR-based approaches we describe here will be equally valuable in studying substrate- or intermediate-bound states of other metalloenzymes where key mechanistic questions remain open, such as nitrogenase, formate dehydrogenase, or carbon monoxide dehydrogenase. We therefore hope that this Account will encourage future studies which unify information from different techniques across bioinorganic chemistry.

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

[NiFe] hydrogenases are enzymes that catalyze the splitting of molecular hydrogen according to the reaction H₂ → 2H⁺ + 2e^-. Most of these enzymes are inhibited even by low traces of O₂. However, a special group of O₂-tolerant hydrogenases exists. A member of this group is the membrane-bound [NiFe] hydrogenase from Ralstonia eutropha ( ReMBH). The ReMBH harbors an unusual iron sulfur cluster with composition 4Fe3S(6Cys) that is able to undergo structural changes triggering the flow of two electrons to the [NiFe] active site. These electrons promote oxygen reduction at the active site, preventing, in this way, aerobic inactivation of the enzyme. In the superoxidized state, the [4Fe3S] cluster binds to a hydroxyl group that originates from either molecular oxygen or water reaching the site. Both reactions, oxygen reduction to water at the [NiFe]- or [4Fe3S]-centers and oxygen evolution from water at the proximal cluster, require the delivery of protons regulated by a subtle communication mechanism between these metal centers. In this work, we sequentially apply multiscale modeling techniques as quantum mechanical/molecular mechanics methods and classical molecular dynamics simulations to investigate the role of two distinct proton transfer pathways connecting the [NiFe] active site and the [4Fe3S] proximal cluster of ReMBH in the protection mechanism against an oxygen attack. Although the "glutamate" pathway is preferred by protons migrating toward the active site to avoid inactivation by O₂, the "histidine" pathway plays an essential role in delivering protons for O₂ reduction at the proximal cluster. The results obtained in this work not only provide new pieces to the puzzling catalytic mechanisms governing O₂-tolerant hydrogenases but also highlight the relevance of dynamics in the proper description of biochemical reactions in general.

Computational analyses, molecular dynamics, and mutagenesis studies of unprocessed form of [NiFe] hydrogenase reveal the role of disorder for efficient enzyme maturation

Biological production and oxidation of hydrogen is mediated by hydrogenases, key enzymes for these energy-relevant reactions. Synthesis of [NiFe] hydrogenases involves a complex series of biochemical reactions to assemble protein subunits and metallic cofactors required for enzyme function. A final step in this biosynthetic pathway is the processing of a C-terminal tail (CTT) from its large subunit, thus allowing proper insertion of nickel in the unique NiFe(CN)₂CO cofactor present in these enzymes. In silico modelling and Molecular Dynamics (MD) analyses of processed vs. unprocessed forms of Rhizobium leguminosarum bv. viciae (Rlv) hydrogenase large subunit HupL showed that its CTT (residues 582-596) is an intrinsically disordered region (IDR) that likely provides the required flexibility to the protein for the final steps of proteolytic maturation. Prediction of pKa values of ionizable side chains in both forms of the enzyme's large subunit also revealed that the presence of the CTT strongly modify the protonation state of some key residues around the active site. Furthermore, MD simulations and mutant analyses revealed that two glutamate residues (E27 in the N-terminal region and E589 inside the CTT) likely contribute to the process of nickel incorporation into the enzyme. Computational analysis also revealed structural details on the interaction of Rlv hydrogenase LSU with the endoprotease HupD responsible for the removal of CTT.

Diffusion network of CO in FeFe-Hydrogenase

FeFe-hydrogenase is an efficient enzyme to produce H₂ under optimal conditions. However, the activity of this enzyme is highly sensitive to the presence of inhibitory gases CO and O₂ that cause irreversible damage to the active site. Therefore, a detailed knowledge of the diffusion pathways of these inhibitory gases is necessary to develop strategies for designing novel enzymes that are tolerant to these gases. In this work, we studied the diffusion pathways of CO in the CpI FeFe-hydrogenase from Clostridium pasteurianum. Specifically, we used several enhanced sampling and free-energy simulation methods to reconstruct a three-dimensional free-energy surface for CO diffusion which revealed 45 free-energy minima forming an interconnected network of pathways. We discovered multiple pathways of minimal free-energy as diffusion portals for CO and found that previously suggested hydrophobic pathways are not thermodynamically favorable for CO diffusion. We also observed that the global minimum in the free-energy surface is located in the vicinity of the active-site metal cluster, the H-cluster, which suggests a high-affinity for CO near the active site. Among 19 potential residues that we propose as candidates for future mutagenesis studies, 11 residues are shared with residues that have been previously proposed to increase the tolerance of this enzyme for O₂. We hypothesize that these shared candidate residues are potentially useful for designing new variants of this enzyme that are tolerant to both inhibitory gases.

Hydrogen bonding effect between active site and protein environment on catalysis performance in H2-producing [NiFe] hydrogenases

The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol^-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H⁺/e^- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S^- in cysteine residues shows a larger influence on the energy profile than that via the CN^- ligands on Fe.

The oxygen reduction reaction on [NiFe] hydrogenases

Oxygen tolerance capacity is critical for hydrogen oxidation/evolution catalysts.

Quantum chemical approaches to [NiFe] hydrogenase

The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.

Proton Transfer in the Catalytic Cycle of [NiFe] Hydrogenases: Insight from Vibrational Spectroscopy

Catalysis of H₂ production and oxidation reactions is critical in renewable energy systems based around H₂ as a clean fuel, but the present reliance on platinum-based catalysts is not sustainable. In nature, H₂ is oxidized at minimal overpotential and high turnover frequencies at [NiFe] catalytic sites in hydrogenase enzymes. Although an outline mechanism has been established for the [NiFe] hydrogenases involving heterolytic cleavage of H₂ followed by a first and then second transfer of a proton and electron away from the active site, details remain vague concerning how the proton transfers are facilitated by the protein environment close to the active site. Furthermore, although [NiFe] hydrogenases from different organisms or cellular environments share a common active site, they exhibit a broad range of catalytic characteristics indicating the importance of subtle changes in the surrounding protein in controlling their behavior. Here we review recent time-resolved infrared (IR) spectroscopic studies and IR spectroelectrochemical studies carried out in situ during electrocatalytic turnover. Additionally, we re-evaluate the significant body of IR spectroscopic data on hydrogenase active site states determined through more conventional solution studies, in order to highlight mechanistic steps that seem to apply generally across the [NiFe] hydrogenases, as well as steps which so far seem limited to specific groups of these enzymes. This analysis is intended to help focus attention on the key open questions where further work is needed to assess important aspects of proton and electron transfer in the mechanism of [NiFe] hydrogenases.

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.

Structure and function of [NiFe] hydrogenases

Hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons via a heterolytic splitting mechanism. The active sites of [NiFe] hydrogenases comprise a dinuclear Ni-Fe center carrying CO and CN^- ligands. The catalytic activity of the standard (O₂-sensitive) [NiFe] hydrogenases vanishes under aerobic conditions. The O₂-tolerant [NiFe] hydrogenases can sustain H₂ oxidation activity under atmospheric conditions. These hydrogenases have very similar active site structures that change the ligand sphere during the activation/catalytic process. An important structural difference between these hydrogenases has been found for the proximal iron-sulphur cluster located in the vicinity of the active site. This unprecedented [4Fe-3S]-6Cys cluster can supply two electrons, which lead to rapid recovery of the O₂ inactivation, to the [NiFe] active site.

Isotopic fractionation associated with [NiFe]- and [FeFe]-hydrogenases

RATIONALE

Hydrogenases catalyze the reversible formation of H2 from electrons and protons with high efficiency. Understanding the relationships between H2 production, H2 uptake, and H2-H2O exchange can provide insight into the metabolism of microbial communities in which H2 is an essential component in energy cycling.

METHODS

We used stable H isotopes (1H and 2H) to probe the isotope effects associated with three [FeFe]-hydrogenases and three [NiFe]-hydrogenases.

RESULTS

All six hydrogenases displayed fractionation factors for H2 formation that were significantly less than 1, producing H2 that was severely depleted in 2H relative to the substrate, water. Consistent with differences in their active site structure, the fractionation factors for each class appear to cluster, with the three [NiFe]-hydrogenases (α = 0.27–0.40) generally having smaller values than the three [FeFe]-hydrogenases (α = 0.41–0.55). We also obtained isotopic fractionation factors associated with H2 uptake and H2-H2O exchange under conditions similar to those utilized for H2 production, providing a more complete picture of the reactions catalyzed by hydrogenases.

CONCLUSIONS

The fractionation factors determined in our studies can be used as signatures for different hydrogenases to probe their activity under different growth conditions and to ascertain which hydrogenases are most responsible for H2 production and/or uptake in complex microbial communities.

Implementation of photobiological H2 production: the O 2 sensitivity of hydrogenases

The search for the ultimate carbon-free fuel has intensified in recent years, with a major focus on photoproduction of H2. Biological sources of H2 include oxygenic photosynthetic green algae and cyanobacteria, both of which contain hydrogenase enzymes. Although algal and cyanobacterial hydrogenases perform the same enzymatic reaction through metallo-clusters, their hydrogenases have evolved separately, are expressed differently (transcription of algal hydrogenases is anaerobically induced, while bacterial hydrogenases are constitutively expressed), and display different sensitivity to O2 inactivation. Among various physiological factors, the sensitivity of hydrogenases to O2 has been one of the major factors preventing implementation of biological systems for commercial production of renewable H2. This review addresses recent strategies aimed at engineering increased O2 tolerance into hydrogenases (as of now mainly unsuccessful), as well as towards the development of methods to bypass the O2 sensitivity of hydrogenases (successful but still yielding low solar conversion efficiencies). The author concludes with a description of current approaches from various laboratories to incorporate multiple genetic traits into either algae or cyanobacteria to jointly address limiting factors other than the hydrogenase O2 sensitivity and achieve more sustained H2 photoproduction activity.

Hydrated Excess Protons Can Create Their Own Water Wires

Grotthuss shuttling of an excess proton charge defect through hydrogen bonded water networks has long been the focus of theoretical and experimental studies. In this work we show that there is a related process in which water molecules move ("shuttle") through a hydrated excess proton charge defect in order to wet the path ahead for subsequent proton charge migration. This process is illustrated through reactive molecular dynamics simulations of proton transport through a hydrophobic nanotube, which penetrates through a hydrophobic region. Surprisingly, before the proton enters the nanotube, it starts "shooting" water molecules into the otherwise dry space via Grotthuss shuttling, effectively creating its own water wire where none existed before. As the proton enters the nanotube (by 2-3 Å), it completes the solvation process, transitioning the nanotube to the fully wet state. By contrast, other monatomic cations (e.g., K(+)) have just the opposite effect, by blocking the wetting process and making the nanotube even drier. As the dry nanotube gradually becomes wet when the proton charge defect enters it, the free energy barrier of proton permeation through the tube via Grotthuss shuttling drops significantly. This finding suggests that an important wetting mechanism may influence proton translocation in biological systems, i.e., one in which protons "create" their own water structures (water "wires") in hydrophobic spaces (e.g., protein pores) before migrating through them. An existing water wire, e.g., one seen in an X-ray crystal structure or MD simulations without an explicit excess proton, is therefore not a requirement for protons to transport through hydrophobic spaces.

A threonine stabilizes the NiC and NiR catalytic intermediates of [NiFe]-hydrogenase

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

Amino acid modified Ni catalyst exhibits reversible H2 oxidation/production over a broad pH range at elevated temperatures

Hydrogenases interconvert H2 and protons at high rates and with high energy efficiencies, providing inspiration for the development of molecular catalysts. Studies designed to determine how the protein scaffold can influence a catalytically active site have led to the synthesis of amino acid derivatives of [Ni(P2(R)N2(R'))2](2+) complexes, [Ni(P2(Cy)N2(Amino acid))2](2+) (CyAA). It is shown that these CyAA derivatives can catalyze fully reversible H2 production/oxidation at rates approaching those of hydrogenase enzymes. The reversibility is achieved in acidic aqueous solutions (pH = 0-6), 1 atm 25% H2/Ar, and elevated temperatures (tested from 298 to 348 K) for the glycine (CyGly), arginine (CyArg), and arginine methyl ester (CyArgOMe) derivatives. As expected for a reversible process, the catalytic activity is dependent upon H2 and proton concentrations. CyArg is significantly faster in both directions (∼300 s(-1) H2 production and 20 s(-1) H2 oxidation; pH = 1, 348 K, 1 atm 25% H2/Ar) than the other two derivatives. The slower turnover frequencies for CyArgOMe (35 s(-1) production and 7 s(-1) oxidation under the same conditions) compared with CyArg suggests an important role for the COOH group during catalysis. That CyArg is faster than CyGly (3 s(-1) production and 4 s(-1) oxidation) suggests that the additional structural features imparted by the guanidinium groups facilitate fast and reversible H2 addition/release. These observations demonstrate that outer coordination sphere amino acids work in synergy with the active site and can play an important role for synthetic molecular electrocatalysts, as has been observed for the protein scaffold of redox active enzymes.

Atomic partitioning of M-H2 bonds in [NiFe] hydrogenase--a test case of concurrent binding

The possibility of simultaneous addition of η(2)-H2 to both the metals (Ni and Fe) in the active site of the as isolated state of the enzyme (Ni-SI) is examined here by an atom-by-atom electronic energy partitioning based on the QTAIM method. Results show that the 4LS state prefers H2 removal than addition. Destabilization of the atomic basins of the thiolate bridges and decrease of the electrophilicity of the Fe and Ni, resulting in poor back donation to the CO ligand, are the bottlenecks that hamper dihydrogen activation simultaneously. The study helps to understand why such states are seldom accessed in the activation of dihydrogen. Moreover, Ni has been found to be the natural choice for the dihydrogen binding.

Molecular dynamics study of the proposed proton transport pathways in [FeFe]-hydrogenase

Possible proton transport pathways in Clostridium pasteurianum (CpI) [FeFe]-hydrogenase were investigated with molecular dynamics simulations. This study was undertaken to evaluate the functional pathway and provide insight into the hydrogen bonding features defining an active proton transport pathway. Three pathways were evaluated, two of which consist of water wires and one of predominantly amino acid residues. Our simulations suggest that protons are not transported through water wires. Instead, the five-residue motif (Glu282, Ser319, Glu279, H2O, Cys299) was found to be the likely pathway, consistent with previously made experimental observations. The pathway was found to have a persistent hydrogen bonded core (residues Cys299 to Ser319), with less persistent hydrogen bonds at the ends of the pathway for both H2 release and H2 uptake. Single site mutations of the four residues have been shown experimentally to deactivate the enzyme. The theoretical evaluation of these mutations demonstrates redistribution of the hydrogen bonds in the pathway, resulting in enzyme deactivation. Finally, coupling between the protein dynamics near the proton transport pathway and the redox partner binding regions was also found as a function of H2 uptake and H2 release states, which may be indicative of a correlation between proton and electron movement within the enzyme.

Influence of the protein structure surrounding the active site on the catalytic activity of [NiFeSe] hydrogenases

A combined experimental and theoretical study of the catalytic activity of a [NiFeSe] hydrogenase has been performed by H/D exchange mass spectrometry and molecular dynamics simulations. Hydrogenases are enzymes that catalyze the heterolytic cleavage or production of H2. The [NiFeSe] hydrogenases belong to a subgroup of the [NiFe] enzymes in which a selenocysteine is a ligand of the nickel atom in the active site instead of cysteine. The aim of this research is to determine how much the specific catalytic properties of this hydrogenase are influenced by the replacement of a sulfur by selenium in the coordination of the bimetallic active site versus the changes in the protein structure surrounding the active site. The pH dependence of the D2/H(+) exchange activity and the high isotope effect observed in the Michaelis constant for the dihydrogen substrate and in the single exchange/double exchange ratio suggest that a "cage effect" due to the protein structure surrounding the active site is modulating the enzymatic catalysis. This "cage effect" is supported by molecular dynamics simulations of the diffusion of H2 and D2 from the outside to the inside of the protein, which show different accumulation of these substrates in a cavity next to the active site.