Neutron Scattering in Biology - Fitter Gutberlet and Katsaras
.pdf16 Quasielastic Neutron Scattering: Applications |
377 |
of the former is appreciably larger than that of the latter. This may be due to one or several of various conceivable reasons. Let us first consider the possibility, that the additional quasielastic scattering may be due to the mobility of small side-groups located on the ligand, which might be higher than that of side-groups in the main part of the complex. Because of the low weight fraction of the ligand, this is not considered as very probable, although it cannot be completely excluded. Further studies employing a deuterated ligand could however clarify this question. Second, multiphonon scattering can in principle lead to a quasielastic contribution due to phonon–phonon annihilation. This e ect should, however, be negligible in the concerned Q-range (below 2.2 ˚A−1), and especially at the low temperature of the experiment, where it is estimated to be of the order of only 2% of the one-phonon scattering. Furthermore such a component is not likely to be very di erent for the two kinds of sample. Similar arguments apply for a multiple scattering contribution, which again is negligibly small, because of the large value of the employed sample transmissions (97.7%).
Finally, we come to what is believed to be the main reason for the observed di erence in the spectra. The complexation may lead to an increased importance of relaxational modes and possibly of very-low-frequency optical vibrations in this whole system and thus to an increased flexibility of the macromolecular ensemble, as compared to the uncomplexed case. The intensity increase also extends further into the inelastic region up to 40 cm−1. It is likely, but cannot be proved at present, that the latter e ect might be due to complexation-caused damping of the vibrations in the latter intermediate frequency region. Damping could lead to a shift of modes from higher to lower frequencies. The deviation from purely harmonic behavior of both systems, (DHFR+NADPH) and (DHFR+NADPH+MTX), is demonstrated by calculating, in the limit Q → 0, a “generalized” density of states g(ω) using Eq. 16.17 in Sect. 16.3.3. This is shown in Fig. 16.12 presenting the g(ω) curves for the complexed and for the uncomplexed protein. 7 It is seen, that in the frequency region below 20 cm−1 the obtained generalized function does not show the well-known ω2-behavior expected for purely harmonic acoustic vibrations at low energies. Furthermore, it seems to be clear that the curves start with finite values at ω = 0, although measured data are of course not available at zero frequency (see Sect.16.3.3). This is of course the phenomenological signature of the quasielastic component. It is also evident from the figure, that the complexed form of DHFR has an appreciably larger g(ω) in
7Note that the low-frequency part (ω below 20 cm−1) of this function must be considered as purely phenomenological: because stochastic fluctuations generally do not obey Bose–Einstein statistics, the division of the scattering function by the Bose phonon-population factor does not have quite the same justification as in the case of harmonic vibrations. But this does not change the main result, namely the observation that complexation of DHFR leads to additional relaxational phenomena causing an increase in the quasielastic scattering component.
378 R.E. Lechner et al.
|
8 |
|
|
|
|
|
|
1.0 |
|
|
|
|
|
) |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
1− |
|
|
|
|
|
|
cm |
|
|
|
|
|
|
(mode |
0 |
40 |
80 |
120 |
|
|
0 |
|
|
||||
0.5 |
|
|
|
|
|
|
g(ω) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.0 |
|
|
|
|
|
|
0 |
5 |
|
10 |
15 |
20 |
|
|
|
|
w (cm−1) |
|
|
Fig. 16.12. Generalized density of states g(ω) for uncomplexed DHFR (lower curve) and for DHFR complexed with methotrexate (upper curve) at 120 K. Data, derived from measurements with the direct-geometry TOF spectrometer IN6 (see Sect. 16.3.1 in Part I, this volume), have been summed over a large range of scattering angles. Inset: enlarged frequency region of the spectra; after [93]
this low-frequency region, suggesting that complexation of DHFR softens the enzyme and makes this macromolecule more flexible, an e ect which could be relevant for its biochemical activity. However, we have to note, that we do not know, whether this behavior is also relevant to the same extent at physiological conditions, where the level of hydration established initially (see above) is e ective.
Finally, it is interesting to note that, using NMR relaxation experiments, other authors have observed an increased conformational flexibility on binding a hydrophobic ligand to mouse major urinary protein [104]. However, we should note that a flexibility increase does not necessarily occur as a rule upon ligand binding, since this depends on the strength of the involved interaction. In other investigations using NMR and crystallography, a flexibility decrease was observed on binding small organic ligands to proteins [105–107].
16.5 Low-Dimensional Systems
16.5.1 Two-Dimensional Long-Range Di usion
of Rotating Molecules
Low-dimensional di usion plays an important role for certain biological objects, such as ion transport channels and membrane surfaces. Obviously, the concept of low-dimensionality should not be understood for instance in the sense of strictly planar motions, since especially on biological membranes there will be appreciable local deviations from this idealized picture. The essential
16 Quasielastic Neutron Scattering: Applications |
379 |
condition is, that the di using particles stay within or near a fictive di usion plane (e.g., representing the membrane surface) at least during the observation time defined by the energy resolution of the experiment (see Sect. 15.3, Eq. 15.30 in Part I, this volume and the related discussion). This also implies that the mean-square displacement parallel to the di usion layer, the particle has acquired during this time, is much larger than perpendicular to it.
The incoherent neutron scattering function in the low-Q limit for longrange translational di usion (TD) was already considered above (see Eqs. 16.6 and 16.7. Accordingly, in three dimensions, when the di usive motion is isotropic, with a di usion coe cient D3D, we have for this function:
S3D |
(Q, ω) = |
1 |
|
D3DQ2 |
, |
(16.30) |
|
π (D3DQ2)2 + ω2 |
|||||||
TD |
|
|
|
||||
|
|
|
|
Certain analytic results [108] concerning the anisotropic case, where the di u- sion process is restricted to the surface of a plane or to the planes of a layered structure, will be discussed in the following. If single-crystalline samples are available, the dimensionality of the di usional motion can of course be determined directly by studying the orientation dependence in a QINS experiment. For a single crystal we have [109]:
S2D |
(Q, ω) = |
1 |
|
D2D(Q sin θ)2 |
, |
(16.31) |
|
π (D2D(Q sin θ)2)2 + ω2 |
|||||||
TD |
|
|
|
||||
|
|
|
|
where D2D is the coe cient of self-di usion in two dimensions, Q sin θ the component of the scattering vector in the di usion plane, and θ the angle between Q and the normal to this plane. If single-crystals are not available, one has to resort to polycrystalline samples requiring orientational averaging of the above expression. It is known, that the resulting integral over all orientations exhibits a logarithmic singularity at zero energy transfer [109]. This is caused by the fact that di usion planes which are perpendicular (or close to perpendicular) to the scattering vector Q contribute elastic (or almost elastic) scattering to the QINS function. Fortunately, the logarithmic singularity is cancelled by finite resolution. It has been shown that for resolution functions which have the shape of a Lorentzian (with HWHM equal to H) or of a sum of Lorentzians (which can be the case in BSC experiments) the orientationally averaged resolution-broadened QINS function for 2D-di usion,
|
|
|
|
|
1 |
0 |
2π |
1 |
0 |
π |
|
||||
< STD2D |
(Q, ω) >orientres |
.= |
dΦ |
sin θ dθ L(θ, ω), |
(16.32) |
||||||||||
|
|
|
|
||||||||||||
2π |
2 |
||||||||||||||
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L(θ, ω) = |
|
1 |
|
|
D2D(Q sin θ)2 + H |
|
(16.33) |
|||||||
|
π [D2D |
(Q sin θ)2 |
+ H]2 + ω2 |
||||||||||||
|
|
|
|
can be written in closed form [108]. The explicit expression is lengthy and will not be given here. The determination of such a characteristic line shape is very
380 R.E. Lechner et al.
di cult. It requires an experimental resolution width H D2DQ2, in order to be able to distinguish it from a simple Lorentzian shape. This problem is discussed in detail in [108].
The study of the anisotropy of translational di usion in strongly anisotropic material, such as biological membranes, is straight-forward, if the membranes are well aligned in the samples. Pertinent investigations with quasielastic neutron scattering techniques have been carried out on purple membrane [110,111] and on the superficial layer of porcine skin (stratum corneum) [112]. As an example we discuss the transport of water molecules on the surface of purple membrane which has a two-dimensional crystalline structure. In biological membranes, proton di usion connected with conduction of protons provides an important mechanism of energy transduction in living organisms. The purple membrane (PM) of Halobacterium salinarum, for instance, contains the protein bacteriorhodopsin (BR) which becomes a one-dimensional stochastically pulsed proton conductor, when activated by light (see Sect. 16.5.2). This is a light-driven proton pump generating an electrochemical gradient across the membrane, which is employed by the bacterium as an energy source. After having been pumped from the cell-interior to the membrane surface, the protons are transported by a mechanism of surface di usion towards other proteins located within the same membrane: The light-generated electrochemical potential across the cell membrane is utilised by the halobacteria to furnish the driving force for ATP synthesis by energizing the rotation of the turbinelike machinery in ATP synthase. Water molecules near the surface are known to be relevant for this biological function, since they have been shown to assist the proton conductivity [113]. Hydration water, its interaction with the surface of biological macromolecules and macromolecular complexes, and its di usion generally play an important role in structure, dynamics, and function of biological systems [114, 115]. It is therefore of considerable interest to study the proton di usion mechanism within and close to the hydration layers of membranes. Similar to the case of bulk water, it is expected that during the di usion process protons are exchanged between water molecules acting either as acceptors (forming for instance (H3O)+) or as donors (producing (OH)− ions). A “solid-like” Grotthuss feature [116–118] is added to the di usion of protons in the liquid water phase by the presence of fixed protonation sites on the surface of purple membrane. It is worthwhile to note, that these sites are arranged in space in a perfectly regular manner, since the bacteria use the most e cient packing of BR: trimers of BR molecules embedded in a lipid bilayer matrix are aligned in a two-dimensional hexagonal single-crystalline structure [119]. However, the pH-value not being very different from 7, the concentration of charged particles is very low. The protons spend most of the time as part of di using neutral water molecules with only rare events of exchange between di erent “vehicles” [116]. Therefore, for the purpose of analyzing neutron scattering experiments, the whole mechanism of di usion may to a good approximation be classified as that of molecular di usion [116].
384 R.E. Lechner et al.
photocycle and the actual vectorial proton transport activity of the proton pump are drastically modified [135–137].
It is interesting to compare the hydration dependence of the photocycle to its variation observed as a function of temperature. While the quantum e ciency of BR is practically temperature independent, the decay of all intermediates relying on thermal activation obviously is not. The photocycle slows down with decreasing temperature and is “frozen in” at some intermediate state, when the time constant of the latter has become practically infinite due to the low value of the temperature. It was found, for instance, that bacteriorhodopsin can be “captured” in the M2-intermediate state at 260 K [138]. Furthermore, according to [139, 140], the photocycle stops at the M-intermediate state below T 220 K (at the M1-intermediate state near 230 K, after [141]), at the L-intermediate below T 180 K (near 155 K, after [141]) and at the K-intermediate below T 150 K (near 90 K, after [141]). The N and O intermediates were not observed at all by low-temperature spectrophotometry. Below we discuss the question, whether the T - and h-dependences of the photocycle are significantly correlated with the T - and h-dependent e ects observed in the static and dynamic structure of multilayered stacks of PM.
In a previous study [100] we have investigated the static structure of multilayered stacks of PM as a function of temperature (T ), after having established well-defined levels of hydration (h) by equilibration at fixed relative humidities (r.h.), at room temperature. The lamellar spacing dl of these systems was measured with neutron di raction as a function of T and h. The observed large T -dependent variations of dl indicate that PM is partially dehydrated, when cooled below a “hydration water freezing point.” This e ect is reversible, but a hysteresis is observed, when PM is rehydrated upon reheating. This phenomenon of dehydration and rehydration, induced by cooling and reheating, respectively, appears to be a general property of biological membranes (see also [142]). It is caused by the presence of hydration forces and by the specific (di erent) temperature dependences of the chemical potentials of interbilayer water and ice. These forces result in a local freezing point depression and are the main reason, why nucleation of ice crystals first occurs outside of the space between bilayers. The temperature variation of the chemical potentials leads to a di erence in vapor pressure, causing water to be successively extracted from the interbilayer space in the presence of these ice crystals, when the temperature is lowered. Vice versa upon heating, the water originating from the melting of the same ice crystals, within the closed system, is “sucked” back into the interbilayer space. This e ect is due to the temperature-dependent change in the balance between the hydration forces and the chemical force resulting from ice formation. It is important to note, that this leads to what amounts to an e ective spatial separation of the crystallized water from the biological surface.
For the temperature dependence of the lamellar spacing the following qualitative behavior of dl was observed [100]: During the cooling cycle of samples
16 Quasielastic Neutron Scattering: Applications |
385 |
equilibrated at a given relative humidity (r.h.), starting at room temperature (295 K), dl stays approximately constant down to a few degrees below the freezing point of bulk water (Tf = 273.15 K for pure H2O and 276.97 K for pure D2O). Then, at a temperature which we denote by Tfh, a discontinuity occurs, which is connected with a large decrease in dl by a step of the order of 20–30% . The hydration water remaining bound to the purple membrane below about 240 K is nonfreezing. Its amount was found to be hnf = 0.24(±0.02) [g D2O/g BR] for all PM samples equilibrated at room temperature in the presence of D2O vapor at ≥ 84 % r.h..
Figure 16.15 shows an example of such a study of a PM specimen [143]. The lamellar lattice constant of the PM multilayer stack is displayed as a function of temperature in a heating cycle; the sample had been equilibrated at room temperature in an atmosphere of 98% r.h. (D2O vapor). The experimental values (triangles) are shown together with a phenomenological model curve. The vertical arrows labeled with letters, indicate the approximate limiting temperatures, below which – in the course of the photocycle – BR does no longer return to the ground state, but stops at the intermediates K, L, M1 and M2, respectively.
|
100 |
120 |
140 |
160 |
180 |
200 |
220 |
240 |
260 |
280 |
300 |
|
90 |
|
|
|
|
|
|
|
Tmh |
|
90 |
(Å) |
85 |
PM 98 % r.h. at 295K |
|
|
|
|
|
85 |
|||
80 |
|
|
|
|
|
|
|
|
|
80 |
|
l |
|
|
|
|
|
|
|
|
|
||
d |
|
|
|
|
|
|
|
|
|
||
spacing |
75 |
|
|
|
|
|
|
|
|
|
75 |
70 |
|
|
|
|
|
|
|
|
|
70 |
|
lattice |
|
|
|
|
|
|
|
|
|
||
65 |
|
|
|
|
|
|
|
|
|
65 |
|
Lamellar |
60 |
|
|
|
|
|
|
|
M2 |
|
60 |
K |
|
|
|
L |
|
M1 |
|
|
|||
55 |
|
|
|
|
|
|
|
|
|
55 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
* |
|
|
50 |
|
|
|
|
|
|
|
Tmh |
50 |
|
|
100 |
120 |
140 |
160 |
180 |
200 |
220 |
240 |
260 |
280 |
300 |
|
|
|
|
|
|
T (K) |
|
|
|
|
|
Fig. 16.15. Lamellar lattice constant of a purple membrane multilayer stack as a function of temperature in a heating cycle; the sample (labeled H1) had been equilibrated at room temperature in an atmosphere of 98% r.h. (D2O vapor). The experimental values (triangles) obtained with the membrane-di ractometer V2 at BENSC in Berlin, are shown together with a phenomenological model curve. The vertical arrows indicate the approximate limiting temperatures, below which – in the course of the photocycle – BR does no longer return to the ground state, but stops at the intermediates K, L, M1, and M2, respectively [141]. Figure taken from [143]
386 R.E. Lechner et al.
Let us now discuss possible correlations between the T -dependent dehydration/rehydration behavior of PM and the variation of its dynamic structure with temperature. It has been shown by studies of purple membrane using quasielastic incoherent neutron scattering (QINS), that the ability of bacteriorhodopsin to functionally relax and complete the photocycle initiated by the absorption of a photon, is strongly correlated with the onset of low-frequency, large-amplitude “anharmonic” molecular motions. This manifests itself as the well-known “dynamical transition” [35, 144–147] starting from a low-temperature harmonic toward a high-temperature “anharmonic” regime [36]. More precisely, the “dynamical transition” announces the onset of overdamped vibrational and/or localized (i.e., spatially restricted) di u- sive (stochastic) molecular motions at temperatures in the neighborhood of 200 K. It is characterized by the appearance of quasielastic neutron scattering due to these motions9. This has also been made visible indirectly by studying the temperature dependence of elastic scattering (“elastic-window scan”, see Sect. 15.3.3 in Part I, this volume), which allows to extract an atomic meansquare displacement <u2>. In this case, for fully hydrated PM, the onset of di usive motions was found near 230 K by an analysis of <u2>. Here, the validity of a Gaussian approximation to the motion of all the hydrogen atoms in PM was assumed, which – for the employed energy resolution (IN13, with FWHM = 10 eV) – includes the frequencies between 1010 s−1 and 1013 s−1 [36]. In this experiment, a large Q-range (up to Q = 4.5 ˚A−1) was used . However at higher energy resolution (FWHM = 1 eV), and in a lower Q-range (up to Q = 1.8 ˚A−1), elastic-window scans with IN16 yielded the onset of a dynamical transition already near 150 K [148, 149]10. From
9Note that it is not possible to give one very precise temperature value for this “transition,” for two reasons: (i) there is a multitude of di erent molecular motions that are gradually activated with rising temperature; the transition is therefore occurring in a continuous way over a certain temperature range; (ii) any newly arising motion can only be detected, when it occurs within the energy window of the experiment; therefore, the temperature dependence of this “transition” phenomenon is also correlated with the resolution-dependent variation of experimental observability; see Sect. 15.3.3, Fig. 15.7, in Part I, this volume
10The elastic-window scan method has also been employed in numerous other biological experiments, in order to determine the temperature dependence of motional amplitudes (“mean-square displacements”), concerning for instance the environmentdependence of confined di usive protein motions (see Sect. 16.4.3), the dynamics of hydrated starch saccharides [150], the e ect of myelin basic protein on the dynamics of oriented lipids [151], proton mobilities in crambin and glutathionine S- transferase [152], the dynamic properties of an oriented lipid/DNA complex (proposed as DNA vector in gene therapy) [153], the influence of solvent composition on global dynamics of human butyrylcholinesterase powders [154], the relation between the glass-forming character and the cryoprotective capability of disaccharide–water mixtures [155, 156], and the comparison of the macromolecular dynamics in psychrophile, mesophile, thermophile, and hyperthermophile bacteria [157]