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.

Some Insights into the Stereochemistry of Inhibition of Macrophage Migration Inhibitory Factor with 2-Fluoro-p-hydroxycinnamate and Its Analogues from Molecular Dynamics Simulations.

A study of the interaction of cinnamate analogues with macrophage migration inhibitory factor (MIF) and P1G mutant from molecular dynamics simulations

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations.

Mechanistic insight into the mode of inhibition of dietary flavonoids; targeting macrophage migration inhibitory factor

Effect of Polarization on HIV-1protease and Fluoro-Substituted Inhibitors Binding Energies by Large Scale Molecular Dynamics Simulations

Discovery of Novel Inhibitors Targeting the Macrophage Migration Inhibitory Factor Via Structure-Based Virtual Screening and Bioassays.

An Analysis of MIF Structural Features that Control Functional Activation of CD74

From docking to dynamics: Unveiling the potential non-peptide and non-covalent inhibitors of M pro from natural products

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

Repurposing old drugs as novel inhibitors of human MIF from structural and functional analysis

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

Molecular Modeling and Molecular Dynamics Simulation Studies on Pyrrolopyrimidine-Based Α-Helix Mimetic As Dual Inhibitors of MDM2 and MDMX

Molecular docking and molecular dynamic simulations of inhibitor ZJ0777 and CIE binding to acetohydroxyacid synthases

Energetic Analysis of the Two Controversial Drug Binding Sites of the M2 Proton Channel in Influenza A Virus.

Binding Mechanism of Inhibitors to SARS-CoV-2 Main Protease Deciphered by Multiple Replica Molecular Dynamics Simulations.

Identification and Structure-Activity Relationships of Dietary Flavonoids As Human Macrophage Migration Inhibitory Factor (MIF) Inhibitors.

Free Energy Calculations on the Two Drug Binding Sites in the M2 Proton Channel

Structural interactions dictate the kinetics of macrophage migration inhibitory factor inhibition by different cancer-preventive isothiocyanates.

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

Computational Study on the Inhibition Mechanism of a Cyclic Peptide MaD5 to PfMATE: Insight from Molecular Dynamics Simulation, Free Energy Calculation and Dynamical Network Analysis

Discovery of Covalent Inhibitors for MIF Tautomerase Via Cocrystal Structures with Phantom Hits from Virtual Screening.