Ultrafast X-ray probing of water structure below the homogeneous ice nucleation temperature.

Water has a number of anomalous physical properties, and some of these become drastically enhanced on supercooling below the freezing point. Particular interest has focused on thermodynamic response functions that can be described using a normal component and an anomalous component that seems to diverge at about 228 kelvin (refs 1-3). This has prompted debate about conflicting theories that aim to explain many of the anomalous thermodynamic properties of water. One popular theory attributes the divergence to a phase transition between two forms of liquid water occurring in the 'no man's land' that lies below the homogeneous ice nucleation temperature (TH) at approximately 232 kelvin and above about 160 kelvin, and where rapid ice crystallization has prevented any measurements of the bulk liquid phase. In fact, the reliable determination of the structure of liquid water typically requires temperatures above about 250 kelvin. Water crystallization has been inhibited by using nanoconfinement, nanodroplets and association with biomolecules to give liquid samples at temperatures below TH, but such measurements rely on nanoscopic volumes of water where the interaction with the confining surfaces makes the relevance to bulk water unclear. Here we demonstrate that femtosecond X-ray laser pulses can be used to probe the structure of liquid water in micrometre-sized droplets that have been evaporatively cooled below TH. We find experimental evidence for the existence of metastable bulk liquid water down to temperatures of 227(-1)(+2) kelvin in the previously largely unexplored no man's land. We observe a continuous and accelerating increase in structural ordering on supercooling to approximately 229 kelvin, where the number of droplets containing ice crystals increases rapidly. But a few droplets remain liquid for about a millisecond even at this temperature. The hope now is that these observations and our detailed structural data will help identify those theories that best describe and explain the behaviour of water.

The inhomogeneous structure of water at ambient conditions.

Small-angle X-ray scattering (SAXS) is used to demonstrate the presence of density fluctuations in ambient water on a physical length-scale of approximately 1 nm; this is retained with decreasing temperature while the magnitude is enhanced. In contrast, the magnitude of fluctuations in a normal liquid, such as CCl(4), exhibits no enhancement with decreasing temperature, as is also the case for water from molecular dynamics simulations under ambient conditions. Based on X-ray emission spectroscopy and X-ray Raman scattering data we propose that the density difference contrast in SAXS is due to fluctuations between tetrahedral-like and hydrogen-bond distorted structures related to, respectively, low and high density water. We combine our experimental observations to propose a model of water as a temperature-dependent, fluctuating equilibrium between the two types of local structures driven by incommensurate requirements for minimizing enthalpy (strong near-tetrahedral hydrogen-bonds) and maximizing entropy (nondirectional H-bonds and disorder). The present results provide experimental evidence that the extreme differences anticipated in the hydrogen-bonding environment in the deeply supercooled regime surprisingly remain in bulk water even at conditions ranging from ambient up to close to the boiling point.

Increasing correlation length in bulk supercooled H2O, D2O, and NaCl solution determined from small angle x-ray scattering.

Using small angle x-ray scattering, we find that the correlation length of bulk liquid water shows a steep increase as temperature decreases at subzero temperatures (supercooling) and that it can, similar to the thermodynamic response functions, be fitted to a power law. This indicates that the anomalous properties of water are attributable to fluctuations between low- and high-density regions with rapidly growing average size upon supercooling. The substitution of H(2)O with D(2)O, as well as the addition of NaCl salt, leads to substantial changes of the power law behavior of the correlation length. Our results are consistent with the proposed existence of a liquid-liquid critical point in the deeply supercooled region but do not exclude a singularity-free model.

Scaling behavior of the reorientation kinetics of block copolymers exposed to electric fields.

We have followed the reorientation kinetics of various block copolymer solutions exposed to an external electric DC field. The characteristic time constants follow a power law indicating that the reorientation is driven by a decrease in electrostatic energy. Moreover, the observed exponent suggests an activated process in line with the expectations for a nucleation and growth process. When properly scaled, the data collapse onto a single master curve spanning several orders of magnitude both in reduced time and in reduced energy. The power law dependence of the rate of reorientation derived from computer simulations based on dynamic density functional theory agrees well with the experimental observations. First experiments in AC electric fields at sufficiently high frequencies confirm the notion that the reorientation process is dominated by differences in the dielectric constants rather than by mobile ions.

Probing fast kinetics in complex fluids by combined rapid mixing and small-angle X-ray scattering.

Stopped-flow mixing coupled to small-angle X-ray scattering (SAXS) is an established technique for investigating structural kinetics in solution down to the millisecond range. More recently, the emphasis has shifted to the sub-millisecond range using continuous flow microfluidic mixing devices. The aim of this article is to review the present status and limitations when applying mixing techniques to a wide range of soft matter and biological systems. In the case of SAXS, special consideration of the mixing quality is necessary for a quantitative description of the scattered intensity. This is demonstrated through two representative examples involving protein refolding and micellar self-assembly.

Supramolecular structures of peptide assemblies in membranes by neutron off-plane scattering: method of analysis.

In a previous paper (Yang et al., Biophys. J. 75:641-645, 1998), we showed a simple, efficient method of recording the diffraction patterns of supramolecular peptide assemblies in membranes where the samples were prepared in the form of oriented multilayers. Here we develop a method of analysis based on the diffraction theory of two-dimensional liquids. Gramicidin was used as a prototype model because its pore structure in membrane in known. At full hydration, the diffraction patterns of alamethicin and magainin are similar to gramicidin except in the scale of q (the momentum transfer of scattering), clearly indicating that both alamethicin and magainin form pores in membranes but of different sizes. When the hydration of the multilayer samples was decreased while the bilayers were still fluid, the in-plane positions of the membrane pores became correlated from one bilayer to the next. We believe that this is a new manifestation of the hydration force. The effect is most prominent in magainin patterns, which are used to demonstrate the method of analysis. When magainin samples were further dehydrated or cooled, the liquid-like diffraction turned into crystal-like patterns. This discovery points to the possibility of investigating the supramolecular structures with high-order diffraction.

Dynamics of the self-assembly of unilamellar vesicles.

We have studied the transient stages in the formation of unilamellar vesicles with millisecond time resolution. The self-assembly was initiated by rapid mixing of equimolar amounts of anionic and zwitterionic micelles and the transient micellar entities were probed by time-resolved small-angle x-ray scattering. Within the mixing time, original micelles transformed to disklike micelles which evolved further to a critical size and then closed to form monodisperse unilamellar vesicles within a second. Subsequent growth led to an unexpected broadening of the vesicle size distribution.

Crystallization of antimicrobial pores in membranes: magainin and protegrin.

Membrane pores spontaneously formed by antimicrobial peptides in membranes were crystallized for the first time by manipulating the sample hydration and temperature. Neutron diffraction shows that magainins and protegrins form stable pores in fully hydrated fluid membranes. At lower hydration levels or low temperature, the membrane multilayers crystallize. In one crystalline phase, the pores in each bilayer arrange in a regular hexagonal array and the bilayers are stacked into a hexagonal ABC lattice, corresponding to the cubic close-packed structure of spheres. In another crystalline phase, the bilayers are modulated into the rippled multilamellae, corresponding to a 2D monoclinic lattice. The phase diagrams are described. Crystallization of the membrane pores provides possibilities for diffraction studies that might provide useful information on the pore structures.

Neutron off-plane scattering of aligned membranes. I. Method Of measurement.

We describe a method of measuring neutron scattering of aligned membranes with the momentum transfer oriented parallel or partly perpendicular to the plane of the membranes. The method obtains the complete information for the structures within fluid membranes obtainable by scattering. Data from alamethicin- and magainin-induced pores are presented. Although the in-plane scattering curves of these two peptides are similar to each other, their off-plane scattering patterns are strikingly distinct. Magainin pores exhibit intermembrane correlations.

Barrel-stave model or toroidal model? A case study on melittin pores.

Transmembrane pores induced by amphiphilic peptides, including melittin, are often modeled with the barrel-stave model after the alamethicin pore. We examine this assumption on melittin by using two methods, oriented circular dichroism (OCD) for detecting the orientation of melittin helix and neutron scattering for detecting transmembrane pores. OCD spectra of melittin were systematically measured. Melittin can orient either perpendicularly or parallel to a lipid bilayer, depending on the physical condition and the composition of the bilayer. Transmembrane pores were detected when the helices oriented perpendicularly to the plane of the bilayers, not when the helices oriented parallel to the bilayers. The evidence that led to the barrel-stave model for alamethicin and that to the toroidal model for magainin were reviewed. The properties of melittin pores are closely similar to that of magainin but unlike that of alamethicin. We conclude that, among naturally produced peptides that we have investigated, only alamethicin conforms to the barrel-stave model. Other peptides, including magainins, melittin and protegrins, all appear to induce transmembrane pores that conform to the toroidal model in which the lipid monolayer bends continuously through the pore so that the water core is lined by both the peptides and the lipid headgroups.

Experimental evidence for hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

Hydrophobic matching, in which transmembrane proteins cause the surrounding lipid bilayer to adjust its hydrocarbon thickness to match the length of the hydrophobic surface of the protein, is a commonly accepted idea in membrane biophysics. To test this idea, gramicidin (gD) was embedded in 1, 2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1, 2-myristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers at the peptide/lipid molar ratio of 1:10. Circular dichroism (CD) was measured to ensure that the gramicidin was in the beta6.3 helix form. The bilayer thickness (the phosphate-to-phosphate distance, or PtP) was measured by x-ray lamellar diffraction. In the Lalpha phase near full hydration, PtP is 30.8 A for pure DLPC, 32.1 A for the DLPC/gD mixture, 35.3 A for pure DMPC, and 32.7 A for the DMPC/gD mixture. Gramicidin apparently stretches DLPC and thins DMPC toward a common thickness as expected by hydrophobic matching. Concurrently, gramicidin-gramicidin correlations were measured by x-ray in-plane scattering. In the fluid phase, the gramicidin-gramicidin nearest-neighbor separation is 26.8 A in DLPC, but shortens to 23.3 A in DMPC. These experiments confirm the conjecture that when proteins are embedded in a membrane, hydrophobic matching creates a strain field in the lipid bilayer that in turn gives rise to a membrane-mediated attractive potential between proteins.

Theoretical analysis of hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

We present a quantitative analysis of the effects of hydrophobic matching and membrane-mediated protein-protein interactions exhibited by gramicidin embedded in dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) bilayers (Harroun et al., 1999. Biophys. J. 76:937-945). Incorporating gramicidin, at 1:10 peptide/lipid molar ratio, decreases the phosphate-to-phosphate (PtP) peak separation in the DMPC bilayer from 35.3 A without gramicidin to 32.7 A. In contrast, the same molar ratio of gramicidin in DLPC increases the PtP from 30.8 A to 32.1 A. Concurrently, x-ray in-plane scattering showed that the most probable nearest-neighbor separation between gramicidin channels was 26.8 A in DLPC, but reduced to 23.3 A in DMPC. In this paper we review the idea of hydrophobic matching in which the lipid bilayer deforms to match the hydrophobic surface of the embedded proteins. We use a simple elasticity theory, including thickness compression, tension, and splay terms to describe the membrane deformation. The energy of membrane deformation is compared with the energy cost of hydrophobic mismatch. We discuss the boundary conditions between a gramicidin channel and the lipid bilayer. We used a numerical method to solve the problem of membrane deformation profile in the presence of a high density of gramicidin channels and ran computer simulations of 81 gramicidin channels to find the equilibrium distributions of the channels in the plane of the bilayer. The simulations contain four parameters: bilayer thickness compressibility 1/B, bilayer bending rigidity Kc, the channel-bilayer mismatch Do, and the slope of the interface at the lipid-protein boundary s. B, Kc, and Do were experimentally measured; the only free parameter is s. The value of s is determined by the requirement that the theory produces the experimental values of bilayer thinning by gramicidin and the shift in the peak position of the in-plane scattering due to membrane-mediated channel-channel interactions. We show that both hydrophobic matching and membrane-mediated interactions can be understood by the simple elasticity theory.

Collective dynamics in fully hydrated phospholipid bilayers studied by inelastic x-ray scattering.

The short wavelength density fluctuation of DLPC (dilaurylphosphatidylcholine) bilayers close to full hydration has been studied by the inelastic x-ray scattering technique below and above the main transition temperature. The analysis based on a generalized three effective eigenmode theory allows us to construct the dispersion relation of the high frequency sound mode for the first time. The marked softening of the excitation near k = 14 nm(-1), corresponding to the lipid chain-chain correlation peak in the structure factor, in the L(alpha) phase implies prevalent occurrences of short-wavelength in-plane motions of lipid chains that might be of importance for transportation of small molecules across membranes.

Membrane thinning effect of the beta-sheet antimicrobial protegrin.

Lipid bilayers containing the antimicrobial peptide protegrin-1 (PG-1) were studied by lamellar X-ray diffraction. Previously, we have shown that the peptide exists in two distinct states when associated with lipid bilayers depending on the peptide concentration [Heller, W. T., Waring, A. J., Lehrer, R. I., and Huang, H. W. (1998) Biochemistry 37, 17331-17338]. For concentrations below a lipid-dependent threshold, PG-1 exhibits a unique oriented circular dichroism spectrum called the S state. X-ray experiments show that in this state PG-1 decreases the thickness of the lipid bilayer in proportion to the peptide concentration, similar to alamethicin's membrane thinning effect. This indicates that the S state is adsorbed in the headgroup region of the lipid bilayer, where the peptide is in an inactive state. For PG-1 above the threshold concentration, X-ray diffraction shows that the interaction between the peptide and the bilayer changes significantly. These results suggest that PG-1 has the same concentration-gated mechanism of action as alamethicin.