Improved 3D-QSAR prediction by multiple-conformational alignment: A case study on PTP1B inhibitors

Three-dimension quantitative structure activity relationship (3D-QSAR) was one of the major statistical techniques to investigate the correlation of biological activity with structural properties of candidate molecules, and the accuracy of statistic greatly depended on molecular alignment methodology. Exhaustive conformational search and successful conformational superposition could extremely improve the predictive accuracy of QSAR modeling. In this work, we proposed a solution to optimize QSAR prediction by multiple-conformational alignment methods, with a set of 40 flexible PTP1B inhibitors as case study. Three different molecular alignment methods were used for the development of 3D-QSAR models listed as following: (1) docking-based alignment (DBA); (2) pharmacophore-based alignment (PBA) and (3) co-crystallized conformer-based alignment (CCBA). Among these three alignments, it was indicated that the CCBA was the best and the fastest strategy in 3D-QSAR development, with the square correlation coefficient (r²) and cross-validated squared correlation coefficient (q²) of comparative molecular field analysis (CoMFA) were 0.992 and 0.694; the r² and q² of comparative molecular similarity indices analysis (CoMSIA) were 0.972 and 0.603, respectively. The alignment methodologies used here not only generated a robust QSAR model with useful molecular field contour maps for designing novel PTP1B inhibitors, but also provided a solution for constructing accurate 3D-QSAR model for various disease targets. Undoubtedly, such attempt in QSAR analysis would greatly help us to understand essential structural features of inhibitors required by its target, and so as to discover more promising chemical derivatives.     

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.     

Comparing the quality and predictiveness between 3D QSAR models obtained from manual and automated alignment

A set of 113 flexible cyclic urea inhibitors of human immunodeficiency virus protease (HIV-1 PR) was used to compare the quality and predictive power of CoMFA and CoMSIA models for manually or automatically aligned inhibitor set. Inhibitors that were aligned automatically with molecular docking were in agreement with information obtained from existing X-ray structures. Both alignment methods produced statistically significant CoMFA and CoMSIA models, with the best q(2) value being 0.649 and the best predictive r(2) being 0.754. The manual alignment gave statistically higher values, whereas the automated alignment gave more robust models for predicting the activities of an external inhibitor set. Both models utilized similar amino acids in the HIV-1 PR active site, supporting the idea that hydrogen bonds form between an inhibitor and the backbone carbonyl oxygens of Gly48 and Gly48' and also the backbone NH group of Asp30, Gly48, Asp29', and Gly48' of the enzyme. These results suggest that an automated inhibitor alignment can yield predictive 3D QSAR models that are well comparable to manual methods. Thus, an automated alignment method in creating 3D QSAR models is encouragable when a well-characterized structure of the target protein is available.     

A 3D QSAR study on a set of dopamine D4 receptor antagonists

The molecular alignments obtained from a previously reported pharmacophore model have been employed in a three-dimensional quantitative structure-activity relationship (3D QSAR) study, to obtain a more detailed insight into the structure-activity relationships for D(2) and D(4) receptor antagonists. The frequently applied CoMFA method and the related CoMSIA method were used. Statistically significant models have been derived with these two methods, based on a set of 32 structurally diverse D(2) and D(4) receptor antagonists. The CoMSIA and the CoMFA methods produced equally good models expressed in terms of q(2) values. The predictive power of the derived models were demonstrated to be high. Graphical interpretation of the results, provided by the CoMSIA method, brings to light important structural features of the compounds related to either low- or high-affinity D(2) or D(4) antagonism. The results of the 3D QSAR studies indicate that bulky N-substituents decrease D(2) binding, whereas D(4) binding is enhanced. Electrostatically favorable and unfavorable regions exclusive to D(2) receptor binding were identified. Likewise, certain hydrogen-bond acceptors can be used to lower D(2) affinity. These observations may be exploited for the design of novel dopamine D(4) selective antagonists.     

Triple-layered QSAR studies on substituted 1,2,4-trioxanes as potential antimalarial agents: superiority of the quantitative pharmacophore-based alignment over common substructure-based alignment

This study reports the utilization of three approaches – pharmacophore, CoMFA/CoMSIA and HQSAR studies – to identify the essential structural requirements in 3D chemical space for the modulation of the antimalarial activity of substituted 1,2,4-trioxanes. The superiority of quantitative pharmacophore-based alignment (QuantitativePBA) over global minima energy conformer-based alignment (GMCBA) has been reported in CoMFA and CoMSIA studies. The developed models showed good statistical significance in internal validation (q ², group cross-validation and bootstrapping) and performed very well in predicting the antimalarial activity of test set compounds. Structural features in terms of their steric, electrostatic and hydrophobic interactions in 3D space have been found to be important for the antimalarial activity of substituted 1,2,4-trioxanes. Further, the HQSAR studies based on the same training and test set acted as an additional tool to find the sub-structural fingerprints of substituted 1,2,4-trioxanes for their antimalarial activity. Together, these studies may facilitate the design and discovery of new substituted 1,2,4-trioxanes with potent antimalarial activity.     

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

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

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

苯并呋喃／噻吩联二苯类PTP1B抑制剂三维构效关系研究

主要采用比较分子力场分析方法（CoMFA）对苯并呋喃／噻吩联二苯类PTP1B (protein tyrosine phosphatase 1B)抑制剂进行了三维构效关系的研究，考察了 静电场、立体场和氢键场对构效关系的影响，交叉系数q^2的值达到0．58，表明 CoMFA得到的构效关系模型比较理想，同时test set中分子的预测活性也表明，模 型具有较好的预测能力，研究还表明，氢键场的加入不一定有利于模型的改善，通 过对分子场等值面图的分析，可以观察到叠合分子周围立体场和静电场对化合物活 性的影响，为改进原有化合物的结构，提高它们的活性提供了指导，还尝试采用比 较分子相似性指数分析方法（CoMFA）对这一系列化合物作了研究，结果表明虽然 CoMFA中加入了疏水场，但是对于研究的体系，CoMFA的模型质量并没有显著提高。     

Molecular electrostatic potentials as input for the alignment of HIV-1 integrase inhibitors in 3D QSAR

Comparative molecular similarity indices analysis (CoMSIA), a three-dimensional quantitative structure activity relationship (3D QSAR) paradigm, was used to examine the correlations between the calculated physicochemical properties and the in vitro activities (3'-processing and 3'-strand transfer inhibition) of a series of human immunodeficiency virus type 1 (HIV-1) integrase inhibitors. The training set consisted of 34 molecules from five structurally diverse classes: salicylpyrazolinones, dioxepinones, coumarins, quinones, and benzoic hydrazides. The data set was aligned using extrema of molecular electrostatic potentials (MEPs). The predictive ability of the resultant model was evaluated using a test set comprised of 7 molecules belonging to a different structural class of thiazepinediones. A CoMSIA model using an MEP-based alignment showed considerable internal as well external predictive ability (r2(cv) = 0.821, r2(pred) = 0.608 for 3'-processing; and r2(cv) = 0.759, r2(pred.) = 0.660 for 3'-strand transfer).     

3D QSAR Modeling and Molecular Docking Studies on a Series of Triazole Analogues as Antibacterial Agents

The 3D QSAR analysis using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques is performed on novel nalidixic acid based 1,2,4-triazole derivatives suggested earlier as antibacterial agents. The CoMFA and CoMSIA models employed for a training set of 28 compounds gives reliable values of Q² (0.53 and 0.52, respectively) and R² (0.79 and 0.85, respectively). The contour maps produced by the CoMFA and CoMSIA models are used to determine a three-dimensional quantitative structure-activity relationship. Based on the 3D QSAR contours new molecules with high predicted activities are designed. In addition, surflex-docking is performed to confirm the stability of predicted molecules in the receptor.     

Synthesis of 1,2,4-oxadiazole derivatives: anticancer and 3D QSAR studies

Monatshefte für Chemie - Chemical Monthly - A 3D QSAR study was performed on 1,2,4-oxadiazole derivatives using the [(SW) kNN MFA], CoMFA, and CoMSIA techniques. On the basis of 3D QSAR...     

QM/MM based 3D QSAR models for potent B-Raf inhibitors

Three dimensional (3D) quantitative structure-activity relationship studies of 37 B-Raf inhibitors, pyrazole-based derivatives, were performed. Based on the co-crystallized compound (PDB ID: 3D4Q), several alignment methods were utilized to derive reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models. Receptor-guided alignment with quantum mechanics/molecular mechanics (QM/MM) minimization led to the best CoMFA model (q ² = 0.624, r ² = 0.959). With the same alignment, a statistically reliable CoMSIA model with steric, H-bond acceptor, and hydrophobic fields was also derived (q ² = 0.590, r ² = 0.922). Both models were validated with an external test set, which gave satisfactory predictive r ² values of 0.926 and 0.878, respectively. Contour maps from CoMFA and CoMSIA models revealed important structural features responsible for increasing biological activity within the active site and explained the correlation between biological activity and receptor-ligand interactions. New fragments were identified as building blocks which can replace R1-3 groups through combinatorial screening methods. By combining these fragments a compound with a high bioactivity level prediction was found. These results can offer useful information for the design of new B-Raf inhibitors.     

QSAR and Pharmacophore Studies of Thiazolidine-4-carboxylic Acid Derivatives as Novel Influenza Neuraminidase Inhibitors Using HQSAR, Topomer CoMFA and CoMSIA

In order to understand the chemical-biological interactions governing their activities toward neuraminidase (NA), QSAR models of 28 thiazolidine-4-carboxylic acid derivatives with inhibitory influenza A virus were developed. The obtained HQSAR (hologram quantitative structure activity relationship), Topomer CoMFA and CoMSIA (comparative molecular similarity indices analysis) models were robust and had good exterior predictive capabilities. Moreover, QSAR modeling results elucidated that hydrogen bonds highly contributed to the inhibitory activity, then electrostatic and hydrophobic factors. Squared multiple correlation coefficients (R2) of HQSAR, Topomer CoMFA and CoMSIA models were 0.994, 0.978 and 0.996, respectively. Squared cross-validated correlation coefficients (Q2) of HQSAR, Topomer CoMFA and CoMSIA models were in turn 0.951, 919 and 0.820. Furthermore, squared multiple correlation coefficients for the test set (R2test) of HQSAR, CoMFA and CoMSIA models were 0.879, 0.912 and 0.953, respectively. Squared cross-validated correlation coefficients for the test set (Q2ext) of HQSAR, Topomer CoMFA and CoMSIA models were 0.867, 0.884 and 0.899, correspondingly.     

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.