Computational modeling of catalytic reactivity