A: All atom molecular simulation systems

Similar documents
Interactions of Polyethylenimines with Zwitterionic and. Anionic Lipid Membranes

Supporting Information

8 Influence of permeation modulators on the behaviour of a SC lipid model mixture

Effects of Epicholesterol on the Phosphatidylcholine Bilayer: A Molecular Simulation Study

Structure of Dipalmitoylphosphatidylcholine/Cholesterol Bilayer at Low and High Cholesterol Concentrations: Molecular Dynamics Simulation

AFM In Liquid: A High Sensitivity Study On Biological Membranes

Under the Influence of Alcohol: The Effect of Ethanol and Methanol on Lipid Bilayers

Gold nanocrystals at DPPC bilayer. Bo Song, Huajun Yuan, Cynthia J. Jameson, Sohail Murad

MOLECULAR DYNAMICS SIMULATION OF MIXED LIPID BILAYER WITH DPPC AND MPPC: EFFECT OF CONFIGURATIONS IN GEL-PHASE

Phase Transition Behaviours of the Supported DPPC Bilayer. Investigated by Sum Frequency Generation (SFG) and Atomic Force

Visualizing Biopolymers and Their Building Blocks

Fluid Mozaic Model of Membranes

Methods and Materials

Emulsification. An Introduction to the Emulsification of Lipids within a Watery Environment. By Noel Ways

Penetration of Gold Nanoparticle through Human Skin: Unraveling Its Mechanisms at the Molecular Scale

The Interaction between Lipid Bilayers and Biological Membranes. Chapter 18

Supplementary Figures

Cholesterol-Sphingomyelin Interactions: A Molecular Dynamics Simulation Study

Water Isotope Effect on the Phosphatidylcholine Bilayer Properties: A Molecular Dynamics Simulation Study

Supplementary Information: A Critical. Comparison of Biomembrane Force Fields: Structure and Dynamics of Model DMPC, POPC, and POPE Bilayers

Biochemistry 1 Recitation1 Cell & Water

arxiv:cond-mat/ v1 [cond-mat.soft] 14 Aug 2002

Supplementary information for Effects of Stretching Speed on. Mechanical Rupture of Phospholipid/Cholesterol Bilayers: Molecular

2. Which is likely to be a nonpolar solvent? A. B. B shows a carboxyl group, while A is only carbons and hydrogens.

The main biological functions of the many varied types of lipids include: energy storage protection insulation regulation of physiological processes

Life Sciences 1a. Practice Problems 4

Neutron reflectivity in biology and medicine. Jayne Lawrence

Biochemistry. 2. Besides carbon, name 3 other elements that make up most organic compounds.

Bilayer Deformation, Pores & Micellation Induced by Oxidized Lipids

Biology Chapter 2 Review

Detergent solubilised 5 TMD binds pregnanolone at the Q245 neurosteroid potentiation site.

Coarse-grained model for phospholipid/cholesterol bilayer employing inverse Monte Carlo with thermodynamic constraints

Lessons of Slicing Membranes: Interplay of Packing, Free Area, and Lateral Diffusion in Phospholipid/Cholesterol Bilayers

Interaction of Functionalized C 60 Nanoparticles with Lipid Membranes

Chapter 2 Part 3: Organic and Inorganic Compounds

SUPPLEMENTARY INFORMATION

Andrey A. Gurtovenko*, and Ilpo Vattulainen,,

Physiologically-relevant levels of sphingomyelin, but not GM1, induces a beta-sheet-rich structure in the amyloid-beta(1-42) monomer

Lipids and Membranes

Supplementary Figures

Supporting Information Identification of Amino Acids with Sensitive Nanoporous MoS 2 : Towards Machine Learning-Based Prediction

Coarse grained simulations of Lipid Bilayer Membranes

Danish Research Institute of Translational Neuroscience DANDRITE, Nordic-EMBL Partnership

Biological Membranes. Lipid Membranes. Bilayer Permeability. Common Features of Biological Membranes. A highly selective permeability barrier

Biochemistry Worksheet

Shape Transformations of Lipid Bilayers Following Rapid Cholesterol

Phospholipid Component Volumes: Determination and Application to Bilayer Structure Calculations

BIOLOGICALLY IMPORTANT MOLECULES

FIGURE A. The phosphate end of the molecule is polar (charged) and hydrophilic (attracted to water).

Simulationen von Lipidmembranen

Probing protein heterogeneity in the plasma membrane using PALM and pair correlation analysis

Will s Pre-Test. (4) A collection of cells that work together to perform a function is termed a(n): a) Organelle b) Organ c) Cell d) Tissue e) Prison

and controllable behavior - Supplementary Information

Molecular modeling of the pathways of vesicle membrane interaction. Tongtao Yue and Xianren Zhang

The Transport and Organization of Cholesterol in Planar Solid-Supported Lipid Bilayer Depend on the Phospholipid Flip-Flop Rate

Lipids are used to store and excess energy from extra carbohydrates in animals

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Interaction Between Amyloid-b (1 42) Peptide and Phospholipid Bilayers: A Molecular Dynamics Study

Biological Molecules B Lipids, Proteins and Enzymes. Triglycerides. Glycerol

Physical Cell Biology Lecture 10: membranes elasticity and geometry. Hydrophobicity as an entropic effect

List of Figures. List of Tables

A student carried out a series of food tests on a sample of milk. The student s observations and conclusions are shown in Table 1.1.

Impact of cholesterol on voids in phospholipid membranes

MARTINI Coarse-Grained Model of Triton TX-100 in Pure DPPC. Monolayer and Bilayer Interfaces. Supporting Information

Improved Coarse-Grained Modeling of Cholesterol-Containing Lipid Bilayers

Simulationen von Lipidmembranen

Chapter 7: Membranes

Good Afternoon! 11/30/18

Nanoparticle Permeation Induces Water Penetration, Ion Transport, and Lipid Flip-Flop

Ch. 5 The S & F of Macromolecules. They may be extremely small but they are still macro.

Irina V. Ionova, Vsevolod A. Livshits, and Derek Marsh

0.5 nm nm acyl tail region (hydrophobic) 1.5 nm. Hydrophobic repulsion organizes amphiphilic molecules: These scales are 5 10xk B T:

2.1.1 Biological Molecules

Condensed Complexes and the Calorimetry of Cholesterol-Phospholipid Bilayers

Chapter Three (Biochemistry)

Introduction to Macromolecules. If you were to look at the nutrition label of whole milk, what main items stick out?

Details of Organic Chem! Date. Carbon & The Molecular Diversity of Life & The Structure & Function of Macromolecules

WHAT IS A LIPID? OBJECTIVE The objective of this worksheet is to understand the structure and function of lipids

Membrane Structure and Membrane Transport of Small Molecules. Assist. Prof. Pinar Tulay Faculty of Medicine

Chapter 2 Transport Systems

The Structure and Function of Biomolecules

CHAPTER 4 EFFECT OF OXALIC ACID ON THE OPTICAL, THERMAL, DIELECTRIC AND MECHANICAL BEHAVIOUR OF ADP CRYSTALS

Lipids and Membranes

Asymmetry of lipid bilayers induced by monovalent salt: Atomistic molecular-dynamics study

Supplementary Figure 1 (previous page). EM analysis of full-length GCGR. (a) Exemplary tilt pair images of the GCGR mab23 complex acquired for Random

Carbohydrates, Lipids, Proteins, and Nucleic Acids

Molecular dynamics simulation studies of lipid bilayer systems

3.9 Carbohydrates. Provide building materials and energy storage. Are molecules that contain carbon, hydrogen and oxygen in a 1:2:1 ratio

Supplementary Materials for

Supporting material. Membrane permeation induced by aggregates of human islet amyloid polypeptides

ALLOSTERIC REGULATION OF GPCR ACTIVITY BY PHOSPHOLIPIDS

Page 2. The pie charts show the results. What conclusions can you make from the data? [Extra space]

Membrane transport. Pharmacy Dr. Szilvia Barkó

Polyoxometalate Macroion Induced Phase and Morphology

WHY IS THIS IMPORTANT?

Topic 3: Molecular Biology

Membrane Structure and Function - 1

Chapter 3. Table of Contents. Section 1 Carbon Compounds. Section 2 Molecules of Life. Biochemistry

Macromolecules Carbohydrates A COMPLEX COLORING EXPERIENCE

Transcription:

Cholesterol level affects surface charge of lipid membranes in physiological environment Aniket Magarkar a, Vivek Dhawan b, Paraskevi Kallinteri a, Tapani Viitala c, Mohammed Elmowafy c, Tomasz Róg d, Alex Bunker a a. Centre for Drug Research, Faculty of Pharmacy, University of Helsinki, Helsinki, Finland; b. Bombay College of Pharmacy, University of Mumbai, Mumbai, India; c. Division of Biopharmaceutics and Pharmacokinetics, Faculty of Pharmacy, University of Helsinki, Helsinki, Finland; d. Department of Physics, Tampere University of Technology, Tampere, Finland; Supplementary Information A: All atom molecular simulation systems Cholesterol Content % 16.66% 2% 25% 33.33% 5% DSPC/POPC 288 24 23 216 192 144 Cholesterol 48 58 72 96 144 Water 126 126 126 126 126 126 Na Ions 3 3 3 3 3 3 Cl Ions 3 3 3 3 3 3 Table S1: Composition of the lipid bilayers used in the molecular dynamics simulations The table shows the simulation system s molecular content used in this study. B: Surface area per lipid calculations.7 Area per lipid (nm2).65.6.55.5.45.4.35 % 1% 2% 3% 4% 5% Cholesterol content in lipid bilayer POPC bilayer DSPC bilayer Figure S2: Area per lipid: The figure shows that as the cholesterol content increases the area per lipid decreases that is the membrane bilayer contracts. This is in accordance with the previous studies. (The area per lipid was calculated with the following formula: Area per lipid = (Total area of membrane bilayer (.39* number of cholesterol molecules))/ number of PC molecules.

!"#$%% #&#$%% Figure S3: The figure shows area per lipid as a function of time of DSPC/POPC lipid bilayer containing cholesterol. The trajectories post 1 ns (marked by dashed line) were considered as equilibrated and thus used for the analysis.

C: Partial mass density profile of ions in simulation systems Density (kg m -3 ) 1 1 Phospate Na - % Cholesterol Cl - % Cholesterol Na - 16% Cholesterol Cl - 16% Cholesterol Na - 2% Cholesterol Cl - 2% Cholesterol Na - 25% Cholesterol Cl - 25% Cholesterol Na - 33% Cholesterol Cl - 33% Cholesterol Na - 5% Cholesterol Cl - 5% Cholesterol 1 1-1 1 2-1 1 2 Box (nm) Box (nm)!"#$!%#$ Figure S4: Detailed mass density profile of the Phosphate head groups of DSPC (A) and POPC (B) with Na + and Cl -. In the absence of cholesterol, the Na + are localized with membrane head groups. But as the cholesterol content is increased, association of Na + ions with membrane head groups decreases. D: Na + interactions with phospholipid bilayers )!" (!"!"#$%% *"+, - "./12"3/"454.6,15".78,956" '!" &!" %!" $!" #!" +, - "./12"3/"454.6,15".78,956"!"!" #!" $!" %!" &!" '!!" '#!" '$!" '%!" '&!" #!!" *!" :745";1<=" )!" (!" '!" &!" %!" $!" #!" #&#$%%!"!" $!" &!" (!" *!" #!!" #$!" #&!" #(!" #*!" $!!" :745";1<=" Figure S5: % Bound Na + to membrane headgroups vs. time. The figure represents the % amount of Na + ions bound to DSPC and POPC membrane containing cholesterol. The trajectories post 1 ns (marked by dashed line) were considered as equilibrated and thus used for the analysis.

Number of Na+ bound 16 14 12 1 8 6 4 2 % 16.66% 2% 25% 33.33% 5% Cholesterol content in DSPC bilayer Number of Na+ bound 8 7 6 5 4 3 2 1 % 16.66% 2% 25% 33.% 5% Cholesterol content in POPC bilayer Figure S6: Identity of oxygens to which Na + ions are bound. For all systems we see that binding (see text for definition of binding) to the phosphate oxygens dominates

456789":;"4< =" #!" '&" '%" '$" '#" '!" &" %" $" #"!" ()*+,()*+" ()*+,+-./" ()*+,()*+" ()*+,+-./" ()*+,()*+" ()*+,+-./" ()*+,()*+" ()*+,+-./" ()*+,()*+" ()*+,+-./" ()*+,()*+" ()*+,+-./"!" '%1%%" #!" #2" 33133" 2!" +>:?8@A89:?"+:BA8BA"CB"()*+"7C?<D89"!"#$%&#'()*+#$',-'+,)(.'/%'1,(2'' (#'! (#&! (#%! (#$! (! "#'! "#&! "#%! "#$! "! "! (&#&&! $"! $1! //#//! 1"! 34,5#26#$,5'7,(6#(6'1('8983'+15%:#$' @! 456789":;"4<="?! >! =! <! ;! :! *.*+,+-./" *.*+,+-./" *.*+,+-./" *.*+,+-./" *.*+,+-./" *.*+,+-./"!"#$%&#'()*+#$',-'/%'+,)(.'1,(2' %#1! %! /#1! /! $#1! $! (#1! (! "#1! "! 2"#1!!" '%1%%" #!" #2" 33" 2!" +>:?8@A89:?"+:BA8BA"CB"*.*+"7C?<D89" "! (&#&&! $"! $1! //#//! 1"! 34,5#26#$,5'7,(6#(6'1(';<83'+15%:#$'! Figure S7 Plots showing stoichiometry of number of sodium ions that bridge different headgroup binding sites.

1 Pure system 16.66% Chol. 2% Chol. 25% Chol. 33% Chol. 5 % Chol. 1 ACF ACF.1!"!#$.1 %&!#$ 5 1 15 2 Time (ns)!"#$ %& %' %( %) %* %+ %, %- %. %/ 123 &4.- (4'56) &4'& &4&-. &4(, '4+* &4&&() ))4+ ',78$9:;4 &4,/ '4-56) &4() &4&*+ &4&), &4(- &4&*' '4(/ &4&&-* '/4( (&78$9:;4 &4)) '4&56+ &4+& &4&'. &4'( &4(. &4&*' '4+* &4&&/- '-4, (+78$9:;4 &4). (4,56) &4). &4&)/ &4'+ &4*- &4&-( (4+( &4&&/. (/4+ ))78$9:;4 &4(+ )4+56) &4*, &4&.* &4( &4,, &4&-- )4&/ &4&'(' )+4+ +&78$9:;4 &4'- '4-56) &4)+ &4&// &4)& &4+/ &4'+ (4,- &4&)) '.4' #<#$ %& %' %( %) %* %+ %, %- %. %/ 123 &4)( '4-56, &4() &4&(, &4)+ &4'( &4&.. '4,' &4&&.& )-4) ',78$9:;4 &4(. '4,56* &4*. &4&*/ &4'- &4+' &4&,( )4*- &4&&-( /'4+ (&78$9:;4 &4)) '4.56* &4() &4&)+ &4)+ &4'. &4&.. (4+- &4&&.& +/4( (+78$9:;4 &4(& *4(56* &4)+ &4&** &4)& &4', &4'* '4-. &4&') )/4, ))78$9:;4 &4(..4'56) &4'- &4&&(- &4*' &4'. &4'( '4.) &4&'+ )*4/ +&78$9:;4 &4'& (4&56* &4), &4&), &4(- &4'. &4() &4., &4&), '&4-+ Figure S8: Auto-correlation functions (ACF) of Na+ binding with the interface in the POPC (upper panel) and DSPC (lower panel) bilayers. The autocorrelation function was defined as: ACF(t )= < hw(t +t ) hw(t)> / < hw(t)>, where the function hw(t) is 1 if the

Na + ion is bonded with lipids, and otherwise (J. Phys. Chem. B 114, 11784 11792). The ACF curves show multiexponential fits in all cases, for the pure DSPC system 4 exponentials, and for the POPC systems 5 exponentials. The fit parameters are included in the above table, where the parameters correspond to a fit of the decay to the function y = a e a1x + a 2 e a3x + a 4 e a 5 x + a 6 e a 7 x + a 8 e a 9 x. In all cases the curves show clear multiexponential fits, for the pure DSPC system four exponentials, and for the POPC systems five exponentials, thus four decay modes for all systems and an extra mode is absent in the gel system. The fit parameters of the multi exponential fit are included. While interpreting this data one can discard the fastest mode that is both small amplitude, and decays over the first few data points, thus is neither well sampled nor important. This fast time scale can be associated with fast local oscillations. The slowest decaying exponential is likely the result of finite system size effects. This leaves us with two decay modes for the gel phase DSPC system three decay modes for all other systems (all DSPC-Chol POPC). Unambiguous mechanistic interpretation of which specific interactions the different decay rates correspond to is not possible. As cholesterol content increases the ion mobility increases and the mechanism for the increase can be seen to be increased strength of the faster decay modes, while the rates of the different decay modes seem to remain approximately constant. The characteristic binding time can be approximated as the time taken for the case of all POPC bilayers the initial decay in the ACF is completed in ~4 ns time, indicating full exchange of ions between interface and water. For the case of DSPC only in the two highest cholesterol concentration curves is the constant value reached. These results indicate that cholesterol increases the rate of ion exchange between water and membrane, and our results for the pure DSPC system are in agreement with our previously published results (Stepniewski et al., J. Chem. Phys. B 21).

Number of bound Na+ ions/pc group.18.17.16.15.14.13.12.11.1 % 1% 2% 3% 4% 5% Cholesterol content in lipid bilayer DSPC POPC Figure S9: Number of bound Ions per PC headgroup present. For both DSPC and POPC these decline with increasing cholesterol content.

E: Cl - interactions with phospholipid bilayers % Ions bound 1 9 8 7 6 5 4 3 2 1 16.66% 2% 25% 33.33% 5% Cholesterol content of DSPC membrane bilayer Bound to Na+ Bound to Membrane Headgroups Free % Ions bound 1 9 8 7 6 5 4 3 2 1 % 16.66% 2% 25% 33.33% 5% Cholesterol content of POPC membrane bilayer Bound to Na+ Bound to Membrane Headgroups Free Figure S1: Cl - ion binding, (see text for definition of binding) in the DSPC and POPC bilayer simulations. We see no significant binding of Cl - ions to lipid headgroups.

1!"#$%&'($)*+,-./$ % Cholesterol 16.66% Cholesteol 2% Cholesterol 25% Cholesterol 33.33% Cholesterol 5% Cholesterol RDF.5 1!#$'1'($)*+,-./$ 1 2 Distance (nm) % Cholesterol 16.66% Cholesterol 2% Cholesterol 25% Cholesterol 33.33% Cholesterol 5% Cholesterol RDF.5 1 2 Distance (nm) Figure S11: Radial distribution function of distance between Cl - ions and choline group of DSPC/POPC. The figure shows no sharp peak indicating there is no strong interaction between them.

F: interactions of cholesterol with phospholipids 7 6 5 4 DSPC POPC 3 2 1 1% 2% 3% 4% 5% Figure S12A: Number of hydrogen bonds between hydroxyl group of cholesterol and phospholipids. This shows as the cholesterol content in the lipid bilayer increases, the number of hydrogen bonds between the hydroxyl group of cholesterol and the phopholipid head group..5.4.3 DSPC.2 POPC.1 1% 2% 3% 4% 5% Figure S12B: Number of hydrogen bonds between hydroxyl group of cholesterol and phospholipids per DSPC/POPC molecules. This shows as the cholesterol content in the lipid bilayer increases, the likelihood that the DSPC/POPC lipid is bonded to choleterol increases.

Cos<P>.4.2 5% Cholesterol 33% Cholesterol 25% Cholesterol 2% Cholesterol 16% Cholesterol % Cholesterol -.2.1 Cos<P> -.1 -.2 1 2 3 Position relative to phosphate peak(nm) Figure S13 Water ordering plots normal to the membrane surface. We see in both the case of POPC (upper panel) and DSPC (lower panel) that the qualitative form of the water ordering is unchanged from previous computational results (J. Phys. Chem. B. 114, 11784-11792 (21)). For the case of DPPC the bilayer water ordering is systematically decreased with increasing cholesterol concentration while for the case of POPC for smaller cholesterol concentrations (- 25mol%) water ordering remains relatively constant and decreases for higher concentrations. Water ordering at the lipid/water interface originates from strong interactions between lipid headgroups and water, e.g. hydrogen bonds. The addition of cholesterol into the ordered DSPC bilayer increases the area of the whole membranes systematically with increasing cholesterol concentration (spacing effect) thus the area per individual headgroup increases and effectively it is seen as a decrease of the water ordering. For the case of POPC in lower concentrations the total area does not increase as a result of the condensing cholesterol effect which rich saturation at 2-3mol% (condensing effect does not occur in DSPC as lipids in the gel phase are already condensed), after this point the total area of the membrane becomes larger which again manifests itself in a decreased degree of water ordering. Thus the whole effect might be attributed to cholesterol spacing vs condensing effect.