CDK4与靛玉红及其衍生物复合物结构的模建

细胞周期蛋白依赖性激酶(Cyclin-dependent Kinases, CDKs)是细胞周期调控的重要因子, 也是治疗癌症的一类重要的药物靶标. 靛玉红是传统中药当归龙荟丸中治疗慢性疾病的有效成分, 靛玉红及其衍生物5-磺酸基-靛玉红对CDKs具有有效的抑制作用. 以获得晶体结构的CDK2与5-磺酸根-靛玉红的复合物为模板, 通过同源模建和分子对接的方法构建出的CDK4与靛玉红及其衍生物的结合模式. 结合CDK4与两种抑制剂的复合物结构, 解释了靛玉红和衍生物5-磺酸基-靛玉红之间抑制活性的差别, 同时也分析了5-磺酸基-靛玉红对CDK2和CDK4不同选择性的原因. 所建CDK4结构为进一步进行基于结构的抗癌药物设计提供了合理的模型.     

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位引入负电性取代基则能提高化合物的生物活性。该模型为进一步设计具有更优选择性和活性的化合物提供了理论依据。     

同源模建和分子对接研究HDAC1/HDAC8的选择性

组蛋白去乙酰化酶(HDACs)是近年来治疗肿瘤的重要靶标之一.由于HDACs包含多种亚型,且各亚型的生理功能存在一定的差异,其选择性抑制剂的开发已成为当前的研发热点.我们通过同源模建的HDAC1结构,与已有的HDAC8晶体结构的活性位点进行比较分析,探讨了对两者选择性有重要影响的残基,为基于受体的选择性抑制剂研究提供重要信息.同时选择了52个HDAC抑制剂,分别建立了HDAC1、HDAC8的活性值与对接打分值的线性回归模型.所建的HDAC1和HDAC8的线性构效关系模型的非交叉验证系数R2分别为0.82和0.80,表明具有一定的统计学意义.利用所建模型对已设计合成的化合物进行了预测,预测结果对HDAC1、HDAC8选择性抑制剂的优化改造提供了一定的指导意义.     

含磷嘧啶类CDK9抑制剂的分子对接、3D-QSAR和分子动力学模拟

选取64个具有潜力的含磷嘧啶类细胞周期依赖性蛋白激酶(CDK9)小分子抑制剂,采用分子对接方法研究了该类小分子与CDK9的结合作用,结果表明,分子构象、氢键形成、疏水性和氨基酸残基Cys106在此类抑制剂与CDK9的结合过程中具有重要作用.在配体叠合的基础上,运用比较分子力场分析(Co MFA)、比较分子相似性指数分析(Co MSIA)和Topomer Co MFA(T-COMFA)研究了分子结构与抑制活性的关系,发现由训练集立体场、静电场和疏水场组合的Co MSIA模型为最优模型,其内部交叉验证相关系数(Q2=0.557)、非交叉验证相关系数(R2=0.959)和外部预测相关系数(r2=0.863)具有统计学意义,该模型的三维等值线图直观显示了化合物的活性与其三维结构的关系.根据这些结果设计了10个具有新结构的含磷嘧啶类化合物,分子对接和分子动力学模拟结果表明,新化合物和CDK9的结合模式与原化合物64相同,自由能分析从理论上证明了新化合物64d的CDK9抑制活性优于化合物64,并且显示含磷基团与残基Asp109的静电场能在化合物与CDK9作用过程中有重要作用.     

3-吡啶基醚类化合物的分子全息QSAR研究

采用分子全息定量构效关系(HQSAR, hologram quantitative structure-activity relationship)方法, 研究了28个3-吡啶基醚类化合物对乙酰胆碱α4β2受体的亲和性与它们的分子结构之间的关系, 讨论了分子碎片大小、分子碎片亚结构类型以及分子全息长度对QSAR的影响, 得到了较好的HQSAR模型, 模型的交叉验证系数平方q²＝0.670, 非交叉相关系数平方r²＝0.965, 偏差S＝0.093. 利用HQSAR的颜色编码, 对化合物中不同基团对亲和活性的影响进行了讨论, 对新配体的合成具有一定的指导作用.     

糖原合成激酶-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.该模型在一定程度上反映了结合部位的性质要求,解释马来酰胺类抑制剂的构效关系,具有较好的预测能力.     

5,8-二甲基乙酰紫草素衍生物QSAR及分子设计

综合运用密度泛函理论、分子力学、统计学及比较分子力场分析等方法研究了抗肿瘤药物5,8-二甲基乙酰紫草素衍生物的二维(2D)、三维(3D)定量构效关系(QSAR).所建最优2D-QSAR方程的交叉验证系数(q2)和拟合相关系数(R2)分别为0.618和0.736.CoMFA模型的q2和非交叉验证系数(r2)分别为0.703和0.982,预测相关系数R2pred为0.746.结果表明,所建立的2D/3D-QSAR模型都具有良好的统计学意义及合理、可信的预报能力,可以预测未知化合物的活性。基于此研究结果,设计了4个具有较高抗肿瘤活性的新化合物.     

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

糖原合成酶激酶-3和细胞周期蛋白依赖性激酶抑制剂paullones的研究进展

对抑制剂paullones的来源、合成,以及与激酶相互作用的分子机理、选择性、细胞效应进行了综述.指出糖原合成酶激酶-3(GSK-3)是参与糖代谢的主要限速酶之一,人类多种疾病均与其活性调节异常有关;细胞周期蛋白依赖性激酶(CDKs)是一类丝氨酸/苏氨酸蛋白激酶,能够促使细胞进行有序的生长、增殖、休眠或进入凋亡.以CDKs为靶点的药物可以阻断细胞周期,控制癌细胞增殖,从而达到抗肿瘤的目的.Paullones对CDKs和GSK-3均具有良好的选择性,是CDKs和GSK-3的有效抑制剂;迄今为止,文献报道的paullones化合物接近120个.     

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抑制剂的设计具有指导意义。     

Understanding and modulating cyclin-dependent kinase inhibitor specificity: molecular modeling and biochemical evaluation of pyrazolopyrimidinones as CDK2/cyclin A and CDK4/cyclin D1 inhibitors

Summary Cyclin-dependent kinases (CDKs) play a key role in regulating the cell cycle. The cyclins, their activating agents, and endogenous CDK inhibitors are frequently mutated in human cancers, making CDKs interesting targets for cancer chemotherapy. Our aim is the discovery of selective CDK4/cyclin D1 inhibitors. An ATP-competitive pyrazolopyrimidinone CDK inhibitor was identified by HTS and docked into a CDK4 homology model. The resulting binding model was consistent with available SAR and was validated by a subsequent CDK2/inhibitor crystal structure. An iterative cycle of chemistry and modeling led to a 70-fold improvement in potency. Small substituent changes resulted in large CDK4/CDK2 selectivity changes. The modeling revealed that selectivity is largely due to hydrogen-bonded interactions with only two kinase residues. This demonstrates that small differences between enzymes can efficiently be exploited in the design of selective inhibitors.     

Differential binding of inhibitors to active and inactive CDK2 provides insights for drug design

The cyclin-dependent kinases (CDKs) have been characterized in complex with a variety of inhibitors, but the majority of structures solved are in the inactive form. We have solved the structures of six inhibitors in both the monomeric CDK2 and binary CDK2/cyclinA complexes and demonstrate that significant differences in ligand binding occur depending on the activation state. The binding mode of two ligands in particular varies substantially in active and inactive CDK2. Furthermore, energetic analysis of CDK2/cyclin/inhibitors demonstrates that a good correlation exists between the in vitro potency and the calculated energies of interaction, but no such relationship exists for CDK2/inhibitor structures. These results confirm that monomeric CDK2 ligand complexes do not fully reflect active conformations, revealing significant implications for inhibitor design while also suggesting that the monomeric CDK2 conformation can be selectively inhibited.     

Virtual screening studies to design potent CDK2-cyclin A inhibitors

The cell division cycle is controlled by cyclin-dependent kinases (CDK), which consist of a catalytic subunit (CDK1-CDK8) and a regulatory subunit (cyclin A-H). Pharmacophore analysis indicates that the best inhibitor model consists of (1) two hydrogen bond acceptors, (2) one hydrogen bond donor, and (3) one hydrophobic feature. The HypoRefine pharmacophore model gave an enrichment factor of 1.31 and goodness of fit score of 0.76. Docking studies were carried out to explore the structural requirements for the CDK2-cyclin A inhibitors and to construct highly predictive models for the design of new inhibitors. Docking studies demonstrate the important role of hydrogen bond and hydrophobic interactions in determining the inhibitor-receptor binding affinity. The validated pharmacophore model is further used for retrieving the most active hits/lead from a virtual library of molecules. Subsequently, docking studies were performed on the hits, and novel series of potent leads were suggested based on the interaction energy between CDK2-cyclin A and the putative inhibitors.     

Combined ligand and structure based approaches for narrowing on the essential physicochemical characteristics for CDK4 inhibition

In the absence of an experimentally determined 3D structure of CDK4 (Cyclin-Dependent Kinase 4), QSARs (Quantitative Structure Activity Relationship) have been explored to rationalize binding affinity in terms of physicochemical and structural parameters. Further, docking on a homology model of CDK4 validated the derived QSARs and predicted the binding mode of this series of inhibitors. Relevant parameters and leave-one-out (LOO) cross-validation (q (2)) as well as an external test set validation (r (2) pred) judged the statistical significance and predictive ability of the models. Docking enabled a better understanding of protein-ligand interaction and provided a mechanistic interpretation in terms of physicochemical characteristics. It identified a unique hydrogen bonding between the imidazole of His-95 and the pyridine nitrogen in the ligand. It rationalized the need for R 2 substituents to be bulky and polar, while the substituent at R 8 to be hydrophobic and comparatively less steric. It also explained why at R 6 a variety of substituents are tolerated and how the presence of methyl at R 5 enhances binding affinity.     

Combined 2D and 3D-QSAR,molecular modelling and docking studies of pyrazolodiazepinones as novel phosphodiesterase 2 inhibitors

Selective inhibition of phosphodiesterase 2 (PDE2) in cells where it is located elevates cyclic guanosine monophosphate (cGMP) and acts as novel analgesic with antinociceptive activity. Three-dimensional quantitative structure–activity relationship (QSAR) studies for pyrazolodiazepinone inhibitors exhibiting PDE2 inhibition were performed using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and Topomer CoMFA, and two-dimensional QSAR study was performed using a Hologram QSAR (HQSAR) method. QSAR models were generated using training set of 23 compounds and were validated using test set of nine compounds. The optimum partial least squares (PLS) for CoMFA-Focusing, CoMSIA-SDH, Topomer CoMFA and HQSAR models exhibited good ‘leave-one-out’ cross validated correlation coefficient (q²) of 0.790, 0.769, 0.840 and 0.787, coefficient of determination (r²) of 0.999, 0.964, 0.979 and 0.980, and high predictive power (r²_pred) of 0.796, 0.833, 0.820 and 0.803 respectively. Docking studies revealed that those inhibitors able to bind to amino acid Gln859 by cGMP binding orientation called ‘glutamine-switch’, and also bind to the hydrophobic clamp of PDE2 isoform, could possess high selectivity for PDE2. From the results of all the studies, structure–activity relationships and structural requirements for binding to active site of PDE2 were established which provide useful guidance for the design and future synthesis of potent PDE2 inhibitors.     

靛玉红类CDK1抑制剂的同源模建、分子对接及3D-QSAR研究

细胞周期蛋白依赖性激酶1的异常表达会导致G2期的停滞及多种肿瘤的发生,故CDK1近年来已成为一个理想的治疗靶点. 本文以细胞分裂调控蛋白2的同源体为模板,同源模建了CDK1的结构,并与靛玉红类小分子抑制剂进行分子对接. 分别运用三种叠合方法进行分子叠合,并在此基础上采用Sybyl 7.1中的比较分子场分析(CoMFA)模块及Discovery Studio 3.0中的三维定量构效关系(3D-QSAR)模块(以下简称为DS)分别建立了3D-QSAR模型. 其中,将分子对接叠合与公共骨架叠合联合运用的叠合方法所得3D-QSAR模型的评价参数是最佳的(CoMFA：q²=0.681,r²=0.909,r_pred.²=0.836; DS：q²=0.579,r²=0.971,r_pred.²=0.795,其中q²为交叉验证系数,r²为非交叉验证系数). 本文的研究结果在对靛玉红类小分子进行结构修饰设计出新的CDK1抑制剂方面,可提供重要的理论基础.     

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

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

Development of CDK4/6 Inhibitors: A Five Years Update

The inhibition of cyclin dependent kinases 4 and 6 plays a role in aromatase inhibitor resistant metastatic breast cancer. Three dual CDK4/6 inhibitors have been approved for the breast cancer treatment that, in combination with the endocrine therapy, dramatically improved the survival outcomes both in first and later line settings. The developments of the last five years in the search for new selective CDK4/6 inhibitors with increased selectivity, treatment efficacy, and reduced adverse effects are reviewed, considering the small-molecule inhibitors and proteolysis-targeting chimeras (PROTACs) approaches, mainly pointing at structure-activity relationships, selectivity against different kinases and antiproliferative activity.     

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

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

Molecular dynamics simulations on the inhibition of Cyclin-Dependent Kinases 2 and 5 in the presence of activators

Interests in CDK2 and CDK5 have stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. Molecular dynamics (MD) simulations have not only become a viable approach to drug design because of advances in computer technology but are increasingly an integral part of drug discovery processes. It is common in MD simulations of inhibitor/CDK complexes to exclude the activator of the CDKs in the structural models to keep computational time tractable. In this paper, we present simulation results of CDK2 and CDK5 with roscovitine using models with and without their activators (cyclinA and p25). While p25 was found to induce slight changes in CDK5, the calculations support that cyclinA leads to significant conformational changes near the active site of CDK2. This suggests that detailed and structure-based inhibitor design targeted at these CDKs should employ activator-included models of the kinases. Comparisons between P/CDK2/cyclinA/roscovitine and CDK5/p25/roscovitine complexes reveal differences in the conformations of the glutamine around the active sites, which may be exploited to find highly selective inhibitors with respect to CDK2 and CDK5.