Dysregulated protein phosphorylation is a primary culprit in multiple physiopathological states. Hence, although analysis of signaling cascades on a proteome-wide scale would provide significant insight into both normal and aberrant cellular function, such studies are simultaneously limited by sheer biological complexity and concentration dynamic range. In principle, immobilized metal affinity chromatography (IMAC) represents an ideal enrichment method for phosphoproteomics. However, anecdotal evidence suggests that this technique is not widely and successfully applied beyond analysis of simple standards, gel bands, and targeted protein immunoprecipitations. Here, we report significant improvements in IMAC-based methodology for enrichment of phosphopeptides from complex biological mixtures. Moreover, we provide detailed explanation for key variables that in our hands most influenced the outcome of these experiments. Our results indicate 5- to 10-fold improvement in recovery of singly- and multiply phosphorylated peptide standards in addition to significant improvement in the number of high-confidence phosphopeptide sequence assignments from global analysis of cellular lysate. In addition, we quantitatively track phosphopeptide recovery as a function of phosphorylation state, and provide guidance for impedance-matching IMAC column capacity with anticipated phosphopeptide content of complex mixtures. Finally, we demonstrate that our improved methodology provides for identification of phosphopeptide distributions that closely mimic physiological conditions.