The infrared spectrum of phosphoenol pyruvate (PEP) in aqueous solution was studied experimentally and theoretically in its fully ionized, singly protonated and doubly protonated form. The density functional theory with the B3LYP functional and with the 6-31G(d,p), 6-31++G(d,p), and 6-311++G(d,p) basis sets were used in the theoretical study. The calculations with the two latter basis sets and the CPCM continuum model for water showed good agreement with the experiments except for vibrations assigned to hydroxyl groups. These needed to be modeled with explicit water molecules. The effects of deuteration and of (13)C(2,3) labeling of PEP were reproduced by the calculations.