The electronic structures of model oxy-, nitrosyl-, and carbon monoxy ferroporphyrin compounds with N-methylimidazole as an axial ligand are calculated for a number of different conformations of the O2, NO, and CO ligands using the iterative extended Huckel method. For the dioxygen ligand a bent end-on (Pauling) geometry with a low-energy off-axis displacement is favored. The “anomalously large” electric field gradient at the iron nucleus can be accounted for by a paired iron(II)-dioxygen configuration. The carbon monoxide ligand is linear in model compounds but bent or tilted in protein. The reduced affinity observed for the intact protein is explained by these results. The nitrosyl ligand is bent both in model compounds and in the protein. The calculated electronic structure is consistent with the spin density observed in hyperfine splitting in electron spin resonance spectra and gives a quadrupole splitting in excellent agreement with experiment.