Phase Behavior of GM1-Containing DMPC ... - ACS Publications

Sep 3, 2018 - (4) Shevchenko, A.; Simons, K. Lipidomics: coming to grips with lipid diversity. Nat. Rev. Mol. Cell Biol. .... Cell Death Differ. 2006,...
0 downloads 0 Views 1MB Size
Subscriber access provided by EKU Libraries

Biological and Environmental Phenomena at the Interface

Phase Behavior of GM1 containing DMPC-Cholesterol Monolayer: Experimental and Theoretical Study Zarrin Shahzadi, Subhasis Das, Tanushree Bala, and Chaitali Mukhopadhyay Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.8b02621 • Publication Date (Web): 03 Sep 2018 Downloaded from http://pubs.acs.org on September 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Phase Behavior of GM1 containing DMPC-Cholesterol Monolayer: Experimental and Theoretical Study Zarrin Shahzadi, Subhasis Das, Tanushree Bala, and Chaitali Mukhopadhyay* Department of Chemistry, University of Calcutta, 92, A. P. C. Road, Kolkata – 700009, India * Address for Correspondence: Chaitali Mukhopadhyay Department of Chemistry, University of Calcutta 92, A. P. C. Road, Kolkata – 700 009, India. Email: [email protected], [email protected]

1 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Organization and distribution of lipids in cellular membranes play an important role in a diverse range of biological processes, such as membrane trafficking and signaling. Here we present the combined experimental and simulated results to elucidate the phase behavioral features of Ganglioside Monosialo 1 (GM1) containing mixed monolayer of the lipids 1, 2-Dimyristoyl-snglycero-3-phosphocholine (DMPC) and Cholesterol (CHOL). Two monolayers having compositions DMPC/CHOL and GM1/DMPC/CHOL are investigated at air-water and air-solid interfaces using Langmuir Blodgett experiments and Scanning Electron Microscope (SEM) respectively to ascertain the phase behavior change of the monolayers. Surface pressure isotherms and SEM imaging of domain formation indicate that addition of GM1 to the monolayer at low surface pressure causes a fluidization of the system but once the system attains the surface pressure corresponding to its liquid-condensed phase, the monolayer becomes more ordered than the system devoid of GM1 and interacts among each other more co-operatively. Besides, condensing effect of Cholesterol on the DMPC monolayer was also verified by our experiments. Apart from these, the effects induced by GM1 on the phase behavior of the binary mixture of DMPC/CHOL were studied with and without applying liquid expanded (LE)-liquid condensed (LC) equilibrium surface pressure using molecular dynamics simulation. Our Molecular dynamics (MD) simulation results give an atomistic level explanation of our experimental findings and furnish a similar conclusion. Introduction Lipids being the major components of the biological membrane have always been under great scrutiny.1-4 Lateral organization of lipids in cellular membranes and their interaction among each other play an important role in a diverse range of processes, such as membrane trafficking and

2 ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

signaling.5 The composition of lipid in the cellular membrane has a great impact on the features and function of the membrane.6 In order to understand the biological processes such as subcellular transport and invasion by infectious agent, a sound understanding of the physical and chemical properties of the membrane lipids is essential. Naturally occurring and mimic membranes have been studied extensively in the past and many physicochemical characteristics have been extrapolated by various physical methods, including differential scanning calorimetry7,

diffraction8,9,

X-ray

neutron

diffraction10,

optical

microscopy11,

IR

spectroscopy12and atomic force microscopy.13 Our current motive was to understand the principles behind phase behavior and dynamics of GM1 containing membrane by determining how these individual lipid molecules interact with each other before and after the liquidexpanded(LE)/liquid-condensed(LC) equilibrium. In order to mimic biological membrane for our study, we chose 1, 2-Dimyristoyl-sn-glycero-3-phosphocholine (DMPC), Cholesterol (CHOL) and Ganglioside Monosialo 1 (GM1) as membrane components.

O

O

CH3 P H3C

N+

O

O O-

O H

DMPC O

CH3

O

OH

H

CHOL

H H

3 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 33

OH CH2OH O OH O

OH O

OH CH2OH

CH3

N H

O

O O

H3C

O CH2OH

HO N H

O

O

OH

O

H

O HO OH

-

OH

OOC

NH

GM1

O

HO O CH2OH

H

OH

OH

Figure 1: Chemical structures of the lipids used in this work. (Prepared using ChemDraw Ultra 8.0)

DMPC was chosen for our experiments because over the years, DMPC membranes have been rigorously studied using a variety of physical techniques14,15 and as a result, they’ve been used to study problems regarding different biological systems like drug-membrane interactions16, cholesterol flip-flop17, membrane properties18 etc. DMPC contains phopsphocholine group and saturated lipid tail; phosphocholine is one of the major components of mammalian membranes and because of having fully saturated acyl chains DMPC monlayers and bilayers are robust and do not require special preparatory environments.19 Besides, from biophysical point of view it is 4 ACS Paragon Plus Environment

Page 5 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

convenient to monitor phase transitions of saturated model lipids and for this fact DMPC has been widely used for Langmuir Blodgett experiments.20,21 Cholesterol influences individual molecular behavior thus regulating the overall mechanical behavior and eventually modulates biological membranes and their functions.22,23Glycolipids are basically lipid molecules which contains sugar groups, found in most animal cell plasma membranes, they regulate various physiological events at the cell surface.24,25 Gangliosides contain one or more negatively charged sialic acid groups and they are the most complex form of glycolipids.26 Gangliosides play a major role in a number of cellular functions, including cell recognition and adhesion, signal transduction and cell growth regulation.27 They constitute 5–10% of the total lipid mass in nerve cells. The external surfaces of certain cells contain10–20mol% gangliosides because they exist mostly on the outer leaflet of the cell membrane.28 One of the most commonly studied gangliosides is GM1, a member of the glycosphingolipid family that contains four neutral sugar groups and one sialic acid residue.29 The organization of gangliosides in model membranes has been broadly investigated with various methods, such as AFM,30 secondary ion mass spectrometry, 31 fluorescence, 32 NMR.33 Lee et al. in their work explain the effect of ganglioside on the phase behavior of dipalmitoylphosphatidlycholine (DPPC) lipid. With varying concentration of GM1 they found that at a lower concentration of GM1, it has a condensing effect on DPPC monolayer, whereas, at higher concentration, it has a fluidizing effect on the monolayer.28 Yuan and Johnston reported about the GM1 rich domain formation in DPPC/CHOL monolayer.34 It is identified that cholesterol promotes the formation of GM1 clusters.35 To begin with; we used Langmuir Blodgett isotherms to monitor the phase transitions in monolayer because pressure/area isotherm is an effective way to monitor the interfacial behavior of lipid molecules. The surface area isotherm provides information about phases of lipid molecules at the

5 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

surface along with the corresponding changes in surface pressure.36 The surface pressure thus obtained corresponding to the LE/LC equilibrium phase of monolayers was used for further study. Besides, there are only a handful of computer simulation investigations to support a large number of existing experimental studies on Langmuir Blodgett (LB) isotherms.37-39 We used Molecular dynamics simulation to study the detailed structure of the monolayer under investigation. The chosen model of the monolayer of interest was observed at atomistic level to complement the experimental findings. Our current report provides a more realistic depiction of phase behavioral change of model membrane containing DMPC and Cholesterol in presence and absence of GM1 by both, experimental and theoretical studies. Experimental Section Materials DMPC and cholesterol were purchased from Sigma Chemical, USA and used without further purification. DEAE–Sephadex was obtained from Pharmacia. Ammonium acetate, KCl, KI, and Na2S2O3 were obtained from Merck and IR-120 cation exchange resin was from Mallinckrodt Chemical Works, New York. All the solvents were of the spectroscopic grade. Each experiment was performed at least in triplicate at 37°C. Methods Ganglioside Monosialo 1(GM1) extraction Ganglioside (GM1) was isolated and purified from goat brain following our published protocol.26 Briefly, brain tissue was homogenized and gangliosides mixtures were obtained from it. After extensive dialysis of the eluted mixture of gangliosides, it was lyophilized. The neuraminidase digestion removed sialic acid from polysialo-gangliosides and in the end gave the GM1. The

6 ACS Paragon Plus Environment

Page 6 of 33

Page 7 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

final solution was again lyophilized and GM1 was dissolved in water. The solution was passed through a cation exchange IR-120 resin column so that any cation present in the mixture gets removed. After extensive dialysis followed by lyophilization, white powdered GM1 was obtained. The purity of GM1 was checked by TLC and 1H NMR spectroscopy. (Figure S1 and S2 of supporting information) Langmuir Blodgett Isotherms DMPC and cholesterol (10 mg/ml) were dissolved in chloroform. GM1 (4 mg/ml) was dissolved in chloroform/methanol (v/v, 80:20). DMPC/CHOL mixture was prepared in 2:1 molar ratio and GM1/DMPC/CHOL mixture was prepared with 10 mol% GM1 in 2:1 DMPC/CHOL. The monolayer of DMPC, GM1, DMPC/CHOL and GM1/DMPC/CHOL mixtures were prepared on Langmuir-Blodgett trough (KSV NIMA medium trough) using deionized water as the subphase. The sample solutions were spread on the sub-phase surface; after solvent evaporation (30 mins), the surface pressure was measured with a precision of 0.1 mN/m using a Wilhelmy balance. The compression rate of isotherm was 25cm2/min. To protect the experimental setup from dust it was placed in a laminar flow hood in which temperature was kept constant at 37°C. Surface pressure was recorded simultaneously as a function of the molecular area. The accuracy of measurements was 0.02 Å2 area per molecule. In order to monitor the LE/LC equilibrium surface pressure of binary (DMPC-CHOL) and ternary (DMPC-CHOL-GM1) mixtures and get a clear picture of phase transition, another set of LB-experiments was performed with lower concentrations of lipids: DMPC and cholesterol (1 mg/ml), GM1 (0.4 mg/ml). DMPC/CHOL mixture was prepared in 2:1 molar ratio and GM1/DMPC/CHOL mixture was prepared with 10 mol% GM1 in 2:1 DMPC/CHOL as used in our preceding experiment. The compression rate of isotherm was also lowered to 1.0 cm2/min. (Figure S3 of supporting information)

7 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Scanning Electron Microscope (SEM) SEM measurements were performed with a Zeiss EVO 18, Special Edition microscope. Triplicate samples were prepared for each monolayer composition at 37°C, lipid monolayers from the Langmuir trough were transferred onto Si (111) wafers at definite surface pressure (2mN/m, 7mN/m) during the compression by vertically dipping the Si (111) wafers into the trough and slowly pulling it up (2 mm/min). The Si (111) wafers with samples deposited onto it, were then thoroughly dried in air before imaging40 Theoretical Section System Setup Monolayer Model. Two monolayers were modeled maintaining the same composition as used in our experiments i.e. (i) 2:1 DMPC/CHOL (40 DMPC and 20 Cholesterol molecules) (ii) In order to closely mimic the 10% GM1 in 2:1 DMPC/CHOL solution as used in our experiments, 5 GM1 molecules were added to system (i), but to avoid any steric clashes 5 DMPC molecules were removed. Therefore, the final ratio of lipid components in the monolayer was 1:7:4 GM1/DMPC/CHOL. In order to stabilize each of the monolayer systems, an additional monolayer was constructed and rotated in such a way that the lipid head groups were facing each other. 4071 TIP3P water models were then added between the gap of the two monolayers to correspond to the system as lipid tail−lipid head−water−lipid head−lipid tail.37-41This was followed by a short energy minimization to release the stress of the system. As the gel-to-fluid transition temperature of DMPC is 23.9°C

42

, hence in order to maintain bilayer fluidity, the system was equilibrated for

50 ns at the target temperature of 310 K. The equilibrated system was used as the starting configuration for each simulation. Both the systems were simulated with and without LE/LC 8 ACS Paragon Plus Environment

Page 8 of 33

Page 9 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

equilibrium surface pressure. The surface pressure Π is calculated as Π=γ0-γ, where γ0 is the surface tension of pure water and γ is the surface tension of the monolayer coated air-water interface.43 In order to attain the LE/LC equilibrium surface pressure obtained from our experimental data i.e. ∼4.5mN/m, surface tension of γ=45 mN/m was applied to both the system. Using γ0=49.5mN/m for TIP3P water at 310K 44, our target surface pressure was achieved and in order to attain zero surface pressure corresponding to the uplift region of experimental isotherm, surface tension of magnitude γ0 i.e. 49.5mN/m was applied to both the system. One thing to keep in mind is that despite the surface pressure of the system being zero there was interfacial surface tension acting between the lipid components which provided stability to the system. A 2fs time step was used to integrate the equation of motion and each simulation was run for 150 ns. Five Na+ ions were added to both systems to maintain electro neutrality. The system was prepared using CHARMM and was equilibrated using NAMD_2.10 package.45 with the standard CHARMM36 force field.46 Simulation protocol. The isobaric−isothermal (NPT) ensemble with imposed 3D periodic boundary conditions was used for all the MD simulations. The Langevin dynamics were used to maintain the temperature for all the simulations, while by using Nose−Hoover−Langevin piston the pressure was kept constant at 1 atm.47 Long-range electrostatic interactions were calculated using smooth Particle Mesh Ewald method calculations.48 Short-range interactions were computed with cut off at 10 Å. All bond lengths involving hydrogen atoms were constrained using the RATTLE49 and SETTLE50 algorithm. The trajectory analysis was performed with CHARMM 51(Chemistry at Harvard Macromolecular Mechanics). Results and Discussion Experimental

9 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In order to obtain the pressure-area isotherms of the monolayers, we performed LangmuirBlodgett experiments and the results are as follow: Monolayers at air/water interface Pure lipids Figure 2 shows the surface pressure-area per molecule isotherm of the DMPC, GM1, 2:1 DMPCCHOL and 10% GM1 in 2:1 DMPC-CHOL at the air/water interface. The GM1 forms a stable monolayer at the air/water interface. The lateral distribution of the GM1 in the monolayer varies between (100-160Å2) in literature reports and depends on the ganglioside source, the content of unsaturated hydrocarbon chains and purity.34,52-54 As seen in inset of Figure 2, the lift-off area of the GM1 molecule in the Langmuir monolayer recorded at the air/water interface at temperature 37°C was close to ∼140.6Å2, which was in agreement with previously published result 140Å2.53 The slight discrepancy in the values reported here could be attributed to different experimental conditions and source of materials used.54-56 Isotherm did not show any plateau as such. Hence there was not any liquid-expanded (LE)/liquid-condensed (LC) phase equilibrium, representing that GM1 remains in LE phase till its collapse that was at 41.4mN/m also supported by the literature.53 Moreover, there was a kink at ∼69Å2 which might be due to the increased interaction between the hydrophobic tails of GM1.This observation was also supported by previous finding. 28

10 ACS Paragon Plus Environment

Page 10 of 33

Page 11 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 2: Monolayer compression isotherms of 2:1 DMPC-CHOL and GM1 in 2:1 DMPCCHOL at 37°C. The inset curves are showing monolayer compression isotherm of pure DMPC and pure GM1 at 37°C. The isotherm for pure DMPC has uplift at ∼112.1 Å2and the isotherm collapsed at 36.4 mN/m. The lift-off area and the collapse pressure was in very good agreement with the literature results with reported values ∼115 Å2

57

and ∼37mN/m respectively.21,58 The isotherm of 2:1 DMPC-

CHOL mixture was shifted to a much lesser molecular area. Liftoff occurred at ∼73.4Å2 and it has a collapse pressure of 39.2mN/m, in agreement with published data for DMPC-CHOL20 and DPPC-CHOL.59 It signified a condensing effect induced by the cholesterol molecules to DMPC. The cholesterol molecule consists of higher hydrophobic structure and by the time it is added to the DMPC system, the molecules rearrange among their selves so that the contact between hydrophobic sterol ring and water gets minimized. This distribution causes energetically favored strong van der Waals cohesive interactions between sterol rings of cholesterol molecules and the

11 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

alkyl chains of DMPC molecules.59 One more noticeable observation was that the collapse surface pressure of 2:1 DMPC/CHOL was 39.2 mN/m which was slightly higher than the collapse surface pressure of pure DMPC i.e. 36.4mN/m, which implied that the introduction of cholesterol molecules to DMPC monolayer increased the stability of the monolayer. Addition of GM1 to 2:1 DMPC-CHOL mixture could not change the nature of the isotherm of the lipid mixture as such but the uplift of the isotherm was shifted towards the right, supported by the literature.60 Addition of GM1 shifted the uplift to higher molecular area i.e. ∼78.5 Å2, while there was no change in the phase behavior. One notable observation in case of binary (DMPC-CHOL) and ternary mixture (DMPC-CHOL-GM1) was at lower surface pressure. At surface pressure∼4.5 mN/m we could see the two isotherms intersecting each other. At the beginning when the surface pressure was low (below 4.5 mN/m), addition of GM1 caused an increment in area per molecule of the system; hinting towards fluidizing effect of GM1 on the lipid monolayer but once the surface pressure was increased, the isotherm of ternary system (DMPC-CHOL-GM1) moved to the left of binary system (DMPC-CHOL) and exhibited lower area per molecule than the binary system. This suggested that after a certain surface pressure, GM1 might have a condensing effect on the lipid monolayer. Still it was too early to conclude anything and in order to further fortify our notion we performed other experiments. Moreover in order to monitor the LE/LC equilibrium plateau of the binary and ternary system, we performed LB-experiments with lower concentrations and lower compression rate (as mentioned in Methods section) and found the LE/LC equilibrium plateau exists at ∼4.5mN/m (Figure S3 in supporting information). SEM results

12 ACS Paragon Plus Environment

Page 12 of 33

Page 13 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

SEM images were taken at two different surface pressures concurrently to image the effect of GM1 on the phase behavior of the monolayer. In order to do this, the binary and ternary lipid monolayers were transferred to the silicon wafer at two different surface pressures, low and high. Low surface pressure (~2mN/m) was corresponding to the LE phase of the monolayer and relatively high surface pressure (∼7mN/m) corresponding to the beginning of LC phase of the monolayer. The surface pressure was obtained from the corresponding isotherm of that respective system. The first image taken was of 2:1 DMPC-CHOL at lower surface pressure (Figure 3a). The image obtained, provided us the picture of the beginning of phase separation and majorly exhibited LE domains along with few LC domains.

Figure3: SEM image of (a) 2:1 DMPC-CHOL monolayer at lower surface pressure.(∼2mN/m) (b) 2:1 monolayer DMPC-CHOL at higher surface pressure.(∼7mN/m) (c) GM1 added to 2:1 DMPC-CHOL monolayer at lower surface pressure.(∼2mN/m) (d) GM1 added to 2:1DMPCCHOL monolayer at higher surface pressure.(∼7mN/m)

13 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

On imaging the same mixture at the surface pressure corresponding to the beginning of the LC phase (Figure 3b), it showed condensed domains of diameter 9.7±0.01µm. A significant increase in the amount of the condensed domain with many of the smaller domains now having coalesced into a larger one was observed. Addition of GM1 to 2:1 DMPC-CHOL changed the morphology of monolayer significantly. At lower surface pressure, the image showed that GM1 preferentially remained in a more ordered phase and it distributed the whole lipid monolayer into small clusters, forming roughly circular smaller domains of average diameter 0.55±0.001 µm. On comparing this image (Figure 3c) with the image of 2:1 DMPC-CHOL without GM1 (Figure 3a), we could say that GM1 chopped the condensed domains and distributed the whole system evenly into smaller domains. On taking the image of the same system at relatively higher surface pressure, i.e. surface pressure corresponding to the beginning of the LC phase of the system, we found that system attained condensed domains of size of 20.7±0.01µm which was much larger than the condensed domains of monolayer without GM1. Up to this, it was reflected from the pressure-area isotherms and SEM images that at the beginning when the DMPC-CHOL monolayer remained mainly in the LE phase, addition ofGM1 to the system had a fluidizing effect on the monolayer but once the monolayer crossed the surface pressure corresponding to the LE/LC equilibrium, the system was more compact than the system without GM1 and GM1 caused the formation of condensed domains of larger dimensions. Simulation Results So far from our experimental observations, we found that GM1 exhibits contrary effects to the phase behavior of the monolayer in the absence and in presence of surface pressure. In order to

14 ACS Paragon Plus Environment

Page 14 of 33

Page 15 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

analyze this at the atomistic level, we performed MD simulation of the concerned monolayer with and without LE/LC equilibrium surface pressure applied to it and from the simulation trajectories we carried out the following analyses. Preferential partitioning of lipid components SEM images indicated that from an initial homogeneous distribution, upon achieving the LE/LC equilibrium, GM1 molecules promoted the formation of large condensed domains. To specify the extent of mixing of monolayer components in presence and absence of LE/LC equilibrium surface pressure, we studied the preferential partitioning of GM1 along with other components in the monolayer. The preferential partitioning of membrane components is defined as the relative number of contacts of a lipid species with each of the other constituents, corrected for the total number of lipids in the system61: PA=

  ∑

Where, PA is the preferential partitioning of membrane component A to another component B, CA the number of contacts of component A and nA the number of molecules of component A. The data were averaged over the last 10 ns of simulation trajectory.

15 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure4: Preferential partitioning matrix of the monolayer. (a) Without GM1 without any surface pressure (b) without GM1 with LE/LC equilibrium surface pressure (c) with GM1 without any surface pressure (d) with GM1with LE/LC equilibrium surface pressure. If a lipid has more than one contact with another lipid this interaction is only counted once. Two lipids are defined as being in contact if the distance between the molecules is less than 6 A˚.

From figure 4(a), it was evident that when there was no surface pressure DMPC has a strong preference for DMPC molecules, whereas cholesterol exhibited the preference for both DMPC and other cholesterol molecules. In presence of LE/LC equilibrium surface pressure, it was found that slowly the preferential partitioning of individual components among themselves was decreased whereas the same was increased for cross components indicating the mixing of monolayer components. From SEM image of DMPC-CHOL monolayer we have seen that after attaining the LC surface pressure, the monolayer forms condensed domains. The increase in preferential partitioning for cross components observed here could be attributed to this condensed domain formation at high surface pressure. Considering the GM1 containing 16 ACS Paragon Plus Environment

Page 16 of 33

Page 17 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

monolayers, there was an increase in preferential partitioning of cholesterol and GM1 for each other at LE/LC equilibrium surface pressure. This was in accordance with the previous reports where it was shown that GM1 has a high affinity for cholesterol molecules.62,63 and GM1 forms micro domains in presence of cholesterol.64 However, GM1 has high preference for other GM1 molecules to form self-aggregates was also reflected in this data.63,65,66Overall, we could also say that after applying the LE/LC equilibrium surface pressure, the mixing of lipid components was increased thus favoring the formation of condensed domains as imaged by SEM. Area per lipid Average area per lipid from the last 10 ns of the trajectory was calculated for lipid molecules and graphs are shown in figure 5. As shown in Table 1, in absence of any surface pressure, the average area per lipid attained by the GM1 containing monolayer at the end of the trajectory was higher than that for the system without GM1. The system with GM1 reached a value of area per lipid of 70.4±0.2 Å2, whereas the same in absence of GM1 was 66.1±0.2 Å2. This observation also supported our findings from experiments where it was found that GM1 has a fluidizing effect on monolayer before the monolayer attains surface pressure corresponding to its LE/LC equilibrium. While on applying the LE/LC equilibrium surface pressure to the monolayer, addition of GM1 increased the compactness of the monolayer. On applying LE/LC equilibrium surface pressure, the area per lipid of the monolayer without GM1 was found to be 65.0±0.2 Å2, whereas the same for the monolayer containing GM1 was 61.2±0.3 Å2. The SEM images showed that once the system attains LE/LC equilibrium, afterward the monolayers forms condensed domains, size of the condensed domain being larger in presence of GM1. Here also the GM1 containing

17 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

monolayerhaving LE/LC equilibrium surface pressure applied to it, exhibited the least area per lipid thus reflecting condensed state of the monolayer Table 1: Area per lipid (Å2) of different systems.

System

Area per lipid (Å2)

Without Surface pressure, without

66.1±0.2

GM1 With

LE/LC

equilibrium

surface

65.0±0.2

pressure, without GM1 Without Surface pressure, with GM1 With

LE/LC

equilibrium

70.4±0.2

surface 61.2±0.3

pressure, with GM1

18 ACS Paragon Plus Environment

Page 18 of 33

Page 19 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure5: Area per lipid of the monolayer. (a) without GM1 without surface pressure (b) without GM1 with LE/LC equilibrium surface pressure (c) with GM1 without surface pressure (d) with GM1 with LE/LC equilibrium surface pressure. Lateral diffusion coefficient The lateral diffusion coefficient of different component of the binary (DMPC-CHOL) and ternary (DMPC-CHOL-GM1) systems in presence and absence of LE/LC equilibrium surface pressure are given in Table 2. Table 2: Diffusion coefficient of different lipid components of different systems.(10-7cm2s-1) Systems (a)Without Surface pressure, without GM1 (b)With LE/LC equilibrium surface pressure, without GM1 (c)Without Surface pressure, with GM1 (d)With LE/LC equilibrium surface pressure, with GM1

DMPC

CHOL

1.77±0.20

1.43±0.22

1.30±0.35

1.26±0.20

2.63±0.31

1.83±0.20

1.21±0.34

0.50±0.18

The lateral diffusion coefficient of components of system (a) was in well agreement with the litrature. Czub et al. in their work report lateral diffusion coefficient of DMPC and CHOL as 0.50±0.73 cm2/s and 0.6±1.79 cm2/s respectively at 300K.67 Our value was slightly higher than the reported value which could be because of the higher temperature applied in our simulations. On considering the binary system first, it was observed that on applying LE/LC equilibrium surface pressure, the lateral diffusion coefficient values for both the components were decreased to some extent signifying slight condensation of the system. However, for the system containing GM1 devoid of any surface pressure, the lateral diffusion coeffiecint values were increased

19 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

convincingly. Persuing our idea of fluidization effect of GM1 previous, we could say that addtion of GM1 increased the lateral motion of other lipid components. For system (d) there was a remarkable decrease in lateral diffusion coefficeint of DMPC and CHOL. This might be due to the condesing effect induced by the GM1 molecules in presence of LE/LC equilibrium surface pressure. The newly formed condensed domains provide very less scope to the lipid molecules to move freely which eventually decreases the lateral diffusion coefficient. Radial distribution

Figure6: Radial distribution function g(r) profile of lipid components in DMPC-CHOL system (a,b,c) and DMPC-CHOL-GM1 system (d,e,f), without any surface pressure. Red (i) and blue (f) lines indicate the beginning and the end of the simulation, respectively. As concluded from our previous results, before the monolayer attained specific surface pressure corresponding to LE/LC equilibrium, the addition of GM1 to the monolayer caused a fluidization of the system; this was further fortified with our radial distribution function (RDF) calculations. The data were averaged over initial and final 10 ns of simulation trajectory. In case of the system

20 ACS Paragon Plus Environment

Page 20 of 33

Page 21 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

where there was no GM1 or surface pressure applied, with the progress of simulation time, it was observed that the peak height for DMPC-CHOL at 6Å was increased significantly i.e. from 1.3 to 1.9 (Figure6c), indicating a mixing between DMPC and CHOL. The same DMPC-CHOL peak for the monolayer containing GM1 was however lowered. Moreover, GM1 also caused the decrement in the peak values for DMPC-DMPC and CHOL-CHOL. As depicted in the figure, Peak for CHOL-CHOL system within the first shell was unaltered but for second shell it was lowered from 2.9 to 2.1 (Figure6d). Also, for the DMPC-DMPC system, the peak within the first shell was lowered from 5.1 to 4.1 (Figure6e). In addition to this, the peak heights of RDF for individual lipids coupled with GM1 did not increase significantly. (Figure S4 of supporting information). These results indicated that addition of GM1 to the system, in absence of any surface

pressure,

caused

a

fluidization

21 ACS Paragon Plus Environment

of

the

monolayer.

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure7: Radial distribution function g(r) profile of lipid components in DMPC-CHOL system (a,b,c) and DMPC-CHOL-GM1 system (d,e,f), with LE/LC equilibrium surface pressure. Red (i) and blue (f) lines indicate the beginning and the end of the simulation, respectively.

On calculating the RDF of monolayer devoid of GM1 having LE/LC equilibrium surface pressure applied to it, it was observed that the peak height was increased for all three pairs i.e. CHOLCHOL, CHOL-DMPC, DMPC-DMPC, increase being maximum for DMPC-CHOL pair (Figure7a,b,c). Now, this was in agreement with our LB isotherm and SEM images. LB isotherm revealed that after a certain surface pressure which is mainly the LE/LC equilibrium surface pressure (∼4.5 mN/m in our experiments) the 2:1 DMPC-CHOL monolayer started attaining LC phase and this result was also reflected in SEM image of 2:1 DMPC-CHOL monolayer taken at ∼7mN/m showing the formation of the condensed domains. Now as we stated above that on applying LE/LC equilibrium surface pressure to the monolayer containing GM1, it was found that peak heights for all three pairs i.e. CHOL-CHOL, CHOL-DMPC and DMPC-DMPC increased remarkably. In addition to this the peaks values for individual components of monolayer paired with GM1, also increased significantly (Figure S5 of supporting information). One noticeable thing here, though there was an increase in the peak values in presence of LE/LC equilibrium surface pressure for all the lipid pairs for both the systems i.e. with and without GM1, this increase was much greater in magnitude in case of the system where GM1 was present in the system. These high increments in peak values were also reflected in our SEM image where it was observed that GM1 containing monolayer at relatively high surface pressure formed larger domain. GM1-GM1 RDF peaks also increased remarkably (from 2.2 to 7.2) near 7Å pointing

22 ACS Paragon Plus Environment

Page 22 of 33

Page 23 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

towards a localized domain of GM1 within the monolayer (Figure S5 of supporting information) and reassured the formation of GM1 aggregates.63,65,66 Therefore, the results obtained from our experiments and simulations reveal that when GM1 is added to DMPC-CHOL membrane, it has a unique mixing behaviour upon the monolayer. At lower surface pressure i.e. surface pressure corresponding to the LE phase, GM1 has a fluidizing effect on the system, whereas as soon as the system reaches LE/LC equilibrium surface pressure, the GM1 imposes a condensing effect to the system and makes it more compact. This could be explained if we have a closer look at the geometry of GM1 (Figure1). GM1 itself has a fluid nature due to large headgroup consisting of four sugar groups coupled with the negatively charged sialic acid. This large head group of GM1 is the reason behind the loose packing of GM1 molecules and thus causes fluidity of the GM1.28 Therefore, the structural feature of GM1 could be the probable reason behind the monolayer fluidization at low surface pressure. On the other hand, when LE/LC equilibrium surface pressure is applied to the system, the GM1 molecules come closer and start to form aggregates with the help of cholesterol molecules. The hydrophobic interaction between the sterol ring of cholesterol molecules and tail region of GM1 promotes the co-operative packing.65,68 Moreover, due to compression, the lipid molecules come closer and intermolecular hydrogen bonding between DMPC molecules and head groups of GM1 forms. As represented in figure 7, one of the terminal glycans i.e. galactose forms the hydrogen bond with the phosphorus attached oxygen of the DMPC molecule. Earlier this kind of intermolecular hydrogen bonding has been reported between 1-palmitoyl-2-oleoyl-sn-glycero-3phosphocholine (POPC) and GM1 molecules and are found to be responsible for the condensed domain formation.35 Altogether these interactions i.e. hydrogen bonding and hydrophobic

23 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

interaction (shown in figure 7) cause the tight packing of lipid molecules and formation of larger LC domains at high surface pressure.

Figure8: Snapshot of a condensed domain showing hydrogen bonding between DMPC-GM1 and hydrophobic interaction between CHOL-GM1, taken at the end of simulation trajectory of the system having LE/LC equilibrium surface pressure applied to it. GM1, DMPC and CHOL are represented in yellow, red and blue color respectively. Conclusions: Effect of GM1 on the phase behavior of monolayer DMPC-CHOL was studied. Our experimental results suggested that initially, addition of GM1 has a fluidizing effect on the monolayer and GM1 distributes the DMPC-CHOL monolayer into small circular domains but eventually when the system attains LE/LC equilibrium, afterward GM1 aids the monolayer to bind more co-operatively. The system becomes more compact with higher stability and forms larger condensed domains. These findings were also supported by our theoretical results. Supporting Information:

24 ACS Paragon Plus Environment

Page 24 of 33

Page 25 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

TLC plate of purified ganglioside GM1 1

H NMR spectra of purified ganglioside GM1

Monolayer compression isotherms of 2:1 DMPC-CHOL and GM1 in 2:1 DMPC-CHOL at 37°C with lower concentrations and lower compression rate. Radial distribution function g(r) profile of GM1 paired with lipid components in DMPC-CHOLGM1 system, without any surface pressure and with LE/LC equilibrium surface pressure. Acknowledgements Z.S. and S.D. thank the University Grants Commission (UGC), New Delhi, for the award of a Junior Research Fellowship. This work is partially funded by the Center for Advanced Studies, Department of Chemistry, University of Calcutta. References: (1)Koldsø, H.; Shorthouse, D.; Hélie, J.;Sansom, M. S. P. Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers. PLoS Comput. Biol .2014, 10,10. (2) van Meer, G.; Voelker D. R.; Feigenson G. W. Membrane lipids: where they are and how they behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112-124. (3) van Meer, G.; de Kroon, A. I. P. M. Lipid map of the mammalian cell. J. Cell Sci. 2011, 124, 5–8. (4) Shevchenko, A.; Simons, K. Lipidomics : coming to grips with lipid diversity. Nat. Rev. Mol. Cell Biol. 2010, 11, 593-598. (5) Simons, K.; Toomre, D. Lipid rafts and signal transduction. Nat. Rev. Mol. Cell Biol. 2000, 1, 31-39.

25 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(6) McHenry, A. J.; Sciacca, M. F. M.; Brender, J. R.; Ramamoorthy, A. Does cholesterol suppress the antimicrobial peptide induced disruption of lipid raft containing membranes? Biochim. Biophys. Acta - Biomembr. 2012, 1818, 3019–3024. (7) Demetzos, C. Differential Scanning Calorimetry (DSC): A Tool to Study the Thermal Behavior of Lipid Bilayers and Liposomal Stability. J. Liposome Res. 2008, 18, 159-173. 8) McIntosh, T. J.; Simon, S. A. Long- and Short-Range Interactions between Phospholipid/Ganglioside GM1 Bilayers. Biochemistry 1994, 33, 10477–10486. (9) Tristram-Nagle, S.; Liu, Y.; Legleiter, J.; Nagle, J. F. Structure of gel phase DMPC determined by x-ray diffraction. Biophys. J. 2002, 83, 3324–3335. (10) Léonard, A.; Escrive, C.; Laguerre, M.; Payroula, E. P.; Neri, W.; Pott, T.; Katsaras, J.; Dufourc, E. J. Location of cholesterol in DMPC membranes. A comparative study by neutron diffraction and molecular mechanics simulation. Langmuir 2001, 17, 2019–2030. (11) Dimova, R.; Aranda, S.; Bezlyepkina, N.; Nikolov, V.; Riske, K. A.; Lipowsky, R. A practical guide to giant vesicles. Probing the membrane nanoregime via optical microscopy. J. Phys. Condens. Matter 2006, 18, 1151-1176. (12) Yan, C.; Thomaz, J. E.; Wang, Y.; Nishida, J.; Yuan, R.; Breen, J. P.; Fayer, M. D. Ultrafast to Ultraslow Dynamics of a Langmuir Monolayer at the Air/Water Interface Observed with Reflection Enhanced 2D IR Spectroscopy. J. Am. Chem. Soc. 2017, 139, 16518-16527. (13) Lin, W.C.; Blanchette, C. D.; Ratto, T. V.; Longo, M. L. Lipid Asymmetry in DLPC/DSPCSupported Lipid Bilayers: A Combined AFM and Fluorescence Microscopy Study. Biophys. J. 2006, 90, 228-237. (14) Aussenac, F.; Laguerre, M.; Schmitter, J. M.; Dufourc, E. J. Detailed Structure and Dynamics of Bicelle Phospholipids Using Selectively Deuterated and Perdeuterated Labels. 2H NMR and Molecular Mechanics Study. Langmuir 2003, 19, 10468-10479. (15) Jämbeck, J. P. M.; Lyubartsev, A. P.; Derivation and Systematic Validation of a Refined All -Atom Force Field for Phosphatidylcholine Lipids. J. Phys. Chem. B. 2012, 116, 3164−3179.

26 ACS Paragon Plus Environment

Page 26 of 33

Page 27 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(16) Komljenovic, I.; Marquardt, D.; Harroun, T. A.; Sternin, E. Location of chlorhexidine in DMPC model membranes: a neutron diffraction study. Chem. Phys. Lipids 2010, 163, 480– 48. (17) Kucerka, N.; Marquardt, D.; Harroun, T. A.; Nieh, M. P.; Wassall, S. R.; Jong, D. H.; Schafer, L. V.; Marrink, S. J.; Katsaras, J. Cholesterol in Bilayers with PUFA Chains: Doping with DMPC or POPC Results in Sterol Reorientation and Membrane-Domain Formation. Biochemistry 2010, 49, 7485–7493. (18) Goñi, F. M.; Alonso, A. Biophysics of sphingolipids I. Membrane properties of sphingosine, ceramides and other simple sphingolipids. Biochim. Biophys. Acta 2006, 1758, 1902–1921. (19) Marquardt, D.; Williams. J. A.; Kinnun, J. J.; Kucerka, N.; Atkinson, J.; Wassall, S. R.; Katsaras, J.; Harroun, T. A. Dimyristoyl Phosphatidylcholine: A Remarkable Exception to α-Tocopherol’s Membrane Presence. J. Am. Chem. Soc. 2014, 136, 203−210. (20) Serro, A. P.; Galante, R; Kozica, A.; Paradiso, P.; Goncalves da Silva, A. M. P. S.; Luzyanina, K. V.; Fernandes, A.C.; Saramago, B. Effect of tetracaine on DMPC and DMPC + cholesterol biomembrane models: Liposomes and monolayers. Colloids Surf. B Biointerfaces 2014, 116, 63–71. (21) Castelli, F.; Sarpietro, M. G.; Rocco, F.; Ceruti, M.; Cattel, L. Interaction of lipophilic gemcitabine prodrugs with biomembrane models studied by Langmuir-Blodgett technique. J. Colloid Interface Sci. 2007, 313, 363–368. (22) Brown, D. A.; London, E. Structure and function of sphingolipid- and cholesterol-rich membrane rafts. J. Biol. Chem. 2000, 275, 17221–17224. (23) Maxfield, F. R.; Tabas, I. Role of cholesterol and lipid organization in disease. Nature 2005, 438, 612–621. (24) Hirai, M.; Iwase, H.; Arai, S.; Takizawa, T.; Hayashi, K. Interaction of gangliosides with proteins depending on oligosaccharide chain and protein surface modification. Biophys. J. 1998, 74, 1380–1387. (25) Azzo, A.; Tessitore, A.; Sano, R. Gangliosides as apoptotic signals in ER stress response. Cell Death Differ. 2006, 13, 404–414.

27 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(26) Chatterjee, C.;Mukhopadhyay, C. Melittin–GM1 Interaction: A Model for a Side-by-Side Complex. Biochem. Biophys. Res. Commun. 2002, 292, 579–585. (27) Chaves, E. P.; Sipione, S. Sphingolipids and gangliosides of the nervous system in membrane function and dysfunction. FEBS Lett. 2010, 584, 1748–1759. (28) Frey, S. L.; Chi, E. Y.; Arratia, C.; Majewski, J.; Kjaer, K.; Lee, K. Y. C. Condensing and fluidizing effects of ganglioside GM1 on phospholipid films. Biophys. J. 2008, 94, 3047– 3064. (29) Schnaar, R. L.; Gerardy-Schahn, R.; Hildebrandt, H. Sialic Acids in the Brain: Gangliosides and Polysialic Acid in Nervous System Development, Stability, Disease, and Regeneration. Physiol. Rev. 2014, 94, 461–518. (30) Yuan, C.; Furlong, J.; Burgos, P.; Johnston, L. J. The size of lipid rafts: An atomic force microscopy study of ganglioside GM1 domains in sphingomyelin/DOPC/cholesterol membranes. Biophys. J. 2002, 82, 2526–2535. (31) Lozano, M. M.; Liu, Z.; Sunnick, E.; Janshoff, A.; Kumar, K.; Boxer, S. G. Colocalization of the Ganglioside GM1 and Cholesterol Detected by Secondary Ion Mass Spectroscopy. J. Am. Chem. Soc. 2013, 135, 5620-5630. (32) McLaurin, J.; Chakrabartty, A. Membrane disruption by Alzheimer ??-amyloid peptides mediated through specific binding to either phospholipids or gangliosides. Implications for neurotoxicity. J. Biol. Chem. 1996, 271, 26482–26489. (33) Demarco, M. I.; Woods, R. J. Atomic-resolution conformational analysis of the GM3 ganglioside in a lipid bilayer and its implications for ganglioside-protein recognition at membrane surfaces. Glycobiology 2009, 19, 344–355. (34) Yuan, C.; Johnston, L. J. Distribution of Ganglioside GM1 in l-αDipalmitoylphosphatidylcholine/Cholesterol Monolayers: A Model for Lipid Rafts. Biophys. J. 2000, 79, 2768–2781. (35) Mori, K.; Mahmood, M. I.; Neya, S.; Matsuzaki, K.;Hoshino, T. Formation of GM1 ganglioside clusters on the lipid membrane containing sphingomyeline and cholesterol. J. Phys. Chem. B. 2012, 116, 5111–5121.

28 ACS Paragon Plus Environment

Page 28 of 33

Page 29 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(36) Yu, Z. W.; Jin, J.; Cao, Y. Characterization of the liquid-expanded to liquid-condensed phase transition of monolayers by means of compressibility. Langmuir 2002, 18, 4530– 4531. (37) Duncan, S. L.; Larson, R. G. Comparing Experimental and Simulated Pressure-Area Isotherms for DPPC. Biophys. J. 2008, 94, 2965–2986. (38) Nielsen, S. O.; Lopez, C. F.; Moore, P. B.; Shelley, J. C.; Klein, M. L. Molecular Dynamics Investigations of Lipid Langmuir Monolayers Using a Coarse-Grain Model. J. Phys. Chem. B. 2003, 107, 13911–13917. (39) Baoukina, S.; Monticelli, L.; Marrink, S. J.; Tieleman, D. P. Pressure-area isotherm of a lipid monolayer from molecular dynamics simulations. Langmuir 2007, 23, 12617–12623. (40) Cote, L. J.; Kim, F.; Huang, J. Langmuir-Blodgett Assembly of Graphite Oxide Single Layers. J. Am. Chem. Soc. 2009, 131, 1043–1049. (41) Basu, I.; Mukhopadhyay, C. Insights into binding of cholera toxin to GM1 containing membrane. Langmuir 2014, 30, 15244–15252. (42) Bahri, M. A.; Heyne, B. J.; Hans, P.; Seret, A. E.; Mickalad, A. A.; Hoebeke, M. D. Quantification of lipid bilayer effective microviscosity and fluidity effect induced by propofol. Biophys. Chem. 2005, 114, 53–61. (43) Kwok, D.Y.; Vollhardt, D.; Miller, R.; Li, D.; Neumann, A.W. Axisymmetric drop shape analysis as a film balance. Colloids Surf. A 1994, 88, 51-58. (44) Ismail, A. E.; Grest, G. S.; Stevens, M. J. Capillary waves at the liquid-vapor interface and the surface tension of water models. J. Chem. Phys. 2006, 125, 1-10. (45) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. (46) Best, R. B.; Zhu, X.; Shim, J.; Lopes, P. E. M.; Mittal, J.; Feig, M.; MacKerell, A. D. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ 1 and χ 2 Dihedral Angles. J. Chem. Theory Comput. 2012, 8, 3257–3273.

29 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(47) Feller, S. E.; Zhang, Y.; Pastor, R. W. Constant pressure molecular dynamics simulation: The Langevin piston method. J. Chem. Phys. 1995, 103, 4613-4621. (48) Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N-log(N) method for ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089-10092. (49) Andersen, H. C. Rattle: A ‘velocity’ version of the shake algorithm for molecular dynamics calculations. J. Comput. Phys. 1983, 52, 24–34. (50) Miyamoto, S.; Kollman, P. A. Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models. J. Comput. Chem. 1992, 13, 952-962. (51) Brooks, B. R.; Bruccoleri, R. E.; Olafson, B. D.; States, D. J.; Swaminathan, S.; Karplus, M. CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4,187-217. (52) Ohta, Y.; Yokoyama, S.; Sakai, H.; Abe, M. Membrane properties of mixed ganglioside GM1/phosphatidylcholine monolayers. Colloids Surfaces B Biointerfaces 2004, 33, 191– 197. (53) Maggio, B.; Cumar, F.;Caputto, R. Surface behaviour of gangliosides and related glycosphingolipids. Biochem. J. 1978, 171, 559–565. (54) Röefzaad, M.; Klüner, T.; Brand, I. Orientation of the GM1 ganglioside in Langmuir– Blodgett monolayers: a PM IRRAS and computational study. Phys. Chem. Chem. Phys. 2009, 11, 10140. (55) Kato, T. What is the Characteristic Time of Measurement of π-A Isotherms? Necessity of a Constant Strain Rate of Compression of Insoluble Monolayers for π-A Measurements. Langmuir 1990, 6, 870–872. (56) Phan, M. D., Lee, J., Shin, K. Collapsed States of Langmuir Monolayers. J. Oleo Sci. 2016, 2015, 385–397. (57) Li, M.; Retter, U.; Lipkowski, J. Kinetic Studies of Spreading DMPC Vesicles at the Air Solution Interface Using Film Pressure Measurements. Langmuir 2005, 21, 4356-4361.

30 ACS Paragon Plus Environment

Page 30 of 33

Page 31 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(58) Castelli, F.; Sarpietro, M. G.; Micieli, D.; Stella, B.; Rocco, F.; Cattel, L. Enhancement of gemcitabine affinity for biomembranes by conjugation with squalene: Differential scanning calorimetry and Langmuir-Blodgett studies using biomembrane models. J. Colloid Interface Sci. 2007, 316, 43–52.

(59) Guzmán, E.; Liggieri, L.; Santini, E.; Ferrari, M.; Ravera, F. Mixed DPPC-cholesterol Langmuir monolayers in presence of hydrophilic silica nanoparticles. Colloids Surf. B Biointerfaces 2013, 105, 284–293. (60) Yuan, C.; Johnston, L. J. Distribution of Ganglioside GM1 in l-αDipalmitoylphosphatidylcholine/Cholesterol Monolayers: A Model for Lipid Rafts1. Biophys. J. 2000, 79, 2768–2781. (61) de Jong, D. H.; Lopez, C. A.; Marrink, S. J. Molecular view on protein sorting into liquidordered membrane domains mediated by gangliosides and lipid anchors. Faraday Discuss. 2013, 161, 347–363. (62) Hung, W. C.; Lee, M. T.; Chen, F. Y.; Huang, H. W. The Condensing Effect of Cholesterol in Lipid Bilayers. Biophys. J. 2007, 92, 3960–3967. (63) Wydro, P.; Knapczyk, S.; Łapczyńska, M. Variations in the condensing effect of cholesterol on saturated versus unsaturated phosphatidylcholines at low and high sterol concentration. Langmuir 2011, 27, 5433–5444. (64) Mori, K.; Mahmood, M. I.; Neya, S.; Matsuzaki, K.; Hoshino, T. Formation of GM1 ganglioside clusters on the lipid membrane containing sphingomyeline and cholesterol. J. Phys. Chem. B 2012, 116, 5111–5121. (65) Basu, I.; Mukhopadhyay, C. In silico phase separation in the presence of GM1 in ternary and quaternary lipid bilayers. Phys. Chem. Chem. Phys. 2015, 17, 17130–17139. (66) Marushchak, D.; Gretskaya, N.; Mikhalyov, I.; Johansson, L. B. Å. Self-aggregation - An intrinsic property of GM1 in lipid bilayers. Mol. Membr. Biol. 2007, 24, 102–112.

31 ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(67) Czub, J.; Baginski, M. Comparative Molecular Dynamics Study of Lipid Membranes Containing Cholesterol and Ergosterol. Biophys. J. 2006, 90, 2368-2382. (68) Mao, Y.; Shang, Z.; Imai, Y.; Hoshino, T.; Tero, R.; Tanaka, M.; Yamamoto, N.; Yanagisawa, K.; Urisu, T. Surface-induced phase separation of a sphingomyelin/cholesterol/ganglioside GM1-planar bilayer on mica surfaces and microdomain molecular conformation that accelerates Aβ oligomerization. Biochim. Biophys. Acta - Biomembr. 2010, 1798, 1090–1099.

32 ACS Paragon Plus Environment

Page 32 of 33

Page 33 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Table of contents only

33 ACS Paragon Plus Environment