基于支持向量机方法的HERG钾离子通道抑制剂分类模型

对human ether-a-gō-gō related genes(HERG)钾离子通道(钾通道)抑制剂,计算了表征分子组成、电荷分布、拓扑、几何结构及物理化学性质等特征的1559个分子描述符.采用Fischer Score(F-Score)排序过滤和Monte Carlo模拟退火法相结合从中筛选与HERG钾通道抑制剂分类相关的分子描述符.采用支持向量机(SVM)方法,分别以IC50=1.0、10.0μmol·L-1为分类标准,建立了三个分类预测模型.对367个训练集分子,用五重交叉验证.得到正、负样本的平均预测精度分别为84.8%-96.6%、80.7%-97.7%,其总的平均预测精度为87.1%-97.2%,优于其它文献报道结果.对97个外部测试集分子,所建三个模型的总样本预测精度在67.0%-90.1%之间,接近或优于其它文献报道结果.     

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.     

磺酰脲类化合物除草活性的QSAR研究

采用密度泛函理论方法, 在B3LYP/6-31G(d)水平下, 计算了23种磺酰类化合物的分子极化率及分子骨架中各原子的Milliken电荷. 提出了一种新的QSAR建模方法, 并据此对其中18种化合物进行多元线性回归分析, 建立了除草活性的预测模型(R=0.96, R2=0.92, r2adj=0.88, F=26.26, q2=0.71, p<0.01, SE=0.36), 对剩余五种化合物进行预测, 结果吻合. 该模型从化合物的亲水性、分子几何特征的角度对如何提高磺酰脲类化合物的除草活性进行了分析, 并对提高化合物除草活性的方法做出预测: 提高苯环和嘧啶环取代基的亲水性, 增加N13周围的电子云密度, 为苯环接入较小的取代基团, 在嘧啶环上接入较大取代基团都可提高化合物的除草活性. 预测结果与3D-QSAR方法的预测结果一致.     

Combined 3D-QSAR modeling and molecular docking study on indolinone derivatives as inhibitors of 3-phosphoinositide-dependent protein kinase-1

3-Phosphoinositide-dependent protein kinase-1 (PDK1) is a promising target for developing novel anticancer drugs. In order to understand the structure-activity correlation of indolinone-based PDK1 inhibitors, we have carried out a combined molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling study. The study has resulted in two types of satisfactory 3D-QSAR models, including the CoMFA model (r(2)=0.907; q(2)=0.737) and CoMSIA model (r(2)=0.991; q(2)=0.824), for predicting the biological activity of new compounds. The detailed microscopic structures of PDK1 binding with inhibitors have been studied by molecular docking. We have also developed docking-based 3D-QSAR models (CoMFA with q(2)=0.729; CoMSIA with q(2)=0.79). The contour maps obtained from the 3D-QSAR models in combination with the docked binding structures help to better interpret the structure-activity relationship. All of the structural insights obtained from both the 3D-QSAR contour maps and molecular docking are consistent with the available experimental activity data. This is the first report on 3D-QSAR modeling of PDK1 inhibitors. The satisfactory results strongly suggest that the developed 3D-QSAR models and the obtained PDK1-inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and in future drug design.     

QSAR modeling on benzo[c]phenanthridine analogues as topoisomerase I inhibitors and anti-cancer agents

Benzo[c]phenanthridine (BCP) derivatives were identified as topoisomerase I (TOP-I) targeting agents with pronounced antitumor activity. In this study, hologram-QSAR, 2D-QSAR and 3D-QSAR models were developed for BCPs on topoisomerase I inbibitory activity and cytotoxicity against seven tumor cell lines including RPMI8402, CPT-K5, P388, CPT45, KB3-1, KBV-1and KBH5.0. The hologram, 2D, and 3D-QSAR models were obtained with the square of correlation coefficient R2 = 0.58-0.77, the square of the crossvalidation coefficient q2 = 0.41-0.60 as well as the external set's square of predictive correlation coefficient r2 = 0.5-0.80. Moreover, the assessment method based on reliability test with confidence level of 95% was used to validate the predictive power of QSAR models and to prevent over-fitting phenomenon of classical QSAR models. Our QSAR model could be applied to design new analogues of BCPs with higher antitumor and topoisomerase I inhibitory activity.