Vladimir A. Daragan

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.

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.

Specific Binding of Ethanol to Cholesterol in Organic Solvents

Although ethanol has been reported to affect cholesterol homeostasis in biological membranes, the molecular mechanism of action is unknown. Here, nuclear magnetic resonance (NMR) spectroscopic techniques have been used to investigate possible direct interactions between ethanol and cholesterol in various low dielectric solvents (acetone, methanol, isopropanol, DMF, DMSO, chloroform, and CCl(4)). Measurement of (13)C chemical shifts, spin-lattice and multiplet relaxation times, as well as self-diffusion coefficients, indicates that ethanol interacts weakly, yet specifically, with the HC-OH moiety and the two flanking methylenes in the cyclohexanol ring of cholesterol. This interaction is most strong in the least polar-solvent carbon tetrachloride where the ethanol-cholesterol equilibrium dissociation constant is estimated to be 2 x 10(-3) M. (13)C-NMR spin-lattice relaxation studies allow insight into the geometry of this complex, which is best modeled with the methyl group of ethanol sandwiched between the two methylenes in the cyclohexanol ring and the hydroxyl group of ethanol hydrogen bonded to the hydroxyl group of cholesterol.

Comparison of 13CαH and 15NH backbone dynamics in protein GB1

This study presents a site‐resolved experimental view of backbone CαH and NH internal motions in the 56‐residue immunoglobulin‐binding domain of streptococcal protein G, GB1. Using 13CαH and 15NH NMR relaxation data [T1, T2, and NOE] acquired at three resonance frequencies (1H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(ω) = 2ωJ(ω) to provide a model‐independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(ω) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all CαH and NH groups. Deconvolution of F(ω) 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αH internal motions are more broadly distributed on the nanosecond time scale, and larger CαH order parameters are related to correlated bond rotations for Cα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 β‐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.

13C multiplet nuclear magnetic resonance relaxation-derived ring puckering and backbone dynamics in proline-containing glycine-based peptides.

13CH2-multiplet nuclear magnetic resonance relaxation studies on proline (P)-containing glycine (G)-based peptides, GP, PG, GPG, PGG, and GPGG, provided numerous dipolar auto- and cross-correlation times for various motional model analyses of backbone and proline-ring bond rotations. Molecular dynamics simulations and bond rotation energy profiles were calculated to assess which motions could contribute most to observed relaxation phenomena. Results indicate that proline restricts backbone psi 1, psi 2, and phi 2 motions by 50% relative to those found for a polyglycine control peptide. psi 1 rotations are more restricted in the trans-proline isomer state than in the cis form. A two-state jump model best approximates proline ring puckering which in water could occur either by the C gamma endo-exo or by the C2 interconversion mechanism. The temperature dependence (5 degrees to 75 degrees C) of C beta, and C gamma, and C delta angular changes is rather flat, suggesting a near zero enthalpic contribution to the ring puckering process. In lower dielectric solvents, dimethylsulfoxide and methanol, which may mimic the hydrophobic environment within a protein, the endo-exo mechanism is preferred.

Lysine side-chain dynamics derived from 13C-multiplet NMR relaxation studies on di- and tripeptides.

Summary13C NMR relaxation data have been used to determine dipolar auto- and cross-correlation times for the di- and tripeptides GK, KG and GKG, primarily to analyze lysine side-chain motional dynamics. In general, correlation times are largest for backbone positions and decrease on going through the lysine side chain, consistent with the idea of increased mobility at Cδ and Cɛ methylenes. Correlation times, however, vary with the peptide ionization state. In the zwitterionic state of GK, for example, both auto-and cross-correlation times are at their maximum values, indicating reduced internal motions probably resulting from intramolecular electrostatic interactions. Modifying the charge state increases motional fluctuations. Activation energies determined from the temperature dependence of CH rotational autocorrelation times at neutral pH are approximately equal for glycine and lysine Cα and lysine Cβ and Cγ positions (4.1±0.2 to 4.5±0.2 kcal/mol) and tend to decrease slightly for lysine Cδ and Cɛ (3.8±0.2 to 4.3±0.2 kcal/mol). The sign of lysine side-chain cross-correlations could not be explained by using any available rotational model, including one parameterized for multiple internally restricted rotations and anisotropic overall tumbling. Molecular and stochastic dynamics calculations were performed to obtain insight into correlated internal rotations and coupled overall tumbling and internal motions. Relatively strong correlations were found for i,i+1 backbone and lysine side-chain internal bond rotations. Stochastic dynamics calculations were more successful at explaining experimentally observed correlation times. In the fully charged state, a preferred conformation was detected with an all-trans lysine side chain.

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.