Our approach unifies these concepts by automatically determining the CV along which anharmonicity is maximal, guiding exploration along the slowest PES ascent. Instead of explicitly evaluating high-dimensional PES directions, the method first samples the local minimum of the reactant state. By analyzing asymmetries in the local distribution-specifically its skewness-it infers the optimal biasing direction. This eliminates reliance on predefined knowledge, enabling automated reaction pathway discovery and efficient sampling of the most probable reaction channels.