3-D QSAR studies on histone deacetylase inhibitors. A GOLPE/GRID approach on different series of compounds

Docking simulation and three-dimensional quantitative structure-activity relationships (3D-QSARs) analyses were conducted on four series of HDAC inhibitors. The studies were performed using the GRID/GOLPE combination using structure-based alignment. Twelve 3-D QSAR models were derived and discussed. Compared to previous studies on similar inhibitors, the present 3-D QSAR investigation proved to be of higher statistical value, displaying for the best global model r2, q2, and cross-validated SDEP values of 0.94, 0.83, and 0.41, respectively. A comparison of the 3-D QSAR maps with the structural features of the binding site showed good correlation. The results of 3D-QSAR and docking studies validated each other and provided insight into the structural requirements for anti-HDAC activity. To our knowledge this is the first 3-D QSAR application on a broad molecular diversity training set of HDACIs.     

毒蕈碱受体激动剂的三维定量构效关系研究

采用比较分子场分析法(CoMFA)研究了55个四氢吡啶类毒蕈碱受体激动剂的三维定量构效关系(3D-QSAR), 建立了具有较强预测能力的3D-QSAR模型. 所得模型的交叉验证相关系数(q²)为0.507, 常规相关系数(R²)为0.982 , 标准方差为0.218, 说明系列化合物分子周围立体场和静电场的分布与生物活性间存在良好的相关性. 模型不仅很好地预测了训练集和测试集化合物的活性, 而且为设计活性更高的受体激动剂提供了理论依据.     

External validation and prediction employing the predictive squared correlation coefficient test set activity mean vs training set activity mean

The external prediction capability of quantitative structure-activity relationship (QSAR) models is often quantified using the predictive squared correlation coefficient, q (2). This index relates the predictive residual sum of squares, PRESS, to the activity sum of squares, SS, without postprocessing of the model output, the latter of which is automatically done when calculating the conventional squared correlation coefficient, r (2). According to the current OECD guidelines, q (2) for external validation should be calculated with SS referring to the training set activity mean. Our present findings including a mathematical proof demonstrate that this approach yields a systematic overestimation of the prediction capability that is triggered by the difference between the training and test set activity means. Example calculations with three regression models and data sets taken from literature show further that for external test sets, q (2) based on the training set activity mean may become even larger than r (2). As a consequence, we suggest to always use the test set activity mean when quantifying the external prediction capability through q (2) and to revise the respective OECD guidance document accordingly. The discussion includes a comparison between r (2) and q (2) value ranges and the q (2) statistics for cross-validation.     

Eigenvalue analysis of peroxisome proliferator-activated receptor gamma agonists

Eigenvalue analysis (EVA) was conducted on a series of potent agonists of peroxisome proliferator-activated receptor gamma (PPARgamma). Predictive EVA quantitative structure-activity relationship (QSAR) models were established using the SYBYL package, which had conventional r2 and cross-validated coefficient (q2) values up to 0.920 and 0.587 for the AM1 method and 0.863 and 0.586 for the PM3 method, respectively. These models were validated by a test set containing 18 compounds. The capability to predict by these two models for PPARgamma agonists, with the best predictive r2pred value of 0.614 for AM1 and 0.822 for PM3 methods, set a successful example for applying a similar approach in building QSAR models for PPARalpha and -delta that could potentially offer a new opportunity in the design of novel PPAR modulators.     

3D-QSAR study on heterocyclic topoisomerase II inhibitors using CoMSIA

Selective topoisomerase II (Topo II) inhibitors have interested to a great extent for the design of new antitumoral compounds in recent years. Comparative molecular similarity indices analysis (CoMSIA) was performed on a series of previously synthesized benzoxazole, benzimidazole, and oxazolo(4,5-b)pyridine derivatives as eukaryotic Topo II inhibitors. A training set of 16 heterocyclic compounds was used to establish the CoMSIA model. They were constructed and geometrically optimized using SYBYL v7.0. The predictive ability of the model was assessed using a test set of 7 compounds. The best model has demonstrated a good fit having r2 value of 0.968 and cross-validated coefficient q2 value as 0.562 including steric and hydrophobic fields. The hydrophobic interactions showed a dominant role for increasing Topo II inhibitor activity and hydrophilic substituent was found more important than hydrophobic one on the 5 or 6 position of benzazole moiety. The model obtained from the present study can be useful for the modification and/or evaluation of the development of new Topo II inhibitors as potential antitumor compounds.     

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.     

(13)C NMR quantitative spectrometric data-activity relationship (QSDAR) models of steroids binding the aromatase enzyme

Five quantitative spectroscopic data-activity relationships (QSDAR) models for 50 steroidal inhibitors binding to aromatase enzyme have been developed based on simulated (13)C nuclear magnetic resonance (NMR) data. Three of the models were based on comparative spectral analysis (CoSA), and the two other models were based on comparative structurally assigned spectral analysis (CoSASA). A CoSA QSDAR model based on five principal components had an explained variance (r(2)) of 0.78 and a leave-one-out (LOO) cross-validated variance (q(2)) of 0.71. A CoSASA model that used the assigned (13)C NMR chemical shifts from a steroidal backbone at five selected positions gave an r(2) of 0.75 and a q(2) of 0.66. The (13)C NMR chemical shifts from atoms in the steroid template position 9, 6, 3, and 7 each had correlations greater than 0.6 with the relative binding activity to the aromatase enzyme. All five QSDAR models had explained and cross-validated variances that were better than the explained and cross-validated variances from a five structural parameter quantitative structure-activity relationship (QSAR) model of the same compounds. QSAR modeling suffers from errors introduced by the assumptions and approximations used in partial charges, dielectric constants, and the molecular alignment process of one structural conformation. One postulated reason that the variances of QSDAR models are better than the QSAR models is that (13)C NMR spectral data, based on quantum mechanical principles, are more reflective of binding than the QSAR model's calculated electrostatic potentials and molecular alignment process. The QSDAR models provide a rapid, simple way to model the steroid inhibitor activity in relation to the aromatase enzyme.     

Three-dimensional quantitative structure-activity relationship of nucleosides acting at the A3 adenosine receptor: analysis of binding and relative efficacy

The binding affinity and relative maximal efficacy of human A3 adenosine receptor (AR) agonists were each subjected to ligand-based three-dimensional quantitative structure-activity relationship analysis. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) used as training sets a series of 91 structurally diverse adenosine analogues with modifications at the N6 and C2 positions of the adenine ring and at the 3', 4', and 5' positions of the ribose moiety. The CoMFA and CoMSIA models yielded significant cross-validated q2 values of 0.53 (r2 = 0.92) and 0.59 (r2 = 0.92), respectively, and were further validated by an external test set (25 adenosine derivatives), resulting in the best predictive r2 values of 0.84 and 0.70 in each model. Both the CoMFA and the CoMSIA maps for steric or hydrophobic, electrostatic, and hydrogen-bonding interactions well reflected the nature of the putative binding site previously obtained by molecular docking. A conformationally restricted bulky group at the N6 or C2 position of the adenine ring and a hydrophilic and/or H-bonding group at the 5' position were predicted to increase A3AR binding affinity. A small hydrophobic group at N6 promotes receptor activation. A hydrophilic and hydrogen-bonding moiety at the 5' position appears to contribute to the receptor activation process, associated with the conformational change of transmembrane domains 5, 6, and 7. The 3D-CoMFA/CoMSIA model correlates well with previous receptor-docking results, current data of A3AR agonists, and the successful conversion of the A3AR agonist into antagonists by substitution (at N6) or conformational constraint (at 5'-N-methyluronamide).     

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.     

A new procedure for improving the predictiveness of CoMFA models and its application to a set of dihydrofolate reductase inhibitors

Summary A new automated procedure to improve the predictive quality of CoMFA models for both training and test sets is described. A model of greater consistency is generated by performing small reorientations of the underlying molecules for which too low activities are calculated. In order to predict activities of test compounds, the most similar molecules in the previously optimized model are identified and used as a basis for the prediction. This method has been applied to two independent sets of dihydrofolate reductase inhibitors (80 compounds each, serving as training sets), resulting in a significant increase of the cross-validated r² value. For both models, the predictive r² value for a test set consisting of 70 compounds was improved substantially.     

Molecular modeling of the intestinal bile acid carrier: A comparative molecular field analysis study

A structure–binding activity relationship for the intestinal bile acidtransporter has been developed using data from a series of bile acid analogsin a comparative molecular field analysis (CoMFA). The studied compoundsconsisted of a series of bile acid–peptide conjugates, withmodifications at the 24 position of the cholic acid sterol nucleus, andcompounds with slight modifications at the 3, 7, and 12 positions. For theCoMFA study, these compounds were divided into a training set and a test set,comprising 25 and 5 molecules, respectively. The best three-dimensionalquantitative structure–activity relationship model found rationalizesthe steric and electrostatic factors which modulate affinity to the bile acidcarrier with a cross-validated, conventional and predictive r²of 0.63, 0.96, and 0.69, respectively, indicating a good predictive model forcarrier affinity. Binding is facilitated by positioning an electronegativemoiety at the 24–27 position, and also by steric bulk at the end of theside chain. The model suggests substitutions at positions 3, 7, 12, and 24that could lead to new substrates with reasonable affinity for the carrier.     

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.     

The comparative molecular surface analysis (CoMSA) with modified uniformative variable elimination-PLS (UVE-PLS) method: application to the steroids binding the aromatase enzyme

The application of the CoMSA method to analyze 3D QSAR of 50 steroid aromatase inhibitors is described. The 3D QSAR model obtained, reaching a value of cross-validated q(2) = 0.96 (s = 0.31), significantly outperforms those reported in the literature for the CoMFA or CoSA (CoSASA). It is shown that the Uniformative Variable Elimination UVE-PLS or modified iterative UVE procedure (IVE-PLS) can be used for indicating the regions contributing to the binding activity. Thus, after separating the series into two groups of the training and test molecules quite correct external predictions result from the processing of the training set. We proved that the procedure of the data elimination provides stable results, if tested in 50 random runs of the IVE-PLS-CoMSA with different training/test sets. Depending upon the procedure used the quality of the predictions for 25 test molecules is given by SDEP = sum(y(pred)-y(obs))(2)/n)(1/2) = 0.321 - 0.782.     

3-氨磺酰苯甲酸类AKR1C3抑制剂的3D-QSAR和分子对接研究

醛酮还原酶1C3(AKR1C3)作为治疗前列腺癌的新靶点已成为研究热点,3-氨磺酰苯甲酸衍生物对其具有高效的选择性和抑制活性。本文采用比较分子场分析(COMFA)和比较分子相似性指数分析(COMSIA)方法,将经分子对接后的34个优势构象组成训练集和11个优势构象组成测试集,构建三维定量构效关系(3D-QSAR)模型。COMFA模型的交叉验证系数(q2),非交叉验证系数(R2),标准偏差(SEE)和F值分别为0.761,0.973,0.122,185.963;自举法回归系数为R2bs=0.98。最佳组合COMSIA模型的q2,R2,SEE,F和R2bs分别为0.734,0.984,0.097,147.850,0.994。COMFA和COMSIA模型的系统外部测试R2pred分别为0.864和0.756,r2m分别为0.8127和0.5377。这些结果表明,所建立的QSAR模型具有较高的可靠性和较强预测能力。经三维等势图分析可知,在2、5或6位适当增加取代基体积,或在5位引入氢键受体,或在7位引入负电性取代基则能提高化合物的生物活性。该模型为进一步设计具有更优选择性和活性的化合物提供了理论依据。