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.

Benchmarking the performance of DFT functionals for absorption and fluorescence spectra of EGFR inhibitor AG-1478 using TD-DFT

Insights into Halogen-Induced Changes in 4-Anilinoquinazoline EGFR Inhibitors: A Computational Spectroscopic Study

Optical spectra and conformation pool of tyrosine kinase inhibitor PD153035 using a robust quantum mechanical conformation search

Design, Synthesis, Anti-Tumor Activity, and Molecular Modeling of Quinazoline and Pyrido[2,3-D]Pyrimidine Derivatives Targeting Epidermal Growth Factor Receptor

Molecular Mechanism Study of EGFR Allosteric Inhibitors Using Molecular Dynamics Simulations and Free Energy Calculations

Deducing the Conformational Properties of a Tyrosine Kinase Inhibitor in Solution by Optical Spectroscopy and Computational Chemistry

Predictive power of long-range corrected functionals on the spectroscopic properties of tetrapyrrole derivatives for photodynamic therapy

Benchmarking Charge-Transfer Excited States in TADF Emitters: ∆DFT outperforms TD-DFT for Emission Energies

Can Range-Separated and Hybrid Dft Functionals Predict Low-Lying Excitations? A Tookad Case Study

Explicit solvation modeling in the accurate TD-DFT prediction of absorption spectra for natural nucleobases and fluorescent nucleobase analogues

Toward an Accurate Description of Thermally Activated Delayed Fluorescence: Equal Importance of Electronic and Geometric Factors

Identification of the structural features of quinazoline derivatives as EGFR inhibitors using 3D-QSAR modeling, molecular docking, molecular dynamics simulations and free energy calculations.

The eXact integral simplified time-dependent density functional theory (XsTD-DFT)

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

Broadband X-Ray Absorption Spectra from Time-Dependent Kohn-Sham Calculations

Assessing and applying DFT approaches for geometries and UV–Vis absorption spectra of tetragonal iron(II) complexes

In vitro and in silico evaluation of 4'-hydroxy-[1,1'-biphenyl]-4-carbohydrazide Schiff base and oxadiazole derivatives targeting EGFR allosteric site

QM/MM simulations of EFGR with afatinib reveal the role of the β-dimethylaminomethyl substitution

Comparative DFT Studies of Optoelectronic Properties of MTPA Derivatives

QM/MM Simulations of Afatinib-EGFR Addition: The Role of β-Dimethylaminomethyl Substitution

Methodological Survey of Simplified TD-DFT Methods for Fast and Accurate Interpretation of UV–Vis–NIR Spectra of Phthalocyanines