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.     

CoMFA modeling of human catechol O-methyltransferase enzyme kinetics

Three-dimensional QSAR models with different charge calculation methods (MOPAC-AM1-ESP, MOPAC-AM1-Coulson and Gasteiger-Hückel) were developed for predicting all three enzyme kinetic parameters Km, Vmax and Vmax/Km for catecholic substrates of human soluble catechol O-methyltransferase (S-COMT). The empirical parameters of 45 substrates were correlated to the steric and electronic molecular fields of the substrates utilizing Comparative Molecular Field Analysis (CoMFA). Alignment rules for CoMFA were developed based on the catalytic mechanism and crystal structure of S-COMT, and the analysis was optimized using an all-space search technique. Leave-one-out and leave-n-out cross-validation (with 5 and 10 cross-validation groups) was carried out, and all developed models proved to be statistically significant with q2 values up to 0.84. The models based on MOPAC charge calculations predicted the empirical values clearly better than the Gasteiger-Hückel method. The derived CoMFA coefficient contour maps of steric and electrostatic interactions correlated clearly with the S-COMT crystallographic structures.     

基于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 （CoMFA） of Curcumin-related Compounds for Anticancer Activity

Three-dimensional （3D） quantitative structure-activity relationship （QSAR） studies of 44 curcumin-related compounds have been carried out based on our previously reported result for their anticancer activity against pancreas cancer Panc-I cells and colon cancer HT-29 cells. The established 3D-QSAR models from the comparative molecular field analysis （CoMFA） in training set showed not only significant statistical quality, but also satisfying predictive ability, with high correlation coefficient values （R12= 0.911, R22= 0.985） and cross-validation coefficient values （q2= 0.580, q22= 0.722）. Based on the CoMFA contour maps, some key structural factors responsible for anticancer activity of these series of compounds were revealed. The results provide some useful theoretical references for understanding the mechanism of action, designing new curcumin-related compounds with anticancer activity and predicting their activities prior to synthesis.     

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.     

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.     

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.     

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

Structure-based rational quest for potential novel inhibitors of human HMG-CoA reductase by combining CoMFA 3D QSAR modeling and virtual screening

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) catalyzes the formation of mevalonate. In many classes of organisms, this is the committed step leading to the synthesis of essential compounds, such as cholesterol. However, a high level of cholesterol is an important risk factor for coronary heart disease, for which an effective clinical treatment is to block HMGR using inhibitors like statins. Recently the structures of catalytic portion of human HMGR complexed with six different statins have been determined by a delicate crystallography study (Istvan and Deisenhofer Science 2001, 292, 1160-1164), which established a solid basis of structure and mechanism for the rational design, optimization, and development of even better HMGR inhibitors. In this study, three-dimensional quantitative structure-activity relationship (3D QSAR) with comparative molecular field analysis (CoMFA) was performed on a training set of up to 35 statins and statin-like compounds. Predictive models were established by using two different ways: (1) Models-fit, obtained by SYBYL conventional fit-atom molecular alignment rule, has cross-validated coefficients (q2) up to 0.652 and regression coefficients (r2) up to 0.977. (2) Models-dock, obtained by FlexE by docking compounds into the HMGR active site, has cross-validated coefficients (q2) up to 0.731 and regression coefficients (r2) up to 0.947. These models were further validated by an external testing set of 12 statins and statin-like compounds. Integrated with CoMFA 3D QSAR predictive models, molecular surface property (electrostatic and steric) mapping and structure-based (both ligand and receptor) virtual screening have been employed to explore potential novel hits for the HMGR inhibitors. A representative set of eight new compounds of non-statin-like structures but with high pIC(50) values were sorted out in the present study.     

QSAR models using a large diverse set of estrogens

Endocrine disruptors (EDs) have a variety of adverse effects in humans and animals. About 58,000 chemicals, most having little safety data, must be tested in a group of tiered assays. As assays will take years, it is important to develop rapid methods to help in priority setting. For application to large data sets, we have developed an integrated system that contains sequential four phases to predict the ability of chemicals to bind to the estrogen receptor (ER), a prevalent mechanism for estrogenic EDs. Here we report the results of evaluating two types of QSAR models for inclusion in phase III to quantitatively predict chemical binding to the ER. Our data set for the relative binding affinities (RBAs) to the ER consists of 130 chemicals covering a wide range of structural diversity and a 6 orders of magnitude spread of RBAs. CoMFA and HQSAR models were constructed and compared for performance. The CoMFA model had a r2 = 0.91 and a q2LOO = 0.66. HQSAR showed reduced performance compared to CoMFA with r2 = 0.76 and q2LOO = 0.59. A number of parameters were examined to improve the CoMFA model. Of these, a phenol indicator increased the q2LOO to 0.71. When up to 50% of the chemicals were left out in the leave-N-out cross-validation, the q2 remained significant. Finally, the models were tested by using two test sets; the q2pred for these were 0.71 and 0.62, a significant result which demonstrates the utility of the CoMFA model for predicting the RBAs of chemicals not included in the training set. If used in conjunction with phases I and II, which reduced the size of the data set dramatically by eliminating most inactive chemicals, the current CoMFA model (phase III) can be used to predict the RBA of chemicals with sufficient accuracy and to provide quantitative information for priority setting.     

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.     

Insight into the structural requirements of narlaprevir-type inhibitors of NS3/NS4A protease based on HQSAR and molecular field analyses

In the life cycle of hepatitis C virus (HCV), NS3/NS4A protease has been proved to play a vital role in the replication of the HCV virus. Narlaprevir and its derivatives, the inhibitors of NS3/NS4A, would be potentially developed as important anti-HCV drugs in the future. In this study, quantitative structure-activity relationship (QSAR) analyses for 190 narlaprevir derivatives were conducted using comparative molecular field analysis (CoMFA), comparative molecular indices analysis (CoMSIA) and hologram quantitative structure-activity relationship (HQSAR) techniques. Both of the best CoMFA and HQSAR models showed statistical significance for the training set and good predictive accuracy for the test set, which strongly manifested the robustness of the CoMFA and HQSAR models. The CoMFA contour maps and the HQSAR contribution maps were both presented. Furthermore, based on the essential factors for ligand binding derived from the QSAR models, sixteen new derivatives were designed and some of them showed higher inhibitory activities confirmed by our models and molecular docking studies. General speaking, this study provides useful suggestions for the design of potential anti-HCV drugs.     

马来酰胺类糖原合成酶激酶-3β抑制剂的分子对接和三维定量构效关系

通过分子对接和三维定量构效关系(3D-QSAR)两种方法来确定两类马来酰胺类的糖原合成酶激酶-3β(GSK-3β)抑制剂的结合方式.首先,用分子对接确定抑制剂与GSK-3β的结合模式及其相互作用;然后用比较分子力场分析法(CoMFA)与比较分子相似性指数分析法(CoMSIA)对48个化合物做三维定量构效关系的分析.两种方法得出的交互验证回归系数分别为0.669(CoMFA)和0.683(CoMSIA),证明该模型具有很好的统计相关性,同时也说明该模型具有较高的预测能力.根据该模型提供的信息,设计出9个预测性较好的分子.     

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.     

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.     

Comparison of 3D quantitative structure-activity relationship methods: Analysis of the in vitro antimalarial activity of 154 artemisinin analogues by hypothetical active-site lattice and comparative molecular field analysis

Two three-dimensional quantitative structure-activity relationship (3D-QSAR) methods, comparative molecular field analysis (CoMFA) and hypothetical active site lattice (HASL), were compared with respect to the analysis of a training set of 154 artemisinin analogues. Five models were created, including a complete HASL and two trimmed versions, as well as two CoMFA models (leave-one-out standard CoMFA and the guided-region selection protocol). Similar r2 and q2 values were obtained by each method, although some striking differences existed between CoMFA contour maps and the HASL output. Each of the four predictive models exhibited a similar ability to predict the activity of a test set of 23 artemisinin analogues, although some differences were noted as to which compounds were described well by either model.