The MACE model's applicability to open-shell structures allowed us to predict bond dissociation energies, of interest because of their relevance in cytochrome P450 metabolism prediction. We have shown that MACE matches the reference DFT values within RMSE of 1.16 -1.59 kcal/mol, with the most appropriate comparison leading to a RMSE of 1.37 kcal/mol. In addition, the order of BDEs within a molecule is predicted better with MACE than it is with a semi-empirical AM1 method, commonly used for high-throughput BDE and site of metabolism predictions, a more recent GFN2-xTB semi-empirical method or an alternative ALFABET model for predicting the bond dissociation enthalpies. Overall, this study showcases the potential of MLIPs in modelling structures beyond near-equilibrium or closed-shell compounds. Going towards a method to replace semi-empirical methods, such as commonly used AM1, in metabolism modelling, the current work would need to be extended to include more chemical elements than carbon, hydrogen and oxygen.