3-氨磺酰苯甲酸类AKR1C3抑制剂的3D-QSAR和分子对接研究

醛酮还原酶1C3(AKR1C3)作为治疗前列腺癌的新靶点已成为研究热点,3-氨磺酰苯甲酸衍生物对其具有高效的选择性和抑制活性。本文采用比较分子场分析(COMFA)和比较分子相似性指数分析(COMSIA)方法,将经分子对接后的34个优势构象组成训练集和11个优势构象组成测试集,构建三维定量构效关系(3D-QSAR)模型。COMFA模型的交叉验证系数(q2),非交叉验证系数(R2),标准偏差(SEE)和F值分别为0.761,0.973,0.122,185.963;自举法回归系数为R2bs=0.98。最佳组合COMSIA模型的q2,R2,SEE,F和R2bs分别为0.734,0.984,0.097,147.850,0.994。COMFA和COMSIA模型的系统外部测试R2pred分别为0.864和0.756,r2m分别为0.8127和0.5377。这些结果表明,所建立的QSAR模型具有较高的可靠性和较强预测能力。经三维等势图分析可知,在2、5或6位适当增加取代基体积,或在5位引入氢键受体,或在7位引入负电性取代基则能提高化合物的生物活性。该模型为进一步设计具有更优选择性和活性的化合物提供了理论依据。     

雌激素β受体喹啉类配体的分子对接及3D-QSAR研究

对33个喹啉衍生物的雌激素β受体活性进行了分子对接以及比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA). 对接结果显示氢键和疏水作用是配体与受体结合的主要因素,同时结果亦显示对接结合能与观测值pIC50具有极显著的线性相关性. 根据对接后各优势构象将33个样本进行叠合并进行CoMFA与CoMSIA研究,均得到了较优的结果,其中以选用立体场、静电场和疏水场建立的CoMSIA模型结果最优,其主成分数,r2,q2（LOO）和r2pred分别为2, 0.894, 0.708和0.802. 构效关系模型分析显示基团的空间位阻、电性及疏水作用是影响活性的主要因素     

酪氨酸蛋白磷酸酯酶1B抑制剂的分子对接和三维定量构效关系研究

用一种柔性分子对接方法(FlexX)将12个2-草酰胺苯甲酸类抑制剂和酪氨酸蛋白磷酸酯酶(PTP1B)活性口袋进行分子对接,对接程序预测的抑制剂和酶之间的相互作用能与抑制活性之间有很好的相关性(非线性相关系数R²达0.859),这说明对接结果可以比较准确地预测抑制剂和PTP1B之间的结合模式.然后,将33个同类抑制剂的骨架叠合在分子对接预测的结合构象上,用比较分子力场分析方法(CoMFA)对其进行三维定量活性构效关系研究,得到的CoMFA模型具有很好的统计相关性(交互验证回归系数q²为0.650),并可以准确地预测测试集6个化合物的活性(平均标准偏差为0.177).同时,由CoMFA模型得出的抑制剂改造信息与用FlexX预测的结合模式是一致的,进一步证明我们预测的结合模式是正确的.为研究这类抑制剂和PTP1B的结合模式及对抑制剂进行结构改造提供了信息.     

二芳基三嗪类HIV-1逆转录酶抑制剂的分子对接及3D-QSAR研究

对40个二芳基三嗪类HIV-1逆转录酶抑制剂进行了分子对接研究,结果表明,2个疏水性芳香取代基团与结合口袋底部形成的疏水和范德华相互作用、三嗪环母核及其R4取代基与结合位点产生的氢键和静电作用以及R4取代基与袋口形成的空间位阻效应是影响该系列化合物活性的重要因素.根据对接优势构象进行分子叠合和比较分子相似性指数分析(C...     

HIV-1整合酶链转移抑制剂的3D-QSAR、分子对接和分子动力学模拟研究

为了获得高活性、结构新颖的整合酶链转移(INST)抑制剂,本文采用Co MFA和Co MSIA两种方法对32个萘啶类INST抑制剂进行了三维定量构效关系研究,并建立了相关模型,其交叉验证系数分别为q~2=0. 809和q~2=0. 816,拟合验证系数分别为r~2=0. 998和r~2=0. 981,表明所建立的模型是可靠的且具有一定的预测能力。利用分子对接探讨小分子化合物与INST蛋白的相互作用模式,结果表明,萘啶类化合物主要通过疏水作用和氢键作用与INSTIs蛋白结合。最后通过分子动力学模拟进一步验证对接结果发现,对接的结合模式与分子动力学模拟得到的结果是一致的。本研究获得的综合模型和推论可以为开发有效的HIV INSTIs提供重要的理论信息。     

脒类KARI酶抑制剂的分子对接和3D-QSAR研究

对设计合成的20个单脒类化合物的水稻KARI酶体外抑制活性和体内除草活性进行了分子对接和三维定量构效关系研究. 前者采用AutoDock3.0方法, 研究发现化合物活性变化趋势与分子对接计算结果基本一致, 通过分析化合物9与KARI酶活性氨基酸残基结合模式发现, 残基Glu319, Asp315, Glu496, Gly253, Met254, Cys517等对氢键和静电相互作用以及疏水作用都有重要贡献; 后者研究采用比较分子力场分析(CoMFA)方法, 结果表明立体场和静电场对活性的贡献分别为67.8%和32.2%, 交叉验证系数r_cv²＝0.774, 非交叉验证r²＝0.999, F＝1593.134, 标准偏差s＝0.036, 所建立的3D-QSAR模型对化合物除草活性具有较好的预测能力. 两种方法研究结果为进一步设计合成更高活性的KARI酶抑制剂提供了指导.     

趋化因子CCR2的同源模建,分子对接和3D-QSAR研究

趋化因子CCR2参与炎症反应、免疫移植排斥和肿瘤的发生,已成为新的研究热点。本文以CCR5的晶体结构为模板,同源模建CCR2的结构,并用CCR2小分子抑制剂与其进行分子对接以得到小分子的最优构象。在对接叠合的基础上建立了QSAR模型,采用比较分子场分析(Co MFA)以及比较分子相似性分析(Co MSIA)研究得到Co MFA和Co MSIA模型最佳评价参数分别为q2=0.743,r2=0.968和q2=0.68,r2=0.978。3D-QSAR模型的等势图分析表明,改造配体R3基团可提高化合物活性。所建模型稳定性好、预测性强,对基于CCR2的小分子抑制剂的设计、优化和改造提供了参考。     

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

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

Design of Novel IRAK4 Inhibitors Using Molecular Docking,Dynamics Simulation and 3D-QSAR Studies

Treatment of several autoimmune diseases and types of cancer has been an intense area of research over the past two decades. Many signaling pathways that regulate innate and/or adaptive immunity, as well as those that induce overexpression or mutation of protein kinases, have been targeted for drug discovery. One of the serine/threonine kinases, Interleukin-1 Receptor Associated Kinase 4 (IRAK4) regulates signaling through various Toll-like receptors (TLRs) and interleukin-1 receptor (IL1R). It controls diverse cellular processes including inflammation, apoptosis, and cellular differentiation. MyD88 gain-of-function mutations or overexpression of IRAK4 has been implicated in various types of malignancies such as Waldenström macroglobulinemia, B cell lymphoma, colorectal cancer, pancreatic ductal adenocarcinoma, breast cancer, etc. Moreover, over activation of IRAK4 is also associated with several autoimmune diseases. The significant role of IRAK4 makes it an interesting target for the discovery and development of potent small molecule inhibitors. A few potent IRAK4 inhibitors such as PF-06650833, RA9 and BAY1834845 have recently entered phase I/II clinical trial studies. Nevertheless, there is still a need of selective inhibitors for the treatment of cancer and various autoimmune diseases. A great need for the same intrigued us to perform molecular modeling studies on 4,6-diaminonicotinamide derivatives as IRAK4 inhibitors. We performed molecular docking and dynamics simulation of 50 ns for one of the most active compounds of the dataset. We also carried out MM-PBSA binding free energy calculation to identify the active site residues, interactions of which are contributing to the total binding energy. The final 50 ns conformation of the most active compound was selected to perform dataset alignment in a 3D-QSAR study. Generated RF-CoMFA (q² = 0.751, ONC = 4, r² = 0.911) model revealed reasonable statistical results. Overall results of molecular dynamics simulation, MM-PBSA binding free energy calculation and RF-CoMFA model revealed important active site residues of IRAK4 and necessary structural properties of ligand to design more potent IRAK4 inhibitors. We designed few IRAK4 inhibitors based on these results, which possessed higher activity (predicted pIC₅₀) than the most active compounds of the dataset selected for this study. Moreover, ADMET properties of these inhibitors revealed promising results and need to be validated using experimental studies.     

HIV-1逆转录酶抑制剂的3D-QSAR研究及分子设计

采用Topomer Co MFA方法对24个二芳基苯胺衍生物进行三维定量构效关系研究,建立了3DQSAR模型,所得优化模型的非交叉相关系数、交互验证系数以及外部验证的复相关系数分别为0.928,0.654和0.940,结果表明该模型具有良好的稳定性和预测能力。采用分子对接技术对药物与受体的作用机制进行了研究,结果显示,药物与HIV-1逆转录酶的LYS172,GLU138,LYS101等位点作用明显。运用这些信息进行分子设计,在理论上获得了一些具有较高活性的新的二芳基苯胺类抗艾滋病药物,该QSAR的研究结果可为新药合成提供理论参考。     

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.     

芳香噻嗪类衍生物作为选择性醛糖还原酶抑制剂的分子对接和基于受体的三维定量构效关系研究

芳香噻嗪类衍生物被证明是一类选择性较好的高活性醛糖还原酶抑制剂(ARIs).本文对44个芳香噻嗪类化合物进行了分子对接(docking)和三维定量构效关系(3D-QSAR)研究,并探索了此类化合物与醛糖还原酶(ALr²)的作用机理.醛糖还原酶与醛还原酶(ALR1)活性位点的叠加结果显示, ALr²中残基Leu 300和Cys298的存在是化合物1m具有高选择性的原因.分别建立了比较分子场分析方法(CoMFA, q² = 0.649, r² =0.934; q²:交叉验证相关系数, r²:非交叉验证相关系数)和比较分子相似性指数分析方法(CoMSIA, q² = 0.746, r² = 0.971)模型,并对影响此类化合物生物活性的结构进行了鉴定.结果显示,两个模型均具有较高预测能力,并通过测试集中的7个化合物进行了验证,其中CoMFA模型和CoMSIA模型的预测相关系数(r_Pred²)分别为0.748和0.828. 3D-QSAR模型中的三维等值线图表明,在化合物1m的苄基环上C3和C4位置以及苯并噻嗪母核上C5和C7位置进行改进可能对生物活性的提高有利,此预测与我们前期报道的苯并噻嗪母核C7位改进结果一致.本文所建3D-QSAR模型能够在理性设计具有更高生物活性的新型ARIs中发挥重要作用.     

DCN1-UBE2M相互作用抑制剂的3D-QSAR、分子对接和分子动力学研究

类泛素化是一种蛋白质翻译后修饰,其异常会导致神经退行性疾病和多种肿瘤的发生,因此它被视为有希望的抗肿瘤靶标.研究表明,抑制DCN1-UBE2M相互作用可选择性阻遏类泛素化.本文基于哌啶基脲类DCN1-UBE2M相互作用抑制剂进行3D-QSAR、分子对接和分子动力学模拟研究.利用3D-QSAR中的CoMFA和CoMSIA...     

α-咔啉类GSK-3β抑制剂的合成、活性评价和分子对接研究

糖原合成酶激酶3β(GSK-3β)是神经退化性疾病如早老性痴呆、II型糖尿病、癌症等多种重复杂性疾病的药物靶标,α-咔啉类化合物具有多种生物活性.以GSK-3β为作用靶点,利用Graebe-Ullmann反应合成了α-咔啉,并采用Buchwald-Hartwig偶联反应合成了9个4位取代的α-咔啉胺类衍生物.体外的GSK-3β激酶酶抑制活性测试表明其中4个化合物具有一定的抑制活性,其中化合物3c的IC50值达到了5.1μmol L-1.同时采用了分子对接方法分别研究了化合物2,3c与大分子蛋白的作用模式,初步探讨了影响抑制活性的原因,为进一步的研究提供了依据.     

Design Two Novel Tetrahydroquinoline Derivatives against Anticancer Target LSD1 with 3D-QSAR Model and Molecular Simulation

Lysine-specific demethylase 1 (LSD1) is a histone-modifying enzyme, which is a significant target for anticancer drug research. In this work, 40 reported tetrahydroquinoline-derivative inhibitors targeting LSD1 were studied to establish the three-dimensional quantitative structure–activity relationship (3D-QSAR). The established models CoMFA (Comparative Molecular Field Analysis (q² = 0.778, Rpred2 = 0.709)) and CoMSIA (Comparative Molecular Similarity Index Analysis (q² = 0.764, Rpred2 = 0.713)) yielded good statistical and predictive properties. Based on the corresponding contour maps, seven novel tetrahydroquinoline derivatives were designed. For more information, three of the compounds (D1, D4, and Z17) and the template molecule 18x were explored with molecular dynamics simulations, binding free energy calculations by MM/PBSA method as well as the ADME (absorption, distribution, metabolism, and excretion) prediction. The results suggested that D1, D4, and Z17 performed better than template molecule 18x due to the introduction of the amino and hydrophobic groups, especially for the D1 and D4, which will provide guidance for the design of LSD1 inhibitors.     

3C-like蛋白酶抑制剂的构效关系、分子对接和分子动力学

3C-like蛋白酶是中东呼吸综合征冠状病毒（MERS-CoV）等其它冠状病毒的繁殖过程中极为重要的蛋白酶。它已成为人类在抗冠状病毒领域中的研究热点。本文基于计算生物学方法对与MERS-CoV同属的蝙蝠冠状病毒HKU4（HKU4-CoV）的43个肽类3C-like蛋白酶抑制剂分子,建立三维定量构效关系（3D-QSAR）模型。在基于配体叠合的基础上,发现比较分子相似性指数分析法（CoMSIA）中的四个场组合（位阻场、静电场、氢键供体场与氢键受体场）为最优的模型（Q²=0.522,R_ncv²=0.996,R_pre²=0.904;Q²：交叉验证相关系数,R_ncv²：非交叉验证相关系数,R_pre²：验证集分子的预测值相关系数）,并借助该模型通过分子对接（docking）与分子动力学（MD）方法阐明了配受体结合作用。实验结果表明：（1）基于最优的CoMSIA模型基础上的三维等势图形象地说明了分子基团的位阻作用、静电作用、氢键供体与氢键受体作用对分子生物活性的影响;（2）分子对接研究结果显示了疏水性以及结晶水、氨基酸His166和Glu169在配体和受体结合过程中产生重要作用;（3）分子动力学模拟进一步验证了分子对接模型的可靠性,并发现了两个新的关键氨基酸Ser24与Gln192,它们与配体产生了两个较强的氢键。此外,根据这些结果,一些新的具有潜在抑制活性的肽类化合物作为3C-like蛋白酶抑制剂被获得。以上结果能够帮助深入了解3C-like蛋白酶与肽类抑制剂的作用机理,并且能够为今后的抗MERS-CoV药物设计提供有价值的参考。     

3D-QSAR和分子对接研究吲哚咔唑类细胞周期蛋白激酶抑制剂的选择性

细胞周期蛋白激酶(cyclin-dependent kinases, CDKs)是近年来治疗肿瘤的重要靶标. 由于大多数激酶ATP结合位点的保守性, CDK选择性激酶抑制剂的开发成为当前的研发难点和热点. 针对吲哚咔唑类CDK抑制剂, 我们采用比较分子力场分析方法(CoMFA)建立了CDK2-QSAR(quantitative structure-activity relationship)和CDK4-QSSR(quantitative structure-selectivity relationship)模型. 所建模型的交叉验证系数q2分别为0.722和0.703; 非交叉验证系数r2分别为0.977和0.946, 表明其具有较好的预测能力. 同时, 用分子对接的方法分析了这类化合物与CDK4同源模建结构的作用模式, 根据这两个模型发现, 吲哚咔唑类化合物的R5和R6位长链取代对CDK4的选择性具有一定的影响, 而且结合其作用模式比较合理地解释了这类抑制剂的选择性原因, 这对CDKs的选择性研究具有一定的指导意义.     

磺胺基羟肟酸类HDAC抑制剂三维定量构效关系

组蛋白去乙酰化酶(HDAC)对染色质分布和基因调节起着重要的作用,也是治疗癌症和其它疾病的新靶点.羟肟酸类抑制剂是目前研究最多的组蛋白去乙酰化酶抑制剂.应用比较分子力场(CoMFA)法对一系列磺胺基羟肟酸类HDAC抑制剂进行了结构活性关系研究,得到的模型具有较高的交叉验证系数(q²=0.704).并在此基础上,建立了非交叉验证的偏最小二乘分析(PLS)模型.用该模型对随机选择的6个化合物组成的测试集进行了预测,得到了令人满意的结果,所建模型具有良好的预测能力.本研究对于设计高活性的HDAC抑制剂及抗癌药物都有指导意义.     

QSAR,Molecular Docking,MD Simulation and MMGBSA Calculations Approaches to Recognize Concealed Pharmacophoric Features Requisite for the Optimization of ALK Tyrosine Kinase Inhibitors as Anticancer Leads

ALK tyrosine kinase ALK TK is an important target in the development of anticancer drugs. In the present work, we have performed a QSAR analysis on a dataset of 224 molecules in order to quickly predict anticancer activity on query compounds. Double cross validation assigns an upward plunge to the genetic algorithm–multi linear regression (GA-MLR) based on robust univariate and multivariate QSAR models with high statistical performance reflected in various parameters like, fitting parameters; R² = 0.69–0.87, F = 403.46–292.11, etc., internal validation parameters; Q²_LOO = 0.69–0.86, Q²_LMO = 0.69–0.86, CCC_cv = 0.82–0.93, etc., or external validation parameters Q²_F1 = 0.64–0.82, Q²_F2 = 0.63–0.82, Q²_F3 = 0.65–0.81, R²_ext = 0.65–0.83 including RMSE_tr < RMSE_cv. The present QSAR evaluation successfully identified certain distinct structural features responsible for ALK TK inhibitory potency, such as planar Nitrogen within four bonds from the Nitrogen atom, Fluorine atom within five bonds beside the non-ring Oxygen atom, lipophilic atoms within two bonds from the ring Carbon atoms. Molecular docking, MD simulation, and MMGBSA computation results are in consensus with and complementary to the QSAR evaluations. As a result, the current study assists medicinal chemists in prioritizing compounds for experimental detection of anticancer activity, as well as their optimization towards more potent ALK tyrosine kinase inhibitor.     

丙型肝炎病毒抑制剂的三维药效团和构效关系

通过ＣＡＴＡＬＹＳＴ软件包得到了两类ＨＣＶ ＮＳ３丝氨酸蛋白酶抑制剂的三维药效团模型。尽管这两类抑制剂具有完全不同的骨架结构,但得到的药效团却具有共同的特性。这当这两类抑制剂和受体发生相互作用时,可能采用了相似的结合模式。根据药效团模型,进行了三维构效关系的研究。结果表明,得到的药效团模型具有良好的预测能力（线性回归系数Ｒ＝０．８９）。