Neuraminidase pharmacophore model derived from diverse classes of inhibitors

A three-dimensional pharmacophore model was developed based on 22 currently available inhibitors, which were carefully selected with great diversity in both molecular structure and bioactivity, for discovering new potent neuraminidase (NA) inhibitors to fight against avian influenza virus. The best hypothesis (Hypo1), consisting of five features, namely, one positive ionizable group, one negative ionizable group, one hydrophobic point, and two hydrogen-bond donors, has a correlation coefficient of 0.902, a root mean square deviation of 1.392, and a cost difference of 72.88, suggesting that a highly predictive pharmacophore model was successfully obtained. The application of the model shows great success in predicting the activities of 88 known NA inhibitors in our test set with a correlation coefficient of 0.818 with a cross-validation of 98% confidence level. Accordingly, our model should be reliable in identifying structurally diverse compounds with desired biological activity.     

Pharmacophore mapping of diverse classes of farnesyltransferase inhibitors

Protein farnesyltransferase (FTase) is a zinc-dependent enzyme that catalyzes the attachment of a farnesyl lipid group to the sulfur atom of a cysteine residue of numerous proteins involved in cell signaling including the oncogenic H-Ras protein. Pharmacophore models were developed by using Catalyst HypoGen program with a training set of 22 farnesyltransferase inhibitors (FTIs), which were carefully selected with great diversity in both molecular structure and bioactivity for discovering new potent FTIs. The best pharmacophore hypothesis (Hypo 1), consisting of four features, namely, one hydrogen-bond acceptor (HBA), one hydrophobic point (HY), and two ring aromatics (RA), has a correlation coefficient of 0.961, a root mean square deviation (RMSD) of 0.885, and a cost difference of 62.436, suggesting that a highly predictive pharmacophore model was successfully obtained. For the test series, a classification scheme was used to distinguish highly active from moderately active and inactive compounds on the basis of activity ranges. Hypo 1 was validated with 181 test set compounds, which has a correlation coefficient of 0.713 between estimated activity and experimentally measured activity. The model was further validated by screening a database spiked with 25 known inhibitors. The model picked up all 25 known inhibitors giving an enrichment factor of 10.892. The results demonstrate that the hypothesis derived in this study can be considered to be a useful and reliable tool in identifying structurally diverse compounds with desired biological activity.     

Pharmacophore identification of Raf-1 kinase inhibitors

A three-dimensional pharmacophore model was developed based on 25 currently available Raf-1 kinase inhibitors. The best pharmacophore hypothesis (Hypo1), consisting of four chemical features (one hydrogen-bond acceptor, one hydrogen-bond donor, and two hydrophobic groups), has a correlation coefficient of 0.972. The results of our study provide a valuable tool in designing new leads with desired biological activity by virtual screening.     

Pharmacophore identification of KSP inhibitors

A three-dimensional pharmacophore model was developed based on 25 currently available KSP (kinesin spindle protein) inhibitors in Catalyst software package. The best pharmacophore hypothesis (Hypo1), consisting of four chemical features (one hydrogen-bond acceptor, one hydrogen-bond donor, one aromatic ring, and one hydrophobic group), has a correlation coefficient of 0.965. The results of our study provide a valuable tool in designing new leads with desired biological activity by virtual screening.     

Pharmacophore modeling and virtual screening to identify potential RET kinase inhibitors

Chemical features based 3D pharmacophore model for REarranged during Transfection (RET) tyrosine kinase were developed by using a training set of 26 structurally diverse known RET inhibitors. The best pharmacophore hypothesis, which identified inhibitors with an associated correlation coefficient of 0.90 between their experimental and estimated anti-RET values, contained one hydrogen-bond acceptor, one hydrogen-bond donor, one hydrophobic, and one ring aromatic features. The model was further validated by a testing set, Fischer’s randomization test, and goodness of hit (GH) test. We applied this pharmacophore model to screen NCI database for potential RET inhibitors. The hits were docked to RET with GOLD and CDOCKER after filtering by Lipinski’s rules. Ultimately, 24 molecules were selected as potential RET inhibitors for further investigation.     

Pharmacophore modeling and virtual screening for designing potential 5-Lipoxygenase inhibitors

Inhibitors of the 5-Lipoxygenase (5-LOX) pathway have a therapeutic potential in a variety of inflammatory disorders such as asthma. In this study, chemical feature based pharmacophore models of inhibitors of 5-LOX have been developed with the aid of HipHop and HypoGen modules within Catalyst program package. The best quantitative pharmacophore model, Hypo1, which has the highest correlation coefficient (0.97), consists of two hydrogen-bond acceptors, one hydrophobic feature and one ring aromatic feature. Hypo1 was further validated by test set and cross validation method. The application of the model shows great success in predicting the activities of 65 known 5-LOX inhibitors in our test set with a correlation coefficient of 0.85 with a cross validation of 95% confidence level, proving that the model is reliable in identifying structurally diverse compounds for inhibitory activity against 5-LOX. Furthermore, Hypo1 was used as a 3D query for screening Maybridge and NCI databases within catalyst and also drug like compounds obtained from Enamine Ltd, which follow Lipinski’s rule of five. The hit compounds were subsequently subjected to filtering by docking and visualization, to identify the potential lead molecules. Finally 5 potential lead compounds, identified in the above process, were evaluated for their inhibitory activities. These studies resulted in the identification of two compounds with potent inhibition of 5-LOX activity with IC₅₀ of 14 μM and 35 μM, respectively. These studies thus validate the pharmacophore model generated and suggest the usefulness of the model in screening of various small molecule libraries and identification of potential lead compounds for 5-LOX inhibition.     

Virtual screening and pharmacophore studies for ftase inhibitors using Indian plant anticancer compounds database

Farnesyl transferase (FTase) is an enzyme responsible for post-translational modification in proteins having a carboxy-terminal CaaX motif in human. It catalyzes the attachment of a lipid group in proteins of RAS superfamily, which is essential in signal transduction. FTase has been recognized as an important target for anti cancer therapeutics. In this work, we performed virtual screening against FTase with entire 125 compounds from Indian Plant Anticancer Database using AutoDock 3.0.5 software. All compounds were docked within binding pocket containing Lys164, Tyr300, His248 and Tyr361 residues in crystal structure of FTase. These complexes were ranked according to their docking score, using methodology that was shown to achieve maximum accuracy. Finally we got three potent compounds with the best Autodock docking Score (Vinorelbine: -21.28 Kcal/mol, Vincristine: -21.74 Kcal/mol and Vinblastine: -22.14 Kcal/mol) and their energy scores were better than the FTase bound co-crystallized ligand (L- 739: -7.9 kcal/mol). These three compounds belong to Vinca alkaloids were analyzed through Python Molecular Viewer for their interaction studies. It predicted similar orientation and binding modes for these compounds with L-739 in FTase.Thus from the complex scoring and binding ability it is concluded that these Vinca alkaloids could be promising inhibitors for FTase. A 2-D pharmacophore was generated for these alkaloids using LigandScout to confirm it. A shared feature pharmacophore was also constructed that shows four common features (one hydogen bond Donar, Two hydrogen bond Acceptor and one ionizable area) help compounds to interact with this enzyme.     

Ligand-based and structure-based approaches in identifying ideal pharmacophore against c-Jun N-terminal kinase-3

Structure and ligand based pharmacophore modeling and docking studies carried out using diversified set of c-Jun N-terminal kinase-3 (JNK3) inhibitors are presented in this paper. Ligand based pharmacophore model (LBPM) was developed for 106 inhibitors of JNK3 using a training set of 21 compounds to reveal structural and chemical features necessary for these molecules to inhibit JNK3. Hypo1 consisted of two hydrogen bond acceptors (HBA), one hydrogen bond donor (HBD), and a hydrophobic (HY) feature with a correlation coefficient (r²) of 0.950. This pharmacophore model was validated using test set containing 85 inhibitors and had a good r² of 0.846. All the molecules were docked using Glide software and interestingly, all the docked conformations showed hydrogen bond interactions with important hinge region amino acids (Gln155 and Met149) and these interactions were compared with Hypo1 features. The results of ligand based pharmacophore model (LBPM) and docking studies are validated each other. The structure based pharmacophore model (SBPM) studies have identified additional features, two hydrogen bond donors and one hydrogen bond acceptor. The combination of these methodologies is useful in designing ideal pharmacophore which provides a powerful tool for the discovery of novel and selective JNK3 inhibitors.     

BCR-ABL tyrosine kinase inhibitor pharmacophore model derived from a series of phenylaminopyrimidine-based (PAP) derivatives

To reveal novel insights into the inhibition of BCR-ABL tyrosine kinase, pharmacophore mapping studies were performed for a series of phenylaminopyrimidine-based (PAP) derivatives, including imatinib (Gleevec). A seven-point pharmacophore model with one hydrophobic group (H), two hydrogen bond donors (D) and four aromatic rings (R) was developed using phase (pharmacophore alignment & scoring engine). The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of 0.886 and a survival score of 4.97 for training set molecules. The model showed excellent predictive power, with a correlation coefficient of Q² = 0.768 for an external test set of ten molecules. The results obtained from our studies provide a valuable tool for designing new lead molecules with potent activity.     

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Pharmacophore modeling and virtual screening studies to identify new c-Met inhibitors

Mesenchymal epithelial transition factor (c-Met) is an attractive target for cancer therapy. Three-dimensional pharmacophore hypotheses were built based on a set of known structurally diverse c-Met inhibitors. The best pharmacophore model, which identified inhibitors with an associated correlation coefficient of 0.983 between their experimental and estimated IC(50) values, consisted of two hydrogen-bond acceptors, one hydrophobic, and one ring aromatic feature. The highly predictive power of the model was rigorously validated by test set prediction and Fischer's randomization method. The high values of enrichment factor and receiver operating characteristic (ROC) score indicated the model performed fairly well at distinguishing active from inactive compounds. The model was then applied to screen compound database for potential c-Met inhibitors. A filtering protocol, including druggability and molecular docking, were also applied in hits selection. The final 38 molecules, which exhibited good estimated activities, desired binding mode and favorable drug likeness were identified as potential c-Met inhibitors. Their novel backbone structures could be served as scaffolds for further study, which may facilitate the discovery and rational design of potent c-Met kinase inhibitors.     

Pharmacophore modeling, virtual screening and 3D-QSAR studies of 5-tetrahydroquinolinylidine aminoguanidine derivatives as sodium hydrogen exchanger inhibitors

Sodium hydrogen exchanger (SHE) inhibitor is one of the most important targets in treatment of myocardial ischemia. In the course of our research into new types of non-acylguanidine, SHE inhibitory activities of 5-tetrahydroquinolinylidine aminoguanidine derivatives were used to build pharmacophore and 3D-QSAR models. Genetic Algorithm Similarity Program (GASP) was used to derive a 3D pharmacophore model which was used in effective alignment of data set. Eight molecules were selected on the basis of structure diversity to build 10 different pharmacophore models. Model 1 was considered as the best model as it has highest fitness score compared to other nine models. The obtained model contained two acceptor sites, two donor atoms and one hydrophobic region. Pharmacophore modeling was followed by substructure searching and virtual screening. The best CoMFA model, representing steric and electrostatic fields, obtained for 30 training set molecules was statistically significant with cross-validated coefficient (q(2)) of 0.673 and conventional coefficient (r(2)) of 0.988. In addition to steric and electrostatic fields observed in CoMFA, CoMSIA also represents hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. CoMSIA model was also significant with cross-validated coefficient (q(2)) and conventional coefficient (r(2)) of 0.636 and 0.986, respectively. Both models were validated by an external test set of eight compounds and gave satisfactory prediction (r(pred)(2)) of 0.772 and 0.701 for CoMFA and CoMSIA models, respectively. This pharmacophore based 3D-QSAR approach provides significant insights that can be used to design novel, potent and selective SHE inhibitors.     

Towards the synthesis of bisubstrate inhibitors of protein farnesyltransferase: Synthesis and biological evaluation of new farnesylpyrophosphate analogues

Protein farnesyltransferase (FTase) has recently appeared as a new target of parasitic diseases, a field poor in drugs in development. With the aim of creating new bisubstrate inhibitors of FTase, new farnesyl pyrophosphate analogues have been studied. Farnesyl analogues with a malonic acid function exhibited the best inhibitory activity on FTase. This group was introduced into our imidazole-containing model leading to new compounds with submicromolar activities. Kinetic experiments have been realized to determine their binding mode to the enzyme.     

Pharmacophore models generation by catalyst and phase consensus-based virtual screening protocol against PI3Kα inhibitors

Phosphoinositide 3-kinases (PI3Ks) family has emerged as promising targets for novel therapeutic agents against neoplastic diseases. Pharmacophore and 3D-quantitative structure–activity relationship modelling were applied to study the structure–activity relationship of PI3K inhibitors. The best HypoGen pharmacophore hypothesis Hypo1 with a correlation coefficient of 0.961 consists of one hydrogen-bond acceptor, one hydrogen-bond donor and two hydrophobic features, whereas the best phase hypothesis AADRRR.378 with favourable statistics (q² = 0.7368, r² = 0.9863) has two hydrogen-bond acceptors, one hydrogen-bond donor and three ring aromatic features. Multiple methods, such as Fischer validation, molecular docking and mapping of test set molecules, were carried out to validate these pharmacophore models. Furthermore, a comparative molecular similarity indices analysis candidate hypothesis model was generated as a supplement of pharmacophore hypothesis. Detailed protein–ligand binding information obtained by Glide was utilised in compound optimisation and virtual screening. A molecular database of 133 known inhibitors and 6179 decoys was built for a screening test to quantitatively analyse various hypotheses and scoring parameters. Finally, we designed a workflow integrating HypoGen pharmacophore searching, phase pharmacophore searching and molecular docking for screening the database. With an improved criterion of enrichment factor (EF = 17.43) and ROC curve (AUC = 0.946), this workflow would provide us an original method for novel PI3K inhibitors.     

Design, synthesis, and activity of 4-quinolone and pyridone compounds as nonthiol-containing farnesyltransferase inhibitors

As a part of our efforts to identify potent inhibitors of farnesyltransferase (FTase), modification of the structure of tipifarnib through structure-based design was undertaken by replacing the 2-quinolones with 4-quinolones and pyridones, and subsequent relocation of the D-ring to the N-methyl group on the imidazole ring. This study has yielded a novel series of potent and selective FTase inhibitors. The X-ray structure of tipifarnib (1) in complex with FTase was described.     

In silico approaches to identify novel myeloid cell leukemia-1 (Mcl-1) inhibitors for treatment of cancer

Myeloid cell leukemia-1 (Mcl-1) has been a validated and attractive target for cancer therapy. Over-expression of Mcl-1 in many cancers allows cancer cells to evade apoptosis and contributes to the resistance to current chemotherapeutics. Here, we identified new Mcl-1 inhibitors using a multi-step virtual screening approach. First, based on two different ligand-receptor complexes, 20 pharmacophore models were established by simultaneously using ‘Receptor-Ligand Pharmacophore Generation’ method and manual build feature method, and then carefully validated by a test database. Then, pharmacophore-based virtual screening (PB-VS) could be performed by using the 20 pharmacophore models. In addition, docking study was used to predict the possible binding poses of compounds, and the docking parameters were optimized before performing docking-based virtual screening (DB-VS). Moreover, a 3D QSAR model was established by applying the 55 aligned Mcl-1 inhibitors. The 55 inhibitors sharing the same scaffold were docked into the Mcl-1 active site before alignment, then the inhibitors with possible binding conformations were aligned. For the training set, the 3D QSAR model gave a correlation coefficient r² of 0.996; for the test set, the correlation coefficient r² was 0.812. Therefore, the developed 3D QSAR model was a good model, which could be applied for carrying out 3D QSAR-based virtual screening (QSARD-VS). After the above three virtual screening methods orderly filtering, 23 potential inhibitors with novel scaffolds were identified. Furthermore, we have discussed in detail the mapping results of two potent compounds onto pharmacophore models, 3D QSAR model, and the interactions between the compounds and active site residues.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Pharmacophore modeling and virtual screening for designing potential PLK1 inhibitors

Pharmacophore models of Polo-like kinase-1 (PLK1) inhibitors have been established by using the HipHop and HypoGen algorithms implemented in the Catalyst software package. The best quantitative pharmacophore model, Hypo1, which has the highest correlation coefficient (0.9895), consists of one hydrogen bond acceptor, one hydrogen bond donor, one hydrophobic feature, and one hydrophobic aliphatic feature. Hypo1 was further validated by test set and cross validation method. Then Hypo1 was used as a 3D query to screen several databases including Specs, NCI, Maybridge, and Chinese Nature Product Database (CNPD). The hit compounds were subsequently subjected to filtering by Lipinski's rule of five and docking study to refine the retrieved hits and as a result to reduce the rate of false positive. Finally, a total of 20 compounds were selected and have been shifted to in vitro and in vivo studies. As far as we know, this is the first report on the pharmacophore modeling even the first publicly reported virtual screening study of PLK1 inhibitors.     

ASK1 pharmacophore model derived from diverse classes of inhibitors

The three-dimensional pharmacophore model of apoptosis signal-regulating kinase 1 (ASK1) inhibitors has been developed with PharmaGist program. The positions of pharmacophore features in the model correspond to conformations of ASK1 highly active inhibitors in which they interact with ATP-binding site of ASK1. The generated pharmacophore model allows accurately predict active and inactive compounds and can be of great use for virtual screening aimed at discovering novel ASK1 inhibitors.     

Receptor-guided 3D-QSAR approach for the discovery of c-kit tyrosine kinase inhibitors

Studies of the the three-dimensional quantitative structure-activity relationships for ninety-five c-kit tyrosine kinase inhibitors were performed. Based on a co-crystallized compound (1 T46), known inhibitors were aligned with c-kit by induced-fit docking, and multiple training/test set splitting was performed to validate the selected pharmacophore model. The best pharmacophore model consisted of five features: one hydrogen-bond donor and four aromatic rings. Reliable statistics were obtained (R(2) = 0.95, R(pred) (2) = 0.75), and the model was validated by using it to select c-kit inhibitors from a database; 82.1% of the hits it retrieved were active. Accordingly, our model can be reliably used to identify new c-kit inhibitors, and can provide useful information when designing new inhibitors.