We have investigated the structure and the vibrational spectrum of peroxynitrite anion in aqueous solution by means of combined quantum-classical (QM/MM) molecular dynamics simulations. In our QM/MM scheme, the reactant was modeled using density functional theory with a Gaussian basis set and the solvent was described using the mean-field TIP4P and the polarizable TIP4P-FQ force fields. The choice of basis sets, functionals and force field parameters has been validated by performing calculations on isolated peroxynitrite and on small peroxynitrite-water complexes. Poor values for isolated peroxynitrite structural properties and vibrational frequencies are found for most ab initio methods, particularly regarding the ON-OO(-) bond distance and the harmonic stretching frequency, probably due to the singlet-triplet instability found in the HF wave function. On the other hand, DFT methods yield results in better agreement with high level CCSD(T) ab initio calculations. We have computed the vibrational spectrum for aqueous peroxynitrite by calculating the Fourier transform of the velocity autocorrelation function extracted from the QM-MM molecular dynamics simulations. Our computational scheme, which allows for the inclusion of both anharmonicity and solvent effects, is able to clarify previous discrepancies regarding the experimental spectra assignments and to shed light on the subtle interplay between solvation and peroxynitrite structure and properties.