Homology modeling of G-protein-coupled seven-transmembrane receptors (GPCRs) remains a challenge despite the increasing number of released GPCR crystal structures. This challenge can be attributed to the low sequence identity and structural diversity of the ligand-binding pocket of GPCRs. We have developed an optimized GPCR structure modeling method based on multiple GPCR crystal structures. This method was designed to be applicable to distantly related receptors of known structural templates. CXC chemokine receptor (CXCR7) is a potential drug target for cancer chemotherapy. Homology modeling, docking, and virtual screening for CXCR7 were carried out using our method. The predicted docking poses of the known antagonists were different from the crystal structure of human CXCR4 with the small-molecule antagonist IT1t. Furthermore, 21 novel CXCR7 ligands with IC50 values of 1.29-11.4 μM with various scaffolds were identified by structure-based virtual screening.