To help better understand the effects of small solutes on protein stability, we carry out atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), with a tri-glycine peptide, which is a good model for protein backbone. Multiple solute concentrations are analyzed and each solute-peptide-water ternary system is studied with ~200–300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Γ23) and experimentally measured values suggests that semi-quantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Γ23) result can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Γ23 into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results on GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.