Molecular dynamics simulation of short-wavelength collective dynamics of phospholipid membranes

Insights on Hydrogen Bond Network of Water in Phospholipid Membranes: An Infrared Study at Varying Hydration

Water in membrane interphases is vital for cellular biological functions, but despite its importance, the structure and function of biological water remain elusive. Here, by studying the OH stretching mode in partially hydrated lipid multilayers by FTIR measurements, relevant information on the water structure near the surface with lipid membranes has been gathered. The water hydrogen bond network is highly perturbed in the first layers that are in contact with the lipid membrane, exhibiting strong deviations from tetrahedral symmetry and a significant number of defects, such as isolated water molecules and a large number of hydrogen-bonded water dimers in the interphase region. These findings support the hypothesis that water chains form in phospholipid membranes, and are involved in the proton transfer across lipid bilayers by phosphate groups of opposing lipids. Furthermore, we have determined that even at very low hydration levels, a small amount of water is embedded within the confined spaces of the hydrocarbon region of phospholipid bilayers, which could potentially contribute to the structural stability of the lipid membrane.

Intra-protein interacting collective modes in the terahertz frequency region

Understanding how proteins work requires a thorough understanding of their internal dynamics. Proteins support a wide range of motions, from the femtoseconds to seconds time scale, relevant to crucial biological functions. In this context, the term "protein collective dynamics" refers to the complex patterns of coordinated motions of numerous atoms throughout the protein in the sub-picosecond time scale (terahertz frequency region). It is hypothesized that these dynamics have a substantial impact on the regulation of functional dynamical mechanisms, including ligand binding and allosteric signalling, charge transport direction, and the regulation of thermodynamic and thermal transport properties. Using the theoretical framework of hydrodynamics, the collective dynamics of proteins had previously been described in a manner akin to that of simple liquids, i.e. in terms of a single acoustic-like excitation, related to intra-protein vibrational motions. Here, we employ an interacting-mode model to analyse the results from molecular dynamics simulations and we unveil that the vibrational landscape of proteins is populated by multiple acoustic-like and low-frequency optic-like modes, with mixed symmetry and interfering with each other. We propose an interpretation at the molecular level of the observed scenario that we relate to the side-chains and the hydrogen-bonded networks dynamics. The present insights provide a perspective for understanding the molecular mechanisms underlying the energy redistribution processes in the interior of proteins.

The Phonon Theory of Liquids and Biological Fluids: Developments and Applications

Among the three basic states of matter (solid, liquid, and gas), the liquid state has always eluded general theoretical approaches for describing liquid energy and heat capacity. In this Viewpoint, we derive the phonon theory of liquids and biological fluids stemming from Frenkel's microscopic picture of the liquid state. Specifically, the theory predicts the existence of phonon gaps in vibrational spectra of liquids and a thermodynamic boundary in the supercritical state. Direct experimental evidence reaffirming these theoretical predictions was achieved through a combination of techniques using static compression X-ray diffraction and inelastic X-ray scattering on deeply supercritical argon in a diamond anvil cell. Furthermore, these findings inspired and then led to the discovery of phonon gaps in liquid crystals (mesogens), block copolymers, and biological membranes. Importantly, phonon gaps define viscoelastic crossovers in cellular membranes responsible for lipid self-diffusion, lateral molecular-level stress propagation, and passive transmembrane transport of small molecules and solutes. Finally, molecular interactions mediated by external stimuli result in synaptic activity controlling biological membranes' plasticity resulting in learning and memory. Therefore, we also discuss learning and memory effects─equally important for neuroscience as well as for the development of neuromorphic devices─facilitated in biological membranes by external stimuli.

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.

Molecular Picture of the Transient Nature of Lipid Rafts

In biological membranes, lipid rafts are now thought to be transient and nanoscopic. However, the mechanism responsible for these nanoscopic assemblies remains poorly understood, even in the case of model membranes. As a result, it has proven extremely challenging to probe the physicochemical properties of lipid rafts at the molecular level. Here, we use all-atom molecular dynamics (MD) simulations and inelastic X-ray scattering (IXS), an intrinsically nanoscale technique, to directly probe the energy transfer and collective short-wavelength dynamics (phonons) of biologically relevant model membranes. We show that the nanoscale propagation of stress in lipid rafts takes place in the form of collective motions made up of longitudinal (compression waves) and transverse (shear waves) molecular vibrations. Importantly, we provide a molecular picture for the so-called van der Waals mediated "force from lipid" [Anishkin, A. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 7898], a key parameter for the ionic channel mechano-transduction and the mechanism for the lipid transfer of molecular level stress [Aponte-Santamarı́a, C. et al. J. Am. Chem. Soc. 2017, 139, 13588]. Specifically, we describe how lipid rafts are formed and maintained through the propagation of molecular stress, lipid raft rattling dynamics, and a relaxation process. Eventually, the rafts dissipate through the self-diffusion of lipids making up the rafts. We also show that the molecular stress and viscoelastic properties of transient lipid rafts can be modulated through the use of hydrophobic biomolecules such as melatonin and tryptophan. Ultimately, the herein proposed mechanism describing the molecular interactions for the formation and dissolution of lipid rafts may offer insights as to how lipid rafts enable biological function.

Functional lipid pairs as building blocks of phase-separated membranes

Biological membranes exhibit a great deal of compositional and phase heterogeneity due to hundreds of chemically distinct components. As a result, phase separation processes in cell membranes are extremely difficult to study, especially at the molecular level. It is currently believed that the lateral membrane heterogeneity and the formation of domains, or rafts, are driven by lipid-lipid and lipid-protein interactions. Nevertheless, the underlying mechanisms regulating membrane heterogeneity remain poorly understood. In the present work, we combine inelastic X-ray scattering with molecular dynamics simulations to provide direct evidence for the existence of strongly coupled transient lipid pairs. These lipid pairs manifest themselves experimentally through optical vibrational (a.k.a. phononic) modes observed in binary (1,2-dipalmitoyl-sn-glycero-3-phosphocholine [DPPC]-cholesterol) and ternary (DPPC-1,2-dioleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-glycero-3-phosphocholine [DOPC/POPC]-cholesterol) systems. The existence of a phononic gap in these vibrational modes is a direct result of the finite size of patches formed by these lipid pairs. The observation of lipid pairs provides a spatial (subnanometer) and temporal (subnanosecond) window into the lipid-lipid interactions in complex mixtures of saturated/unsaturated lipids and cholesterol. Our findings represent a step toward understanding the lateral organization and dynamics of membrane domains using a well-validated probe with a high spatial and temporal resolution.

Dynamic structure factor of a lipid bilayer in the presence of a high electric field

The influence of a high average electric field (∼1 V/nm) in the hydrophobic interior of a bilayer lipid membrane on short-wavelength in-plane phononic motions of lipid chains is considered. The average electric field is assumed to be nearly constant on a picosecond time scale and a nanometer length scale. This field may be induced, for instance, by externally applied subnanosecond electric pulses or the membrane dipole potential. Using a generalized hydrodynamic approach, we derive a corresponding electrohydrodynamic model generalized to high wave numbers. In the considered approximation, all electric field effects are reduced only to a constant contribution to the generalized isothermal compressibility modulus. The corresponding dynamic structure factor for a lipid bilayer is derived. We show that due to polarization effects, the high field can critically impact the dynamics of longitudinal acousticlike modes at wave numbers near the major peak of the static structure factor. We estimate quantitatively that for typical lipid bilayers, transverse high electric fields can cause strong phonon energy softening, enhancement of phonon population, and formation of a gap in the dispersion of excitation frequency. The results obtained agree with simulations of the initiation of lipid bilayer electropores, suggesting that the proposed model reproduces the essential features of the field's impact on atomic density fluctuations. The proposed mechanism may have significant implications for the understanding of electroporation, passive molecular transport, and spontaneous pore formation in lipid bilayers.

Nonequilibrium dynamic structure factor of a lipid bilayer in the presence of an in-plane temperature gradient

There is rapidly increasing evidence that nanoscale temperature heterogeneities are involved in important biological processes. Combining nanoheating and nanoscale thermosensors forms the basis of emerging unique methods of cell therapy, tissue engineering, and regenerative medicine. Understanding corresponding phenomena seems to require a mesoscopic nonequilibrium hydrodynamic theory. In this paper, a Langevin-type model of dynamics of phonon modes propagating along a bilayer lipid membrane in the presence of an in-plane temperature gradient is proposed. Corresponding quantitative estimates for the Brillouin components of the nonequilibrium dynamic structure factor and the equal-time longitudinal momentum-density correlation function for a lipid bilayer are obtained. The analysis reveals that for typical values of parameters of lipid bilayer, the longitudinal temperature gradient of the order of 5qK for wave numbers q from 0.01 to 1nm^{-1} induces significant asymmetry of the Brillouin components of the dynamic structure factor and long-range spatial correlations in the plane of the bilayer. The corresponding membrane temperature gradients seem to be typical or achievable for cellular processes responsible for intracellular temperature variations and such external physical impacts as high-intensity electromagnetic pulses or heating of membrane-associated nanoparticles.

Multiple Interacting Collective Modes and Phonon Gap in Phospholipid Membranes

We combine Brillouin neutron scattering measurements with recent inelastic X-ray scattering [ Zhernenkov et al. Nat. Commun. 2016 , 7 , 11575 ] to propose a model for the collective dynamics of phospholipid bilayers. Neutron and X-ray spectra were fitted by the model response function associated with the Hamiltonian of an interacting-phonon system. This approach allows for a comprehensive and unprecedented picture of the vibrational collective features of phospholipids. At low wavevectors Q, the dispersion relations can be interpreted in terms of two acoustic-like modes, one longitudinal and one transverse, plus a dispersionless optic-like mode. The transverse mode of the liquid phase shows a phonon gap that can be linked to a passive transport mechanism through membranes, an interpretation that was proposed in Zhernenkov et al. At higher Q values, the interaction of the longitudinal acoustic excitation with the dispersionless mode gives rise to a pattern that is consistent with avoided-crossing behavior. Evidence is found for a slow- to fast-sound transition, similar to bulk water and other biomolecules.

Crossover from picosecond collective to single particle dynamics defines the mechanism of lateral lipid diffusion

It has been widely accepted that the thermally excited motions of the molecules in a cell membrane is the prerequisite for a cell to carry its biological functions. On the other hand, the detailed mapping of the ultrafast picosecond single-molecule and the collective lipid dynamics in a cell membrane remains rather elusive. Here, we report all-atom molecular dynamics simulations of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer over a wide range of temperature. We elucidate a molecular mechanism underlying the lateral lipid diffusion in a cell membrane across the gel, rippled, and liquid phases using an analysis of the longitudinal and transverse current correlation spectra, the velocity auto-correlation functions, and the molecules mean square displacements. The molecular mechanism is based on the anomalous ultrafast vibrational properties of lipid molecules at the viscous-to-elastic crossover. The macroscopic lipid diffusion coefficients predicted by the proposed diffusion model are in a good agreement with experimentally observed values. Furthermore, we unveil the role of water confined at the water-lipid interface in triggering collective vibrations in a lipid bilayer.

Anomalous Nanoscale Optoacoustic Phonon Mixing in Nematic Mesogens

Recent inelastic X-ray scattering (IXS) experiments on mesogens have revealed entirely new capabilities with regards to their nanoscale phonon-assisted heat management. Mesogens such as nematic liquid crystals (LCs) are appealing systems for study because their structure and morphology can easily be tuned. We report on Q-resolved ultra-high-resolution IXS, X-ray diffraction, and THz time-domain spectroscopy experiments combined with large-scale all-atom molecular dynamics simulations on the dynamic properties of 5CB LCs. For the first time, we observe a strong mixing of phonon excitations originating from independent in-phase and out-of-phase van-der-Waals-mediated displacement patterns. The coexistence of transverse acoustic and optical modes of 5CB LCs at near room temperature is revealed through the emergent transverse phonon gap and THz light-phonon coupling taking place within the same energy range. Furthermore, our experimental observations are supported by analysis showing correlations of spontaneous fluctuations of LCs on picosecond time scales. These findings are significant for the design of a new generation of soft molecular vibration-sensitive nanoacoustic and optomechanical applications.

Aspirin locally disrupts the liquid-ordered phase

Local structure and dynamics of lipid membranes play an important role in membrane function. The diffusion of small molecules, the curvature of lipids around a protein and the existence of cholesterol-rich lipid domains (rafts) are examples for the membrane to serve as a functional interface. The collective fluctuations of lipid tails, in particular, are relevant for diffusion of membrane constituents and small molecules in and across membranes, and for structure and formation of membrane domains. We studied the effect of aspirin (acetylsalicylic acid, ASA) on local structure and dynamics of membranes composed of dimyristoylphosphocholine (DMPC) and cholesterol. Aspirin is a common analgesic, but is also used in the treatment of cholesterol. Using coherent inelastic neutron scattering experiments and molecular dynamics (MD) simulations, we present evidence that ASA binds to liquid-ordered, raft-like domains and disturbs domain organization and dampens collective fluctuations. By hydrogen-bonding to lipid molecules, ASA forms 'superfluid' complexes with lipid molecules that can organize laterally in superlattices and suppress cholesterol's ordering effect.

Role of Terahertz (THz) Fluctuations in the Allosteric Properties of the PDZ Domains

With the aim of investigating the relationship between the fast fluctuations of proteins and their allosteric behavior, we perform molecular dynamics simulations of two model PDZ domains with differential allosteric responses. We focus on protein dynamics in the THz regime (0.1–3 THz) as opposed to lower frequencies. By characterizing the dynamic modulation of the protein backbone induced by ligand binding in terms of single residue and pairwise distance fluctuations, we identify a response nucleus modulated by the ligand that is visible only at THz frequencies. The residues of this nucleus undergo a significant stiffening and an increase in mutual coordination upon binding. Additionally, we find that the dynamic modulation is significantly more intense for the side chains, where it is also redistributed to distal regions not immediately in contact with the ligand allowing us to better define the response nucleus at THz frequencies. The overlap between the known allosterically responding residues of the investigated PDZ domains and the modulated region highlighted here suggests that fast THz dynamics could play a role in allosteric mechanisms.

Probing Intermolecular Interactions in Phospholipid Bilayers by Far-Infrared Spectroscopy

Fast thermal fluctuations and low frequency phonon modes are thought to play a part in the dynamic mechanisms of many important biological functions in cell membranes. Here we report a detailed far-infrared study of the molecular subpicosecond motions of phospholipid bilayers at various hydrations. We show that these systems sustain several low frequency collective modes and deduce that they arise from vibrations of different lipids interacting through intermolecular van der Waals forces. Furthermore, we observe that the low frequency vibrations of lipid membrane have strong similarities with the subpicosecond motions of liquid water and suggest that resonance mechanisms are an important element to the dynamics coupling between membranes and their hydration water.

Molecular Dynamics of POPC Phospholipid Bilayers through the Gel to Fluid Phase Transition: An Incoherent Quasi-Elastic Neutron Scattering Study

The microscopic dynamics for the gel and liquid-crystalline phase of highly aligned D2O-hydrated bilayers of 1-palmitoyl-oleoyl-sn-glycero-phosphocholine (POPC) were investigated in the temperature range from 248 to 273 K by using incoherent quasi-elastic neutrons scattering (QENS). We develop a model for describing the molecular motions of the liquid phase occurring in the 0.3 to 350 ps time range. Accordingly, the complex dynamics of hydrogen are described in terms of simple dynamical processes involving different parts of the phospholipid chain. The analysis of the data evidences the existence of three different motions: the fast motion of hydrogen vibrating around the carbon atoms, the intermediate motion of carbon atoms in the acyl chains, and the slower translational motion of the entire phospholipid molecule. The influence of the temperature on these dynamical processes is investigated. In particular, by going from gel to liquid-crystalline phase, we reveal an increase of the segmental motion mainly affecting the terminal part of the acyl chains and a change of the diffusional dynamics from a localized rattling-like motion to a confined diffusion.

Revealing the mechanism of passive transport in lipid bilayers via phonon-mediated nanometre-scale density fluctuations

The passive transport of molecules through a cell membrane relies on thermal motions of the lipids. However, the nature of transmembrane transport and the precise mechanism remain elusive and call for a comprehensive study of phonon excitations. Here we report a high resolution inelastic X-ray scattering study of the in-plane phonon excitations in 1,2-dipalmitoyl-sn-glycero-3-phosphocholine above and below the main transition temperature. In the gel phase, for the first time, we observe low-frequency transverse modes, which exhibit a phonon gap when the lipid transitions into the fluid phase. We argue that the phonon gap signifies the formation of short-lived nanometre-scale lipid clusters and transient pores, which facilitate the passive molecular transport across the bilayer plane. Our findings suggest that the phononic motion of the hydrocarbon tails provides an effective mechanism of passive transport, and illustrate the importance of the collective dynamics of biomembranes.

The molecular transport through bio-membranes of cells heavily relies on the dynamics of lipids, but the related mechanism remains unknown. Here, Zhernenkov et al. observe the propagating transverse phonon mode with a finite band gap and suggest its connection to short-lived local lipid clustering.