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.     

The low-temperature heat capacity of solid proteins.

Several harmonic models of protein fluctuations are used to calculate the heat capacity. They get the spectral density of conformational modes from inelastic neutron scattering, normal mode calculations, or macroscopic elasticity (Debye model). It is assumed that the low-frequency spectral density depends only weakly on temperature and protein species. The Debye model predicts temperatures below which modes are primarily in their ground states: 10 and 80 K for the lattice and conformational modes, respectively. The models differ most below 100 K. The mode calculations yield the most accurate predictions, though all three models are within twofold of the data. The heat capacity has the power law form aTb for T less than 30 K. The experimental b's of proteins are 1.6-1.8, and the theoretical, 1.1-1.3. One possible explanation for the discrepancy is the occurrence of transitions between discrete conformations. All of the models approach the measured data in the range 100-200 K. They are very similar above 200 K, where the heat capacity includes significant contributions from bond stretching and bending. This masks the possible anharmonic behavior of the conformational modes. Hydration substantially increases the heat capacity above 200 K. This effect seems to be a consequence of conformational transitions that have higher energy than the ones seen with low hydration. The analysis also predicts that denaturation with constant hydration produces a negligible increase of heat capacity. The larger increment in solution arises from the different hydration of the folded and unfolded states, and is responsible for the existence of cold denaturation. This phenomenon is thus predicted not to occur when the hydration is constant.     

Metallodrug-protein interaction probed by synchrotron terahertz and neutron scattering spectroscopy

This experimental work applied coherent synchrotron-radiation terahertz spectroscopy and inelastic neutron scattering to address two processes directly associated with the mode of action of metal-based anticancer agents that can severely undermine chemotherapeutic treatment: drug binding to human serum albumin, occurring during intravenous drug transport, and intracellular coordination to thiol-containing biomolecules (such as metallothioneins) associated with acquired drug resistance. Cisplatin and two dinuclear platinum (Pt)- and palladium (Pd)-polyamine agents developed by this research group, which have yielded promising results toward some types of human cancers, were investigated. Complementary synchrotron-radiation-terahertz and inelastic neutron scattering data revealed protein metalation, through S- and N-donor ligands from cysteine, methionine, and histidine residues. A clear impact of the Pt and Pd agents was evidenced, drug binding to albumin and metallothionein having been responsible for significant changes in the overall protein conformation, as well as for an increased flexibility and possible aggregation.     

Vibrational mode analysis of guanine by neutron inelastic scattering

The low-temperature neutron inelastic spectrum of guanine has been measured. In order to assign the intense peaks observed in this spectrum, a normal mode analysis has been performed, using the Wilson GF-method. The theoretical treatment is based on a non-redundant set of internal coordinates, and a simplified valence force-field approximation. Only the fundamentals have been considered for simulating the internal vibrational mode spectrum. The calculations account for the spectral shape as well as the main observed peaks. Offprint requests to: M. Ghomi     

Direct measurement of hydration-related dynamic changes in lysozyme using inelastic neutron scattering spectroscopy

Inelastic neutron scattering spectroscopy is used to investigate dynamic changes in lysozyme powder at two different low D2O hydrations (0.07g D2O/g protein and 0.20 g D2O/g protein). In the higher hydration sample, the inelastic scattering between 0.8 and 4.0 cm-1 energy transfer is increased and the elastic scattering is decreased. The decreased elastic scattering suggests increased atomic amplitudes of motion and the increased 0.8 to 4.0 cm-1 scattering suggests increased motions in this frequency range. Comparison with normal mode models of lysozyme dynamics shows that the inelastic difference occurs in the frequency region predicted for the lowest frequency, largest amplitude, global modes of the molecular [M. Levitt, C. Sander and P.S. Stern, J. Mol. Biol. 181, 423 (1985). B. Brooks and M. Karplus, Proc. Natl. Acad. Sci (U.S.A) 82, 4995 (1985), R.E. Bruccoleri, M. Karplus and J.A. McCammon, Biopolymers 25 1767 (1986)]. Our results are consistent with a model in which an increased number of low frequency global modes are present in the higher hydrated sample.     

Effect of Urey-Bradley-Shimanouchi force field on the harmonic dynamics of proteins.

A normal mode analysis of bovine pancreatic trypsin inhibitor is carried out by using a Urey-Bradley-Shimanouchi potential energy function. The density of vibrational states, the magnitudes, and time scales of the atomic fluctuations are compared with experimental and theoretical results obtained by the more commonly used potential energy functions. The atomic fluctuations of Lys-15 are subject to extensive considerations as this residue is buried in the trypsin specificity pocket. It is found that Arg-17 is likely to be of importance in order to understand the way BPTI binds on trypsin.     

Effects of macromolecular crowding on an intrinsically disordered protein characterized by small-angle neutron scattering with contrast matching

Small-angle neutron scattering was used to examine the effects of molecular crowding on an intrinsically disordered protein, the N protein of bacteriophage λ, in the presence of high concentrations of a small globular protein, bovine pancreatic trypsin inhibitor (BPTI). The N protein was labeled with deuterium, and the D₂O concentration of the solvent was adjusted to eliminate the scattering contrast between the solvent and unlabeled BPTI, leaving only the scattering signal from the unfolded protein. The scattering profile observed in the absence of BPTI closely matched that predicted for an ensemble of random conformations. With BPTI added to a concentration of 65 mg/mL, there was a clear change in the scattering profile representing an increase in the mass fractal dimension of the unfolded protein, from 1.7 to 1.9, as expected if crowding favors more compact conformations. The crowding protein also inhibited aggregation of the unfolded protein. At 130 mg/mL BPTI, however, the fractal dimension was not significantly different from that measured at the lower concentration, contrary to the predictions of models that treat the unfolded conformations as convex particles. These results are reminiscent of the behavior of polymers in concentrated melts, suggesting that these synthetic mixtures may provide useful insights into the properties of unfolded proteins under crowding conditions.     

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.     

Quantum Behavior of Water Protons in Protein Hydration Shell

Quantum effects on the water proton dynamics over the surface of a hydrated protein are measured by means of broadband dielectric spectroscopy and deep inelastic neutron scattering. Dielectric spectroscopy indicates a reduced energy barrier for a hydrogenated protein sample compared to a deuterated one, along with a large and temperature-dependent isotopic ratio, in good agreement with theoretical studies. Recent deep inelastic neutron scattering data have been reanalyzed, and now show that the momentum distribution of water protons reflects a characteristic delocalization at ambient temperatures. These experimental findings might have far-reaching implications for enzymatic catalysis involving proton transfer processes, as in the case of the lysozyme protein studied in this report.     

Influence of pressure on structure and dynamics of bovine pancreatic trypsin inhibitor (BPTI): small angle and quasi-elastic neutron scattering studies

We have studied the influence of pressure on structure and dynamics of a small protein belonging to the enzymatic catalysis: the bovine pancreatic trypsin inhibitor (BPTI). Using a copper-beryllium high-pressure cell, we have performed small angle neutron scattering (SANS) experiment on NEAT spectrometer at HMI (Berlin, Germany). In the SANS configuration, the evolution of the radius of gyration and of the shape of the protein under pressures up to 6,000 bar has been studied. When increasing pressure from atmospheric pressure up to 6,000 bar, the pressure effects on the global structure of BPTI result on a reduction of the radius of gyration from 13.4 A down to 12.0 A. Between 5,000 and 6,000 bar, some transition already detected by FTIR [N. Takeda, K. Nakano, M. Kato, Y. Taniguchi, Biospectroscopy, 4, 1998, pp. 209-216] is observed. The pressure effect is not reversible because the initial value of the radius of gyration is not recovered after pressure release. By extending the range of wave-vectors to high q, we have observed a change of the form factor (shape) of the BPTI under pressure. At atmospheric pressure BPTI exhibits an ellipsoidal form factor that is characteristic of the native state. When the pressure is increased from atmospheric pressure up to 6,000 bar, the protein keeps its ellipsoidal shape. The parameters of the ellipsoid vary and the transition detected between 5,000 and 6,000 bar in the form factor of BPTI is in agreement with the FTIR results. After pressure release, the form factor of BPTI is characteristic of an ellipsoid of revolution with a semi-axis a, slightly elongated with respect to that of the native one, indicating that the pressure-induced structural changes on the protein are not reversible. The global motions and the internal dynamics of BPTI protein have been investigated in the same pressure range by quasi-elastic neutron scattering experiments on IN5 time-of-flight spectrometer at ILL (Grenoble, France). The diffusion coefficients D and the internal relaxation times of BPTI deduced from the analysis of the intermediate scattering functions show a slowing down of protein dynamics when increasing pressure.     

Collective chain dynamics in lipid bilayers by inelastic x-ray scattering

We investigated the application of inelastic x-ray scattering (IXS) to lipid bilayers. This technique directly measures the dynamic structure factor S(q,omega) which is the space-time Fourier transform of the electron density correlation function of the measured system. For a multiatomic system, the analysis of S(q,omega) is usually complicated. But for multiple bilayers of lipid, S(q,omega) is dominated by chain-chain correlations within individual bilayers. Thus IXS provides a unique probe for the collective dynamics of lipid chains in a bilayer that cannot be obtained by any other method. IXS of dimyristoyl phosphatidylcholine and dimyristoyl phosphatidylcholine + cholesterol at two different concentrations were measured. S(q,omega) was analyzed by three-mode hydrodynamic equations, including a thermal diffusive mode and two propagating acoustic modes. We obtained the dispersion curves for the phonons that represent the collective in-plane excitations of lipid chains. The effect of cholesterol on chain dynamics was detected. Our analysis shows the importance of having a high instrument resolution as well as the requirement of sufficient signal-to-noise ratio to obtain meaningful results from such an IXS experiment. The requirement on signal-to-noise also applies to molecular dynamics simulations.     

The BPTI decamer observed in acidic pH crystal forms pre-exists as a stable species in solution

Bovine pancreatic trypsin inhibitor (BPTI) crystallizes under acidic pH conditions in the presence of thiocyanate, chloride and sulfate ions, yielding three different polymorphs in P2(1), P6(4)22 and P6(3)22 space groups, respectively. In all three crystal forms, the same decamer is found in the packing (ten BPTI molecules organized through two perpendicular 2-fold and 5-fold axes as a well-defined and compact object) in contrast to the monomeric crystal forms observed at basic pH conditions. The crystallization of BPTI under acidic conditions (pH 4.5) was investigated by small angle X-ray scattering with both under- and supersaturated BPTI solutions. Data showed the oligomerization of BPTI molecules under all investigated conditions. Accordingly, various mixtures of discrete oligomers (n=1 to 10) were considered. Calculated scattering curves were obtained using models based on the crystallographic structures, and the experimental patterns were analyzed as a linear combination of the model curves using a non-linear curve fitting procedure. The results, confirmed by gel filtration experiments, unambiguously demonstrate the co-existence of two different BPTI particles in solution: a monomer and a decamer, with no evidence of any other intermediates. Moreover, using both approaches, the fraction of decamers was found to increase with increasing salt concentration, even beyond the solubility curve. We therefore propose that at acidic pH, BPTI crystallizes following a two step process: decamers are first built in under- and supersaturated solutions, upon which crystal growth proceeds by decamer stacking. Indeed, those BPTI crystals should best be described as "BPTI decamer" crystals.     

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.     

Direct Measurement of Hydration-Related Dynamic Changes in Lysozyme using Inelastic Neutron Scattering Spectroscopy

Abstract

Inelastic neutron scattering spectroscopy is used to investigate dynamic changes in lysozyme powder at two different low D,0 hydrations (0.07g D²,O/g protein and 0,20 g D²,O/g protein). In the higher hydration sample, the inelastic scattering between 0.8 and 4.0 cm^?1 energy transfer is increased and the elastic scattering is decreased. The decreased elastic scattering suggests increased atomic amplitudes of motion and the increased 0.8 to 4.0 cm^?1 scattering suggests increased motions in this frequency range. Comparison with normal mode models of lysozyme dynamics shows that the inelastic difference occurs in the frequency region predicted for the lowest frequency, largest amplitude, global modes of the molccule[M. Levitt, C. Sanderand P. S. Stern, J. Mol. Biol. 181. 423 (1985). B.Brooks and M.Karplus.Prot. Natl Acad. Sci (U.S.A) 82. 4995 (1985), R.E. Bruccoleri, M. Karplus and J.A. McCammon, Biopolymers 25 1767 (1986)]. Our results are consistent with a model in which an increased number of low frequency global modes are present in the higher hydrated sample.     

Iron normal mode dynamics in (nitrosyl)iron(II)tetraphenylporphyrin from X-ray nuclear resonance data

The complete iron atom vibrational spectrum has been obtained by refinement of normal mode calculations to nuclear inelastic x-ray absorption data from (nitrosyl)iron(II)tetraphenylporphyrin, FeTPP(NO), a useful model for heme dynamics in myoglobin and other heme proteins. Nuclear resonance vibrational spectroscopy (NRVS) provides a direct measurement of the frequency and iron amplitude for all normal modes involving significant displacement of (57)Fe. The NRVS measurements on isotopically enriched single crystals permit determination of heme in-plane and out-of-plane modes. Excellent agreement between the calculated and experimental values of frequency and iron amplitude for each mode is achieved by a force-field refinement. Significantly, we find that the presence of the phenyl groups and the NO ligand leads to substantial mixing of the porphyrin core modes. This first picture of the entire iron vibrational density of states for a porphyrin compound provides an improved model for the role of iron atom dynamics in the biological functioning of heme proteins.     

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.     

Coupling of laser excitation and inelastic neutron scattering: attempt to probe the dynamics of light-induced C-phycocyanin dynamics

Excitation energy transfer (EET) in light-harvesting antennae is a highly efficient key event in photosynthesis, where light-induced dynamics of the antenna pigment-protein complexes may play a functional role. So far, however, the relationship between EET and protein dynamics remains unknown. C-phycocyanin (C-PC) is the main pigment/protein complex present in the cyanobacterial antenna, called "phycobilisome". The aim of the present study was to investigate light-induced C-PC internal thermal motions (ps timescale) measured by inelastic neutron scattering. To synchronize the beginning of the laser flash (6 ns duration) with that of the neutron test pulse ( approximately 87 mus duration), we developed a novel type of "time-resolved" experimental setup on MIBEMOL time-of-flight neutron spectrometer (LLB, France). Data acquisition has been modified to get quasi-simultaneously "light" and "dark" measurements (with and without laser, respectively) and eliminate many spurious effects that could occur on the sample during the experiment. The study was carried out on concentrated C-PC ( approximately 135 g/L protein in D(2)O phosphate buffer), contained in an aluminium/sapphire sample holder (almost "transparent" for neutrons) and homogeneously illuminated inside an "integrating sphere". We observed very similar incoherent dynamical structure factors of C-PC with or without light. The vibrational density of states showed two very slightly increased vibrational modes with light, at approximately 30 and approximately 50 meV ( approximately 240 and approximately 400 cm(-1), respectively). These effects have to be verified by further experiments before probing any temporal evolution, by introducing a time delay between the laser flash and the neutron test pulse.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

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.     

Temperature dependence of dynamics of hydrated myoglobin. Comparison of force field calculations with neutron scattering data

Molecular dynamics is used to probe the atomic motions of the carboxy-myoglobin protein as a function of temperature. Simulations of 150 picoseconds in length are carried out on the protein at 20, 60, 100, 180, 220, 240, 260, 280, 300, 320 and 340 K. The simulations attempt to mimic neutron scattering experiments very closely by including a partial hydration shell around the protein. Theoretical elastic, quasielastic and inelastic neutron scattering data are derived from the trajectories and directly compared with experiment. Compared to experiment, the simulation-derived elastic scattering curves show a decrease in intensity as a function of the scattering wavevector, q2. The inelastic and quasielastic spectra show that the inelastic peak is shifted to lower frequency than the experimental value, while quasielastic behavior is in good agreement with experiment. This suggests that the theoretical model is too flexible in the harmonic limit (low temperature), but accurately reproduces high-temperature behavior. Time correlation functions of the intermediate scattering function are determined. At low temperature there is one fast decay process, and at high temperatures there is an additional slow relaxation process that is due to quasielastic scattering. The average atomic fluctuations show that the protein behaves harmonically at low temperatures. At approximately 210 K, a glass-like transition in atomic fluctuations is seen. Above the transition temperature, the atomic fluctuations exhibit both harmonic and anharmonic behavior. Comparison of protein mobility behavior with experiment indicate the fluctuations derived from simulations are larger in the harmonic region. However, the anharmonic region agrees very well with experiment. The anharmonicity is large at all temperatures, with a gradual monotonic increase from 0.5 at 20 K to greater than 0.7 at 340 K without a noticeable change at the glass transition temperature. Heavy-atom dihedral transitions are monitored as a function of temperature. Trends in the type of dihedral transitions that occur with temperature are clearly visible. Dihedral transitions involving backbone atoms occur only above the glass transition temperature. The overall protein behavior results suggest that at low temperatures there is purely vibrational motion with one fast decay process, and above the glass transition temperature there is more anharmonic motion with a fast and a slower relaxation process occurring simultaneously.(ABSTRACT TRUNCATED AT 400 WORDS)