Several studies using IR-UV double-resonance spectroscopy have reported νOH wavenumbers for phenol in small phenol-water clusters, showing red shifts of ∼130, ∼270 and ∼310 cm -1 for clusters of 1, 2 and 3 water molecules, respectively. Based on these experimental results, previous theoretical studies on phenol vibrational spectra have concluded that the νOH wavenumber is mainly dependent on the curvature of the hydrogen-bonding interaction. This experimental trend highlights the sensitivity of the νOH wavenumber to hydrogen-bond interactions, and we explore this feature with our ML/MM implementation to assess how well it captures the impact of solvation on the potential energy surface and vibrational modes.