On the combined use of discrete solvent models and continuum descriptions of solvent effects in ligand exchange reactions: a case study of the uranyl(VI) aquo ion