Carl Diehl

Protein Flexibility and Conformational Entropy in Ligand Design Targeting the Carbohydrate Recognition Domain of Galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein−ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. 15N and 2H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein−carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

The Carbohydrate-Binding Site in Galectin-3 Is Preorganized To Recognize a Sugarlike Framework of Oxygens: Ultra-High-Resolution Structures and Water Dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate–protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultra-high-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 Å at 100 K, 1.25 Å at 298 K) and in complex with lactose (0.86 Å) or glycerol (0.9 Å). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of β-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design.

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3.

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used 15N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the 15N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.

Starting-Condition Dependence of Order Parameters Derived from Molecular Dynamics Simulations

We have studied how backbone N-H S(2) order parameters calculated from molecular dynamics simulations depend on the method used to calculate them, the starting conditions, and the length of the simulations. Using the carbohydrate binding domain of galectin-3 in the free and lactose-bound states as a test case, we compared the calculated order parameters with experimental data from NMR relaxation. The results indicate that the sampling can be improved by using several starting structures, taking into account conformational heterogeneity reported in crystal structures. However, the improvement is rather limited, and for 93% of the dihedrals that have alternative conformations in the crystal structures, the conformational space is well sampled even if a single conformation is used as the starting structure. Moreover, the agreement with experimental data is improved when using several short simulations, rather than a single long simulation. In the present case, we find that ∼10 independent simulations provide sufficient sampling, and the ideal length of the simulations is ∼10 ns, which is ∼25% longer than the global correlation time for rotational diffusion. On the other hand, the equilibration time appears to be less important, and our results suggest that an equilibration time of 0.25 ns is sufficient. We have also compared four different methods to extract the order parameters from the simulations, namely, the autocorrelation function and isotropic reorientational eigenmode dynamics using three different window sizes. Overall, the four methods yield comparable results, but large differences between the methods may serve to pinpoint cases for which the calculated parameters are unreliable.

Structural analysis of a complex between small ubiquitin-like modifier 1 (SUMO1) and the ZZ domain of CREB-binding protein (CBP/p300) reveals a new interaction surface on SUMO

We have recently discovered that the ZZ zinc finger domain represents a novel small ubiquitin-like modifier (SUMO) binding motif. In this study we identify the binding epitopes in the ZZ domain of CBP (CREB-binding protein) and SUMO1 using NMR spectroscopy. The binding site on SUMO1 represents a unique epitope for SUMO interaction spatially opposite to that observed for canonical SUMO interaction motifs (SIMs). HADDOCK docking simulations using chemical shift perturbations and residual dipolar couplings was employed to obtain a structural model for the ZZ domain-SUMO1 complex. Isothermal titration calorimetry experiments support this model by showing that the mutation of key residues in the binding site abolishes binding and that SUMO1 can simultaneously and non-cooperatively bind both the ZZ domain and a canonical SIM motif. The binding dynamics of SUMO1 was further characterized using 15N Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersions, which define the off rates for the ZZ domain and SIM motif and show that the dynamic binding process has different characteristics for the two cases. Furthermore, in the absence of bound ligands SUMO1 transiently samples a high energy conformation, which might be involved in ligand binding.

pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK(a) values. However, pK(a) values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK(a) values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK(a) values were determined for carboxyl groups by monitoring their pH-dependent (13)C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK(a) values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK(a) values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK(a) values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.

Conformational Entropy of FK506 Binding to FKBP12 Determined by Nuclear Magnetic Resonance Relaxation and Molecular Dynamics Simulations

FKBP12 (FK506 binding protein 12 kDa) is an important drug target. Nuclear magnetic resonance (NMR) order parameters, describing amplitudes of motion on the pico- to nanosecond time scale, can provide estimates of changes in conformational entropy upon ligand binding. Here we report backbone and methyl-axis order parameters of the apo and FK506-bound forms of FKBP12, based on 15N and 2H NMR relaxation. Binding of FK506 to FKBP12 results in localized changes in order parameters, notably for the backbone of residues E54 and I56 and the side chains of I56, I90, and I91, all positioned in the binding site. The order parameters increase slightly upon FK506 binding, indicating an unfavorable entropic contribution to binding of TΔ S = -18 ± 2 kJ/mol at 293 K. Molecular dynamics simulations indicate a change in conformational entropy, associated with all dihedral angles, of TΔ S = -26 ± 9 kJ/mol. Both these values are significant compared to the total entropy of binding determined by isothermal titration calorimetry and referenced to a reactant concentration of 1 mM ( TΔ S = -29 ± 1 kJ/mol). Our results reveal subtle differences in the response to ligand binding compared to that of the previously studied rapamycin-FKBP12 complex, despite the high degree of structural homology between the two complexes and their nearly identical ligand-FKBP12 interactions. These results highlight the delicate dependence of protein dynamics on drug interactions, which goes beyond the view provided by static structures, and reinforce the notion that protein conformational entropy can make important contributions to the free energy of ligand binding.

Structure and Interactions of a Dimeric Variant of sHIP, a Novel Virulence Determinant of Streptococcus pyogenes

Streptococcus pyogenes is one of the most significant bacterial pathogens in the human population mostly causing superficial and uncomplicated infections (pharyngitis and impetigo) but also invasive and life-threatening disease. We have previously identified a virulence determinant, protein sHIP, which is secreted at higher levels by an invasive compared to a non-invasive strain of S. pyogenes. The present work presents a further characterization of the structural and functional properties of this bacterial protein. Biophysical and structural studies have shown that protein sHIP forms stable tetramers both in the crystal and in solution. The tetramers are composed of four helix-loop-helix motifs with the loop regions connecting the helices displaying a high degree of flexibility. Owing to interactions at the tetramer interface, the observed tetramer can be described as a dimer of dimers. We identified three residues at the tetramer interface (Leu84, Leu88, Tyr95), which due to largely non-polar side-chains, could be important determinants for protein oligomerization. Based on these observations, we produced a sHIP variant in which these residues were mutated to alanines. Biophysical experiments clearly indicated that the sHIP mutant appear only as dimers in solution confirming the importance of the interfacial residues for protein oligomerisation. Furthermore, we could show that the sHIP mutant interacts with intact histidine-rich glycoprotein (HRG) and the histidine-rich repeats in HRG, and inhibits their antibacterial activity to the same or even higher extent as compared to the wild type protein sHIP. We determined the crystal structure of the sHIP mutant, which, as a result of the high quality of the data, allowed us to improve the existing structural model of the protein. Finally, by employing NMR spectroscopy in solution, we generated a model for the complex between the sHIP mutant and an HRG-derived heparin-binding peptide, providing further molecular details into the interactions involving protein sHIP.