Ligand Binding Thermodynamics in Drug Discovery: Still a Hot Tip?

Solvation Structure and Thermodynamic Mapping (SSTMap): An Open-Source, Flexible Package for the Analysis of Water in Molecular Dynamics Trajectories

We have developed SSTMap, a software package for mapping structural and thermodynamic water properties in molecular dynamics trajectories. The package introduces automated analysis and mapping of local measures of frustration and enhancement of water structure. The thermodynamic calculations are based on Inhomogeneous Fluid Solvation Theory (IST), which is implemented using both site-based and grid-based approaches. The package also extends the applicability of solvation analysis calculations to multiple molecular dynamics (MD) simulation programs by using existing cross-platform tools for parsing MD parameter and trajectory files. SSTMap is implemented in Python and contains both command-line tools and a Python module to facilitate flexibility in setting up calculations and for automated generation of large data sets involving analysis of multiple solutes. Output is generated in formats compatible with popular Python data science packages. This tool will be used by the molecular modeling community for computational analysis of water in problems of biophysical interest such as ligand binding and protein function.

Exploiting a water network to achieve enthalpy-driven, bromodomain-selective BET inhibitors

Within the last decade, the Bromodomain and Extra-Terminal domain family (BET) of proteins have emerged as promising drug targets in diverse clinical indications including oncology, auto-immune disease, heart failure, and male contraception. The BET family consists of four isoforms (BRD2, BRD3, BRD4, and BRDT/BRDT6) which are distinguished by the presence of two tandem bromodomains (BD1 and BD2) that independently recognize acetylated-lysine (KAc) residues and appear to have distinct biological roles. BET BD1 and BD2 bromodomains differ at five positions near the substrate binding pocket: the variation in the ZA channel induces different water networks nearby. We designed a set of congeneric 2- and 3-heteroaryl substituted tetrahydroquinolines (THQ) to differentially engage bound waters in the ZA channel with the goal of achieving bromodomain selectivity. SJ830599 (9) showed modest, but consistent, selectivity for BRD2-BD2. Using isothermal titration calorimetry, we showed that the binding of all THQ analogs in our study to either of the two bromodomains was enthalpy driven. Remarkably, the binding of 9 to BRD2-BD2 was marked by negative entropy and was entirely driven by enthalpy, consistent with significant restriction of conformational flexibility and/or engagement with bound waters. Co-crystallography studies confirmed that 9 did indeed stabilize a water-mediated hydrogen bond network. Finally, we report that 9 retained cytotoxicity against several pediatric cancer cell lines with EC50 values comparable to BET inhibitor (BETi) clinical candidates.

Enthalpy-Based Screening of Focused Combinatorial Libraries for the Identification of Potent and Selective Ligands

In modern drug discovery, the ability of biophysical methods, including nuclear magnetic resonance spectroscopy or surface plasmon resonance, to detect and characterize ligand-protein interactions accurately and unambiguously makes these approaches preferred versus conventional biochemical high-throughput screening of large collections of compounds. Nonetheless, ligand screening strategies that address simultaneously potency and selectivity have not yet been fully developed. In this work, we propose a novel method for screening large collections of combinatorial libraries using enthalpy measurements as a primary screening technique. We demonstrate that selecting binders that are driven by enthalpy (ΔH) results in agents that are not only potent but also more selective for a given target. This general and novel approach, we termed ΔH screening of fPOS (enthalpy screening of focused positional scanning library), combines the principles of focused combinatorial chemistry with rapid calorimetry measurements to efficiently identify potent and selective inhibitors.

An update on anticancer drug development and delivery targeting carbonic anhydrase IX

The expression of carbonic anhydrase (CA) IX is up-regulated in many types of solid tumors in humans under hypoxic and acidic microenvironment. Inhibition of CA IX enzymatic activity with selective inhibitors, antibodies or labeled probes has been shown to reverse the acidic environment of solid tumors and reduce the tumor growth establishing the significant role of CA IX in tumorigenesis. Thus, the development of potent antitumor drugs targeting CA IX with minimal toxic effects is important for the target-specific tumor therapy. Recently, several promising antitumor agents against CA IX have been developed to treat certain types of cancers in combination with radiation and chemotherapy. Here we review the inhibition of CA IX by small molecule compounds and monoclonal antibodies. The methods of enzymatic assays, biophysical methods, animal models including zebrafish and Xenopus oocytes, and techniques of diagnostic imaging to detect hypoxic tumors using CA IX-targeted conjugates are discussed with the aim to overview the recent progress related to novel therapeutic agents that target CA IX in hypoxic tumors.

X-ray Characterization and Structure-Based Optimization of Striatal-Enriched Protein Tyrosine Phosphatase Inhibitors

Excessive activity of striatal-enriched protein tyrosine phosphatase (STEP) in the brain has been detected in numerous neuropsychiatric disorders including Alzheimer's disease. Notably, knockdown of STEP in an Alzheimer mouse model effected an increase in the phosphorylation levels of downstream STEP substrates and a significant reversal in the observed cognitive and memory deficits. These data point to the promising potential of STEP as a target for drug discovery in Alzheimer's treatment. We previously reported a substrate-based approach to the development of low molecular weight STEP inhibitors with K_i values as low as 7.8 μM. Herein, we disclose the first X-ray crystal structures of inhibitors bound to STEP and the surprising finding that they occupy noncoincident binding sites. Moreover, we utilize this structural information to optimize the inhibitor structure to achieve a K_i of 110 nM, with 15-60-fold selectivity across a series of phosphatases.

Ligand efficiency metrics in drug discovery: the pros and cons from a practical perspective

INTRODUCTION

Ligand efficiency metrics are almost universally accepted as a valuable indicator of compound quality and an aid to reduce attrition. Areas covered: In this review, the authors describe ligand efficiency metrics giving a balanced overview on their merits and points of weakness in order to enable the readers to gain an informed opinion. Relevant theoretical breakthroughs and drug-like properties are also illustrated. Several recent exemplary case studies are discussed in order to illustrate the main fields of application of ligand efficiency metrics. Expert opinion: As a medicinal chemist guide, ligand efficiency metrics perform in a context- and chemotype-dependent manner; thus, they should not be used as a magic box. Since the 'big bang' of efficiency metrics occurred more or less ten years ago and the average time to develop a new drug is over the same period, the next few years will give a clearer outlook on the increased rate of success, if any, gained by means of these new intriguing tools.

Determination of the endothelin-1 recognition sites of endothelin receptor type A by the directed-degeneration method

G-protein coupled receptors (GPCRs) play indispensable physiological roles in cell proliferation, differentiation, and migration; therefore, identifying the mechanisms by which ligands bind to GPCRs is crucial for developing GPCR-targeting pharmaceutics and for understanding critical biological functions. Although some structural information is available regarding the interactions between GPCRs and their small molecule ligands, knowledge of how GPCRs interact with their corresponding macromolecule ligands, such as peptides and proteins, remains elusive. In this study, we have developed a novel strategy to investigate the precise ligand recognition mechanisms involved in the interaction of endothelin receptor type A (ET_A) with its ligand, endothelin-1 (ET-1); we call this method “directed degeneration” method. Through flow cytometric screening of a randomized ET_A library, statistical analysis of the identified sequences, and biochemical studies, the ligand interaction map was successfully obtained.

Hot-Spotting with Thermal Scanning: A Ligand- and Structure-Independent Assessment of Target Ligandability

Evaluating the ligandability of a protein target is a key component when defining hit-finding strategies or when prioritize among drug targets. Computational as well as biophysical approaches based on nuclear magnetic resonance (NMR) fragment screening are powerful approaches but suffer from specific constraints that limit their usage. Here, we demonstrate the applicability of high-throughput thermal scanning (HTTS) as a simple and generic biophysical fragment screening method to reproduce assessments from NMR-based screening. By applying this method to a large set of proteins we can furthermore show that the assessment is predictive of the success of high-throughput screening (HTS). The few divergences for targets of low ligandability originate from the sensitivity differences of the orthogonal biophysical methods. We thus applied a new strategy making use of modulations in the solvent structure to improve assay sensitivity. This novel approach enables improved ligandability assessments in accordance with NMR-based assessments and more importantly positions the methodology as a valuable option for biophysical fragment screening.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Structure-based virtual screening and optimization of modulators targeting Hsp90-Cdc37 interaction

Identification of novel Hsp90 inhibitors to disrupt Hsp90-Cdc37 protein-protein interaction (PPI) could be an alternative strategy to achieve Hsp90 inhibition. In this paper, a series of small molecules targeting Hsp90-Cdc37 complex are addressed and characterized. The molecules' key characters are determined by utilizing a structure-based virtual screening workflow, derivatives synthesis, and biological evaluation. Structural optimization and structure-activity relationship (SAR) analysis were then carried out on the virtual hit of VS-8 with potent activity, which resulted in the discovery of compound 10 as a more potent regulator of Hsp90-Cdc37 interaction with a promising inhibitory effect (IC₅₀ = 27 μM), a moderate binding capacity (K_D = 40 μM) and a preferable antiproliferative activity against several cancer lines including MCF-7, SKBR3 and A549 cell lines (IC₅₀ = 26 μM, 15 μM and 38 μM respectively). All the data suggest that compound 10 exhibits moderate inhibitory effect on Hsp90-Cdc37 and could be regard as a first evidence of a non-natural compound targeting Hsp90-Cdc37 PPI.

Estimation of kinetic and thermodynamic ligand-binding parameters using computational strategies

Kinetic and thermodynamic ligand–protein binding parameters are gaining growing importance as key information to consider in drug discovery. The determination of the molecular structures, using particularly x-ray and NMR techniques, is crucial for understanding how a ligand recognizes its target in the final binding complex. However, for a better understanding of the recognition processes, experimental studies of ligand–protein interactions are needed. Even though several techniques can be used to investigate both thermodynamic and kinetic profiles for a ligand–protein complex, these procedures are very often laborious, time consuming and expensive. In the last 10 years, computational approaches have enormous potential in providing insights into each of the above effects and in parsing their contributions to the changes in both kinetic and thermodynamic binding parameters. The main purpose of this review is to summarize the state of the art of computational strategies for estimating the kinetic and thermodynamic parameters of a ligand–protein binding.

Harnessing the Versatility of Optical Biosensors for Target-Based Small-Molecule Drug Discovery

Optical biosensors entered target-based small-molecule drug discovery more than two decades ago and have since transformed into a value-adding component in the decision-making process. Here, we briefly highlight the major application areas of optical biosensors and focus on desirable profiles of such platforms in order to ensure their effective use in small molecule drug discovery. Furthermore, we will emphasize current technology-based constraints and discuss experimental strategies to address these limitations as well as provide a view of necessary technology improvements for next generation platforms.

Intrinsic Thermodynamics and Structures of 2,4- and 3,4-Substituted Fluorinated Benzenesulfonamides Binding to Carbonic Anhydrases

The goal of rational drug design is to understand structure-thermodynamics correlations in order to predict the chemical structure of a drug that would exhibit excellent affinity and selectivity for a target protein. In this study we explored the contribution of added functionalities of benzenesulfonamide inhibitors to the intrinsic binding affinity, enthalpy, and entropy for recombinant human carbonic anhydrases (CA) CA I, CA II, CA VII, CA IX, CA XII, and CA XIII. The binding enthalpies of compounds possessing similar chemical structures and affinities were found to be very different, spanning a range from -90 to +10 kJ mol-1 , and are compensated by a similar opposing entropy contribution. The intrinsic parameters of binding were determined by subtracting the linked protonation reactions. The sulfonamide group pKa values of the compounds were measured spectrophotometrically, and the protonation enthalpies were measured by isothermal titration calorimetry (ITC). Herein we describe the development of meta- or ortho-substituted fluorinated benzenesulfonamides toward the highly potent compound 10 h, which exhibits an observed dissociation constant value of 43 pm and an intrinsic dissociation constant value of 1.1 pm toward CA IX, an anticancer target that is highly overexpressed in various tumors. Fluorescence thermal shift assays, ITC, and X-ray crystallography were all applied in this work.

Interaction of human serum albumin with uremic toxins: a thermodynamic study

Interaction of uremic toxins with HSA is studied by ITC and understood in terms of thermodynamic driving forces.

Binding thermodynamics discriminates fragments from druglike compounds: a thermodynamic description of fragment-based drug discovery

Small is beautiful - reducing the size and complexity of chemical starting points for drug design allows better sampling of chemical space, reveals the most energetically important interactions within protein-binding sites and can lead to improvements in the physicochemical properties of the final drug. The impact of fragment-based drug discovery (FBDD) on recent drug discovery projects and our improved knowledge of the structural and thermodynamic details of ligand binding has prompted us to explore the relationships between ligand-binding thermodynamics and FBDD. Information on binding thermodynamics can give insights into the contributions to protein-ligand interactions and could therefore be used to prioritise compounds with a high degree of specificity in forming key interactions.

Optimization and bioevaluation of Cdc37-derived peptides: An insight into Hsp90-Cdc37 protein-protein interaction modulators

Targeting Hsp90-Cdc37 protein-protein interaction (PPI) is becoming an alternative approach for future anti-cancer drug development. We previously reported the discovery of an eleven-residue peptide (Pep-1) with micromolar activity for the disruption of Hsp90-Cdc37 PPI. Efforts to improve upon the Pep-1 led to the discovery of more potent modulators for Hsp90-Cdc37 PPI. Through the analysis of peptides binding patterns, more peptides were designed for further verification which resulted in Pep-5, the shortest peptide targeting Hsp90-Cdc37, exerting the optimal structure and the most efficient binding mode. Subsequent MD simulation analysis also confirmed that Pep-5 could perform more stable binding ability and better ligand properties than Pep-1. Under the premise of retentive binding capacity, Pep-5 exhibited lower molecular weight and higher ligand efficiency with a K_d value of 5.99μM (Pep-1 K_d=6.90μM) in both direct binding determination and biological evaluation. The optimal and shortest Pep-5 might provide a breakthrough and a better model for the future design of small molecule inhibitors targeting Hsp90-Cdc37 PPI.

Biophysics in drug discovery: impact, challenges and opportunities

Over the past 25 years, biophysical technologies such as X-ray crystallography, nuclear magnetic resonance spectroscopy, surface plasmon resonance spectroscopy and isothermal titration calorimetry have become key components of drug discovery platforms in many pharmaceutical companies and academic laboratories. There have been great improvements in the speed, sensitivity and range of possible measurements, providing high-resolution mechanistic, kinetic, thermodynamic and structural information on compound-target interactions. This Review provides a framework to understand this evolution by describing the key biophysical methods, the information they can provide and the ways in which they can be applied at different stages of the drug discovery process. We also discuss the challenges for current technologies and future opportunities to use biophysical methods to solve drug discovery problems.

Allosteric modulation of AURKA kinase activity by a small-molecule inhibitor of its protein-protein interaction with TPX2

The essential mitotic kinase Aurora A (AURKA) is controlled during cell cycle progression via two distinct mechanisms. Following activation loop autophosphorylation early in mitosis when it localizes to centrosomes, AURKA is allosterically activated on the mitotic spindle via binding to the microtubule-associated protein, TPX2. Here, we report the discovery of AurkinA, a novel chemical inhibitor of the AURKA-TPX2 interaction, which acts via an unexpected structural mechanism to inhibit AURKA activity and mitotic localization. In crystal structures, AurkinA binds to a hydrophobic pocket (the 'Y pocket') that normally accommodates a conserved Tyr-Ser-Tyr motif from TPX2, blocking the AURKA-TPX2 interaction. AurkinA binding to the Y- pocket induces structural changes in AURKA that inhibit catalytic activity in vitro and in cells, without affecting ATP binding to the active site, defining a novel mechanism of allosteric inhibition. Consistent with this mechanism, cells exposed to AurkinA mislocalise AURKA from mitotic spindle microtubules. Thus, our findings provide fresh insight into the catalytic mechanism of AURKA, and identify a key structural feature as the target for a new class of dual-mode AURKA inhibitors, with implications for the chemical biology and selective therapeutic targeting of structurally related kinases.

The Role of Conformational Changes in Molecular Recognition

Conformational changes of molecules are crucial elements in many biochemical processes, and also in molecular recognition. Here, we present a novel exact mathematical equation for the binding free energy of a receptor-ligand pair. It shows that the energetic contribution due to conformational changes upon molecular recognition is defined by the so-called Kullback-Leibler (KL) divergence between the probability distributions of the conformational ensemble of the biomolecule in the bound and free states. We show that conformational changes always contribute positively to the change in free energy and therefore disfavor the association process. Using the example of ligands binding to a flexible cavity of T4 lysozyme, we illustrate that, due to enthalpy-entropy compensation, the conformational entropy is a misleading quantity for assessing the conformational contribution to the binding free energy, in contrast to the KL divergence, which is the correct quantity to use in this context.

Integrating biophysics with HTS-driven drug discovery projects

Over the past decade biophysics has become an established discipline in HTS hit triaging, owing to its high fidelity in detecting protein-ligand interactions. Many pharma companies are using biophysical techniques to filter HTS output for false positives, as will be discussed in this review. Moreover, I will demonstrate how the earlier application of biophysics, already at the HTS assay development stage, is potentially even more impactful. Two key areas here are early mode-of-action studies and ensuring that the HTS assay and subsequent cascade are fit for purpose. Top-level results from 20 in-house projects are shown to underpin the impact of these studies.

Discovery and identification of Cdc37-derived peptides targeting the Hsp90–Cdc37 protein–protein interaction

In order to explore the key residues of the Hsp90–Cdc37 binding interface for further design of peptide inhibitors, a combined strategy of molecular dynamics simulation and MM-PBSA analysis was performed.