Docking-based 3D-QSAR study of pyridyl aminothiazole derivatives as checkpoint kinase 1 inhibitors

Checkpoint kinase 1 (Chk1) is a promising target for the design of novel anticancer agents. In the present work, molecular docking simulations and three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on pyridyl aminothiazole derivatives as Chk1 inhibitors. AutoDock was used to determine the probable binding conformations of all the compounds inside the active site of Chk1. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were developed based on the docking conformations and alignments. The CoMFA model produced statistically significant results with a cross-validated correlation coefficient (q²) of 0.608 and a coefficient of determination (r²) of 0.972. The reliable CoMSIA model with q² of 0.662 and r² of 0.970 was obtained from the combination of steric, electrostatic and hydrogen bond acceptor fields. The predictive power of the models were assessed using an external test set of 14 compounds and showed reasonable external predictabilities (r²_pred) of 0.668 and 0.641 for CoMFA and CoMSIA models, respectively. The models were further evaluated by leave-ten-out cross-validation, bootstrapping and progressive scrambling analyses. The study provides valuable information about the key structural elements that are required in the rational design of potential drug candidates of this class of Chk1 inhibitors.     

2D and 3D‐QSAR studies on antiproliferative thiazolidine analogs

Two‐dimensional (2D) and three‐dimensional (3D) quantitative structure–activity relationships (QSARs) of 22 thiazolidine analogs with antiproliferative activity expressed as pIC₅₀, which is defined as the negative value of the logarithm of necessary molar concentration of these compounds to cause 50% growth inhibition against melanoma cell lines WM‐164, have been studied by using a combined method of the DFT, MM2 and statistics for 2D, as well as the comparative molecular field analysis (CoMFA) method for 3D. The established 2D‐QSAR model in training set comprised of random 18 compounds shows not only significant statistical quality, but also predictive ability, with the square of adjusted correlation coefficient (R = 0.832) and the square of the cross‐validation coefficient (q² = 0.803). The same model was further applied to predict pIC₅₀ values of the four compounds in the test set, and the resulting R reaching 0.784, further confirms that this 2D‐QSAR model has high predictive ability. The 3D‐QSAR model also shows good correlative and predictive capabilities in terms of R² (0.956) and q² (0.615) obtained from CoMFA model. Further, the robustness of the CoMFA model was verified by bootstrapping analysis (100 runs) with R (0.979) and SD_bs (0.056). It is very interesting to find that the results from 2D‐ and 3D‐QSAR analyses accord with each other, and they all show that the steric interaction plays a crucial role in determining the cytotoxicities of the compounds, and that selecting a moderate‐size or appropriate‐hydrophobicity substituent R as well as increasing the negative charges of C₄ on phenyl ring at the same time are advantageous to improving the cytotoxicity. Such results can offer some useful theoretical references for directing the molecular design and understanding the action mechanism of this kind of compound with antiproliferative activity. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

3D-QSAR study of hallucinogenic phenylalkylamines by using CoMFA approach

The three-dimensional quantitative structure–activity relationship (3D-QSAR) has been studied on 90 hallucinogenic phenylalkylamines by the comparative molecular field analysis (CoMFA). Two conformations were compared during the modeling. Conformation I referred to the amino group close to ring position 6 and conformation II related to the amino group trans to the phenyl ring. Satisfactory results were obtained by using both conformations. There were still differences between the two models. The model based on conformation I got better statistical results than the one about conformation II. And this may suggest that conformation I be preponderant when the hallucinogenic phenylalkylamines interact with the receptor. To further confirm the predictive capability of the CoMFA model, 18 compounds with conformation I were randomly selected as a test set and the remaining ones as training set. The best CoMFA model based on the training set had a cross-validation coefficient q ² of 0.549 at five components and non cross-validation coefficient R ² of 0.835, the standard error of estimation was 0.219. The model showed good predictive ability in the external test with a coefficient R _pre² of 0.611. The CoMFA coefficient contour maps suggested that both steric and electrostatic interactions play an important role. The contributions from the steric and electrostatic fields were 0.450 and 0.550, respectively.     

Molecular description of pyrimidine-based inhibitors with activity against FAK combining 3D-QSAR analysis,molecular docking and molecular dynamics

Focal adhesion kinase (FAK) is a promising target for developing more effective anticancer drugs. To better understand the structure-activity relationships and mechanism of actions of FAK inhibitors, a molecular modeling study using 3D-QSAR, molecular docking, molecular dynamics simulations, and binding free energy analysis were conducted. Two types of satisfactory 3D-QSAR models were generated, comprising the CoMFA model (R²_cv = 0.528, R²_pred = 0.7557) and CoMSIA model (R²_cv = 0.757, R²_pred = 0.8362), for predicting the inhibitory activities of novel inhibitors. The derived contour maps indicate structural characteristics for substituents on the template. Molecular docking, molecular dynamic simulations and binding free energy calculations further reveal that the binding of inhibitors to FAK is mainly contributed from hydrophobic, electrostatic and hydrogen bonding interactions. In addition, some key residues (Arg14, Glu88, Cys90, Arg138, Asn139, Leu141, and Leu155) responsible for ligand-receptor binding are highlighted. All structural information obtained from 3D-QSAR models and molecular dynamics is consist with the available experimental activities. All the results will facilitate the optimization of this series of FAK inhibitors with higher inhibitory activities.     

Studies on molecular mechanism between SHP2 and pyridine derivatives by 3D-QSAR,molecular docking and MD simulations

BackgroundSrc homology 2 (SH2)-containing protein tyrosine phosphatase 2 (SHP2) as a major phosphatase would affect the development of tumors by regulating several cellular processes, and is a significant potential target for cancer treatment.MethodsIn the present work, a series of pyridine derivatives possessing a wide range of inhibitory activity was employed to investigate the structural requirements by developing three dimensional quantitative structure–activity relationship (3D-QSAR) models using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. The results show that CoMFA (R²_cv = 0.646, R²_pred = 0.5587) and CoMSIA (R²_cv = 0.777, R²_pred = 0.7131) have excellent stability and predictability. The relationship between the inhibitory activity and structure of the inhibitors was analyzed by the derived contour maps. Furthermore, the QSAR models were validated by molecular docking and molecular dynamics simulations, which were also applied to reveal the potential molecular mechanism of these inhibitors.FindingsIt was found that Arg110, Asn216, Thr218, Thr252 and Pro490 play a crucial role in stabilizing the inhibitors. Additionally, MM/PBSA calculations provided the binding free energy were also conducted to explain the discrepancy of binding activities. Overall, the outcomes of this work could provide useful information and theoretical guidance for the development of novel and potent SHP2 inhibitors.     

含磷嘧啶类CDK9抑制剂的分子对接、3D-QSAR和分子动力学模拟

选取64个具有潜力的含磷嘧啶类细胞周期依赖性蛋白激酶(CDK9)小分子抑制剂,采用分子对接方法研究了该类小分子与CDK9的结合作用,结果表明,分子构象、氢键形成、疏水性和氨基酸残基Cys106在此类抑制剂与CDK9的结合过程中具有重要作用.在配体叠合的基础上,运用比较分子力场分析(Co MFA)、比较分子相似性指数分析(Co MSIA)和Topomer Co MFA(T-COMFA)研究了分子结构与抑制活性的关系,发现由训练集立体场、静电场和疏水场组合的Co MSIA模型为最优模型,其内部交叉验证相关系数(Q2=0.557)、非交叉验证相关系数(R2=0.959)和外部预测相关系数(r2=0.863)具有统计学意义,该模型的三维等值线图直观显示了化合物的活性与其三维结构的关系.根据这些结果设计了10个具有新结构的含磷嘧啶类化合物,分子对接和分子动力学模拟结果表明,新化合物和CDK9的结合模式与原化合物64相同,自由能分析从理论上证明了新化合物64d的CDK9抑制活性优于化合物64,并且显示含磷基团与残基Asp109的静电场能在化合物与CDK9作用过程中有重要作用.