Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes

Neutron spin echo shows pHLIP is capable of retarding membrane thickness fluctuations

Cell membranes are responsible for a range of biological processes that require interactions between lipids and proteins. While the effects of lipids on proteins are becoming better understood, our knowledge of how protein conformational changes influence membrane dynamics remains rudimentary. Here, we performed experiments and computer simulations to study the dynamic response of a lipid membrane to changes in the conformational state of pH-low insertion peptide (pHLIP), which transitions from a surface-associated (SA) state at neutral or basic pH to a transmembrane (TM) α-helix under acidic conditions. Our results show that TM-pHLIP significantly slows down membrane thickness fluctuations due to an increase in effective membrane viscosity. Our findings suggest a possible membrane regulatory mechanism, where the TM helix affects lipid chain conformations, and subsequently alters membrane fluctuations and viscosity.

Significance of in situ quantitative membrane property-morphology relation (QmPMR) analysis

Deformation of the cell membrane is well understood from the viewpoint of protein interactions and free energy balance. However, the various dynamic properties of the membrane, such as lipid packing and hydrophobicity, and their relationship with cell membrane deformation are unknown. Therefore, the deformation of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and oleic acid (OA) giant unilamellar vesicles (GUVs) was induced by heating and cooling cycles, and time-lapse analysis was conducted based on the membrane hydrophobicity and physical parameters of "single-parent" and "daughter" vesicles. Fluorescence ratiometric analysis by simultaneous dual-wavelength detection revealed the variation of different hydrophilic GUVs and enabled inferences of the "daughter" vesicle composition and the "parent" membrane's local composition during deformation; the "daughter" vesicle composition of OA was lower than that of the "parents", and lateral movement of OA was the primary contributor to the formation of the "daughter" vesicles. Thus, our findings and the newly developed methodology, named in situ quantitative membrane property-morphology relation (QmPMR) analysis, would provide new insights into cell deformation and accelerate research on both deformation and its related events, such as budding and birthing.

Kinetic control of liposome size by direct lipid transfer

Spontaneous lipid vesiculation and related size distribution are traditionally studied in the framework of equilibrium thermodynamics and continuum mechanics, overlooking the kinetic aspects of the process. In the scenario of liposomes consisting of different lipid molecules dispersed in the same medium - a non-equilibrium situation -, the system evolves driven by lipid monomer transfer among the different liposomes. This process encompasses time-dependent changes in liposome size and size distribution, thus predicting size and composition at a given time would entail the control of the size of liposomes by kinetic means, an asset in the framework of diagnostics and synthetic biology. We introduce a direct transfer model, based on the fact that monomers are highly reactive species and apply it to saturated phospholipid molecules differing in hydrophobic chain length. Considering a well-defined gamma-type liposome size distribution, we demonstrate a clear liposome size-composition correlation and are able to predict liposome size and size distribution at any time in the transfer process. The size-composition correlation opens up new prospects for the control of the self-assembling properties of lipids and thereby the control of the liposome size.

Free Energy Surface and Molecular Mechanism of Slow Structural Transitions in Lipid Bilayers

Lipid membrane remodeling, crucial for many cellular processes, is governed by the coupling of membrane structure and shape fluctuations. Given the importance of the ∼ nm length scale, details of the transition intermediates for conformational change are not fully captured by a continuum-mechanical description. Slow dynamics and the lack of knowledge of reaction coordinates (RCs) for biasing methods pose a challenge for all-atom (AA) simulations. Here, we map system dynamics on Langevin dynamics in a normal mode space determined from an elastic network model representation for the lipid-water Hamiltonian. AA molecular dynamics (MD) simulations are used to determine model parameters, and Langevin dynamics predictions for bilayer structural, mechanical, and dynamic properties are validated against MD simulations and experiments. Transferability to describe the dynamics of a larger lipid bilayer and a heterogeneous membrane-protein system is assessed. A set of generic RCs for pore formation in two tensionless bilayers is obtained by coupling Langevin dynamics to the underlying energy landscape for membrane deformations. Structure evolution is carried out by AA MD, wherein the generic RCs are used in a path metadynamics or an umbrella sampling simulation to determine the thermodynamics of pore formation and its molecular determinants, such as the role of distinct bilayer motions, lipid solvation, and lipid packing.

Rotation of the c-Ring Promotes the Curvature Sorting of Monomeric ATP Synthases

ATP synthases are proteins that catalyse the formation of ATP through the rotatory movement of their membrane-spanning subunit. In mitochondria, ATP synthases are found to arrange as dimers at the high-curved edges of cristae. Here, a direct link is explored between the rotatory movement of ATP synthases and their preference for curved membranes. An active curvature sorting of ATP synthases in lipid nanotubes pulled from giant vesicles is found. Coarse-grained simulations confirm the curvature-seeking behaviour of rotating ATP synthases, promoting reversible and frequent protein-protein contacts. The formation of transient protein dimers relies on the membrane-mediated attractive interaction of the order of 1.5 kB T produced by a hydrophobic mismatch upon protein rotation. Transient dimers are sustained by a conic-like arrangement characterized by a wedge angle of θ ≈ 50°, producing a dynamic coupling between protein shape and membrane curvature. The results suggest a new role of the rotational movement of ATP synthases for their dynamic self-assembly in biological membranes.

Origin of the nonlinear structural and mechanical properties in oppositely curved lipid mixtures

Structural and mechanical properties of membranes such as thickness, tail order, bending modulus and curvature energetics play crucial role in controlling various cellular functions that depend on the local lipid organization and membrane reshaping. While behavior of these biophysical properties are well understood in single component membranes, very little is known about how do they change in the mixed lipid membranes. Often various properties of the mixed lipid bilayers are assumed to change linearly with the mole fractions of the constituent lipids which, however, is true for "ideal" mixing only. In this study, using molecular dynamics simulations, we show that structural and mechanical properties of binary lipid mixture change nonlinearly with the lipid mole fractions, and the strength of the nonlinearity depends on two factors - spontaneous curvature difference and locally inhomogeneous interactions between the lipid components.

Mechanical Responses of a Single Myelin Layer: A Molecular Simulation Study

The myelin sheath provides insulation to the brain's neuron cells, which aids in signal transmission and communication with the body. Degenerated myelin hampers the connection between the glial cells, which are the front row responders during traumatic brain injury mitigation. Thus, the structural integrity of the myelin layer is critical for protecting the brain tissue from traumatic injury. At the molecular level, myelin consists of a lipid bilayer, myelin basic proteins (MBP), proteolipid proteins (PLP), water and ions. Structurally, the myelin sheath is formed by repeatedly wrapping forty or more myelin layers around an axon. Here, we have used molecular dynamic simulations to model and capture the tensile response of a single myelin layer. An openly available molecular dynamic solver, LAMMPS, was used to conduct the simulations. The interatomic potentials for the interacting atoms and molecules were defined using CHARMM force fields. Following a standard equilibration process, the molecular model was stretched uniaxially at a deformation rate of 5 Å/ps. We observed that, at around 10% applied strain, the myelin started to cohesively fail via flaw formation inside the bilayers. Further stretching led to a continued expansion of the defect inside the bilayer, both radially and transversely. This study provides the cellular-level mechanisms of myelin damage due to mechanical load.

In-plane and out-of-plane gigahertz sound velocities of saturated and unsaturated phospholipid bilayers from cryogenic to room temperatures

Here, we examined the gigahertz sound velocities of hydrated multibilayers of saturated (1,2-dimyristoyl-sn-glycero-3-phosphocholine, DMPC) and unsaturated (1,2-dioleoyl-sn-glycero-3-phosphocholine, DOPC) phospholipids by Brillouin spectroscopy. Out-of-plane and in-plane (lateral) phonons were studied independently of each other. Similar strong temperature dependences of the sound velocities were found for phonons of both types. The sound velocities in the low-temperature limit were two-fold higher than that at physiological temperatures; a significant part of the changes in sound velocity occurs in the solid-like gel phase. The factors that may be involved in the peculiar behavior of sound velocity include changes in the chain conformational state, relaxation susceptibility, changes in the elastic modulus at infinite frequencies, and lateral packing of molecules.

Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes

Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.

Revealing the Impacts of Chemical Complexity on Submicrometer Sea Spray Aerosol Morphology

Sea spray aerosol (SSA) ejected through bursting bubbles at the ocean surface is a complex mixture of salts and organic species. Submicrometer SSA particles have long atmospheric lifetimes and play a critical role in the climate system. Composition impacts their ability to form marine clouds, yet their cloud-forming potential is difficult to study due to their small size. Here, we use large-scale molecular dynamics (MD) simulations as a "computational microscope" to provide never-before-seen views of 40 nm model aerosol particles and their molecular morphologies. We investigate how increasing chemical complexity impacts the distribution of organic material throughout individual particles for a range of organic constituents with varying chemical properties. Our simulations show that common organic marine surfactants readily partition between both the surface and interior of the aerosol, indicating that nascent SSA may be more heterogeneous than traditional morphological models suggest. We support our computational observations of SSA surface heterogeneity with Brewster angle microscopy on model interfaces. These observations indicate that increased chemical complexity in submicrometer SSA leads to a reduced surface coverage by marine organics, which may facilitate water uptake in the atmosphere. Our work thus establishes large-scale MD simulations as a novel technique for interrogating aerosols at the single-particle level.

Neutron scattering studies on dynamics of lipid membranes

Neutron scattering methods are powerful tools for the study of the structure and dynamics of lipid bilayers in length scales from sub Å to tens to hundreds nm and the time scales from sub ps to μs. These techniques also are nondestructive and, perhaps most importantly, require no additives to label samples. Because the neutron scattering intensities are very different for hydrogen- and deuterium-containing molecules, one can replace the hydrogen atoms in a molecule with deuterium to prepare on demand neutron scattering contrast without significantly altering the physical properties of the samples. Moreover, recent advances in neutron scattering techniques, membrane dynamics theories, analysis tools, and sample preparation technologies allow researchers to study various aspects of lipid bilayer dynamics. In this review, we focus on the dynamics of individual lipids and collective membrane dynamics as well as the dynamics of hydration water.

Compressible viscoelasticity of cell membranes determined by gigahertz-frequency acoustic vibrations

Membrane viscosity is an important property of cell biology, which determines cellular function, development and disease progression. Various experimental and computational methods have been developed to investigate the mechanics of cells. However, there have been no experimental measurements of the membrane viscosity at high-frequencies in live cells. High frequency measurements are important because they can probe viscoelastic effects. Here, we investigate the membrane viscosity at gigahertz-frequencies through the damping of the acoustic vibrations of gold nanoplates. The experiments are modeled using a continuum mechanics theory which reveals that the membranes display viscoelasticity, with an estimated relaxation time of ca. 5.7 + 2.4 / - 2.7 ps. We further demonstrate that membrane viscoelasticity can be used to differentiate a cancerous cell line (the human glioblastoma cells LN-18) from a normal cell line (the mouse brain microvascular endothelial cells bEnd.3). The viscosity of cancerous cells LN-18 is lower than that of healthy cells bEnd.3 by a factor of three. The results indicate promising applications of characterizing membrane viscoelasticity at gigahertz-frequency in cell diagnosis.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Membrane free-energy landscapes derived from atomistic dynamics explain nonuniversal cholesterol-induced stiffening

All lipid membranes have inherent morphological preferences and resist deformation. Yet adaptations in membrane shape can and do occur at multiple length scales. While this plasticity is crucial for cellular physiology, the factors controlling the morphological energetics of lipid bilayers and the dominant mechanisms of membrane remodeling remain unclear. An ongoing debate regarding the universality of the stiffening effect of cholesterol underscores the challenges facing this field, both experimentally and theoretically, even for simple lipid mixtures. On the computational side, we have argued that enhanced- sampling all-atom molecular dynamics simulations are uniquely suited for quantification of membrane conformational energetics, not only because they minimize a-priori assumptions, but also because they permit analysis of bilayers in deformed states. To showcase this approach, we examine reported inconsistencies between alternative experimental measurements of bending moduli for cholesterol-enriched membranes. Specifically, we analyze lipid bilayers with different chain saturation, and compute free-energy landscapes for curvature deformations distributed over areas from ∼5 to ∼60 nm 2 . These enhanced simulations, totaling over 100 microseconds of sampling time, enable us to directly quantify both bending and tilt moduli, and to dissect the contributing factors and molecular mechanisms of curvature generation at each length scale. Our results show that cholesterol effects are lipid-specific, in agreement with giantvesicle measurements, and explain why experiments probing nanometer scale lipid dynamics diverge. In summary, we demonstrate that quantitative structure-mechanics relationships can now be established for heterogenous membranes, paving the way for addressing open fundamental questions in cell membrane mechanics.

Significance

Elucidating the energetics and mechanisms of membrane remodeling is an essential step towards understanding cell physiology. This problem is challenging, however, because membrane bending involves both large-scale and atomic-level dynamics, which are difficult to measure simultaneously. A recent controversy regarding the stiffening effect of cholesterol, which is ubiquitous in animal cells, illustrates this challenge. We show how enhanced molecular-dynamics simulations can bridge this length-scale gap and reconcile seemingly incongruent observations. This approach facilitates a conceptual connection between lipid chemistry and membrane mechanics, thereby providing a solid basis for future research on remodeling phenomena, such as in membrane trafficking or viral infection.

Modelling the role of membrane mechanics in cell adhesion on titanium oxide nanotubes

Titanium surface treated with titanium oxide nanotubes was used in many studies to quantify the effect of surface topography on cell fate. However, the predicted optimal diameter of nanotubes considerably differs among studies. We propose a model that explains cell adhesion to a nanostructured surface by considering the deformation energy of cell protrusions into titanium nanotubes and the adhesion to the surface. The optimal surface topology is defined as a geometry that gives the membrane a minimum energy shape. A dimensionless parameter, the cell interaction index, was proposed to describe the interplay between the cell membrane bending, the intrinsic curvature, and the strength of cell adhesion. Model simulation shows that an optimal nanotube diameter ranging from 20 nm to 100 nm (cell interaction index between 0.2 and 1, respectively) is feasible within a certain range of parameters describing cell membrane adhesion and bending. The results indicate a possibility to tune the topology of a nanostructural surface in order to enhance the proliferation and differentiation of cells mechanically compatible with the given surface geometry while suppressing the growth of other mechanically incompatible cells.

The impact of lipid polyunsaturation on the physical and mechanical properties of lipid membranes

The lipid composition of cellular membranes and the balance between the different lipid components can be impacted by aging, certain pathologies, specific diets and other factors. This is the case in a subgroup of individuals with psychiatric disorders, such as schizophrenia, where cell membranes of patients have been shown to be deprived in polyunsaturated fatty acids (PUFAs), not only in brain areas where the target receptors are expressed but also in peripheral tissues. This PUFA deprivation thus represents a biomarker of such disorders that might impact not only the interaction of antipsychotic medications with these membranes but also the activation and signaling of the targeted receptors embedded in the lipid membrane. Therefore, it is crucial to understand how PUFAs levels alterations modulate the different physical properties of membranes. In this paper, several biophysical approaches were combined (Laurdan fluorescence spectroscopy, atomic force microscopy, differential scanning calorimetry, molecular modeling) to characterize membrane properties such as fluidity, elasticity and thickness in PUFA-enriched cell membranes and lipid model systems reflecting the PUFA imbalance observed in some diseases. The impact of both the number of unsaturations and their position along the chain on the above properties was investigated. Briefly, data revealed that PUFA presence in membranes increases membrane fluidity, elasticity and flexibility and decreases its thickness and order parameter. Both the level of unsaturation and their position affect these membrane properties.

Capturing Lipid Nanodisc Shape and Properties Using a Continuum Elastic Theory

Lipid nanodiscs are nanometric bilayer patches enveloped by confining structures, commonly composed of membrane scaffolding proteins (MSPs). To resolve the interplay between MSP geometry, lipid confinement, and membrane material properties on the nanodisc shape, we apply a continuum elastic theory accounting for lipid bending, tilting, and area deformations. The equilibrium nanodisc shape is then determined by minimizing the elastic free energy functional. Analytic expressions derived under simplifying assumptions demonstrate that the nanodisc shape is sensitive to its size, lipid density, and the lipid tilt and thickness imposed at the contact with the MSP. Under matching physical parameters, these expressions quantitatively reproduce the shape of nanodiscs seen in molecular dynamics simulations, but only if lipid tilt is explicitly considered. We further demonstrate how the bending rigidity can be extracted from the membrane shape profile by fitting the numerically minimized full elastic functional to the membrane shape found in simulations. This fitting procedure faithfully informs on the bending rigidity of nanodiscs larger than ca. 5 nm in radius. The fitted profiles accurately reproduce the increase in bending modulus found using real-space fluctuation analysis of simulated nanodiscs and, for large nanodiscs, also accurately resolve its spatial variations. Our study shows how deformations in lipid patches confined in nanodiscs can be well described by a continuum elastic theory and how this fit can be used to determine local material properties from shape analysis of nanodiscs in simulations. This methodology could potentially allow direct determination of lipid properties from experiments, for example cryo-electron microscopy images of lipid nanodiscs, thereby allowing to guide the development of future nanodisc formulations with desired properties.

Nanoscale Bending Dynamics in Mixed-Chain Lipid Membranes

Lipids that have two tails of different lengths are found throughout biomembranes in nature, yet the effects of this asymmetry on the membrane properties are not well understood, especially when it comes to the membrane dynamics. Here we study the nanoscale bending fluctuations in model mixed-chain 14:0–18:0 PC (MSPC) and 18:0–14:0 PC (SMPC) lipid bilayers using neutron spin echo (NSE) spectroscopy. We find that despite the partial interdigitation that is known to persist in the fluid phase of these membranes, the collective fluctuations are enhanced on timescales of tens of nanoseconds, and the chain-asymmetric lipid bilayers are softer than an analogous chain-symmetric lipid bilayer with the same average number of carbons in the acyl tails, di-16:0 PC (DPPC). Quantitative comparison of the NSE results suggests that the enhanced bending fluctuations at the nanosecond timescales are consistent with experimental and computational studies that showed the compressibility moduli of chain-asymmetric lipid membranes are 20% to 40% lower than chain-symmetric lipid membranes. These studies add to growing evidence that the partial interdigitation in mixed-chain lipid membranes is highly dynamic in the fluid phase and impacts membrane dynamic processes from the molecular to mesoscopic length scales without significantly changing the bilayer thickness or area per lipid.

Biophysics of Membrane Stiffening by Cholesterol and Phosphatidylinositol 4,5-bisphosphate (PIP2)

Cell membranes regulate a wide range of phenomena that are implicated in key cellular functions. Cholesterol, a critical component of eukaryotic cell membranes, is responsible for cellular organization, membrane elasticity, and other critical physicochemical parameters. Besides cholesterol, other lipid components such as phosphatidylinositol 4,5-bisphosphate (PIP2) are found in minor concentrations in cell membranes yet can also play a major regulatory role in various cell functions. In this chapter, we describe how solid-state deuterium nuclear magnetic resonance (²H NMR) spectroscopy together with neutron spin-echo (NSE) spectroscopy can inform synergetic changes to lipid molecular packing due to cholesterol and PIP2 that modulate the bending rigidity of lipid membranes. Fundamental structure-property relations of molecular self-assembly are illuminated and point toward a length and time-scale dependence of cell membrane mechanics, with significant implications for biological activity and membrane lipid-protein interactions.

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

Changes in Mechanical Properties of Vesicles by Mucin in Aqueous Solution

The mechanical properties of vesicles were investigated as they were prepared, according to the ratio of mucin to dipalmitoylphosphatidylcholine (DPPC), using an atomic force microscope (AFM). After the confirmation of the vesicle adsorption on a mica surface, an AFM-tip deflection, caused by the interaction between the tip and the vesicle, was measured. The deflection showed that the tip broke through into the vesicle twice. Each break meant a tip-penetration into the upper and lower portion of the vesicle. Only the first penetration allowed the Hertzian model available to estimate the vesicle mechanical moduli. Two moduli reduced as the ratio of mucin to DPPC increased to 0.5, but the moduli were little changed above the 0.5 ratio. These results seem to be a platform for the effect of the mucin on the plasma-membrane anchoring and cellular signaling.

Polymyxins induce lipid scrambling and disrupt the homeostasis of Gram-negative bacteria membrane

Polymyxins are increasingly used as the last-line therapeutic option for the treatment of infections caused by multidrug-resistant Gram-negative bacteria. However, efforts to address the resistance in superbugs are compromised by a poor understanding of the bactericidal modes because high-resolution detection of the cell structure is still lacking. By performing molecular dynamics simulations at a coarse-grained level, here we show that polymyxin B (PmB) disrupts Gram-negative bacterial membranes by altering lipid homeostasis and asymmetry. We found that the binding of PmBs onto the asymmetric outer membrane (OM) loosens the packing of lipopolysaccharides (LPS) and induces unbalanced bending torque between the inner and outer leaflets, which in turn triggers phospholipids to flip from the inner leaflet to the outer leaflet to compensate for the stress deformation. Meanwhile, some LPSs may be detained on the inner membrane (IM). Then, the lipid-scrambled OM undergoes phase separation. Defects are created at the boundaries between LPS-rich domains and phospholipid-rich domains, which consequently facilitate the uptake of PmB across the OM. Finally, PmBs target LPSs detained on the IM and similarly perturb the IM. This lipid Scramble, membrane phase Separation, and peptide Translocation model depicts a novel mechanism by which polymyxins kill bacteria and sheds light on developing a new generation of polymyxins or antibiotic adjuvants with improved killing activities and higher therapeutic indices.

Structural and mechanical properties of the red blood cell's cytoplasmic membrane seen through the lens of biophysics

Red blood cells (RBCs) are the most abundant cell type in the human body and critical suppliers of oxygen. The cells are characterized by a simple structure with no internal organelles. Their two-layered outer shell is composed of a cytoplasmic membrane (RBC cm ) tethered to a spectrin cytoskeleton allowing the cell to be both flexible yet resistant against shear stress. These mechanical properties are intrinsically linked to the molecular composition and organization of their shell. The cytoplasmic membrane is expected to dominate the elastic behavior on small, nanometer length scales, which are most relevant for cellular processes that take place between the fibrils of the cytoskeleton. Several pathologies have been linked to structural and compositional changes within the RBC cm and the cell's mechanical properties. We review current findings in terms of RBC lipidomics, lipid organization and elastic properties with a focus on biophysical techniques, such as X-ray and neutron scattering, and Molecular Dynamics simulations, and their biological relevance. In our current understanding, the RBC cm 's structure is patchy, with nanometer sized liquid ordered and disordered lipid, and peptide domains. At the same time, it is surprisingly soft, with bending rigidities κ of 2-4 kBT. This is in strong contrast to the current belief that a high concentration of cholesterol results in stiff membranes. This extreme softness is likely the result of an interaction between polyunsaturated lipids and cholesterol, which may also occur in other biological membranes. There is strong evidence in the literature that there is no length scale dependence of κ of whole RBCs.

Partial Volumes of Phosphatidylcholines and Vitamin E: α-Tocopherol Prefers Disordered Membranes

Despite its discovery over 95 years ago, the biological and nutritional roles of vitamin E remain subjects of much controversy. Though it is known to possess antioxidant properties, recent assertions have implied that vitamin E may not be limited to this function in living systems. Through densitometry measurements and small-angle X-ray scattering we observe favorable interactions between α-tocopherol and unsaturated phospholipids, with more favorable interactions correlating to an increase in lipid chain unsaturation. Our data provide evidence that vitamin E may preferentially associate with oxygen sensitive lipids─an association that is considered innate for a viable membrane antioxidant.

Development of resistance to 5-fluorouracil affects membrane viscosity and lipid composition of cancer cells

The investigations reported here were designed to determine whether the bulk plasma membrane is involved in mechanisms of acquired resistance of colorectal cancer cells to 5-fluorouracil (5-FU). Fluorescence lifetime imaging microscopy (FLIM) of live cultured cells stained with viscosity-sensitive probe BODIPY 2 was exploited to non-invasively assess viscosity in the course of treatment and adaptation to the drug. In parallel, lipid composition of membranes was examined with the time-of-flight secondary ion mass spectrometry (ToF-SIMS). Our results showed that a single treatment with 5-FU induced only temporal changes of viscosity in 5-FU sensitive cells immediately after adding the drug. Acquisition of chemoresistance was accompanied by persistent increase of viscosity, which was preserved upon treatment without any changes. Lipidomic analysis revealed that the resistant cells had a lower level of monounsaturated fatty acids and increased sphingomyelin or decreased phosphatidylcholine in their membranes, which partly explain increase of the viscosity. Thus, we propose that a high membrane viscosity mediates the acquisition of resistance to 5-FU.

The effect of cholesterol on highly curved membranes measured by nanosecond Fluorescence Correlation Spectroscopy

Small lipid vesicles (with diameter ≤ 100nm) with their highly curved membranes comprise a special class of biological lipid bilayers. The mechanical properties of such membranes are critical for their function, e.g. exocytosis. Cholesterol is a near-universal regulator of membrane properties in animal cells. Yet measurements of the effect of cholesterol on the mechanical properties of membranes have remained challenging, and the interpretation of such measurements has remained a matter of debate. Here we show that nanosecond fluorescence correlation spectroscopy (FCS) can directly measure the ns-microsecond rotational correlation time (τr) of a lipid probe in high curvature vesicles with extraordinary sensitivity. Using a home-built 4-Pi fluorescence cross-correlation spectrometer containing polarization-modulating elements, we measure the rotational correlation time (τr) of Nile Red in neurotransmitter vesicle mimics. As the cholesterol mole fraction increases from 0 to 50 %, τr increases from 17 ± 1 to 112 ± 12 ns, indicating a viscosity change of nearly a factor of 7. These measurements are corroborated by solid-state NMR results, which show that the order parameter of the lipid acyl chains increases by about 50% for the same change in cholesterol concentration. Additionally, we measured the spectral parameters of polarity-sensitive fluorescence dyes, which provide an indirect measure of viscosity. The green/red ratio of Nile Red and the generalized polarization of Laurdan show consistent increases of 1.3x and 2.6x, respectively. Our results demonstrate that rotational FCS can directly measure the viscosity of highly curved membranes with higher sensitivity and wider dynamic range compared to other conventional techniques. Significantly, we observe that the viscosity of neurotransmitter vesicle mimics is remarkably sensitive to their cholesterol content.

The bending rigidity of the red blood cell cytoplasmic membrane

An important mechanical property of cells is the membrane bending modulus, κ. In the case of red blood cells (RBCs) there is a composite membrane consisting of a cytoplasmic membrane and an underlying spectrin network. Literature values of κ are puzzling, as they are reported over a wide range, from 5 kBT to 230 kBT. To disentangle the contribution of the cytoplasmic membrane from the spectrin network, we investigated the bending of red blood cell cytoplasmic membranes (RBCcm) in the absence of spectrin and adenosine triphosphate (ATP). We used a combination of X-ray diffuse scattering (XDS), neutron spin-echo (NSE) spectrometry and Molecular Dynamics (MD) simulations. Our results indicate values of κ of order 4 kBT to 6 kBT, relatively small compared to literature values for most single component lipid bilayers. We suggest two ways this relative softness might confer biological advantage.

Multiscale lipid membrane dynamics as revealed by neutron spectroscopy

The plasma membrane is one of the principal structural components of the cell and, therefore, one of the key components of the cellular life. Because the membrane's dynamics links the membrane's structure and function, the complexity and the broad range of the membrane's motions are essential for the enormously diverse functionality of the cell membrane. Even for the main membrane component, the lipid bilayer, considered alone, the range and complexity of the lipid motions are remarkable. Spanning the time scale from sub-picosecond to minutes and hours, the lipid motion in a bilayer is challenging to study even when a broad array of dynamic measurement techniques is employed. Neutron scattering plays a special role among such dynamic measurement techniques, particularly, because it involves the energy transfers commensurate with the typical intra- and inter- molecular dynamics and the momentum transfers commensurate with intra- and inter-molecular distances. Thus, using neutron scattering-based techniques, the spatial and temporal information on the lipid motion can be obtained and analysed simultaneously. Protium vs. deuterium sensitivity and non-destructive character of the neutron probe add to the remarkable prowess of neutron scattering for elucidating the lipid dynamics. Herein we present an overview of the neutron scattering-based studies of lipid dynamics in model membranes, with a discussion of the direct relevance and implications to the real-life cell membranes. The latter are much more complex systems than simple model membranes, consisting of heterogeneous non-stationary domains composed of lipids, proteins, and other small molecules, such as carbohydrates. Yet many fundamental aspects of the membrane behavior and membrane interactions with other molecules can be understood from neutron scattering measurements of the model membranes. For example, such studies can provide a great deal of information on the interactions of antimicrobial compounds with the lipid matrix of a pathogen membrane, or the interactions of drug molecules with the plasma membrane. Finally, we briefly discuss the recently emerging field of neutron scattering membrane studies with a reach far beyond the model membrane systems.

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.

Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale

Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.

A vesicle microrheometer for high-throughput viscosity measurements of lipid and polymer membranes

Viscosity is a key property of cell membranes that controls mobility of embedded proteins and membrane remodeling. Measuring it is challenging because existing approaches involve complex experimental designs and/or models, and the applicability of some methods is limited to specific systems and membrane compositions. As a result there is scarcity of systematic data, and the reported values for membrane viscosity vary by orders of magnitude for the same system. Here, we show how viscosity of membranes can be easily obtained from the transient deformation of giant unilamellar vesicles. The approach enables a noninvasive, probe-independent, and high-throughput measurement of the viscosity of membranes made of lipids or polymers with a wide range of compositions and phase state. Using this novel method, we have collected a significant amount of data that provides insights into the relation between membrane viscosity, composition, and structure.

Relationship between Viscosity and Acyl Tail Dynamics in Lipid Bilayers

Membrane viscosity is a fundamental property that controls molecular transport and structural rearrangements in lipid membranes. Given its importance in many cell processes, various experimental and computational methods have been developed to measure the membrane viscosity, yet the estimated values depend highly on the method and vary by orders of magnitude. Here we investigate the molecular origins of membrane viscosity by measuring the nanoscale dynamics of the lipid acyl tails using x-ray and neutron spectroscopy techniques. The results show that the membrane viscosity can be estimated from the structural relaxation times of the lipid tails.

The dynamic face of lipid membranes

Cell membranes - primarily composed of lipids, sterols, and proteins - form a dynamic interface between living cells and their environment. They act as a mechanical barrier around the cell while selectively facilitating material transport, signal transduction, and various other functions necessary for the cell viability. The complex functionality of cell membranes and the hierarchical motions and responses they exhibit demand a thorough understanding of the origin of different membrane dynamics and how they are influenced by molecular additives and environmental cues. These dynamic modes include single-molecule diffusion, thermal fluctuations, and large-scale membrane deformations, to name a few. This review highlights advances in investigating structure-driven dynamics associated with model cell membranes, with a particular focus on insights gained from neutron scattering and spectroscopy experiments. We discuss the uniqueness of neutron contrast variation and its remarkable potential in probing selective membrane structure and dynamics on spatial and temporal scales over which key biological functions occur. We also present a summary of current and future opportunities in synergistic combinations of neutron scattering with molecular dynamics (MD) simulations to gain further understanding of the molecular mechanisms underlying complex membrane functions.

Polystyrene perturbs the structure, dynamics, and mechanical properties of DPPC membranes: An experimental and computational study

Synthetic plastic oligomers can interact with the cells of living organisms by different ways. They can be intentionally administered to the human body as part of nanosized biomedical devices. They can be inhaled by exposed workers, during the production of multicomponent, polymer-based nanocomposites. They can leak out of food packaging. Most importantly, they can result from the degradation of plastic waste, and enter the food chain. A physicochemical characterization of the effects of synthetic polymers on the structure and dynamics of cell components is still lacking. Here, we combine a wide spectrum of experimental techniques (calorimetry, x-ray, and neutron scattering) with atomistic Molecular Dynamics simulations to study the interactions between short chains of polystyrene (25 monomers) and model lipid membranes (DPPC, in both gel and fluid phase). We find that doping doses of polystyrene oligomers alter the thermal properties of DPPC, stabilizing the fluid lipid phase. They perturb the membrane structure and dynamics, in a concentration-dependent fashion. Eventually, they modify the mechanical properties of DPPC, reducing its bending modulus in the fluid phase. Our results call for a systematic, interdisciplinary assessment of the mechanisms of interaction of synthetic, everyday use polymers with cell membranes.

Experimental demonstration of the novel "van-Hove integral method (vHI)" for measuring diffusive dynamics by elastic neutron scattering

Quasi-elastic neutron scattering (QENS)-based on the seminal work of Nobel Laureate Brockhouse-has been one of the major methods for studying pico-second to nano-second diffusive dynamics over the past 70 years. This is regarded as an "inelastic" method for dynamics. In contrast, we recently proposed a new neutron-scattering method for dynamics, which uses the elastic line of the scattering to access system dynamics directly in the time domain (Benedetto and Kearley in Sci Rep 9:11284, 2019). This new method has been denoted "vHI" that stands for "van Hove Integral". The reason is that, under certain conditions, the measured elastic intensity corresponds to the running-time integral of the intermediate scattering function, [Formula: see text], up to a time that is inversely proportional to the energy band-width incident on the sample. As a result, [Formula: see text] is accessed from the time derivative of the measured vHI profile. vHI has been supported by numerical and Monte-Carlo simulations, but has been difficult to validate experimentally due to the lack of a suitable instrument. Here we show that vHI works in practice, which we achieved by using a simple modification to the standard QENS backscattering spectrometer methodology. Basically, we varied the neutron-energy band-widths incident at the sample via a step-wise variation of the frequency of the monochromator Doppler-drive. This provides a measurement of the vHI profile at the detectors. The same instrument and sample were also used in standard QENS mode for comparison. The intermediate scattering functions, [Formula: see text], obtained by the two methods-vHI and QENS-are strikingly similar providing a direct experimental validation of the vHI method. Perhaps surprisingly, the counting statistics of the two methods are comparable even though the instrument used was expressly designed for QENS. This shows that the methodology modification adopted here can be used in practice to access vHI profiles at many of the backscattering spectrometers worldwide. We also show that partial integrations of the measured QENS spectrum cannot provide the vHI profile, which clarifies a common misconception. At the same time, we show a novel approach which does access [Formula: see text] from QENS spectra.

Stiffening Effect of the [Bmim][Cl] Ionic Liquid on the Bending Dynamics of DMPC Lipid Vesicles

The elastic properties of the cellular lipid membrane play a crucial role for life. Their alteration can lead to cell malfunction, and in turn, being able to control them holds the promise of effective therapeutic and diagnostic approaches. In this context, due to their proven strong interaction with lipid bilayers, ionic liquids (ILs)-a vast class of organic electrolytes-may play an important role. This work focuses on the effect of the model imidazolium-IL [bmim][Cl] on the bending modulus of DMPC lipid vesicles, a basic model of cellular lipid membranes. Here, by combining small-angle neutron scattering and neutron spin-echo spectroscopy, we show that the IL, dispersed at low concentrations at the bilayer-water interface, (i) diffuses into the lipid region, accounting for five IL-cations for every 11 lipids, and (ii) causes an increase of the lipid bilayer bending modulus, up to 60% compared to the neat lipid bilayer at 40 °C.

Collective dynamics in lipid membranes containing transmembrane peptides

Biological membranes are composed of complex mixtures of lipids and proteins that influence each other's structure and function. The biological activities of many channel-forming peptides and proteins are known to depend on the material properties of the surrounding lipid bilayer. However, less is known about how membrane-spanning channels affect the lipid bilayer properties, and in particular, their collective fluctuation dynamics. Here we use neutron spin echo spectroscopy (NSE) to measure the collective bending and thickness fluctuation dynamics in dimyristoylphosphatidylcholine (di 14 : 0 PC, DMPC) lipid membranes containing two different antimicrobial peptides, alamethicin (Ala) and gramicidin (gD). Ala and gD are both well-studied antimicrobial peptides that form oligomeric membrane-spanning channels with different structures. At low concentrations, the peptides did not have a measurable effect on the average bilayer structure, yet significantly changed the collective membrane dynamics. Despite both peptides forming transmembrane channels, they had opposite effects on the relaxation time of the collective bending fluctuations and associated effective bending modulus, where gD addition stiffened the membrane while Ala addition softened the membrane. Meanwhile, the lowest gD concentrations enhanced the collective thickness fluctuation dynamics, while the higher gD concentrations and all studied Ala concentrations dampened these dynamics. The results highlight the synergy between lipids and proteins in determining the collective membrane dynamics and that not all peptides can be universally treated as rigid bodies when considering their effects on the lipid bilayer fluctuations.

Biomembrane Structure and Material Properties Studied With Neutron Scattering

Cell membranes and their associated structures are dynamical supramolecular structures where different physiological processes take place. Detailed knowledge of their static and dynamic structures is therefore needed, to better understand membrane biology. The structure–function relationship is a basic tenet in biology and has been pursued using a range of different experimental approaches. In this review, we will discuss one approach, namely the use of neutron scattering techniques as applied, primarily, to model membrane systems composed of lipid bilayers. An advantage of neutron scattering, compared to other scattering techniques, is the differential sensitivity of neutrons to isotopes of hydrogen and, as a result, the relative ease of altering sample contrast by substituting protium for deuterium. This property makes neutrons an ideal probe for the study of hydrogen-rich materials, such as biomembranes. In this review article, we describe isotopic labeling studies of model and viable membranes, and discuss novel applications of neutron contrast variation in order to gain unique insights into the structure, dynamics, and molecular interactions of biological membranes. We specifically focus on how small-angle neutron scattering data is modeled using different contrast data and molecular dynamics simulations. We also briefly discuss neutron reflectometry and present a few recent advances that have taken place in neutron spin echo spectroscopy studies and the unique membrane mechanical data that can be derived from them, primarily due to new models used to fit the data.

Value of models for membrane budding

The budding of membranes and curvature generation is common to many forms of trafficking in cells. Clathrin-mediated endocytosis, as a prototypical example of trafficking, has been studied in great detail using a variety of experimental systems and methods. Recently, advances in experimental methods have led to great strides in insights on the molecular mechanisms and the spatiotemporal dynamics of the protein machinery associated with membrane curvature generation. These advances have been ably supported by computational models, which have given us insights into the underlying mechanical principles of clathrin-mediated endocytosis. On the other hand, targeted experimental perturbation of membranes has lagged behind that of proteins in cells. In this area, modeling is especially critical to interpret experimental measurements in a mechanistic context. Here, we discuss the contributions made by these models to our understanding of endocytosis and identify opportunities to strengthen the connections between models and experiments.

Effect of gold nanoparticle incorporation into oil-swollen surfactant lamellar membranes

An oil-swollen surfactant membrane is employed to measure the effects of incorporated hydrophobically functionalized gold nanoparticles (AuNPs) on the structure and dynamics of the membranes. While maintaining an average AuNP diameter of approximately 5 nm, the membrane thickness was varied from 5 nm to 7.5 nm by changing the amount of oil in the membrane. The membranes become softer as the proportion of oil is increased, while the thickness fluctuations become slower. We attribute this to an increased fluctuation wavelength. Incorporation of AuNPs in the membrane induces membrane thinning and softening. Oil molecules surround the nanoparticles in the membrane and help their relatively homogeneous distribution. AuNPs significantly alter the membrane's structure and dynamics through thinning of the membrane, increased compressibility, and possible diffusion of AuNPs inside the membrane.