19F Fast Magic-Angle Spinning NMR Spectroscopy on Microcrystalline Complexes of Fluorinated Ligands and the Carbohydrate Recognition Domain of Galectin-3.

Structural characterization of protein-ligand binding interfaces at atomic resolution is essential for improving the design of specific and potent inhibitors. Herein, we explored fast 19F- and 1H-detected magic angle spinning NMR spectroscopy to investigate the interaction between two fluorinated ligand diastereomers with the microcrystalline galectin-3 carbohydrate recognition domain. The detailed environment around the fluorine atoms was mapped by 2D 13C-19F and 1H-19F dipolar correlation experiments and permitted characterization of the binding interface. Our results demonstrate that 19F MAS NMR is a powerful tool for detailed characterization of protein-ligand interfaces and protein interactions at the atomic level.

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.

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.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins.

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

Cosolvent Dimethyl Sulfoxide Influences Protein-Ligand Binding Kinetics via Solvent Viscosity Effects: Revealing the Success Rate of Complex Formation Following Diffusive Protein-Ligand Encounter.

Protein-ligand-exchange kinetics determines the duration of biochemical signals and consequently plays an important role in drug design. Binding studies commonly require solubilization of designed ligands in solvents such as dimethyl sulfoxide (DMSO), resulting in residual amounts of DMSO following titration of solubilized ligands into aqueous protein samples. Therefore, it is critical to establish whether DMSO influences protein-ligand binding. Here, we address the general and indirect effect of DMSO on protein-ligand binding caused by solvent viscosity, which is strongly dependent on the relative concentrations of DMSO and water. As a model system, we studied the binding of a drug-like ligand to the carbohydrate recognition domain of galectin-3 in the presence of variable amounts of DMSO. We used isothermal titration calorimetry to characterize binding thermodynamics and 15N NMR relaxation to monitor kinetics. The binding enthalpy is not affected, but we observe a subtle trend of increasingly unfavorable entropy of binding, and consequently decreased affinity, with increasing DMSO concentration. The increasing concentration of DMSO results in a reduced association rate of binding, while the dissociation rate is less affected. The observed association rate is inversely proportional to the viscosity of the DMSO-water mixture, as expected from theory, but significantly reduced from the diffusion-controlled limit. By comparing the viscosity dependence of the observed association rate with that of the theoretical diffusion-controlled association rate, we estimate the success rate of productive complex formation following an initial encounter of proteins and ligands, showing that only one out of several hundred binding "attempts" are successful.

Identification of Distinct Soluble States During Fibril Formation Using Multilinear Analysis of NMR Diffusion Data.

Protein misfolding and self-assembling into amyloid structures are associated with a number of diseases. Characterization of protein amyloid formation reactions is a challenging task as transient populations of multiple species are involved. Here we outline a method for identification and characterization of the individual soluble states during protein amyloid formation. The method combines NMR translational diffusion measurements with multilinear data analysis.

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility.

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using 15N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

Characterizing Fast Conformational Exchange of Aromatic Rings Using Residual Dipolar Couplings: Distinguishing Jumplike Flips from Other Exchange Mechanisms.

Aromatic ring flips are a hallmark of protein dynamics. They are experimentally studied by NMR spectroscopy, where recent advances have led to improved characterization across a wide range of time scales. Results on different proteins have been interpreted as continuous diffusive ring rotations or jumplike flips, leading to diverging views of the protein interior as being fluidlike or solidlike, respectively. It is challenging to distinguish between these mechanisms and other types of conformational exchange because chemical-shift-mediated line broadening provides only conclusive evidence for ring flips only if the system can be moved from the slow- to intermediate/fast-exchange regime. Moreover, whenever the chemical shift difference between the two symmetry-related sites is close to zero, it is not generally possible to determine the exchange time scale. Here we resolve these issues by measuring residual dipolar coupling (RDC)-mediated exchange contributions using NMR relaxation dispersion experiments on proteins dissolved in dilute liquid crystalline media. Excellent agreement is found between the experimental difference in RDC between the two symmetry-related sites and the value calculated from high-resolution X-ray structures, demonstrating that dynamics measured for F52 in the B1 domain of protein G reports on distinct, jumplike flips rather than other types of conformational exchange.

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 19F 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.

Transition-State Compressibility and Activation Volume of Transient Protein Conformational Fluctuations.

Proteins are dynamic entities that intermittently depart from their ground-state structures and undergo conformational transitions as a critical part of their functions. Central to understanding such transitions are the structural rearrangements along the connecting pathway, where the transition state plays a special role. Using NMR relaxation at variable temperature and pressure to measure aromatic ring flips inside a protein core, we obtain information on the structure and thermodynamics of the transition state. We show that the isothermal compressibility coefficient of the transition state is similar to that of short-chain hydrocarbon liquids, implying extensive local unfolding of the protein. Our results further indicate that the required local volume expansions of the protein can occur not only with a net positive activation volume of the protein, as expected from previous studies, but also with zero activation volume by compaction of remote void volume, when averaged over the ensemble of states.

1H R1ρ relaxation dispersion experiments in aromatic side chains.

Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic 1H R1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated 1H-13C spin pairs with vicinal 2H-12C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the 3JHH couplings and limited chemical shift differences among 1H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, kflip, agree well with previous results from 13C R1ρ relaxation dispersion experiments, and yield 1H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic1H R1ρ relaxation dispersion experiment in combination with the site-selective 1H-13C/2H-12C labeling scheme enable measurement of exchange rates up to kex = 2kflip = 80,000 s-1, and serve as a useful complement to previously developed 13C-based methods.

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.

Rapid measurement of heteronuclear transverse relaxation rates using non-uniformly sampled R1ρ accordion experiments.

Multidimensional, heteronuclear NMR relaxation methods are used extensively to characterize the dynamics of biological macromolecules. Acquisition of relaxation datasets on proteins typically requires significant measurement time, often several days. Accordion spectroscopy offers a powerful means to shorten relaxation rate measurements by encoding the "relaxation dimension" into the indirect evolution period in multidimensional experiments. Time savings can also be achieved by non-uniform sampling (NUS) of multidimensional NMR data, which is used increasingly to improve spectral resolution or increase sensitivity per unit time. However, NUS is not commonly implemented in relaxation experiments, because most reconstruction algorithms are inherently nonlinear, leading to problems when estimating signal intensities, relaxation rate constants and their error bounds. We have previously shown how to avoid these shortcomings by combining accordion spectroscopy with NUS, followed by data reconstruction using sparse exponential mode analysis, thereby achieving a dramatic decrease in the total length of longitudinal relaxation experiments. Here, we present the corresponding transverse relaxation experiment, taking into account the special considerations required for its successful implementation in the framework of the accordion-NUS approach. We attain the highest possible precision in the relaxation rate constants by optimizing the NUS scheme with respect to the Cramér-Rao lower bound of the variance of the estimated parameter, given the total number of sampling points and the spectrum-specific signal characteristics. The resulting accordion-NUS R1ρ relaxation experiment achieves comparable precision in the parameter estimates compared to conventional CPMG (Carr-Purcell-Meiboom-Gill) R2 or spin-lock R1ρ experiments while saving an order of magnitude in experiment time.

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.

Ultrastructural evidence for self-replication of Alzheimer-associated Aß42 amyloid along the sides of fibrils.

The nucleation of Alzheimer-associated Aß peptide monomers can be catalyzed by preexisting Aß fibrils. This leads to autocatalytic amplification of aggregate mass and underlies self-replication and generation of toxic oligomers associated with several neurodegenerative diseases. However, the nature of the interactions between the monomeric species and the fibrils during this key process, and indeed the ultrastructural localization of the interaction sites have remained elusive. Here we used NMR and optical spectroscopy to identify conditions that enable the capture of transient species during the aggregation and secondary nucleation of the Aß42 peptide. Cryo-electron microscopy (cryo-EM) images show that new aggregates protrude from the entire length of the progenitor fibril. These protrusions are morphologically distinct from the well-ordered fibrils dominating at the end of the aggregation process. The data provide direct evidence that self-replication through secondary nucleation occurs along the sides of fibrils, which become heavily decorated under the current solution conditions (14 µM Aß42, 20 mM sodium phosphate, 200 µM EDTA, pH 6.8).

Adsorption of unfolded Cu/Zn superoxide dismutase onto hydrophobic surfaces catalyzes its formation of amyloid fibrils.

Intracellular aggregates of superoxide dismutase 1 (SOD1) are associated with amyotrophic lateral sclerosis. In vivo, aggregation occurs in a complex and dense molecular environment with chemically heterogeneous surfaces. To investigate how SOD1 fibril formation is affected by surfaces, we used an in vitro model system enabling us to vary the molecular features of both SOD1 and the surfaces, as well as the surface area. We compared fibril formation in hydrophilic and hydrophobic sample wells, as a function of denaturant concentration and extraneous hydrophobic surface area. In the presence of hydrophobic surfaces, SOD1 unfolding promotes fibril nucleation. By contrast, in the presence of hydrophilic surfaces, increasing denaturant concentration retards the onset of fibril formation. We conclude that the mechanism of fibril formation depends on the surrounding surfaces and that the nucleating species might correspond to different conformational states of SOD1 depending on the nature of these surfaces.

Revealing Well-Defined Soluble States during Amyloid Fibril Formation by Multilinear Analysis of NMR Diffusion Data.

Amyloid fibril formation is a hallmark of neurodegenerative disease caused by protein aggregation. Oligomeric protein states that arise during the process of fibril formation often coexist with mature fibrils and are known to cause cell death in disease model systems. Progress in this field depends critically on development of analytical methods that can provide information about the mechanisms and species involved in oligomerization and fibril formation. Here, we demonstrate how the powerful combination of diffusion NMR and multilinear data analysis can efficiently disentangle the number of involved species, their kinetic rates of formation or disappearance, spectral contributions, and diffusion coefficients, even without prior knowledge of the time evolution of the process or chemical shift assignments of the various species. Using this method we identify oligomeric species that form transiently during aggregation of human superoxide dismutase 1 (SOD1), which is known to form misfolded aggregates in patients with amyotrophic lateral sclerosis. Specifically, over a time course of 42 days, during which SOD1 fibrils form, we detect the disappearance of the native monomeric species, formation of a partially unfolded intermediate in the dimer to tetramer size range, subsequent formation of a distinct similarly sized species that dominates the final spectrum detected by solution NMR, and concomitant appearance of small peptide fragments.

Rotamer Jumps, Proton Exchange, and Amine Inversion Dynamics of Dimethylated Lysine Residues in Proteins Resolved by pH-Dependent 1H and 13C NMR Relaxation Dispersion.

Post-translational methylation of lysine side chains is of great importance for protein regulation, including epigenetic control. Here, we present specific 13CHD2 labeling of dimethylated lysines as a sensitive probe of the structure, interactions, and dynamics of these groups, and outline a theoretical and experimental framework for analyzing their conformational dynamics using 1H and 13C CPMG relaxation dispersion experiments. Dimethylated lysine side chains in calcium-loaded calmodulin show a marked pH dependence of their Carr-Purcell-Meiboom-Gill (CPMG) dispersion profiles, indicating complex exchange behavior. Combined analysis of 1H and 13C CPMG relaxation dispersions requires consideration of 12-state correlated exchange of the two methyl groups due to circular three-state rotamer jumps around the Cε-Nζ axis combined with proton exchange and amine inversion. Taking into account a number of fundamental constraints, the exchange model can be reduced to include only three fitted parameters, namely, the geometric average of the rotamer-jump rate constants, the rate constant of deprotonation of Nζ, and the chemical shift difference between the trans and gauge positions of the 13C or 1H nuclei. The pH dependence indicates that protonation of the end group dramatically slows down rotamer exchange for some lysine residues, whereas deprotonation leads to rapid amine inversion coupled with rotamer scrambling. The observed variation among residues in their exchange behavior appears to depend on the structural environment of the side chain. Understanding this type of exchange process is critical to correctly interpreting NMR spectra of methylated lysine side chains. The exchange model presented here forms the basis for studying the structure and dynamics of epigenetically modified lysine side chains and perturbations caused by changes in pH or interactions with target proteins.