Structure and dynamics of water at the interface with phospholipid bilayers

Influence of cholesterol on hydrogen-bond dynamics of water molecules in lipid-bilayer systems at varying temperatures

Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting a considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoylphosphatidylcholine (DPPC) and palmitoyl sphingomyelin; they are distinct with respect to the structures of the hydrogen-bond (H-bond) acceptors. Our investigation focuses on the dynamic properties and H-bonds of water molecules in the lipid-membrane systems, with a particular emphasis on the influence of Chol at varying temperatures. Notably, in the gel phase at 303 K, the presence of Chol extends the lifetimes of H-bonds of the oxygen atoms acting as H-bond acceptors within DPPC with water molecules by a factor of 1.5-2.5. In the liquid-crystalline phase at 323 K, on the other hand, H-bonding dynamics with lipid membranes remain largely unaffected by Chol. This observed shift in H-bonding states serves as a crucial key to unraveling the subtle control mechanisms governing water dynamics in lipid-membrane systems.

Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.

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 Design and Optimization of Ceramide NP-Loaded Liposomes to Restore the Skin Barrier

The impairment of skin integrity derived from derangement of the orthorhombic lateral organization is mainly caused by dysregulation of ceramide amounts in the skin barrier. Ceramides, fatty acids, and cholesterol-containing nano-based formulations have been used to impair the skin barrier. However, there is still a challenge to formulate novel formulations consisting of ceramides due to their chemical structure, poor aqueous solubility, and high molecular weight. In this study, the design and optimization of Ceramide 3 (CER-NP)-loaded liposomes are implemented based on response surface methodology (RSM). The optimum CER-NP-loaded liposome was selected based on its particle size (PS) and polydispersity index (PDI). The optimum CER-NP-loaded liposome was imagined by observing the encapsulation by using a confocal laser scanning microscope (CLSM) within fluorescently labeled CER-NP. The characteristic liquid crystalline phase and lipid chain conformation of CER-NP-loaded liposomes were determined using attenuated total reflectance infrared spectroscopy (ATR-IR). The CER-NP-loaded liposomes were imagined using a field emission scanning electron microscope (FE-SEM). Finally, the in vitro release of CER-NP from liposomes was examined using modified Franz Cells. The experimental and predicted results were well correlated. The CLSM images of optimized liposomes were conformable with the other studies, and the encapsulation efficiency of CER-NP was 93.84 ± 0.87%. ATR-IR analysis supported the characteristics of the CER-NP-loaded liposome. In addition, the lipid chain conformation shows similarity with skin barrier lipid organization. The release pattern of CER-NP liposomes was fitted with the Korsmeyer-Peppas model. The cytotoxicity studies carried out on HaCaT keratinocytes supported the idea that the liposomes for topical administration of CER-NP could be considered relatively safe. In conclusion, the optimized CER-NP-loaded liposomes could have the potential to restore the skin barrier function.

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.

Interdependence of cholesterol distribution and conformational order in lipid bilayers

We show, via molecular simulations, that not only does cholesterol induce a lipid order, but the lipid order also enhances cholesterol localization within the lipid leaflets. Therefore, there is a strong interdependence between these two phenomena. In the ordered phase, cholesterol molecules are predominantly present in the bilayer leaflets and orient themselves parallel to the bilayer normal. In the disordered phase, cholesterol molecules are mainly present near the center of the bilayer at the midplane region and are oriented orthogonal to the bilayer normal. At the melting temperature of the lipid bilayers, cholesterol concentration in the leaflets and the bilayer midplane is equal. This result suggests that the localization of cholesterol in the lipid bilayers is mainly dictated by the degree of ordering of the lipid bilayer. We validate our findings on 18 different lipid bilayer systems, obtained from three different phospholipid bilayers with varying concentrations of cholesterol. To cover a large temperature range in simulations, we employ the Dry Martini force field. We demonstrate that the Dry and the Wet Martini (with polarizable water) force fields produce comparable results.

Integrating In Vitro and In Silico Analysis of a Cationic Antimicrobial Peptide Interaction with Model Membranes of Colistin-Resistant Pseudomonas aeruginosa Strains

Bacterial antibiotic resistance is a serious global public health concern. Infections caused by colistin-resistant Pseudomonas aeruginosa (CRPa) strains represent a serious threat due to their considerable morbidity and mortality rates, since most of the current empirical antibiotic therapies are ineffective against these strains. Accordingly, cationic antimicrobial peptides (CAMPs) have emerged as promising alternatives to control resistant bacteria. In this study, the interaction of a CAMP derived from cecropin D-like (∆M2) with model membranes mimicking bacterial biomembranes of wild-type (WTPa) strains of P. aeruginosa and CRPa was evaluated through in vitro and in silico approaches. In vitro interaction was determined by infrared spectroscopy, whereas in silico molecular dynamics was performed to predict specific interactions between amino acids of ∆M2 and lipids of model membrane systems. Experimental analysis showed this peptide interacted with the lipids of bacterial-like model membranes of WTPa and CRPa. In both cases, an increase in the concentration of peptides induced an increase in the phase transition temperature of the lipid systems. On the other hand, the peptides in solution underwent a transition from a random to a helical secondary structure after interacting with the membranes mostly favored in the CRPa system. The α-helix structure percentage for ΔM2 interacting with WTPa and CRPa lipid systems was 6.4 and 33.2%, respectively. Finally, molecular dynamics showed ∆M2 to have the most affinities toward the phospholipids palmitoyl-oleyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) that mimic membranes of WTPa and CRPa, respectively. This work provides clues for elucidating the membrane-associated mechanism of action of ∆M2 against colistin-susceptible and -resistant strains of Pseudomonas aeruginosa.

Algorithm to catalogue topologies of dynamic lipid hydrogen-bond networks

Lipid membrane interfaces host reactions essential for the functioning of cells. The hydrogen-bonding environment at the membrane interface is particularly important for binding of proteins, drug molecules, and ions. We present here the implementation and applications of a depth-first search algorithm that analyzes dynamic lipid interaction networks. Lipid hydrogen-bond networks sampled transiently during simulations of lipid bilayers are clustered according to main types of topologies that characterize three-dimensional arrangements of lipids connected to each other via short water bridges. We characterize the dynamics of hydrogen-bonded lipid clusters in simulations of model POPE and POPE:POPG membranes that are often used for bacterial membrane proteins, in a model of the Escherichia coli membrane with six different lipid types, and in POPS membranes. We find that all lipids sample dynamic hydrogen-bonded networks with linear, star, or circular arrangements of the lipid headgroups, and larger networks with combinations of these three types of topologies. Overall, linear lipid-water bridges tend to be short. Water-mediated lipid clusters in all membranes with PE lipids tend to be somewhat small, with about four lipids in all membranes studied here. POPS membranes allow circular arrangements of three POPS lipids to be sampled frequently, and complex arrangements of linear, star, and circular paths may also be sampled. These findings suggest a molecular picture of the membrane interface whereby lipid molecules transiently connect in clusters with somewhat small spatial extension.

Influence of the extent of hydrophobicity on water organization and dynamics on 2D graphene oxide surfaces

Graphene oxide (GO) nanomaterials are being extensively explored for a wide spectrum of applications, ranging from water desalination to fuel cell applications due to their tunable mechanical, thermal, and electrical...

Influence of water models on water movement through AQP1

Water diffusion through membrane proteins is a key aspect of cellular function. Essential processes of cellular metabolism are driven by osmotic pressure, which depends on water channels. Membrane proteins such as aquaporins (AQPs) are responsible for enabling water permeation through the cell membrane. AQPs are highly selective, allowing only water and relatively small polar molecules to cross the membrane. Experimentally, estimation of water flux through membrane proteins is still a challenge, and hence, accurate simulations of water permeation are of particular importance. We present a numerical study of water diffusion through AQP1 comparing three water models: TIP3P, OPC, and TIP4P/2005. Bulk diffusion, diffusion permeability, and osmotic permeability are computed and compared among all models. The results show that there are significant differences between TIP3P (a particularly widespread model for simulations of biological systems) and the more recently developed TIP4P/2005 and OPC models. We demonstrate that OPC and TIP4P/2005 reproduce protein-water interactions and dynamics in very good agreement with experimental data. From this study, we find that the choice of the water model has a significant effect on the computed water dynamics as well as its molecular behavior within a biological nanopore.

Phenolic compounds alter the ion permeability of phospholipid bilayers via specific lipid interactions

This study aims to understand the role of specific phenolic-lipid interactions in the membrane-altering properties of phenolic compounds. We combine tethered lipid bilayer (tBLM) electrical impedance spectroscopy (EIS) with all-atom molecular dynamics (MD) simulations to study the membrane interactions of six phenolic compounds: caffeic acid methyl ester, caffeic acid, 3,4 dihydroxybenzoic acid, chlorogenic acid, syringic acid and p-coumaric acid. tBLM/EIS experiments showed that caffeic acid methyl ester, caffeic acid and 3,4 dihydroxybenzoic acid significantly increase the permeability of phospholipid bilayers to Na⁺ ions. In contrast, chlorogenic acid, syringic acid and p-coumaric acid showed no effect. Experiments with lipids lacking the phosphate group show a significant decrease in the membrane-altering effects indicating that specific phenolic-lipid interactions are critical in altering ion permeability. MD simulations confirm that compounds that alter ion permeability form stable interactions with the phosphate oxygen. In contrast, inactive phenolic compounds are superficially bound to the membrane surface and primarily interact with interfacial water. Our combined results show that compounds with similar structures can have very different effects on ion permeability in membranes. These effects are governed by specific interactions at the water-lipid interface and show no correlation with lipophilicity. Furthermore, none of the compounds alter the overall structure of the phospholipid bilayer as determined by area per lipid and order parameters. Based on data from this study and previous findings, we propose that phenolic compounds can alter membrane ion permeability by causing local changes in lipid packing that subsequently reduce the energy barrier for ion-induced pores.

Hydration Layer of Only a Few Molecules Controls Lipid Mobility in Biomimetic Membranes

Self-assembly of biomembranes results from the intricate interactions between water and the lipids’ hydrophilic head groups. Therefore, the lipid–water interplay strongly contributes to modulating membrane architecture, lipid diffusion, and chemical activity. Here, we introduce a new method of obtaining dehydrated, phase-separated, supported lipid bilayers (SLBs) solely by controlling the decrease of their environment’s relative humidity. This facilitates the study of the structure and dynamics of SLBs over a wide range of hydration states. We show that the lipid domain structure of phase-separated SLBs is largely insensitive to the presence of the hydration layer. In stark contrast, lipid mobility is drastically affected by dehydration, showing a 6-fold decrease in lateral diffusion. At the same time, the diffusion activation energy increases approximately 2-fold for the dehydrated membrane. The obtained results, correlated with the hydration structure of a lipid molecule, revealed that about six to seven water molecules directly hydrating the phosphocholine moiety play a pivotal role in modulating lipid diffusion. These findings could provide deeper insights into the fundamental reactions where local dehydration occurs, for instance during cell–cell fusion, and help us better understand the survivability of anhydrobiotic organisms. Finally, the strong dependence of lipid mobility on the number of hydrating water molecules opens up an application potential for SLBs as very precise, nanoscale hydration sensors.

Abnormal Properties of Low-Dimensional Confined Water

Water molecules confined to low-dimensional spaces exhibit unusual properties compared to bulk water. For example, the alternating hydrophilic and hydrophobic nanodomains on flat silicon wafer can induce the abnormal spreading of water (contact angles near 0°) which is caused by the 2D capillary effect. Hence, exploring the physicochemical properties of confined water from the nanoscale is of great value for understanding the challenges in material science and promoting the applications of nanomaterials in the fields of mass transport, nanofluidic designing, and fuel cell. The knowledge framework of confined water can also help to better understand the complex functions of the hydration layer of biomolecules, and even trace the origin of life. In this review, the physical properties, abnormal behaviors, and functions of the confined water are mainly summarized through several common low-dimensional water formats in the fields of solid/air-water interface, nanochannel confinement, and biological hydration layer. These researches indicate that the unusual behaviors of the confined water depend strongly on the confinement size and the interaction between the molecules and confining surface. These diverse properties of confined water open a new door to materials science and may play an important role in the future development of biology.

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.

Ultrafast Dynamics at the Lipid-Water Interface: DMSO Modulates H-Bond Lifetimes

Dimethyl sulfoxide (DMSO) is a common cosolvent and cryopreservation agent used to freeze cells and tissues. DMSO alters the H-bond structure of water, but its interactions with biomolecules and, specifically, with biological interfaces remain poorly understood. Here we investigate the effects of DMSO on the H-bond dynamics at the lipid-water interface using a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations. Ester carbonyl absorption spectra show that DMSO dehydrates the interface, and simulations show that the area per lipid is decreased. Ultrafast 2D IR spectra measure the time scales of frequency fluctuations at the ester carbonyl positions located precisely between the hydrophobic and hydrophilic regions of the membrane. 2D IR measurements show that low DMSO concentrations (<10 mol %) induce ∼40% faster H-bond dynamics compared with pure water, whereas increased concentrations (>10-20 mol %) once again slow down the dynamics. This slow-fast-slow trend is described in terms of two different solvation regimes. Below 10 mol %, DMSO weakens the interfacial H bond, leading to faster "bulk-like" dynamics, whereas above 10 mol %, water molecules become "relatively immobilized" as the H-bond networks becoming disrupted by the H-bond donor/acceptor imbalance at the interface. These studies are an important step toward characterizing the environments around lipid membranes, which are essential to numerous biological processes.

Interfacial structure in the liquid-liquid extraction of rare earth elements by phosphoric acid ligands: a molecular dynamics study

Solvent extraction (SX), wherein two immiscible liquids, one containing the extractant molecules and the other containing the solute to be extracted are brought in contact to effect the phase transfer of the solute, underpins metal extraction and recovery processes. The interfacial region is of utmost importance in the SX process, since besides thermodynamics, the physical and chemical heterogeneity at the interface governs the kinetics of the process. Yet, a fundamental understanding of this heterogeneity and its implications for the extraction mechanism are currently lacking. We use molecular dynamics (MD) simulations to study the liquid-liquid interface under conditions relevant to the SX of Rare Earth Elements (REEs) by a phosphoric acid ligand. Simulations revealed that the extractant molecules and varying amounts of acid and metal ions partitioned to the interface. The presence of these species had a significant effect on the interfacial thickness, hydrogen bond life times and orientations of the water molecules at the interface. Deprotonation of the ligands was essential for the adsorption of the metal ions at the interface, with these ions forming a number of different complexes at the interface involving one to three extractant molecules and four to eight water molecules. Although the interface itself was rough, no obvious 'finger-like' water protrusions penetrating the organic phase were seen in our simulations. While the results of our work help us gain fundamental insights into the sequence of events leading to the formation of a variety of interfacial complexes, they also emphasize the need to carry out a more detailed atomic level study to understand the full mechanism of extraction of REEs from the aqueous to organic phases by phosphoric acid ligands.

Perturbation of Hydrogen-Bonding Networks over Supported Lipid Bilayers by Poly(allylamine hydrochloride)

Water is vital to many biochemical processes and is necessary for driving fundamental interactions of cell membranes with their external environments, yet it is difficult to probe the membrane/water interface directly and without the use of external labels. Here, we employ vibrational sum frequency generation spectroscopy to understand the role of interfacial water molecules above bilayers formed from zwitterionic (phosphatidylcholine) and anionic (phosphatidylglycerol, PG, and phosphatidylserine, PS) lipids as they are exposed to the common polycation poly(allylamine hydrochloride) (PAH) in 100 mM NaCl. We show that as the concentration of PAH is increased, the interfacial water molecules are irreversibly displaced and find that it requires 10 times more PAH to displace interfacial water molecules from membranes formed from purely zwitterionic lipids when compared to membranes that contain the anionic PG and PS lipids. This outcome is likely due to the difference in (1) the energy with which water molecules are bound to the lipid headgroups, (2) the number of water molecules bound to the headgroups, which is related to the headgroup area, and (3) the electrostatic interactions between the PAH molecules and the negatively charged lipids that are favored when compared to the zwitterionic lipid headgroups. The findings presented here contribute to establishing causal relationships in nanotoxicology and to understanding, controlling, and predicting the initial steps that lead to the lysis of cells exposed to membrane-disrupting polycations or to transfection.

Unique Structures and Vibrational Spectra of Protic Ionic Liquids Confined in TiO2 Slits: The Role of Interfacial Hydrogen Bonds

The ionic liquid (IL)/titanium dioxide (TiO₂) interface exists in many application systems, such as nanomaterial synthesis, catalysis, and electrochemistry systems. The nanoscale interfacial properties in the above systems are a common issue. However, directly detecting the interfacial properties of nanoconfined ILs by experimental methods is still challenging. To help better learn about the interfacial issue, molecular dynamics simulations have been performed to explore the structures, vibration spectra, and hydrogen bond (HB) properties at the IL/TiO₂ interface. Ethylammonium nitrate (EAN) ILs confined in TiO₂ slit pores with different pore widths were studied. A unique vibrational spectrum appeared for EAN ILs confined in a 0.7 nm TiO₂ slit, and this phenomenon is related to interfacial hydrogen bonds (HBs). An analysis of the HB types indicated that the interfacial NH₃⁺ group of the cations was in an asymmetric HB environment in the 0.7 nm TiO₂ slit, which led to the disappearance of the symmetric N-H stretching mode. In addition, the significant increase in the HB strength between NH₃⁺ groups and the TiO₂ surface slowed down the stretching vibration of the N-H bond, resulting in one peak in the vibrational spectra at a lower frequency. For the first time, our simulation work establishes a molecular-level relationship between the vibrational spectrum and the local HB environment of nanoconfined ILs at the IL/TiO₂ interface, and this relationship is helpful for interface design in related systems.

Characterization of the interface of binary mixed DOPC:DOPS liposomes in water: The impact of charge condensation

Solutions of liposomes composed of binary mixtures of anionic dioleoylphosphatidylserine (DOPS) and zwitterionic dioleoylphosphatidylcholine (DOPC) are investigated with label-free angle-resolved (AR) second harmonic scattering (SHS) and electrophoretic mobility measurements. The membrane surface potential is extracted from the AR-SHS response. The surface potential changes from -10 to -145 mV with varying DOPS content ( from 0% to 100%) and levels off already at ∼ 10 % DOPS content. The ζ-potential shows the same trend but with a drastically lower saturation value (-44 mV). This difference is explained by the formation of a condensed layer of Na⁺ counterions around the outer leaflet of the liposome as predicted by charge condensation theories for polyelectrolyte systems.

Unique orientations and rotational dynamics of a 1-butyl-3-methyl-imidazolium hexafluorophosphate ionic liquid at the gas–liquid interface: the effects of the hydrogen bond and hydrophobic interactions

Two-stage rotational motions of the interfacial [BMIM]⁺ cations are essentially determined by both hydrophobic and hydrogen-bonding interactions.

Effects of Dimethyl Sulfoxide on Surface Water near Phospholipid Bilayers

Despite much effort to probe the properties of dimethyl sulfoxide (DMSO) solution, the effects of DMSO on water, especially near plasma membrane surfaces, still remain elusive. By performing molecular dynamics simulations at varying DMSO concentrations (X_DMSO), we study how DMSO affects structural and dynamical properties of water in the vicinity of phospholipid bilayers. As proposed by a number of experiments, our simulations confirm that DMSO induces dehydration from bilayer surfaces and disrupts the H-bond structure of water. However, DMSO-enhanced water diffusivity at solvent-bilayer interfaces, an intriguing discovery reported by a spin-label measurement, is not confirmed in our simulations. To resolve this discrepancy, we examine the location of the spin label (Tempo) relative to the solvent-bilayer interface. In accord with the evidence in the literature, our simulations, which explicitly model Tempo-phosphatidylcholine, find that the Tempo moiety is equilibrated at ∼8-10 Å below the bilayer surface. Furthermore, the DMSO-enhanced surface-water diffusion is confirmed only when water diffusion is analyzed around the Tempo moiety that is immersed below the bilayer surface, which implies that the experimentally detected signal of water using Tempo stems from the interior of bilayers, not from the interface. Our analysis finds that the increase of water diffusion below the bilayer surface is coupled to the increase of area per lipid with an increasing X_DMSO(≲10mol%). Underscoring the hydrophobic nature of the Tempo moiety, our study calls for careful re-evaluation of the use of Tempo in measurements on lipid bilayer surfaces.

Single Particle Dynamics at the Intrinsic Surface of Various Apolar, Aprotic Dipolar, and Hydrogen Bonding Liquids As Seen from Computer Simulations

We investigate the single molecule dynamics at the intrinsic liquid/vapor interface of five different molecular liquids (carbon tetrachloride, acetone, acetonitrile, methanol, and water). After assessing that the characteristic residence times in the surface layer are long enough for a meaningful definition of several transport properties within the layer itself, we characterize the dynamics of the individual molecules at the liquid surface by analyzing their normal and lateral mean-square displacements and lateral velocity autocorrelation functions and, in the case of the hydrogen bonding liquids (i.e., water and methanol), also the properties of the hydrogen bonds. Further, dynamical properties as well as the clustering of the molecules residing unusually long in the surface layer are also investigated. The global picture emerging from this analysis is that of a noticeably enhanced dynamics of the molecules at the liquid surface, with diffusion coefficients up to 4 times larger than in the bulk, and the disappearance of the caging effect at the surface of all liquids but water. The dynamics of water is dominated by the strong hydrogen bonding structure also at the liquid surface.

Water populations in restricted environments of lipid membrane interphases

We employ molecular dynamics simulations to study the hydration properties of Dipalmitoylphosphatidylcholine (DPPC) bilayers, both in the gel and the liquid crystalline states. We show that while the tight hydration centers (PO and CO moieties) are significantly hydrated in both phases, the gel-fluid transition involves significant changes at the second hydration shell, particularly at the buried region between the hydrocarbon tails. Thus, while almost no buried water population exists in the gel state below the carbonyls, this hydrophobic region becomes partially water accesible in the liquid crystalline state. We shall also show that such water molecules present a lower H-bond coordination as compared to the molecules at the primary hydration shell. This means that, while the latter are arranged in relatively compact nanoclusters (as already proposed), the buried water molecules tend to organize themselves in less compact structures, typically strings or branched strings, with a scarce population of isolated molecules. This behavior is similar to that observed in other hydration contexts, like water penetrating carbon nanotubes or model hydrophobic channels or pores, and reflects the reluctance of water to sacrifice HB coordination.

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.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Origin of 1/f noise in hydration dynamics on lipid membrane surfaces

Water molecules on lipid membrane surfaces are known to contribute to membrane stability by connecting lipid molecules and acting as a water bridge. Although water structures and diffusivities near the membrane surfaces have been extensively studied, hydration dynamics on the surfaces has remained an open question. Here we investigate residence time statistics of water molecules on the surface of lipid membranes using all-atom molecular dynamics simulations. We show that hydration dynamics on the lipid membranes exhibits 1/f noise. Constructing a dichotomous process for the hydration dynamics, we find that residence times in each state follow a power-law with exponential cutoff and that the process can be regarded as a correlated renewal process where interoccurrence times are correlated. The results imply that the origin of the 1/f noise in hydration dynamics on the membrane surfaces is a combination of a power-law distribution with cutoff of interoccurrence times of switching events and a long-term correlation between the interoccurrence times. These results suggest that the 1/f noise attributed to the correlated renewal process may contribute to the stability of the hydration layers and lipid membranes.

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.

Membrane Hydration: A Hint to a New Model for Biomembranes

The classical view of a biological membrane is based on the Singer-Nicholson mosaic fluid model in which the lipid bilayer is the structural backbone. Under this paradigm, many studies of biological processes such as, permeability, active transport, enzyme activity and adhesion and fusion processes have been rationalized considering the lipid membrane as a low dielectric slab of hydrocarbon chains with polar head groups exposed to water at each side in which oil/water partition prevails. In spite of several analyses and evidence available in relation to membrane hydration, water is not taken into account as a functional component. For this purpose, new insights in the water organization in restricted environments and the thermodynamical and mechanical properties emerging from them are specifically analysed and correlated.This chapter summarizes the progress of the studies of water in membranes along the book in order to give a more realistic structural and dynamical picture accounting for the membrane functional properties.

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.

Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes

Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10(-5) cm(2)/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10(-8) cm(2)/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2-8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a "wagging of the tails" frequency around 30 cm(-1), which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.

Structural and thermodynamic properties of water-membrane interphases: significance for peptide/membrane interactions

Water appears as a common intermediary in the mechanisms of interaction of proteins and polypeptides with membranes of different lipid composition. In this review, how water modulates the interaction of peptides and proteins with lipid membranes is discussed by correlating the thermodynamic response and the structural changes of water at the membrane interphases. The thermodynamic properties of the lipid-protein interaction are governed by changes in the water activity of monolayers of different lipid composition according to the lateral surface pressure. In this context, different water populations can be characterized below and above the phase transition temperature in relation to the CH₂ conformers' states in the acyl chains. According to water species present at the interphase, lipid membrane acts as a water state regulator, which determines the interfacial water domains in the surface. It is proposed that those domains are formed by the contact between lipids themselves and between lipids and the water phase, which are needed to trigger adsorption-insertion processes. The water domains are essential to maintain functional dynamical properties and are formed by water beyond the hydration shell of the lipid head groups. These confined water domains probably carries information in local units in relation to the lipid composition thus accounting for the link between lipidomics and aquaomics. The analysis of these results contributes to a new insight of the lipid bilayer as a non-autonomous, responsive (reactive) structure that correlates with the dynamical properties of a living system.

Interaction of neurotransmitters with a phospholipid bilayer: a molecular dynamics study

We have performed a series of molecular dynamics simulations to study the interactions between the neurotransmitters (NTs) γ-aminobutyrate (GABA), glycine (GLY), acetylcholine (ACH) and glutamate (GLU) as well as the amidated/acetylated γ-aminobutyrate (GABA(neu)) and the osmolyte molecule glycerol (GOL) with a dipalmitoylphosphatidylcholine (DPPC) bilayer. In agreement with previously published experimental data, we found the lowest membrane affinity for the charged molecules and a moderate affinity for zwitterionic and polar molecules. The affinity can be ranked as follows: ACH-GLU<<GABA<GLY<<GABA(neu)<<GOL. The latter three penetrated the bilayer at most with the deepest location being close to the glycerol backbone of the phospholipids. Even at that position, these solutes were noticeably hydrated and carried ∼30-80% of the bulk water along. The mobility of hydration water at the solute is also affected by the penetration into the bilayer. Two time scales of exchanging water molecules could be determined. In the bulk phase, the hydration layer contains ∼20% slow exchanging water molecules which increases 2-3 times as the solutes entered the bilayer. Our results indicate that there is no intermediate exchange of slow moving water molecules from the solutes to the lipid atoms and vice versa. Instead, the exchange relies on the reservoir of unbounded ("free") water molecules in the interfacial bilayer region. Results from the equilibrium simulations are in good agreement with the results from umbrella sampling simulations, which were conducted for the four naturally occurring NTs. Free energy profiles for ACH and GLU show a minimum of ∼2-3 kJ/mol close to the bilayer interface, while for GABA and GLY, a minimum of respectively ∼2 kJ/mol and ∼5 kJ/mol is observed when these NTs are located in the vicinity of the lipid glycerol backbone. The most important interaction of NTs with the bilayer is the charged amino group of NTs with the lipid phosphate group.

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.

Hydration repulsion between membranes and polar surfaces: simulation approaches versus continuum theories

A review of various computer simulation approaches for the study of the hydration repulsion between lipid membranes and polar surfaces is presented. We discuss different methods and compare their advantages and limitations. We consider interaction pressures, interaction thermodynamics, and interaction mechanisms. We take a close look at the influence of the experimental boundary conditions and at repulsion mechanisms due to the unfavorable overlap of interfacial water layers. To this end, we analyze several distinct water order parameters in simulations of interacting polar surfaces and compare the results to the predictions of simple continuum theories.

Origin of subdiffusion of water molecules on cell membrane surfaces

Water molecules play an important role in providing unique environments for biological reactions on cell membranes. It is widely believed that water molecules form bridges that connect lipid molecules and stabilize cell membranes. Using all-atom molecular dynamics simulations, we show that translational and rotational diffusion of water molecules on lipid membrane surfaces exhibit subdiffusion and aging. Moreover, we provide evidence that both divergent mean trapping time (continuous-time random walk) and long-correlated noise (fractional Brownian motion) contribute to this subdiffusion. These results suggest that subdiffusion on cell membranes causes the water retardation, an enhancement of cell membrane stability, and a higher reaction efficiency.

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.

Nontrivial behavior of water in the vicinity and inside lipid bilayers as probed by molecular dynamics simulations

The atomic-scale diffusion of water in the presence of several lipid bilayers mimicking biomembranes is characterized via unconstrained molecular dynamics (MD) simulations. Although the overall water dynamics corresponds well to literature data, namely, the efficient braking near polar head groups of lipids, a number of interesting and biologically relevant details observed in this work have not been sufficiently discussed so far; for instance, the fact that waters "sense" the membrane unexpectedly early, before water density begins to decrease. In this "transitional zone" the velocity distributions of water and their H-bonding patterns deviate from those in the bulk solution. The boundaries of this zone are well preserved even despite the local (<1 nm size) perturbation of the lipid bilayer, thus indicating a decoupling of the surface and bulk dynamics of water. This is in excellent agreement with recent experimental data. Near the membrane surface, water movement becomes anisotropic, that is, solvent molecules preferentially move outward the bilayer. Deep in the membrane interior, the velocities can even exceed those in the bulk solvent and undergo large-scale fluctuations. The analysis of MD trajectories of individual waters in the middle part of the acyl chain region of lipids reveals a number of interesting rare phenomena, such as the fast (ca. 50 ps) breakthrough across the membrane or long-time (up to 750 ps) "roaming" between lipid leaflets. The analysis of these events was accomplished to delineate the mechanisms of spontaneous water permeation inside the hydrophobic membrane core. It was shown that such nontrivial dynamics of water in an "alien" environment is driven by the dynamic heterogeneities of the local bilayer structure and the formation of transient atomic-scale "defects" in it. The picture observed in lipid bilayers is drastically different from that in a primitive membrane mimic, a hydrated cyclohexane slab. The possible biological impact of such phenomena in equilibrated lipid bilayers is discussed.