Molecular docking guided 3D-QSAR CoMFA analysis of N-4-Pyrimidinyl-1H-indazol-4-amine inhibitors of leukocyte-specific protein tyrosine kinase

Inhibition of leukocyte-specific protein tyrosine kinase (Lck) activity offers one of the approaches for the treatment of T-cell mediated inflammatory disorders including rheumatoid arthritis, transplant rejection and inflammatory bowel disease. To explore the relationship between the structures of the N-4 Pyrimidinyl-1H-indazol-4-amines and their Lck inhibition, 3D-QSAR study using CoMFA analysis have been performed on a dataset of 42 molecules. The bioactive conformation of the template molecule, selected as the most potent molecule 23 from the series was obtained by performing molecular docking at the ATP binding site of Lck, which is then used to build the rest of the molecules in the series. The constructed CoMFA model is robust with of 0.603 and conventional r² of 0.983. The predictive power of the developed model was obtained using a test set of 10 molecules, giving predictive correlation coefficient of 0.921. CoMFA contour analysis was performed to obtain useful information about the structural requirements for the Lck inhibitors which could be utilized in its future design. Figure CoMFA steric contour map. Sterically favored areas (contribution level 80%) are represented by green polyhedra. Sterically disfavored areas (contribution level 20%) are represented by yellow polyhedra. The active molecule 23 shown in capped sticks. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

3D-QSAR and molecular docking studies of 2-pyrimidinecarbonitrile derivatives as inhibitors against falcipain-3

The three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on a series of falcipain-3 inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. A training set containing 42 molecules served to establish the QSAR models. The optimum CoMFA and CoMSIA models obtained for the training set were statistically significant with cross-validated correlation coefficients r(cv)(2) (q(2)) of 0.549 and 0.608, and conventional correlation coefficients (r(2)) of 0.976 and 0.932, respectively. An independent test set of 12 molecules validated the external predictive power of both models with predicted correlation coefficients (r(pred)(2)) for CoMFA and CoMSIA as 0.697 and 0.509, respectively. The docking of inhibitors into falcipain-3 active site using GOLD software revealed the vital interactions and binding conformation of the inhibitors. The CoMFA and CoMSIA field contour maps agree well with the structural characteristics of the binding pocket of falcipain-3 active site, which suggests that the information rendered by 3D-QSAR models and the docking interactions can provide guidelines for the development of improved falcipain-3 inhibitors.     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

All-Orientation Search and All-Placement Search in Comparative Molecular Field Analysis

Three sets of molecules have been used to study the conventional CoMFA procedure. For all the three test sets, the resulting q² values were observed to vary simply because of the change in the orientation or placement of the aligned molecules. The reason is believed to root in the imperfect sampling of the molecular field. We have introduced two new strategies, all-orientation search (AOS) and all-placement search (APS), to optimize the sampling process. By rotating and translating the molecular aggregate within the grid systematically, all the possible samplings of the molecular field are tested and subsequently the one with the highest q² value can be picked out. We have also demonstrated that the combined application of AOS/APS with GOLPE procedure can yield results better than the ones by using them respectively.Electronic Supplementary Material available.     

Hologram and 3D-quantitative structure toxicity relationship studies of azo dyes

Amino azobenzenes are important dyes in the food and textile industry but their application is limited due to their mutagenicity. Computational modeling techniques were used to help understand the factors responsible for mutagenicity, and several quantitative structure toxicity relationship (QSTR) models have been derived. HQSTR (hologram QSTR) analyses indicated that different substituents at sites on both rings contribute to mutagenicity. Fragment parameters such as bond (B) and connectivity(C), as well as donor-acceptor (DA)-based model provide significant results (q² = 0.59, r² = 0.92, ) explaining these harmful effect. HQSTR results indicated that a bulky group at ring “Y” and small group at ring “X” might help to decrease mutagenicity. 3D-QSTR based on comparative molecular field analyses (CoMFA) and comparative molecular similarity index analyses (CoMSIA) are also in agreement with HQSTR. The 3D QSTR studies reveal that steric and electrostatic field effects have a strong relationship with mutagenicity (for CoMFA: q² = 0.51, r² = 0.95, and for CoMSIA: q² = 0.51, r² = 0.93 and ). In summary, negative groups and steric bulk at ring “Y” and small groups at carbon-3 of ring “X” might be helpful in reducing the mutagenicity of azo dyes.     

CoMFA based de novo design of Pyrrolidine Carboxamides as Inhibitors of Enoyl Acyl Carrier Protein Reductase from Mycobacterium tuberculosis

InhA, the enoyl acyl carrier protein reductase (EACP reductase) from Mycobacterium tuberculosis, is one of the key enzymes involved in the mycobacterial fatty acid elongation cycle and has been validated as an effective target for the development of anti-microbial agents. We report here, comparative molecular field analysis (CoMFA) studies and subsequent de novo ligand design using the LeapFrog program on pyrrolidine carboxamides, which have been reported as selective inhibitors of EACP reductase from Mycobacterium tuberculosis. The CoMFA model, constructed from the inhibitors used in this study has been successfully used to rationalize the structure-activity relationship of pyrrolidine carboxamides. The CoMFA model produced statistically significant results with cross-validated and conventional correlation coefficients of 0.626 and 0.953 respectively. Further, the predictive ability of CoMFA model was determined using a test set which gave predictive correlation coefficient r ² _pred of 0.880, indicating good predictive power. Finally, Leapfrog was used to propose 13 new pyrrolidine carboxamide analogues, based on the information derived from the CoMFA contour maps. The designed molecules showed better predicted activity using the CoMFA model with respect to the already reported systems; hence suggesting that newly proposed molecules in this series of compounds may be more potent and selective toward EACP reductase inhibition.     

Synthesis of novel hetero ring fused pyridine derivatives; Their anticancer activity,CoMFA and CoMSIA studies

A series of novel furo[2,3-b]pyridine-2-carboxamide 4a–h/pyrido[3′,2′:4,5]furo[3,2-d] pyrimidin-4(3H)-one derivatives 5a–p were prepared from pyridin 2(1H) one 1 via selective O-alkylation with α-bromoethylester followed by cyclization, then reaction with different aliphatic primary amines to obtain 4 and further reaction with triethyl orthoacetate/triethyl orthoformate. Also prepared novel furo[2,3-b]pyridine-2-carbohydrazide Schiff’s bases 7a–h and pyrido [3′,2′:4,5]furo[3,2-d]pyrimidin-4(3H)-one derivatives 8a–h starting from furo[2,3-b]pyridine carboxylate derivatives 3 by reaction with hydrazine hydrate to form 6 and reaction with diverse substituted aldehydes and cyclization. Products 4a–h, 5a–p, 7a–h and 8a–h were screened against four human cancer cell lines (HeLa, COLO205, Hep G2 and MCF 7) and one normal cell line (HEK 293). Compounds 4e, 4f, 4g, 5h, 7c, 7d, 7e and 7f showed significant anticancer activity against all the cell lines at micro molar concentration and found to be non-toxic to normal cell line. Studies for HeLa, COLO205 and MCF-7 using CoMFA and CoMSIA. Models from 3D-QSAR provided a strong basis for future rational design of more active and selective HeLa, COLO205 and MCF-7 cell line inhibitors.