Photon statistics in single molecule
orientational imaging

doi:10.1364/OE.15.013597

« Show journal navigation

Photon statistics in single molecule orientational imaging

Matthew R. Foreman, Sherif S. Sherif, and Peter Török »View Author Affiliations

Optics Express, Vol. 15, Issue 21, pp. 13597-13606 (2007)
http://dx.doi.org/10.1364/OE.15.013597

View Full Text Article

Acrobat PDF (242 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Signal-to-noise considerations
3. Probability density function of the number of detected photons
4. Probability density function of time averaged intensity
5. Discussion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (22)
Cited By
Figures (3)
Metrics

Abstract

Optical techniques in single molecule imaging rely heavily on photon counting for data acquisition. Extraction of information from the recorded readings is often done by means of statistical signal processing, however this requires a full knowledge of the photoelectron statistics. In addition to counting statistics we include a specific form of random signal variations namely reorientational dynamics, or wobble to derive the general probability density function of the number of detected photons. The relative importance of the two factors is dependent upon the total number of photons in the system and results are given in all regimes.

1. Introduction

Single molecule detection (SMD) has become an important technique in recent years for studying dynamic processes such as chemical reactions and molecular motions at a fundamental level [1

K. D. Weston and L. S. Goldner “Orientation imaging and reorientation dynamics of single dye molecules,” J. Phys. Chem. B 105 3453-3462 (2001) [CrossRef]

]. Historically these processes are usually studied using methods based on ensemble averaging of a sample of molecules, however these are based on the fundamental assumption that the sample is homogeneous. This ergodic hypothesis is invalid for many biological and chemical systems. Studies on single molecules are thus advantageous as information, such as statistical distributions of particular quantities, is not lost by averaging.

Single molecule imaging techniques, such as fluorescence microscopy, can also be used to track bio-molecular motions. This has applications in the pharmaceutical industry where a good understanding of processes such as protein folding [2

D.M. Warshaw, E. Hayes, D. Gaffney, A.M. Lauzon, J.R. Wu, G. Kennedy, K. Trybus, S. Lowey, and C. Berger, “Myosin conformational states determined by single fluorophore polarization,” Proc. Natl. Acad. Sci. U.S.A. 95 8034-8039 (1998) [CrossRef] [PubMed]

] and molecule motions [3

H. P. Lu, L. Y. Xun, and X. S. Xie, “Single molecule enzymatic dynamics,” Science 282 1877–1882 (1998) [CrossRef] [PubMed]

] is vital to new drug development. Relevant tracking techniques often employ fixed-site fluorescent probe molecules [4

R. E. Dale and S. C. Hopkins “Model-Independent analysis of the orientation of fluorescent probes with restricted mobility in muscle fibers, Biophys. J. 76 1606–1618 (1999) [CrossRef] [PubMed]

] whose properties, such as their orientation and position, vary according to their local environment. These changes can be detected with a suitable experimental setup [5

T. M. Jovin, M. Bartholdi, W. L. C. Vaz, and R. H. Austin “Rotational diffusion of biological macromolecules by time-resolved delayed luminescence (phosphorescence, fluorescence) anisotropy,” Ann. N.Y. Acad. Sci. 366 176–196 (1981) [CrossRef] [PubMed]

, 6

T. Ha, T. Enderle, D. F. Ogletree, D. S. Chemla, P. R. Selvin, and S. Weiss, “Probing the interaction between two single molecules: Fluorescence resonance energy transfer between a single donor and a single acceptor” Proc. Natl. Acad. Sci. U.S.A. 93 6264–6268 (1996). [CrossRef] [PubMed]

Optical techniquesin single molecule imaging almost always require the use of photon counting since individual fluorescent molecules are very weak light sources. Under these conditions the accuracy of measurements are limited by random variations in the measured signal and statistical processing must thus often be used to extract the desired information. This however requires a good understanding of the random processes present. In this work we thus determine the probability density function (PDF) of the number of photons detected during a measurement of finite length as a means of describing the full statistical nature of the detection process.

In the next section we consider the importance of different contributions to photon statistics and find that both noise and statistical signal variations must be considered. As such in Section 3 we derive a general expression for the PDF of photon numbers allowing for both influences. Furthermore in Section 4 we consider one form that signal fluctuations may take, namely those arising from random rotations which fluorophores may undergo; a phenomenon we refer to as wobbling. Wobble of fluorophores, which within the framework of classical electrodynamics can be considered as electric dipole emitters, can be either a continuous angular variation or discrete orientational jumps [7

T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss “Hindered rotational diffusion and rotational jumps of single molecules,” Phys. Rev. Lett. 80 2093–2096 (1998). [CrossRef]

]. We use the Laplace transform to obtain the PDF of the time averaged intensity when changes occur discretely. Successive jumps may depend on previous dipole orientations and we consider the case when this is true and when it is not. Limiting forms for both slow and fast wobble are also given. Further results are then derived for continuous variation in dipole orientation, which is found to be of a differing functional form to the discrete case. A short discussion of the results is finally given in Section 5.

2. Signal-to-noise considerations

Statistical fluctuations in the number of detected photons can derive either from noise present in a system or from random variations in the signal itself. The relative importance of these sources is conveniently parameterised by the signal-to-noise ratio (SNR). In the context of single molecule imaging Basché [8

Th. Basché, W. P. Ambrose, and W. E. Moerner, “Optical spectra and kinetics of single impurity molecules in a polymer: spectral diffusion and persistent spectral hole burning,” J. Opt. Soc. Am. B. 9 829–836 (1992). [CrossRef]

] states that the practically obtainable SNR can be approximated by

SNR = \frac{\frac{DqσP t_{0}}{{AE}_{p}}}{\sqrt{(\frac{DqσP t_{0}}{{AE}_{p}}) + C_{b} P t_{0} + N_{d} t_{0}}}

(1)

where D is an instrument dependent collection factor typically ranging from 1-8%, q is the fluorescence quantum yield, σ is the peak absorption crosssection, P is the laser power, t ₀ is the integration time, A is the beam area, E_p is the energy of a photon in the beam, C_b is the background count per watt of excitation power (typically around 2×10⁸photons/Ws in confocal experiments) and N_d is the dark count of the detector. Figure 1 shows the behaviour of the SNR over a range of experimental conditions from which it can be seen that a value no better than around 15dB is to be expected. Consequently noise properties of the detection process play an important role in determining the statistical behaviour of the detected signal.

Fig. 1. Contour plot of SNR (dB) vs laser power and focused spot size (as parameterised by the NA of a focusing lens) assuming a wavelength of 395nm and the following parameter values based on use of GFP: D = 7%, q = 0.79 [9

G. H. Patterson, S. N. Knobel, W. D. Sharif, S. R. Kain, and D. W. Piston, “Use of the Green Flurorescent Protein and its mutants in quantitative fluorescence microscopy,” Biophys. J. 73 2782–2790 (1997). [CrossRef] [PubMed]

], t ₀ = 0.01s, C_b = 2×10⁸ photons/Ws, N_d = 50. For numerical apertures greater than unity we have assumed an oil immersion lens of refractive index 1.5. Saturation effects are included such that σ = σ₀/(1 +I/I_s ), where σ₀ was taken to be 4×10^-16 cm², I = P/A and the saturation intensity I_s was assumed to be 10³ W/cm² [10

D. J. Pikas, S. M. Kirkpatrick, E. Tewksbury, L. L. Brott, R. R. Naik, M. O. Stone, and W. M. Dennis, “Nonlinear Saturation and Lasing Characteristics of Green Fluorescent Protein,” J. Phys. Chem. B 106 4831–4837 (2002) [CrossRef]

]. The inset shows the variation of the SNR versus integration time for a 5mW laser focused through a 0.95 NA lens.

3. Probability density function of the number of detected photons

Data acquisition in single molecule experiments is invariably done by means of photon counting in which the predominant source of noise is quantisation noise. Denoting the number of photons arriving at the detector during a measurement of duration t ₀ by N(t ₀), the output reading is of the form I_out = GN, where G is some gain factor. The arrival of photons at the detector is a Poisson random process [11

A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).

] and we can hence write the probability mass function of N

p_{N} (n) = \frac{{(I (t_{0}))}^{n}}{n!} \exp (- I (t_{0}))

(2)

where p_N (n) denotes the probability that N(t ₀) = n. We note here that we use the convention whereby an upper case letter denotes a random process and/or variable, whilst the lower case equivalent denotes a particular outcome. I(t ₀) is the average rate of arrival of photons (intensity) or equivalently the time average of ℐ(t) where ℐ(t) is the instantaneous rate of arrival of photons at the detector at time t i.e.

I (t_{0}) = \int_{0}^{t_{0}} ℐ (t) dt

(3)

As an example many experimental setups use polarisation sensitive methods e.g. [12

T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss “Single molecule dynamics studied by polarization modulation,” Phys. Rev. Lett. 77 3979–3982 (1996). [CrossRef] [PubMed]

] whereby the intensity of the detected signal is proportional to the square of the dot product of the illuminating field and the electric dipole moment giving

ℐ (t) = A \cos^{2} (θ (t) - β)

(4)

where θ(t) is the transverse orientation of the dipole at time t, β is the transverse angle of the plane of polarisation of incident light and A is a constant.

For a stationary dipole Eq. (2) fully describes the photon statistics at the detector, however a change in dipole orientation will cause a change in ℐ . If this change is random the arrival of photons at the detector and hence their subsequent detection is termed a doubly stochastic process. Possible sources of such randomness include fluctuations in the illuminating light source and/or movementof the molecule. It is this latter factor that we concentrate on here since for tracking applications the molecule’s environment is unlikely to be static. Furthermore we consider only orientational changes since probe molecules are normally rigidly fixed to targets. Under these circumstances I(t ₀) is a random variable and the probabilities as given by Eq. (2) differ for each possible value. As such we recast Eq. (2) by conditioning the probabilities on a particular outcome i(t ₀) i.e.

p_{N} (n | i) = \frac{i^{n}}{n!} e^{- i}

(5)

where p_N is now a conditional probability and we have dropped the functional dependence on t ₀ for clarity.

Assuming knowledge of the random nature of the time average, as characterised by its PDF fI(i) (see Section 4), we can use the identity [11

A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).

]

f_{N, I} (n, i) = p_{N} (n | i) f_{I} (i)

(6)

to find the joint PDF of N and I i.e. the probability that N = n and I = i. Integrating over the joint PDF gives the PDF of the number of detected photons

p_{N} (n) = \int_{0}^{\infty} f_{I} (i) \frac{{(ηi)}^{n}}{n!} e^{- ηi} di .

(7)

where we have also included the non-ideal nature of the detector by introduction of the quantum efficiency η. Eq. (7

T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss “Hindered rotational diffusion and rotational jumps of single molecules,” Phys. Rev. Lett. 80 2093–2096 (1998). [CrossRef]

) is equivalent to averaging the conditional probability with respect to the average intensity and requires knowledge of f_I (i) which is discussed in the following section.

4. Probability density function of time averaged intensity

4.1. Discrete reorientational jumps

In this section we turn our attention to determining the PDF of the time averaged intensity f_I (i). We first consider the case when changes in the orientation of a dipole occur discretely. This could for example be associated with the desorption and readsorption of fluorophores from and onto a glass surface [12

T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss “Single molecule dynamics studied by polarization modulation,” Phys. Rev. Lett. 77 3979–3982 (1996). [CrossRef] [PubMed]

]. In what follows we shall talk of an electric dipole in an orientational state, by which it is meant that the dipole makes an angle θ to the x-axis in the x-y plane as illustrated in Fig. 2. The dipole then remains fixed at this angle for a time τ before moving to a new state. It is this transverse angle that the signal ℐ in many experimental techniques is dependent on (c.f. Eq. 4). Techniques based on structured illumination and total internal reflection [13

B. Sick, B. Hecht, and L. Novotny “Orientational imaging of single molecules by annular illumination,” Phys. Rev. Lett. 85 4482–4485 (2000). [CrossRef] [PubMed]

, 14

R. M. Dickson, D. J. Norris, and W. E. Moerner, “Simultaneous imaging of individual molecules aligned both parallel and perpendicular to the optic axis,” Phys. Rev. Lett. 81 5322–5325 (1998) [CrossRef]

] do however exist in which the signal depends upon the full three dimensional orientation of a dipole, although these are often restricted to very specific circumstances. Here we restrict our discussion to a two dimensional system for simplicity. Conceptially the full three dimensional situation is identical and requires only minor mathematical modifications as is discussed in the appendix.

Fig. 2. An electric dipole p has a transverse orientation described by the angle θ. Here z describes the optical axis.

Assuming that M different orientational states are occupied during a single measurement the time averaged intensity is given by:

i = A (\cos^{2} (θ_{1} - β) τ_{1} + \cos^{2} (θ_{2} - β) τ_{2} + \dots + \cos^{2} (θ_{M} - β) τ_{M})

(8)

where θ _j and τ _j are the parameters corresponding to the j ^th occupied angular state. Without loss of generality the dipole is assumed to be initially orientated parallel to the x-axis. It should be further noted that changes in the dipole angle are assumed to occur instantaneously.

From the law of rare events [15

K. Itô, Introduction to probability theory (Cambridge University Press, 1984).

] we can assert that M be a Poisson random variable or equivalently that the length of time a dipole remains in each state is distributed according to an exponential law i.e.

f_{τ} (τ) = v \exp (- vτ)

(9)

where f _τ(τ) denotes the PDF of τ and v is the average rate at which dipole jump events occur.

Since a measured intensity is always positive we may use the Laplace transform to find the PDF of I. The Laplace transform X ^*(s) of a random variable X is defined as

X^{*} (s) = ℒ [f_{X} (x)] = \int_{0}^{\infty} f_{X} (x) e^{- sx} dx = E [e^{- sX}]

(10)

where ℒ[…] denotes the Laplace transform and E[…] denotes the expected value.

From [11

A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).

] and Eq. (8) the Laplace transform for fI(i) is given by

I_{M}^{*} (s) = Z_{1}^{*} (s) Z_{2}^{*} (s) \dots Z_{M}^{*} (s)

(11)

where we have let Z_j = Acos²(Θ_j -β)τ_j. The PDF of the average intensity I is then given by:

f_{I} (i) = \sum_{m = 0}^{\infty} p_{M} (m) f_{I} (i | m)

(12)

where f_I (i|m) = ℒ^-1 (I ^* _M=m(s)) and the weighted summation over the possible values of M=m is required since the number of reorientations during a measurement is random.

With this knowledge in hand it remains to find an explicit expression for Z ^* _j(s). From Eqs. (4) and (10) we can write

Z_{j}^{*} (s) = \int_{- π}^{π} \int_{0}^{\infty} f_{Θ, τ} (θ, τ) \exp (- sA \cos^{2} (θ - β) τ) dτdθ

(13)

where f _Θ,τ(θ,τ) is the joint probability distribution of Θ and τ. Since dipole angle and state occupancy time are independent this is given by the product of the marginal probability distributions f _Θ(θ) and f _τ(τ). Using Eq. (9) we can then write

Z_{j}^{*} (s) = \int_{- π}^{π} \int_{0}^{\infty} f_{Θ} (θ) v \exp (- vτ - sA \cos^{2} (θ - β) τ) dτdθ

= \int_{- π}^{π} \frac{v f_{Θ} (θ)}{v + sA \cos^{2} (θ - β)} dθ

(14)

The physical process governing the random wobble of the electric dipole will dictate the form of the probability distribution for Θ. For example, rebinding of a fluorophore to a probe site may be modeled using a uniform PDF

f_{Θ}^{uni} (θ) = {\begin{matrix} 1 / 2 Δ for - Δ \leq θ < Δ \\ 0 otherwise \end{matrix}

(15)

Standard integration tables [16

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, London, 1980).

] then give the analytic result

Z_{uni}^{*} (s) = \frac{1}{2 Δ \sqrt{sAv + v^{2}}} \times

[\arctan (\sqrt{\frac{v}{sA + v}} \tan (Δ - β)) + \arctan (\sqrt{\frac{v}{sA + v}} \tan (Δ + β))]

(16)

Finding a full analytical result for f_I (i) would be complicated however in the limits of small and large v we can find simpler results. These limits correspond to only a few, and to many events per measurement respectively. As the rate at which events occur decreases the contribution from later terms in Eq. (12) becomes negligible. In the limit of v ≪ 1 only the first term produces a significant contribution and we can consider the dipole as fixed during a single measurement and hence

f_{I} (i) = f_{ℐ} (ℐ)

(17)

i.e. the PDF of the average intensity is the same as the PDF for the instantaneous intensity. Fortunately this agrees with our intuitive expectations.

When dipole wobble is on a time scale much shorter than the duration of a measurement we must consider many terms in our summation i.e. large v. Since each value of τ_j is independent each Z_j term is also independent. There are then two cases to consider; that when each subsequent value of θ is independent and that when they are not. In the former case we can invoke the Central Limit Theorem which states the PDF of a sum of independent, identically distributed random variables tends to a Gaussian distribution as the number of terms increases. As such the PDF of the average intensity in the limit of large v is given by

f_{I} (i) = \frac{1}{\sqrt{2 π σ_{I}^{2}}} \exp (- \frac{i^{2}}{2 σ_{I}^{2}}) .

(18)

Assuming dependence of consecutive terms means the PDF of the dipole angle θ_j is centered on its previous outcome, θ _j-1. For a particular realisation of θ, that is to say one possible outcome of the sequence of dipole orientations, we can write

f_{Θ_{j}} (θ_{j}) = f_{Θ} (θ - θ_{j - 1})

(19)

When averaged over all possible realisations the result is similar to Eq. (15

K. Itô, Introduction to probability theory (Cambridge University Press, 1984).

) except now the width of the distribution increases with each subsequent jump. Consequentially the condition of identical distributions required for validity of the Central Limit Theorem is not satisfied. If, however, the Lyapunov condition [17

W. Feller, Probability Theory and its Applications (John Wiley and Sons Inc., New York, 1950).

] is satisfied then the Central Limit Theorem still applies. Numerical simulations show that this is the case.

4.2. Continuous angular variation

Changes in dipole orientation may occur continuously and it is here that we consider how this affects the PDF of the time averaged intensity. It can be shown [18

P. Debye, Polar Molecules (Dover Publications, New York, 1945)

] that the probability distribution function of the orientation of the dipole at a time t satisfies the differential equation:

\frac{\partial f}{\partial t} = α \frac{\partial^{2} f}{\partial θ^{2}}

(20)

subject to the initial condition f(t = 0) = δ(θ-θ ₀), where δ represents the Dirac delta function. This diffusion equation holds when subsequent orientations are dependent on the previous orientation. A solution to Eq. (20

J. Yguerabide, H. F. Epstein, and L. Stryer “Segmental flexibility in an antibody molecule,” J. Mol. Biol. 51 573–590 (1970). [CrossRef] [PubMed]

) is

f_{Θ} (θ, t) = \frac{1}{\sqrt{4 παt}} \exp (- \frac{{(θ - θ_{0})}^{2}}{4 αt}) .

(21)

This result applies to all time scales over which θ may vary as is set by the diffusion coefficient α.

To find the PDF of the average intensity we first use a transformation of variables to find the PDF of the instantaneous intensity fℐ(ℐ) which we then need only integrate over the length of a measurement to give our desired result. Thus

f_{I} (i) = \frac{1}{t_{0}} \sum_{k} \int_{0}^{t_{0}} \frac{f_{Θ} (θ_{k}, t)}{\sqrt{i (A - i)}} dt .

(22)

where θ _k are the solutions to the equation i = Acos²(θ - β) and the 1/t ₀ factor is to ensure correct normalisation of the PDF. The integral can be evaluated using the substitution x ² = t ^-1 and integration by parts which yields:

f_{I} (i) = \frac{1}{\sqrt{πα t_{0} i} (A - i)} \sum_{k} [\exp (- \frac{θ_{k}^{2}}{4 α t_{0}}) - \frac{∣ θ_{k} ∣}{2 α t_{0}} erfc (\sqrt{\frac{θ_{k}^{2}}{4 α t_{0}}})]

(23)

Fig. 3. a) Histogram of the time averaged intensity for a dipole undergoing continuous angular diffusion with α = 5, β = π/4, t ₀ = 10^-3s and A = 10⁵ photons/s shown with theoretical fits for differing diffusion coefficients. b) Variance of the number of detected photons as a function of the peak signal strength A and the cumulative probability functions for wobbly and stationary dipoles (inset) for the same parameter values as a).

where erfc(…) denotes the complimentary error function.

For the independent case f _Θ(θ,t) can not depend on time (assuming the physical cause of the wobble does not vary in time) and as such Eq. (22) reduces to f_I (i) = f_ℐ (ℐ).

Figure 3a) shows a histogram of the result of Monte-Carlo simulations with 10⁴ realisations for continuous variation and a diffusion coefficient of α = 5. Various theoretical fits, as based on Eq. (23), are also drawn from which it can be seen that for α = 0 (no dipole wobble) the PDF is identical to that of a Poisson distribution as would be expected. Good agreement can also be seen between the simulated and theoretical results.

Furthermore, using these PDFs it is possible to calculate the total cumulative probability of N taking any value below n as plotted in the inset of Fig. 3b) as a function of n. Confidence levels including or neglecting dipole wobble can then be calculated. Assuming the values

β = \frac{π}{4}

, θ₀ = 0, A = 10⁵ photons/s, t ₀ = 10^-3s and α = 5 we calculated that when neglecting dipole wobble an experimental measurement can determine the orientation of a dipole within a range of 1.78° with 90% confidence. Inclusion of dipole wobble causes this to increase to 2.43°. Such a discrepency further highlights the need to include dipole wobble in statistical processing and error analysis.

5. Discussion

It has been shown that the variation of the orientation of a dipole over the course of a finite duration measurement can alter the statistical properties of the number of photoelectrons induced in a photon counting detector. Although analytic evaluation of Eqs. (7) and (12) will in general not be possible we can make some general observations regarding the probability functions involved.

Considering first the PDF of the time averaged intensity we can expect different forms and behaviour for differing dependence conditions and time scales. More specifically, although the distribution will always be peaked around the initial angle, when subsequent orientations of the dipole are dependent on earlier positions the distribution is narrower for slower changes, whilst the converse is true when independence holds. This can be understood since the dependent situation is essentially a diffusion problem and so the larger the ratio of spreading rate (as given by the diffusion coefficient) to integration time the larger the range of angles the dipole can cover during a measurement. On the other hand the distribution focuses when successive orientations are independent since the central probability peak for each Z term is reinforced with each additional term in the average.

For wobbling on faster time scales the PDF of the average intensity has been shown to tend to that of a bell-shaped Gaussian distribution (for discrete variations). Slower time scales will exhibit a sharper more centralised distribution, since for small n i.e. slow variation only a few terms significantly contribute to the average performed by the detector. In this case the peaked nature of both the exponentially distributed state occupancy times and the Poisson PDF for the number of events per measurement dominate. For larger v the Poisson PDF becomes smoother and the position of the peak moves to larger m. Low m terms of Eq. (12

T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss “Single molecule dynamics studied by polarization modulation,” Phys. Rev. Lett. 77 3979–3982 (1996). [CrossRef] [PubMed]

) are then negligible and the peaked nature of the exponential PDF is less dominant. Eventually the Poisson PDF tends to a Gaussian itself whereby it acts as an envelope for the PDF of the average dipole angle.

In terms of the photoelectron statistics we can say that for smaller angular ranges of dipole oscillation one would expect less deviation from conventional Poissonian behaviour. Furthermore if the variations are on a timescale much longer than the integration time then the additional random behaviour will be unobservable. On the other hand if fluctuations are much faster than the detector response the effects are likely to again go unnoticed. That said dipole wobble has been seen at many different time scales ranging from the subnanosecond level [19

P. Wahl, K. Tawada, and J.C. Auchet “Study of tropomyosin labelled with a fluorescent probe by pulse fluorimetry in polarized light. Interaction of that protein with troponin and actin,” Eur. J. Biochem. 88 421–424 (1978) [CrossRef] [PubMed]

], through the millisecond regime [5

] and higher [7

T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss “Hindered rotational diffusion and rotational jumps of single molecules,” Phys. Rev. Lett. 80 2093–2096 (1998). [CrossRef]

, 20

J. Yguerabide, H. F. Epstein, and L. Stryer “Segmental flexibility in an antibody molecule,” J. Mol. Biol. 51 573–590 (1970). [CrossRef] [PubMed]

]. In conjunction with the varying time resolution of different experimental setups [21

W. E. Moerner and D. P. Fromm, “Methods of single molecule fluorescence spectroscopy and microscopy,” Rev. Sci. Instrum. 74 3597–3619 (2003) [CrossRef]

] and the large angular ranges over which fluorophores can oscillate e.g. 26° has been observed [22

I. Munro, I. Pecht, and L. Stryer “Subnanosecond motions of Tryptophan residues in proteins,” Proc. Natl. Acad. Sci. USA 76 56–60 (1979) [CrossRef] [PubMed]

], it is likely that non-Poisson behaviour will be frequently encountered.

Finally we draw attention to Fig. 3b) which shows a plot of the expected variance of experimental data as a function of the number of photons in the system (as parameterised by A), when signal variations from photon counting and dipole wobble are considered separately. Quadratic behaviour can be seen for the case of dipole wobble only, whilst for photon counting the linear behaviour expected from a pure Poisson random variable is evident. The relative importance of the two factors can be seen. At very low light intensities, where it is likely to be impractical to conduct experiments, photon counting dominates. For the intermediate regime both influences are comparable until eventually at higher intensities the molecular wobble dominates.

Our work has detailed the statistical behaviour in each of these regimes (Eqs. (2), (7) and (12) or (23) respectively) and allows further development of suitable signal processing algorithms and estimation procedures so as to maximise the precision of angular information extracted from experimental data and thus help to fulfil the potential of single molecule techniques.

Appendices

Appendix - Three dimensional dipole wobble

Earlier discussion has been restricted to the case of two dimensional dipole wobble, however we here explicitly state the mathematical modifications required to accommodate rotation in all three dimensions. Obviously such three dimensional variation is only of significance if the measured signal is sensitive to the full three dimensional orientation as described by the two angles θ and ϕ as shown in Fig. 2 i.e.

ℐ = ℐ (θ, ϕ)

(24)

Random variation in the orientation of the dipole is described by the joint PDF of Θ and Φ. Since in most physical situations Θ and Φ will be independent we can write the PDF in the form

f_{Θ, Φ} (θ, ϕ) = f_{Θ} (θ) f_{Φ} (ϕ)

(25)

When considering the discrete case this means Eq. (13

B. Sick, B. Hecht, and L. Novotny “Orientational imaging of single molecules by annular illumination,” Phys. Rev. Lett. 85 4482–4485 (2000). [CrossRef] [PubMed]

) becomes a triple integral

Z_{j}^{*} (s) = \int_{- π}^{π} \int_{0}^{π} \int_{0}^{\infty} f_{Θ} (θ) f_{Φ} (ϕ) f_{τ} (τ) \exp (- s ℐ (θ, ϕ) τ) dτdϕdθ

(26)

however all of the subsequent working remains unchanged. In the continuous case we must solve the two dimensional diffusion equation

\frac{\partial f}{\partial t} = α (\frac{\partial^{2} f}{\partial θ^{2}} + \frac{\partial^{2} f}{\partial ϕ^{2}})

(27)

to give the joint PDF

f_{Θ, Φ} (θ, ϕ) = \frac{1}{4 παt} \exp (- \frac{{(θ - θ_{0})}^{3}}{4 αt}) \exp (- \frac{{(ϕ - ϕ_{0})}^{2}}{4 αt})

(28)

which when integrated according to the three dimensional analogue of Eq. (22

I. Munro, I. Pecht, and L. Stryer “Subnanosecond motions of Tryptophan residues in proteins,” Proc. Natl. Acad. Sci. USA 76 56–60 (1979) [CrossRef] [PubMed]

) yields

f_{I} (i) = \frac{1}{8 πα t_{0} \sqrt{i (A - i)}} \sum_{k} [Γ (0, \frac{{(θ_{k} - θ_{0})}^{2} + {(ϕ_{k} - ϕ_{0})}^{2}}{4 t_{0} α})]

(29)

where Γ(a,z) = ∫^∞ _z x ^a-1 e ^-x dx is the incomplete Gamma function and θ_k and ϕ_k are the solutions to the equation i = ℐ(θ,ϕ).

References and links

1.	K. D. Weston and L. S. Goldner “Orientation imaging and reorientation dynamics of single dye molecules,” J. Phys. Chem. B 105 3453-3462 (2001) [CrossRef]
2.	D.M. Warshaw, E. Hayes, D. Gaffney, A.M. Lauzon, J.R. Wu, G. Kennedy, K. Trybus, S. Lowey, and C. Berger, “Myosin conformational states determined by single fluorophore polarization,” Proc. Natl. Acad. Sci. U.S.A. 95 8034-8039 (1998) [CrossRef] [PubMed]
3.	H. P. Lu, L. Y. Xun, and X. S. Xie, “Single molecule enzymatic dynamics,” Science 282 1877–1882 (1998) [CrossRef] [PubMed]
4.	R. E. Dale and S. C. Hopkins “Model-Independent analysis of the orientation of fluorescent probes with restricted mobility in muscle fibers, Biophys. J. 76 1606–1618 (1999) [CrossRef] [PubMed]
5.	T. M. Jovin, M. Bartholdi, W. L. C. Vaz, and R. H. Austin “Rotational diffusion of biological macromolecules by time-resolved delayed luminescence (phosphorescence, fluorescence) anisotropy,” Ann. N.Y. Acad. Sci. 366 176–196 (1981) [CrossRef] [PubMed]
6.	T. Ha, T. Enderle, D. F. Ogletree, D. S. Chemla, P. R. Selvin, and S. Weiss, “Probing the interaction between two single molecules: Fluorescence resonance energy transfer between a single donor and a single acceptor” Proc. Natl. Acad. Sci. U.S.A. 93 6264–6268 (1996). [CrossRef] [PubMed]
7.	T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss “Hindered rotational diffusion and rotational jumps of single molecules,” Phys. Rev. Lett. 80 2093–2096 (1998). [CrossRef]
8.	Th. Basché, W. P. Ambrose, and W. E. Moerner, “Optical spectra and kinetics of single impurity molecules in a polymer: spectral diffusion and persistent spectral hole burning,” J. Opt. Soc. Am. B. 9 829–836 (1992). [CrossRef]
9.	G. H. Patterson, S. N. Knobel, W. D. Sharif, S. R. Kain, and D. W. Piston, “Use of the Green Flurorescent Protein and its mutants in quantitative fluorescence microscopy,” Biophys. J. 73 2782–2790 (1997). [CrossRef] [PubMed]
10.	D. J. Pikas, S. M. Kirkpatrick, E. Tewksbury, L. L. Brott, R. R. Naik, M. O. Stone, and W. M. Dennis, “Nonlinear Saturation and Lasing Characteristics of Green Fluorescent Protein,” J. Phys. Chem. B 106 4831–4837 (2002) [CrossRef]
11.	A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).
12.	T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss “Single molecule dynamics studied by polarization modulation,” Phys. Rev. Lett. 77 3979–3982 (1996). [CrossRef] [PubMed]
13.	B. Sick, B. Hecht, and L. Novotny “Orientational imaging of single molecules by annular illumination,” Phys. Rev. Lett. 85 4482–4485 (2000). [CrossRef] [PubMed]
14.	R. M. Dickson, D. J. Norris, and W. E. Moerner, “Simultaneous imaging of individual molecules aligned both parallel and perpendicular to the optic axis,” Phys. Rev. Lett. 81 5322–5325 (1998) [CrossRef]
15.	K. Itô, Introduction to probability theory (Cambridge University Press, 1984).
16.	I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, London, 1980).
17.	W. Feller, Probability Theory and its Applications (John Wiley and Sons Inc., New York, 1950).
18.	P. Debye, Polar Molecules (Dover Publications, New York, 1945)
19.	P. Wahl, K. Tawada, and J.C. Auchet “Study of tropomyosin labelled with a fluorescent probe by pulse fluorimetry in polarized light. Interaction of that protein with troponin and actin,” Eur. J. Biochem. 88 421–424 (1978) [CrossRef] [PubMed]
20.	J. Yguerabide, H. F. Epstein, and L. Stryer “Segmental flexibility in an antibody molecule,” J. Mol. Biol. 51 573–590 (1970). [CrossRef] [PubMed]
21.	W. E. Moerner and D. P. Fromm, “Methods of single molecule fluorescence spectroscopy and microscopy,” Rev. Sci. Instrum. 74 3597–3619 (2003) [CrossRef]
22.	I. Munro, I. Pecht, and L. Stryer “Subnanosecond motions of Tryptophan residues in proteins,” Proc. Natl. Acad. Sci. USA 76 56–60 (1979) [CrossRef] [PubMed]

OCIS Codes
(000.5490) General : Probability theory, stochastic processes, and statistics
(110.4280) Imaging systems : Noise in imaging systems
(180.2520) Microscopy : Fluorescence microscopy

ToC Category:
Imaging Systems

History
Original Manuscript: May 17, 2007
Revised Manuscript: September 28, 2007
Manuscript Accepted: September 28, 2007
Published: October 2, 2007

Virtual Issues
Vol. 2, Iss. 11 Virtual Journal for Biomedical Optics

Citation
Matthew R. Foreman, Sherif S. Sherif, and Peter Török, "Photon statistics in single molecule orientational imaging," Opt. Express 15, 13597-13606 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-21-13597

Sort: Author | Year | Journal | Reset

References

K. D. Weston and L. S. Goldner "Orientation imaging and reorientation dynamics of single dye molecules," J. Phys. Chem. B 105, 3453-3462 (2001). [CrossRef]
D. M. Warshaw, E. Hayes, D. Gaffney, A. M. Lauzon, J. R. Wu, G. Kennedy, K. Trybus, S. Lowey, and C. Berger, "Myosin conformational states determined by single fluorophore polarization," Proc. Natl. Acad. Sci. U.S.A. 95, 8034-8039 (1998). [CrossRef] [PubMed]
H. P. Lu, L. Y. Xun, and X. S. Xie, "Single molecule enzymatic dynamics," Science 282, 1877-1882 (1998). [CrossRef] [PubMed]
R. E. Dale and S. C. Hopkins "Model-Independent analysis of the orientation of fluorescent probes with restricted mobility in muscle fibers," Biophys. J. 76, 1606-1618 (1999). [CrossRef] [PubMed]
T. M. Jovin, M. Bartholdi, W. L. C. Vaz, and R. H. Austin, "Rotational diffusion of biological macromolecules by time-resolved delayed luminescence (phosphorescence, fluorescence) anisotropy," Ann. N. Y. Acad. Sci. 366, 176-196 (1981). [CrossRef] [PubMed]
T. Ha, T. Enderle, D. F. Ogletree, D. S. Chemla, P. R. Selvin and S. Weiss, "Probing the interaction between two single molecules: Fluorescence resonance energy transfer between a single donor and a single acceptor," Proc. Natl. Acad. Sci. U.S.A. 93, 6264-6268 (1996). [CrossRef] [PubMed]
T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss, "Hindered rotational diffusion and rotational jumps of single molecules," Phys. Rev. Lett. 80, 2093-2096 (1998). [CrossRef]
Th. Basché, W. P. Ambrose, and W. E. Moerner, "Optical spectra and kinetics of single impurity molecules in a polymer: spectral diffusion and persistent spectral hole burning," J. Opt. Soc. Am. B. 9, 829-836 (1992). [CrossRef]
G. H. Patterson, S. N. Knobel, W. D. Sharif, S. R. Kain, and D. W. Piston, "Use of the green flurorescent protein and its mutants in quantitative fluorescence microscopy," Biophys. J. 73, 2782-2790 (1997). [CrossRef] [PubMed]
D. J. Pikas, S. M. Kirkpatrick, E. Tewksbury, L. L. Brott, R. R. Naik, M. O. Stone, and W.M. Dennis, "Nonlinear saturation and lasing characteristics of green fluorescent protein," J. Phys. Chem. B 106, 4831-4837 (2002). [CrossRef]
A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).
T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss, "Single molecule dynamics studied by polarization modulation," Phys. Rev. Lett. 77, 3979-3982 (1996). [CrossRef] [PubMed]
B. Sick, B. Hecht, and L. Novotny "Orientational imaging of single molecules by annular illumination," Phys. Rev. Lett. 85, 4482-4485 (2000). [CrossRef] [PubMed]
R. M. Dickson, D. J. Norris, and W. E. Moerner, "Simultaneous imaging of individual molecules aligned both parallel and perpendicular to the optic axis," Phys. Rev. Lett. 81, 5322-5325 (1998). [CrossRef]
K. Itô, Introduction to Probability Theory (Cambridge University Press, 1984).
I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, London, 1980).
W. Feller, Probability Theory and its Applications (John Wiley and Sons Inc., New York, 1950).
P. Debye, Polar Molecules (Dover Publications, New York, 1945).
P. Wahl, K. Tawada, and J. C. Auchet, "Study of tropomyosin labelled with a fluorescent probe by pulse fluorimetry in polarized light. Interaction of that protein with troponin and actin," Eur. J. Biochem. 88, 421-424 (1978). [CrossRef] [PubMed]
J. Yguerabide, H. F. Epstein, and L. Stryer, "Segmental flexibility in an antibody molecule," J. Mol. Biol. 51, 573-590 (1970). [CrossRef] [PubMed]
W. E. Moerner and D. P. Fromm, "Methods of single molecule fluorescence spectroscopy and microscopy," Rev. Sci. Instrum. 74, 3597-3619 (2003). [CrossRef]
I. Munro, I. Pecht, and L. Stryer "Subnanosecond motions of Tryptophan residues in proteins," Proc. Natl. Acad. Sci. USA 76, 56-60 (1979). [CrossRef] [PubMed]

Ambrose, W. P.
Auchet, J. C.
Austin, R. H.
Bartholdi, M.
Basché, Th.
Berger, C.
Brott, L. L.
Chemla, D. S.
Dale, R. E.
Dennis, W.M.
Dickson, R. M.
Enderle, T.
Epstein, H. F.
Fromm, D. P.
Gaffney, D.
Glass, J.
Goldner, L. S.
Ha, T.
Hayes, E.
Hecht, B.
Hopkins, S. C.
Jovin, T. M.
Kain, S. R.
Kennedy, G.
Kirkpatrick, S. M.
Knobel, S. N.
Lauzon, A. M.
Lowey, S.
Lu, H. P.
Moerner, W. E.
Munro, I.
Naik, R. R.
Norris, D. J.
Novotny, L.
Ogletree, D. F.
Patterson, G. H.
Pecht, I.
Pikas, D. J.
Piston, D.W.
Selvin, P. R.
Sharif, W. D.
Sick, B.
Stone, M. O.
Stryer, L.
Tawada, K.
Tewksbury, E.
Trybus, K.
Vaz, W. L. C.
Wahl, P.
Warshaw, D. M.
Weiss, S.
Weston, K. D.
Wu, J. R.
Xie, X. S.
Xun, L. Y.
Yguerabide, J.

Ann. N. Y. Acad. Sci.

T. M. Jovin, M. Bartholdi, W. L. C. Vaz, and R. H. Austin, "Rotational diffusion of biological macromolecules by time-resolved delayed luminescence (phosphorescence, fluorescence) anisotropy," Ann. N. Y. Acad. Sci. 366, 176-196 (1981). [CrossRef] [PubMed]

Biophys. J.

R. E. Dale and S. C. Hopkins "Model-Independent analysis of the orientation of fluorescent probes with restricted mobility in muscle fibers," Biophys. J. 76, 1606-1618 (1999). [CrossRef] [PubMed]
G. H. Patterson, S. N. Knobel, W. D. Sharif, S. R. Kain, and D. W. Piston, "Use of the green flurorescent protein and its mutants in quantitative fluorescence microscopy," Biophys. J. 73, 2782-2790 (1997). [CrossRef] [PubMed]

Eur. J. Biochem.

P. Wahl, K. Tawada, and J. C. Auchet, "Study of tropomyosin labelled with a fluorescent probe by pulse fluorimetry in polarized light. Interaction of that protein with troponin and actin," Eur. J. Biochem. 88, 421-424 (1978). [CrossRef] [PubMed]

J. Mol. Biol.

J. Yguerabide, H. F. Epstein, and L. Stryer, "Segmental flexibility in an antibody molecule," J. Mol. Biol. 51, 573-590 (1970). [CrossRef] [PubMed]

J. Opt. Soc. Am. B.

Th. Basché, W. P. Ambrose, and W. E. Moerner, "Optical spectra and kinetics of single impurity molecules in a polymer: spectral diffusion and persistent spectral hole burning," J. Opt. Soc. Am. B. 9, 829-836 (1992). [CrossRef]

J. Phys. Chem. B

D. J. Pikas, S. M. Kirkpatrick, E. Tewksbury, L. L. Brott, R. R. Naik, M. O. Stone, and W.M. Dennis, "Nonlinear saturation and lasing characteristics of green fluorescent protein," J. Phys. Chem. B 106, 4831-4837 (2002). [CrossRef]
K. D. Weston and L. S. Goldner "Orientation imaging and reorientation dynamics of single dye molecules," J. Phys. Chem. B 105, 3453-3462 (2001). [CrossRef]

Phys. Rev. Lett.

T. Ha, J. Glass, T. Enderle, D. S. Chemla, and S. Weiss, "Hindered rotational diffusion and rotational jumps of single molecules," Phys. Rev. Lett. 80, 2093-2096 (1998). [CrossRef]
T. Ha, T. Enderle, D. S. Chemla, P. R. Selvin, and S. Weiss, "Single molecule dynamics studied by polarization modulation," Phys. Rev. Lett. 77, 3979-3982 (1996). [CrossRef] [PubMed]
B. Sick, B. Hecht, and L. Novotny "Orientational imaging of single molecules by annular illumination," Phys. Rev. Lett. 85, 4482-4485 (2000). [CrossRef] [PubMed]
R. M. Dickson, D. J. Norris, and W. E. Moerner, "Simultaneous imaging of individual molecules aligned both parallel and perpendicular to the optic axis," Phys. Rev. Lett. 81, 5322-5325 (1998). [CrossRef]

Proc. Natl. Acad. Sci. U.S.A.

D. M. Warshaw, E. Hayes, D. Gaffney, A. M. Lauzon, J. R. Wu, G. Kennedy, K. Trybus, S. Lowey, and C. Berger, "Myosin conformational states determined by single fluorophore polarization," Proc. Natl. Acad. Sci. U.S.A. 95, 8034-8039 (1998). [CrossRef] [PubMed]
T. Ha, T. Enderle, D. F. Ogletree, D. S. Chemla, P. R. Selvin and S. Weiss, "Probing the interaction between two single molecules: Fluorescence resonance energy transfer between a single donor and a single acceptor," Proc. Natl. Acad. Sci. U.S.A. 93, 6264-6268 (1996). [CrossRef] [PubMed]

Proc. Natl. Acad. Sci. USA

I. Munro, I. Pecht, and L. Stryer "Subnanosecond motions of Tryptophan residues in proteins," Proc. Natl. Acad. Sci. USA 76, 56-60 (1979). [CrossRef] [PubMed]

Rev. Sci. Instrum.

W. E. Moerner and D. P. Fromm, "Methods of single molecule fluorescence spectroscopy and microscopy," Rev. Sci. Instrum. 74, 3597-3619 (2003). [CrossRef]

Science

H. P. Lu, L. Y. Xun, and X. S. Xie, "Single molecule enzymatic dynamics," Science 282, 1877-1882 (1998). [CrossRef] [PubMed]

Other

A. Leon-Garcia, Probability and Random Processes for Electrical Engineering (Addison-Wesley Publishing Company Inc., 1994).
K. Itô, Introduction to Probability Theory (Cambridge University Press, 1984).
I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, London, 1980).
W. Feller, Probability Theory and its Applications (John Wiley and Sons Inc., New York, 1950).
P. Debye, Polar Molecules (Dover Publications, New York, 1945).

2003, Moerner, Rev. Sci. Instrum.
2002, Pikas, J. Phys. Chem. B
2001, Weston, J. Phys. Chem. B
2000, Sick, Phys. Rev. Lett.
1999, Dale, Biophys. J.
1998, Warshaw, Proc. Natl. Acad. Sci. U.S.A.
1998, Lu, Science
1998, Ha, Phys. Rev. Lett.
1998, Dickson, Phys. Rev. Lett.
1997, Patterson, Biophys. J.
1996, Ha, Proc. Natl. Acad. Sci. U.S.A.
1996, Ha, Phys. Rev. Lett.
1992, Basché, J. Opt. Soc. Am. B.
1981, Jovin, Ann. N. Y. Acad. Sci.
1979, Munro, Proc. Natl. Acad. Sci. USA
1978, Wahl, Eur. J. Biochem.
1970, Yguerabide, J. Mol. Biol.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.

« Previous Article | Next Article »

OSA is a member of CrossRef.

Select a Conference
Session or Paper #
Year

OpticsInfoBase

Optics Express

Photon statistics in single molecule orientational imaging

Article Tools

Article Navigation

Abstract

1. Introduction

2. Signal-to-noise considerations

3. Probability density function of the number of detected photons

4. Probability density function of time averaged intensity

4.1. Discrete reorientational jumps

4.2. Continuous angular variation

5. Discussion

Appendices

Appendix - Three dimensional dipole wobble

References and links

References

Ann. N. Y. Acad. Sci.

Biophys. J.

Eur. J. Biochem.

J. Mol. Biol.

J. Opt. Soc. Am. B.

J. Phys. Chem. B

Phys. Rev. Lett.

Proc. Natl. Acad. Sci. U.S.A.

Proc. Natl. Acad. Sci. USA

Rev. Sci. Instrum.

Science

Other

Cited By

Figures