Yuyang Wang1, Peter Zijlstra1. 1. Molecular Biosensing for Medical Diagnostics, Faculty of Applied Physics, and Institute for Complex Molecular Systems, Eindhoven University of Technology, PO Box 513, 5600 MB, Eindhoven, The Netherlands.
Abstract
We present a numerical study on plasmon-enhanced single-molecule enzymology. We combine Brownian dynamics and electromagnetic simulations to calculate the enhancement of fluorescence signals of fluorogenic substrate converted by an enzyme conjugated to a plasmonic particle. We simulate the Brownian motion of a fluorescent product away from the active site of the enzyme, and calculate the photon detection rate taking into account modifications of the excitation and emission processes by coupling to the plasmon. We show that plasmon enhancement can boost the signal-to-noise ratio (SNR) of single turnovers by up to 100 fold compared to confocal microscopy. This enhancement factor is a trade-off between the reduced residence time in the near-field of the particle, and the enhanced emission intensity due to coupling to the plasmon. The enhancement depends on the size, shape and material of the particle and the photophysical properties of the fluorescent product. Our study provides guidelines on how to enhance the SNR of single-molecule enzyme studies and may aid in further understanding and quantifying static and dynamic heterogeneity.
We present a numerical study on plasmon-enhanced single-molecule enzymology. We combine Brownian dynamics and electromagnetic simulations to calculate the enhancement of fluorescence signals of fluorogenic substrate converted by an enzyme conjugated to a plasmonic particle. We simulate the Brownian motion of a fluorescent product away from the active site of the enzyme, and calculate the photon detection rate taking into account modifications of the excitation and emission processes by coupling to the plasmon. We show that plasmon enhancement can boost the signal-to-noise ratio (SNR) of single turnovers by up to 100 fold compared to confocal microscopy. This enhancement factor is a trade-off between the reduced residence time in the near-field of the particle, and the enhanced emission intensity due to coupling to the plasmon. The enhancement depends on the size, shape and material of the particle and the photophysical properties of the fluorescent product. Our study provides guidelines on how to enhance the SNR of single-molecule enzyme studies and may aid in further understanding and quantifying static and dynamic heterogeneity.
Enzymes are
vital biomolecules
that catalyze thousands of biochemical reactions.[1] The study of enzyme kinetics started with the work of Michaelis
and Menten,[2] which plays an essential role
in analyzing the kinetics of an ensemble of enzymes. Such ensemble
studies assume that all enzymes in solution are identical; however,
recent single-molecule studies revealed that the catalytic rates in
a population of enzymes are heterogeneous (static heterogeneity).[3,4] Furthermore, it was observed that the turnover rate of individual
enzymes exhibits temporal fluctuations (dynamic heterogeneity).[5−8] It may well be that such heterogeneity has far-reaching biological
consequences, slowly fluctuating enzymes are, for example, implicated
in variations in antibiotic resistance of genetically identical bacteria.[9]Single-molecule confocal microscopy has
been used to record the
kinetic behavior of single immobilized enzymes,[5−8,10−14] and often fluorogenic substrates are used that are converted by
the enzyme from a dark to an emitting state.[15,16] The activity of the enzyme is then visualized as a series of fluorescence
bursts caused by the diffusion of the fluorescent product out of the
confocal volume.[16] This approach has been
used to study the activity of enzymes such as cholesterol oxidase,[5] β-galactosidase,[7] chymotrypsin,[14] lipase B,[6,11] and λ-exonuclease.[12] These pioneering
studies reported temporal fluctuations in activity on time scales
of milliseconds to seconds, which were attributed to conformational
fluctuations.[5,7,11]However, the average signal-to-noise ratio (SNR) of these bursts
ranges from 1 to 4 in practice.[13,14] The SNR is limited
by the fluorescent background caused by, for example, the presence
of a small fraction of autohydrolized product in the large confocal
volume. In addition, the SNR of individual bursts is broadly distributed
around its mean since each fluorescent product diffuses out of the
confocal volume through a random path determined by Brownian motion.
This complicates the interpretation of single-enzyme data and the
quantification of the degree of heterogeneity because a significant
fraction of the bursts exhibits a SNR < 1 and is missed in the
analysis.[14,15,17−19]Recently, zero-mode waveguides (ZMW) have been used to increase
the SNR in the detection of single-molecule fluorescence.[20,21] ZMWs efficiently reduce the fluorescent background from solution-phase
product by confining light to zeptoliter volumes in small apertures
that are milled in a metal film. Using ZMWs, single turnovers of DNA
polymerase were recorded with increased SNR and higher background
fluorophore concentrations.[20] However,
signal enhancements are modest due to quenching by the large volume
of nearby metal.[22,23] The plasmon resonance of single
metal nanoparticles on the other hand provides both confinement of
the excitation to zeptoliter volumes, and can additionally strongly
enhance the emitted intensity because the particle acts as a nanoantenna.[24,25]The origin of the plasmon enhancement of single-molecule fluorescence
is the modification of the excitation and decay rates of a molecular
dipole close to a particle. The locally enhanced field leads to a
dramatic increase of the excitation rate, whereas the nanoantenna
modifies the emission of the fluorophore due to modulations of the
radiative and nonradiative decay rate.[26−29] Fluorescence intensity enhancements
in excess of 103 to 105 have been reported for
single gold nanorods,[30,31] metallic particles on a film,[32] lithographically fabricated nanogaps,[33] and self-assembled gold nanoparticle dimers.[34−36] Plasmonic nanoantennae have already been used to extend fluorescence
correlation spectroscopy to micromolar concentations,[21,37−42] to study single light-harvesting complexes,[43] and to study diffusion of membrane-proteins.[44]Here we propose the use of single plasmonic nanoparticles
to enhance
the SNR of single-molecule enzymology. We combine Brownian dynamics
and electromagnetic simulations to model the number of photons generated
by a single fluorescent product that diffuses away from the active
site of the enzyme. The simulations account for modifications of the
excitation, radiative, and nonradiative rates of the fluorophore by
the presence of the nanoantenna. Our approach differs from previous
ensemble-averaged approaches of plasmon-enhanced fluorescence where
molecular excitation and detection functions were used to simulate
fluorescence correlation spectroscopy (FCS) curves.[37,38] We explicitly simulate the trajectories of each single molecule
and reconstruct the time-dependent intensity and SNR on a molecule-by-molecule
basis. In contrast to fluorescence-correlation spectroscopy models,
this yields information on molecule-to-molecule variations in the
detected intensity and SNR. Establishing these variations is especially
critical for single-molecule enzymology where no event should be missed.
The Brownian dynamics simulations reveal that a fluorescent product
resides in the particle’s near-field 100–1000×
shorter compared to a confocal excitation volume. Despite this shorter
residence time the resulting SNR is increased up to 100-fold due to
enhancements by the nanoantenna. We show that the SNR enhancement
depends on the size, shape and material of the plasmonic particle
and on the photophysical properties of the fluorescent product.
Computational
Approach
We consider a particle–enzyme complex in
the presence of
a solution of fluorogenic substrate molecules, placed in a diffraction-limited
confocal volume of a laser beam (Figure ). In a typical confocal fluorescence measurement,
a focused laser beam illuminates a confined volume of fluorophores
that undergo Brownian motion. The enzyme converts a nonfluorescent
substrate to a fluorescent product, and this generated fluorophore
will follow a Brownian motion trajectory (a series of 3D coordinates
as a function of time) to eventually diffuse out of the focus.
Figure 1
Confocal excitation
of an enzyme-gold nanoparticle complex. An
enzyme is coupled to a gold nanoparticle that is immobilized on a
glass surface (blue) and located in the center of a confocal laser
beam (green). The active site of the enzyme converts nonfluorescent
substrate to an emitting product (red) that follows a Brownian motion
trajectory once released from the active site.
Confocal excitation
of an enzyme-gold nanoparticle complex. An
enzyme is coupled to a gold nanoparticle that is immobilized on a
glass surface (blue) and located in the center of a confocal laser
beam (green). The active site of the enzyme converts nonfluorescent
substrate to an emitting product (red) that follows a Brownian motion
trajectory once released from the active site.To assess the effect of the nanoantenna on the SNR we compare
the
SNR in the absence and the presence of the particle. We focus our
analysis on spherical and rod-shaped gold particles due to their widespread
use and ease of synthesis.[45] The simulation
of the signal and background signals consists of the following steps:Generation
of Brownian motion trajectories
of single product molecules that diffuse away from the active site
of the enzyme. The trajectories consist of a 3D random walk with a
total length of 1 ms and a given time step of 10 ns. The mean squared
displacement (MSD) of the molecule derives from the Stokes–Einstein
equation: with D being the diffusion
coefficient (m–2 s–1), kB is the Boltzmann’s constant (J K–1), T is the absolute temperature
(K), η is the dynamic viscosity of water (Pa s), and r is the hydrodynamic radius of the molecule (m). The MSD
for 3D free diffusion is then given by (see Supporting Information, Figure S4):Calculation of the position-dependent
emission rate, taking into account the plasmon-induced modifications
of the excitation rate, radiative, and nonradiative decay rate of
the fluorophore. The photon count rate PCR of a fluorophore can be
expressed as where Iexc is
the position-dependent excitation intensity (W m–2), σabs is the absorption cross section of the fluorophore
(m2), ϕ is the fluorophore’s quantum yield
defined as , with γr and γnr being the radiative
and nonradiative decay rates (s–1) and the
photon energy (J).Equation is valid in the weak excitation regime, in which the
photon count rate of a fluorophore is linear in incident intensity.
To take into account saturation effects we introduce the saturation
intensity with γtot = γr + γnr, where we ignored triplet state dynamics.[43] This yields the following expression for the
PCR of a single product (see Supporting Information):In the absence
of a nanoparticle, the only
position-dependent term in the above equations is Iexc, whose position dependence is given by the 3D Gaussian
approximation of the confocal volume. The other terms are assumed
position-independent. In the presence of a plasmonic particle, however,
the decay rates γr and γnr will
also be modified depending on the location of the fluorophore with
respect to the particle. These rate-modifications are calculated numerically
using the Boundary Element Method. In all simulations we used the
typical photophysical parameters of a red fluorescence dye in free
space: the absorption cross section σabs of 10–20 m2, the fluorescence lifetime τlifetime of 10 ns, and corresponding total decay rate γtotof 108 s–1.Calculation of the background signal
PCRbg for a certain fraction of autohydrolized substrate
using eq , but without
plasmon-enhancement. Herein we ignore the possible background and
noise from autohydrolyzed product diffusing through the plasmon-enhanced
near field because the near field is >104× smaller
than the confocal volume. We then find PCRbg bywhere here the PCR for a fluorophore in the
center of a 3D Gaussian beam is used, nA is Avogadro’s constant, csubs is the substrate concentration (M), Vconf is the confocal volume (L), and fh is
the fraction of autohydrolyzed substrate. In order to acquire a realistic
signal and noise level as reported in existing experimental work for
the confocal-only case,[7] we set fh = 0.002‰ and the collection efficiency
of the setup ηcol = 6%.We then calculate the
PCR of single diffusing products for each
position of the fluorophore in the Brownian trajectory, yielding the
time-dependent PCR(τ). Time traces of multiple subsequent turnovers
are generated by randomly assigning a time stamp to every turnover
and stitching the corresponding PCR(τ). The average frequency
of events was set to 400 Hz to match the typical turnover rate of
an enzyme with reasonable activity.[1,46] We then generate
the expected experimental timetrace by taking into account the collection
efficiency of the setup ηcol and the integration
time tint. This yields the number of detected
photons per integration time, Ndet(t):Based on the above calculations, we
determine the signal-to-noise ratio (SNR) of each event by dividing
the peak photon counts from a single turnover by the shotnoise originating
from both the signal and background counts.
Results and Discussion
The electromagnetic field around
the particle is enhanced compared
to the incident illumination, resulting in an enhanced excitation
rate of fluorophores located in this near field. The excitation enhancement
can therefore be defined as , where I0 is
the maximum intensity at the center of an incident 3D Gaussian excitation
beam. In Figure a,
we show the numerically calculated near field intensity around a nanorod
of 82 × 30 nm2, evaluated on resonance with the longitudinal
plasmon at 705 nm. We observe an enhanced intensity around the two
tips of the gold nanorod of nearly , in good agreement with previous studies.[31] This enhancement decays rapidly away from the
particle surface on length-scales of ∼5 nm. The excitation
is thus strongly confined to the particle surface with a volume nearly
104× smaller than the diffraction limit.
Figure 2
Numerical calculations
of the local intensity enhancement and rate-enhancements
around an 82 × 30 nm2 gold nanorod. (a) Excitation
enhancement shown by normalized Iexc containing
both distributions from the 3D Gaussian point-spread-function of an
excitation beam with a maximum intensity of I0, and the gold nanorod near field. The near field enhancement
around the particle was calculated for excitation with a plane wave
at 705 nm (resonant with the longitudinal plasmon). (b) Zoom of the
gold nanorod in (a). (c, d) Quantum yield enhancement for a dipolar
emitter with an intrinsic quantum yield ϕ0 of 1%
and 25%, respectively (orientation averaged). The emitter is modeled
by a narrow absorption line at 705 nm, with a narrow emission line
with a Stokes shift of 25 nm. (e, f) Total fluorescence enhancement
for excitation at 705 nm and emission at 730 nm. Note the difference
in the color scale.
Numerical calculations
of the local intensity enhancement and rate-enhancements
around an 82 × 30 nm2 gold nanorod. (a) Excitation
enhancement shown by normalized Iexc containing
both distributions from the 3D Gaussian point-spread-function of an
excitation beam with a maximum intensity of I0, and the gold nanorod near field. The near field enhancement
around the particle was calculated for excitation with a plane wave
at 705 nm (resonant with the longitudinal plasmon). (b) Zoom of the
gold nanorod in (a). (c, d) Quantum yield enhancement for a dipolar
emitter with an intrinsic quantum yield ϕ0 of 1%
and 25%, respectively (orientation averaged). The emitter is modeled
by a narrow absorption line at 705 nm, with a narrow emission line
with a Stokes shift of 25 nm. (e, f) Total fluorescence enhancement
for excitation at 705 nm and emission at 730 nm. Note the difference
in the color scale.In addition to the enhanced
excitation rate, also the decay rates
of the fluorophore are modified around the particle, resulting in
a modified quantum yield. We quantify the emission enhancement as , which was calculated for a single-wavelength
dipole emitting at 730 nm (Stokes shifted by 25 nm from the excitation)
and was orientation averaged (see Supporting Information). Note here that we assume the dipole has rotational correlation
times shorter than the fluorescence lifetime and that the dipole orientation
is fast tumbling within a fluorescent decay (see the influence of
a slowly tumbling molecule in Figure S3). In Figure c we
show for ϕ0 = 1%. The emission
enhancement depends strongly on the molecule–surface separation
and exhibits a maximum value of 28 at a distance of 4 nm along the
center axis of the nanorod. The significant increase in quantum yield
especially for positions around the tips of the nanorod is due to
the enhancement of the radiative decay rate due to coupling of the
dipole emission to the dipolar plasmon of the nanorod. The drastic
decrease of the quantum yield at a distance less than 4 nm is attributed
to enhanced nonradiative decay through higher order plasmon modes,[47] which accounts for quenching by the nearby particle.
Usually single-molecule enzyme studies are performed with higher quantum
yield dyes, for which we show the enhancement factors for ϕ0 = 25% in Figure d. We find a maximum of
only 2 at a distance of 15 nm from the
surface because of the already high ϕ0, here the
excited-state decay is limited by the radiation efficiency of plasmons.[47] Note here that the distance of maximum quantum
yield enhancement is different for dyes with different intrinsic quantum
yield. This is due to the fact that the enhancement of the nonradiative
rate decays faster with separation than the radiative rate enhancement,
resulting in a maximal quantum yield enhancement at a certain distance
from the surface that depends on the intrinsic quantum yield. The
total fluorescence enhancement is
then a combination of excitation enhancement
and emission enhancement, given byThis total fluorescence enhancement is shown in Figure e,f, again for two
values of the intrinsic
quantum yield. For ϕ0 = 1%, the maximum is as high as 20000 at
a distance of 2.5
nm from the tip, whereas for ϕ0 = 25%, decreases to a maximum
of 1000 at a distance
of 2.8 nm. These results imply that the majority of the enhancement
is due to excitation enhancement, whereas the emission enhancement
is significant only for fluorophores with a low ϕ0. The positions with maximum fluorescence enhancement will be used
as the exact location of the enzyme active site, which is the starting
point of the Brownian motion trajectories.In addition to enhancement
of the excitation rate and quantum yield,
also the saturation intensity Isat is
modified in the presence of a plasmonic particle.[43] Since higher Isat directly
leads to higher photon emission rates, an increased Isat is beneficial to the detection of a fluorescent product.
Experimentally 50-fold higher photon count rate has been reported,
combined with a 10-fold higher total photon count.[43] In our calculation, where , the increase of saturation level is determined
by the modification of distance-dependent total decay rate γtot and is found to be 3 orders of magnitude enhanced on the
particle surface, and decreases rapidly within 10 nm.To quantify
the effect of plasmon-enhancement, we consider the
enzyme reaction in the presence and absence of the plasmonic particle.
In Figure a we show
a typical trajectory of a single fluorophore diffusing away from the
center of a Gaussian focus, and Figure b shows the time-dependent PCR(τ)of a single
turnover event according to eq . The residence time in the confocal volume varies from 100
μs to 1 ms for individual fluorophores due to the random nature
of the trajectories. In Figure c,d we show the trajectory and the plasmon-enhanced photon
counts for a fluorophore (ϕ0 = 25%) starting at 2.8
nm from the tip of a 82 × 30 nm2 nanorod in a 3D Gaussian
excitation beam. The plasmon-enhanced time trace shows a 1000-fold
enhancement in maximum photon detection rate compared to the confocal-only
case, in agreement with Figure . Although the intensity is enhanced 1000-fold, the fluorophore
spends on average 100× shorter in the near-field compared to
a confocal volume.
Figure 3
Simulated Brownian motion trajectories of a fluorophore
out of
a 3D Gaussian focus (a, c) without and (b, d) with a 82 × 30
nm2 AuNR. The blue half-ellipses illustrate the approximate
1/e2 beam size of the 3D Gaussian focus.
(b, d) Corresponding PCR(τ) as a function of time for the trajectories
shown in (a) and (c). The inset in (d) shows the first few points
in the trajectory projected on the x–z plane. The dotted line in the inset illustrates the planar
infinite substrate, acting as a hard boundary for the diffusion. The
simulation is performed for ϕ0 = 25% and I0 = 108 W m–2 (corresponding
to a 705 nm laser with a power of 20 μW focused with an objective
lens with NA = 1.2) and a time step of 10 ns in the Brownian motion
calculations.
Simulated Brownian motion trajectories of a fluorophore
out of
a 3D Gaussian focus (a, c) without and (b, d) with a 82 × 30
nm2 AuNR. The blue half-ellipses illustrate the approximate
1/e2 beam size of the 3D Gaussian focus.
(b, d) Corresponding PCR(τ) as a function of time for the trajectories
shown in (a) and (c). The inset in (d) shows the first few points
in the trajectory projected on the x–z plane. The dotted line in the inset illustrates the planar
infinite substrate, acting as a hard boundary for the diffusion. The
simulation is performed for ϕ0 = 25% and I0 = 108 W m–2 (corresponding
to a 705 nm laser with a power of 20 μW focused with an objective
lens with NA = 1.2) and a time step of 10 ns in the Brownian motion
calculations.In Figure a,b we
show the comparison of the simulated time trace of single enzyme turnovers
without and with plasmonic particle respectively for tint = 1 ms. Despite this much shorter residence time,
the plasmon enhancement wins and we find a significant enhancement
inNdet(t) for experimentally
relevant binning times of tint = 1 ms.
For the confocal excitation, we used the SNR reported in previous
experimental results as a benchmark, and by varying the autohydrolyzed
fraction fh in eq and collection efficiency ηcol in eq , we find an
average SNR = 3 for I0 = 2 × 109 W m–2 that closely resembles literature
reports.[7,14] This indicates that our simulation generates
the correct signal and background intensities. As can be seen, the
plasmon-enhanced time trace (b) exhibits a significantly higher SNR
where all events are clearly visible above the background.
Figure 4
Simulated time
traces of single enzyme turnovers. The number of
detected photons Ndet(t) with shot-noise superposed, without (a) and with (b) a AuNR, and
(c) the histogram showing the SNR distribution evaluated from the
time traces shown in (a) and (b). Simulation is based on an 82 ×
30 nm2 AuNR, which is under 705 nm laser excitation, with I0 = 2 × 109 W m–2. For both simulations we used tint =
1 ms and ϕ0 = 25%.
Simulated time
traces of single enzyme turnovers. The number of
detected photons Ndet(t) with shot-noise superposed, without (a) and with (b) a AuNR, and
(c) the histogram showing the SNR distribution evaluated from the
time traces shown in (a) and (b). Simulation is based on an 82 ×
30 nm2 AuNR, which is under 705 nm laser excitation, with I0 = 2 × 109 W m–2. For both simulations we used tint =
1 ms and ϕ0 = 25%.In Figure c we
show the histogram of the SNR of each event, which shows a 7-fold
increase in average SNR due to plasmon enhancement. More importantly,
in the plasmon-enhanced time trace, all simulated events exhibit SNR
> 1, whereas ∼30% of events exhibit SNR < 1 for confocal-only
excitation. This is crucial for single enzyme measurements because
a high SNR will significantly increase the accuracy of waiting time
determination and will alleviate the uncertainties now experienced
in various threshold-finding algorithms.[14]We show in Figure the SNR with (red curve) and without (blue curve) plasmon
enhancement,
as a function of incident laser intensityI0. We evaluate the SNR under 705 nm laser intensity I0 from 5 × 106 to 2 × 1010 W m–2 corresponding to incident laser powers of
1 μW to 2 mW focused by an objective lens with NA = 1.2. The
errors bars indicate the full width of the distribution of SNR per
event, which spans 1–2 decades due to the broadly distributed
residence time of diffusing fluorophores caused by Brownian motion.
Figure 5
SNR in
single enzyme turnover detection as a function of laser
intensity.The simulation is based on an 82 × 30 nm2 AuNR, with 705 nm laser excitation, tint = 1 ms, and ϕ0 = 25%.
SNR in
single enzyme turnover detection as a function of laser
intensity.The simulation is based on an 82 × 30 nm2 AuNR, with 705 nm laser excitation, tint = 1 ms, and ϕ0 = 25%.For the confocal-only simulations, we find an average SNR
> 1 for I0 > 2 × 108 W m–2, and the SNR does not increase anymore beyond I0 = 5 × 109 W m–2 due
to fluorescence saturation (see eq ). This results in a maximum average SNR in confocal-only
simulations of ∼4. With plasmon enhancement we show that the
SNR can be enhanced by more than an order of magnitude under these
conditions. The increase in SNR under plasmon enhancement is attributed
to the strong fluorescence enhancement around the gold nanorod, implying
that although the residence time is shorter by ∼100×,
a significant increase in PCR(τ) is still expected due to fluorescence
enhancement. However, although the saturation level of fluorophores
close to the particle is increased due to strong increase of γtot, the SNR still saturates at a value of ∼25 due to
the increased (but still limited) saturation intensity.The
structure and activity of proteins can be impaired when they
are heated at a certain temperature for extended periods of time.
From heat conduction calculations we estimate the surface temperature
rise and find that for I0 = 1 × 108 W m–2, the average temperature rise is
19 K (see Supporting Information, Figure S1), well below the denaturing point of many enzymes.[48] In the remainder of our calculations we therefore use I0 = 1 × 108 W m–2.The overall enhancement is a trade-off between fluorescence
enhancement
and the residence time of the fluorophore in the near-field. The latter
depends on the size of the particle, which is illustrated in Figure a. There we show
averaged and normalized time traces of single fluorophores generated
by the enzyme in the absence and presence of a particle, averaged
over ∼400 events. For the confocal-only case we observe an
average residence time of ∼10 μs, whereas the residence
time in the near-field of a nanorod is typically >100× shorter.
As expected, we find that the residence time increases with nanoparticle
width, a direct consequence of the greater extent of the near-field
for larger particles.
Figure 6
(a) Normalized PCR time traces of single fluorophores
diffusing
away from the active site of an enzyme coupled to a AuNR with different
widths. Time traces are averaged over ∼400 trajectories. The
particles have a fixed plasmon wavelength λsp = 705
nm. (b) Average SNR as a function of binning time tint for a nanorod with a width of 30 nm and the confocal-only
case. The excitation wavelength is 705 nm with I0 = 1 × 108 W m–2 and ϕ0 = 25% for all calculations to facilitate comparison. In order
to keep λsp = 705 nm for all particles the aspect
ratios vary slightly.
(a) Normalized PCR time traces of single fluorophores
diffusing
away from the active site of an enzyme coupled to a AuNR with different
widths. Time traces are averaged over ∼400 trajectories. The
particles have a fixed plasmon wavelength λsp = 705
nm. (b) Average SNR as a function of binning time tint for a nanorod with a width of 30 nm and the confocal-only
case. The excitation wavelength is 705 nm with I0 = 1 × 108 W m–2 and ϕ0 = 25% for all calculations to facilitate comparison. In order
to keep λsp = 705 nm for all particles the aspect
ratios vary slightly.The short residence times we find in Figure a imply that the SNR will depend strongly
on the binning time used in the experiments. This effect is shown
in Figure b, where
we compare the confocal-only case with a 30 nm wide nanorod. For all
evaluated binning times tint the average
SNR is significantly increased (by a factor of 10 × 100) due
to plasmon enhancement. For the confocal-only case, we find an optimum
binning time of about 10 μs, which matches with the average
residence time of the product in the laser focus. For shorter integration
time, not all signal is collected within one integration time, resulting
in a lower SNR of the bursts. In the presence of a 30 nm wide nanorod
we find that the SNR increases monotonously with decreasing binning
time, with 100-fold enhancement compared to the confocal SNR for tint = 1 μs. Again we expect an optimum
in the SNR when tint matches the residence
time in the near-field, but this now occurs for tint < 1 μs, which is difficult to access experimentally.
For experimentally feasible binning times of 100 μs we find
that the SNR in the presence of a 30 nm wide particle is approximately
enhanced 5× compared to the confocal-only case.In Figure a, we
show the SNR as a function of particle width. To compare between different
plasmonic materials, we also calculate the SNR enhancement for silver
nanorods (AgNRs). As expected we find that the SNR enhancement with
AgNRs is higher than AuNRs by a factor up to 10, caused by a reduced
absorption by the bulk metal which results in a higher local field
enhancement and dipolar radiation efficiency (see Figure S5 in Supporting Information). For widths from 10 to
20 nm the SNR increases for both AuNRs and AgNRs, which is the result
of an increased residence time in the near-field combined with a higher
emission enhancement for larger particles. For increasing widths radiation
damping causes broadening of the plasmon resonance[49] and concurrent reduction in fluorescence enhancement.[50] For AgNRs, significant radiation damping starts
for ∼20 nm width, while for AuNRs this onset is shifted to
widths of 30–40 nm (see Figure S6 in Supporting Information) This implies that for AuNRs, even though emission
enhancement is stronger for larger particles, this effect is similar
in magnitude as the reduced excitation enhancement due to radiation
damping. The result is a SNR that does not strongly depend on particle
size for particles with a width larger than 30 nm. For AgNRs, the
stronger radiation damping for larger particles results in a decreasing
SNR for particles with a width larger than 20 nm. This trade-off between
emission enhancement and radiation damping results in an optimum particle
size that is different for gold versus silver due to a different onset
of radiation damping.
Figure 7
(a) SNR as a function of the width of the particle for
both gold
and silver nanorods with λSP = 705 nm, the dotted
blue line indicates the confocal-only SNR. Simulation parameters: I0 = 2 × 109 W m–2 for confocal-only and I0 = 1 ×
108 W m–2 for the plasmon-enhanced case,
ϕ0 = 25%, and tint =
100 μs. Note that the calculations for confocal-only were performed
with a higher power density than the enzyme-nanoparticle conjugates,
see the text for further details. (b) Dependence of the average SNR
on ϕ0 for particles with λSP = 705
nm. The AuNR has size 82 × 30 nm2, the AgNR is 72
× 20 nm2.
(a) SNR as a function of the width of the particle for
both gold
and silver nanorods with λSP = 705 nm, the dotted
blue line indicates the confocal-only SNR. Simulation parameters: I0 = 2 × 109 W m–2 for confocal-only and I0 = 1 ×
108 W m–2 for the plasmon-enhanced case,
ϕ0 = 25%, and tint =
100 μs. Note that the calculations for confocal-only were performed
with a higher power density than the enzyme-nanoparticle conjugates,
see the text for further details. (b) Dependence of the average SNR
on ϕ0 for particles with λSP = 705
nm. The AuNR has size 82 × 30 nm2, the AgNR is 72
× 20 nm2.It is known that fluorophores with a low intrinsic quantum
yield
are enhanced more by coupling to the plasmon.[47] We therefore also investigate the influence of ϕ0on the average SNR, see Figure b. For the confocal-only case, we observe an increasing
SNR with ϕ0simply due to the fact that more signal
is generated for a higher ϕ0. The SNR scales approximately
with the square-root of ϕ0, as expected for shot-noise
limited detection with low (but here nonzero) background. However,
due to the strong decrease of fluorescence enhancement as a function
of ϕ0, the plasmon-enhanced SNR decreases with ϕ0 (see Supporting Information, Figure S7). In addition, lowering ϕ0 reduces the background
due to autohydrolyzed substrate, resulting in optimum SNR for low
ϕ0. This implies that simply by lowering ϕ0 (by, e.g., using quenchers such as methyl viologen[51]) one can improve the SNR in single-molecule
enzymology further by a factor of 2–3.Our results indicate
that the maximum experimentally feasible enhancement
for AuNRs is ∼100× when a fluorogenic substrate with low
ϕ0 is used in combination with 30 nm wide particles.
For silver particles, the enhancement is increased to 300× for
20 nm particles. Note that we assumed here that the fluorescent product
immediately dissociates from the enzyme. In reality a finite dissociation
constant will be associated with the enzyme–product interaction,[14] which results in the product remaining bound
for a short time. This will further increase the SNR due to the longer
residence time in the near-field of the particle. The enhanced SNR
provided by coupling to a plasmonic particle will aid in the quantification
of the kinetic heterogeneity of enzymes by allowing for more reliable
thresholding and burst-detection.[14] Moreover,
it may allow for the first direct measurement of the enzyme–product
dissociation process by relating the measured enhancement factor to
the residence time of the product in the near-field of the particle.
In addition, it has been reported that the enzymatic activity and
long-term stability are modified for nanoparticle-bound enzymes.[52] Our method may therefore allow for the quantification
of the interfacial activation of enzymes and the elucidation of the
underlying mechanisms at the single-molecule level.
Conclusions
In conclusion, we studied plasmon-enhanced single-molecule enzymology
by combining Brownian dynamics and electromagnetic simulations. We
find dramatically enhanced average SNR caused by the enhancement of
the fluorescence intensity of a product near a plasmonic nanorod.
We show that, although the average residence time is greatly reduced
compared to a confocal-only excitation, the SNR can still be largely
improved, and a 100–300-fold increase in average SNR has been
found for a low quantum yield fluorophore with a binning time of 100 μs.
These enhancements are achieved for enzymes conjugated to the tip
of a particle, which can be experimentally achieved using recently
reported tip-specific functionalization protocols.[53,54] Our results provide a practical guideline for choosing the optimum
particle size and experimental parameters such as intrinsic quantum
yield, excitation power, and binning time. The increase in SNR has
important implications in the study of single-molecule enzymology,
since we show that in the presence of a plasmonic nanoparticle all
bursts exhibit SNR > 1, which may tremendously help the data analysis
in single enzyme kinetics measurement. This may aid in further quantification
and understanding of dynamic and static heterogeneity.
Authors: Valentin Flauraud; Raju Regmi; Pamina M Winkler; Duncan T L Alexander; Hervé Rigneault; Niek F van Hulst; María F García-Parajo; Jérôme Wenger; Jürgen Brugger Journal: Nano Lett Date: 2017-02-14 Impact factor: 11.189