Equilibrium structures for oxygen and carbon monoxide complexes of iron porphyrin (FeP), computed with Car-Parrinello molecular dynamics, are in agreement with experimental data of heme models, In addition, we provide information on binding energies and spin-structure relationships, helping to understand the chemistry of the heme prosthetic group. (C) 1997 Published by Elsevier Science B.V.