Quasiparticle (QP) wave functions, also known as Dyson orbitals, extend the concept of single-particle states to interacting electron systems. Here we employ many-body perturbation theory in the GW approximation to calculate the QP wave functions for a semiempirical model describing a pi-conjugated molecular wire in contact with a metal surface. We find that image charge effects pull the frontier molecular orbitals toward the metal surface, while orbitals with higher or lower energy are pushed away. This affects both the size of the energetic image charge shifts and the coupling of the individual orbitals to the metal substrate. Full diagonalization of the QP equation and, to some extent, self-consistency in the GW self-energy, is important to describe the effect, which is not captured by standard density functional theory or Hartree-Fock. These results should be important for the understanding and theoretical modeling of electron transport across metal-molecule interfaces.