The three acid dissociation constants (pK a s) of the 8-hydroxyquinoline-5-sulfonic acid were computed using a computational protocol based on Density Functional Theory. A hybrid exchange correlation functional was applied and bulk solvent effects were treated within the framework of the Polarizable Continuum Model. Direct solute–solvent interactions were taken into account adding explicit water molecules. The computed pK a s are in line with the experimental data and allow better defining the first pK a , confirmed to be negative. From the calculated pK a s, ‘ab initio’ distribution diagrams of the relative concentration of the different species in solution as a function of pH were drawn.