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.

The medically important drug target galectin-3 binds galactose-containing moieties on glycoproteins through an intricate pattern of hydrogen bonds to a largely polar surface-exposed binding site. All successful inhibitors of galectin-3 to date have been based on mono- or disaccharide cores closely resembling natural ligands. A detailed understanding of the H-bonding networks in these natural ligands will provide an improved foundation for the design of novel inhibitors. Neutron crystallography is an ideal technique to reveal the geometry of hydrogen bonds because the positions of hydrogen atoms are directly detected rather than being inferred from the positions of heavier atoms as in X-ray crystallography. We present three neutron crystal structures of the C-terminal carbohydrate recognition domain of galectin-3: the ligand-free form and the complexes with the natural substrate lactose and with glycerol, which mimics important interactions made by lactose. The neutron crystal structures reveal unambiguously the exquisite fine-tuning of the hydrogen bonding pattern in the binding site to the natural disaccharide ligand. The ligand-free structure shows that most of these hydrogen bonds are preserved even when the polar groups of the ligand are replaced by water molecules. The protonation states of all histidine residues in the protein are also revealed and correlate well with NMR observations. The structures give a solid starting point for molecular dynamics simulations and computational estimates of ligand binding affinity that will inform future drug design.

Hydration free energy calculations are often used to validate molecular simulation methodologies and molecular mechanics force fields. We use the free-energy perturbation method together with the AMOEBA polarizable force field and the Poltype parametrization protocol to predict the hydration free energies of 52 molecules as part of the SAMPL4 blind challenge. For comparison, similar calculations are performed using the non-polarizable General Amber force field. Against our expectations, the latter force field gives the better results compared to experiment. One possible explanation is the sensitivity of the AMOEBA results to the conformation used for parametrization.

We have calculated relative binding affinities for eight tetrafluorophenyl-triazole-thiogalactoside inhibitors of galectin-3 with the alchemical free-energy perturbation approach. We obtain a mean absolute deviation from experimental estimates of only 2–3 kJ/mol and a correlation coefficient (R2) of 0.5–0.8 for seven relative affinities spanning a range of up to 11 kJ/mol. We also studied the effect of using different methods to calculate the charges of the inhibitor and different sizes of the perturbed group (the atoms that are described by soft-core potentials and are allowed to have differing coordinates). However, the various approaches gave rather similar results and it is not possible to point out one approach as consistently and significantly better than the others. Instead, we suggest that such small and reasonable variations in the computational method can be used to check how stable the calculated results are and to obtain a more accurate estimate of the uncertainty than if performing only one calculation with a single computational setup.Electronic supplementary materialThe online version of this article (10.1007/s10822-018-0110-5) contains supplementary material, which is available to authorized users.

Galectin-3 is an important protein in molecular signalling events involving carbohydrate recognition, and an understanding of the hydrogen-bonding patterns in the carbohydrate-binding site of its C-terminal domain (galectin-3C) is important for the development of new potent inhibitors. The authors are studying these patterns using neutron crystallography. Here, the production of perdeuterated human galectin-3C and successive improvement in crystal size by the development of a crystal-growth protocol involving feeding of the crystallization drops are described. The larger crystals resulted in improved data quality and reduced data-collection times. Furthermore, protocols for complete removal of the lactose that is necessary for the production of large crystals of apo galectin-3C suitable for neutron diffraction are described. Five data sets have been collected at three different neutron sources from galectin-3C crystals of various volumes. It was possible to merge two of these to generate an almost complete neutron data set for the galectin-3C-lactose complex. These data sets provide insights into the crystal volumes and data-collection times necessary for the same system at sources with different technologies and datacollection strategies, and these insights are applicable to other systems.