Temperature coefficients of peptides dissolved in hexafluoroisopropanol monitor distortions of helices

Summary Temperature coefficients are widely used as an indication of solvent accessibility to amide protons. Low temperature coefficients are related to low accessibility and are often interpreted as evidence for intramolecular hydrogen bonding. Conformational shifts, i.e. the difference between chemical shifts of a particular residue in a structured and in a random-coil conformation, provide information on secondary structure. In particular, negative CH^α conformational shifts are often used to delineate the extent of helical stretches. NH conformational shifts show large oscillations within a helix that have been interpreted as the result of helix distortions affecting hydrogen bond lengths. In the ocurse of the study of differnet peptides that adopt a helical structure in the presence of the structure-inducing solvent hexafluoroisopropanol (HFIP), we have found a strong correlation between temperature coefficients and amide conformational shifts. However, contrary to the initial expectations, lower temperature coefficients were associated to amide protons involved in longer, and presumably weaker, hydrogen bonds. The correlation can be explained, however, assuming that, in helical peptides dissolved in HFIP, temperature affects the chemical shift of amide protons mainly by changing the average length of intramolecular hydrogen bonds and changes in solvent accessibility play only a secondary role under these experimental conditions. The pattern of temperature coefficients in helical peptides can therefore be used to identify short or long hydragen bonds causing bending of the helix axis.     

Temperature coefficients of peptides dissolved in hexafluoroisopropanol monitor distortions of helices

Temperature coefficients are widely used as an indication of solvent accessibility to amide protons. Low temperature coefficients are related to low accessibility and are often interpreted as evidence for intramolecular hydrogen bonding. Conformational shifts, i.e. the difference between chemical shifts of a particular residue in a structured and in a random-coil conformation, provide information on secondary structure. In particular, negative CH conformational shifts are often used to delineate the extent of helical stretches. NH conformational shifts show large oscillations within a helix that have been interpreted as the result of helix distortions affecting hydrogen bond lengths. In the course of the study of different peptides that adopt a helical structure in the presence of the structure-inducing solvent hexafluoroisopropanol (HFIP), we have found a strong correlation between temperature coefficients and amide conformational shifts. However, contrary to the initial expectations, lower temperature coefficients were associated to amide protons involved in longer, and presumably weaker, hydrogen bonds. The correlation can be explained, however, assuming that, in helical peptides dissolved in HFIP, temperature affects the chemical shift of amide protons mainly by changing the average length of intramolecular hydrogen bonds and changes in solvent accessibility play only a secondary role under these experimental conditions. The pattern of temperature coefficients in helical peptides can therefore be used to identify short or long hydrogen bonds causing bending of the helix axis.     

Short hydrogen bonds in proteins

Short hydrogen bonds are present in many chemical and biological systems. It is well known that these short hydrogen bonds are found in the active site of enzymes and aid enzyme catalysis. This study aims to systematically characterize all short hydrogen bonds from a nonredundant dataset of protein structures. The study has revealed that short hydrogen bonds are commonly found in proteins and are widely present in different regions of the protein chain, such as the backbone or side chain, and in different secondary structural regions such as helices, strands and turns. The frequency of occurrence of donors and acceptors from the charged side chains as well as from the neutral backbone atoms is equally high. This suggests that short hydrogen bonds in proteins occur either due to increased strength or due to geometrical constraints and this has been illustrated from several examples.     

Vibrational frequency shifts as a probe of hydrogen bonds: thermal expansion and glass transition of myoglobin in mixed solvents

The contribution of hydrogen bonds to protein-solvent interactions and their impact on structural flexibility and dynamics of myoglobin are discussed. The shift of vibrational peak frequencies with the temperature of myoglobin in sucrose/water and glycerol/water solutions is used to probe the expansion of the hydrogen bond network. We observe a characteristic change in the temperature slope of the O–H stretching frequency at the glass transition which correlates with the discontinuity of the thermal expansion coefficient. The temperature-difference spectra of the amide bands show the same tendency, indicating that stronger hydrogen bonding in the bulk affects the main-chain solvent interactions in parallel. However, the hydrogen bond strength decreases relative to the bulk solvent with increasing cosolvent concentration near the protein surface, which suggests preferential hydration. Weaker and/or fewer hydrogen bonds are observed at low degrees of hydration. The central O–H stretching frequency of protein hydration water is red-shifted by 40 cm^–1 relative to the bulk. The shift increases towards lower temperatures, consistent with contraction and increasing strength of the protein-water bonds. The temperature slope shows a discontinuity near 180 K. The contraction of the network has reached a critical limit which leads to frozen-in structures. This effect may represent the molecular mechanism underlying the dynamic transition observed for the mean square displacements of the protein atoms and the heme iron of myoglobin. Received: 10 July 1996 / Accepted: 10 April 1997     

The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

Temperature influence on the energy of nonspecific and specific interactions taking place between 4-aminophthalimide (4-AP) and homogeneous solvents.

Temperature induced spectral shifts of the 4-aminophthalimide (4-AP) emission spectra have been measured and compared to the predictions of the McRae solvent induced shift theory (J. Phys. Chem., 1957, 61, 562-572). Three moderately polar chloroalkanes selected as nonspecifically interacting media, and six hydrogen accepting or/and electron pair donating solvents have been used as the media in which the temperature influence on 4-AP-solvent interactions has been studied in the range of 180-320 K. Using the ab initio determined 4-AP ground state dipole moment and fitting appropriate expression originating from the mentioned theory to the shifts found in the chloroalkanes it has been possible to estimate the 4-AP excited state dipole moment, the probe excited state Onsager radius and its gas phase emission spectrum position. Using these values the thermochromic shifts of 4-AP emission spectra in hydrogen bond forming solvents have been predicted and compared to the experimental one. Temperature has been found to have different impact on the changes, upon excitation of the probe, in the mean values of the energies of different hydrogen bonds formed by 4-AP with solvents molecules.     

High-precision measurement of hydrogen bond lengths in proteins by nuclear magnetic resonance methods.

We have compared hydrogen bond lengths on enzymes derived with high precision (< or = +/- 0.05 A) from both the proton chemical shifts (delta) and the fractionation factors (phi) of the proton involved with those obtained from protein X-ray crystallography. Hydrogen bond distances derived from proton chemical shifts were obtained from a correlation of 59 O--H....O hydrogen bond lengths, measured by small molecule high-resolution X-ray crystallography, with chemical shifts determined by solid-state nuclear magnetic resonance (NMR) in the same crystals (McDermott A, Ridenour CF, Encyclopedia of NMR, Sussex, U.K.: Wiley, 1996:3820-3825). Hydrogen bond distances were independently obtained from fractionation factors that yield distances between the two proton wells in quartic double minimum potential functions (Kreevoy MM, Liang TM, J Am Chem Soc, 1980;102:3315-3322). The high-precision hydrogen bond distances derived from their corresponding NMR-measured proton chemical shifts and fractionation factors agree well with each other and with those reported in protein X-ray structures within the larger errors (+/-0.2-0.8 A) in distances obtained by protein X-ray crystallography. The increased precision in measurements of hydrogen bond lengths by NMR has provided insight into the contributions of short, strong hydrogen bonds to catalysis for several enzymatic reactions.     

Cold-induced changes in the protein ubiquitin

Conformational changes are essential for protein-protein and protein-ligand recognition. Here we probed changes in the structure of the protein ubiquitin at low temperatures in supercooled water using NMR spectroscopy. We demonstrate that ubiquitin is well folded down to 263 K, although slight rearrangements in the hydrophobic core occur. However, amide proton chemical shifts show non-linear temperature dependence in supercooled solution and backbone hydrogen bonds become weaker in the region that is most prone to cold-denaturation. Our data suggest that the weakening of the hydrogen bonds in the β-sheet of ubiquitin might be one of the first events that occur during cold-denaturation of ubiquitin. Interestingly, the same region is strongly involved in ubiquitin-protein complexes suggesting that this part of ubiquitin more easily adjusts to conformational changes required for complex formation.     

Hydrogen bonds in rubredoxins from mesophilic and hyperthermophilic organisms

The extent and strength of the hydrogen bond networks in rubredoxins from the hyperthermophile Pyrococcus furiosus (PfRd), and its mesophilic analogue Clostridium pasteurianum (CpRd), are examined and compared using NMR spectroscopy. NMR parameters examined in this study include through-hydrogen bond (h3)J(NC)(') scalar couplings and (1)H, (13)C, and (15)N chemical shifts, as well as covalent (1)J(NH) and (1)J(NC)(') scalar couplings. These parameters have allowed the characterization in solution of 12 hydrogen bonds in each protein. Despite a 83% sequence homology and a low RMSD for the backbone heavy atoms (0.648 A) in the crystalline state, subtle, but definite, changes have been identified in the detailed hydrogen-bonding patterns. CpRd shows an increased number of hydrogen bonds in the triple-stranded beta-sheet and an additional hydrogen bond in the multiple-turn segment including residues 14-32. On the other hand, PfRd exhibits an overall strengthening of N-H...O=C hydrogen bonds in the loops involved at the metal binding site as well as evidence for an additional NH...S(Cys) hydrogen bond involving the alanine residue 44. These data, as well as temperature dependence of the NMR parameters, suggest that the particular NMR hydrogen bond pattern found in the hyperthermophile rubredoxin leads to an increased stabilization at the metal binding pocket. It seems to result from a subtle redistribution of hydrogen-bonding interactions between the triple-stranded beta-sheet and the actual metal binding site.     

The structure of a glycogen phosphorylase glucopyranose spirohydantoin complex at 1.8 A resolution and 100 K: the role of the water structure and its contribution to binding.

A glucopyranose spirohydantoin (a pyranose analogue of the potent herbicide, hydantocidin) has been identified as the highest affinity glucose analogue inhibitor of glycogen phosphorylase b (GPb). In order to elucidate the structural features that contribute to the binding, the structures of GPb in the native T state conformation and in complex with glucopyranose spirohydantoin have been determined at 100 K to 2.0 A and 1.8 A resolution, respectively, and refined to crystallographic R values of 0.197 (R[free] 0.248) and 0.182 (R[free] 0.229), respectively. The low temperature structure of GPb is almost identical to that of the previously determined room temperature structure, apart from a decrease in overall atomic temperature factors ((B) room temperature GPb = 34.9 A2; (B) 100 K GPb = 23.4 A2). The glucopyranose spirohydantoin inhibitor (Ki = 3.0 microM) binds at the catalytic site and induces small changes in two key regions of the protein: the 280s loop (residues 281-286) that results in a decrease in mobility of this region, and the 380s loop (residues 377-385) that undergoes more significant shifts in order to optimize contact to the ligand. The hydantoin group, that is responsible for increasing the affinity of the glucose compound by a factor of 10(3), makes only one hydrogen bond to the protein, from one of its NH groups to the main chain oxygen of His377. The other polar groups of the hydantoin group form hydrogen bonds to five water molecules. These waters are involved in extensive networks of hydrogen bonds and appear to be an integral part of the protein structure. Analysis of the water structure at the catalytic site of the native enzyme, shows that five waters are displaced by ligand binding and that there is a significant decrease in mobility of the remaining waters on formation of the GPb-hydantoin complex. The ability of the inhibitor to exploit existing waters, to displace waters and to recruit new waters appears to be important for the high affinity of the inhibitor.     

Proton nuclear magnetic resonance study of an active pentapeptide fragment of ubiquitin

The aqueous solution conformation of Tyr-Asn-Ile-Gln-Lys (UB5) corresponding to positions 59-63 of the polypeptide, ubiquitin, has been investigated by proton NMR. Like the parent protein, UB5 induces nonspecifically both T and B lymphocyte differentiation. The various NH and CH resonances of this pentapeptide have been assigned, and its solution conformation has been probed through a study of chemical shift variations with pH, temperature dependence of amide hydrogen chemical shifts, vicinal NH--C alpha H and C alpha H--C beta H2 coupling constant data, and amide hydrogen-exchange rates. The latter were measured in H2O by using a combination of transfer of solvent saturation and saturation recovery NMR experiments. The data are compatible with the assumption of a highly motile dynamic equilibrium among different conformations for this peptide. The various secondary amide hydrogens remain essentially exposed to the solvent. The temperature-dependence study of the amide hydrogen chemical shifts also did not reveal any strong internal hydrogen bonds. A rotamer population analysis of tyrosine and asparagine side chains suggests that two of the rotomers are predominantly populated for each of these residues. From these results, a picture emerges of the dynamic conformation of UB5 in aqueous solution.     

1H-NMR investigation of the interaction of the amino terminal domain of the LexA repressor with a synthetic half-operator.

A synthetic half-operator DNA-duplex, d(GCTACTGTATGT), containing a portion of the proposed recognition sequence (CTGT) of several "SOS" genes, has been synthesized. The dodecamer has been characterized through 1H-NMR spectroscopy. Complete assignment of exchangeable hydrogen bonded imino protons has been achieved by applying 1D NOE techniques and an analysis of the temperature dependence of the chemical shifts. In order to determine the specific role of the CTGT consensus sequence in the overall recognition process, the oligonucleotide duplex has been titrated with the amino terminal DNA binding domain of the LexA repressor. The observation of substantial changes of 1H-NMR chemical shifts in the imino proton region upon interaction with the protein strongly suggests that the protein binds specifically to the operator DNA. The largest deviations of 1H-NMR chemical shifts upon protein binding have been observed for protons assigned to the CTGT segment, thus strongly suggesting a direct involvement of this sequence in the binding process. At high potassium chloride concentrations the 1H-NMR chemical shift deviations are reverted which is consistent with the known drop in the affinity constant of LexA for operator DNA at high salt concentrations.     

Comparative (1)H NMR and molecular modeling study of hydroxy protons of beta-D-Galp-(1-->4)-beta-D-GlcpNAc-(1-->2)-alpha-D-Manp-(1-->O)(CH(2))(7)CH(3) analogues in aqueous solution

The (1)H chemical shifts, coupling constants, temperature coefficients, exchange rates, and inter-residual ROEs have been measured, in aqueous solution, for the hydroxy and amine/amide proton resonances of a set of beta-D-Galp-(1-->4)-beta-D-GlcpNAc-(1-->2)-alpha-D-Manp-(1-->O)(CH(2))(7)CH(3) analogues. From the structural data, a few significant structural features could be ascertained, such as a preferential anti-conformation for the amide protons of the N-acetyl and N-propionyl groups. The introduction of systematic modifications at Gal 2-C and Gal 6-C resulted in alterations of the Gal 4-OH, Gal 3-OH, and GlcNAc 3-OH areas, since variations in chemical shifts and temperature coefficient were observed. In order to verify the possibility of hydrogen bonds, molecular dynamics simulations in the gas phase and explicit solvent were performed and correlated with the experimental data. A network of hydrogen bonds to solvent molecules was observed, but no strong intramolecular hydrogen bonding was observed.     

Hydrogen bonding on the ice-binding face of a beta-helical antifreeze protein indicated by amide proton NMR chemical shifts

The dependence of amide proton chemical shifts on temperature is used as an indication of the hydrogen bonding properties in a protein. The amide proton temperature coefficients of the beta-helical antifreeze protein from Tenebrio molitor are examined to determine their hydrogen bonding state in solution. The temperature-dependent chemical shift behavior of the amides in T. molitor antifreeze protein varies widely throughout the protein backbone; however, very subtle effects of hydrogen bonding can be distinguished using a plot of chemical shift deviation (CSD) versus the backbone amide chemical shift temperature gradient (Deltadelta/DeltaT). We show that differences between the two ranks of ice-binding threonine residues on the surface of the protein indicate that threonine residues in the left-hand rank participate in intrastrand hydrogen bonds that stabilize the flat surface required for optimal ice binding.     

Molecular dynamics simulations of helix denaturation.

An understanding of the structural transitions that an alpha-helix undergoes will help to elucidate such motions in proteins and their role in protein folding. We present the results of molecular dynamics simulations to investigate these transitions in a short polyalanine peptide (13 residues) both in vacuo and in the presence of solvent. The denaturation of this peptide was monitored as a function of temperature (ranging from 5 to 200 degrees C). In vacuo, the helical state predominated at all temperatures, whereas in solution the helix melted with increasing temperature. The peptide was predominantly helical at low temperature in solution, while at intermediate temperatures the peptide spent the bulk of the time fluctuating between different conformations with intermediate amounts of helix, e.g. not completely helical nor entirely non-helical. Many of these conformations consisted of short helical segments with intervening non-helical residues. At high temperature the peptide unfolded and adopted various collapsed unstructured states. The intrahelical hydrogen bonds that break at high temperature were not fully compensated by hydrogen bonds with water molecules in the partially unfolded forms of the peptide. Increases in temperature disrupted both the helical structure and the peptide-water interactions. Water played a major but indirect role in facilitating unfolding, as opposed to specifically competing for the intrapeptide hydrogen bonds. The implications of our results to protein folding are discussed.     

Similarities between intra- and intermolecular hydrogen bonds in RNA kissing complexes found by means of cross-correlated relaxation

The bond lengths and dynamics of intra- and intermolecular hydrogen bonds in an RNA kissing complex have been characterized by determining the NMR relaxation rates of various double- and triple-quantum coherences that involve an imino proton and two neighboring nitrogen-15 nuclei belonging to opposite bases. New experiments allow one to determine the chemical shift anisotropy of the imino protons. The bond lengths derived from dipolar relaxation and the lack of modulations of the nitrogen chemical shifts indicate that the intermolecular hydrogen bonds which hold the kissing complex together are very similar to the intramolecular hydrogen bonds in the double-stranded stem of the RNA.     

Refined crystal structure of the triphosphate conformation of H-ras p21 at 1.35 A resolution: implications for the mechanism of GTP hydrolysis.

The crystal structure of the H-ras oncogene protein p21 complexed to the slowly hydrolysing GTP analogue GppNp has been determined at 1.35 A resolution. 211 water molecules have been built into the electron density. The structure has been refined to a final R-factor of 19.8% for all data between 6 A and 1.35 A. The binding sites of the nucleotide and the magnesium ion are revealed in high detail. For the stretch of amino acid residues 61-65, the temperature factors of backbone atoms are four times the average value of 16.1 A2 due to the multiple conformations. In one of these conformations, the side chain of Gln61 makes contact with a water molecule, which is perfectly placed to be the nucleophile attacking the gamma-phosphate of GTP. Based on this observation, we propose a mechanism for GTP hydrolysis involving mainly Gln61 and Glu63 as activating species for in-line attack of water. Nucleophilic displacement is facilitated by hydrogen bonds from residues Thr35, Gly60 and Lys16. A mechanism for rate enhancement by GAP is also proposed.     

Membrane channel forming polypeptides. 270-MHz proton magnetic resonance studies of the aggregation of the 11-21 fragment of suzukacillin in organic solvents

270-MHz 1H NMR studies of the 11-21 suzukacillin fragment Boc-Gln-Aib-Leu-Aib-Gly-Leu-Aib-Pro-Val-Aib-Aib-OMe (11-G) and its analogue Boc-Ala-Aib-Leu-Aib-Gly-Leu-Aib-Pro-Val-Aib-Aib-OMe (11-A) have been carried out in CDCl3 and (CD3)2SO. The NH chemical shifts and their temperature coefficients have been measured as a function of peptide concentration in both solvents. It is established that replacement of Gln by Ala is without effect on backbone conformation. Both peptides adopt highly folded 310 helical conformations stabilized by seven intramolecular 4 leads to hydrogen bonds. Nonlinear temperature dependences are demonstrated for free NH groups in the Gln(1) peptide. Aggregation is mediated by intermolecular hydrogen bonds formed by solvent-exposed NH groups. A major role for the Gln side chain in peptide association is suggested by differences in the NMR behavior of the Gln(1) and Ala(1) peptides. For the Gln(1) peptide in CDCl3, the carboxamide side chain carbonyl group forms an intramolecular hydrogen bond to the peptide backbone, while the trans side chain NH shows evidence for intermolecular interactions. In (CD3)2SO, the cis carboxamide NH is involved in intermolecular hydrogen bonding. The possible role of the central Gln residue in stabilizing aggregates of peptide channel formers is discussed, and a model for hexameric association is postulated.     

1H-NMR Investigation of the Interaction of the Amino Terminal Domain of the LexA Repressor with a Synthetic Half-Operator

Abstract

A synthetic half-operator DNA-duplex, d(GCTACTGTATGT), containing a portion of the proposed recognition sequence (CTGT) of serveral “SOS” genes, has been synthesized. The dodecamer has been characterized through ¹H-NMR spectroscopy. Complete assignment of exchangeable hydrogen bonded imino protons has been acheived by applying 1D NOE techniques and an analysis of the temperature dependence of the chemical shifts. In order to determine the specific role of the CTGT consensus sequence in the overall recognition process, the oligonucleotide duplex has been titrated with the amino terminal DNA binding domain of the LexA repressor. The observation of substantial changes of ¹H-NMR chemical shifts in the imino proton region upon interaction with the protein strongly suggests that the protein binds specifically to the operator DNA. The largest deviations of ¹H-NMR chemical shifts upon protein binding have been observed for protons assigned to the CTGT segment, thus strongly suggesting a direct involvement of this sequence in the binding process. At high potassium chloride concentrations the ¹H-NMR chemical shift deviations are reverted which is consistent with the known drop in the affinity constant of LexA for operator DNA at high salt concentrations.