Discovery of dual binding site acetylcholinesterase inhibitors identified by pharmacophore modeling and sequential virtual screening techniques

Dual binding site acetylcholinesterase (AChE) inhibitors are promising for the treatment of Alzheimer’s disease (AD). They alleviate the cognitive deficits and AD-modifying agents, by inhibiting the β-amyloid (Aβ) peptide aggregation, through binding to both the catalytic and peripheral anionic sites, the so called dual binding site of the AChE enzyme. In this Letter, chemical features based 3D-pharmacophore models were developed based on the eight potent and structurally diverse AChE inhibitors (I-VIII) obtained from high-throughput in vitro screening technique. The best 3D-pharmacophore model, Hypo1, consists of two hydrogen-bond acceptor lipid, one hydrophobe, and two hydrophobic aliphatic features obtained by Catalyst/HIPHOP algorithm adopted in Discovery studio program. Hypo1 was used as a 3D query in sequential virtual screening study to filter three small compound databases. Further, a total of nine compounds were selected and followed on in vitro analysis. Finally, we identified two leads—Specs1 (IC₅₀ = 3.279 μM) and Spec2 (IC₅₀ = 5.986 μM) dual binding site compounds from Specs database, having good AChE enzyme inhibitory activity.     

Ligand based virtual screening and biological evaluation of inhibitors of chorismate mutase (Rv1885c) from Mycobacterium tuberculosis H37Rv

We have identified new lead candidates that possess inhibitory activity against Mycobacterium tuberculosis H37Rv chorismate mutase by a ligand-based virtual screening optimized for lead evaluation in combination with in vitro enzymatic assay. The initial virtual screening using a ligand-based pharmacophore model identified 95 compounds from an in-house small molecule database of 15,452 compounds. The obtained hits were further evaluated by molecular docking and 15 compounds were short listed based on docking scores and the other scoring functions and subjected to biological assay. Chorismate mutase activity assays identified four compounds as inhibitors of M. tuberculosis chorismate mutase (MtCM) with low K(i) values. The structural models for these ligands in the chorismate mutase binding site will facilitate medicinal chemistry efforts for lead optimization against this protein.     

Discovery of novel SecA inhibitors of Candidatus Liberibacter asiaticus by structure based design

Candidatus Liberibacter asiaticus is the causal agent of Huanglongbing (HLB) disease of citrus. Current management practices have not been able to control HLB and stop the spread of HLB. The current study is focused on screening small molecule inhibitors against SecA protein of Ca. L. asiaticus. Homology modeling, structure based virtual screening and molecular docking methods have been used to find the novel inhibitory compounds against SecA activity at ATP binding region. At 20 μm 17 compounds showed >50% inhibition and four compounds had more than 65% inhibition. The most active compound has IC₅₀ value of 2.5 μM. The differences between the activities of the compounds are explained by their inter-molecular interactions at ATP binding site.     

Identification of novel inhibitors of HCV RNA-dependent RNA polymerase by pharmacophore-based virtual screening and in vitro evaluation

Hepatitis C virus (HCV) is the major etiological agent of non-A, non-B hepatitis where no effective treatment is available. The HCV NS5B with RNA-dependent RNA polymerase (RdRp) activity is a key target for the treatment of HCV infection. Here we report novel NS5B polymerase inhibitors identified by virtual screening and in vitro evaluation of their inhibitory activities. On the basis of a newly identified binding pocket of NS5B, distinct from the nucleotide binding site but highly conserved among various HCV isolates, we performed virtual screening of compounds that fit this binding pocket from the available chemical database of 3.5 million compounds. The inhibitory activities of the in silico selected 119 compounds were estimated with in vitro RdRp assay. Three compounds with IC50 values of about 20 μM were identified, and their kinetic analyses suggest that these compounds are noncompetitive inhibitors with respect to the ribonucleotide substrate. Furthermore, the single-point mutations of the conserved residues in the binding pocket of NS5B resulted in the significant decrease of the RdRp activity, indicating that the binding pocket presented here might be important for the therapeutic intervention of HCV. These novel inhibitors would be useful for the development of effective anti-HCV agents.     

A novel class of Hsp90 inhibitors isolated by structure-based virtual screening

A novel class of 3-phenyl-2-styryl-3H-quinazolin-4-one Hsp90 inhibitors with in vitro anti-tumor activity are identified by structure-based virtual screening of a chemical database with docking simulations in the N-terminal ATP-binding site, in vitro ATPase assay using yeast Hsp90, and cell-based Her2 degradation assay in a consecutive fashion. These results exemplify the usefulness of the structure-based virtual screening with molecular docking in drug discovery. The structural features responsible for a tight binding of the inhibitors in the active site of Hsp90 are discussed in detail.     

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.     

Receptor based 3D-QSAR to identify putative binders of Mycobacterium tuberculosis Enoyl acyl carrier protein reductase

In the current study, the applicability and scope of 3D-QSAR models (CoMFA and CoMSIA) to complement virtual screening using 3D pharmacophore and molecular docking is examined and applied to identify potential hits against Mycobacterium tuberculosis Enoyl acyl carrier protein reductase (MtENR). Initially CoMFA and CoMSIA models were developed using series of structurally related arylamides as MtENR inhibitors. Docking studies were employed to position the inhibitors into MtENR active site to derive receptor based 3D-QSAR models. Both CoMFA and CoMSIA yielded significant cross validated q² values of 0.663 and 0.639 and r² values of 0.989 and 0.963, respectively. The statistically significant models were validated by a test set of eight compounds with predictive r² value of 0.882 and 0.875 for CoMFA and CoMSIA. The contour maps from 3D-QSAR models in combination with docked binding structures help to better interpret the structure activity relationship. Integrated with CoMFA and CoMSIA predictive models structure based (3D-pharmacophore and molecular docking) virtual screening have been employed to explore potential hits against MtENR. A representative set of 20 compounds with high predicted IC₅₀ values were sorted out in the present study.     

Design, synthesis, and biological evaluation of substituted hydrazone and pyrazole derivatives as selective COX-2 inhibitors: Molecular docking study

New arylhydrazone derivatives and a series of 1,5-diphenyl pyrazoles were designed and synthesized from 1-(4-chlorophenyl)-4,4,4-trifuorobutane-1,3-dione 1. The newly synthesized compounds were investigated in vivo for their anti-inflammatory activities using carrageenan-induced rat paw oedema model. Moreover, they were tested for their inhibitory activity against ovine COX-1 and COX-2 using an in vitro cyclooxygenase (COX) inhibition assay. Some of the new compounds (2f, 6a and 6d) showed a reasonable in vitro COX-2 inhibitory activity, with IC?? value of 0.45 μM and selectivity index of 111.1. A virtual screening was carried out through docking the designed compounds into the COX-2 binding site to predict if these compounds have analogous binding mode to the COX-2 inhibitors. Docking study of the synthesized compounds 2f, 6a and 6d into the active site of COX-2 revealed a similar binding mode to SC-558, a selective COX-2 inhibitor.     

Identification of low-molecular inhibitors of proteinase ASK1

Protein kinase ASK1 (Apoptosis signal-regulating kinase 1) plays a key role in cell differentiation, aging and apoptosis. High activity of the kinase is associated with several pathologies. The ASK1 inhibitors might be therapeutic for patients with neurodegenerative, cardiovascular diseases and fibrous histiocytoma. In this work the identification of ASK1 inhibitors was performed by the methods of computer modeling and biochemical testing in vitro. The virtual screening experiments were carried out targeting the ATP binding site of ASK1 by browsing the database which contained 164 840 compounds of diverse chemical classes. The best-scored 300 ligands have been taken for the kinase assay analysis. In vitro tests revealed that derivatives of 2-thioxo-thiazolidin-4-one exhibited inhibitory activity against ASK1. The most active compound was 5-bromo-3-(4-oxo-2-thioxo-thiazolidin-5-ylidene)-1,3-dihydro-indol-2-one (IC50 = 2 microM). Binding mode for inhibitors of this class with ASK1 ATP-binding site was proposed. Our results can be used for further optimization and developing more potent and selective inhibitors of ASK1.     

Discovery of novel HCV polymerase inhibitors using pharmacophore-based virtual screening

We report the use of pharmacophore-based virtual screening as an efficient tool for the discovery of novel HCV polymerase inhibitors. A three-dimensional pharmacophore model for the HCV-796 binding site, NNI site IV inhibitor, to the enzyme was built by means of the structure-based focusing module in Cerius2 program. Using these models as a query for virtual screening, we produced a successful example of using pharmacophore-based virtual screening to identify novel compounds with HCV replicon assay through inhibition of HCV polymerization. Among the hit compounds, compounds 1 and 2 showed 56% and 48% inhibition of NS5B polymerization activity at 20 μM, respectively. In addition, compound 1 also exhibited replicon activity with EC₅₀ value of 2.16 μM. Following up the initial hit, we obtained derivatives of compound 1 and evaluated polymerization inhibition activity and HCV replicon assay. These results provide information necessary for the development of more potent NS5B inhibitors.     

Ca-asp bound X-ray structure and inhibition of Bacillus anthracis dihydroorotase (DHOase)

Dihydroorotase (DHOase) is the third enzyme in the de novo pyrimidine synthesis pathway and is responsible for the reversible cyclization of carbamyl-aspartate (Ca-asp) to dihydroorotate (DHO). DHOase is further divided into two classes based on several structural characteristics, one of which is the length of the flexible catalytic loop that interacts with the substrate, Ca-asp, regulating the enzyme activity. Here, we present the crystal structure of Class I Bacillus anthracis DHOase with Ca-asp in the active site, which shows the peptide backbone of glycine in the shorter loop forming the necessary hydrogen bonds with the substrate, in place of the two threonines found in Class II DHOases. Despite the differences in the catalytic loop, the structure confirms that the key interactions between the substrate and active site residues are similar between Class I and Class II DHOase enzymes, which we further validated by mutagenesis studies. B. anthracis DHOase is also a potential antibacterial drug target. In order to identify prospective inhibitors, we performed high-throughput screening against several libraries using a colorimetric enzymatic assay and an orthogonal fluorescence thermal binding assay. Surface plasmon resonance was used for determining binding affinity (K_D) and competition analysis with Ca-asp. Our results highlight that the primary difference between Class I and Class II DHOase is the catalytic loop. We also identify several compounds that can potentially be further optimized as potential B. anthracis inhibitors.     

Novel GSK-3beta inhibitors from sequential virtual screening

Glycogen synthase kinase-3 (GSK-3beta) has been emerging as a key therapeutic target for type-2 diabetics, Alzheimer's disease, cancer, and chronic inflammation. For the purpose of finding biologically active and novel compounds and providing new idea for drug-design, we performed virtual screening using commercially available database. Three-dimensional common feature pharmacophore model was developed by using HipHop program provided in Catalyst software and it was used as a query for screening database. Recursive partitioning (RP) model was developed as a filtering system, which was able to classify active and inactive compounds. Eventually, a sequential virtual screening procedure (SQSP) was conducted by applying the common feature pharmacophore and RP model in succession to discover novel potent GSK-3beta inhibitors. The final 56 hit compounds were carefully selected considering predicted docking mode in crystal structures. Subsequent enzyme assay for human GSK-3beta protein confirmed that three compounds of these hit compounds exhibit micromolar inhibitory activity. Here, we report novel hit compounds and their binding mode in the active site of GSK-3beta crystal structure.     

Discovery of isatin and 1H-indazol-3-ol derivatives as d-amino acid oxidase (DAAO) inhibitors

d-Amino acid oxidase (DAAO) is a potential target in the treatment of schizophrenia as its inhibition increases brain d-serine level and thus contributes to NMDA receptor activation. Inhibitors of DAAO were sought testing [6+5] type heterocycles and identified isatin derivatives as micromolar DAAO inhibitors. A pharmacophore and structure-activity relationship analysis of isatins and reported DAAO inhibitors led us to investigate 1H-indazol-3-ol derivatives and nanomolar inhibitors were identified. The series was further characterized by pK_a and isothermal titration calorimetry measurements. Representative compounds exhibited beneficial properties in in vitro metabolic stability and PAMPA assays. 6-fluoro-1H-indazol-3-ol (37) significantly increased plasma d-serine level in an in vivo study on mice. These results show that the 1H-indazol-3-ol series represents a novel class of DAAO inhibitors with the potential to develop drug candidates.     

Discovery of potential pancreatic cholesterol esterase inhibitors using pharmacophore modelling, virtual screening, and optimization studies

Pancreatic cholesterol esterase (CEase) is a serine hydrolase involved in the hydrolysis of variety of lipids and transport of free cholesterol. In this study, pharmacophore hypotheses based on known inhibitors were generated using common feature pharmacophore generation protocol available in Discovery Studio program. The best pharmacophore model containing two hydrogen bond acceptor and three hydrophobic features was selected and validated. It was further used in screening three diverse chemical databases. Hit compounds were subjected to drug-likeness and molecular docking studies. Four hits, namely SEW00846, NCI0040784, GK03167, and CD10645, were selected based on the GOLD fitness score and interaction with active site amino acids. All hit compounds were further optimized to improve their binding in the active site. The optimized compounds were found to have improved binding at the active site. Strongly binding optimized hits at the active site can act as virtual leads in potent CEase inhibitor designing.     

Is rat an appropriate animal model to study the involvement of D-serine catabolism in schizophrenia? Insights from characterization of D-amino acid oxidase

D-Amino acid oxidase (DAAO; EC1.4.3.3) has been proposed to play a main role in the degradation of D-serine, an allosteric activator of the N-methyl-D-aspartate-type glutamate receptor in the human brain, and to be associated with the onset of schizophrenia. To prevent excessive D-serine degradation, novel drugs for schizophrenia treatment based on DAAO inhibition were designed and tested on rats. However, the properties of rat DAAO are unknown and various in vivo trials have demonstrated the effects of DAAO inhibitors on d-serine concentration in rats. In the present study, rat DAAO was efficiently expressed in Escherichia coli. The recombinant enzyme was purified as an active, 40 kDa monomeric flavoenzyme showing the basic properties of the dehydrogenase-oxidase class of flavoproteins. Rat DAAO differs significantly from the human counterpart because: (a) it possesses a different substrate specificity; (b) it shows a lower kinetic efficiency, mainly as a result of a low substrate affinity; (c) it differs in affinity for the binding of classical inhibitors; (d) it is a stable monomer in the absence of an active site ligand; and (e) it interacts with the mammalian protein modulator pLG72 yielding a ~100 kDa complex in addition to the ~200 kDa one, as formed by the human DAAO. Furthermore, the concentration of endogenous D-serine in U87 glioblastoma cells was not affected by transfection with rat DAAO, whereas it was significantly decreased when expressing the human homologue. These results raise doubt on the use of the rat as a model system for testing new drugs against schizophrenia and indicate a different physiological function of DAAO in rodents and humans.     

Computational approach to the identification of novel Aurora-A inhibitors

Aurora kinase A has been emerging as a key therapeutic target for the design of anticancer drugs. For the purpose of finding biologically active and novel compounds and providing new ideas for drug-design, we performed virtual screening using commercially available databases. A three-dimensional common feature pharmacophore model was developed with the HipHop program provided in the Catalyst software package, and this model was used as a query for screening the databases. A recursive partitioning (RP) model was developed as a filtering system, which was able to classify active and inactive compounds. Eventually, a step-wise virtual screening procedure was conducted by applying the common feature pharmacophore and the RP model in succession to discover novel potent Aurora-A inhibitors. A total of 68 compounds were selected for testing of their in vitro anticancer activities against various human cancer cell lines. Based on the activity data, we have identified fifteen compounds that warrant further investigation. Several compounds have a high inhibition rate (above 80% at 10 ??M) and a GI₅₀ lower than 5 ??M for the cell lines DU145 and HT29. Enzyme assay for these compounds identified hits with micro molar activity. Compound C11 has the highest activity (IC₅₀ = 5.09 ??M). The hits obtained from this screening scheme could be potential drug candidates after further optimization.     

Exploration of Virtual Candidates for Human HMG-CoA Reductase Inhibitors Using Pharmacophore Modeling and Molecular Dynamics Simulations

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.     

Discovery of Novel New Delhi Metallo-β-Lactamases-1 Inhibitors by Multistep Virtual Screening

The emergence of NDM-1 containing multi-antibiotic resistant "Superbugs" necessitates the needs of developing of novel NDM-1inhibitors. In this study, we report the discovery of novel NDM-1 inhibitors by multi-step virtual screening. From a 2,800,000 virtual drug-like compound library selected from the ZINC database, we generated a focused NDM-1 inhibitor library containing 298 compounds of which 44 chemical compounds were purchased and evaluated experimentally for their ability to inhibit NDM-1 in vitro. Three novel NDM-1 inhibitors with micromolar IC₅₀ values were validated. The most potent inhibitor, VNI-41, inhibited NDM-1 with an IC₅₀ of 29.6 ± 1.3 μM. Molecular dynamic simulation revealed that VNI-41 interacted extensively with the active site. In particular, the sulfonamide group of VNI-41 interacts directly with the metal ion Zn1 that is critical for the catalysis. These results demonstrate the feasibility of applying virtual screening methodologies in identifying novel inhibitors for NDM-1, a metallo-β-lactamase with a malleable active site and provide a mechanism base for rational design of NDM-1 inhibitors using sulfonamide as a functional scaffold.     

Effect of varying chain length between P1 and P1′ position of tripeptidomimics on activity of angiotensin-converting enzyme inhibitors

One of the efficient modes of treatments of chronic hypertension and cardiovascular disorders has been to restrain the formation of angiotensin-II by inhibiting the action of angiotensin-converting enzyme (ACE) on angiotensin-I. The efforts continue towards achieving superior molecules or drugs with improved affinities, better bioavailability and thus to achieve long duration of action with minimum side effects. Previously, we reported a library of tripeptidomimics of Ornithyl–Proline (Orn–Pro) conjugated with various unnatural amino acids and carboxylic acid derived heterocyclics based on the SAR studies of existing ACE inhibitors. Their synthesis and screening for possible inhibitors of angiotensin-converting enzyme (ACE) revealed that increase in the backbone chain length by one carbon atom results in a sudden decrease in their activity. Therefore, in the present study heterocycles with different chain length were introduced to interact with the hydrophobic S₂ sub-site of ACE and screened for their in vitro ACE inhibition activity. Further, their binding interaction with C-domain of somatic ACE was also determined. Docking and consequent LUDI scores showed good correlation with binding of these molecules in the active site of ACE. Results suggest that heterocycles with C₃ chain length are more appropriate for the effective binding of the tripeptidomimics within the active site of ACE.     

Discovery of a natural PI3Kδ inhibitor through virtual screening and biological assay study

Phosphoinositide-3-kinase-δ (PI3Kδ) is a key regulator in the process of IgE mediated mast cell degranulation, which directly induces allergic diseases, such as asthma. This study is aimed at discovery of natural PI3Kδ inhibitors from Chinese medicine and evaluating their anti-mast cell degranulation activity. A combined virtual screening based on 3D pharmacophore model and molecular docking was used to screen for bioactive ingredients directly targeting PI3Kδ. Then, an in vitro kinase inhibition assay was conducted to evaluate the PI3Kδ inhibitory activity of the virtual screening hits. Subsequently, a β-hexosaminidase release assay was performed to verify the anti-mast cell degranulation activity of the active compounds. Finally, ginkgoneolic acid was identified as a PI3Kδ inhibitor (IC₅₀?=?2.49?μM) and exhibited anti-mast cell degranulation activity in vitro (IC₅₀?=?2.40?μM). Docking studies showed that Glu826, Val827 and Val828 were key amino acid residues for PI3Kδ inhibitory activity. Ginkgoneolic acid may be a potential lead compound for developing effective and safe PI3Kδ-inhibiting drugs.