Virtual screening and bioassay study of novel inhibitors for dengue virus mRNA cap (nucleoside-2'O)-methyltransferase

We report high-throughput structure-based virtual screening of putative Flavivirus 2′-O-methyltransferase inhibitors together with results from subsequent bioassay tests of selected compounds. Potential inhibitors for the S-adenosylmethionine binding site were explored using 2D similarity searching, pharmacophore filtering and docking. The inhibitory activities of 15 top-ranking compounds from the docking calculations were tested on a recombinant methyltransferase with the RNA substrate ^7MeGpppAC₅. Local and global docking simulations were combined to estimate the ligand selectivity for the target site. The results of the combined computational and experimental screening identified a novel inhibitor, with a previously unknown scaffold, that has an IC₅₀ value of 60 μM.     

Structure-based virtual screening and biological evaluation of potent and selective ADAM12 inhibitors

We describe a series of potent and selective inhibitors of ADAM12 that were discovered using computational screening of a focused virtual library. The initial structure-based virtual screening selected 64 compounds from a 3D database of 67,062 molecules. Being evaluated by a cell-based ADAM12 activity assay, compounds 5, 11, 14, 16 were further identified as the potent and selective inhibitors of ADAM12 with low nanomolar IC₅₀ values. The mechanism underlying the potency and selectivity of a representative compound, 5, was investigated through molecular docking studies.     

Lead generation of heat shock protein 90 inhibitors by a combination of fragment-based approach, virtual screening, and structure-based drug design

Heat shock protein 90 (Hsp90) is a molecular chaperone which regulates maturation and stabilization of its substrate proteins, known as client proteins. Many client proteins of Hsp90 are involved in tumor progression and survival and therefore Hsp90 can be a good target for developing anticancer drugs. With the aim of efficiently identifying a new class of orally available inhibitors of the ATP binding site of this protein, we conducted fragment screening and virtual screening in parallel against Hsp90. This approach quickly identified 2-aminotriazine and 2-aminopyrimidine derivatives as specific ligands to Hsp90 with high ligand efficiency. In silico evaluation of the 3D X-ray Hsp90 complex structures of the identified hits allowed us to promptly design CH5015765, which showed high affinity for Hsp90 and antitumor activity in human cancer xenograft mouse models.     

Discovery and development of heat shock protein 90 inhibitors

Heat shock protein 90 (Hsp90) is an important target in cancer because of its role in maintaining transformation and has recently become the focus of several drug discovery and development efforts. While compounds with different modes of action are known, the focus of this review is on those classes of compounds which inhibit Hsp90 by binding to the N-terminal ATP pocket. These include natural product inhibitors such as geldanamycin and radicicol and synthetic inhibitors comprised of purines, pyrazoles, isoxazoles and other scaffolds. The synthetic inhibitors have been discovered either by structure-based design, high throughput screening and more recently using fragment-based design and virtual screening techniques. This review will discuss the discovery of these different classes, as well as their development as potential clinical agents.     

A fluorescence polarization assay for inhibitors of Hsp90

Hsp90 encodes a ubiquitous molecular chaperone protein conserved among species which acts on multiple substrates, many of which are important cell-signaling proteins. Inhibition of Hsp90 function has been promoted as a mechanism to degrade client proteins involved in tumorigenesis and disease progression. Several assays to monitor inhibition of Hsp90 function currently exist but are limited in their use for a drug discovery campaign. Using data from the crystal structure of an initial hit compound, we have developed a fluorescence polarization assay to monitor binding of compounds to the ATP-binding site of Hsp90. This assay is very robust (Z' > 0.9) and can detect affinity of compounds with IC50s to 40 nM. We have used this assay in conjunction with cocrystal structures of small molecules to drive a structure-based design program aimed at the discovery and optimization of a novel class of potent Hsp90 inhibitors.     

Computational aided mechanistic understanding of Camellia sinensis bioactive compounds against co-chaperone p23 as potential anticancer agent

Co-chaperon p23 has been well established as molecular chaperon for the heat shock protein 90 (Hsp90) that further leads to immorality in cancer cells by providing defense against Hsp90 inhibitors, and as stimulating agent for generating overexpressed antiapoptotic proteins, that is, Hsp70 and Hsp27. The natural compounds such as catechins from Camellia sinensis (green tea) are also well known for inhibition activity against various cancer. However, molecular interaction profile and potential lead bioactive compounds against co-chaperon p23 from green tea are not yet reported. To this context, we study the various secondary metabolites of green tea against co-chaperon p23 using structure-based virtual screening from Traditional Chinese Medicine (TCM) database. Following 26 compounds were obtained from TCM database and further studied for extra precision molecular docking that showed binding score between −10.221 and −2.276 kcal/mol with co-chaperon p23. However, relative docking score to known inhibitors, that is, ailanthone (−4.54 kcal/mol) and gedunin ( 3.60 kcal/mol) along with ADME profile analysis concluded epicatechin (−7.013 kcal/mol) and cis-theaspirone (−4.495 kcal/mol) as potential lead inhibitors from green tea against co-chaperone p23. Furthermore, molecular dynamics simulation and molecular mechanics generalized born surface area calculations validated that epicatechin and cis-theaspirone have significantly occupied the active region of co-chaperone p23 by hydrogen and hydrophobic interactions with various residues including most substantial amino acids, that is, Thr90, Ala94, and Lys95. Hence, these results supported the fact that green tea contained potential compounds with an ability to inhibit the cancer by disrupting the co-chaperon p23 activity.     

In silico discovery of 2-amino-4-(2,4-dihydroxyphenyl)thiazoles as novel inhibitors of DNA gyrase B

Cyclothialidines are a class of bacterial DNA gyrase B (GyrB) subunit inhibitors, targeting its ATP-binding site. Starting from the available structural information on cyclothialidine GR122222X (2), an in silico virtual screening campaign was designed combining molecular docking calculations with three-dimensional structure-based pharmacophore information. A novel class of 2-amino-4-(2,4-dihydroxyphenyl)thiazole based inhibitors (5–9) with low micromolar antigyrase activity was discovered.     

Discovery of the inhibitors of tumor necrosis factor alpha with structure-based virtual screening

Tumor necrosis factor alpha (TNF-α) has been considered as one of the attractive drug targets for allergic diseases including asthma. We have been able to identify five novel TNF-α inhibitors with a drug-design protocol involving the structure-based virtual screening and in vitro cell-based assay for antagonistic activity. Because the newly discovered inhibitors are structurally diverse and have the desirable physicochemical properties as a drug candidate, they deserve a further investigation as anti-asthmatic drugs. The interactions of the identified inhibitors in the binding site of TNF-α dimer are addressed in detail to understand the mechanisms for the stabilization of the inactive dimeric form of TNF-α.     

Dynamic and Multi-Pharmacophore Modeling for Designing Polo-Box Domain Inhibitors

The polo-like kinase 1 (Plk1) is a critical regulator of cell division that is overexpressed in many types of tumors. Thus, a strategy in the treatment of cancer has been to target the kinase activity (ATPase domain) or substrate-binding domain (Polo-box Domain, PBD) of Plk1. However, only few synthetic small molecules have been identified that target the Plk1-PBD. Here, we have applied an integrative approach that combines pharmacophore modeling, molecular docking, virtual screening, and in vitro testing to discover novel Plk1-PBD inhibitors. Nine Plk1-PBD crystal structures were used to generate structure-based hypotheses. A common pharmacophore model (Hypo1) composed of five chemical features was selected from the 9 structure-based hypotheses and used for virtual screening of a drug-like database consisting of 159,757 compounds to identify novel Plk1-PBD inhibitors. The virtual screening technique revealed 9,327 compounds with a maximum fit value of 3 or greater, which were selected and subjected to molecular docking analyses. This approach yielded 93 compounds that made good interactions with critical residues within the Plk1-PBD active site. The testing of these 93 compounds in vitro for their ability to inhibit the Plk1-PBD, showed that many of these compounds had Plk1-PBD inhibitory activity and that compound Chemistry_28272 was the most potent Plk1-PBD inhibitor. Thus Chemistry_28272 and the other top compounds are novel Plk1-PBD inhibitors and could be used for the development of cancer therapeutics.     

Targeting conserved water molecules: Design of 4-aryl-5-cyanopyrrolo[2,3-d]pyrimidine Hsp90 inhibitors using fragment-based screening and structure-based optimization

Inhibitors of the Hsp90 molecular chaperone are showing promise as anti-cancer agents. Here we describe a series of 4-aryl-5-cyanopyrrolo[2,3-d]pyrimidine ATP competitive Hsp90 inhibitors that were identified following structure-driven optimization of purine hits revealed by NMR based screening of a proprietary fragment library. Ligand-Hsp90 X-ray structures combined with molecular modeling led to the rational displacement of a conserved water molecule leading to enhanced affinity for Hsp90 as measured by fluorescence polarization, isothermal titration calorimetry and surface plasmon resonance assays. This displacement was achieved with a nitrile group, presenting an example of efficient gain in binding affinity with minimal increase in molecular weight. Some compounds in this chemical series inhibit the proliferation of human cancer cell lines in vitro and cause depletion of oncogenic Hsp90 client proteins and concomitant elevation of the co-chaperone Hsp70. In addition, one compound was demonstrated to be orally bioavailable in the mouse. This work demonstrates the power of structure-based design for the rapid evolution of potent Hsp90 inhibitors and the importance of considering conserved water molecules in drug design.     

2,4-dihydroxy benzaldehyde derived Schiff bases as small molecule Hsp90 inhibitors: Rational identification of a new anticancer lead

Hsp90 is a molecular chaperone that heals diverse array of biomolecules ranging from multiple oncogenic proteins to the ones responsible for development of resistance to chemotherapeutic agents. Moreover they are over-expressed in cancer cells as a complex with co-chaperones and under-expressed in normal cells as a single free entity. Hence inhibitors of Hsp90 will be more effective and selective in destroying cancer cells with minimum chances of acquiring resistance to them. In continuation of our goal to rationally develop effective small molecule azomethines against Hsp90, we designed few more compounds belonging to the class of 2,4-dihydroxy benzaldehyde derived imines (1–13) with our validated docking protocol. The molecules exhibiting good docking score were synthesized and their structures were confirmed by IR, ¹H NMR and mass spectral analysis. Subsequently, they were evaluated for their potential to suppress Hsp90 ATPase activity by Malachite green assay. The antiproliferative effect of the molecules were examined on PC3 prostate cancer cell lines by adopting 3-(4,5-dimethythiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay methodology. Finally, schiff base 13 emerged as the lead molecule for future design and development of Hsp90 inhibitors as anticancer agents.     

Structure-based discovery of a new class of Hsp90 inhibitors

Docking-based virtual screening identified 1-(2-phenol)-2-naphthol compounds as a new class of Hsp90 inhibitors of low to sub-micromolar potency. Here we report the binding affinities and cellular activities of several members of this class. A high resolution crystal structure of the most potent compound reveals its binding mode in the ATP binding site of Hsp90, providing a rationale for the observed activity of the series and suggesting strategies for developing compounds with improved properties.     

Novel human mPGES-1 inhibitors identified through structure-based virtual screening

Microsomal prostaglandin E synthase-1 (mPGES-1) is an inducible prostaglandin E synthase after exposure to pro-inflammatory stimuli and, therefore, represents a novel target for therapeutic treatment of acute and chronic inflammatory disorders. It is essential to identify mPGES-1 inhibitors with novel scaffolds as new leads or hits for the purpose of drug design and discovery that aim to develop the next-generation anti-inflammatory drugs. Herein we report novel mPGES-1 inhibitors identified through a combination of large-scale structure-based virtual screening, flexible docking, molecular dynamics simulations, binding free energy calculations, and in vitro assays on the actual inhibitory activity of the computationally selected compounds. The computational studies are based on our recently developed three-dimensional (3D) structural model of mPGES-1 in its open state. The combined computational and experimental studies have led to identification of new mPGES-1 inhibitors with new scaffolds. In particular, (Z)-5-benzylidene-2-iminothiazolidin-4-one is a promising novel scaffold for the further rational design and discovery of new mPGES-1 inhibitors. To our best knowledge, this is the first time a 3D structural model of the open state mPGES-1 is used in structure-based virtual screening of a large library of available compounds for the mPGES-1 inhibitor identification. The positive experimental results suggest that our recently modeled trimeric structure of mPGES-1 in its open state is ready for the structure-based drug design and discovery.     

Discovery and biological evaluation of novel alpha-glucosidase inhibitors with in vivo antidiabetic effect

Discovery of alpha-glucosidase inhibitors has been actively pursued with the aim to develop therapeutics for the treatment of diabetes and the other carbohydrate-mediated diseases. We have identified four novel alpha-glucosidase inhibitors by means of a drug design protocol involving the structure-based virtual screening under consideration of the effects of ligand solvation in the scoring function and in vitro enzyme assay. Because the newly identified inhibitors reveal in vivo antidiabetic activity as well as a significant potency with more than 70% inhibition of the catalytic activity of alpha-glucosidase at 50 microM, all of them seem to deserve further development to discover new drugs for diabetes. Structural features relevant to the interactions of the newly identified inhibitors with the active site residues of alpha-glucosidase are discussed in detail.     

Identification of novel inhibitors of extracellular signal-regulated kinase 2 based on the structure-based virtual screening

Extracellular signal-regulated kinase 2 (ERK2) has become an attractive target for the development of therapeutics for the treatment of cancer. We have been able to identify eight new inhibitors of ERK2 by means of a drug design protocol involving the virtual screening with docking simulations and in vitro enzyme assay. The newly discovered inhibitors can be categorized into three structural classes and reveal a significant potency with IC(50) values ranging from 1 to 30 microM. Therefore, all of the three inhibitor scaffolds deserve further development by structure-activity relationship or de novo design methods. Structural features relevant to the stabilizations of the newly identified inhibitors in the ATP-binding site of ERK2 are discussed in detail.     

Discovery of novel EGFR tyrosine kinase inhibitors by structure-based virtual screening

By using of structure-based virtual screening, 13 novel epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors were discovered from 197,116 compounds in the SPECS database here. Among them, 8 compounds significantly inhibited EGFR kinase activity with IC(50) values lower than 10 μM. 3-{[1-(3-Chloro-4-fluorophenyl)-3,5-dioxo-4-pyrazolidinylidene]methyl}phenyl 2-thiophenecarboxylate (13), particularly, was the most potent inhibitor possessing the IC(50) value of 3.5 μM. The docking studies also provide some useful information that the docking models of the 13 compounds are beneficial to find a new path for designing novel EGFR inhibitors.     

Application of docking-based comparative intermolecular contacts analysis to validate Hsp90α docking studies and subsequent in silico screening for inhibitors

Heat shock protein (Hsp90α) has been recently implicated in cancer, prompting several attempts to discover and optimize new Hsp90α inhibitors. Towards this end, we docked 83 diverse Hsp90α inhibitors into the ATP-binding site of this chaperone using several docking-scoring settings. Subsequently, we applied our newly developed computational tool-docking-based comparative intramolecular contacts analysis (dbCICA)-to assess the different docking conditions and select the best settings. dbCICA is based on the number and quality of contacts between docked ligands and amino acid residues within the binding pocket. It assesses a particular docking configuration based on its ability to align a set of ligands within a corresponding binding pocket in such a way that potent ligands come into contact with binding site spots distinct from those approached by low-affinity ligands, and vice versa. The optimal dbCICA models were translated into valid pharmacophore models that were used as 3D search queries to mine the National Cancer Institute's structural database for new inhibitors of Hsp90α that could potentially be used as anticancer agents. The process culminated in 15 micromolar Hsp90α ATPase inhibitors.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

Targeting the hydrophobic region of Hsp90’s ATP binding pocket with novel 1,3,5-triazines

Heat shock protein 90 (Hsp90) is a molecular chaperone that plays an important role in regulating the maturation and stabilization of many oncogenic proteins. In an attempt to discover a new class of Hsp90 inhibitors, a series of 1,3,5-triazine compounds were rationally designed, synthesized, and their biological activities were evaluated. Compound 3b was found to degrade Hsp90’s client proteins of Her2, Met and Akt and to induce the expression level of Hsp70. The binding mode of 3b in the ATP-binding site of Hsp90 was predicted by the molecular docking.     

3D-QSAR,molecular docking,and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.