Pro4 prolyl peptide bond isomerization in human galectin-7 modulates the monomer-dimer equilibrum to affect function.

Human galectin-7 (Gal-7; also termed p53-induced gene 1 product) is a multifunctional effector by productive pairing with distinct glycoconjugates and protein counter-receptors in the cytoplasm and nucleus, as well as on the cell surface. Its structural analysis by NMR spectroscopy detected doubling of a set of particular resonances, an indicator of Gal-7 existing in two conformational states in slow exchange on the chemical shift time scale. Structural positioning of this set of amino acids around the P4 residue and loss of this phenomenon in the bioactive P4L mutant indicated cis-trans isomerization at this site. Respective resonance assignments confirmed our proposal of two Gal-7 conformers. Mapping hydrogen bonds and considering van der Waals interactions in molecular dynamics simulations revealed a structural difference for the N-terminal peptide, with the trans-state being more exposed to solvent and more mobile than the cis-state. Affinity for lactose or glycan-inhibitable neuroblastoma cell surface contact formation was not affected, because both conformers associated with an overall increase in order parameters (S2). At low µM concentrations, homodimer dissociation is more favored for the cis-state of the protein than its trans-state. These findings give direction to mapping binding sites for protein counter-receptors of Gal-7, such as Bcl-2, JNK1, p53 or Smad3, and to run functional assays at low concentration to test the hypothesis that this isomerization process provides a (patho)physiologically important molecular switch for Gal-7.

Lactose binding to human galectin-7 (p53-induced gene 1) induces long-range effects through the protein resulting in increased dimer stability and evidence for positive cooperativity.

The product of p53-induced gene 1 is a member of the galectin family, i.e., galectin-7 (Gal-7). To move beyond structural data by X-ray diffraction, we initiated the study of the lectin by nuclear magnetic resonance (NMR) and circular dichroism spectroscopies, and molecular dynamics (MD) simulations. In concert, our results indicate that lactose binding to human Gal-7 induces long-range effects (minor conformational shifts and changes in structural dynamics) throughout the protein that result in stabilization of the dimer state, with evidence for positive cooperativity. Monte Carlo fits of (15)N-Gal-7 HSQC titrations with lactose using a two-site model yield K1 = 0.9 ± 0.6 × 10(3) M(-1) and K2 = 3.4 ± 0.8 × 10(3) M(-1). Ligand binding-induced stabilization of the Gal-7 dimer was supported by several lines of evidence: MD-based calculations of interaction energies between ligand-loaded and ligand-free states, gel filtration data and hetero-FRET spectroscopy that indicate a highly reduced tendency for dimer dissociation in the presence of lactose, CD-based thermal denaturation showing that the transition temperature of the lectin is significantly increased in the presence of lactose, and saturation transfer difference (STD) NMR using a molecular probe of the monomer state whose presence is diminished in the presence of lactose. MD simulations with the half-loaded ligand-bound state also provided insight into how allosteric signaling may occur. Overall, our results reveal long-range effects on Gal-7 structure and dynamics, which factor into entropic contributions to ligand binding and allow further comparisons with other members of the galectin family.

Heat capacities and a snapshot of the energy landscape in protein GB1 from the pre-denaturation temperature dependence of backbone NH nanosecond fluctuations.

Protein stability is usually characterized calorimetrically by a melting temperature and related thermodynamic parameters. Despite its importance, the microscopic origin of the melting transition and the relationship between thermodynamic stability and dynamics remains a mystery. Here, NMR relaxation parameters were acquired for backbone 15NH groups of the 56 residue immunoglobulin-binding domain of streptococcal protein G over a pre-denaturation temperature range of 5-50 degrees C. Relaxation data were analyzed using three methods: the standard three-Lorentzian model free approach; the F(omega)=2omegaJ(omega) spectral density approach that yields motional correlation time distributions, and a new approach that determines frequency-dependent order parameters. Regardless of the method of analysis, the temperature dependence of internal motional correlation times and order parameters is essentially the same. Nanosecond time-scale internal motions are found for all NHs in the protein, and their temperature dependence yields activation energies ranging up to about 33kJ/mol residue. NH motional barrier heights are structurally correlated, with the largest energy barriers being found for residues in the most "rigid" segments of the fold: beta-strands 1 and 4 and the alpha-helix. Trends in this landscape also parallel the free energy of folding-unfolding derived from hydrogen-deuterium (H-D) exchange measurements, indicating that the energetics for internal motions occurring on the nanosecond time-scale mirror those occurring on the much slower time-scale of H-D exchange. Residual heat capacities, derived from the temperature dependence of order parameters, range from near zero to near 100J/mol K residue and correlate with this energy landscape. These results provide a unique picture of this protein's energy landscape and a relationship between thermodynamic stability and dynamics that suggests thermosensitive regions in the fold that could initiate the melting process.

Comparison of (13)C(alpha)H and (15)NH backbone dynamics in protein GB1.

This study presents a site-resolved experimental view of backbone C(alpha)H and NH internal motions in the 56-residue immunoglobulin-binding domain of streptococcal protein G, GB1. Using (13)C(alpha)H and (15)NH NMR relaxation data [T(1), T(2), and NOE] acquired at three resonance frequencies ((1)H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(omega) = 2omegaJ(omega) to provide a model-independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(omega) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all C(alpha)H and NH groups. Deconvolution of F(omega) curves effectively separates overall tumbling and internal motional correlation time distributions to yield more accurate order parameters than determined by using standard model free approaches. Compared to NH groups, C(alpha)H internal motions are more broadly distributed on the nanosecond time scale, and larger C(alpha)H order parameters are related to correlated bond rotations for C(alpha)H fluctuations. Motional parameters for NH groups are more structurally correlated, with NH order parameters, for example, being larger for residues in more structured regions of beta-sheet and helix and generally smaller for residues in the loop and turns. This is most likely related to the observation that NH order parameters are correlated to hydrogen bonding. This study contributes to the general understanding of protein dynamics and exemplifies an alternative and easier way to analyze NMR relaxation data.

Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding.

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.

A simple method to measure 13CH2 heteronuclear dipolar cross-correlation spectral densities.

Here, we report a method to simultaneously determine CH2 cross-correlation spectral densities and T1 relaxation times in the laboratory and rotating frames. To accomplish this, we have employed an indirect approach that is based on measurement of differences in relaxation rates acquired with and without cross-correlation terms. The new method, which can be employed using multidimensional NMR and standard relaxation pulse sequences, is validated experimentally by investigation of a selectively 13C-enriched hexadecapeptide and the uniformly 13C-enriched immunoglobulin-binding domain of streptococcal protein G (GB1). Use of this approach makes determination of CH2 cross-correlation spectral densities in uniformly 13C-enriched proteins now routine and provides novel information concerning their internal motions.

Protein dynamics using frequency-dependent order parameters from analysis of NMR relaxation data.

A novel approach is described to analyze NMR relaxation data on proteins. This method introduces the frequency-dependent order parameter, S(2)(omega), in order to estimate contributions to the generalized order parameter S(2) from different motional frequencies occurring on the picosecond to nanosecond time scales. S(2)(omega) is defined as the sum of a specified set of weighting coefficients from the Lorentzian expansion of the spectral density function. 15N NMR relaxation data (500, 600, and 800 MHz) on protein GB1 exemplify the method. Using this approach provides information on motional restrictions over specific frequency or time scale ranges and provides a normalized comparison of motional restrictions between proteins having different overall tumbling correlation times.

Lactose binding to galectin-1 modulates structural dynamics, increases conformational entropy, and occurs with apparent negative cooperativity.

Galectins are a family of lectins with a conserved carbohydrate recognition domain that interacts with beta-galactosides. By binding cell surface glycoconjugates, galectin-1 (gal-1) is involved in cell adhesion and migration processes and is an important regulator of tumor angiogenesis. Here, we used heteronuclear NMR spectroscopy and molecular modeling to investigate lactose binding to gal-1 and to derive solution NMR structures of gal-1 in the lactose-bound and unbound states. Structure analysis shows that the beta-strands and loops around the lactose binding site, which are more open and dynamic in the unbound state, fold in around the bound lactose molecule, dampening internal motions at that site and increasing motions elsewhere throughout the protein to contribute entropically to the binding free energy. CD data support the view of an overall more open structure in the lactose-bound state. Analysis of heteronuclear single quantum coherence titration binding data indicates that lactose binds the two carbohydrate recognition domains of the gal-1 dimer with negative cooperativity, in that the first lactose molecule binds more strongly (K(1)=21+/-6 x 10(3) M(-1)) than the second (K(2)=4+/-2 x 10(3) M(-1)). Isothermal calorimetry data fit using a sequential binding model present a similar picture, yielding K(1)=20+/-10 x 10(3) M(-1) and K(2)=1.67+/-0.07 x 10(3) M(-1). Molecular dynamics simulations provide insight into structural dynamics of the half-loaded lactose state and, together with NMR data, suggest that lactose binding at one site transmits a signal through the beta-sandwich and loops to the second binding site. Overall, our results provide new insight into gal-1 structure-function relationships and to protein-carbohydrate interactions in general.