Due to dopamine's chemical structure and the fact that it has three pKa values, its deprotonation process, in aqueous solution, may involve different chemical species. For instance, the first deprotonation step, from the fully protonated dopamine molecule (H3DA+) to the neutral one (H2DA), will result in zwitterionic species if a proton from one of the OH groups in the catechol ring is lost or into a neutral species if the proton is lost from the amino group. Given that the interaction of such a product with its environment will be quite different depending on its nature, it is very important, therefore, to have an accurate knowledge of which is the dopamine chemical species that results after each deprotonation step. In order to gain a better understanding of dopamine chemistry and to establish a plausible dopamine deprotonation pathway, the optimized geometries of the aforementioned species were calculated in this work by means of the density functionals theory (B3LYP/6-311+G(d,p)) in both cases: in vacuo and with solvent effect, to assess, among other theoretical criteria, the proton affinities of the different dopamine species. This permitted us to propose the following reaction pathway: [reaction in text]. Moreover, the calculations of the chemical shift (NMR-GIAO) modeling the effect of the solvent with a continuum method (PCM) was in agreement with the 13C NMR experimental spectra, which confirmed even further the proposed deprotonation pathway.