Theoretical study of the methylprolithospermate's pKa in aqueous solution

Protonation and deprotonation reactions play an important role in the biological processes. Such processes are fully thermodynamic and occur in complex media where functional groups are surrounded by other chemicals such as molecules, ions, water, etc. In this paper, the compound of interest is the methylprolithospermate (MPL) molecule, which was recently isolated from salvia yunnanensis CH Wright. The intended functional groups for deprotonation are -OH and -COOH. These are the potential groups depicted at four sites located in MPL. Three thermodynamic cycles were required for the estimation of MPL's pKa by the means of DFT at the levels of B3LYP and PW6B95 with the dual basis set 6-311++G(d,p)//6-31+G(d). Cycle 1 is the direct pKa calculation, Cycle 2 is also the direct pKa calculation with proton hydration, and Cycle 3, the proton exchange process with hydroxyl anion hydration. In these cycles, water molecules as, reagents as well as products, were held as monomers in one hand (monomer cycle approaches), and cluster on the other hand (cluster cycle approaches). The isodesmic scheme was used as benchmark method for pKa calculation and yielded: pKa1 = 3.6 ± 0.3, pKa2 = 6.6 ± 0.2, pKa3 = 9.3 ± 0.2, and pKa4 = 14.4. The inclusion of binding energies of MPL anions in the pKa calculation has improved the results. The best cycles depend on the level of waters and the level of acidity. The more the number of waters clustered to MPL anions, the more the deviation of pKa values.