In this study, the temperature-dependent solubility of nicotinamide (niacin) was measured in six neat solvents and five aqueous-organic binary mixtures (methanol, 1,4-dioxane, acetonitrile, DMSO and DMF). It was discovered that the selected set of organic solvents offer all sorts of solvent effects, including co-solvent, synergistic, and anti-solvent features, enabling flexible tuning of niacin solubility. In addition, differential scanning calorimetry was used to characterize the fusion thermodynamics of nicotinamide. In particular, the heat capacity change upon melting was measured. The experimental data were interpreted by means of COSMO-RS-DARE (conductor-like screening model for realistic solvation–dimerization, aggregation, and reaction extension) for concentration dependent reactions. The solute–solute and solute–solvent intermolecular interactions were found to be significant in all of the studied systems, which was proven by the computed mutual affinity of the components at the saturated conditions. The values of the Gibbs free energies of pair formation were derived at an advanced level of theory (MP2), including corrections for electron correlation and zero point vibrational energy (ZPE). In all of the studied systems the self-association of nicotinamide was found to be a predominant intermolecular complex, irrespective of the temperature and composition of the binary system. The application of the COSMO-RS-DARE approach led to a perfect match between the computed and measured solubility data, by optimizing the parameter of intermolecular interactions.