Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.     

Development of novel 3D-QSAR combination approach for screening and optimizing B-Raf inhibitors in silico

B-Raf is a member of the RAF family of serine/threonine kinases: it mediates cell division, differentiation, and apoptosis signals through the RAS-RAF-MAPK pathway. Thus, B-Raf is of keen interest in cancer therapy, such as melanoma. In this study, we propose the first combination approach to integrate the pharmacophore (PhModel), CoMFA, and CoMSIA models for B-Raf, and this approach could be used for screening and optimizing potential B-Raf inhibitors in silico. Ten PhModels were generated based on the HypoGen BEST algorithm with the flexible fit method and diverse inhibitor structures. Each PhModel was designated to the alignment rule and screening interface for CoMFA and CoMSIA models. Therefore, CoMFA and CoMSIA models could align and recognize diverse inhibitor structures. We used two quality validation methods to test the predication accuracy of these combination models. In the previously proposed combination approaches, they have a common factor in that the number of training set inhibitors is greater than that of testing set inhibitors. In our study, the 189 known diverse series B-Raf inhibitors, which are 7-fold the number of training set inhibitors, were used as a testing set in the partial least-squares validation. The best validation results were made by the CoMFA09 and CoMSIA09 models based on the Hypo09 alignment model. The predictive r(2)(pred) values of 0.56 and 0.56 were derived from the CoMFA09 and CoMSIA09 models, respectively. The CoMFA09 and CoMSIA09 models also had a satisfied predication accuracy of 77.78% and 80%, and the goodness of hit test score of 0.675 and 0.699, respectively. These results indicate that our combination approach could effectively identify diverse B-Raf inhibitors and predict the activity.     

Pharmacophore and Docking-based 3D-QSAR Studies on HIV-1 Integrase Inhibitors

Integrase（IN） plays an essential role in the process of HIV-1 replication.IN inhibitors of diketo acid derivatives（DKAs） were analysed by the Comparative Molecular Field Analysis（CoMFA） and Comparative Molecular Similarity Induces Analysis（CoMSIA） methods.A set of 42 compounds were randomly selected as the training set（35） and test set（7）.Firstly,a good pharmacophore（goodness of hit=0.787） was obtained and used to align ligands.Then,predictive models were constructed with the CoMFA and CoMSIA methods based on the pharmacophore alignment.As a result,the CoMS1A method yielded the best model with an r2 of 0.955 and a q2 of 0.665,which can predict the activities of the tested DKAs very well（r2=0.559）.Finally,DKAs were docked into IN,and the predicit modes were superimposed on the contour maps obtained from the best CoMSIA model.The superimposed maps gave a visualized and meaningful insight into the inhibitory behaviors,providing significantly useful information for the rational drug design of anti-IN agents.     

3D QSAR analyses of novel tyrosine kinase inhibitors based on pharmacophore alignment

In an effort to develop a quantitative ligand-binding model for the receptor tyrosine kinases, a pharmacophore search was first used to identify structural features that are common in two novel sets of 12 molecules of the 3-substituted indolin-2-ones and 19 compounds of the benzylidene malononitriles with low-to-high affinity for HER2, a kind of receptor tyrosine kinase. The common pharmacophore model based on these 31 compounds was used as a template to obtain the aligned molecular aggregate, which provided a good starting point for 3D-QSAR analysis of only the 19 benzylidene malononitriles. Two molecular field analysis (MFA) techniques, including CoMFA and CoMSIA, were used to derive the quantitative structure-activity relationships of the studied molecules. From the studied results, it was obvious that the 3D-QSAR models based on the pharmacophore alignment were superior to those based on the simple atom-by-atom fits. Considering the flexibility of the studied molecules and the difference between the active conformers and the energy-lowest conformers, the pharmacophore model can usually provide the common features for the flexible regions. Moreover, the best CoMSIA model based on the pharmacophore hypothesis gave good statistical measure from partial least-squares analysis (PLS) (q(2) = 0.71), which was slightly better than the CoMFA one. Our study demonstrated that pharmacophore modeling and CoMSIA research could be effectively combined. Results obtained from both methods helped with understanding the specific activity of some compounds and designing new specific HER2 inhibitors.     

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.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

Pharmacophore elucidation and 3D-QSAR analysis of a new class of highly potent inhibitors of acid ceramidase based on maximum common substructure and field fit alignment methods

Comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA) studies have been carried on a series of 2,4-dioxopyrimidine-1-carboxamides as acid ceramidase inhibitors. Two alignment rules for the compounds were defined using maximum common substructure and field fit. The best orientation was then searched by all-orientation search strategy, to minimize the effect of the initial orientation of the structures. The Kennard Stone algorithm was used to divide the entire set into training (25 compounds) and test (7 compounds) sets. Pharmacophore model identification was also performed using DISCOtech algorithm and refinement was carried out using GASP, to highlight important structural features that could be responsible for the inhibitory activity. All constructed models showed appropriate statistical parameters in terms of q ² and r _pred ² . Based upon the information obtained from CoMFA, CoMSIA, and developed pharmacophore pattern, some key features that may be used to design new inhibitors for acid ceramidase have been identified.     

CoMFA and CoMSIA studies on C-aryl glucoside SGLT2 inhibitors as potential anti-diabetic agents

SGLT2 has become a target of therapeutic interest in diabetes research. CoMFA and CoMSIA studies were performed on C-aryl glucoside SGLT2 inhibitors (180 analogues) as potential anti-diabetic agents. Three different alignment strategies were used for the compounds. The best CoMFA and CoMSIA models were obtained by means of Distill rigid body alignment of training and test sets, and found statistically significant with cross-validated coefficients (q ²) of 0.602 and 0.618, respectively, and conventional coefficients (r ²) of 0.905 and 0.902, respectively. Both models were validated by a test set of 36 compounds giving satisfactory predicted correlation coefficients (r ² _pred) of 0.622 and 0.584 for CoMFA and CoMSIA models, respectively. A comparison was made with earlier 3D QSAR study on SGLT2 inhibitors, which shows that our 3D QSAR models are better than earlier models to predict good inhibitory activity. CoMFA and CoMSIA models generated in this work can provide useful information to design new compounds and helped in prediction of activity 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.     

Design of novel FLT-3 inhibitors based on dual-layer 3D-QSAR model and fragment-based compounds in silico

FMS-like tyrosine kinase 3 (FLT-3) is strongly correlated with acute myeloid leukemia, but no FLT-3-inhibitor cocomplex structure is available to assist the design of therapeutic inhibitors. Hence, we propose a dual-layer 3D-QSAR model for FLT-3 that integrates the pharmacophore, CoMFA, and CoMSIA. We then coupled the model with the fragment-based design strategy to identify novel FLT-3 inhibitors. In the first layer, the previously established model, Hypo02, was evaluated in terms of its correlation coefficient (r), RMS, cost difference, and configuration cost, with values of 0.930, 1.24, 106.45, and 16.44, respectively. Moreover, Fischer's cross-validation test of data generated by Hypo02 yielded a 98% confidence level, and the validation of the testing set yielded a best r value of 0.87. The features of Hypo02 were separated into two parts and then used to screen the MiniMaybridge fragment compound database. Nine novel FLT-3 inhibitors were generated in this layer. In the second layer, Hypo02 was subjected to an alignment rule to generate CoMFA- and CoMSIA-based models, for which the partial least-squares validation method was utilized. The values of q(2), r(2), and predictive r(2) were 0.58, 0.98, and 0.76, respectively, derived from the CoMFA model with steric and electrostatic fields. The CoMSIA model with five different fields yielded values of 0.54, 0.97, and 0.76 for q(2), r(2), and predictive r(2), respectively. The CoMFA and CoMSIA models were used to constrain 3D structures of the nine novel FLT-3 inhibitors. This dual-layer 3D-QSAR model constitutes a valuable tool to easily and quickly screen and optimize novel potential FLT-3 inhibitors for the treatment of acute myeloid leukemia.     

CBP/P300溴结构域联芳基类抑制剂三维定量构效关系研究

CREB结合蛋白(CBP)和与其高度同源的P300蛋白是组蛋白乙酰化酶的两个亚型,两者通过它们的溴结构域(bromodomain,BRD)与染色质结合,目前,CBP/P300已经成为人类在肿瘤靶点领域中的研究热点。本研究基于CBP/P300溴结构域联芳基类抑制剂建立三维定量构效关系,采用比较分子力场分析法(Co MFA)和比较分子相似性指数分析法(Co MSIA)分别建立35个已知活性抑制剂的3D-QSAR模型,以确定CBP/P300溴结构域联芳基类抑制剂分子结构与生物活性之间的定量关系。Co MFA和Co MSIA模型活性数据p IC50的预测值与实验值基本一致,说明这两个模型具有较高的预测能力和统计学意义。根据Co MFA和Co MSIA模型所提供的立体场、静电场、疏水场、氢键给体场、氢键供体场等信息提出了改善此类抑制剂活性的药物设计思路,为指导设计具有更高活性的新分子和预测更加有效的CBP/P300溴结构域抑制剂提供理论依据。     

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies

In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

新型三唑类抗真菌化合物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 系统研究了40个新型三唑类化合物抗真菌活性的三维定量构效关系. 在CoMFA研究中, 研究了两种药效构象对模型的影响, 并考察了网格点步长对统计结果的影响. 在CoMSIA研究中, 系统考察了各种分子场组合、网格点步长和衰减因子对模型统计结果的影响, 发现立体场、静电场、疏水场和氢键受体场的组合得到最佳模型. 所建立CoMFA和CoMSIA模型的交叉相关系数q²值分别为0.718和0.655, 并都具有较强的预测能力. CoMFA和CoMSIA模型的三维等值线图直观地解释了化合物的构效关系, 阐明了化合物结构中苯环上各位置取代基对抗真菌活性的影响, 为进一步结构优化提供了重要依据.     

Receptor-based modeling and 3D-QSAR for a quantitative production of the butyrylcholinesterase inhibitors based on genetic algorithm

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for a series of structurally related steroidal alkaloids as butyrylcholinesterase (BuChE) inhibitors. Docking studies were employed to position the inhibitors into the BuChE active site to determine the most probable binding mode. The strategy was to explore multiple inhibitor conformations in producing a more reliable 3D-QSAR model. These multiple conformations were derived using the FlexS program. The conformation selection step for CoMFA was done by genetic algorithm. The genetic algorithm based CoMFA approach was found to be the best. Both CoMFA and CoMSIA yielded significant cross-validated q(2) values of 0.701 and 0.627 and the r(2) values of 0.979 and 0.982, respectively. These statistically significant models were validated by a test set of five compounds. Comparison of CoMFA and CoMSIA contour maps helped to identify structural requirements for the inhibitors and serves as a basis for the design of the next generation of the inhibitor analogues. The results demonstrate that the combination of ligand-based and receptor-based modeling with use of a genetic algorithm is a powerful approach to build 3D-QSAR models. These data can be used for the lead optimization process with respect to inhibition enhancement which is important for the drug discovery and development for Alzheimer's disease.     

Three-dimensional Quantitative Structure-Activity Relationship Analyses of a Series of Butenolide ETA Antagonists

Two kinds of Three-dimensional Quantitative Structure-activity Relationship(3D-QSAR) methods,comparative molecular filed analysis(CoMFA) and comparative molecular similarity indices analysis (CoMSIA) ,were applied to analyze the structure-activity relationship of a series of 63 butenolide ETA selective antagonists with respect to their inhibition against human ETA receptor,The CoMFA and CoMSIA models were developed for the conceivable alignment of the molecules based on a template structure from the crystallized data.The statistical results from the initial orientation of the aligned molecules show that the 3D-QSAR model from CoMFA(q^2=0.543) is obviously superior to that from the conventional CoMSIA(q^2=0.407).In order to refine the model,all-space search (ASS) was applied to minimize the field sampling process.By rotating and translating the molecular aggregate within the grid systematically,all the possible samplings of the molecular fields were tested and subsequently the one with the highest q^2 was picked out .The comparison of the sensitivity of CoMFA and CoMSIA to different space orientation shows that the CoMFA q^2 values are more sensitive to the translations and rotations of the aligned molecules with respect to the lattice than those of CoMSIA.The best CoMFA model from ASS was further refined by the region focused technique.The high quality of the best model is indicated by the high corss-validated correlation and the prediction on the external test set.The CoMFA coefficient contour plots identify several key features that explain the wide range of activities,which may help us to design new effective ETA selective antagonists.     

3D-QSAR Analysis of a Series of 1,2,3-Triazole-chromenone Derivatives as an Acetylcholinesterase Inhibitor against Alzheimer's Disease

Chromenones have attracted much attention since they are excellent acetylcholinesterase inhibitor(AChEi). The 1,2,3-triazoles are multifunctional anti-acetylcholinesterase(AChE) agents. In this paper, we report the three-dimensional quantitative structure-activity relationship(3D-QSAR) study of 25 1,2,3-triazolechromenone derivatives based comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA). To construct CoMFA and CoMSIA models, the 25 active molecules were randomly divided into the training and test sets. The obtained cross-validation Q~2 of the CoMFA model, the coefficient of non-cross-validation R~2, and the test value F are 0.597, 0.994, and 396.726, respectively. The cross-validation Q~2 of the CoMSIA model, the coefficient of the non-cross-validation R~2, and the test value F are 0.721, 0.979, and 131.107, respectively. The predictive correlation coefficient(r_(pred)~2) is 0.728 for CoMFA and 0.805 for CoMSIA, which verifies that the model is predictable. Based on the potential maps of CoMFA and CoMSIA, a library containing a set of potent AChEi was designed. The inhibitory potential of the compounds in this library was found to be greater than the inhibitory potential of the most active compounds in the data set. The results obtained from this study laid the foundation for the development of effective drugs for AChEi.     

3D-QSAR Studies on a Series of Indoleamide Derivatives as Antiplasmodial Drugs

In this study, three-dimensional quantitative structure-activity relationship(3D-QSAR) was studied for the antiplasmodial activity of a series of novel indoleamide derivatives by comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(Co MSIA). 3D-QSAR model was established by a training set of 20 compounds and was externally validated by a test set of 4 compounds. The best prediction(Q~2 = 0.593 and 0.527, R~2 = 0.990 and 0.953, r_(pred)~2 = 0.967 and 0.962 for CoMFA and CoMSIA) was obtained according to CoMFA and CoMSIA. Those parameters indicated the model was reliable and predictable. We designed several molecules with high activities according to the contour maps produced by the CoMFA and CoMSIA models.     

Strategies to design pyrazolyl urea derivatives for p38 kinase inhibition: a molecular modeling study

The p38 protein kinase is a serine–threonine mitogen activated protein kinase, which plays an important role in inflammation and arthritis. A combined study of 3D-QSAR and molecular docking has been undertaken to explore the structural insights of pyrazolyl urea p38 kinase inhibitors. The 3D-QSAR studies involved comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA). The best CoMFA model was derived from the atom fit alignment with a cross-validated r ² (q ²) value of 0.516 and conventional r ² of 0.950, while the best CoMSIA model yielded a q ² of 0.455 and r ² of 0.979 (39 molecules in training set, 9 molecules in test set). The CoMFA and CoMSIA contour maps generated from these models provided inklings about the influence of interactive molecular fields in the space on the activity. GOLD, Sybyl (FlexX) and AutoDock docking protocols were exercised to explore the protein–inhibitor interactions. The integration of 3D-QSAR and molecular docking has proffered essential structural features of pyrazolyl urea inhibitors and also strategies to design new potent analogues with enhanced activity.