The Mo-site and its ligand environment of the FeMo-cofactor (FeMo-co) were studied using the hybrid density functional method B3LYP. The structure and stability of the model complex (S-ligand)3(N-ligand)Mo[(S)-OCH(CH3)C(O)O-] along with its various protonated and reduced/oxidized forms were calculated. Several hypotheses were tested: (i) ligand environment of the Mo-site, (ii) monodentate vs bidentate coordination of the Mo-bound homocitrate ligand, (iii) substrate coordination to the Mo center, and (iv) Mo-His interaction. It was found that the decoordination of one of the homocitrate (lactate in the model) "legs", the bidentate-->monodentate rearrangement, does not occur spontaneously upon either single/double protonation or one-electron reduction. However, it could occur only upon substrate coordination to the Mo-center of the single-protonated forms of the complex. It was shown that one-electron reduction, single-protonation, and substrate coordination facilitate the bidentate<-->monodentate rearrangement of the homocitrate (lactate) ligand of FeMo-co. It was demonstrated that the smallest acceptable model of His ligand in FeMo-co is methylimidazolate (MeIm-). Our studies suggest that the epsilon-N of the FeMo-co-bound His residue is not protonated, and as a consequence the cluster is tightly bound to the protein matrix via a strong Mo-N delta bond.