Eigen value analysis of HIV-1 integrase inhibitors.

A three-dimensional quantitative structure activity relationship using the eigen value analysis (EVA) paradigm applied to 41 HIV-1 integrase inhibitors that inhibit integrase mediated cleavage (3'-processing step) and integration (3'-strand transfer step) in vitro was performed. The training set consisted of 35 molecules from five structurally diverse classes: salicylhydrazines, lichen acids, coumarins, quinones, and thiazolothiazepines. Models derived using semiempirical (MOPAC AM1 and PM3) calculated normal-mode frequencies were compared. The predictive ability of each resultant model was evaluated using a test set comprised of six molecules belonging to a different structural class: hydrazides. Models derived using AM1 method showed considerable internal as well as external predictivity (r(2)(cv) = 0.806, r(2)(pred) = 0.761 for 3'-processing and r(2)(cv) = 0.677, r(2)(pred) = 0.591 for 3'-strand transfer).     

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表位肽进行结构改造并基于此进行治疗性疫苗分子设计提供理论基础.     

Molecular modeling studies on 11H-dibenz[b,e]azepine and dibenz[b,f][1,4]oxazepine derivatives as potent agonists of the human TRPA1 receptor

A computational strategy based on comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) was performed on a series of the 11H-dibenz[b,e]azepine and dibenz[b,f][1,4]oxazepine derivatives as potent agonists of the human TRPA1 receptor. The CoMFA and CoMSIA models resulting from a 21 molecule training set gave r2(cv) values of 0.631 and 0.542 and r2 values of 0.986 and 0.981, respectively. The statistically significant models were validated by a test set of five compounds with predictive r2(pred). values of 0.967 and 0.981 for CoMFA and CoMSIA, respectively. A systemic external validation was also performed on the established models. The information obtained from 3D counter maps could facilitate the design of more potent human TRPA1 receptor agonists.     

CoMFA and CoMSIA studies on 1,3-bis(benzylidene)-3,4-dihydro-1H-naphthalen-2-one, 2,6-bis(benzylidene)cyclohexanone, and 3,5-bis(benzylidene)-4-piperidone series of cytotoxic compounds

A number of 1,3-bis(benzylidene)-3,4-dihydro-1H-naphthalen-2-ones, 2,6-bis(benzylidene)cyclohexanones, and 3,5-bis(benzylidene)-4-piperidones possess significant potencies toward L1210, Molt 4/C8, and CEM cell lines. The objective of the current 3D QSAR study is to discover some of the structural parameters which govern cytotoxic potencies. The CoMFA models with steric and electrostatic fields provided satisfactory statistical data [(r2cv = 0.485, r2ncv = 0.834, r2pred = 0.591), (r2cv = 0.532, r2ncv = 0.850, r2pred = 0.729), and (r2cv = 0.561, r2ncv = 0.864, r2pred = 0.666)] in regard to the cytotoxic potencies observed toward L1210, Molt 4/C8, and CEM cell lines, respectively. The CoMSIA model with steric, electrostatic, hydrophobic, and H-bond donor fields exhibited r2cv = 0.513, r2ncv = 0.833, and r2pred = 0.562 for cytotoxic activity toward L1210 cells, while the best CoMSIA models were obtained by a combination of steric, electrostatic, and hydrophobic fields which yielded statistically significant data [(r2cv = 0.531, r2ncv = 0.828, r2pred = 0.652) and (r2cv = 0.560, r2ncv = 0.841, r2pred = 0.729)] to explain the cytotoxicity toward Molt 4/C8 and CEM cells, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps can be used in the design of more potent cytotoxins.     

3D QSAR and molecular docking studies of benzimidazole derivatives as hepatitis C virus NS5B polymerase inhibitors

The urgent need for novel HCV antiviral agents has provided an impetus for understanding the structural requisites of NS5B polymerase inhibitors at the molecular level. Toward this objective, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of 67 HCV NS5B polymerase inhibitors were performed using two methods. First, ligand-based 3D QSAR studies were performed based on the lowest energy conformations employing the atom fit alignment method. Second, receptor-based 3D QSAR models were derived from the predicted binding conformations obtained by docking all NS5B inhibitors at the allosteric binding site of NS5B (PDB ID: 2dxs). Results generated from the ligand-based model were found superior (r2cv values of 0.630 for CoMFA and 0.668 for CoMSIA) to those obtained by the receptor-based model (r2cv values of 0.536 and 0.561 for CoMFA and CoMSIA, respectively). The predictive ability of the models was validated using a structurally diversified test set of 22 compounds that had not been included in a preliminary training set of 45 compounds. The predictive r2 values for the ligand-based CoMFA and CoMSIA models were 0.734 and 0.800, respectively, while the corresponding predictive r2 values for the receptor-based CoMFA and CoMSIA models were 0.538 and 0.639, respectively. The greater potency of the tryptophan derivatives over that of the tyrosine derivatives was interpreted based on CoMFA steric and electrostatic contour maps. The CoMSIA results revealed that for a NS5B inhibitor to have appreciable inhibitory activity it requires hydrogen bond donor and acceptor groups at the 5-position of the indole ring and an R substituent at the chiral carbon, respectively. Interpretation of the CoMFA and CoMSIA contour maps in context of the topology of the allosteric binding site of NS5B provided insight into NS5B-inhibitor interactions. Taken together, the present 3D QSAR models were found to accurately predict the HCV NS5B polymerase inhibitory activity of structurally diverse test set compounds and to yield reliable clues for further optimization of the benzimidazole derivatives in the data set.     

3D-QSAR Analysis of a Series of Dihydroquinolizinone Derivatives as a Hepatitis B Virus Expression Inhibitor

In this study, we explored a three-dimensional quantitative structure-activity relationship(3D-QSAR) model of 63 HBV viral gene expression inhibitors containing dihydroquinolizinones. Two high predictive QSAR models have been built, including comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA). The internal validation parameter(CoMFA, q~2 = 0.701, r~2 = 0.999; CoMSIA, q~2 = 0.721, r~2 = 0.998) and external validation parameter(CoMFA, r~2_(pred = 0.999); CoMSIA, r~2_(pred = 0.999)) indicated that the models have good predictive abilities and significant statistical reliability. We designed several molecules with potentially higher predicted activity on the basis of the result of the models. This work might provide useful information to design novel HBV viral gene expression inhibitors.     

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.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

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 and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q² = 0.589 and r² = 0.927 and q² = 0.473 and r² = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.     

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.     

A 3D-QSAR CoMSIA study on 3-azolylmethylindoles as anti-leishmanial agents

A three-dimensional quantitative structure-activity relationship (3D-QSAR) study using Comparative Molecular Similarity Indices Analysis (CoMSIA) was conducted on a series of 3-azolylmethylindoles as anti-leishmanial agents. Evaluation of 24 compounds synthesized in our laboratory served to establish the model. A random search was performed on the library of compounds, and molecules of the training set were aligned on common elements of template molecule 13, one of the most active compounds. The best predictions were obtained from multifit procedure with a CoMSIA model combining steric, electrostatic, hydrophobic and hydrogen bond acceptor fields (q2 = 0.594, r2 = 0.897). The model was validated using an external test set of 7 compounds giving a satisfactory predictive r2 value of 0.649. Information obtained from CoMSIA contour maps could be used for further design of more promising inhibitors.     

3D-QSAR Analysis of a Series of 1,2,3-Triazole-chromenone Derivatives as an Acetylcholinesterase Inhibitor against Alzheimer's Disease

Chromenones have attracted much attention since they are excellent acetylcholinesterase inhibitor(AChEi). The 1,2,3-triazoles are multifunctional anti-acetylcholinesterase(AChE) agents. In this paper, we report the three-dimensional quantitative structure-activity relationship(3D-QSAR) study of 25 1,2,3-triazolechromenone derivatives based comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA). To construct CoMFA and CoMSIA models, the 25 active molecules were randomly divided into the training and test sets. The obtained cross-validation Q~2 of the CoMFA model, the coefficient of non-cross-validation R~2, and the test value F are 0.597, 0.994, and 396.726, respectively. The cross-validation Q~2 of the CoMSIA model, the coefficient of the non-cross-validation R~2, and the test value F are 0.721, 0.979, and 131.107, respectively. The predictive correlation coefficient(r_(pred)~2) is 0.728 for CoMFA and 0.805 for CoMSIA, which verifies that the model is predictable. Based on the potential maps of CoMFA and CoMSIA, a library containing a set of potent AChEi was designed. The inhibitory potential of the compounds in this library was found to be greater than the inhibitory potential of the most active compounds in the data set. The results obtained from this study laid the foundation for the development of effective drugs for AChEi.     

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

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

3D QSAR studies on protein tyrosine phosphatase 1B inhibitors: comparison of the quality and predictivity among 3D QSAR models obtained from different conformer-based alignments

A set of 65 flexible peptidomimetic competitive inhibitors (52 in the training set and 13 in the test set) of protein tyrosine phosphatase 1B (PTP1B) has been used to compare the quality and predictive power of 3D quantitative structure-activity relationship (QSAR) comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for the three most commonly used conformer-based alignments, namely, cocrystallized conformer-based alignment (CCBA), docked conformer-based alignment (DCBA), and global minima energy conformer-based alignment (GMCBA). These three conformers of 5-[(2S)-2-({(2S)-2-[(tert-butoxycarbonyl)amino]-3-phenylpropanoyl}amino)3-oxo-3-pentylamino)propyl]-2-(carboxymethoxy)benzoic acid (compound number 66) were obtained from the X-ray structure of its cocrystallized complex with PTP1B (PDB ID: 1JF7), its docking studies, and its global minima by simulated annealing. Among the 3D QSAR models developed using the above three alignments, the CCBA provided the optimal predictive CoMFA model for the training set with cross-validated r2 (q2)=0.708, non-cross-validated r2=0.902, standard error of estimate (s)=0.165, and F=202.553 and the optimal CoMSIA model with q2=0.440, r2=0.799, s=0.192, and F=117.782. These models also showed the best test set prediction for the 13 compounds with predictive r2 values of 0.706 and 0.683, respectively. Though the QSAR models derived using the other two alignments also produced statistically acceptable models in the order DCBA>GMCBA in terms of the values of q2, r2, and predictive r2, they were inferior to the corresponding models derived using CCBA. Thus, the order of preference for the alignment selection for 3D QSAR model development may be CCBA>DCBA>GMCBA, and the information obtained from the CoMFA and CoMSIA contour maps may be useful in designing specific PTP1B inhibitors.     

3D-QSAR and molecular modeling of HIV-1 integrase inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) methods were applied on a series of inhibitors of HIV-1 integrase with respect to their inhibition of 3-processing and 3-end joining steps in vitro.The training set consisted of 27 compounds belonging to the class of thiazolothiazepines. The predictive ability of each model was evaluated using test set I consisting of four thiazolothiazepines and test set II comprised of seven compounds belonging to an entirely different structural class of coumarins. Maximum Common Substructure (MCS) based method was used to align the molecules and this was compared with other known methods of alignment. Two methods of 3D QSAR: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were analyzed in terms of their predictive abilities. CoMSIA produced significantly better results for all correlations. The results indicate a strong correlation between the inhibitory activity of these compounds and the steric and electrostatic fields around them. CoMSIA models with considerable internal as well as external predictive ability were obtained. A poor correlation obtained with hydrophobic field indicates that the binding of thiazolothiazepines to HIV-1 integrase is mainly enthalpic in nature. Further the most active compound of the series was docked into the active site using the crystal structure of integrase. The binding site was formed by the amino acid residues 64-67, 116, 148, 151-152, 155-156, and 159. The comparison of coefficient contour maps with the steric and electrostatic properties of the receptor shows high level of compatibility.     

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

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

3D-QSAR Studies on 3-Substituted N-Benzhydryl-nortropane Analogs as Nociception Receptor Agonists

The nociceptin receptor(NOP) has been involved in multiple biological functions, including pain, anxiety, cough, substance abuse, cardiovascular control, and immunity. Thus, selective NOP agonists might have clinical potential for the treatment of related diseases. In the present work, three-dimensional quantitative structure-activity relationship(3D-QSAR) studies were performed on a series of 3-substituted N-benzhydryl-nortropane analogs as NOP agonists using comparative molecular field analysis(Co MFA) and comparative molecular similarity indices analysis(CoM SIA) techniques. The statistically significant models were obtained with 54 compounds in training set by ligand-based atom-by-atom matching alignment. The CoM FA model gave cross-validated coefficient(q2) value of 0.530 using 6 components, non-cross-validated(r2) value of 0.921 with estimated F value of 93.668, and standard error of estimate(SEE) of 0.185. The best Co MSIA model resulted in q2 = 0.592, r2 = 0.945, N = 10, SEE = 0.162, and F = 75.654, based on steric, electrostatic, hydrophobic and hydrogen bond acceptor fields. The predictive ability of the Co MFA and CoM SIA models was further validated using a test set of 18 molecules that were not included in the training set, which resulted in predictive correlation coefficients(r2pred) of 0.551 and 0.637, respectively. Moreover, the CoM FA and CoM SIA contour maps identified the features important for exhibiting potent binding affinities on NOP, and can thus serve as a useful guide for the design of potential NOP agonists.