skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on April 10, 2026

Title: Molecular Insights into Cholesterol Concentration Effects on Planar and Curved Lipid Bilayers for Liposomal Drug Delivery
Abstract Liposomal carriers provide a flexible and effective strategy for delivering therapeutics across a broad spectrum of diseases. Cholesterol is frequently included in these systems to improve membrane rigidity and limit permeability. Despite its widespread use, the optimal cholesterol-to-lipid proportion for achieving stable and efficient liposome performance remains to be fully determined. In this work, we apply all-atom molecular dynamics simulations to explore how different cholesterol concentrations influence the structural and dynamic characteristics of 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) bilayers, considering both planar and curved membrane geometries. Bilayers with cholesterol molar ratios of 0%, 10%, 20%, 30%, 40%, and 50% were simulated, and key biophysical parameters including area per lipid (APL), membrane thickness, leaflet interdigitation, and deuterium order parameters (SCD) were analyzed. In planar bilayers, increasing cholesterol concentration led to a progressive decrease in APL from approximately 60°A2to 40°A2, accompanied by increased membrane thickness and lipid ordering, consistent with cholesterol’s classical condensing effect. In contrast, curved bilayers exhibited a cholesterol-induced expansion effect, particularly in the inner leaflet, where APL increased from approximately 60°A2to 90°A2with rising cholesterol levels. SCD profiles showed that cholesterol enhanced tail ordering up to 40% concentration, beyond which the effect plateaued or slightly declined, suggesting structural saturation or packing frustration. Membrane thickness displayed a monotonic increase in planar bilayers but followed a nonlinear trend in curved systems due to curvature-induced stress. These findings highlight that cholesterol’s influence on membrane properties is highly dependent on bilayer geometry and asymmetry. While planar bilayers exhibit predictable responses, curved systems reveal nonclassical behaviors that challenge traditional models of cholesterol-lipid interactions. This work provides molecular-level insights and establishes a computational framework for the rational design of liposomal systems, emphasizing the need to account for curvature and asymmetry in membrane engineering.  more » « less
Award ID(s):
1945465
PAR ID:
10600747
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
bioRxiv
Date Published:
Format(s):
Medium: X
Institution:
bioRxiv
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The plasma membrane of eukaryotic cells is known to be compositionally asymmetric. Certain phospholipids, such as sphingomyelin and phosphatidylcholine species, are predominantly localized in the outer leaflet, while phosphatidylethanolamine and phosphatidylserine species primarily reside in the inner leaflet. While phospholipid asymmetry between the membrane leaflets is well established, there is no consensus about cholesterol distribution between the two leaflets. We have performed a systematic study, via molecular simulations, of how the spatial distribution of cholesterol molecules in different “asymmetric” lipid bilayers are affected by the lipids’ backbone, head-type, unsaturation, and chain-length by considering an asymmetric bilayer mimicking the plasma membrane lipids of red blood cells, as well as seventeen other asymmetric bilayers comprising of different lipid types. Our results reveal that the distribution of cholesterol in the leaflets is solely a function of the extent of ordering of the lipids within the leaflets. The ratio of the amount of cholesterol matches the ratio of lipid order in the two leaflets, thus providing a quantitative relationship between the two. These results are understood by the observation that asymmetric bilayers with equimolar amount of lipids in the two leaflets develop tensile and compressive stresses due to differences in the extent of lipid order. These stresses are alleviated by the transfer of cholesterol from the leaflet in compressive stress to the one in tensile stress. These findings are important in understanding the biology of the cell membrane, especially with regard to the composition of the membrane leaflets. 
    more » « less
  2. Many cellular lipid bilayers consist of leaflets that differ in their lipid composition — a non-equilibrium state actively maintained by cellular sorting processes that counter passive lipid flip-flop. While this lipidomic aspect of membrane asymmetry has been known for half a century, its elastic and thermodynamic ramifications have garnered attention only fairly recently. Notably, the torque arising when lipids of different spontaneous curvature reside in the two leaflets can be counterbalanced by a difference in lateral mechanical stress between them. Such membranes can be essentially flat in their relaxed state, despite being compositionally strongly asymmetric, but they harbor a surprisingly large but macroscopically invisible differential stress. This hidden stress can affect a wide range of other membrane properties, such as the resistance to bending, the nature of phase transitions in its leaflets, and the distribution of flippable species, most notably sterols. In this short note we offer a concise overview of our recently proposed basic framework for capturing the interplay between curvature, lateral stress, leaflet phase behavior, and cholesterol distribution in generally asymmetric membranes, and how its implied signatures might be used to learn more about the hidden but physically consequential differential stress. 
    more » « less
  3. Phosphatidic acid (PA) is a signaling lipid that is produced enzymatically from phosphatidylcholine (PC), lysophosphatidic acid, or diacylglycerol. Compared to PC, PA lacks a choline moiety on the headgroup, making the headgroup smaller than that of PC and PA, and PA has a net negative charge. Unlike the cylindrical geometry of PC, PA, with its small headgroup relative to the two fatty acid tails, is proposed to support negatively curved membranes. Thus, PA is thought to play a role in a variety of biological processes that involve bending membranes, such as the formation of intraluminal vesicles in multivesicular bodies and membrane fusion. Using supported tubulated lipid bilayers (STuBs), the extent to which PA localizes to curved membranes was determined. STuBs were created via liposome deposition with varying concentrations of NaCl (500 mM to 1 M) on glass to form supported bilayers with connected tubules. The location of fluorescently labeled lipids relative to tubules was determined by imaging with total internal reflection or confocal fluorescence microscopy. The accumulation of various forms of PA (with acyl chains of 16:0-6:0, 16:0-12:0, 18:1-12:0) were compared to PC and the headgroup labeled phosphatidylethanolamine (PE), a lipid that has been shown to accumulate at regions of curvature. PA and PE accumulated more at tubules and led to the formation of more tubules than PC. Using large unilamellar liposomes in a dye-quenching assay, the location of the headgroup labeled PE was determined to be mostly on the outer, positively curved leaflet, whereas the tail labeled PA was located more on the inner, negatively curved leaflet. This study demonstrates that PA localizes to regions of negative curvature in liposomes and supports the formation of curved, tubulated membranes. This is one way that PA could be involved with curvature formation during a variety of cell processes. 
    more » « less
  4. null (Ed.)
    It has become increasingly apparent that the lipid composition of cell membranes affects the function of transmembrane proteins such as ion channels. Here, we leverage the structural and functional diversity of small viral K+ channels to systematically examine the impact of bilayer composition on the pore module of single K+ channels. In vitro–synthesized channels were reconstituted into phosphatidylcholine bilayers ± cholesterol or anionic phospholipids (aPLs). Single-channel recordings revealed that a saturating concentration of 30% cholesterol had only minor and protein-specific effects on unitary conductance and gating. This indicates that channels have effective strategies for avoiding structural impacts of hydrophobic mismatches between proteins and the surrounding bilayer. In all seven channels tested, aPLs augmented the unitary conductance, suggesting that this is a general effect of negatively charged phospholipids on channel function. For one channel, we determined an effective half-maximal concentration of 15% phosphatidylserine, a value within the physiological range of aPL concentrations. The different sensitivity of two channel proteins to aPLs could be explained by the presence/absence of cationic amino acids at the interface between the lipid headgroups and the transmembrane domains. aPLs also affected gating in some channels, indicating that conductance and gating are uncoupled phenomena and that the impact of aPLs on gating is protein specific. In two channels, the latter can be explained by the altered orientation of the pore-lining transmembrane helix that prevents flipping of a phenylalanine side chain into the ion permeation pathway for long channel closings. Experiments with asymmetrical bilayers showed that this effect is leaflet specific and most effective in the inner leaflet, in which aPLs are normally present in plasma membranes. The data underscore a general positive effect of aPLs on the conductance of K+ channels and a potential interaction of their negative headgroup with cationic amino acids in their vicinity. 
    more » « less
  5. Membrane undulations play a vital role in many biological processes, including the regulation of membrane protein activity. The asymmetric lipid composition of most biological membranes complicates theoretical description of these bending fluctuations, yet experimental data that would inform any such a theory is scarce. Here, we used neutron spin-echo (NSE) spectroscopy to measure the bending fluctuations of large unilamellar vesicles (LUV) having an asymmetric transbilayer distribution of high- and low-melting lipids. The asymmetric vesicles were prepared using cyclodextrin-mediated lipid exchange, and were composed of an outer leaflet enriched in egg sphingomyelin (ESM) and an inner leaflet enriched in 1-palmitoyl-2-oleoyl-phosphoethanolamine (POPE), which have main transition temperatures of 37 °C and 25 °C, respectively. The overall membrane bending rigidity was measured at three temperatures: 15 °C, where both lipids are in a gel state; 45 °C, where both lipids are in a fluid state; and 30 °C, where there is gel-fluid co-existence. Remarkably, the dynamics for the fluid asymmetric LUVs (aLUVs) at 30 °C and 45 °C do not follow trends predicted by their symmetric counterparts. At 30 °C, compositional asymmetry suppressed the bending fluctuations, with the asymmetric bilayer exhibiting a larger bending modulus than that of symmetric bilayers corresponding to either the outer or inner leaflet. We conclude that the compositional asymmetry and leaflet coupling influence the internal dissipation within the bilayer and result in membrane properties that cannot be directly predicted from corresponding symmetric bilayers. 
    more » « less