We use the meta-MD method developed by Grimme to automatically determine reaction networks for possible catalyst candidates. From the reaction networks we extract the presence/absence of catalytic activity as well as a mechanism for the catalytic cycle. We find that the TS of the rate determining step for the found catalyst candidate can be used for catalyst optimization (lead optimization).