Discovery of potential negative allosteric modulators of mGluR5 from natural products using pharmacophore modeling, molecular docking, and molecular dynamics simulation studies

mGluR5, which belongs to the G-protein-coupled receptor superfamily, is believed to be associated with many human diseases, such as a wide range of neurological disorders, gastroesophageal reflux disease, and cancer. Comparing with compounds that target on the orthosteric binding site, significant roles have been established for mGluR5 negative allosteric modulators (NAMs) due to their higher subtype selectivity and more suitable pharmacokinetic profiles. Nevertheless, to date, none of them have come to market for various reasons. In this study, a 3D quantitative pharmacophore model was generated by using the HypoGen module in Discovery Studio 4.0. With several validation methods ultilized, the optimal pharmacophore model Hypo2 was selected to discover potential mGluR5 NAMs from natural products. Two hundred and seventeen potential NAMs were obtained after being filtered by Lipinski's rule (>= 4). Then, molecular docking was used to refine the pharmacophore-based screening results and analyze the binding mode of NAMs and mGluR5. Three compounds, aglaiduline, 5-Omicron-ethyl-hirsutanonol, and yakuchinone A, with good ADMET properties, acceptable Fit value and estimated value, and high docking score, were reserved for a molecular dynamics 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.