Fragment-Based Drug Discovery Applied to Hsp90. Discovery of Two Lead Series with High Ligand Efficiency

Design, synthesis and biological evaluation of indole-2-carboxylic acid derivatives as novel HIV-1 integrase strand transfer inhibitors

Integrase plays an important role in the life cycle of HIV-1, and integrase strand transfer inhibitors (INSTIs) can effectively impair the viral replication. However, drug resistance mutations have been confirmed to decrease the efficacy of INSTI during the antiviral therapy. Herein, indole-2-carboxylic acid (1) was found to inhibit the strand transfer of integrase, and the indole nucleus of compound 1 was observed to chelate with two Mg2+ ions within the active site of integrase. Through optimization of compound 1, a series of indole-2-carboxylic acid derivatives were designed and synthesized, and compound 17a was proved to markedly inhibit the effect of integrase, with IC50 value of 3.11 μM. Binding mode analysis of 17a demonstrated that the introduced C6 halogenated benzene ring could effectively bind with the viral DNA (dC20) through π-π stacking interaction. These results indicated that indole-2-carboxylic acid is a promising scaffold for the development of integrase inhibitors.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Combined In-Solution Fragment Screening and Crystallographic Binding-Mode Analysis with a Two-Domain Hsp70 Construct

Heat shock protein 70 (Hsp70) isoforms are key players in the regulation of protein homeostasis and cell death pathways and are therefore attractive targets in cancer research. Developing nucleotide-competitive inhibitors or allosteric modulators, however, has turned out to be very challenging for this protein family, and no Hsp70-directed therapeutics have so far become available. As the field could profit from alternative starting points for inhibitor development, we present the results of a fragment-based screening approach on a two-domain Hsp70 construct using in-solution NMR methods, together with X-ray-crystallographic investigations and mixed-solvent molecular dynamics simulations. The screening protocol resulted in hits on both domains. In particular, fragment binding in a deeply buried pocket at the substrate-binding domain could be detected. The corresponding site is known to be important for communication between the nucleotide-binding and substrate-binding domains of Hsp70 proteins. The main fragment identified at this position also offers an interesting starting point for the development of a dual Hsp70/Hsp90 inhibitor.

Tuning Potential Functions to Host-Guest Binding Data

Software to more rapidly and accurately predict protein-ligand binding affinities is of high interest for early-stage drug discovery, and physics-based methods are among the most widely used technologies for this purpose. The accuracy of these methods depends critically on the accuracy of the potential functions that they use. Potential functions are typically trained against a combination of quantum chemical and experimental data. However, although binding affinities are among the most important quantities to predict, experimental binding affinities have not to date been integrated into the experimental data set used to train potential functions. In recent years, the use of host-guest complexes as simple and tractable models of binding thermodynamics has gained popularity due to their small size and simplicity, relative to protein-ligand systems. Host-guest complexes can also avoid ambiguities that arise in protein-ligand systems such as uncertain protonation states. Thus, experimental host-guest binding data are an appealing additional data type to integrate into the experimental data set used to optimize potential functions. Here, we report the extension of the Open Force Field Evaluator framework to enable the systematic calculation of host-guest binding free energies and their gradients with respect to force field parameters, coupled with the curation of 126 host-guest complexes with available experimental binding free energies. As an initial application of this novel infrastructure, we optimized generalized Born (GB) cavity radii for the OBC2 GB implicit solvent model against experimental data for 36 host-guest systems. This refitting led to a dramatic improvement in accuracy for both the training set and a separate test set with 90 additional host-guest systems. The optimized radii also showed encouraging transferability from host-guest systems to 59 protein-ligand systems. However, the new radii are significantly smaller than the baseline radii and lead to excessively favorable hydration free energies (HFEs). Thus, users of the OBC2 GB model currently may choose between GB cavity radii that yield more accurate binding affinities and GB cavity radii that yield more accurate HFEs. We suspect that achieving good accuracy on both will require more far-reaching adjustments to the GB model. We note that binding free-energy calculations using the OBC2 model in OpenMM gain about a 10× speedup relative to corresponding explicit solvent calculations, suggesting a future role for implicit solvent absolute binding free-energy (ABFE) calculations in virtual compound screening. This study proves the principle of using host-guest systems to train potential functions that are transferrable to protein-ligand systems and provides an infrastructure that enables a range of applications.

Casuarina glauca branchlets' extract as a potential treatment for ulcerative colitis: chemical composition, in silico and in vivo studies

Ulcerative colitis (UC) is an inflammatory bowel disease that is often resistant to current treatment options, leading to a need for alternative therapies. Herbal products have shown promise in managing various conditions, including UC. However, the potential of Casuarina glauca branchlets ethanolic extract (CGBRE) in treating UC has not been explored. This study aimed to analyze the chemical composition of CGBRE and evaluate its efficacy in UC treatment through in silico and in vivo experiments. LC-ESI-MS/MS was used to identify 86 compounds in CGBRE, with 21 potential bioactive compounds determined through pharmacokinetic analysis. Network pharmacology analysis revealed 171 potential UC targets for the bioactive compounds, including EGFR, LRRK2, and HSP90 as top targets, which were found to bind to key CGBRE compounds through molecular docking. Molecular docking findings suggested that CGBRE may be effective in the prevention or treatment of ulcerative colitis mediated by these proteins, where key CGBRE compounds exhibited good binding affinities through formation of numerous interactions. In vivo studies in rats with acetic acid-induced UC demonstrated that oral administration of 300 mg/kg CGBRE for 6 days reduced UC symptoms and colonic expression of EGFR, LRRK2, and HSP90. These findings supported the therapeutic potential of CGBRE in UC and suggested the need for further preclinical and clinical investigation.

Small-molecule dual inhibitors targeting heat shock protein 90 for cancer targeted therapy

Heat shock protein 90, also known as Hsp90, is an extensively preserved molecular chaperone that performs a critical function in organizing various biological pathways and cellular operations. As a potential drug target, Hsp90 is closely linked to cancer. Hsp90 inhibitors are a class of drugs that have been extensively studied in preclinical models and have shown promise in a variety of diseases, especially cancer. However, Hsp90 inhibitors have encountered several challenges in clinical development, such as low efficacy, toxicity, or drug resistance, few Hsp90 small molecule inhibitors have been approved worldwide. Nonetheless, combining Hsp90 inhibitors with other tumor inhibitors, such as HDAC inhibitors, tubulin inhibitors, and Topo II inhibitors, has been shown to have synergistic antitumor effects. Consequently, the development of Hsp90 dual-target inhibitors is an effective strategy in cancer treatment, as it enhances potency while reducing drug resistance. This article provides an overview of Hsp90's domain structure and biological functions, as well as a discussion of the design, discovery, and structure-activity relationships of Hsp90 dual inhibitors, aiming to provide insights into clinical drug research from a medicinal chemistry perspective and discover novel Hsp90 dual inhibitors.

WaterKit: Thermodynamic Profiling of Protein Hydration Sites

Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.

In silico identification of potential Hsp90 inhibitors via ensemble docking, DFT and molecular dynamics simulations

The molecular chaperone heat shock protein 90 (Hsp90) has emerged as one of the most exciting targets for anticancer drug development and Hsp90 inhibitors are potentially useful chemotherapeutic agents in cancer. Within the current study, Hsp90 inhibitors that entered different phases of clinical trials were subjected to Zinc15 structure query to find similar compounds (≥ 78%). Obtained small molecules (1-29) with defined similarity cut-off were docked into ensemble of Hsp90-α NTDs. Docked complexes were ranked on the basis of binding modes and Gibbs free energies as Hsp90 binders (cut-off point; ΔG_b ≤ -12 kcal/mol). Top-ranked compounds were subjected to energy decomposition analysis per residue of binding pocket via density functional theory (DFT) calculations in B3LYP level of theory. Subsequent MD simulations of the top-ranked complexes were performed for 100 ns to explore the stable binding modes during a reasonable period in explicit water. Results of molecular docking and intermolecular binding analysis indicated that H-bond, hydrophobic and salt bridge interactions were determinant forces in complex formation. Compounds 19 and 20 were well accommodated in binding pocket of Hsp90 via relatively varied conformations. It was revealed that Asn51 and Phe138 were key residues that interacted stably to 19 and 20. Although primary mechanism of action for proposed molecules are unknown and yet to be explored, results of the present study revealed key structural features for future structure-guided optimization toward potent inhibitors of Hsp90-α NTD. HighlightsHsp90 inhibitors that entered different phases of clinical trials were subjected to Zinc15 based structure query to afford potential enzyme inhibitors 19 and 20.Quantum chemical calculations confirmed docking results and verified pivotal role of a conserved residues (Asn51, Leu103, Phe138 and Tyr139) in making effective hydrogen bonds.MD simulations of top-ranked docked derivatives revealed the achievement of stable binding modes with less conformational variation of 20 than 19 in the active site of Hsp90-α NTD.H-bond, hydrophobic contacts and salt bridge interactions were determinant forces in binding interactions of in silico hits.Resorcinol and isoxazole were important structural motifs of in silico hits in binding to the active site of Hsp90-α NTD.Communicated by Ramaswamy H. Sarma.

Agathis robusta Bark Extract Protects from Renal Ischemia-Reperfusion Injury: Phytochemical, In Silico and In Vivo Studies

Background: Acute kidney injury (AKI) induced by renal ischemia-reperfusion injury (RIRI) is associated with a high incidence of mortality. Existing therapies are mainly supportive, with no available nephroprotective agent. The purpose of this study is to examine the potential protective effect of Agathis robusta Bark Extract (ARBE) in RIRI. Methods: The chemical composition of ARBE was examined by LC-ESI-MS/MS. Network pharmacology was utilized to identify the RIRI molecular targets that could be aimed at by the identified major components of ARBE. Experimentally validated protein–protein interactions (PPIs) and compound-target networks were constructed using the STRING database and Cytoscape software. Molecular docking studies were employed to assess the interaction of the most relevant ARBE compounds with the hub RIRI-related targets. Furthermore, ARBE was tested in a rat model of RIRI. Results: The phytochemical analysis identified 95 components in ARBE, 37 of which were majors. Network analysis identified 312 molecular targets of RIRI that were associated with ARBE major compounds. Of these 312, the top targets in the experimentally validated PPI network were HSP90, EGFR, and P53. The most relevant compounds based on their peak area and network degree value included narcissoside, isorhamnetin-3-O-glucoside, and syringetin-3-O-glucoside, among others. Docking studies of the most relevant compounds revealed significant interactions with the top RIRI-related targets. In the in vivo RIRI experiments, pretreatment of ARBE improved kidney function and structural changes. ARBE reduced the renal expression of p-NfkB and cleaved caspase-3 by downregulating HSP90 and P53 in rats exposed to RIRI. Conclusion: Taken together, this study revealed the chemical composition of ARBE, depicted the interrelationship of the bioactive ingredients of ARBE with the RIRI-related molecular targets, and validated a nephroprotective effect of ARBE in RIRI.

Evaluating the use of absolute binding free energy in the fragment optimisation process

Key to the fragment optimisation process within drug design is the need to accurately capture the changes in affinity that are associated with a given set of chemical modifications. Due to the weakly binding nature of fragments, this has proven to be a challenging task, despite recent advancements in leveraging experimental and computational methods. In this work, we evaluate the use of Absolute Binding Free Energy (ABFE) calculations in guiding fragment optimisation decisions, retrospectively calculating binding free energies for 59 ligands across 4 fragment elaboration campaigns. We first demonstrate that ABFEs can be used to accurately rank fragment-sized binders with an overall Spearman's r of 0.89 and a Kendall τ of 0.67, although often deviating from experiment in absolute free energy values with an RMSE of 2.75 kcal/mol. We then also show that in several cases, retrospective fragment optimisation decisions can be supported by the ABFE calculations. Comparing against cheaper endpoint methods, namely N_wat-MM/GBSA, we find that ABFEs offer better ranking power and correlation metrics. Our results indicate that ABFE calculations can usefully guide fragment elaborations to maximise affinity.

Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix)

BACKGROUND

Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules.

METHOD

To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions.

RESULT

We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential.

CONCLUSION

Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.

Selective Inhibition of the Hsp90α Isoform

The 90 kDa heat shock protein (Hsp90) is a molecular chaperone that processes nascent polypeptides into their biologically active conformations. Many of these proteins contribute to the progression of cancer, and consequently, inhibition of the Hsp90 protein folding machinery represents an innovative approach toward cancer chemotherapy. However, clinical trials with Hsp90 N-terminal inhibitors have encountered deleterious side effects and toxicities, which appear to result from the pan-inhibition of all four Hsp90 isoforms. Therefore, the development of isoform-selective Hsp90 inhibitors is sought to delineate the pathological role played by each isoform. Herein, we describe a structure-based approach that was used to design the first Hsp90α-selective inhibitors, which exhibit >50-fold selectivity versus other Hsp90 isoforms.

Assay design and development strategies for finding Hsp90 inhibitors and their role in human diseases

Heat shock protein 90 (Hsp90) is a molecular chaperone that facilitates the maturation of its client proteins including protein kinases, transcription factors, and steroid hormone receptors which are structurally and functionally diverse. These client proteins are involved in various cellular signaling pathways, and Hsp90 is implicated in various human diseases including cancer, inflammation, and diseases associated with protein misfolding; thus making Hsp90 a promising target for drug discovery. Some of its client proteins are well-known cancer targets. Instead of targeting these client proteins individually, however, targeting Hsp90 is more practical for cancer drug development. Efforts have been invested in recognizing potential drugs for clinical use that inhibit Hsp90 activity and result in the prevention of Hsp90 client maturation and dampening of subsequent signaling cascades. Here, we discuss current assays and technologies used to find and characterize Hsp90 inhibitors that include biophysical, biochemical, cell-based assays and computational modeling. This review highlights recent discoveries that N-terminal isoform-selective compounds and inhibitors that target the Hsp90 C-terminus that may offer the potential to overcome some of the detriments observed with pan Hsp90 inhibitors. The tools and assays summarized in this review should be used to develop Hsp90-targeting drugs with high specificity, potency, and drug-like properties that may prove immensely useful in the clinic.

Targeting the HSP90-CDC37-kinase chaperone cycle: A promising therapeutic strategy for cancer

Heat shock protein 90 (HSP90) is an indispensable molecular chaperone that facilitates the maturation of numerous oncoproteins in cancer cells, including protein kinases, ribonucleoproteins, steroid hormone receptors, and transcription factors. Although over 30 HSP90 inhibitors have steadily entered clinical trials, further clinical advancement has been restricted by their limited efficacy, inevitable heat shock response, and multiple side-effects, likely induced via an ATP inhibition mechanism. Since both ATP and various co-chaperones play essential roles in the HSP90 chaperone cycle to achieve integrated function, optimal therapeutics require an understanding of the dynamic interactions among HSP90, ATP, and cochaperones. To date, continuous research has promoted the exploration of the cochaperone cell division cycle 37 (CDC37) as a kinase-specific recognizer and has shown that the HSP90-CDC37-kinase complex is particularly relevant in cancers. Indeed, disrupting the HSP90-CDC37-kinase complex, rather than totally blocking the ATP function of HSP90, is emerging as an alternative way to avoid the limitations of current inhibitors. In this review, we first briefly introduce the HSP90-CDC37-kinase cycle and present the currently available approaches for inhibitor development targeting this cycle and provide insights into selective regulation of the kinase clients of HSP90 by more directional ways.

Structure-Based Design of Potent and Orally Active Isoindolinone Inhibitors of MDM2-p53 Protein-Protein Interaction

Inhibition of murine double minute 2 (MDM2)-p53 protein-protein interaction with small molecules has been shown to reactivate p53 and inhibit tumor growth. Here, we describe rational, structure-guided, design of novel isoindolinone-based MDM2 inhibitors. MDM2 X-ray crystallography, quantum mechanics ligand-based design, and metabolite identification all contributed toward the discovery of potent in vitro and in vivo inhibitors of the MDM2-p53 interaction with representative compounds inducing cytostasis in an SJSA-1 osteosarcoma xenograft model following once-daily oral administration.

Fragment-based drug discovery: opportunities for organic synthesis

This Review describes the increasing demand for organic synthesis to facilitate fragment-based drug discovery (FBDD), focusing on polar, unprotected fragments. In FBDD, X-ray crystal structures are used to design target molecules for synthesis with new groups added onto a fragment via specific growth vectors. This requires challenging synthesis which slows down drug discovery, and some fragments are not progressed into optimisation due to synthetic intractability. We have evaluated the output from Astex's fragment screenings for a number of programs, including urokinase-type plasminogen activator, hematopoietic prostaglandin D2 synthase, and hepatitis C virus NS3 protease-helicase, and identified fragments that were not elaborated due, in part, to a lack of commercially available analogues and/or suitable synthetic methodology. This represents an opportunity for the development of new synthetic research to enable rapid access to novel chemical space and fragment optimisation.

Inhibition of the NLRP3 inflammasome by HSP90 inhibitors

Excessive and dysregulated inflammation is known to contribute to disease progression. HSP90 is an intracellular chaperone known to regulate inflammatory processes including the NLRP3 inflammasome and secretion of the pro‐inflammatory cytokine interleukin(IL)‐1β. Here, primarily using an in vitro inflammasome ASC speck assay, and an in vivo model of murine peritonitis, we tested the utility of HSP90 inhibitors as anti‐inflammatory molecules. We report that the HSP90 inhibitor EC144 effectively inhibited inflammatory processes including priming and activation of NLRP3 in vitro and in vivo. A specific inhibitor of the β HSP90 isoform was ineffective suggesting the importance of the α isoform in inflammatory signalling. EC144 inhibited IL‐1β and IL‐6 in vivo when administered orally, and was brain‐penetrant. These data suggest that HSP90 inhibitors may be useful for targeting inflammation in diverse diseases that are worsened by the presence of inflammation.

Old and New Approaches to Target the Hsp90 Chaperone

The 90-kDa heat shock protein (Hsp90) is a molecular chaperone that ensures cellular proteostasis by maintaining the folding, stabilization, activation, and degradation of over 400 client proteins. Hsp90 is not only critical for routine protein maintenance in healthy cells, but also during states of cellular stress, such as cancer and neurodegenerative diseases. Due to its ability to affect phosphorylation of numerous client proteins, inhibition of Hsp90 has been an attractive anticancer approach since the early 1990's, when researchers identified a druggable target on the amino terminus of Hsp90 for a variety of cancers. Since then, 17 Hsp90 inhibitors that target the chaperone's Nterminal domain, have entered clinical trials. None, however, have been approved thus far by the FDA as a cancer monotherapy. In these trials, a major limitation observed with Hsp90 inhibition at the N-terminal domain was dose-limiting toxicities and relatively poor pharmacokinetic profiles. Despite this, preclinical and clinical research continues to show that Hsp90 inhibitors effectively target cancer cell death and decrease tumor progression supporting the rationale for the development of novel Hsp90 inhibitors. Here, we present an in-depth overview of the Hsp90 inhibitors used in clinical trials. Finally, we present current shifts in the field related to targeting the carboxy-terminal domain of Hsp90 as well as to the development of isoform-selective inhibitors as a means to bypass the pitfalls of current Hsp90 inhibitors and improve clinical trial outcomes.

NMR as a "Gold Standard" Method in Drug Design and Discovery

Studying disease models at the molecular level is vital for drug development in order to improve treatment and prevent a wide range of human pathologies. Microbial infections are still a major challenge because pathogens rapidly and continually evolve developing drug resistance. Cancer cells also change genetically, and current therapeutic techniques may be (or may become) ineffective in many cases. The pathology of many neurological diseases remains an enigma, and the exact etiology and underlying mechanisms are still largely unknown. Viral infections spread and develop much more quickly than does the corresponding research needed to prevent and combat these infections; the present and most relevant outbreak of SARS-CoV-2, which originated in Wuhan, China, illustrates the critical and immediate need to improve drug design and development techniques. Modern day drug discovery is a time-consuming, expensive process. Each new drug takes in excess of 10 years to develop and costs on average more than a billion US dollars. This demonstrates the need of a complete redesign or novel strategies. Nuclear Magnetic Resonance (NMR) has played a critical role in drug discovery ever since its introduction several decades ago. In just three decades, NMR has become a "gold standard" platform technology in medical and pharmacology studies. In this review, we present the major applications of NMR spectroscopy in medical drug discovery and development. The basic concepts, theories, and applications of the most commonly used NMR techniques are presented. We also summarize the advantages and limitations of the primary NMR methods in drug development.

Hsp90 chaperones have an energetic hot-spot for binding inhibitors

Although Hsp90-family chaperones have been extensively targeted with ATP-competitive inhibitors, it is unknown whether high affinity is achieved from a few highly stabilizing contacts or from many weaker contacts within the ATP-binding pocket. A large-scale analysis of Hsp90α:inhibitor structures shows that inhibitor hydrogen-bonding to a conserved aspartate (D93 in Hsp90α) stands out as most universal among Hsp90 inhibitors. Here we show that the D93 region makes a dominant energetic contribution to inhibitor binding for both cytosolic and organelle-specific Hsp90 paralogs. For inhibitors in the resorcinol family, the D93:inhibitor hydrogen-bond is pH-dependent because the associated inhibitor hydroxyl group is titratable, rationalizing a linked-protonation event previously observed by the Matulis group. The inhibitor hydroxyl group pKa associated with the D93 hydrogen-bond is therefore critical for optimizing the affinity of resorcinol derivatives, and we demonstrate that spectrophotometric measurements can determine this pKa value. Quantifying the energetic contribution of the D93 hotspot is best achieved with the mitochondrial Hsp90 paralog, yielding 3-6 kcal/mol of stabilization (35-60% of the total binding energy) for a diverse set of inhibitors. The Hsp90 Asp93➔Asn substitution has long been known to abolish nucleotide binding, yet puzzlingly, native sequences of structurally similar ATPases, such as Topoisomerasese II, have an asparagine at this same crucial site. While aspartate and asparagine sidechains can both act as hydrogen bond acceptors, we show that a steric clash prevents the Hsp90 Asp93➔Asn sidechain from adopting the necessary rotamer, whereas this steric restriction is absent in Topoisomerasese II.

Design and Synthesis of Fungal-Selective Resorcylate Aminopyrazole Hsp90 Inhibitors

The molecular chaperone Hsp90, essential in all eukaryotes, plays a multifaceted role in promoting survival, virulence, and drug resistance across diverse pathogenic fungal species. The chaperone is also critically important, however, to the pathogen's human host, preventing the use of known clinical Hsp90 inhibitors in antifungal applications due to concomitant host toxicity issues. With the goal of developing Hsp90 inhibitors with acceptable therapeutic indices for the treatment of invasive fungal infections, we initiated a program to design and synthesize potent inhibitors with selective activity against fungal Hsp90 isoforms over their human counterparts. Building on our previously reported derivatization of resorcylate natural products to produce fungal-selective compounds, we have developed a series of synthetic aminopyrazole-substituted resorcylate amides with broad, potent, and fungal-selective Hsp90 inhibitory activity. Herein we describe the synthesis of this series, as well as biochemical structure-activity relationships driving selectivity for the Hsp90 isoforms expressed by Cryptococcus neoformans and Candida albicans, two pathogenic fungi of major clinical importance.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Thermodynamic profiling for fragment-based lead discovery and optimization

Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.

1D NMR WaterLOGSY as an efficient method for fragment-based lead discovery

WaterLOGSY is a sensitive ligand-observed NMR experiment for detection of interaction between a ligand and a protein and is now well-established as a screening technique for fragment-based lead discovery. Here we develop and assess a protocol to derive ligand epitope mapping from WaterLOGSY data and demonstrate its general applicability in studies of fragment-sized ligands binding to six different proteins (glycogen phosphorylase, protein peroxiredoxin 5, Bcl-x_L, Mcl-1, HSP90, and human serum albumin). We compare the WaterLOGSY results to those obtained from the more widely used saturation transfer difference experiments and to the 3D structures of the complexes when available. In addition, we evaluate the impact of ligand labile protons on the WaterLOGSY data. Our results demonstrate that the WaterLOGSY experiment can be used as an additional confirmation of the binding mode of a ligand to a protein.

MDA-9/Syntenin: An emerging global molecular target regulating cancer invasion and metastasis

With few exceptions, metastasis is the terminal stage of cancer with limited therapeutic options. Metastasis consists of numerous phenotypic and genotypic alterations of cells that are directly and indirectly induced by multiple intrinsic (cellular) and extrinsic (micro-environmental) factors. To metastasize, a cancer cell often transitions from an epithelial to mesenchymal morphology (EMT), modifies the extracellular matrix, forms emboli and survives in the circulation, escapes immune surveillance, adheres to sites distant from the initial tumor and finally develops a blood supply (angiogenesis) and colonizes in a secondary niche (a micrometastasis). Scientific advances have greatly enhanced our understanding of the precise molecular and genetic changes, operating independently or collectively, that lead to metastasis. This review focuses on a unique gene, melanoma differentiation associated gene-9 (also known as Syntenin-1; Syndecan Binding Protein (sdcbp); mda-9/syntenin), initially cloned and characterized from metastatic human melanoma and shown to be a pro-metastatic gene. In the last two decades, our comprehension of the diversity of actions of MDA-9/Syntenin on cellular phenotype has emerged. MDA-9/Sytenin plays pivotal regulatory roles in multiple signaling cascades and orchestrates both metastatic and non-metastatic events. Considering the relevance of this gene in controlling cancer invasion and metastasis, approaches have been developed to uniquely and selectively target this gene. We also provide recent updates on strategies that have been successfully employed in targeting MDA-9/Syntenin resulting in profound pre-clinical anti-cancer activity.

Fragment Hits: What do They Look Like and How do They Bind?

A “fragment hit”, a molecule of low molecular weight that has been validated to bind to a target protein, can be an effective chemical starting point for a drug discovery project. Our ability to find and progress fragment hits could potentially be improved by enhancing our understanding of their binding properties, which to date has largely been based on tacit knowledge and reports from individual projects. In the work reported here, we systematically analyzed the molecular and binding properties of fragment hits using 489 published protein–fragment complexes. We identified a number of notable features that these hits tend to have in common, including preferences in buried surface area upon binding, hydrogen bonding and other directional interactions with the protein targets, structural topology, functional-group occurrence, and degree of carbon saturation. In the future, taking account of these preferences in designing and selecting fragments to screen against protein targets may increase the chances of success in fragment screening campaigns.

The nature of ligand efficiency

Ligand efficiency is a widely used design parameter in drug discovery. It is calculated by scaling affinity by molecular size and has a nontrivial dependency on the concentration unit used to express affinity that stems from the inability of the logarithm function to take dimensioned arguments. Consequently, perception of efficiency varies with the choice of concentration unit and it is argued that the ligand efficiency metric is not physically meaningful nor should it be considered to be a metric. The dependence of ligand efficiency on the concentration unit can be eliminated by defining efficiency in terms of sensitivity of affinity to molecular size and this is illustrated with reference to fragment-to-lead optimizations. Group efficiency and fit quality are also examined in detail from a physicochemical perspective. The importance of examining relationships between affinity and molecular size directly is stressed throughout this study and an alternative to ligand efficiency for normalization of affinity with respect to molecular size is presented.

The succinct synthesis of AT13387, a clinically relevant Hsp90 inhibitor

AT13387 is an orally bioavailable clinical candidate developed to inhibit theheat shock protein 90 (Hsp90). This article describes a modified synthetic route for the multi-gram production of AT13387 in 46% overall yield. The modified synthetic route is short, avoids stringent reaction conditions and difficult purifications, which led to increase in an overall yield.

Chaperome heterogeneity and its implications for cancer study and treatment

The chaperome is the collection of proteins in the cell that carry out molecular chaperoning functions. Changes in the interaction strength between chaperome proteins lead to an assembly that is functionally and structurally distinct from each constituent member. In this review, we discuss the epichaperome, the cellular network that forms when the chaperome components of distinct chaperome machineries come together as stable, functionally integrated, multimeric complexes. In tumors, maintenance of the epichaperome network is vital for tumor survival, rendering them vulnerable to therapeutic interventions that target critical epichaperome network components. We discuss how the epichaperome empowers an approach for precision medicine cancer trials where a new target, biomarker, and relevant drug candidates can be correlated and integrated. We introduce chemical biology methods to investigate the heterogeneity of the chaperome in a given cellular context. Lastly, we discuss how ligand-protein binding kinetics are more appropriate than equilibrium binding parameters to characterize and unravel chaperome targeting in cancer and to gauge the selectivity of ligands for specific tumor-associated chaperome pools.

The role of small-angle scattering in structure-based screening applications

In many biomolecular interactions, changes in the assembly states and structural conformations of participants can act as a complementary reporter of binding to functional and thermodynamic assays. This structural information is captured by a number of structural biology and biophysical techniques that are viable either as primary screens in small-scale applications or as secondary screens to complement higher throughput methods. In particular, small-angle X-ray scattering (SAXS) reports the average distance distribution between all atoms after orientational averaging. Such information is important when for example investigating conformational changes involved in inhibitory and regulatory mechanisms where binding events do not necessarily cause functional changes. Thus, we summarise here the current and prospective capabilities of SAXS-based screening in the context of other methods that yield structural information. Broad guidelines are also provided to assist readers in preparing screening protocols that are tailored to available X-ray sources.

Probing the role of Arg97 in Heat shock protein 90 N-terminal domain from the parasite Leishmania braziliensis through site-directed mutagenesis on the human counterpart

In Brazil, the mucocutaneous form of leishmaniasis, caused by the parasite Leishmania braziliensis, is a widespread and very challenging disease responsible for disfiguration and, in the most severe cases, death. Heat shock protein 90 (Hsp90) is a ubiquitous molecular chaperone playing a pivotal role in the folding process of client proteins, and therefore its activity is fundamental for cell survival and proliferation. Since the chaperone activity requires ATP hydrolysis, molecules able to occupy the ATP binding pocket in the protein N-terminal domain (NTD) act as Hsp90 inhibitors. The development of selective molecules targeting the ATPase site of protozoan Hsp90 is tricky for the high homology with the human Hsp90 NTD (hNTD). Notably, only the human Lys112 is replaced by Arg97 in the L. braziliensis enzyme. Recently, this difference has been probed to design selective inhibitors targeting parasite Hsp90s. Here, a reliable protocol for expression and purification of LbHsp90-NTD (LbNTD) was developed but its structural characterization was unsuccessful. The role of Arg97 in LbNTD was hence probed by means of the "leishmanized" K112R variant of hNTDα. To deeply investigate the role of this residue, also the hNTDα K112A variant was generated. Structural studies performed on hNTDα and its variants using various ADP and ATP analogues and cAMP revealed that this residue is not crucial for nucleotide binding. This finding strongly suggests that Arg97 in LbNTD and more generally the conserved arginine residue in parasite Hsp90s are not exploitable for the development of selective inhibitors.

Methyl group reorientation under ligand binding probed by pseudocontact shifts

Liquid-state NMR spectroscopy is a powerful technique to elucidate binding properties of ligands on proteins. Ligands binding in hydrophobic pockets are often in close proximity to methyl groups and binding can lead to subtle displacements of methyl containing side chains to accommodate the ligand. To establish whether pseudocontact shifts can be used to characterize ligand binding and the effects on methyl groups, the N-terminal domain of HSP90 was tagged with caged lanthanoid NMR probe 5 at three positions and titrated with a ligand. Binding was monitored using the resonances of leucine and valine methyl groups. The pseudocontact shifts (PCS) caused by ytterbium result in enhanced dispersion of the methyl spectrum, allowing more resonances to be observed. The effects of tag attachment on the spectrum and ligand binding are small. Significant changes in PCS were observed upon ligand binding, indicating displacements of several methyl groups. By determining the cross-section of PCS iso-surfaces generated by two or three paramagnetic centers, the new position of a methyl group can be estimated, showing displacements in the range of 1-3 Å for methyl groups in the binding site. The information about such subtle but significant changes may be used to improve docking studies and can find application in fragment-based drug discovery.

Polypharmacology of Berberine Based on Multi-Target Binding Motifs

Background: Polypharmacology is emerging as the next paradigm in drug discovery. However, considerable challenges still exist for polypharmacology modeling. In this study, we developed a rational design to identify highly potential targets (HPTs) for polypharmacological drugs, such as berberine. Methods and Results: All the proven co-crystal structures locate berberine in the active cavities of a redundancy of aromatic, aliphatic, and acidic residues. The side chains from residues provide hydrophobic and electronic interactions to aid in neutralization for the positive charge of berberine. Accordingly, we generated multi-target binding motifs (MBM) for berberine, and established a new mathematical model to identify HPTs based on MBM. Remarkably, the berberine MBM was embodied in 13 HPTs, including beta-secretase 1 (BACE1) and amyloid-β_1-42 (Aβ_1-42). Further study indicated that berberine acted as a high-affinity BACE1 inhibitor and prevented Aβ_1-42 aggregation to delay the pathological process of Alzheimer's disease. Conclusion: Here, we proposed a MBM-based drug-target space model to analyze the underlying mechanism of multi-target drugs against polypharmacological profiles, and demonstrated the role of berberine in Alzheimer's disease. This approach can be useful in derivation of rules, which will illuminate our understanding of drug action in diseases.

Structure-guided design of an Hsp90β N-terminal isoform-selective inhibitor

The 90 kDa heat shock protein (Hsp90) is a molecular chaperone responsible for folding proteins that are directly associated with cancer progression. Consequently, inhibition of the Hsp90 protein folding machinery results in a combinatorial attack on numerous oncogenic pathways. Seventeen small-molecule inhibitors of Hsp90 have entered clinical trials, all of which bind the Hsp90 N-terminus and exhibit pan-inhibitory activity against all four Hsp90 isoforms. pan-Inhibition of Hsp90 appears to be detrimental as toxicities have been reported alongside induction of the pro-survival heat shock response. The development of Hsp90 isoform-selective inhibitors represents an alternative approach towards the treatment of cancer that may limit some of the detriments. Described herein is a structure-based approach to design isoform-selective inhibitors of Hsp90β, which induces the degradation of select Hsp90 clients without concomitant induction of Hsp90 levels. Together, these initial studies support the development of Hsp90β-selective inhibitors as a method to overcome the detriments associated with pan-inhibition.

The molecular chaperone Hsp90 oversees the folding of many proteins associated with cancer progression but existing small-molecule inhibitors of this pathway are not isoform-selective. Here, the authors rationally design an Hsp90 inhibitor that displays high selectivity for the Hsp90β isoform.

Dynamics revelation of conformational changes and binding modes of heat shock protein 90 induced by inhibitor associations

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment. The current work integrates molecular dynamics (MD) simulations, calculations of binding free energy, and principal component (PC) analysis with scanning of inhibitor–residue interaction to probe the binding modes of inhibitors YK9, YKJ and YKI to Hsp90 and identify the hot spot of the inhibitor–Hsp90 binding. The results suggest that the introductions of two groups G1 and G2 into YKJ and YKI strengthen the binding ability of YKJ and YKI to Hsp90 compared to YK9. PC analysis based MD trajectories prove that inhibitor bindings exert significant effects on the conformational changes, internal dynamics and motion modes of Hsp90, especially for the helix α2 and the loops L1 and L2. The calculations of residue-based free energy decomposition and scanning of the inhibitor–Hsp90 interaction suggest that six residues L107, G108, F138, Y139, W162 and F170 construct the common hot spot of the inhibitor–residue interactions. Moreover the substitutions of the groups G1 and G2 in YKJ and YKI lead to two additional hydrogen bonding interactions and multiple hydrophobic interactions for bindings of YKJ and YKI to Hsp90. This work is also expected to contribute theoretical hints for the design of potent inhibitors toward Hsp90.

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment.

Fragment-based drug discovery and its application to challenging drug targets

Fragment-based drug discovery (FBDD) is a technique for identifying low molecular weight chemical starting points for drug discovery. Since its inception 20 years ago, FBDD has grown in popularity to the point where it is now an established technique in industry and academia. The approach involves the biophysical screening of proteins against collections of low molecular weight compounds (fragments). Although fragments bind to proteins with relatively low affinity, they form efficient, high quality binding interactions with the protein architecture as they have to overcome a significant entropy barrier to bind. Of the biophysical methods available for fragment screening, X-ray protein crystallography is one of the most sensitive and least prone to false positives. It also provides detailed structural information of the protein-fragment complex at the atomic level. Fragment-based screening using X-ray crystallography is therefore an efficient method for identifying binding hotspots on proteins, which can then be exploited by chemists and biologists for the discovery of new drugs. The use of FBDD is illustrated here with a recently published case study of a drug discovery programme targeting the challenging protein-protein interaction Kelch-like ECH-associated protein 1:nuclear factor erythroid 2-related factor 2.

Resorcinol-Based Grp94-Selective Inhibitors

Glucose regulated protein 94 (Grp94) is the endoplasmic reticulum resident of the 90 kDa heat shock protein (Hsp90) family and represents a promising therapeutic target for the treatment of several diseases. Grp94 is the most unique member of the 90 kDa heat shock protein family due to a five amino acid insertion into its primary sequence, which creates hydrophobic subpockets exclusive to Grp94 that can be utilized for selective inhibition. The first resorcinol-based Grp94-selective inhibitor to take advantage of the hydrophobic S2 subpocket has been developed and shown to manifest low nanomolar affinity and ∼10-fold selectivity for Grp94. Furthermore, these Grp94-selective inhibitors manifest low micromolar GI₅₀ values against multiple myeloma cells, supporting Grp94 as an emerging target for the treatment of this disease.

POVME 3.0: Software for Mapping Binding Pocket Flexibility

We present a substantial update to the open-source POVME binding pocket analysis software. New capabilities of POVME 3.0 include a flexible chemical coloring scheme for feature identification, postanalysis tools for comparing large ensembles of pockets (e.g., from molecular dynamics simulations), and the introduction of scripts and methods that facilitate binding pocket comparison and analysis. We envision the use of this software for visualization of binding pocket dynamics, selection of representative structures for ensemble docking, and incorporation of molecular dynamics results into ligand design efforts.