From section -The Quantum Chemical Workflow, we can determine the empirical values a and b in Equation ??. For each set of deprotonated sites in a molecule, we extract computed ∆G • min and fit it against the experimental pK a values. Hereafter, we convert the computed ∆G • min to computed QM pK a values using Equation . We then inspect outliers that exceed an absolute pK a unit difference of 5 pK a units between the experimental pK a value and the QM computed pK a value. We choose an absolute pK a unit difference of 5 pK a units to ensure the QM computed pK a is well above the error that is to be expected on the level of theory we are using (CAM-B3LYP). The observed outliers typically result from one of the following reasons: calculation errors concerning the expected minimum pK a site, discrepancies between literature structures and database structures, mislabeled experimental pK a values, or extrapolated pK a values. Notably, the extrapolated pK a values correspond to compounds beyond the scale (pK a ≥ 35) measurable in DMSO due to DMSO's autoprotolysis (pK a =35). For more information regarding finding and removing outliers, see SI -Finding outliers section. After multiple iterations, we identified 695 molecules to have reliable experimental pK a values and computed ∆G