The molecular structure of a phosphatidylserine bilayer determined by scattering and molecular dynamics simulations

Small-angle X-ray and neutron scattering applied to lipid-based nanoparticles: Recent advancements across different length scales

Lipid-based nanoparticles (LNPs), ranging from nanovesicles to non-lamellar assemblies, have gained significant attention in recent years, as versatile carriers for delivering drugs, vaccines, and nutrients. Small-angle scattering methods, employing X-rays (SAXS) or neutrons (SANS), represent unique tools to unveil structure, dynamics, and interactions of such particles on different length scales, spanning from the nano to the molecular scale. This review explores the state-of-the-art on scattering methods applied to unveil the structure of lipid-based nanoparticles and their interactions with drugs and bioactive molecules, to inform their rational design and formulation for medical applications. We will focus on complementary information accessible with X-rays or neutrons, ranging from insights on the structure and colloidal processes at a nanoscale level (SAXS) to details on the lipid organization and molecular interactions of LNPs (SANS). In addition, we will review new opportunities offered by Time-resolved (TR)-SAXS and -SANS for the investigation of dynamic processes involving LNPs. These span from real-time monitoring of LNPs structural evolution in response to endogenous or external stimuli (TR-SANS), to the investigation of the kinetics of lipid diffusion and exchange upon interaction with biomolecules (TR-SANS). Finally, we will spotlight novel combinations of SAXS and SANS with complementary on-line techniques, recently enabled at Large Scale Facilities for X-rays and neutrons. This emerging technology enables synchronized multi-method investigation, offering exciting opportunities for the simultaneous characterization of the structure and chemical or mechanical properties of LNPs.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

A closer look at calcium-induced interactions between phosphatidylserine-(PS) doped liposomes and the structural effects caused by inclusion of gangliosides or polyethylene glycol- (PEG) modified lipids

The effects of polyethylene glycol- (PEG) modified lipids and gangliosides on the Ca2+ induced interaction between liposomes composed of palmitoyl-oleoyl phosphatidylethanolamine (POPE) and palmitoyl-oleoyl phosphatidylserine (POPS) was investigated at physiological ionic strength. Förster resonance energy transfer (FRET) studies complemented with dynamic light scattering (DLS) and cryo-transmission electron microscopy (Cryo-EM) show that naked liposomes tend to adhere, rupture, and collapse on each other's surfaces upon addition of Ca2+, eventually resulting in the formation of large multilamellar aggregates and bilayer sheets. Noteworthy, the presence of gangliosides or PEGylated lipids does not prevent the adhesion-rupture process, but leads to the formation of small, long-lived bilayer fragments/disks. PEGylated lipids seem to be more effective than gangliosides at stabilizing these structures. Attractive interactions arising from ion correlation are proposed to be a driving force for the liposome-liposome adhesion and rupture processes. The results suggest that, in contrast with the conclusions drawn from previous solely FRET-based studies, direct liposome-liposome fusion is not the dominating process triggered by Ca2+ in the systems studied.

Deciphering Lipid Arrangement in Phosphatidylserine/Phosphatidylcholine Mixed Membranes: Simulations and Experiments

Phosphatidylserine (PS) exposure on the plasma membrane is crucial for many cellular processes including apoptotic cell recognition, blood clotting regulation, cellular signaling, and intercellular interactions. In this study, we investigated the arrangement of PS headgroups in mixed PS/phosphatidylcholine (PC) bilayers, serving as a simplified model of the outer leaflets of mammalian cell plasma membranes. Combining atomistic-scale molecular dynamics (MD) simulations with Langmuir monolayer experiments, we unraveled the mutual miscibility of POPC and POPS lipids and the intricate intermolecular interactions inherent to these membranes as well as the disparities in position and orientation of PC and PS headgroups. Our experiments revealed micrometer-scale miscibility at all mole fractions of POPC and POPS, marked by modest deviations from ideal mixing with no apparent microscale phase separation. The MD simulations, meanwhile, demonstrated that these deviations were due to strong electrostatic interactions between like-lipid pairs (POPC-POPC and POPS-POPS), culminating in lateral segregation and nanoscale clustering. Notably, PS headgroups profoundly affect the ordering of the lipid acyl chains, leading to lipid elongation and subtle PS protrusion above the zwitterionic membrane. In addition, PC headgroups are more tilted with respect to the membrane normal, while PS headgroups align at a smaller angle, making them more exposed to the surface of the mixed PC/PS membranes. These findings provide a detailed molecular-level account of the organization of mixed PC/PS membranes, corroborated by experimental data. The insights gained here extend our comprehension of the physiological role of PSs.

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.

Lipid Corona Formation on Micro- and Nanoplastic Particles Modulates Uptake and Toxicity in A549 Cells

Plastic waste is a global issue leaving no continents unaffected. In the environment, ultraviolet radiation and shear forces in water and land contribute to generating micro- and nanoplastic particles (MNPP), which organisms can easily take up. Plastic particles enter the human food chain, and the accumulation of particles within the human body is expected. Crossing epithelial barriers and cellular uptake of MNPP involves the interaction of plastic particles with lipids. To this end, we generated unilamellar vesicles from POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine) and POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) and incubated them with pristine, carboxylated, or aminated polystyrene spheres (about 1 µm in diameter) to generate lipid coronas around the particles. Lipid coronas enhanced the average particle sizes and partially changed the MNPP zeta potential and polydispersity. In addition, lipid coronas led to significantly enhanced uptake of MNPP particles but not their cytotoxicity, as determined by flow cytometry. Finally, adding proteins to lipid corona nanoparticles further modified MNPP uptake by reducing the uptake kinetics, especially in pristine and carboxylated plastic samples. In conclusion, our study demonstrates for the first time the impact of different types of lipids on differently charged MNPP particles and the biological consequences of such modifications to better understand the potential hazards of plastic exposure.

Gold-Supported Lipid Membranes Formed by Redox-Triggered Vesicle Fusion on Binary Self-Assembled Monolayers: Ion-Pairing Association and Surface Hydrophilicity

The assembly of biomimetic, planar supported lipid bilayers (SLBs) by the popular vesicle fusion method, which relies on the spontaneous adsorption and rupture of small unilamellar vesicles from aqueous solution on a solid surface, typically works with a limited range of support materials and lipid systems. We previously reported a conceptual advance in the formation of SLBs from vesicles in the gel or fluid phase using the interfacial ion-pairing association of charged phospholipid headgroups with electrochemically generated cationic ferroceniums bound to a self-assembled monolayer (SAM) chemisorbed to gold. This redox-driven approach lays down a single bilayer membrane on the SAM-modified gold surface at room temperature within minutes and is compatible with both anionic and zwitterionic phospholipids. The present work explores the effects of the surface ferrocene concentration and hydrophobicity/hydrophilicity on the formation of continuous SLBs of dialkyl phosphatidylserine, dialkyl phosphatidylglycerol, and dialkyl phosphatidylcholine using binary SAMs of ferrocenylundecanethiolate (FcC₁₁S) and dodecanethiolate (CH₃C₁₁S) or hydroxylundecanethiolate (HOC₁₁S) comprising different surface mole fractions of ferrocene (χ_Fc^surf). An increase in the surface hydrophilicity and surface free energy of the FcC₁₁S/HOC₁₁S SAM mitigates the decrease in the attractive ion-pairing interactions resulting from a reduced χ_Fc^surf. SLBs of ≳80% area coverage form on the FcC₁₁S/HOC₁₁S SAM for all the phospholipid types down to χ_Fc^surf of at least 0.2, composition yielding a water contact angle (θ_W) of 44 ± 4°. By contrast, a greater number of ion-pairing interactions is required on the hydrophobic FcC₁₁S/CH₃C₁₁S surface to drive the vesicle fusion process; bilayers or bilayer patches form at χ_Fc^surf ≳ 0.6 (θ_W = 97 ± 3°). These findings will aid in tailoring the surface chemistry of redox-active modified surfaces to widen the conditions that yield supported lipid membranes.

Effect of Anionic Lipids on Mammalian Plasma Cell Membrane Properties

The effect of lipid composition on models of the inner leaflet of mammalian cell membranes has been investigated. Grazing incidence X-ray diffraction and X-ray and neutron reflectivity have been used to characterize lipid packing and solvation, while electrochemical and infrared spectroscopic methods have been employed to probe phase behavior in an applied electric field. Introducing a small quantity of the anionic lipid dimyristoylphosphatidylserine (DMPS) into bilayers of zwitterionic dimyristoylphosphatidylethanolamine (DMPE) results in a significant change in the bilayer response to an applied field: the tilt of the hydrocarbon chains increases before returning to the original tilt angle on detachment of the bilayer. Equimolar mixtures, with slightly closer chain packing, exhibit a similar but weaker response. The latter also tend to incorporate more solvent during this electrochemical phase transition, at levels similar to those of pure DMPS. Reflectivity measurements reveal greater solvation of lipid layers for DMPS > 30 mol %, matching the greater propensity for DMPS-rich bilayers to incorporate water. Taken together, the data indicate that the range of 10-35 mol % DMPS provides optimum bilayer properties (in flexibility and function as a barrier), which may explain why the DMPS content of cell membranes tends to be found within this range.

Lamellarity-Driven Differences in Surface Structural Features of DPPS Lipids: Spectroscopic, Calorimetric and Computational Study

Although single-lipid bilayers are usually considered models of eukaryotic plasma membranes, their research drops drastically when it comes to exclusively anionic lipid membranes. Being a major anionic phospholipid in the inner leaflet of eukaryote membranes, phosphatidylserine-constituted lipid membranes were occasionally explored in the form of multilamellar liposomes (MLV), but their inherent instability caused a serious lack of efforts undertaken on large unilamellar liposomes (LUVs) as more realistic model membrane systems. In order to compensate the existing shortcomings, we performed a comprehensive calorimetric, spectroscopic and MD simulation study of time-varying structural features of LUV made from 1,2-dipalmitoyl-sn-glycero-3-phospho-L-serine (DPPS), whereas the corresponding MLV were examined as a reference. A substantial uncertainty of UV/Vis data of LUV from which only T_m was unambiguously determined (53.9 ± 0.8 °C), along with rather high uncertainty on the high-temperature range of DPPS melting profile obtained from DSC (≈50-59 °C), presumably reflect distinguished surface structural features in LUV. The FTIR signatures of glycerol moiety and those originated from carboxyl group serve as a strong support that in LUV, unlike in MLV, highly curved surfaces occur continuously, whereas the details on the attenuation of surface features in MLV were unraveled by molecular dynamics.

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.

Development of OPLS-AA/M Parameters for Simulations of G Protein-Coupled Receptors and Other Membrane Proteins

G protein-coupled receptors (GPCRs) and other membrane proteins are valuable drug targets, and their dynamic nature makes them attractive systems for study with molecular dynamics (MD) simulations and free energy approaches. Here, we report the development, implementation, and validation of OPLS-AA/M force field parameters to enable simulations of these systems. These efforts include the introduction of post-translational modifications including lipidations and phosphorylation. We also modify previously reported parameters for lipids to be more consistent with the OPLS-AA force field standard and extend their coverage. These new parameters are validated on a variety of test systems, with the results compared to high-level quantum mechanics calculations, experimental data, and simulations with other force fields. The results demonstrate that the new parameters reliably reproduce the behavior of membrane protein systems.

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.

New Insights into the Interaction of Class II Dihydroorotate Dehydrogenases with Ubiquinone in Lipid Bilayers as a Function of Lipid Composition

The fourth enzymatic reaction in the de novo pyrimidine biosynthesis, the oxidation of dihydroorotate to orotate, is catalyzed by dihydroorotate dehydrogenase (DHODH). Enzymes belonging to the DHODH Class II are membrane-bound proteins that use ubiquinones as their electron acceptors. We have designed this study to understand the interaction of an N-terminally truncated human DHODH (HsΔ29DHODH) and the DHODH from Escherichia coli (EcDHODH) with ubiquinone (Q₁₀) in supported lipid membranes using neutron reflectometry (NR). NR has allowed us to determine in situ, under solution conditions, how the enzymes bind to lipid membranes and to unambiguously resolve the location of Q₁₀. Q₁₀ is exclusively located at the center of all of the lipid bilayers investigated, and upon binding, both of the DHODHs penetrate into the hydrophobic region of the outer lipid leaflet towards the Q₁₀. We therefore show that the interaction between the soluble enzymes and the membrane-embedded Q₁₀ is mediated by enzyme penetration. We can also show that EcDHODH binds more efficiently to the surface of simple bilayers consisting of 1-palmitoyl, 2-oleoyl phosphatidylcholine, and tetraoleoyl cardiolipin than HsΔ29DHODH, but does not penetrate into the lipids to the same degree. Our results also highlight the importance of Q₁₀, as well as lipid composition, on enzyme binding.

Lipid21: Complex Lipid Membrane Simulations with AMBER

We extend the modular AMBER lipid force field to include anionic lipids, polyunsaturated fatty acid (PUFA) lipids, and sphingomyelin, allowing the simulation of realistic cell membrane lipid compositions, including raft-like domains. Head group torsion parameters are revised, resulting in improved agreement with NMR order parameters, and hydrocarbon chain parameters are updated, providing a better match with phase transition temperature. Extensive validation runs (0.9 μs per lipid type) show good agreement with experimental measurements. Furthermore, the simulation of raft-like bilayers demonstrates the perturbing effect of increasing PUFA concentrations on cholesterol molecules. The force field derivation is consistent with the AMBER philosophy, meaning it can be easily mixed with protein, small molecule, nucleic acid, and carbohydrate force fields.

Identifying Membrane Lateral Organization by Contrast-Matched Small Angle Neutron Scattering

Lipid domains in model membranes are routinely studied to provide insight into the physical interactions that drive raft formation in cellular membranes. Using small angle neutron scattering, contrast-matching techniques enable the detection of lipid domains ranging from tens to hundreds of nanometers which are not accessible to other techniques without the use of extrinsic probes. Here, we describe a probe-free experimental approach and model-free analysis to identify lipid domains in freely floating vesicles of ternary phase separating lipid mixtures.

Influence of Charge Lipid Head Group Structures on Electric Double Layer Properties

In this study we derived a model for a multicomponent lipid monolayer in contact with an aqueous solution by means of a generalized classical density functional theory and Monte Carlo simulations. Some of the important biological lipid systems were studied as monolayers composed of head groups with different shapes and charge distributions. Starting from the free energy of the system, which includes the electrostatic interactions, additional internal degrees of freedom are included as positional and orientational entropic contributions to the free energy functional. The calculus of variation was used to derive Euler–Lagrange equations, which were solved numerically by the finite element method. The theory and Monte Carlo simulations predict that there are mainly two distinct regions of the electric double layer: (1) the interfacial region, with thickness less than or equal to the length of the fully stretched conformation of the lipid head group, and (2) the outside region, which follows the usual screening of the interface. In the interfacial region, the electric double layer is strongly perturbed, and electrostatic profiles and ion distributions have functionality distinct to classical mean-field theories. Based purely on Coulomb interactions, the theory suggests that the dominant effect on the lipid head group conformation is from the charge density of the interface and the structured lipid mole fraction in the monolayer, rather than the salt concentration in the system.

Vesicle Viewer: Online visualization and analysis of small-angle scattering from lipid vesicles

Small-angle X-ray and neutron scattering are among the most powerful experimental techniques for investigating the structure of biological membranes. Much of the critical information contained in small-angle scattering (SAS) data is not easily accessible to researchers who have limited time to analyze results by hand or to nonexperts who may lack the necessary scientific background to process such data. Easy-to-use data visualization software can allow them to take full advantage of their SAS data and maximize the use of limited resources. To this end, we developed an internet-based application called Vesicle Viewer to visualize and analyze SAS data from unilamellar lipid bilayer vesicles. Vesicle Viewer utilizes a modified scattering density profile (SDP) analysis called EZ-SDP in which key bilayer structural parameters, such as area per lipid and bilayer thickness, are easily and robustly determined. Notably, we introduce a bilayer model that is able to describe an asymmetric bilayer, whether it be chemically or isotopically asymmetric. The application primarily uses Django, a Python package specialized for the development of robust web applications. In addition, several other libraries are used to support the more technical aspects of the project; notable examples are Matplotlib (for graphs) and NumPy (for calculations). By eliminating the barrier of downloading and installing software, this web-based application will allow scientists to analyze their own vesicle scattering data using their preferred operating system. The web-based application can be found at https://vesicleviewer.dmarquardt.ca/.

Lipid Nanodiscs for High-Resolution NMR Studies of Membrane Proteins

Membrane proteins (MPs) play essential roles in numerous cellular processes. Because around 70% of the currently marketed drugs target MPs, a detailed understanding of their structure, binding properties, and functional dynamics in a physiologically relevant environment is crucial for a more detailed understanding of this important protein class. We here summarize the benefits of using lipid nanodiscs for NMR structural investigations and provide a detailed overview of the currently used lipid nanodisc systems as well as their applications in solution-state NMR. Despite the increasing use of other structural methods for the structure determination of MPs in lipid nanodiscs, solution NMR turns out to be a versatile tool to probe a wide range of MP features, ranging from the structure determination of small to medium-sized MPs to probing ligand and partner protein binding as well as functionally relevant dynamical signatures in a lipid nanodisc setting. We will expand on these topics by discussing recent NMR studies with lipid nanodiscs and work out a key workflow for optimizing the nanodisc incorporation of an MP for subsequent NMR investigations. With this, we hope to provide a comprehensive background to enable an informed assessment of the applicability of lipid nanodiscs for NMR studies of a particular MP of interest.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

Structure and Interdigitation of Chain-Asymmetric Phosphatidylcholines and Milk Sphingomyelin in the Fluid Phase

We addressed the frequent occurrence of mixed-chain lipids in biological membranes and their impact on membrane structure by studying several chain-asymmetric phosphatidylcholines and the highly asymmetric milk sphingomyelin. Specifically, we report trans-membrane structures of the corresponding fluid lamellar phases using small-angle X-ray and neutron scattering, which were jointly analyzed in terms of a membrane composition-specific model, including a headgroup hydration shell. Focusing on terminal methyl groups at the bilayer center, we found a linear relation between hydrocarbon chain length mismatch and the methyl-overlap for phosphatidylcholines, and a non-negligible impact of the glycerol backbone-tilting, letting the sn1-chain penetrate deeper into the opposing leaflet by half a CH₂ group. That is, penetration-depth differences due to the ester-linked hydrocarbons at the glycerol backbone, previously reported for gel phase structures, also extend to the more relevant physiological fluid phase, but are significantly reduced. Moreover, milk sphingomyelin was found to follow the same linear relationship suggesting a similar tilt of the sphingosine backbone. Complementarily performed molecular dynamics simulations revealed that there is always a part of the lipid tails bending back, even if there is a high interdigitation with the opposing chains. The extent of this back-bending was similar to that in chain symmetric bilayers. For both cases of adaptation to chain length mismatch, chain-asymmetry has a large impact on hydrocarbon chain ordering, inducing disorder in the longer of the two hydrocarbons.

pH dependent electrical properties of the inner- and outer- leaflets of biomimetic cell membranes

Composition and asymmetry of lipid membranes provide a means for regulation of trans-membrane permeability of ions and small molecules. The pH dependence of these processes plays an important role in the functioning and survival of cells. In this work, we study the pH dependence of membrane electrical resistance and capacitance using electrochemical impedance spectroscopy (EIS), surface plasmon resonance (SPR) and neutron reflectometry (NR) measurements of biomimetic tethered bilayer lipid membranes (tBLMs). tBLMs were prepared with single-component phospholipid compositions, as well as mixtures of phospholipids (phosphatidylcholine, phosphatidylserine, phosphatidylethanolamine, sphingomyelin and cholesterol) that mimic the inner- and outer- leaflets of plasma cell membranes. We found that all studied tBLMs have a resistance maximum at pHs near the pK_as of the phospholipids. SPR and NR indicated that surface concentration of phospholipids and the thickness of the hydrophobic part of the membrane did not change versus pH. We postulate that these maxima are the result of protonation of the phosphate oxygen of the phospholipids and that hydronium ions play a major role in the conductance at pHs < pK_as while sodium ions play the major role at pHs > pK_as. An additional sharp resistance maximum of the PE tBLMs found at pH 5.9 and most likely represents the phosphatidylethanolamine's isoelectric point. The data show the key roles of the characteristic parts of phospholipid molecules: terminal group (choline, carboxyl, amine), phosphate, glycerol and ester oxygens on the permeability and selectivity of ions through the membrane. The interactions between these groups lead to significant differences in the electrical properties of biomimetic models of inner- and outer- leaflets of the plasma cell membranes.

Development of a blocker of the universal phosphatidylserine- and phosphatidylethanolamine-dependent viral entry pathways

Envelope phosphatidylserine (PtdSer) and phosphatidylethanolamine (PtdEtr) have been shown to mediate binding of enveloped viruses. However, commonly used PtdSer binding molecules such as Annexin V cannot block PtdSer-mediated viral infection. Lack of reagents that can conceal envelope PtdSer and PtdEtr and subsequently inhibit infection hinders elucidation of the roles of the envelope phospholipids in viral infection. Here, we developed sTIM1dMLDR801, a reagent capable of blocking PtdSer- and PtdEtr-dependent infection of enveloped viruses. Using sTIM1dMLDR801, we found that envelope PtdSer and/or PtdEtr can support ZIKV infection of not only human but also mosquito cells. In a mouse model for ZIKV infection, sTIM1dMLDR801 reduced ZIKV load in serum and the spleen, indicating envelope PtdSer and/or PtdEtr support in viral infection in vivo. sTIM1dMLDR801 will enable elucidation of the roles of envelope PtdSer and PtdEtr in infection of various virus species, thereby facilitating identification of their receptors and transmission mechanisms.

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.

Expanding MPEx Hydropathy Analysis to Account for Electrostatic Contributions to Protein Interactions with Anionic Membranes

Hydropathy plots are a crucial tool to guide experimental design, as they generate predictions of protein-membrane interactions and their bilayer topology. The predictions are based on experimentally determined hydrophobicity scales, which provide an estimate for the propensity and stability of these interactions. A significant improvement to the accuracy of hydropathy analyses was provided by the development of the popular Wimley-White interfacial and octanol hydrophobicity scales. These scales have been previously incorporated into the freely available MPEx (Membrane Protein Explorer) online application. Here, we introduce a substantial update to MPEx that allows for the consideration of electrostatic contributions to the bilayer partitioning free energy. This component originates from the Coulombic attraction or repulsion of charges between proteins and membranes. Its inclusion in hydropathy calculations increases the accuracy of hydropathy plot predictions and extends their use to more complex systems (i.e., anionic membranes). We illustrate the application of this analysis to studies on the membrane selectivity of antimicrobial peptides, the membrane partitioning of ion-channel gating modifiers, and the amyloid proteins α-synuclein and Tau, as well as pH-dependent bilayer interactions of diphtheria toxin and apoptotic inhibitor Bcl-xL.

Understanding the Structure and Dynamics of Complex Biomembrane Interactions by Neutron Scattering Techniques

The membrane is one of the key structural materials of biology at the cellular level. Composed predominantly of a bilayer of lipids with embedded and bound proteins, it defines the boundaries of the cell and many organelles essential to life and therefore is involved in almost all biological processes. Membrane-specific interactions, such as drug binding to a membrane receptor or the interactions of an antimicrobial compound with the lipid matrix of a pathogen membrane, are of interest across the scientific disciplines. Herein we present a review, aimed at nonexperts, of the major neutron scattering techniques used in membrane studies: small-angle neutron scattering, neutron membrane diffraction, neutron reflectometry, quasielastic neutron scattering, and neutron spin echo. Neutron scattering techniques are well suited to studying biological membranes. The nondestructive nature of cold neutrons means that samples can be measured for long periods without fear of beam damage from ultraviolet, electron, or X-ray radiation, and neutron beams are highly penetrating, thus offering flexibility in samples and sample environments. Most important is the strong difference in neutron scattering lengths between the two most abundant forms of hydrogen, protium and deuterium. Changing the relative amounts of protium/deuterium in a sample allows the production of a series of neutron scattering data sets, enabling the observation of differing components within complex membrane architectures. This approach can be as simple as using the naturally occurring neutron contrast between different biomolecules to study components in a complex by changing the solution H₂O/D₂O ratio or as complex as selectively labeling individual components with hydrogen isotopes. This review presents an overview of each experimental technique with the neutron instrument configuration, related sample preparation and sample environment, and data analysis, highlighted by a special emphasis on using prominent neutron contrast to understand structure and dynamics. This review gives researchers a practical introduction to the often enigmatic suite of neutron beamlines, thereby lowering the barrier to taking advantage of these large-facility techniques to achieve new understandings of membranes and their interactions with other molecules.

Structural Characterization of Biomaterials by Means of Small Angle X-rays and Neutron Scattering (SAXS and SANS), and Light Scattering Experiments

Scattering techniques represent non-invasive experimental approaches and powerful tools for the investigation of structure and conformation of biomaterial systems in a wide range of distances, ranging from the nanometric to micrometric scale. More specifically, small-angle X-rays and neutron scattering and light scattering techniques represent well-established experimental techniques for the investigation of the structural properties of biomaterials and, through the use of suitable models, they allow to study and mimic various biological systems under physiologically relevant conditions. They provide the ensemble averaged (and then statistically relevant) information under in situ and operando conditions, and represent useful tools complementary to the various traditional imaging techniques that, on the contrary, reveal more local structural information. Together with the classical structure characterization approaches, we introduce the basic concepts that make it possible to examine inter-particles interactions, and to study the growth processes and conformational changes in nanostructures, which have become increasingly relevant for an accurate understanding and prediction of various mechanisms in the fields of biotechnology and nanotechnology. The upgrade of the various scattering techniques, such as the contrast variation or time resolved experiments, offers unique opportunities to study the nano- and mesoscopic structure and their evolution with time in a way not accessible by other techniques. For this reason, highly performant instruments are installed at most of the facility research centers worldwide. These new insights allow to largely ameliorate the control of (chemico-physical and biologic) processes of complex (bio-)materials at the molecular length scales, and open a full potential for the development and engineering of a variety of nano-scale biomaterials for advanced applications.

Lipid headgroups alter huntingtin aggregation on membranes

Huntington's Disease is a fatal neurodegenerative disorder caused by expansion of a glutamine repeat region (polyQ) beyond a critical threshold within exon1 of the huntingtin protein (htt). As a consequence of polyQ expansion, htt associates into a variety of aggregate species that are thought to underlie cellular toxicity. Within cells, htt associates with numerous membranous organelles and surfaces that exert influence on the aggregation process. In particular, the first 17 amino acids at the N-terminus of htt (Nt17) serve as a lipid-binding domain that is intrinsically disordered in bulk solution but adopts an amphipathic α-helical structure upon binding membranes. Beyond this, Nt17 is implicated in initiating htt fibrillization. As the interaction between Nt17 and lipid membranes is likely influenced by lipid properties, the impact of lipid headgroups on htt-exon1 aggregation, membrane activity, and the ability to form protein:lipid complexes was determined. Htt-exon1 with a disease-length polyQ domain (46Q) was exposed to lipid vesicles comprised of lipids with either zwitterionic (POPC and POPE) or anionic (POPG and POPS) headgroups. With zwitterionic head groups, large lipid to peptide ratios were required to have a statistically significant impact on htt aggregation. Anionic lipids enhanced htt fibrillization, even at low lipid:protein ratios, and this was accompanied by changes in aggregate morphology. Despite the larger impact of anionic lipids, htt-exon1(46Q) was more membrane active with zwitterionic lipid systems. The ability of Nt17 to form complexes with lipids was also mediated by lipid headgroups as zwitterionic ionic lipids more readily associated with multimeric forms of Nt17 in comparison with anionic lipids. Collectively, these results highlight the complexity of htt/membrane interactions and the resulting impact on the aggregation process.

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

We have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Ų, a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Ų). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S_CD) profiles calculated from the simulations were in relatively good agreement with S_CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Ų resulted in good agreement between the simulation and experimental scattering form factors, but not with S_CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

Membrane Interactions of Cy3 and Cy5 Fluorophores and Their Effects on Membrane-Protein Dynamics

Organic fluorophores, such as Cy3 and Cy5, have been widely used as chemical labels to probe the structure and dynamics of membrane proteins. Although a number of previous studies have reported on the possibility of some of the water-soluble fluorophores to interact with lipid bilayers, detailed fluorophore-lipid interactions and, more importantly, the potential effect of such interactions on the natural dynamics of the labeled membrane proteins have not been well studied. We have performed a large set of all-atom molecular dynamics simulations employing the highly mobile membrane mimetic model to describe spontaneous partitioning of the fluorophores into lipid bilayers with different lipid compositions. Spontaneous membrane partitioning of Cy3 and Cy5 fluorophores captured in these simulations proceeds in two steps. Electrostatic interaction between the fluorophores and the lipid headgroups facilitates the initial, fast membrane association of the fluorophores, followed by slow insertion of hydrophobic moieties into the lipid bilayer core. After the conversion of the resulting membrane-bound systems to full-membrane representations, biased-exchange umbrella sampling simulations are performed for free energy calculations, revealing a higher energy barrier for partitioning into negatively charged (phosphatidylserine or phosphatidylcholine) membranes than purely zwitterionic (phosphatidylcholine or phosphatidylethanolamine) ones. Furthermore, the potential effect of fluorophore-lipid interactions on membrane proteins has been examined by covalently linking Cy5 to single- and multipass transmembrane helical proteins. Equilibrium simulations show strong position-dependent effects of Cy5-tagging on the structure and natural dynamics of membrane proteins. Interactions between the tagged protein and Cy5 were also observed. Our results suggest that fluorophore-lipid interactions can affect the structure and dynamics of membrane proteins to various extents, especially in systems with higher structural flexibility.

The structures of polyunsaturated lipid bilayers by joint refinement of neutron and X-ray scattering data

We present the detailed structural analysis of polyunsaturated fatty acid-containing phospholipids namely, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC). A newly developed molecular dynamics (MD) simulation parsing scheme for lipids containing fatty acids with multiple double bonds was implemented into the scattering density profile (SDP) model to simultaneously refine differently contrasted neutron and X-ray scattering data. SDP analyses of scattering data at 30 °C yielded lipid areas of 71.1 Ų and 70.4 Ų for PDPC and SDPC bilayers, respectively, and a model free analysis of PDPC at 30 °C resulted in a lipid area of 72 Ų. In addition to bilayer structural parameters, using area-constrained MD simulations we determined the area compressibility modulus, K_A, to be 246.4 mN/m, a value similar to other neutral phospholipids.

The Fats of Life: Using Computational Chemistry to Characterise the Eukaryotic Cell Membrane

Our current knowledge of the structural dynamics and complexity of lipid bilayers is still developing. Computational techniques, especially molecular dynamics simulations, have increased our understanding significantly as they allow us to model functions that cannot currently be experimentally resolved. Here we review available computational tools and techniques, the role of the major lipid species, insights gained into lipid bilayer structure and function from molecular dynamics simulations, and recent progress towards the computational modelling of the physiological complexity of eukaryotic lipid bilayers.

Improved Cation Binding to Lipid Bilayers with Negatively Charged POPS by Effective Inclusion of Electronic Polarization

Phosphatidylserine (PS) lipids are important signaling molecules and the most common negatively charged lipids in eukaryotic membranes. The signaling can be often regulated by calcium, but its interactions with PS headgroups are not fully understood. Classical molecular dynamics (MD) simulations can potentially give detailed description of lipid-ion interactions, but the results strongly depend on the used force field. Here, we apply the electronic continuum correction (ECC) to the Amber Lipid17 parameters of 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipid to improve its interactions with K⁺, Na⁺, and Ca²⁺ ions. The partial charges of the headgroup, glycerol backbone, and carbonyls of POPS, bearing a unit negative charge, were scaled with a factor of 0.75, derived for monovalent ions, and the Lennard-Jones σ parameters of the same segments were scaled with a factor of 0.89. The resulting ECC-POPS model gives more realistic interactions with Na⁺ and Ca²⁺ cations than the original Amber Lipid17 parameters when validated using headgroup order parameters and the "electrometer concept". In ECC-lipids simulations, populations of complexes of Ca²⁺ cations with more than two PS lipids are negligible, and interactions of Ca²⁺ cations with only carboxylate groups are twice more likely than with only phosphate groups, while interactions with carbonyls almost entirely involve other groups as well. Our results pave the way for more realistic MD simulations of biomolecular systems with anionic membranes, allowing signaling processes involving PS and Ca²⁺ to be elucidated.

On the Mechanism of Bilayer Separation by Extrusion, or Why Your LUVs Are Not Really Unilamellar

Extrusion through porous filters is a widely used method for preparing biomimetic model membranes. Of primary importance in this approach is the efficient production of single bilayer (unilamellar) vesicles that eliminate the influence of interlamellar interactions and strictly define the bilayer surface area available to external reagents such as proteins. Submicroscopic vesicles produced using extrusion are widely assumed to be unilamellar, and large deviations from this assumption would impact interpretations from many model membrane experiments. Using three probe-free methods-small angle X-ray and neutron scattering and cryogenic electron microscopy-we report unambiguous evidence of extensive multilamellarity in extruded vesicles composed of neutral phosphatidylcholine lipids, including for the common case of neutral lipids dispersed in physiological buffer and extruded through 100-nm diameter pores. In such preparations, only ∼35% of lipids are externally accessible and this fraction is highly dependent on preparation conditions. Charged lipids promote unilamellarity as does decreasing solvent ionic strength, indicating the importance of electrostatic interactions in determining the lamellarity of extruded vesicles. Smaller extrusion pore sizes also robustly increase the fraction of unilamellar vesicles, suggesting a role for membrane bending. Taken together, these observations suggest a mechanistic model for extrusion, wherein the formation of unilamellar vesicles involves competition between bilayer bending and adhesion energies. The findings presented here have wide-ranging implications for the design and interpretation of model membrane studies, especially ensemble-averaged observations relying on the assumption of unilamellarity.

Intrinsic Curvature-Mediated Transbilayer Coupling in Asymmetric Lipid Vesicles

We measured the effect of intrinsic lipid curvature, J0, on structural properties of asymmetric vesicles made of palmitoyl-oleoyl-phosphatidylethanolamine (POPE; J0<0) and palmitoyl-oleoyl-phosphatidylcholine (POPC; J0∼0). Electron microscopy and dynamic light scattering were used to determine vesicle size and morphology, and x-ray and neutron scattering, combined with calorimetric experiments and solution NMR, yielded insights into leaflet-specific lipid packing and melting processes. Below the lipid melting temperature we observed strong interleaflet coupling in asymmetric vesicles with POPE inner bilayer leaflets and outer leaflets enriched in POPC. This lipid arrangement manifested itself by lipids melting cooperatively in both leaflets, and a rearrangement of lipid packing in both monolayers. On the other hand, no coupling was observed in vesicles with POPC inner bilayer leaflets and outer leaflets enriched in POPE. In this case, the leaflets melted independently and did not affect each other’s acyl chain packing. Furthermore, we found no evidence for transbilayer structural coupling above the melting temperature of either sample preparation. Our results are consistent with the energetically preferred location of POPE residing in the inner leaflet, where it also resides in natural membranes, most likely causing the coupling of both leaflets. The loss of this coupling in the fluid bilayers is most likely the result of entropic contributions.

Lipid-modulation of membrane insertion and refolding of the apoptotic inhibitor Bcl-xL

Bcl-xL is a member of the Bcl-2 family of apoptotic regulators, responsible for inhibiting the permeabilization of the mitochondrial outer membrane, and a promising anti-cancer target. Bcl-xL exists in the following conformations, each believed to play a role in the inhibition of apoptosis: (a) a soluble folded conformation, (b) a membrane-anchored (by its C-terminal α8 helix) form, which retains the same fold as in solution and (c) refolded membrane-inserted conformations, for which no structural data are available. Previous studies established that in the cell Bcl-xL exists in a dynamic equilibrium between soluble and membranous states, however, no direct evidence exists in support of either anchored or inserted conformation of the membranous state in vivo. In this in vitro study, we employed a combination of fluorescence and EPR spectroscopy to characterize structural features of the bilayer-inserted conformation of Bcl-xL and the lipid modulation of its membrane insertion transition. Our results indicate that the core hydrophobic helix α6 inserts into the bilayer without adopting a transmembrane orientation. This insertion disrupts the packing of Bcl-xL and releases the regulatory N-terminal BH4 domain (α1) from the rest of the protein structure. Our data demonstrate that both insertion and refolding of Bcl-xL are modulated by lipid composition, which brings the apparent pK_a of insertion to the threshold of physiological pH. We hypothesize that conformational rearrangements associated with the bilayer insertion of Bcl-xL result in its switching to a so-called non-canonical mode of apoptotic inhibition. Presented results suggest that the alteration in lipid composition before and during apoptosis can serve as an additional factor regulating the permeabilization of the mitochondrial outer membrane.

Phosphatidylserine Asymmetry Promotes the Membrane Insertion of a Transmembrane Helix

The plasma membrane (PM) contains an asymmetric distribution of lipids between the inner and outer bilayer leaflets. A lipid of special interest in eukaryotic membranes is the negatively charged phosphatidylserine (PS). In healthy cells, PS is actively sequestered to the inner leaflet of the PM, but PS redistributes to the outer leaflet when the cell is damaged or at the onset of apoptosis. However, the influence of PS asymmetry on membrane protein structure and folding are poorly understood. The pH low insertion peptide (pHLIP) adsorbs to the membrane surface at a neutral pH, but it inserts into the membrane at an acidic pH. We have previously observed that in symmetric vesicles, PS affects the membrane insertion of pHLIP by lowering the pH midpoint of insertion. Here, we studied the effect of PS asymmetry on the membrane interaction of pHLIP. We developed a modified protocol to create asymmetric vesicles containing PS and employed Annexin V labeled with an Alexa Fluor 568 fluorophore as a new probe to quantify PS asymmetry. We observed that the membrane insertion of pHLIP was promoted by the asymmetric distribution of negatively charged PS, which causes a surface charge difference between bilayer leaflets. Our results indicate that lipid asymmetry can modulate the formation of an α-helix on the membrane. A corollary is that model studies using symmetric bilayers to mimic the PM may fail to capture important aspects of protein-membrane interactions.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Pathways of amyloid-beta absorption and aggregation in a membranous environment

Aggregation of misfolded oligomeric amyloid-beta (Aβ) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aβ aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aβ 16-22 (K-L-V-F-F-A-E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge - zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine: POPS) - on Aβ 16-22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aβ peptides that contribute towards faster emergence of ordered β-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide-lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aβ aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process.

Calcium-Induced Lipid Nanocluster Structures: Sculpturing of the Plasma Membrane

The plasma membrane of the cell is a complex, tightly regulated, heterogeneous environment shaped by proteins, lipids, and small molecules. Ca²⁺ ions are important cellular messengers, spatially separated from anionic lipids. After cell injury, disease, or apoptotic events, anionic lipids are externalized to the outer leaflet of the plasma membrane and encounter Ca²⁺, resulting in dramatic changes in the plasma membrane structure and initiation of signaling cascades. Despite the high chemical and biological significance, the structures of lipid-Ca²⁺ nanoclusters are still not known. Previously, we demonstrated by solid-state nuclear magnetic resonance (NMR) spectroscopy that upon binding to Ca²⁺, individual phosphatidylserine lipids populate two distinct yet-to-be-characterized structural environments. Here, we concurrently employ extensive all-atom molecular dynamics (MD) simulations with our accelerated membrane mimetic and detailed NMR measurements to identify lipid-Ca²⁺ nanocluster conformations. We find that major structural characteristics of these nanoclusters, including interlipid pair distances and chemical shifts, agree with observable NMR parameters. Simulations reveal that lipid-ion nanoclusters are shaped by two characteristic, long-lived lipid structures induced by divalent Ca²⁺. Using ab initio quantum mechanical calculations of chemical shifts on MD-captured lipid-ion complexes, we show that computationally observed conformations are validated by experimental NMR data. Both NMR measurements of diluted specifically labeled lipids and MD simulations reveal that the basic structural unit that reshapes the membrane is a Ca²⁺-coordinated phosphatidylserine tetramer. Our combined computational and experimental approach presented here can be applied to other complex systems in which charged membrane-active molecular agents leave structural signatures on lipids.

Effects of phosphatidylserine-containing supported lipid bilayers on the polarization of macrophages

Placement of dental implants initiates inflammatory foreign body response, in which macrophages play a central role and affect the subsequent tissue healing process such as bone formation. The purpose of this study was to fabricate phosphatidylserine (PS)-containing supported lipid bilayers (SLBs) on a titanium surface to regulate the polarization of macrophages, a critical factor that affects following tissue healing and regeneration. The fluorescent recovery after photobleaching images showed that the percentage of PS had a critical influence on the fluidity, and 20% PS had the highest fluidity. Furthermore, more expanded and elongated cells were observed in the SLB-coated groups. transforming growth factor-β1 and vascular endothelial growth factor, the key cytokine markers of M2 macrophage polarization, were increased in the SLB-coated groups, especially in the 20% PS group. Consistently, cells cultured on the SLB-coated titanium exhibited the distribution of CD206⁺ , which is a M2 macrophage specific maker. The results of this study demonstrated M2 polarization of macrophages on PS-SLB-coated titanium discs, which suggests the application of PS-SLB as an immune-regulating coating material to improve tissue reactions to dental implants. © 2018 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 106A: 2625-2633, 2018.

The impact of deuteration on natural and synthetic lipids: A neutron diffraction study

The structural investigation of cellular membranes requires access to model systems where the molecular complexity is representative of the cellular environment and that allow for the exploitation of structural techniques. Neutron scattering, and in particular neutron diffraction can provide unique and detailed information on the structure of lipid membranes. However, deuterated samples are desirable to fully exploit this powerful method. Recently, the extraction of lipids from microorganisms grown in deuterated media was demonstrated to be both an attracting route to obtain complex lipid mixtures resembling the composition of natural membranes, and to producing deuterated molecules in a very convenient way. A full characterization of these deuterated extracts is hence pivotal for their use in building up model membrane systems. Here we report the structural characterization of lipid extracts obtained from Pichia pastoris by means of neutron diffraction measurements. In particular, we compare the structure of membranes extracted from yeast cells grown in a standard culture medium and in a corresponding deuterated culture medium. The results show that the different molecular composition of the deuterated and protiated lipid extracts induce different structural organization of the lipid membranes. In addition, we compare these membranes composed of extracted yeast lipids with stacked bilayers prepared from synthetic lipid mixtures.

Role of Phosphorylation in Moesin Interactions with PIP2-Containing Biomimetic Membranes

Moesin, a protein of the ezrin, radixin, and moesin family, which links the plasma membrane to the cytoskeleton, is involved in multiple physiological and pathological processes, including viral budding and infection. Its interaction with the plasma membrane occurs via a key phosphoinositide, the phosphatidyl(4,5)inositol-bisphosphate (PIP2), and phosphorylation of residue T558, which has been shown to contribute, in cellulo, to a conformationally open protein. We study the impact of a double phosphomimetic mutation of moesin (T235D, T558D), which mimics the phosphorylation state of the protein, on protein/PIP2/microtubule interactions. Analytical ultracentrifugation in the micromolar range showed moesin in the monomer and dimer forms, with wild-type (WT) moesin containing a slightly larger fraction (∼30%) of dimers than DD moesin (10-20%). Only DD moesin was responsive to PIP2 in its micellar form. Quantitative cosedimentation assays using large unilamellar vesicles and quartz crystal microbalance on supported lipid bilayers containing PIP2 reveal a specific cooperative interaction for DD moesin with an ability to bind two PIP2 molecules simultaneously, whereas WT moesin was able to bind only one. In addition, DD moesin could subsequently interact with microtubules, whereas WT moesin was unable to do so. Altogether, our results point to an important role of these two phosphorylation sites in the opening of moesin: since DD moesin is intrinsically in a more open conformation than WT moesin, this intermolecular interaction is reinforced by its binding to PIP2. We also highlight important differences between moesin and ezrin, which appear to be finely regulated and to exhibit distinct molecular behaviors.

Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties

Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.

Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties

Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.

High-Affinity Interactions of Beryllium(2+) with Phosphatidylserine Result in a Cross-Linking Effect Reducing Surface Recognition of the Lipid

Beryllium has multiple industrial applications, but its manufacture is associated with a serious occupational risk of developing chronic inflammation in the lungs known as berylliosis, or chronic beryllium disease. Although the Be²⁺-induced abnormal immune responses have recently been linked to a specific MHC-II allele, the nature of long-lasting granulomas is not fully understood. Here we show that Be²⁺ binds with a micromolar affinity to phosphatidylserine (PS), the major surface marker of apoptotic cells. Isothermal titration calorimetry indicates that, like that of Ca²⁺, binding of Be²⁺ to PS liposomes is largely entropically driven, likely by massive desolvation. Be²⁺ exerts a compacting effect on PS monolayers, suggesting cross-linking through coordination by both phosphates and carboxyls in multiple configurations, which were visualized in molecular dynamics simulations. Electrostatic modification of PS membranes by Be²⁺ includes complete neutralization of surface charges at ∼30 μM, accompanied by an increase in the boundary dipole potential. The data suggest that Be²⁺ can displace Ca²⁺ from the surface of PS, and being coordinated in a tight shell of four oxygens, it can mask headgroups from Ca²⁺-mediated recognition by PS receptors. Indeed, 48 μM Be²⁺ added to IC-21 cultured macrophages specifically suppresses binding and engulfment of PS-coated silica beads or aged erythrocytes. We propose that Be²⁺ adsorption at the surface of apoptotic cells may potentially prevent normal phagocytosis, thus causing accumulation of secondary necrotic foci and the resulting chronic inflammation.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.