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/nm
2 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 μm
2/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 10
5 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
with a
w = 250 nm and r
i 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
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/nm
2 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
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
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
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
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 (
) 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
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.
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
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/(6
tframe) ≤ 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,
is proportional to
. Assuming minimal image readout time, minimal background noise, and
texp ≈
tframe,
and
σr are approximately inversely related to each other, and an imaging frame rate could be found that results in
≈
σ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
being decreased by 90% (from 20 ms to 2 ms),
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
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
and there are no longer benefits of decreasing
tframe.
Increasing the photon flux can improve both
σr and
. 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.