Hydration dynamics of a lipid membrane: Hydrogen bond networks and lipid-lipid associations

Optical and molecular features of negatively curved surfaces created by POPE lipids: A crucial role of the initial conditions

Membrane fusion is closely related to plasma membrane domains rich in cone-shaped phosphatidylethanolamine (PE) lipids that can reverse membrane curvature under certain conditions. The phase transition of PE-based lipid membranes from the lamellar fluid phase (Lα) to the inverse hexagonal phase (HII) is commonly taken as a general model in reconstructing the membrane fusion pathway, and whose structural features have been mostly described so far using structural and microscopic techniques. The aim of this paper is to decipher the optical and molecular features of Lβ → Lα and especially of Lα → HII transition of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) lipids at pH = 7.0 when they are initially prepared in the form of both multi- and unilamellar liposomes (MLVs and LUVs). The distinction between optical properties of MLS- and LUVs-derived HII phase, provided from turbidity-sensitive temperature-dependent UV-Vis spectra, was attributed to different formation mechanisms of HII phase. Most importantly, from FTIR spectroscopic data of POPE lipids in Lβ (15 °C), Lα (50 °C) and HII (85 °C) phases we identified the changes in molecular features of POPE lipids during phase transitions. Among the latter, by far the most significant is different hydration pattern of POPE lipids in MLVs- and LUVs-derived HII phase which extends from the polar-apolar interface all the way to the terminal amino group of the POPE lipid, along with the changes in the conformation of glycerol backbone as evidenced by the signature of α-methylene groups. Molecular dynamics simulations confirmed higher water penetration in HII phase and provided insight into hydrogen bonding patterns.

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.

Rotational Dynamics of Water near Osmolytes by Molecular Dynamics Simulations

The behavior of water molecules around organic molecules has attracted considerable attention as a crucial factor influencing the properties and functions of soft matter and biomolecules. Recently, it has been suggested that the change in protein stability upon the addition of small organic molecules (osmolytes) is dominated by the change in the water dynamics caused by the osmolyte, where the dynamics of not only the directly interacting water molecules but also the long-range hydration layer affect the protein stability. However, the relation between the long-range structure of hydration water in various solutions and the water dynamics remains unclear at the molecular level. We performed density-functional tight-binding molecular dynamics simulations to elucidate the varying rotational dynamics of water molecules in 15 osmolyte solutions. A positive correlation was observed between the rotational relaxation time and our proposed normalized parameter obtained by dividing the number of hydrogen bonds between water molecules by the number of nearest-neighbor water molecules. For the 15 osmolyte solutions, an increase or a decrease in the value of the normalized parameter for the second hydration shell tended to result in water molecules with slow and fast rotational dynamics, respectively, thus illustrating the importance of the second hydration shell for the rotational dynamics of water molecules. Our simulation results are anticipated to advance the current understanding of water dynamics around organic molecules and the long-range structure of water molecules.

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.

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Coarse-grained molecular dynamics simulation of cation distribution profiles on negatively charged lipid membranes during phase separation

Revealing the ion distributions on a charged lipid membrane in aqueous solution under the influence of long-range interactions is essential for understanding the origin of the stability of the bilayer structure and the interaction between biomembranes and various electrolytes. However, the ion distributions and their dynamics associated with the phase separation process of the lipid bilayer membrane are still unclear. We perform coarse-grained molecular dynamics simulations to reveal the Na⁺ and Cl^- distributions on charged phospholipid bilayer membranes during phase separation. During the phase separation, cations closely follow the position of negatively charged lipids on a microsecond timescale and are rapidly redistributed parallel to the lipid bilayer. In the homogenous mixture of zwitterionic and negatively charged lipids, cations weakly follow negatively charged lipids, indicating the strong interaction between cations and negatively charged lipids. We also compare cation concentrations as a function of surface charge density obtained by our simulation with those obtained by a modified Poisson-Boltzmann theory. Including the ion finite size makes the statistical results consistent, suggesting the importance of the ion-ion interactions in aqueous solution. Our simulation results advance our understanding of ion distribution during phase separation.

Real Space and Time Imaging of Collective Headgroup Dipole Motions in Zwitterionic Lipid Bilayers

Lipid bilayers are supramolecular structures responsible for a range of processes, such as transmembrane transport of ions and solutes, and sorting and replication of genetic materials, to name just a few. Some of these processes are transient and currently, cannot be visualized in real space and time. Here, we developed an approach using 1D, 2D, and 3D Van Hove correlation functions to image collective headgroup dipole motions in zwitterionic phospholipid bilayers. We show that both 2D and 3D spatiotemporal images of headgroup dipoles are consistent with commonly understood dynamic features of fluids. However, analysis of the 1D Van Hove function reveals lateral transient and re-emergent collective dynamics of the headgroup dipoles-occurring at picosecond time scales-that transmit and dissipate heat at longer times, due to relaxation processes. At the same time, the headgroup dipoles also generate membrane surface undulations due a collective tilting of the headgroup dipoles. A continuous intensity band of headgroup dipole spatiotemporal correlations-at nanometer length and nanosecond time scales-indicates that dipoles undergo stretching and squeezing elastic deformations. Importantly, the above mentioned intrinsic headgroup dipole motions can be externally stimulated at GHz-frequency scale, enhancing their flexoelectric and piezoelectric capabilities (i.e., increased conversion efficiency of mechanical energy into electric energy). In conclusion, we discuss how lipid membranes can provide molecular-level insights about biological learning and memory, and as platforms for the development of the next generation of neuromorphic computers.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

Structural and Dynamical Studies of a Lipid-Nanoclay Composite Layer at the Air-Water Interface

We modulate the adsorption affinities of nanoclay particles for the air-water interface by changing the cationic surface charge composition of the lipid monolayer and thereby tune the attractive electrostatic interaction between the positively charged lipid layer and the zwitterionic nanoclay particles in the water subphase. Our findings emphasize the significance of electrostatic interaction between lipids and the nanoclay, as well as its impact on the structural and viscoelastic features of the composite layer. We use surface pressure (Π)-mean molecular area (A) isotherms, atomic force microscope (AFM), Brewster angle microscopy (BAM), and energy dispersive X-ray spectrsocopy (EDXS) measurements to analyze the structure phases of lipid and lipid-nanoclay composite interfacial layer. The Π-A isotherm curve shows that the lipid-nanoclay composite layer has a larger lift-off area than the neat lipid layer, indicating that nanoparticles adsorb at the lipid layer via electrostatic interaction between lipid and nanoclay molecules. The surface density of the adsorbed nanoclay particles increases with an increase in the composition of the cationic lipid molecules. The stress relaxation response of the composite layer, measured using step compression measurements, exhibits exponential decay and ubiquitous dependence on the cationic dimyristoy-trimethylammonium propane (DMTAP) composition in the lipid layer with crossover to faster relaxation dynamics at DMTAP > 0.75. The power-law study of the frequency-dependent dynamic viscoelastic responses of the interfacial layer, measured using the barrier oscillation method, reveals a transition from glass-like response from neat lipid layer to gel-like dynamic response for the lipid-nanoclay composite layer. A solid-like behavior is evident for all the interface layers with dilation elastic modulus (E') > dilational viscous modulus (E″); however, the dynamic response of the neat layer is largely frequency-independent, whereas lipid-nanoclay composite layers with DMTAP > 0.75 reveal a frequency-dependent dynamic responses. The frequency-dependent power-law exponent of E', E″ increases on increasing the fractional composition of cationic DMTAP from 0.1 to 1.0, which forms a saturated interface of laponite particles and behaves as a viscoelastic gel in 2D.

Lipid/water interface of galactolipid bilayers in different lyotropic liquid-crystalline phases

In this study, carried out using computational methods, the organisation of the lipid/water interface of bilayers composed of galactolipids with both α-linolenoyl acyl chains is analysed and compared in three different lyotropic liquid-crystalline phases. These systems include the monogalactosyldiglyceride (MGDG) and digalactosyldiglyceride (DGDG) bilayers in the lamellar phase, the MGDG double bilayer during stalk phase formation and the inverse hexagonal MGDG phase. For each system, lipid-water and direct and water-mediated lipid-lipid interactions between the lipids of one bilayer leaflet and those of two apposing leaflets at the onset of new phase (stalk) formation, are identified. A network of interactions between DGDG molecules and its topological properties are derived and compared to those for the MGDG bilayer.

Effect of encapsulated protein on the dynamics of lipid sponge phase: a neutron spin echo and molecular dynamics simulation study

Lipid membranes are highly mobile systems with hierarchical, time and length scale dependent, collective motions including thickness fluctuations, undulations, and topological membrane changes, which play an important role in membrane interactions. In this work we have characterised the effect of encapsulating two industrially important enzymes, β-galactosidase and aspartic protease, in lipid sponge phase nanoparticles on the dynamics of the lipid membrane using neutron spin echo (NSE) spectroscopy and molecular dynamics (MD) simulations. From NSE, reduced membrane dynamics were observed upon enzyme encapsulation, which were dependent on the enzyme concentration and type. By fitting the intermediate scattering functions (ISFs) with a modified Zilman and Granek model including nanoparticle diffusion, an increase in membrane bending rigidity was observed, with a larger effect for β-galactosidase than aspartic protease at the same concentration. MD simulations for the system with and without aspartic protease showed that the lipids relax more slowly in the system with protein due to the replacement of the lipid carbonyl-water hydrogen bonds with lipid-protein hydrogen bonds. This indicates that the most likely cause of the increase in membrane rigidity observed in the NSE measurements was dehydration of the lipid head groups. The dynamics of the protein itself were also studied, which showed a stable secondary structure of protein over the simulation, indicating no unfolding events occurred.

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.

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".

New insights on the mechanism of polyethylenimine transfection and their implications on gene therapy and DNA vaccines

Polyethylenimine (PEI) has been demonstrated as an efficient DNA delivery vehicle both in vitro and in vivo. There is a consensus that PEI-DNA complexes enter the cells by endocytosis and escape from endosomes by the so-called "proton sponge" effect. However, little is known on how and where the polyplexes are de-complexed for DNA transcription and replication to occur inside the cell nucleus. To better understand this issue, we (i) tracked the cell internalization of PEI upon transfection to human epithelial cells and (ii) studied the interaction of PEI with phospholipidic layers mimicking nuclear membranes. Both the biological and physicochemical experiments provided evidence of a strong binding affinity between PEI and the lipidic bilayer. Firstly, confocal microscopy revealed that PEI alone could not penetrate the cell nucleus; instead, it arranged throughout the cytoplasm and formed a sort of aureole surrounding the nuclei periphery. Secondly, surface tension measurements, fluorescence dye leakage assays, and differential scanning calorimetry demonstrated that a combination of hydrophobic and electrostatic interactions between PEI and the phospholipidic monolayers/bilayers led to the formation of stable defects along the model membranes, allowing the intercalation of PEI through the monolayer/bilayer structure. Results are also supported by molecular dynamics simulation of the pore formation in PEI-lipidic bilayers. As discussed throughout the text, these results might shed light on a the mechanism in which the interaction between PEI and the nucleus membrane might play an active role on the DNA release: on the one hand, the PEI-membrane interaction is anticipated to facilitate the DNA disassembly from the polyplex by establishing a competition with DNA for the PEI binding and on the other hand, the forming defects are expected to serve as channels for the entrance of de-complexed DNA into the cell nucleus. A better understanding of the mechanism of transfection of cationic polymers opens paths to development of more efficiency vectors to improve gene therapy treatment and the new generation of DNA vaccines.

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.

Interfacial hydration determines orientational and functional dimorphism of sterol-derived Raman tags in lipid-coated nanoparticles

Lipid-coated noble metal nanoparticles (L-NPs) combine the biomimetic surface properties of a self-assembled lipid membrane with the plasmonic properties of a nanoparticle (NP) core. In this work, we investigate derivatives of cholesterol, which can be found in high concentrations in biological membranes, and other terpenoids, as tunable, synthetic platforms to functionalize L-NPs. Side chains of different length and polarity, with a terminal alkyne group as Raman label, are introduced into cholesterol and betulin frameworks. The synthesized tags are shown to coexist in two conformations in the lipid layer of the L-NPs, identified as "head-out" and "head-in" orientations, whose relative ratio is determined by their interactions with the lipid-water hydrogen-bonding network. The orientational dimorphism of the tags introduces orthogonal functionalities into the NP surface for selective targeting and plasmon-enhanced Raman sensing, which is utilized for the identification and Raman imaging of epidermal growth factor receptor-overexpressing cancer cells.

Rotational Dynamics of Water at the Phospholipid Bilayer Depending on the Head Groups Studied by Molecular Dynamics Simulations

Hydration states are a crucial factor that affect the self-assembly and properties of soft materials and biomolecules. Although previous experiments have revealed that the hydration state strongly depends on the chemical structure of lipid molecules, the mechanisms at the molecular level remain unknown. Classical and density-functional tight-binding (DFTB) molecular dynamics (MD) simulations are employed to determine the mechanisms underlying dissimilar water dynamics between lipid membranes with phosphatidylcholine (PC) and phosphatidylethanolamine (PE) head groups. The classical MD simulation shows that rotational relaxations of water are faster on the PE lipid than on the PC lipid, which is consistent with a previous experimental study using terahertz spectroscopy. Furthermore, DFTB-MD simulation of N(CH₃)₄⁺ and NH₄⁺ ions, which correspond to the different head groups in PC and PE, respectively, shows qualitative agreement with the classical MD simulation. Remarkably, the PE lipids and the NH₄⁺ ions break hydrogen bonds between water molecules in the secondary hydration shell. In contrast, the PC lipids and the N(CH₃)₄⁺ ions bind water molecules weakly in both the primary and secondary hydration shells and increase hydrogen bonds between water. Our atomistic simulations show that these changes in the hydrogen-bond network of water molecules cause the different rotational relaxation of water molecules between the two lipids.

Modulation of Membrane Fluidity to Control Interfacial Water Structure and Dynamics in Saturated and Unsaturated Phospholipid Vesicles

The structure and dynamics of interfacial water in biological systems regulate the biochemical reactions. But, it is still enigmatic how the behavior of the interfacial water molecule is controlled. Here, we have investigated the effect of membrane fluidity on the structure and dynamics of interfacial water molecules in biologically relevant phopholipid vesicles. This study delineates that modulation of membrane fluidity through interlipid separation and unsaturation not only mitigate membrane rigidity but also disrupt the strong hydrogen bond (H-bond) network around the lipid bilayer interface. As a result, a disorder in H-bonding between water molecules arises several layers beyond the first hydration shell of the polar headgroup, which essentially modifies the interfacial water structure and dynamics. Furthermore, we have also provided evidence of increasing transportation through these modulated membranes, which enhance the membrane mediated isomerization reaction rate.

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.

Unusual confinement properties of a water insoluble small peptide hydrogel

Unlike polymeric hydrogels, in the case of supramolecular hydrogels, the cross-linked network formation is governed by non-covalent forces. Hence, in these cases, the gelator molecules inside the network retain their characteristic physicochemical properties as no covalent modification is involved. Supramolecular hydrogels thus get dissolved easily in aqueous medium as the dissolution leads to a gain in entropy. Thus, any supramolecular hydrogel, insoluble in bulk water, is beyond the present understanding and hitherto not reported as well. Herein, we present a peptide-based (PyKC) hydrogel which remained insoluble in water for more than a year. Moreover, in the gel state, any movement of solvent or solute to and from the hydrogel is highly restricted resulting in a high degree of compartmentalization. The hydrogel could be re-dissolved in the presence of some biomolecules which makes it a prospective material for in vivo applications. Experimental studies and all atom molecular dynamics simulations revealed that a cysteine containing gelator forms dimers through disulfide linkage which self-assemble into PyKC layers with a distinct PyKC-water interface. The hydrogel is stabilized by intra-molecular hydrogen bonds within the peptide-conjugates and the π-π stacking of the pyrene rings. The unique confinement ability of the hydrogel is attributed to the slow dynamics of water which remains confined in the core region of PyKC via hydrogen bonds. The hydrogen bonds present in the confined water need activation energies to move through the water depleted hydrophobic environment of pyrene rings which significantly reduces water transport across the hydrogel. The compartmentalizing ability is effectively used to protect enzymes for a long time from denaturing agents like urea, heat or methanol. Overall, the presented system shows unique insolubility and confinement properties that could be a milestone in the research of soft-materials.

Hydration-mediated stiffening of collective membrane dynamics by cholesterol

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

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.