Adding calorimetric data to decision making in lead discovery: a hot tip

Surprising Non-Additivity of Methyl Groups in Drug-Kinase Interaction

Drug optimization is guided by biophysical methods with increasing popularity. In the context of lead structure modifications, the introduction of methyl groups is a simple but potentially powerful approach. Hence, it is crucial to systematically investigate the influence of ligand methylation on biophysical characteristics such as thermodynamics. Here, we investigate the influence of ligand methylation in different positions and combinations on the drug-kinase interaction. Binding modes and complex structures were analyzed using protein crystallography. Thermodynamic signatures were measured via isothermal titration calorimetry (ITC). An extensive computational analysis supported the understanding of the underlying mechanisms. We found that not only position but also stereochemistry of the methyl group has an influence on binding potency as well as the thermodynamic signature of ligand binding to the protein. Strikingly, the combination of single methyl groups does not lead to additive effects. In our case, the merger of two methyl groups in one ligand leads to an entirely new alternative ligand binding mode in the protein ligand complex. Moreover, the combination of the two methyl groups also resulted in a nonadditive thermodynamic profile of ligand binding. Molecular dynamics (MD) simulations revealed distinguished characteristic motions of the ligands in solution explaining the pronounced thermodynamic changes. The unexpected drastic change in protein ligand interaction highlights the importance of crystallographic control even for minor modifications such as the introduction of a methyl group. For an in-depth understanding of ligand binding behavior, MD simulations have shown to be a powerful tool.

Flexible Fragment Growing Boosts Potency of Quorum-Sensing Inhibitors against Pseudomonas aeruginosa Virulence

Hit-to-lead optimization is a critical phase in drug discovery. Herein, we report on the fragment-based discovery and optimization of 2-aminopyridine derivatives as a novel lead-like structure for the treatment of the dangerous opportunistic pathogen Pseudomonas aeruginosa. We pursue an innovative treatment strategy by interfering with the Pseudomonas quinolone signal (PQS) quorum sensing (QS) system leading to an abolishment of bacterial pathogenicity. Our compounds act on the PQS receptor (PqsR), a key transcription factor controlling the expression of various pathogenicity determinants. In this target-driven approach, we made use of biophysical screening via surface plasmon resonance (SPR) followed by isothermal titration calorimetry (ITC)-enabled enthalpic efficiency (EE) evaluation. Hit optimization then involved growth vector identification and exploitation. Astonishingly, the latter was successfully achieved by introducing flexible linkers rather than rigid motifs leading to a boost in activity on the target receptor and anti-virulence potency.

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.

Binding site plasticity in viral PPxY Late domain recognition by the third WW domain of human NEDD4

The recognition of PPxY viral Late domains by the third WW domain of the HECT-E3 ubiquitin ligase NEDD4 (hNEDD4-WW3) is essential for the completion of the budding process of numerous enveloped viruses, including Ebola, Marburg, HTLV1 or Rabies. hNEDD4-WW3 has been validated as a promising target for the development of novel host-oriented broad spectrum antivirals. Nonetheless, finding inhibitors with good properties as therapeutic agents remains a challenge since the key determinants of binding affinity and specificity are still poorly understood. We present here a detailed structural and thermodynamic study of the interactions of hNEDD4-WW3 with viral Late domains combining isothermal titration calorimetry, NMR structural determination and molecular dynamics simulations. Structural and energetic differences in Late domain recognition reveal a highly plastic hNEDD4-WW3 binding site that can accommodate PPxY-containing ligands with varying orientations. These orientations are mostly determined by specific conformations adopted by residues I859 and T866. Our results suggest a conformational selection mechanism, extensive to other WW domains, and highlight the functional relevance of hNEDD4-WW3 domain conformational flexibility at the binding interface, which emerges as a key element to consider in the search for potent and selective inhibitors of therapeutic interest.

Machine learning models for hydrogen bond donor and acceptor strengths using large and diverse training data generated by first-principles interaction free energies

We present machine learning (ML) models for hydrogen bond acceptor (HBA) and hydrogen bond donor (HBD) strengths. Quantum chemical (QC) free energies in solution for 1:1 hydrogen-bonded complex formation to the reference molecules 4-fluorophenol and acetone serve as our target values. Our acceptor and donor databases are the largest on record with 4426 and 1036 data points, respectively. After scanning over radial atomic descriptors and ML methods, our final trained HBA and HBD ML models achieve RMSEs of 3.8 kJ mol-1 (acceptors), and 2.3 kJ mol-1 (donors) on experimental test sets, respectively. This performance is comparable with previous models that are trained on experimental hydrogen bonding free energies, indicating that molecular QC data can serve as substitute for experiment. The potential ramifications thereof could lead to a full replacement of wetlab chemistry for HBA/HBD strength determination by QC. As a possible chemical application of our ML models, we highlight our predicted HBA and HBD strengths as possible descriptors in two case studies on trends in intramolecular hydrogen bonding.

Determination of the binding mechanism of histone deacetylase inhibitors

This article places its focus on methods and tools enabling the elucidation of the mechanism by which ligands, small-molecule inhibitors, or substrates interact with zinc-containing bacterial or human members of the histone deacetylase family (HDACs). These methods include biochemical and biophysical approaches and can be subdivided into equilibrium and kinetic methods. More information about the exact mode of action can be obtained by combining these methods with specific mutant variants of the enzymes and/or series of structural similar ligands. All available equilibrium and kinetic data including additional information from 3D structures of HDAC-ligand complexes can be beneficially combined in a data analysis procedure called Integrated Global-Fit analysis eventually providing the most likely binding mechanism.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Thermodynamic analysis of interactions of the Hsp90 with adenosine nucleotides: A comparative perspective

Hsp90s are key proteins in cellular homeostasis since they interact with many client proteins. Several studies indicated that Hsp90s are potential targets for treating diseases, such as cancer or malaria. It has been shown that Hsp90s from different organisms have peculiarities despite their high sequence identity. Therefore, a detailed comparative analysis of several Hsp90 proteins is relevant to the overall understanding of their activity. Accordingly, the goal of this work was to evaluate the interaction of either ADP or ATP with recombinant Hsp90s from different organisms (human α and β isoforms, Plasmodium falciparum, Leishmania braziliensis, yeast and sugarcane) by isothermal titration calorimetry. The measured thermodynamic signatures of those interactions indicated that despite the high identity among all Hsp90s, they have specific thermodynamic characteristics. Specifically, the interactions with ADP are driven by enthalpy but are opposed by entropy, whereas the interaction with ATP is driven by both enthalpy and entropy. Complimentary structural and molecular dynamics studies suggested that specific interactions with ADP that differ from those with ATP may contribute to the observed enthalpies and entropies. Altogether, the data suggest that selective inhibition may be more easily achieved using analogues of the Hsp90-ADP bound state than those of Hsp90-ATP bound state.

2-Aminothiazole Derivatives as Selective Allosteric Modulators of the Protein Kinase CK2. 2. Structure-Based Optimization and Investigation of Effects Specific to the Allosteric Mode of Action

Protein CK2 has gained much interest as an anticancer drug target in the past decade. We had previously described the identification of a new allosteric site on the catalytic α-subunit, along with first small molecule ligands based on the 4-(4-phenylthiazol-2-ylamino)benzoic acid scaffold. In the present work, structure optimizations guided by a binding model led to the identification of the lead compound 2-hydroxy-4-((4-(naphthalen-2-yl)thiazol-2-yl)amino)benzoic acid (27), showing a submicromolar potency against purified CK2α (IC₅₀ = 0.6 μM). Furthermore, 27 induced apoptosis and cell death in 786-O renal cell carcinoma cells (EC₅₀ = 5 μM) and inhibited STAT3 activation even more potently than the ATP-competitive drug candidate CX-4945 (EC₅₀ of 1.6 μM vs 5.3 μM). Notably, the potencies of our allosteric ligands to inhibit CK2 varied depending on the individual substrate. Altogether, the novel allosteric pocket was proved a druggable site, offering an excellent perspective to develop efficient and selective allosteric CK2 inhibitors.

Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option?

Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.

Structures, Activities and Drug-Likeness of Anti-Infective Xanthone Derivatives Isolated from the Marine Environment: A Review

Marine organisms represent almost half of total biodiversity and are a very important source of new bioactive substances. Within the varied biological activities found in marine products, their antimicrobial activity is one of the most relevant. Infectious diseases are responsible for high levels of morbidity and mortality and many antimicrobials lose their effectiveness with time due to the development of resistance. These facts justify the high importance of finding new, effective and safe anti-infective agents. Among the variety of biological activities of marine xanthone derivatives, one that must be highlighted is their anti-infective properties. In this work, a literature review of marine xanthones with anti-infective activity, namely antibacterial, antifungal, antiparasitic and antiviral, is presented. Their structures, biological activity, sources and the methods used for bioactivity evaluation are described. The xanthone derivatives are grouped in three sets: xanthones, hydroxanthones and glycosylated derivatives. Moreover, molecular descriptors, biophysico-chemical properties, and pharmacokinetic parameters were calculated, and the chemical space occupied by marine xanthone derivatives is recognized. The chemical space was compared with marketed drugs and framed accordingly to the drug-likeness concept in order to profile the pharmacokinetic of anti-infective marine xanthone derivatives.

Intrinsic Thermodynamics of Protein-Ligand Binding by Isothermal Titration Calorimetry as Aid to Drug Design

Isothermal titration calorimetry (ITC) is one of the main techniques to determine specific interactions between molecules dissolved in aqueous solution. This technique is commonly used in drug development programs when low-molecular-weight molecules are sought that bind tightly and specifically to a protein (disease target) molecule. The method allows a complete thermodynamic characterization of an interaction, i.e., ITC enables direct determination of the model-independent observed interaction change in enthalpy (ΔH) and a model-dependent observed interaction affinity (change in Gibbs free energy, ΔG) in a single experiment. The product of temperature and change in entropy (TΔS) can be obtained by the subtraction of ΔG from ΔH, and the change in heat capacity (ΔC_p) can be determined as a slope of the temperature dependence of the binding ΔH.Despite the apparent value of ITC in characterization of interactions, it is often forgotten that many protein-ligand binding reactions are linked to protonation-deprotonation reactions or various conformational changes. In such cases, it is important to determine the linked-reaction contributions and obtain the intrinsic values of the changes in Gibbs energy (affinity), enthalpy, and entropy. These energy values can then be used in various SAR-type structure-thermodynamics and combined with structure-kinetics correlations in drug design, when searching for small molecules that would bind the protein target molecule. This manuscript provides a detailed protocol on how to determine the intrinsic values of protein-ligand binding thermodynamics by ITC.

Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design

A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.

Structural basis, stoichiometry, and thermodynamics of binding of the chemokines KC and MIP2 to the glycosaminoglycan heparin

Keratinocyte-derived chemokine (KC or mCXCL1) and macrophage inflammatory protein 2 (MIP2 or mCXCL2) play nonredundant roles in trafficking blood neutrophils to sites of infection and injury. The functional responses of KC and MIP2 are intimately coupled to their interactions with glycosaminoglycans (GAGs). GAG interactions orchestrate chemokine concentration gradients and modulate receptor activity, which together regulate neutrophil trafficking. Here, using NMR, molecular dynamics (MD) simulations, and isothermal titration calorimetry (ITC), we characterized the molecular basis of KC and MIP2 binding to the GAG heparin. Both chemokines reversibly exist as monomers and dimers, and the NMR analysis indicates that the dimer binds heparin with higher affinity. The ITC experiments indicate a stoichiometry of two GAGs per KC or MIP2 dimer and that the enthalpic and entropic contributions vary significantly between the two chemokine-heparin complexes. NMR-based structural models of heparin-KC and heparin-MIP2 complexes reveal that different combinations of residues from the N-loop, 40s turn, β₃-strand, and C-terminal helix form a binding surface within a monomer and that both conserved residues and residues unique to a particular chemokine mediate the binding interactions. MD simulations indicate significant residue-specific differences in their contribution to binding and affinity for a given chemokine and between chemokines. On the basis of our observations that KC and MIP2 bind to GAG via distinct molecular interactions, we propose that the differences in these GAG interactions lead to differences in neutrophil recruitment and play nonoverlapping roles in resolution of inflammation.

Discovery and Optimization of Inhibitors of the Parkinson's Disease Associated Protein DJ-1

DJ-1 is a Parkinson's disease associated protein endowed with enzymatic, redox sensing, regulatory, chaperoning, and neuroprotective activities. Although DJ-1 has been vigorously studied for the past decade and a half, its exact role in the progression of the disease remains uncertain. In addition, little is known about the spatiotemporal regulation of DJ-1, or the biochemical basis explaining its numerous biological functions. Progress has been hampered by the lack of inhibitors with precisely known mechanisms of action. Herein, we have employed biophysical methodologies and X-ray crystallography to identify and to optimize a family of compounds inactivating the critical Cys106 residue of human DJ-1. We demonstrate these compounds are potent inhibitors of various activities of DJ-1 in vitro and in cell-based assays. This study reports a new family of DJ-1 inhibitors with a defined mechanism of action, and contributes toward the understanding of the biological function of DJ-1.

Bayesian analysis of isothermal titration calorimetry for binding thermodynamics

Isothermal titration calorimetry (ITC) is the only technique able to determine both the enthalpy and entropy of noncovalent association in a single experiment. The standard data analysis method based on nonlinear regression, however, provides unrealistically small uncertainty estimates due to its neglect of dominant sources of error. Here, we present a Bayesian framework for sampling from the posterior distribution of all thermodynamic parameters and other quantities of interest from one or more ITC experiments, allowing uncertainties and correlations to be quantitatively assessed. For a series of ITC measurements on metal:chelator and protein:ligand systems, the Bayesian approach yields uncertainties which represent the variability from experiment to experiment more accurately than the standard data analysis. In some datasets, the median enthalpy of binding is shifted by as much as 1.5 kcal/mol. A Python implementation suitable for analysis of data generated by MicroCal instruments (and adaptable to other calorimeters) is freely available online.

Protein-ligand binding enthalpies from near-millisecond simulations: Analysis of a preorganization paradox

Calorimetric studies of protein-ligand binding sometimes yield thermodynamic data that are difficult to understand. Today, molecular simulations can be used to seek insight into such calorimetric puzzles, and, when simulations and experiments diverge, the results can usefully motivate further improvements in computational methods. Here, we apply near-millisecond duration simulations to estimate the relative binding enthalpies of four peptidic ligands with the Grb2 SH2 domain. The ligands fall into matched pairs, where one member of each pair has an added bond that preorganizes the ligand for binding and thus may be expected to favor binding entropically, due to a smaller loss in configurational entropy. Calorimetric studies have shown that the constrained ligands do in fact bind the SH2 domain more tightly than the flexible ones, but, paradoxically, the improvement in affinity for the constrained ligands is enthalpic, rather than entropic. The present enthalpy calculations yield the opposite trend, as they suggest that the flexible ligands bind more exothermically. Additionally, the small relative binding enthalpies are found to be balances of large differences in the energies of structural components such as ligand and the binding site residues. As a consequence, the deviations from experiment in the relative binding enthalpies represent small differences between these large numbers and hence may be particularly susceptible to error, due, for example, to approximations in the force field. We also computed first-order estimates of changes in configurational entropy on binding. These too are, arguably, paradoxical, as they tend to favor binding of the flexible ligands. The paradox is explained in part by the fact that the more rigid constrained ligands reduce the entropy of binding site residues more than their flexible analogs do, at least in the simulations. This result offers a rather general counterargument to the expectation that preorganized ligands should be associated with more favorable binding entropies, other things being equal.

Thermodynamic compatibility between cyclodextrin supramolecular complexes and surfactant

Supramolecular structures based on cyclodextrins have been extensively used for drug delivery systems over decades. However, combining host and guest molecules in a pharmaceutical formulation is not a trivial process, being one of the majors concern the inclusion complex compatibility with other excipients presented in the final formulation. Herein, experimental and theoretical calculations were used to investigate the competition of sodium dodecyl sulfate (SDS) with atenolol (ATE) or losartan (LOS), antihypertensive drugs widely used in the treatment of hypertension. Our findings, using nuclear magnetic resonance and isothermal titrations calorimetry experiments and molecular dynamic simulations demonstrated that LOS remain included into CD cavity after excipient (SDS) addition, which was not verified for ATE ternary system, being the drug displaced by SDS molecule.

In vitro interactions of histones and α-crystallin

The aggregation of crystallins in lenses is associated with cataract formation. We previously reported that mutant crystallins are associated with an increased abundance of histones in knock-in and knockout mouse models. However, very little is known about the specific interactions between lens crystallins and histones. Here, we performed in vitro analyses to determine whether α-crystallin interacts with histones directly. Isothermal titration calorimetry revealed a strong histone–α-crystallin binding with a K_d of 4 × 10⁻⁷ M, and the thermodynamic parameters suggested that the interaction was both entropy and enthalpy driven. Size-exclusion chromatography further showed that histone–α-crystallin complexes are water soluble but become water insoluble as the concentration of histones is increased. Right-angle light scattering measurements of the water-soluble fractions of histone–α-crystallin mixtures showed a decrease in the oligomeric molecular weight of α-crystallin, indicating that histones alter the oligomerization of α-crystallin. Taken together, these findings reveal for the first time that histones interact with and affect the solubility and aggregation of α-crystallin, indicating that the interaction between α-crystallin and histones in the lens is functionally important.

•

Histones are upregulated in α-crystallin mutant mice with hereditary cataracts.

•

Isothermal titration calorimetry reveals histones and α-crystallin interact.

•

These interactions are enthalpy and entropy driven.

•

Histones affect the solubility and aggregation behavior of α-crystallin in vitro.

Catalase and its mysteries

Catalase is one of the firsts in every realm of biological sciences. At the same time it also has a number of unusual features. It has one of the highest turnover numbers of all enzymes. It is essential for neutralizing the noxious hydrogen peroxide both in the nature and the various industries such as dairy, textile and pharmaceutics. It also has the merit of being one of the first protein crystals to be isolated. Ironically its three-dimensional structure was discerned some forty years later. However through the times this senile enzyme has continued to intrigue the scientists by surprising facts and phenomena, such as peculiar interweaving of subunits and remarkable thermal stability. It is also known for suicide inactivation by its own substrate. Catalase is known to be implicated in various medical scenarios and its levels have served as a marker in that capacity. It has even been incorporated into several pharmaceuticals. This review strives to clarify these perspectives. It also draws attention to the biophysical contributions offered by thermodynamics and kinetics in these discoveries. The ultimate aim of this review, however, is to state that the venerable catalase will continue to bewilder us with its mysteries well into the twenty-first century.

Challenges in the Development of a Thiol-Based Broad-Spectrum Inhibitor for Metallo-β-Lactamases

Pathogens, expressing metallo-β-lactamases (MBLs), become resistant against most β-lactam antibiotics. Besides the dragging search for new antibiotics, development of MBL inhibitors would be an alternative weapon against resistant bacterial pathogens. Inhibition of resistance enzymes could restore the antibacterial activity of β-lactams. Various approaches to MBL inhibitors are described; among others, the promising motif of a zinc coordinating thiol moiety is very popular. Nevertheless, since the first report of a thiol-based MBL inhibitor (thiomandelic acid) in 2001, no steps in development of thiol based MBL inhibitors were reported that go beyond clinical isolate testing. In this study, we report on the synthesis and biochemical characterization of thiol-based MBL inhibitors and highlight the challenges behind the development of thiol-based compounds, which exhibit good in vitro activity toward a broad spectrum of MBLs, selectivity against human off-targets, and reasonable activity against clinical isolates.

Crystal structure correlations with the intrinsic thermodynamics of human carbonic anhydrase inhibitor binding

The structure-thermodynamics correlation analysis was performed for a series of fluorine- and chlorine-substituted benzenesulfonamide inhibitors binding to several human carbonic anhydrase (CA) isoforms. The total of 24 crystal structures of 16 inhibitors bound to isoforms CA I, CA II, CA XII, and CA XIII provided the structural information of selective recognition between a compound and CA isoform. The binding thermodynamics of all structures was determined by the analysis of binding-linked protonation events, yielding the intrinsic parameters, i.e., the enthalpy, entropy, and Gibbs energy of binding. Inhibitor binding was compared within structurally similar pairs that differ by para- or meta-substituents enabling to obtain the contributing energies of ligand fragments. The pairs were divided into two groups. First, similar binders—the pairs that keep the same orientation of the benzene ring exhibited classical hydrophobic effect, a less exothermic enthalpy and a more favorable entropy upon addition of the hydrophobic fragments. Second, dissimilar binders—the pairs of binders that demonstrated altered positions of the benzene rings exhibited the non-classical hydrophobic effect, a more favorable enthalpy and variable entropy contribution. A deeper understanding of the energies contributing to the protein-ligand recognition should lead toward the eventual goal of rational drug design where chemical structures of ligands could be designed based on the target protein structure.

NMR-Fragment Based Virtual Screening: A Brief Overview

Fragment-based drug discovery (FBDD) using NMR has become a central approach over the last twenty years for development of small molecule inhibitors against biological macromolecules, to control a variety of cellular processes. Yet, several considerations should be taken into account for obtaining a therapeutically relevant agent. In this review, we aim to list the considerations that make NMR fragment screening a successful process for yielding potent inhibitors. Factors that may govern the competence of NMR in fragment based drug discovery are discussed, as well as later steps that involve optimization of hits obtained by NMR-FBDD.

Discovery of antimicrobial compounds targeting bacterial type FAD synthetases

The increase of bacterial strains resistant to most of the available antibiotics shows a need to explore novel antibacterial targets to discover antimicrobial drugs. Bifunctional bacterial FAD synthetases (FADSs) synthesise the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD). These cofactors act in vital processes as part of flavoproteins, making FADS an essential enzyme. Bacterial FADSs are potential antibacterial targets because of differences to mammalian enzymes, particularly at the FAD producing site. We have optimised an activity-based high throughput screening assay targeting Corynebacterium ammoniagenes FADS (CaFADS) that identifies inhibitors of its different activities. We selected the three best high-performing inhibitors of the FMN:adenylyltransferase activity (FMNAT) and studied their inhibition mechanisms and binding properties. The specificity of the CaFADS hits was evaluated by studying also their effect on the Streptococcus pneumoniae FADS activities, envisaging differences that can be used to discover species-specific antibacterial drugs. The antimicrobial effect of these compounds was also evaluated on C. ammoniagenes, S. pneumoniae, and Mycobacterium tuberculosis cultures, finding hits with favourable antimicrobial properties.

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.

Group Additivity in Ligand Binding Affinity: An Alternative Approach to Ligand Efficiency

Group additivity is a concept that has been successfully applied to a variety of thermochemical and kinetic properties. This includes drug discovery, where functional group additivity is often assumed in ligand binding. Ligand efficiency can be recast as a special case of group additivity where ΔG/HA is the group equivalent (HA is the number of non-hydrogen atoms in a ligand). Analysis of a large data set of protein-ligand binding affinities (K_i) for diverse targets shows that in general ligand binding is distinctly nonlinear. It is possible to create a group equivalent scheme for ligand binding, but only in the context of closely related proteins, at least with regard to size. This finding has broad implications for drug design from both experimental and computational points of view. It also offers a path forward for a more general scheme to assess the efficiency of ligand binding.

Systematic screening of viral entry inhibitors using surface plasmon resonance

Viral binding and entry into host cells for various viruses have been studied extensively, yielding a detailed understanding of the overall viral entry process. As cell entry is an essential and requisite process by which a virus initiates infection, it is an attractive target for therapeutic intervention. The advantages of targeting viral entry are an extracellular target site, relatively easy access for biological interventions, and lower toxicity. Several cell-based strategies and biophysical techniques have been used to screen compounds that block viral entry. These studies led to the discovery of inhibitors against HIV, HCV, influenza, Ebola, and RSV. In recent years, several compounds screened by fragment-based drug discovery have been approved as drugs or are in the final stages of clinical trials. Among fragment screening technologies, surface plasmon resonance has been widely used because it provides accurate information on binding kinetics, allows real-time monitoring of ligand-drug interactions, requires very small sample amounts to perform analyses, and requires no modifications to or labeling of ligands. This review focuses on surface plasmon resonance-based schemes for screening viral entry inhibitors.

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.

Impact of binding mechanism on selective inhibition of histone deacetylase isoforms

Industrialized drug screening campaigns usually deliver hundreds of compounds that are active on a particular pharmaceutical target. In light of high failure rates of drug candidates due to unforeseeable off-target toxicity, the early identification of the most promising compounds with high potential for target selectivity is an urgent need to improve the quality of lead compounds and lower attrition rates in the drug development process. The reliable prediction of the selectivity of active substances for a target protein is a challenging task. A comprehensive study of the binding kinetics, thermodynamics, and selectivity of chemically related ligands of histone deacetylase (HDAC) like amidohydrolase from Pseudomonas aeruginosa (HDAH_pa ) reveals one general binding mechanism for all analyzed compounds consisting of a preceding conformational selection step followed by an optional subsequent induced fit. Depending on the chemical structure, the ligands bind to one or two of at least three protein conformations with different rate constants. Although these kinetic and mechanistic differences hamper the predictability of selectivity for the HDAC inhibitors, we demonstrate that the enthalpy-weighted binding constant KdΔH is a useful metric to predict isoform selectivity of inhibitors against HDAC enzymes and relatively robust toward different but related binding mechanisms.

A Methylidene Group in the Phosphonic Acid Analogue of Phenylalanine Reverses the Enantiopreference of Binding to Phenylalanine Ammonia-Lyases

Aromatic amino acid ammonia‐lyases and aromatic amino acid 2,3‐aminomutases contain the post‐translationally formed prosthetic 3,5‐dihydro‐4‐methylidene‐5H‐imidazol‐5‐one (MIO) group. MIO enzymes catalyze the stereoselective synthesis of α‐ or β‐amino acid enantiomers, making these chemical processes environmentally friendly and affordable. Characterization of novel inhibitors enables structural understanding of enzyme mechanism and recognizes promising herbicide candidates as well. The present study found that both enantiomers of the aminophosphonic acid analogue of the natural substrate phenylalanine and a novel derivative bearing a methylidene at the β‐position inhibited phenylalanine ammonia‐lyases (PAL), representing MIO enzymes. X‐ray methods unambiguously determined the absolute configuration of all tested enantiomers during their synthesis. Enzyme kinetic measurements revealed the enantiomer of the methylidene‐substituted substrate analogue as being a mirror image relation to the natural l‐phenylalanine as the strongest inhibitor. Isothermal titration calorimetry (ITC) confirmed the binding constants and provided a detailed analysis of the thermodynamic driving forces of ligand binding. Molecular docking suggested that binding of the (R)‐ and (S)‐enantiomers is possible by a mirror image packing.

Functional, thermodynamics, structural and biological studies of in silico-identified inhibitors of Mycobacterium tuberculosis enoyl-ACP(CoA) reductase enzyme

Novel chemotherapeutics agents are needed to kill Mycobacterium tuberculosis, the main causative agent of tuberculosis (TB). The M. tuberculosis 2-trans-enoyl-ACP(CoA) reductase enzyme (MtInhA) is the druggable bona fide target of isoniazid. New chemotypes were previously identified by two in silico approaches as potential ligands to MtInhA. The inhibition mode was determined by steady-state kinetics for seven compounds that inhibited MtInhA activity. Dissociation constant values at different temperatures were determined by protein fluorescence spectroscopy. van't Hoff analyses of ligand binding to MtInhA:NADH provided the thermodynamic signatures of non-covalent interactions (ΔH°, ΔS°, ΔG°). Phenotypic screening showed that five compounds inhibited in vitro growth of M. tuberculosis H37Rv strain. Labio_16 and Labio_17 compounds also inhibited the in vitro growth of PE-003 multidrug-resistant strain. Cytotoxic effects on Hacat, Vero and RAW 264.7 cell lines were assessed for the latter two compounds. The Labio_16 was bacteriostatic and Labio_17 bactericidal in an M. tuberculosis-infected macrophage model. In Zebrafish model, Labio_16 showed no cardiotoxicity whereas Labio_17 showed dose-dependent cardiotoxicity. Accordingly, a model was built for the MtInhA:NADH:Labio_16 ternary complex. The results show that the Labio_16 compound is a direct inhibitor of MtInhA, and it may represent a hit for the development of chemotherapeutic agents to treat TB.

The low binding affinity of D-serine at the ionotropic glutamate receptor GluD2 can be attributed to the hinge region

Ionotropic glutamate receptors (iGluRs) are responsible for most of the fast excitatory communication between neurons in our brain. The GluD2 receptor is a puzzling member of the iGluR family: It is involved in synaptic plasticity, plays a role in human diseases, e.g. ataxia, binds glycine and D-serine with low affinity, yet no ligand has been discovered so far that can activate its ion channel. In this study, we show that the hinge region connecting the two subdomains of the GluD2 ligand-binding domain is responsible for the low affinity of D-serine, by analysing GluD2 mutants with electrophysiology, isothermal titration calorimetry and molecular dynamics calculations. The hinge region is highly variable among iGluRs and fine-tunes gating activity, suggesting that in GluD2 this region has evolved to only respond to micromolar concentrations of D-serine.

De Novo Modular Development of a Foldameric Protein-Protein Interaction Inhibitor for Separate Hot Spots: A Dynamic Covalent Assembly Approach

Protein-protein interactions stabilized by multiple separate hot spots are highly challenging targets for synthetic scaffolds. Surface-mimetic foldamers bearing multiple recognition segments are promising candidate inhibitors. In this work, a modular bottom-up approach is implemented by identifying short foldameric recognition segments that interact with the independent hot spots, and connecting them through dynamic covalent library (DCL) optimization. The independent hot spots of a model target (calmodulin) are mapped with hexameric β-peptide helices using a pull-down assay. Recognition segment hits are subjected to a target-templated DCL ligation through thiol-disulfide exchange. The most potent derivative displays low nanomolar affinity towards calmodulin and effectively inhibits the calmodulin-TRPV1 interaction. The DCL assembly of the folded segments offers an efficient approach towards the de novo development of a high-affinity inhibitor of protein-protein interactions.

Opportunities and Challenges for Fatty Acid Mimetics in Drug Discovery

Fatty acids beyond their role as an endogenous energy source and storage are increasingly considered as signaling molecules regulating various physiological effects in metabolism and inflammation. Accordingly, the molecular targets involved in formation and physiological activities of fatty acids hold significant therapeutic potential. A number of these fatty acid targets are addressed by some of the oldest and most widely used drugs such as cyclooxygenase inhibiting NSAIDs, whereas others remain unexploited. Compounds orthosterically binding to proteins that endogenously bind fatty acids are considered as fatty acid mimetics. On the basis of their structural resemblance, fatty acid mimetics constitute a family of bioactive compounds showing specific binding thermodynamics and following similar pharmacokinetic mechanisms. This perspective systematically evaluates targets for fatty acid mimetics, investigates their common structural characteristics, and highlights demands in their discovery and design. In summary, fatty acid mimetics share particularly favorable characteristics justifying the conclusion that their therapeutic potential vastly outweighs the challenges in their design.

Water-Restructuring Mutations Can Reverse the Thermodynamic Signature of Ligand Binding to Human Carbonic Anhydrase

This study uses mutants of human carbonic anhydrase (HCAII) to examine how changes in the organization of water within a binding pocket can alter the thermodynamics of protein-ligand association. Results from calorimetric, crystallographic, and theoretical analyses suggest that most mutations strengthen networks of water-mediated hydrogen bonds and reduce binding affinity by increasing the enthalpic cost and, to a lesser extent, the entropic benefit of rearranging those networks during binding. The organization of water within a binding pocket can thus determine whether the hydrophobic interactions in which it engages are enthalpy-driven or entropy-driven. Our findings highlight a possible asymmetry in protein-ligand association by suggesting that, within the confines of the binding pocket of HCAII, binding events associated with enthalpically favorable rearrangements of water are stronger than those associated with entropically favorable ones.