How does the push/pull effect of the axial ligand influence the catalytic properties of Compound I of catalase and cytochrome P450?

Density functional theory (DFT) calculations on the chemoselective epoxidation versus hydroxylation reactions of propene by oxoiron porphyrin models mimicking the active sites of catalase, cytochrome P450 (P450) and horseradish peroxidase Compound I (CpdI) are presented. The catalase reactions are concerted and proceed via two-state reactivity patterns on competing doublet and quartet spin state surfaces, but the lowest barrier is the one leading to epoxide products on the doublet spin surface. The results are compared with earlier DFT studies of models of cytochrome P450, horseradish peroxide (HRP), taurine/alpha-ketoglutarate dioxygenase and some synthetic oxoiron catalysts.