Identification of novel inhibitors of p-hydroxyphenylpyruvate dioxygenase using receptor-based virtual screening

p-Hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD) is a potential target for herbicide discovery, which catalyzes the conversion of p-hydroxyphenylpyruvate to homogentisate. Triketone HPPD inhibitors have been widely applied to weed control for many years. To discover potential HPPD inhibitors with novel scaffold, a virtual screening procedure combining pharmacophore model, molecular docking and molecular dynamics (MD) simulation was carried out. The crystal complex (PDB: 6ISD ) with sulcotrione was selected to establish the receptor-ligand pharmacophore model (CBP). The verified CBP model with key features and interactions was employed to seek structures with potential inhibitory from ZINC database. The best-fit compounds from pharmacophore screening (2050 compounds) were subsequently subjected to molecular docking with CDOCKER method to identify the binding affinity of inhibitors with amino acid residues in active pocket. Six compounds were obtained as the final hits with the Fe(Ш) coordination and π-π interaction with Phe381 and Phe424. MD results demonstrated that Phe424 made great contributions to binding affinities of all the systems. This work provided powerful insights into the design and development of novel HPPD inhibitor in silico techniques. ZINC000255329531 was identified as the most promising candidate due to its high HPPD inhibition effect (ΔGbind = 25.66 kcal mol−1), in comparison with the commercial triketone HPPD herbicide sulcotrione (ΔGbind = 18.94 kcal mol−1).