LC combined with MS/MS analysis of complex mixtures of protein digests is a reliable and sensitive method for characterization of protein phosphorylation. Peptide retention times (RTs) measured during an LC-MS/MS run depend on both the peptide sequence and the location of modified amino acids. These RTs can be predicted using the LC of biomacromolecules at critical conditions model (BioLCCC). Comparing the observed RTs to those obtained from the BioLCCC model can provide additional validation of MS/MS-based peptide identifications to reduce the false discovery rate and to improve the reliability of phosphoproteome profiling. In this study, energies of interaction between phosphorylated residues and the surface of RP separation media for both "classic" alkyl C18 and polar-embedded C18 stationary phases were experimentally determined and included in the BioLCCC model extended for phosphopeptide analysis. The RTs for phosphorylated peptides and their nonphosphorylated analogs were predicted using the extended BioLCCC model and compared with their experimental RTs. The extended model was evaluated using literary data and a complex phosphoproteome data set distributed through the Association of Biomolecular Resource Facilities Proteome Informatics Research Group 2010 study. The reported results demonstrate the capability of the extended BioLCCC model to predict RTs which may lead to improved sensitivity and reliability of LC-MS/MS-based phosphoproteome profiling.