Two sites of a Pro-Pro diamide were subjected to individual Pro -> Thr point mutations. The parent diamide Pro-Pro as well as selected conformers of the Pro-Thr and Thr-Pro mutant models were subjected to molecular computations at the B3LYP/6-31 G(d) level of theory. At the optimized geometries, thermodynamic functions IN (S, H, and G) were computed. In order to assess relative stabilities of the mutant models, isodesmic reactions were constructed to calculate Delta S, Delta H, and Delta G, relative to the initial Pro-Pro state. The importance of intramolecular hydrogen bonds, involving the -OH group of the Thr side chain, which emerged after the point mutations were also examined. Our findings suggest a novel approach to analyzing the stability of point mutants in peptide models through the analysis of thermodynamic functions.