We introduce an extension of a recently published method\cite{Mendels2018} to
obtain low-dimensional collective variables for studying multiple states free
energy processes in chemical reactions. The only information needed is a
collection of simple statistics of the equilibrium properties of the reactants
and product states. No information on the reaction mechanism has to be given.
The method allows studying a large variety of chemical reactivity problems
including multiple reaction pathways, isomerization, stereo- and
regiospecificity. We applied the method to two fundamental organic chemical
reactions. First we study the \ce{S_N2} nucleophilic substitution reaction of a
\ce{Cl} in \ce{CH_2 Cl_2} leading to an understanding of the kinetic origin of
the chirality inversion in such processes. Subsequently, we tackle the problem
of regioselectivity in the hydrobromination of propene revealing that the
nature of empirical observations such as the Markovinikov's rules lies in the
chemical kinetics rather than the thermodynamic stability of the products