Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Water-coupled low-frequency modes of myoglobin and lysozyme observed by inelastic neutron scattering.

Conformational changes of proteins often involve the relative motion of rigid structural domains. Normal mode analysis and molecular dynamics simulations of small globular proteins predict delocalized vibrations with frequencies below 20 cm(-1), which may be overdamped in solution due to solvent friction. In search of these modes, we have studied deuterium-exchanged myoglobin and lysozyme using inelastic neutron scattering in the low-frequency range at full and low hydration to modify the degree of damping. At room temperature, the hydrated samples exhibit a more pronounced quasielastic spectrum due to diffusive motions than the dehydrated samples. The analysis of the corresponding lineshapes suggests that water modifies mainly the amplitude, but not the characteristic time of fast protein motions. At low temperatures, in contrast, the dehydrated samples exhibit larger motional amplitudes than the hydrated ones. The excess scattering, culminating at 16 cm(-1), is suggested to reflect water-coupled librations of polar side chains that are depressed in the hydrated system by strong intermolecular hydrogen bonding. Both myoglobin and lysozyme exhibit ultra-low-frequency modes below 10 cm(-1) in the dry state, possibly related to the breathing modes predicted by harmonic analysis.     

Fourier transform infrared spectrometric analysis of protein conformation: effect of sampling method and stress factors.

Changes in the amide bands in Fourier transform infrared spectra of proteins are generally attributed to alterations in protein secondary structure. In this study spectra of five different globular proteins were compared in the solid and solution states recorded with several sampling techniques. Spectral differences for each protein were observed between the various sampling techniques and physical states, which could not all be explained by a change in protein secondary structure. For example, lyophilization in the absence of lyoprotectants caused spectral changes that could (partially) have been caused by the removal of hydrating water molecules rather than secondary structural changes. Moreover, attenuated total reflectance spectra of proteins in H2O were not directly comparable to transmission spectra due to the anomalous dispersion effect. Our study also revealed that the amide I, II, and III bands differ in their sensitivities to changes in protein conformation: For example, strong bands in the region 1620-1630 and 1685-1695 cm(-1) were seen in the amide I region of aggregated protein spectra. Surprisingly, absorbance of such magnitudes was not observed in the amide II and III region. It appears, therefore, that only the amide I can be used to distinguish between intra- and intermolecular beta-sheet formation. Considering the differing sensitivity of the different amide modes to structural changes, it is advisable to utilize not only the amide I band, but also the amide II and III bands, to determine changes in protein secondary structure. Finally, it is important to realize that changes in these bands may not always correspond to secondary structural changes of the proteins.     

Quantitative IR spectrophotometry of peptide compounds in water (H2O) solutions. III. Estimation of the protein secondary structure

Infrared spectra of 13 globular proteins have been obtained in the 1800-1480-cm-1 region for H2O solutions. A method for estimating protein secondary structure from the ir spectrum has been developed. The method can also be used for estimating polypeptide and fibrous protein conformation. For the globular and fibrous proteins and polypeptides analyzed, the correlation coefficients between the ir and x-ray estimates of ordered helix, disordered helix, ordered beta-structure, disordered beta-structure, turns, and remainder were 0.98, 0.80, 0.99, 0.87, 0.90, and 0.92 respectively.     

Temperature dependence of the low frequency dynamics of myoglobin. Measurement of the vibrational frequency distribution by inelastic neutron scattering.

Inelastic neutron scattering spectra of myoglobin hydrated to 0.33 g water (D2O)/g protein have been measured in the low frequency range (1-150 cm-1) at various temperatures between 100 and 350 K. The spectra at low temperatures show a well-resolved maximum in the incoherent dynamic structure factor Sinc(q, omega) at approximately 25 cm-1 and no elastic broadening. This maximum becomes gradually less distinct above 180 K due to the increasing amplitude of quasielastic scattering which extends out to 30 cm-1. The vibrational frequency distribution derived independently at 100 and 180 K are very similar, suggesting harmonic behavior at these temperatures. This result has been used to separate the vibrational motion from the quasielastic motion at temperatures above 180 K. The form of the density of states of myoglobin is discussed in relation to that of other amorphous systems, to theoretical calculations of low frequency modes in proteins, and to previous observations by electron-spin relaxation of fractal-like spectral properties of proteins. The onset of quasielastic scattering above 180 K is indicative of a dynamic transition of the system and correlates with an anomalous increase in the atomic mean-squared displacements observed by M?ssbauer spectroscopy (Parak, F., E. W. Knapp, and D. Kucheida. 1982. J. Mol. Biol. 161: 177-194.) and inelastic neutron scattering (Doster, W., S. Cusack, and W. Petry, 1989. Nature [Lond.]. 337: 754-756.) Similar behavior is observed for a hydrated powder of lysozyme suggesting that the low frequency dynamics of globular proteins have common features.     

Raman scattering of chymotrypsinogen A, ribonuclease, and ovalbumin in the aqueous solution and solid state

The Raman spectra of three globular proteins, beef pancreas chymotrypsinogen A, beef pancreas ribonuclease, and hen egg white ovalbumin have been obtained in the solid state and aqueous solution. X-ray diffraction and circular dichroism evidence have indicated that these proteins have a low α-helical content and a large fraction of the residues in the unordered and β-sheet conformation. The frequencies and intensities of the amide I and amide III lines are consistent with assignments based on the Raman spectra of polypeptides. The intense amide III lines observed in all the spectra would be expected for proteins with a low fraction of the residues in the α-helical conformation. Several spectra changes occur upon dissolution of the proteins in water and may be associated with further hydration of the proteins. The spectrum of thermally denatured chymotrypsinogen is presented. A 3 cm^–1 decrease in the frequency of the amide I line of the protein dissolved in D₂O upon heating was observed. This observation is consistent with a denaturation mechanism allowing only slight changes in the secondary structure but an increase in solvent penetration upon going from the native to the reversibly denatured state.     

Protein secondary structures in water from second-derivative amide I infrared spectra

Infrared spectra have been obtained for 12 globular proteins in aqueous solution at 20 degrees C. The proteins studied, which vary widely in the relative amounts of different secondary structures present, include myoglobin, hemoglobin, immunoglobulin G, concanavalin A, lysozyme, cytochrome c, alpha-chymotrypsin, trypsin, ribonuclease A, alcohol dehydrogenase, beta 2-microglobulin, and human class I major histocompatibility complex antigen A2. Criteria for evaluating how successfully the spectra due to liquid and gaseous water are subtracted from the observed spectrum in the amide I region were developed. Comparisons of second-derivative amide I spectra with available crystal structure data provide both qualitative and quantitative support for assignments of infrared bands to secondary structures. Band frequency assignments assigned to alpha-helix, beta-sheet, unordered, and turn structures are highly consistent among all proteins and agree closely with predictions from theory. alpha-Helix and unordered structures can each be assigned to only one band whereas multiple bands are associated with beta-sheets and turns. These findings demonstrate a method of analysis of second-derivative amide I spectra whereby the frequencies of bands due to different secondary structures can be obtained. Furthermore, the band intensities obtained provide a useful method for estimating the relative amounts of different structures.     

Native-like secondary structure of molten globules

The most common evidence for the existence of secondary structure in a globular protein is the presence of a strongly pronounced far-UV circular dichroism (CD) spectrum. Although CD spectra of native proteins are well described and their quantitative analysis is widely used, similar studies for denatured proteins have still to be done. Far-UV CD spectra of nine proteins in the native and the pH-induced molten globule states were acquired and analyzed. Singular value decomposition showed that the spectra of molten globules could be described as a superposition of at least three independent components (most likely alpha-, beta- and irregular structure). A self-consistent procedure of CD spectra analysis revealed the existence of a clear correlation between the shape of the molten globule spectra and the content of secondary structure elements in the corresponding native proteins, as determined from X-ray data. A mathematical expression of this correlation in terms of the Pierson coefficient amounts to the value of 0.9 for both the alpha-helix and the beta-structure. Thus, the secondary structure of proteins in the molten globule state is close to that in the native state.     

Frequency analysis of infrared absorption and vibrational circular dichroism of proteins in D2O solution.

The IR absorption frequencies as derived from second derivatives of the Fourier transform IR spectra of the amide I' bands of globular proteins in D2O are compared to those obtained from band fitting of the vibrational circular dichroism (VCD) spectra. The two sets of frequencies are in very good agreement, yielding consistent ranges where amide I' VCD and IR features occur. Use of VCD to complement the IR allows one to add sign information to the frequency information so that features occurring in the overlapping frequency ranges that might arise from different secondary structures can be better discriminated. From this comparison, it is clear that correlation just of the frequency of a given IR transition to secondary structure can lead to a nonunique solution. Different sign patterns were identified for correlated groups of globular proteins in restricted frequency ranges that have been previously assigned to defined secondary structural elements. Hence, different secondary structural elements must contribute band components to a given frequency range.     

Neutron scattering measurements of intact cells show changes after heat shock consistent with an increase in molecular crowding

Molecular crowding has been shown to be important in many cellular processes. The crowded environment in the cell results in a significant proportion of the cellular water being in contact with macromolecules such as proteins and DNA. These interfacial water molecules show a reduced dynamic motion that has been observed with isolated macromolecules using several biophysical techniques. Previously we investigated the inelastic neutron scattering properties of water closely associated with isolated biomolecules, and showed that interfacial water is strongly perturbed, as judged by its energy transfer spectrum. Here we have probed living cells using inelastic and quasielastic neutron scattering. We have found that mild heat stress ('heat shock'), which causes some proteins to become unfolded in the cell, results in changes in the inelastic neutron scattering in the librational region (45-130 meV). Heat shock also causes a narrowing of the quasielastic scattering peak. These changes can be understood in terms of an increase in the proportion of interfacial water molecules, and a net reduction in proton dynamics.     

Normal vibrations of -glycylglycine

A complete normal coordinate analysis, infrared absorption, and inelastic neutron scattering studies on α-glycylglycine are reported. A Urey-Bradley force field which is also valid for N-deuterated sample is obtained. Assignments in the low-frequency region are based on inelastic neutron spectrum obtained in one phonon and cubic approximation.     

Overestimated accuracy of circular dichroism in determining protein secondary structure

Circular dichroism (CD) is a spectroscopic technique widely used for estimating protein secondary structures in aqueous solution, but its accuracy has been doubted in recent work. In the present paper, the contents of nine globular proteins with known secondary structures were determined by CD spectroscopy and Fourier transform infrared spectroscopy (FTIR) in aqueous solution. A large deviation was found between the CD spectra and X-ray data, even when the experimental conditions were optimized. The content determined by FTIR was in good agreement with the X-ray crystallography data. Therefore, CD spectra are not recommended for directly calculating the content of a protein’s secondary structure.     

Dynamics of hydrogen atoms in superoxide dismutase by quasielastic neutron scattering.

The low energy dynamic of the enzyme Cu,Zn superoxide dismutase have been investigated by means of quasielastic neutron scattering in the temperature range 4-320 K. Below 200 K the scattering is purely elastic, while above this temperature a pronounced decrease in the elastic intensity is observed, together with the onset of a small quasielastic component. This behavior is similar to that previously observed in other more flexible globular proteins, and can be attributed to transitions between slightly different conformational substates of the protein tertiary structure. The presence of only a small quasielastic component, whose intensity is < or = 25% of the total spectrum, is related to the high structural rigidity of this protein.     

Plasticity of secondary structure in the N-terminal region of beta-dystroglycan

The secondary structure content of the N-terminal extracellular domain of beta-dystroglycan (a recombinant fragment extending from positions 654 to 750) has been quantitatively determined by means of CD and FTIR spectroscopies. The elements of secondary structure, namely an 8-10 residue long alpha-helix (10%) and two beta-strands (24%) have been assigned to specific amino acid sequences by means of a GOR constrained prediction method. The remaining 66% of the whole sequence is classified as turns or unordered. The temperature dependence of CD and FTIR spectra has been investigated in detail. A reversible, non-cooperative thermal transition is observed with both CD and FTIR spectroscopies up to 95 degrees C. The profile of the transition is typical of the unfolding of isolated peptides and corresponds to the progressive loss of the secondary structure elements of the protein with no evidence for collapsing phenomena, typical of globular proteins, upon heating.     

A CD determination of the alpha-helix contents of the coat proteins of four filamentous bacteriophages: fd, IKe, Pf1, and Pf3

The CD spectra of four filamentous bacteriophages--fd, IKe, Pf1, and Pf3--were analyzed to determine the alpha-helix contents of their major coat proteins. Measured spectra included the 192-nm band so that analyses could be carried out over the full wavelength range of the reference spectra for protein secondary structures available (a) from globular proteins [J.T. Yang, C.S.C. Wu, and H.M. Martinez (1986) Methods in Enzymology 130, 208-269] and (b) from poly(L-lysine) [N. Greenfield and G.D. Fasman (1960) Biochemistry 8, 4108-4116]. Extended analyses were also performed with the addition of the spectrum of a model beta-turn to the Greenfield and Fasman reference set, with the spectrum of a short alpha-helix in the Yang et al. reference set, and with an estimate of the spectrum of Trp added to both reference sets. The reference set based on the simple poly(L-lysine) polypeptide, plus a spectrum of a model beta-turn or of Trp, gave reasonably good fits to the measured spectra for all four phages and yielded the largest percentages of alpha-helix. The class I phages--fd and IKe--had large percentages of alpha-helix of 98 +/- 2 and 97 +/- 5%, respectively, while the two class II phages--Pf1 and Pf3--had similar but smaller alpha-helix contents of 83 +/- 6 and 84 +/- 2, respectively. While these alpha-helix contents were within the ranges previously reported from CD spectra of these phages in solution, they were more precise, and they indicated that the coat proteins of the intact phages have CD spectra that are probably modeled better by the reference spectra of polypeptides than by those of globular proteins.     

Estimation of protein secondary structure from circular dichroism spectra: inclusion of denatured proteins with native proteins in the analysis

We have expanded our reference set of proteins used in the estimation of protein secondary structure by CD spectroscopy from 29 to 37 proteins by including 3 additional globular proteins with known X-ray structure and 5 denatured proteins. We have also modified the self-consistent method for analyzing protein CD spectra, SELCON3, by including a new selection criterion developed by W. C. Johnson, Jr. (Proteins Struct. Funct. Genet. 35, 307-312, 1999). The secondary structure corresponding to the denatured proteins was approximated to be 90% unordered, owing to the spectral similarity of the denatured proteins and unordered structures. We examined the thermal denaturation of ribonuclease T1 by CD using both the original and expanded sets of reference proteins and obtained more consistent results with the expanded set. The expanded set of reference proteins will be helpful for the determination of protein secondary structure from protein CD spectra with higher reliability, especially of proteins with significant unordered structure content and/or in the course of denaturation.     

Conformation of the ribosomal protein S1 of Thermus thermophilus in solution under different ionic conditions

The structure of protein SI of Thermus thermophilus (M = 61 kDa) in solution at low and moderate ionic strengths (0 M and 100 mM NaCl, respectively) has been studied by small-angle X-ray and neutron scattering. It was found that protein S1 has a globular conformation under both ionic conditions. The modelling of different packing of six homologous domains of S1 on the basis of the NMR-resolved structure of one domain showed that the best fit of calculated scattering patterns from such complexes to experimental ones is observed at a compact package of the domains. The calculated value of the radius of gyration of the models is 28-29 angtroms, which is characteristic for globular proteins with a molecular mass of about 60 kDa. It was found that protein S1 has a tendency to form associates, and the type of the associate depends on ionic strength. These associates have, in general, two or three monomers at a moderate ionic strength, while at a low ionic strength the number of monomers exceeds three and they are packed in a compact manner. Strongly elongated associates were observed in neutron experiments at a moderate ionic strength in heavy water. The association of protein molecules was also confirmed by the data of dynamic light scattering. From these data, the translational diffusion coefficient of protein S1 at a moderate ionic strength was calculated to be (D20,w = (2.7 +/- 0.1) x 10(-7)cm2/s). This value is essentially smaller than the expected value (D20,w = (5.8 - 6.0) x 10(-7)cm2/s) for the S1 monomer in the globular conformation, indicating the association of protein molecules under equilibrium conditions.     

Determination of the secondary structure of proteins from the amide I band of the laser Raman spectrum

The amide I band in the laser Raman spectrum of proteins has been resolved into six components, each representing residues in a different type of secondary structure. These structure types are ordered or bihydrogen-bonded helix (believed to be located in the center of helical segments), disordered or monohydrogen-bonded helix (believed to be located at the ends of helical segments), antiparallel beta sheet, parallel beta sheet, reverse turn, and undefined. The Raman spectrum representing 100% of each type of residue conformation has been computed from the solvent-subtracted Raman spectra of ten proteins with known secondary structure, plus poly-l-lysine using a least-squares solution of the overdetermined system of equations. Linear combinations of these reference spectra were then fitted to the experimental amide I spectra of these and other proteins to estimate the fractions of residues in these conformations. Statistical tests suggest that the discrimination between bihydrogen-bonded helix and monohydrogen-bonded helix is significant as is the discrimination between parallel and antiparallel β-sheet. However, the discrimination between random structure and turns has not yet been accomplished by these studies. The absolute difference between X-ray and Raman estimates of structure for 17 protein samples is generally less than 6%. We conclude that detailed and reasonably accurate estimates of secondary structure can be derived from the amide I spectra of proteins.     

Relationship between1H and13C NMR chemical shifts and the secondary and tertiary structure of a zinc finger peptide

Summary Essentially complete assignments have been obtained for the¹H and protonated¹³C NMR spectra of the zinc finger peptide Xfin-31 in the presence and absence of zinc. The patterns observed for the¹H and¹³C chemical shifts of the peptide in the presence of zinc, relative to the shifts in the absence of zinc, reflect the zinc-mediated folding of the unstructured peptide into a compact globular structure with distinct elements of secondary structure. Chemical shifts calculated from the 3D solution structure of the peptide in the presence of zinc and the observed shifts agree to within ca. 0.2 and 0.6 ppm for the backbone C^aH and NH protons, respectively. In addition, homologous zinc finger proteins exhibit similar correlations between their¹H chemical shifts and secondary structure.     

Low-frequency Raman spectra of lysozyme crystals and oriented DNA films: dynamics of crystal water.

We observed low-frequency Raman spectra of tetragonal lysozyme crystals and DNA films, with varying water content of the samples. The spectra are fitted well by sums of relaxation modes and damped harmonic oscillators in the region from approximately 1 cm(-1) to 250 cm(-1). The relaxation modes are due to crystal water, and the distribution of relaxation times is determined. In wet samples, the relaxation time of a small part of the water molecules is a little longer than that of bulk water. The relaxation time of a considerable part of the crystal water, which belongs mainly to the secondary hydration shell, is an order of magnitude longer than that of bulk water. Furthermore, the relaxation time of some water molecules in the primary hydration shell of semidry samples is shorter than we expected. Thus we have shown that low-frequency Raman measurements combined with properly oriented samples can give specific information on the dynamics of hydration water in the ps range. On the other hand, we concluded, based on polarized Raman spectra of lysozyme crystals, that the damped oscillators correspond to essentially intramolecular vibrational modes.