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