Pharmacophore-based virtual screening, molecular docking and molecular dynamics simulation for identification of potential ERK inhibitors

As the downstream component of the mitogen-activated protein kinases (MAPK) pathway, the extracellular signal-regulated kinase (ERK) is responsible for phosphorylating a broad range of substrates in cell proliferation, differentiation, and survival. Direct targeting the ERK proteins by the piperidinopyrimidine urea-based inhibitors has been demonstrated to be an effective way to block the MAPK signaling pathway in inhibiting tumor growth. In order to discover better inhibitors, a computer-aided drug design (CADD) approach was employed to reveal the pharmacological characteristics and mechanisms of action. The pharmacophore model was generated on the basis of the compounds with eight features, i.e., four hydrogen bond acceptor atoms, one hydrogen bond donor atom, and three hydrophobic centers. A total of 14 hit compounds were obtained through virtual screening. Two potential inhibitors, namely VS01 and VS02, have been identified by molecular docking and molecular dynamics simulations. Both compounds are capable of attaching to the ERK pocket precisely. The binding free energies of VS01 and VS02 are about 15 kJ/mol and 4 kJ/mol stronger than that of the clinic Ulixertinib because of the characteristic hydrogen bonding, electrostatic, and hydrophilic interactions. The present theoretical investigations shed new light on the rational design of the potential ERK inhibitors to stimulate further experimental tests.Communicated by Ramaswamy H. Sarma