In silico screening and study of novel ERK2 inhibitors using 3D QSAR, docking and molecular dynamics

•A pharmacophore model based on known ERK2 inhibitors was proposed.•Five chemical databases were screened using the model.•Docking study was conducted for top predicted ERK2 inhibitors.•Virtual counter screening docking has been conducted.•Molecular dynamics study confirmed stability of the ligand–kinase complex.

ERK2 is a dual specificity protein kinase, part of the Ras/Raf/MEK/ERK signal transduction cascade. It forms an interesting target for inhibition based on its relationship with cell proliferation and oncogenesis. A 3D QSAR pharmacophore model (Hypo1) with high correlation (r = 0.938) was developed for ERK2 ATP site on the basis of experimentally known inhibitors. The model included three hydrogen bonds, and one hydrophobic site. Assessment of Hypo1 through Fisher randomization, cost analysis, leave one out method and decoy test suggested that the model can reliably detect ERK2 inhibitors. Hypo1 has been used for virtual screening of potential inhibitors from ZINC, Drug Bank, NCI, Maybridge and Chembank databases. Using Hypo1 as a query, databases have been interrogated for compounds who meet the pharmacophore features. The resulting hit compounds were subject to docking and analysis. Docking and molecular dynamics analysis showed that in order to achieve a higher potency compounds have to interact with catalytic site, glycine rich loop, Hinge region, Gatekeeper region and ATP site entrance residues. We also identified catalytic site and Glycine rich loop as important regions to bind by molecules for better potency and selectivity.

Graphical abstractFigure optionsDownload full-size imageDownload high-quality image (259 K)Download as PowerPoint slide