We present a method for evaluating second-order Moller-Plesset (MP2) energy and gradients for solvated molecules described within the polarizable continuum model (PCM). The explicit inclusion of solvent effects into the evaluation of the relaxed MP2 density through the Z-vector technique is reported and analyzed. Applications to some one-electron response properties (dipoles, electrostatic molecular potentials, electric field gradients) as well as nuclear gradients are presented.