P-glycoprotein (P-gp) is a member of the ATP-Binding Cassette transporter superfamily and mediates transmembrane efflux of many drugs. Since it is involved in multi-drug resistance activity in various cancer cells, the development of P-gp inhibitor is one of the major concerns in anticancer therapy. Human P-gp protein has at least two "functional" drug binding sites that are called "H" site and "R" site, hence it has multi-bindingspecificities. Though the amino acid residues that constitute in drug binding pockets have been proposed by previous experimental evidences, the shapes and the binding poses are not revealed clearly yet. In this study, human P-gp structure was built by homology modeling with available crystal structure of mouse P-gp as a template and docking simulations were performed with inhibitors such as verapamil, hoechst33342, and rhodamine123 to construct the interaction between human P-gp and its inhibitors. The docking simulations were performed 500 times for each inhibitor, and then the interaction frequency of the amino acids at the binding poses was analyzed. With the analysis results, we proposed highly contributing residues that constitute binding pockets of the human P-gp for the inhibitors. Using the highly contributing residues, we proposed the locations and the shapes of verapamil binding site and "R" site, and suggested the possible position of "H" site.