Paradoxically, Most Flexible Ligand Binds Most Entropy-Favored: Intriguing Impact of Ligand Flexibility and Solvation on Drug–Kinase Binding

Calculations of Absolute Free Energies, Enthalpies, and Entropies for Drug Binding

Computer simulation methods can aid in the rational design of drugs aimed at a specific target, typically a protein. The affinity of a drug for its target is given by the free energy of binding. Binding can be further characterized by the enthalpy and entropy changes in the process. Methods exist to determine exact free energies, enthalpies, and entropies that are dependent only on the quality of the potential model and adequate sampling of conformational degrees of freedom. Entropy and enthalpy are roughly an order of magnitude more difficult to calculate than the free energy. This project combines a replica exchange method for enhanced sampling, designed to be efficient for protein-sized systems, with free energy calculations. This approach, replica exchange with dynamical scaling (REDS), uses two conventional simulations at different temperatures so that the entropy can be found from the temperature dependence of the free energy. A third replica is placed between them, with a modified Hamiltonian that allows it to span the temperature range of the conventional replicas. REDS provides temperature-dependent data and aids in sampling. It is applied to the bromodomain-containing protein 4 (BRD4) system. We find that for the force fields used, the free energies are accurate but the entropies and enthalpies are not, with the entropic contribution being too positive. Reproducing the entropy and enthalpy of binding appears to be a more stringent test of the force fields than reproducing the free energy.

Insights into Modeling Approaches in Chemistry: Assessing Ligand-Protein Binding Thermodynamics Based on Rigid-Flexible Model Molecules

In the field of chemistry, model compounds find extensive use for investigating complex objects. One prime example of such object is the protein-ligand supramolecular interaction. Prediction the enthalpic and entropic contribution to the free energy associated with this process, as well as the structural and dynamic characteristics of protein-ligand complexes poses considerable challenges. This review exemplifies modeling approaches used to study protein-ligand binding (PLB) thermodynamics by employing pairs of conformationally constrained/flexible model molecules. Strategically designing the model molecules can reduce the number of variables that influence thermodynamic parameters. This enables scientists to gain deeper insights into the enthalpy and entropy of PLB, which is relevant for medicinal chemistry and drug design. The model studies reviewed here demonstrate that rigidifying ligands may induce compensating changes in the enthalpy and entropy of binding. Some "rules of thumb" have started to emerge on how to minimize entropy-enthalpy compensation and design efficient rigidified or flexible ligands.

Analysis of the Structural Dynamics of Proteins in the Ligand-Unbound and -Bound States by Diffracted X-ray Tracking

Although many protein structures have been determined at atomic resolution, the majority of them are static and represent only the most stable or averaged structures in solution. When a protein binds to its ligand, it usually undergoes fluctuation and changes its conformation. One attractive method for obtaining an accurate view of proteins in solution, which is required for applications such as the rational design of proteins and structure-based drug design, is diffracted X-ray tracking (DXT). DXT can detect the protein structural dynamics on a timeline via gold nanocrystals attached to the protein. Here, the structure dynamics of single-chain Fv antibodies, helix bundle-forming de novo designed proteins, and DNA-binding proteins in both ligand-unbound and ligand-bound states were analyzed using the DXT method. The resultant mean square angular displacements (MSD) curves in both the tilting and twisting directions clearly demonstrated that structural fluctuations were suppressed upon ligand binding, and the binding energies determined using the angular diffusion coefficients from the MSD agreed well with the binding thermodynamics determined using isothermal titration calorimetry. In addition, the size of gold nanocrystals is discussed, which is one of the technical concerns of DXT.

Thermodynamic characterization of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its acyclic analogs

Cyclization of small molecules is a widely applied strategy in drug design for ligand optimization to improve affinity, as it eliminates the putative need for structural preorganization of the ligand before binding, or to improve pharmacokinetic properties. In this work, we provide a deeper insight into the binding thermodynamics of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its linear analogs. Characterization of the thermodynamic binding profiles by isothermal titration calorimetry experiments revealed an unfavorable entropy of the macrocycle compared to the open linear reference ligands. Molecular dynamic simulations and X-ray crystal structure analysis indicated only minor benefits from macrocyclization to fixate a favorable conformation, while linear ligands retained some flexibility even in the protein-bound complex structure, possibly explaining the initially surprising effect of a higher entropic penalty for the macrocyclic ligand.

Pharmacophore Modeling of Targets Infested with Activity Cliffs via Molecular Dynamics Simulation Coupled with QSAR and Comparison with other Pharmacophore Generation Methods: KDR as Case Study

Activity cliffs (ACs) are defined as pairs of structurally similar compounds with large difference in their potencies against certain biotarget. We recently proposed that potent AC members induce significant entropically-driven conformational modifications of the target that unveil additional binding interactions, while their weakly-potent counterparts are enthalpically-driven binders with little influence on the protein target. We herein propose to extract pharmacophores for ACs-infested target(s) from molecular dynamics (MD) frames of purely "enthalpic" potent binder(s) complexed within the particular target. Genetic function algorithm/machine learning (GFA/ML) can then be employed to search for the best possible combination of MD pharmacophore(s) capable of explaining bioactivity variations within a list of inhibitors. We compared the performance of this approach with established ligand-based and structure-based methods. Kinase inserts domain receptor (KDR) was used as a case study. KDR plays a crucial role in angiogenic signaling and its inhibitors have been approved in cancer treatment. Interestingly, GFA/ML selected, MD-based, pharmacophores were of comparable performances to ligand-based and structure-based pharmacophores. The resulting pharmacophores and QSAR models were used to capture hits from the national cancer institute list of compounds. The most active hit showed anti-KDR IC50 of 2.76 µM.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Prediction of Water Distributions and Displacement at Protein-Ligand Interfaces

The retention and displacement of water molecules during formation of ligand-protein interfaces play a major role in determining ligand binding. Understanding these effects requires a method for positioning of water molecules in the bound and unbound proteins and for defining water displacement upon ligand binding. We describe an algorithm for water placement and a calculation of ligand-driven water displacement in >9000 protein-ligand complexes. The algorithm predicts approximately 38% of experimental water positions within 1.0 Å and about 83% within 1.5 Å. We further show that the predicted water molecules can complete water networks not detected in crystallographic structures of the protein-ligand complexes. The algorithm was also applied to solvation of the corresponding unbound proteins, and this allowed calculation of water displacement upon ligand binding based on differences in the water network between the bound and unbound structures. We illustrate use of this approach through comparison of water displacement by structurally related ligands at the same binding site. This method for evaluation of water displacement upon ligand binding may be of value for prediction of the effects of ligand modification in drug design.

Unveiling protein dynamics in solution with field-cycling NMR relaxometry

Field-cycling NMR relaxometry is a well-established technique that can give information on molecular structure and dynamics of biological systems. It provides the nuclear relaxation rates as a function of the applied magnetic field, starting from fields as low as ~ 10^-4 T up to about 1-3 T. The profiles so collected, called nuclear magnetic relaxation dispersion (NMRD) profiles, can be extended to include the relaxation rates at the largest fields achievable with high resolution NMR spectrometers. By exploiting this wide range of frequencies, the NMRD profiles can provide information on motions occurring on time scales from 10^-6 to 10^-9 s. ¹H NMRD measurements have proved very useful also for the characterization of paramagnetic proteins, because they can help characterise a number of parameters including the number, distance and residence time of water molecules coordinated to the paramagnetic center, the reorientation correlation times and the electron spin relaxation time, and the electronic structure at the metal site.

Applications of isothermal titration calorimetry in pure and applied research from 2016 to 2020

The last 5 years have seen a series of advances in the application of isothermal titration microcalorimetry (ITC) and interpretation of ITC data. ITC has played an invaluable role in understanding multiprotein complex formation including proteolysis-targeting chimeras (PROTACS), and mitochondrial autophagy receptor Nix interaction with LC3 and GABARAP. It has also helped elucidate complex allosteric communication in protein complexes like trp RNA-binding attenuation protein (TRAP) complex. Advances in kinetics analysis have enabled the calculation of kinetic rate constants from pre-existing ITC data sets. Diverse strategies have also been developed to study enzyme kinetics and enzyme-inhibitor interactions. ITC has also been applied to study small molecule solvent and solute interactions involved in extraction, separation, and purification applications including liquid-liquid separation and extractive distillation. Diverse applications of ITC have been developed from the analysis of protein instability at different temperatures, determination of enzyme kinetics in suspensions of living cells to the adsorption of uremic toxins from aqueous streams.

Mapping Water Thermodynamics on Drug Candidates via Molecular Building Blocks: a Strategy to Improve Ligand Design and Rationalize SAR

The consideration of interactions involving water molecules in protein-ligand binding is widely appreciated in drug discovery nowadays. However, it is not ultimately clear how insights about these interactions translate into molecular design concepts. In this work, we introduce a computational strategy that, trained with high-precision experimental data, allows for the decomposition of water-related thermodynamic properties into chemically relevant building blocks (BBs) of a given ligand scaffold. For each of these BBs, a score based on solvation energy and entropy is computed, thus enabling the analysis of solvent-related affinity contributions for individual BBs. We find the nonvariable BB in a congeneric ligand pair to have a larger impact on the binding affinity than the variable part thus suggesting strong cooperative effects. Furthermore, we find enhanced solute-solvent interactions for a BB due to the presence of a C-F bond. Our investigation may be used to design drug molecules with tailored solvent thermodynamic properties.

Polytorsional-amide/carboxylates-directed Cd(ii) coordination polymers exhibiting multi-functional sensing behaviors

By rational assembly of polytorsional-amide [N,N′-bis(4-methylenepyridin-4-yl)-1,4-naphthalene dicarboxamide (L)] and polytorsional-carboxylates [H₂ADI = adipic acid, H₂PIM = pimelic acid, H₂SUB = suberic acid], three new Cd-based coordination polymers (CPs) C₃₀H₃₀CdN₄O₇ (1), C₃₁H₃₂CdN₄O₇ (2) and C_31.03H_30.55CdCl_0.24N₄O_5.52 (3) were successfully synthesized. CPs 1–2 and 3 are 2D networks and a 3D framework, which all display 3,5-connected topologies with different structural details. The effects of carboxylates with different carbon chains on the structure of the complexes were studied. Fluorescence experiments show that CPs 1–3 have good multi-functional sensing ability for metal cations (Fe³⁺), anions (MnO₄⁻, CrO₄²⁻, Cr₂O₇²⁻) and organochlorine pesticides (2,6-dichloro-4-nitroamine) with good anti-interference and recyclable characteristics. The possible sensing mechanism is also investigated in detail.

Three (3,5)-connected Cd(ii) coordination polymers induced by polytorsional-amide/carboxylates exhibiting controllable multifunctional fluorescent sensing activities.

Some thermodynamic effects of varying nonpolar surfaces in protein-ligand interactions

Understanding how making structural changes in small molecules affects their binding affinities for targeted proteins is central to improving strategies for rational drug design. To assess the effects of varying the nature of nonpolar groups upon binding entropies and enthalpies, we designed and prepared a set of Grb2-SH2 domain ligands, Ac-pTyr-Ac₆c-Asn-(CH₂)_n-R, in which the size and electrostatic nature of R groups at the pTyr+3 site were varied. The complexes of these ligands with the Grb2-SH2 domain were evaluated in a series of studies in which the binding thermodynamics were determined using isothermal titration calorimetry, and binding interactions were examined in crystallographic studies of two different complexes. Notably, adding nonpolar groups to the pTyr+3 site leads to higher binding affinities, but the magnitude and energetic origins of these effects vary with the nature of the R substituent. For example, enhancements to binding affinities using aliphatic R groups are driven by more favorable changes in binding entropies, whereas aryl R groups improve binding free energies through a combination of more favorable changes in binding enthalpies and entropies. However, enthalpy/entropy compensation plays a significant role in these associations and mitigates against any significant variation in binding free energies, which vary by only 0.8 kcal•mol^-1, with changes in the electrostatic nature and size of the R group. Crystallographic studies show that differences in ΔG° or ΔH° correlate with buried nonpolar surface area, but they do not correlate with the total number of polar or van der Waals contacts. The relative number of ordered water molecules and relative order in the side chains at pTyr+3 correlate with differences in -TΔS°. Overall, these studies show that burial of nonpolar surface can lead to enhanced binding affinities arising from dominating entropy- or enthalpy-driven hydrophobic effects, depending upon the electrostatic nature of the apolar R group.

Hydration-Induced Disorder Lowers the Energy Barriers for Methyl Rotation in Drug Molecules

The thermally activated dynamics of methyl groups are important for biochemical activity as they allow for a more efficient sampling of the energy landscape. Here, we compare methyl rotations in the dry and variously hydrated states of three primary drugs under consideration to treat the recent coronavirus disease (COVID-19), namely, hydroxychloroquine and its sulfate, dexamethasone and its sodium diphosphate, and remdesivir. We find that the main driving force behind the considerable reduction in the activation energy for methyl rotations in the hydrated state is the hydration-induced disorder in the methyl group local environments. Furthermore, the activation energy for methyl rotations in the hydration-induced disordered state is much lower than that in an isolated drug molecule, indicating that neither isolated molecules nor periodic crystalline structures can be used to analyze the potential landscape governing the side group dynamics in drug molecules. Instead, only the explicitly considered disordered structures can provide insight.

Two Methods, One Goal: Structural Differences between Cocrystallization and Crystal Soaking to Discover Ligand Binding Poses

In lead optimization, protein crystallography is an indispensable tool to analyze drug binding. Binding modes and non-covalent interaction inventories are essential to design follow-up synthesis candidates. Two protocols are commonly applied to produce protein-ligand complexes: cocrystallization and soaking. Because of its time and cost effectiveness, soaking is the more popular method. Taking eight ligand hinge binders of protein kinase A, we demonstrate that cocrystallization is superior. Particularly for flexible proteins, such as kinases, and larger ligands cocrystallization captures more reliable the correct binding pose and induced protein adaptations. The geometrical discrepancies between soaking and cocrystallization appear smaller for fragment-sized ligands. For larger flexible ligands that trigger conformational changes of the protein, soaking can be misleading and underestimates the number of possible polar interactions due to inadequate, highly impaired positions of protein amino-acid side and main chain atoms. Thus, if applicable cocrystallization should be the gold standard to study protein-ligand complexes.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

STACKED - Solvation Theory of Aromatic Complexes as Key for Estimating Drug Binding

The use of fragments to biophysically characterize a protein binding pocket and determine the strengths of certain interactions is a computationally and experimentally commonly applied approach. Almost all drug like molecules contain at least one aromatic moiety forming stacking interactions in the binding pocket. In computational drug design, the strength of stacking and the resulting optimization of the aromatic core or moiety is usually calculated using high level quantum mechanical approaches. However, as these calculations are performed in a vacuum, solvation properties are neglected. We close this gap by using Grid Inhomogeneous Solvation Theory (GIST) to describe the properties of individual heteroaromatics and complexes and thereby estimate the desolvation penalty. In our study, we investigated the solvation free energies of heteroaromatics frequently occurring in drug design projects in complex with truncated side chains of phenylalanine, tyrosine, and tryptophan. Furthermore, we investigated the properties of drug-fragments crystallized in a fragment-based lead optimization approach investigating PDE-10-A. We do not only find good correlation for the estimated desolvation penalty and the experimental binding free energy, but our calculations also allow us to predict prominent interaction sites. We highlight the importance of including the desolvation penalty of the respective heteroaromatics in stacked complexes to explain the gain or loss in affinity of potential lead compounds.

Structural basis of subtype-selective competitive antagonism for GluN2C/2D-containing NMDA receptors

N-Methyl-D-aspartate receptors (NMDARs) play critical roles in the central nervous system. Their heterotetrameric composition generates subtypes with distinct functional properties and spatio-temporal distribution in the brain, raising the possibility for subtype-specific targeting by pharmacological means for treatment of neurological diseases. While specific compounds for GluN2A and GluN2B-containing NMDARs are well established, those that target GluN2C and GluN2D are currently underdeveloped with low potency and uncharacterized binding modes. Here, using electrophysiology and X-ray crystallography, we show that UBP791 ((2S*,3R*)-1-(7-(2-carboxyethyl)phenanthrene-2-carbonyl)piperazine-2,3-dicarboxylic acid) inhibits GluN2C/2D with 40-fold selectivity over GluN2A-containing receptors, and that a methionine and a lysine residue in the ligand binding pocket (GluN2D-Met763/Lys766, GluN2C-Met736/Lys739) are the critical molecular elements for the subtype-specific binding. These findings led to development of UBP1700 ((2S*,3R*)-1-(7-(2-carboxyvinyl)phenanthrene-2-carbonyl)piperazine-2,3-dicarboxylic acid) which shows over 50-fold GluN2C/2D-selectivity over GluN2A with potencies in the low nanomolar range. Our study shows that the L-glutamate binding site can be targeted for GluN2C/2D-specific inhibition.

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.

Strong Adverse Contribution of Conformational Dynamics to Streptavidin-Biotin Binding

Molecular dynamics plays an important role for the biological function of proteins. For protein ligand interactions, changes of conformational entropy of protein and hydration layer are relevant for the binding process. Quasielastic neutron scattering (QENS) was used to investigate differences in protein dynamics and conformational entropy of ligand-bound and ligand-free streptavidin. Protein dynamics were probed both on the fast picosecond time scale using neutron time-of-flight spectroscopy and on the slower nanosecond time scale using high-resolution neutron backscattering spectroscopy. We found the internal equilibrium motions of streptavidin and the corresponding mean square displacements (MSDs) to be greatly reduced upon biotin binding. On the basis of the observed MSDs, we calculated the difference of conformational entropy ΔS_conf of the protein component between ligand-bound and ligand-free streptavidin. The rather large negative ΔS_conf value (-2 kJ mol^-1 K^-1 on the nanosecond time scale) obtained for the streptavidin tetramer seems to be counterintuitive, given the exceptionally high affinity of streptavidin-biotin binding. Literature data on the total entropy change ΔS observed upon biotin binding to streptavidin, which includes contributions from both the protein and the hydration water, suggest partial compensation of the unfavorable ΔS_conf by a large positive entropy gain of the surrounding hydration layer and water molecules that are displaced during ligand binding.

Strategies for Late-Stage Optimization: Profiling Thermodynamics by Preorganization and Salt Bridge Shielding

Structural fixation of a ligand in its bioactive conformation may, due to entropic reasons, improve affinity. We present a congeneric series of thrombin ligands with a variety of functional groups triggering preorganization prior to binding. Fixation in solution and complex formation have been characterized by crystallography, isothermal titration calorimetry (ITC), and molecular dynamics (MD) simulations. First, we show why these preorganizing modifications do not affect the overall binding mode and how key interactions are preserved. Next, we demonstrate how preorganization thermodynamics can be largely dominated by enthalpy rather than entropy because of the significant population of low-energy conformations. Furthermore, a salt bridge is shielded by actively reducing its surface exposure, thus leading to an enhanced enthalpic binding profile. Our results suggest that the consideration of the ligand solution ensemble by MD simulation is necessary to predict preorganizing modifications that enhance the binding behavior of already promising binders.

Hydration of Aromatic Heterocycles as an Adversary of π-Stacking

Hydration is one of the key players in the protein-ligand binding process. It not only influences the binding process per se, but also the drug's absorption, distribution, metabolism, and excretion properties. To gain insights into the hydration of aromatic cores, the solvation thermodynamics of 40 aromatic mono- and bicyclic systems, frequently occurring in medicinal chemistry, are investigated. Thermodynamics is analyzed with two different methods: grid inhomogeneous solvation theory (GIST) and thermodynamic integration (TI). Our results agree well with previously published experimental and computational solvation free energies. The influence of adding heteroatoms to aromatic systems and how the position of these heteroatoms impacts the compound's interactions with water is studied. The solvation free energies of these heteroaromatics are highly correlated to their gas phase interaction energies with benzene: compounds showing a high interaction energy also have a high solvation free energy value. Therefore, replacing a compound with one having a higher gas phase interaction energy might not result in the expected improvement in affinity. The desolvation costs counteract the higher stacking interactions, hence weakening or even inverting the expected gain in binding free energy.

Diamondoid Amino Acid-Based Peptide Kinase A Inhibitor Analogues

The incorporation of diamondoid amino acids (DAAs) into peptide-like drugs is a general strategy to improve lipophilicity, membrane permeability, and metabolic stability of peptidomimetic pharmaceuticals. We designed and synthesized five novel peptidic DAA-containing kinase inhibitors of protein kinase A using a sophisticated molecular dynamics protocol and solid-phase peptide synthesis. By means of a thermophoresis binding assay, NMR, and crystal structure analysis, we determined the influence of the DAAs on the secondary structure and binding affinity in comparison to the native protein kinase inhibitor, which is purely composed of proteinogenic amino acids. Affinity and binding pose are largely conserved. One variant showed 6.5-fold potency improvement, most likely related to its increased side chain lipophilicity. A second variant exhibited slightly decreased affinity presumably due to loss of hydrogen-bond contacts to surrounding water molecules of the first solvation shell.

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.

Interplay between Conformational Entropy and Solvation Entropy in Protein-Ligand Binding

Understanding the driving forces underlying molecular recognition is of fundamental importance in chemistry and biology. The challenge is to unravel the binding thermodynamics into separate contributions and to interpret these in molecular terms. Entropic contributions to the free energy of binding are particularly difficult to assess in this regard. Here we pinpoint the molecular determinants underlying differences in ligand affinity to the carbohydrate recognition domain of galectin-3, using a combination of isothermal titration calorimetry, X-ray crystallography, NMR relaxation, and molecular dynamics simulations followed by conformational entropy and grid inhomogeneous solvation theory (GIST) analyses. Using a pair of diastereomeric ligands that have essentially identical chemical potential in the unbound state, we reduced the problem of dissecting the thermodynamics to a comparison of the two protein-ligand complexes. While the free energies of binding are nearly equal for the R and S diastereomers, greater differences are observed for the enthalpy and entropy, which consequently exhibit compensatory behavior, ΔΔ H°(R - S) = -5 ± 1 kJ/mol and - TΔΔ S°(R - S) = 3 ± 1 kJ/mol. NMR relaxation experiments and molecular dynamics simulations indicate that the protein in complex with the S-stereoisomer has greater conformational entropy than in the R-complex. GIST calculations reveal additional, but smaller, contributions from solvation entropy, again in favor of the S-complex. Thus, conformational entropy apparently dominates over solvation entropy in dictating the difference in the overall entropy of binding. This case highlights an interplay between conformational entropy and solvation entropy, pointing to both opportunities and challenges in drug design.

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.