Next Article in Journal
Short Review on Predicting Fouling in RO Desalination
Previous Article in Journal
Characterizing the Performance of Gas-Permeable Membranes as an Ammonia Recovery Strategy from Anaerobically Digested Dairy Manure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Revealing the Effects of Nanoscale Membrane Curvature on Lipid Mobility

by
Abir Maarouf Kabbani
,
Xinxin Woodward
and
Christopher V. Kelly
*
Department of Physics and Astronomy, Wayne State University, Detroit, MI 48201, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Membranes 2017, 7(4), 60; https://doi.org/10.3390/membranes7040060
Submission received: 12 September 2017 / Revised: 29 September 2017 / Accepted: 12 October 2017 / Published: 18 October 2017
(This article belongs to the Section Membrane Surfaces and Interfaces)

Abstract

:
Recent advances in nanoengineering and super-resolution microscopy have enabled new capabilities for creating and observing membrane curvature. However, the effects of curvature on single-lipid diffusion have yet to be revealed. The simulations presented here describe the capabilities of varying experimental methods for revealing the effects of nanoscale curvature on single-molecule mobility. Traditionally, lipid mobility is revealed through fluorescence recovery after photobleaching (FRAP), fluorescence correlation spectroscopy (FCS), and single particle tracking (SPT). However, these techniques vary greatly in their ability to detect the effects of nanoscale curvature on lipid behavior. Traditionally, FRAP and FCS depend on diffraction-limited illumination and detection. A simulation of FRAP shows minimal effects on lipids diffusion due to a 50 nm radius membrane bud. Throughout the stages of the budding process, FRAP detected minimal changes in lipid recovery time due to the curvature versus flat membrane. Simulated FCS demonstrated small effects due to a 50 nm radius membrane bud that was more apparent with curvature-dependent lipid mobility changes. However, SPT achieves a sub-diffraction-limited resolution of membrane budding and lipid mobility through the identification of the single-lipid positions with ≤15 nm spatial and ≤20 ms temporal resolution. By mapping the single-lipid step lengths to locations on the membrane, the effects of membrane topography and curvature could be correlated to the effective membrane viscosity. Single-fluorophore localization techniques, such SPT, can detect membrane curvature and its effects on lipid behavior. These simulations and discussion provide a guideline for optimizing the experimental procedures in revealing the effects of curvature on lipid mobility and effective local membrane viscosity.

1. Introduction

The shape of biological membranes is precisely controlled for diverse, essential, cellular processes, such as regulating organelle morphology, exocytosis/endocytosis, pathogen vulnerability/protection, and effective therapeutic targeting [1,2,3]. Accordingly, the dysregulation of membrane curvature is broadly implicated in cardiovascular disease, viral infections, cancer, Alzheimer's disease, Huntington disease, diabetes, and other diseases [4,5,6]. Each of these processes requires the fusion or fission of <50 nm radius vesicles with otherwise near planar membranes via precise regulation of the local curvature-generating forces [7]. Cellular membrane shape regulation incorporates a wide variety of proteins that can bend the membrane. For example, Bin/Amphiphysin/Rvs (BAR) domains have an intrinsic molecular shape [8], clathrin proteins create a scaffold [9], and intrinsically disordered proteins apply steric repulsion to induce membrane curvature [10]. The underlying non-specific, lipid-based influences remain relatively unknown in complex cellular membranes, although the importance of some key lipids has been previously demonstrated [11,12,13,14].
Experimental studies on the effects of curvature on membrane behavior are becoming possible with nanoengineering and super-resolution microscopy. For example, naturally occurring plasma membrane tubules [15] and engineered membrane tubules [16] have demonstrated protein and lipid sorting dependent on the membrane curvature. Nanoscale tubules have been created with model membrane via microbead pulling [17], protein crowding [18], or molecular motor pulling [19]. Prior experimental attempts to reveal diffusion differences on membrane tubules of varying radii were complicated by coupled variations in membrane composition and tension. However, slower diffusion is consistently observed on membrane tubules of smaller radii [20,21,22], as expected theoretically [23,24,25,26].
By engineering curvature on a solid substrate, modeled or living membranes may assume the substrate topography if the substrate radius of curvature and the membrane-substrate adhesion are sufficiently large. These engineered buds represent a membrane shape similar to endocytic pits preceding vesiculation and the post-fusion state of exocytosis, although the engineered structures are static, while endocytosis and exocytosis require dynamic membrane changes. Nanoscale membrane buds have been formed over substrates patterned via electron-beam lithography [16,27,28,29,30] and polystyrene nanoparticles [31,32,33,34]. These studies have revealed the effects of membrane curvature on protein sorting [16,27,34], lipid phases [28,29,30], and single-lipid dynamics [31,32,33].
Substrate nanoengineering has enabled the creation of the same membrane topographies simulated in this manuscript. The limited studies to date regarding the effective membrane viscosity on a nanoscale buds have reported that curvature can have no effect [32] or can slow the lipid mobility to 4–10% of the planar system [33,34,35,36]; however, these reports vary widely in the applied experimental methods. A focus of this manuscript is to demonstrate how varying observational methods can yield varying results on the effects of nanoscale curvature. Further, the collection of membrane morphologies simulated here represent a subset of the diverse membrane shapes created during endocytosis and exocytosis.
Understanding the effects of curvature on lipid dynamics will require separately resolving the two leaflets and the variation across ≤100 nm diameter membrane buds. Within the membrane bud, there is positive and negative principal curvature and positive and negative Gaussian curvature. The molecular structure of the constituent lipids and proteins can lead to curvature-induced nanoscale molecular sorting with compositional variation laterally across the membrane or between leaflets [21,37]. No known experiments to date have been able to distinguish the effects of nanoscale positive versus negative curvature on the effective membrane viscosity. However, membrane curvature generally seems to have the net effect of increasing the local effective membrane viscosity and slowing of the lipid and protein mobility [20,21,33,34,35,36]. Accordingly, sub-diffraction-limited spatial resolution will be necessary to experimentally measure how membrane curvature influences lipid and protein dynamics.
Optical techniques are traditionally limited by diffraction to a spatial resolution of >200 nm. Nanoscopic optical methods such as direct stochastic optical reconstruction microscopy [(d)STORM], fluorescence photoactivated localization microscopy [(f)PALM], and stimulated depletion emission microscopy (STED) [38,39,40,41,42] have improved the resolution of optical microscopy to >10 nm. Many variations on these techniques have been adapted to yield fluorophore height or orientation [43,44,45,46,47,48,49,50]. In particular, polarized localization microscopy (PLM) was designed to reveal nanoscale, membrane curvature with sub-diffraction-limited resolution [33]. PLM combines single-molecule localization microscopy (SMLM) and polarized total internal reflection fluorescence microscopy (TIRFM) [11,33]. Polarized TIRFM is sensitive to membrane orientation by selectively coupling linearly polarized fluorescence excitation with lipidated indocarbocyanine dyes (i.e., DiI, DiO, DiD) that maintain their fluorescence dipole moment in the plane of the membrane [51]. Pointillist SMLM methods, such as (d)STORM, (f)PALM, and PLM provide raw data that can be interpreted for high-throughput single-particle tracking dependent on membrane curvature. Tracking individual fluorophores that stay on for multiple sequential frames enables the observation of single-molecule diffusion rates versus membrane topology. For example, DiI molecules diffuse on curved membranes at <10% of the speed at which they diffuse on flat membranes [33,35]. Analysis of single-molecule diffusion rates relative to membrane topology may reveal information on the local environment (i.e., lipid phases or molecular crowding) cause by membrane bending.
In this manuscript, we demonstrate the capabilities of various fluorescence techniques to reveal lipid dynamics relative to membrane curvature. We focus on the three most common methods of measuring lipid mobility: fluorescence recovery after photobleaching (FRAP), fluorescence correlation spectroscopy (FCS), and single-particle tracking (SPT). Through Monte Carlo simulations of Brownian diffusing lipids over membrane buds of varying heights, we demonstrate the ability of each of these techniques in revealing the presence of the membrane bud, the single-molecule dynamics on the bud, and the effects of curvature on lipid mobility. Our simulations demonstrate how FRAP was not sufficiently sensitive to reveal that a bud was present under any of our simulation conditions. FCS revealed the bud’s presence, but FCS is typically limited to diffraction-limited length scales. SPT, however, measured the effects of membrane topography change with and without curvature-induced alteration to lipid mobility on each part of the membrane bud. By mapping the single-lipid steps over space, buds of varying heights and membranes of laterally varying viscosity could be distinguished. Within these simulations, we consider the effects of lipid diffusion variations with membrane curvature could have on the collected data, but we do not advocate for any particular function of curvature dependence on the diffusion rates or distinguish between the different types of lipids. Through carefully chosen methods, SPT data can reveal spatial information across the sample with <20 nm resolution. Guidelines are provided for designing SPT experiments to optimize the resolution of membrane curvature and its effects on molecular mobility.

2. Methods

The diffusion of lipids through membrane buds was simulated and analyzed to mimic the expected experimental results that would be obtained by a variety of fluorescence-based methods. All of the simulations were performed with custom MATLAB (MathWorks, Inc., Natick, MA, USA) programming, which is available in the Supplemental Material. Membrane buds were modeled with a radius of curvature equal to 50 nm and varying heights above a surrounding planar membrane (hbud). The bud membrane was smoothly connected to the surrounding planar membrane with a radius of curvature equal to 20 nm along the principal plane radial from the bud center (Figure 1), as done previously [33,34]. hbud = 0 represents the case of a planar membrane with no bud protrusion. When hbud = 140 nm, the bud had detached from planar membrane such that there was no diffusion between the vesicle and the planar membrane, and the vesicle was assumed to not contribute to the observed lipid diffusions as if the newly formed endosome had quickly left the proximity of the plasma membrane. Simulations of diffusion on a vesicle disconnected from a surrounding SLB have been recently published [35]. This is a minimalistic model of membrane shape during endocytosis or exocytosis that qualitatively matches electron microscopy images of cell plasma membranes [3,16] and allows for the single parameter hbud to report the stage of progression as the membrane transitions between a plane and the highly bent moment of fission or fussion.
The trajectories of the individual lipids were simulated upon the budding topography via a Monte-Carlo method for a discrete set of randomly distributed points. The discrete points were created at a density of 4 points/nm2 across the bud top, the bud-to-planar membrane neck, and the surrounding planar membrane. At each time step, the lipid moved to one of the 110 ± 10 random points within 3 nm. This resulted in an average single step distance of 2 nm. To mimic the diffusion coefficient (D) of 1 μm2/s over many steps, each time step would correspond to 1.1 μs. Each trajectory of each lipid started 1 μm away from the bud center, then diffused randomly upon the simulated membrane until it was >1 μm away from the bud center. More than 105 different trajectories were simulated for each condition, and 1300 ± 100 of those trajectories made it onto the membrane bud for each hbud. Example trajectories over a bud of hbud = 120 nm are shown in Figure 1C.
To mimic curvature-induced slowing of lipid diffusion, the effective time per simulation step was changed to be 11 μs or 28 μs for each 2 nm step whenever the simulated lipid was on the bud to mimic DBud = 0.1 or 0.04 μm2/s, respectively. Simulations were performed with the varying values of DBud to reveal how the various observation methods would report curvature-induced lipid slowing.
The analyses performed in this manuscript were limited to the z-projection of the fluorescence signal into the imaging xy-plane. The fluorescence emission was assumed to have no z-dependence or polarization dependence. These two assumptions apply well to topographies that vary in thickness by less than the diffraction-limit and for fluorophores that tumble randomly through all of the orientations, respectively.

2.1. Mimicking Fluorescence Recovery after Photobleaching

Fluorescence recovery after photobleaching (FRAP) measures the recovery of a fluorescence signal after the fluorophores within a region of the sample were bleached. Here, two bleaching conditions were used and compared. The first bleaching method, “complete bleaching”, was identical to the processes of coming to equilibrium, as described above, in which all the fluorophores within 1 μm from the bud center were bleached. The second method, “Gaussian bleaching”, took an equilibrated distribution of fluorophores and imposed a bleaching probability upon the ith fluorophore equal to P G B ( i ) =   exp ( 2 r i 2 / w 2 ) with a w = 250 nm and ri equal to the lateral distance of the fluorophore from the bud center. In both cases, simulated trajectories were analyzed with a Gaussian illumination (I) centered on the membrane bud according to
I ( t ) = Σ i   exp ( 2 r i 2 / w 2 )
with an illumination width (w) set equal to 250 nm, as would be expected for a diffraction-limited illumination. Upon Gaussian bleaching, I(t) was reduced to 50 ± 5% of the equilibrated intensity.
New single-lipid trajectories were started at the perimeter of the 2 μm diameter system every 0.08 μs for an equilibrium density of 0.00146 lipids/nm2 or approximately 10 mol % fluorescent lipids. Since an area larger than the illumination spot was initially bleached with the complete bleaching method, only the I(t) after the 25% of the recovery had occurred was analyzed. The increasing I(t) was fit to
I F R A P   F i t ( t ) =   A ( 1 e t / τ F R A P )
The fitting variable A represents the steady-state magnitude of I and is proportional to the steady-state fluorophore density, the membrane area, the illumination intensity, and the fluorescence emission collection efficiency. The fitting variable τFRAP represents the characteristic FRAP recovery time. This fitting model is an approximation to the expected recovery for planar systems with uniform [52] or Gaussian illumination [53]; this approximation was used since the exact recovery shape over the non-planar membrane topography would with hbud, which would not be typically known a priori.

2.2. Mimicking Fluorescence Correlation Spectroscopy

Fluorescence correlation spectroscopy (FCS) examines the fluctuations in the steady state I versus t signal. This is performed by calculating the autocorrelation (G) as a function of lag time (τ) and finding the characteristic fluctuation time (τFCS). In these simulations, I(t) for FCS was calculated from the single-molecule trajectories through a Gaussian illumination profile according to Equation (1), as would be expected for typical confocal FCS. G was calculated from I(t) according to
G ( τ ) =   < δ I ( t ) δ I ( t τ ) > / < I ( t ) > 2
The angle brackets (< >) represent the average over t and δI(t) = I(t) − <I(t)>. The correlation time (τFCS) in I(t) was found by fitting G(τ) according to
G F i t ( τ ) = G 0 ( 1 ( τ / τ F C S ) 2 ) 1
as is expected for two-dimensional (2D) Brownian diffusion. The fitting variable, G0, is inversely proportional to the number of diffusers simultaneously observed and other experimental conditions not relevant in these simulations. With a membrane bud present, the autocorrelation is not expected to fit perfectly to Equation (4); however, the inherent averaging incorporated into an autocorrelation analysis makes finding minor populations difficult, and complex fitting functions are typically unwarranted [54]. For this analysis, Equation (1) was assumed to also be the spatial detection sensitivity, which is the standard approximation of the Airy point spread function in a diffraction-limited system [55] and especially accurate when a confocal collection pinhole is used with a diameter of one Airy unit.

2.3. Mimicking Single Particle Tracking

Single-particle tracking (SPT) includes identifying the center of each single-fluorophore image via computational analysis of a movie of sparse, dynamic fluorophores. From the motion of the single-molecules between sequential frames, single-molecules trajectories were observed. The single-steps lengths (s) observed in a region of the sample may be fit to a 2D Maxwell-Boltzmann or Rayleigh Distribution to determine the local D, where the probability distribution (P) of step lengths over a time step of Δt for a single Brownian diffuser in a uniform membrane is expected to be
P ( s ) = s 2 D Δ t e s 2 4 D Δ t
The observed single-molecule steps were grouped and fitted according to their location in the sample (x, y) or distance from the bud center (r) such that D could be measured at different locations in the sample. The ability to gain spatial resolution for variations in D across a sample is the key differences between a single-step length analysis used here versus the more traditional mean squared displacement (MSD) analysis, as discussed below. Optimization of MSD analysis has been recently described [56].
However, the single-step analysis in experimental systems is affected by the 2D localization uncertainty ( σ r 2 = σ x 2 + σ y 2 ) and the camera exposure duration (texp) to yield a systematic difference between the D found from fitting Equation (5) to experimental data (DFit) and the D that would be found from an idealized system (DReal), according to
D R e a l = ( D F i t σ r 2 2 Δ t ) / ( 1 t e x p 3 Δ t ) .
Frequently, Δt is equal to the time between adjacent frames (tframe) used in this analysis, however, any Δt that is a multiple of tframe are permitted. tframe is equal to the sum of texp, and the frame read time such that the inverse of 1/tframe equal the imaging frequency.
When the membrane is not parallel to the coverslip, then the z-component of the lipid diffusion within the membrane, results in a slowing of the lipid through the xy-plane. It was not possible to extract the in-plane diffusion rate from the observed diffusion through the xy-plane (Dxy), when both the membrane topography and the influence of curvature on membrane viscosity are unknown. In the below analysis, DReal was calculated under the approximation that the membrane was parallel to the coverslip, and this value was reported as Dxy to be explicit in representing that the membrane topography was contributing to the measurements.

2.4. Assumptions Used in the Simulations

Super-resolution methods may incorporate a light polarization and z-dependence in the localization probability versus membrane height or orientation. The localization probability can vary as a fluorophore diffuses away from the objective focal plane in the z-direction. However, this probability change occurs only when variations in z are greater than the diffraction limit are allowable. Since the nanoscale buds simulated here vary in z by <140 nm, no variation in the probability of fluorophore localization is expected. Polarization was important in FCS of vesicles containing DiD [35]; however, no differences in s- versus p-polarized illumination was observed in SPT of DiI [33]. The simulations presented here mimic randomly tumbling fluorophores with no polarization dependence.
Simulations performed here do not incorporate any differences in the lipid diffusion due to the membrane-substrate adhesion. The diffusion of transmembrane proteins and >100 nm diameter lipid domains can be greatly slowed in SLBs versus freestanding membranes; however, the diffusion of individual lipids are significantly less affected by the substrate [57]. We would expect to see significant substrate effects for fluorescent lipids with larger head-group labeling [58] or for systems that may form >100 nm phase-separated domains.

3. Results

3.1. Membrane Bending and FRAP

Experimental FRAP conditions, such as the laser-based illumination with complete or Gaussian bleaching and small observation areas, were mimicked in the simulations presented here. The presence of the membrane bud and the extra membrane area associated with it were not sufficient to consistently influence the FRAP recovery time. Upon incorporating a curvature-induced slowing of the lipid diffusion on the bud, where DBud was reduced to 0.1 or 0.04 μm2/s, while maintaining DPlane = 1 μm2/s, the variation between repeated simulations proved to provide insignificant differences in τFRAP versus hbud for both complete bleaching (Figure 2A,B) and Gaussian bleaching (Figure 2C,D). Further, these results were created from a perfectly centered Gaussian illumination on the bud center (Equation (1)). Upon slight deviations of the illumination from the bud center, the effects of the bud on τFRAP were reduced further. Within the range of DBud explored here, the identification of the bud and the curvature-induced changes in lipid diffusion were undetectable by FRAP.

3.2. Membrane Bending and FCS

The presence of a single nanoscale membrane bud could be detected with FCS with or without curvature-induced lipid slowing (Figure 3). With diffraction-limited illumination (Equation (1)) and imprecise centering of the observation spot over the membrane bud, simulated FCS detected the effects of the bud on the lipid dwell time. For example, with DPlane and DBud = 1 μm2/s, the addition of the extra membrane area for a bud of hbud = 100 results in a 60% slowing of the intensity fluctuations; τFCS = 12.3 ± 0.6 ms and 20 ± 1 ms when hbud = 0 and 100 nm, respectively. This change is consistent with the 60% extra membrane area created by the bud when hbud = 100 nm versus 0 nm (Figure 1B).
When the membrane bud incorporated slower lipid diffusion, the effects of the membrane bud on the FCS result were more pronounced. For example, a 1700% slowing τFCS was observed between hbud = 0 and 100 nm when DBud = 0.1 μm2/s and DPlane = 1 μm2/s, with τFCS = 12.3 ± 0.6 ms and 210 ± 20 ms, respectively; the autocorrelations for varying hbud are shown for this condition in Figure 3A.
The effects of the membrane bud on the FCS results were apparent but reduced if the fluorescence illumination was not well centered over the membrane bud. For example, the effects on τFCS was decreased by an average 15% when the center of the membrane bud was 100 nm offset away from the center of the diffraction-limited FCS illumination spot, as calculated by the average ratio of the centered versus offset τFCS values for hbud from 5 nm to 135 nm (Figure 3B).

3.3. Membrane Bending and SPT

Mapping the locations of single molecule steps enables measuring Dxy at specific locations within the sample (Figure 4). When a membrane with consistent viscosity and consistent lipid mobility was simulated, the variations in the measured Dxy across the sample were due to the membrane topology change, and the measured shorter lipid step lengths on tilted membranes as the lipid trajectories are projected on the imaging xy-plane (Figure 4A). The effects of the membrane bud are enhanced when the curved membrane reduces the lipid mobility (Figure 4B,C). When the membrane bending occurs in a rotationally symmetrical way, as would be expected for an endocytic pit and the engineered curvature shown here, a radial averaging provides greater clarity in the effects of membrane budding on the SPT results (Figure 5).
In an SPT experiment, a balance between acquisition frame rate, localization precision, laser power, and on fluorophore density is required to achieve meaningful data. The simulations performed here demonstrate that a faster frame rate at the cost of worse localization precision yields.
A finer resolution of the lipid mobility for most experimental conditions. The most common experimental SMLM and SPT parameters result in σr = 15 nm and tframe = 20 ms. The effects of the membrane bud for the membrane topography simulated here, assuming DPlane = DBud, resulted in a negligible change in Dxy versus r until hbud ≥ 60 nm (Figure 5A). The greatest loss of signal with these conditions came from the long distance the fluorophore moved between adjacent frames, and individual lipids sampling both the curved and planar portions of the membrane within a single step. For tframe = 20 ms and D = 1 μm2/s, the expected step length of a single fluorophore between adjacent frames ( s ¯ 1 ) would be (4Dtframe)1/2 = 280 nm. This is much farther than noise added by the single-fluorophore localization imprecision.
When the membrane curvature caused the local lipid diffusion coefficient to slow, the membrane bud was readily apparent via SPT with σr = 15 nm and tframe = 20 ms for any hbud (Figure 5A–C). SPT revealed Dxy for the top of the bud separate from the surrounding planar SLB; however, the diffusion around the bud neck was challenging to interpret when considering the combined effects of membrane tilt and data blur.
Significant benefits can be obtained from SPT results by employing faster imaging frame rates. Keeping all other parameters constant, the effects of decreasing tframe from 20 ms to 2 ms resulted in a 3x improvement in blur induced by diffusion between frames (Figure 5B,E). Additionally, the shorter tframe provided an improvement in data statistics for a given total acquisition time, this is mainly due to the increased number of independent steps observed per time. Even if decreasing tframe comes at the cost of increasing σr to 45 nm, there was still a clear improvement in the obtained results (Figure 5B,D). This relative importance of tframe and σr is further demonstrated by showing how marginal improvements in the resolution of Dxy were obtained if tframe = 2 ms and σr was decreased from 15 nm to 3 nm, which is feasible when fluorescent beads or metal nanoparticles are used for tracking rather than single fluorophores (Figure 5E,F).

4. Discussion

Resolving the nanoscale biophysical effects of membrane curvature remains experimentally challenging. Optical techniques such as PLM enable the detection of curvature with higher sensitivity and resolution than conventional optical techniques, while providing a biologically friendly imaging environment. However, decoupling the effects of varying membrane area, orientation, and curvature can be challenging when the primary data collected is the z-projection of sample topography changing on sub-diffraction-limited length scales. The focus of this manuscript is to provide a framework for protocols and capabilities of experimental measurements of lipid diffusion as influenced by nanoscale membrane curvature. The simulation results presented here guide the experimental procedures that will improve the understanding of membrane curvature.

4.1. Capabilities of FRAP, FCS, and SPT

The method of draping a supported lipid bilayer over nanoparticles has been used in prior experimental studies to reveal the influence of curvature on lipid dynamics and protein sorting. FRAP [30,31] and SPT [32,33] have been used to confirm the continuity of the membrane over the nanoparticles. FRAP was used to confirm membrane integrity by matching the final fluorescence emission intensity after bleaching and recovery to the before-bleaching intensity. SPT was used to confirm membrane integrity through the observation of single-particle trajectories connecting the bud to the surrounding SLB. However, FRAP was unable to reveal any difference in the recovery rate of the fluorescence due to the presence of curvature in experiments [31,32], as expected by the simulations performed here (Figure 2). The setup in these simulations matches most prior experimental approaches where a single membrane bud is being observed at a time. Should the bud density become high enough for multiple buds were present within the FRAP observation region, FRAP may prove to provide bud-dependent results.
FCS can detect the effects of the bud and curvature-induced slowing on lipid diffusion in these simulations with diffraction-limited illumination (Equation (1)) and imprecise centering of the observation spot over the membrane bud (Figure 3). However, the ability for FCS to detect membrane buds without curvature-induced lipid slowing was limited to mature buds (hbud > 55 nm) for a 20% change in τFCS. With curvature induced slowing (DBud ≤ DPlane/10), significant changes to τFCS are apparent as soon as the bud is formed. Additionally, it would be feasible to provide sub-diffraction-limited resolution of the bud’s location in the sample without a complementary signal (e.g., nanoparticle, atomic force microscopy, or clathrin co-localization) by finding where τFCS was most slowed when scanning the illumination beam over the sample. Without a complementary signal co-localized to the slowed τFCS, however, it would be difficult to confirm that τFCS slowing was caused to a membrane budding rather than another membrane defect, such as membrane-substrate interaction or sorting to domains of lipid phase separation.
FCS is a valuable technique for identifying the presence of a membrane bud and/or the magnitude of slowing induced by the curvature; however, the ability of FCS to reveal the lipid mobility on different parts of the bud are prohibited by the relatively large size of the diffraction-limited illumination (w = 250 nm) as compared to the bud radius (50 nm). Potentially, sub-diffraction-limited STED-FCS could yield a greater resolution of the lipid mobility on distinct parts of the bud as well as increase sensitivity to the membrane topography itself [59]; however, STED-FCS is technically challenging, expensive, and rare. Sub-bud resolution may be achievable with FCS by exploiting rotationally-limited diffusion of fluorophores with polarized excitation or emission, as has been done on lipid vesicles of varying size [35].
SPT is unique among these techniques in that it provides sub-diffraction-limited spatial resolution with the capability to reveal the locations on the bud that most affect the local lipid diffusion. By binning the single-lipid step size versus distance from the center of the membrane bud and fitting the resulting histogram of step lengths to the Rayleigh distribution, the effects of lipid topography can be revealed directly even without curvature-affected lipid mobility (Figure 4A and Figure 5A). The effect of the membrane bending was observed via Dxy in these simulations without any curvature-induced changes to D. Similarly, observation of Dxy across a sample could reveal previously unknown membrane topography if there were curvature induced change to lipid mobility were known. This is particularly shown in Figure 4C, in which the vertical edge of the membrane bud, the horizontal top of the bud, and the surrounding SLB are each individual identifiable.

4.2. Comparative Ease of Performing FRAP, FCS, and SPT

The decision of which technique to use requires considering both the type of information needed to be obtained and the associated technical challenges. FRAP, FCS, and SPT each provide benefits regarding the specific membrane processes that they reveal and the ease by which they are experimentally performed. FRAP is the easiest to these techniques in both the execution and analysis of the experiment. FRAP can be carried out on large observation regions with a conventional epifluorescence microscope through the opening and closing of a field diaphragm in the conjugate image plane. Alternatively, greater spatial resolution can be gained with diffraction-limited bleaching and illumination, as simulated here. Although these methods of performing FRAP can provide a coarse analysis of membrane integrity and lipid mobility, it can be difficult to achieve a high intensity of the fluorescence emission and small enough observation areas to provide precise measurements. Even if performed with a relatively weak and slow bleaching procedure, FRAP can reveal the fraction of the diffusers that are immobile and the average diffusion coefficient of the mobile diffusers in a large observation area. This is especially valuable for demonstrating the continuity of a model membrane.
The correlation of intensity versus time for FCS only reports the diffusers that move through the diffraction-limited observation spot over the ~30 s of data collection; FCS does not incorporate any information from immobile particles on a membrane. However, FCS is more sensitive to sub-populations of diffusers than FRAP and provides greater accuracy in the measured diffusion coefficients of the mobile diffusers. Commercial FCS setups require laser illumination, expensive detectors, hardware correlators, and software, but can provide analysis of the results in real-time. Homebuilt FCS setups may use high-frame rate EMCCD or sCMOS cameras and custom software for correlation calculation and fitting. Since FCS can reveal late-stage bud formation and the effects of bending on lipid mobility, it is feasible that future incorporations of FCS will be used to report the biophysical ramifications of membrane bending.
As shown above, SPT provides the best spatial resolution of membrane bending and the effects of bending on lipid mobility. However, SPT requires significantly more effort in data collection and analysis. SPT requires a high photon flux of the emission (>5000 photons/s) for precise single-fluorophore fitting (σr < 50 nm) with fast frame rates (≥50 Hz), which often requires oxygen-scavenging buffers to reduce fluorophore oxidization to provide more photons per fluorophore ‘on’ state and greater conversion from the fluorophore ‘off’ to ‘on’ states. The raw SPT data typically comprises ≥5000 individual camera images, from which the single-molecule locations are calculated. The locations are linked for trajectory analysis, MSD analysis, and/or single-step length analysis. A complicating factor of SPT is that the data quality and signal-to-noise ratio can vary between experiments, such that user confirmation is needed for the analysis of each experiment. Despite these experimental challenges, the resolution benefits of SPT commonly justify its implementation.

4.3. SPT without Long Trajectories

Typically, the single-particle trajectories are analyzed by calculating the MSD versus Δt such that 2D Brownian diffusion results in a linear relationship of
M S D ( Δ t ) =   4 D Δ t +   2 σ r 2 8 D R t f r a m e
The localization uncertainty and camera blur contribute to these last two terms of Equation (7), respectively. Camera blur is a spreading of the acquired image due to the motion of the subject during the finite single-frame acquisition time. Camera blur depends on the motion blur coefficient (R), such that R = texp/(6tframe) ≤ 1/6 for continuous exposures within each frame [60]. MSD analysis has the potential to provide a precise D for a single diffuser and prior analysis has optimized the experimental parameters for MSD analysis [56]. However, MSD versus Δt analysis spatially averages each trajectory by providing a single D for all the space explored during the single trajectory. For trajectories that include dozens of sequential localizations, this analysis may provide an average over many square microns of the sample, well beyond the extent of a membrane bud. Conventional MSD analysis is unable to provide the curvature-dependent diffusion of lipids when only a small fraction of each trajectory is on the membrane bud.

4.4. Effects of Frame Rate on SPT

Individual fluorophore localizations for SMLM or SPT are frequently made to ≤20 nm certainty with precisely optimized imaging buffers, illumination intensities, and detection optics. However, the spatial resolution of the mobility measurement is often determined by the distance between sequential localizations rather than the precision of single localizations. SMLM methods yield σr that is approximately proportional to the inverse of the square root of the number of photons collected per single fluorophore image (N), while N is proportional to texp [61]. Similarly, s ¯ 1   is proportional to t f r a m e . Assuming minimal image readout time, minimal background noise, and texptframe, s ¯ 1   and σr are approximately inversely related to each other, and an imaging frame rate could be found that results in s ¯ 1   σr, for which the best resolution of D across a sample could be achieved. This is demonstrated by comparing panels B and D in Figure 5 for which a substantial improvement in the resolution of D results from t f r a m e   being decreased by 90% (from 20 ms to 2 ms), s ¯ 1 decreasing by 70% (from 280 nm to 90nm), and σr increasing by 200% (from 15 nm to 45 nm). Presumably, decreasing tframe further would have resulted in even greater resolution benefit until s ¯ 1 equaled σr; however, experimental realities, such as image background noise, becomes significant at shorter tframe when N < 100 such that σr increases significantly faster than 1 / N   and there are no longer benefits of decreasing tframe.
Increasing the photon flux can improve both σr and s ¯ 1 . For example, single biomolecules labeled with metal nanoparticles have yielded more fluorescence emission to demonstrate membrane hop diffusion with σr = 17 nm and tframe = 0.03 ms [62]. Metal nanoparticles detected in a non-fluorescence, interferometric microscope have yielded σr = 1.7 nm and tframe = 1 ms [63]. However, these experiments with ≤1 ms frame rates depend on >20 nm diameter gold nanoparticle labels and are associated experimental uncertainties that are not present with single-fluorophore labels. The uncertainties associated specifically with gold nanoparticle labels include the reduced specificity of the number of lipids per nanoparticle, the effects of drag on the nanoparticle, the non-specific binding between the nanoparticle and the other membrane components, and the local heating that may be caused by the gold absorption of the illumination.

5. Conclusions

There are numerous challenges for observing the effects of membrane bending at physiological length scales. Diverse super-resolution optical techniques are providing resolution of nanoscale membrane topography; however, the dynamical effects of curvature remain largely unknown. The ability for super-resolution optical techniques such as PLM to reveal nanoscale membrane bending is expanding the experimental capabilities for membrane curvature detection. The capacity to engineer membrane bending through the creation of supported lipid bilayers draped over nanoengineered substrates allow for the experimental creation of membrane topographies that are analogous to endocytosis and exocytosis. With SPT, the spatial resolution of lipid mobility and membrane bending can be observed with a higher precision than detectable with FCS or FRAP. The fitting of the histogram of single-step sizes distribution enables the calculation of the lipid diffusion coefficients that are corrected for the localization uncertainty and the camera blur. Future experimental implementations of PLM and SPT will reveal the effects of membrane bending on the membrane viscosity and lipid mobility. Through asymmetric labeling of model membranes, the specific contribution of each bilayer leaflet will be determined in the nanoscale budding membrane. The sum of these results will contribute to the greater understanding of membrane biophysics and the mechanisms of cellular regulation of membrane topography.

Supplementary Materials

The custom MATLAB scripts used in this work are available online at www.mdpi.com/2077-0375/7/4/60/s1.

Acknowledgments

The authors thank the Wayne State University High Performance Computing Grid. A.M.K. was funded by Thomas C. Rumble Fellowship Award. Financial support was provided by Wayne State University laboratory startup funds and Richard J. Barber. This material is based upon work supported by the National Science Foundation under Grant No. DMR1652316.

Author Contributions

A.M.K., X.W. and C.V.K. conceived of and designed the simulations; A.M.K., X.W. and C.V.K. analyzed the data; A.M.K., X.W. and C.V.K. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. 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. [Google Scholar] [CrossRef] [PubMed]
  2. McMahon, H.T.; Boucrot, E. Membrane curvature at a glance. J. Cell Sci. 2015, 128, 1065–1070. [Google Scholar] [CrossRef] [PubMed]
  3. Bonifacino, J.S.; Glick, B.S. The mechanisms of vesicle budding and fusion. Cell 2004, 116, 153–166. [Google Scholar] [CrossRef]
  4. Chaibva, M.; Burke, K.A.; Legleiter, J. Curvature enhances binding and aggregation of huntingtin at lipid membranes. Biochemistry (Moscow) 2014, 53, 2355–2365. [Google Scholar] [CrossRef] [PubMed]
  5. Williams, T.L.; Serpell, L.C. Membrane and surface interactions of Alzheimer’s Aβ peptide—Insights into the mechanism of cytotoxicity. FEBS J. 2011, 278, 3905–3917. [Google Scholar] [CrossRef] [PubMed]
  6. Michel, V.; Bakovic, M. Lipid rafts in health and disease. Biol. Cell Auspices Eur. Cell Biol. Organ. 2007, 99, 129–140. [Google Scholar] [CrossRef] [PubMed]
  7. Stachowiak, J.C.; Brodsky, F.M.; Miller, E.A. A cost-benefit analysis of the physical mechanisms of membrane curvature. Nat. Cell Biol. 2013, 15, 1019–1027. [Google Scholar] [CrossRef] [PubMed]
  8. Mim, C.; Unger, V.M. Membrane curvature and its generation by BAR proteins. Trends Biochem. Sci. 2012, 37, 526–533. [Google Scholar] [CrossRef] [PubMed]
  9. Fotin, A.; Cheng, Y.; Sliz, P.; Grigorieff, N.; Harrison, S.C.; Kirchhausen, T.; Walz, T. Molecular model for a complete clathrin lattice from electron cryomicroscopy. Nature 2004, 432, 573–579. [Google Scholar] [CrossRef] [PubMed]
  10. Busch, D.J.; Houser, J.R.; Hayden, C.C.; Sherman, M.B.; Lafer, E.M.; Stachowiak, J.C. Intrinsically disordered proteins drive membrane curvature. Nat. Commun. 2015, 6, 7875. [Google Scholar] [CrossRef] [PubMed]
  11. Baumgart, T.; Capraro, B.R.; Zhu, C.; Das, S.L. Thermodynamics and Mechanics of Membrane Curvature Generation and Sensing by Proteins and Lipids. Annu. Rev. Phys. Chem. 2011, 62, 483–506. [Google Scholar] [CrossRef] [PubMed]
  12. Sadeghi, S.; Müller, M.; Vink, R.L.C. Raft formation in lipid bilayers coupled to curvature. Biophys. J. 2014, 107, 1591–1600. [Google Scholar] [CrossRef] [PubMed]
  13. Larsen, J.B.; Jensen, M.B.; Bhatia, V.K.; Pedersen, S.L.; Bjørnholm, T.; Iversen, L.; Uline, M.; Szleifer, I.; Jensen, K.J.; Hatzakis, N.S.; et al. Membrane curvature enables N-Ras lipid anchor sorting to liquid-ordered membrane phases. Nat. Chem. Biol. 2015, 11, 192–194. [Google Scholar] [CrossRef] [PubMed]
  14. Parthasarathy, R.; Yu, C.; Groves, J.T. Curvature-modulated phase separation in lipid bilayer membranes. Langmuir 2006, 22, 5095–5099. [Google Scholar] [CrossRef] [PubMed]
  15. Rosholm, K.R.; Leijnse, N.; Mantsiou, A.; Tkach, V.; Pedersen, S.L.; Wirth, V.F.; Oddershede, L.B.; Jensen, K.J.; Martinez, K.L.; Hatzakis, N.S.; et al. Membrane curvature regulates ligand-specific membrane sorting of GPCRs in living cells. Nat. Chem. Biol. 2017, 13, 724–729. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, W.; Hanson, L.; Lou, H.-Y.; Akamatsu, M.; Chowdary, P.D.; Santoro, F.; Marks, J.R.; Grassart, A.; Drubin, D.G.; Cui, Y.; et al. Nanoscale manipulation of membrane curvature for probing endocytosis in live cells. Nat. Nanotechnol. 2017, 12, 750–756. [Google Scholar] [CrossRef] [PubMed]
  17. Roux, A.; Cuvelier, D.; Nassoy, P.; Prost, J.; Bassereau, P.; Goud, B. Role of curvature and phase transition in lipid sorting and fission of membrane tubules. EMBO J. 2005, 24, 1537–1545. [Google Scholar] [CrossRef] [PubMed]
  18. Stachowiak, J.C.; Hayden, C.C.; Sasaki, D.Y. Steric confinement of proteins on lipid membranes can drive curvature and tubulation. Proc. Natl. Acad. Sci. USA 2010, 107, 7781–7786. [Google Scholar] [CrossRef] [PubMed]
  19. Roux, A.; Cappello, G.; Cartaud, J.; Prost, J.; Goud, B.; Bassereau, P. A minimal system allowing tubulation with molecular motors pulling on giant liposomes. Proc. Natl. Acad. Sci. USA 2002, 99, 5394–5399. [Google Scholar] [CrossRef] [PubMed]
  20. Domanov, Y.A.; Aimon, S.; Toombes, G.E.S.; Renner, M.; Quemeneur, F.; Triller, A.; Turner, M.S.; Bassereau, P. Mobility in geometrically confined membranes. Proc. Natl. Acad. Sci. USA 2011, 108, 12605–12610. [Google Scholar] [CrossRef] [PubMed]
  21. Zhu, C.; Das, S.L.; Baumgart, T. Nonlinear Sorting, Curvature Generation, and Crowding of Endophilin N-BAR on Tubular Membranes. Biophys. J. 2012, 102, 1837–1845. [Google Scholar] [CrossRef] [PubMed]
  22. Tian, A.; Baumgart, T. Sorting of Lipids and Proteins in Membrane Curvature Gradients. Biophys. J. 2009, 96, 2676–2688. [Google Scholar] [CrossRef] [PubMed]
  23. Daniels, D.R.; Turner, M.S. Diffusion on membrane tubes: A highly discriminatory test of the Saffman-Delbruck theory. Langmuir ACS J. Surf. Colloids 2007, 23, 6667–6670. [Google Scholar] [CrossRef] [PubMed]
  24. Klaus, C.J.S.; Raghunathan, K.; DiBenedetto, E.; Kenworthy, A.K. Analysis of diffusion in curved surfaces and its application to tubular membranes. Mol. Biol. Cell 2016, 27, 3937–3946. [Google Scholar] [CrossRef] [PubMed]
  25. Henle, M.L.; Levine, A.J. Hydrodynamics in curved membranes: The effect of geometry on particulate mobility. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2010, 81, 011905. [Google Scholar] [CrossRef] [PubMed]
  26. Renner, M.; Domanov, Y.; Sandrin, F.; Izeddin, I.; Bassereau, P.; Triller, A. Lateral Diffusion on Tubular Membranes: Quantification of Measurements Bias. PLoS ONE 2011, 6. [Google Scholar] [CrossRef] [PubMed]
  27. Xie, C.; Hanson, L.; Cui, Y.; Cui, B.X. Vertical nanopillars for highly localized fluorescence imaging. Proc. Natl. Acad. Sci. USA 2011, 108, 3894–3899. [Google Scholar] [CrossRef] [PubMed]
  28. Ogunyankin, M.O.; Longo, M.L. Metastability in pixelation patterns of coexisting fluid lipid bilayer phases imposed by e-beam patterned substrates. Soft Matter 2013, 9, 2037–2046. [Google Scholar] [CrossRef] [PubMed]
  29. Ogunyankin, M.O.; Huber, D.L.; Sasaki, D.Y.; Longo, M.L. Nanoscale Patterning of Membrane-Bound Proteins Formed through Curvature-Induced Partitioning of Phase-Specific Receptor Lipids. Langmuir 2013, 29. [Google Scholar] [CrossRef] [PubMed]
  30. Ogunyankin, M.O.; Torres, A.; Yaghmaie, F.; Longo, M.L. Lipid domain pixelation patterns imposed by e-beam fabricated substrates. Langmuir ACS J. Surf. Colloids 2012, 28, 7107–7113. [Google Scholar] [CrossRef] [PubMed]
  31. Black, J.C.; Cheney, P.P.; Campbell, T.; Knowles, M.K. Membrane curvature based lipid sorting using a nanoparticle patterned substrate. Soft Matter 2014, 10, 2016–2023. [Google Scholar] [CrossRef] [PubMed]
  32. Cheney, P.P.; Weisgerber, A.W.; Feuerbach, A.M.; Knowles, M.K. Single Lipid Molecule Dynamics on Supported Lipid Bilayers with Membrane Curvature. Membranes 2017, 7, 15. [Google Scholar] [CrossRef] [PubMed]
  33. Kabbani, A.M.; Kelly, C.V. The Detection of Nanoscale Membrane Bending with Polarized Localization Microscopy. Biophys. J. 2017, 113, 1782–1794. [Google Scholar] [CrossRef]
  34. Kabbani, A.M.; Kelly, C.V. Nanoscale membrane budding induced by CTxB on quasi-one component lipid bilayers detected by polarized localization microscopy. Biophys. J. 2017, 113, 1795–1806. [Google Scholar] [CrossRef]
  35. Chao, M.-H.; Lin, Y.-T.; Dhenadhayalan, N.; Lee, H.-L.; Lee, H.-Y.; Lin, K.-C. 3D Probed Lipid Dynamics in Small Unilamellar Vesicles. Small Weinh. Bergstr. Ger. 2017, 13. [Google Scholar] [CrossRef] [PubMed]
  36. Risselada, H.J.; Marrink, S.J. Curvature effects on lipid packing and dynamics in liposomes revealed by coarse grained molecular dynamics simulations. Phys. Chem. Chem. Phys. PCCP 2009, 11, 2056–2067. [Google Scholar] [CrossRef] [PubMed]
  37. Cooke, I.R.; Deserno, M. Coupling between Lipid Shape and Membrane Curvature. Biophys. J. 2006, 91, 487–495. [Google Scholar] [CrossRef] [PubMed]
  38. Betzig, E.; Patterson, G.H.; Sougrat, R.; Lindwasser, O.W.; Olenych, S.; Bonifacino, J.S.; Davidson, M.W.; Lippincott-Schwartz, J.; Hess, H.F. Imaging Intracellular Fluorescent Proteins at Nanometer Resolution. Science 2006, 313, 1642–1645. [Google Scholar] [CrossRef] [PubMed]
  39. Hess, S.T.; Girirajan, T.P.K.; Mason, M.D. Ultra-High Resolution Imaging by Fluorescence Photoactivation Localization Microscopy. Biophys. J. 2006, 91, 4258–4272. [Google Scholar] [CrossRef] [PubMed]
  40. Rust, M.J.; Bates, M.; Zhuang, X. Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM). Nat. Methods 2006, 3, 793–796. [Google Scholar] [CrossRef] [PubMed]
  41. Willig, K.I.; Kellner, R.R.; Medda, R.; Hein, B.; Jakobs, S.; Hell, S.W. Nanoscale resolution in GFP-based microscopy. Nat. Methods 2006, 3, 721–723. [Google Scholar] [CrossRef] [PubMed]
  42. Van de Linde, S.; Löschberger, A.; Klein, T.; Heidbreder, M.; Wolter, S.; Heilemann, M.; Sauer, M. Direct stochastic optical reconstruction microscopy with standard fluorescent probes. Nat. Protoc. 2011, 6, 991–1009. [Google Scholar] [CrossRef] [PubMed]
  43. Deschout, H.; Zanacchi, F.C.; Mlodzianoski, M.; Diaspro, A.; Bewersdorf, J.; Hess, S.T.; Braeckmans, K. Precisely and accurately localizing single emitters in fluorescence microscopy. Nat. Methods 2014, 11, 253–266. [Google Scholar] [CrossRef] [PubMed]
  44. Huang, B.; Jones, S.A.; Brandenburg, B.; Zhuang, X. Whole-cell 3D STORM reveals interactions between cellular structures with nanometer-scale resolution. Nat. Methods 2008, 5, 1047–1052. [Google Scholar] [CrossRef] [PubMed]
  45. Shtengel, G.; Galbraith, J.A.; Galbraith, C.G.; Lippincott-Schwartz, J.; Gillette, J.M.; Manley, S.; Sougrat, R.; Waterman, C.M.; Kanchanawong, P.; Davidson, M.W.; et al. Interferometric fluorescent super-resolution microscopy resolves 3D cellular ultrastructure. Proc. Natl. Acad. Sci. USA 2009, 106, 3125–3130. [Google Scholar] [CrossRef] [PubMed]
  46. Schnitzbauer, J.; McGorty, R.; Huang, B. 4Pi fluorescence detection and 3D particle localization with a single objective. Opt. Express 2013, 21, 19701–19708. [Google Scholar] [CrossRef] [PubMed]
  47. Pavani, S.R.P.; DeLuca, J.G.; Piestun, R. Polarization sensitive, three-dimensional, single-molecule imaging of cells with a double-helix system. Opt. Express 2009, 17, 19644–19655. [Google Scholar] [CrossRef] [PubMed]
  48. Aguet, F.; Geissbühler, S.; Märki, I.; Lasser, T.; Unser, M. Super-resolution orientation estimation and localization of fluorescent dipoles using 3-D steerable filters. Opt. Express 2009, 17, 6829–6848. [Google Scholar] [CrossRef] [PubMed]
  49. Backer, A.S.; Backlund, M.P.; Lew, M.D.; Moerner, W.E. Single-molecule orientation measurements with a quadrated pupil. Opt. Lett. 2013, 38, 1521–1523. [Google Scholar] [CrossRef] [PubMed]
  50. Backlund, M.P.; Lew, M.D.; Backer, A.S.; Sahl, S.J.; Grover, G.; Agrawal, A.; Piestun, R.; Moerner, W.E. Simultaneous, accurate measurement of the 3D position and orientation of single molecules. Proc. Natl. Acad. Sci. USA 2012, 109, 19087–19092. [Google Scholar] [CrossRef] [PubMed]
  51. Sund, S.E.; Swanson, J.A.; Axelrod, D. Cell Membrane Orientation Visualized by Polarized Total Internal Reflection Fluorescence. Biophys. J. 1999, 77, 2266–2283. [Google Scholar] [CrossRef]
  52. Soumpasis, D.M. Theoretical analysis of fluorescence photobleaching recovery experiments. Biophys. J. 1983, 41, 95–97. [Google Scholar] [CrossRef]
  53. Axelrod, D.; Koppel, D.E.; Schlessinger, J.; Elson, E.; Webb, W.W. Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. Biophys. J. 1976, 16, 1055–1069. [Google Scholar] [CrossRef]
  54. Guo, S.-M.; He, J.; Monnier, N.; Sun, G.; Wohland, T.; Bathe, M. Bayesian approach to the analysis of fluorescence correlation spectroscopy data II: Application to simulated and in vitro data. Anal. Chem. 2012, 84, 3880–3888. [Google Scholar] [CrossRef] [PubMed]
  55. Ries, J.; Schwille, P. Fluorescence correlation spectroscopy. BioEssays News Rev. Mol. Cell. Dev. Biol. 2012, 34, 361–368. [Google Scholar] [CrossRef] [PubMed]
  56. Vestergaard, C.L. Optimizing experimental parameters for tracking of diffusing particles. Phys. Rev. E 2016, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Seu, K.J.; Pandey, A.P.; Haque, F.; Proctor, E.A.; Ribbe, A.E.; Hovis, J.S. Effect of Surface Treatment on Diffusion and Domain Formation in Supported Lipid Bilayers. Biophys. J. 2007, 92, 2445–2450. [Google Scholar] [CrossRef] [PubMed]
  58. Shi, X.; Kohram, M.; Zhuang, X.; Smith, A.W. Interactions and Translational Dynamics of Phosphatidylinositol Bisphosphate (PIP2) Lipids in Asymmetric Lipid Bilayers. Langmuir ACS J. Surf. Colloids 2016, 32, 1732–1741. [Google Scholar] [CrossRef] [PubMed]
  59. Eggeling, C.; Ringemann, C.; Medda, R.; Schwarzmann, G.; Sandhoff, K.; Polyakova, S.; Belov, V.N.; Hein, B.; von Middendorff, C.; Schonle, A.; et al. Direct observation of the nanoscale dynamics of membrane lipids in a living cell. Nature 2009, 457, 1159–1162. [Google Scholar] [CrossRef] [PubMed]
  60. Berglund, A.J. Statistics of camera-based single-particle tracking. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2010, 82. [Google Scholar] [CrossRef] [PubMed]
  61. Thompson, R.E.; Larson, D.R.; Webb, W.W. Precise nanometer localization analysis for individual fluorescent probes. Biophys. J. 2002, 82, 2775–2783. [Google Scholar] [CrossRef]
  62. Fujiwara, T.; Ritchie, K.; Murakoshi, H.; Jacobson, K.; Kusumi, A. Phospholipids undergo hop diffusion in compartmentalized cell membrane. J. Cell Biol. 2002, 157, 1071–1081. [Google Scholar] [CrossRef] [PubMed]
  63. Hsieh, C.-L.; Spindler, S.; Ehrig, J.; Sandoghdar, V. Tracking single particles on supported lipid membranes: Multimobility diffusion and nanoscopic confinement. J. Phys. Chem. B 2014, 118, 1545–1554. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Simulations were performed by mimicking the membrane topography throughout the fission and fusion of a single 50 nm radius vesicle with a planar membrane. (A) By maintaining a 50 nm radius of curvature on top and a smooth, 20-nm radius of curvature connection to the surrounding membrane, the stage of membrane budding was tracked via the maximum bud height above the surrounding membrane (hbud); (B) Bud growth results in an increased area of total membrane within a 250 nm of the bud center. When hbud > 137 nm, the bud undergoes fission, a separate vesicle is formed, which is assumed to diffuse away and not contribute to these simulations; (C) Simulated trajectories are shown in assorted colors over a membrane topography with hbud = 120 nm.
Figure 1. Simulations were performed by mimicking the membrane topography throughout the fission and fusion of a single 50 nm radius vesicle with a planar membrane. (A) By maintaining a 50 nm radius of curvature on top and a smooth, 20-nm radius of curvature connection to the surrounding membrane, the stage of membrane budding was tracked via the maximum bud height above the surrounding membrane (hbud); (B) Bud growth results in an increased area of total membrane within a 250 nm of the bud center. When hbud > 137 nm, the bud undergoes fission, a separate vesicle is formed, which is assumed to diffuse away and not contribute to these simulations; (C) Simulated trajectories are shown in assorted colors over a membrane topography with hbud = 120 nm.
Membranes 07 00060 g001
Figure 2. Increasing the bud height does not result in significant changes to the fluorescence recovery after photobleaching (FRAP) recovery time, including when the lipid diffusion was slowed to DBud = 0.1 or 0.04 μm2/s on the curved membrane for both (A,B) complete bleaching and (C,D) Gaussian Bleaching; (A,C) Simulated I(t) traces while DPlane = 1 μm2/s and DBud = 0.1 μm2/s shows the recovery of I(t) after bleaching. There was no apparent trend in the recovery rate changing with bud height; (B,D) The recovery rate was quantified by fitting Equation (2) to find τFRAP from I(t) of each condition. Error bars represent the standard error of the mean between four repeated simulations.
Figure 2. Increasing the bud height does not result in significant changes to the fluorescence recovery after photobleaching (FRAP) recovery time, including when the lipid diffusion was slowed to DBud = 0.1 or 0.04 μm2/s on the curved membrane for both (A,B) complete bleaching and (C,D) Gaussian Bleaching; (A,C) Simulated I(t) traces while DPlane = 1 μm2/s and DBud = 0.1 μm2/s shows the recovery of I(t) after bleaching. There was no apparent trend in the recovery rate changing with bud height; (B,D) The recovery rate was quantified by fitting Equation (2) to find τFRAP from I(t) of each condition. Error bars represent the standard error of the mean between four repeated simulations.
Membranes 07 00060 g002
Figure 3. A membrane bud may be detected by fluorescence correlation spectroscopy (FCS) if hbud > 80 nm or the membrane curvature slowed the local lipid diffusion. While a lipid on the planar membrane experiences DPlane = 1 μm2/s, varying the diffusion rate of the lipid on the curved membrane and the location of the bud within the excitation spot affects the observed FCS results. (A) A shifting of G(τ) to longer lag times was apparent when DBud = 0.1 μm2/s and the bud was centered on the 250 nm wide Gaussian illumination; (B) If DBud = 1 μm2/s and hbud ≤ 50 nm, a minimal effect of the bud was observed on the measured τFCS; the extra membrane area and varying membrane orientation are not sufficient to affect the FCS results within the uncertainty of finding τFCS. However, if DBud = 0.1 μm2/s, the presence of the bud becomes clear for hbud ≥ 15 nm. The bud can be offset by 200 nm from the center of the illumination and still yield a clear change in the τFCS. If the observation spot was off-centered from the bud by 100 nm or 200 nm, the effects of the bud are decreased by an average of 15% or 66%, respectively.
Figure 3. A membrane bud may be detected by fluorescence correlation spectroscopy (FCS) if hbud > 80 nm or the membrane curvature slowed the local lipid diffusion. While a lipid on the planar membrane experiences DPlane = 1 μm2/s, varying the diffusion rate of the lipid on the curved membrane and the location of the bud within the excitation spot affects the observed FCS results. (A) A shifting of G(τ) to longer lag times was apparent when DBud = 0.1 μm2/s and the bud was centered on the 250 nm wide Gaussian illumination; (B) If DBud = 1 μm2/s and hbud ≤ 50 nm, a minimal effect of the bud was observed on the measured τFCS; the extra membrane area and varying membrane orientation are not sufficient to affect the FCS results within the uncertainty of finding τFCS. However, if DBud = 0.1 μm2/s, the presence of the bud becomes clear for hbud ≥ 15 nm. The bud can be offset by 200 nm from the center of the illumination and still yield a clear change in the τFCS. If the observation spot was off-centered from the bud by 100 nm or 200 nm, the effects of the bud are decreased by an average of 15% or 66%, respectively.
Membranes 07 00060 g003
Figure 4. Dxy mapped over the sample through single particle tracking (SPT). The single steps from all trajectories are binned according to the average x position and the average y position of the two localizations. This two-dimensional (2D) binning of the x and y positions allows for the analysis of all the step lengths in each region on the sample, as represented by a single pixel in these images. Here, hbud = 100 nm, σr = 15 nm, DPlane = 1 μm2/s, and tframe = 2 ms while (A) DBud = 1 μm2/s; (B) DBud = 0.1 μm2/s; and (C) DBud = 0.04 μm2/s. The bud induced slowing in (A) was due to the membrane topography causing the lipid to move slower through the xy-plane with constant local diffusion in the membrane of varying orientation; this effect also contributes to (B,C). The distinction between the bud and the planar membrane is enhanced when there is a greater difference between DPlane and DBud.
Figure 4. Dxy mapped over the sample through single particle tracking (SPT). The single steps from all trajectories are binned according to the average x position and the average y position of the two localizations. This two-dimensional (2D) binning of the x and y positions allows for the analysis of all the step lengths in each region on the sample, as represented by a single pixel in these images. Here, hbud = 100 nm, σr = 15 nm, DPlane = 1 μm2/s, and tframe = 2 ms while (A) DBud = 1 μm2/s; (B) DBud = 0.1 μm2/s; and (C) DBud = 0.04 μm2/s. The bud induced slowing in (A) was due to the membrane topography causing the lipid to move slower through the xy-plane with constant local diffusion in the membrane of varying orientation; this effect also contributes to (B,C). The distinction between the bud and the planar membrane is enhanced when there is a greater difference between DPlane and DBud.
Membranes 07 00060 g004
Figure 5. Azimuthal averaging of the spatial mapping of Dxy around a bud (i.e., Figure 4) improved the data statistics and enabled easier comparison of Dxy versus distance from the bud center (r) and hbud. Simulated SPT results for Dxy are presented with varying σr, tframe, and DBud. (A) When DPlane = DBud = 1 μm2/s and, the single molecules maintain a uniform average local speed through the membrane, and the observed variations in Dxy were due to the tilt of the membrane; (B,C) As DBud decreased to 0.1 or 0.04 μm2/s, the bud became increasingly apparent at all hbud; The (DF) significant improvements in the resolution of Dxy were obtained by decreasing tframe, even at the expense of increasing σr; (F) Differences in Dxy may be observed across the bud, including slowing at the vertical edge of the bud.
Figure 5. Azimuthal averaging of the spatial mapping of Dxy around a bud (i.e., Figure 4) improved the data statistics and enabled easier comparison of Dxy versus distance from the bud center (r) and hbud. Simulated SPT results for Dxy are presented with varying σr, tframe, and DBud. (A) When DPlane = DBud = 1 μm2/s and, the single molecules maintain a uniform average local speed through the membrane, and the observed variations in Dxy were due to the tilt of the membrane; (B,C) As DBud decreased to 0.1 or 0.04 μm2/s, the bud became increasingly apparent at all hbud; The (DF) significant improvements in the resolution of Dxy were obtained by decreasing tframe, even at the expense of increasing σr; (F) Differences in Dxy may be observed across the bud, including slowing at the vertical edge of the bud.
Membranes 07 00060 g005

Share and Cite

MDPI and ACS Style

Kabbani, A.M.; Woodward, X.; Kelly, C.V. Revealing the Effects of Nanoscale Membrane Curvature on Lipid Mobility. Membranes 2017, 7, 60. https://doi.org/10.3390/membranes7040060

AMA Style

Kabbani AM, Woodward X, Kelly CV. Revealing the Effects of Nanoscale Membrane Curvature on Lipid Mobility. Membranes. 2017; 7(4):60. https://doi.org/10.3390/membranes7040060

Chicago/Turabian Style

Kabbani, Abir Maarouf, Xinxin Woodward, and Christopher V. Kelly. 2017. "Revealing the Effects of Nanoscale Membrane Curvature on Lipid Mobility" Membranes 7, no. 4: 60. https://doi.org/10.3390/membranes7040060

APA Style

Kabbani, A. M., Woodward, X., & Kelly, C. V. (2017). Revealing the Effects of Nanoscale Membrane Curvature on Lipid Mobility. Membranes, 7(4), 60. https://doi.org/10.3390/membranes7040060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop