Insight into the structural requirements of proton pump inhibitors based on CoMFA and CoMSIA studies

In the present study, a series of 179 quinoline and quinazoline heterocyclic analogues exhibiting inhibitory activity against Gastric (H+/K+)-ATPase were investigated using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA) methods. Both the models exhibited good correlation between the calculated 3D-QSAR fields and the observed biological activity for the respective training set compounds. The most optimal CoMFA and CoMSIA models yielded significant leave-one-out cross-validation coefficient, q(2) of 0.777, 0.744 and conventional cross-validation coefficient, r(2) of 0.927, 0.914 respectively. The predictive ability of generated models was tested on a set of 52 compounds having broad range of activity. CoMFA and CoMSIA yielded predicted activities for test set compounds with r(pred)(2) of 0.893 and 0.917 respectively. These validation tests not only revealed the robustness of the models but also demonstrated that for our models r(pred)(2) based on the mean activity of test set compounds can accurately estimate external predictivity. The factors affecting activity were analyzed carefully according to standard coefficient contour maps of steric, electrostatic, hydrophobic, acceptor and donor fields derived from the CoMFA and CoMSIA. These contour plots identified several key features which explain the wide range of activities. The results obtained from models offer important structural insight into designing novel peptic-ulcer inhibitors prior to their synthesis.     

An investigation of structurally diverse carbamates for acetylcholinesterase (AChE) inhibition using 3D-QSAR analysis

In order to identify the essential structural features and physicochemical properties for acetylcholinesterase (AChE) inhibitory activity in some carbamate derivatives, the systematic QSAR (Quantitative Structure Activity Relationship) studies (CoMFA, advance CoMFA and CoMSIA) have been carried out on a series of (total 78 molecules) taking 52 and 26 molecules in training and test set, respectively. Statistically significant 3D-QSAR (three-dimensional Quantitative Structure Activity Relationship) models were developed on training set molecules using CoMFA and CoMSIA and validated against test set compounds. The highly predictive models (CoMFA q(2)=0.733, r(2)=0.967, predictive r(2)=0.732, CoMSIA q(2)=0.641, r(2)=0.936, predictive r(2)=0.812) well explained the variance in binding affinities both for the training and the test set compounds. The generated models suggest that steric, electrostatic and hydrophobic interactions play an important role in describing the variation in binding affinity. In particular the carbamoyl nitrogen should be more electropositive; substitutions on this nitrogen should have high steric bulk and hydrophobicity while the amino nitrogen should be electronegative in order to have better activity. These studies may provide important insights into structural variations leading to the development of novel AChE inhibitors which may be useful in the development of novel molecules for the treatment of Alzheimer's disease.     

黄酮类化合物的3D-QSAR研究

从NCI数据库中,筛选出67个与矢车菊黄素类似的天然黄酮化合物.采用CoMFA方法研究其构效关系,构建CoMFA模型,其模型相关系数为q2=0.599,r2=0.919,验证模型的预测能力和拟合能力较好.通过分子场等势图,可直观分子周围立体和静电特征对化合物活性的影响,为设计高活性黄酮衍生物提供理论依据.     

N,N-二甲基-2-溴苯乙胺类衍生物活性与CoMFA、CoMSIA的研究

用比较分子力场分析(CoMFA)法和比较分子相似性指数分析(CoMSIA)法,建立N,N-二甲基-2-溴苯乙胺类化合物的3D-QSAR模型。CoMFA模型中,其交叉验证系数q2=0．792,传统的相关系数R2=0．955(R=0．978),相应立体场贡献为77．4％、静电场贡献为22．6％,优于文献的报导。CoMSIA研究中,其交叉验证系数q2=0．757,传统的相关系数R2=0．917 (R=0．958),其疏水场、立体场、静电场贡献依次为:42．9％、39．5％、17．6％。用两种模型分别预测检测集分子的活性,结果与实验值较吻合。说明所建的模型具有较好的预测能力。通过分析CoMFA分子场等值线图,可为优化N,N-二甲基-2-溴苯乙胺类衍生物的结构提供理论指导。     

类青蒿素抗疟活性3D-QSAR研究

采用比较分子力场分析(CoMFA)法,将40个类青蒿素的抗疟活性作了定量构效关系的研究。所得模型交叉验证系数q2为0．601,非交叉验证系数r2为0．982,标准偏差SE=0．140,F=302．246,其中立体场与静电场的贡献分别为33．0％和67．0％。CoMFA离散图表明,增大C13上取代基的体积以及增强O15上取代基的电负性等,都有利于提高化合物的抗疟活性。该模型的离散图为改造此类化合物的结构提供理论依据和研究方向。     

甾类化合物的结构与抗炎活性相关性3D—QSAR研究

研究了甾族化合物的结构与抗炎活性相关的三维构效关系，采用比较分子力场分析法寻找结构与活性的关系，建立了一个可靠的构效关系模型，为活性分子的进一步修和改造提供了信息。     

Investigation of structural requirements of anticancer activity at the paclitaxel/tubulin binding site using CoMFA and CoMSIA

CoMFA and CoMSIA analysis were utilized in this investigation to define the important interacting regions in paclitaxel/tubulin binding site and to develop selective paclitaxel-like active compounds. The starting geometry of paclitaxel analogs was taken from the crystal structure of docetaxel. A total of 28 derivatives of paclitaxel were divided into two groups-a training set comprising of 19 compounds and a test set comprising of nine compounds. They were constructed and geometrically optimized using SYBYL v6.6. CoMFA studies provided a good predictability (q(2)=0.699, r(2)=0.991, PC=6, S.E.E.=0.343 and F=185.910). They showed the steric and electrostatic properties as the major interacting forces whilst the lipophilic property contribution was a minor factor for recognition forces of the binding site. These results were in agreement with the experimental data of the binding activities of these compounds. Five fields in CoMSIA analysis (steric, electrostatic, hydrophobic, hydrogen-bond acceptor and donor properties) were considered contributors in the ligand-receptor interactions. The results obtained from the CoMSIA studies were: q(2)=0.535, r(2)=0.983, PC=5, S.E.E.=0.452 and F=127.884. The data obtained from both CoMFA and CoMSIA studies were interpreted with respect to the paclitaxel/tubulin binding site. This intuitively suggested where the most significant anchoring points for binding affinity are located. This information could be used for the development of new compounds having paclitaxel-like activity with new chemical entities to overcome the existing pharmaceutical barriers and the economical problem associated with the synthesis of the paclitaxel analogs. These will boost the wide use of this useful class of compounds, i.e. in brain tumors as the most of the present active compounds have poor blood-brain barrier crossing ratios and also, various tubulin isotypes has shown resistance to taxanes and other antimitotic agents.     

3D-QSAR and molecular docking studies on pyrazolopyrimidine derivatives as glycogen synthase kinase-3beta inhibitors

Glycogen synthase kinase-3 (GSK-3), a serine/threonine kinase, is a fascinating enzyme with diverse biological actions in intracellular signaling systems, making it an emerging target for diseases such as diabetes mellitus, cancer, chronic inflammation, bipolar disorders and Alzheimer's disease. It is important to inhibit GSK-3 selectively and the net effect of the GSK-3 inhibitors thus should be target specific, over other phylogenetically related kinases such as CDK-2. In the present work, we have carried out three-dimensional quantitative structure activity relationship (3D-QSAR) studies on novel class of pyrazolopyrimidine derivatives as GSK-3 inhibitors reported to have improved cellular activity. Docked conformation of the most active molecule in the series, which shows desirable interactions in the receptor, was taken as template for alignment of the molecules. Statistically significant CoMFA and CoMSIA models were generated using 49 molecules in training set. By applying leave-one-out (LOO) cross-validation study, r(cv)2 values of 0.53 and 0.48 for CoMFA and CoMSIA, respectively and non-cross-validated (r(ncv)2) values of 0.98 and 0.92 were obtained for CoMFA and CoMSIA models, respectively. The predictive ability of CoMFA and CoMSIA models was determined using a test set of 12 molecules which gave predictive correlation coefficients (r(pred)2) of 0.47 and 0.48, respectively, indicating good predictive power. Based upon the information derived from CoMFA and CoMSIA contour maps, we have identified some key features that explain the observed variance in the activity and have been used to design new pyrazolopyrimidine derivatives. The designed molecules showed better binding affinity in terms of estimated docking scores with respect to the already reported systems; hence suggesting that newly designed molecules can be more potent and selective towards GSK-3beta inhibition.     

Three-dimensional QSAR of HPPD inhibitors, PSA inhibitors, and anxiolytic agents: effect of tautomerism on the CoMFA models

The present study was design to examine the effect of tautomerism upon the CoMFA results. Three selected data sets involving protropic tautomerism, which are 21 p-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitors, 35 inhibitors of puromycin-sensitive aminopeptidase (PSA), and 67 anxiolytic agents, were used for this purpose. Atom-by-atom alignment technique was adopted to superimpose the molecules in the data sets onto a template. The structural alignments using different tautomeric forms had no significant difference except the atoms involved in tautomerism, which ensures, to a great extent, that the differences of the CoMFA results result primarily from the tautomerism. All-orientation and all-placement search (AOS-APS) based CoMFA models, in addition to the conventional ones, were derived for each system and proved to be capable of yielding much improved statistical results. In the cases of the data sets of HPPD inhibitors and PSA inhibitors, excellent AOS-APS CoMFA models (q2>0.8 with four components for the former and q2>0.7 with seven components for the latter) were obtained, and almost no significant difference in statistical quality was observed when using different tautomeric forms to derive the models. However, it was not the case when treating the data set of anxiolytic agents. The keto tautomer, which was the active form of the PBI type inhibitors, produced measurably better results (q2=0.54 with eight components) than that the enol one (q2=0.37 with five components), indicating the importance of selecting proper tautomer in the CoMFA studies. Furthermore, there existed some substantial differences of the electrostatic field contours between the two different tautomeric forms for all of the three systems considered, whereas the differences in the steric field contour maps were limited. This implies that the resulting new potent ligands may be quite different if one utilizes the CoMFA models of different tautomeric forms for guiding further structural refinements.     

新型肟醚拟除虫菊酯化合物杀粘虫活性3D-QSAR研究

用比较分子力场分析（CoMFA）方法对50个自主设计并合成的新型芳基烷基肟醚化合物杀粘虫活性的定量构效关系进行了研究。所得模型交叉验证参数q^2为0．598，非交叉验证相关系数r^2为0．949，标准偏差SE=0．136，F=92．419，影响药效的立体场与静电场的贡献分别50．6％和49．4％。CoMFA等值线图表明化合物肟苯环对位取代基是影响其杀粘虫活性的主要因素，体积和电性增加有利于提高化合物的杀粘虫活性。该模型三维等值线图为该类化合物的结构改造提供了理论依据。     

真蛸胺类杀虫剂药效团模型的确定

用比较分子场分析方法（ＣｏＭＦＡ）研究了２１种真蛸胺类杀虫剂的三维定量构效关系，确定了药效团模型，得到了具有较好预测能力的ＣｏＭＦＡ模型。研究表明，苯环与支链Ｎ原子及它们间的距离对于这些化合物的活性有重要影响，而且立体与静电场对其药效均起作用。     

The comparative molecular surface analysis (COMSA): a novel tool for molecular design

A new method allowing for 3-D QSAR analysis and the prediction of biological activity is presented. Unlike comparative molecular field analysis (CoMFA)-like techniques, it is based not on a comparison of the properties characterizing a discrete set of points but on the mean electrostatic potential (MEP) calculated and labeling specific areas defined on the molecular surface. A Kohonen self-organizing neural network and partial least square (PLS) analysis have been used for performing such an operation. The series of steroids complexing the corticosteroid (CBG) and testosterone (TBG) globulins, which forms a benchmark measuring the performance of the methods in molecular design, and a series of benzoic acids described by the Hammett sigma constants is used for testing the method. It is demonstrated that a method can be used efficiently to evaluate the responses determined both by the combination of electrostatic and steric effects or by electrostatic effects alone, therefore, two different schemes were developed. The first one, which involves PLS analysis of the full comparative networks, covers both steric and electrostatic effects. This scheme works well for both the CBG and TBG data. The second scheme takes into account only the properties (MEP) of these regions within molecules that can be superimposed with the template molecule. This scheme provides the best predictive power for the benzoic acids series. Comparison of the results from a CoMFA analysis proves that method is at least as effective for the responses limited by electrostatic effects, although it significantly outperforms CoMFA for CBG affinity which is dominated by steric effects.     

氰基丙烯酸酯类化合物的3D-QSAR研究

3D-Quantitative Structure—Activity Relationships of a series of cyanoacrylate inhibitors which are known to act by blocking photosynthetic electron transport close to photosystem II reaction center in the thylakoid membranes of plant chloroplasts have been investigated using comparative molecular field analysis(CoMFA). The active conformations of the title compounds were obtained by constrained systematic search program. The resulting CoMFA model with considerable predictive ability shows that the contribution of steric effects for activities of cyanoacrylate inhibitors is more important than electrostatic effects.     

3D QSAR studies of 1,3,4-benzotriazepine derivatives as CCK2 receptor antagonists

A number of CCK(2) antagonists have been reported to play an important role in controlling gastric acid-related conditions, nervous system related disorders and certain types of cancer. To obtain the helpful information for designing potent antagonists with novel structures and to investigate the quantitative structure-activity relationship of a group of 62 different CCK(2) receptor antagonists with varying structures and potencies, CoMFA, CoMSIA, and HQSAR studies were carried out on a series of 1,3,4-benzotriazepine-based CCK(2) receptor antagonists. QSAR models were derived from a training set of 49 compounds. By applying leave-one-out (LOO) cross-validation study, cross-validated (r(cv)(2)) values of 0.673 and 0.608 and non-cross-validated (r(ncv)(2)) values of 0.966 and 0.969 were obtained for the CoMFA and CoMSIA models, respectively. The predictive ability of the CoMFA and CoMSIA models was determined using a test set of 13 compounds, which gave predictive correlation coefficients (r(pred)(2)) of 0.793 and 0.786, respectively. HQSAR was also carried out as a complementary study, and the best HQSAR model was generated using atoms, bonds, hydrogen atoms, and chirality as fragment distinction with fragment size (2-5) and six components showing r(cv)(2) and r(ncv)(2) values of 0.744 and 0.918, respectively. CoMFA steric and electrostatic, CoMSIA hydrophobic and hydrogen bond acceptor fields, and HQSAR atomic contribution maps were used to analyze the structural features of the datasets that govern their antagonistic potency.     

Structural investigation of PAP derivatives by CoMFA and CoMSIA reveals novel insight towards inhibition of Bcr-Abl oncoprotein

Molecular modeling by 3D-QSAR comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were employed on a series of phenylaminopyrimidine-based (PAP) Bcr-Abl inhibitors. The chemical structures of 63 PAP analogues were aligned using a template extracted from the crystal structure of STI571 bound to Abl kinase. Subsequently, the structures built were divided into training and test sets that include 53 and 10 compounds, respectively. Statistical results showed that the 3D-QSAR models generated from CoMSIA were superior to CoMFA (CoMSIA; q2=0.66, r2=0.94, N=3, F=139.09, r2pred=0.64 while CoMFA; q2=0.53, r2=0.73, N=3, F=43.53, r2pred=0.61). Based on the contour interpretation, the attachment of hydrophobic and bulky groups to the phenyl and pyrrolidine (D- and E-ring of NS-187, respectively) along with highly electronegative groups around the D-ring are important structural features for the design of second-generation Bcr-Abl inhibitors. The generated models are predictive based on reproducible values of the predicted compared with experimental activities in the test set. Further, the complementary analysis of contour maps to the Bcr-Abl binding site suggested the anchor points for binding affinity.     

3D-QSAR studies of D3R antagonists and 5-HT1AR agonists

Combination of dopamine D₃ antagonism and serotonin 5-HT_1A agonism leads to an effective way to atypical antipsychotics. In this work, two predictive 3D-QSAR models were bulit for D₃R antagonists and 5-HT_1AR agonists, respectively. Based on the steric and electrostatic information of contour maps, four compounds with improved predicted activities were newly designed. In addition, molecular docking and ADMET properties suggested that designed molecules had strong interactions with receptors and low hepatotoxicity. This work sheds light on the design of bifunctional novel antipsychotic drugs for D₃R antagonists and 5HT_1AR agonists.     

3D QSAR studies on a series of potent and high selective inhibitors for three kinases of RTK family

For targets belonging to the same family of receptors, inhibitors often act at more than one biological target and produce a synergistic effect. Separate CoMFA and CoMSIA models were developed from our data set for the KDR, cKit and Tie-2 inhibitors. These models showed excellent internal predictability and consistency, and validation using test-set compounds yielded a good predictive power for the pIC(50) value. The field contour maps (CoMFA and CoMSIA) corresponding to the KDR, cKit and Tie-2 kinase subtypes reflected the characteristic similarities and differences between these types. These maps provided valuable information to facilitate structural modifications of the inhibitor to increase selectivity for the KDR over cKit and Tie-2.