Hydration dynamics near a model protein surface

The evolution of water dynamics from dilute to very high concentration solutions of a prototypical hydrophobic amino acid with its polar backbone, N-acetyl-leucine-methylamide (NALMA), is studied by quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulation for both the completely deuterated and completely hydrogenated leucine monomer. The NALMA-water system and the QENS data together provide a unique study for characterizing the dynamics of different hydration layers near a prototypical hydrophobic side chain and the backbone of which it is attached. We observe several unexpected features in the dynamics of these biological solutions under ambient conditions. The NALMA dynamics shows evidence of de Gennes narrowing, an indication of coherent long timescale structural relaxation dynamics. The translational and rotational water dynamics at the highest solute concentrations are found to be highly suppressed as characterized by long residential time and slow diffusion coefficients. The analysis of the more dilute concentration solutions models the first hydration shell with the 2.0 M spectra. We find that for outer layer hydration dynamics that the translational diffusion dynamics is still suppressed, although the rotational relaxation time and residential time are converged to bulk-water values. Molecular dynamics analysis of the first hydration shell water dynamics shows spatially heterogeneous water dynamics, with fast water motions near the hydrophobic side chain, and much slower water motions near the hydrophilic backbone. We discuss the hydration dynamics results of this model protein system in the context of protein function and protein-protein recognition.     

The observation of structural transitions of a single protein molecule

Coherent neutron scattering measurements of an amorphous, in vivo deuterated C-phycocyanin are compared with a calculation of the individual protein molecule's coherent static structure factor. Both show the significant features associated with known structure factors of several amorphous materials, most notably, an unusually sharp first diffraction peak occurring near 1.4 A(-1). We show that in the protein, such a peak results from the product of a form factor associated with correlations of atoms within individual amino acids and a structural term expressing inter-amino-acid correlations. The measurement, interpreted through behavior of the first diffraction peak, indicates that inter-amino-acid correlations - a measure of the protein's medium-range structure - undergo transitions which are primarily related to hydration rather than to temperature.     

The molecular basis of the temperature- and pH-induced conformational transitions in elastin-based peptides

Elastin undergoes an inverse temperature transition and collapses at high temperatures in both simulation and experiment. We investigated a pH-dependent modification of this transition by simulating a glutamic acid (Glu)-substituted elastin at varying pHs and temperatures. The Glu-substituted peptide collapsed at higher temperature than the unsubstituted elastin when Glu was charged. The charge effects could be reversed by neutralization of the Glu carboxyl groups at low pH, and in that case the peptide collapsed at a lower temperature. The collapse was accompanied by the formation of beta-turns and short distorted beta-sheets. Formation of contacts between hydrophobic side chains drives the collapse at high temperature, but interactions between water and polar groups (Glu and main chain) can attenuate this effect at high pH. The overall competition and balance of the polar and nonpolar groups determined the conformational states of the peptide. Water hydration contributed to the conformational transition, and the peptide and its hydration shell must be considered. Structurally, waters near polar residues mainly formed hydrogen bonds with the protein atoms, while waters around the hydrophobic side chains tended to be parallel to the peptide groups to maximize water-water interactions.     

Alteration of Water Structure by Peptide Clusters Revealed by Neutron Scattering in the Small-Angle Region (below 1???1)

Solution scattering of neutrons and x-rays can provide direct information on local interactions of importance for biomolecular folding and structure. Here, neutron scattering experiments are combined with molecular-dynamics simulation to interpret the scattering signal of a series of dipeptides with varying degrees of hydrophobicity (GlyAla, GlyPro, and AlaPro) in concentrated aqueous solution (1:20 solute/water ratio) in which the peptides form large segregates (up to 50–60 amino acids). Two main results are found: 1), the shift to lower Q of the so-called water-ring peak (Q ≈ 2 Å⁻¹) arises mainly from an overlap of water-peptide and peptide-peptide correlations in the region of 1.3 < Q < 2 Å⁻¹, rather than from a shift of the water signal induced by the presence of the clusters; and 2), in the low-Q region (Q ≈ 0.6 Å⁻¹) a positive peak is observed originating from both the solute-solute correlations and changes in the water structure induced by the formation of the clusters. In particular, the water molecules are found to be more connected than in the bulk with hydrogen-bonding directions tangential to the exposed hydrophobic surfaces, and this effect increases with increasing peptide hydrophobicity. This work demonstrates that important information on the (hydrophobic) hydration of biomolecules can be obtained in the very-small-angle region.     

Neutron scattering from water in DNA and muscle tissue

Neutron scattering studies of rotational modes of water in DNA and muscle tissue are reported. Spectra from wet DNA and rat muscle are compared with those from ice and water. No evidence is found for an ice-like water component in biological systems. The relevance of neutron scattering for such a study is examined and possible discrepancies with other measurements are discussed. The similarity of the neutron spectra for wet DNA and muscle tissue is taken as evidence that water-nucleotide and water-peptide interactions are similar to water-water interactions when probed on a short time scale.     

Vibrational modes of hemoglobin in red blood cells.

Equine red blood cells were washed in saline heavy water (2H2O) to exchange the hydrogen atoms of the non-hemoglobin components with deuterons. This led to novel neutron scattering measurements of protein vibrations within a cellular system and permitted a comparison with inelastic neutron scattering measurements on purified horse hemoglobin, either dry or wetted with 2H2O. As a function of wavevector transfer Q and the frequency transfer v the neutron response typified by the dynamic structure factor S(Q, v) was found to be similar for extracted and cellular hemoglobin at low and high temperatures. At 77 K, in the cells, a peak in S(Q, v) due to the protein was found near 0.7 THz, approximately half the frequency of a strong peak in the aqueous medium. Measurements at higher temperatures (170 and 230 K) indicated similar small shifts downwards in the peak frequencies of both components. At 260 K the low frequency component became predominantly quasielastic, but a significant inelastic component could still be ascribed to the aqueous scattering. Near 295 K the frequency responses of both components were similar and centered near zero. When scattering due to water is taken into account it appears that the protein neutron response in, or out of, red blood cells is little affected by hydration in the low frequency regime where Van der Waals forces are thought to be effective.     

Water structure around dipeptides in aqueous solutions

The bulk water structure around small peptide fragments-glycyl-L: -alanine, glycyl-L: -proline and L: -alanyl-L: -proline-has been determined by a combination of neutron diffraction with isotopic substitution and empirical potential structural refinement techniques. The addition of each of the dipeptides to water gives rise to decreased water-water coordination in the surrounding water solvent. Additionally, both the O(w)-O(w) radial distribution functions and the water-water spatial density functions in all of the solutions indicate an electrostrictive effect in the second water coordination shell of the bulk water network. This effect is not observed in similar experiments on the amino acid L: -proline alone in solution, which is one component of two of the peptides measured here.     

Short-range order and collective dynamics of DMPC bilayers: a comparison between molecular dynamics simulations, X-ray, and neutron scattering experiments

We present an extensive comparison of short-range order and short wavelength dynamics of a hydrated phospholipid bilayer derived by molecular dynamics simulations, elastic x-ray, and inelastic neutron scattering experiments. The quantities that are compared between simulation and experiment include static and dynamic structure factors, reciprocal space mappings, and electron density profiles. We show that the simultaneous use of molecular dynamics and diffraction data can help to extract real space properties like the area per lipid and the lipid chain ordering from experimental data. In addition, we assert that the interchain distance can be computed to high accuracy from the interchain correlation peak of the structure factor. Moreover, it is found that the position of the interchain correlation peak is not affected by the area per lipid, while its correlation length decreases linearly with the area per lipid. This finding allows us to relate a property of the structure factor quantitatively to the area per lipid. Finally, the short wavelength dynamics obtained from the simulations and from inelastic neutron scattering are analyzed and compared. The conventional interpretation in terms of the three-effective-eigenmode model is found to be only partly suitable to describe the complex fluid dynamics of lipid chains.     

Pressure denaturation of staphylococcal nuclease studied by neutron small-angle scattering and molecular simulation

We studied the pressure-induced folding/unfolding transition of staphylococcal nuclease (SN) over a pressure range of approximately 1-3 kilobars at 25 degrees C by small-angle neutron scattering and molecular dynamics simulations. We find that applying pressure leads to a twofold increase in the radius of gyration derived from the small-angle neutron scattering spectra, and P(r), the pair distance distribution function, broadens and shows a transition from a unimodal to a bimodal distribution as the protein unfolds. The results indicate that the globular structure of SN is retained across the folding/unfolding transition although this structure is less compact and elongated relative to the native structure. Pressure-induced unfolding is initiated in the molecular dynamics simulations by inserting water molecules into the protein interior and applying pressure. The P(r) calculated from these simulations likewise broadens and shows a similar unimodal-to-bimodal transition with increasing pressure. The simulations also reveal that the bimodal P(r) for the pressure-unfolded state arises as the protein expands and forms two subdomains that effectively diffuse apart during initial stages of unfolding. Hydrophobic contact maps derived from the simulations show that water insertions into the protein interior and the application of pressure together destabilize hydrophobic contacts between these two subdomains. The findings support a mechanism for the pressure-induced unfolding of SN in which water penetration into the hydrophobic core plays a central role.     

Organotin compounds alter the physical organization of phosphatidylcholine membranes

Organotin compounds have a broad range of biological activities and are ubiquitous contaminants in the environment. Their toxicity mainly lies in their action on the membrane. In this contribution we study the interaction of tributyltin and triphenyltin with model membranes composed of phosphatidylcholines of different acyl chain lengths using differential scanning calorimetry, (31)P-nuclear magnetic resonance, X-ray diffraction and infrared spectroscopy. Organotin compounds broaden the main gel to liquid-crystalline phase transition, shift the transition temperature to lower values and induce the appearance of a new peak below the main transition peak. These effects are more pronounced in the case of tributyltin and are quantitatively larger as the phosphatidylcholine acyl chain length decreases. Both tributyltin and triphenyltin increase the enthalpy change of the transition in all the phosphatidylcholine systems studied except in dilauroylphosphatidylcholine. Organotin compounds do not affect the macroscopic bilayer organization of the phospholipid but do affect the degree of hydration of its carbonyl moiety. The above evidence supports the idea that organotin compounds are located in the upper part of the phospholipid palisade near the lipid/water interface.     

Definition of lipid membrane structural parameters from neutronographic experiments with the help of the strip function model.

Neutron diffraction is an effective method for investigating model and biological membranes. Yet, to obtain accurate structural information it is necessary to use deuterium labels and much time is needed to acquire experimental data as there are a large number of diffraction reflections to register. This paper offers a way to define the hydrophobic boundary position in lipid membranes with high accuracy and for this purpose it is sufficient to take into consideration three structural factors. The method is based on modeling the density of the neutron diffraction amplitude rho(x) in the direction of the bilayer plane normal by means of a strip function, but it also takes into consideration the fact that the multiplication of the strip function amplitude rho i by the step width zi-zi-1 makes the sum of neutron scattering amplitudes of the atoms included in the step region. On the basis of the analysis of a large number of experimental data for different membranes, the effectiveness of this method in the determination of the position of hydrophilic/hydrophobic boundary is demonstrated, including the case of various rho(x) modifications in the region of polar heads and also the different phase states of membranes. However, it is shown in the present paper that the strip function model is not an adequate instrument for the determination of other structural parameters of membranes.     

The nature of the hydrophobic binding of small peptides at the bilayer interface: implications for the insertion of transbilayer helices

One method of obtaining useful information about the physical chemistry of peptide/bilayer interactions is to relate thermodynamic parameters of the interactions to structural parameters obtained by diffraction methods. We report here the results of the application of this approach to interactions of hydrophobic tripeptides of the form Ala-X-Ala-O-tert-butyl with lipid bilayers. The thermodynamic constants (delta Gt, delta Ht, and delta St) for the transfer of the tripeptides from water into DMPC vesicles were determined for X = Leu, Phe, and Trp and found to be consistent with those expected for hydrophobic interactions above the phase transition of DMPC. Combining these results with the earlier ones of Jacobs and White [(1986) Biochemistry 25, 2605-2612], the favorable free energies of transfer with different amino acids in the -X- position increase in the order Gly less than Ala less than Leu less than Phe less than Trp in agreement with the Nozaki and Tanford [(1971) J. Biol. Chem. 246, 2211-2217] hydrophobicity scale. Determination of the location of Ala-[2H5]Trp-Ala-O-tert-butyl in oriented DOPC bilayers by neutron diffraction shows that the most hydrophobic peptide of the series is confined to the bilayer headgroup/water region. Refinement of the diffraction measurements shows that only 13% of the tryptophan is associated with the hydrocarbon core. The distribution of the water tends to mirror that of the peptide. Unlike peptide-free bilayers, 5% of the water penetrates the hydrocarbon, which is about 100-fold greater than expected. A quantitative thermodynamic analysis of the interfacial binding of the peptides suggests that (1) the hydrophobic interactions are 60-70% complete upon binding at the bilayer interface, (2) the interface is likely to play an important role in helix formation and insertion, (3) the hydrogen bond status of amino acid side chains is crucial to insertion, and (4) an a priori lack of knowledge of the status of such bonds could limit the precision of hydrophobicity plots. We introduce an interfacial hydrophobicity scale, IFH(h), with a variable hydrogen bond parameter (h) that permits one to consider explicitly hydrogen bonding in transbilayer helix searches.     

Lateral interactions among phosphatidylcholine and phosphatidylethanolamine head groups in phospholipid monolayers and bilayers

We develop theory for the lateral interactions among the zwitterionic head groups of phospholipids in monolayers and bilayers, particularly phosphatidylcholine (PC) and phosphatidylethanolamine (PE). With the P- end of the head group anchored at the water/hydrocarbon interface, a balance of two effects dictates the angle that the P--N+ dipole makes with respect to the plane of the bilayer: N+ is driven toward water due to the (Born) electrostatic free energy, but the hydrophobic effect drives the methyl and methylene groups around the N+ charge toward the hydrocarbon. The only adjustable parameter of the model is the average fluctuation of the oil/water interface or, alternatively, the dielectric constant of the hydrocarbon phase. The model predicts that at 5 degrees C the head group dipole should lie largely in the bilayer plane, in accord with X-ray, neutron diffraction, and NMR studies. The theory makes the novel prediction that the N+ end of the dipole becomes increasingly submerged in hydrocarbon with increasing temperature, leading to strongly enhanced lateral repulsion between PC head groups. This prediction is in good agreement with second and third viral coefficients of monolayer lateral pressures, and with the temperature dependence of the former. The theoretical model is consistent with head group fluctuations measured by neutron diffraction of PC and PE bilayers. Because PE has a smaller hydrophobic cluster near N+, its lateral repulsion should be much smaller and less temperature dependent than for PC, also in agreement with equation-of-state measurements. This suggests why at high density PE monolayers have higher melting temperatures than PC monolayers and more propensity for reversed curvature.     

New insights into the structure and hydration of a stratum corneum lipid model membrane by neutron diffraction

The structure and hydration of a stratum corneum (SC) lipid model membrane composed of N-(-hydroxyoctadecanoyl)-phytosphingosine (CER6)/cholesterol (Ch)/palmitic acid (PA)/cholesterol sulfate (ChS) were characterized by neutron diffraction. The neutron scattering length density across the SC lipid model membrane was calculated from measured diffraction peak intensities. The internal membrane structure and water distribution function across the bilayer were determined. The low hydration of the intermembrane space is a major feature of the SC lipid model membrane. The thickness of the water layer in the SC lipid model membrane is about 1 Å at full hydration. For the composition 55% CER6/25% Ch/15% PA/5% ChS, in a partly dehydrated state (60% humidity) and at 32°C, the lamellar repeat distance and the membrane thickness have the same value of 45.6 Å . The hydrophobic region of the membrane has a thickness of 31.2 Å . A decrease of the Ch content increases the membrane thickness. The water diffusion through the SC lipid model multilamellar membrane is a considerably slow process relative to that through phospholipid membranes. In excess water, the membrane hydration follows an exponential law with two characteristic times of 93 and 44 min. At 81°C and 97% humidity, the membrane separates into two phases with repeat distances of 45.8 and 40.5 Å . Possible conformations of CER6 molecules in the dry and hydrated multilayers are discussed.     

Site-specific deuterium order parameters and membrane-bound behavior of a peptide fragment from the intracellular domain of HIV-1 gp41.

The behavior of the cytolytic peptide fragment 828-848 (P828) from the carboxy-terminus of the envelope glycoprotein gp41 of HIV-1 in membranes was investigated by solid-state 2H NMR on P828 with the selectively deuterated isoleucines I3, I13, I16, and I20. The quadrupole splittings of the I3 side chain show significant sensitivity to the main phase-transition temperature of the lipid, consistent with partial penetration of the N-terminal peptide region into the hydrophobic core of the membrane. In contrast, the quadrupole splittings of I13, I16, and I20 are in agreement with a location of the C-terminal portion of the peptide near the lipid/water interface. The perturbation of the bilayer by the peptide was studied by 2H NMR on sn-1 chain deuterated 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoserine membranes. Peptide incorporation results in a significant reduction of lipid chain order toward the bilayer center, but only a modest reduction near the lipid glycerol. These observations suggest a penetration of the partially structured peptide backbone into the membrane/water interface region that reduces lateral packing density and decreases order in the hydrophobic core. In addition, the structure of the peptide was investigated free in water and bound to SDS micelles by high-resolution NMR. P828 is unstructured in water but exists in a flexible partially helical conformation when bound to negatively charged liposomes or micelles. The flexible helix covers the first 14 residues of the peptide, whereas the C-terminus of the peptide, where three of the six positively charged arginine residues are located, appears to be unstructured. The peptide-induced changes in lipid chain order profiles indicate that membrane curvature stress is the driving force for the cytolytic behavior of P828.     

Single-channel water permeabilities of Escherichia coli aquaporins AqpZ and GlpF

From equilibrium molecular dynamics simulations we have determined single-channel water permeabilities for Escherichia coli aquaporin Z (AqpZ) and aquaglyceroporin GlpF with the channels embedded in lipid bilayers. GlpF's osmotic water permeability constant pf exceeds by 2-3 times that of AqpZ and the diffusive permeability constant (pd) of GlpF is found to exceed that of AqpZ 2-9-fold. Achieving complete water selectivity in AqpZ consequently implies lower transport rates overall relative to the less selective, wider channel of GlpF. For AqpZ, the ratio pf/pd congruent with 12 is close to the average number of water molecules in the channel lumen, whereas for GlpF, pf/pd congruent with 4. This implies that single-file structure of the luminal water is more pronounced for AqpZ, the narrower channel of the two. Electrostatics profiles across the pore lumens reveal that AqpZ significantly reinforces water-channel interactions, and weaker water-water interactions in turn suppress water-water correlations relative to GlpF. Consequently, suppressed water-water correlations across the narrow selectivity filter become a key structural determinant for water permeation causing luminal water to permeate slower across AqpZ.     

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.     

The structure of polyunsaturated lipid bilayers important for rhodopsin function: a neutron diffraction study

The structure of oriented 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine bilayers with perdeuterated stearoyl- or docosahexaenoyl hydrocarbon chains was investigated by neutron diffraction. Experiments were conducted at two different relative humidities, 66 and 86%. At both humidities we observed that the polyunsaturated docosahexaenoyl chain has a preference to reside near the lipid water interface. That leaves voids in the bilayer center that are occupied by saturated stearoyl chain segments. This uneven distribution of saturated- and polyunsaturated chain densities is likely to result in membrane elastic stress that modulates function of integral receptor proteins like rhodopsin.     

Effect of lesions on the dynamics of DNA on the picosecond and nanosecond timescales using a polarity sensitive probe

This paper explores the effects of structural modifications on the fast dynamics of DNA and the ability of time-resolved Stokes shift spectroscopy to measure those changes. The time-resolved Stokes shift of a synthetic coumarin base-pair replacement within an oligomer is measured between 40 ps and 40 ns. Comparisons are made between 17mers without modification, with a deleted base near the coumarin and with the coumarin placed near the end of the oligomer. The deletion of a next-to-nearest-neighbor base pair does not change the subnanosecond dynamics, but does cause an additional motion with a time constant of ~20 ns. A candidate for this motion is the flipping of the abasic sugar out of the helix and the concomitant intrusion of water into the interior of the helix. A nearby chain end causes little change in the dynamics after 1 ns but leads to a reduction in the amplitude of the dynamics between 40 ps and 1 ns. We suggest that at the chain end, where DNA on one side of the probe has been replaced by water, the charge- stabilizing dynamics have the same overall amplitude, but that much of the relaxation occurs before the start of the measurement time window.     

Protein folding in membranes: insights from neutron diffraction studies of a membrane beta-sheet oligomer

Studies of the assembly of the hexapeptide Acetyl-Trp-Leu₅ (AcWL₅) into β-sheets in membranes have provided insights into membrane protein folding. Yet, the exact structure of the oligomer in the lipid bilayer is unknown. Here we use neutron diffraction to study the disposition of the peptides in bilayers. We find that pairs of adjacent deuterium-labeled leucines have no well-defined peak or dip in the transmembrane distribution profiles, indicative of heterogeneity in the depth of membrane insertion. At the same time, the monomeric homolog AcWL₄ exhibits a homogeneous, well-defined, interfacial location in neutron diffraction experiments. Thus, although the bilayer location of monomeric AcWL₄ is determined by hydrophobicity matching or complementarity within the bilayer, the AcWL₅ molecules in the oligomer are positioned at different depths within the bilayer because they assemble into a staggered transmembrane β-sheet. The AcWL₅ assembly is dominated by protein-protein interactions rather than hydrophobic complementarity. These results have implications for the structure and folding of proteins in their native membrane environment and highlight the importance of the interplay between hydrophobic complementarity and protein-protein interactions in determining the structure of membrane proteins.