Modeling of the Electrostatic Interaction and Catalytic Activity of [NiFe] Hydrogenases on a Planar Electrode.

Hydrogenases are a group of enzymes that have caught the interest of researchers in renewable energies, due to their ability to catalyze the redox reaction of hydrogen. The exploitation of hydrogenases in electrochemical devices requires their immobilization on the surface of suitable electrodes, such as graphite. The orientation of the enzyme on the electrode is important to ensure a good flux of electrons to the catalytic center, through an array of iron-sulfur clusters. Here we present a computational approach to determine the possible orientations of a [NiFe] hydrogenase (PDB 1e3d) on a planar electrode, as a function of pH, salinity, and electrode potential. The calculations are based on the solution of the linearized Poisson-Boltzmann equation, using the PyGBe software. The results reveal that electrostatic interactions do not truly immobilize the enzyme on the surface of the electrode, but there is instead a dynamic equilibrium between different orientations. Nonetheless, after averaging over all thermally accessible orientations, we find significant differences related to the solution's salinity and pH, while the effect of the electrode potential is relatively weak. We also combine models for the protein adsoption-desorption equilibria and for the electron transfer between the proteins and the electrode to arrive at a prediction of the electrode's activity as a function of the enzyme concentration.

Theoretical insights into [NiFe]-hydrogenases oxidation resulting in a slowly reactivating inactive state.

[NiFe]-hydrogenases catalyse the relevant H2 â†’ 2H+ + 2e- reaction. Aerobic oxidation or anaerobic oxidation of this enzyme yields two inactive states called Ni-A and Ni-B. These states differ for the reactivation kinetics which are slower for Ni-A than Ni-B. While there is a general consensus on the structure of Ni-B, the nature of Ni-A is still controversial. Indeed, several crystallographic structures assigned to the Ni-A state have been proposed, which, however, differ for the nature of the bridging ligand and for the presence of modified cysteine residues. The spectroscopic characterization of Ni-A has been of little help due to small differences of calculated spectroscopic parameters, which does not allow to discriminate among the various forms proposed for Ni-A. Here, we report a DFT investigation on the nature of the Ni-A state, based on systematic explorations of conformational and configurational space relying on accurate energy calculations, and on comparisons of theoretical geometries with the X-ray structures currently available. The results presented in this work show that, among all plausible isomers featuring various protonation patterns and oxygenic ligands, the one corresponding to the crystallographic structure recently reported by Volbeda et al. (J Biol Inorg Chem 20:11-22, 19)-featuring a bridging hydroxide ligand and the sulphur atom of Cys64 oxidized to bridging sulfenate-is the most stable. However, isomers with cysteine residues oxidized to terminal sulfenate are very close in energy, and modifications in the network of H-bond with neighbouring residues may alter the stability order of such species.