As aforementioned, we assumed that the proportions of molecules with low, medium and high solubilities in the test sets were identical to those in the training set. Therefore, following model development we sorted all molecules according to their predicted solubility scores and assigned respective classes using the same percentages of compounds as observed in the training set.