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.

Revealing the binding modes and the unbinding of 14-3-3σ proteins and inhibitors by computational methods

Revealing The Binding Modes And The Unbinding Of 14-3-3 Sigma Proteins And Inhibitors By Computational Methods

Binding Modes and Conformational Changes of FK506-binding Protein 51 Induced by Inhibitor Bindings: Insight into Molecular Mechanisms Based on Multiple Simulation Technologies

Exploring the Binding Mechanism Between Human Profilin (PFN1) and Polyproline-10 Through Binding Mode Screening

Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods

A Computational Insight into Binding Modes of Inhibitors XD29, XD35, and XD28 to Bromodomain-Containing Protein 4 Based on Molecular Dynamics Simulations.

Understanding the Selectivity of Inhibitors Toward PI4KIIIα and PI4KIIIβ Based Molecular Modeling

Multi-Scaled Explorations of Binding-Induced Folding of Intrinsically Disordered Protein Inhibitor IA3 to Its Target Enzyme

Molecular Dynamics Investigations Suggest a Non-specific Recognition Strategy of 14-3-3σ Protein by Tweezer: Implication for the Inhibition Mechanism.

A Computational Analysis of the Binding Model of MDM2 with Inhibitors

Role of Bridging Water Molecules in GSK3β‐inhibitor Complexes: Insights from QM/MM, MD, and Molecular Docking Studies

Binding Affinity and Dissociation Pathway Predictions for a Series of USP7 Inhibitors with Pyrimidinone Scaffold by Multiple Computational Methods.

Binding Mechanism of Inhibitors to P38α MAP Kinase Deciphered by Using Multiple Replica Gaussian Accelerated Molecular Dynamics and Calculations of Binding Free Energies.

Theoretical Studies on the Selectivity Mechanisms of Glycogen Synthase Kinase 3β (GSK3β) with Pyrazine ATP-competitive Inhibitors by 3DQSAR, Molecular Docking, Molecular Dynamics Simulation and Free Energy Calculations

Interpretation of the Binding Affinities of PTP1B Inhibitors with the MM-GB/SA Method and the X-score Scoring Function.

Molecular Dynamic Simulations Give Insight into the Mechanism of Binding Between 2-Aminothiazole Inhibitors and CDK5

Molecular mechanism of inhibitor bindings to bromodomain-containing protein 9 explored based on molecular dynamics simulations and calculations of binding free energies

Exploring the Binding Mechanism of a Supramolecular Tweezer CLR01 to 14-3-3σ Protein Via Well-Tempered Metadynamics

Molecular Dynamics Simulations Studies and Free Energy Analysis on Inhibitors of Mdm2-P53 Interaction

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

Structure-based and multiple potential three-dimensional quantitative structure-activity relationship (SB-MP-3D-QSAR) for inhibitor design.