R基团搜索技术用于硼酸三肽蛋白酶体抑制剂的分子设计

本文采用Topomer CoMFA对44个Tyropeptin硼酸三肽类蛋白酶体抑制剂进行三维定量构效关系(3D-QSAR)分析。所得最优模型的拟合、交互验证、及外部验证的复相关系数分别为0.983、0.651、0.963。采用Topomer search对ZINC数据库进行R基团的虚拟筛选,得到具有特定活性贡献的R基团,以活性最高的分子为模板过滤,得到7个R1和5个R2基团。并以此设计得到活性优于模板分子的20个新化合物。结果表明,所建立的Topomer CoMFA模型具有良好的稳定性和预测能力,基于R集团的Topomer search技术可以有效筛选,并为设计出新的蛋白酶体抑制剂提供理论依据。     

含1,6-二氢-6-哒嗪酮的1,3,4-噁二唑啉和噻唑烷-4-酮类化合物的合成

以1,6-二氢-6-哒嗪酮-3-甲酰肼(1)与芳香醛反应得到相应的1,6-二氢-6-哒嗪酮-3-羰基芳香醛腙(2a～2e). 再将2a～2e与乙酸酐作用, 合环得到一系列含有1,3,4-噁二唑啉环的衍生物3a～3e; 在DMF中用HSCH₂COOH合环, 得到含有4-噻唑烷酮的衍生物4a～4e. 所有化合物的结构均经元素分析, IR, ¹H NMR 和MS谱得以证实.     

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

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

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

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

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.     

3D-QSAR Studies on the Anti-tumor Activity of N-Aryl-salicylamide Derivatives

In the present work, comparative molecular field analysis(CoMFA) techniques were used to perform three-dimensional quantitative structure-activity relationship(3 D-QSAR) studies on the anti-tumor activity(pHi, i = 1, 2, 3, 4) of N-aryl-salicylamide derivatives against four cancer cell lines, including A549, MCF-7, SGC-7901, and Bel-7402. 12 compounds were randomly selected as the training set to establish the prediction models, which were verified by the test set of 5 compounds containing template molecule. The contributions of steric and electrostatic fields to pH1, pH2, pH_3, and pH_4were 23.8% and 76.2%, 20.1% and 79.9%, 18.7% and 81.3%, and 14.3% and 85.7%, respectively. The cross-validation(Rcv2) and non-cross-validation coefficients(R2) were 0.826 and 0.963 for pH1, 0.867 and 0.974 for pH2, 0.941 and 0.989 for pH_3, and 0.797 and 0.961 for pH_4, respectively. The CoMFA models were then used to predict the activities of the compounds, and it was found that the models had strong stability and good predictability. Based on the CoMFA contour maps, some key structural factors responsible for the anticancer activity of the series of compounds were revealed. The results provide some useful theoretical references for understanding the mechanism of action, designing new N-aryl-salicylamide derivatives with high anti-tumor activity, and predicting their activities.     

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.     

苯并噁嗪酮衍生物的3D-QSAR分析

苯并噁嗪酮衍生物是近年来发现的一类抗血小板聚集化合物,在前人研究的基础上利用比较分子场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)对23个苯并噁嗪酮衍生物进行了三维定量构效关系(3D-QSAR)研究.其中CoMFA模型交叉验证系数Q2=0.703,回归系数R2=0.994,计算值与实验值的平均方差SEE=0.053,统计方差比F=184.773;CoMSIA模型Q2=0.847,R2=0.992,SEE=0.058,F=171.670.两种方法得到的模型都具有较好的预测能力.结果表明,标题化合物中8-位取代基R1静电效应起主要作用;2-位取代基R2立体效应占主导作用,但官能团大小要适中.根据研究结果设计了六种活性较高的化合物.     

含呋喃环双酰肼类衍生物的合成、杀虫活性及3D-QSAR研究

为发现新的昆虫生长调节剂,经单取代苯基呋喃甲酰氯与取代苯甲酰肼反应得到22个未见文献报道的含呋喃环双酰肼类化合物,其结构均通过了IR,1HNMR和元素分析确认.初步生测结果表明,所有目标化合物对豆蚜(Aphisfabae)均有活性,部分目标化合物表现出较好或中等的杀幼虫活性.化合物Ia,Ib和Ic在药剂浓度为0.05%时,对豆蚜的死亡抑制率分别为81.8%,58.4%和52.2%,其中化合物Ia对若蚜的蜕皮和成蚜产雌能力具有一定的抑制作用.而大部分目标化合物在药剂浓度为0.1%,0.05%和0.001%时,对3龄粘虫(Mythimna separate)、棉红蜘蛛(Tetranchus urticae)和尖音淡色库蚊(Culex pipiens pallens)幼虫杀虫活性不明显.采用比较分子力场分析(CoMFA)方法,对22个化合物的杀蚜虫活性进行三维定量构效关系(3D-QSAR)研究.在CoMFA研究中,考察了不同力场和电荷下网格点步长对统计结果的影响.建立了三维定量构效关系CoMFA模型(q2=0.518,r2=0.936).CoMFA模型的立体场、静电场三维等值线图不仅直观地解释了结构与活性的关系,而且为后续优化该系列化合物提供了理论依据.     

Molecular Docking and 3D-QSAR Studies on a Series of Fused Heterocyclic Amides as B-Raf Inhibitors

B-Raf has been identified as promising targets for novel anticancer agents. To further explore the interactions between small molecules and B-Raf, and to elucidate structural characteristics that influence the B-Raf kinase activity, molecular docking and three-dimensional quantitative structure-activity relationship(3D-QSAR) studies were performed on a dataset of 75 Type Ⅱ inhibitors. Molecular docking was applied to explore the detailed binding process between the inhibitors and B-Raf kinase in its DFG-out inactive conformation. Based on the conformations obtained by molecular docking strategy, 3D-QSAR models, including comparative molecular field analysis(CoMFA) and comparative molecular similarity indexes analysis(CoMSIA), were constructed. The established 3D-QSAR models show significant statistical quality and satisfactory predictive ability, with high q~2 and r~2 values: CoMFA model(q~2= 0.759, r~2 = 0.922), and CoMSIA model(q~2 = 0.685, r~2 = 0.945). The systemic external validation indicated that both CoMFA and CoMSIA models were quite robust and possess high predictive abilities with r~2 pred values of 0.633 and 0.708, respectively. Several key structural features accounting for the inhibitory activities of these compounds were discussed based on the 3D contour maps generated by the CoMFA and CoMSIA models, which were in good accordance with the docking results. These theoretical results rendered by 3D-QSAR models along with the docking may provide a useful reference for understanding the action mechanism and designing novel potential B-Raf inhibitors.     

脒类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酶抑制剂提供了指导.     

2D-QSAR Using MLR and 3D-QSAR Using CoMSIA Studies on the Toxicity of Aromatic Hydrocarbons on Larval Sinonvaculina Constricta

Aromatic hydrocarbons,one of the persistent organic pollutants(POPs),has been usually found in mussels,accumulated for their hard mobility and activities in harbours and estuaries.In this study,based on the 96 hr-LC50 of 12 aromatic hydrocarbons with larval sinonvaculina constricta,three-dimensional quantitative structure-activity relationship(3D-QSAR) technique:comparative molecular similarity indices analysis(CoMSIA) and 2D-QSAR technique:multiple linear regression(MLR) were described to obtain more detailed insight into the structure-activity relationships between the molecular structure and bio-activity.The results show the MLR model based on density functional theory(DFT) calculation carried out at the B3LYP/6-311** level with Gaussian 03 program yielded a very good correlation with a coefficient squared R2 of 0.716 and a cross-validated Q2 of 0.874.The dipole moment and enthalpy,as the thermodynamic parameters,were two important factors influencing pLC50.Correspondingly,CoMSIA based on the partial least-squares(PLS) methodology with steric,electrostatic,hydrophobic,H-bond donor and acceptor fields contributing simultaneously were employed and the values of R2 and the cross validation with leave-One-Out(LOO) Q2LOO were 0.585 and 0.990,respectively,which reveals the structure features,such as the electronegative substituent(nitro-group),hydrophobic groups(the benzene ring) and H-bond(nitro-group),related to the toxicity.The results of 2D-QSAR employing MLR model and 3D-QSAR employing CoMSIA model provide the useful information for predicting the toxicity of other aromatic hydrocarbons by comparing the molecular structures of similar compounds.     

Molecular Docking and 3D-QSAR Research of Biphenyl Carboxylic Acid MMP-3 Inhibitors

The molecular docking by LigandFit docking of Discovery Studios 2.5 was employed to the three-dimensional quantitative structure-activity relationship(3D-QSAR) studies of biphenyl carboxylic acid MMP3 inhibitors.A significant correlation coefficient was obtained between dock scores and biological activities.Based on the optimal docking conformations,3D-HoVAIF was employed to the QSAR studies of 51 biphenyl carboxylic acid MMP-3 inhibitors.R2 and Q_CV2(leave-one-out,LOO) of the optimal 3D-HoVAIF-PLS model were 0.873 and 0.841 respectively.The conclusions obtained from the PLS analysis were in agreement with the docking results.     

Docking and 3D-QSAR Studies on the Imidazo[1,5-c]pyrimidine Derivative as EED Inhibitors

Embryonic ectoderm development(EED) has become a novel target for cancer treatment. In this study, a series of EED inhibitors was subjected to a three-dimensional quantitative structure-activity relationship(3D-QSAR) and molecular docking. Accordingly, this is the first of such 3D-QSAR studies in a series of EED inhibitors displaying anti-cancer pharmacological profiles. The CoMFA(q~2 = 0.792, r~2 = 0.994, r~2_(pred) = 0.74) and CoMSIA(q~2 = 0.873, r~2 = 0.994, r~2_(pred) = 0.81) models demonstrated good robustness and predictive ability. Moreover, molecular docking suggested that cation-π,π-π stacking and hydrogen bonding interactions were the main factors affecting the activity of these inhibitors. Five new small molecules were designed based on the CoMFA and CoMSIA contour maps. These molecules were then submitted to further ADME studies, in which the ADME properties of the five designed molecules were found to be within a reasonable range. In view of the corresponding findings, this study may provide theoretical guidance for the rational design of novel EED inhibitors.     

Design of 5H-pyrrolo[2,3-b]pyrazine-2-phenyl Ether Derivatives as Potential JAK3 Inhibitor Based on Surflex Docking,3D-QSAR Modeling and Reverse Docking

Tyrosine protein kinase JAK3 has a very important significance on organ transplantation and the treatment of autoimmune diseases, which has been a potential therapeutic target. In recent years, a large number of JAK3 inhibitors have been reported. However, the poor selectivity and side effects have limited their widespread use in clinical practice. In order to solve this problem, 52 potential small-molecule inhibitors were combined with JAK1, JAK2 and JAK3 respectively to obtain the optimal conformation of small molecules. On the basis of that we established 3D quantitative structure-activity relationships(3D-QSAR) model. Comparative molecular field analysis(Co MFA) and molecular similarity analysis(Co MSIA) were used to evaluate the model. We took advantage of reverse docking to explore the underlying toxicity and side effects. Combining 3D quantitative structure-activity relationships, surflex-dock and reverse docking results, ten 5 H-pyrrolo[2,3-b]pyrazine-2-phenyl ether derivatives based on the most optimal selectivity and activity compound 39 were designed. It can be seen from Co MFA and Co MSIA predicted active values of designed molecules that the selectivity of designed small molecules was improved obviously. Among them, compounds 61 and 62 could become the potential small molecule compounds.     

Molecular Docking, 3D-QSAR and Molecular Dynamics Simulation Studies of Substituted Pyrimidines as Selective Covalent Janus Kinase 3 Inhibitors

Janus kinase 3(JAK3) is a member of Janus kinase(JAK) family, and it represents a promising target for the treatment of immune diseases and cancers. However, no highly selective inhibitors of JAK3 have been developed. For discovering the binding mechanism of JAK3 and these inhibitors, a molecular modeling study combining molecular docking, three-dimensional quantitative structure-activity relationships(3D-QSAR), molecular dynamics and binding free energy calculations was performed on a series of pyrimidine-based compounds which could bind with the unique residue Cys909 of JAK3 kinase as the selective inhibitors of JAK3 in this work. The optimum Co MFA and Co MSIA models were generated based on the conformations obtained by molecular docking. The results showed that the models have satisfactory predicted capacity in both internal and external validation. Furthermore, a 50 ns molecular dynamics simulation was carried out to determine the detailed binding process of inhibitors with different activities. It was demonstrated that hydrogen bond interactions with Leu828, Glu903, Tyr904, Leu905 and Leu956 of JAK3 are significant for activity increase, and the Van der Waals interaction is mainly responsible for stable complex.     

CoMFA 3D-QSAR Analysis of Epothilones Based on Docking Conformation and Alignment

Epothilones belong to a class of novel microtubule stabilizing and anti-mitotic agents, which have a paclitaxel-like mechanism of action. A three-dimensional quantitative structure-activity relationship （3D-QSAR） model was built for epothilones by the method of comparative molecular field analysis （CoMFA） combined with the flexible docking technology. The docking CoMFA model gave a good cross-validated value of q2=0.784 with an optimized component of 6 and the conventional correlation coefficient of r^2=0.985. The statistical results show that the model has good ability to predict the activity of the studied compounds. At last, the docking CoMFA model was analyzed through contour maps complemented with MOLCAD-generated active site potential surface in the α,β-tubulin receptor, which can provide important information for the structure-based drug design.     

雪花胺类化合物的三维构效关系研究

雪花胺是一类重要的乙酰胆碱脂酶抑制剂，有可能发展成治疗阿尔茨海默病的药物. 利用分子力学计算得到了这类化合物的优势构象，并对这些化合物进行了CoMFA研究．发现，对于雪花胺类化合物，影响其药效的主要因素是空间结构，电荷作用力的影响较小．对空间因素的进一步分析发现，对于该类分子，不宜用空阻较大的基团进行取代．由电荷影响的分析得到了在不同位置上取代基所应有的电荷性质．三维定量构效关系研究为基于雪花胺的抑制剂设计提供了方案．