A transferable MLIP trained on neutral molecules and sp 3 carbon radicals, as explained above, allows us to calculate the homolytic bond dissociation energies (BDEs) and radical reaction energies. Accurate BDE estimation is of interest in the context of CYP metabolism prediction, where the relative reactivity of different sites is crucial in determining the site(s)