Effective simulation of biological systems: Choice of density functional and basis set for heme-containing complexes