Protein glycosylation and phosphorylation play essential roles in biological systems. The crosstalk of both glycosylation and phosphorylation on one protein represents an unveiled biological function. To realize the analyses of both glycopeptides and phosphopeptides, a simultaneous enrichment method of N-glycopeptides, mono-phosphopeptides and multi-phosphopeptides was developed based on a multi-functional dual-metal centered zirconium metal-organic framework that provided multiple interactions for HILIC, IMAC, and MOAC for glycopeptides and phosphopeptides. Based on a careful optimization of sample loading and elution conditions for the simultaneous enrichment of glycopeptides and phosphopeptides with the zirconium metal-organic framework, a total of 1011 N-glycopeptides derived from 410 glycoproteins and 1996 phosphopeptides including 741 multi-phosphopeptides derived from 1189 phosphoproteins could be identified from a HeLa cell digest. The simultaneous enrichment approach for glycopeptides and mono-/multi-phosphopeptides demonstrates the great potential of the combined interactions for HILIC, IMAC, and MOAC in integrated post-translational modification proteomics research.