A combined broken-symmetry density functional and electrostatics approach has been used to model the one-electron reduced and protonated state of the iron-molybdenum cofactor active site of nitrogenase. The active site of the protein contains Fe, Mo, S, N, and O atoms, and many possible sites for protonation have been examined. A novel hydridic proton asymmetrically located in the central cavity created by six Fe sites is most favored from the calculations. Under physiological turnover conditions of low electron flux, the formation of this iron-hydride intermediate may represent a first step towards cofactor liberation of dihydrogen in the absence of dinitrogen.
