Membrane tubulation in lipid vesicles triggered by the local application

Sep 14, 2017 - Experimental and theoretical studies on ion-lipid interactions, predict that binding of calcium ions to cell membranes leads to macrosc...
0 downloads 8 Views 4MB Size
Download PDF
Article pubs.acs.org/Langmuir

Membrane Tubulation in Lipid Vesicles Triggered by the Local Application of Calcium Ions Baharan Ali Doosti,# Weria Pezeshkian,⊥ Dennis S. Bruhn,⊥ John H. Ipsen,⊥ Himanshu Khandelia,⊥ Gavin D. M. Jeffries,# and Tatsiana Lobovkina*,# #

Department of Chemistry and Chemical Engineering, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden Center for Biomembrane Physics (MEMPHYS), Department of Physics, Chemistry and Pharmacy (FKF), University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark

⊥

S Supporting Information *

ABSTRACT: Experimental and theoretical studies on ion−lipid interactions predict that binding of calcium ions to cell membranes leads to macroscopic mechanical effects and membrane remodeling. Herein, we provide experimental evidence that a point source of Ca2+ acting upon a negatively charged membrane generates spontaneous curvature and triggers the formation of tubular protrusions that point away from the ion source. This behavior is rationalized by strong binding of the divalent cations to the surface of the charged bilayer, which effectively neutralizes the surface charge density of outer leaflet of the bilayer. The mismatch in the surface charge density of the two leaflets leads to nonzero spontaneous curvature. We probe this mismatch through the use of molecular dynamics simulations and validate that calcium ion binding to a lipid membrane is sufficient to generate inward spontaneous curvature, bending the membrane. Additionally, we demonstrate that the formed tubular protrusions can be translated along the vesicle surface in a controlled manner by repositioning the site of localized Ca2+ exposure. The findings demonstrate lipid membrane remodeling in response to local chemical gradients and offer potential insights into the cell membrane behavior under conditions of varying calcium ion concentrations.

■

and rigidity13,15,17,18 and triggers vesicular aggregation,19,20 in addition to membrane fusion mediation18,20,21 and spontaneous curvature induction.22,23 To date, the majority of experimental studies introduce calcium ions to lipid membranes in bulk solution, which results in a homogeneous distribution of ions over the membrane surface. However, the utilization of a localized concentration gradient has been shown to alter membrane behavior compared to bulk exposure, leading to local perturbation. For example, crista-like membrane protrusions were observed during local proton stimulation in giant vesicles, which were not present in bulk pH changes.24 Additionally, it was shown that localized exposure of Ca2+ to flat negatively charged double lipid bilayer leads to protrusive spreading in lipid membranes.25 Herein we demonstrate that the membrane of a negatively charged GUV, locally exposed to Ca2+, forms tubular protrusions, a process which is not observed when the bulk calcium concentration is increased. This membrane tubulation is described as the response to the trans-bilayer mismatch of surface charge density upon calcium ion binding from one side of the bilayer. Molecular dynamics simulations were performed on a lipid bilayer model, which demonstrated that heterogeneous calcium ion binding can be sufficient to induce inward

INTRODUCTION Calcium is an important biological agent involved in many cellular processes, ranging from cell division to apoptosis,1−3 including molecular transport,1 endo- and exocytosis,4 and communication between cells.5 The principles of calcium ion interaction with the cell membrane can be elucidated using artificial lipid bilayers to mimic multicomponent biological membranes. Examples of these artificial systems include giant unilamellar lipid vesicles (GUVs), lipid nanotubes, and supported bilayers. The lipid composition of these model membranes is engineered in such a manner to contain selected amounts of important biological lipids, tuning the ion−lipid interactions. Often, these lipids are phosphatidylcholine (PC), phosphatidylserine (PS), phosphatidylinositol (PI) or phosphatidic acid (PA). The PS lipids are of particular interest due to their role in signaling, coagulation, and cell behavior under stress conditions and apoptosis.6 Due to the ubiquitous occurrence of calcium in many cellular processes, the response of lipid membranes to ion binding has been studied extensively using complementary techniques such as 1H NMR,7 X-ray diffraction,8 thermodynamic techniques,9 and spectroscopic studies,10 in addition to molecular dynamics simulations.11,12 Experimental and theoretical studies have shown that binding of calcium ions to the lipid membrane results in lipid dehydration and tighter packing of the head groups and the hydrophobic regions of the lipids.11,13−16 This lateral membrane reorganization increases membrane tension © XXXX American Chemical Society

Received: April 28, 2017 Revised: September 7, 2017 Published: September 14, 2017 A

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir

Figure 1. Experimental setup and typical examples of spherical and undulating vesicles. (A) A schematic representation of the experimental setup. The giant unilamellar vesicle (GUV) is coupled to a multilamellar vesicle (MLV) and deposited on a glass surface in a 1 cm sized buffer droplet (not drawn). The ion gradient is established at a distance of 3−4 μm from the membrane surface by applying a positive pressure (p) through the tapered glass micropipet filled with the ion solution. The illustration is not drawn to scale. (B−D) Typical example of a spherical GUV at different time points. (E−G) Typical example of a GUV with visible thermal undulations in the membrane at different time points. The black arrows highlight the visibly deformed areas showing fluctuations in the membrane. The fluorescence images are enhanced and inverted to improve visualization of the GUV membranes. KCl (5 mM). The micropipet was used to serve as a local point-source for the ion injection (injection pressure 20−25 hPa). Upon application of a positive pressure, the content of the micropipet was slowly released into the solution creating an ion gradient. The positioning of glass micropipets to vesicles was controlled by a PatchStar micromanipulator (Scientifica, Uckfield, UK). A microinjection system (Eppendorf Femtojet) was used for point-wise application of the ions. The micropipet was positioned at a distance of 3−4 μm from the vesicle surface. Microscopy Observation and Data Analysis. An inverted fluorescence microscopy system (Leica DM IRB, Wetzlar, Germany), equipped with a 100× oil immersion objective (Leica, 100×, 1.4 NA, oil immersion), was used for observation during investigation, and a 488 nm laser line (Cobolt MLD-488 nm) enabled excitation of the ATTO488 fluorophore. A camera (Prosilica Ex 1920, Allied Vision Technologies GmbH, Thuringia, Germany) and a custom-made script in Labview 2009 (National Instruments) was used to collect the data. Fluorescence images and movies were edited and improved using NIH ImageJ software and VirtualDub 1.10.4. See Supporting Information, section S1A, for the details. The simulation figures were produced using Matlab, Inkscape 0.91, and VMD 1.9.3. Simulation Methods. Coarse grained molecular dynamics simulations were performed using GROMACS 4.6.730−32 software and the Martini force field.33,34 The size of the systems and the simulation time limited us to use a coarse grained model. Each system consists of two identical and symmetric bilayers to create two solvent compartments (Figure 4A,B). They each consist of 646 lipids, with 80% DOPC (dioleoylphosphatidylcholine) and 20% DOPS or 60% DOPC and 40% DOPS (Supporting Information, section S2B, Table S1). The negative charges from DOPS are balanced by adding Na+ equally between the two water compartments. Every tenth water particle is replaced by an “antifreeze particle” in accordance with Marrink et al.34 Additionally, in one system Na+ and Cl− are added in the middle compartment (Figure 4B) to reach a concentration of around 0.5 M while in the other Ca2+ and Cl− is added (Figure 4A) to reach around 0.25 M Ca2+ (or 0.5 M Cl−). This guarantees the same amount of positive charge in the central compartment. For each system, 20 simulations were performed to minimize statistical uncertainties in the estimated quantities. Each simulation was run for 2 μs with a time step of 30 fs. The final 1 μs was used for

spontaneous curvature and bend the membrane. Additionally, we reveal that formed tubular protrusions can be translated around the GUV membrane in the direction of increasing calcium ion concentration, by moving the source of ions along the vesicle surface.

■

EXPERIMENTAL SECTION

Chemicals. Soybean Polar Lipid Extract (SPE) and 1,2-dioleoyl-snglycero-3-phospho-L-serine (DOPS) were purchased from Avanti Polar Lipids, Inc. (Alabaster, USA). Chloroform, KCl, NaCl, and 4-(2hydroxyethyl)piperazine-1-ethanesulfonic acid (HEPES) solution (1 M) were obtained from Sigma-Aldrich (Missouri, USA). CaCl2 was purchased from KEBO lab (Sweden). MgCl2 (MgCl2·6H20) was purchased from Sharlau Chemie S.A. (Spain). 1,2-Dioleoyl-sn-glycero3-phosphoethanolamine with the fluorophore ATTO-488 (ATTO488DOPE) was purchased from ATTO-TEC (Germany). Giant Vesicle Preparation. The giant lipid vesicles were composed of SPE consisting of PE (phosphatidylethanolamine) 22.1%, PI 18.4%, PC 45.7%, PA 6.9%, other lipids 6.9% with or without addition of 20 wt % of DOPS. Both lipid mixtures included the fluorescently labeled lipid ATTO488-DOPE (1 wt %). The DOPS molecules were added to the lipid mixture to mimic the cytoplasmic side of the cell plasma membrane.26,27 Giant vesicles were prepared using a dehydration−rehydration method28 with modifications.29 First, small unilamellar vesicles (SUVs) were produced as previously described,25 after which, a droplet of the SUV solution (5 μL) was placed on a glass coverslip and dehydrated in a vacuum desiccator for 20 min. The dry lipid film was placed for 4 min at room temperature. Then, 50 μL of 10 mM HEPES buffer (pH 7.4) was used to rehydrate the dry lipid film for 5 min to form vesicles. The vesicle solution was thereafter transferred onto a new glass coverslip holding a 300 μL droplet of HEPES buffer, and vesicles were allowed settle and adhere to the glass surface for 30 min. Local Injection of Ions. Tapered glass micropipets (GC100TF-10, Clark Electromedical Instruments, Reading, UK) were pulled on a P2000 CO2 laser-puller instrument (Sutter Instruments, Novato, USA) resulting in a glass capillary with a tip diameter of 0.2−0.3 μm. The tapered glass micropipets were filled with HEPES buffer supplemented with CaCl2 (1, 2, or 5 mM), MgCl2 (5 mM), NaCl (7.5 or 75 mM), or B

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir analysis. The temperature and pressure were kept constant at 310 K and 1 bar, respectively, using the velocity-rescale thermostat35 and the Parrinello−Rahman barostat.36 All other settings are in line with the Martini standards.33,34 For a translational invariant, symmetric and continuous bilayer in a simulation box with periodic boundary conditions, the bilayer curvature does not appear as visible geometric bending. The nonzero spontaneous curvature will appear as a lateral average force asymmetry that can be observed by asymmetry in the lateral pressure profile. Therefore, we have calculated the spontaneous curvature from the simulations using the κc0 = ∫ zπ(z) dz 12,37,38 where κ is the bending rigidity of the bilayer, c0 is the spontaneous curvature, z is the axis perpendicular to the plane of the bilayer, starting from the bilayer center, and π(z) is the so-called lateral pressure profile, defined as π(z) = PL(z) − PN(z), where PL(z) = (Pxx(z) + Pyy(z))/2 is the in-plane pressure and PN(z) is the normal pressure. The lateral pressure profile is calculated using the method of Vanegas et al.43,44 The bending rigidity of the bilayers was determined from the spectrum of the bilayer shape fluctuations that takes the form k T ⟨|u(q)|2 ⟩ = 4 B 2 for large wavelengths, where u(q) is the Fourier

We utilized two standard membrane compositions, SPE/ DOPS and SPE, in order to address the effect of local calcium application to the GUV membranes at various ionic concentrations and at two different morphological regimes, (i) spherical GUVs and (ii) thermally undulating GUVs. The interaction of Ca2+ with spherical SPE/DOPS vesicles results in bending of the membrane, and the formation of dense and closely packed tubular-like structures that protrude into the lumen of the GUVs. In a typical experiment, the membrane tubulation begins with the accumulation of membrane material at the application site, which can be observed by an increase of fluorescence and small membrane indentations (Figure 2A,B).

κq + σq

transform of the bilayer midsurface in the Monge representation, kB is the Boltzmann constant, T is the temperature, and σ is the membrane tension.39,40 In order to capture the large wavelength membrane fluctuations, for each composition of lipids and ions concentration, we have simulated a system with a single large bilayer patch (Supporting Information, section S2B, Table S1, last 4 rows. For more details on the bending rigidity calculations, see Supporting Information, section S2A, Figure S3).

■

RESULTS AND DISCUSSION GUVs are an appropriate model system for studying the properties of cell membranes due to their close resemblance to mammalian plasma membranes, including size and lipid composition, and their ability to be integrated with established interrogation schemes, such as micromanipulation and electroporation.41,42 Within our study, GUVs were prepared as described previously28 in such a manner to be connected to multilamellar vesicles (MLVs), required as a source of lipids for the tube formation.43 A schematic illustration of the MLV− GUV system is shown in Figure 1A. This vesicle platform has been used in previous studies for mimicking cell tubular structures,41 exocytosis, and membrane interaction with polymers.44,45 The morphology of the prepared GUVs (typically 5−15 μm in diameter) was studied using fluorescence microscopy; the majority (∼75%) were found to be spherical, as illustrated in Figure 1B−D (see also Supporting Information, Movie 1). However, a fraction of the prepared GUVs (∼25%) were undergoing visible thermal undulations, further referred to as undulating, thus having a lower membrane tension, likely as a result of an excess of membrane material (Figure 1E−G, Supporting Information, Movie 2). For both vesicle populations, we assumed homogeneous lipid distribution with zero spontaneous curvature. Injection of one solution into the buffer leads to diffusional spreading, resulting in a concentration gradient from tip to the vesicle. To probe for this diffusional spread of ions, we used 0.1 mM fluorescein solution and monitored the fluorescence intensity profile away from the pipet tip. By measuring the normalized intensity values between the micropipet tip and the vesicle surface, we estimated a concentration decrease of approximately 75% from the initial value (i.e., from 0.1 mM to 0.025 mM) at the membrane surface from a distance of ∼3 μm. The concentration of the ions mentioned further in the text refers to the ion concentration within the micropipet, unless it is otherwise stated.

Figure 2. Formation of tubular protrusions in SPE/DOPS GUVs triggered by the local action of Ca2+ and Mg2+ ions. (A) The MLV− GUV system prior to Ca2+ (5 mM) gradient introduction. (B) Increase in fluorescence intensity of the membrane and formation of small tubular protrusions directed inside the GUV upon Ca2+ injection. (C) After an additional 3 s of exposure to Ca2+ ions, the tubulation area increases in size, and the tubular protrusions become longer. The schematic drawings on the right from the fluorescence images are simplified illustrations of formation and growth of the tubular protrusions. The images are not drawn to the scale. (D) An undulating GUV exposed to Ca2+ gradient (1 mM). Thin, long, and sparse tubules are formed at the site of Ca2+ injection. The dashed black lines outline the GUV membrane. (E) An undulating GUV exposed to a Mg2+ gradient (5 mM), which forms sparse tubular protrusions directed inside the GUV. The pipet tip in panels D and E is out of frame of the image. The fluorescence images are enhanced and inverted to improve visualization of the tubular protrusions.

Holding the micropipet close to the membrane surface while injecting Ca2+ promotes the growth of tubules (Figure 2C). These formed tubular protrusions are directed away from the membrane, inside the vesicles, that is, the membrane bends away from the ion binding site. The protrusions can reach a length of 3−5 μm after ∼5 s and continue to grow upon Ca2+ C

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir

Figure 3. Translation of membrane tubular structures along the GUV surface by moving the point source of calcium ions. (A) The MLV−GUV system prior to calcium exposure (5 mM). (B) Tubular protrusions are formed by the injection of the calcium ions. (C,D) Translating the micropipet around the vesicles surface, resulting in the entire tubulation area moving in the direction of the pipet tip. The fluorescence images are enhanced and inverted to aid in tubular protrusions visualization.

supply. They appear closely packed and dynamic and fluctuate within the vesicle, a concise video of which can be found in the Supporting Information, Movie 3. The spherical SPE/DOPS vesicles were exposed to various local Ca2+ concentrations, 1, 2, and 5 mM. Tubular protrusions formed using both 2 and 5 mM concentrations, but not at 1 mM. This indicates a threshold concentration at which protrusions are formed. We further examined undulating SPE/DOPS vesicles using the same concentration range (1, 2, and 5 mM Ca2+), and formation of tubular protrusions was observed at all three concentrations. The difference in membrane behavior at lower calcium ion concentration, when compared to the spherical GUVs, highlights the importance of the excess membrane material, which likely aids in the formation of the protrusions. The observed protrusions were generated in a shorter period of time, within 1−2 s, while appearing sparse in comparison with spherical vesicles. These lower densities of protrusions allowed the individual tubes to be resolved using fluorescence microscopy (Figure 2D, Supporting Information, Movie 4). The length of the individual tubules was measured to be 3−5 μm, similar to those produced in spherical vesicles. To elucidate the influence of DOPS molecules on the process of membrane tubulation, we excluded them from the lipid mixture and tested the SPE-only vesicles with local application of 5 mM of Ca2+. The undulating vesicles responded to calcium ion stimulation in a manner similar to the SPE/DOPS composition, forming loosely packed tubular protrusions (Supporting Information, Movie 7). However, the local application of calcium ions did not result in tubules for spherical vesicles, indicating that formation of the tubular protrusions for such vesicles is DOPS-dependent. To further investigate the role of calcium in the formation of tubular protrusions, another divalent cation, Mg2+, was tested. We observed that the spherical GUVs (both SPE/DOPS and SPE) exposed to 5 mM of Mg2+ did not produce tubular protrusions. In contrast, for undulating vesicles, the tubular protrusions were similar to those generated by Ca2+ ions, Figure 2E, Supporting Information, Movies 5 and 8. This indicates that the ion size plays a role in the force required to bend the membrane. The role of the exposure gradient was probed by conducting experiments where Ca2+ or Mg2+ was introduced to the vesicle solution in the bulk preparation, to a final concentration of 5 mM for each cation. In these instances, tubular protrusions were not observed. The vesicles (both SPE/DOPS and SPE) instead ruptured resulting in lipid patches adhered to the glass

surface (Supporting Information, section S1B, Figure S1). However, it should be noted that if GUV rupture could be prevented, it would allow investigation of the formation of tubular protrusions at bulk exposure. An example of such an experiment would include placing the MLV−GUV vesicle system on a microfabricated pillar, thereby limiting the adhesion of the lipid membrane to a surface. The effect of valency on the formation of the tubular protrusions was also probed, exposing the vesicles to monovalent cations (Na+, K+). No observable tubulation was present at concentrations of 7.5 or 75 mM of Na+ and 5 mM of K+, thus excluding the role of osmotic stress. The influence of hydrodynamic flow was excluded by performing control experiments where the micropipet was filled with pure HEPES buffer, which resulted in no observable effect on the membrane. Intrigued by the fact that tubular protrusions were induced in lipid vesicles by the local action of cations, we explored the possibility to translocate these tubules within the membrane by controlling the Ca2+ ion concentration profile at the membrane surface. By carefully translating the micropipet around the spherical SPE/DOPS vesicle surface (taking care to maintain the same distance between the vesicle and the pipet tip), we were able to guide the tubular protrusions in a contactless manner toward the pipet tip, that is, in the direction of increasing Ca2+ concentration (Figure 3, Supporting Information, Movie 6). Such a translational movement of the formed tubular protrusions is likely accompanied by the formation of new tubules at the new location of the calcium ion source, as we demonstrated in Supporting Information, section S1C, Figure S2. The observed controlled translation of tubular protrusions can have important implications for further understanding the mechanisms of lipid mobility and distribution within cell membranes, which can be guided by spatiotemporal calcium ion oscillations in the vicinity of the membrane. The tubular protrusions that form within the vesicles are reminiscent of the ones that are produced in model membranes upon localized proton flow,24 protein−protein crowding effect,46 binding of bacterial toxic proteins, for example, Shiga and cholera toxin47,48 or pore-forming protein perforin49 (at homogeneous concentrations). Interestingly, a generic mechanism is often shared among these phenomena, that is, a small increment of local curvature is induced by the inclusions (proteins)50,51 and the curvature inducing inclusions cluster and colocalize.52,53 This leads to a buildup of spontaneous D

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir

Figure 4. Molecular dynamics simulations of DOPC/DOPS bilayers under Ca2+ and Na+ ions gradients. (A, B) Snapshots for the 20%/Na+ and 20%/Ca2+ systems, respectively. The two bilayers create two solvent compartments to impose different ionic concentrations on each side of a bilayer. (Blue sphere) Na+ ions, (green sphere) Ca2+ ions, and lipids are shown in stick representation. Water particles and Cl− ions are not shown for clarity. (C) Averaged pressure profiles for each system. Asymmetry in the profile indicates spontaneous curvature and is more pronounced for systems containing Ca2+, as indicated in the inset. (D) Total charge density of the systems. The negative surface charge density of one side of the bilayer is effectively neutralized for systems containing Ca2+. (E) Charge densities of Na+ and Ca2+ for each system (Note, Ca2+ exists in only two systems). Notice, panels C, D, and E present the results for only one bilayer in the system. Ca2+ ions are localized on the surface of the bilayer, effectively neutralizing the surface charge density of the bilayer.

away from the Ca2+ ions; consistent with the experimental results (Figure 4C). However, for the 20%/Na+ system, the curvature was solely due to thermally driven deformation of the membrane surface, fluctuating around zero. This is also consistent with the experimental findings. Upon increase of the DOPS concentration to 40%, the radius of curvature of 40%/Ca2+ system decreased to 22 ± 4 nm (corresponding to 255 pN/μm spontaneous tension), while for the 40%/Na+ system, no spontaneous curvature was observed (Figure 4C). Ca2+ tends to bind to the surface of the bilayer and effectively neutralize the negative charge density of the bilayer (see Figure 4D, E). The adsorbed charge between two monolayers is 0.61 e/nm2 and 0.67 e/nm2 for 20%/Ca2+ and 40%/Ca2+ systems respectively, while this value is around 0.25 e/nm2 for both Na+ systems (Supporting Information, section S2C, Table S2). The resulting large asymmetry in the Ca2+ systems causes spontaneous membrane bending.55 Although some Na+ ions also bind to the membrane surface (Figure 4E), they only make a small reduction in the bilayer surface charge density (Figure 4D). Therefore, Na+ fails to induce large membrane curvature. Nevertheless, we should notice that this phenomenon described here is generic and is not limited only to the case Ca2+ and Na+. In general, multivalent counterions tend to bind strongly to the surface56 of the charged membrane and can neutralize the surface charge density,57 while the distribution of the monovalent counterions is more extended around charged surfaces (large Gouy−Chapman length).58 Our results suggest

tension (σs = 2κCm 2 , where Cm is the spontaneous mean curvature) that subsequently drives the formation of membrane tubular protrusions. Therefore, a sufficiently large domain with a small spontaneous curvature is enough for tubulation formation.54 In order to understand how Ca2+ binding induces spontaneous curvature to a charged bilayer, we have performed coarse grained molecular dynamics simulations of bilayers containing DOPC and two different concentrations of DOPS lipids (20% or 40%). These model systems contain two identical and symmetric bilayers to impose different ionic concentrations on each side of a bilayer (Figure 4A,B). The negative charge from the DOPS lipids is countered by adding Na+ evenly between the two solvent compartments. Additionally, in the middle solvent compartment, we insert Ca2+ and Cl− in one system (we refer to this system as 20%/Ca2+) and Na+ and Cl− in another (referred to as 20%/Na+). Spontaneous curvature of the bilayers was found by calculating the first moment of the pressure profile (see Simulation Methods) that is nonzero only for asymmetric profiles (Figure 4C). Figure 4C clearly shows that for all systems containing Ca2+ the pressure profile is asymmetric while for Na+ systems the profile is relatively symmetric. The results show that for the Ca2+ system with 20% DOPS concentration (20%/Ca2+), a local membrane curvature is induced with a curvature radius (inverse of the spontaneous mean curvature, 1/Cm) of 52 ± 12 nm (corresponding to 45 pN/μm spontaneous tension) directed E

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir

fluctuations of intracellular calcium ion concentration can vary by orders of magnitude, these concentration bursts may be sufficient to generate spontaneous curvature and lead to membrane deformations. Calcium ion variations may also directly regulate lipid mobility and lateral membrane reorganization, whereby membrane components sense the Ca2+ profile within the microenvironment. Biological cells contain a wide variety of charged molecular species, which may screen or magnify the Ca2+ ion effect, the contribution of each could be dissected by using our approach.

that the spontaneous curvature induced by divalent ion adsorption is directly related to the amount of adsorbed ions (Figure 4D,E). Therefore, the larger ions that have lower adsorption19 will induce a smaller spontaneous curvature and subsequently a lower spontaneous tension. This explains the different behaviors of Mg2+ and Ca2+ in the formation of the tubular protrusions. Notice, if there is no asymmetry in the Ca2+ exposure, the binding of Ca2+ ions induces lipid packing, instead of membrane bending (Supporting Information, section S2D, Table S3). Additionally, we ranked the effect of Ca2+, Mg2+, and Na+ ions according to their ability to generate tubular protrusions in the SPE/DOPS vesicles, such that Ca2+ > Mg2+ > Na+, where no membrane tubulation by Na+ was observed. We compared this ranking to the Hofmeister series, that is, a series in which the ions are empirically ranked with regards to the effect of these ions on, for example, the solubility of proteins in aqueous solution. According to the Hofmeister series, these cations are ranked as Na+ > Mg2+ > Ca2+.59 Interestingly, in our study on formation of tubular protrusions, the ranking is pointing in the direction of the reverse Hofmeister series. Screening of other cations according to their ability to induce formation of tubular protrusions would allow us to address the complete Hofmeister series.

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.langmuir.7b01461. Images of lipid vesicles exposed to Ca2+ and Mg2+, details on the microscopy image adjustments, and simulation details (PDF) Example of a SPE/DOPS giant lipid vesicle with spherical morphology (AVI) Example of a SPE/DOPS giant lipid vesicle with undulating morphology (AVI) Spherical SPE/DOPS giant lipid vesicle exposed to a Ca2+ gradient (5 mM), showing formation and growth of tubular protrusions directed inward in the GUV (AVI) Undulating SPE/DOPS giant lipid vesicle exposed to Ca2+ gradient (1 mM) (AVI) Undulating SPE/DOPS giant lipid vesicle exposed to Mg2+gradient (5 mM), showing formation of long and sparse tubular protrusions upon local Mg2+ stimulation (AVI) Formation and directed movement of tubular protrusions along the GUV surface in SPE/DOPS giant lipid vesicle exposed to Ca2+ ions (5 mM) (AVI) Undulating SPE giant lipid vesicle exposed to Ca2+ gradient (5 mM), showing formation of long and sparse tubular protrusions at the site of Ca2+ injection (AVI) Undulating SPE giant lipid vesicle exposed to Mg2+ gradient (5 mM), showing formation of long and sparse tubular protrusions at the site of Mg2+injection (AVI)

■

SUMMARY AND CONCLUSIONS Our results demonstrate that application of calcium ions to the membrane of giant lipid vesicles triggers membrane bending and tubulation, which can only be observed under conditions of local membrane stimulation. In particular, we show that the localized binding of Ca2+ to the outer membrane of negatively charged giant lipid vesicles generates inward spontaneous curvature sufficient to induce membrane bending, leading to the formation of tubular protrusions directed inside the vesicles. Lipid vesicles of two different compositions (SPE/DOPS and SPE) and two different morphological regimes (spherical and thermally undulating GUVs) were tested. We observed that within vesicles having spherical morphology, the tubular protrusions can only be observed in DOPS-enriched membranes. However, for undulating vesicles, the membrane protrusions were formed for both lipid mixtures, as well as in response to another divalent cation, Mg2+. The observed inward spontaneous curvature and membrane bending can best be described as a response of the bilayer to a large reduction in the surface charge density of one membrane leaflet upon binding of divalent cations. A charged bilayer has a strong repulsive coulomb interaction (a positive Maxwell stress) in both monolayers. Divalent cations bind strongly to the surface of the charged bilayer and neutralize the surface charge density of the bound monolayer. This asymmetry in the stress causes the membrane to bend. Molecular dynamics simulations, examining the effect of Ca2+ binding, demonstrated that sufficient force can be generated to induce spontaneous curvature, bending the membrane, leading to the formation of tubular protrusions, which agreed well with our experimental findings. The demonstrated lipid membrane behavior upon local application of Ca2+ offers unique insights into possible mechanisms of cell membrane remodeling, including that ion−lipid interaction alone can lead to the membrane deformations and generate highly curved tubular protrusions. These curved membrane regions are of fundamental importance in many cellular processes such as cell division, endo- and exocytosis, and trafficking of proteins and lipids between cellular compartments. Considering that cytosolic

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Baharan Ali Doosti: 0000-0001-9149-0751 Weria Pezeshkian: 0000-0001-5509-0996 Dennis S. Bruhn: 0000-0003-1362-4747 Tatsiana Lobovkina: 0000-0002-8923-1328 Author Contributions

B.A.D. and T.L. conceived the study. B.A.D., G.D.M.J., and T.L. participated in the experimental study. W.P., D.S.B., J.H.I., and H.K. performed the theoretical study. B.A.D., W.P., and D.S.B. are the major contributors to the experimental (B.A.D.) and theoretical work (W.P., D.S.B.). The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes

The authors declare no competing financial interest. F

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir

■

(19) Martín-Molina, A.; Rodríguez-Beas, C.; Faraudo, J. Effect of Calcium and Magnesium on Phosphatidylserine Membranes: Experiments and All-Atomic Simulations. Biophys. J. 2012, 102, 2095−2103. (20) Leckband, D. E.; Helm, C. A.; Israelachvili, J. Role of Calcium in the Adhesion and Fusion of Bilayers. Biochemistry 1993, 32, 1127− 1140. (21) Wilschut, J.; Duezguenes, N.; Papahadjopoulos, D. Calcium/ Magnesium Specificity in Membrane Fusion: Kinetics of Aggregation and Fusion of Phosphatidylserine Vesicles and the Role of Bilayer Curvature. Biochemistry 1981, 20, 3126−3133. (22) Simunovic, M.; Lee, K. Y. C.; Bassereau, P. Celebrating Soft Matter’s 10th Anniversary: Screening of the Calcium-Induced Spontaneous Curvature of Lipid Membranes. Soft Matter 2015, 11, 5030−5036. (23) Yaghmur, A.; Sartori, B.; Rappolt, M. The Role of Calcium in Membrane Condensation and Spontaneous Curvature Variations in Model Lipidic Systems. Phys. Chem. Chem. Phys. 2011, 13, 3115−3125. (24) Khalifat, N.; Puff, N.; Bonneau, S.; Fournier, J.-B.; Angelova, M. I. Membrane Deformation under Local pH Gradient: Mimicking Mitochondrial Cristae Dynamics. Biophys. J. 2008, 95, 4924−4933. (25) Lobovkina, T.; Gözen, I.; Erkan, Y.; Olofsson, J.; Weber, S. G.; Orwar, O. Protrusive Growth and Periodic Contractile Motion in Surface-Adhered Vesicles Induced by Ca2+-Gradients. Soft Matter 2010, 6, 268−272. (26) Kay, J. G.; Koivusalo, M.; Ma, X.; Wohland, T.; Grinstein, S. Phosphatidylserine Dynamics in Cellular Membranes. Mol. Biol. Cell 2012, 23, 2198−2212. (27) Vernier, P. T.; Ziegler, M. J.; Sun, Y.; Chang, W. V.; Gundersen, M. A.; Tieleman, D. P. Nanopore Formation and Phosphatidylserine Externalization in a Phospholipid Bilayer at High Transmembrane Potential. J. Am. Chem. Soc. 2006, 128, 6288−6289. (28) Criado, M.; Keller, B. U. A Membrane Fusion Strategy for Single-Channel Recordings of Membranes Usually Non-Accessible to Patch-Clamp Pipette Electrodes. FEBS Lett. 1987, 224, 172−176. (29) Karlsson, M.; Nolkrantz, K.; Davidson, M. J.; Stromberg, A.; Ryttsen, F.; Akerman, B.; Orwar, O. Electroinjection of Colloid Particles and Biopolymers into Single Unilamellar Liposomes and Cells for Bioanalytical Applications. Anal. Chem. 2000, 72, 5857−5862. (30) Berendsen, H. J. C.; van der Spoel, D.; van Drunen, R. GROMACS: A Message-Passing Parallel Molecular Dynamics Implementation. Comput. Phys. Commun. 1995, 91, 43−56. (31) Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435−447. (32) Pronk, S.; Páll, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M. R.; Smith, J. C.; Kasson, P. M.; van der Spoel, D.; et al. GROMACS 4.5: A High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 2013, 29, 845−854. (33) Marrink, S. J.; de Vries, A. H.; Mark, A. E. Coarse Grained Model for Semiquantitative Lipid Simulations. J. Phys. Chem. B 2004, 108, 750−760. (34) Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; de Vries, A. H. The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. B 2007, 111, 7812−7824. (35) Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. (36) Rahman, M.; Parrinello, M. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182−7190. (37) Safran, S. A. Statistical Thermodynamics of Surfaces, Interfaces, and Membranes; Addison-Wesley: Reading, MA, 1994. (38) Rózẏ cki, B.; Lipowsky, R. Spontaneous Curvature of Bilayer Membranes from Molecular Simulations: Asymmetric Lipid Densities and Asymmetric Adsorption. J. Chem. Phys. 2015, 142, 054101. (39) Brochard, F.; Lennon, J. F. Frequency Spectrum of the Flicker Phenomenon in Erythrocytes. J. Phys. 1975, 36, 1035−1047.

ACKNOWLEDGMENTS The authors acknowledge A-S. Cans for informative discussions. We also thank the Swedish Research Council (Grant 21220109) and Area of Advance Nano at Chalmers University of Technology (Grant 21420010) for funding. Simulations for the work described in this paper were performed using resources provided by the DeiC National HPC Centre, University of Southern Denmark.

■

REFERENCES

(1) Monteith, G. R.; McAndrew, D.; Faddy, H. M.; RobertsThomson, S. J. Calcium and Cancer: Targeting Ca2+ Transport. Nat. Rev. Cancer 2007, 7, 519−530. (2) Berridge, M. J.; Bootman, M. D.; Lipp, P. Calcium - a Life and Death Signal. Nature 1998, 395, 645−648. (3) Prevarskaya, N.; Skryma, R.; Shuba, Y. Calcium in Tumour Metastasis: New Roles for Known Actors. Nat. Rev. Cancer 2011, 11, 609−618. (4) Wollman, R.; Meyer, T. Coordinated Oscillations in Cortical Actin and Ca2+ Correlate with Cycles of Vesicle Secretion. Nat. Cell Biol. 2012, 14, 1261−1269. (5) Ross, W. N. Understanding Calcium Waves and Sparks in Central Neurons. Nat. Rev. Neurosci. 2012, 13, 157−168. (6) Birge, R. B.; Boeltz, S.; Kumar, S.; Carlson, J.; Wanderley, J.; Calianese, D.; Barcinski, M.; Brekken, R. A.; Huang, X.; et al. Phosphatidylserine Is a Global Immunosuppressive Signal in Efferocytosis, Infectious Disease, and Cancer. Cell Death Differ. 2016, 23, 962−978. (7) Roux, M.; Bloom, M. Ca2 + , Mg2 + , Li + , Na + , and K + Distributions in the Headgroup Region of Binary Membranes of Phosphatidylcholine and Phosphatidylserine As Seen by Deuterium. Biochemistry 1990, 29, 7077−7089. (8) Waldbillig, R. C.; Robertson, J. D.; McIntosh, T. J. Images of Divalent Cations in Unstained Symmetric and Asymmetric Lipid Bilayers. Biochim. Biophys. Acta, Biomembr. 1976, 448, 1−14. (9) Tarafdar, P. K.; Chakraborty, H.; Dennison, S. M.; Lentz, B. R. Phosphatidylserine Inhibits and Calcium Promotes Model Membrane Fusion. Biophys. J. 2012, 103, 1880−1889. (10) Melcrová, A.; Pokorna, S.; Pullanchery, S.; Kohagen, M.; et al. The Complex Nature of Calcium Cation Interactions with Phospholipid Bilayers. Sci. Rep. 2016, 6, 38035. (11) Pedersen, U. R.; Leidy, C.; Westh, P.; Peters, G. H. The Effect of Calcium on the Properties of Charged Phospholipid Bilayers. Biochim. Biophys. Acta, Biomembr. 2006, 1758, 573−582. (12) Bruhn, D. S.; Lomholt, M. A.; Khandelia, H. Quantifying the Relationship between Curvature and Electric Potential in Lipid Bilayers. J. Phys. Chem. B 2016, 120, 4812−4817. (13) Sinn, C. G.; Antonietti, M.; Dimova, R. Binding of Calcium to Phosphatidylcholine−phosphatidylserine Membranes. Colloids Surf., A 2006, 282−283, 410−419. (14) Binder, H.; Zschörnig, O. The Effect of Metal Cations on the Phase Behavior and Hydration Characteristics of Phospholipid Membranes. Chem. Phys. Lipids 2002, 115, 39−61. (15) Boettcher, J. M.; Davis-Harrison, R. L.; Clay, M. C.; Nieuwkoop, A. J.; Ohkubo, Y. Z.; Tajkhorshid, E.; Morrissey, J. H.; Rienstra, C. M. Atomic View of Calcium-Induced Clustering of Phosphatidylserine in Mixed Lipid Bilayers. Biochemistry 2011, 50, 2264−2273. (16) Porasso, R. D.; López Cascales, J. J. Study of the Effect of Na+ and Ca2+ Ion Concentration on the Structure of an Asymmetric DPPC/DPPC + DPPS Lipid Bilayer by Molecular Dynamics Simulation. Colloids Surf., B 2009, 73, 42−50. (17) Dluhy, R. A.; Cameron, D. G.; Mantsch, H. H.; Mendelsohn, R. Fourier Transform Infrared Spectroscopic Studies of the Effect of Calcium Ions on Phosphatidylserine. Biochemistry 1983, 22, 6318− 6325. (18) Ohki, S.; Zschörnig, O. Ion-Induced Fusion of Phosphatidic Acid Vesicles and Correlation between Surface Hydrophobicity and Membrane Fusion. Chem. Phys. Lipids 1993, 65, 193−204. G

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX

Article

Langmuir (40) Loubet, B.; Lomholt, M. A.; Khandelia, H. Tension Moderation and Fluctuation Spectrum in Simulated Lipid Membranes under an Applied Electric Potential. J. Chem. Phys. 2013, 139, 164902. (41) Karlsson, M.; Sott, K.; Cans, A.-S.; Karlsson, A.; Karlsson, R.; Orwar, O. Micropipet-Assisted Formation of Microscopic Networks of Unilamellar Lipid Bilayer Nanotubes and Containers. Langmuir 2001, 17, 6754−6758. (42) Sott, K.; Karlsson, M.; Pihl, J.; Hurtig, J.; Lobovkina, T.; Orwar, O. Micropipet Writing Technique for Production of Two-Dimensional Lipid Bilayer Nanotube-Vesicle Networks on Functionalized and Patterned Surfaces. Langmuir 2003, 19, 3904−3910. (43) Vanegas, J. M.; Torres-Sánchez, A.; Arroyo, M. Importance of Force Decomposition for Local Stress Calculations in Biomembrane Molecular Simulations. J. Chem. Theory Comput. 2014, 10, 691−702. (44) Ollila, O. H. S.; Risselada, H. J.; Louhivuori, M.; Lindahl, E.; Vattulainen, I.; Marrink, S. J. 3D Pressure Field in Lipid Membranes and Membrane-Protein Complexes. Phys. Rev. Lett. 2009, 78101, 1−4. (45) Wegrzyn, I.; Jeffries, G. D. M.; Nagel, B.; Katterle, M.; Gerrard, S. R.; Brown, T.; Orwar, O.; Jesorka, A. Membrane Protrusion Coarsening and Nanotubulation within Giant Unilamellar Vesicles. J. Am. Chem. Soc. 2011, 133, 18046−18049. (46) Stachowiak, J. C.; Schmid, E. M.; Ryan, C. J.; Ann, H. S.; Sasaki, D. Y.; Sherman, M. B.; Geissler, P. L.; Fletcher, D. a.; Hayden, C. C. Membrane Bending by Protein-Protein Crowding. Nat. Cell Biol. 2012, 14, 944−949. (47) Römer, W.; Berland, L.; Chambon, V.; Gaus, K.; Windschiegl, B.; Tenza, D.; Aly, M. R. E.; Fraisier, V.; Florent, J.-C.; Perrais, D.; et al. Shiga Toxin Induces Tubular Membrane Invaginations for Its Uptake into Cells. Nature 2007, 450, 670−675. (48) Ewers, H.; Römer, W.; Smith, A. E.; Bacia, K.; Dmitrieff, S.; Chai, W.; Mancini, R.; Kartenbeck, J.; Chambon, V.; Berland, L.; et al. GM1 Structure Determines SV40-Induced Membrane Invagination and Infection. Nat. Cell Biol. 2010, 12, 11−20. (49) Praper, T.; Sonnen, A. F.; Andrighetti, A. O.; Viero, G.; Morris, K. J.; Volpi, E.; Lunelli, L.; Serra, M. D.; Froelich, C. J.; Gilbert, R. J. C.; et al. Perforin Activity at Membranes Leads to Invaginations and Vesicle Formation. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 21016− 21021. (50) Pezeshkian, W.; Hansen, A. G.; Johannes, L.; Khandelia, H.; Shillcock, J. C.; Kumar, P. B. S.; Ipsen, J. H. Membrane Invagination Induced by Shiga Toxin B-Subunit: From Molecular Structure to Tube Formation. Soft Matter 2016, 12, 5164−5171. (51) Jarsch, I. K.; Daste, F.; Gallop, J. L. Membrane Curvature in Cell Biology: An Integration of Molecular Mechanisms. J. Cell Biol. 2016, 214, 375−387. (52) Pezeshkian, W.; Gao, H.; Arumugam, S.; Becken, U.; Bassereau, P.; Florent, J.-C.; Ipsen, J. H.; Johannes, L.; Shillcock, J. C. Mechanism of Shiga Toxin Clustering on Membranes. ACS Nano 2017, 11, 314− 324. (53) Liu, Y.; Agudo-Canalejo, J.; Grafmueller, A.; Dimova, R.; Lipowsky, R. Patterns of Flexible Nanotubes Formed by LiquidOrdered and Liquid-Disordered Membranes. ACS Nano 2016, 10, 463−474. (54) Lipowsky, R. Spontaneous Tubulation of Membranes and Vesicles Reveals Membrane Tension Generated by Spontaneous Curvature. Faraday Discuss. 2013, 161, 305−331. (55) Chou, T.; Jarić, M. V.; Siggia, E. D. Electrostatics of Lipid Bilayer Bending. Biophys. J. 1997, 72, 2042−2055. (56) Martín-Molina, A.; Rodríguez-Beas, C.; Faraudo, J. Charge Reversal in Anionic Liposomes: Experimental Demonstration and Molecular Origin. Phys. Rev. Lett. 2010, 104, 168103. (57) Pezeshkian, W.; Nikoofard, N.; Norouzi, D.; Mohammad-Rafiee, F.; Fazli, H. Distribution of Counterions and Interaction between Two Similarly Charged Dielectric Slabs: Roles of Charge Discreteness and Dielectric Inhomogeneity. Phys. Rev. Lett. E 2012, 85, 061925. (58) Naji, A.; Jungblut, S.; Netz, R. R.; Moreira, A. G. Electrostatic Interactions in Strongly Coupled Soft Matter. Phys. A 2005, 352, 131− 170.

(59) Okur, H. I.; Hladílková, J.; Rembert, K. B.; Cho, Y.; Heyda, J.; Dzubiella, J.; Cremer, P. S.; Jungwirth, P. Beyond the Hofmeister Series: Ion-Specific Effects on Proteins and Their Biological Functions. J. Phys. Chem. B 2017, 121, 1997−2014.

H

DOI: 10.1021/acs.langmuir.7b01461 Langmuir XXXX, XXX, XXX−XXX