A density functional theory study of the influence of the various functional groups of the molybdopterin ligand on electronic and geometric properties of active-site models for the molybdenum and tungsten cofactors has been undertaken. We used analogous molybdenum and tungsten complexes with increasingly accurate representation of the molybdopterin ligands and compared bond lengths, angles, charge distribution, composition of the binding orbitals, as well as the redox potentials in relation to each other. On the basis of our findings, we suggest using ligand systems including the pyrane and the pyrazine rings, besides the dithiolene function, to obtain sufficiently reliable computational, but also synthetic, models for the molybdenum and tungsten cofactors, whereas the second ring of the pterin might be neglected for efficiency reasons.Electronic supplementary materialThe online version of this article (doi:10.1007/s00775-009-0548-y) contains supplementary material, which is available to authorized users.