Line Activity of Ganglioside GM1 Regulates the ... - ACS Publications

Mar 21, 2017 - National University of Science and Technology “MISiS”, 4 Leninskii Prospekt, Moscow, 119049 Russia. §. Moscow Institute of Physics...
0 downloads 0 Views 2MB Size
Subscriber access provided by Indiana University

Article

LINE ACTIVITY OF GANGLIOSIDE GM1 REGULATES RAFT SIZE DISTRIBUTION IN A CHOLESTEROL-DEPENDENT MANNER Timur R. Galimzyanov, Anna S. Lyushnyak, Veronika V. Aleksandrova, Liudmila A. Shilova, Ilya I. Mikhalyov, Irina M. Molotkovskaya, Sergey A. Akimov, and Oleg V. Batishchev Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b00404 • Publication Date (Web): 21 Mar 2017 Downloaded from http://pubs.acs.org on March 24, 2017

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 free 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 accessible to all readers and 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.

Langmuir 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 30

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

Abstract 45x29mm (300 x 300 DPI)

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

Figure 1. Schematic cross-section of the raft and the surrounding membrane by the plane, which is perpendicular to the phase boundary. The raft is shown in grey color. Cartesian coordinate system is introduced: Ox axis is perpendicular to the raft boundary, Oy axis is parallel to the boundary, Oz axis is perpendicular to the membrane plane. L is the distance between the boundaries of the domain monolayer in the top and bottom membrane leaflets. The values related to the bottom monolayer we denote by the index “b”, and those for the top monolayer by the index “t”. The thickness of undeformed Lo monolayer is hr, the thickness of undeformed Ld monolayer is hs. In the zone I (surrounding membrane), thicknesses of monolayers coincide and are equal to hs. In the raft zone III, they are equal to hr. In the intermediate zone II, thicknesses of top and bottom monolayers are different: hr for the top monolayer and hs for the bottom one. 82x58mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

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. AFM topography images of membranes prepared from DOPC:eSM:Chol (1:1:1, by mole) mixture containing ganglioside GM1 with indicated concentrations in mol.%. 87x129mm (300 x 300 DPI)

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

Figure 3. Dependence of average Lo domain size on GM1 content (in mole %) in the membrane from DOPC:eSM:Chol (1:1:1, by mole). 201x141mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30

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 4. The dependence of the raft boundary line tension on the position of ganglioside-rich stripe. The following parameters were used: the equilibrium thickness of Lo monolayer hr = 1.8 nm and Ld monolayer hs = 1.3 nm; spontaneous curvatures of the raft and the surrounding membrane were assumed to be zero (Jr = 0, Js = 0); splay modulus of Lo and Ld monolayers were BR = 30 kBT and BS = 10 kBT, respectively; tilt modulus Kt = 10 kBT/nm2 = 40 mN/m; ganglioside-rich stripe width l = 1 nm, its spontaneous curvature Jp = 0.1 nm–1. 108x156mm (300 x 300 DPI)

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

Figure 5. Dependence of the line tension of the raft boundary on the total ganglioside content for two sets of parameters: (A) our system, DOPC:eSM:Chol 1:1:1; (B) the system described in ref49, DOPC:eSM:Chol 2:2:1. 56x20mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 6 of 30

Page 7 of 30

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

LINE ACTIVITY OF GANGLIOSIDE GM1 REGULATES RAFT SIZE DISTRIBUTION IN A CHOLESTEROL-DEPENDENT MANNER

T.R. Galimzyanov†,‡, A.S. Lyushnyak†,ǁ, V.V. Aleksandrova‡, L.A. Shilova†,ǁ, I.I. Mikhalyov§, I.M. Molotkovskaya§, S.A. Akimov†,‡, O.V. Batishchev†,ǁ*



A.N. Frumkin Institute of Physical Chemistry and Electrochemistry, Russian Academy of

Sciences, 31/4, Leninskii prospekt, Moscow, 119071 Russia ‡

National University of Science and Technology “MISiS”, 4, Leninskii prospekt, Moscow, 119049

Russia ǁ

Moscow Institute of Physics and Technology, 9, Institutskii per., Dolgoprudnyi, Moscow region,

141700 Russia §

M.M. Shemyakin and Yu.A. Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of

Sciences, 16/10, Miklukho-Maklaya str., Moscow, 117997 Russia

* Corresponding author. E-mail: [email protected]; Fax: +7-495-952-55-82; Tel: +7-495-95486-73 (O.V.B.)

ABSTRACT

Liquid-ordered lipid domains, also called rafts, are assumed to be important players in different cellular processes, mainly signal transduction and membrane trafficking. They are thicker than the disordered part of the membrane and are thought to form for compensating the hydrophobic mismatch between transmembrane proteins and lipid environment. Despite the existence of such structures in vivo is still an open question, they are observed in model systems of multicomponent lipid bilayers. Moreover, the predictions obtained from model experiments allow explaining 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

different physiological processes possibly involving rafts. Here we present the results of the study of regulation of raft size distribution by ganglioside GM1. Combining atomic force microscopy with theoretical considerations basing on the theory of membrane elasticity, we predict that this glycolipid should change the line tension of raft boundaries in two different ways, mainly depending on cholesterol content. These results explain the shedding of gangliosides from the surface of tumor cells and the following ganglioside-induced apoptosis of T-lymphocytes in a raftdependent manner. Moreover, the generality of the model allows predicting line activity of different membrane components based on their molecular geometry.

CARTOON FOR ABSTRACT

INTRODUCTION

Lipid matrix of cell membranes is highly heterogenic due to a large variety of lipid species1. Although the fluidity of biological membranes is one of the main properties determining their normal functionality, there is a possibility of formation of ordered lipid domains differing by structure and composition from the rest of the membrane2-4. These domains, also called rafts, are responsible for proper biological functions of certain membrane proteins5-8, involved in signal transduction in cells9,10, and play a significant role in processes of membrane trafficking11,12, cytoskeletal organization13,14, and even viral infection15. They are thicker than the other (liquid2 ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30

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

disordered) part of the membrane16-18. Deformations arising to compensate the thickness mismatch determine an energy penalty per unit length of the boundary between raft and surrounding membrane called line tension19,20. Acting against entropy, line tension regulates the equilibrium size distribution of domains17,21. It is commonly believed that the mission of ordered lipid domains is to compensate hydrophobic mismatch between thick transmembrane part of the proteins and surrounding membrane22. That is why many authors contribute raft-forming function to the proteins23-25. Nevertheless, in the past two decades the role of lipids in emerging the domain formation has been established26-28. As Simons and Ikonen state5, rafts should be enriched in cholesterol and sphingolipids. There are also many evidences of formation of liquid-ordered lipid domains in ternary mixtures without sphingolipids, but containing other lipids with saturated hydrocarbon tails18,29. This means that lipids can form a large variety of domains, which, in turn, may influence protein clustering and cellular processes30. The regulating role of lipids in these events is still under debates, mainly due to difficulties in study of rafts in different systems. The well-known discrepancy exists between the results obtained in model systems and in cell membranes, since nobody sees micron-sized lipid domains in cell membranes, while they commonly exist in model ones31-33. The first problem is impossibility of observation of nanoscopic domains by optical techniques. Other methods, such as atomic force microscopy, requires specific model systems, e.g., supported lipid monolayers or bilayers16,34, while biochemical techniques give contradictory results35,36. Some authors describe rafts as ensemble of nanoscale domains, merging into micrometer-sized assemblies under certain stimulus19,24,37. This possibility of domain merging could be finely regulated by lipid chemical structure38 or membrane composition39. On the other hand, the size of domains is dictated by the line tension of their boundaries. This parameter and corresponding raft morphology strongly depend on the lipid composition38,39. Several lipid species, called line-active components40,41, can powerfully influence the line tension at raft boundaries in a concentration-dependent way. For example, gangliosides, having high affinity to the liquid-ordered phase of the membrane, vary the size of lipid domains42. They can also form gel-phase domains in 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

the membrane by its own42,43, or in combination with other membrane components44. These glycosphingolipids are involved in a large variety of cellular processes, such as reception, neuronal differentiation, and the formation of axons, dendrites, and synapses45. Moreover, in refs46,47 authors show that these lipids shedded from cancer cells may trigger the apoptosis of T-lymphocytes. Monosialoganglioside GM1, the most studied representative of ganglioside family, demonstrates activity on domain behavior in different model systems, from lipid monolayers to giant unilamellar vesicles34,48-50. However, the question remains unanswered is the ability of GM1 to change size distribution of lipid domains dramatically, being in very small concentration, fractions of molar percent49. It means that this lipid should accumulate in a narrow zone near domain boundary, significantly influencing the line tension. Since GM1 hydrocarbon tails are fully saturated, we cannot associate this effect with the model of line activity suggested in ref41, which attributed the line activity to “hybrid” lipids possessing one saturated and one unsaturated tail only. Recently we have demonstrated that two leaflets of bilayer lipid domain should not be completely in register, so there is a narrow transition zone around the domain necessary to compensate elastic deformations arising from height mismatch at domain boundary51. Numerous molecular dynamics simulations52-54 support this result. In the present study, we conducted series of experiments on supported bilayer lipid membranes using atomic force microscopy and combined them with theoretical calculations to show the effect of small fractions of GM1 on size distribution of liquid-ordered domains depending on lipid composition of the membrane.

MATERIALS AND METHODS

A. Experimental section Materials. 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC), egg sphingomyelin (eSM), Npalmitoyl-D-erythro-sphingosylphosphorylcholine (pSM) and cholesterol (Chol) were purchased from Avanti Polar Lipids (USA). Monosialoganglioside (GM1), methanol (> 99.0%), chloroform (> 4 ACS Paragon Plus Environment

Page 10 of 30

Page 11 of 30

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

99.0%), sodium chloride (NaCl) and HEPES were purchased from Sigma (USA). All chemicals were used without further purification. In experiments, we used Milli-Q water with a resistivity of 18.2 MΩ⋅cm. Preparation of supported lipid bilayer (sBLM). Liposomes were prepared as described in ref55. In brief, lipid solution containing DOPC, eSM (pSM in control experiments) and Chol dissolved in chloroform in a molar ratio 1:1:1 or 2:2:1 plus 0.25-2 mol. % of GM1 dissolved in chloroform : methanol (3:1 by volume) with a final lipid concentration of 1 mg/ml was rotary evaporated on the bottom of a glass flask and kept under vacuum for one hour. The resulting dried lipid film was resuspended in the working buffer solution (100 mM NaCl, 10 mM HEPES, pH 7.1). After 10 freeze-thaw cycles for 5 min each and vigorous vortexing the solution was extruded 19 times through two polycarbonate membranes with pore diameter of 100 nm (Avestin, Canada) to give unilamellar liposomes. The supported bilayer lipid membrane (sBLM) was created by deposing the liposomes on freshly cleaved mica (200 µl/~1.75 cm2) followed by incubation for 1 h at room temperature (RT) to allow lipid bilayer formation on the support. Then the sample was rinsed five times with working buffer solution and leaved covered by it to prevent bilayer dehydration. Atomic force microscopy (AFM). AFM experiments were performed using Multimode Nanoscope IV setup (Veeco Digital Instruments, USA) equipped with J type scanner and electrochemical fluid cell. All experiments were performed in tapping mode in working buffer solution. For scanning, we used Ultrasharp SiN3 cantilevers with nominal spring constant of 0.06 N/m and tip radius of approximately 2 nm (Bruker, USA). Z-axis resolution was 0.1 nm. After placing of a sample into the cell and filling of the cell with working buffer solution the sample was heated to 50°C with the rate of 2°C per 1 min using Thermal Applications Controller (Bruker, USA) and then incubated for 10 min at this temperature to allow all lipids to melt to liquid-disordered state. After that, the sample was cooled down to 18°C with the rate of 1°C per min, again incubated for 10 min and scanned. All images were scanned with dimensions of 3×3 µm2 and processed using WSxM software56 after scanning. 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

B. Theoretical section We consider the bilayer lipid membrane as a continuous liquid-crystalline elastic medium in which liquid-ordered domains (rafts) form because of the phase separation19,20. For simplicity, we assume that the size of domains is large enough, so we can consider their boundaries locally as straight lines. Thus, the system possesses translational symmetry along the raft boundary, and we can treat it as effectively one-dimensional. This assumption is valid if the radius of the raft is considerably greater than the characteristic length of elastic deformations, which is about several nanometers51. Numerous studies report that thickness of raft is larger than the thickness of the surrounding membrane16-18. If raft monolayers are homogeneous and flat up to the boundary, jump in thickness of high energy will exist at the boundary. To prevent exposure of hydrophobic membrane interior to water raft monolayers should deform at the boundary51,57. We analyze the energy of membrane deformations at the domain boundary in the framework of theory of elasticity of liquid crystals adapted to lipid membranes58. We take into account two fundamental membrane deformations: splay and tilt. We minimize the elastic energy functional with respect to spatial distribution of deformations, obtain and solve Euler-Lagrange equations under appropriate boundary conditions. Substituting the deformations back into the energy functional and integrating over the monolayer surfaces, we obtain the energy of the boundary (per unit length), or line tension. As we show in refs51,57, the line tension of liquid ordered domain boundary is minimal when opposed monolayer domains are not completely in register, but relatively shifted by some equilibrium distance of about 2-4 nm (Fig. 1). In the framework of “elastic” approach, we consider the ganglioside, having bulk polar part and relatively small hydrophobic region, as the component possessing positive spontaneous curvature. Ganglioside can laterally distribute into liquid-ordered, Lo, and liquid disordered, Ld, phases, as well as into the intermediate zone between them (zone II in Fig. 1). However, GM1 is known to be the marker of Lo phase34. Thus, we assume that its content in Ld phase is negligible. It seems reasonable that to provide the line activity the gangliosides should be 6 ACS Paragon Plus Environment

Page 12 of 30

Page 13 of 30

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

enriched in a narrow stripe-like zone close to the raft boundary. Only in this case boundary structure can be changed significantly by small amounts of GM1. The exact distribution of the GM1 between the stripe and the bulk raft phase, the width of the stripe, and its lateral position we find from the condition of minimal elastic energy of the system.

Figure 1. Schematic cross-section of the raft and the surrounding membrane by the plane, which is perpendicular to the phase boundary. The raft is shown in grey color. Cartesian coordinate system is introduced: Ox axis is perpendicular to the raft boundary, Oy axis is parallel to the boundary, Oz axis is perpendicular to the membrane plane. L is the distance between the boundaries of the domain monolayer in the top and bottom membrane leaflets. The values related to the bottom monolayer we denote by the index “b”, and those for the top monolayer by the index “t”. The thickness of undeformed Lo monolayer is hr, the thickness of undeformed Ld monolayer is hs. In the zone I (surrounding membrane), thicknesses of monolayers coincide and are equal to hs. In the raft zone III, they are equal to hr. In the intermediate zone II, thicknesses of top and bottom monolayers are different: hr for the top monolayer and hs for the bottom one.

The spontaneous curvature of homogeneous monolayer of non-raft phase, JS, does not depend on the coordinates and concentration of ganglioside. In Lo phase, the spontaneous curvature, JR, is assumed to be proportional to the concentration of ganglioside, cR:

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

J R = J R 0 (1 − cR ) + J G cR ,

Page 14 of 30

(1)

where JG is the spontaneous curvature of pure GM1 monolayer, JR0 is the spontaneous curvature of Lo phase monolayer without ganglioside. The spontaneous curvature of the stripe with high content of ganglioside GM1 we determine similarly:

J st = J R 0 (1 − cst ) + J G cst ,

(2)

where cst is the ganglioside concentration in the stripe. GM1 concentration in the bulk raft phase and raft stripe in each monolayer should obey the condition of matter conservation:

S R cR + Lst d st cst = S0cG ,

(3)

where SR is the raft phase area, Lst is the length of the boundary of the raft phase, dst is the width of the stripe enriched by ganglioside GM1, S0 is the total area of the lipid monolayer surface, cG is the total concentration of ganglioside in the system. The location, width and composition of the stripes containing ganglioside are determined by minimizing the total elastic energy under the condition of Eq. (3). The resulting expressions are analytical, but very bulky, so we will illustrate the results graphically in the following section.

RESULTS AND DISCUSSION

AFM experiments. We performed AFM experiments for the membrane from “canonical” raftforming mixture of DOPC:SM:Chol (1:1:1) with and without ganglioside GM1. Resulting topography images are shown in Figure 2. In absence of GM1 we observed two clearly distinguishable lipid phases with the height difference of about 1 nm, meaning the formation of ordered domains in our system22-24. The average size of the higher phase domains and their fraction were in in good agreement with reported values for this system39. Although shape of domains was not perfectly round, it was also far from dendrite-like or fractal, seen in monolayers47. This fact, together with large average size of domains, allows concluding that: (i) domains were in liquid state; (ii) line tension of domain boundaries was rather high25. In ref39 authors demonstrate that Lo 8 ACS Paragon Plus Environment

Page 15 of 30

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

domains merge over time. It means that the area of the individual domains may change in time with the preserved conditions of the system. After such process, size of initial domains could be distinguished only near the merged domain boundary. Therefore, we overlapped circles of different diameter over domains of liquid-ordered phase near its boundary to find the best match (see Supporting Information, Fig. S1), and average diameter of such circles was set as the raft diameter for given lipid composition (Fig. 3). This phase pattern changed dramatically when the system contained GM1. For 0.25 mol. % of GM1 we observed fragmentation of the domains to smaller ones, with average diameter near 60 nm, and the formation of circular domains with the size of approximately 40 nm in diameter for 0.5 mol. % of GM1 (Fig. 3). Such abrupt reduction of average size of Lo domains indicates the sharp decrease of line tension of phase boundary21. It is generally accepted that GM1 preferentially redistributes to the liquid-ordered part of the membrane34. Note, that the distance between rafts was of the same order of magnitude as their size (tens of nanometers), so such pattern of rafts in membrane cannot be observed using conventional optical fluorescent microscopy50, resulting in artificial treatment of such pattern as homogeneous lipid membrane. This observation supports the idea that only part of real phase diagram for multicomponent lipid system can be detected optically50. At 0.75 mol. % of GM1 rafts became bigger again, and for 1 mol. % of GM1 in the membrane we observed the Lo domains similar in shape and size to the case of GM1-free system (Fig. 3). Some clusters of height near 2.5 nm and average diameter of about several nanometers appeared for 0.75 mol. % and higher content of GM1 (see Fig. 2). According to ref49, such clusters could be GM1 domains since this glycolipid forms domains by its own, or gel-phase clusters of GM1 with SM and Chol44. Another possibility is that the clusters resulted from long-tail fraction of egg SM, since this lipid extract from Avanti contains 3 mol. % of 22:0 fraction and 3 mol. % of 24:1 fraction. To exclude the case we performed the same experiments with synthetic pSM and obtained similar pattern of phase separation (see Supporting Information, Fig. S2), so long-tail fractions of egg SM are not responsible for formation of the above mentioned clusters. The increased average size of domains indicates the growth of line 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

tension at their boundaries. For 2 mol. % of GM1 the average domain size and the morphology of domains almost totally coincide with that for the ganglioside-free raft-forming mixture (Fig. 2), with exception of some higher clusters staying near the domain boundaries. Therefore, line tension, being decreased at 0.25 mol. % and 0.5 mol. % of GM1, started to grow up for higher amounts of GM1 in the system. Increase of line tension for high GM1 content (up to 10 mol. %) is observed by fluorescent microscopy in giant unilamellar vesicles48. At the same time, this fact contradicts the results of the ref49, where authors show monotonic decrease of domain size with increasing molar fraction of GM1 in the membrane. However, the authors used different raft-forming mixture of lipids (DOPC:eSM:Chol ratio is 2:2:1 by mole). For clarity reasons we repeated their experiments for 0, 0.5 and 1 mol. % of GM1 and obtained the same results as in the ref. 49 (see Supporting Information, Fig. S3). In the next section, we analyzed the dependence of domain line tension on GM1 concentration basing on thickness mismatch model of raft boundary structure51,57.

10 ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30

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. AFM topography images of membranes prepared from DOPC:eSM:Chol (1:1:1, by mole) mixture containing ganglioside GM1 with indicated concentrations in mol.%.

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

Figure 3. Dependence of average Lo domain size on GM1 content (in mole %) in the membrane from DOPC:eSM:Chol (1:1:1, by mole).

Theoretical model. To clarify the non-linear effect of GM1 on the line tension of raft boundary observed in AFM experiments, we studied the distribution of ganglioside in the raft phase using the continuum model51,57. In this model, we assumed the energy of the raft boundary to be determined by deformations arising in order to compensate the thickness mismatch between raft and surrounding membrane. We calculated the energy of the raft boundary as a function of the ganglioside position in the membrane. A typical dependence of the line tension, γ, on the location of GM1 is shown in Fig. 4. We obtained that the optimal location of ganglioside is a narrow stripe is the top (Lo) monolayer of the intermediate zone II (see Fig. 1). It correlates with the fact that the transition zone of raft boundary is the most deformed region of lipid bilayer and ganglioside allows partial relaxing of the elastic stress. Thus, in further calculations we considered only one ganglioside-rich stripe located in the top monolayer of the transitional zone II (Fig. 1).

12 ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30

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 4. The dependence of the raft boundary line tension on the position of ganglioside-rich stripe. The following parameters were used: the equilibrium thickness of Lo monolayer hr = 1.8 nm and Ld monolayer hs = 1.3 nm; spontaneous curvatures of the raft and the surrounding membrane were assumed to be zero (Jr = 0, Js = 0); splay modulus of Lo and Ld monolayers were BR = 30 kBT and BS = 10 kBT, respectively; tilt modulus Kt = 10 kBT/nm2 = 40 mN/m; ganglioside-rich stripe width l = 1 nm, its spontaneous curvature Jp = 0.1 nm–1.

We calculated the dependence of line tension of the raft boundary on the overall content of GM1 in the membrane (Fig. 5) for two sets of parameters. The first set (Fig. 5A) corresponded to our experimental system (DOPC:eSM:Chol 1:1:1): BR = 30 kBT (kBT ∼ 4×10–21 J), BS = 10 kBT, Kt = 10 kBT/nm2 = 40 mN/m, JG = 0.25 nm–1, hs = 1.3 nm, hr = 1.8 nm, JS = –0.14 nm–1, JR = –0.15 nm–1. The second set of parameters (Fig. 5B) corresponded to the experimental conditions of the work by R. Bao et al.49 for the DOPC:eSM:Chol 2:2:1 system: BR = 30 kBT, BS = 10 kBT, Kt = 10 kBT/nm2, JG 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

Page 20 of 30

= 0.25 nm–1, hS = 1.3 nm, hR = 1.8 nm, JS = –0.09 nm–1, JR = –0.06 nm–1. Values of monolayer elastic moduli were taken from the refs59,60, spontaneous curvature of GM1 JG = 0.25 nm–1 was taken from the ref61, spontaneous curvature of all other components were taken from refs62,63: Js(Chol) = –0.44 nm-1, Js(DOPC) = –0.11 nm-1, Js(SM) = 0.12 nm-1. Spontaneous curvatures of the phases were set as concentration-weighted sum of spontaneous curvatures of their components. Phase composition is taken from DOPC:eSM:Chol ternary phase diagram measured in ref39 and summarized in table 1.

Table 1. Phase compositions and spontaneous curvatures. Component ratio (DOPC:SM:Chol) 1:1:1 2:2:1

Liquid-ordered phase Composition Spontaneous curvature 5% DOPC, –0.15 nm–1 48% eSM, 47% Chol 2% DOPC, –0.06 nm–1 67% eSM, 31% Chol

Liquid-disordered phase Composition Spontaneous curvature 72% DOPC, –0.14 nm–1 11% eSM, 17% Chol 75% DOPC, –0.09 nm–1 18% eSM, 7% Chol

Figure 5. Dependence of the line tension of the raft boundary on the total ganglioside content for two sets of parameters: (A) our system, DOPC:eSM:Chol 1:1:1; (B) the system described in ref49, DOPC:eSM:Chol 2:2:1.

In our experiments, small fraction of ganglioside significantly decreased the line tension of the raft boundary (see Fig. 2). It means that up to the concentration of about 0.5 mol. % this line-active 14 ACS Paragon Plus Environment

Page 21 of 30

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

component distributed exclusively to the raft boundary region. It resulted in increase of the width of the ganglioside-rich stripe and its spontaneous curvature. Although the total concentration of GM1 in the membrane is relatively small, its local concentration in the stripe near the domain boundary may reach tens of molar percent. This results in the formation of “gel”-phase clusters of ganglioside with eSM and Chol44 near the domain boundary already at 0.75 mol. % of GM1. Since the total area of such clusters is rather small (see Fig. 2), we assumed that GM1 mostly distributes between raft boundary and domain bulk region. However, the line tension trend for GM1 content higher than 0.5 mol. % depends on the initial spontaneous curvature of monolayer of raft and non-raft phases, i.e. on composition of the membrane (Fig. 5). Composition of phases in 1:1:1 and 2:2:1 systems differs mainly in SM and Chol amount (see Table 1). Absolute value of cholesterol spontaneous curvature (Js = –0.44 nm-1) is about four times larger than for SM (Js = 0.12 nm-1). Thus, cholesterol content generally determines the overall spontaneous curvature of the phase. After the threshold concentration of GM1 in the stripe, larger negative spontaneous curvature of phases in 1:1:1 system leads to upward trend of line tension (Fig. 5A). Relatively low spontaneous curvature in 2:2:1 membrane leads to reverse trend of line tension (Fig. 5B) described in the ref49 as decrease of average size of Lo domains.

CONCLUSIONS

In the present work, we developed the analytical model describing the behavior of liquid-ordered domains in the presence of line-active components such as ganglioside GM1. In the framework of the model, the major property responsible for line activity of the component was its non-zero spontaneous curvature, i.e. conical or inverted-conical effective shape of the molecule. Results of calculations predict that the influence of GM1 on the pattern of domain size distribution in lipid membrane greatly depends on the initial composition of the system described in terms of spontaneous curvature of its components. Cholesterol possesses significantly more pronounced 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

spontaneous curvature of all the components used in raft-forming lipid mixtures, thus its content generally determines the overall spontaneous curvature of a phase. Therefore, we can conclude that in two cases presented here, the main factor determining the behavior of the system was amount of cholesterol: depending on its content, GM1 at certain concentration can either decrease or increase the line tension of the Lo domain boundary. The former case is observed together with the high content of sphingomyelin typical for outer leaflet of plasma membrane. This lipid forms gel-like domains even at physiological temperature (37°C), and this effect is the most pronounced in the absence of macroscopic Lo domains, as it is shown by raft disruption by cholesterol depletion64. The coexistence of gel and liquid-disordered lipid phases is deadly for cells since this leads to increase of ion permeability of plasma membrane65. On the other hand, for cells with high content of cholesterol, such as tumor cells66, the effect of gangliosides is reverse, as for the case of our DOPC:eSM:Chol 1:1:1. For concentration of ganglioside higher than 0.5 mol. % we should obtain almost linear increase of the line tension of the raft boundary leading to budding of raft vesicles at concentrations of GM1 higher than 2 mol. %48. The interplay of the observed two distinct dependencies of raft line tension on GM1 content one can see for the tumor-induced apoptosis in cytotoxic T-lymphocytes67. It is well established that tumor cells shed ganglioside and gangliosidecontaining microvesicles from their surface68, and we could explain this effect in the frameworks of the proposed model by high line tension of Lo domains. Incorporation of shedded gangliosides into normal cell membranes, such as plasma membranes of activated T-lymphocyte, should reduce the line tension of raft boundaries, for DOPC:eSM:Chol 2:2:1 leading to vanishing of Lo domains and following formation of gel-like domains of sphingomyelin64 or GM1-eSM heterodimers49. This process, acting similar to cholesterol extraction from cell membranes, will increase ion permeability of plasma membrane65, calcium mobilization in T-lymphocytes69, and disrupt raft-associated signal pathways necessary for normal functioning of T-cell70. All of the observed results can lead to cell death.

16 ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30

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

The last thing noteworthy is the fact that we treat the ganglioside in our system just as a component with given spontaneous curvature. Thus, the proposed model is not limited exclusively to GM1, and applicable to the wide range of non-bilayer components of lipid membrane, which could also influence the raft size distribution and interaction.

ACKNOWLEDGEMENTS

Authors thank Professor Sarah L. Keller for fruitful suggestions for the work and Dr. V.I. Zolotarevskiy for great technical assistance. The work was supported by the Russian Foundation for Basic Research (projects #16-34-01203, #15-54-15006), by the Ministry of Education and Science of the Russian Federation in the framework of Increase of Competitiveness Program of “MISiS”, and by grant of the President of the Russian Federation MK-6058.2016.4.

SUPPORTING INFORMATION AVAILABLE Figures S1-S3. The Supporting Information is available free of charge via the Internet at http://pubs.acs.org/.

REFERENCES

(1) 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. (2) Simons, K.; van Meer, G. Lipid sorting in epithelial cells. Biochemistry 1988, 27, 6197–6202. (3) Lagerholm, B.C.; Weinreb, G.E.; Jacobson, K.; Thompson, N.L. Detecting microdomains in intact cell membranes. Annu. Rev. Phys. Chem. 2005, 56, 309–336. (4) Ayuyan, A.G.; Cohen, F.S. Raft composition at physiological temperature and pH in the absence of detergents. Biophys. J. 2008, 94, 2654-2666. 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

(5) Simons, K.; Ikonen, E. Functional rafts in cell membranes. Nature 1997 387, 569–572. (6) Anderson, R.G.; Jacobson, K. A role for lipid shells in targeting proteins to caveolae, rafts, and other lipid domains. Science 2002, 296, 1821–1825. (7) Pike, L.J. Rafts defined: a report on the Keystone Symposium on Lipid Rafts and Cell Function. J. Lipid Res. 2006, 47, 1597-1598. (8) Epand, R.M. Proteins and cholesterol-rich domains. Biochim. Biophys. Acta 2008, 1778, 1576– 1582. (9) Manes, S.; Mira, E.; Gomez-Moulton, C.; Lacalle, R.A.; Keller, P.; Labrador, J.P.; Martinez, A.C. Membrane raft microdomains mediate front-rear polarity in migrating cells. EMBO J. 1999, 18, 6211–6220. (10) Aman, M.J.; Ravichandran, K.S. A requirement for lipid rafts in B cell receptor induced Ca(2+) flux. Curr. Biol. 2001, 10, 393–396. (11) Brown, D.A. Rose, J. Sorting of GPI-anchored proteins to glycolipid enriched membrane subdomains during transport to the apical cell surface. Cell 1992, 6, 533–544. (12) Gonnord, P.; Blouin, C.M.; Lamaze, C. Membrane trafficking and signaling: two sides of the same coin. Semin. Cell Dev. Biol. 2012, 23, 154–164. (13) Head, B.P.; Patel, H.H.; Roth, D.M.; Murray, F.; Swaney, J.S.; Niesman, I.R.; Farquhar, M.G.; Insel, P.A. Microtubules and actin microfilaments regulate lipid raft/caveolae localization of adenylyl cyclase signaling components. J. Biol. Chem. 2006, 281, 26391–26399. (14) Suzuki, T.; Zhang, J.; Miyazawa, S.; Liu, Q.; Farzan, M.R.; Yao, W.D. Association of membrane rafts and postsynaptic density: proteomics, biochemical, and ultrastructural analyses. J. Neurochem. 2011, 119, 64–77. (15) Scheiffele, P.; Rietveld, A.; Wilk, T.; Simons, K.. Influenza viruses select ordered lipid domains during budding from the plasma membrane. J. Biol. Chem. 1999, 274, 2038–2044. (16) Rinia, H.A.; Snel, M.M.; van der Eerden, J.P.; de Kruijff, B. Visualizing detergent resistant domains in model membranes with atomic force microscopy. Febs Lett. 2001, 501, 92–96. 18 ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30

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

(17) Garcia-Saez, A.J.; Chiantia, S.; Schwille, P. Effect of line tension on the lateral organization of lipid membranes. J. Biol. Chem. 2007, 282, 33537–33544. (18) Heftberger, P.; Kollmitzer, B.; Rieder, A.A.; Amenitsch, H.; Pabst, G. In situ determination of structure and fluctuations of coexisting fluid membrane domains. Biophys. J. 2015, 108, 854-862. (19) Samsonov, A.V.; Mihalyov, I.; Cohen, F.S. Characterization of cholesterol– sphingomyelin domains and their dynamics in bilayer membranes. Biophys. J. 2001, 81, 1486–1500. (20) Baumgart, T.; Hess, S.T.; Webb, W.W. Imaging coexisting fluid domains in biomembrane models coupling curvature and line tension. Nature 2003, 425, 821–824. (21) Frolov, V.A.J.; Chizmadzhev, Y.A.; Cohen, F.S.; Zimmerberg, J. “Entropic traps” in the kinetics of phase separation in multicomponent membranes stabilize nanodomains. Biophys. J. 2006, 91, 189–205. (22) Killian, J.A. Hydrophobic mismatch between proteins and lipids in membranes. Biochim. Biophys. Acta 1998, 1376, 401–415. (23) Akimov, S.A.; Frolov, Vl.A.; Kuzmin, P.I.; Zimmerberg, J.; Chizmadzhev, Yu.A.; Cohen, F.S. Domain formation in membranes caused by lipid wetting of protein. Phys. Rev. E 2008, 77, 0519011–051901-17. (24) Fan, J.; Sammalkorpi, M.; Haataja, M. Formation and regulation of lipid microdomains in cell membranes: theory, modeling, and speculation. FEBS Lett. 2010, 584, 1678–1684. (25) Lillemeier, B.F.; Pfeiffer, J.R.; Surviladze, Z.; Wilson, B.S.; Davis, M.M. Plasma membraneassociated proteins are clustered into islands attached to the cytoskeleton. Proc. Natl. Acad. Sci. USA 2006, 103, 18992–18997. (26) Brown, D.A., London, E. Structure and function of sphingolipid-and cholesterol-rich membrane rafts. J. Biol. Chem. 2000, 275, 17221–17224. (27) Dietrich, C.; Bagatolli, L.A.; Volovyk, Z.N.; Thompson, N.L.; Levi, M.; Jacobson, K.; Gratton, E. Lipid rafts reconstituted in model membranes. Biophys. J. 2001, 80, 1417–1428.

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

Page 26 of 30

(28) Heberle, F.A.; Feigenson, G.W. Phase separation in lipid membranes. Cold Spring Harb. Perspect. Biol. 2011, 3, a004630. (29) Marsh, D. Cholesterol-induced fluid membrane domains: a compendium of lipid-raft ternary phase diagrams. Biochim. Biophys. Acta 2009, 1788, 2114–2123. (30) Mollinedo, F.; Gajate, C. Lipid rafts and clusters of apoptotic signaling molecule-enriched rafts in cancer therapy. Future Oncol. 2010, 6, 811e21. (31) Munro, S. Lipid rafts: elusive or illusive? Cell, 2003, 115, 377–388. (32) Harder, T.; Scheiffele, P.; Verkade, P.; Simons, K. Lipid domain structure of the plasma membrane revealed by patching of membrane components. J. Cell Biol. 1998, 141, 929–942. (33) Sengupta, P.; Baird, B.; Holowka, D. Lipid rafts, fluid/fluid phase separation, and their relevance to plasma membrane structure and function. Semin. Cell Dev. Biol. 2007, 18, 583–590. (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) Lichtenberg, D.; Goni, F.M.; Heerklotz, H. Detergent-resistant membranes should not be identified with membrane rafts. Trends Biochem. Sci. 2005, 30, 430–436. (36) Shogomori, H.; Brown, D.A. Use of detergents to study membrane rafts: the good, the bad, and the ugly. Biol. Chem. 2003, 384, 1259–1263. (37) Eggeling, C.; Ringemann, C.; Medda, R.; Schwarzmann, G.; Sandhoff, K.; Polyakova, S.; Belov, V.N.; Hein, B.; von Middendorff, C.; Schönle, A.; Hell, S.W. Direct observation of the nanoscale dynamics of membrane lipids in a living cell. Nature 2009, 457, 1159–1163. (38) Staneva, G.; Osipenko, D.S.; Galimzyanov, T.R.; Pavlov, K.V.; Akimov, S.A. Metabolic Precursor of Cholesterol Causes Formation of Chained Aggregates of Liquid-Ordered Domains. Langmuir, 2016, 32, 1591–1600.

20 ACS Paragon Plus Environment

Page 27 of 30

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

(39) Khadka, N.K.; Ho, C.S.; Pan, J. Macroscopic and Nanoscopic Heterogeneous Structures in a Three-Component Lipid Bilayer Mixtures Determined by Atomic Force Microscopy. Langmuir, 2015, 31, 12417–12425. (40) Trabelsi, S.; Zhang, S.; Lee, T.R.; Schwartz, D.K. Linactants: surfactant analogues in two dimensions. Phys. Rev. Lett. 2008, 100, 037802. (41) Brewster, R.; Safran, S.A. Line active hybrid lipids determine domain size in phase separation of saturated and unsaturated lipids. Biophys. J. 2010, 98, L21–L23. (42) Westerlund, B.; Slotte, J.P. How the molecular features of glycosphingolipids affect domain formation in fluid membranes. Biochim. Biophys. Acta 2009, 1788, 194–201. (43) Šachl, R.; Amaro, M.; Aydogan, G.; Koukalová, A.; Mikhalyov, I.I.; Boldyrev, I.A.; Humpolíčková J.; Hof, M. On multivalent receptor activity of GM1 in cholesterol containing membranes. Biochim. Biophys. Acta, 2015, 1853, 850–857. (44) 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 GM1planar bilayer on mica surfaces and microdomain molecular conformation that accelerates Aβ oligomerization. Biochim. Biophys. Acta, 2010, 1798, 1090–1099. (45) de Chaves, E.P.; Sipione, S. Sphingolipids and gangliosides of the nervous system in membrane function and dysfunction. FEBS Lett. 2010, 584, 1748–1759. (46) Ladisch, S.; Kitada, S.; Hays, E.F. Gangliosides Shed by Tumor Cells Enhance Tumour Formation in Mice. J. Clin. Invest. 1987, 79, 1879–1882. (47) Heitger, A.; Ladisch, S. Gangliosides Block Antigen Presentation by Human Monocytes, Biochim. Biophys. Acta 1996, 1303, 161–168. (48) Akimov, S.A.; Hlaponin, E.A.; Bashkirov, P.V.; Boldyrev, I.A.; Mikhalyov, I.I.; Telford, W.G.; Molotkovskaya, I.M. Ganglioside GM1 increases line tension at raft boundary in model membranes. Biol. Membrany, 2009, 26, 234–239.

21 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

(49) Bao, R.; Li, L.; Qiu, F.; Yang, Y. Atomic force microscopy study of ganglioside GM1 concentration effect on lateral phase separation of sphingomyelin / dioleoylphosphatidylcholine / cholesterol bilayers. J. Phys. Chem. B 2011, 115, 5923–5929. (50) Puff, N.; Watanabe, C.; Seigneuret, M.; Angelova, M.I.; Staneva, G. Lo/Ld phase coexistence modulation induced by GM1. Biochim. Biophys. Acta 2014, 1838, 2105–2114. (51) Galimzyanov, T.R.; Molotkovsky, R.J.; Bozdaganyan, M.E.; Cohen, F.S.; Pohl, P.; Akimov, S.A. Elastic membrane deformations govern interleaflet coupling of lipid-ordered domains. Phys. Rev. Lett. 2015, 115, 088101. (52) Pantano, D.A.; Moore, P.B.; Klein, M.L.; Discher, D.E. Raft registration across bilayers in a molecularly detailed model. Soft Matter 2011, 7, 8182–8191. (53) Risselada, H.J.; Marrink, S.J. The molecular face of lipid rafts in model membranes. Proc. Natl. Acad. Sci. USA 2008, 105, 17367–17372. (54) Perlmutter, J.D., Sachs, J.N. Interleaflet interaction and asymmetry in phase separated lipid bilayers: molecular dynamics simulations. J. Am. Chem. Soc. 2011, 133, 6563–6577. (55) Krenn, B.M.; Egorov, A.; Romanovskaya-Romanko, E.; Wolschek, M.; Nakowitsch, S.; Ruthsatz, T.; Kiefmann, B.; Morokutti, A.; Wressnigg, N.; Sturlan, S.; Ferko, B.; Batishchev, O.V.; Indenbom, A.V.; Zhu, R.; Kastner, M.; Hinterdorfer, P.; Humer, J.; Kiselev, O.; Muster, T.; Romanova, J. Single HA2 mutation increases the infectivity and immunogenicity of live attenuated H5N1 intranasal influenza vaccine candidate lacking NS1. PLoS ONE 2011, 6, e18577. (56) Horcas, I.; Fernandez, R.; Gomez-Rodriguez, J.M.; Colchero, J.; Gomez-Herrero J.; Baro, A.M. WSxM: A software for scanning probe microscopy and a tool for nanotechnology. Rev. Sci. Instrum. 2007, 78, 013705. (57) Galimzyanov, T.R.; Molotkovsky, R.J.; Kuzmin, P.I.; Akimov, S.A. Stabilization of bilayer structure of raft due to elastic deformations of membrane. Biol. Membrany 2011, 28, 307–314. (58) Hamm, M.; Kozlov, M.M. Elastic energy of tilt and bending of fluid membranes. Eur. Phys. J. E 2000, 3, 323–335. 22 ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30

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

(59) Rawicz, W.; Olbrich, K.C.; McIntosh, T.; Needham, D.; Evans, E. Effect of chain length and unsaturation on elasticity of lipid bilayers. Biophys. J. 2000, 79, 328–339. (60) Baumgart, T.; Das, S.; Webb, W.W.; Jenkins, J.T. Membrane elasticity in giant vesicles with fluid phase coexistence. Biophys. J. 2005, 89, 1067–1080. (61) Dimova, R.; Dasgupta, R.; Fricke, N.; Liu, Y.; Agudo-Canalejo, J.; Grafmuller, A.; Lipowsky, R. Spontaneous Tubulation in Giant Vesicles Induced by GM1 or PEG Adsorption. Biophys. J. 2016, 110, 244a. (62) Chen, Z.; Rand, R.P. The influence of cholesterol on phospholipid membrane curvature and bending elasticity. Biophys. J. 1997, 73, 267–276. (63) Venable, R. M.; Sodt, A.J.; Rogaski, B.; Rui, H.; Hatcher, E.; MacKerell, A.D., Pastor, R.W.; Klauda, J.B. CHARMM all-atom additive force field for sphingomyelin: elucidation of hydrogen bonding and of positive curvature. Biophys. J. 2014, 107, 134–145. (64) Nishimura, S.Y.; Vrljic, M.; Klein, L.O.; McConnell, H.M.; Moerner, E. Cholesterol depletion induces solid-like regions in the plasma membrane. Biophys. J. 2006, 90, 927–38. (65) Clerc, S.G.; Thompson, T.E. Permeability of dimyristoyl phosphatidylcholine/dipalmitoyl phosphatidylcholine bilayer membranes with coexisting gel and liquid-crystalline phases. Biophys. J. 1995, 68, 2333–41. (66) Krycer, J.R.; Brown, A.J. Cholesterol accumulation in prostate cancer: a classic observation from a modern perspective. Biochim. Biophys. Acta 2013, 1835, 219–229. (67) Molotkovskaya, I.M.; Kholodenko, R.V.; Molotkovsky, J.G. Influence of gangliosides on the IL-2- and IL-4-dependent cell proliferation. Neurochem. Res. 2002, 27, 761–770. (68) Kong, Y.; Li R.; Ladisch, S. Natural forms of shed tumor gangliosides. Biochim Biophys Acta 1998, 1394, 43–56. (69) Pizzo, P.; Giurisato, E.; Tassi, M.; Benedetti, A.; Pozzan, T.; Viola, A. Lipid rafts and T cell receptor signaling: a critical re‐evaluation. Eur. J. Immunol. 2002, 32, 3082–91.

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

(70) Kabouridis, P.S.; Jury, E.C. Lipid rafts and T-lymphocyte function: implications for autoimmunity FEBS Lett. 2008, 582, 3711–3718.

24 ACS Paragon Plus Environment

Page 30 of 30