A method based on statistical mechanical fluctuation solution theory for composition derivatives of activity coefficients is employed for estimating dilute solubilities of 11 solid pharmaceutical solutes in nearly 70 mixed aqueous and nonaqueous solvent systems. The solvent mixtures range from nearly ideal to strongly nonideal. The database covers a temperature range from 293 to 323 K. Comparisons with available data and other existing solubility methods show that the method successfully describes a variety of observed mixed solvent solubility behaviors using solute-solvent parameters from global regression of ternary data as well as predictions based on pure solvent solubilities with an average error of about 10% on mole fractions,