Complex Lipid Bilayers Research Articles

Model Mechanism for Lipid Uptake by the Human STARD2/PC-TP Phosphatidylcholine Transfer Protein.

The human StAR-related lipid transfer domain protein 2 (STARD2), also known as phosphatidylcholine (PC) transfer protein, is a single-domain lipid transfer protein thought to transfer PC lipids between intracellular membranes. We performed extensive μs-long molecular dynamics simulations of STARD2 of its apo and holo forms in the presence or absence of complex lipid bilayers. The simulations in water reveal ligand-dependent conformational changes. In the 2 μs-long simulations of apo STARD2 in the presence of a lipid bilayer, we observed spontaneous reproducible PC lipid uptake into the protein hydrophobic cavity. We propose that the lipid extraction mechanism involves one to two metastable states stabilized by choline-tyrosine or choline-tryptophane cation-π interactions. Using free energy perturbation, we evaluate that PC-tyrosine cation-π interactions contribute 1.8 and 2.5 kcal/mol to the affinity of a PC-STARD2 metastable state, thus potentially providing a significant decrease of the energy barrier required for lipid desorption.

Membrane specificity of the human cholesterol transfer protein STARD4

STARD4 regulates cholesterol homeostasis by transferring cholesterol between the plasma membrane and endoplasmic reticulum. The STARD4 structure features a helix-grip fold surrounding a large hydrophobic cavity holding the sterol. Its access is controlled by a gate formed by the Ω1 and Ω4 loops and the C-terminal α-helix. Little is known about the mechanisms by which STARD4 binds to membranes and extracts/releases cholesterol. All available structures of STARD4 are without a bound sterol and display the same closed conformation of the gate. The cholesterol transfer activity of the mouse STARD4 is enhanced in the presence of anionic lipids, and in particular of phosphatidylinositol biphosphates (PIP2) for which two binding sites were proposed on the mouse STARD4 surface. Yet only one of these sites is conserved in human STARD4. We here report the results of a liposome microarray-based assay and microseconds-long molecular dynamics simulations of human STARD4 with complex lipid bilayers mimicking the composition of the donor and acceptor membranes. We show that the binding of apo form of human STARD4 is sensitive to the presence of PIP2 through two specific binding sites, one of which was not identified on mouse STARD4. We report two novel conformations of the gate in holo-STARD4: a yet-unobserved close conformation and an open conformation of Ω4 shedding light on the opening/closure mechanism needed for cholesterol uptake/release. Overall, the modulation of human STARD4 membrane-binding by lipid composition, and by the presence of the cargo supports the capacity of human STARD4 to achieve directed transfer between specific organelle membranes.

Interaction of detergent with complex mimics of bacterial membranes

Detergents are valuable tools to extract membrane proteins for biophysical, biochemical, and structural scrutiny. The detergent-driven solubilization of bilayers made from a single lipid species is commonly described in terms of pseudo-phase diagrams and a three-stage model accounting for three ranges comprising (i) intact vesicles, (ii) vesicle/micelle co-existence, or (iii) mixed micelles. Moreover, the pseudo-phase boundaries thus determined can often be quantitatively rationalized in terms of the molecular shapes of the lipid and the detergent used. Yet, it has remained unclear to what extent this approach can be applied to multi-component lipid membranes that more closely mimic the compositional complexity of cellular membranes. Here, we studied how lipid mixtures composed of palmitoyl oleoyl phosphatidylethanolamine (POPE), palmitoyl oleoyl phosphatidylglycerol (POPG), and tetraoleoyl cardiolipin (TOCL) are solubilized by the commonly used zwitterionic detergent lauryldimethylamine N-oxide using isothermal titration calorimetry. While phase diagrams of the diverse lipid mixtures showed the typical ranges of the three-stage model, we found that POPG-rich POPE/POPG bilayers are more difficult to solubilize than POPG-poor POPE/POPG bilayers. In turn, POPE/POPG/TOCL bilayers became increasingly resistant to detergent with increasing TOCL content. Since POPG is nearly cylindrically shaped and TOCL adopts inverted cone-like shapes under current buffer conditions, our solubilization data do not align with shape-based arguments. Instead, additional electrostatic interactions between lipids and detergents lead to non-additive mixing behavior affecting the resilience of complex lipid bilayers against solubilization.

Modelling the effects of E/Z photoisomerization of a cyclocurcumin analogue on the properties of cellular lipid membranes.

The use of photosensitive molecules capable of isomerizing under light stimuli, and thus induce perturbation in biological systems, is becoming increasingly popular for potential light-activated chemotherapeutic purposes. We recently show that a cyclocurcumin derivative (CCBu), may be suitable for light-activated chemotherapy and may constitute a valuable alternative to traditional photodynamic therapy, due to its oxygen-independent mechanism of action, which allows the treatment of hypoxic solid tumors. In particular, we have shown that the E/Z photoisomerization of CCBu correlates with strong perturbations of model lipid bilayers. In this work, we perform all-atom classical molecular dynamics for a more complex bilayer, whose composition is, thus, much closer to eukaryotic outer cell membranes. We have evidenced important differences in the interaction pathway between CCBu and the complex lipid bilayer as compared to previous models, concerning both the membrane penetration capacity and the isomerization-induced perturbations. While we confirm that structural perturbations of the lipid membrane are induced by isomerization, we also show how the use of a simplified membrane model can result in an oversimplification of the system and hinder key physical and biological phenomena. Although, CCBu may be considered as a suitable candidate for light-activated chemotherapy, we also underline how the inclusion of bulkier substituents, inducing larger perturbations upon photoisomerization, may enhance its efficiency.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure–function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck’s activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck’s conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures

Membrane reshaping is an essential biological process. The chemical composition of lipid membranes determines their mechanical properties and thus the energetics of their shape. Hundreds of distinct lipid species make up native bilayers, and this diversity complicates efforts to uncover what compositional factors drive membrane stability in cells. Simplifying assumptions, therefore, are used to generate quantitative predictions of bilayer dynamics based on lipid composition. One assumption commonly used is that “per lipid” mechanical properties are both additive and constant—that they are an intrinsic property of lipids independent of the surrounding composition. Related to this is the assumption that lipid bulkiness, or “shape,” determines its curvature preference, independently of context. In this study, all-atom molecular dynamics simulations on three separate multilipid systems were used to explicitly test these assumptions, applying methodology recently developed to isolate properties of single lipids or nanometer-scale patches of lipids. The curvature preference experienced by populations of lipid conformations were inferred from their redistribution on a dynamically fluctuating bilayer. Representative populations were extracted by both structural similarity and semi-automated hidden Markov model analysis. The curvature preferences of lipid dimers were then determined and compared with an additive model that combines the monomer curvature preference of both the individual lipids. In all three systems, we identified conformational subpopulations of lipid dimers that showed non-additive curvature preference, in each case mediated by a special chemical interaction (e.g., hydrogen bonding). Our study highlights the importance of specific chemical interactions between lipids in multicomponent bilayers and the impact of interactions on bilayer stiffness. We identify two mechanisms of bilayer softening: diffusional softening, driven by the dynamic coupling between lipid distributions and membrane undulations, and conformational softening, driven by the inter-conversion between distinct dimeric conformations.

Nanoscale Features of Tunable Bacterial Outer Membrane Models Revealed by Correlative Microscopy.

The rise of antibiotic resistance is a growing worldwide human health issue, with major socioeconomic implications. An understanding of the interactions occurring at the bacterial membrane is crucial for the generation of new antibiotics. Supported lipid bilayers (SLBs) made from reconstituted lipid vesicles have been used to mimic these membranes, but their utility has been restricted by the simplistic nature of these systems. A breakthrough in the field has come with the use of outer membrane vesicles derived from Gram-negative bacteria to form SLBs, thus providing a more physiologically relevant system. These complex bilayer systems hold promise but have not yet been fully characterized in terms of their composition, ratio of natural to synthetic components, and membrane protein content. Here, we use correlative atomic force microscopy (AFM) with structured illumination microscopy (SIM) for the accurate mapping of complex lipid bilayers that consist of a synthetic fraction and a fraction of lipids derived from Escherichia coli outer membrane vesicles (OMVs). We exploit the high resolution and molecular specificity that SIM can offer to identify areas of interest in these bilayers and the enhanced resolution that AFM provides to create detailed topography maps of the bilayers. We are thus able to understand the way in which the two different lipid fractions (natural and synthetic) mix within the bilayers, and we can quantify the amount of bacterial membrane incorporated into the bilayer. We prove the system’s tunability by generating bilayers made using OMVs engineered to contain a green fluorescent protein (GFP) binding nanobody fused with the porin OmpA. We are able to directly visualize protein–protein interactions between GFP and the nanobody complex. Our work sets the foundation for accurately understanding the composition and properties of OMV-derived SLBs to generate a high-resolution platform for investigating bacterial membrane interactions for the development of next-generation antibiotics.

Serotonergic drugs modulate the phase behavior of complex lipid bilayers

Serotonin is an endogenous neurotransmitter involved in both physiological and pathophysiological processes. Traditionally, serotonin acts as a ligand for G protein-coupled receptors (GPCRs) leading to subsequent cell signaling. However, serotonin can also bind to lipid membranes with high affinity and modulate the phase behavior in 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC)/N-palmitoyl-D-erythro-sphingosylphosphorylcholine (PSM)/cholesterol model membranes mimicking the outer leaflet of the plasma membrane. Here, we investigated if serotonergic drugs containing the pharmacophore from serotonin would also modulate phase behavior in lipid membranes in a similar fashion. We used 2H NMR spectroscopy to explore the phase behavior of POPC/PSM/cholesterol (4/4/2 molar ratio) mixtures in the presence of the serotonergic drugs aripiprazole, BRL-54443, BW-723C86, and CP-135807 at a lipid to drug molar ratio of 10:1. POPC and PSM were perdeuterated in the palmitoyl chain, respectively, and prepared in individual samples. Numerical lineshape simulations of the 2H NMR spectra were used to calculate the order parameter profiles and projected lengths of the saturated acyl chains. All serotonergic drugs induce two components in 2H NMR spectra, indicating that they increased the hydrophobic mismatch between the thickness of the coexisting lipid phases leading to larger domain sizes, relatively similarly to serotonin. AFM force indentation and Raman spectral studies, which interrogate membrane mechanical properties, also indicate changes in membrane order in the presence of these drugs. These findings highlight how serotonergic drugs alter membrane phase behavior and could modulate both target and other membrane proteins, possibly explaining the side effects observed for serotonergic and other clinically relevant drugs.

Ag nanoparticle interaction with complex lipidbilayer: an atomistic and coarse-grained study

Understanding nanoparticle (NP) translocation through complex cellular membranes is required for pre-evaluating their toxicity. Silver (Ag) is a common material, extensively applied in NP form for its apparent anti-bacterial properties. While so widely available, data about their interaction with lipid membranes, particularly the role of size and surface functionalization in membrane binding, is relatively scarce. As a first step, we performed all-atom biased and unbiased molecular dynamic simulations for bare Ag NPs for three sizes (1/3/5 nm diameter) with a complex membrane of DPPC:DOPC:POPC:Chol in 5:2:2:1 representing model lung membrane.

Formulation Development and Evaluation of Transfero-somes Nano Gel Loaded with Levocetirizine Transfero-somes

A new drug delivery technology called transferosomes were came into existence which is an artificial vesicle designed to show the characteristics of a cell vesicle suitable for controlled and potentially directed drug delivery. Transferosome is a highly flexible and stress- responsive compound, complex compound and highly deformed vesicle with an aqueous core surrounded by the complex lipid bilayer thus enabling it to deliver both hydrophilic and hydrophobic drug. Urticaria is a general condition distinguished as brief erythematous and oedematous plaques or papules with defined erythematous borders and central clearing, identified as hives/wheals. Levocetirizine which is used for the treatment is a H1 anti histamine that is actively used in treatment of urticaria. The aim of the present research work was to investigate the potential of transferosome formulations for transdermal delivery of levocetirizine. The transferosomes were formulated by lipid film hydration technique using Rotary vacuum Evaporator. In the present work levocetirizine vesicle was efficaciously formulated using an appropriate ratio of tween 80 and soya lecithin and was incorporated into gel since gel formulations are easy to administer and patient compliance. Levocetirizine transferosomes was evaluated for vesicle characteristicslike zeta potential, poly-dispersibility index, TEM, and stability study. Theaverage sizes of transferosome were found to be 566.6nm and poly dispersity index (PI) was found to be 0.532, zeta potential of the transferosome was found to be -2.3 mV which indicates that transferosome formulation is stable. Levocetirizine gel was prepared by using various concentrations of Carbopol 934 and is evaluated for their gel characteristics like pH, viscosity, Spreadability, extrudability, homogeneity, drug content, diffusion etc. Gel containing 2% Carbopol show best and promising results. Levocetirizine transferosomal gel were efficaciously formulated by utilizing levocetirizine transferosome which was prepared by thin film hydration method by using soya lecithin and tween 80 in the ratio 85:15. Other Oral formulation used in treatment for urticaria have disadvantage like poor bioavailability, 1st pass metabolism, patient non-compliance etc. Transferosome nano gel loaded with levocetirizine transferosome was found to be more effective than other oral formulations used in treatment for urticaria since transferosome are capable of passing through lipid layer and delivering drug into systemic circulation with maximum bioavailability.

Molecular Simulations Reveal the Dynamics of the T-Cell Receptor in a T-Cell Model Membrane

The T-cell antigen receptor (TCR) initiates immune responses by recognising a wide variety of foreign peptides presented by Major Histocompatibility Complex (MHC) proteins. The TCR consists of non-covalently associated TCRαβ, CD3εγ, CD3εδ and ζζ dimers. The CD3 and ζζ subunits have long intracellular tails that regulate TCR signalling. Despite significant knowledge of the TCR structure and topology, lack of information on the arrangement of its cytoplasmic tails limits our understanding of the molecular mechanism of intracellular signal propagation in a T-cell. In this study, we used molecular modelling to model the entire TCR structure and performed multi-scale molecular dynamics simulations of the TCR embedded in a complex lipid bilayer which mimics the T-cell membrane. Our simulations revealed conformational changes within the transmembrane region of the resting TCR, and also showed that the intracellular signalling tails of CD3 and ζζ subunits fold into an intertwined conformation forming contacts with each other. Additionally, the tails strongly interacted with the inner leaflet of the membrane via cationic stretches, whilst intracellular tyrosine side-chains transiently penetrated into the hydrophobic region of the membrane. Phosphatidyl inositol phosphate (PIP) lipids and cholesterol formed preferential protein-lipid interactions creating a unique local lipid environment around the TCR. Our study provides insights into the TCR's complete structure and dynamics in a membrane in atomic resolution. Further, this potentially forms the basis to study the molecular mechanism of initiation of TCR signalling.

Targeting Extracellular Vesicles to Dendritic Cells and Macrophages.

Targeting protein therapeutics to specific cells and tissues is a major challenge in modern medicine. Improving the specificity of protein therapeutic delivery will significantly enhance efficiency in drug development. One of the promising tools for protein delivery is extracellular vesicles (EVs) that are enveloped by a complex lipid bilayer. EVs are secreted by almost all cell types and possess significant advantages: biocompatibility, stability, and the ability to penetrate the blood–brain barrier. Overexpression of the vesicular stomatitis virus protein G (VSV-G) was shown to promote EV formation by the producer cell. We have developed an EV-based system for targeted delivery of protein cargoes to antigen-presenting cells (APCs). In this study, we show that attachment of a recombinant llama nanobody α-CD206 to the N-terminus of a truncated VSV-G increases the selectivity of EV cargo delivery mainly to APCs. These results highlight the outstanding technological and biomedical potential of EV-based delivery systems for correcting the immune response in patients with autoimmune, viral, and oncological diseases.

Molecular Simulations of Intact Anion Exchanger 1 Reveal Specific Domain and Lipid Interactions

Anion exchanger 1 (AE1) is responsible for the exchange of bicarbonate and chloride across the erythrocyte plasma membrane. Human AE1 consists of a cytoplasmic and a membrane domain joined by a 33-residue flexible linker. Crystal structures of the individual domains have been determined, but the intact AE1 structure remains elusive. In this study, we use molecular dynamics simulations and modeling to build intact AE1 structures in a complex lipid bilayer that resembles the native erythrocyte plasma membrane. AE1 models were evaluated using available experimental data to provide an atomistic view of the interaction and dynamics of the cytoplasmic domain, the membrane domain, and the connecting linker in a complete model of AE1 in a lipid bilayer. Anionic lipids were found to interact strongly with AE1 at specific amino acid residues that are linked to diseases and blood group antigens. Cholesterol was found in the dimeric interface of AE1, suggesting that it may regulate subunit interactions and anion transport.

Interplay of Curvature, Lipid Segregation and Stability Modulation in Complex Lipid Bilayers

Interplay of Curvature, Lipid Segregation and Stability Modulation in Complex Lipid Bilayers

Molecular diffusion and nano-mechanical properties of multi-phase supported lipid bilayers.

Understanding the properties of cell membranes is important in the fields of fundamental and applied biology. While the characterization of simplified biological membrane mimics comprising liquid phase lipids has been routinely performed due to the ease of fabrication, the characterization of more realistic membrane mimics comprising multi-phase lipids remains challenging due to more complicated fabrication requirements. Herein, we report a convenient approach to fabricate and characterize multi-phase supported lipid bilayers (SLBs). We employed the solvent-assisted lipid bilayer (SALB) formation method to fabricate mixed lipid bilayers comprising liquid phase 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and gel phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipids at room temperature. The fabrication procedure was performed inside a newly designed microfluidic chamber, which facilitated the subsequent characterization of the SLBs without exposure to air. The SLBs were then characterized via fluorescence microscopy, fluorescence recovery after photobleaching (FRAP), atomic force microscopy (AFM) and AFM-based force-distance measurements. Interestingly, results from these characterization techniques revealed that regardless of the gel phase composition, the SALB formation method consistently yielded uniform SLBs at room temperature, even though the transition temperature of DPPC is considerably higher. Furthermore, the composition ratio of DOPC and DPPC in the precursor solution is well reproduced in the fabricated SLBs. We also identified from diffusivity measurements that a high ratio of gel phase lipid revitalizes lipid-lipid interactions, which led to reduced molecular fluidity and the suppression of thermal undulation within the SLBs. Taken together, our results highlight the robustness of the SALB formation method that allows the fabrication of complex lipid bilayers with a high degree of precision, which is suitable for functional studies of biological membranes.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name “cholesterolic gel” domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

Label-free and charge-sensitive dynamic imaging of lipid membrane hydration on millisecond time scales

Biological membranes are highly dynamic and complex lipid bilayers, responsible for the fate of living cells. To achieve this function, the hydrating environment is crucial. However, membrane imaging typically neglects water, focusing on the insertion of probes, resonant responses of lipids, or the hydrophobic core. Owing to a recent improvement of second-harmonic (SH) imaging throughput by three orders of magnitude, we show here that we can use SH microscopy to follow membrane hydration of freestanding lipid bilayers on millisecond time scales. Instead of using the UV/VIS resonant response of specific membrane-inserted fluorophores to record static SH images over time scales of >1,000 s, we SH imaged symmetric and asymmetric lipid membranes, while varying the ionic strength and pH of the adjacent solutions. We show that the nonresonant SH response of water molecules aligned by charge-dipole interactions with charged lipids can be used as a label-free probe of membrane structure and dynamics. Lipid domain diffusion is imaged label-free by means of the hydration of charged domains. The orientational ordering of water is used to construct electrostatic membrane potential maps. The average membrane potential depends quadratically on an applied external bias, which is modeled by nonlinear optical theory. Spatiotemporal fluctuations on the order of 100-mV changes in the membrane potential are seen. These changes imply that membranes are very dynamic, not only in their structure but also in their membrane potential landscape. This may have important consequences for membrane function, mechanical stability, and protein/pore distributions.

Formation of nanopores in DiynePC-DPPC complex lipid bilayers triggered by on-demand photo-polymerization.

Vesicles have unique characteristics that enable the release of drugs as well as encapsulation while maintaining biocompatibility. A photo-polymerizable liposome composed of 1,2-bis(10,12-tricosadiynoyl)-sn-glycero-3-phosphocholine (23:2 DiynePC) has been investigated as vehicles for triggered delivery of drugs to cells. In this study, we confirmed for the first time that supported lipid bilayers (SLBs) prepared with a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC)/DiynePC mixture generated pores ca. 100–300 nm in size on the membrane after UV polymerization. This direct observation was done by analyzing the SLBs formed with the DPPC/DiynePC mixture by employing atomic force microscopy (AFM) in a liquid environment. However, photo-polymerization did not occur in the 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/DiynePC mixed bilayer and pores were not formed. A theoretical study was performed to explore the phase behavior of the lipid mixtures. A coarse-grained model of DiynePC was developed that is comparable with the Martini force field; the parameters were validated against atomistic simulations. Transition from fluidic to gel phase was observed only when DiynePC was mixed with DPPC, whereas the DOPC mixture remained fluidic over the entire domain. This implies a correlation between the formation of DiynePC-rich gel phase domains and the generation of pores after polymerization. The size of the pores were found to be controlled by the amount of polymerizable lipid which results in higher release rate of encapsulated calcein from the vesicles with larger pores.

Multiscale molecular dynamics simulations of lipid interactions with P-glycoprotein in a complex membrane

P-glycoprotein (P-gp) can transport a wide range of very different hydrophobic organic molecules across the membrane. Its ability to extrude molecules from the cell creates delivery problems for drugs that target proteins in the central nervous system (CNS) and also causes drug-resistance in many forms of cancer. Whether a drug will be susceptible to export by P-gp is difficult to predict and currently this is usually assessed with empirical and/or animal models. Thus, there is a need to better understand how P-gp works at the molecular level in order to fulfil the 3Rs: Refinement, reduction and replacement of animals in research. As structural information increasingly becomes available, our understanding at the molecular level improves. Proteins like P-gp are however very dynamic entities and thus one of the most appropriate ways to study them is with molecular dynamics simulations, especially as this can capture the influence of the surrounding environment. Recent parameterization developments have meant that it is now possible to simulate lipid bilayers that more closely resemble in vivo membranes in terms of their composition. In this report we construct a complex lipid bilayer that mimics the composition of brain epithelial cells and examine the interactions of it with P-gp. We find that the negatively charged phosphatidylserine lipids in the inner leaflet of the membrane tend to form an annulus around P-gp. We also observed the interaction of cholesterol with three distinct areas of the P-gp. Potential of mean force (PMF) calculations suggest that a crevice between transmembrane helices 10 and 12 has particularly favourable interaction energy for cholesterol.

Publisher’s note

P-glycoprotein (P-gp) can transport a wide range of very different hydrophobic organic molecules across the membrane. Its ability to extrude molecules from the cell creates delivery problems for drugs that target proteins in the central nervous system (CNS) and also causes drug-resistance in many forms of cancer. Whether a drug will be susceptible to export by P-gp is difficult to predict and currently this is usually assessed with empirical and/or animal models. Thus, there is a need to better understand how P-gp works at the molecular level in order to fulfil the 3Rs: Refinement, reduction and replacement of animals in research. As structural information increasingly becomes available, our understanding at the molecular level improves. Proteins like P-gp are however very dynamic entities and thus one of the most appropriate ways to study them is with molecular dynamics simulations, especially as this can capture the influence of the surrounding environment. Recent parameterization developments have meant that it is now possible to simulate lipid bilayers that more closely resemble in vivo membranes in terms of their composition. In this report we construct a complex lipid bilayer that mimics the composition of brain epithelial cells and examine the interactions of it with P-gp. We find that the negatively charged phosphatidylserine lipids in the inner leaflet of the membrane tend to form an annulus around P-gp. We also observed the interaction of cholesterol with three distinct areas of the P-gp. Potential of mean force (PMF) calculations suggest that a crevice between transmembrane helices 10 and 12 has particularly favourable interaction energy for cholesterol.