Combined 3D-QSAR modeling and molecular docking study on indolinone derivatives as inhibitors of 3-phosphoinositide-dependent protein kinase-1

3-Phosphoinositide-dependent protein kinase-1 (PDK1) is a promising target for developing novel anticancer drugs. In order to understand the structure-activity correlation of indolinone-based PDK1 inhibitors, we have carried out a combined molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling study. The study has resulted in two types of satisfactory 3D-QSAR models, including the CoMFA model (r(2)=0.907; q(2)=0.737) and CoMSIA model (r(2)=0.991; q(2)=0.824), for predicting the biological activity of new compounds. The detailed microscopic structures of PDK1 binding with inhibitors have been studied by molecular docking. We have also developed docking-based 3D-QSAR models (CoMFA with q(2)=0.729; CoMSIA with q(2)=0.79). The contour maps obtained from the 3D-QSAR models in combination with the docked binding structures help to better interpret the structure-activity relationship. All of the structural insights obtained from both the 3D-QSAR contour maps and molecular docking are consistent with the available experimental activity data. This is the first report on 3D-QSAR modeling of PDK1 inhibitors. The satisfactory results strongly suggest that the developed 3D-QSAR models and the obtained PDK1-inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and in future drug design.     

3 D-QSAR Analysis of Agonists of nAChRs： Epibatidine Analogues

A 3 D-QSAR about nAChRs agonists epibatidine analogues was performed using theCoMFA and CoMSIA. The correlation coefficients were R2cv = 0.546, R2cv = 0.907 in CoMFA andR2cv = 0.655, R2,~ = 0.962 in CoMSIA of the final model. The prediction using the final models tothe test set was r2 = 0.675 in CoMFA and r2 = 0.462 in CoMSIA. This model will be useful in thedesign of novel compounds with high affinity.     

Receptor-based modeling and 3D-QSAR for a quantitative production of the butyrylcholinesterase inhibitors based on genetic algorithm

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for a series of structurally related steroidal alkaloids as butyrylcholinesterase (BuChE) inhibitors. Docking studies were employed to position the inhibitors into the BuChE active site to determine the most probable binding mode. The strategy was to explore multiple inhibitor conformations in producing a more reliable 3D-QSAR model. These multiple conformations were derived using the FlexS program. The conformation selection step for CoMFA was done by genetic algorithm. The genetic algorithm based CoMFA approach was found to be the best. Both CoMFA and CoMSIA yielded significant cross-validated q(2) values of 0.701 and 0.627 and the r(2) values of 0.979 and 0.982, respectively. These statistically significant models were validated by a test set of five compounds. Comparison of CoMFA and CoMSIA contour maps helped to identify structural requirements for the inhibitors and serves as a basis for the design of the next generation of the inhibitor analogues. The results demonstrate that the combination of ligand-based and receptor-based modeling with use of a genetic algorithm is a powerful approach to build 3D-QSAR models. These data can be used for the lead optimization process with respect to inhibition enhancement which is important for the drug discovery and development for Alzheimer's disease.     

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溴结构域抑制剂提供理论依据。     

3D-QSAR and Surflex Docking Studies of a Series of Alkaline Phosphatase Inhibitors

Alkaline phosphatases(APs) include the placental AP(PLAP), germ cell AP(GCAP), intestinal AP(IAP) and tissue nonspecific AP(TNAP). Over expression of TNAP in smooth muscle cells of kidney and vessels provokes the progress of such serious diseases as end-stage renal disease, idiopathic infantile arterial calcification, ankylosis, osteoarthritis and diabetes. In order to design and optimize the potent TNAP inhibitors, comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA) were used to analyze 3D structure-activity relationships(3D-QSAR) of TNAP inhibitors. The 3D-QSAR model(CoMFA with q~2 = 0.521, r~2 = 0.930; CoMSIA with q~2 = 0.529, r~2 = 0.933) had a good predictability. Surflex-dock was used to reveal the binding mode between the inhibitors and TNAP protein. CoMFA, CoMSIA and docking results provide guidance for the discovery of TNAP inhibitors. Finally, eight new compounds as potential TNAP inhibitors were designed.     

含呋喃环双酰脲类衍生物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 对27个新型双酰基脲类化合物的杀蚊幼虫(Aedes aegypti L.)活性进行三维定量构效关系(3D-QSAR)研究. 在CoMFA研究中, 考察了网格点步长对统计结果的影响. 在CoMSIA研究中, 系统考察了各种分子场组合、网格点步长和衰减因子对模型统计结果的影响, 发现立体场和氢键供体场的组合得到最佳模型. 所建立的CoMFA和CoMSIA模型的非交叉验证相关系数r2值分别为0.828和0.841, 并都具有较强的预测能力. CoMFA和CoMSIA模型的三维等值图不仅直观地解释了结构与活性的关系, 而且为后续优化该系列化合物提供了理论依据.     

STAT3环烷烃并噻吩类衍生物抑制剂的三维定量构效关系研究

本文对STAT3抑制剂的化学结构与生物活性之间的关系进行研究。采用三维定量构效关系(3D-QSAR)中的比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)方法针对52个STAT3抑制剂建立3D-QSAR模型,阐明了抑制剂化学结构与其生物活性之间的关系。所构建的CoMFA模型交叉验证系数为0.548,非交叉验证系数为0.754,标准偏差为0.278,显著系数为58.297;所构建的CoMSIA模型交叉验证系数为0.892,非交叉验证系数为0.597,标准偏差为0.192,显著系数为57.794。结果显示CoMFA和CoMSIA模型具有良好的稳定性和预测能力。3D-QSAR模型等势图提供的相关场信息对新型STAT3抑制剂的设计具有指导意义。     

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

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

苯并呋喃类N-肉豆蔻酰基转移酶抑制剂的三维定量构效关系研究

摘要采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 系统地研究了40个苯并呋喃类N-肉豆蔻酰基转移酶(NMT)抑制剂的三维定量构效关系. 在CoMFA研究中, 考察了网格点步长对模型统计结果的影响. 在CoMSIA研究中, 研究了各种分子场组合、 网格点步长和衰减因子对模型统计结果的影响, 发现立体场、 静电场、 疏水场和氢键受体场的组合可得到最佳模型. 所建立的CoMFA和CoMSIA模型的交叉相关系数q2值分别为0.759和0.730, 均具有较强的预测能力. 利用CoMFA和CoMSIA模型的三维等值线图直观地解释了化合物的构效关系, 阐明了化合物结构中苯并呋喃环上各位置取代基对抑酶活性的影响, 为进一步结构优化提供了重要依据.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

3D-QSAR Analysis of a Series of Dihydroquinolizinone Derivatives as a Hepatitis B Virus Expression Inhibitor

In this study, we explored a three-dimensional quantitative structure-activity relationship(3D-QSAR) model of 63 HBV viral gene expression inhibitors containing dihydroquinolizinones. Two high predictive QSAR models have been built, including comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA). The internal validation parameter(CoMFA, q~2 = 0.701, r~2 = 0.999; CoMSIA, q~2 = 0.721, r~2 = 0.998) and external validation parameter(CoMFA, r~2_(pred = 0.999); CoMSIA, r~2_(pred = 0.999)) indicated that the models have good predictive abilities and significant statistical reliability. We designed several molecules with potentially higher predicted activity on the basis of the result of the models. This work might provide useful information to design novel HBV viral gene expression inhibitors.     

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.     

3D-QSAR Studies on a Series of Indoleamide Derivatives as Antiplasmodial Drugs

In this study, three-dimensional quantitative structure-activity relationship(3D-QSAR) was studied for the antiplasmodial activity of a series of novel indoleamide derivatives by comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(Co MSIA). 3D-QSAR model was established by a training set of 20 compounds and was externally validated by a test set of 4 compounds. The best prediction(Q~2 = 0.593 and 0.527, R~2 = 0.990 and 0.953, r_(pred)~2 = 0.967 and 0.962 for CoMFA and CoMSIA) was obtained according to CoMFA and CoMSIA. Those parameters indicated the model was reliable and predictable. We designed several molecules with high activities according to the contour maps produced by the CoMFA and CoMSIA models.     

Molecular Modeling Studies of Vascular Endothelial Growth Factor Receptor Tyrosine Kinase Inhibitors Combining Molecular Docking and 3D-QSAR Methods

The vascular endothelial growth factor (VEGF) and its receptor tyrosine kinases VEGFR-2 or kinase insertdomain receptor (KDR) have emerged as attractive targets for the design of novel anticancer agents. In the present work, molecular docking method combined with three dimensional quantitative structure-activity relationships (comparative molecular field analysis (CoMFA) and comparative molecular similarity indice analysis (CoMSIA)) to analyze the possible interactions between KDR and those derivatives which acted as selective inhibitors. The CoMFA and CoMSIA models gave a cross-validated coefficient Q2 of 0.713 and 0.549, non-cross-validated R2 values of 0.974 and 0.878, and predicted R2 values of 0.966 and 0.823, respectively. The 3D contour maps generated by the CoMFA and CoMSIA models were used to identify the key structural requirements responsible for the biological activity. The information obtained from 3D-QSAR and docking studies were very helpful to design novel selective inhibitors of KDR with desired activity and good chemical property.     

苯并咪唑类缓蚀剂的3D-QSAR研究及分子设计

采用比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 对苯并咪唑衍生物抗盐酸腐蚀的缓蚀性能进行了三维定量构效关系研究, 并使用留一法交叉验证手段对3D-QSAR模型的稳定性及预测能力进行了分析. 结果表明, 立体场、静电场和氢键供体场(电子给体)是影响苯并咪唑缓蚀剂缓蚀性能的主要因素; 所构建的CoMFA模型(q2=0.541, R2=0.996)和CoMSIA模型(q2=0.581, R2=0.987)均具有较好的统计学稳定性和预测能力. 基于3D-QSAR等势图设计出了几种具有较好缓蚀性能的苯并咪唑化合物, 为油气田新型缓蚀剂的研发提供了一种新思路.     

Development of novel 3D-QSAR combination approach for screening and optimizing B-Raf inhibitors in silico

B-Raf is a member of the RAF family of serine/threonine kinases: it mediates cell division, differentiation, and apoptosis signals through the RAS-RAF-MAPK pathway. Thus, B-Raf is of keen interest in cancer therapy, such as melanoma. In this study, we propose the first combination approach to integrate the pharmacophore (PhModel), CoMFA, and CoMSIA models for B-Raf, and this approach could be used for screening and optimizing potential B-Raf inhibitors in silico. Ten PhModels were generated based on the HypoGen BEST algorithm with the flexible fit method and diverse inhibitor structures. Each PhModel was designated to the alignment rule and screening interface for CoMFA and CoMSIA models. Therefore, CoMFA and CoMSIA models could align and recognize diverse inhibitor structures. We used two quality validation methods to test the predication accuracy of these combination models. In the previously proposed combination approaches, they have a common factor in that the number of training set inhibitors is greater than that of testing set inhibitors. In our study, the 189 known diverse series B-Raf inhibitors, which are 7-fold the number of training set inhibitors, were used as a testing set in the partial least-squares validation. The best validation results were made by the CoMFA09 and CoMSIA09 models based on the Hypo09 alignment model. The predictive r(2)(pred) values of 0.56 and 0.56 were derived from the CoMFA09 and CoMSIA09 models, respectively. The CoMFA09 and CoMSIA09 models also had a satisfied predication accuracy of 77.78% and 80%, and the goodness of hit test score of 0.675 and 0.699, respectively. These results indicate that our combination approach could effectively identify diverse B-Raf inhibitors and predict the activity.     

3D-QSAR and Docking Studies of Pyrido[2,3-d]pyrimidine Derivatives as Wee1 Inhibitors

In order to investigate the inhibiting mechanism and obtain some helpful information for de-signing functional inhibitors against Wee1, three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies have been performed on 45 pyrido[2,3-d] pyrim-idine derivatives acting as Wee1 inhibitors. Two optimal 3D-QSAR models with significant statistical quality and satisfactory predictive ability were established, including the CoMFA model (q²=0.707, R²=0.964) and CoMSIA model (q²=0.645, R²=0.972). The external val-idation indicated that both CoMFA and CoMSIA models were quite robust and had high predictive power with the predictive correlation coefficient values of 0.707 and 0.794, essen-tial parameter r²_m values of 0.792 and 0.826, the leave-one-out r²_m(LOO) values of 0.781 and 0.809, r²_m(overall) values of 0.787 and 0.810, respectively. Moreover, the appropriate binding orientations and conformations of these compounds interacting with Wee1 were revealed by the docking studies. Based on the CoMFA and CoMSIA contour maps and docking analyses, several key structural requirements of these compounds responsible for inhibitory activity were identified as follows: simultaneously introducing high electropositive groups to the sub-stituents R1 and R5 may increase the activity, the substituent R2 should be smaller bulky and higher electronegative, moderate-size and strong electron-withdrawing groups for the substituent R3 is advantageous to the activity, but the substituent X should be medium-size and hydrophilic. These theoretical results help to understand the action mechanism and design novel potential Wee1 inhibitors.     

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

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

Homology model-guided 3D-QSAR studies of HIV-1 integrase inhibitors

In the present study, we report the exploration of binding modes of potent HIV-1 integrase (IN) inhibitors MK-0518 (raltegravir) and GS-9137 (elvitegravir) as well as chalcone and related amide IN inhibitors we recently synthesized and the development of 3D-QSAR models for integrase inhibition. Homology models of DNA-bound HIV-1 IN were constructed on the basis of the X-ray crystal structure of the foamy virus IN-DNA complex (PDB ID: 3L2T ) and used for docking. The binding modes of raltegravir and elvitegravir in our homology models are in accordance with those in the foamy virus structure revealing interactions important for inhibitor-IN binding. To gain further insights into the structural requirements for IN inhibition, three-dimensional quantitative structure activity relationship (3D-QSAR) studies were conducted using raltegravir, elvitegravir, and their analogs; our synthesized 3-keto salicylic acid IN inhibitor series; as well as other structurally related HIV-1 IN inhibitors. In the first part of the study with 103 compounds, atom-fit alignments, I and II, and docking-based alignment, III, were used to develop 3D-QSAR models 1, 2, and 3, respectively, each comprising comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSARs. This initial analysis indicated that the docking-based (structure-based) model 3 performed better than the atom-fit (ligand-based) models 1 and 2, in terms of statistical significance and robustness. Thus, the docking-based alignment was then subsequently used with an expanded data set of 296 compounds for building a more comprehensive 3D-QSAR, model 4. Model 4 afforded good q2 values of 0.70 and 0.75 for CoMFA and CoMSIA 3D-QSARs, respectively, and showed good predictive performance on an external validation test set of 59 compounds with predictive r2 values up to 0.71. The HIV IN-DNA homology model of biological relevance and the comprehensive 3D-QSAR models developed in the present study provide insights and new predictive tools for structure-based design and optimization of IN inhibitors.