A detailed QSPR investigation of the water solubility (logS) of 1063 solid neutral chemicals, agrochemicals, drugs and prodrugs has been carried out. The application of the ™General Solubility Equation∫ of Yalkowsky et al. resulted in a correlation between experimental and calculated solubility with rather modest statistic criteria. It was found that 191 compounds (18%) have calculated logS with deviations above one logarithmic unit. A comparison of experimental values with those calculated by an equation previously derived for liquid chemicals demonstrates that a part of the total set of compounds containing certain substructures such as chloroalkyls, phenyls, biphenyls, nitrogen (acyclic and in cycles), benzanthracenes, phenols, and chemicals with ether, ester, and N,N-disubstituted carboxamide groups obey to the equation for liquids. Not unexpectedly, the other part of compounds containing both strong H-bond acceptor and donor groups is essentially less soluble than calculated for compounds in the liquid state. It is proposed that this low solubility is connected with a specific H-bond association of these compounds in their crystal lattice.Two approaches were tested for the derivation of quantitative models for prediction of the solubility of solid chemicals and drugs. The first approach is based on the QSPR for liquids extended by several indicator variables for different functional groups. The second approach is based on the combination of chemical similarity and traditional QSAR techniques. In the framework of this approach different numbers of structural and physicochemical neighbors of a compound-of-interest were considered together with the corresponding HYBOT descriptor data. This method enables to calculate the solubility of a compound-of-interest by using the solubility of nearest neighbor compounds and the difference between descriptor data of those neighbors and the corresponding data of the compound-of-interest. It was found that the use of nearest neighbor compounds with high Tanimoto index value (i.e. good structural similarity) and close H-bond acceptor and donor factors (good physicochemical similarity) ensures good prediction of water solubility with average absolute errors on the level of the error of experimental logS determination.