Dynamics of hydration water in deuterated purple membranes explored by neutron scattering

The function and dynamics of proteins depend on their direct environment, and much evidence has pointed to a strong coupling between water and protein motions. Recently however, neutron scattering measurements on deuterated and natural-abundance purple membrane (PM), hydrated in H(2)O and D(2)O, respectively, revealed that membrane and water motions on the ns-ps time scale are not directly coupled below 260 K (Wood et al. in Proc Natl Acad Sci USA 104:18049-18054, 2007). In the initial study, samples with a high level of hydration were measured. Here, we have measured the dynamics of PM and water separately, at a low-hydration level corresponding to the first layer of hydration water only. As in the case of the higher hydration samples previously studied, the dynamics of PM and water display different temperature dependencies, with a transition in the hydration water at 200 K not triggering a transition in the membrane at the same temperature. Furthermore, neutron diffraction experiments were carried out to monitor the lamellar spacing of a flash-cooled deuterated PM stack hydrated in H(2)O as a function of temperature. At 200 K, a sudden decrease in lamellar spacing indicated the onset of long-range translational water diffusion in the second hydration layer as has already been observed on flash-cooled natural-abundance PM stacks hydrated in D(2)O (Weik et al. in J Mol Biol 275:632-634, 2005), excluding thus a notable isotope effect. Our results reinforce the notion that membrane-protein dynamics may be less strongly coupled to hydration water motions than the dynamics of soluble proteins.     

Protein dynamics in solution and powder measured by incoherent elastic neutron scattering: the influence of Q-range and energy resolution

Incoherent elastic neutron scattering (IENS) has been widely used to measure intramolecular atomic mean square displacements (MSDs) of proteins in powder and in solution. The instrumental energy resolution and the wave vector transfer (Q-range) determine, respectively, the time and length scales of observable motions. In order to investigate contributions of diffusive motions to MSDs measured by this method, we calculated the elastic intensity for several simple scattering functions. We showed that continuous translational diffusion contributes to MSDs in a Q-range where the energy width of the scattering function is of the order of the instrumental energy resolution. We discuss the choice of instruments adapted to focus on intramolecular motions in the presence of solvent or global macromolecular diffusion. The concepts developed are applied to interpret experimental data from H₂O- and D₂O-hydrated proteins. Finally, analogies between the Gaussian approximation in IENS and the Guinier approximation in small-angle scattering are discussed.     

Incoherent elastic and quasi-elastic neutron scattering investigation of hemoglobin dynamics

In this work we investigate the dynamic properties of hemoglobin in glycerolD(8)/D(2)O solution using incoherent elastic (ENS) and quasi-elastic (QENS) neutron scattering. Taking advantage of complementary energy resolutions of backscattering spectrometers at ILL (Grenoble), we explore motions in a large space-time window, up to 1 ns and 14 A; moreover, in order to cover the harmonic and anharmonic protein dynamics regimes, the elastic experiments have been performed over the wide temperature interval of 20-300 K. To study the dependence of the measured dynamics upon the protein quaternary structure, both deoxyhemoglobin (in T quaternary conformation) and carbonmonoxyhemoglobin (in R quaternary conformation) have been investigated. From the ENS data the mean square displacements of the non-exchangeable hydrogen atoms of the protein and their temperature dependence are obtained. In agreement with previous results on hydrated powders, a dynamical transition at about 220 K is detected. The results show interesting differences between the two hemoglobin quaternary conformations, the T-state protein appearing more rigid and performing faster motions than the R-state one; however, these differences involve motions occurring in the nanosecond time scale and are not detected when only faster atomic motions in the time scale up to 100 ps are investigated. The QENS results put in evidence a relevant Lorentzian quasi-elastic contribution. Analysis of the dependence of the Elastic Incoherent Structure Factor (EISF) and of the Lorentzian halfwidth upon the momentum transfer suggests that the above quasi-elastic contribution arises from the diffusion inside a confined space, values of confinement radius and local diffusion coefficient being compatible with motions of hydrogen atoms of the amino acid side chains. When averaged over the whole range of momentum transfer the QENS data put in evidence differences between deoxy and carbonmonoxy hemoglobin and confirm the quaternary structure dependence of the protein dynamics in the nanosecond time scale.     

Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from M?ssbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.     

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)     

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.     

Neutron Laue macromolecular crystallography

Recent progress in neutron protein crystallography such as the use of the Laue technique and improved neutron optics and detector technologies have dramatically improved the speed and precision with which neutron protein structures can now be determined. These studies are providing unique and complementary insights on hydrogen and hydration in protein crystal structures that are not available from X-ray structures alone. Parallel improvements in modern molecular biology now allow fully (per)deuterated protein samples to be produced for neutron scattering that essentially eradicate the large-and ultimately limiting-hydrogen incoherent scattering background that has hampered such studies in the past. High quality neutron data can now be collected to near atomic resolution (approximately 2.0 A) for proteins of up to approximately 50 kDa molecular weight using crystals of volume approximately 0.1 mm3 on the Laue diffractometer at ILL. The ability to flash-cool and collect high resolution neutron data from protein crystals at cryogenic temperature (15 K) has opened the way for kinetic crystallography on freeze trapped systems. Current instrument developments now promise to reduce crystal volume requirements by a further order of magnitude, making neutron protein crystallography a more accessible and routine technique.     

C-phycocyanin hydration water dynamics in the presence of trehalose: an incoherent elastic neutron scattering study at different energy resolutions

We present a study of C-phycocyanin hydration water dynamics in the presence of trehalose by incoherent elastic neutron scattering. By combining data from two backscattering spectrometers with a 10-fold difference in energy resolution we extract a scattering law S(Q,omega) from the Q-dependence of the elastic intensities without sampling the quasielastic range. The hydration water is described by two dynamically different populations--one diffusing inside a sphere and the other diffusing quasifreely--with a population ratio that depends on temperature. The scattering law derived describes the experimental data from both instruments excellently over a large temperature range (235-320 K). The effective diffusion coefficient extracted is reduced by a factor of 10-15 with respect to bulk water at corresponding temperatures. Our approach demonstrates the benefits and the efficiency of using different energy resolutions in incoherent elastic neutron scattering over a large angular range for the study of biological macromolecules and hydration water.     

Harmonic behavior of trehalose-coated carbon-monoxy-myoglobin at high temperature.

Embedding biostructures in saccharide glasses protects them against extreme dehydration and/or exposure to very high temperature. Among the saccharides, trehalose appears to be the most effective bioprotectant. In this paper we report on the low-frequency dynamics of carbon monoxy myoglobin in an extremely dry trehalose glass measured by neutron spectroscopy. Under these conditions, the mean square displacements and the density of state function are those of a harmonic solid, up to room temperature, in contrast to D2O-hydrated myoglobin, in which a dynamical transition to a nonharmonic regime has been observed at approximately 180 K (Doster et al., 1989. Nature. 337:754-756). The protective effect of trehalose is correlated, therefore, with a trapping of the protein in a harmonic potential, even at relatively high temperature.     

Molecular dynamics decomposition of temperature-dependent elastic neutron scattering by a protein solution

Molecular dynamics simulations are performed of bovine pancreatic trypsin inhibitor in a cryosolution over a range of temperatures from 80 to 300 K and the origins identified of elastic dynamic neutron scattering from the solution. The elastic scattering and mean-square displacement calculated from the molecular dynamics trajectories are in reasonable agreement with experiments on a larger protein in the same solvent. The solvent and protein contributions to the scattering from the simulation model are determined. At lower temperatures (< approximately 200 K) or on shorter timescales ( approximately 10 ps) the scattering contributions are proportional to the isotopic nuclear scattering cross-sections of each component. However, for T > 200 K marked deviations from these cross-sections are seen due to differences in the dynamics of the components of the solution. Rapid activation of solvent diffusion leads to the variation with temperature of the total elastic intensity being determined largely by that of the solvent. At higher temperatures (>240 K) and longer times ( approximately 100 ps) the protein makes the only significant contribution to the scattering, the solvent scattering having moved out of the accessible time-space window. Decomposition of the protein mean-square displacement shows that the observed dynamical transition in the solution at 200-220 K involves activation of both internal motions and external whole-molecule rotational and translational diffusion. The proportion that the external dynamics contributes to the protein mean-square displacement increases to approximately 30 and 60% at 300 K on the 10- and 100-ps timescales, respectively.     

Evolution of the internal dynamics of two globular proteins from dry powder to solution.

Myoglobin and lysozyme picosecond internal dynamics in solution is compared to that in hydrated powders by quasielastic incoherent neutron scattering. This technique is sensitive to the motions of the nonexchangeable hydrogen atoms in a sample. Because these are homogeneously distributed throughout the protein structure, the average dynamics of the protein is described. We first propose an original data treatment to deal with the protein global motions in the case of solution samples. The validity of this treatment is checked by comparison with classical measurements of the diffusion constants. The evolution with the scattering vector of the width and relative contribution of the quasielastic component was then used to derive information on the amount of local diffusive motions and their characteristic average relaxation time. From dry powder to coverage by one water layer, the surface side chains progressively acquire the possibility to diffuse locally. On subsequent hydration, the main effect of water is to improve the rate of these diffusive motions. Motions with higher average amplitude occur in solution, about three times more than for a hydrated powder at complete coverage, with a shorter average relaxation time, approximately 4.5 ps compared to 9.4 ps for one water monolayer.     

Temperature- and hydration-dependent protein dynamics in photosystem II of green plants studied by quasielastic neutron scattering

Protein dynamics in hydrated and vacuum-dried photosystem II (PS II) membrane fragments from spinach has been investigated by quasielastic neutron scattering (QENS) in the temperature range between 5 and 300 K. Three distinct temperature ranges can be clearly distinguished by active type(s) of protein dynamics: (A) At low temperatures (T < 120 K), the protein dynamics of both dry and hydrated PS II is characterized by harmonic vibrational motions. (B) In the intermediate temperature range (120 < T < 240 K), the total mean square displacement total slightly deviates from the predicted linear behavior. The QENS data indicate that this deviation, which is virtually independent of the extent of hydration, is due to a partial onset of diffusive protein motions. (C) At temperatures above 240 K, the protein flexibility drastically changes because of the onset of diffusive (large-amplitude) protein motions. This dynamical transition is clearly hydration-dependent since it is strongly suppressed in dry PS II. The thermally activated onset of protein flexibility as monitored by QENS is found to be strictly correlated with the temperature-dependent increase of the electron transport efficiency from Q(A)(-) to QB (Garbers et al. (1998) Biochemistry 37, 11399-11404). Analogously, the freezing of protein mobility by dehydration in dry PS II appears to be responsible for the blockage of Q(A)(-) reoxidation by Q(B) at hydration values lower than 45% r.h. (Kaminskaya et al. (2003) Biochemistry 42, 8119-8132). Similar effects were observed for reactions of the water-oxidizing complex as outlined in the Discussion section.     

Picosecond molecular motions in bacteriorhodopsin from neutron scattering.

The characteristics of internal molecular motions of bacteriorhodopsin in the purple membrane have been studied by quasielastic incoherent neutron scattering. Because of the quasihomogeneous distribution of hydrogen atoms in biological molecules, this technique enables one to study a wide variety of intramolecular motions, especially those occurring in the picosecond to nanosecond time scale. We performed measurements at different energy resolutions with samples at various hydration levels within a temperature range of 10-300 K. The analysis of the data revealed a dynamical transition at temperatures Td between 180 K and 220 K for all motions resolved at time scales ranging from 0.1 to a few hundred picoseconds. Whereas below Td the motions are purely vibrational, they are predominantly diffusive above Td, characterized by an enormously broad distribution of correlation times. The variation of the hydration level, on the other hand, mainly affects motions slower than a few picoseconds.     

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.     

Thermal fluctuations of DNA enclosed by glycerol–water glassy matrices: an elastic neutron scattering investigation

Through elastic neutron scattering measurements, we investigated the thermal fluctuations of DNA enclosed by glycerol–water glassy matrices, at different levels of hydration, over the wide temperature range from 20 to 300 K. For all the samples, the extracted hydrogen mean square displacements (MSD) show a purely vibrational harmonic trend at very low temperatures, and a first onset of anharmonic dynamics above ∼100 K. Such onset is consistent with the activation of DNA methyl group rotational motions. Then, at a certain temperature T _d, the MSD show a second onset of anharmonicity, which corresponds to the DNA dynamical transition. The T _d values vary as a function of the hydration degree of the environment. The crucial role of the solvent mobility to activate the DNA thermal fluctuations is proposed, together with a preferential hydration effect of the DNA phosphate groups. Finally, a comparison between the average mobility of homologous samples of DNA and the lysozyme protein is considered. Advanced neutron scattering and complementary techniques to study biological systems. Contributions from the meetings, “Neutrons in Biology”, STFC Rutherford Appleton Laboratory, Didcot, UK, 11–13 July and “Proteins At Work 2007”, Perugia, Italy, 28–30 May 2007.     

Selective deuteration of tryptophan and methionine residues in maltose binding protein: a model system for neutron scattering

We describe methods that have been developed within the ILL-EMBL Deuteration Laboratory for the production of maltose binding protein (MBP) that has been selectively labelled either with deuterated tryptophan or deuterated methionine (single labelling), or both (double labelling). MBP is used as an important model system for biophysical studies, and selective labelling can be helpful in the analysis of small-angle neutron scattering (SANS) data, neutron reflection (NR) data, and high-resolution neutron diffraction data. The selective labelling was carried out in E. coli high-cell density cultures using auxotrophic mutants in minimal medium containing the required deuterated precursors. Five types of sample were prepared and studied: (1) unmodified hydrogenated MBP (H-MBP), (2) perdeuterated MBP (D-MBP), (3) singly labelled MBP with the tryptophan residues deuterated (D-trp MBP), (4) singly labelled MBP with methionine residues deuterated (D-met MBP) and (5) doubly labelled MBP with both tryptophan and methionine residues deuterated (D-trp/met MBP). Labelled samples were characterised by size exclusion chromatography, gel electrophoresis, light scattering and mass spectroscopy. Preliminary small-angle neutron scattering (SANS) experiments have also been carried out and show measurable differences between the SANS data recorded for the various labelled analogues. More detailed SANS experiments using these labelled MBP analogues are planned; the degree to which such data could enhance structure determination by SANS is discussed.     

Combination of Neutron Scattering and Molecular Dynamics to Determine Internal Motions in Biomolecules

Abstract

The forms and frequencies of atomic dynamics on the pico- and nanosecond timescales are accessible experimentally using incoherent neutron scattering. Molecular dynamics simulations cover the same space and time domains and neutron scattering intensities can be calculated from the simulations for direct comparison with experiment. To illustrate the complementarity of neutron scattering and molecular dynamics we examine measured and simulation-derived elastic incoherent scattering profiles from myoglobin and from the crystalline alanine dipeptide. Elastic incoherent scattering gives information on the geometry of the volume accessible to the atoms in the samples. The simulation-derived dipeptide elastic scattering profiles are in reasonable accord with experiment, deviations being due to the sampling limitations in the simulations and experimental detector normalisation procedures. The simulated dynamics is decomposed, revealing characteristic profiles due to rotational diffusional and translational vibrational motions of the methyl groups. In myoglobin, for which the timescale of the simulation matches more closely that accessible to the experiment, good agreement is seen for the elastic incoherent structure factor. This indicates that the space sampled by the hydrogen atoms in the protein on the timescale <100 ps is well represented by the simulation. Part of the helix atom fluctuations can be described in terms of rigid helix motions.     

Frozen spin targets in ribosomal structure research.

Polarized neutron scattering strongly depends on nuclear spin polarisation, particularly on proton spin polarisation. A single proton in a deuterated environment then is as efficient as 10 electrons in X-ray anomalous diffraction. Neutron scattering from the nuclear spin label is controlled by the polarisation of neutron spins and nuclear spins. Pure deuteron spin labels and proton spin labels are created by NMR saturation. We report on results obtained from the large subunit of E. coli ribosomes which have been obtained at the research reactor of GKSS using the polarized target facility developed by CERN. The nuclear spins were oriented with respect to an external field by dynamic nuclear polarisation. Proton spin polarisations of more than 80% were obtained in ribosomes at temperatures below 0.5 K. At T = 130 mK the relaxation time of the polarized target is one month (frozen spin target). Polarized small-angle neutron scattering of the in situ structure of rRNA and the total ribosomal protein (TP) has been determined from the frozen spin targets of the large ribosomal subunit, which has been deuterated in the TP and rRNA respectively. The results agree with those from neutron scattering in H2O/D2O mixtures obtained at room temperature. This is a necessary prerequisite for the planned determination of the in situ structure of individual ribosomal proteins and especially of that of ribosome bound mRNA and tRNAs.     

Measurement of coherent Debye-Waller factor in in vivo deuterated C-phycocyanin by inelastic neutron scattering.

Quasielastic neutron scattering measurements of dry and 35% D2O hydrated amorphous protein powder of C-phycocyanin were made as a function of temperature ranging from 313K down to 100K. The protein is grown from blue-green algae cultured in D2O and is deuterated up to 99%. The scattering is thus dominated by coherent scattering. Within the best energy resolution of the time-of-flight instrument, which is 28 mueV FWHM, the scattering appears entirely elastic. For this reason we are able to extract a coherent Debye-Waller factor by making an independent measurement of the static structure factor. We observe a considerable difference in the q dependence of the Debye-Waller factor between the dry and hydrated proteins; furthermore, there is an interesting temperature dependence of the Debye-Waller factor that is quite different from that predicted for dense hard-sphere liquids.     

Dynamics of Protein and its Hydration Water: Neutron Scattering Studies on Fully Deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<∼200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <∼180–190 K and is not affected by the proteins’ methyl group rotations. The dynamics of the hydration water exhibits changes at ∼180–190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <∼3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.