In this work a density functional theory study of the CO addition reaction to FeIFeI and FeIFeII models of the active site of [FeFe] hydrogenases is presented. A series of model complexes, ranging from simple diiron model complexes of the binuclear [2Fe]H subcluster to the full H-cluster, have been investigated. For each system, the thermodynamic parameters for the CO adduct formation, a reaction that mimics the enzyme CO inhibition, were computed. Parallel to the investigation of the CO addition reaction, the structural features of the various FeIFeI and FeIFeII species have been evaluated, with particular attention to the issue of the ligand arrangement as a function of the redox state. CO affinity depends on the redox state of the model and the chemical nature of its ligands. FeIFeII species are more favorable to form the CO adducts than the reduced FeIFeI species. According to the computed free energies and enthalpies for the CO adduct formation from Fe2(pdt)(CO)5L models, the CO affinity follows the ...