P2Y1 receptor (P2Y1R), which belongs to G protein-coupled receptors (GPCRs), is an important target in ADP-induced platelet aggregation. The crystal structure of P2Y1R has been solved recently, which revealed orthosteric and allosteric ligand-binding sites with the details of ligand-protein binding modes. And it suggests that P2Y1R antagonists, which recognize two distinct sites, could potentially provide an efficacious and safe antithrombotic profile. In present paper, 2D similarity search, pharmacophore based screening, and molecular docking were used to explore the potential natural P2Y1R antagonists. 2D similarity search was used to classify orthosteric and allosteric antagonists of P2Y1R. Based on the result, pharmacophore models were constructed and validated by the test set. Optimal models were selected to discover potential P2Y1R antagonists of orthosteric and allosteric sites from Traditional Chinese Medicine (TCM). And the hits were filtered by Lipinski's rule. Then molecular docking was used to refine the results of pharmacophore based screening and analyze the binding mode of the hits and P2Y1R. Finally, two orthosteric and one allosteric potential compounds were obtained, which might be used in future P2Y1R antagonists design. This work provides a reliable guide for discovering natural P2Y1R antagonists acting on two distinct sites from TCM.