Cholesterol Flip-Flop in Heterogeneous Membranes

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

Time-resolved small-angle neutron scattering for characterization of molecular exchange in lipid nanoparticle therapeutics

HYPOTHESIS

Nano-scale dynamics of self-assembled therapeutics play a large role in their biological function. However, assessment of such dynamics remains absent from conventional pharmaceutical characterization. We hypothesize that time-resolved small-angle neutron scattering (TR-SANS) can reveal their kinetic properties. For lipid nanoparticles (LNP), limited molecular motion is important for avoiding degradation prior to entering cells while, intracellularly, enhanced molecular motion is then vital for effective endosomal escape. We propose TR-SANS for quantifying molecular exchange in LNPs and, therefore, enabling optimization of opposing molecular behaviors of a pharmaceutical in two distinct environments.

EXPERIMENTS

We use TR-SANS in combination with traditional SANS and small-angle x-ray scattering (SAXS) to experimentally quantify nano-scale dynamics and provided unprecedented insight to molecular behavior of LNPs.

FINDINGS

LNPs have molecular exchange dynamics relevant to storage and delivery which can be captured using TR-SANS. Cholesterol exchanges on the time-scale of hours even at neutral pH. As pH drops below the effective pKa of the ionizable lipid, molecular exchange occurs faster. The results give insight into behavior enabling delivery and provide a quantifiable metric by which to compare formulations. Successful analysis of this multi-component system also expands the opportunities for using TR-SANS to characterize complex therapeutics.

Measuring the Transmembrane Registration of Lipid Domains in Droplet Interface Bilayers through Tensiometry

Diverse collections of lipids self-assemble into domains within biological membranes, and these domains are typically organized in both the transverse and lateral directions of the membrane. The ability of the membrane to link these domains across the membrane's interior grants cells control over features on the external cellular surface. Numerous hypothesized factors drive the cross-membrane (or transverse) coupling of lipid domains. In this work we seek to isolate these transverse lipid-lipid influences in a simple model system using droplet interface bilayers (DIBs) to better understand the associated mechanics. DIBs enable symmetric and asymmetric combinations of domain-forming lipid mixtures within a model bilayer, and the evolving energetics of the membrane may be tracked using drop-shape analysis. We find that symmetric distributions of domain-forming lipids produce long-lasting, gradual shifts in the DIB membrane energetics that are not observed in asymmetric distributions of the lipids where the domain-forming lipids are only within one leaflet. The approach selected for this work provides experimental measurement of the mismatch penalty associated with antiregistered lipid domains as well as measurements of the influence of rafts on DIB behaviors with suggestions for their future use as a model platform.

Strategies for Improved pDNA Loading and Protection Using Cationic and Neutral LNPs with Industrial Scalability Potential Using Microfluidic Technology

Purpose

In recent years, microfluidic technologies have become mainstream in producing gene therapy nanomedicines (NMeds) following the Covid-19 vaccine; however, extensive optimizations are needed for each NMed type and genetic material. This article strives to improve LNPs for pDNA loading, protection, and delivery, while minimizing toxicity.

Methods

The microfluidic technique was optimized to form cationic or neutral LNPs to load pDNA. Classical "post-formulation" DNA addition vs "pre" addition in the aqueous phase were compared. All formulations were characterized (size, homogeneity, zeta potential, morphology, weight yield, and stability), then tested for loading efficiency, nuclease protection, toxicity, and cell uptake.

Results

Optimized LNPs formulated with DPPC: Chol:DOTAP 1:1:0.1 molar ratio and 10 µg of DOPE-Rhod, had a size of 160 nm and good homogeneity. The chemico-physical characteristics of cationic LNPs worsened when adding 15 µg/mL of pDNA with the "post" method, while maintaining their characteristics up to 100 µg/mL of pDNA with the "pre" addition remaining stable for 30 days. Interestingly, neutral LNPs formulated with the same method loaded up to 50% of the DNA. Both particles could protect the DNA from nucleases even after one month of storage, and low cell toxicity was found up to 40 µg/mL LNPs. Cell uptake occurred within 2 hours for both formulations with the DNA intact in the cytoplasm, outside of the lysosomes.

Conclusion

In this study, the upcoming microfluidic technique was applied to two strategies to generate pDNA-LNPs. Cationic LNPs could load 10x the amount of DNA as the classical approach, while neutral LNPs, which also loaded and protected DNA, showed lower toxicity and good DNA protection. This is a big step forward at minimizing doses and toxicity of LNP-based gene therapy.

Simulating asymmetric membranes using P21 periodic boundary conditions

Molecular dynamics (MD) simulations of symmetric lipid bilayers are now well established, while those of asymmetric ones are considerably less developed. This disjunction arises in part because the surface tensions of leaflets in asymmetric bilayers can differ (unlike those of symmetric ones), and there is no simple way to determine them without assumptions. This chapter describes the use of P21 periodic boundary conditions (PBC), which allow lipids to switch leaflets, to generate asymmetric bilayers under the assumption of equal chemical potentials of lipids in opposing leaflets. A series of examples, ranging from bilayers with one lipid type to those with peptides and proteins, provides a guide for the use of P21 PBC. Critical properties of asymmetric membranes, such as spontaneous curvature, are highly sensitive to differences in the leaflet surface tensions (or differential stress), and equilibration with P21 PBC substantially reduces differential stress of asymmetric bilayers assembled with surface area-based methods. Limitations of the method are discussed. Technically, the nonstandard unit cell is difficult to parallelize and to incorporate restraints. Inherently, the assumption of equal chemical potentials, and therefore the method itself, is not applicable to all target systems. Despite these limitations, it is argued that P21 simulations should be considered when designing equilibration protocols for MD studies of most asymmetric membranes.

Exploring membrane asymmetry and its effects on membrane proteins

Plasma membranes utilize free energy to maintain highly asymmetric, non-equilibrium distributions of lipids and proteins between their two leaflets. In this review we discuss recent progress in quantitative research enabled by using compositionally controlled asymmetric model membranes. Both experimental and computational studies have shed light on the nuanced mechanisms that govern the structural and dynamic coupling between compositionally distinct bilayer leaflets. This coupling can increase the membrane bending rigidity and induce order - or lipid domains - across the membrane. Furthermore, emerging evidence indicates that integral membrane proteins not only respond to asymmetric lipid distributions but also exhibit intriguing asymmetric properties themselves. We propose strategies to advance experimental research, aiming for a deeper, quantitative understanding of membrane asymmetry, which carries profound implications for cellular physiology.

Studying lipid flip-flop in asymmetric liposomes using 1H NMR and TR-SANS

The specific spatial and temporal distribution of lipids in membranes play a crucial role in determining the biochemical and biophysical properties of the system. In nature, the asymmetric distribution of lipids is a dynamic process with ATP-dependent lipid transporters maintaining asymmetry, and passive transbilayer diffusion, that is, flip-flop, counteracting it. In this chapter, two probe-free techniques, 1H NMR and time-resolved small angle neutron scattering, are described in detail as methods of investigating lipid flip-flop rates in synthetic liposomes that have been generated with an asymmetric bilayer composition.

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.

Molecular mechanism of calcitriol enhances membrane water permeability

Helicobacter pylori (H. pylori) exhibits a unique membrane lipid composition, including dimyristoyl phosphatidylethanolamine (DMPE) and cholesterol, unlike other Gram-negative bacteria. Calcitriol has antimicrobial activity against H. pylori, but cholesterol enhances antibiotics resistance in H. pylori. This study explored the changes in membrane structure and the molecular mechanisms of cholesterol/calcitriol translocation using well-tempered metadynamics (WT-MetaD) simulations and microsecond conventional molecular dynamics (CMD) simulations. Calcitriol facilitated water transport across the membrane, while cholesterol had the opposite effect. The differing effects might result from the tail 25-hydroxyl group and a wider range of orientations of calcitriol in the DMPE/dimyristoyl phosphatidylglycerol (DMPG) (3:1) membrane. Calcitriol moves across the bilayer center without changing its orientation along the membrane Z-axis, becomes parallel to the membrane surface at the membrane-water interface, and then rotates approximately 90° in this interface. The translocation mechanism of calcitriol is quite different from the flip-flop of cholesterol. Moreover, calcitriol crossed from one layer to another more easily than cholesterol, causing successive perturbations to the hydrophobic core and increasing water permeation. These results improve our understanding of the relationship between cholesterol/calcitriol concentrations and the lipid bilayer structure and the role of lipid composition in water permeation.

The Structure of Oxysterols Determines Their Behavior at Phase Boundaries: Implications for Model Membranes and Structure-Activity Relationships

The presence of an additional polar group in the cholesterol backbone increases the hydrophilicity of resulting compounds (oxysterols), determines their arrangement at the phase boundary, and interactions with other lipids and proteins. As a result, physicochemical properties of biomembranes (i.e., elasticity, permeability, and ability to bind proteins) are modified, which in turn may affect their functioning. The observed effect depends on the type of oxysterol and its concentration and can be both positive (e.g., antiviral activity) or negative (disturbance of cholesterol homeostasis, signal transduction, and protein segregation). The membrane activity of oxysterols has been successfully studied using membrane models (vesicles, monolayers, and solid supported films). Membrane models, in contrast to the natural systems, provide the possibility to selectively examine the specific aspect of biomolecule-membrane interactions. Moreover, the gradual increase in the complexity of the used model allows to understand the molecular phenomena occurring at the membrane level. The interest in research on artificial membranes has increased significantly in recent years, mainly due to the development of modern and sophisticated physicochemical methods (static and dynamic) in both the micro- and nanoscale, which are applied with the assistance of powerful theoretical calculations. This review provides an overview of the most important findings on this topic in the current literature.

Understanding the free-energy landscape of phase separation in lipid bilayers using molecular dynamics

Liquid-liquid phase separation inside the cell often results in biological condensates that can critically affect cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For these reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation-essentially measuring the sign (but not the amplitude) of the free-energy change upon separation-rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from MD simulations.

Understanding the Free Energy Landscape of Phase Separation in Lipid Bilayers using Molecular Dynamics

Liquid-liquid phase separation (LLPS) inside the cell often results in biological condensates that can critically impact cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the so-called "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For this reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation - essentially measuring the sign (but not the amplitude) of the free energy change upon separation - rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from molecular dynamics simulations.

SIGNIFICANCE

Standard molecular dynamics simulations can determine the sign the free energy change upon phase separation, but not the amplitude. We present a new method to determine the phase separation free energy for lipid membranes, based on a enhanced sampling using the weighted ensemble method combined with a novel collective variable, validated using coarse-grained simulations applied to several simple systems. The new method will be valuable as a way to develop models that connect molecular-level structural features to the thermodynamics of phase separation.

Simulated dynamic cholesterol redistribution favors membrane fusion pore constriction

Endo- and exocytosis proceed through a highly strained membrane fusion pore topology regardless of the aiding protein machinery. The membrane's lipid components bias fusion pores toward expansion or closure, modifying the necessary work done by proteins. Cholesterol, a key component of plasma membranes, promotes both inverted lipid phases with concave leaflets (i.e., negative total curvature, which thins the leaflet) and flat bilayer phases with thick, ordered hydrophobic interiors. We demonstrate by theory and simulation that both leaflets of nascent catenoidal fusion pores have negative total curvature. Furthermore, the hydrophobic core of bilayers with strong negative Gaussian curvature is thinned. Therefore, it is an open question whether cholesterol will be enriched in these regions because of the negative total curvature or depleted because of the membrane thinning. Here, we compare all-atom molecular dynamics simulations (built using a procedure to create specific fusion pore geometries) and theory to understand the underlying reasons for lipid redistribution on fusion pores. Our all-atom molecular dynamics simulations resolve this question by showing that cholesterol is strongly excluded from the thinned neck of fusion and fission pores, revealing that thickness (and/or lipid order) influences cholesterol distributions more than curvature. The results imply that cholesterol exclusion can drive fusion pore closure by creating a small, cholesterol-depleted zone in the neck. This model agrees with literature evidence that membrane reshaping is connected to cholesterol-dependent lateral phase separation.

Preventing the amyloid-beta peptides accumulation on the cell membrane by applying GHz electric fields: A molecular dynamic simulation

Alzheimer's disease is associated with accumulating different amyloid peptides on the nerve cell membranes. The non-thermal effects of the GHz electric fields in this topic have yet to be well recognized. Hence, in this study, the impacts of 1 and 5 GHz electric fields on the amyloid peptide proteins accumulation on the cell membrane have been investigated, utilizing molecular dynamics (MD) simulation. The obtained results indicated that this range of electric fields did not significantly affect the peptide structure. Moreover, it was found that the peptide penetration into the membrane was increased as the field frequency was increased when the system was exposed to a 20 mv/nm oscillating electric field. In addition, it was observed that the protein-membrane interaction is reduced significantly in the presence of the 70 mv/nm electric field. The molecular level results reported in this study could be helpful in better understanding Alzheimer's disease.

Coarse-grained modeling of polystyrene-modified CNTs and their interactions with lipid bilayers

In the present work, we describe Martini3 coarse-grained models of polystyrene and carboxyl-terminated polystyrene functionalized carbon nanotubes (CNTs) and investigate their interactions with lipid bilayers with and without cholesterol (CHOL) using molecular dynamics simulations. By changing the polystyrene chain length and grafting density at the end ring of the CNTs at two different nanotube concentrations, we observe the translocation of nanoparticles as well as changes in the lipid bilayer properties. Our results show that all developed models passively diffuse into the membranes without causing any damage to the membrane integrity, although high concentrations of CNTs induce structural and elastic changes in lipid bilayers. In the presence of CHOL, increasing CNT concentration results in decreased rates of CHOL transmembrane motions. On the other hand, CNTs are prone to lipid and polystyrene blockage, which affects their equilibrated configurations, and tilting behavior within the membranes. Hence, we demonstrate that polystyrene-functionalized CNTs are promising drug-carrier agents. However, polystyrene chain length and grafting density are important factors to consider to enhance the efficiency of drug delivery.

Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-Hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.

IFITM3 blocks influenza virus entry by sorting lipids and stabilizing hemifusion

Interferon-induced transmembrane protein 3 (IFITM3) inhibits the entry of numerous viruses through undefined molecular mechanisms. IFITM3 localizes in the endosomal-lysosomal system and specifically affects virus fusion with target cell membranes. We found that IFITM3 induces local lipid sorting, resulting in an increased concentration of lipids disfavoring viral fusion at the hemifusion site. This increases the energy barrier for fusion pore formation and the hemifusion dwell time, promoting viral degradation in lysosomes. In situ cryo-electron tomography captured IFITM3-mediated arrest of influenza A virus membrane fusion. Observation of hemifusion diaphragms between viral particles and late endosomal membranes confirmed hemifusion stabilization as a molecular mechanism of IFITM3. The presence of the influenza fusion protein hemagglutinin in post-fusion conformation close to hemifusion sites further indicated that IFITM3 does not interfere with the viral fusion machinery. Collectively, these findings show that IFITM3 induces lipid sorting to stabilize hemifusion and prevent virus entry into target cells.

Structure and function meet at the nicotinic acetylcholine receptor-lipid interface

The nicotinic acetylcholine receptor (nAChR) is a transmembrane protein that mediates fast intercellular communication in response to the endogenous neurotransmitter acetylcholine. It is the best characterized and archetypal molecule in the superfamily of pentameric ligand-gated ion channels (pLGICs). As a typical transmembrane macromolecule, it interacts extensively with its vicinal lipid microenvironment. Experimental evidence provides a wealth of information on receptor-lipid crosstalk: the nAChR exerts influence on its immediate membrane environment and conversely, the lipid moiety modulates ligand binding, affinity state transitions and gating of ion translocation functions of the receptor protein. Recent cryogenic electron microscopy (cryo-EM) studies have unveiled the occurrence of sites for phospholipids and cholesterol on the lipid-exposed regions of neuronal and electroplax nAChRs, confirming early spectroscopic and affinity labeling studies demonstrating the close contact of lipid molecules with the receptor transmembrane segments. This new data provides structural support to the postulated "lipid sensor" ability displayed by the outer ring of M4 transmembrane domains and their modulatory role on nAChR function, as we postulated a decade ago. Borrowing from the best characterized nAChR, the electroplax (muscle-type) receptor, and exploiting new structural information on the neuronal nAChR, it is now possible to achieve an improved depiction of these sites. In combination with site-directed mutagenesis, single-channel electrophysiology, and molecular dynamics studies, the new structural information delivers a more comprehensive portrayal of these lipid-sensitive loci, providing mechanistic explanations for their ability to modulate nAChR properties and raising the possibility of targetting them in disease.

The Relationship of Cholesterol Responses to Mitochondrial Dysfunction and Lung Inflammation in Chronic Obstructive Pulmonary Disease

Hyperlipidemia is frequently reported in chronic obstructive pulmonary disease (COPD) patients and is linked to the progression of the disease and its comorbidities. Hypercholesterolemia leads to cholesterol accumulation in many cell types, especially immune cells, and some recent studies suggest that cholesterol impacts lung epithelial cells' inflammatory responses and mitochondrial responses. Several studies also indicate that targeting cholesterol responses with either statins or liver X receptor (LXR) agonists may be plausible means of improving pulmonary outcomes. Equally, cholesterol metabolism and signaling are linked to mitochondrial dysfunction and inflammation attributed to COPD progression. Here, we review the current literature focusing on the impact of cigarette smoke on cholesterol levels, cholesterol efflux, and the influence of cholesterol on immune and mitochondrial responses within the lungs.

The effect of rhamnolipids on fungal membrane models as described by their interactions with phospholipids and sterols: An in silico study

Introduction: Rhamnolipids (RLs) are secondary metabolites naturally produced by bacteria of the genera Pseudomonas and Burkholderia with biosurfactant properties. A specific interest raised from their potential as biocontrol agents for crop culture protection in regard to direct antifungal and elicitor activities. As for other amphiphilic compounds, a direct interaction with membrane lipids has been suggested as the key feature for the perception and subsequent activity of RLs. Methods: Molecular Dynamics (MD) simulations are used in this work to provide an atomistic description of their interactions with different membranous lipids and focusing on their antifungal properties. Results and discussion: Our results suggest the insertion of RLs into the modelled bilayers just below the plane drawn by lipid phosphate groups, a placement that is effective in promoting significant membrane fluidification of the hydrophobic core. This localization is promoted by the formation of ionic bonds between the carboxylate group of RLs and the amino group of the phosphatidylethanolamine (PE) or phosphatidylserine (PS) headgroups. Moreover, RL acyl chains adhere to the ergosterol structure, forming a significantly higher number of van der Waals contact with respect to what is observed for phospholipid acyl chains. All these interactions might be essential for the membranotropic-driven biological actions of RLs.

PI(4,5)P2 and Cholesterol: Synthesis, Regulation, and Functions

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) is the most abundant membrane phosphoinositide and cholesterol is an essential component of the plasma membrane (PM). Both lipids play key roles in a variety of cellular functions including as signaling molecules and major regulators of protein function. This chapter provides an overview of these two important lipids. Starting from a brief description of their structure, synthesis, and regulation, the chapter continues to describe the primary functions and signaling processes in which PI(4,5)P₂ and cholesterol are involved. While PI(4,5)P₂ and cholesterol can act independently, they often act in concert or affect each other's impact. The chapters in this volume on "Cholesterol and PI(4,5)P₂ in Vital Biological Functions: From Coexistence to Crosstalk" focus on the emerging relationship between cholesterol and PI(4,5)P₂ in a variety of biological systems and processes. In this chapter, the next section provides examples from the ion channel field demonstrating that PI(4,5)P₂ and cholesterol can act via common mechanisms. The chapter ends with a discussion of future directions.

Snapshots of ABCG1 and ABCG5/G8: A Sterol's Journey to Cross the Cellular Membranes

The subfamily-G ATP-binding cassette (ABCG) transporters play important roles in regulating cholesterol homeostasis. Recent progress in the structural data of ABCG1 and ABCG5/G8 disclose putative sterol binding sites that suggest the possible cholesterol translocation pathway. ABCG1 and ABCG5/G8 share high similarity in the overall molecular architecture, and both transporters appear to use several unique structural motifs to facilitate cholesterol transport along this pathway, including the phenylalanine highway and the hydrophobic valve. Interestingly, ABCG5/G8 is known to transport cholesterol and phytosterols, whereas ABCG1 seems to exclusively transport cholesterol. Ligand docking analysis indeed suggests a difference in recruiting sterol molecules to the known sterol-binding sites. Here, we further discuss how the different and shared structural features are relevant to their physiological functions, and finally provide our perspective on future studies in ABCG cholesterol transporters.

Orientation of Cholesterol in Polyunsaturated Lipid Bilayers

The local normal to the fluid liquid crystalline phase of the lipid membrane is an axis of motional symmetry for the molecules that make up the bilayer. The presence of cholesterol in the membrane increases not only the lipid hydrocarbon chain order but also the strength of the membrane's orienting potential. Cholesterol undergoes rapid reorientation about a diffusion axis that is roughly aligned with the long molecular axis, but there is also a slower reorientation of the diffusion axis, or "wobble", relative to the local bilayer normal. The extent of this second, slower motion depends on the degree of order of the lipids that make up the bilayer. We use ²H nuclear magnetic resonance of deuterium-labeled cholesterol to investigate quantitatively the effect of lipid chain unsaturation on cholesterol orientation in a series of phospholipid bilayers. We find that the hydrocarbon chains in membranes composed of polyunsaturated lipids are much more highly disordered than those in membranes composed of saturated lipids but that cholesterol remains aligned roughly along the bilayer normal.

Conformations of Three Types of Ultra-Long-Chain Fatty Acids in Multicomponent Lipid Bilayers

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in certain types of tissues, but their biological roles remain unknown. Here, we report how the conformation of ULCFAs depends on the length and unsaturated-bond ratio of the ultra-long chains and the composition of the host bilayer membrane using molecular dynamics simulations. The ultra-long chain of ULCFAs flips between the two leaflets and fluctuates among three conformations: elongated, L-shaped, and turned. Furthermore, we found that the saturated ultra-long chain exhibited an elongated conformation more frequently than the unsaturated chain. In addition, the truncation of the ultra-long chain at C26 had little effect on the remaining ULCFAs. ULCFAs respond to lipid-density differences in the two leaflets, and the ratio of the elongated and turned conformations changed to reduce this difference. However, in cholesterol-containing membranes, ULCFAs exhibit no density difference after the flip-flop of cholesterol removes the difference.

Distribution of cholesterol in asymmetric membranes driven by composition and differential stress

Many lipid membranes of eukaryotic cells are asymmetric, which means the two leaflets differ in at least one physical property, such as lipid composition or lateral stress. Maintaining this asymmetry is helped by the fact that ordinary phospholipids rarely transition between leaflets, but cholesterol is an exception: its flip-flop times are in the microsecond range, so that its distribution between leaflets is determined by a chemical equilibrium. In particular, preferential partitioning can draw cholesterol into a more saturated leaflet, and phospholipid number asymmetry can force it out of a compressed leaflet. Combining highly coarse-grained membrane simulations with theoretical modeling, we investigate how these two driving forces play against each other until cholesterol's chemical potential is equilibrated. The theory includes two coupled elastic sheets and a Flory-Huggins mixing free energy with a χ parameter. We obtain a relationship between χ and the interaction strength between cholesterol and lipids in either of the two leaflets, and we find that it depends, albeit weakly, on lipid number asymmetry. The differential stress measurements under various asymmetry conditions agree with our theoretical predictions. Using the two kinds of asymmetries in combination, we find that it is possible to counteract the phospholipid number bias, and the resultant stress in the membrane, via the control of cholesterol mixing in the leaflets.

Improving our picture of the plasma membrane: Rafts induce ordered domains in a simplified model cytoplasmic leaflet

By study of asymmetric membranes, models of the cell plasma membrane (PM) have improved, with more realistic properties of the asymmetric lipid composition of the membrane being explored. We used hemifusion of symmetric giant unilamellar vesicles (GUVs) with a supported lipid bilayer (SLB) to engineer bilayer leaflets of different composition. During hemifusion, only the outer leaflets of GUV and SLB are connected, exchanging lipids by simple diffusion. aGUVs were detached from the SLB for study. In general these aGUVs are formed with one leaflet that phase-separates into Ld (liquid disordered) + Lo (liquid ordered) phases, and another leaflet with lipid composition that would form a single fluid phase in a symmetric bilayer. We observed that ordered phases of either Lo or Lβ (gel phase) induce an ordered domain in the apposed fluid leaflet that lacks high melting lipids. Results suggest both an inter-leaflet and an intra-leaflet redistribution of cholesterol. We used C-Laurdan spectral images to investigate the lipid packing/order of aGUVs, finding that cholesterol partitions into the induced ordered domains. We suggest this behavior to be commonplace, that when Ld + Lo phase separation occurs in a cell PM exoplasmic leaflet, an induced order domain forms in the cytoplasmic leaflet.

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.

Full scale structural, mechanical and dynamical properties of HIV-1 liposomes

Enveloped viruses are enclosed by a lipid membrane inside of which are all of the components necessary for the virus life cycle; viral proteins, the viral genome and metabolites. Viral envelopes are lipid bilayers that adopt morphologies ranging from spheres to tubes. The envelope is derived from the host cell during viral replication. Thus, the composition of the bilayer depends on the complex constitution of lipids from the host-cell's organelle(s) where assembly and/or budding of the viral particle occurs. Here, molecular dynamics (MD) simulations of authentic, asymmetric HIV-1 liposomes are used to derive a unique level of resolution of its full-scale structure, mechanics and dynamics. Analysis of the structural properties reveal the distribution of thicknesses of the bilayers over the entire liposome as well as its global fluctuations. Moreover, full-scale mechanical analyses are employed to derive the global bending rigidity of HIV-1 liposomes. Finally, dynamical properties of the lipid molecules reveal important relationships between their 3D diffusion, the location of lipid-rafts and the asymmetrical composition of the envelope. Overall, our simulations reveal complex relationships between the rich lipid composition of the HIV-1 liposome and its structural, mechanical and dynamical properties with critical consequences to different stages of HIV-1's life cycle.

Molecular Organization of a Raft-like Domain in a Polyunsaturated Phospholipid Bilayer: A Supervised Machine Learning Analysis of Molecular Dynamics Simulations

Numerous health benefits are associated with omega-3 polyunsaturated fatty acids (n-3 PUFA) consumed in fish oils. An understanding of the mechanism remains elusive. The plasma membrane as a site of action is the focus in this study. With large-scale all-atom MD simulations run on a model membrane (1050 lipid molecules), we observed the evolution over time (6 μs) of a circular (raft-like) domain composed of N-palmitoylsphingomyelin (PSM) and cholesterol embedded into a surrounding (non-raft) patch composed of polyunsaturated 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) (1:1:1 mol). A supervised machine learning algorithm was developed to characterize the migration of each lipid based on molecular conformation and the local environment. PDPC molecules were seen to infiltrate the ordered raft-like domain in a small amount, while a small concentration of PSM and cholesterol molecules was seen to migrate into the disordered non-raft region. Enclosing the raft-like domain, a narrow (∼2 nm in width) interfacial zone composed of PDPC, PSM, and cholesterol that buffers the substantial difference in order (ΔS_CD ≈ 0.12) between raft-like and non-raft environments was seen to form. Our results suggest that n-3 PUFA regulate the architecture of lipid rafts enriched in sphingolipids and cholesterol with a minimal effect on order within their interior in membranes.

Deciphering lipid transfer between and within membranes with time-resolved small-angle neutron scattering

This review focuses on time-resolved neutron scattering, particularly time-resolved small angle neutron scattering (TR-SANS), as a powerful in situ noninvasive technique to investigate intra- and intermembrane transport and distribution of lipids and sterols in lipid membranes. In contrast to using molecular analogues with potentially large chemical tags that can significantly alter transport properties, small angle neutron scattering relies on the relative amounts of the two most abundant isotope forms of hydrogen: protium and deuterium to detect complex membrane architectures and transport processes unambiguously. This review discusses advances in our understanding of the mechanisms that sustain lipid asymmetry in membranes-a key feature of the plasma membrane of cells-as well as the transport of lipids between membranes, which is an essential metabolic process.

The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study

Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO₂/TiO₂) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO₂ or TiO₂ experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.

Interdomain exchange and the flip-flop of cholesterol in ternary component lipid membranes and their effects on heterogeneous cholesterol diffusion

Cell membranes are heterogeneous with a variety of lipids, cholesterol, and proteins and are composed of domains of different compositions. Such heterogeneous environments make the transport of cholesterol complicated: cholesterol not only diffuses within a particular domain but also travels between domains. Cholesterol also flip-flops between upper and lower leaflets such that cholesterol may reside both within leaflets and in the central region between two leaflets. How the presence of multiple domains and the interdomain exchange of cholesterol would affect the cholesterol transport, however, remains elusive. In this study, therefore, we perform molecular dynamics simulations up to 100μs for ternary component lipid membranes, which consist of saturated lipids (dipalmitoylphosphatidylcholine, DPPC), unsaturated lipids (dilinoleylphosphatidylcholine, DIPC), and cholesterol. The ternary component membranes in our simulations form two domains readily: DPPC and DIPC domains. We find that the diffusion of cholesterol molecules is much more heterogeneous and non-Gaussian than expected for binary component lipid membranes of lipids and cholesterol. The non-Gaussian parameter of the cholesterol molecules is about four times larger in the ternary component lipid membranes than in the binary component lipid membranes. Such non-Gaussian and heterogeneous transport of cholesterol arises from the interplay among the interdomain kinetics, the different diffusivity of cholesterol in different domains, and the flip-flop of cholesterol. This suggests that in cell membranes that consist of various domains and proteins, the cholesterol transport can be very heterogeneous. We also find that the mechanism of the interdomain exchange differs for different domains: cholesterol tends to exit the DIPC domain along the central region of the membrane for the DIPC-to-DPPC transition, while the cholesterol is likely to exit the DPPC domain within the membrane leaflet for the DPPC-to-DIPC transition. Also, the interdomain exchange kinetics of cholesterol for the DPPC-to-DIPC transition is up to 7.9 times slower than the DIPC-to-DPPC transition.

Role of cholesterol flip-flop in oxidized lipid bilayers

We performed a series of molecular dynamics simulations of cholesterol (Chol) in nonoxidized 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and in binary mixtures of PLPC-oxidized-lipid-bilayers with 0-50% Chol concentration and oxidized lipids with hydroperoxide and aldehyde oxidized functional groups. From the 60 unbiased molecular dynamics simulations (total of 161 μs), we found that Chol inhibited pore formation in the aldehyde-containing oxidized lipid bilayers at concentrations greater than 11%. For both pure PLPC bilayer and bilayers with hydroperoxide lipids, no pores were observed at any Chol concentration. Furthermore, increasing cholesterol concentration led to a change of phase state from the liquid-disordered to the liquid-ordered phase. This condensing effect of Chol was observed in all systems. Data analysis shows that the addition of Chol results in an increase in bilayer thickness. Interestingly, we observed Chol flip-flop only in the aldehyde-containing lipid bilayer but neither in the PLPC nor the hydroperoxide bilayers. Umbrella-sampling simulations were performed to calculate the translocation free energies and the Chol flip-flop rates. The results show that Chol's flip-flop rate depends on the lipid bilayer type, and the highest rate are found in aldehyde bilayers. As the main finding, we shown that Chol stabilizes the oxidized lipid bilayer by confining the distribution of the oxidized functional groups.

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial β-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

LiPyphilic: A Python Toolkit for the Analysis of Lipid Membrane Simulations

Molecular dynamics simulations are now widely used to study emergent phenomena in lipid membranes with complex compositions. Here, we present LiPyphilic-a fast, fully tested, and easy-to-install Python package for analyzing such simulations. Analysis tools in LiPyphilic include the identification of cholesterol flip-flop events, the classification of local lipid environments, and the degree of interleaflet registration. LiPyphilic is both force field- and resolution-agnostic, and by using the powerful atom selection language of MDAnalysis, it can handle membranes with highly complex compositions. LiPyphilic also offers two on-the-fly trajectory transformations to (i) fix membranes split across periodic boundaries and (ii) perform nojump coordinate unwrapping. Our implementation of nojump unwrapping accounts for fluctuations in the box volume under the NPT ensemble-an issue that most current implementations have overlooked. The full documentation of LiPyphilic, including installation instructions and links to interactive online tutorials, is available at https://lipyphilic.readthedocs.io/en/latest.

Nonconverged Constraints Cause Artificial Temperature Gradients in Lipid Bilayer Simulations

Molecular dynamics (MD) simulations have become an indispensable tool to investigate phase separation in model membrane systems. In particular, simulations based on coarse-grained (CG) models have found widespread use due to their increased computational efficiency, allowing for simulations of multicomponent lipid bilayers undergoing phase separation into liquid-ordered and liquid-disordered domains. Here, we show that a significant temperature difference between molecule types can artificially arise in CG MD membrane simulations with the standard Martini simulation parameters in GROMACS. In particular, the linear constraint solver (LINCS) algorithm does not converge with its default settings, resulting in serious temperature differences between molecules in a time step-dependent manner. We demonstrate that the underlying reason for this behavior is the presence of highly constrained moieties, such as cholesterol. Their presence can critically impact numerous structural and dynamic membrane properties obtained from such simulations. Furthermore, any preference of these molecules toward a certain membrane phase can lead to spatial temperature gradients, which can amplify the degree of phase separation or even induce it in compositions that would otherwise mix well. We systematically investigated the effect of the integration time step and LINCS settings on membrane properties. Our data show that for cholesterol-containing membranes, a time step of 20 fs should be combined with at least lincs_iter = 2 and lincs_order = 12, while using a time step of 30 fs requires at least lincs_iter = 3 and lincs_order = 12 to bring the temperature differences to a level where they do not perturb central membrane properties. Moreover, we show that in cases where stricter LINCS settings are computationally too demanding, coupling the lipids in multiple groups to the temperature bath offers a practical workaround to the problem, although the validity of this approach should be further verified. Finally, we show that similar temperature gradients can also emerge in atomistic simulations using the CHARMM force field in combination with settings that allow for a 5 fs integration step.

The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane

A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

Cholesterol activates BK channels by increasing KCNMB1 protein levels in the plasmalemma

Calcium-/voltage-gated, large-conductance potassium channels (BKs) control critical physiological processes, including smooth muscle contraction. Numerous observations concur that elevated membrane cholesterol (CLR) inhibits the activity of homomeric BKs consisting of channel-forming alpha subunits. In mammalian smooth muscle, however, native BKs include accessory KCNMB1 (β1) subunits, which enable BK activation at physiological intracellular calcium. Here, we studied the effect of CLR enrichment on BK currents from rat cerebral artery myocytes. Using inside-out patches from middle cerebral artery (MCA) myocytes at [Ca2+]free=30 μM, we detected BK activation in response to in vivo and in vitro CLR enrichment of myocytes. While a significant increase in myocyte CLR was achieved within 5 min of CLR in vitro loading, this brief CLR enrichment of membrane patches decreased BK currents, indicating that BK activation by CLR requires a protracted cellular process. Indeed, blocking intracellular protein trafficking with brefeldin A (BFA) not only prevented BK activation but led to channel inhibition upon CLR enrichment. Surface protein biotinylation followed by Western blotting showed that BFA blocked the increase in plasmalemmal KCNMB1 levels achieved via CLR enrichment. Moreover, CLR enrichment of arteries with naturally high KCNMB1 levels, such as basilar and coronary arteries, failed to activate BK currents. Finally, CLR enrichment failed to activate BK channels in MCA myocytes from KCNMB1-/- mouse while activation was detected in their wild-type (C57BL/6) counterparts. In conclusion, the switch in CLR regulation of BK from inhibition to activation is determined by a trafficking-dependent increase in membrane levels of KCNMB1 subunits.

Effect of alcohol on the phase separation in model membranes

The discovery of coexisting liquid-ordered and liquid-disordered phases in multicomponent lipid bilayers has received widespread attention due to its potential relevance for biological systems. One of the many open questions is how the presence of additional components affects the nature of the coexisting phases. Of particular interest is the addition of alcohols because their anesthetic properties may arise from modulating bilayer behavior. We use coarse-grained Molecular Dynamics simulations to gain insight into the partitioning preferences of linear n-alcohols into ordered and disordered bilayers alongside their effects on local membrane structure. We find that alcohols cause only small changes to membrane composition alongside a lack of significant effects on membrane thickness and lipid tail order. Cholesterol and n-alcohol trans-bilayer motion is measured and found to be near or within the range of previous atomistic results. The cholesterol flip-flop rates increase with both n-alcohol length and concentration for octanol, dodecanol, and hexadecanol, indicating a decrease in lipid order. Umbrella sampling simulations of removing cholesterol from tertiary membranes find no significant difference with or without n-alcohols at various concentrations. Simulations of a phase-separated bilayer show that octanol preferentially partitions into the liquid-disordered phase in a ratio of approximately 3:1 over the liquid-ordered phase. Furthermore, partition coefficients of alcohol in single-phase membranes show a preference for longer alcohols (dodecanol and hexadecanol) to partition preferentially into the liquid-ordered phase, while decreasing the length of the alcohol reverses this trend. Our work tests experimental results while also investigating the ability for coarse-grained MARTINI simulations to capture minute differences in model membrane spatial arrangements on the nanoscale level.

Active cholesterol 20 years on

This review considers the following hypotheses, some well-supported and some speculative. Almost all of the sterol molecules in plasma membranes are associated with bilayer phospholipids in complexes of varied strength and stoichiometry. These complexes underlie many of the material properties of the bilayer. The small fraction of cholesterol molecules exceeding the binding capacity of the phospholipids is thermodynamically active and serves diverse functions. It circulates briskly among the cell membranes, particularly through contact sites linking the organelles. Active cholesterol provides the upstream feedback signal to multiple mechanisms governing plasma membrane homeostasis, pegging the sterol level to a threshold set by its phospholipids. Active cholesterol could also be the cargo for various inter-organelle transporters and the form excreted from cells by reverse transport. Furthermore, it is integral to the function of caveolae; a mediator of Hedgehog regulation; and a ligand for the binding of cytolytic toxins to membranes. Active cholesterol modulates a variety of plasma membrane proteins-receptors, channels and transporters-at least in vitro.

Toward a Membrane-Centric Biology

With advancements of modern biophysical tools and superresolution imaging, cell biology is entering a new phase of research with technological power fitting for membrane dynamics analyses. However, our current knowledge base of cellular signaling events is mostly built on a network of protein interactions, which is incompatible with the essential roles of membrane activities in those events. The lack of a theoretical platform is rendering biophysical analyses of membrane biology supplementary to the protein-centric paradigm. We hypothesize a framework of signaling events mediated by lipid dynamics and argue that this is the evolutionarily obligatory developmental path of cellular complexity buildup. In this framework, receptors are the late comers, integrating into the pre-existing membrane based signaling events using their lipid interface as the point of entry. We further suggest that the reason for cell surface receptors to remain silent at the resting state is via the suppression effects of their surrounding lipids. The avoidance of such a suppression, via ligand binding or lipid domain disruption, enables the receptors to autonomously integrate themselves into the preexisting networks of signaling cascades.

Composition dependence of cholesterol flip-flop rates in physiological mixtures

Unbiased, all-atom simulations of mixtures representative of the inner and outer leaflets of a mammalian red blood cell and a synaptic vesicle reveal many cholesterol flip-flop events over the 5 μsec duration of the simulations. Enough events are observed for a direct estimate of the flip-flop rate. Slower rates are found in more ordered membranes, and faster rates in more disordered membranes, consistent with earlier reports in the literature. However, the rates found here are neither as fast as the fastest nor as slow as the slowest rates obtained by previous simulations. The difference likely stems from the compositions studied here, which unlike previous work include exclusively lipids with differing acyl chains, as observed in mammalian lipidomes.

On the lipid flip-flop and phase transition coupling

We measured the passive lipid flip-flop of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in solid supported lipid bilayers across their main gel to fluid (Lβ → Lα) phase transition. By performing time and temperature resolved neutron reflectometry experiments, we demonstrated that asymmetric systems prepared in the gel phase are stable for at least 24 hours. Lipid flip-flop was found to be intrinsically linked to the amount of lipid molecules in the fluid phase. Moreover, the  increase of this amount during the broad phase transition was found to be the main key factor for the timing of the flip-flop process. By measuring different temperature scan rate, we could demonstrate that, in the case of supported bilayers and for the temperature investigated, the lipid flip flop is characterised by an activation energy of 50 kJ mol-1 and a timescale on the order of few hours. Our results demonstrate the origin on the discrepancies between passive flip-flop in bulk systems and at interfaces.

Flip-flop asymmetry of cholesterol in model membranes induced by thermal gradients

Lipid asymmetry is a crucial property of biological membranes and significantly influences their physical and mechanical properties. It is responsible for maintaining different chemical environments on the external and internal surfaces of cells and organelles and plays a vital role in many biological processes such as cell signalling and budding. In this work we show, using non-equilibrium molecular dynamics (NEMD) simulations, that thermal fields can induce lipid asymmetry in biological membranes. We focus our investigation on cholesterol, an abundant lipid in the plasma membrane, with a rapid flip-flop rate, significantly influencing membrane properties. We demonstrate that thermal fields induce membrane asymmetry with cholesterol showing thermophobic behaviour and therefore accumulating on the cold side of the membrane. This work highlights a possible experimental route to preparing and controlling asymmetry in synthetic membranes.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.