Investigation of ionization equilibria of L-aspartic and L-glutamic acids in water at a microscopic level by DFT method

The pK_a values of different dissociable groups of L-Aspartic and L-Glutamic acids in vacuo and in aqueous medium over a wide pH range have been estimated by DFT/B3LYP/6-31G(d) and 6-31G++(d,p) methods. For both the amino acids discrete water molecules (n?=?0, 3 and 6) have been used to get the first hydration sphere. Starting from a low pH, all possible tautomers resulting from each dissociation step are assumed to exist in a cyclic equilibrium. The structures of the species involved in the dissociation and tautomerization processes have been optimized in vacuo and also in aqueous medium considering H-bonded water molecules under the PCM formalism. For obtaining pK_a in aqueous medium the difference in Gibbs energy of the clusters H₃O⁺.mH₂O and (H₂O) _m+1 (m?=?an integer) is necessary and has been evaluated from computed literature data. Calculations reveal that in vacuo the neutral or less charged species predominate but in aqueous medium the zwitterionic or more chargeable forms contribute appreciably. The Gibbs energy changes for the microscopic cyclic equilibria have been estimated theoretically. These lead to overall (macroscopic) pK_a values for the ionization steps which are in good agreement with available experimental data for both the amino acids.