XUV photoelectron spectroscopy (XPS) is a powerful method for investigating the electronic structures of molecules. However, the correct interpretation of results in the condensed phase requires theoretical models that account for solvation. Here we present experimental aqueous-phase XPS of two organic biomimetic molecular switches, NAIP and p-HDIOP. These switches are structurally similar, but have opposite charges and thus present a stringent benchmark for solvation models which need to reproduce the observed ΔeBE = 1.1 eV difference in electron binding energy compared to the 8 eV difference predicted in the gas phase. We present calculations using implicit and explicit solvent models. The latter employs the average solvent electrostatic configuration and free energy gradient (ASEC-FEG) approach. Both nonequilibrium polarizable continuum models and ASEC-FEG calculations give vertical binding energies in good agreement with the experiment for three different computational protocols. Counterions, explicitly accounted for in ASEC-FEG, contribute to the stabilization of molecular states and reduction of ΔeBE upon solvation.