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.

Quantitative clustering and matching of conformers using average electron densities

Average Electron Density: A Quantitative Tool for Evaluating Non-Classical Bioisosteres of Amides

Clustering of atoms relative to vector space in the Z-matrix coordinate system and 'graphical fingerprint' analysis of 3D pharmacophore structure

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

Modern Semiempirical Electronic Structure Methods and Machine Learning Potentials for Drug Discovery: Conformers, Tautomers, and Protonation States

Sorting drug conformers in enzyme active sites: the XTB way

Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design

Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes.

A Machine Learning-Based High-Precision Density Functional Method for Drug-Like Molecules

First-Principles Molecular Structure Search with a Genetic Algorithm

An artificial intelligence approach for tackling conformational energy uncertainties in chiroptical spectroscopies

Discrimination of Aminobiphenyl Isomers in the Gas Phase and Investigation of Their Complex Conformations

Abstract Wrk2-04: Virtual screening of ultra-large chemical spaces for novel chemotype discovery

Quantitatively integrating molecular structure and bioactivity profile evidence into drug-target relationship analysis

Molecular Gas Phase Conformational Ensembles

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

CASTELO: Clustered Atom Subtypes aidEd Lead Optimization -- a combined machine learning and molecular modeling method

An optimization algorithm for conformer generation based on the bond contribution ranking

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

Drug Delivery Process Simulation - Quantifying the Conformation Dynamics of Paclitaxel and Cremophor EL