ABSTRACT Accurate binding energies of molecules to water clusters are relevant for understanding intermolecular interactions and various chemical applications. They enter models of interstellar chemical processes, as binding to icy grains influences surface reactions and thus affects calculated gas-phase abundances. Unfortunately, many astrochemical molecules (especially radicals and ions) are incompletely characterised in these models. To address this, we report computational searches for optimal structures and benchmark binding and condensation energies for sets of neutral, radical, cationic, and anionic molecules of astrochemical relevance with clusters of N = 1 − 4 water molecules. These calculations utilised reliable density functionals for geometry optimisation, and coupled cluster (CCSD(T)) single point calculations with large basis sets. Four energetic binding motifs (weak, intermediate, strong or covalently bonded) were observed depending on the chemical nature of the guest molecule. Neutral closed and open-shell molecules with strong dipoles and a greater potential for hydrogen bonding are more tightly bound to water clusters compared to non-polar ones. For closed-shell cationic and anionic species, barrier-less reactions with water clusters occur, which reveals radical-free routes to molecular processing in the gas phase and on amorphous ice surfaces.