We present a computational methodology based on a polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM) approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems. This approach is applied for the calculation of Infrared (IR), Vibrational Circular Dichroism (VCD), Raman and Raman Optical Activity (ROA) spectra of aqueous solutions of (l)-methyl lactate and (S)-glycidol. Remarkable agreement between calculations and experiments is reported, showing the reliability and accuracy of the methodology, especially with respect to standard continuum solvation approaches.