The most rigorous way is performing molecular dynamics simulation of the substance in a given solvent environment, using such as GROMACS program.

From pure quantum chemistry perspective, you can calculate solvation free energy using SMD implicit solvation model, which is a key factor determining solubility, it is expected that the more negative the value, the higher the solubility. However, you should also consider vaporization free energy of the substance, which reflects the difficulty of scattering the substance from solid state to ideal gas state.