Discovery of new inhibitors of aldose reductase from molecular docking and database screening

Aldose reductase (ALR2) is a target enzyme for the treatment of diabetic complications. Owing to the limited number of currently available drugs for the treatment of diabetic complications, the discovery of new inhibitors of ALR2 that can potentially be optimized as drugs appears highly desirable. In this study, a molecular docking analysis of the structures of more than 127,000 organic compounds contained in the National Cancer Institute database was performed to find and score molecules that are complementary to ALR2. Besides retrieving several carboxylic acid derivatives, which are known to generally inhibit aldose reductase, docking proposed other families of putative inhibitors such as sulfonic acids, nitro-derivatives, sulfonamides and carbonyl derivatives. Twenty-five compounds, chosen as the highest-scoring representatives of each of these families, were tested as aldose reductase inhibitors. Five of them were found to inhibit aldose reductase in the micromolar range. For these active compounds, selectivity with respect to the closely-related aldehyde reductase was determined by measuring the corresponding inhibitory activities. The structures of the complexes between the new lead inhibitors and aldose reductase, here refined with molecular mechanics and molecular dynamics calculations, suggest that new pharmacophoric groups can bind aldose reductase very efficiently. In the case of the family of the nitro-derivative inhibitors, a class of particularly interesting compounds, a round of optimizations was performed with the synthesis and biological evaluation of a series of derivatives aimed at testing the proposed binding mode and at improving interaction with active site residues. Starting from a hit compound having an IC(50) of 42 microM, the most potent compound synthesized showed a 10-fold increase in inhibitory activity and 10-fold selectivity with respect to ALR1, and structure--activity relationships of the designed compounds were in agreement with the proposed mode of binding at the active site.     

Structure-based design of inhibitors of human L-xylulose reductase modelled into the active site of the enzyme

The program GRID was used to design potential inhibitors of human L-xylulose reductase based on a model of the holoenzyme in complex with n-butyric acid. The inclusion of phosphate or carboxylate functional groups in the ligand suggested an increase in the net binding energy of the complex up to 2.8- and 4.0-fold, respectively. This study may be useful in the development of potent and specific inhibitors of the enzyme.     

Structure-based virtual screening,molecular docking,ADMET and molecular simulations to develop benzoxaborole analogs as potential inhibitor against Leishmania donovani trypanothione reductase

Visceral leishmaniasis (VL) is the most fatal form of leishmaniasis and it affects 70 countries worldwide. Increasing drug resistant for antileishmanial drugs such as miltefosine, sodium stibogluconate and pentamidine has been reported in the VL endemic region. Amphotericin B has shown potential antileishmanial activity in different formulations but its cost of treatment and associated nephrotoxicity have limited its use by affected people living in the endemic zone. To control the VL infection in the affected countries, it is necessary to develop new antileishmanial compounds with high efficacy and negligible toxicity. Computer aided programs such as binding free energy estimation; ADMET prediction and molecular dynamics simulation can be used to investigate novel antileishmanial molecules in shorter duration. To develop antileishmanial lead molecule, we performed standard precision (SP) docking for 1160 benzoxaborole analogs along with reference inhibitors against trypanothione reductase of Leishmania parasite. Furthermore, extra precision (XP) docking, ADMET prediction, prime MM-GBSA was conducted over 115 ligands, showing better docking score than reference inhibitors to get potential antileishmanial compounds. Simultaneously, area under the curve (AUC) was estimated using ROC plot to validate the SP and XP docking protocol. Later on, two benzoxaborole analogs with best MM-GBSA ΔG-bind were subjected to molecular simulation and docking confirmation to ensure the ligand interaction with TR. The presented drug discovery based on computational study confirms that BOB27 can be used as a potential drug candidate and warrants further experimental investigation to fight against VL in endemic areas.     

Structure-based virtual screening approach to the discovery of novel PTPMT1 phosphatase inhibitors

Dual-specificity protein tyrosine phosphatase localized to mitochondrion 1 (PTPMT1) has recently proved to be a promising therapeutic target for the treatment of type II diabetes. Herein we report the first example for a successful application of the structure-based virtual screening to identify the novel inhibitors of human PTPMT1. These inhibitors were computationally screened for having desirable physicochemical properties as a drug candidate and reveal a high potency with IC(50) values ranging from 0.7 to 17.3μM. Therefore, they deserve consideration for further development by structure-activity relationship studies to optimize the antidiabetic activities. Structural features relevant to the stabilization of the newly identified inhibitors in the active site of PTPMT1 are addressed in detail.     

Structure-based virtual screening approach to the discovery of phosphoinositide 3-kinase alpha inhibitors

Phosphoinositide 3-kinase alpha (PI3Kα) has proved to be an attractive target for the development of therapeutics for the treatment of cancer. Herein we report a successful application of the structure-based virtual screening to identify the novel inhibitors of PI3Kα. These inhibitors have desirable physicochemical properties as a drug candidate and reveal a moderate potency with IC₅₀ values ranging from 20 to 40 μM. Therefore, they deserve a consideration for further development by structure-activity relationship (SAR) studies to optimize the inhibitory activities. Structural features relevant to the stabilization of the newly identified inhibitors in the ATP-binding site of PI3Kα are addressed in detail.     

Structure-based virtual screening approach to the discovery of p38 MAP kinase inhibitors

p38 Mitogen-activated protein kinase (MAPK) has been considered to be a promising target for the development of therapeutics for various immunologic diseases. Herein we report an example for a successful application of the virtual screening with protein-ligand docking to identify the novel inhibitors of p38α MAPK. These inhibitors were screened for having desirable physicochemical properties as a drug candidate and compound 1-3 revealed a moderate inhibitory activity with IC(50) values ranging from 0.7 to 20 μM. Therefore, they deserve a consideration for further development by structure-activity relationship (SAR) studies to optimize the inhibitory activities. Structural features relevant to the stabilization of the newly identified inhibitors in the ATP-binding site of p38 MAPK are addressed in detail.     

Structure-based drug discovery of carbonic anhydrase inhibitors

Inhibition of the metalloenzyme carbonic anhydrase (CA; EC 4.2.1.1) has pharmacologic applications in the field of anti-glaucoma, anti-convulsant and anti-cancer agents. But recently, it has also emerged that these enzymes have the potential for designing anti-infective drugs (anti-fungal and anti-bacterial agents) with a novel mechanism of action. Sulphonamides and their isosteres (sulphamates/sulphamides) constitute the main class of CA inhibitors (CAIs), which bind to the metal ion from the enzyme active site. Recently, the dithiocarbamates (DTCs), possessing a similar mechanism of action, were reported as a new class of inhibitors. These types of CAIs will be discussed in detail in this review. Novel drug design strategies have been reported ultimately based on the tail approach for obtaining sulphonamides/DTCs, which exploit more external binding regions within the enzyme active site (in addition to coordination to the metal ion), leading thus to isoform-selective compounds. Most of the promising data have been obtained by combining x-ray crystallography of enzyme-inhibitor adducts with novel synthetic approaches for generating chemical diversity. Whereas sulphonamide – NO donating hybrid drugs were reported as effective anti-glaucoma agents, most of the interesting new inhibitors were designed for inhibiting specifically the tumour-associated isoforms CA IX and XII, validated targets for imaging and treatment of hypoxic tumours. Promising compounds that inhibit CAs from bacterial and fungal pathogens, of the DTC and carboxylate types, will be also reviewed.     

Structure-based discovery of novel flavonol inhibitors of human protein kinase CK2

Serine/threonine protein kinase CK2 controls vast variety of fundamental processes in cell life; however, despite long period of study, its functional role is not completely determined. CK2 has a significant pathogenic potential and its activity is strictly associated with the development of various kinds of disorders. There are a growing number of facts that inhibitors of CK2 could be used as pharmaceutical agents for the cancer treatment, viral infections, and inflammatory diseases. In this article, we report structural and biological data on the novel synthetic flavonol derivatives, 3-hydroxy-4'-carboxyflavones, possessing a high inhibitory activity toward CK2. With the aid of combinatorial organic synthesis, molecular modeling techniques and biochemical in vitro tests, we studied the structure-activity relationships of flavonol derivatives and developed binding model describing their key intermolecular interactions with the CK2 ATP-binding site. Obtained data show that the synthetic 3-hydroxy-4'-carboxyflavones possess the highest activity among flavonol inhibitors of CK2 known till date.     

Structure-based discovery of cellular-active allosteric inhibitors of FAK

In order to develop potent and selective focal adhesion kinase (FAK) inhibitors, synthetic studies on pyrazolo[4,3-c][2,1]benzothiazines targeted for the FAK allosteric site were carried out. Based on the X-ray structural analysis of the co-crystal of the lead compound, 8-(4-ethylphenyl)-5-methyl-1,5-dihydropyrazolo[4,3-c][2,1]benzothiazine 4,4-dioxide 1 with FAK, we designed and prepared 1,5-dimethyl-1,5-dihydropyrazolo[4,3-c][2,1]benzothiazin derivatives which selectively inhibited kinase activity of FAK without affecting seven other kinases. The optimized compound, N-(4-tert-butylbenzyl)-1,5-dimethyl-1,5-dihydropyrazolo[4,3-c][2,1]benzothiazin-8-amine 4,4-dioxide 30 possessed significant FAK kinase inhibitory activities both in cell-free (IC₅₀ = 0.64 μM) and in cellular assays (IC₅₀ = 7.1 μM). These results clearly demonstrated a potential of FAK allosteric inhibitors as antitumor agents.     

Virtual screening,docking, and dynamics of potential new inhibitors of dihydrofolate reductase from Yersinia pestis

In the present work, we propose to design drugs that target the enzyme dihydrofolate redutase (DHFR) as a means of a novel drug therapy against plague. Potential inhibitors of DHFR from Yersinia pestis (YpDHFR) were selected by virtual screening and subjected to docking, molecular dynamics (MD) simulations, and Poisson–Boltzmann surface area method, in order to evaluate their interactions in the active sites of YpDHFR and human DHFR (HssDHFR). The results suggested selectivity for three compounds that were further used to propose the structures of six new potential selective inhibitors for YpDHFR.     

Virtual screening for inhibitors of human aldose reductase

The inhibition of aldose reductase (AR) provides an interesting strategy to prevent the complications of chronic diabetes. Although a large number of different AR inhibitors are known, very few of these compounds exhibit sufficient efficacy in clinical trials. We performed a virtual screening based on the ultrahigh resolution crystal structure of the inhibitor IDD594 in complex with human AR. AR operates on a large scale of structurally different substrates. To achieve this pronounced promiscuity, the enzyme can adapt rather flexibly to its substrates. Likewise, it has a similar adaptability for the binding of inhibitors. We applied a protocol of consecutive hierarchical filters to search the Available Chemicals Directory. In the first selection step, putative ligands were chosen that exhibit functional groups to anchor the anion-binding pocket of AR. Subsequently, a pharmacophore model based on the binding geometry of IDD594 and the mapping of the binding pocket in terms of putative "hot spots" of binding was applied as a second consecutive filter. In a third and final filtering step, the remaining candidate molecules were flexibly docked into the binding pocket of IDD594 with FlexX and ranked according to their estimated DrugScore values. Out of 206 compounds selected by this search and complemented by a cluster analysis and visual inspection, 9 compounds were selected and subjected to biological testing. Of these, 6 compounds showed IC50 values in the micromolar range. According to the proposed binding mode, the two inhibitors BTB02809 (IC50 = 2.4 +/- 0.5 microM) and JFD00882 (IC50 = 4.1 +/- 1.0 microM) both place a nitro group into the hydrophobic specificity pocket of human AR in an orientation coinciding with the position of the bromine atom of IDD594. The interaction of this Br with Thr113 has been identified as a key feature that is responsible for selectivity enhancement.     

The discovery of potential acetylcholinesterase inhibitors: A combination of pharmacophore modeling,virtual screening,and molecular docking studies

Background  

Alzheimer's disease (AD) is the most common cause of dementia characterized by progressive cognitive impairment in the elderly people. The most dramatic abnormalities are those of the cholinergic system. Acetylcholinesterase (AChE) plays a key role in the regulation of the cholinergic system, and hence, inhibition of AChE has emerged as one of the most promising strategies for the treatment of AD.     

Lead discovery using molecular docking

As the structures of more and more proteins and nucleic acids become available, molecular docking is increasingly considered for lead discovery. Recent studies consider the hit-rate enhancement of docking screens and the accuracy of docking structure predictions. As more structures are determined experimentally, docking against homology-modeled targets also becomes possible for more proteins. With more docking studies being undertaken, the 'drug-likeness' and specificity of docking hits is also being examined.     

Structure-based virtual screening of Src kinase inhibitors

Src is an important target in multiple processes associated with tumor growth and development, including proliferation, neovascularization, and metastasis. In this study, hit identification was performed by virtual screening of commercial and in-house compound libraries. Docking studies for the hits were performed, and scoring functions were used to evaluate the docking results and to rank ligand-binding affinities. Subsequently, hit optimization for potent and selective candidate Src inhibitors was performed through focused library design and docking analyses. Consequently, we report that a novel compound ‘43’ with an IC₅₀ value of 89 nM, representing (S)-N-(4-(5-chlorobenzo[d][1,3]dioxol-4-ylamino)-7-(2-methoxyethoxy)quinazolin-6-yl)pyrrolidine-2-carboxamide, is highly selective for Src in comparison to EGFR (IC₅₀ ratio > 80-fold) and VEGFR-2 (IC₅₀ ratio > 110-fold). Compound 43 exerted anti-proliferative effects on Src-expressing PC3 human prostate cancer and A431 human epidermoid carcinoma cells, with calculated IC₅₀ values of 1.52 and 0.78 μM, respectively. Moreover, compound 43 (0.1 μM) suppressed the phosphorylation of extracellular signal-regulated kinases and p90 ribosomal S6 kinase, downstream molecules of Src, in a time-dependent manner, in both PC3 and A431 cell lines. The docking structure of compound 43 with Src disclosed that the chlorobenzodioxole moiety and pyrrolidine ring of C-6 quinazoline appeared to fit tightly into the hydrophobic pocket of Src. Additionally, the pyrrolidine NH forms a hydrogen bond with the carboxyl group of Asp348. These results confirm the successful application of virtual screening studies in the lead discovery process, and suggest that our novel compound 43 can be an effective Src inhibitor candidate for further lead optimization.     

The discovery of steroids and other novel FKBP inhibitors using a molecular docking program

The molecular docking computer program SANDOCK was used to screen small molecule three-dimensional databases in the hunt for novel FKBP inhibitors. Spectroscopic measurements confirmed binding of over 20 compounds to the target protein, some with dissociation constants in the low micromolar range. The discovery that FK506 binding protein is a steroid binding protein may be of wider biological significance. Two-dimensional NMR was used to determine the steroid binding mode and confirmed the interactions predicted by the docking program.     

Computational screening of dipeptidyl peptidase IV inhibitors from micoroalgal metabolites by pharmacophore modeling and molecular docking

Dipeptidyl peptidase IV (DPP‐IV) catalyzes conversion of GLP‐1 (glucagon like peptide 1) to inert structure, which results in insufficient secretion of insulin and increase in postprandial blood glucose level. The present study attempts to identify novel inhibitors from bioactive metabolites present in microalgae against DPP‐IV through virtual screening, molecular docking, and pharmacophore modeling for the active target. Possible binding modes of all 60 ligands against DPP‐IV receptor were constructed using MTiOpenScreen virtual screening server. Pharmacophore model was built based on identified 38 DPP‐IV test ligands by using the web‐based PharmaGist program which encompasses hydrogen‐bond acceptors, hydrophobic groups, spatial features, and aromatic rings. The pharmacophore model having highest scores was selected to screen active DPP‐IV ligands. Highest scoring model was used as a query in ZincPharmer screening. All identified ligands were filtered, based on the Lipinski's rule‐of‐five and were subjected to docking studies. In the process of docking analyses, we considered different bonding modes of one ligand with multiple active cavities of DPP‐IV with the help of AutoDock 4.0. The docking analyses indicate that the bioactive constituents, namely, β‐stigmasterol, barbamide, docosahexaenoic acid, arachidonic acid, and harman showed the best binding energies on DPP‐IV receptor and hydrogen bonding with ASP545, GLY741, TYR754, TYR666, ARG125, TYR547, SER630, and LYS554 residues. This study concludes that docosahexaenoic acid, arachidonic acid, β‐stigmasterol, barbamide, harman, ZINC58564986, ZINC56907325, ZINC69432950, ZINC69431828, ZINC73533041, ZINC84287073, ZINC69849395, and ZINC10508406 act as possible DPP‐IV inhibitors.     

Structure-based tailoring of compound libraries for high-throughput screening: discovery of novel EphB4 kinase inhibitors

High-throughput docking is a computational tool frequently used to discover small-molecule inhibitors of enzymes or receptors of known three-dimensional structure. Because of the large number of molecules in chemical libraries, automatic procedures to prune multimillion compound collections are useful for high-throughput docking and necessary for in vitro screening. Here, we propose an anchor-based library tailoring approach (termed ALTA) to focus a chemical library by docking and prioritizing molecular fragments according to their binding energy which includes continuum electrostatics solvation. In principle, ALTA does not require prior knowledge of known inhibitors, but receptor-based pharmacophore information (hydrogen bonds with the hinge region) is additionally used here to identify molecules with optimal anchor fragments for the ATP-binding site of the EphB4 receptor tyrosine kinase. The 21,418 molecules of the focused library (from an initial collection of about 730,000) are docked into EphB4 and ranked by force-field-based energy including electrostatic solvation. Among the 43 compounds tested in vitro, eight molecules originating from two different anchors show low-micromolar activity in a fluorescence-based enzymatic assay. Four of them are active in a cell-based assay and are potential anti-angiogenic compounds.     

Drug discovery for DNA break repair system by screening from TCM database and molecular dynamics approach

In search for novel anti-tumour agents targeting non-homologous end joining pathway, we conducted a virtual screening using traditional Chinese medicine (TCM) on Ku heterodimer which plays an important role on DNA break repair system. Docking results gave glycyrrhizic acid and macedonoside C as potential TCM candidates. Both glycyrrhizic acid and macedonoside C were docked onto the Ku heterodimer with high-predicted binding affinities as evidenced in their Dock Score. These two ligands also interact with key residues which have been previously shown to influence Ku DNA-binding affinity. Both compounds show consistent hydrogen bonding interactions with key residues throughout the 10- ns dynamics simulation run. Furthermore, glycyrrhizic acid and macedonoside C are able to form additional hydrogen bond interactions with positively charged Ku heterodimer surface. Such interactions strengthen the binding interaction of the top two TCM compounds with Ku heterodimer. Glycyrrhizic acid and macedonoside C are products of licorice (Glycyrrhiza glabra), which is a commonly used herb in TCM formulations. Overall, we report glycyrrhizic acid and macedonoside C as potential anti-tumour agents.     

Structure-based virtual screening of chemical libraries for drug discovery

One of the main goals in drug discovery is to identify new chemical entities that have a high likelihood of binding to the target protein to elicit the desired biological response. To this end, virtual screening is being increasingly used as a complement to high-throughput screening to improve the speed and efficiency of the drug discovery and development process. The availability of inexpensive high-performance computing platforms in recent years has transformed this field into one that is highly diverse and rapidly evolving, where large chemical databases have been successfully screened to identify hits for a wide range of targets such as Bcl-2 family proteins, G protein-coupled receptors, kinases, metalloproteins, nuclear hormone receptors, proteases and many more.