螺环吲哚类MDM2抑制剂的分子对接、定量构效关系和分子动力学模拟

采用分子对接、三维定量构效关系(3D-QSAR)和分子动力学方法研究了21个螺环吲哚类化合物与MDM2蛋白的相互作用, 并建立了相关预测模型. 比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)模型的交互验证相关系数q²分别为0.573 和0.651, 非交互验证相关系数r²分别为0.948和0.980. 分子对接得到的结合模式与分子动力学模拟得到的结果一致, 结合模式表明该类螺环吲哚化合物主要通过疏水相互作用和氢键与MDM2结合. 基于上述相互作用模型设计并合成了6个新结构螺环吲哚化合物, 并在MDM2高表达的前列腺癌LNCaP细胞株上测定其活性, 结果表明化合物5, 6的半数抑制浓度均低于1μg·mL^-1, 可作为新的抗肿瘤药物先导化合物进一步深入研究. 本研究对以MDM2为靶点的新结构螺环吲哚类抑制剂的开发提供了理论和实验依据.     

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

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

4,5,6-三取代嘧啶苯磺酰脲类化合物的生物活性、分子对接与3D-QSAR关系研究

用柔性分子对接方法(FlexX)将15个4,5,6-三取代嘧啶苯磺酰脲化合物以及3个不含5-位取代嘧啶苯磺酰脲化合物(分别为4,6-双取代嘧啶和4-取代嘧啶)和乙酰羟酸合成酶(AHAS)活性口袋进行了对接, 对接程序预测的抑制剂和酶之间的相互作用能与抑制活性之间有一定的相关性, 相关系数为0.660. 然后采用比较分子相似性指数分析(CoMSIA)对27个新型4,5,6-三取代嘧啶苯磺酰脲类化合物的除草活性进行三维定量构效关系(3D-QSAR)研究. 建立了三维定量构效关系CoMSIA模型, 立体场、静电场和氢键的贡献分别为47.3%, 32.8%, 19.9%. 交叉验证系数q2值为0.520. 根据CoMSIA模型的立体场、静电场、氢键给体场三维等值线图不仅直观地解释了结构与活性的关系, 并且与用FlexX预测的结合模式相一致, 证明了我们预测的结合模式是可靠的, 为进一步设计高活性的标题化合物提供较好的理论指导.     

肿瘤血管阻断剂黄酮-8-乙酸类衍生物的三维定量关系研究

选择了肿瘤血管阻断剂黄酮-8-乙酸类衍生物, 采用比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)方法进行三维定量构效关系的研究. 33个化合物建立了预测模型, 6个化合物作为训练集进行模型验证. 其中CoMFA模型的交叉验证系数q²＝0.621, 最佳主成分数为4, 标准偏差s_press＝0.345, 非交叉相关系数r²＝0.945, 标准偏差s＝0.131, F＝120.455. CoMSIA模型的交叉验证系数q²＝0.700, 最佳主成分数为5, 标准偏差s_press＝0.312, 非交叉相关系数r²＝0.946, 标准偏差s＝0.133, F＝94.193. 计算结果表明, 构建的CoMFA和CoMSIA模型具有良好的预测能力, 可用于指导该类化合物的设计.     

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

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

橙酮类DRAK2抑制剂的3D-QSAR及分子对接研究

本文对橙酮类DRAK2抑制剂的化学结构与生物活性之间的关系进行研究。采用三维定量构效关系（3D-QSAR）中的比较分子力场分析（CoMFA）和比较分子相似性指数分析（CoMSIA）方法针对59个DRAK2抑制剂建立3D-QSAR模型,阐明了抑制剂化学结构与其生物活性之间的关系。所构建的CoMFA模型交叉验证系数（q²）为0.625,非交叉验证系数（r²）为0.811,标准偏差（S）为0.365,Fisher检验值（F）为59.971;所构建的CoMSIA模型q²为0.62,r²为0.846,S为0.333,F值为56.453。内部和外部验证参数表明,生成的3D-QSAR模型均具有良好的预测能力和显著的统计学可靠性。分子对接实验与等势图的一致性,进一步表明本次分子模拟是可靠的。本研究对发现新型的潜在的更高活性的橙酮类DRAK2抑制剂具有指导意义。     

喹啉酮类BRD4抑制剂的3D-QSAR研究

使用比较分子力场分析法(CoMFA)和比较分子相似性指数法(Co MSIA)对33个已报道的喹啉酮类BRD4抑制剂进行3D-QSAR模型建立,研究了其化学结构和生物活性间的关系,并用计算机辅助药物设计(computer-aided drug design,CADD)设计出7个喹啉酮类抑制剂。结果表明,建立的CoMFA(q²=0.926,r²=0.997,r²_pred=0.744)和Co MSIA(q²=0.939,r²=0.991,r²_pred=0.786)模型具有较好的预测能力,基于这些模型设计的7个新喹啉酮类BRD4抑制剂具有高活性,并对其进行ADMET性质评价和类药性分析。以上研究结果有助于改造和开发更加有效的喹啉酮类BRD4抑制剂。     

新型二氟甲基磷酸类酪氨酸蛋白磷酸酯酶1B抑制剂的分子动力学模拟和结合自由能计算

通过分子对接建立了一系列含二氟甲基磷酸基团(DFMP)或二氟甲基硫酸基团(DFMS)的抑制剂与酪氨酸蛋白磷酸酯酶1B(PTP1B)的相互作用模式, 并通过1 ns的分子动力学模拟和molecular mechanics/generalized Born surface area (MM/GBSA)方法计算了其结合自由能. 计算获得的结合自由能排序和抑制剂与靶酶间结合能力排序一致; 通过基于主方程的自由能计算方法, 获得了抑制剂与靶酶残基间相互作用的信息, 这些信息显示DFMP/DFMS基团的负电荷中心与PTP1B的221位精氨酸正电荷中心之间的静电相互作用强弱决定了此类抑制剂的活性, 进一步的分析还显示位于DFMP/DFMS基团中的氟原子或其他具有适当原子半径的氢键供体原子会增进此类抑制剂与PTP1B活性位点的结合能力.     

HIV-1整合酶与抑制剂金精三羧酸复合物的分子动力学模拟

用分子对接程序(Autodock)将含有一个Mg²⁺的HIV-1整合酶核心区(以下简称IN-A)与抑制剂小分子金精三羧酸(简称Aurin)进行对接,预测其未知的复合物结构,然后用分子动力学(MD)方法对IN-A与Aurin的对接结果进行了950 ps的模拟.MD模拟结果发现,IN-A与Aurin形成了两个稳定的氢键,Mg²⁺也与Aurin上的氧原子形成了稳定的配键,IN-A与Aurin之间的静电相互作用能和范德华相互作用能的平均值分别为-205.8和-162.7 kJ/mol.根据MD模拟得到的IN-A与Aurin相互作用后的构象变化信息,我们对对接复合物结构进行了修正,给出了更加合理和稳定的复合物预测结构.本工作得到的HIV-1整合酶与抑制剂Aurin的结合模式信息将有助于设计和改造出效果更好的抗HIV-1整合酶的先导化合物.     

6-氮杂吲唑类Pim-1激酶抑制剂的3D-QSAR、药效团模型研究和分子设计

Pim-1激酶通过作用于多种信号通路或靶点影响肿瘤的发生发展,近年来被认为是肿瘤治疗的良好靶标。本文采用SYBYL-X2.1.1软件中的Topomer CoMFA、GALAHAD模块建立计算机模型,研究39个基于6-氮杂吲唑环的Pim-1激酶抑制剂的三维定量构效关系及药效团特征元素。结果显示,Topomer CoMFA建模所得交叉验证系数(q²)和相关系数(r²)分别为0.756和0.951,结合外部验证表明此3D-QSAR模型具有较高预测能力及较好的统计学稳定性,同时,用等势图描述了R1、R2基团处立体场、静电场对活性的具体影响。药效团研究结果表明,含氢键受体的芳香杂环母核结构,以及侧链取代基中含有芳香杂环结构对化合物的活性贡献较大。最后根据上述模型信息新设计了15个Pim-1激酶抑制剂分子并完成活性预测及分子对接模式研究,其中4个分子的预测pIC₅₀高于建模分子中活性最好的化合物17,Surflex-Dock分析显示新设计分子均与Pim-1激酶形成较强氢键相互作用。基于6-氮杂吲唑环的Pim-1激酶抑制剂的3D-QS...     

CoMFA 3D-QSAR analysis of HIV-1 RT nonnucleoside inhibitors, TIBO derivatives based on docking conformation and alignment

HIV-1 RT is one of the key enzymes in the duplication of HIV-1. Inhibitors of HIV-1 RT are classified as nonnucleoside RT inhibitors (NNRTIs) and nucleoside analogues. NNRTIs bind in a region not associated with the active site of the enzyme. Within the NNRTI category, there is a set of inhibitors commonly referred to as TIBO inhibitors. Fifty TIBO inhibitors were used in the work to build 3-D QSAR models. The two known crystal structures of complexes are used to investigate and validate the docking protocol. The results show that the docking simulations reproduce the crystal complexes very well with RMSDs of approximately 1 A and approximately 0.6 A for 1REV and 1COU, respectively. The alignment of molecules and "active" conformation selection are the key to a successful 3D-QSAR model by CoMFA. The flexible docking (Autodock3) was used on determination of "active" conformation and molecular alignment, and CoMFA and CoMSIA were used to develop 3D-QSAR models of 50 TIBOs in the work. The 3D-QSAR models demonstrate a good ability to predict the activity of studied compounds (r2 = 0.972, 0.944, q2 = 0.704, 0.776). It is shown that the steric and electrostatic properties predicted by CoMFA contours can be related to the binding structure of the complex. The results demonstrate that the combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models.     

3D-QSAR and docking studies on 2-arylbenzoxazole and linker-Y transthyretin amyloidogenesis inhibitors

Transthyretin (TTR), a plasma protein with a tetramer structure, could form amyloid fibril associated with several human diseases through the dissociation of tetramer and the misfolding of monomer. These amyloidogenesis can be inhibited by small molecules which bind to the central channel of TTR. A number of small molecules like 2-arylbenzoxazoles (ABZ) analogues are proposed as promising therapeutic strategy to treat amyloidosis. In this work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies were performed on series of 2-arylbenzoxazoles (ABZ) and linker-Y analogues to investigate the inhibitory activities of TTR amyloidogenesis at atomic level. Significant correlation coefficients for ABZ series (CoMFA, r ² = 0.877, q ² = 0.431; CoMSIA, r ² = 0.836, q ² = 0.447) and those for linker-Y series (CoMFA, r ² = 0.828, q ² = 0.522; CoMSIA, r ² = 0.800, q ² = 0.493) were obtained, and the generated models were validated using test sets. In addition, docking studies on 6 compounds binding to TTR were performed to analyze the forward or reverse binding mode and interactions between molecules and TTR. These results from 3D-QSAR and docking studies have great significance for designing novel TTR amyloidogenesis inhibitors in the future.     

Docking and 3D-QSAR Analysis on a Series of Pyridone-based EZH2 Inhibitors

Enhancer of Zeste homolog 2(EZH2) is closely correlated with malignant tumor and regarded as a promising target to treat B-cell lymphoma. In our research, the molecular docking and three-dimensional quantitative structure-activity relationships(3D-QSAR) studies were performed on a series of pyridone-based EZH2 compounds. Molecular docking allowed us to study the critical interactions at the binding site of EZH2 protein with inhibitors and identify the practical conformations of ligands in binding pocket. Moreover, the docking-based alignment was applied to derive the reliable 3D-QSAR models. Comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA) provided available ability of visualization. All the derived 3D-QSAR models were considered to be statistically significant with respect to the internal and external validation parameters. For the CoMFA model, q~2 = 0.649, r~2 = 0.961 and r~2 pred = 0.877. For the CoMSIA model, q~2 = 0.733, r~2 = 0.980 and r~2 pred = 0.848. With the above arguments, we extracted the correlation between the biological activity and structure. Based on the binding interaction and 3D contour maps, several new potential inhibitors with higher biological activity predicted were designed, which still awaited experimental validation. These theoretical conclusions could be helpful for further research and exploring potential EZH2 inhibitors.     

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies

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.     

3D QSAR studies on protein tyrosine phosphatase 1B inhibitors: comparison of the quality and predictivity among 3D QSAR models obtained from different conformer-based alignments

A set of 65 flexible peptidomimetic competitive inhibitors (52 in the training set and 13 in the test set) of protein tyrosine phosphatase 1B (PTP1B) has been used to compare the quality and predictive power of 3D quantitative structure-activity relationship (QSAR) comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for the three most commonly used conformer-based alignments, namely, cocrystallized conformer-based alignment (CCBA), docked conformer-based alignment (DCBA), and global minima energy conformer-based alignment (GMCBA). These three conformers of 5-[(2S)-2-({(2S)-2-[(tert-butoxycarbonyl)amino]-3-phenylpropanoyl}amino)3-oxo-3-pentylamino)propyl]-2-(carboxymethoxy)benzoic acid (compound number 66) were obtained from the X-ray structure of its cocrystallized complex with PTP1B (PDB ID: 1JF7), its docking studies, and its global minima by simulated annealing. Among the 3D QSAR models developed using the above three alignments, the CCBA provided the optimal predictive CoMFA model for the training set with cross-validated r2 (q2)=0.708, non-cross-validated r2=0.902, standard error of estimate (s)=0.165, and F=202.553 and the optimal CoMSIA model with q2=0.440, r2=0.799, s=0.192, and F=117.782. These models also showed the best test set prediction for the 13 compounds with predictive r2 values of 0.706 and 0.683, respectively. Though the QSAR models derived using the other two alignments also produced statistically acceptable models in the order DCBA>GMCBA in terms of the values of q2, r2, and predictive r2, they were inferior to the corresponding models derived using CCBA. Thus, the order of preference for the alignment selection for 3D QSAR model development may be CCBA>DCBA>GMCBA, and the information obtained from the CoMFA and CoMSIA contour maps may be useful in designing specific PTP1B inhibitors.     

4-苯基-1,3-硒唑衍生物的合成及其对蛋白酪氨酸磷酸酯酶-1B抑制活性

蛋白酪氨酸磷酸酯酶-1B(PTP1B)是抗糖尿病治疗的重要靶点,因此创制活性优良的PTP1B抑制剂具有重要意义。 本文设计并合成了11个含1,3-硒唑和1,2,4-三唑活性组块新型结构目标分子(ZLXZ1-ZLXZ11),并利用傅里叶变换红外光谱仪(FTIR)、核磁共振波谱仪(NMR)和高分辨质谱(HRMS)等对其进行了结构表征。 首先选择ZLXZ1和ZLXZ11在MOE 2015.10程序上,与PTP1B进行分子对接模拟,结果表明,在ZLXZ1分子中硒唑环上的硒原子与PTP1B中副催化位点Tyr46、Ala217、Lys120和Asp 48分别形成了π-H作用和氢键作用。 在ZLXZ11分子中硒唑上的硒原子与PTP1B中Asp181、Arg221和Asp48形成了氢键作用。 在分子对接模拟的基础上,测试了11个目标分子的抑制活性,结果表明,所有目标分子的抑制率均在87.02%以上,其中3个目标分子PTP1B抑制活性高于阳性参照物齐墩果酸,抑制活性优良,有望成为潜在的PTP1B抑制剂。     

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提供重要的理论信息。