NMR spin grouping and correlation exchange analysis. Application to low hydration NaDNA paracrystals.

The NMR spin-grouping technique is applied to low hydration oriented fibers of NaDNA to study the role of exchange in determining the apparent (observed) spin relaxation of the system. The analysis proceeds in three steps: first, the apparent proton relaxation is measured at high fields, with both selective and nonselective inversion pulse sequences, and in the rotating frame. The spin-grouping technique is used in all spin-lattice relaxation measurements to provide the optimum apparent relaxation characterization of the sample. Next, all apparent results are analyzed for exchange. In this analysis the results from the high field and rotating frame experiments (which probe the exchange at two different time scales) are correlated to determine the inherent (or true) spin relaxation parameters of each of the proton groups in the system. The results of selective inversion T1 measurements are also incorporated into the exchange analysis. Finally, the dynamics of each spin group are inferred from the inherent relaxation characterization. The low hydration NaDNA structure is such that the exchange between the protons on the water and those on the NaDNA is limited, a priori, to dipolar mixing. The results of the exchange analysis indicate that the dipolar mixing between water and NaDNA protons is faster than the spin diffusion within the NaDNA proton group itself. The spin-diffusion on the macromolecule is the bottleneck for the exchange between the water protons and the NaDNA protons. The water protons serve as the relaxation sink both at high fields and in the rotating frame for the total NaDNA-water spin bath. The inherent relaxation of the water is characteristic of water undergoing anisotropic motion with a fast reorientational correlation time about one axis (5 X 10(-10) less than or equal to tau r less than or equal to 8 X 10(-9)S) which is about three orders of magnitude slower than that of water in the bulk; and a slow tumbling correlation time for this axis (1.5 x 10(-7) less than or equal to tau t less than or equal to 8 x 10(-7)S) which is two orders of magnitude slower yet.     

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.     

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.     

Ion solvation and water structure in potassium halide aqueous solutions

The structure of water and the nature of ionic hydration is explored in aqueous solutions of potassium fluoride, chloride, bromide and iodide over a range of concentrations up to 4.8 ion pairs per 100 water molecules, using the combined techniques of neutron diffraction with hydrogen isotope substitution. The diffraction data are interpreted using the method of empirical potential structure refinement, which attempts to build a three-dimensional model of the scattering system consistent with the diffraction data. The water structure is strongly perturbed in the first hydration shells of both anion and cation, but is found to be only mildly perturbed outside of this region, with the largest effects occurring with the smallest anion and highest concentrations. For the potassium ion there are strong orientational correlations in the first hydration shell, with the water molecules lying with their dipole moments pointing almost directly away from the cation on average, but with an angular spread of approximately +/-60 degrees which is mildly dependent on the anion type present. For all the anions the water molecules in the first shell are strongly oriented with one O-H vector pointing directly towards the anion on average, with an angular spread of approximately +/-10 degrees for F(-), increasing to approximately +/-22 degrees for I(-). For both anions and cations the second hydration shell is much more disordered than the first, but there is a weak pattern of orientational correlation which becomes more pronounced with the larger anions. There is some evidence that the fluoride ion structures water significantly in its first hydration shell, but not beyond. The findings throw further light on recent findings that the orientational relaxation time for water outside the first shell of dissolved ions is the same as in the bulk liquid.     

Determination of the number of water molecules in the proton pathway of bacteriorhodopsin using neutron diffraction data

It has been shown that water molecules participate in the proton pathway of bacteriorhodopsin. Large efforts have been made to determine with various biophysical methods the number of water molecules involved. Neutron diffraction H2O/D2O exchange experiments have been often used to reveal the position of water even with low-resolution diffraction data. With this technique, care must be taken with the limitations of the difference Fourier method which are commonly applied to analyze the data. In this paper we compare the results of the difference Fourier method applied to measured diffraction data (not presented here) and models with those from alternative methods introduced here: (1) a computer model calculation procedure to determine a label's scattering length density based on a comparison of intensity differences derived from models and intensity differences from our measurements; (2) a method based on the Parseval formula. Both alternative methods have been evaluated and tested using results of neutron diffraction experiments on purple membranes (Hauss et al. 1994). Our findings indicate that the difference Fourier method applied to low-resolution diffraction data can successfully determine the position of localized water molecules but underestimates their integrated scattering length density in the presence of labels in other positions. Furthermore, we present the results of neutron diffraction experiments on purple membranes performed to determine the number of water molecules in the projected area of the Schiff base at 86%, 75% and 57% relative humidity (r.h.). We found 19 +/- 2 exchangeable protons at 75% r.h., which means at least 8-9 water molecules are indispensable for normal pump function.     

Bilayer thickness and lipid interface area in unilamellar extruded 1,2-diacylphosphatidylcholine liposomes: a small-angle neutron scattering study

Small-angle neutron scattering (SANS) experiments have been performed on large unilamellar liposomes prepared from 1,2-dilauroylphosphatidylcholine (DLPC), 1,2-dimyristoyl-phosphatidylcholine (DMPC) and 1,2-distearoylphosphatidylcholine (DSPC) in heavy water by extrusion through polycarbonate filters with 500 A pores. The neutron scattering intensity I(Q) in the region of scattering vectors Q corresponding to 0.0015 A(-2) < or = Q(2) < or = 0.0115 A(-2) was fitted using a step function model of bilayer neutron scattering length density and supposing that the liposomes are spherical and have a Gaussian distribution of radii. Using the lipid volumetric data, and supposing that the thickness of bilayer polar region equals to d(H) = 9+/-1 A and the water molecular volume intercalated in the bilayer polar region is the same as in the aqueous bulk aqueous phase, the steric bilayer thickness d(L), the lipid surface area A(L) and the number of water molecules per lipid molecule N intercalated in the bilayer polar region were obtained: d(L) = 41.58+/-1.93 A, A(L) = 57.18+/-1.00 A(2) and N = 6.53+/-1.93 in DLPC at 20 degrees C, d(L) = 44.26+/-1.42 A, A(L) = 60.01+/-0.75 A(2) and N = 7.37+/-1.94 in DMPC at 36 degrees C, and d(L) = 49.77+/-1.52 A, A(L) = 64.78+/-0.46 A(2) and N = 8.67+/-1.97 in DSPC at 60 degrees C. After correcting for area thermal expansivity alpha approximately 0.00417 K(-1), the lipid surface area shows a decrease with the lipid acyl chain length at 60 degrees C: A(L) = 67.56+/-1.18 A(2) in DLPC, A(L) = 66.33+/-0.83 A(2) in DMPC and A(L) = 64.78+/-0.46 A(2) in DSPC. It is also shown that a joint evaluation of SANS and small-angle X-ray scattering on unilamellar liposomes can be used to obtain the value of d(H) and the distance of the lipid phosphate group from the bilayer hydrocarbon region d(H1).     

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.     

Hydration of ds-DNA and ss-DNA by neutron quasielastic scattering

Quasielastic neutron scattering measurements were performed in hydrated samples of ds-DNA and ss-DNA. The samples were hydrated in a high relative humidity atmosphere, and their final water content was 0.559 g H(2)O/g ds-DNA and 0.434 g H(2)O/g ss-DNA. The measurements were performed at 8 and 5.2 A for the ds-DNA sample, and at 5.2 A for the ss-DNA sample. The temperature was in both cases 298 K. Analysis of the obtained data indicates that in the ds-DNA sample we can distinguish two types of protons-those belonging to water molecules strongly attached to the ds-DNA surface and another fraction belonging to water that diffuses isotropically in a sphere of radius 2.8 A, with a local diffusion coefficient of 2.2 x 10(-5) cm(2) s(-1). For ss-DNA, on the other hand, no indication was found of motionally restricted or confined water. Further, the fraction of protons strongly attached to the ds-DNA surface corresponds to 0.16 g H(2)O/g ds-DNA, which equals the amount of water that is released by ds-DNA upon thermal denaturation, as studied by one of us (G.M.) by differential scanning calorimetry. This value also equals the difference between the critical hydration values of ds-DNA and ss-DNA, also determined by DSC. These results represent, thus, a completely independent measurement of water characteristics and behavior in ds- and ss-DNA at critical hydration values, and therefore substantiate the previous suggestions/conclusions of the results obtained by calorimetry.     

Index to Authors,Volume 7

Abstract

The thermally stimulated depolarization current (TSDC) measurements in frozen aqueous solutions, gels and solid layers of NaDNA show typically up to three dipolar overlapping peaks in the low-temperature range of 80—;150 K. Up to four discrete relaxation peaks have been observed at higher temperatures above 150 K. The low-temperature TSDC peaks are due to the dipolar relaxations of free and loosely bound water which crystallizes. Part of bound water especially in the first hydration shell of DNA molecule is at low temperatures in the form of glass. The transition of this glass from solidlike behavior to liquidlike behavior observed mainly in gels and solid samples is associated with a previously founded TSDC relaxation peak. The peak is at its maximum at 165- 250 K depending on the sample humidity. Existence of this relaxation in the samples with water contents in a broad range confirms, that the slowly relaxing shell (minimally 5–7 water molecules/nucleotide) closely associated with DNA double helix retains its characteristics. Also another peak of the high-temperature band at 180–205 K which was observed in the samples at hydration 2–1800 g H₂O/g dry NaDNA is due to a relaxation in the sample volume. At the highest temperatures relax the space charges trapped at the electrodes.     

Protein structure and hydration probed by SANS and osmotic stress

Interactions governing protein folding, stability, recognition, and activity are mediated by hydration. Here, we use small-angle neutron scattering coupled with osmotic stress to investigate the hydration of two proteins, lysozyme and guanylate kinase (GK), in the presence of solutes. By taking advantage of the neutron contrast variation that occurs upon addition of these solutes, the number of protein-associated (solute-excluded) water molecules can be estimated from changes in both the zero-angle scattering intensity and the radius of gyration. Poly(ethylene glycol) exclusion varies with molecular weight. This sensitivity can be exploited to probe structural features such as the large internal GK cavity. For GK, small-angle neutron scattering is complemented by isothermal titration calorimetry with osmotic stress to also measure hydration changes accompanying ligand binding. These results provide a framework for studying other biomolecular systems and assemblies using neutron scattering together with osmotic stress.     

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction.

The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Neutron scattering study of the binding of tRNAPhe to Escherichia coli phenylalanyl-tRNA synthetase

Escherichia coli phenylalanyl-tRNA synthetase has been characterized by small-angle neutron scattering. In solution (20 mM imidazole hydrochloride, pH 7.6, 10 mM 2-mercaptoethanol, and 0.1 mM ethylenediaminetetraacetic acid), this enzyme has a molecular weight of 227K +/- 20K with a radius of gyration of 48.3 +/- 0.6 A, independent of the presence of MgCl2 up to 50 mM. The change of the scattering upon adding tRNAPhe to the enzyme has been followed with 10 mM MgCl2 present in the buffer. One enzyme molecule is capable of binding two tRNAPhe molecules with affinity constants larger than 10(6) M-1. Parallel titration experiments in 73% 2H2O, close to the matching point of tRNA, show that the RG of the enzyme is not changed by the binding of one or two tRNAPhe molecules. These results are compared with quasi-electric light scattering studies [Holler, E., Wang, C. C., & Ford, N.C., Jr. (1981) Biochemistry 20, 861-867] where the addition of either MgCl2 or tRNAPhe was shown to cause dramatic changes of the apparent translational diffusion constant of phenylalanyl-tRNA synthetase.     

A model for water motion in crystals of lysozyme based on an incoherent quasielastic neutron-scattering study

This paper reports an incoherent quasielastic neutron scattering study of the single particle, diffusive motions of water molecules surrounding a globular protein, the hen egg-white lysozyme. For the first time such an analysis has been done on protein crystals. It can thus be directly related and compared with a recent structural study of the same sample. The measurement temperature ranged from 100 to 300 K, but focus was on the room temperature analysis. The very good agreement between the structural and dynamical studies suggested a model for the dynamics of water in triclinic crystals of lysozyme in the time range approximately 330 ps and at 300 K. Herein, the dynamics of all water molecules is affected by the presence of the protein, and the water molecules can be divided into two populations. The first mainly corresponds to the first hydration shell, in which water molecules reorient themselves fivefold to 10-fold slower than in bulk solvent, and diffuse by jumps from hydration site to hydration site. The long-range diffusion coefficient is five to sixfold less than for bulk solvent. The second group corresponds to water molecules further away from the surface of the protein, in a second incomplete hydration layer, confined between hydrated macromolecules. Within the time scale probed they undergo a translational diffusion with a self-diffusion coefficient reduced approximately 50-fold compared with bulk solvent. As protein crystals have a highly crowded arrangement close to the packing of macromolecules in cells, our conclusion can be discussed with respect to solvent behavior in intracellular media: as the mobility is highest next to the surface, it suggests that under some crowding conditions, a two-dimensional motion for the transport of metabolites can be dominant.     

Thermotropic properties of bipolar lipids of Sulfolobus solfataricus and of their mixtures with dipalmitoylphosphatidylcholine

The thermotropic properties of the bipolar lipids, glycerol dialkylglycerol tetraether (GDGT) and glycerol dialkylnonitol tetraether (GDNT), were determined at different degrees of hydration and in mixtures with dipalmitoylphosphatidylcholine (DPPC). The number of water molecules rendered unfreezable by the GDNT molecule is 10+/-1.5 and that by the GDGT molecule 2.8+/-0.7 or about 1.1-1.5 H2O molecules per OH group. Binding of water molecules causes randomization of the two polar heads from the oriented form prevailing in the dry state. The hydration seems to be a cooperative process extending over a whole lipid domain. DPPC added in small amounts to GDNT interacts preferentially with the nonitol halves of the molecules separating them from the glycerol half molecules. In the cooperative interaction domain each DPPC molecule is surrounded by up to six GDNT molecules. Cooperative domains formed during the interaction of DPPC with GDGT are less pronounced. In both cases they affect the thermotropic properties of the system.     

Dynamical properties of phospholipid bilayers from computer simulation

We present the results of a 10-ns molecular dynamics simulation of a dipalmitoylphosphatidylcholine/water system. The main emphasis of the present study is on the investigation of the stability over a long time and the dynamic properties of the water/membrane system. The motion of the lipid molecules is characterized by the center of mass movement and the displacement of individual atom groups. Because of the slow movement of the headgroup atoms, their contributions to the dipole potential vary slowly and with a large amplitude. Nevertheless, the water molecules compensate the strong fluctuations and maintain an almost constant total dipole potential. From the lateral displacement of the center of masses, we calculate the lateral diffusion coefficient to be Dlat = (3 +/- 0.6) x 10(-7) cm2/s, in agreement with neutron scattering results. The rotational motion is also investigated in our simulations. The calculated value for the rotational diffusion coefficient parallel to the molecular long axis, D = (1.6 +/- 0.1) x 10(8) s-1, is in good agreement with the experiment.     

Hydration of DNA: a neutron scattering study of oriented NaDNA

A preliminary neutron-scattering study of wet-spun oriented NaDNA samples with various contents of H₂O or D₂O has been performed to obtain information about the hydration of NaDNA. The preferred orientation was in the direction of the crystallographic c axis. Perpendicular to this direction the samples could be regarded as poly-crystalline. It has been found that the preferred crystallite orientation was also reflected in the water structure. This fact, taken together with the small energy widths of the measured near-elastic peaks indicates that water molecules in the sample are relatively strongly bound to the DNA helix either directly or via other water molecules. Angular distribution measurements of elastically scattered neutrons were performed with the momentum transfer vector directed along and perpendicular to the DNA helix axis. Different Debye-Waller factors were obtained, indicating different water dynamics in the two directions. Comparison between the H₂O and D₂O measurements suggests that the difference is caused by motions of rotational nature. The shape of the frequency distribution function for the H₂O molecules in the NaDNA sample was somewhat similar to that which has been measured in ice.     

Hydrophobic thickness, lipid surface area and polar region hydration in monounsaturated diacylphosphatidylcholine bilayers: SANS study of effects of cholesterol and beta-sitosterol in unilamellar vesicles

The influence of a mammalian sterol cholesterol and a plant sterol beta-sitosterol on the structural parameters and hydration of bilayers in unilamellar vesicles made of monounsaturated diacylphosphatidylcholines (diCn:1PC, n=14-22 is the even number of acyl chain carbons) was studied at 30 degrees C using small-angle neutron scattering (SANS). Recently published advanced model of lipid bilayer as a three-strip structure was used with a triangular shape of polar head group probability distribution (Kucerka et al., Models to analyze small-angle neutron scattering from unilamellar lipid vesicles, Physical Review E 69 (2004) Art. No. 051903). It was found that 33 mol% of both sterols increased the thickness of diCn:1PC bilayers with n=18-22 similarly. beta-sitosterol increased the thickness of diC14:1PC and diC16:1PC bilayers a little more than cholesterol. Both sterols increased the surface area per unit cell by cca 12 A(2) and the number of water molecules located in the head group region by cca 4 molecules, irrespective to the acyl chain length of diCn:1PC. The structural difference in the side chain between cholesterol and beta-sitosterol plays a negligible role in influencing the structural parameters of bilayers studied.