Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit.

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Conformational Dynamics of Phytoglobin BvPgb1.2 from Beta vulgaris ssp. vulgaris.

Plant hemoglobins, often referred to as phytoglobins, play important roles in abiotic stress tolerance. Several essential small physiological metabolites can be bound to these heme proteins. In addition, phytoglobins can catalyze a range of different oxidative reactions in vivo. These proteins are often oligomeric, but the degree and relevance of subunit interactions are largely unknown. In this study, we delineate which residues are involved in dimer formation of a sugar beet phytoglobin type 1.2 (BvPgb1.2) using NMR relaxation experiments. E. coli cells harboring a phytoglobin expression vector were cultivated in isotope-labeled (2H, 13C and 15N) M9 medium. The triple-labeled protein was purified to homogeneity using two chromatographic steps. Two forms of BvPgb1.2 were examined, the oxy-form and the more stable cyanide-form. Using three-dimensional triple-resonance NMR experiments, sequence-specific assignments for CN-bound BvPgb1.2 were achieved for 137 backbone amide cross-peaks in the 1H-15N TROSY spectrum, which amounts to 83% of the total number of 165 expected cross-peaks. A large proportion of the non-assigned residues are located in α-helixes G and H, which are proposed to be involved in protein dimerization. Such knowledge around dimer formation will be instrumental for developing a better understanding of phytoglobins' roles in planta.

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.

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation.

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from 15N relaxation and side-chain dynamics from 13C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

Mapping the energy landscape of protein-ligand binding via linear free energy relationships determined by protein NMR relaxation dispersion.

Biochemical signaling is mediated by complexes between macromolecular receptors and their ligands, with the duration of the signal being directly related to the lifetime of the ligand-receptor complex. In the field of drug design, the recognition that drug efficacy in vivo depends on the lifetime of the drug-protein complex has spawned the concept of designing drugs with particular binding kinetics. To advance this field it is critical to investigate how the molecular details of designed ligands might affect the binding kinetics, as well as the equilibrium binding constant. Here we use protein NMR relaxation dispersion to determine linear free energy relationships involving the on- and off-rates and the affinity for a series of congeneric ligands targeting the carbohydrate recognition domain of galectin-3. Using this approach we determine the energy landscape and the position of the transition state along the reaction coordinate of protein-ligand binding. The results show that ligands exhibiting reduced off-rates achieve this by primarily stabilizing the bound state, but do not affect the transition state to any greater extent. The transition state forms early, that is, it is located significantly closer to the free state than to the bound state, suggesting a critical role of desolvation. Furthermore, the data suggest that different subclasses of ligands show different behavior with respect to these characteristics.

Structural characterization of the microbial enzyme urocanate reductase mediating imidazole propionate production.

The human microbiome can produce metabolites that modulate insulin signaling. Type 2 diabetes patients have increased circulating concentrations of the microbially produced histidine metabolite, imidazole propionate (ImP) and administration of ImP in mice resulted in impaired glucose tolerance. Interestingly, the fecal microbiota of the patients had increased capacity to produce ImP, which is mediated by the bacterial enzyme urocanate reductase (UrdA). Here, we describe the X-ray structures of the ligand-binding domains of UrdA in four different states, representing the structural transitions along the catalytic reaction pathway of this unexplored enzyme linked to disease in humans. The structures in combination with functional data provide key insights into the mechanism of action of UrdA that open new possibilities for drug development strategies targeting type 2 diabetes.

Systematic Tuning of Fluoro-galectin-3 Interactions Provides Thiodigalactoside Derivatives with Single-Digit nM Affinity and High Selectivity.

Symmetrical and asymmetrical fluorinated phenyltriazolyl-thiodigalactoside derivatives have been synthesized and evaluated as inhibitors of galectin-1 and galectin-3. Systematic tuning of the phenyltriazolyl-thiodigalactosides' fluoro-interactions with galectin-3 led to the discovery of inhibitors with exceptional affinities (Kd down to 1-2 nM) in symmetrically substituted thiodigalactosides as well as unsurpassed combination of high affinity (Kd 7.5 nM) and selectivity (46-fold) over galectin-1 for asymmetrical thiodigalactosides by carrying one trifluorphenyltriazole and one coumaryl moiety. Studies of the inhibitor-galectin complexes with isothermal titration calorimetry and X-ray crystallography revealed the importance of fluoro-amide interaction for affinity and for selectivity. Finally, the high affinity of the discovered inhibitors required two competitive titration assay tools to be developed: a new high affinity fluorescent probe for competitive fluorescent polarization and a competitive ligand optimal for analyzing high affinity galectin-3 inhibitors with competitive isothermal titration calorimetry.

A De Novo Designed Coiled-Coil Peptide with a Reversible pH-Induced Oligomerization Switch.

Protein conformational switches have many useful applications but are difficult to design rationally. Here we demonstrate how the isoenergetic energy landscape of higher-order coiled coils can enable the formation of an oligomerization switch by insertion of a single destabilizing element into an otherwise stable computationally designed scaffold. We describe a de novo designed peptide that was discovered to switch between a parallel symmetric pentamer at pH 8 and a trimer of antiparallel dimers at pH 6. The transition between pentamer and hexamer is caused by changes in the protonation states of glutamatic acid residues with highly upshifted pKa values in both oligomer forms. The drastic conformational change coupled with the narrow pH range makes the peptide sequence an attractive candidate for introduction of pH sensing into other proteins. The results highlight the remarkable ability of simple-α helices to self-assemble into a vast range of structural states.