The aim of the present study was to investigate the extent to which the Henderson-Hasselbalch (HH) relationship can be used to predict the pH-dependent aqueous solubility of cationic drugs. The pH-dependent solubility for 25 amines, carrying a single positive charge, was determined with a small-scale shake flask method. Each sample was prepared as a suspension in 150 mM phosphate buffer. The pH-dependent solubility curves were obtained using at least 10 different pH values. The intrinsic solubility, the solubility at the pKa and the solubility at pH values reflecting the pH of the bulk and acid microclimate in the human small intestine (pH 7.4 and 6.5, respectively) were determined for all compounds. The experimental study revealed a large diversity in slope, from -0.5 (celiprolol) to -8.6 (hydralazine) in the linear pH-dependent solubility interval, which is in sharp contrast to the slope of -1 assumed by the HH equation. In addition, a large variation in the range of solubility between the completely uncharged and completely charged drug species was observed. The range for disopyramide was only 1.1 log units, whereas that for amiodarone was greater than 6.3 log units, pointing at the compound specific response to counter-ion effects. In conclusion, the investigated cationic drugs displayed compound specific pH-dependent solubility profiles, indicating that that the HH equation in many cases will only give rough estimations of the pH-dependent solubility of drugs in divalent buffer systems.