Orientational Dynamics of Water in Phospholipid Bilayers with Different Hydration Levels

Dynamical Crossover of Interfacial Water upon Melting of a Lipid Bilayer: Influence of Different Parts of the Headgroups

All-atom molecular dynamics simulations of a 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer in contact with liquid water were performed at different temperatures ranging from 285 to 320 K. We have investigated the heterogeneity and dynamical transitions in interfacial water as the lipid bilayer undergoes a melting transition. Results are obtained for water at the outer surface of the bilayer and for those buried more deeply in the lipid chains of the bilayer. It is found that lipid bilayer melting influences both the structure and dynamics of interfacial water. The number of interfacial water molecules shows a jump in the melting of the bilayer. The temperature dependence of the diffusivity and orientational relaxation of interfacial water molecules exhibits a dynamical crossover upon melting of the bilayer. The extent of dynamical crossover is found to be rather strong with significant changes in activation barriers for interfacial water around the carbonyl groups, which are deeply buried toward the lipid chains of the bilayer. The dynamical crossover gradually decreases as one moves further away from the outer surface, and it essentially vanishes for water in the region of 5-10 Å from the outer surface. It is found that the lipid melting-induced dynamical crossover of interfacial water is significant only for water that is in close proximity to the bilayer surface or deeply buried into it. The current results reveal that water molecules in different parts of the interface respond differently on melting of the bilayer. The current study also shows that the carbonyl-bound water molecules can play an important role in the phase transition of the bilayer as the temperature is raised through its melting point.

The Nanometer-Scale Proximity of Bilayers Facilitates Intermembrane Lipid Transfer

Biological membranes approach one another in various biological phenomena, such as lipid transport at membrane contact sites and membrane fusion. The proximity of two bilayers may cause environmental changes in the interbilayer space and alter the dynamics of lipid molecules. Here, we investigate the structure and dynamics of vesicles aggregated due to the depletion attraction caused by polyethylene glycol (PEG) through static and dynamic small-angle neutron scattering. Manipulation of the interbilayer distance using PEG-conjugated lipids reveals that lipid molecules rapidly transfer between vesicles when the opposing bilayers are within ∼2 nm of each other. This distance corresponds to a region in which water molecules are more structured than in bulk water. Kinetic analysis suggests that the decrease in water entropy is responsible for the progression of lipid transfer. These results provide a basis for understanding the dynamic function of biomembranes in confined regions.

How Far Is "Bulk Water" from Interfaces? Depends on the Nature of the Surface and What We Measure

Using systematic molecular dynamics (MD) simulations, we revisit the question: At what distance from an interface do the properties of "bulk water" get recovered? We have considered three different kinds of interfaces: nonpolar (hydrophobic; isooctane-water interface), charged (negative; AOT bilayer), and polar (zwitterionic; POPC bilayer). In order to interrogate the extent of perturbation of the interfacial water molecules as a function of the distance from the interface, we utilize a diverse range of structural and dynamical parameters. To capture the structural perturbations, we look into local density (translational order), local tetrahedral order parameter, and dipolar orientation of the water molecules. We also explore the anisotropic diffusion of the water molecules in the direction perpendicular to the interface as well as the planar diffusion parallel to the interface in a distance dependent manner. In addition, the orientational time correlation functions have been computed to understand the extent of slowdown in the rotational dynamics. As expected, the electrostatic field emanating from the charged AOT interface seems to have the highest long-range effect on the orientational order and dynamics of the water molecules, whereas specific interactions like hydrogen bonding and electrostatic interaction lead to significant trapping and kinetic slowdown for both AOT and POPC (zwitterionic) very close to the interface. Our analysis highlights that not only the length-scale of perturbation depends on the nature of the interfaces and specific interactions but also the type of water property that we measure/calculate. Different water properties seem to have widely different length-scale of perturbation. Orientational order parameters seem to be perturbed to a much longer length-scale as compared to translational order parameters. The global orientational order of water can be perturbed even up to ∼4-5 nm near the negatively charged AOT surface in the absence of any extra electrolyte. This observation has significant implication toward the interpretation of experimental measurements as well since different spectroscopic techniques would probe different parameters or water properties with possible mutual disagreement and inconsistency between different types of measurements. Thus, our study provides a broader and unifying perspective toward the aspect of "context dependent" structural and dynamical perturbation of "interfacial water".

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

How stereochemistry of lipid components can affect lipid organization and the route of liposome internalization into cells

Though liposome-based drugs are in clinical use, the mechanism of cell internalization of liposomes is yet an object of controversy. The present experimental investigation, carried out on human glioblastoma cells, indicated different internalization routes for two diastereomeric liposomes. Molecular dynamics simulations of the lipid bilayers of the two formulations indicated that the different stereochemistry of a lipid component controls some parameters such as area per lipid molecule and fluidity of lipid membranes, surface potential and water organization at the lipid/water interface, all of which affect the interaction with biomolecules and cell components.

Redefining the concept of hydration water near soft interfaces

Water determines the properties of biological systems. Therefore, understanding the nature of the mutual interaction between water and biosystems is of primary importance for a proper assessment of any biological activity, e.g., the efficacy of new drugs or vaccines. A convenient way to characterize the interactions between biosystems and water is to analyze their impact on water density and dynamics in the proximity of the interfaces. It is commonly accepted that water bulk density and dynamical properties are recovered at distances of the order of 1 nm away from the surface of biological systems. This notion leads to the definition of hydration or biological water as the nanoscopic layer of water covering the surface of biosystems and to the expectation that all the effects of the water-interface interaction are limited to this thin region. Here, we review some of our latest contributions, showing that phospholipid membranes affect the water dynamics, structural properties, and hydrogen bond network at a distance that is more than twice as large as the commonly evoked ∼1nm thick layer and of the order of 2.4 nm. Furthermore, we unveil that at a shorter distance ∼0.5nm from the membrane, instead, there is an additional interface between lipid-bound and unbound water. Bound water has a structural role in the stability of the membrane. Our results imply that the concept of hydration water should be revised or extended and pave the way to a deeper understanding of the mutual interactions between water and biological systems.

Network Topology in Water Nanoconfined between Phospholipid Membranes

Water provides the driving force for the assembly and stability of many cellular components. Despite its impact on biological functions, a nanoscale understanding of the relationship between its structure and dynamics under soft confinement has remained elusive. As expected, water in contact with biological membranes recovers its bulk density and dynamics at ∼1 nm from phospholipid headgroups but surprisingly enhances its intermediate range order (IRO) over a distance, at least, twice as large. Here, we explore how the IRO is related to the water's hydrogen-bond network (HBN) and its coordination defects. We characterize the increased IRO by an alteration of the HBN up to more than eight coordination shells of hydration water. The HBN analysis emphasizes the existence of a bound-unbound water interface at ∼0.8 nm from the membrane. The unbound water has a distribution of defects intermediate between bound and bulk water, but with density and dynamics similar to bulk, while bound water has reduced thermal energy and many more HBN defects than low-temperature water. This observation could be fundamental for developing nanoscale models of biological interactions and for understanding how alteration of the water structure and topology, for example, due to changes in extracellular ions concentration, could affect diseases and signaling. More generally, it gives us a different perspective to study nanoconfined water.

Composition Dependence of Water Permeation Across Multicomponent Gel-Phase Bilayers

The permeability of multicomponent phospholipid bilayers in the gel phase is investigated via molecular dynamics simulation. The physical role of the different molecules is probed by comparing multiple mixed-component bilayers containing distearylphosphatidylcholine (DSPC) with varying amounts of either the emollient isostearyl isostearate or long-chain alcohol (dodecanol, octadecanol, or tetracosanol) molecules. Permeability is found to depend on both the tail packing density and hydrogen bonding between lipid headgroups and water. Whereas the addition of emollient or alcohol molecules to a gel-phase DSPC bilayer can increase the tail packing density, it also disturbed the hydrogen-bonding network, which in turn can increase interfacial water dynamics. These phenomena have opposing effects on bilayer permeability, which is found to depend on the balance between enhanced tail packing and decreased hydrogen bonding.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Perspective: Structure and ultrafast dynamics of biomolecular hydration shells

The structure and function of biomolecules can be strongly influenced by their hydration shells. A key challenge is thus to determine the extent to which these shells differ from bulk water, since the structural fluctuations and molecular excitations of hydrating water molecules within these shells can cover a broad range in both space and time. Recent progress in theory, molecular dynamics simulations, and ultrafast vibrational spectroscopy has led to new and detailed insight into the fluctuations of water structure, elementary water motions, and electric fields at hydrated biointerfaces. Here, we discuss some central aspects of these advances, focusing on elementary molecular mechanisms and processes of hydration on a femto- to picosecond time scale, with some special attention given to several issues subject to debate.

Dynamics of Water in Gemini Surfactant-Based Lyotropic Liquid Crystals

The dynamics of water confined to nanometer-sized domains is important in a variety of applications ranging from proton exchange membranes to crowding effects in biophysics. In this work, we study the dynamics of water in gemini surfactant-based lyotropic liquid crystals (LLCs) using molecular dynamics simulations. These systems have well characterized morphologies, for example, hexagonal, gyroid, and lamellar, and the surfaces of the confining regions can be controlled by modifying the headgroup of the surfactants. This allows one to study the effect of topology, functionalization, and interfacial curvature on the dynamics of confined water. Through analysis of the translational diffusion and rotational relaxation, we conclude that the hydration level and resulting confinement length scale is the predominate determiner of the rates of water dynamics, and other effects, namely, surface functionality and curvature, are largely secondary. This novel analysis of the water dynamics in these LLC systems provides an important comparison for previous studies of water dynamics in lipid bilayers and reverse micelles.

Structural Interpretation of the Large Slowdown of Water Dynamics at Stacked Phospholipid Membranes for Decreasing Hydration Level: All-Atom Molecular Dynamics

Hydration water determines the stability and function of phospholipid membranes as well as the interaction of membranes with other molecules. Experiments and simulations have shown that water dynamics slows down dramatically as the hydration decreases, suggesting that the interfacial water that dominates the average dynamics at low hydration is slower than water away from the membrane. Here, based on all-atom molecular dynamics simulations, we provide an interpretation of the slowdown of interfacial water in terms of the structure and dynamics of water–water and water–lipid hydrogen bonds (HBs). We calculate the rotational and translational slowdown of the dynamics of water confined in stacked phospholipid membranes at different levels of hydration, from completely hydrated to poorly hydrated membranes. For all hydrations, we analyze the distribution of HBs and find that water–lipids HBs last longer than water–water HBs and that at low hydration most of the water is in the interior of the membrane. We also show that water–water HBs become more persistent as the hydration is lowered. We attribute this effect (i) to HBs between water molecules that form, in turn, persistent HBs with lipids; (ii) to the hindering of the H-bonding switching between water molecules due to the lower water density at the interface; and (iii) to the higher probability of water–lipid HBs as the hydration decreases. Our interpretation of the large dynamic slowdown in water under dehydration is potentially relevant in understanding membrane biophysics at different hydration levels.

Anomalous surface diffusion of protons on lipid membranes

The cellular energy machinery depends on the presence and properties of protons at or in the vicinity of lipid membranes. To asses the energetics and mobility of a proton near a membrane, we simulated an excess proton near a solvated DMPC bilayer at 323 K, using a recently developed method to include the Grotthuss proton shuttling mechanism in classical molecular dynamics simulations. We obtained a proton surface affinity of -13.0 ± 0.5 kJ mol(-1). The proton interacted strongly with both lipid headgroup and linker carbonyl oxygens. Furthermore, the surface diffusion of the proton was anomalous, with a subdiffusive regime over the first few nanoseconds, followed by a superdiffusive regime. The time- and distance dependence of the proton surface diffusion coefficient within these regimes may also resolve discrepancies between previously reported diffusion coefficients. Our simulations show that the proton anomalous surface diffusion originates from restricted diffusion in two different surface-bound states, interrupted by the occasional bulk-mediated long-range surface diffusion. Although only a DMPC membrane was considered in this work, we speculate that the restrictive character of the on-surface diffusion is highly sensitive to the specific membrane conditions, which can alter the relative contributions of the surface and bulk pathways to the overall diffusion process. Finally, we discuss the implications of our findings for the energy machinery.

Effects of surface pressure on the properties of Langmuir monolayers and interfacial water at the air-water interface

The effects of surface pressure on the physical properties of Langmuir monolayers of palmitic acid (PA) and dipalmitoylphosphatidic acid (DPPA) at the air/water interface are investigated through molecular dynamics simulations with atomistic force fields. The structure and dynamics of both monolayers and interfacial water are compared across the range of surface pressures at which stable monolayers can form. For PA monolayers at T = 300 K, the untilted condensed phase with a hexagonal lattice structure is found at high surface pressure, while the uniformly tilted condensed phase with a centered rectangular lattice structure is observed at low surface pressure, in agreement with the available experimental data. A state with uniform chain tilt but no periodic spatial ordering is observed for DPPA monolayers on a Na(+)/water subphase at both high and low surface pressures. The hydrophobic acyl chains of both monolayers pack efficiently at all surface pressures, resulting in a very small number of gauche defects. The analysis of the hydrogen-bonding structure/dynamics at the monolayer/water interface indicates that water molecules hydrogen-bonded to the DPPA head groups reorient more slowly than those hydrogen-bonded to the PA head groups, with the orientational dynamics becoming significantly slower at high surface pressure. Possible implications for physicochemical processes taking place on marine aerosols in the atmosphere are discussed.

Coupled diffusion in lipid bilayers upon close approach

Biomembrane interfaces create regions of slowed water dynamics in their vicinity. When two lipid bilayers come together, this effect is further accentuated, and the associated slowdown can affect the dynamics of larger-scale processes such as membrane fusion. We have used molecular dynamics simulations to examine how lipid and water dynamics are affected as two lipid bilayers approach each other. These two interacting fluid systems, lipid and water, both slow and become coupled when the lipid membranes are separated by a thin water layer. We show in particular that the water dynamics become glassy, and diffusion of lipids in the apposed leaflets becomes coupled across the water layer, while the "outer" leaflets remain unaffected. This dynamic coupling between bilayers appears mediated by lipid-water-lipid hydrogen bonding, as it occurs at bilayer separations where water-lipid hydrogen bonds become more common than water-water hydrogen bonds. We further show that such coupling occurs in simulations of vesicle-vesicle fusion prior to the fusion event itself. Such altered dynamics at membrane-membrane interfaces may both stabilize the interfacial contact and slow fusion stalk formation within the interface region.

Long-Range Lipid-Water Interaction as Observed by ATR-FTIR Spectroscopy

It is commonly assumed that the structure of water at a lipid-water interface is influenced mostly in the first hydration layer. However, recent results from different experimental methods show that perturbation extends through several hydration layers. Due to its low light penetration depth, attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy is specifically suited to study interlamellar water structure in multibilayers. Results obtained by this technique confirm the long-range water structure disturbance. Consequently, in confined membrane environments nearly all water molecules can be perturbed. It is important to note that the behavior of confined water molecules differs significantly in samples prepared in excess water and in partially hydrated samples. We show in what manner the interlamellar water perturbation is influenced by the hydration level and how it is sequentially modified with a step-by-step dehydration of samples either by water evaporation or by osmotic pressure. Our results also indicate that besides different levels of hydration the lipid-water interaction is modulated by different lipid headgroups and different lipid phases as well. Therefore, modification of interlamellar water properties may clarify the role of water-mediated effects in biological processes.

Mosaic of Water Orientation Structures at a Neutral Zwitterionic Lipid/Water Interface Revealed by Molecular Dynamics Simulations

Ordering of water structures near the surface of biological membranes has been recently extensively studied using interface-selective techniques like vibrational sum frequency generation (VSFG) spectroscopy. The detailed structures of interface water have emerged for charged lipids, but those for neutral zwitterionic lipids remain obscure. We analyze an all-atom molecular dynamics (MD) trajectory of a hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer to characterize the orientation of interface waters in different chemical environments. The structure and dynamics of interfacial waters strongly depend on both their vertical position along the bilayer normal as well as vicinal lipid charged groups. Water orientation in the vicinity of phosphate groups is opposite to that around choline groups. The results are consistent with observed VSFG spectra and demonstrate that a mosaic of water orientation structures exists on the surface of a neutral zwitterionic phospholipid bilayer, reflecting rapid water exchange and the influence of local chemical environments.

Effects of millimeter wave irradiation and equivalent thermal heating on the activity of individual neurons in the leech ganglion

Many of today's radiofrequency-emitting devices in telecommunication, telemedicine, transportation safety, and security/military applications use the millimeter wave (MMW) band (30–300 GHz). To evaluate the biological safety and possible applications of this radiofrequency band for neuroscience and neurology, we have investigated the physiological effects of low-intensity 60-GHz electromagnetic irradiation on individual neurons in the leech midbody ganglia. We applied incident power densities of 1, 2, and 4 mW/cm² to the whole ganglion for a period of 1 min while recording the action potential with a standard sharp electrode electrophysiology setup. For comparison, the recognized U.S. safe exposure limit is 1 mW/cm² for 6 min. During the exposure to MMWs and gradual bath heating at a rate of 0.04°C/s (2.4°C/min), the ganglionic neurons exhibited similar dose-dependent hyperpolarization of the plasma membrane and decrease in the action potential amplitude. However, narrowing of the action potential half-width during MMW irradiation at 4 mW/cm² was 5 times more pronounced compared with that during equivalent bath heating of 0.6°C. Even more dramatic difference in the effects of MMW irradiation and bath heating was noted in the firing rate, which was suppressed at all applied MMW power densities and increased in a dose-dependent manner during gradual bath heating. The mechanism of enhanced narrowing of action potentials and suppressed firing by MMW irradiation, compared with that by gradual bath heating, is hypothesized to involve specific coupling of MMW energy with the neuronal plasma membrane.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. III. Water clustering and vibrational energy transfer

Water clustering and connectivity around lipid bilayers strongly influences the properties of membranes and is important for functions such as proton and ion transport. Vibrational anisotropic pump-probe spectroscopy is a powerful tool for understanding such clustering, as the measured anisotropy depends upon the time-scale and degree of intra- and intermolecular vibrational energy transfer. In this article, we use molecular dynamics simulations and theoretical vibrational spectroscopy to help interpret recent experimental measurements of the anisotropy of water in lipid multi-bilayers as a function of both lipid hydration level and isotopic substitution. Our calculations are in satisfactory agreement with the experiments of Piatkowski, Heij, and Bakker, and from our simulations we can directly probe water clustering and connectivity. We find that at low hydration levels, many water molecules are in fact isolated, although up to 70% of hydration water forms small water clusters or chains. At intermediate hydration levels, water forms a wide range of cluster sizes, while at higher hydration levels, the majority of water molecules are part of a large, percolating water cluster. Therefore, the size, number, and nature of water clusters are strongly dependent on lipid hydration level, and the measured anisotropy reflects this through its dependence on intermolecular energy transfer.

Structural Dynamics of Hydrated Phospholipid Surfaces Probed by Ultrafast 2D Spectroscopy of Phosphate Vibrations

The properties of biomembranes depend in a decisive way on interactions of phospholipids with hydrating water molecules. To map structural dynamics of a phospholipid-water interface on the length and time scale of molecular motions, we introduce the phospholipid symmetric and asymmetric phosphate stretch vibrations as probes of interfacial hydrogen bonds and electrostatic interactions. The first two-dimensional infrared spectra of such modes and a line shape analysis by density matrix theory reveal two distinct structural dynamics components; the first 300 fs contribution is related to spatial fluctuations of charged phospholipid head groups with additional water contributions at high hydration levels; the second accounts for water-phosphate hydrogen bonds persisting longer than 10 ps. Our results reveal a relatively rigid hydration shell around phosphate groups, a behavior relevant for numerous biomolecular systems.

Diffusion of water and selected atoms in DMPC lipid bilayer membranes

Molecular dynamics simulations have been used to determine the diffusion of water molecules as a function of their position in a fully hydrated freestanding 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) bilayer membrane at 303 K and 1 atm. The diffusion rate of water in a ∼10 Å thick layer just outside the membrane surface is reduced on average by a factor of ∼2 relative to bulk. For water molecules penetrating deeper into the membrane, there is an increasing reduction in the average diffusion rate with up to one order of magnitude decrease for those deepest in the membrane. A comparison with the diffusion rate of selected atoms in the lipid molecules shows that ∼6 water molecules per lipid molecule move on the same time scale as the lipids and may therefore be considered to be tightly bound to them. The quasielastic neutron scattering functions for water and selected atoms in the lipid molecule have been simulated and compared to observed quasielastic neutron scattering spectra from single-supported bilayer DMPC membranes.

Power-law trapping of water molecules on the lipid-membrane surface induces water retardation

Cell membranes provide unique local environments for biological reactions, where the diffusion of biomolecules as well as water molecules plays critical roles. Translational and rotational motions of water molecules near membranes are known to be slower than those in bulk. Using all-atom molecular dynamics simulations of a membrane, we show that the temperature dependence of the water molecular motions on the membrane surface is different from that in bulk. Decreasing temperature enhances the water retardation on the membrane surface, and the lateral motions of water molecules are correlated with the vertical motions. We find that trapping times of water molecules onto membrane surfaces are distributed according to a power-law distribution and that the power-law exponents depend linearly on temperature, suggesting a random energy landscape picture. Moreover, we find that water molecules on the membrane surfaces exhibit subdiffusions in translational motions.

Transition of the hydration state of a surfactant accompanying structural transitions of self-assembled aggregates

What role does water play in the self-assembly of soft materials? To understand the correlation between the hydration state and the various self-assembled structures of a nonionic surfactant, terahertz time-domain spectroscopy has been performed for a C(12)E(5) solution with complementary use of small-angle x-ray scattering. Precise observations of the hydration state show clearly that transitions of the hydration state are accompanied by structural phase transitions of the surfactant from hexagonal to micelle to lamellae. These transitions of hydration state suggest that water is not a homogeneous solvent, and the interaction between water and the soft material is important for self-assembly.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. II. Two-dimensional infrared and peak shift observables within different theoretical approximations

In a previous report, we calculated the infrared absorption spectrum and both the isotropic and anisotropic pump-probe signals for the OD stretch of isotopically dilute water in dilauroylphosphatidylcholine (DLPC) multi-bilayers as a function of the lipid hydration level. These results were then compared to recent experimental measurements and are in generally good agreement. In this paper, we will further investigate the structure and dynamics of hydration water using molecular dynamics simulations and calculations of the two-dimensional infrared and vibrational echo peak shift observables for hydration water in DLPC membranes. These observables have not yet been measured experimentally, but future comparisons may provide insight into spectral diffusion processes and hydration water heterogeneity. We find that at low hydration levels the motion of water molecules inside the lipid membrane is significantly arrested, resulting in very slow spectral diffusion. At higher hydration levels, spectral diffusion is more rapid, but still slower than in bulk water. We also investigate the effects of several common approximations on the calculation of spectroscopic observables by computing these observables within multiple levels of theory. The impact of these approximations on the resulting spectra affects our interpretation of these measurements and reveals that, for example, the cumulant approximation, which may be valid for certain systems, is not a good approximation for a highly heterogeneous environment such as hydration water in lipid multi-bilayers.

Vibrational spectroscopy of water at interfaces

Understanding liquid water's behavior at the molecular level is essential to progress in fields as disparate as biology and atmospheric sciences. Moreover, the properties of water in bulk and water at interfaces can be very different, making the study of the hydrogen-bonding networks therein very important. With recent experimental advances in vibrational spectroscopy, such as ultrafast pulses and heterodyne detection, it is now possible to probe the structure and dynamics of bulk and interfacial water in unprecedented detail. We consider here three aqueous interfaces: the water liquid-vapor interface, the interface between water and the surfactant headgroups of reverse micelles, and the interface between water and the lipid headgroups of aligned multi-bilayers. In the first case, sum-frequency spectroscopy is used to probe the interface. In the second and third cases, the confined water pools are sufficiently small that techniques of bulk spectroscopy (such as FTIR, pump-probe, two-dimensional IR, and the like) can be used to probe the interfacial water. In this Account, we discuss our attempts to model these three systems and interpret the existing experiments. For the water liquid-vapor interface, we find that three-body interactions are essential for reproducing the experimental sum-frequency spectrum, and presumably for the structure of the interface as well. The observed spectrum is interpreted as arising from overlapping and canceling positive and negative contributions from molecules in different hydrogen-bonding environments. For the reverse micelles, our theoretical models confirm that the experimentally observed blue shift of the water OD stretch (for dilute HOD in H(2)O) arises from weaker hydrogen bonding to sulfonate oxygens. We interpret the observed slow-down in water rotational dynamics as arising from curvature-induced frustration. For the water confined between lipid bilayers, our theoretical models confirm that the experimentally observed red shift of the water OD stretch arises from stronger hydrogen bonding to phosphate oxygens. We develop a model for heterogeneous vibrational lifetime distributions, and we implement the model to calculate isotropic and anisotropic pump-probe decays. We then compare these results with experimental data. Clearly, recent experimental advances in vibrational spectroscopy have led to beautiful new results, providing information about the structure and dynamics of water at interfaces. These experimental and concomitant theoretical advances (particularly the unified theoretical framework of non-linear response functions) have greatly contributed to our understanding of this unique and important substance.

Aqueous solutions at the interface with phospholipid bilayers

In a sense, life is defined by membranes, because they delineate the barrier between the living cell and its surroundings. Membranes are also essential for regulating the machinery of life throughout many interfaces within the cell's interior. A large number of experimental, computational, and theoretical studies have demonstrated how the properties of water and ionic aqueous solutions change due to the vicinity of membranes and, in turn, how the properties of membranes depend on the presence of aqueous solutions. Consequently, understanding the character of aqueous solutions at their interface with biological membranes is critical to research progress on many fronts. The importance of incorporating a molecular-level description of water into the study of biomembrane surfaces was demonstrated by an examination of the interaction between phospholipid bilayers that can serve as model biological membranes. The results showed that, in addition to well-known forces, such as van der Waals and screened Coulomb, one has to consider a repulsion force due to the removal of water between surfaces. It was also known that physicochemical properties of biological membranes are strongly influenced by the specific character of the ions in the surrounding aqueous solutions because of the observation that different anions produce different effects on muscle twitch tension. In this Account, we describe the interaction of pure water, and also of aqueous ionic solutions, with model membranes. We show that a symbiosis of experimental and computational work over the past few years has resulted in substantial progress in the field. We now better understand the origin of the hydration force, the structural properties of water at the interface with phospholipid bilayers, and the influence of phospholipid headgroups on the dynamics of water. We also improved our knowledge of the ion-specific effect, which is observed at the interface of the phospholipid bilayer and aqueous solution, and its connection with the Hofmeister series. Nevertheless, despite substantial progress, many issues remain unresolved. Thus, for example, we still cannot satisfactorily explain the force of interaction between phospholipid bilayers immersed in aqueous solutions of NaI. Although we try to address many issues here, the scope of the discussion is limited and does not cover such important topics as the influence of ionic solutions on phases of bilayers, the influence of salts on the properties of Langmuir monolayers containing lipid molecules, or the influence of aqueous solutions on bilayers containing mixtures of lipids. We anticipate that the future application of more powerful experimental techniques, in combination with more advanced computational hardware, software, and theory, will produce molecular-level information about these important topics and, more broadly, will further illuminate our understanding of interfaces between aqueous solutions and biological membranes.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. I. Infrared spectra and ultrafast pump-probe observables

The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.

Long-range hydration effect of lipid membrane studied by terahertz time-domain spectroscopy

The hydration state of biomolecules is believed to affect their self-assembly. The hydration state of phospholipid bilayers is studied precisely by terahertz spectroscopy, by which water perturbed by a lipid membrane is detected sensitively from the observation of the relaxation dynamics of water molecules in the subpicosecond time scale. Combined with x-ray observation of the lamellar structure of the lipid, a long-range hydration effect on up to 4-5 layers of water is confirmed. Most water molecules in the lamellae fall into the hydration water, and condensation of them is also indicated.

Orientational dynamics for an amphiphilic-solvent solution

In this work, we performed Monte Carlo simulations on a lattice model for spontaneous amphiphilic aggregation, in order to study the orientational and hydrogen-bonding dynamics of water on different regions inside the micellar solution. We employed an associating lattice gas model that mimics the aqueous solvent, which presents a rich phase diagram with first- and second-order transition lines. Even though this is a simplified model, it makes possible to investigate the orientational dynamics of water in an equilibrium solution of amphiphiles, as well as the influence of the different phases of the solvent in the interfacial and bulk water dynamics. By means of extensive simulations, we showed that, at high temperatures, the behavior of the orientational relaxation and hydrogen bonding of water molecules in the bulk, first, and second hydration shells are considerable different. We observe the appearance of a very slow component for water molecules in the first hydration shell of micelles when the system reaches a high-density phase, consistent with previous theoretical and experimental studies concerning biological water. Also, at high temperatures, we find that water molecules in the second hydration shell of micelles have an orientational decay similar to that of bulk water, but with a generally slower dynamics. Otherwise, at low temperatures, we have two components for the orientational relaxation of bulk water in the low density liquid phase, and only a single component in the high density liquid (HDL) phase, which reflect the symmetry properties of the different phases of the solvent model. In the very dense region of water molecules in the first hydration shell of micelles at low temperatures, we find two components for the orientational relaxation on both liquid phases, one of them much slower than that in the single component of bulk water in the HDL phase. This happens even though our model does not present any hindrance to the water rotational freedom caused by the presence of the amphiphiles.