mGluR5, which belongs to G-protein-coupled receptor (GPCR) superfamily, is believed to be associated with many human diseases, such as a wide range of neurological disorders, gastroesophageal reflux disease (GERD) and cancer. Comparing with compounds target on orthosteric binding site, significant roles have been established for mGluR5 negative allosteric modulators (NAMs) due to its higher subtype-selectivity and more suitable pharmacokinetic profiles. Nevertheless, to date, none of them were come into market for various reasons.In this study, 3D quantitative pharmacophore model was generated by using HypoGen module in DS. With several validation methods ultilized, the optimal pharmacophore model Hypo2 was selected to discover potential mGluR5 NAM s from natural products. 217 potential NAMs were obtained after filtered by Lipinski's rule (≥4). Then, molecular docking was used to refine the pharmacophore-based screening results and analyze the binding mode of NAM s and mGluR5. 3 compounds, Aglaiduline, 5-O-Ethyl-hirsutanonol, and Yakuchinone A, which with good ADM ET properties, acceptable Fit value and estimate value, and high docking score, were reserved for molecular dynamics (M D) simulation study. All of them have stability of ligand binding. From our computational results, there might exhibit drug-like negative allosteric moderating effects on mGluR5 in these natural products. This work provides a reliable method for discovering mGluR5 NAMs from natural products.