The dopamine D3 receptor (D3R) is an important central nervous system target for treating various neurological diseases. D3R antagonists modulate the improvement of psychostimulant addiction and relapse, while D3R agonists can enhance the response to dopaminergic stimulation and have potential applications in treating Parkinson’s disease, which highlights the importance of identifying novel D3R ligands. Therefore, we performed auto dock Vina-based virtual screening and D3R-binding-affinity assays to identify human D3R ligands with diverse structures. All molecules in the ChemDiv library (>1,500,000) were narrowed down to a final set of 37 molecules for the binding assays. Twenty-seven compounds exhibited over 50% inhibition of D3R at a concentration of 10 μM, and 23 compounds exhibited over 70% D3R inhibition at a concentration of 10 μM. Thirteen compounds exhibited over 80% inhibition of D3R at a concentration of 10 μM and the IC₅₀ values were measured. The IC₅₀ values of the five compounds with the highest D3R-inhibition rates ranged from 0.97 μM to 1.49 μM. These hit compounds exhibited good structural diversity, which prompted us to investigate their D3R-binding modes. After trial and error, we combined unbiased molecular dynamics simulation (MD) and molecular mechanics generalized Born surface area (MM/GBSA) binding free-energy calculations with the reported protein–ligand-binding pose prediction method using induced-fit docking (IFD) and binding pose metadynamics (BPMD) simulations into a self-consistent and computationally efficient method for predicting and verifying the binding poses of the hit ligands to D3R. Using this IFD-BPMD-MD-MM/GBSA method, we obtained more accurate and reliable D3R–ligand-binding poses than were obtained using the reported IFD-BPMD method. This IFD-BPMD-MD-MM/GBSA method provides a novel paradigm and reference for predicting and validating other protein–ligand binding poses.