Anomalous diffusion of water molecules in hydrated lipid bilayers

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.

Simulations reveal that antimicrobial BP100 induces local membrane thinning, slows lipid dynamics and favors water penetration

MD simulations reveal that BP100 peptide induces local membrane thinning and negative curvature, slows lipid dynamics and increases the water life time in the lipid hydrophobic core and transmembrane water transport in the direction of the peptide.

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.

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and X-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts (Bolmatov et al., 2020). Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager's regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

Dynamic coupling of a hydration layer to a fluid phospholipid membrane: intermittency and multiple time-scale relaxations

Understanding the coupling of a hydration layer and a lipid membrane is crucial to gaining access to membrane dynamics and understanding its functionality towards various biological processes.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Anomalous Dynamics of Water in Polyamide Matrix

Water in polymer matrixes is likely to show anomalous dynamics, a problem that has not been well understood yet. Here, we performed atomistic molecular dynamics simulations to study the water dynamics in a polyamide (PA) matrix, the bulk phase of well-known reverse osmosis membranes. For time-dependent ensemble average, water molecules experienced ballistic diffusion at a shorter time scale, followed by a crossover from subdiffusion to Brownian diffusion at a time scale ∼10 ns, and non-Gaussian diffusion, an indication of anomalous dynamics, sticks on even in the Brownian diffusion region. The anomalous dynamics mainly originates from two distinct motions including small-step continuous diffusion and jumping diffusion. The jumping motion has a mean length of 3.08 ± 0.31 Å and characteristic relaxation time of 0.218 ± 0.040 ns, which dominates the water diffusion in a fully hydrated PA matrix. It comprised low- and high-frequency jumps; the former is almost unchanged, and the latter remarkably increases with the increase of the hydration level. Surrounding neighbors of water strongly affect the jumping frequency, which exponentially or linearly decays with the increase in the number of atoms from the PA matrix. Although the PA matrix is flexible, associated with the water dynamics, the translocation of water is mainly through either tracing the position of neighboring water or jumping into the adjacent accommodation space.

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.

Hydration-mediated stiffening of collective membrane dynamics by cholesterol

Hydration water governs the cholesterol-induced changes in collective headgroup dynamics in lipid bilayers.

Influence of Cholesterol on the Dynamics of Hydration in Phospholipid Bilayers

We investigate the dynamics of interfacial waters in dipalmitoylphosphatidylcholine (DPPC) bilayers upon the addition of cholesterol, by molecular dynamics simulations. Our data reveal that the inclusion of cholesterol modifies the membrane aqueous interfacial dynamics: waters diffuse faster, their rotational decay time is shorter, and the DPPC/water hydrogen bond dynamics relaxes faster than in the pure DPPC membrane. The observed acceleration of the translational water dynamics agrees with recent experimental results, in which, by means of NMR techniques, an increment of the surface water diffusivity is measured upon the addition of cholesterol. A microscopic analysis of the lipid/water hydrogen bond network at the interfacial region suggests that the mechanism underlying the observed water mobility enhancement is given by the rupture of a fraction of interlipid water bridge hydrogen bonds connecting two different DPPC molecules, concomitant to the formation of new lipid/solvent bonds, whose dynamics is faster than that of the former. The consideration of a simple two-state model for the decay of the hydrogen bond correlation function yielded excellent results, obtaining two well-separated characteristic time scales: a slow one (∼250 ps) associated with bonds linking two DPPC molecules, and a fast one (∼15 ps), related to DPPC/solvent bonds.

Hydration water dynamics around a protein surface: a first passage time approach

A stochastic noise-driven dynamic model is proposed to study the diffusion of water molecules around a protein surface, under the effect of thermal fluctuations that arise due to the collision of water molecules with the surrounding environment. The underlying dynamics of such a system may be described in the framework of the generalized Langevin equation, where the thermal fluctuations are assumed to be algebraically correlated in time, which governs the non-Markovian behavior of the system. Results of the calculations of mean-square displacement and the velocity autocorrelation function reveal that the hydration water around the protein surface follows subdiffusive dynamics at long times. Analytical expressions for the first passage time distribution, survival probability, mean residence time and mean first passage time of water molecules are derived for different boundary conditions, to analyze hydration water dynamics under the effect of thermally correlated noise. The results depict a unimodal distribution of the first passage time unlike Brownian motion. The survival probability of hydration water follows a stretched exponential decay for both boundary conditions. The mean residence time of the hydration water molecule for different initial positions increases with increase in the complexity/heterogeneity of the surrounding environment for both boundary conditions. The mean first passage time of the water molecule to reach the absorbing/reflecting boundary follows an asymptotic power law with respect to the thickness of the hydration layer, and increases with increase in the complexity/heterogeneity of the environment.

Hydration Dynamics of a Peripheral Membrane Protein

Water dynamics in the hydration shell of the peripheral membrane protein annexin B12 were studied using MD simulations and Overhauser DNP-enhanced NMR. We show that retardation of water motions near phospholipid bilayers is extended by the presence of a membrane-bound protein, up to around 10 Å above that protein. Near the membrane surface, electrostatic interactions with the lipid head groups strongly slow down water dynamics, whereas protein-induced water retardation is weaker and dominates only at distances beyond 10 Å from the membrane surface. The results can be understood from a simple model based on additive contributions from the membrane and the protein to the activation free energy barriers of water diffusion next to the biomolecular surfaces. Furthermore, analysis of the intermolecular vibrations of the water network reveals that retarded water motions near the membrane shift the vibrational modes to higher frequencies, which we used to identify an entropy gradient from the membrane surface toward the bulk water. Our results have implications for processes that take place at lipid membrane surfaces, including molecular recognition, binding, and protein-protein interactions.

Anomalous Dynamics of a Lipid Recognition Protein on a Membrane Surface

Pleckstrin homology (PH) domains are lipid-binding modules present in peripheral membrane proteins which interact with phosphatidyl-inositol phosphates (PIPs) in cell membranes. We use multiscale molecular dynamics simulations to characterize the localization and anomalous dynamics of the DAPP1 PH domain on the surface of a PIP-containing lipid bilayer. Both translational and rotational diffusion of the PH domain on the lipid membrane surface exhibit transient subdiffusion, with an exponent α ≈ 0.5 for times of less than 10 ns. In addition to a PIP₃ molecule at the canonical binding site of the PH domain, we observe additional PIP molecules in contact with the protein. Fluctuations in the number of PIPs associated with the PH domain exhibit 1/f noise. We suggest that the anomalous diffusion and long-term correlated interaction of the PH domain with the membrane may contribute to an enhanced probability of encounter with target complexes on cell membrane surfaces.

Theoretical vibrational sum-frequency generation spectroscopy of water near lipid and surfactant monolayer interfaces. II. Two-dimensional spectra

The structural stability and function of biomolecules is strongly influenced by the dynamics and hydrogen bonding of interfacial water. Understanding and characterizing the dynamics of these water molecules require a surface-sensitive technique such as two-dimensional vibrational sum-frequency generation (2DSFG) spectroscopy. We have combined theoretical 2DSFG calculations with molecular dynamics simulations in order to investigate the dynamics of water near different lipid and surfactant monolayer surfaces. We show that 2DSFG can distinguish the dynamics of interfacial water as a function of the lipid charge and headgroup chemistry. The dynamics of water is slow compared to the bulk near water-zwitterionic and water-anionic interfaces due to conformational constraints on interfacial water imposed by strong phosphate-water hydrogen bonding. The dynamics of water is somewhat faster near water-cationic lipid interfaces as no such constraint is present. Using hydrogen bonding and rotational correlation functions, we characterize the dynamics of water as a function of the distance from the interface between water and zwitterionic lipids. We find that there is a transition from bulk-like to interface-like dynamics approximately 7 Å away from a zwitterionic phosphatidylcholine monolayer surface.

Structure and dynamics of ionic micelles: MD simulation and neutron scattering study

Fully atomistic molecular dynamics (MD) simulations have been carried out on sodium dodecyl sulfate (SDS), an anionic micelle, and three cationic (CnTAB; n = 12, 14, 16) micelles, investigating the effects of size, the form of the headgroup, and chain length. They have been used to analyze neutron scattering data. MD simulations confirm the dynamical model of global motion of the whole micelle, segmental motion (headgroup and alkyl chain), and fast torsional motion associated with the surfactants that is used to analyze the experimental data. It is found that the solvent surrounding the headgroups results in their significant mobility, which exceeds that of the tails on the nanosecond time scale. The middle of the chain is found to be least mobile, consolidating the micellar configuration. This dynamical feature is similar for all the ionic micelles investigated and therefore independent of headgroup form and charge and chain length. Diffusion constants for global and segmental motion of the different micelles are consistent with experimentally obtained values as well as known structural features. This work provides a more realistic model of micelle dynamics and offers new insight into the strongly fluctuating surface of micelles which is important in understanding micelle dispersion and related functionality, like drug delivery.

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.

Atomically detailed lipid bilayer models for the interpretation of small angle neutron and X-ray scattering data

We present a new atom density profile (ADP) model and a statistical approach for extracting structural characteristics of lipid bilayers from X-ray and neutron scattering data. Models for five lipids with varying head and tail chemical composition in the fluid phase, 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), are optimized using a simplex based method to simultaneously reproduce both neutron and X-ray scattering data. Structural properties are determined using statistical analysis of multiple optimal model structures. The method and models presented make minimal assumptions regarding the atomic configuration, while taking into account the underlying physical properties of the system. The more general model and statistical approach yield data with well defined uncertainties, indicating the precision in determining density profiles, atomic locations, and bilayer structural characteristics. Resulting bilayer structures include regions exhibiting large conformational variation. Due to the increased detail in the model, the results demonstrate the possibility of a distinct hydration layer within the interfacial (backbone) region.

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.

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.