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.     

3D-QSAR study,docking molecular and simulation dynamic on series of benzimidazole derivatives as anti-cancer agents

Benzimidazole is an important heterocyclic organic compound which has a structural analogy to nucleotides found in human body and hence is an important pharmacophore in medicinal chemistry. The anti-cancer activities for a diverse set of benzimidazole as anti-cancer agents against breast cancer cell line (MCF7) assay have been subjected to 3D-QSAR (3-Dimensional Quantitative Structural-Activity Relationship) studies. Both CoMFA and CoMSIA models exhibit significant results in terms of statistical parameters as determination coefficients R² > 0.9 and Leave One Out cross-validation determination coefficients Q²> 6. The predictive quality of both 3D QSAR models have been assessed by external validation and Y-randomization test. Five new compounds have been designed and predicted by in silico ADMET method. In the second part, we have used the docking molecular and simulation dynamics (MD) to investigate the bonding interactions and stability of the designed compounds into the Pin1. Then, we have compared them to Trastuzumab and Tamoxifen as a standard inhibitors drug of breast cancer. The designed compounds form stable hydrogen and hydrophobic bonding interactions with the residues Lys63, Gln131, Ser154, Arg 68 and Arg69 of Pin1 receptor during 100 ns as a time of the simulation. The obtained results showed that the new benzimidazole are useful as a template for future design of more potent inhibitors against breast cancer cell lines (MCF7).     

3D-QSAR and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q² = 0.589 and r² = 0.927 and q² = 0.473 and r² = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.     

Molecular modelling studies on cinnoline-based BTK inhibitors using docking and structure-based 3D-QSAR

ABSTRACT

BTK inhibitors have been proved as an effective target for B-cell malignancies. Ibrutinib is the most advanced irreversible BTK inhibitor for treating mantle cell lymphoma/chronic lymphocytic leukaemia but with existing drug resistance and adverse effects. To design novel effective and safety reversible BTK inhibitors, 115 newly cinnoline analogues were selected to perform molecular docking and 3D-QSAR study because of the main scaffold similarity to Ibrutinib. Both established CoMFA and CoMSIA models obtained high predictive and satisfactory value. CoMFA/CoMSIA contour maps demonstrated that bulky substitutions are preferred at R₁ and R₃ positions, and introducing hydrophilic and negative electrostatic substitutions at R₁ positions is important for improving BTK inhibitory activities. These results will be useful to provide clues for rationally designing novel and high potency BTK inhibitors.     

苯并嗪酮衍生物的3D-QSAR分析

苯并嗪酮衍生物是近年来发现的一类抗血小板聚集化合物, 在前人研究的基础上利用比较分子场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)对23个苯并嗪酮衍生物进行了三维定量构效关系(3D-QSAR)研究. 其中CoMFA模型交叉验证系数Q2=0.703, 回归系数R2=0.994, 计算值与实验值的平均方差SEE=0.053, 统计方差比F=184.773; CoMSIA模型Q2=0.847, R2=0.992, SEE=0.058, F=171.670. 两种方法得到的模型都具有较好的预测能力. 结果表明, 标题化合物中8-位取代基R1静电效应起主要作用; 2-位取代基R2立体效应占主导作用, 但官能团大小要适中. 根据研究结果设计了六种活性较高的化合物.     

A 3D-QSAR CoMSIA study on 3-azolylmethylindoles as anti-leishmanial agents

A three-dimensional quantitative structure-activity relationship (3D-QSAR) study using Comparative Molecular Similarity Indices Analysis (CoMSIA) was conducted on a series of 3-azolylmethylindoles as anti-leishmanial agents. Evaluation of 24 compounds synthesized in our laboratory served to establish the model. A random search was performed on the library of compounds, and molecules of the training set were aligned on common elements of template molecule 13, one of the most active compounds. The best predictions were obtained from multifit procedure with a CoMSIA model combining steric, electrostatic, hydrophobic and hydrogen bond acceptor fields (q ²?=?0.594, r ²?=?0.897). The model was validated using an external test set of 7 compounds giving a satisfactory predictive r ² value of 0.649. Information obtained from CoMSIA contour maps could be used for further design of more promising inhibitors.     

马来酰胺类糖原合成酶激酶-3β抑制剂的分子对接和三维定量构效关系

通过分子对接和三维定量构效关系(3D-QSAR)两种方法来确定两类马来酰胺类的糖原合成酶激酶-3β(GSK-3β)抑制剂的结合方式. 首先, 用分子对接确定抑制剂与GSK-3β结合模式及其相互作用; 然后用比较分子力场分析法(CoMFA)与比较分子相似性指数分析法(CoMSIA)对48个化合物做三维定量构效关系的分析. 两种方法得出的交互验证回归系数分别为0.669(CoMFA)和0.683(CoMSIA), 证明该模型具有很好的统计相关性, 同时也说明该模型具有较高的预测能力.根据该模型提供的信息, 设计出9个预测活性较好的分子.     

Computational Analysis of CRTh2 receptor antagonist: A Ligand-based CoMFA and CoMSIA approach

CRTh2 receptor is an important mediator of inflammatory effects and has attracted much attention as a therapeutic target for the treatment of conditions such as asthma, COPD, allergic rhinitis and atopic dermatitis. In pursuit of better CRTh2 receptor antagonist agents, 3D-QSAR studies were performed on a series of 2-(2-(benzylthio)-1H-benzo[d]imidazol-1-yl) acetic acids. There is no crystal structure information available on this protein; hence in this work, ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed by atom by atom matching alignment using systematic search and simulated annealing methods. The 3D-QSAR models were generated with 10 different combinations of test and training set molecules, since the robustness and predictive ability of the model is very important. We have generated 20 models for CoMFA and 100 models for CoMSIA based on two different alignments. Each model was validated with statistical cut off values such as q² > 0.4, r² > 0.5 and r²_pred > 0.5. Based on better q² and r²_pred values, the best predictions were obtained for the CoMFA (model 5 q² = 0.488, r²_pred = 0.732), and CoMSIA (model 45 q² = 0.525, r²_pred = 0.883) from systematic search conformation alignment. The high correlation between the cross-validated/predicted and experimental activities of a test set revealed that the CoMFA and CoMSIA models were robust. Statistical parameters from the generated QSAR models indicated the data is well fitted and have high predictive ability. The generated models suggest that steric, electrostatic, hydrophobic, H-bond donor and acceptor parameters are important for activity. Our study serves as a guide for further experimental investigations on the synthesis of new CRTh2 antagonist.     

皂甙的三维定量构效关系研究

针对目标分子柔性大的特点,在比较分子场分析(CoMFA)方法中采用交叉验证相关系数平方R^2引导的构象选择法。对12个皂甙分子的生物活性进行了三维定量构效关系研究。探讨了几种探针对构效关系结果的影响,并选择了一种较合理的“复合”探针方案。应用该复合探针构建CoMFA模型,发现影响药效的立体场与静电场的贡献分别为40%和40%,其它能量项的贡献为20%。该模型交叉验证的相关系数平方R^2为0.653,非交叉验证的R^2为0.991,方差比F(4,7)值130.195(即置信度99%以上),活性预计的标准偏差与极差比(s/△γ)为4.2%,表明模型具有较好的预测能力。根据该模型,预计在指定位置添加位阻较大的基团活性值提高将会比较明显。     

Additional synthesis on thiophene-containing trisubstituted methanes (TRSMs) as inhibitors of M. tuberculosis and 3D-QSAR studies

Abstract

We earlier reported thiophene-containing trisubstituted methanes (TRSMs) as novel cores carrying anti-tubercular activity, and identified S006-830 as the phenotypic lead with potent bactericidal activity against single- and multi-drug resistant clinical isolates of Mycobacterium tuberculosis (M. tb). In this work, we carried out additional synthesis of several TRSMs. The reaction scheme essentially followed the Grignard reaction and Friedel–Crafts alkylation, followed by insertion of a dialkylaminoethyl chain. We also performed microbiological evaluations including in vitro screening against the virulent strain M. tb H37Rv, cytotoxicity assessment in the Vero C-1008 cell line, and 3D-QSAR studies with comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). CoMFA and CoMSIA models yielded good statistical results in terms of q² and r² values, suggesting the validity of the models. It was concluded that a para-substituted benzene ring with bulkier electron-donating groups and aminoalkyl chains are required for higher inhibitory capacity against M. tuberculosis. We believe that these insights will rationally guide the design of newer, optimal, TRSMs.     

Glycogen synthase kinase-3 inhibition by 3-anilino-4-phenylmaleimides: insights from 3D-QSAR and docking

Glycogen synthase kinase-3, a serine/threonine kinase, has been implicated in a wide variety of pathological conditions such as diabetes, Alzheimer’s disease, stroke, bipolar disorder, malaria and cancer. Herein we report 3D-QSAR analyses using CoMFA and CoMSIA and molecular docking studies on 3-anilino-4-phenylmaleimides as GSK-3α inhibitors, in order to better understand the mechanism of action and structure-activity relationship of these compounds. Comparison of the active site residues of GSK-3α and GSK-3β isoforms shows that all the key amino acids involved in polar interactions with the maleimides for the β isoform are the same in the α isoform, except that Asp133 in the β isoform is replaced by Glu196 in the α isoform. We prepared a homology model for GSK-3α, and showed that the change from Asp to Glu should not affect maleimide binding significantly. Docking studies revealed the binding poses of three subclasses of these ligands, namely anilino, N-methylanilino and indoline derivatives, within the active site of the β isoform, and helped to explain the difference in their inhibitory activity.     

新型三唑类抗真菌化合物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 系统研究了40个新型三唑类化合物抗真菌活性的三维定量构效关系. 在CoMFA研究中, 研究了两种药效构象对模型的影响, 并考察了网格点步长对统计结果的影响. 在CoMSIA研究中, 系统考察了各种分子场组合、网格点步长和衰减因子对模型统计结果的影响, 发现立体场、静电场、疏水场和氢键受体场的组合得到最佳模型. 所建立CoMFA和CoMSIA模型的交叉相关系数q²值分别为0.718和0.655, 并都具有较强的预测能力. CoMFA和CoMSIA模型的三维等值线图直观地解释了化合物的构效关系, 阐明了化合物结构中苯环上各位置取代基对抗真菌活性的影响, 为进一步结构优化提供了重要依据.     

CBP/P300溴结构域联芳基类抑制剂三维定量构效关系研究

CREB结合蛋白(CBP)和与其高度同源的P300蛋白是组蛋白乙酰化酶的两个亚型,两者通过它们的溴结构域(bromodomain,BRD)与染色质结合,目前,CBP/P300已经成为人类在肿瘤靶点领域中的研究热点。本研究基于CBP/P300溴结构域联芳基类抑制剂建立三维定量构效关系,采用比较分子力场分析法(Co MFA)和比较分子相似性指数分析法(Co MSIA)分别建立35个已知活性抑制剂的3D-QSAR模型,以确定CBP/P300溴结构域联芳基类抑制剂分子结构与生物活性之间的定量关系。Co MFA和Co MSIA模型活性数据p IC50的预测值与实验值基本一致,说明这两个模型具有较高的预测能力和统计学意义。根据Co MFA和Co MSIA模型所提供的立体场、静电场、疏水场、氢键给体场、氢键供体场等信息提出了改善此类抑制剂活性的药物设计思路,为指导设计具有更高活性的新分子和预测更加有效的CBP/P300溴结构域抑制剂提供理论依据。     

A combined CoMFA and CoMSIA 3D-QSAR study of benzamide type antibacterial inhibitors of the FtsZ protein in drug-resistant Staphylococcus aureus

A major problem today is bacterial resistance to antibiotics and the small number of new therapeutic agents approved in recent years. The development of new antibiotics capable of acting on new targets is urgently required. The filamenting temperature-sensitive Z (FtsZ) bacterial protein is a key biomolecule for bacterial division and survival. This makes FtsZ an attractive new pharmacological target for the development of antibacterial agents. There have been several attempts to develop ligands able to inhibit FtsZ. Despite the large number of synthesized compounds that inhibit the FtsZ protein, there are no quantitative structure–activity relationships (QSAR) that allow for the rational design and synthesis of promising new molecules. We present the first 3D-QSAR study of a large and diverse set of molecules that are able to inhibit the FtsZ bacterial protein. We summarize a set of chemical changes that can be made in the steric, electrostatic, hydrophobic and donor/acceptor hydrogen-bonding properties of the pharmacophore, to generate new bioactive molecules against FtsZ. These results provide a rational guide for the design and synthesis of promising new antibacterial agents, supported by the strong statistical parameters obtained from CoMFA (r²_pred = 0.974) and CoMSIA (r²_pred = 0.980) analyses.     

CoMFA and CoMSIA 3D-quantitative structure-activity relationship model on benzodiazepine derivatives, inhibitors of phosphodiesterase IV

Recently, we reported structurally novel PDE4 inhibitors based on 1,4-benzodiazepine derivatives. The main interest in developing bezodiazepine-based PDE4 inhibitors is in their lack of adverse effects of emesis with respect to rolipram-like compounds. A large effort has thus been made toward the structural optimization of this series. In the absence of structural information on the inhibitor binding mode into the PDE4 active site, 2D-QSAR (H-QSAR) and two 3D-QSAR (CoMFA and CoMSIA) methods were applied to improve our understanding of the molecular mechanism controlling the PDE4 affinity of the benzodiazepine derivatives. As expected, the CoMSIA 3D contour maps have provided more information on the benzodiazepine interaction mode with the PDE4 active site whereas CoMFA has built the best tool for activity prediction. The 2D pharmacophoric model derived from CoMSIA fields is consistent with the crystal structure of the PDE4 active site reported recently. The combination of the 2D and 3D-QSAR models was used not only to predict new compounds from the structural optimization process, but also to screen a large library of bezodiazepine derivatives.