The dopamine D3 receptor is an important CNS target for the treatment of a variety of neurological diseases. Selective dopamine D3 receptor antagonists modulate the improvement of psychostimulant addiction and relapse. In this study, five and six featured pharmacophore models of D3R antagonists were generated and evaluated with the post-hoc score combining two survival scores of active and inactive. Among Top 10 models, APRRR215 and AHPRRR104 were chosen based on the coefficient of determination (APRRR215: R2training = 0.80; AHPRRR104: R2training = 0.82) and predictability (APRRR215: Q2test = 0.73, R2predictive = 0.82; AHPRRR104: Q2test = 0.86, R2predictive = 0.74) of their 3D-quantitative structure–activity relationship models. Pharmacophore-based virtual screening of a large compound library from eMolecules (> 3 million compounds) using two optimal models expedited the search process 100-fold speed increase compared to the docking-based screening (HTVS scoring function in Glide) and identified a series of hit compounds having promising novel scaffolds. After the screening, docking scores, as an adjuvant predictor, were added to two fitness scores (from the pharmacophore models) and predicted Ki (from PLSs of the QSAR models) to improve accuracy. Final selection of the most promising hit compounds were also evaluated for CNS-like properties as well as expected D3R antagonism.