楝酰胺(Rocaglamide)类化合物的三维定量构效关系

采用比较分子力场分析(CoMFA)和比较分子相似因子分析(CoMSIA)方法,对训练集中的26个楝酰胺(Rocaglamide)类化合物进行了三维定量构效关系(3D-QSAR)研究,最终建立的CoMFA模型和CoMSlA模型的q<'2>分别为0.593和0.656.并对测试集中的5个化合物的生物活性进行了预测,结果表明...     

基于5-芳基乙内酰脲类化合物的CoMFA和HQSAR研究

用CoMFA和HQSAR两种QSAR方法研究了50个乙内酰脲类分子的定量构效关系.本研究从构象搜索所得的低能结构出发构建化合物分子的构象, 建立CoMFA模型,并进行了全空间搜索. HQSAR本质上是一种二维的QSAR方法,与CoMFA方法相比,该方法在数据处理方面,比CoMFA方法快捷,并且可重复性好.两种方法均得到了较好分析结果, CoMFA的交叉验证相关系数q2 值为0.815, HQSAR的q2值为0.893.这些方程有力地说明了该类分子在(R,R)-N-3,5-dinitrobenzoyl-1,2-diamine型手性固定相上拆分过程中的影响因素,对今后类似拆分的实验研究提供了理论支持.     

Comparative molecular field analysis of artemisinin derivatives: Ab initio versus semiempirical optimized structures

Based on the belief that structural optimization methods, producing structures more closely to the experimental ones, should give better, i.e. more relevant, steric fields and hence more predictive CoMFA models, comparative molecular field analyses of artemisinin derivatives were performed based on semiempirical AM1 and HF/3-21G optimized geometries. Using these optimized geometries, the CoMFA results derived from the HF/3-21G method are found to be usually but not drastically better than those from AM1. Additional calculations were performed to investigate the electrostatic field difference using the Gasteiger and Marsili charges, the electrostatic potential fit charges at the AM1 level, and the natural population analysis charges at the HF/3-21G level of theory. For the HF/3-21G optimized structures no difference in predictability was observed, whereas for AM1 optimized structures such differences were found. Interestingly, if ionic compounds are omitted, differences between the various HF/3-21G optimized structure models using these electrostatic fields were found.     

基于分子对接的6-萘甲基取代HEPT类HIV-1逆转录酶抑制剂构效关系研究

采用分子对接方法得到了一系列6-萘甲基取代HEPT类逆转录酶抑制剂分子与HIV-1逆转录酶复合物模型,从中抽取出抑制剂分子的活性构象,进一步应用CoMFA和CoMSIA方法建立了具有较好预测能力的3D-QSAR模型,深入探讨了这些化合物的定量构效关系,为进一步的药物设计奠定了良好的基础.另外,以化合物13及其相应的β异构体24为代表,结合量子化学从头算分子轨道理论方法考察了它们的前线轨道,为阐明α和β系列化合物的活性差异提供了理论依据.     

二氢吡啶类化合物的三维定量构效关系

通过分子力学和量子化学计算,得出两种二氢吡啶衍生物的低能构象,再应用比较分子力场分析方法(CoMFA)和比较分子相似性指数分析方法(CoMSIA)分别对两种构象的43个二氢吡啶衍生物进行3D-QSAR研究. 计算结果表明,用两种方法建立的两种构象的构效关系模型均有较好的预测能力.通过分析CoMFA和CoMSIA的系数等势图,直观地了解二氢吡啶衍生物的结构对生物活性的影响,为进一步设计高活性的二氢吡啶衍生物提供一定的理论依据.     

Receptor-based modeling and 3D-QSAR for a quantitative production of the butyrylcholinesterase inhibitors based on genetic algorithm

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for a series of structurally related steroidal alkaloids as butyrylcholinesterase (BuChE) inhibitors. Docking studies were employed to position the inhibitors into the BuChE active site to determine the most probable binding mode. The strategy was to explore multiple inhibitor conformations in producing a more reliable 3D-QSAR model. These multiple conformations were derived using the FlexS program. The conformation selection step for CoMFA was done by genetic algorithm. The genetic algorithm based CoMFA approach was found to be the best. Both CoMFA and CoMSIA yielded significant cross-validated q(2) values of 0.701 and 0.627 and the r(2) values of 0.979 and 0.982, respectively. These statistically significant models were validated by a test set of five compounds. Comparison of CoMFA and CoMSIA contour maps helped to identify structural requirements for the inhibitors and serves as a basis for the design of the next generation of the inhibitor analogues. The results demonstrate that the combination of ligand-based and receptor-based modeling with use of a genetic algorithm is a powerful approach to build 3D-QSAR models. These data can be used for the lead optimization process with respect to inhibition enhancement which is important for the drug discovery and development for Alzheimer's disease.     

Exploring the molecular basis for selective cytotoxicity of lamellarins against human hormone-dependent T47D and hormone-independent MDA-MB-231 breast cancer cells

Abstract  

The common structural requirements for cytotoxicity of lamellarins against two human breast cancer cell lines were determined using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. Twenty lamellarins were selected to serve as the training set, whereas another group of six compounds were used as the test set. The best CoMFA and CoMSIA models for both cell lines yielded satisfactory predictive ability with r _cv² values in the range of 0.659–0.728. Additionally, the contour maps obtained from both the CoMFA and CoMSIA models agreed well with the experimental results and may be used in the design of more potent cytotoxic compounds for human breast cancers. Both analyses not only suggested structural requirements of various substituents around the lamellarin skeleton for their cytotoxic activity against both human breast cancer cell lines but also revealed the molecular basis for the differences between the saturated and unsaturated D-rings of the lamellarins.     

几种改进的CoMFA方法比较研究血小板活化因子拮抗剂

由于传统的比较分子场分析(CoMFA)方法本身存在一些缺陷，使得分子的叠合 规则以及叠合分子的空间取向和空间位置等因素对q~2的影响很大，因此相继提出 了几种改进的CoMFA方法。为了优化CoMFA结果，应用传统的CoMFA方法和交叉验证 的R~2引导的区域选择法(q~2-GRS)、全取向搜索法（AOS）、全空间搜索法（APS） 以及比较分子相似性指数(CoMSIA)等四种改进的CoMFA方法，对18个pinusolide类 衍生物这类新发现的血小板活化因子（PAF）拮抗剂进行了比较研究。结果表明四 种改进的CoMFA方法得到的q~2值均比传统CoMFA的高。q~2-GRS方法得到的q~2值有 所提高，但综合结果并不理想，AOS与APS得到的q~2较为理想，而在CoMSIA中， q~2几乎不受空间取向或空间位置的影响。同时我们引人基于样本的偏最小二乘法 （SAMPLS）取代原AOS／APS程序中的传统PLS进行统计分析，明显提高了其运行速 度。最后，根据q~2最高的CoMFA模型和CoMSIA模型设计了几个预测活性更高的 pinusolide类似物。     

香水百合香气成分的气相色谱保留指数三维定量构效关系研究

采用比较分子场分析（CoMFA）和比较分子相似性指数分析（CoMSIA）方法,研究了香水百合中38种香气成分分子结构与气相色谱保留指数值之间的定量构效关系。用外部测试集验证法和留一交叉验证法对模型的稳健性和预测能力进行了检验,并通过CoMSIA模型和CoMFA模型的分子场三维等势图研究了这些化合物分子中不同化学结构对保留指数值的影响。检验结果表明,所建立的CoMSIA模型和CoMFA模型都具有较好的稳健性和预测能力,且能够合理解释结构对保留指数值的影响,可应用于对香水百合香气成分的色谱保留指数值的预测。与CoMFA模型相比,CoMSIA模型的预测准确度更高,在香水百合香气成分的色谱定量构效关系研究中,显然有更好的应用前景。     

HLA-A*0201限制性CTL表位肽的三维定量构效关系的研究

运用比较分子力场(CoMFA)和比较分子相似性指数分析(CoMFA)方法研究了50个HLA-A^*0201限制性CTL表位九肽结构与亲和性间的关系,另外15个表位九肽作为预测集用于检验模型的预测能力.结果表明采用CoMSIA得到的构效关系模型(q^2=0.628,r^2=0.997,F=840.419)要明显优于采用CoMFA得到的构效关系模型.在CoMSIA计算中,当引入疏水场时,三维构效关系模型得到明显改善,通过该三维构效关系模型,可较精确地估算预测集中15个CTL表位肽与HLA-A^*0201间的亲和力(r^2pred=0.743).通过分析分子场等值面图在空间的分布,可以观察到表位肽分子周围的立体及疏水特征对表位肽与HLA-A^*0201间结合亲和力的影响,从而为进一步对CTL表位肽进行结构改造并基于此进行治疗性疫苗分子设计提供理论基础.     

3D-QSAR Study on the Inhibitory Activity of Flavonoids on PIM-1 Kinase

20 Typical flavonoids were selected for study on the interaction between them and PIM-1 kinase with the comparative molecular field analysis method(CoMFA) as well as the comparative molecular similarity index analysis method(CoMSIA) based on molecule docking.3D-QSAR models between these flavonoids and receptor PIM-1 kinase were established.The obtained optimal cross-validation correlation coefficient Q2 for CoMFA model was 0.582,and the non-cross-validation correlation coefficient R2 was 0.955;the corresponding values for CoMSIA model were 0.790 and 0.974,respectively.These two models showed fairly fine stability and predictive ability.In addition,molecule docking results revealed the key residues in the receptor cavity and their specific action ways with flavonoids.     

Computer-aided molecular design of (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule targeted agents in prostate cancer

Microtubules are tube-shaped, filamentous and cytoskeletal proteins that are essential in all eukaryotic cells. Microtubule is an attractive and promising target for anticancer agents. In this study, three-dimensional quantitative structure activity relationships (3D-QSAR) including comparative molecular field analysis, CoMFA, and comparative molecular similarity indices analysis, CoMSIA, were performed on a set of 45 (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule-targeted anti-prostate cancer agents. Automated grid potential analysis, AutoGPA module in Molecular Operating Environment 2009.10 (MOE) as a new 3D-QSAR approach with the pharmacophore-based alignment was carried out on the same dataset. AutoGPA-based 3D-QSAR model yielded better prediction parameters than CoMFA and CoMSIA. Based on the contour maps generated from the models, some key features were identified in (E)-N-Aryl-2-arylethene-sulfonamide analogues that were responsible for the anti-cancer activity. Virtual screening was performed based on pharmacophore modeling and molecular docking to identify the new inhibitors from ZINC database. Seven top ranked compounds were found based on Gold score fitness function. In silico ADMET studies were performed on compounds retrieved from virtual screening in compliance with the standard ranges.     

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.     

Pharmacophore Screening on Piperidinecarboxamides Derivatives Based on GALAHAD and CoMFA Models

In order to discover the novel anticonvulsant drugs, pharmacophore screening of the anticonvulsant inhibitors was enforced. Genetic Algorithm with Linear Assignment for Hypermolecular Alignment of Datasets (GALAHAD) and Comparative Molecular Field Analysis (CoMFA) studies were combined to implement our research. Firstly, multiple models were generated using GALAHAG based on high active molecules. Secondly, several of them were validated using the CoMFA study. Finally, a good values of q² from training set and promising predictive power from test set were obtained based on one model simutaneously. One model had been selected as the most reasonable pharmacophore model. The results of the CoMFA study based on the model 1 suggested that both steric and electrostatic interactions played important roles.     

Combretastatins类微管蛋白抑制剂的定量构效关系与结合模式

以Combretastatins的B环改造化合物为研究对象, 采用遗传函数分析方法进行了二维定量构效关系研究. 研究结果表明, Apol, PMI-mag, Dipole-mag, Hbond donor和RadOfGyration等描述符对该系列抑制剂活性的贡献最大. 采用比较分子场分析方法(CoMFA)和比较分子相似因子分析方法(CoMSIA)进行了三维定量构效关系研究, 建立的CoMFA和CoMSIA模型的交叉验证相关系数q2分别为0.630和0.634, 具有较强的预测能力. 利用CoMFA和CoMSIA模型的三维等势图解析了Combretastatins类化合物的构效关系, 阐明了B环上各取代基对抑制微管蛋白聚合活性的影响, 同时应用分子对接方法分析并验证了定量构效关系模型.     

Improved CoMFA modeling by optimization of settings

The possibility of improving the predictive ability of comparative molecular field analysis (CoMFA) by settings optimization has been evaluated to show that CoMFA predictive ability can be improved. Ten different CoMFA settings are evaluated, producing a total of 6120 models. This method has been applied to nine different data sets, including the widely used benchmark steroid data set, as well as eight other data sets proposed as QSAR benchmarking data sets by Sutherland et al. (J. Med. Chem. 2004, 47, 5541-5554). All data sets have been studied using training and test sets to allow for both internal (q(2)) and external (r(2)(pred)) predictive ability assessment. CoMFA settings optimization was successful in developing models with improved q(2) and r(2)(pred) as compared to default CoMFA modeling. Optimized CoMFA is compared with comparative molecular similarity indices analysis (CoMSIA) and holographic quantitative structure-activity relationship (HQSAR) models and found to consistently produce models with improved or equivalent q(2) and r(2)(pred). The ability of settings optimization to improve model predictive ability has been validated using both internal and external predictions, and the risk of chance correlation has been evaluated using response variable randomization tests.     

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.     

3 D-QSAR Analysis of Agonists of nAChRs： Epibatidine Analogues

A 3 D-QSAR about nAChRs agonists epibatidine analogues was performed using theCoMFA and CoMSIA. The correlation coefficients were R2cv = 0.546, R2cv = 0.907 in CoMFA andR2cv = 0.655, R2,~ = 0.962 in CoMSIA of the final model. The prediction using the final models tothe test set was r2 = 0.675 in CoMFA and r2 = 0.462 in CoMSIA. This model will be useful in thedesign of novel compounds with high affinity.