This dataset includes the computational workflows of a density functional theory based study of the oxygen evolution reaction (OER) on a manganese oxide catalyst in presence and absence of iron dopant. The thermodynamic OER overpotential has been computed by using a surface slab model based on a layered birnessite bulk structure of MnO₂ considering supercells with two and four MnO₂ units and by varying the intercalation with KOH, the amount of Fe dopant and the dopant positions. In addition, the dependence of the oxidation state of the active site atoms (either Mn or Fe) on the directly bound OER intermediate species, has been investigated. The results suggest a decrease of up to 310 mV in the thermodynamic OER overpotential upon doping the considered model structures with Fe that is consistent with the experimentally measured total overpotential decrease of 190 mV.