Elucidation of Hydrogen Bonding Patterns in Ligand-Free, Lactose- and Glycerol-Bound Galectin-3C by Neutron Crystallography to Guide Drug Design

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Synthesis and Binding Mode Predictions of Novel Siglec-7 Ligands

Siglec-7 regulates immune cell activity and is a promising target for immunomodulation. Here, we report the discovery of novel sialic acid derivatives binding to Siglec-7. Synthesis and affinity measurements are complemented by high-quality models of sialoside-Siglec-7 complexes based on molecular dynamics (MD) simulations on the microsecond time scale. We provide details for the predicted binding modes for the new ligands, e.g., that an extension of the carbon backbone leads to a different molecular interaction pattern with the receptor and the nearby water structure than found for known Siglec-7 ligands. Further on, we uncover some shortcomings of the GLYCAM06 and GAFF2 force fields when used for the simulation of sialoside-based glycomimetics. Our results open new opportunities for the rational design of Siglec-7 inhibitors. In addition, we provide strategies on how to use and visualize MD simulations to describe and investigate sialoside-Siglec complexes in general.

Targeting galectin-3 in inflammatory and fibrotic diseases

Galectin (Gal)-3 is a β-galactoside-binding lectin emerging as a key player in cardiac, hepatic, renal, and pulmonary fibrosis and inflammation, respiratory infections caused by COVID-19, and neuroinflammatory disorders. Here, we review recent information highlighting Gal-3 as a relevant therapeutic target in these specific disease conditions. While a causal link was difficult to establish until now, we discuss how recent strategic breakthroughs allowed us to identify new-generation Gal-3 inhibitors with improved potency, selectivity, and bioavailability, and report their usefulness as valuable tools for proof-of-concept studies in various preclinical models of the aforementioned diseases, with emphasis on those actually in clinical stages. We also address critical views and suggestions intended to expand the therapeutic opportunities provided by this complex target.

Biosensor Characterization from Cratylia mollis Seed Lectin (Cramoll)-MOF and Specific Carbohydrate Interactions in an Electrochemical Model

Biosensors are small devices known for their selectivity, high specificity and sensitivity to the respective analyte, at low concentrations. We developed an electrochemical biosensor using the crystalline polymer MOF-[Cu3 (BTC)2 (H2 O)2 ]n to characterize Cratylia mollis seed lectin (Cramoll) and its interaction with free carbohydrate (glucose) and carbohydrates on the surface of rabbit erythrocytes. The electrochemical potentials presented by the exponential curves that vary from 96 to 142 mV in relation to concentrations of 10 to 20 mM of glucose are decisive for the use of the system containing gold electrode/MOF/Cramoll for the characterization of biological models due to its high sensitivity. As well as the kinetic behavior presented in the cyclic voltammograms, with a cathodic current response of 0.000 3 A for a glucose concentration of 15 mM. These results were due to the high specificity of Cramoll under these conditions, promoting stability of surface charges at the Cramoll/electrode interface. This phenomenon facilitates the monitoring of the interaction with free glucose present in the electrolyte medium by potentiometric and amperometric methods and with carbohydrates present on the surface of rabbit erythrocytes through the potentiometric method. Through scanning electron microscopy (SEM) it was possible to observe Cramoll immobilized on the MOF surface, proving the specificity of the ligand (glucose-lectin) through the morphological lectin changes in this process. This electrochemical model, Cramoll/MOF biosensor, is effective for evaluating free lectin/carbohydrate or in the erythrocyte membrane.

Structural Characterization of Rat Galectin-5, an N-Tailed Monomeric Proto-Type-like Galectin

Galectins are multi-purpose effectors acting via interactions with distinct counterreceptors based on protein-glycan/protein recognition. These processes are emerging to involve several regions on the protein so that the availability of a detailed structural characterization of a full-length galectin is essential. We report here the first crystallographic information on the N-terminal extension of the carbohydrate recognition domain of rat galectin-5, which is precisely described as an N-tailed proto-type-like galectin. In the ligand-free protein, the three amino-acid stretch from Ser2 to Ser5 is revealed to form an extra β-strand (F0), and the residues from Thr6 to Asn12 are part of a loop protruding from strands S1 and F0. In the ligand-bound structure, amino acids Ser2–Tyr10 switch position and are aligned to the edge of the β-sandwich. Interestingly, the signal profile in our glycan array screening shows the sugar-binding site to preferentially accommodate the histo-blood-group B (type 2) tetrasaccharide and N-acetyllactosamine-based di- and oligomers. The crystal structures revealed the characteristically preformed structural organization around the central Trp77 of the CRD with involvement of the sequence signature’s amino acids in binding. Ligand binding was also characterized calorimetrically. The presented data shows that the N-terminal extension can adopt an ordered structure and shapes the hypothesis that a ligand-induced shift in the equilibrium between flexible and ordered conformers potentially acts as a molecular switch, enabling new contacts in this region.

The impact of folding modes and deuteration on the atomic resolution structure of hen egg-white lysozyme

The biological function of a protein is intimately related to its structure and dynamics, which in turn are determined by the way in which it has been folded. In vitro refolding is commonly used for the recovery of recombinant proteins that are expressed in the form of inclusion bodies and is of central interest in terms of the folding pathways that occur in vivo. Here, biophysical data are reported for in vitro-refolded hydrogenated hen egg-white lysozyme, in combination with atomic resolution X-ray diffraction analyses, which allowed detailed comparisons with native hydrogenated and refolded perdeuterated lysozyme. Distinct folding modes are observed for the hydrogenated and perdeuterated refolded variants, which are determined by conformational changes to the backbone structure of the Lys97-Gly104 flexible loop. Surprisingly, the structure of the refolded perdeuterated protein is closer to that of native lysozyme than that of the refolded hydrogenated protein. These structural differences suggest that the observed decreases in thermal stability and enzymatic activity in the refolded perdeuterated and hydrogenated proteins are consequences of the macromolecular deuteration effect and of distinct folding dynamics, respectively. These results are discussed in the context of both in vitro and in vivo folding, as well as of lysozyme amyloidogenesis.

Molecular mechanism of interspecies differences in the binding affinity of TD139 to Galectin-3

Galectin-3 (Gal-3), a β-galactoside-binding lectin, has been implicated in a plethora of pathological disorders including fibrosis, inflammation, cancer and metabolic diseases. TD139-a thio-digalactoside inhibitor developed by Galecto Biotech as a potential therapeutic for idiopathic pulmonary fibrosis-is the most advanced small-molecule Gal-3 inhibitor in clinical studies. It binds to human Gal-3 with high affinity but has lower affinity towards mouse and rat homologs, which is also manifested in the differential inhibition of Gal-3 function. Using biophysical methods and high-resolution X-ray co-crystal structures of TD139 and Gal-3 proteins, we demonstrate that a single amino acid change corresponding to A146 in human Gal-3 is sufficient for the observed reduction in the binding affinity of TD139 in rodents. Site-directed mutagenesis of A146V (in human Gal-3) and V160A (in mouse Gal-3) was sufficient to interchange the affinities, mainly by affecting the off rates of the inhibitor binding. In addition, molecular dynamics simulations of both wild-type and mutant structures revealed the sustained favorable noncovalent interactions between the fluorophenyl ring and the active site A146 (human Gal-3 and mouse V160A) that corroborate the finding from biophysical studies. Current findings have ramifications in the context of optimization of drug candidates against Gal-3.

Computational Study of Potential Galectin-3 Inhibitors in the Treatment of COVID-19

Galectin-3 is a carbohydrate-binding protein and the most studied member of the galectin family. It regulates several functions throughout the body, among which are inflammation and post-injury remodelling. Recent studies have highlighted the similarity between Galectin-3's carbohydrate recognition domain and the so-called "galectin fold" present on the N-terminal domain of the S1 sub-unit of the SARS-CoV-2 spike protein. Sialic acids binding to the N-terminal domain of the Spike protein are known to be crucial for viral entry into humans, and the role of Galectin-3 as a mediator of lung fibrosis has long been the object of study since its levels have been found to be abnormally high in alveolar macrophages following lung injury. In this context, the discovery of a double inhibitor may both prevent viral entry and reduce post-infection pulmonary fibrosis. In this study, we use a database of 56 compounds, among which 37 have known experimental affinity with Galectin-3. We carry out virtual screening of this database with respect to Galectin-3 and Spike protein. Several ligands are found to exhibit promising binding affinity and interaction with the Spike protein's N-terminal domain as well as with Galectin-3. This finding strongly suggests that existing Galectin-3 inhibitors possess dual-binding capabilities to disrupt Spike-ACE2 interactions. Herein we identify the most promising inhibitors of Galectin-3 and Spike proteins, of which five emerge as potential dual effective inhibitors. Our preliminary results warrant further in vitro and in vivo testing of these putative inhibitors against SARS-CoV-2 with the hope of being able to halt the spread of the virus in the future.

In Structural Glycobiology, Deuterium provides the Details

In this issue of Structure, Gadjos et al. (2021b) determine the structure of a bacterial lectin in complex with L-fucose by neutron diffraction of both perdeuterated protein and carbohydrate ligand. The structure provides insight into lectin-ligand interactions, opening avenues for drug design targeting bacterial lectins for intervention in infectious disease.

Exploring ligand dynamics in protein crystal structures with ensemble refinement

Understanding the dynamics of ligands bound to proteins is an important task in medicinal chemistry and drug design. However, the dominant technique for determining protein-ligand structures, X-ray crystallography, does not fully account for dynamics and cannot accurately describe the movements of ligands in protein binding sites. In this article, an alternative method, ensemble refinement, is used on six protein-ligand complexes with the aim of understanding the conformational diversity of ligands in protein crystal structures. The results show that ensemble refinement sometimes indicates that the flexibility of parts of the ligand and some protein side chains is larger than that which can be described by a single conformation and atomic displacement parameters. However, since the electron-density maps are comparable and Rfree values are slightly increased, the original crystal structure is still a better model from a statistical point of view. On the other hand, it is shown that molecular-dynamics simulations and automatic generation of alternative conformations in crystallographic refinement confirm that the flexibility of these groups is larger than is observed in standard refinement. Moreover, the flexible groups in ensemble refinement coincide with groups that give high atomic displacement parameters or non-unity occupancy if optimized in standard refinement. Therefore, the conformational diversity indicated by ensemble refinement seems to be qualitatively correct, indicating that ensemble refinement can be an important complement to standard crystallographic refinement as a tool to discover which parts of crystal structures may show extensive flexibility and therefore are poorly described by a single conformation. However, the diversity of the ensembles is often exaggerated (probably partly owing to the rather poor force field employed) and the ensembles should not be trusted in detail.

Entropy-Entropy Compensation between the Protein, Ligand, and Solvent Degrees of Freedom Fine-Tunes Affinity in Ligand Binding to Galectin-3C

Molecular recognition is fundamental to biological signaling. A central question is how individual interactions between molecular moieties affect the thermodynamics of ligand binding to proteins and how these effects might propagate beyond the immediate neighborhood of the binding site. Here, we investigate this question by introducing minor changes in ligand structure and characterizing the effects of these on ligand affinity to the carbohydrate recognition domain of galectin-3, using a combination of isothermal titration calorimetry, X-ray crystallography, NMR relaxation, and computational approaches including molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). We studied a congeneric series of ligands with a fluorophenyl-triazole moiety, where the fluorine substituent varies between the ortho, meta, and para positions (denoted O, M, and P). The M and P ligands have similar affinities, whereas the O ligand has 3-fold lower affinity, reflecting differences in binding enthalpy and entropy. The results reveal surprising differences in conformational and solvation entropy among the three complexes. NMR backbone order parameters show that the O-bound protein has reduced conformational entropy compared to the M and P complexes. By contrast, the bound ligand is more flexible in the O complex, as determined by ¹⁹F NMR relaxation, ensemble-refined X-ray diffraction data, and MD simulations. Furthermore, GIST calculations indicate that the O-bound complex has less unfavorable solvation entropy compared to the other two complexes. Thus, the results indicate compensatory effects from ligand conformational entropy and water entropy, on the one hand, and protein conformational entropy, on the other hand. Taken together, these different contributions amount to entropy-entropy compensation among the system components involved in ligand binding to a target protein.

Synthesis of an Exhaustive Library of Naturally Occurring Galf-Manp and Galp-Manp Disaccharides. Toward Fingerprinting According to Ring Size by Advanced Mass Spectrometry-Based IM-MS and IRMPD

Nature offers a huge diversity of glycosidic derivatives. Among numerous structural modulations, the nature of the ring size of hexosides may induce significant differences on both biological and physicochemical properties of the glycoconjugate of interest. On this assumption, we expect that small disaccharides bearing either a furanosyl entity or a pyranosyl residue would give a specific signature, even in the gas phase. On the basis of the scope of mass spectrometry, two analytical techniques to register those signatures were considered, i.e., the ion mobility (IM) and the infrared multiple photon dissociation (IRMPD), in order to build up cross-linked databases. d-Galactose occurs in natural products in both tautomeric forms and presents all possible regioisomers when linked to d-mannose. Consequently, the four reducing Galf-Manp disaccharides as well as the four Galp-Manp counterparts were first synthesized according to a highly convergent approach, and IM-MS and IRMPD-MS data were second collected. Both techniques used afforded signatures, specific to the nature of the connectivity between the two glycosyl entities.

Visualization of hydrogen atoms in a perdeuterated lectin-fucose complex reveals key details of protein-carbohydrate interactions

Carbohydrate-binding proteins from pathogenic bacteria and fungi have been shown to be implicated in various pathological processes, where they interact with glycans present on the surface of the host cells. These interactions are part of the initial processes of infection of the host and are very important to study at the atomic level. Here, we report the room temperature neutron structures of PLL lectin from Photorhabdus laumondii in its apo form and in complex with deuterated L-fucose, which is, to our knowledge, the first neutron structure of a carbohydrate-binding protein in complex with a fully deuterated carbohydrate ligand. A detailed structural analysis of the lectin-carbohydrate interactions provides information on the hydrogen bond network, the role of water molecules, and the extent of the CH-π stacking interactions between fucose and the aromatic amino acids in the binding site.

Time-Dependent Fluorescence Spectroscopy to Quantify Complex Binding Interactions

Measuring the binding affinity for proteins that can aggregate or undergo complex binding motifs presents a variety of challenges. In this study, fluorescence lifetime measurements using intrinsic tryptophan fluorescence were performed to address these challenges and to quantify the binding of a series of carbohydrates and carbohydrate-functionalized dendrimers to recombinant human galectin-3. Collectively, galectins represent an important target for study; in particular, galectin-3 plays a variety of roles in cancer biology. Galectin-3 binding dissociation constants (K D) were quantified: lactoside (73 ± 4 μM), methyllactoside (54 ± 10 μM), and lactoside-functionalized G(2), G(4), and G(6)-PAMAM dendrimers (120 ± 58 μM, 100 ± 45 μM, and 130 ± 25 μM, respectively). The chosen examples showcase the widespread utility of time-dependent fluorescence spectroscopy for determining binding constants, including interactions for which standard methods have significant limitations.

Automated orientation of water molecules in neutron crystallographic structures of proteins

The structure and function of proteins are strongly affected by the surrounding solvent water, for example through hydrogen bonds and the hydrophobic effect. These interactions depend not only on the position, but also on the orientation, of the water molecules around the protein. Therefore, it is often vital to know the detailed orientations of the surrounding ordered water molecules. Such information can be obtained by neutron crystallography. However, it is tedious and time-consuming to determine the correct orientation of every water molecule in a structure (there are typically several hundred of them), which is presently performed by manual evaluation. Here, a method has been developed that reliably automates the orientation of a water molecules in a simple and relatively fast way. Firstly, a quantitative quality measure, the real-space correlation coefficient, was selected, together with a threshold that allows the identification of water molecules that are oriented. Secondly, the refinement procedure was optimized by varying the refinement method and parameters, thus finding settings that yielded the best results in terms of time and performance. It turned out to be favourable to employ only the neutron data and a fixed protein structure when reorienting the water molecules. Thirdly, a method has been developed that identifies and reorients inadequately oriented water molecules systematically and automatically. The method has been tested on three proteins, galectin-3C, rubredoxin and inorganic pyrophosphatase, and it is shown that it yields improved orientations of the water molecules for all three proteins in a shorter time than manual model building. It also led to an increased number of hydrogen bonds involving water molecules for all proteins.

Current status and future opportunities for serial crystallography at MAX IV Laboratory

Over the last decade, serial crystallography, a method to collect complete diffraction datasets from a large number of microcrystals delivered and exposed to an X-ray beam in random orientations at room temperature, has been successfully implemented at X-ray free-electron lasers and synchrotron radiation facility beamlines. This development relies on a growing variety of sample presentation methods, including different fixed target supports, injection methods using gas-dynamic virtual-nozzle injectors and high-viscosity extrusion injectors, and acoustic levitation of droplets, each with unique requirements. In comparison with X-ray free-electron lasers, increased beam time availability makes synchrotron facilities very attractive to perform serial synchrotron X-ray crystallography (SSX) experiments. Within this work, the possibilities to perform SSX at BioMAX, the first macromolecular crystallography beamline at  MAX IV Laboratory in Lund, Sweden, are described, together with case studies from the SSX user program: an implementation of a high-viscosity extrusion injector to perform room temperature serial crystallography at BioMAX using two solid supports - silicon nitride membranes (Silson, UK) and XtalTool (Jena Bioscience, Germany). Future perspectives for the dedicated serial crystallography beamline MicroMAX at MAX IV Laboratory, which will provide parallel and intense micrometre-sized X-ray beams, are discussed.

Accurate Backbone 13 C and 15 N Chemical Shift Tensors in Galectin-3 Determined by MAS NMR and QM/MM: Details of Structure and Environment Matter

Chemical shift tensors obtained from solid-state NMR spectroscopy are very sensitive reporters of structure and dynamics in proteins. While accurate 13 C and 15 N chemical shift tensors are accessible by magic angle spinning (MAS) NMR, their quantum mechanical calculations remain challenging, particularly for 15 N atoms. Here we compare experimentally determined backbone 13 Cα and 15 NH chemical shift tensors by MAS NMR with hybrid quantum mechanics/molecular mechanics/molecular dynamics (MD-QM/MM) calculations for the carbohydrate-binding domain of galectin-3. Excellent agreement between experimental and computed 15 NH chemical shift anisotropy values was obtained using the Amber ff15ipq force field when solvent dynamics was taken into account in the calculation. Our results establish important benchmark conditions for improving the accuracy of chemical shift calculations in proteins and may aid in the validation of protein structure models derived by MAS NMR.

Interactive model building in neutron macromolecular crystallography

This chapter aims to give an overview of the process of interactive model building in macromolecular neutron crystallography for the researcher transitioning from X-ray crystallography alone. The two most popular programs for refinement and model building, phenix.refine and Coot, respectively, are used as examples, and familiarity with the programs is assumed. Some work-arounds currently required for proper communication between the programs are described. We also discuss the appearance of nuclear density maps and how this differs from that of electron density maps. Advice is given to facilitate deposition of jointly refined neutron/X-ray structures in the Protein Data Bank.

Protein X-ray Crystallography and Drug Discovery

With the advent of structural biology in the drug discovery process, medicinal chemists gained the opportunity to use detailed structural information in order to progress screening hits into leads or drug candidates. X-ray crystallography has proven to be an invaluable tool in this respect, as it is able to provide exquisitely comprehensive structural information about the interaction of a ligand with a pharmacological target. As fragment-based drug discovery emerged in the recent years, X-ray crystallography has also become a powerful screening technology, able to provide structural information on complexes involving low-molecular weight compounds, despite weak binding affinities. Given the low numbers of compounds needed in a fragment library, compared to the hundreds of thousand usually present in drug-like compound libraries, it now becomes feasible to screen a whole fragment library using X-ray crystallography, providing a wealth of structural details that will fuel the fragment to drug process. Here, we review theoretical and practical aspects as well as the pros and cons of using X-ray crystallography in the drug discovery process.

Water structure in solution and crystal molecular dynamics simulations compared to protein crystal structures

The function of proteins is influenced not only by the atomic structure but also by the detailed structure of the solvent surrounding it. Computational studies of protein structure also critically depend on the water structure around the protein. Herein we compare the water structure obtained from molecular dynamics (MD) simulations of galectin-3 in complex with two ligands to crystallographic water molecules observed in the corresponding crystal structures. We computed MD trajectories both in a water box, which mimics a protein in solution, and in a crystallographic unit cell, which mimics a protein in a crystal. The calculations were compared to crystal structures obtained at both cryogenic and room temperature. Two types of analyses of the MD simulations were performed. First, the positions of the crystallographic water molecules were compared to peaks in the MD density after alignment of the protein in each snapshot. The results of this analysis indicate that all simulations reproduce the crystallographic water structure rather poorly. However, if we define the crystallographic water sites based on their distances to nearby protein atoms and follow these sites throughout the simulations, the MD simulations reproduce the crystallographic water sites much better. This shows that the failure of MD simulations to reproduce the water structure around proteins in crystal structures observed both in this and previous studies is caused by the problem of identifying water sites for a flexible and dynamic protein (traditionally done by overlaying the structures). Our local clustering approach solves the problem and shows that the MD simulations reasonably reproduce the water structure observed in crystals. Furthermore, analysis of the crystal MD simulations indicates a few water molecules that are close to unmodeled electron density peaks in the crystal structures, suggesting that crystal MD could be used as a complementary tool for identifying and modelling water in protein crystallography.

Shikonin induces colorectal carcinoma cells apoptosis and autophagy by targeting galectin-1/JNK signaling axis

Colorectal carcinoma (CRC) is the third most common malignant tumor pathology worldwide. Despite progress in surgical procedures and therapy options, CRC is still a considerable cause of cancer-related mortality. In this study, we tested the antitumor effects of shikonin in CRC and tried to identify its potential mechanism. The potential target, molecular mechanism as well as in vitro and in vivo antitumor effects of shikonin in CRC cells were determined by an integrative protocol including quantitative proteomics, RT-PCR, western blotting, RNA interference and overexpression, apoptosis and autophagy assays, etc. Galectin-1 was a potential target of shikonin from the iTRAQ-based proteomic analysis in shikonin-treated SW620 cell. The overexpression and RNA silencing of galectin-1 in two CRC cells suggested that the shikonin sensitivity was correlation to galectin-1 levels. The ROS accumulation induced by shikonin was important to the formation of galectin-1 dimers. Dimer galectin-1 was found to be associated with the activation of JNK and downstream apoptosis or autophagy. Moreover, through functional in vitro studies, we showed that differences in galectin-1 level affected tumor cell proliferation, migration, and invasion. In summary, shikonin induced CRC cells apoptosis and autophagy by targeting galectin-1 and JNK signaling pathway both in vitro and in vivo, which suggested a potential novel therapy target for CRC.

Glycans in drug discovery

Glycans are key players in many biological processes. They are essential for protein folding and stability and act as recognition elements in cell-cell and cell-matrix interactions. Thus, being at the heart of medically relevant biological processes, glycans have come onto the scene and are considered hot spots for biomedical intervention. The progress in biophysical techniques allowing access to an increasing molecular and structural understanding of these processes has led to the development of effective therapeutics. Indeed, strategies aimed at designing glycomimetics able to block specific lectin-carbohydrate interactions, carbohydrate-based vaccines mimicking self- and non-self-antigens as well as the exploitation of the therapeutic potential of glycosylated antibodies are being pursued. In this mini-review the most prominent contributions concerning recurrent diseases are highlighted, including bacterial and viral infections, cancer or immune-related pathologies, which certainly show the great promise of carbohydrates in drug discovery.

Aminopyrimidine-galactose hybrids are highly selective galectin-3 inhibitors

Galactopyranosides with aryl-aminopyrimidine moieties at O3 inhibit the tumor and immunity-related galectin-3 with high selectivity over other galectins.

Galectins are a family of carbohydrate recognition proteins involved in, among other things, modulating cell signalling and cell–environment interactions, giving them roles in several pathologies like cancer and idiopathic lung fibrosis. Hence, developing new galectin inhibitors with high affinity and high selectivity is important to be able to target such diseases. Most existing galectin inhibitors have a disaccharide scaffold, but there has been success as of late in developing monogalactoside inhibitors such as α-arylthioglycosides. Here, we report aminopyrimidine-derivatised galactosides as good galectin-3 inhibitors with affinities down to 1.7 μM and a more than 300-fold selectivity over galectin-1. Mutant studies replacing Arg144 in galectin-3 with lysine and serine support the hypothesis that the binding of the derivatives involves interactions with Arg144. Molecular dynamics simulations converged to stable poses of the inhibitor aminopyrimidine moiety with polar interactions with Asp148 and Ser237, while the aryl-aminopyrimidine ring stacked onto the side chain of Arg144. Hence, combining an aminopyrimidine motif with a phenyl α-thiogalactoside motif offers an attractive route towards highly selective galectin-3 inhibitors.

Refinement of protein structures using a combination of quantum-mechanical calculations with neutron and X-ray crystallographic data

Neutron crystallography is a powerful method to determine the positions of H atoms in macromolecular structures. However, it is sometimes hard to judge what would constitute a chemically reasonable model, and the geometry of H atoms depends more on the surroundings (for example the formation of hydrogen bonds) than heavy atoms, so that the empirical geometry information for the H atoms used to supplement the experimental data is often less accurate. These problems may be reduced by using quantum-mechanical calculations. A method has therefore been developed to combine quantum-mechanical calculations with joint crystallographic refinement against X-ray and neutron data. A first validation of this method is provided by re-refining the structure of the galectin-3 carbohydrate-recognition domain in complex with lactose. The geometry is improved, in particular for water molecules, for which the method leads to better-resolved hydrogen-bonding interactions. The method has also been applied to the active copper site of lytic polysaccharide monooxygenase and shows that the protonation state of the amino-terminal histidine residue can be determined.

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.

Are crystallographic B-factors suitable for calculating protein conformational entropy?

It is not possible to obtain reliable entropy estimates from crystallographic B-factors even with re-refined or room-temperature crystal structures.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

The Neutron Macromolecular Crystallography Instruments at Oak Ridge National Laboratory: Advances, Challenges, and Opportunities

The IMAGINE and MaNDi instruments, located at Oak Ridge National Laboratory High Flux Isotope Reactor and Spallation Neutron Source, respectively, are powerful tools for determining the positions of hydrogen atoms in biological macromolecules and their ligands, orienting water molecules, and for differentiating chemical states in macromolecular structures. The possibility to model hydrogen and deuterium atoms in neutron structures arises from the strong interaction of neutrons with the nuclei of these isotopes. Positions can be unambiguously assigned from diffraction studies at the 1.5–2.5 Å resolutions, which are typical for protein crystals. Neutrons have the additional benefit for structural biology of not inducing radiation damage to protein crystals, which can be critical in the study of metalloproteins. Here we review the specifications of the IMAGINE and MaNDi beamlines and illustrate their complementarity. IMAGINE is suitable for crystals with unit cell edges up to 150 Å using a quasi-Laue technique, whereas MaNDi provides neutron crystallography resources for large unit cell samples with unit cell edges up to 300 Å using the time of flight (TOF) Laue technique. The microbial culture and crystal growth facilities which support the IMAGINE and MaNDi user programs are also described.