A computational approach has been developed to automatically generate and analyse the structures of the intermediates of palladium-catalysed carbon-hydrogen (C-H) activation reactions as well as to predict the final products. Implemented as a high-performance computing cluster tool, it has been shown to correctly choose the mechanism and rationalise regioselectivity of chosen examples from open literature reports. The developed methodology is capable of predicting reactivity of various substrates by differentiation between two major mechanisms - proton abstraction and electrophilic aromatic substitution. An attempt has been made to predict new C-H activation reactions. This methodology can also be used for the automated reaction planning, as well as a starting point for microkinetic modelling.