Pharmacophore modeling, multiple docking, and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G 2 checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhi...

Full description

Saved in:
Bibliographic Details
Published in:Journal of biomolecular structure & dynamics 2019-07, Vol.37 (10), p.2703-2715
Main Authors: Hu, Yanqiu, Zhou, Lu, Zhu, Xiaohong, Dai, Duoqian, Bao, Yinfeng, Qiu, Yaping
Format: Article
Language:English
Subjects:
binding free energy calculations
dynamics
pharmacophore model
quantum polarized ligand docking
Wee1 inhibitors
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G 2 checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R 2  = 0.9212, Q 2  = 0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski's rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50 ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand-protein complex. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified. Communicated by Ramaswamy H. Sarma
ISSN:0739-1102
1538-0254
DOI:10.1080/07391102.2018.1495576