Infrared and Raman spectra of three drug molecules, aspirin, caffeine and ibuprofen, in gas phase and in aqueous solution have been simulated using hybrid density functional theory. The long range solvent effect is modelled by the polarizable continuum model, while the short range hydrogen bonding effects are taken care of by the super-molecular approach with explicit inclusion of water molecules. The calculated spectra are found to compare well with available experimental results. The agreement obtained make grounds for proposing theoretical modeling as a tool for characterizing changes in the bonding environments of drug molecules in terms of particular variations in their IR and Raman spectra.