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

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

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

2,4-二氨基嘧啶类黏着斑激酶抑制剂的3D-QSAR研究

本文选取42个2,4-二氨基嘧啶类黏着斑激酶(FAK)小分子抑制剂,分别以比较分子场分析法(CoMFA)与相似性指数分析法(CoMSIA)构建3D-QSAR模型,评价模板分子、公共骨架点、最小能量优化参数、分子构象等因素对模型优化的影响。分析最优模型中立体场、静电场以及氢键等因素对抑制剂活性的影响,并应用分子对接分析该类抑制剂与FAK蛋白的相互作用。结果表明,选择化合物16作为模板分子,骨架A作为公共骨架点,最小能量优化参数中电荷、最大迭代系数、最低能量限定值分别为MMFF94、1000、0.01kcal/mol时所构建的模型最优。以CoMFA和CoMSIA构建的3D-QSAR模型的交叉验证系数(q~2)分别为0.666和0.736,非交叉验证系数(R~2)分别为0.990和0.989,表明此模型具有良好的预测能力。分子对接分析显示,其与FAK的氨基酸残基CYS502、ASP564形成了重要的氢键作用,并与周围残基形成了较强的疏水作用。通过3D-QSAR的构建与分子对接分析,可指导2,4-二氨基嘧啶类FAK小分子抑制剂的进一步结构优化设计。     

糖原合成激酶-3β马来酰胺类抑制剂构效关系的理论研究

建立糖原合成激酶-3β(GSK-3β)抑制剂的三维构效关系,可预测新的糖原合成激酶-3β抑制剂.通过确定分子的药效构象,与选定的模板分子进行叠合,采用比较分子力场分析法(Co MFA)和比较分子相似性指数分析法(Co MSIA)分别建立38个糖原合成激酶-3β抑制剂的3D-QSAR模型.比较分子力场分析法所建立的模型的决定系数q2=0.711,非交叉验证系数r2=0.887,标准偏差ES=0.411,显著系数F=109.856;比较分子相似性指数分析法所建立模型的决定系数q2=0.605,非交叉验证系数r2=0.931,标准偏差ES=0.326,显著系数F=122.122.该模型在一定程度上反映了结合部位的性质要求,解释马来酰胺类抑制剂的构效关系,具有较好的预测能力.     

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

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

胆固醇酯转运蛋白抑制剂的3D-QSAR模型

利用VolSurf参数和比较分子场分析(CoMFA)方法对N,N-二取代三氟-3-氨基-2-丙醇衍生物类胆固醇酯转运蛋白抑制剂进行了三维定量构效关系(3D-QSAR)模型研究, 均得到较好的结果, 训练集模型具有良好的预测能力. VolSurf参数分析表明抑制活性高的分子必须具有合适的亲水性、多的氢键给体和少的氢键受体; 在一定范围内, 分子量大、表面光滑且非球性高的分子抑制活性高; 高疏水性以及质量中心与疏水区中心的高不平衡性对活性是不利因素. CoMFA结果表明, 立体作用对活性的影响较静电作用稍强, N-苯基取代基苯氧基的间位体积大且正电性强的基团对活性有利, N-苄基取代基的间位体积大且合适的电负性对活性有利, 而苄基的对位立体位阻的增加则对活性不利. VolSurf参数提供了分子整体性质信息, CoMFA提供了取代基信息, 两者互为补充, 对该类抑制剂新化合物的设计具有指导意义.     

人类免疫缺陷病毒整合酶二酮酸类抑制剂的三维药效团构建

应用遗传算法相似性程序(GASP), 以作用于I型人类免疫缺陷病毒(human immun-odeficiency virus type 1, HIV-1)整合酶(IN)的二酮酸类(diketoacids, DKAs)抑制剂构建药效团模型. 所选训练集分子均具有可靠的类药性特征及DKAs药效团特征. 尝试将抑制剂与药效团叠合后的构象和抑制剂与IN的对接构象进行叠合, 得到药效团模型与分子对接构象中IN残基的相对位置, 并基于抑制剂的药效团模型特征与周围IN氨基酸残基位置的匹配情况进行药效团特征的修改. 所得最优药效团由1个疏水特征、3对氢键特征和1个氢键供体特征组成. 该药效团的命中物质量(goodness of hit, GH)为0.56, 产出率(Y)达63.6%, 假阳性率(FP)为0.41%. 该药效团具有较好的置信度, 产出率较高而假阳性率较低, 可用于数据库搜索发现新的具有DKAs药效团特征的活性化合物, 也可为先导化合物的改造提供帮助.     

柔性比较分子场分析方法在糖原磷酸化酶抑制剂三维定量构效关系研究中的应用

通过取消比较分子场分析法(CoMFA)中作为探针安置的网格格点, 并辅以群体智能算法来优化不同类型探针在药物分子周围最佳分布情况, 以此把柔性机制引入到CoMFA当中, 从而形成了一种新的3D-QSAR研究方法: 柔性比较分子场分析法(FCoMFA). 初步尝试使用FCoMFA对47个糖原磷酸化酶抑制剂研究结果表明: FCoMFA能够有效找寻配基与受体的活性位点作用方式并给出具有较强稳定性和预测能力的QSAR模型, 同时可直观通过探针分布模式图观察不同取代基对药物活性的影响情况.     

2,4-二氨基嘧啶类FAK抑制剂的3D-QSAR研究

摘要 本文选取42个2,4-二氨基嘧啶类FAK小分子抑制剂,分别以比较分子场分析法（CoMFA）与相似性指数分析法（CoMSIA）构建3D-QSAR模型,评价模板分子、公共骨架点、最小能量优化参数、分子构象等因素对模型优化的影响。分析最优模型中立体场、静电场以及氢键等因素对抑制剂活性的影响,并应用分子对接分析该类抑制剂与FAK蛋白的相互作用。结果表明选择16号化合物作为模板分子,骨架A作为公共骨架点,最小能量优化参数中电荷、最大迭代系数、最低能量限定值分别为MMFF94、1000、0.01 Kcal/mol时所构建的模型最优。以CoMFA和CoMSIA构建的3D-QSAR模型的交叉验证系数（q2）分别为0.666和0.736,非交叉验证系数（R2）分别为0.990和0.989,表明此模型具有良好的预测能力。分子对接分析显示,其与FAK的氨基酸残基CYS502、ASP564形成了重要的氢键作用,并与周围残基形成了较强的疏水作用。通过3D-QSAR的构建与分子对接分析,可指导2,4-二氨基嘧啶类FAK小分子抑制剂的进一步结构优化设计。     

DATA类逆转录酶抑制剂的三维定量构效关系

采用对接方法得到HIV-1抑制剂DATA(二芳基三嗪类)分子的活性构象, 进一步用比较分子场分析(CoMFA)和比较分子相似性分析(CoMSIA)法对DATA类逆转录酶抑制剂(RTIs)的三维定量构效关系(3D-QSAR)进行了研究, 建立3D-QSAR模型, 以指导进一步结构修饰. 用此模型预测了5个DATA类似物, 预测偏差较小, 表明了所建立的模型具有较强的预测能力.     

Three-dimensional Quantitative Structure-activity Relationship Models of HIV-1 Integrase Inhibitors of DKAs

As one of the three viral encoded enzymes of HIV-1 infection, HIV-1 integrase has become an attractive drug target for the treatment. Diketoacid compounds (DKAs) are one kind of potent and selective inhibitors of HIV-1 IN. In the present work, two three-dimensional QSAR techniques (CoMFA and CoMSIA) were employed to correlate the molecular structure with the activity of inhibiting the strand transfer for 147 DKAs. The all-oritation search (AOS) and all-placement search (APS) were used to optimize the CoMFA model. The diketo and keto-enol tautomers of DKAs were also used to establish the CoMFA models. The results indicated that the enol was the dominant conformation in the HIV-1 IN and DKAs complexes. It can provide a new method and reference to identify the bioactive conformation of drugs by using QSAR analysis. The best CoMSIA model, with five fields combined, implied that the hydrophobic field is very important as well as the steric and electrostatic fields. All models indicated favorable internal validation. A comparative analysis with the three models demonstrated that the CoMFA model seems to be more predictive. The contour maps could afford steric, electrostatic, hydrophobic and H-bond information about the interaction of ligand-receptor complex visually. The models would give some useful guidelines for designing novel and potent HIV-1 integrase inhibitors.     

靶向整合过程不同步骤的HIV-1整合酶抑制剂的结构类型和研究进展

HIV-1整合酶是病毒复制所必需的三个基本酶之一, 为病毒所特有, 人体无对应的酶, 因此整合酶是理想的抗HIV药物设计的新靶标. HIV-1整合酶催化病毒DNA插入宿主染色体的过程涉及到整合酶与前病毒DNA形成整合前复合物、病毒DNA的3’末端切断和DNA链转移等步骤, 目前研究得最多的HIV-1整合酶抑制剂是抑制链转移反应的芳基二酮酸化合物, 其中的电子等排体衍生物Raltegravir (MK-0518)于2007年10月被美国食品药品管理局(FDA)批准上市, 而GS-9137处于三期临床试验, 此外还有多个处于临床前研究和临床阶段的药物. 根据抑制剂的不同作用机理, 本综述介绍了近年来所报道的HIV-1整合酶抑制剂的结构类型、药效团模型、研究进展及化学合成, 将整合酶抑制剂分为链转移反应抑制剂、整合酶-DNA结合抑制剂、整合酶3’端切除反应抑制剂、非专一性整合酶抑制剂以及多肽类抑制剂等几大类. 其中链转移反应抑制剂结构类型最丰富、发展最快. 整合酶抑制剂的出现丰富了高效抗逆转录病毒疗法(HAART), 为多重抗药性艾滋病患者提供了新的有效的治疗方案.     

用分子模拟方法研究HIV-1整合酶与咖啡酰基类抑制剂的相互作用

用分子对接和分子动力学(MD)模拟方法研究了一类咖啡酰基和没食子酰基类HIV-1整合酶抑制剂与整合酶之间的相互作用模式, 结果表明该类抑制剂分子上的两个侧链基团(咖啡酰基或没食子酰基)与整合酶的DDE基序之间的相互作用对抑制整合酶活性起到关键作用. 当侧链基团为没食子酰基时, 可以提高该类抑制剂与整合酶的结合能力. 采用线性相互作用能方法(LIE)计算了该类抑制剂与整合酶之间的结合自由能, 预测值与实验值相吻合, 均方根偏差RMSD为1.39 kJ•mol^-1, 以上结果可为基于结构的HIV-1整合酶抑制剂设计提供有用的信息.     

Pharmacophore and Docking-based 3D-QSAR Studies on HIV-1 Integrase Inhibitors

Integrase（IN） plays an essential role in the process of HIV-1 replication.IN inhibitors of diketo acid derivatives（DKAs） were analysed by the Comparative Molecular Field Analysis（CoMFA） and Comparative Molecular Similarity Induces Analysis（CoMSIA） methods.A set of 42 compounds were randomly selected as the training set（35） and test set（7）.Firstly,a good pharmacophore（goodness of hit=0.787） was obtained and used to align ligands.Then,predictive models were constructed with the CoMFA and CoMSIA methods based on the pharmacophore alignment.As a result,the CoMS1A method yielded the best model with an r2 of 0.955 and a q2 of 0.665,which can predict the activities of the tested DKAs very well（r2=0.559）.Finally,DKAs were docked into IN,and the predicit modes were superimposed on the contour maps obtained from the best CoMSIA model.The superimposed maps gave a visualized and meaningful insight into the inhibitory behaviors,providing significantly useful information for the rational drug design of anti-IN agents.     

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.     

De novo design based identification of potential HIV-1 integrase inhibitors: A pharmacoinformatics study

In the present study, pharmacoinformatics paradigms include receptor-based de novo design, virtual screening through molecular docking and molecular dynamics (MD) simulation are implemented to identify novel and promising HIV-1 integrase inhibitors. The de novodrug/ligand/molecule design is a powerful and effective approach to design a large number of novel and structurally diverse compounds with the required pharmacological profiles. A crystal structure of HIV-1 integrase bound with standard inhibitor BI-224436 is used and a set of 80,000 compounds through the de novo approach in LigBuilder is designed. Initially, a number of criteria including molecular docking, in-silico toxicity and pharmacokinetics profile assessments are implied to reduce the chemical space. Finally, four de novo designed molecules are proposed as potential HIV-1 integrase inhibitors based on comparative analyses. Notably, strong binding interactions have been identified between a few newly identified catalytic amino acid residues and proposed HIV-1 integrase inhibitors. For evaluation of the dynamic stability of the protein-ligand complexes, a number of parameters are explored from the 100 ns MD simulation study. The MD simulation study suggested that proposed molecules efficiently retained their molecular interaction and structural integrity inside the HIV-1 integrase. The binding free energy is calculated through the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach for all complexes and it also explains their thermodynamic stability. Hence, proposed molecules through de novo design might be critical to inhibiting the HIV-1 integrase.     

Efficient 3D database screening for novel HIV-1 IN inhibitors

We describe the use of pharmacophore modeling as an efficient tool in the discovery of novel HIV-1 integrase (IN) inhibitors. A three-dimensional hypothetical model for the binding of diketo acid analogues to the enzyme was built by means of the Catalyst program. Using these models as a query for virtual screening, we found several compounds that contain the specified 3D patterns of chemical functions. Biological testing shows that our strategy was successful in searching for new structural leads as HIV-1 IN inhibitors.     

Synthesis and HIV-1 integrase inhibition activity of some N-arylindoles

Eight simple N-arylindoles were designed, synthesized and evaluated as human immunodeficiency virus (HIV)-1 integrase inhibitors in vitro for the first time. Among these compounds, 3b, 3e and 3g demonstrated significant anti-HIV-1 integrase activity. Especially 3b showed the highest anti-HIV-1 integrase activity with EC50 value of 7.88 microg/ml and TI value of 24.61. Meantime, some structure-activity relationships were also observed and will provide a new lead for design and discovery of more potent N-arylindoles as HIV-1 integrase inhibitors.