What measure is used to characterize the peak tall thin nature of distribution?

Introduction

Efforts to develop new drugs are not limited to the physicochemical properties of pharmaceuticals. They also include explorations of effective ways to deliver those drugs without compromising efficacy or safety1,2,3,4,5,6,7,8,9,10. Despite advances in molecular biology research, many drugs still have serious side effects due to the lack of a specific target and correct control release profile, and these side effects limit our ability to design optimal medications for many diseases, including cancer, neurodegenerative diseases and infectious diseases11,12,13,14,15. To address this issue, researchers have developed several new modes of drug delivery system (DDS) that have entered clinical practice, including nanoparticles based on polymers, noble metals and lipid based carriers3,16,17,18,19. The interactions and stability of such materials are strongly dependent on carrier size, whose characterization is crucial in assessing the quality and determining the efficiency of the DDS20,21,22. In particular, chemical modification of nanoparticles is necessary to make them suitable for physiological conditions, and accurate measurement of their size is necessary for quality control23,24,25. This requirement has become more significant as nanoparticles, and their chemical modifications, have been developed for more specific purposes1,26,27. Likewise, lipid vesicles, due to their versatile engineering capabilities, have been combined with various therapeutic agents to achieve desired pharmaceutical properties28,29. Due to the inherent self-assembly of lipids, validation of their size and distribution is essential to understand the physical properties that directly correlate with drug efficacy. All of these factors highlight the importance of using accurate and precise measurement techniques to characterize the size distribution of biological and synthetic nanoparticle suspensions.

In the analysis of macromolecular assemblies, various techniques are used to measure the physical properties of samples, including imaging30,31, separation of particles32,33, scattered light34,35 and those measurements are related to the size by conversions relying on various physical principles36. Direct imaging techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM) and atomic force microscopy (AFM), are some of the most popular methods to obtain the topographical size of particles, as well as their shape and texture31,36,37,38,39 Imaging has been preferred due to its intuitive high-resolution visualization of particles and the minimal influence of artifacts in size determination36. However, imaging methods require laborious sample preparation steps, and the sample must be removed from its native or working environment, often resulting in a deformation to the samples. In addition, throughput is limited and limited sampling may result in biased information36.

Another strategy is to separate the particles in the sample, creating a spatial macromolecular redistribution in a medium, in which the degree of separation is determined by the mass or volume of the macromolecules and can be converted into their size36. This approach is a feature of various techniques, including size exclusion chromatography (SEC), asymmetrical flow field-flow fractionation (AF4) and analytical ultracentrifugation (AUC), which measure differences in the elution, sedimentation or diffusion of particles33,40,41. As the particles in the sample are spatially separated depending on their differences in the course of measurement, these techniques can be combined with other size measurement techniques, such as multi-angle light scattering (MALS), to improve the size resolution or measure additional properties, such as molecular weight42,43. As the techniques involve separation of the measured sample, they provide more useful size measurements of polydisperse samples but introduce distortion of the sample condition due to the medium used36.

Among non-destructive measurement techniques44, dynamic light scattering (DLS) is the most widely used due to its simplicity. Upon laser illumination, the intensity of the light scattered by the particles in suspension changes both temporally and spatially depending on the size and weight of the particles and can be converted into size information34,36. Despite being a powerful and accessible tool, DLS has several drawbacks due to the inherent limitation of intensity-biased detection34,45. DLS determines a particle’s size from fluctuations of the scattered light resulted from the Brownian motion of the particles. The intensity of the scattered light is proportional to the square of the volume of the particle, which makes DLS very sensitive to the presence of large particles45,46. Small amounts of large aggregates or dust particles can disturb the size determination if the main population is significantly smaller in size.

To address this problem, other techniques including AFM, SEM, TEM, AF4 and AUC can be used to verify the size determined by DLS38,47,48,49. However, the definition of size measured by one technique can be different from that of the others, which makes the comparison complicated as it requires careful interpretation of the data obtained. In particular, DLS measures the hydrodynamic diameter from the diffusion coefficient of particles in suspension, and this diffusion coefficient is converted into the diameter of an assumed hypothetical sphere that has the same diffusion coefficient34. In contrast, SEM, for instance, can obtain a geometric size of particles given by measuring the width of individual particles from the image36,50. Recently, a size characterization tool called nanoparticle tracking analysis (NTA) was introduced to acquire the size of particles by determining their diffusion coefficient, meaning that the definition of the size measured by NTA is identical to that of DLS. NTA can be a good method to verify the results of DLS because they measure the same physical property51,52. Whereas DLS reads the intensity change of scattered light to find the diffusion coefficient of particles, NTA calculates the diffusion coefficient based on the movements of individual particles in successive optical video images51,52. This difference in the detection principles of DLS and NTA results in a difference in the way that they report size, i.e., the quantities of the particles measured by DLS and NTA are weighted by intensity and number, respectively, which makes NTA an excellent technique for verifying DLS results.

Initial studies on NTA focused on the validation of NTA measurements of mono- and multimodal nanoparticle samples and compared the results with DLS47,49,52,53. These studies confirmed that NTA provides comparable results in determining the size of mono- and multimodal nanoparticles compared with DLS. The validation was extended to comparisons with AFM47, TEM38 and AUC49 techniques, where the distribution reported by NTA was a good match for those from the techniques on monomodal samples. When the comparison was extended to polydisperse samples such as proteins and vesicles52,54,55,56, the results showed better size resolution in NTA compared with DLS.

While NTA showed better size resolution than DLS on polydisperse samples, the studies did not observe a narrow distribution from NTA on monodisperse samples. The analysis software for NTA used in the previous studies acquired the size distribution by directly converting the displacement of tracks into the size, which is prone to a stochastic error in determining the size from the particle tracking. Although the studies recognized the uncertainty resulting from the stochastic error47,49,56, they used the size distributions from the direct conversion due to the lack of methods to mitigate the stochastic error. To solve this issue, Saveyn et al. and Walker suggested size estimation methods based on the maximum likelihood estimation (MLE) principle57,58. These methods search a narrower distribution that maximizes the likelihood on the acquired track data either by assuming a certain size distribution with a few adjustable parameters or by taking iterative steps, and can successfully recover a narrower size distribution of monomodal samples from NTA track data. Using this strategy to recover a better size distribution from NTA, monomodal reference samples were measured to compare the results with TEM30. Further, Kestens et al. compared versions 2.3 and 3.0 of NTA software, which use the direct conversion and the iterative approach, respectively59, and found that the iterative approach has a better size resolution on mono- and multi-modal reference samples.

Inspired by previous studies52,57,58,59, this study aimed to evaluate NTA results with respect to the size estimation methods for NTA. First, monodisperse reference standards of known size values were measured using NTA and compared with DLS. The size analysis of acquired particle tracks of NTA was processed by direct conversion of the tracks and with the two MLE-based methods, namely finite track length adjustment (FTLA) and the iterative method. The comparison confirmed that both of the two MLE-based methods can recover a narrow size distribution of the monodisperse reference standards. The estimation methods were also applied to estimate the size distribution of polydisperse vesicle samples prepared in various polydispersity conditions to verify if FTLA and the iterative method are applicable for polydisperse samples of unknown size distribution. While the iterative method achieved results comparable with DLS, the results from FTLA deviated from those of DLS as it assumes an arbitrary size distribution that may not be appropriate for polydisperse samples. In addition, the size distribution of NTA acquired with the estimation methods was converted into an intensity-weighted size distribution from its number-weighted distribution to investigate if the conversion gives a better comparison with the DLS results.

Results and Discussion

Size Measurement of Polystyrene Latex Nanoparticle Standards

To validate NTA and its size distribution analysis methods, monodisperse polystyrene (PS) latex standards were measured by DLS and NTA. The mean size and the polydispersity index (PI) of the DLS results were acquired by cumulant analysis of the measured auto-correlation function and are shown in Fig. 1. The mean size obtained by DLS is in good agreement with the nominal values, with 91 ± 1, 278 ± 1 and 352 ± 3 nm acquired for the 92, 269 and 343 nm particle size standards, respectively. The PIs of the samples were 0.037 ± 0.029, 0.016 ± 0.009 and 0.035 ± 0.018, respectively, indicating that the samples were monodisperse.

Figure 1

(a) Mean and (b) relative variance of the size distributions of PS latex nanoparticle samples measured by DLS and NTA. Mean size is displayed as a value relative to its respective nominal sample size. The size information from DLS was acquired using the cumulant method, while the results from NTA were acquired by applying three different methods, i.e., direct conversion of acquired tracks (DC), FTLA and the iterative method (IC).

Full size image

Size information from NTA, shown in Fig. 1, was extracted from the track data with segment length greater than 5, and processed by the direct conversion method. The mean sizes of the 92, 269 and 343 nm PS standards were 94 ± 4, 261 ± 3 and 320 ± 3 nm, respectively, close to their respective nominal values. The RVs of the 92, 269 and 343 nm PS standards were 0.075 ± 0.009, 0.033 ± 0.007 and 0.044 ± 0.001, confirming that the size distributions were monodisperse.

To determine if the RV of the NTA results from the direct conversion could be reduced by considering a limited track segment length, the two MLE-based size estimation methods, i.e., FTLA and the iterative method, were applied to the particle tracks acquired from the standard samples. The size distributions of the 92 nm standard sample acquired by direct conversion, FTLA and the iterative method are shown in Fig. 2. The distributions of the other standard samples are presented in the Supporting Information. As shown in Fig. 1, mean sizes of 88 ± 4, 250 ± 3 and 294 ± 7 nm acquired by FTLA and 89 ± 4, 250 ± 4 and 295 ± 6 nm acquired by the iterative method were determined for the 92, 269 and 343 nm standard samples, respectively. RVs of 0.019 ± 0.003, 0.010 ± 0.008 and 0.052 ± 0.022 by FTLA and 0.038 ± 0.003, 0.041 ± 0.010 and 0.058 ± 0.002 by the iterative method were acquired for the 92, 269 and 343 nm standard samples, respectively. The RVs of the 256 and 343 nm standards acquired by the FTLA and iterative methods are not very different from those derived by direct conversion, whereas that of the 92 nm standard is significantly reduced. This suggests that the MLE-based methods are effective for size distribution estimation especially when the particle size is small, which tends to result in a limitation in acquiring tracks of high enough track segment length due to the large diffusion coefficient. For small particle samples, the error in size estimation by direct conversion is enlarged as the error is inversely proportional to the square root of the track segment length, making the observed size distribution broad57. For large particle samples, the particles in the sample are tracked for long enough that the error in size estimation by direct conversion is as small as that of the MLE-based methods.

Figure 2

The size distributions of 92-nm PS latex standard measured by DLS and NTA. The size distribution from DLS assumes a log-normal distribution acquired from the cumulant method, while those from NTA are from the three different analysis methods, i.e., direct conversion (DC), FTLA and the iterative method (IC).

Full size image

In the comparison, the mean sizes of the standard samples determined by the FTLA and iterative methods match each other well despite their different approaches to finding the maximum likelihood. However, the RVs from FTLA are smaller than the corresponding values from the iterative method. The assumed size distribution of FTLA has the advantage of a smoother size distribution but at the expense of reduced likelihood60, resulting in the smaller RVs of these PS standard sample measurements. Given that the samples are monodisperse, assuming a log-normal distribution for the size distribution estimation does not decrease the likelihood compared with the corresponding result from the iterative method (see Supporting Information).

POPC Vesicles

POPC vesicles prepared in various polydispersity conditions were measured by DLS. The results are presented in Fig. 3. The mean size of the vesicle samples increased with increasing pore size of the extrusion filter. For vesicle samples extruded through filters with pore sizes of 50 nm (V50), 100 nm (V100), 200 nm (V200) and 400 nm (V400), DLS revealed mean sizes of 91 ± 4, 132 ± 3, 213 ± 4 and 458 ± 13 nm, respectively, and PIs of 0.089 ± 0.032, 0.064 ± 0.035, 0.077 ± 0.048 and 0.291 ± 0.013, respectively. Despite the increase in mean size, the PI does not vary much for pore sizes of 50, 100 and 200 nm, indicating that these samples are relatively monodisperse. However, it jumps to about 0.3 for the 400 nm pore size, showing that the vesicle sample extruded through 400 nm pores is highly polydisperse.

Figure 3

(a) Mean and relative variance of the size distributions of the POPC vesicle samples measured by DLS and NTA and (b) relative mean size from NTA with respect to their corresponding mean size from DLS. DLS values are derived from the cumulant method, while the NTA results are obtained by applying the three different analysis methods, i.e., direct conversion (DC), FTLA and the iterative method (IC).

Full size image

For the vesicles processed by freeze and thaw (FT) treatment before being extruded through a 400-nm pore filter, the mean size shows a gradual decrease with increasing number of cycles. As shown in Fig. 3, for 3, 5, 7, 9, 11, 13, 15 and 17 FT cycles (FT3 to FT17), the mean sizes were 282 ± 4, 223 ± 2, 185 ± 2, 161 ± 1, 148 ± 1, 131 ± 1, 117 ± 2 and 109 ± 1 nm, respectively. The PI results were 0.273 ± 0.021, 0.224 ± 0.022, 0.197 ± 0.029, 0.149 ± 0.023, 0.181 ± 0.024, 0.168 ± 0.015, 0.189 ± 0.025 and 0.197 ± 0.021, respectively, indicating that the polydispersity of the samples reduced upon successive FT cycles but did not decrease enough for the sample to be regarded as monodisperse.

The NTA measurements of the same samples were analyzed by using direct conversion, which determined mean sizes of 100 ± 2, 119 ± 9, 211 ± 6 and 329 ± 17 nm and RVs of 0.149 ± 0.059, 0.151 ± 0.082, 0.097 ± 0.002 and 0.136 ± 0.037 for the vesicle samples extruded through pores of 50, 100, 200 and 400 nm, respectively, and mean sizes of 243 ± 12, 225 ± 9, 201 ± 7, 184 ± 8, 171 ± 13, 158 ± 15, 135 ± 4 and 137 ± 1 nm and RVs of 0.171 ± 0.024, 0.288 ± 0.029, 0.223 ± 0.016, 0.206 ± 0.009, 0.273 ± 0.024, 0.305 ± 0.059, 0.245 ± 0.010 and 0.274 ± 0.019 for the vesicle samples treated by 3, 5, 7, 9, 11, 13, 15 and 17 FT cycles, respectively, before being extruded through a 400 nm filter, as shown in Fig. 3.

For the different extrusion filter pore sizes, the mean sizes measured by NTA were similar to those from DLS, increasing with increasing filter pore size. However, the RVs acquired from NTA indicated that the samples were polydisperse and did not significantly vary with filter pore size, unlike the PI values obtained by DLS. For the various numbers of FT cycles, the trend of mean size measured by NTA also matched that obtained by DLS, decreasing gradually as the number of cycles increased. However, the RVs of the samples from NTA show a gradual increase as the number of FT cycles increases, which is opposite to the trend of the PI measured by DLS, probably because the FT treatment led to homogenization of the vesicle samples61. Remarkably, for some of the samples the mean size acquired from NTA was larger than its corresponding mean size from DLS, in contrast to the belief that the mean size of an intensity-weighted distribution is larger than that of its corresponding number-weighted distribution45,48. As the polydispersity of the samples is very large, it is more unlikely that a larger mean size will be obtained from NTA.

As direct conversion is not appropriate in determining standard deviation in the measurements of the PS latex standards, the two MLE-based methods were applied for comparison. The results are presented in Fig. 3. Using the FTLA method, NTA obtained mean sizes of 87 ± 5, 106 ± 5, 193 ± 5 and 275 ± 14 nm and RVs of 0.083 ± 0.049, 0.068 ± 0.045, 0.051 ± 0.001 and 0.283 ± 0.033 for the filter pore sizes of 50, 100, 200 and 400 nm, respectively, and mean sizes of 212 ± 13, 195 ± 6, 173 ± 5, 160 ± 5, 138 ± 8, 131 ± 8, 113 ± 3 and 113 ± 1 nm and RVs of 0.147 ± 0.018, 0.203 ± 0.013, 0.135 ± 0.007, 0.112 ± 0.006, 0.178 ± 0.021, 0.170 ± 0.040, 0.128 ± 0.005 and 0.137 ± 0.006 for 3, 5, 7, 9, 11, 13, 15 and 17 FT cycles, respectively. The iterative method determined mean sizes of 88 ± 5, 108 ± 4, 193 ± 5 and 279 ± 14 nm and RVs of 0.103 ± 0.055, 0.105 ± 0.050, 0.064 ± 0.002 and 0.212 ± 0.031 for the filter pore size of 50, 100, 200 and 400 nm, respectively, and mean sizes of 216 ± 11, 196 ± 6, 173 ± 5, 160 ± 5, 138 ± 8, 131 ± 8, 113 ± 4 and 112 ± 1 nm and RVs of 0.193 ± 0.012, 0.249 ± 0.019, 0.176 ± 0.010, 0.145 ± 0.009, 0.228 ± 0.024, 0.243 ± 0.043, 0.178 ± 0.011 and 0.180 ± 0.011 for 3, 5, 7, 9, 11, 13, 15 and 17 FT cycles, respectively. The size distribution profiles of the vesicle samples analyzed by the different methods for NTA measurements are presented in the Supporting Information.

Regardless of the condition of the vesicle samples, the mean sizes of the standard sample measurements estimated by the two MLE-based methods match very well despite the different approaches of the two methods. However, the RVs from the two methods are different, with the RVs from the FTLA method being smaller than those from the iterative method except for the 400-nm filter pore size condition. Although FTLA finds a smaller RV than the iterative method, as noted above, it is at the cost of allowing a smaller likelihood value by assuming a certain shape for the size distribution, i.e., a log-normal distribution, for parameter optimization60. In fact, the likelihood value of the size distribution acquired through FTLA is smaller than that from the iterative method except in one case (see Supporting Information). Moreover, some of the likelihood values are even smaller than their corresponding values from the direct conversion, i.e., without any statistical processing. This implies that FTLA is inappropriate for the estimation of size distributions of NTA measurements on polydisperse samples without a proper assumption of the sample size distribution and that the iterative method is preferable to FTLA for such samples54.

By applying the iterative method as the preferred method for the size distribution estimation of NTA measurements, the size distributions of the vesicle samples acquired by DLS and NTA are compared in Fig. 4. A comparison of the distribution profiles shows that the size distribution of the vesicle samples is monomodal for most of the samples except V400 and the distribution from NTA is shifted toward smaller sizes than its corresponding distribution from DLS. This shift can be attributed to the different physical quantity that the two measurement techniques acquire, i.e., the intensity of scattered light and number of particles for DLS and NTA, respectively, indicating a larger shift for higher polydispersity45,48. However, the shift (i.e., the difference in the mean sizes measured by DLS and NTA), is very small compared with the acquired PI and RV of the samples even after applying the iterative method for better size estimation of NTA.

Figure 4

Size distribution of the POPC vesicle samples prepared in various polydispersity conditions measured by DLS and NTA. DLS results (solid red line) are from the cumulant method while NTA results (solid blue line) are acquired by applying the iterative method. For comparison, the size distribution of NTA is reconstructed into an intensity-weighted distribution (dotted black line) by introducing the thin-shell sphere model for the conversion. The vesicle samples were prepared by extrusion through (a) 50-nm, (b) 100-nm, (c) 200-nm and (d) 400-nm pore filters or by freeze-thaw treatment of (e-l) 3 to 17 cycles before extrusion through a 400-nm pore filter.

Full size image

Filipe et al. proposed that the small shift can be explained by the lower detection limit of DLS, which detects small particles below 30 nm, which also explains the relatively high PI value measured by DLS52. It expects a smaller RV for a size distribution acquired by NTA than the corresponding PI by DLS since NTA would neglect the small size range from the true size distribution and report a narrower distribution. However, the results of FT11 and FT13 show that their respective RV from NTA is larger than their corresponding PI from DLS while the mean size of each sample measured by NTA and DLS is similar, which is contradictory to the expectation.

Theoretically, the small difference in the mean sizes measured by DLS and NTA can be due to the complex nature of the light scattering intensity in relation to the scattering form factor of the sample48. Therefore, the number-weighted size distribution acquired by NTA was reconstructed into an intensity-weighted distribution (as shown in Fig. 4) by assuming an isotropic thin-shell hollow sphere model for the vesicle samples and introducing the RGD approximation to construct the form factor48,55,62. The reconstruction revealed mean sizes of 102 ± 1, 118 ± 9, 191 ± 4 and 279 ± 46 nm and RVs of 0.065, 0.056, 0.033 and 0.271 for V50, V100, V200 and V400, respectively, and mean sizes of 209 ± 17, 201 ± 7, 173 ± 4, 165 ± 3, 147 ± 6, 148 ± 5, 131 ± 2 and 130 ± 3 nm and RVs of 0.116 ± 0.045, 0.169 ± 0.021, 0.090 ± 0.020, 0.076 ± 0.002, 0.113 ± 0.019, 0.128 ± 0.002, 0.110 ± 0.003 and 0.102 ± 0.013 for the vesicle samples treated with 3 to 17 FT cycles, respectively. Figure 5 compares the mean sizes and RVs of the original number-weighted and the converted intensity-weighted distributions from NTA with those from DLS. The mean sizes measured by NTA and those acquired from its reconstructed intensity-weighted size distribution are very close to each other. However, the RVs from the reconstructed NTA size distributions are significantly low for some samples, but it is questionable if those samples can be considered monodisperse based only on the reconstructed values. This result indicates that converting the NTA size distribution to an intensity-weighted distribution does not effectively mitigate the difference in the mean sizes observed by DLS and NTA and that it is better to compare the measurements by DLS and NTA when the size distributions are expressed in their original weightings.

Figure 5

(a) Mean size and (b) relative variance of the POPC vesicle samples measured by NTA compared with DLS (open red square). NTA results were analyzed by the iterative method (solid black triangle) and reconstructed into an intensity-weighted distribution (open black triangle) using the thin-shell sphere model for the conversion.

Full size image

Conclusion

This paper introduces and compares two particle size measurement techniques, DLS and NTA, for the size characterization of polydisperse macromolecular assemblies. While both techniques acquire size information by detecting the diffusion coefficient of the measured particles, they differ in the quantity of the size distribution they report, with NTA reporting number and DLS reporting intensity. Three size distribution estimation methods for NTA were tested to compare their performance using monomodal PS latex standards. The two MLE-based methods, i.e., the FTLA and iterative methods, produced results comparable to each other and in line with those from DLS while the direct conversion resulted in a larger variance, especially for small particles, due to the limitation of obtaining sufficiently long particle tracks. This result indicates that an MLE-based approach should be applied for the accurate measurement of small sized particles with NTA.

The two MLE-based size distribution estimation methods for NTA were further tested on measurements of polydisperse vesicle samples prepared in various polydispersities, which obtained the same mean sizes despite their different strategies to finding the optimal size distribution from the given track data. However, the calculated likelihood of the acquired size distributions obtained by the two methods indicated that FTLA sacrificed likelihood at the expense of a smoother size distribution, further suggesting that the iterative method is preferable for polydisperse samples.

The results for the vesicle samples obtained by NTA using the iterative method were then compared with those from DLS. The mean sizes were comparable except for those samples with a very high polydispersity index, i.e., V400 and FT3. However, the mean sizes of some samples, e.g., FT9 to FT17, measured by NTA were larger than or very close to those from DLS, which seems contradictory to the fact that the mean size from an intensity-weighted size distribution is larger than that from its corresponding number-weighted size distribution. Although the low size detection limit of NTA was pointed out as a source of the contradiction by reducing the variance measured by NTA, it does not fully explain why the mean sizes of V50 were less different between DLS and NTA compared with V100 and V200 despite its higher PI value.

Additionally, the number-weighted size distribution of the vesicle samples from NTA was converted into an intensity-weighted distribution by assuming the thin-shell hollow sphere model with the RGD approximation to verify the influence of the different quantities that DLS and NTA produce. While the mean size given by NTA after conversion was very close to the value before conversion, the RV after the conversion was much reduced compared to that from DLS and NTA. Considering the nature of the vesicle samples, it is questionable whether the small relative variance after reconstruction is reliable. Therefore, the conversion of a number-weighted size distribution of NTA into an intensity-weighted one does not seem to effectively explain the difference in the mean sizes measured by DLS and NTA, and the size information from DLS and NTA is better compared in their original weightings.

Theory

Nanoparticle Tracking Analysis

Similar to DLS, NTA extracts the size information of particles in suspension by measuring their diffusion coefficient51. As illustrated in Fig. 6, by taking sequential images of illuminated particles in suspension on a periodic time interval, the displacement of a particle can be identified from successive images and constructed into a track. To determine the displacement of particles, NTA compares the two-dimensional location of particles in an image frame with the subsequent frame. In doing so, NTA sets a certain threshold distance, also known as the maximum jump distance, to properly identify if the two particles in the two adjacent frames are the same particle. If any single particle is found in the successive frame within the threshold distance from the location of the particle in the previous image, the two particles are recognized as the same particle and make a track segment. In the same manner, this recognition process is performed for subsequent frames and track segments are combined to construct a particle’s track. The track terminates if there is no particle in the following frame or if there are more than two particles within the threshold distance. Because of the nature of the tracking process, the track segment length, and the number of track segments, are finite and variable.

Figure 6

(a) Instrumental schematics of NTA detection. The particles suspended in liquid are illuminated by laser beam. Scattered light by the particles is monitored by the microscope, from which the position of the particles is identified by the software. (b) The software tracks the trajectory of the particles and constructs particle tracks. The average squared displacement z of each track is then converted into the corresponding hydrodynamic diameter d according to the Stokes-Einstein equation.

Full size image

For a track of a track segment length n, the mean squared displacement of the track z, expressed as

$${z}=\frac{1}{n}\sum _{i=1}^{n}\,{r}_{i}^{2}$$

(1)

where ri is the two-dimensional displacement of ith track segment of the track, is translated into the diffusion coefficient assuming a 2D Brownian motion, which is related by

$${z}=4{D}{\rm{\Delta }}{t}$$

(2)

where D is the diffusion coefficient of the tracked particle and Δt is the time interval of the image frames51,52,57,58. Then the acquired diffusion coefficient D of the track is converted to the hydrodynamic diameter relying on the Stokes-Einstein equation,

$$D=\frac{{k}_{B}T}{3\pi \eta d}$$

(3)

where kB is the Boltzmann constant, T is the temperature, η is the viscosity of the medium and d is the hydrodynamic diameter of the tracked particle, and comprises the size distribution of the measured sample.

Size Distribution Estimation Methods for NTA

Uncertainty in the Mean Squared Displacement Measurement in NTA

Measurement of the mean squared displacement z in NTA assumes that a particle is tracked for long enough that the measured mean squared displacement is close enough to the ideal mean squared displacement. However, as the track segment length is finite, the acquired mean squared displacement z of a track has statistical uncertainty that is inversely proportional to the square root of the track segment length n54,63. Therefore, small track segment lengths would make the measured size distribution significantly broader than its true size distribution and they should be excluded to have a narrower size distribution.

Maximum Likelihood Estimation (MLE) with an Assumed Arbitrary Distribution with Parameters

As particles observed by NTA undergo Brownian motion, the probability distribution of the mean squared displacement z of a tracked particle whose size is d with a track segment length n can be described by a gamma distribution as follows58:

$$P(z|n,d)=\frac{{(ns)}^{n}}{(n-1)!}{z}^{n-1}\,{e}^{-nsz}$$

(4)

where s is defined as \(1/4D{\rm{\Delta }}t=3\pi \eta d/4{\rm{\Delta }}t{k}_{B}T\). An observed set of tracks Φ of a sample measured by NTA can be related to a particle size distribution f(d) by the integration

$${\rm{\Phi }}={\int }_{0}^{\infty }\,f(x)P(z|n,x)\,dx$$

(5)

Although the inversion of the integration would enable recovery of f(d) by finding Q(d|n, z) that satisfies

$$f(d)=\sum _{{\rm{\Phi }}}\,{\rm{\Phi }}\cdot Q(d|n,z)$$

(6)

the inversion relation cannot be used because Q(d|n, z) depends on f(d)60.

Due to the difficulty of finding Q for the inversion, finite track length adjustment (FTLA) can be used to determine f(d) by maximizing the likelihood of the observed list of tracks Φ with an assumed size distribution for f(d) as illustrated in Fig. 7 57. For an assumed size distribution with a few parameters, for example a log-normal distribution with its geometric mean and geometric standard deviation set as parameters, an “ideal” set of particles di (with j = 1–999) is drawn, which comprises 1000-quantiles of the size distribution, i.e., those of particle size di that satisfy C(dj) = j/1000 (j = 1–999) for the cumulative size distribution C(d) of f(d). For each particle size dj, the likelihood of producing a track whose mean squared displacement z and track segment length n is given by

$$L(z,n;{d}_{j})=P(z|n,{d}_{j})$$

(7)

Figure 7

Schematic representation of the size determination methods based on maximum likelihood estimation. By assuming a size distribution described by adjustable parameters, the FTLA method looks for the parameter values that maximize the likelihood of the given NTA tracks. The iterative correction method refines the size distribution that approaches to the maximum likelihood over each iteration.

Full size image

Hence, the likelihood of producing such a track with a given particle set dj (j = 1–999) is given as

$$L(z,n;\{d\})=\sum _{j=1}^{999}P(z|n,{d}_{j})$$

(8)

Likewise, the likelihood for all the observed list of tracks Φ becomes

$$L({\rm{\Phi }};\{d\})=\prod _{k=1}^{K}(\sum _{j=1}^{999}P({z}_{k}|{n}_{k},{d}_{j}))$$

(9)

where zk and nk are the mean squared displacement and the track segment length of the k th track of Φ, respectively, and K is the total number of the observed tracks. For practical reasons, the maximization is performed on the logarithm of the likelihood, as maximizing a positive function also maximizes its logarithm58:

$$LL({\rm{\Phi }};\{d\})=\,\mathrm{log}(L)=\sum _{k=1}^{K}\,\mathrm{log}(\sum _{j=1}^{999}\,P({z}_{k}|{n}_{k},{d}_{j}))$$

(10)

With the log-likelihood, an optimal parameter set for the size distribution is sought for its maximum value, which determines the size distribution for the observed tracks. As this approach assumes an arbitrary size distribution, we can expect a decrease in the likelihood for the benefit of a simpler solution60.

Maximum Likelihood Estimation by Iterative Correction

For a size distribution f(d), the log-likelihood of producing the observed list of tracks Φ can be expressed as

$$LL({\rm{\Phi }})=\sum _{k=1}^{K}\,\mathrm{log}(\frac{1}{{\sum }_{b=1}^{B}\,f({d}_{b})}\sum _{b=1}^{B}\,P({z}_{k}|{n}_{k},{d}_{b})f({d}_{b}))$$

(11)

where b is the bin number for the segmented range of particle size d58.

Then its differentiation with respect to the size distribution is given as

$$\frac{\partial LL}{\partial f}=\sum _{k=1}^{K}\frac{P({z}_{k}|{n}_{k},{d}_{b})}{{\sum }_{b=1}^{B}\,P({z}_{k}|{n}_{k},{d}_{b})f({d}_{b})}-\frac{K}{{\sum }_{b=1}^{B}\,f({d}_{b})}$$

(12)

At the maximum of the likelihood, the differential becomes zero, leading to

$$\sum _{k=1}^{K}\frac{P({z}_{k}|{n}_{k},{d}_{b})}{{\sum }_{b=1}^{B}\,P({z}_{k}|{n}_{k},{d}_{b})f({d}_{b})}=\frac{K}{{\sum }_{b=1}^{B}\,f({d}_{b})}$$

(13)

From this relation, it is possible to find the size distribution by an iterative procedure, given as

$${f}^{(r+1)}({d}_{b})={f}^{(r)}({d}_{b})\cdot \frac{1}{K}\sum _{k=1}^{K}\frac{P({z}_{k}|{n}_{k},{d}_{b})}{{\sum }_{b=1}^{B}\,P({z}_{k}|{n}_{k},{d}_{b})f({d}_{b})/{\sum }_{b=1}^{B}\,f({d}_{b})}$$

(14)

where f (r)(d) is the r th estimate of the size distribution from the iterations as illustrated in Fig. 7 60. For the initial size distribution f (0)(d), a uniform distribution is chosen58,64,65.

For the termination criteria of the iteration, the chi-squared statistic of the error between a histogram of the mean squared displacement Hz and that calculated from the rth iteration solution H(r) can be used, where

$${\chi }^{2}=\sum _{m=1}^{M}\frac{{({H}_{z}({z}_{m})-{H}^{(r)}({z}_{m}))}^{2}}{{H}^{(r)}({z}_{m})}$$

(15)

$${H}^{(r)}({z}_{m})=\sum _{n={N}_{{\min }}}^{{N}_{\max }}{N}_{n}\sum _{b=1}^{B}\frac{P({z}_{m}|n,{d}_{b}){f}^{(r)}({d}_{b}){\rm{\Delta }}{z}_{m}}{{\sum }_{b=1}^{B}\,{f}^{(r)}({d}_{b})}$$

(16)

where m is the bin number for the mean squared displacement and Nn is the number of tracks with track segment length n, and in this study the iteration terminates when the change in χ2 becomes smaller than 1% of the previous value58.

Conversion of Number-Weighted Distribution of NTA

The size distribution measured by DLS is weighted by the intensity of the scattered light, which is dependent on the particle size, whereas that from NTA is weighted by the number. In comparing the two quantities, one must be converted to match the weighting of the other.

The number-weighted size distribution f(d) from NTA can be converted into an intensity-weighted distribution fI(d) given by

$${f}_{I}(d)=\frac{f(d)I(q,d)}{{\int }_{0}^{\infty }f(d)I(q,d)dd}$$

(17)

where q is the scattering vector and I(q,d) is the form factor of the measured particles66,67. For vesicles, we assume a thin-shell hollow sphere model48,68,69, so that

$$I(q,d)={(\pi {d}^{2}t)}^{2}P(q,d)$$

(18)

where P(q,d) is the structural factor for the vesicle approximated by the Rayleigh-Gans-Debye (RGD) approximation given by

$$P(q,d)={(\frac{\sin (q\cdot d/2)}{q\cdot d/2})}^{2}$$

(19)

Materials and Methods

Reagents

Polystyrene latex standards (92 ± 3, 269 ± 7 and 343 ± 9 nm) were purchased from Thermo Scientific (Rockford, IL., USA). 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was purchased from Avanti Polar Lipids Inc., (Alabaster, AL, USA).

Vesicle preparation

Vesicles composed of egg phosphatidylcholine (PC) (Avanti Polar Lipids Inc., Alabaster, AL) were prepared at a lipid concentration of ∼5 mg/mL and then diluted before the experiment. Briefly, dried lipid films were hydrated with 10 mM Tris (pH 7.5) buffer solution with 150 mM NaCl, and the sample was then vortexed periodically for 5 min. Vesicle samples were extruded through polycarbonate membranes with either 50, 100, 200 or 400 nm pores sized by a miniextruder (Avanti Polar Lipids). For those samples pretreated with a range of freeze-thaw cycles, freeze-thaw treatment was performed before extrusion on newly hydrated lipid films by using a previously described methodology based on freeze-thaw pretreatment and then extrusion61. Specifically, in each treatment cycle, the vesicle suspension was frozen in liquid nitrogen for 30 s, before thawing in an 80 °C water bath for 90 s, and then finally being vortexed to complete each cycle. After the freeze-thaw cycles of 3, 5, 7, 9, 11, 13, 15 or 17 repetitions, vesicles were sized by an extruder (Avanti Polar Lipids) through 400 nm polycarbonate membrane pores. All aqueous solutions and buffers were prepared in Milli-Q water with a minimum resistivity of 18.2 MΩ·cm (Millipore, Billerica, MA).

Experimental setup

For the DLS measurements, a ZetaPals particle size analyzer (Brookhaven Instruments, Holtsville, NY) with a 658.0 nm monochromatic laser was used. For each sample, three independent runs of 1 min were performed. To avoid unnecessary reflection, all measurements were taken at a scattering angle of 90°, and the measured intensity autocorrelation function was fitted to yield the intensity-weighted size distribution of particles in solution. The deconvolution of the autocorrelation function was done using the cumulants method, which was applied to calculate the intensity-weighted log-normal profile of the size distribution expressed by the average effective diameter and its polydispersity.

NTA measurements were made with an LM10HS (Nanosight Limited, Amesbury, UK) equipped with a scientific CMOS camera, a 20x objective lens, a blue laser module (405 nm, LM12 version C) and NTA software version 3.1. A 1-ml disposable syringe was used to inject the samples into the instrument chamber. The video data for NTA measurements were collected for 30 seconds, repeated three times for each sample. The detection threshold of the NTA software was set to 5 and the maximum jump distance and the minimum track segment length were both set to auto.

Detected tracks were then translated into a size distribution using three different methods, i.e., direct conversion of the detected particle tracks, maximum likelihood estimation with an assumed distribution (the FTLA method) and maximum likelihood estimation by iterative correction. For the conversion, valid tracks were acquired from the detected tracks by the software, with tracks of a small track segment length excluded from the selection.

FTLA assumes a certain shape for the size distribution to maximize likelihood, and we assumed a log-normal distribution in this study so that the mean and standard deviation parameters could be optimized to produce maximum likelihood57. For the calculation, 1000-quantiles are generated from a log-normal distribution while varying its mean size and standard deviation, which represent an “ideal” set of particles for the assumed size distribution57.

In the iterative method, the size distribution f(d) is refined by the iteration

$$\begin{array}{rcl}{f}^{(j+1)}({d}_{i}) & = & {f}^{(j)}({d}_{i})\cdot \frac{1}{N}\sum _{k}\,\frac{P({z}_{k};{n}_{k},{d}_{i})}{{Q}^{(j)}}\\ {Q}^{(j)} & = & \sum _{i}\,P({z}_{k};{n}_{k},{d}_{i}){f}^{(j)}({d}_{i})\end{array}$$

(20)

where f (j)(di) is a normalized fraction of those particles in the size distribution of the jth iteration whose size is between di and di+1, and N is the number of tracks in the given set of tracks. The initial size distribution f (0)(d) is set to a uniform distribution so that f (0)(d) = 1/M, where M is the number of bins for the size distribution. The iterations are performed until the change in the chi-squared value between the histogram of the displacement of the observed tracks and that of the estimated size distribution from the iteration is less than 1% of the previous value58.

In comparing the results from DLS and NTA, the mean size of NTA from the acquired size distribution is compared with the mean hydrodynamic size Rz of DLS while the relative variance (RV, the variance divided by the square of the mean size) of NTA is compared with the PI of DLS.

Data Availability

The data that support the findings of this study are present in the paper and the Supplementary Information. Additional information is available from the authors upon reasonable request.

References

1. Shi, J., Votruba, A. R., Farokhzad, O. C. & Langer, R. Nanotechnology in drug delivery and tissue engineering: from discovery to applications. Nano Lett 10, 3223–3230, https://doi.org/10.1021/nl102184c (2010).

   Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

2. Parveen, S., Misra, R. & Sahoo, S. K. Nanoparticles: a boon to drug delivery, therapeutics, diagnostics and imaging. Nanomedicine 8, 147–166, https://doi.org/10.1016/j.nano.2011.05.016 (2012).

   Article  CAS  PubMed  Google Scholar 

3. Blanco, E., Shen, H. & Ferrari, M. Principles of nanoparticle design for overcoming biological barriers to drug delivery. Nat Biotechnol 33, 941–951, https://doi.org/10.1038/nbt.3330 (2015).

   Article  CAS  PubMed  PubMed Central  Google Scholar 

4. Chang, H. I. & Yeh, M. K. Clinical development of liposome-based drugs: formulation, characterization, and therapeutic efficacy. Int J Nanomedicine 7, 49–60, https://doi.org/10.2147/IJN.S26766 (2012).

   Article  CAS  PubMed  Google Scholar 

5. Tan, M. L., Choong, P. F. & Dass, C. R. Recent developments in liposomes, microparticles and nanoparticles for protein and peptide drug delivery. Peptides 31, 184–193, https://doi.org/10.1016/j.peptides.2009.10.002 (2010).

   Article  CAS  PubMed  Google Scholar 

6. Bamrungsap, S. et al. Nanotechnology in therapeutics: a focus on nanoparticles as a drug delivery system. Nanomedicine (Lond) 7, 1253–1271, https://doi.org/10.2217/nnm.12.87 (2012).

   Article  CAS  Google Scholar 

7. Warne, N. W. Development of high concentration protein biopharmaceuticals: the use of platform approaches in formulation development. Eur J Pharm Biopharm 78, 208–212, https://doi.org/10.1016/j.ejpb.2011.03.004 (2011).

   Article  CAS  PubMed  Google Scholar 

8. Tsuchikama, K. & An, Z. Antibody-drug conjugates: recent advances in conjugation and linker chemistries. Protein Cell 9, 33–46, https://doi.org/10.1007/s13238-016-0323-0 (2018).

   Article  CAS  PubMed  Google Scholar 

9. Li, B. L. et al. Directing Assembly and Disassembly of 2D MoS2 Nanosheets with DNA for Drug Delivery. ACS Appl Mater Interfaces 9, 15286–15296, https://doi.org/10.1021/acsami.7b02529 (2017).

   Article  CAS  PubMed  Google Scholar 

10. Komiyama, M., Yoshimoto, K., Sisido, M. & Ariga, K. Chemistry Can Make Strict and Fuzzy Controls for Bio-Systems: DNA Nanoarchitectonics and Cell-Macromolecular Nanoarchitectonics. Bull Chem Soc Jpn 90, 967–1004, https://doi.org/10.1246/bcsj.20170156 (2017).

    Article  Google Scholar 

11. Muller, R. H. & Keck, C. M. Challenges and solutions for the delivery of biotech drugs–a review of drug nanocrystal technology and lipid nanoparticles. J Biotechnol 113, 151–170, https://doi.org/10.1016/j.jbiotec.2004.06.007 (2004).

    Article  CAS  PubMed  Google Scholar 

12. Brigger, I., Dubernet, C. & Couvreur, P. Nanoparticles in cancer therapy and diagnosis. Adv Drug Deliv Rev 64, 24–36, https://doi.org/10.1016/j.addr.2012.09.006 (2012).

    Article  Google Scholar 

13. Barenholz, Y. Doxil®, the first FDA-approved nano-drug: Lessons learned. J Control Release 160, 117–134, https://doi.org/10.1016/j.jconrel.2012.03.020 (2012).

    Article  CAS  PubMed  Google Scholar 

14. Cavalli, A. et al. Multi-target-directed ligands to combat neurodegenerative diseases. J Med Chem 51, 347–372, https://doi.org/10.1021/jm7009364 (2008).

    Article  CAS  PubMed  Google Scholar 

15. Sams-Dodd, F. Target-based drug discovery: is something wrong? Drug Discov Today 10, 139–147, https://doi.org/10.1016/s1359-6446(04)03316-1 (2005).

    Article  CAS  PubMed  Google Scholar 

16. Gregoriadis, G. Engineering liposomes for drug delivery: Progress and problems. Trends Biotechnol 13, 527–537, https://doi.org/10.1016/s0167-7799(00)89017-4 (1995).

    Article  CAS  PubMed  Google Scholar 

17. Rawat, M., Singh, D., Saraf, S. & Saraf, S. Nanocarriers: promising vehicle for bioactive drugs. Biol Pharm Bull 29, 1790–1798 (2006).

    Article  CAS  Google Scholar 

18. Al-Jamal, W. T. & Kostarelos, K. Liposomes: from a clinically established drug delivery system to a nanoparticle platform for theranostic nanomedicine. Acc Chem Res 44, 1094–1104, https://doi.org/10.1021/ar200105p (2011).

    Article  CAS  PubMed  Google Scholar 

19. Metselaar, J. M. & Storm, G. Liposomes in the treatment of inflammatory disorders. Expert Opin Drug Deliv 2, 465–476, https://doi.org/10.1517/17425247.2.3.465 (2005).

    Article  CAS  PubMed  Google Scholar 

20. He, C. B., Yin, L. C., Tang, C. & Yin, C. H. Size-dependent absorption mechanism of polymeric nanoparticles for oral delivery of protein drugs. Biomaterials 33, 8569–8578, https://doi.org/10.1016/j.biomaterials.2012.07.063 (2012).

    Article  CAS  PubMed  Google Scholar 

21. Banerjee, A., Qi, J. P., Gogoi, R., Wong, J. & Mitragotri, S. Role of nanoparticle size, shape and surface chemistry in oral drug delivery. J Control Release 238, 176–185, https://doi.org/10.1016/j.jconrel.2016.07.051 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

22. Coradeghini, R. et al. Size-dependent toxicity and cell interaction mechanisms of gold nanoparticles on mouse fibroblasts. Toxicol Lett 217, 205–216, https://doi.org/10.1016/j.toxlet.2012.11.022 (2013).

    Article  CAS  PubMed  Google Scholar 

23. Braun, N. J., DeBrosse, M. C., Hussain, S. M. & Comfort, K. K. Modification of the protein corona-nanoparticle complex by physiological factors. Mater Sci Eng C Mater Biol Appl 64, 34–42, https://doi.org/10.1016/j.msec.2016.03.059 (2016).

    Article  CAS  PubMed  Google Scholar 

24. Hama, S. et al. Overcoming the polyethylene glycol dilemma via pathological environment-sensitive change of the surface property of nanoparticles for cellular entry. J Control Release 206, 67–74, https://doi.org/10.1016/j.jconrel.2015.03.011 (2015).

    Article  CAS  PubMed  Google Scholar 

25. Natarajan, S. K. & Selvaraj, S. Mesoporous silica nanoparticles: importance of surface modifications and its role in drug delivery. RSC Adv. 4, 14328–14334, https://doi.org/10.1039/c4ra00781f (2014).

    Article  CAS  Google Scholar 

26. Soliman, G. M. et al. Tailoring the efficacy of nimodipine drug delivery using nanocarriers based on A2B miktoarm star polymers. Biomaterials 31, 8382–8392, https://doi.org/10.1016/j.biomaterials.2010.07.039 (2010).

    Article  CAS  PubMed  Google Scholar 

27. Coelho, J. F. et al. Drug delivery systems: Advanced technologies potentially applicable in personalized treatments. EPMA J 1, 164–209, https://doi.org/10.1007/s13167-010-0001-x (2010).

    Article  PubMed  PubMed Central  Google Scholar 

28. Love, W. G., Amos, N., Kellaway, I. W. & Williams, B. D. Specific accumulation of cholesterol-rich liposomes in the inflammatory tissue of rats with adjuvant arthritis. Ann Rheum Dis 49, 611–614 (1990).

    Article  CAS  Google Scholar 

29. Hua, S. & Wu, S. Y. The use of lipid-based nanocarriers for targeted pain therapies. Front Pharmacol 4, 7, https://doi.org/10.3389/fphar.2013.00143 (2013).

    Article  CAS  Google Scholar 

30. de Temmerman, P. J., Verleysen, E., Lammertyn, J. & Mast, J. Size measurement uncertainties of near-monodisperse, near-spherical nanoparticles using transmission electron microscopy and particle-tracking analysis. J Nanopart Res 16, 17, https://doi.org/10.1007/s11051-014-2628-3 (2014).

    Article  Google Scholar 

31. Troiber, C. et al. Comparison of four different particle sizing methods for siRNA polyplex characterization. Eur J Pharm Biopharm 84, 255–264, https://doi.org/10.1016/j.ejpb.2012.08.014 (2013).

    Article  CAS  PubMed  Google Scholar 

32. Fraunhofer, W., Winter, G. & Coester, C. Asymmetrical flow field-flow fractionation and multiangle light scattering for analysis of gelatin nanoparticle drug carrier systems. Anal Chem 76, 1909–1920, https://doi.org/10.1021/ac0353031 (2004).

    Article  CAS  PubMed  Google Scholar 

33. Liu, J., Andya, J. D. & Shire, S. J. A critical review of analytical ultracentrifugation and field flow fractionation methods for measuring protein aggregation. AAPS J 8, E580–E589, https://doi.org/10.1208/aapsj080367 (2006).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

34. Berne, B. J. & Pecora, R. Dynamic light scattering: with applications to chemistry, biology, and physics. Dover edn, (Dover Publications, 2000).

35. Maulucci, G. et al. Particle size distribution in DMPC vesicles solutions undergoing different sonication times. Biophys J 88, 3545–3550, https://doi.org/10.1529/biophysj.104.048876 (2005).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

36. Merkus, H. G. Particle size measurements: fundamentals, practice, quality. (Springer, 2009).

37. Sokolova, V. et al. Characterisation of exosomes derived from human cells by nanoparticle tracking analysis and scanning electron microscopy. Colloids Surf B Biointerfaces 87, 146–150, https://doi.org/10.1016/j.colsurfb.2011.05.013 (2011).

    Article  CAS  PubMed  Google Scholar 

38. Anderson, W., Kozak, D., Coleman, V. A., Jamting, A. K. & Trau, M. A comparative study of submicron particle sizing platforms: accuracy, precision and resolution analysis of polydisperse particle size distributions. J Colloid Interface Sci 405, 322–330, https://doi.org/10.1016/j.jcis.2013.02.030 (2013).

    Article  ADS  CAS  PubMed  Google Scholar 

39. Bootz, A., Vogel, V., Schubert, D. & Kreuter, J. Comparison of scanning electron microscopy, dynamic light scattering and analytical ultracentrifugation for the sizing of poly(butyl cyanoacrylate) nanoparticles. Eur J Pharm Biopharm 57, 369–375, https://doi.org/10.1016/S0939-6411(03)00193-0 (2004).

    Article  CAS  PubMed  Google Scholar 

40. Striegel, A. M. Modern size-exclusion liquid chromatography: practice of gel permeation and gel filtration chromatography. 2nd edn, (Hoboken, N. J.: Wiley, 2009).

41. Yohannes, G., Jussila, M., Hartonen, K. & Riekkola, M. L. Asymmetrical flow field-flow fractionation technique for separation and characterization of biopolymers and bioparticles. J Chromatogr A 1218, 4104–4116, https://doi.org/10.1016/j.chroma.2010.12.110 (2011).

    Article  CAS  PubMed  Google Scholar 

42. Wyatt, P. J. Light scattering and the absolute characterization of macromolecules. Anal Chim Acta 272, 1–40, https://doi.org/10.1016/0003-2670(93)80373-S (1993).

    Article  CAS  Google Scholar 

43. Podzimek, S. Light scattering, size exclusion chromatography, and asymmetric flow field flow fractionation: powerful tools for the characterization of polymers, proteins, and nanoparticles. (Hoboken, N. J.: Wiley, 2011).

44. van der Pol, E., Coumans, F., Varga, Z., Krumrey, M. & Nieuwland, R. Innovation in detection of microparticles and exosomes. J Thromb Haemost 11(Suppl 1), 36–45, https://doi.org/10.1111/jth.12254 (2013).

    Article  PubMed  Google Scholar 

45. Hassan, P. A., Rana, S. & Verma, G. Making sense of Brownian motion: colloid characterization by dynamic light scattering. Langmuir 31, 3–12, https://doi.org/10.1021/la501789z (2015).

    Article  CAS  PubMed  Google Scholar 

46. Stetefeld, J., McKenna, S. A. & Patel, T. R. Dynamic light scattering: a practical guide and applications in biomedical sciences. Biophys Rev 8, 409–427, https://doi.org/10.1007/s12551-016-0218-6 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

47. Boyd, R. D., Pichaimuthu, S. K. & Cuenat, A. New approach to inter-technique comparisons for nanoparticle size measurements: using atomic force microscopy, nanoparticle tracking analysis and dynamic light scattering. Colloids Surf A Physicochem Eng Asp 387, 35–42, https://doi.org/10.1016/j.colsurfa.2011.07.020 (2011).

    Article  CAS  Google Scholar 

48. Hallett, F. R., Watton, J. & Krygsman, P. Vesicle sizing: Number distributions by dynamic light scattering. Biophys J 59, 357–362, https://doi.org/10.1016/S0006-3495(91)82229-9 (1991).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

49. Bell, N. C., Minelli, C., Tompkins, J., Stevens, M. M. & Shard, A. G. Emerging techniques for submicrometer particle sizing applied to Stober silica. Langmuir 28, 10860–10872, https://doi.org/10.1021/la301351k (2012).

    Article  CAS  PubMed  Google Scholar 

50. Amini, R., Brar, S. K., Cledon, M. & Surampalli, R. Y. Intertechnique comparisons for nanoparticle size measurements and shape distribution. J Hazard Toxic Radioact Waste 20, B4015004, https://doi.org/10.1061/(asce)hz.2153-5515.0000286 (2016).

    Article  CAS  Google Scholar 

51. Malloy, A. & Carr, B. NanoParticle tracking analysis - The Halo™ system. Part Part Syst Charact 23, 197–204, https://doi.org/10.1002/ppsc.200601031 (2006).

    Article  Google Scholar 

52. Filipe, V., Hawe, A. & Jiskoot, W. Critical evaluation of nanoparticle tracking analysis (NTA) by NanoSight for the measurement of nanoparticles and protein aggregates. Pharm Res 27, 796–810, https://doi.org/10.1007/s11095-010-0073-2 (2010).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

53. James, A. E. & Driskell, J. D. Monitoring gold nanoparticle conjugation and analysis of biomolecular binding with nanoparticle tracking analysis (NTA) and dynamic light scattering (DLS). Analyst 138, 1212–1218, https://doi.org/10.1039/c2an36467k (2013).

    Article  ADS  CAS  PubMed  Google Scholar 

54. Dragovic, R. A. et al. Sizing and phenotyping of cellular vesicles using nanoparticle tracking analysis. Nanomedicine 7, 780–788, https://doi.org/10.1016/j.nano.2011.04.003 (2011).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

55. Yang, D. T., Lu, X., Fan, Y. & Murphy, R. M. Evaluation of Nanoparticle Tracking for Characterization of Fibrillar Protein Aggregates. AIChE J 60, 1236–1244, https://doi.org/10.1002/aic.14349 (2014).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

56. Gross, J., Sayle, S., Karow, A. R., Bakowsky, U. & Garidel, P. Nanoparticle tracking analysis of particle size and concentration detection in suspensions of polymer and protein samples: Influence of experimental and data evaluation parameters. Eur J Pharm Biopharm 104, 30–41, https://doi.org/10.1016/j.ejpb.2016.04.013 (2016).

    Article  CAS  PubMed  Google Scholar 

57. Saveyn, H. et al. Accurate particle size distribution determination by nanoparticle tracking analysis based on 2-D Brownian dynamics simulation. J Colloid Interface Sci 352, 593–600, https://doi.org/10.1016/j.jcis.2010.09.006 (2010).

    Article  ADS  CAS  PubMed  Google Scholar 

58. Walker, J. G. Improved nano-particle tracking analysis. Meas Sci Technol 23, 065605, https://doi.org/10.1088/0957-0233/23/6/065605 (2012).

    Article  ADS  CAS  Google Scholar 

59. Kestens, V., Bozatzidis, V., De Temmerman, P. J., Ramaye, Y. & Roebben, G. Validation of a particle tracking analysis method for the size determination of nano- and microparticles. J Nanopart Res 19, 271, https://doi.org/10.1007/s11051-017-3966-8 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

60. Lucy, L. B. Iterative technique for rectification of observed distributions. Astronom J 79, 745–754, https://doi.org/10.1086/111605 (1974).

    Article  ADS  Google Scholar 

61. Mayer, L. D., Hope, M. J., Cullis, P. R. & Janoff, A. S. Solute distributions and trapping efficiencies observed in freeze-thawed multilamellar vesicles. Biochim Biophys Acta 817, 193–196 (1985).

    Article  CAS  Google Scholar 

62. Matsuzaki, K. et al. Optical characterization of liposomes by right angle light scattering and turbidity measurement. Biochim Biophys Acta 1467, 219–226 (2000).

    Article  CAS  Google Scholar 

63. Gardiner, C., Ferreira, Y. J., Dragovic, R. A., Redman, C. W. & Sargent, I. L. Extracellular vesicle sizing and enumeration by nanoparticle tracking analysis. J Extracell Vesicles 2, https://doi.org/10.3402/jev.v2i0.19671 (2013).

64. Veklerov, E. & Llacer, J. Stopping rule for the MLE algorithm based on statistical hypothesis testing. IEEE Trans Med Imag 6, 313–319, https://doi.org/10.1109/tmi.1987.4307849 (1987).

    Article  CAS  Google Scholar 

65. Hebert, T. J. Statistical stopping criteria for iterative maximum-likelihood reconstruction of emission images. Phys Med Biol 35, 1221–1232 (1990).

    Article  Google Scholar 

66. Hanus, L. H. & Ploehn, H. J. Conversion of intensity-averaged photon correlation spectroscopy measurements to number-averaged particle size distributions. 1. Theoretical development. Langmuir 15, 3091–3100, https://doi.org/10.1021/la980958w (1999).

    Article  CAS  Google Scholar 

67. Patty, P. J. & Frisken, B. J. Direct determination of the number-weighted mean radius and polydispersity from dynamic light-scattering data. Appl Opt 45, 2209–2216 (2006).

    Article  ADS  Google Scholar 

68. Pecora, R. & Aragon, S. R. Theory of light scattering from hollow spheres. Chem Phys Lipids 13, 1–10 (1974).

    Article  CAS  Google Scholar 

69. Pencer, J. & Hallett, F. R. Effects of vesicle size and shape on static and dynamic light scattering measurements. Langmuir 19, 7488–7497, https://doi.org/10.1021/la0345439 (2003).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Research Foundation of Singapore through a Proof-of-Concept grant (NRF2015NRF-POC0001-19) and as well as through the Center for Precision Biology at Nanyang Technological University.

Author information

Authors and Affiliations

1. School of Materials Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, 639798, Singapore, Singapore

   Ahram Kim, Wei Beng Ng & Nam-Joon Cho

2. Centre for Biomimetic Sensor Science, Nanyang Technological University, 50 Nanyang Drive, 637553, Singapore, Singapore

   Wei Beng Ng & Nam-Joon Cho

3. Particle Characterization Laboratories, Inc. 845 Olive Ave, Suite A, Novato, CA, 94945, USA

   William Bernt

4. School of Chemical and Biomedical Engineering, Nanyang Technological University, 62 Nanyang Drive, 637459, Singapore, Singapore

   Nam-Joon Cho

Authors

1. Ahram Kim

   You can also search for this author in PubMed Google Scholar

2. Wei Beng Ng

   You can also search for this author in PubMed Google Scholar

3. William Bernt

   You can also search for this author in PubMed Google Scholar

4. Nam-Joon Cho

   You can also search for this author in PubMed Google Scholar

Contributions

A.K. and N.J.C. made an initial conceptual idea and wrote the manuscript. A.K. analyzed the results and prepared the figures. A.K., W.B.N. and W.B. carried out the measurements. All the authors reviewed the manuscript.

Corresponding author

Correspondence to Nam-Joon Cho.

Ethics declarations

Competing Interests

The authors declare no competing interests.

Additional information

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Cite this article

Kim, A., Ng, W.B., Bernt, W. et al. Validation of Size Estimation of Nanoparticle Tracking Analysis on Polydisperse Macromolecule Assembly. Sci Rep 9, 2639 (2019). https://doi.org/10.1038/s41598-019-38915-x

Download citation

• Received: 29 October 2018

• Accepted: 11 January 2019

• Published: 25 February 2019

• DOI: https://doi.org/10.1038/s41598-019-38915-x

What measure is used to characterize the peak tall skinny nature of the distribution?

What are the three layers of information required in a strong performance dashboard design?

Are the mathematical methods used to estimate or describe the degree of variation in a given variable of interest?

What information can you find in an analyze dashboard quizlet?