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.     

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.     

3D-QSAR and molecular docking studies of 1,3,5-triazene-2,4-diamine derivatives against r-RNA: novel bacterial translation inhibitors

Aminoglycoside mimetics inhibit bacterial translation by interfering with the ribosomal decoding site. To elucidate the structural properties of these compounds important for antibacterial activity, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were applied to a set of 56 aminoglycosides mimetics. The successful CoMFA model yielded the leave-one-out (LOO) cross-validated correlation coefficient (q(2)) of 0.708 and a non-cross-validated correlation coefficient (r(2)) of 0.967. CoMSIA model gave q(2)=0.556 and r(2)=0.935. The CoMFA and CoMSIA models were validated with 36 test set compounds and showed a good r(pred)(2) of 0.624 and 0.640, respectively. Contour maps of the two QSAR approaches show that electronic effects dominantly determine the binding affinities. These obtained results were agreed well with the experimental observations and docking studies. The results not only lead to a better understanding of structural requirements of bacterial translation inhibitors but also can help in the design of novel bacterial translation inhibitors.     

3D-QSAR and molecular docking studies of selective agonists for the thyroid hormone receptor beta

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) on a series of agonists of thyroid hormone receptor beta (TRbeta), which may lead to safe therapies for non-thyroid disorders while avoiding the cardiac side effects. The reasonable q(2) (cross-validated) values 0.600 and 0.616 and non-cross-validated r(2) values of 0.974 and 0.974 were obtained for CoMFA and CoMSIA models for the training set compounds, respectively. The predictive ability of two models was validated using a test set of 12 molecules which gave predictive correlation coefficients (r(pred)(2)) of 0.688 and 0.674, respectively. The Lamarckian Genetic Algorithm (LGA) of AutoDock 4.0 was employed to explore the binding mode of the compound at the active site of TRbeta. The results not only lead to a better understanding of interactions between these agonists and the thyroid hormone receptor beta but also can provide us some useful information about the influence of structures on the activity which will be very useful for designing some new agonist with desired activity.     

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.     

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.     

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.     

Structural and chemical basis for enhanced affinity and potency for a large series of estrogen receptor ligands: 2D and 3D QSAR studies

The estrogen receptor (ER) is an important drug target for the development of novel therapeutic agents for the treatment of breast cancer. Progress towards the design of more potent and selective ER modulators requires the optimization of multiple ligand-receptor interactions. Comparative molecular field analyses (CoMFA) and hologram quantitative structure-activity relationships (HQSAR) were conducted on a large set of ERalpha modulators. Two training sets containing either 127 or 69 compounds were used to generate QSAR models for in vitro binding affinity and potency, respectively. Significant correlation coefficients (affinity models, CoMFA, r(2)=0.93 and q(2)=0.79; HQSAR, r(2)=0.92 and q(2)=0.71; potency models, CoMFA, r(2)=0.94 and q(2)=0.72; HQSAR, r(2)=0.92 and q(2)=0.74) were obtained, indicating the potential of the models for untested compounds. The generated models were validated using external test sets, and the predicted values were in good agreement with the experimental results. The final QSAR models as well as the information gathered from 3D contour maps should be useful for the design of novel ERalpha modulators having improved affinity and potency.     

Three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses of choline acetyltransferase inhibitors

As a basis for predicting structural features that may lead to the design of more potent and selective inhibitors of choline acetyltransferase (ChAT), the three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were carried out on a series of trans-1-methyl-4-(1-naphthylvinyl)pyridinium (MNVP+) analogs, which are known ChAT inhibitors. 3D-QSAR studies were carried out using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. Since these inhibitors have extremely shallow potential energy minimum energy wells and low barriers to rotation, two dihedral angles unique to these inhibitors were systematically modified to reflect the energetically preferred conformations as determined by force field calculations. An optimum alignment rule was devised based on the conformations obtained from the molecular mechanics studies, using a common substructure alignment method. The studies involve a set of 21 compounds and experimentally determined molar IC50 values were used as the dependent variable in the analysis. The 3D-QSAR models have conventional r2-values of 0.953 and 0.954 for CoMFA and CoMSIA, respectively; similarly, cross-validated coefficient q2-values of 0.755 and 0.834 for CoMFA and CoMSIA, respectively, were obtained. On the basis of these predictive r2-values the model was tested using previously determined IC50 values. CoMSIA 3D-QSAR yielded better results than CoMFA.     

A CoMFA study of COX-2 inhibitors with receptor based alignment

A diverse set of 53 cyclooxygenase-2 (COX-2) inhibitors which were aligned in two different ways were subjected to CoMFA analysis. The first method of alignment of the molecules was based on the binding information sourced from the crystallographic study, from which CoMFA Model 1 was derived. The second mode of alignment was generated by docking the inhibitors in the binding pocket using the DOCK and AFFINITY suite of programs; this gave a second model. The CoMFA Model 2 was slightly better than Model 1 in terms of the statistical parameters r(2) and q(2). The two models could predict very well the activity of a test set of diverse molecules, with a predictive r(2) of 0.593 and 0.768, respectively. Besides the QSAR results, the docking studies give a deep insight into the H-bonding interactions between the inhibitors and residues in the active site of the enzyme, which can be exploited in designing better inhibitors. Useful ideas on activity improvement could be gleaned from these models.     

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研究

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 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.     

Three-dimensional quantitative structure-activity relationship for several bioactive peptides searched by a convex hull-comparative molecular field analysis approach

Three-dimensional (3D) convex hulls are computed for theoretically generated structures of a group of 18 bioactive tachykinin peptides. The number of peptides treated as a training set is 14, whereas that treated as a test set is four. The frequency of atoms of the same atomic type lying at the vertices of all the hulls computed for all the structures in a structural set is counted. Vertex atoms with non-zero frequency counted are collected together as a set of commonly exposed groups. These commonly exposed atoms are then treated as a set of correspondences for aligning all the other 13 structures in a structural set against a common template, which is the structure of the most potent peptide in the set using the FIT module of the SYBYL 6.6 program. Each aligned structural set is then analyzed by the comparative molecular field analysis (CoMFA) module using the C.3 probe having a charge of +1.0. The corresponding cross-validated r2 values range from -0.99 to 0.57 for a number of 73 structural sets studied. The comparative molecular similarity indices analysis (CoMSIA) module within the SYBYL 6.6 package is also used to analyze some of these aligned structural sets. Although the CoMSIA results are in accord with those of CoMFA, it is also found that the CoMFA results of several structural sets can be improved somewhat for conformations of the structures in the sets that are adjusted by constraint energy minimization and then constraint molecular dynamics simulation runs using distance constraints derived from some commonly exposed groups determined for them. This result further implies that the convex hull-CoMFA is a feasible approach to screen the bioactive conformations for molecules of this class.     

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.     

Combining docking, scoring and molecular field analyses to probe influenza neuraminidase-ligand interactions

In this project, several docking conditions, scoring functions and corresponding protein-aligned molecular field analysis (CoMFA) models were evaluated for a diverse set of neuraminidase (NA) inhibitors. To this end, a group of inhibitors were docked into the active site of NA. The docked structures were utilized to construct a corresponding protein-aligned CoMFA models by employing probe-based (H+, OH, CH3) energy grids and genetic partial least squares (G/PLS) statistical analysis. A total of 16 different docking configurations were evaluated, of which some succeeded in producing self-consistent and predictive CoMFA models. However, the best model coincided with docking the ionized ligands into the hydrated form of the binding site via PLP1 scoring function (r2LOO=0.735, r2PRESS against 24 test compounds=0.828). The highest-ranking CoMFA models were employed to probe NA-ligand interactions. Further validation by comparison with a co-crystallized ligand-NA crystallographic structure was performed. This combination of docking/scoring/CoMFA modeling provided interesting insights into the binding of different NA inhibitors.