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.     

基于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型手性固定相上拆分过程中的影响因素,对今后类似拆分的实验研究提供了理论支持.     

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.     

A comparative study of quantitative structure-activity relationship methods based on gallic acid derivatives

By using hologram quantitative structure-activity relationship (HQSAR) and comparative molecular field analysis (CoMFA) methods, the relationships between the structures of 49 gallic acid derivatives and their analgesic activity have been investigated to yield statistically reliable models with considerable predictive power. The best HQSAR model was generated using atoms, bond and connectivity as fragment distinction parameters and fragment size 5-7 from a hologram length of 307 with 3 components. High conventional r2 (r2 = 0.825) and cross-validation r2 (r2(cv) = 0.726) values were obtained. CoMFA analyses varying lattice size and location, grid spacing, probe charges and using, Tripos standard and Indicator force field were performed. The best model was developed with 4 components using sp3-hybridized carbon atom with +1.0 charge as probe, grid spacing (2 A), lattice offset (1.0, 3.0, -2.5). The CoMFA model showed a conventional correlation coefficient r2 of 0.889 and across-validation r2(cv) equals to 0.633. The robustness and predictive ability of the HQSAR and CoMFA models have been validated by means of an external test set. The results indicate that both models possess high statistical quality in the prediction of analgesic potency of novel gallic acid analogs.     

Quantitative structure-activity relationships for a series of selective estrogen receptor-beta modulators

The estrogen receptor-beta subtype (ERbeta) is an attractive drug target for the development of novel therapeutic agents for hormone replacement therapy. Hologram quantitative structure-activity relationships (HQSAR) were conducted on a series of 6-phenylnaphthalene and 2-phenylquinoline derivatives, employing values of ERbeta binding affinity. A training set of 65 compounds served to derive the models. The best statistical HQSAR model (q(2) = 0.73 and r(2) = 0.91) was generated using atoms, bonds, connections and donor and acceptor as fragment distinction parameters, and fragment size default (4-7) with hologram length of 199. The model was used to predict the binding affinity of an external test set of 16 compounds, and the predicted values were in good agreement with the experimental results. The final HQSAR model and the information obtained from 2D contribution maps should be useful for the design of novel ERbeta modulators having improved affinity.     

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.     

Quantitative structure-activity relationships for a series of selective estrogen receptor-beta modulators

The estrogen receptor-beta subtype (ERβ) is an attractive drug target for the development of novel therapeutic agents for hormone replacement therapy. Hologram quantitative structure-activity relationships (HQSAR) were conducted on a series of 6-phenylnaphthalene and 2-phenylquinoline derivatives, employing values of ERβ binding affinity. A training set of 65 compounds served to derive the models. The best statistical HQSAR model (q ²?=?0.73 and r ²?=?0.91) was generated using atoms, bonds, connections and donor and acceptor as fragment distinction parameters, and fragment size default (4–7) with hologram length of 199. The model was used to predict the binding affinity of an external test set of 16 compounds, and the predicted values were in good agreement with the experimental results. The final HQSAR model and the information obtained from 2D contribution maps should be useful for the design of novel ERβ modulators having improved affinity.     

Comparative molecular field analysis (CoMFA) model using a large diverse set of natural, synthetic and environmental chemicals for binding to the androgen receptor

A large number of natural, synthetic and environmental chemicals are capable of disrupting the endocrine systems of experimental animals, wildlife and humans. These so-called endocrine disrupting chemicals (EDCs), some mimic the functions of the endogenous androgens, have become a concern to the public health. Androgens play an important role in many physiological processes, including the development and maintenance of male sexual characteristics. A common mechanism for androgen to produce both normal and adverse effects is binding to the androgen receptor (AR). In this study, we used Comparative Molecular Field Analysis (CoMFA), a three-dimensional quantitative structure-activity relationship (3D-QSAR) technique, to examine AR-ligand binding affinities. A CoMFA model with r2 = 0.902 and q2 = 0.571 was developed using a large training data set containing 146 structurally diverse natural, synthetic, and environmental chemicals with a 10(6)-fold range of relative binding affinity (RBA). By comparing the binding characteristics derived from the CoMFA contour map with these observed in a human AR crystal structure, we found that the steric and electrostatic properties encoded in this training data set are necessary and sufficient to describe the RBA of AR ligands. Finally, the CoMFA model was challenged with an external test data set; the predicted results were close to the actual values with average difference of 0.637 logRBA. This study demonstrates the utility of this CoMFA model for real-world use in predicting the AR binding affinities of structurally diverse chemicals over a wide RBA range.     

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

QSAR prediction of antagonistic activity of PCBs towards human PXR by using heuristic method and best subset modeling

Polychlorinated biphenyls(PCBs) can antagonize human pregnane X receptor(hPXR) activation.Such chemicals could pose a serious threat to the reproductive and developmental ability of humans.The quantitative structure activity relationship(QSAR) provides a promising method for the estimation of PCBs’ antagonistic activity.In this investigation,a QSAR model was developed by using heuristic method and best subset modeling(r2 = 0.873,q2LOO=0.742).The built model was validated externally by splitting the original data set into training and prediction sets.The results of the model derived are as follows:r2 = 0.907,q2LOO=0.709,r2pred=0.676,suggesting developed QSAR model had good robustness and predictive ability.The applicability domain(AD) of the model was assessed by Williams plot.The antagonistic activity(?logKi) of 108 PCBs,which are unavailable by experiment at present,was predicted within the applicability domain of the model.The critical structural features related to the activity of PCBs were identified.     

CoMFA modeling of enzyme kinetics: K(m) values for sulfation of diverse phenolic substrates by human catecholamine sulfotransferase SULT1A3

Three-dimensional QSAR models were developed for predicting kinetic Michaelis constant (K(m)) values for phenolic substrates of human catecholamine sulfating sulfotransferase (SULT1A3). The K(m) values were correlated to the steric and electronic molecular fields of the substrates utilizing Comparative Molecular Field Analysis (CoMFA). The evaluated SULT1A3 substrate data set consisted of 95 different substituted phenols, catechols, catecholamines, steroids, and related structures for which the K(m) values were available. The data set was divided in three different subgroups in the initial analysis: (1). for the first CoMFA model substrates with only one reacting hydroxyl group were selected (n = 51), (2).the second model was build with structurally rigid substrates (n = 59), and (3). finally all substrates of the data set were included in the analysis (n = 95). Substrate molecules were aligned using the aromatic ring and the reacting hydroxyl group as a template. After the initial analysis different substrate alignment rules based on the existing knowledge of the SULT1A3 active site structure were evaluated. After this optimization a final CoMFA model was built including all 95 substrates of the data set. Cross-validated q(2) values (leave-one-out and leave-n-out) and coefficient contour maps were calculated for all derived CoMFA models. All four CoMFA models were statistically significant with q(2) values up to 0.624. These predictive QSAR models will provide us information about the factors that affect substrate binding at the active site of human catecholamine sulfotransferase SULT1A3.     

苯并咪唑类缓蚀剂的3D-QSAR研究及分子设计

采用比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 对苯并咪唑衍生物抗盐酸腐蚀的缓蚀性能进行了三维定量构效关系研究, 并使用留一法交叉验证手段对3D-QSAR模型的稳定性及预测能力进行了分析. 结果表明, 立体场、静电场和氢键供体场(电子给体)是影响苯并咪唑缓蚀剂缓蚀性能的主要因素; 所构建的CoMFA模型(q2=0.541, R2=0.996)和CoMSIA模型(q2=0.581, R2=0.987)均具有较好的统计学稳定性和预测能力. 基于3D-QSAR等势图设计出了几种具有较好缓蚀性能的苯并咪唑化合物, 为油气田新型缓蚀剂的研发提供了一种新思路.     

Predicting toxicity of benzene derivatives by molecular hologram derived quantitative structure-activity relationships (QSARS)

Holographic quantitative structure-activity relationship (HQSAR) is an emerging QSAR technique with the combined application of molecular hologram, which encodes the frequency of occurrence of various molecular fragment types, and the subsequent partial least squares (PLS) regression analysis. Based on molecular hologram, alignment-free QSAR models could be rapidly and easily developed with highly statistical significance and predictive ability. In this paper, the toxicity data for a series of 83 benzene derivatives to the autotrophic Chlorella vulgaris (IGC50, negative logarithmic form of 6-h 50% population growth inhibition concentration in mmol/l) were subjected to HQSAR analysis and this resulted in a model with a high predictive ability. The robustness and predictive ability of the model were validated by "leave-one-out" (LOO) cross-validation procedure and an external testing set. The influence of fragment distinction parameters and fragment size on the quality of the HQSAR model have been also discussed.     

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.     

基于启发式变量筛选方法研究与雌激素β受体结合的化合物结构与活性之间的关系

雌激素类化合物由于其对人和野生动物健康的负面影响而受到广泛关注.雌激素受体存在两种亚型(ERα和ERβ),化合物与两种受体亚型在结合活性和化合物结构特征方面存在差异.以31种与雌激素β受体亚型(ERβ)结合的化合物为研究对象,采用启发式变量筛选方法,从1524个变量中筛选出5个与化合物活性(lgRBA)最相关的变量,然后采用多元线性回归(MLR)建立最佳预测模型.模型相关性显著,而且具有良好的稳健性和预测能力(r2=0.829,q2LOO=0.742,r2pred=0.772,q2ext=0.724,RMSEE=0.395).同时揭示了影响化合物与ERβ受体结合的配体化合物分子的结构特征,并对模型的应用域进行了研究.     

Predicting toxicity of benzene derivatives by molecular hologram derived quantitative structure–activity relationships (QSARS)

Holographic quantitative structure–activity relationship (HQSAR) is an emerging QSAR technique with the combined application of molecular hologram, which encodes the frequency of occurrence of various molecular fragment types, and the subsequent partial least squares (PLS) regression analysis. Based on molecular hologram, alignment-free QSAR models could be rapidly and easily developed with highly statistical significance and predictive ability. In this paper, the toxicity data for a series of 83 benzene derivatives to the autotrophic Chlorella vulgaris (IGC⁵⁰, negative logarithmic form of 6-h 50% population growth inhibition concentration in mmol/l) were subjected to HQSAR analysis and this resulted in a model with a high predictive ability. The robustness and predictive ability of the model were validated by “leave-one-out” (LOO) cross-validation procedure and an external testing set. The influence of fragment distinction parameters and fragment size on the quality of the HQSAR model have been also discussed.     

QSAR analyses on avian influenza virus neuraminidase inhibitors using CoMFA, CoMSIA, and HQSAR

The recent wide spreading of the H5N1 avian influenza virus (AIV) in Asia, Europe and Africa and its ability to cause fatal infections in human has raised serious concerns about a pending global flu pandemic. Neuraminidase (NA) inhibitors are currently the only option for treatment or prophylaxis in humans infected with this strain. However, drugs currently on the market often meet with rapidly emerging resistant mutants and only have limited application as inadequate supply of synthetic material. To dig out helpful information for designing potent inhibitors with novel structures against the NA, we used automated docking, CoMFA, CoMSIA, and HQSAR methods to investigate the quantitative structure-activity relationship for 126 NA inhibitors (NIs) with great structural diversities and wide range of bioactivities against influenza A virus. Based on the binding conformations discovered via molecular docking into the crystal structure of NA, CoMFA and CoMSIA models were successfully built with the cross-validated q (2) of 0.813 and 0.771, respectively. HQSAR was also carried out as a complementary study in that HQSAR technique does not require 3D information of these compounds and could provide a detailed molecular fragment contribution to the inhibitory activity. These models also show clearly how steric, electrostatic, hydrophobicity, and individual fragments affect the potency of NA inhibitors. In addition, CoMFA and CoMSIA field distributions are found to be in well agreement with the structural characteristics of the corresponding binding sites. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing novel potent NA inhibitors.     

3-吡啶基醚类化合物的分子全息QSAR研究

采用分子全息定量构效关系(HQSAR, hologram quantitative structure-activity relationship)方法, 研究了28个3-吡啶基醚类化合物对乙酰胆碱α4β2受体的亲和性与它们的分子结构之间的关系, 讨论了分子碎片大小、分子碎片亚结构类型以及分子全息长度对QSAR的影响, 得到了较好的HQSAR模型, 模型的交叉验证系数平方q²＝0.670, 非交叉相关系数平方r²＝0.965, 偏差S＝0.093. 利用HQSAR的颜色编码, 对化合物中不同基团对亲和活性的影响进行了讨论, 对新配体的合成具有一定的指导作用.