Ligand-based and e-pharmacophore modeling, 3D-QSAR and hierarchical virtual screening to identify dual inhibitors of spleen tyrosine kinase (Syk) and janus kinase 3 (JAK3)

The clinical efficacy of multiple kinase inhibitors has caught the interest of Pharmaceutical and Biotech researchers to develop potential drugs with multi-kinase inhibitory activity for complex diseases. In the present work, we attempted to identify dual inhibitors of spleen tyrosine kinase (Syk) and janus kinase 3 (JAK3), keys players in immune signaling, by developing ideal pharmacophores integrating Ligand-based pharmacophore models (LBPMs) and Structure-based pharmacophore models (SBPMs), thereby projecting the optimum pharmacophoric required for inhibition of both the kinases. The four point LBPM; ADPR.14 suggested the presence of one hydrogen bond acceptor, one hydrogen bond donor, one positive ionizable, and one ring aromatic feature for Syk inhibitory activity and AADH.54 proposed the necessity of two hydrogen bond acceptor, one hydrogen bond donor, and one hydrophobic feature for JAK3 inhibitory activity. To our interest, SBPMs identified additional ring aromatic features required for inhibition of both the kinases. For Syk inhibitory activity, the hydrogen bond acceptor feature indicated by LBPM was devoid of forming hydrogen bonding interaction with the hinge region amino acid residue (Ala451). Thus merging the information revealed by both LBPMs and SBPMs, ideal pharmacophore models i.e. ADPRR.14 (Syk) and AADHR.54 (JAK3) were generated. These models after rigorous statistical validation were used for screening of Asinex database. The systematic virtual screening protocol, including pharmacophore and docking-based screening, ADME property, and MM-GBSA energy calculations, retrieved final 10 hits as dual inhibitors of Syk and JAK3. Final 10 hits thus obtained can aid in the development of potential therapeutic agents for autoimmune disorders. Also the top two hits were evaluated against both the enzymes.     

Pharmacophore mapping and in silico screening to identify new potent leads for A(2A) adenosine receptor as antagonists

A(2A) adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson's disease. A 3D-QSAR study of A(2A) AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A(2A) AR are analysed using docking studies and novel potent ligands of A(2A) AR are projected.     

Pharmacophore mapping and in silico screening to identify new potent leads for A2A adenosine receptor as antagonists

A_2A adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson’s disease. A 3D-QSAR study of A_2A AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A_2A AR are analysed using docking studies and novel potent ligands of A_2A AR are projected.     

Identification of Novel HIV 1- Protease Inhibitors: Application of Ligand and Structure Based Pharmacophore Mapping and Virtual Screening

A combined ligand and structure-based drug design approach provides a synergistic advantage over either methods performed individually. Present work bestows a good assembly of ligand and structure-based pharmacophore generation concept. Ligand-oriented study was accomplished by employing the HypoGen module of Catalyst in which we have translated the experimental findings into 3-D pharmacophore models by identifying key features (four point pharmacophore) necessary for interaction of the inhibitors with the active site of HIV-1 protease enzyme using a training set of 33 compounds belonging to the cyclic cyanoguanidines and cyclic urea derivatives. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond acceptors and two hydrophobic, showed a correlation (r) of 0.90 and a root mean square of 0.71 and cost difference of 56.59 bits between null cost and fixed cost. The model was validated using CatScramble technique, internal and external test set prediction. In the second phase of our study, a structure-based five feature pharmacophore hypothesis was generated which signifies the importance of hydrogen bond donor, hydrogen bond acceptors and hydrophobic interaction between the HIV-1 protease enzyme and its inhibitors. This work has taken a significant step towards the full integration of ligand and structure-based drug design methodologies as pharmacophoric features retrieved from structure-based strategy complemented the features from ligand-based study hence proving the accuracy of the developed models. The ligand-based pharmacophore model was used in virtual screening of Maybridge and NCI compound database resulting in the identification of four structurally diverse druggable compounds with nM activities.     

Construction of a three-dimensional pharmacophore for Bcl-2 inhibitors by flexible docking and the multiple copy simultaneous search method

B-Cell lymphoma-2 (Bcl-2) protein is a new promising target for anticancer drugs. A number of anticancer Bcl-2 inhibitors with diverse chemical structures have been discovered in recent years. In this paper, the flexible docking was performed to determine the binding modes of the representative inhibitors from different structural types. Subsequently, the binding modes of inhibitor were used to construct a primary three- dimensional (3D) pharmacophore model. It proved that this model can effectively disrupt the binding of the BH3 domain of proapoptotic Bcl-2 family members to Bcl-2, and match the structural requirement of a new type of Bcl-2 inhibitors. However, these distances between pharmacophoric points are not optimal due to the fact that not all of individual functional groups are located in the ideal position when inhibitors bind to its receptor. In this paper, we proposed a new idea to improve the quality of the pharmacophore model: the multiple copy simultaneous search (MCSS) method was performed to determine the energetically favorable distribution of functional groups with similar features to these pharmacophoric points in the active site of Bcl-2 first. Then their most energetically favorable minima in the positions near the pharmacophoric points were used to optimize the distances between pharmacophoric points. By examining the binding modes of several inhibitors from the same structural type, it was found that the more potent the inhibitor was, the closer it was to the optimized distances between pharmacophoric points. The optimized 3D pharmacophore model obtained in this paper may provide a good starting point for further rational design of Bcl-2 inhibitors.     

Pharmacoinformatics analysis to identify inhibitors of Mtb-ASADH

Aspartate-semialdehyde dehydrogenase (ASADH; EC 1.2.1.11) is a key enzyme in the biosynthesis of essential amino acids in prokaryotes and fungi, inhibition of ASADH leads to the development of novel antitubercular agents. In the present work, a combined structure and ligand-based pharmacophore modeling, molecular docking, and molecular dynamics (MD) approaches were employed to identify potent inhibitors of mycobacterium tuberculosis (Mtb)-ASADH. The structure-based pharmacophore hypothesis consists of three hydrogen bond acceptor (HBA), two negatively ionizable, and one positively ionizable center, while ligand-based pharmacophore consists of additional one HBA and one hydrogen bond donor features. The validated pharmacophore models were used to screen the chemical databases (ZINC and NCI). The screened hits were subjected to ADME and toxicity filters, and subsequently to the molecular docking analysis. Best-docked 25 compounds carry the characteristics of highly electronegative functional groups (–COOH and –NO₂) on both sides and exhibited the H-bonding interactions with highly conserved residues Arg99, Arg249, and His256. For further validation of docking results, MD simulation studies were carried out on two representative compounds NSC51108 and ZINC04203124. Both the compounds remain bound to the key active residues of Mtb-ASADH during the MD simulations. These identified hits can be further used for lead optimization and in the design more potent inhibitors against Mtb-ASADH.     

Pharmacophoric features of Pseudomonas aeruginosa deacetylase LpxC inhibitors: an electronic and structural analysis

Various electronic properties of structurally diverse synthetic LpxC inhibitors containing oxazoline, aroylserine and thiazoline rings were calculated and correlated with biological activity. These electronic features include the magnitude and locations of 3-dimensional molecular electrostatic potentials, hydrogen bond acceptor/donor density, lowest unoccupied molecular orbital, and highest occupied molecular orbital. Strong correlation of these stereo-electronic properties with LpxC inhibitory potency reveals the potential pharmacophoric features of specific LpxC inhibitors. Thus, these pharmacophoric features of LpxC inhibitors based on electronic and surface analysis could be successfully exploited for designing more potent LpxC inhibitors.     

Green synthesis,biological evaluation,molecular docking studies and 3D-QSAR analysis of novel phenylalanine linked quinazoline-4(3H)-one-sulphonamide hybrid entities distorting the malarial reductase activity in folate pathway

A modified Grimmel's method for N-heterocyclization of phenylalanine linked sulphonamide side arm at position-2 was optimized leading to 2,3-disustituted-4-quinazolin-(3H)-ones. Further, [Bmim][BF₄]-H2O (IL) was used as green solvent as well as catalyst for the synthesis of twenty two hybrid quinazolinone motifs (4a-4v) by N-heterocyclization reaction using microwave irradiation technique. The in vitro screening of the hybrid entities against the malarial species Plasmodium falciparum yielded five potent molecules 4l, 4n, 4r, 4t & 4u owing comparable antimalarial activity to the reference drugs. In continuation, an in silico study was carried out to obtain a pharmacophoric model and quantitative structure activity relationship. We also built a 3D-QSAR model to procure more information that could be applied to design new molecules with more potent Pf-DHFR inhibitory activity. The designed pharmacophore was recognized to be more potent for the selected molecules, exhibiting five pharmacophoric features. The active scaffolds were further evaluated for enzyme inhibition efficacy against alleged receptor Pf-DHFR computationally and in vitro, proving their candidature as lead dihydrofolate reductase inhibitors as well as the selectivity of the test candidates was ascertained by toxicity study against vero cells. The perception of good oral bioavailability was also proved by study of pharmacokinetic properties.     

Structure of lactate dehydrogenase from Plasmodium vivax: complexes with NADH and APADH

Malaria caused by Plasmodium vivax is a major cause of global morbidity and, in rare cases, mortality. Lactate dehydrogenase is an essential Plasmodium protein and, therefore, a potential antimalarial drug target. Ideally, drugs directed against this target would be effective against both major species of Plasmodium, P. falciparum and P. vivax. In this study, the crystal structure of the lactate dehydrogenase protein from P. vivax has been solved and is compared to the equivalent structure from the P. falciparum enzyme. The active sites and cofactor binding pockets of both enzymes are found to be highly similar and differentiate these enzymes from their human counterparts. These structures suggest effective inhibition of both enzymes should be readily achievable with a common inhibitor. The crystal structures of both enzymes have also been solved in complex with the synthetic cofactor APADH. The unusual cofactor binding site in these Plasmodium enzymes is found to readily accommodate both NADH and APADH, explaining why the Plasmodium enzymes retain enzymatic activity in the presence of this synthetic cofactor.     

Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method

BackgroundChagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives.Conclusions/SignificanceFMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs.     

Pharmacophoric pattern matching in files of 3D chemical structures: evaluation of search performance

A computerized retrieval system is described which allows the interactive identification of pharmacophoric patterns in 3D chemical structures. A query pharmacophore in the system is matched against a file of structures using a three-stage retrieval mechanism. The initial search involves a rapid bit-string matching operation using screens denoting the presence or absence of interatomic distances that lie within certain predefined ranges. These ranges are selected by an analysis of the frequencies of occurrence of interatomic distances in the data file that is being searched. The second stage involves the generation of all interatomic distances in potentially matching structures in the file to ascertain whether all of the query distances are present. If this is the case, a final search is carried out to determine whether these distances are in the correct geometric orientation to each other; this search involves a novel matching algorithm based upon the set reduction techniques used for substructure searching in files of 2D chemical structures. The evaluation of the system has involved the searching of several pharmacophoric patterns against a 12 728-compound subset of the Cambridge Crystallographic Data Bank. Measurements of the numbers of structures matching at each of the three levels of search show that the system provides a highly efficient means of identifying structures containing a query pharmacophore.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

Substrate specificity of the ileal and the hepatic Na(+)/bile acid cotransporters of the rabbit. II. A reliable 3D QSAR pharmacophore model for the ileal Na(+)/bile acid cotransporter

To design a reliable 3D QSAR model of the intestinal Na(+)/bile acid cotransporter, we have used a training set of 17 inhibitors of the rabbit ileal Na(+)/bile acid cotransporter. The IC(50) values of the training set of compounds covered a range of four orders of magnitude for inhibition of [(3)H]cholyltaurine uptake by CHO cells expressing the rabbit ileal Na(+)/bile acid cotransporter allowing the generation of a pharmacophore using the CATALYST algorithm. After thorough conformational analysis of each molecule, CATALYST generated a pharmacophore model characterized by five chemical features: one hydrogen bond donor, one hydrogen bond acceptor, and three hydrophobic features. The 3D pharmacophore was enantiospecific and correctly estimated the activities of the members of the training set. The predicted interactions of natural bile acids with the pharmacophore model of the ileal Na(+)/bile acid cotransporter explain exactly the experimentally found structure;-activity relationships for the interaction of bile acids with the ileal Na(+)/bile acid cotransporter (Kramer et al. 1999. J. Lipid. Res. 40: 1604;-1617). The natural bile acid analogues cholyltaurine, chenodeoxycholyltaurine, or deoxycholyltaurine were able to map four of the five features of the pharmacophore model: a) the five-membered ring D and the methyl group at position 18 map one hydrophobic site and the 21-methyl group of the side chain maps a second hydrophobic site; b) one of the alpha-oriented hydroxyl groups at position 7 or 12 fits the hydrogen bond donor feature; c) the negatively charged side chain acts as hydrogen bond acceptor; and d) the hydroxy group at position 3 does not specifically map any of the five binding features of the pharmacophore model. The 3-hydroxy group of natural bile acids is not essential for interactions with ileal or hepatic Na(+)/bile acid cotransporters. A modification of the 3-position of a natural bile acid molecule is therefore the preferred position for drug targeting strategies using bile acid transport pathways.     

Toll-Like Receptor 7 Agonists: Chemical Feature Based Pharmacophore Identification and Molecular Docking Studies

Chemical feature based pharmacophore models were generated for Toll-like receptors 7 (TLR7) agonists using HypoGen algorithm, which is implemented in the Discovery Studio software. Several methods tools used in validation of pharmacophore model were presented. The first hypothesis Hypo1 was considered to be the best pharmacophore model, which consists of four features: one hydrogen bond acceptor, one hydrogen bond donor, and two hydrophobic features. In addition, homology modeling and molecular docking studies were employed to probe the intermolecular interactions between TLR7 and its agonists. The results further confirmed the reliability of the pharmacophore model. The obtained pharmacophore model (Hypo1) was then employed as a query to screen the Traditional Chinese Medicine Database (TCMD) for other potential lead compounds. One hit was identified as a potent TLR7 agonist, which has antiviral activity against hepatitis virus in vitro. Therefore, our current work provides confidence for the utility of the selected chemical feature based pharmacophore model to design novel TLR7 agonists with desired biological activity.     

Conformational and electrostatic analysis of SN1 donor analogue glycomimetic inhibitors of ST3Gal-I mammalian sialyltransferase

Mammalian sialyltransferases play a role in the metastasis of various cancers in humans. Inhibitors of these enzymes will in principle be able to directly inhibit aberrant sialylation in cancer. Inhibitors of ST3Gal-I resembling the donor component of S_N1 Transition State structures were previously evaluated as part of a kinetics study. Here, using classical dynamics simulations and free energy perturbation calculations, we rationalize the performance of three of these donor analogue ST3Gal-I enzyme inhibitors. We find to inhibit the mammalian ST3Gal-I enzyme a donor analogue requires configurationally limited functionality. This is mediated by the binding of the inhibitor to the enzyme. The inhibitor’s ability to interact with Y194 and T272 through a charged group such as a carboxylate is especially important. Furthermore, a conformational rigid form approximating the donor substrate is central. Here this is achieved by an intramolecular hydrogen bond formed between the carboxylate group and one of the ribose hydroxyl groups of the cytidine monophosphate (CMP) leaving group. This intramolecular interaction results in the donor substrate conformer complimenting the form of the catalytic binding site. Finally the carboxylate charge is essential for electrostatic pairing with the binding site. Substituting this group for an alcohol or amide results in severe weakening of the ligand binding. The carboxylate thus proves an to be an irreplaceable functional group and an essential pharmacophore.     

In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a negative regulator of insulin and leptin signalling pathway; hence, it can be considered as a new therapeutic target of intervention for the treatment of type 2 diabetes. Inhibition of this molecular target takes care of both diabetes and obesity, i.e. diabestiy. In order to get more information on identification and optimization of lead, pharmacophore modelling, atom-based 3D QSAR, docking and molecular dynamics studies were carried out on a set of ligands containing thiazolidine scaffold. A six-point pharmacophore model consisting of three hydrogen bond acceptor (A), one negative ionic (N) and two aromatic rings (R) with discrete geometries as pharmacophoric features were developed for a predictive 3D QSAR model. The probable binding conformation of the ligands within the active site was studied through molecular docking. The molecular interactions and the structural features responsible for PTP1B inhibition and selectivity were further supplemented by molecular dynamics simulation study for a time scale of 30 ns. The present investigation has identified some of the indispensible structural features of thiazolidine analogues which can further be explored to optimize PTP1B inhibitors.     

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.     

Crystal structure of gamma-chymotrypsin in complex with 7-hydroxycoumarin

The 1.8 A crystal structure of 7-hydroxycoumarin (7-HC) bound to chymotrypsin reveals that this inhibitor forms a planar cinnamate acyl-enzyme complex. The phenyl ring of the bound inhibitor forms numerous van der Waals contacts in the S1 pocket of the enzyme, with the p-hydroxyl group donating a hydrogen bond to the main-chain oxygen atom of Ser217, and the o-hydroxyl group forming a water-mediated hydrogen bond with the carbonyl oxygen of Val227. The structure of the acyl-enzyme complex suggests that the mechanism of inhibition of 7-HC involves nucleophilic attack by the Ser195 O(gamma) atom on the carbonyl carbon atom of the inhibitor, accompanied by the breaking of the 2-pyrone ring of the inhibitor, and leading to the formation of a cinnamate acyl-enzyme derivative via a tetrahedral transition state. Comparisons with structures of photoreversible cinnamates bound to chymotrypsin reveal that although 7-HC interacts with the enzyme in a similar fashion, the binding of 7-HC to chymotrypsin takes place in a productive conformation in contrast to the photoreversible cinnamates. In summary, the 7-HC-chymotrypsin complex provides basic insight into the inhibition of chymotrypsin by natural coumarins and provides a structural basis for the design of more potent mechanism-based inhibitors against a wide range of biologically important chymotrypsin-like enzymes.     

Design of antimalarial transmission blocking agents: Pharmacophore mapping of ligands active against stage-V mature gametocytes of Plasmodium falciparum

Abstract

The discovery of transmission-blocking (T-B) agents is crucial for preventing and complete removal of malaria infection. However, most of the existing antimalarials are only active against the asexual stages of Plasmodium parasite, but ineffective against the sexual stage (gametocytes). In this background, we have developed pharmacophore models against the stage-V mature gametocytes of P. falciparum parasites. The pharmacophore model (Hypo-1) showed five pharmacophoric features namely, one hydrogen bond donor (HBD), one hydrophobic aliphatic (HYAl), one ring aromatic (RA), and two hydrophobic aromatic (HYAr) essential for the anti-gametocytic activity. The amino, methyl, fused phenyl ring of the quinazoline heterocycle, two phenyl rings of biphenyl moiety (HBD, HYAl, HYAr1, HYAr2 and RA) are the crucial features responsible for the non-specific anti-gametocytic activity (PfG). Subsequently, the model (Hypo-2) developed against the stage-V female gametocytes (PffG) showed the contribution of three pharmacophoric features namely, two hydrogen bond acceptor (HYA) and one RA required for the anti-gametocytic activity. The sulfhydryl, imine and pyridyl groups are observed to be essential for anti-gametocytic activity against female gametocytes. Both the models (PfG and PfGG) showed the classification accuracies of 78.26 and 71.64% for training set compounds and 60.80 and 60.18% for the test set compounds, respectively, for classification of compounds into higher and lower active classes. Also, both the models were found to retain the higher active compounds (IC₅₀ <100?nM) in top 1% of total compounds (actives and decoys) as observed after screening the decoy set compounds.

Communicated by Ramaswamy H Sarma