Inside the Hsp90 inhibitors binding mode through induced fit docking

Molecular docking and dynamics simulation study of quinones and pyrones from Alternaria solani and Alternaria alternata with HSP90: an important therapeutic target of cancer

Although cancer continues to be one of the world's major causes of death, current cancer drugs have many serious side effects. There remains a need for new anticancer agents to overcome these shortcomings. Alternaria is one of the most widespread fungal genera, many species of which produce several classes of metabolites with potential polypharmacological activities. A few quinones and pyrones from Alternaria spp. have proven to exert cytotoxic effects against certain cancer cell lines, but their molecular mode of action is not known. The current study aimed to investigate the potential mechanisms that underlie the anticancer activity of a few selected quinones and pyrones from Alternaria solani and Alternaria alternata by molecular docking and dynamic simulation approaches. The selected metabolites were screened for their binding affinity to Heat shock protein 90 (HSP90), which is a known anticancer drug target. Molecular docking studies have revealed that Macrosporin, Altersolanol B, Fonsecin, and Neoaltenuene have good binding affinities with the target protein and the stabilities of the formed complexes were evaluated through molecular dynamics simulations. By analyzing the Root Mean Square Distance (RMSD), Root Mean Square Fluctuation (RMSF), and Principal Component Analysis (PCA) plots obtained from molecular dynamics simulations, this study shows that the complexes of all 4 lead molecules with target protein are stable over a 100 ns period. Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were used to compute the binding free energies. The lead molecules were studied using in-silico analysis to determine their drug-likeness based on their Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) and physicochemical properties. The results demonstrate that Macrosporin, Fonsecin, and Neoaltenuene could become promising anticancer molecules that target HSP90.Communicated by Ramaswamy H. Sarma.

Mebendazole's Conformational Space and Its Predicted Binding to Human Heat-Shock Protein 90

Recent experimental evidence suggests that mebendazole, a popular antiparasitic drug, binds to heat shock protein 90 (Hsp90) and inhibits acute myeloid leukemia cell growth. In this study we use quantum mechanics (QM), molecular similarity, and molecular dynamics (MD) calculations to predict possible binding poses of mebendazole to the adenosine triphosphate (ATP) binding site of Hsp90. Extensive conformational searches and minimization of the five mebendazole tautomers using the MP2/aug-cc-pVTZ theory level resulted in 152 minima. Mebendazole-Hsp90 complex models were subsequently created using the QM optimized conformations and protein coordinates obtained from experimental crystal structures that were chosen through similarity calculations. Nine different poses were identified from a total of 600 ns of explicit solvent, all-atom MD simulations using two different force fields. All simulations support the hypothesis that mebendazole is able to bind to the ATP binding site of Hsp90.

Synthesis, biological evaluation, and in silico studies of novel chalcone- and pyrazoline-based 1,3,5-triazines as potential anticancer agents

A novel series of triazin-chalcones (7,8)a-g and triazin-N-(3,5-dichlorophenyl)pyrazolines (9,10)a-g were synthesized and evaluated for their anticancer activity against nine different cancer strains. Triazine ketones 5 and 6 were synthesized from the cyanuric chloride 1 by using stepwise nucleophilic substitution of the chlorine atom. These ketones were subsequently subjected to a Claisen-Schmidt condensation reaction with aromatic aldehydes affording chalcones (7,8)a-g. Then, N-(3,5-dichlorophenyl)pyrazolines (9,10)a-g were obtained by cyclocondensation reactions of the respective chalcones (7,8)a-g with 3,5-dichlorophenylhydrazine. Among all the evaluated compounds, chalcones 7d,g and 8g exhibited more potent in vitro anticancer activity, with outstanding GI50 values ranging from 0.422 to 14.9 μM and LC50 values ranging from 5.08 μM to >100 μM. In silico studies, for both ligand- and structure-based, were executed to explore the inhibitory nature of chalcones and triazine derivatives. The results suggested that the evaluated compounds could act as modulators of the human thymidylate synthase enzyme.

Molecular Docking: Challenges, Advances and its Use in Drug Discovery Perspective

Molecular docking is a process through which small molecules are docked into the macromolecular structures for scoring its complementary values at the binding sites. It is a vibrant research area with dynamic utility in structure-based drug-designing, lead optimization, biochemical pathway and for drug designing being the most attractive tools. Two pillars for a successful docking experiment are correct pose and affinity prediction. Each program has its own advantages and drawbacks with respect to their docking accuracy, ranking accuracy and time consumption so a general conclusion cannot be drawn. Moreover, users don't always consider sufficient diversity in their test sets which results in certain programs to outperform others. In this review, the prime focus has been laid on the challenges of docking and troubleshooters in existing programs, underlying algorithmic background of docking, preferences regarding the use of docking programs for best results illustrated with examples, comparison of performance for existing tools and algorithms, state of art in docking, recent trends of diseases and current drug industries, evidence from clinical trials and post-marketing surveillance are discussed. These aspects of the molecular drug designing paradigm are quite controversial and challenging and this review would be an asset to the bioinformatics and drug designing communities.

The in silico identification of potent anti-cancer agents by targeting the ATP binding site of the N-domain of HSP90

To identify new HSP90 inhibitors, the ATP binding site of the N-domain of HSP90 was targeted by molecular docking of a library of 23,129,083 compounds (from the ZINC database) to the ATP binding site of the N-domain of HSP90. Structure-based virtual screen (SBVS) was performed using idock software on the istar web platform. Based on idock binding energies, 40 molecules were considered as HSP90 inhibitors. In the next step, the 40 molecules and the compound AT13387 (Onalespib) were docked to the XJX binding site using AutoDock Vina software. By comparing the binding energies of the 40 molecules selected with compound AT13387, 26 molecules were selected. By applying the rule of five, eight molecules were selected as hit compounds. The interactions of these eight compounds with the XJX binding site were obtained and investigated, and two-dimensional interaction maps were provided for the others. Finally, computing the toxicity of these compounds with the ProTox-II webserver shows that three compounds, namely ZINC89453765, ZINC23918431 and ZINC12414793, can be considered as good HSP90 inhibitors. These compounds are inactive for nuclear receptor signalling and stress response pathways including heat shock response, so do not have the limitations of common HSP90 inhibitors. They are also inactive for hepatotoxicity, carcinogenicity, immunotoxicity, mutagenicity and cytotoxicity.

Optimal strategies for virtual screening of induced-fit and flexible target in the 2015 D3R Grand Challenge

Induced fit or protein flexibility can make a given structure less useful for docking and/or scoring. The 2015 Drug Design Data Resource (D3R) Grand Challenge provided a unique opportunity to prospectively test optimal strategies for virtual screening in these type of targets: heat shock protein 90 (HSP90), a protein with multiple ligand-induced binding modes; and mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4), a kinase with a large flexible pocket. Using previously known co-crystal structures, we tested predictions from methods that keep the receptor structure fixed and used (a) multiple receptor/ligand co-crystals as binding templates for minimization or docking ("close"), (b) methods that align or dock to a single receptor ("cross"), and (c) a hybrid approach that chose from multiple bound ligands as initial templates for minimization to a single receptor ("min-cross"). Pose prediction using our "close" models resulted in average ligand RMSDs of 0.32 and 1.6 Å for HSP90 and MAP4K4, respectively, the most accurate models of the community-wide challenge. On the other hand, affinity ranking using our "cross" methods performed well overall despite the fact that a fixed receptor cannot model ligand-induced structural changes,. In addition, "close" methods that leverage the co-crystals of the different binding modes of HSP90 also predicted the best affinity ranking. Our studies suggest that analysis of changes on the receptor structure upon ligand binding can help select an optimal virtual screening strategy.

Heat-shock protein 90 (Hsp90) as anticancer target for drug discovery: an ample computational perspective

There are over 100 different types of cancer, and each is classified based on the type of cell that is initially affected. If left untreated, cancer can result in serious health problems and eventually death. Recently, the paradigm of cancer chemotherapy has evolved to use a combination approach, which involves the use of multiple drugs each of which targets an individual protein. Inhibition of heat-shock protein 90 (Hsp90) is one of the novel key cancer targets. Because of its ability to target several signaling pathways, Hsp90 inhibition emerged as a useful strategy to treat a wide variety of cancers. Molecular modeling approaches and methodologies have become 'close counterparts' to experiments in drug design and discovery workflows. A wide range of molecular modeling approaches have been developed, each of which has different objectives and outcomes. In this review, we provide an up-to-date systematic overview on the different computational models implemented toward the design of Hsp90 inhibitors as anticancer agents. Although this is the main emphasis of this review, different topics such as background and current statistics of cancer, different anticancer targets including Hsp90, and the structure and function of Hsp90 from an experimental perspective, for example, X-ray and NMR, are also addressed in this report. To the best of our knowledge, this review is the first account, which comprehensively outlines various molecular modeling efforts directed toward identification of anticancer drugs targeting Hsp90. We believe that the information, methods, and perspectives highlighted in this report would assist researchers in the discovery of potential anticancer agents.

Heat shock proteins in multiple myeloma

Heat shock proteins are molecular chaperones with a central role in protein folding and cellular protein homeostasis. They also play major roles in the development of cancer and in recent years have emerged as promising therapeutic targets. In this review, we discuss the known molecular mechanisms of various heat shock protein families and their involvement in cancer and in particular, multiple myeloma. In addition, we address the current progress and challenges in pharmacologically targeting these proteins as anti-cancer therapeutic strategies.

Small molecule inhibitors in acute myeloid leukemia: from the bench to the clinic

Many patients with acute myeloid leukemia will eventually develop refractory or relapsed disease. In the absence of standard therapy for this population, there is currently an urgent unmet need for novel therapeutic agents. Targeted therapy with small molecule inhibitors represents a new therapeutic intervention that has been successful for the treatment of multiple tumors (e.g., gastrointestinal stromal tumors, chronic myelogenous leukemia). Hence, there has been great interest in generating selective small molecule inhibitors targeting critical pathways of proliferation and survival in acute myeloid leukemia. This review highlights a selective group of intriguing therapeutic agents and their presumed targets in both preclinical models and in early human clinical trials.

New benzothieno[3,2-d]-1,2,3-triazines with antiproliferative activity: synthesis, spectroscopic studies, and biological activity

New benzothieno[3,2-d]-1,2,3-triazines, together with precursors triazenylbenzo[b]thiophenes, were designed, synthesized and screened as anticancer agents. The structural features of these compounds prompted us to investigate their DNA binding capability through UV-vis absorption titrations, circular dichroism, and viscometry, pointing out the occurrence of groove-binding. The derivative 3-(4-methoxy-phenyl)benzothieno[3,2-d]-1,2,3-triazin-4(3H)-one showed the highest antiproliferative effect against HeLa cells and was also tested in cell cycle perturbation experiments. The obtained results assessed for the first time the anticancer activity of benzothieno[3,2-d]-1,2,3-triazine nucleus, and we related it to its DNA-binding properties.

Multivariate analysis in the identification of biological targets for designed molecular structures: the BIOTA protocol

In this work the new protocol BIOlogical Target Assignation (BIOTA) for the prediction of the biological target from molecular structures is proposed. BIOTA is based on the Principal Components Analysis (PCA) application on a matrix of ligands versus molecular descriptors. The application of BIOTA could allow to hypothesize the mechanism of action of a candidate drug prior to its biological evaluation or to repurpose old drugs. The protocol can be fine-tuned by choosing opportune targets (biological or not) and molecular descriptors, and it can be useful in every fields in with it is possible to collect set of compounds with known properties. The robustness of the protocol depends from different factors: the correctness of biological data, the optimization of the molecular structures and their molecular descriptors calculation, the selection of the biological targets. The application of BIOTA to a new class of Hsp90 inhibitors was able to predict quite correctly their affinity for Hsp90.

Exploring the anticancer potential of pyrazolo[1,2-a]benzo[1,2,3,4]tetrazin-3-one derivatives: the effect on apoptosis induction, cell cycle and proliferation

In order to investigate their anticancer potential, four new pyrazolo[1,2-a]benzo[1,2,3,4]tetrazinone derivatives, designed through the chemometric protocol VLAK, and three of the most active compounds of the previous series have been evaluated on some cellular events including proliferation, apoptosis induction, and cell cycle. The NCI one dose (10 μM) screening revealed that the 8,9-di-methyl derivative showed activity against Leukemia (CCRF-CEM) and Colon cancer cell line (COLO 205), reaching 81% and 45% of growth inhibition (GI), respectively. Replacement of the two methyl groups with two chlorine atoms maintained the activity toward Leukemia cell (CCRF-CEM, GI 77%) and selectively enhanced the activity against COLO 205 attaining a LD50 in the μM range and against SW-620 a GI of 77%. Interestingly, an appreciable growth inhibition of 47% against therapeutically "refractory" Non-Small Cell Lung Cancer (NCI-H522) was observed. Moreover, the apoptosis induction, based on mitochondrial membrane depolarization, was found in the range EC50 3-5 μM on HeLa cell, evidencing a well defined relationship with the related in vitro cell growth inhibitory assays (MTT) performed against other selected tumor cell lines not included in the NCI tumor panel (HeLa, cervix; H292, lung; LAN-5, CNS; CaCo-2, colon; 16HBE, normal human cell lung) and against MCF-7 tumor cell line (breast). Only for the most active compounds, further cell cycle tests on HeLa displayed a cell arrest on S phase. Thus, a promising new class of anticancer candidates, acting as valuable apoptotic inductors, is proposed.

Synthesis and biological activities of a new class of heat shock protein 90 inhibitors, designed by energy-based pharmacophore virtual screening

The design through energy-based pharmacophore virtual screening has led to aminocyanopyridine derivatives as efficacious new inhibitors of Hsp90. The synthesized compounds showed a good affinity for the Hsp90 ATP binding site in the competitive binding assay. Moreover, they showed an excellent antiproliferative activity against a large number of human tumor cell lines. Further biological studies on the derivative with the higher EC50 confirmed its specific influence on the cellular pathways involving Hsp90.

Homology modeling, molecular dynamics, e-pharmacophore mapping and docking study of Chikungunya virus nsP2 protease

To date, no suitable vaccine or specific antiviral drug is available to treat Chikungunya viral (CHIKV) fever. Hence, it is essential to identify drug candidates that could potentially impede CHIKV infection. Here, we present the development of a homology model of nsP2 protein based on the crystal structure of the nsP2 protein of Venezuelan equine encephalitis virus (VEEV). The protein modeled was optimized using molecular dynamics simulation; the junction peptides of a nonstructural protein complex were then docked in order to investigate the possible protein-protein interactions between nsP2 and the proteins cleaved by nsP2. The modeling studies conducted shed light on the binding modes, and the critical interactions with the peptides provide insight into the chemical features needed to inhibit the CHIK virus infection. Energy-optimized pharmacophore mapping was performed using the junction peptides. Based on the results, we propose the pharmacophore features that must be present in an inhibitor of nsP2 protease. The resulting pharmacophore model contained an aromatic ring, a hydrophobic and three hydrogen-bond donor sites. Using these pharmacophore features, we screened a large public library of compounds (Asinex, Maybridge, TOSLab, Binding Database) to find a potential ligand that could inhibit the nsP2 protein. The compounds that yielded a fitness score of more than 1.0 were further subjected to Glide HTVS and Glide XP. Here, we report the best four compounds based on their docking scores; these compounds have IDs of 27943, 21362, ASN 01107557 and ASN 01541696. We propose that these compounds could bind to the active site of nsP2 protease and inhibit this enzyme. Furthermore, the backbone structural scaffolds of these four lead compounds could serve as building blocks when designing drug-like molecules for the treatment of Chikungunya viral fever.

IKK-beta inhibitors: an analysis of drug-receptor interaction by using molecular docking and pharmacophore 3D-QSAR approaches

The IKK kinases family represents a thrilling area of research because of its importance in regulating the activity of NF-kB transcription factors. The discovery of the central role played by IKK-beta in the activation of transcription in response to apoptotic or inflammatory stimuli allowed to considerate its modulation as a promising tool for the treatment of chronic inflammation and cancer. To date, several IKK-beta inhibitors have been discovered and tested. In this work, an analysis of the interactions between different classes of inhibitors and their biological target was performed, through the application of Molecular Docking and Pharmacophore/3D-QSAR approaches to a set of 141 inhibitors included in the Binding Database. In order to overcome the difficulty due to the lack of crystallographic data for IKK-beta, a homology model of this protein has been built and validated. The results allowed to study in depth the structural bases for the interaction of each family of inhibitors and provided clues for further modifications, with the aim of improving the activity and selectivity of designed drugs targeting this enzyme.

The HSP90 binding mode of a radicicol-like E-oxime determined by docking, binding free energy estimations, and NMR 15N chemical shifts

We determine the binding mode of a macrocyclic radicicol-like oxime to yeast HSP90 by combining computer simulations and experimental measurements. We sample the macrocyclic scaffold of the unbound ligand by parallel tempering simulations and dock the most populated conformations to yeast HSP90. Docking poses are then evaluated by the use of binding free energy estimations with the linear interaction energy method. Comparison of QM/MM-calculated NMR chemical shifts with experimental shift data for a selective subset of backbone (15)N provides an additional evaluation criteria. As a final test we check the binding modes against available structure-activity-relationships. We find that the most likely binding mode of the oxime to yeast HSP90 is very similar to the known structure of the radicicol-HSP90 complex.

Principal component analysis on molecular descriptors as an alternative point of view in the search of new Hsp90 inhibitors

Inhibiting a protein that regulates multiple signal transduction pathways in cancer cells is an attractive goal for cancer therapy. Heat shock protein 90 (Hsp90) is one of the most promising molecular targets for such an approach. In fact, Hsp90 is a ubiquitous molecular chaperone protein that is involved in folding, activating and assembling of many key mediators of signal transduction, cellular growth, differentiation, stress-response and apoptothic pathways. With the aim to analyze which molecular descriptors have the higher importance in the binding interactions of these classes, we first performed molecular docking experiments on the 187 Hsp90 inhibitors included in the BindingDB, a public database of measured binding affinities. Further, for each frozen conformation obtained from the docking, a set of 250 molecular descriptors was calculated, and the resulting Structure/Descriptors matrix was submitted to Principal Component Analysis. From the factor scores it emerged a good clusterization among similar compounds both in terms of structural class and activity spectrum, while examination of the loadings of the first two factors also allowed to study the classes of descriptors which mainly contribute to each one.