In this work, the applicability of this method was extended from atomic liquids (previous work) to water, by introducing a molecule-based definition of the scaling, along with the implementation of additional interactions relevant for rigid solvent molecules, such as reaction-field electrostatics and the handling of bond constraints. Furthermore, a spherical system geometry was defined in addition to the slab geometry, which is more appropriate for solvated globular macromolecules.