Abstract:In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.

Pharmacophore Guided 3D‐QSAR CoMFA Analysis of Amino Substituted Nitrogen Heterocycle Ureas As KDR Inhibitors

Pyrrolo[3,2-d]pyrimidine Derivatives As Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies

3D-QSAR Study on 2,3-Dihydroimidazo[4,5]-pyridin-2-one Derivatives with a Meta Substitution Pattern As V600E BRAF Inhibitors

Structure-based Rational Quest for Potential Novel Inhibitors of Human HMG-CoA Reductase by Combining CoMFA 3D QSAR Modeling and Virtual Screening.

3D-QSAR Study for Inhibitors of Cyclin-Dependent Kinase 4

QSAR analysis of VEGFR-2 inhibitors based on machine learning, Topomer CoMFA and molecule docking

3D-QSAR and 3D-QSSR Studies of Thieno[2,3-D]pyrimidin-4-yl Hydrazone Analogues As CDK4 Inhibitors by CoMFA Analysis

Molecular description of pyrimidine-based inhibitors with activity against FAK combining 3D-QSAR analysis, molecular docking and molecular dynamics

Mapping the Binding Site of a Large Set of Quinazoline Type EGF-R Inhibitors Using Molecular Field Analyses and Molecular Docking Studies.

3D-QSAR modeling and molecular docking studies on a series of 2, 4, 5-trisubstituted imidazole derivatives as CK2 inhibitors

Effect of Substituents in Different Positions of Aminothiazole Hinge-Binding Scaffolds on Inhibitor–CDK2 Association Probed by Interaction Entropy Method

Docking-based CoMFA and CoMSIA Studies on Naphthyl-Substituted Diarylpyrimidines As NNRTIs.

Three-dimensional quantitative structure–activity relationship study for analogues of TQXs using CoMFA and CoMSIA

3D-QSAR Studies on 4-Hydroxyphenylpyruvate Dioxygenase Inhibitors by Comparative Molecular Field Analysis (Comfa)

Design of Novel FLT-3 Inhibitors Based on Dual-Layer 3D-QSAR Model and Fragment-Based Compounds .

Three-dimensional Quantitative Structure–activity Relationship Study on Paullones As CDK Inhibitors Using CoMSIA and CoMFA

Three-Dimensional Quantitative Structure-Activity Relationship Analyses of A Series of Butenolide Eta Antagonists

Quantitative structure activity relationship and molecular simulations for the exploration of natural potent VEGFR-2 inhibitors: an in silico anti-angiogenic study

Pharmacophore modeling and in silico screening for new KDR kinase inhibitors.

The Discovery of Novel Vascular Endothelial Growth Factor Receptor Tyrosine Kinases Inhibitors: Pharmacophore Modeling, Virtual Screening and Docking Studies

In silico design novel dihydropyrimio[4, 5-d]pyrimidine derivatives as inhibitors for colony-stimulating factor-1 receptor based on 3D-QSAR, molecular docking and molecular dynamics simulation