Improving Multicolor Colocalization in Single-Vesicle Flow Cytometry with Vesicle Transit Time.

Immunophenotyping of vesicles, such as extracellular vesicles (EVs), is essential to understanding their origin and biological role. We previously described a custom-built flow analyzer that utilizes a gravity-driven flow, high numerical aperture objective, and micrometer-sized flow channels to reach the sensitivity needed for fast multidimensional analysis of the surface proteins of EVs, even down to the smallest EVs (e.g., ∼30-40 nm). It is difficult to flow focus small EVs, and thus, the transiting EVs exhibit a distribution in particle velocities due to the laminar flow. This distribution of vesicle velocities leads to potentially incorrect results when immunophenotyping nanometer-sized vesicles using cross-correlation analysis (Xcorr), as the order of appearance of the vesicles might not be the same at different spatially offset laser excitation regions. Here, we describe an alternative cross-correlation analysis strategy (Scorr), which uses information on particle transit time across the laser excitation beam width to improve multicolor colocalization in single-vesicle immunoprofiling. We tested the performance of the algorithm for colocalization analysis of multicolor nanobeads and EVs experimentally and via simulations and found that Scorr improved both the efficiency and accuracy of colocalization versus Xcorr. As shown from Monte Carlo simulations, Scorr provided an ∼1.2-4.7-fold increase in the number of colocalized peaks and ensured negligible colocalization of peaks. In silico results were in good agreement with experimental data, which showed an increase in colocalized peaks of ∼1.3-2.5-fold and ∼1.2-2-fold for multicolor beads and EVs, respectively.

Sizing Extracellular Vesicles Using Membrane Dyes and a Single Molecule-Sensitive Flow Analyzer.

Extracellular vesicles (EVs) are membranous particles released by most cells in our body, which are involved in many cell-to-cell signaling processes. Given the nanometer sizes and heterogeneity of EVs, highly sensitive methods with single-molecule resolution are fundamental to investigating their biophysical properties. Here, we demonstrate the sizing of EVs using a fluorescence-based flow analyzer with single-molecule sensitivity. Using a dye that selectively partitions into the vesicle's membrane, we show that the fluorescence intensity of a vesicle is proportional to its diameter. We discuss the constraints in sample preparation which are inherent to sizing nanoscale vesicles with a fluorescent membrane dye and propose several guidelines to improve data consistency. After optimizing staining conditions, we were able to measure the size of vesicles in the range ∼35-300 nm, covering the spectrum of EV sizes. Lastly, we developed a method to correct the signal intensity from each vesicle based on its traveling speed inside the microfluidic channel, by operating at a high sampling rate (10 kHz) and measuring the time required for the particle to cross the laser beam. Using this correction, we obtained a threefold greater accuracy in EV sizing, with a precision of ±15-25%.

Detection of 14 High-Risk Human Papillomaviruses Using Digital LAMP Assays on a Self-Digitization Chip.

Cervical cancer is the fourth-leading cause of cancer deaths among women worldwide and most cases occur in developing countries. Detection of high-risk (HR) HPV, the etiologic agent of cervical cancer, is a primary screening method for cervical cancer. However, the current gold standard for HPV detection, real-time PCR, is expensive, time-consuming, and instrumentation-intensive. A rapid, low-cost HPV detection method is needed for cervical cancer screening in low-resource settings. We previously developed a digital loop-mediated isothermal amplification (dLAMP) assay for rapid, quantitative detection of nucleic acids without the need for thermocycling. This assay employs a microfluidic self-digitization chip to automatically digitize a sample into an array of nanoliter wells in a simple assay format. Here we evaluate the dLAMP assay and self-digitization chip for detection of the commonly tested 14 high-risk HPVs in clinical samples. The dLAMP platform provided reliable genotyping and quantitative detection of the 14 high-risk HPVs with high sensitivity, demonstrating its potential for simple, rapid, and low-cost diagnosis of HPV infection.

Statistical Analysis of Nonuniform Volume Distributions for Droplet-Based Digital PCR Assays.

We present a method to determine the concentration of nucleic acids in a sample by partitioning it into droplets with a nonuniform volume distribution. This digital PCR method requires no special equipment for partitioning, unlike other methods that require nearly identical volumes. Droplets are generated by vortexing a sample in an immiscible oil to create an emulsion. PCR is performed, and droplets in the emulsion are imaged. Droplets with one or more copies of a nucleic acid are identified, and the nucleic acid concentration of the sample is determined. Numerical simulations of droplet distributions were used to estimate measurement error and dynamic range and to examine the effects of the total volume of droplets imaged and the shape of the droplet size distribution on measurement accuracy. The ability of the method to resolve 1.5- and 3-fold differences in concentration was assessed by using simulations of statistical power. The method was validated experimentally; droplet shrinkage and fusion during amplification were also assessed experimentally and showed negligible effects on measured concentration.

Single-Molecule Flow Platform for the Quantification of Biomolecules Attached to Single Nanoparticles.

We describe here a flow platform for quantifying the number of biomolecules on individual fluorescent nanoparticles. The platform combines line-confocal fluorescence detection with near nanoscale channels (1-2 µm in width and height) to achieve high single-molecule detection sensitivity and throughput. The number of biomolecules present on each nanoparticle was determined by deconvolving the fluorescence intensity distribution of single-nanoparticle-biomolecule complexes with the intensity distribution of single biomolecules. We demonstrate this approach by quantifying the number of streptavidins on individual semiconducting polymer dots (Pdots); streptavidin was rendered fluorescent using biotin-Alexa647. This flow platform has high-throughput (hundreds to thousands of nanoparticles detected per second) and requires minute amounts of sample (∼5 µL at a dilute concentration of 10 pM). This measurement method is an additional tool for characterizing synthetic or biological nanoparticles.

Squaraine-based polymer dots with narrow, bright near-infrared fluorescence for biological applications.

This article describes the design and development of squaraine-based semiconducting polymer dots (Pdots) that show large Stokes shifts and narrow-band emissions in the near-infrared (NIR) region. Fluorescent copolymers containing fluorene and squaraine units were synthesized and used as precursors for preparing the Pdots, where exciton diffusion and likely through-bond energy transfer led to highly bright and narrow-band NIR emissions. The resulting Pdots exhibit the emission full width at half-maximum of ∼36 nm, which is ∼2 times narrower than those of inorganic quantum dots in the same wavelength region (∼66 nm for Qdot705). The squaraine-based Pdots show a high fluorescence quantum yield (QY) of 0.30 and a large Stokes shift of ∼340 nm. Single-particle analysis indicates that the average per-particle brightness of the Pdots is ∼6 times higher than that of Qdot705. We demonstrate bioconjugation of the squaraine Pdots and employ the Pdot bioconjugates in flow cytometry and cellular imaging applications. Our results suggest that the narrow bandwidth, high QY, and large Stokes shift are promising for multiplexed biological detections.

Supercoiled pseudocircular domains in single-twisted DNAs under tension: Elastic constants and unwinding dynamics in complexes with Topo I.

Extension versus twist data of Koster et al. (Nature 2005, 434, 671-674) are analyzed to obtain C for the main-chain segments and the twist energy parameter (ET ) for the supercoiled pseudocircular (sp) domain(s) from which C is estimated via simulations. The torsional rigidity in the tension-free sp domain(s) (C = 163 fJ fm) is typical of the unstrained DNA and is less than half the value in the main-chain segments under tension (C = 350-410 fJ fm). Tension is suggested to induce a structural transition to a torsionally stiffer state. Data of Koster et al. for the rate of extension owing to unwinding of a covalent complex of DNA with human Topoisomerase Ib (H Topo I) are analyzed to determine the torque and rate of rotation from which an effective friction coefficient is obtained. A Langevin equation for the unwinding motion in a supercoiled DNA:H Topo I complex is solved to obtain the temporal trajectory of the average winding angle and the time-dependent distribution of winding angles. The mean rate constant for the religation reaction is estimated from the measured probability of reaction per turn. We predict that unwinding proceeds rather far during a single-cleavage and religation cycle, and is effectively completely equilibrated during the 3.2 cleavage and religation cycles that occur during each noncovalent binding and dissociation event. H Topo I is predicted to be completely processive as in accord with observations on calf-thymus Topo I (Brewood et al., Biochemistry 2010, 49, 3367-3380).

Protein quantification at the single vesicle level reveals that a subset of synaptic vesicle proteins are trafficked with high precision.

Protein sorting represents a potential point of regulation in neurotransmission because it dictates the protein composition of synaptic vesicles, the organelle that mediates transmitter release. Although the average number of most vesicle proteins has been estimated using bulk biochemical approaches (Takamori et al., 2006), no information exists on the intervesicle variability of protein number, and thus on the precision with which proteins are sorted to vesicles. To address this, we adapted a single molecule quantification approach (Mutch et al., 2007) and used it to quantify both the average number and variance of seven integral membrane proteins in brain synaptic vesicles. We report that four vesicle proteins, SV2, the proton ATPase, Vglut1, and synaptotagmin 1, showed little intervesicle variation in number, indicating they are sorted to vesicles with high precision. In contrast, the apparent number of VAMP2/synaptobrevin 2, synaptophysin, and synaptogyrin demonstrated significant intervesicle variability. These findings place constraints on models of protein function at the synapse and raise the possibility that changes in vesicle protein expression affect vesicle composition and functioning.

Single-molecule fluorescence quantification with a photobleached internal standard.

In cellular and molecular biology, fluorophores are employed to aid in tracking and quantifying molecules involved in cellular function. We previously developed a sensitive single-molecule quantification technique to count the number of proteins and the variation of the protein number over the population of individual subcellular organelles. However, environmental effects on the fluorescent intensity of fluorophores can make it difficult to accurately quantify proteins using these sensitive techniques. In this letter, we demonstrate the use of photobleaching to extract an accurate single-molecule calibration intensity distribution from the sample directly to avoid any differences in environment that may alter the count. Using this technique, we were able to show that goat antimouse IgG antibody labeled with Alexa Fluor 488, an environmentally insensitive fluorophore, exhibited an average fluorescence equivalent to 4.6 single fluorophores. SynaptopHluorin vesicles, which contain the environmentally sensitive green fluorescent protein, exhibited an average of 4.4 single green fluorescent proteins per vesicle.

Measurements of the acidification kinetics of single SynaptopHluorin vesicles.

Uptake of neurotransmitters into synaptic vesicles is driven by the proton gradient established across the vesicle membrane. The acidification of synaptic vesicles, therefore, is a crucial component of vesicle function. Here we present measurements of acidification rate constants from isolated, single synaptic vesicles. Vesicles were purified from mice expressing a fusion protein termed SynaptopHluorin created by the fusion of VAMP/synaptobrevin to the pH-sensitive super-ecliptic green fluorescent protein. We calibrated SynaptopHluorin fluorescence to determine the relationship between fluorescence intensity and internal vesicle pH, and used these values to measure the rate constant of vesicle acidification. We also measured the effects of ATP, glutamate, and chloride on acidification. We report acidification time constants of 500 ms to 1 s. The rate of acidification increased with increasing extravesicular concentrations of ATP and glutamate. These data provide an upper and a lower bound for vesicle acidification and indicate that vesicle readiness can be regulated by changes in energy and transmitter availability.

Large structural change in isolated synaptic vesicles upon loading with neurotransmitter.

The size of a synaptic vesicle (SV) is generally thought to be determined by the amount of lipid and membrane protein it contains. Once formed, it is thought to remain constant in size. Using fluorescence correlation spectroscopy and cryogenic electron microscopy, we show that glutamatergic vesicles reversibly increase their size upon filling with glutamate. The increase ( approximately 25% in diameter) corresponds to an increase in surface area of approximately 50% and in volume of approximately 100%. This large size increase implies a large structural change in the SV upon loading with neurotransmitters. Vesicles lacking SV protein 2A (SV2A) did not manifest a change in size after loading with glutamate, indicating that SV2A is required for this phenomenon.

Sizing subcellular organelles and nanoparticles confined within aqueous droplets.

This article describes two complementary techniques, single-particle tracking and correlation spectroscopy, for accurately sizing nanoparticles confined within picoliter volume aqueous droplets. Single-particle tracking works well with bright particles that can be continuously illuminated and imaged, and we demonstrated this approach for sizing single fluorescent beads. Fluorescence correlation spectroscopy detects small intensity bursts from particles or molecules diffusing through the confocal probe volume, which works well with dim and rapidly diffusing particles or molecules; we demonstrated FCS for sizing synaptic vesicles confined in aqueous droplets. In combination with recent advances in droplet manipulations and analysis, we anticipate this capability to size single nanoparticles and molecules in free solution will complement existing tools for probing cellular systems, subcellular organelles, and nanoparticles.

Estimation of the persistence length of DNA from the torsion elastic constant and supercoiling free energy: effect of ethylene glycol.

Two different methods are proposed to estimate the persistence length ( P) of DNA from the measured torsion elastic constant (alpha) and the twist energy parameter ( E T ) that governs the supercoiling free energy. The first method involves Monte Carlo simulations and reversible-work calculations of E T for model DNAs that possess the measured alpha and selected trial values of P. Comparison of the computed E T values with the experimental value allows estimation of P (or equivalently the bending elastic constant (kappa beta)) by interpolation. A far simpler, though less accurate, alternative is to solve a previously conjectured analytical relation connecting E T , alpha, kappa beta (or P), and an unknown "constant" ( B). The present simulations are used to ascertain the optimum value of B and to assess the validity and accuracy of that relation. Within the simulation errors, P values obtained from the measured alpha and E T via this analytical expression agree with those determined from the simulations and E T values reckoned from the input alpha and kappa beta by this analytical expression agree with the corresponding simulated values. Although B is found to be insensitive to variation in alpha, it appears to decline slightly with increasing kappa beta. The original analytical expression is modified to take this apparent variation of B with kappa beta into account. By using this modified analytical relation to estimate P (from the measured alpha and E T ) or E T (from the input alpha and kappa beta), much closer agreement is obtained respectively with the values of P or E T obtained from the simulations. As specific examples, these methods are applied to determine P in 0 and 20 w/v % ethylene glycol, which has been shown to induce a structural transition in duplex DNA.

Self-digitization chip for single-cell genotyping of cancer-related mutations.

Cancer is a heterogeneous disease, and patient-level genetic assessments can guide therapy choice and impact prognosis. However, little is known about the impact of genetic variability within a tumor, intratumoral heterogeneity (ITH), on disease progression or outcome. Current approaches using bulk tumor specimens can suggest the presence of ITH, but only single-cell genetic methods have the resolution to describe the underlying clonal structures themselves. Current techniques tend to be labor and resource intensive and challenging to characterize with respect to sources of biological and technical variability. We have developed a platform using a microfluidic self-digitization chip to partition cells in stationary volumes for cell imaging and allele-specific PCR. Genotyping data from only confirmed single-cell volumes is obtained and subject to a variety of relevant quality control assessments such as allele dropout, false positive, and false negative rates. We demonstrate single-cell genotyping of the NPM1 type A mutation, an important prognostic indicator in acute myeloid leukemia, on single cells of the cell line OCI-AML3, describing a more complex zygosity distribution than would be predicted via bulk analysis.

Optical gradient flow focusing.

This paper describes a new method for carrying out flow cytometry, which employs optical gradient forces to guide and focus particles in the fluid flow. An elliptically shaped Gaussian beam was focused at the center of a microchannel to exert radiation pressure on suspended nanoparticles that are passing through the channel, such that these particles are guided to the center of the channel for efficient detection and sorting. To verify the efficiency of this optical-gradient-flow-focusing method, we present numerical simulations of the trajectories of the nanoparticles in both electroosmotic flow (EOF) and pressure-driven flow (PDF).

Quantitative microscopy based on single-molecule fluorescence.

Quantitative microscopy is needed to understand reactions or phenomena carried out by biological molecules such as enzymes, receptors, and membrane-localized proteins. Counting the biomolecules of interest in single organelles or cellular compartments is critical in these approaches. In this brief perspective, we focus on the development of quantitative fluorescence microscopies that measure the precise copy numbers of proteins in cellular organelles or purified samples. We introduce recent improvements in quantitative microscopies to overcome undercounting or overcounting errors in certain conditions. We conclude by discussing biological applications.

Accurate sizing of nanoparticles using confocal correlation spectroscopy.

The ability to accurately size low concentrations of nanoscale particles in small volumes is useful for a broad range of disciplines. Here, we characterize confocal correlation spectroscopy (CCS), which is capable of measuring the sizes of both fluorescent and nonfluorescent particles, such as quantum dots, gold colloids, latex spheres, and fluorescent beads. We accurately measured particles ranging in diameter from 11 to 300 nm, a size range that had been difficult to probe, owing to a phenomenon coined biased diffusion that causes diffusion times, or particle size, to deviate as a function of laser power. At low powers, artifacts mimicking biased diffusion are caused by saturation of the detector, which is especially problematic when probing highly fluorescent or highly scattering nanoparticles. However, at higher powers (>1 mW), autocorrelation curves in both resonant and nonresonant conditions show a structure indicative of an increased contribution from longer correlation times coupled with a decrease in shorter correlation times. We propose that this change in the autocorrelation curve is due to the partial trapping of the particles as they transit the probe volume. Furthermore, we found only a slight difference in the effect of biased diffusion when comparing resonant and nonresonant conditions. Simulations suggest the depth of trapping potential necessary for biased diffusion is > 1 k(B)T. Overcoming artifacts from detector saturation and biased diffusion, CCS is particularly advantageous due to its ability to size particles in the small volumes characteristic of microfluidic channels and aqueous microdroplets. We believe the method will find increasing use in a wide range of applications in measuring nanoparticles and macromolecular systems.

Can reliable torsion elastic constants be determined from FPA data on 24 and 27 base-pair DNAs?

Torsion elastic constants obtained from fluorescence polarization anisotropy (FPA) measurements on fifty-three 24 and 27 base-pair (bp) DNAs were recently reported [F. Pedone, F. Mazzei, D. Santoni, Sequence-dependent DNA torsional rigidity: a tetranucleotide code, Biophys. Chem. 112 (2004) 77-88; F. Pedone, F. Mazzei, M. Matzeu, F. Barone, Torsional constant of 27-mer DNA oligomers of different sequences, Biophys. Chem. 94 (2001) 175-184]. The problem of extracting reliable torsion elastic constants (alpha) from FPA measurements on such short DNAs is examined in detail. The difficulty is illustrated by two (fictitious) 24 bp DNAs with approximately 5-fold different torsion elastic constants and 10% different initial anisotropies (r(0)), which exhibit practically indistinguishable anisotropy decays for all t>1 ns. FPA data were simulated for 24 bp DNAs with different input values of alpha and r(0) in the presence and absence of Poisson noise, and were fitted using different choices of the adjustable and fixed parameters. Experimental data for a 24 bp DNA were fitted in a similar manner. For either the simulated or experimental FPA data, it was not possible to determine both the initial anisotropy, r(0), and the torsion elastic constant, alpha, in a reliable (i.e. statistically significant) manner in the presence of Poisson noise. When r(0) is assumed to be fixed at any particular value in the fitting protocol, a unique best-fit value of alpha is obtained, but that best-fit alpha is extremely sensitive to small deviations of the assumed fixed value of r(0) away from the input r(0)-value of the simulated data. Pedone et al. fitted their FPA data by assuming that r(0)=0.360, and adjusting alpha, the hydrodynamic radius (R(H)), and effective length (L). In fact, the reported best-fit values of R(H) and L lay significantly outside their expected ranges. When this same fitting protocol is applied to simulated data for 27 bp DNAs, better overall agreement with the reported experimental values (alpha, R(H), and L) is obtained for a model, wherein all DNAs have the same typical input alpha=5.9 x 10(-12) dyn cm, R(H)=10.0 A, and L=27 (3.4)+2.7=94.5 A, but a 1.00- to 1.13-fold range of r(0)-values, than for the model of Pedone et al., wherein all DNAs have the same input r(0)=0.360, R(H)=10.0 A, and L=94.5 A, but a approximately 3-fold range of alpha-values. It is concluded that, in the absence of reliable independent estimates of r(0) for every DNA, the alpha-values reported for 24 and 27 bp DNAs cannot be regarded as experimentally justified. The reliability of the torsion elastic constants reported for the 136 distinct tetranucleotide steps, which are inferred from the values reported for the fifty-three 24 and 27 bp DNAs, is also briefly discussed.

Dynamic twisting correlations in a model DNA with uniform torsion elastic constant.

The extent to which the twisting motions of two separate subunits (base pairs) in an elastic filament are correlated is discussed in terms of the two-point correlation function of their azimuthal angular displacements C(Delta,t) identical with <(phi(m)(t) - phi(m)(0))(phi(n)(t) - phi(n)(0))>, where m,n are the sequential subunit indices and Delta = |m - n| is their absolute difference. An approximate expression is derived for C(Delta,t) for an infinitely long model DNA from the analytical theory developed previously to treat the decay of the fluorescence polarization anisotropy of intercalated ethidium. C(Delta,t) is numerically evaluated as a function of Delta for a range of times (20, 40, 60, 80, 100, and 120 ns) for a model DNA with a typical torsion elastic constant. By t = 120 ns, significant dynamic correlations are observed to extend over a domain (Delta) several hundred base pairs. Copyright 1999 John Wiley & Sons, Inc.

Extensions of counterion condensation theory. I. Alternative geometries and finite salt concentration.

The counterion condensation theory originally proposed by Manning is extended to take account of both finite counterion concentration (m(C)) and the actual structure of the array of discrete changes. Counterion condensation is treated here as a binding isotherm problem, in which the unknown free volume is replaced by an unknown local binding constant beta', which is expected to vary with m(C) and polyion structure. The relation between the condensed fraction of counterion charge, r, beta' and m(C) is obtained from the relevant grand partition function via the maximum term method. In the case of the single polyion in a large salt reservoir, the result is practically identical to Manning's equation. In order to determine the values of beta' and r at arbitrary m(C), a second relation between r, beta' and m(C) is required. We propose an alternative auxiliary relation that is equivalent to previous assumptions near m(C) = 0, but which yields qualitatively correct and quantitatively useful results at finite m(C). Simple expressions for r vs. m(C) and beta' vs. m(C) are obtained by simultaneously solving the binding isotherm and auxiliary equations. Then r and beta' are evaluated for five different linear arrays of infinite extent with different geometries: (1) a straight line of charges with uniform axial spacing; (2) two parallel lines of in-phase uniformly spaced charges; (3) a single-helix of discrete charges with uniform axial spacing; (4) a double-helix of discrete charges with uniform axial spacing of pairs of charges; (5) a cylindrical array of many parallel charged lines, chosen to simulate a uniformly charged cylinder. In all cases, the computed binding isotherms exhibit qualitatively correct behavior. As m(C) approaches zero, r approaches the Manning limit, r = 1-1/(L(B)/b) where b is the average axial spacing of electronic charges in the array and L(B) is the Bjerrum length. However, beta' varies with polyion geometry, even in the zero salt limit, and matches the Manning value only in the case of a single straight charged line. With increasing m(C), r declines significantly below its limiting value whenever lambda(b) > or approximately equal 0.3, where lambda is the Debye screening parameter. In the case of cylindrical arrays containing either 2 or 100 parallel charged lines, r also decreases, whenever lambda(d) > or approximately equal 2.0, where d is the diameter of the array. In the case of two parallel charged lines, each with axial charge spacing b=3.4 A, which are separated by d = 200 A, r exhibits a plateau value, 0.76, characteristic of the two combined lines, when lambda(d)<<2.0, and declines with increasing m(C) to a shelf value, 0.52, characteristic of either single line, when lambda(d) > or approximately equal 2.0 and the lines become effectively screened from one another. beta' behaves in a roughly similar fashion. In the case of a cylindrical array of charged lines with the diameter and linear charge density of DNA, the r-values predicted by the present theory agree fairly well with those predicted by non-linear Poisson-Boltzmann theory up to 0.15 M uni-univalent salt.