1,2-dihydroxybenzenes (DHBs) are organic compounds which are widely studied as they are applied to advanced oxidation processes (AOPs). These compounds are also related to the development of oxidative stress, wood biodegradation, and neuronal disease in humans. DHBs are metal ligands with pro-oxidant and antioxidant properties. These activities are related to their chelation properties and a consequence of the deprotonation of their hydroxyl groups. In literature, there are several pKa values for the hydroxyl groups of DHBs. These values vary depending on the experimental conditions or the algorithm used for calculation. In this work, an experimentally validated computational method was implemented in aqueous solution for pKa determination of 24 DHBs. The deprotonation order of the hydroxyl groups in DHB was determined observing a selective deprotonation, which depended on the ability of the substituent to donate or withdraw electron density over the ring.