Study of the stepwise deprotonation reactions of glyphosate and the corresponding pKa values in aqueous solution

Glyphosate (N-(phosphonomethyl)glycine) (Gph) is a herbicide that is broadly used in several countries. Its application to eliminate weeds may have the undesired effect of diminishing the metallic cations found in the soil (e.g., Ni(2+) and Zn(2+)), due to a complexation reaction that depends on the soil's pH. To better understand the molecular structures of glyphosate that are involved in such a complexation reaction, we have studied all possible glyphosate conformations in aqueous solution that may be involved in deprotonation reactions in the pH range from 2 to 11 using the polarizable continuum method (PCM). We have also compared direct (or absolute) methods to calculate pKa values, the cluster-continuum model and the proton-exchange scheme, using different thermodynamic cycles. The best result was achieved when using a proton-exchange scheme, which was able to properly reproduce three glyphosate experimental pKa values predicted for the glyphosate structures and conformations previously determined.