Stabilities and conformational properties of two Pro --> Thr point mutation models were computed at the B3LYP/6-31G(d) level of theory for the parent triamino acid diamide Pro-Pro-Pro (HCO-Pro-Pro-Pro-NH2). Geometrical parameters for the amino acid sequences, used in the molecular orbital computations for Pro-Pro-Thr and Pro-Thr-Pro, were retrieved from the Protein Data Bank. Thermodynamic functions (S, H, G) were computed for the fully optimized geometries. To assess the stabilization energetics of these mutant models, relative to the parent Pro-Pro-Pro reference conformer epsilon(L) epsilon(L) gamma(L), isodesmic reactions were constructed to calculate DeltaS, DeltaH, and DeltaG. The importance of intramolecular hydrogen bonds involving the -OH group of the Thr side chain, which emerged after the point mutations, was also examined to determine the internal stabilization of these peptide models. This study describes an approach to analyzing a point mutation at the center of a peptide chain and compares its stability to that of a point mutation at a terminal end in a small peptide model.