Previous Article |
Table of Contents
| Next Article
Biophysics
Single-molecule high-resolution imaging with photobleaching


*Department of Physics and
Center for Biophysics and Computational Biology, University of Illinois at UrbanaChampaign, Urbana, IL 61801
Communicated by Gordon A. Baym, University of Illinois at UrbanaChampaign, Urbana, IL, March 8, 2004 (received for review January 29, 2004)
| Abstract |
|---|
|
|
|---|
200 nm. Some recent techniques have been developed to circumvent this limitation (3, 4), but they are technically demanding and at this point of limited applicability to biological systems. Fluorescence resonance energy transfer (FRET) can be used to measure distances much smaller, on the order of a few nanometers. However, because of the strong distance dependence of energy transfer, FRET is limited to measuring distances less than
10 nm (5). As a result, there is a "gap" in the resolution attainable by single-color optical spectroscopy, making it difficult to measure separations between 10 and 200 nm. Many biological objects of interest are on this scale, including DNA structures, macromolecular complexes, and motor proteins. Although it is possible to measure distances on these scales by using two or more dyes of different colors (6, 7), this technique presents its own problems, such as how to achieve heterogeneous labeling, and calibration of the distance registration between different color wide-field images. Recently, with the introduction of low-noise high-quantum-yield charge-coupled device (CCD) cameras, it has become practical to localize individual fluorescent dyes at subwavelength scales (810). We recently demonstrated the ability to determine the positions of single molecules at room temperature with a precision of 1.5 nm (11) by using total internal reflection (TIR) microscopy and a cooled CCD. This technique has been used to probe the motility of Myosin V and kinesin (11, 12). Here we present a technique, single-molecule high-resolution imaging with photobleaching (SHRImP), that takes advantage of the quantal photobleaching of fluorescent molecules to resolve two identical fluorophores in the xy plane separated by distances as small as 10 nm, with a precision for single measurements on the order of 5 nm. In SHRImP, two molecules, m1 and m2, with overlapping point-spread functions (PSFs) are imaged continuously, creating a series of images that forms a "movie." A plot of the integrated intensity of the combined PSF shows a two-step signal reduction (Fig. 1a), each step corresponding to the photobleaching of an individual molecule. After the first photobleaching event, only m2 is left, and its PSF can be fit to a Gaussian distribution, allowing us to localize m2 to within 1.5 nm. Now, by selecting one image from before the photobleaching event (Ipre) and one from after (Ipost), we may compute the distribution of m1 alone, Ipre', by taking Ipre - Ipost (Fig. 1b). We may then use this distribution to localize m1, and hence compute the separation between m1 and m2. Although this is a conceptually simple procedure, a more technically involved procedure involves localizing the PSFs simultaneously, by performing a global analysis involving both Ipre and Ipost, fitting Ipost to a single Gaussian, and fitting Ipre to two Gaussians, as discussed below. We used the latter method for the data presented but found that, in practice, both procedures produced nearly identical results.
|
| Materials and Methods |
|---|
|
|
|---|
Sample Preparation. Biotinylated Cy3-DNA samples were prepared as in ref. 11. Briefly, coverslips were cleaned by sonication in acetone, then in KOH, and rinsed between baths. Coverslips were incorporated into flow chambers, and Cy3-DNA was attached via a biotin-streptavidin linkage. The singly labeled DNA sequence was the same as reported in ref. 11. Flow chambers for doubly labeled DNA samples (Integrated DNA Technologies, Coralville, IA) were prepared similarly, but coverslips were silanized by using a procedure derived from ref. 14. Coverslips were cleaned as above and then incubated in 6 mM 3-aminopropyltriethoxysilane (APTES, Sigma-Aldrich), pH
3, at 60°C for 17.5 h. They were then rinsed with water and stored in distilled ethanol. Doubly labeled DNA samples were diluted to
100 pM in T50 buffer (10 mM Tris/1 mM EDTA/50 mM NaCl, pH 7) and flowed directly over the silanized surface. They were allowed to sit for 5 min then flushed with T50 buffer, followed by a flush with deoxygenating imaging buffer. The imaging buffer consisted of 1% by volume gloxy, 1% by volume 2-mercaptoethanol, 0.4% by mass glucose, and 15 mM MgCl2 in T50. Gloxy consisted of 30 mg of glucose oxidase (G-7016, Sigma-Aldrich) and 60 µl of catalase solution (106810, Roche Diagnostics) in 200 µl of T50 buffer, spun down and filtered twice with syringe filters. Doubly labeled DNA consisted of two hybridized complementary DNA strands, each covalently labeled on the 5' end with a Cy3 phosphoramidite dye (Amersham Biosciences). The doubly labeled 51-mer DNA sequence was 5'-TGG ACT CAG AGC ATT CAA GAT GGT CGG TGG ACG GTT GAC CTA CGG CTG CCA, hybridized with its complement. The 40- and 30-mer sequences were the same but truncated from the 3' end and modified so that the last three bases were CCA.
Data Analysis. Biotinylated Cy3-DNA images were analyzed by performing a global fit with SIGMAPLOT 8.0 (SPSS, Chicago), using single pre- and postbleach images selected by hand. Pre- and postbleach images Ipre and Ipost were analyzed by globally fitting both distributions simultaneously to two-dimensional Gaussians. The underlying distribution is expected to be an Airy disk, but it is common to model it with a Gaussian, because it is more mathematically tractable, and because the differences are expected to be minor in practice (10, 11). For fitting, we used the following functions:
![]() | [1a] |
and
![]() | [1b] |
where A and B are the heights of the two Gaussians, x0 and y0 are the coordinates of the centroid of the first Gaussian,
x and
y are the x and y components of the vector separating the centroid of the second Gaussian from the first, wx1, wy1, wx2, and wy2 are the widths of the x and y components of the two Gaussians, and z0 and z1 are the z-axis offsets of the distributions, arising from the sum of the background noise, dark noise, and readout noise from the CCD, as well as any background fluorescence. The separation of the two molecules was computed from
.
To improve throughput, doubly labeled DNA images were analyzed by using custom-written software in IDL (Research Systems, Boulder, CO). Each image frame consisted of several DNA molecules. The software picked individual punctuate objects, and from these, molecules showing a clean, two-step intensity profile similar to Fig. 1a and good fit to a Gaussian were chosen for further analysis. Fluorescence intensity stability was variable among molecules, so it was not possible to determine where photobleaching events occurred for some observed molecules. Some spots also showed many more than two photobleaching events, indicating some "clumping" of the sample. However, overall, 3040% of the observed spots showed clear two-step photobleaching and could be analyzed further.
Each molecule was analyzed in the following way. For values of N ranging from 1 to 4, 2N images from before the first bleach step and 2N images from after this step were chosen. N images from each set were summed, resulting in four composite images, and each of the two postbleach images were paired with each of the two prebleach images to generate a set of four different pairs of images. We refer to this set as a "quartet." Four separations (d) were computed from the quartet by performing a global fit of both the measured prebleach distribution (Ipre) and the measured postbleach distribution (Ipost) simultaneously. The pre- and postbleach fit functions were similar to those in Eq. 1, but each point in the distribution was integrated in a square region over the pixel size of the camera, reflecting the fact that each measured datum corresponds to all of the photons falling on a particular pixel of finite size.
The data quality was evaluated according to the following three criteria. The standard deviation of d within each quartet must be <10 nm, the ellipticity of each PSF must be <20%, and each composite image must have enough photons to theoretically allow localization with at least 1.7-nm precision. The ellipticity of the PSF is
= 1 - wminor/wmajor, where wminor and wmajor are the major and minor axes that form the waist of the Gaussian distribution. From the quartets that passed these tests for each value of N, the one with the smallest standard deviation was chosen, and the mean of the quartet was used as a data point for computing the final histogram.
Of the above cited criteria, the ellipticity requirement is intended to filter out dyes with a high degree of immobility. Our observations of Cy3 dried on coverslips shows that highly immobile dyes have PSFs with much higher ellipticity. It has been demonstrated that very small deviations in the image focus can lead to large systematic errors in the localization of dyes with highly immobile emission dipoles (15). Dyes that change angle rapidly on the time scale of imaging (16), however, are not subject to this systematic error. In practice, however, this criterion was almost never violated when observing dyes in solution.
Dye intensity is variable because of variations in illumination, blinking, and local environment. Allowing up to four sequential images to be summed increases the likelihood of being able to localize the molecule to within the required accuracy. However, the more images that are binned together, the worse the time resolution becomes. Also, the further in time that we are from the photobleaching point, the greater the likelihood is that the dye intensity has changed. But, the SHRImP analysis depends on the assumption that the intensity remains constant over time, because we assume the postphotobleach image of m2 is a good approximation of its contribution to the prephotobleach image. Hence, using larger bins can decrease the accuracy if the intensity fluctuations are large and occur on the time scale of the camera integration time. Taking two sets of images from both before and after the photobleaching event and using them to compute four different distances allows us to use the standard deviation of the quartet distances as a measure of the degree to which intensity fluctuations impact the measurement. Quartets showing a large standard deviation may be due to intensity fluctuations which are too large to yield good measurements.
After all molecules were analyzed, the computed separations were used to generate a histogram. The histogram was fit to a Gaussian plus a constant offset, of the form
![]() | [2] |
The centroid of this Gaussian is the measured value of the separation, with a margin of error given by the standard error of the mean as determined by the fit. The constant offset allows for "background" due to closely spaced but disconnected molecules.
| Results and Discussion |
|---|
|
|
|---|
µ) of 13.5 nm.
µ was determined by SIGMAPLOT based on the computed fit to the fit functions.
|
150 bp (17)], and hence are expected to be relatively rigid and straight.
We found that, in practice, a large percentage of the observed molecules (
70%) did not show good two-step photobleaching behavior. In some cases, this was because they showed only a single photobleaching event or showed more than two photobleaching events, or because the photobleaching transitions were obscured by intensity fluctuations. However, because wide-field illumination and CCD imaging allows highly paralleled data collection, we were easily able to generate enough data for analysis.
Of the molecules showing good two-step behavior, more than two thirds were found to pass all three of the filtering criteria. Almost all of the data passed the ellipticity criterion, and the quartet standard deviation was
3 nm for those that passed the criteria.
Fig. 3 shows histograms of the distance measurements for the three DNA lengths. A single peak is seen in each histogram, and the peak distance increases for increasing DNA length. The histograms were each fit to a one-dimensional Gaussian with an offset (Eq. 2). The offset allows for a "background" due to closely spaced but disconnected molecules, because of the fact that there is no a priori way to differentiate between two covalently linked dye molecules and two that simply lie close to each other. The density of observed molecules was
0.3 µm-2, low enough that we would expect a small incidence of such observations. That such measurements were an observable fraction may be due in part to the clumping discussed in Materials and Methods. It also may be due to individual molecules that were deformed because of interaction with the surface, or to data that were not fit well. The Gaussian fits to the histograms shown in Fig. 3 have widths ranging from 4 to 6 nm, where the "width" refers to the parameter w in Eq. 2. This parameter would be the standard deviation of the measurement in the absence of the "background." Our data filtering is intended to discard molecules that cannot be localized to at least 1.7 nm, and therefore, ideally, we would expect to see a width close to
, or 2.4 nm. The larger observed width may reflect the impact of two factors: first, the experimentally achieved accuracy is often worse than 1.7 nm because of aberrations in the optical system, and is sometimes as high as 2.5 nm. Second, there is a contribution to the width due to the uncertainty in the intensity of the m1 dye because of shot noise and "blinking" (18). As noted above, the average quartet standard deviation was
3 nm. This standard deviation is probably a reasonable measure of the contribution of dye intensity fluctuation to the overall histogram width, and these numbers taken together may account for the magnitude of the observed width.
|
|
| Acknowledgements |
|---|
| Footnotes |
|---|
To whom correspondence should be addressed. E-mail: selvin{at}uiuc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
A. Agrawal, R. Deo, G. D. Wang, M. D. Wang, and S. Nie Nanometer-scale mapping and single-molecule detection with color-coded nanoparticle probes PNAS, March 4, 2008; 105(9): 3298 - 3303. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Zhang, D. Shu, F. Huang, and P. Guo Instrumentation and metrology for single RNA counting in biological complexes or nanoparticles by a single-molecule dual-view system RNA, October 1, 2007; 13(10): 1793 - 1802. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. R. Selvin, T. Lougheed, M. Tonks Hoffman, H. Park, H. Balci, B. H. Blehm, and E. Toprak Fluorescence Imaging with One-Nanometer Accuracy (FIONA) CSH Protocols, October 1, 2007; 2007(20): pdb.top27 - pdb.top27. [Abstract] [Full Text] |
||||
![]() |
M. Bates, B. Huang, G. T. Dempsey, and X. Zhuang Multicolor Super-Resolution Imaging with Photo-Switchable Fluorescent Probes Science, September 21, 2007; 317(5845): 1749 - 1753. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. E. Moerner Single-molecule Chemistry and Biology Special Feature: New directions in single-molecule imaging and analysis PNAS, July 31, 2007; 104(31): 12596 - 12602. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Mutch, B. S. Fujimoto, C. L. Kuyper, J. S. Kuo, S. M. Bajjalieh, and D. T. Chiu Deconvolving Single-Molecule Intensity Distributions for Quantitative Microscopy Measurements Biophys. J., April 15, 2007; 92(8): 2926 - 2943. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sharonov and R. M. Hochstrasser Wide-field subdiffraction imaging by accumulated binding of diffusing probes PNAS, December 12, 2006; 103(50): 18911 - 18916. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Heintzmann and G. Ficz Breaking the resolution limit in light microscopy Brief Funct Genomic Proteomic, December 1, 2006; 5(4): 289 - 301. [Abstract] [Full Text] [PDF] |
||||
![]() |
T.-F. Chan, C. Ha, A. Phong, D. Cai, E. Wan, L. Leung, P.-Y. Kwok, and M. Xiao A simple DNA stretching method for fluorescence imaging of single DNA molecules Nucleic Acids Res., October 18, 2006; 34(17): e113 - e113. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. C. Lagerholm, L. Averett, G. E. Weinreb, K. Jacobson, and N. L. Thompson Analysis Method for Measuring Submicroscopic Distances with Blinking Quantum Dots Biophys. J., October 15, 2006; 91(8): 3050 - 3060. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Betzig, G. H. Patterson, R. Sougrat, O. W. Lindwasser, S. Olenych, J. S. Bonifacino, M. W. Davidson, J. Lippincott-Schwartz, and H. F. Hess Imaging Intracellular Fluorescent Proteins at Nanometer Resolution Science, September 15, 2006; 313(5793): 1642 - 1645. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. R. Brau, P. B. Tarsa, J. M. Ferrer, P. Lee, and M. J. Lang Interlaced Optical Force-Fluorescence Measurements for Single Molecule Biophysics Biophys. J., August 1, 2006; 91(3): 1069 - 1077. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Michalet and S. Weiss Using photon statistics to boost microscopy resolution PNAS, March 28, 2006; 103(13): 4797 - 4798. [Full Text] [PDF] |
||||
![]() |
S. Ram, E. S. Ward, and R. J. Ober Beyond Rayleigh's criterion: A resolution measure with application to single-molecule microscopy PNAS, March 21, 2006; 103(12): 4457 - 4462. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. S. Churchman, H. Flyvbjerg, and J. A. Spudich A Non-Gaussian Distribution Quantifies Distances Measured with Fluorescence Localization Techniques Biophys. J., January 15, 2006; 90(2): 668 - 671. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Robinson, J. P.A. Wood, S. A. Capaldi, A. J. Baron, C. Gell, D. A. Smith, and N. J. Stonehouse Affinity of molecular interactions in the bacteriophage {varphi}29 DNA packaging motor. Nucleic Acids Res., January 1, 2006; 34(9): 2698 - 2709. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Balci, T. Ha, H. L. Sweeney, and P. R. Selvin Interhead Distance Measurements in Myosin VI via SHRImP Support a Simplified Hand-Over-Hand Model Biophys. J., July 1, 2005; 89(1): 413 - 417. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. J. SMITH, T. R. SOSNICK, N. F. SCHERER, and T. PAN Efficient fluorescence labeling of a large RNA through oligonucleotide hybridization RNA, February 1, 2005; 11(2): 234 - 239. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. S. Churchman, Z. Okten, R. S. Rock, J. F. Dawson, and J. A. Spudich Single molecule high-resolution colocalization of Cy3 and Cy5 attached to macromolecules measures intramolecular distances through time PNAS, February 1, 2005; 102(5): 1419 - 1423. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Qu, D. Wu, L. Mets, and N. F. Scherer Nanometer-localized multiple single-molecule fluorescence microscopy PNAS, August 3, 2004; 101(31): 11298 - 11303. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||