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.

Molecular Docking of 3-Methylindole-containing Drugs Binding into Cyp3a4

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

Ensemble Docking Approach to Mitigate Pregnane X Receptor-Mediated CYP3A4 Induction Risk

Computational Investigation of Ligand Binding to the Peripheral Site in CYP3A4: Conformational Dynamics and Inhibitor Discovery.

Investigation of Indazole Unbinding Pathways in CYP2E1 by Molecular Dynamics Simulations

A Comprehensive Analysis of Human CYP3A4 Crystal Structures as a Potential Tool for Molecular Docking-Based Site of Metabolism and Enzyme Inhibition Studies

Docking and 3D-QSAR study on p38α MAP kinase inhibitors

X-ray Crystallography, DFT Calculations and Molecular Docking of Indole-Arylpiperazine Derivatives As Α1a-Adrenoceptor Antagonists

Molecular Docking, 3D-QSAR and Structural Optimization of Indole Biphenylcarboxylic Acids as PPARγ Antagonists

Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism

Molecular Docking, QSAR and Molecular Dynamics Simulation on Spiro-oxindoles As MDM2 Inhibitors

Structure Determinants of Indolin-2-on-3-spirothiazolidinones As MptpB Inhibitors: an in Silico Study

MD Investigation on the Differences in the Dynamic Interactions Between the Specific Ligand Azamulin and Two CYP3A Isoforms, 3A4 and 3A5.

Differences of Atomic-Level Interactions between Midazolam and Two CYP Isoforms 3A4 and 3A5

Combined 3d-Qsar Modeling And Molecular Docking Studies On Pyrrole-Indolin-2-Ones As Aurora A Kinase Inhibitors

Molecular Docking and 3D-QSAR Studies of 2-Substituted 1-Indanone Derivatives As Acetylcholinesterase Inhibitors.

A combination of pharmacophore modeling, atom-based 3D-QSAR, molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Pharmacophore, Qsar, and Binding Mode Studies of Substrates of Human Cytochrome P450 2d6 (cyp2d6) Using Molecular Docking and Virtual Mutations and an Application to Chinese Herbal Medicine Screening

Molecular Modeling of Human Cytochrome P450 2W1 and Its Interactions with Substrates.

Design of 5H-pyrrolo[2,3-b] pyrazine-2-phenyl Ether Derivatives as Potential JAK3 Inhibitor Based on Surflex Docking, 3D-QSAR Modeling and Reverse Docking

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