R.
Maugi
a,
P.
Hauer
b,
J.
Bowen
c,
E.
Ashman
d,
E.
Hunsicker
*e and
M.
Platt
*a
aDepartment of Chemistry, Loughborough University, Loughborough, Leicestershire LE11 3TU, UK. E-mail: m.platt@lboro.ac.uk
bSchool of Chemical and Physical Sciences, Victoria University of Wellington, Wellington 6140, New Zealand
cSchool of Pharmacy and Pharmaceutical Sciences, Cardiff University, Cardiff, CF10 3NB, UK
dCMD Ltd, Cardiff, CF15 7AB, UK
eDepartment of Mathematical Sciences, Centre for Imaging Science, Loughborough University, Loughborough, Leicestershire LE11 3TU, UK. E-mail: e.hunsicker@lboro.ac.uk
First published on 26th November 2019
The discovery and characterisation of nanomaterials represents a multidisciplinary problem. Their properties and applications within biological, physical and medicinal sciences depend on their size, shape, concentration and surface charge. No single technology can currently measure all characteristics. Here we combine resistive pulse sensing with predictive logistic regression models, termed RPS-LRM, to rapidly characterise a nanomaterial's size, aspect ratio, shape and concentration when mixtures of nanorods and nanospheres are present in the same solution. We demonstrate that RPS-LRM can be applied to the characterisation of nanoparticles over a wide size range, and varying aspect ratios, and can distinguish between nanorods over nanospheres when they possess an aspect ratio grater then two. The RPS-LRM can rapidly measure the ratios of nanospheres to nanorods in solution within mixtures, regardless of their relative sizes and ratios i.e. many large nanospherical particles do not interfere with the characterisation of smaller nanorods. This was done with a 91% correct classification of nanospherical particles and 72% correct classification of nanorods even when the fraction of nanorods in solution is as low as 20%. The methodology here will enable the classification of nanomedicines, new nanomaterials and biological analytes in solution.
Over the last two decades, Resistive Pulse Sensing (RPS) has emerged as a promising technology for the analysis of individual nanoparticles. The translocation of a nanoparticle or analyte through a nanopore can be monitored by measuring ionic current.3–7 Each translocation event causes a change to the conductance of the channel known as a Resistive Pulse, Fig. 1a, which is related to the physical properties of the analyte.8,9 RPS is a highly attractive sensing platform with single particle/analyte resolution, providing information on the analyte size, concentration and charge.5,10 They perform the analysis quickly, reliably and to a high level of sensitivity.11–14 Indeed, RPS has exhibited high resolution for the monitoring of various analytes, from metal ions to whole cells15–19
Recent advances in RPS methodologies have focused upon improving the analysis, fabrication and signal to noise ratio,3,20–22 and have started to advance RPS as a method to identify particle shape.23–28 In these, the dimensions of the nanopores used have been similar those of the particles studied. Although these results are promising, the requirement for pore sizes to be similar in size as the particles being analysed has drawbacks. First, requiring knowledge about expected particle sizes limits the technique for the study of unknown or complex mixtures of particles. Second, the low pore-to-particle size ratios are prone to blockages. By contrast, larger pores are less likely to block and have the added benefit of a higher flow rate of liquid through the pore, which reduces the time required to scan particles.29–31 In previous studies, we have used RPS sensors to study aggregates of nanorods and nanospheres,31,32 illustrating that the signal can be used to follow the mean behaviour of the sample under a change in solution properties. We have not, however, previously been able to attribute a specific shape to each individual nanoparticle recorded.
Here we show how high aspect ratio polyurethane (PU) nanopores can be used to identify the shape of individual nanoparticles in solution. We illustrate that the signals recorded by nanospheres and nanorods are sufficiently different to permit their classification on an individual basis. The analysis is quick, analysing hundreds of particles within seconds. Further, the differences we observe experimentally between nanorod and nanosphere pulse shape corresponds to computational models of predicted pulse shapes, opening the possibility of characterising nanoparticle shape more broadly using large pore RPS coupled to statistical methods. Calibrating the response of the nanopore method makes the process reproducible across several nanopores and days and allows us to accurately measure the ratio of nanospheres and nanorods in solution.
000 g mol−1, analytical standard, 50 wt%, P3143, poly(acrylic acid-co-maleic acid) (PAAMA), Mw ∼ 3000 g mol−1 50 wt%, 416053, purchased from SigmaAldrich, UK. Reagents were prepared in purified water with a resistance of 18.2 MΩ cm. Particles were taken from stock (50 μL) and suspended in PEI (1 mL, 5% in H2O). The solution was placed on the rotary wheel for 30 min. The solution was centrifuged for 5 min at 10
000 rpm, and the PEI solution removed from the particles and replaced with water. The sample was vortexed and sonicated until the particles were fully dispersed. This wash step was repeated twice to ensure all excess PEI had been removed. The PEI-coated particles were suspended in PAAMA (5% in 50 mM NaCl) for 30 min and placed on the rotary wheel. The same process for removing the excess PEI was used. The particles were then stored at 2–4 °C in water. The carboxyl polystyrene particles (158 nm) from Bangs Laboratories, US and the CPC200 (200 nm) from Izon Science, NZ were used without modification. The particles were diluted using 50 mM potassium chloride solution (KCl, >99%, P/4240/60 purchased from Fisher Scientific, UK).
Within the PU setup, two main modes of transport exist. The first is electrophoresis. As all particles have a carboxyl surface chemistry they travel towards the anode on the side of the membrane opposite to the sample. The second is convection, caused by fluid flow under the influence of gravity. Here the pores have a vertical orientation and the sample is placed at the top of the membrane.
34 using the RStudio interface. First a method was developed to reproducibly extract each pulse. This involved extracting 1001 time points containing a detected blockade event, identifying the blockade minimum, then extracting 301 time points with the pulse minimum at point 151. This resulted in 150 data points before and after the pulse magnitude. To account for any base line drift, and variation in noise, each pulse was detrended using the first 50 and last 50 time points.
A file consisting of all extracted pulses was created from each run, as well as a file of extracted pulses normalised to depth 1. Finally, the extracted pulses and normalised extracted pulses were each approximated by quadratic b-splines with fixed knots (using the R package cobs) and the coefficients saved in two additional datasets.
Models were built in both the R base package (glm) and in the package glmnet (for lasso penalised regression). Two calibration runs were used to train the models using solutions of either nanospheres or nanorods alone. Following the calibration data two test solutions were run, again containing just nanorods or just nanospheres. This process was repeated for each nanopore, and in total 5 PU nanopores were tested and in terms of their ability to predict the shape of each nanoparticle recorded within the run using the R package pROC.
A computational model was built to determine the theoretical pulse shapes for the nanospheres and the nanorods as a function of particle position within the pore. Fig. 2a shows the simulated normalised pulse shape for a nanorod and nanosphere as they translocate the nanopore. A clear difference in current-position relationship for nanorods and nanospheres is observed as the particle approaches the sensing zone, left hand side. Whilst a smaller difference is predicted as the particle exists the sensing zone, right hand side. Note that whilst we have modelled the nanorod aligned to the axis of the nanopore (see Experimental section), there is no way to predict that the nanorods have a preferred orientation as they approach the opening. For comparison, Fig. 2b shows the average measured pulse shape for >500 nanospherical and nanorod particles. Note that the difference in pulse sharpness between the predicted and observed pulses is due to the difference between horizontal axes—and number of data points. However in terms of a comparison between pulses for the two particles, the average measured pulse show a similar trend to the simulations, which suggests that it should be possible to determine which type of particle has been detected by using information about the pulse shape. To further illustrate the sensitivity of this technique nanorods of different aspect ratios were run through the setup. Fig. 2c shows the average measured pulse shape for >500 particles for varying aspect ratios. In the current setup nanorods with an aspect ratio below 2 are indistinguishable from nanospheres.
Shape data in the form of spline coefficients for the normalised pulses were obtained using quadratic splines with 25 fixed knots and three fixed points (24 coefficients). The spline coefficients provide a mathematical description of the shape of the denoised pulse, standard Fourier methods of denoising are not suitable as pulse lengths are on the same length scale as noise oscillations. A schematic of this is shown in Fig. 3, but it should be noted that the knots were not evenly spaced as represented, more knots were used along the areas of rapid change in current than in the mostly stable regions at the ends.
Five different types of logistic regression models were then fit to the data. Each used one pair of nanorod and nanosphere runs (the training data) and tested for particle classification on data from the remaining pair of runs (the testing data). The first model, A, used only the maximal depth of the pulse, which is has been shown to be approximately proportional to particle volume.29 The second model, B, used the vector (pulse width) of the blockade data described elsewhere.6 The third model, C, used spline coefficients for the extracted signals. This also includes the size data from model A. The fourth model, D, used spline coefficients for b-splines fit to the extracted signals after they were all normalised to have depth 1. This model only examines the shape of the resistive pulse as the size (depth) data has been removed. Models B, C and D were fit for each pore using lasso penalised regression (glmnet package in R). This means that the variables selected (i.e., with non-zero coefficients) for these models were different for different pores. The final model, E, used a fixed set of spline coefficients from the normalised signal data based on the results of the multiple t-tests (coefficients 9, 10, 16, 17, 18 and 19), and fit the model using standard logistic regression. The signal processing and modelling methods are illustrated in Fig. 3.
The performance of the models was compared using the maximum of sensitivity + specificity along the ROC curve on the testing data for each model. Here sensitivity is the percent of correctly classified nanospheres and specificity is the percent of correctly classified nanorods, and thus a perfect score would be 200. For each set of data (two runs of each particle type), there were four choices (two choices of nanosphere run and two of nanorod run) of training and testing data. We refer to the models built and tested in these four different ways as replicates. The model scores for these four replicates across 5 nanopores are recorded in Table 1. Important to note that for all models, each pore had to be calibrated before it could characterise particle shape, i.e. a sample of pure nanorods or nanospheres were passed through each pore first. It should also be noted that the calibration on one pore would not predict particle shape on another, data not shown. We attribute this to differences in pore structure and reproducibility of the manufacturing process, as well as its stability over the course of the experiment (Fig. S3†).44
| Model A | Model B | Model C | Model D | Model E | |
|---|---|---|---|---|---|
| Pore 1 | |||||
| Replicate 1 | 162 | 108 | 168 | 135 | 161 |
| Replicate 2 | 167 | 112 | 174 | 151 | 162 |
| Replicate 3 | 140 | 106 | 144 | 120 | 149 |
| Replicate 4 | 139 | 104 | 141 | 115 | 152 |
| Pore 2 | |||||
| Replicate 1 | 148 | 121 | 145 | 119 | 150 |
| Replicate 2 | 140 | 119 | 137 | 118 | 147 |
| Replicate 3 | 148 | 116 | 149 | 127 | 144 |
| Replicate 4 | 137 | 127 | 143 | 124 | 137 |
| Pore 3 | |||||
| Replicate 1 | 163 | 164 | 159 | 142 | 135 |
| Replicate 2 | 159 | 142 | 155 | 121 | 114 |
| Replicate 3 | 160 | 158 | 161 | 137 | 127 |
| Replicate 4 | 161 | 162 | 153 | 148 | 143 |
| Pore 4 | |||||
| Replicate 1 | 128 | 136 | 167 | 163 | 162 |
| Replicate 2 | 117 | 146 | 155 | 148 | 153 |
| Replicate 3 | 142 | 156 | 156 | 142 | 136 |
| Replicate 4 | 149 | 135 | 145 | 128 | 131 |
| Pore 5 | |||||
| Replicate 1 | 122 | 131 | 134 | 135 | 133 |
| Replicate 2 | 129 | 108 | 125 | 123 | 124 |
| Replicate 3 | 129 | 129 | 136 | 134 | 136 |
| Replicate 4 | 135 | 107 | 127 | 123 | 121 |
| Mean pores 1–5 | 144 | 129 | 149 | 133 | 141 |
| Standard deviation | 15 | 20 | 13 | 13 | 14 |
For model A (depth data only), the overall mean score was 144 ± 15 with a standard deviation of 15, with approximately 74% of nanospherical particles identified as nanospheres, and 70% of nanorod particles identified as nanorods. As mentioned above, particle volume is approximately proportional to pulse magnitude, so particles of different sizes can be easy identified. However (see Fig. 1S†), whilst the nanospherical particles are on average bigger, there is considerable overlap between the sizes, which we see reflected in scores less than 200.
Model B had an average score of 129 ± 20. Interestingly, this performed worse on average than the model using only depth, indicating that pulse width data is not an effective way to capture shape data, and that the extra width data can in fact lead to model overfitting rather than improved predictive performance. This is also indicated by the high standard deviation, which shows that instrument collected width data is not a reliable measure of shape.
The average score of Model C at the optimal classification point was 149 ± 13. Thus we see that on average the inclusion of the pulse shape yields a somewhat better classification model for nanomaterial shape then depth alone, and also with a marginally lower standard deviation, indicating a somewhat more reliable method of model building. Model C also has the powerful advantage that it removes the need for any pre-knowledge of the particle size for the materials analysis.
Model D had an average score of 133 ± 13. Whilst this score is lower than for models A and C (both of which contain depth data), the reproducibility, as measured by standard deviation, is the same as that for the full spline model, and better than for the model of depth only or the model with depth plus width.
To ascertain which parts of the spline model were key in identifying the particle shape, we performed multiple t-tests with a Bonferroni correction on each spline coefficient over all pairs of nanorod and nanosphere data taken for each pore (4 pairs per pore). Fig. 4 visualises the key spline segments i.e. the parts of the curve whose shape shows a significant difference between nanospheres and nanorods. For all pores, the fall in current, i.e. as the particle enters the pore, has splines that show significant features for particle shape (green segments). This corresponds to the theoretical pulse shapes shown in Fig. 2. In some cases, pores 2, 3 and 4, the second part of the pulse i.e. the part that is recorded as it traverses the pore also shows a significant ability to measure particle shape. This most likely reflects the variation in pore-to-pore structure and lack of reproducibility of the manufacturing process. It also suggests that pores of different shapes may have overall better ability than conical pores to distinguish particle shape.
For a model that used shape only data, the best performance was for the model with fixed normalised spline coefficients, model E, which had an average score of 141 ± 14, corresponding to about 75% sensitivity (correct nanosphere classification) and 66% specificity (correct nanorod classification). In general, specificity was below sensitivity at the maximum of the sum in all cases, which we attribute to the effect of variation in nanorod orientation creating a wider range of possible pulse shapes, which are thus more difficult than nanospheres to characterise with a dichotomous model.
Whilst the identification of particle shape for a pure sample is powerful, a final demonstration of the power of the RPS-LRM was to mix different ratios of nanospheres and nanorods together and ask the model to predict the ratio. For this test we used Model D, tested on mixtures of large nanospherical particles 235 nm in diameter and nanorods shown in Fig. 1. Model D was used here instead of model E, to illustrate that with no knowledge of the key splines vectors any user can download the R code here and apply it to their RPS experiments. Particles of different sizes were chosen for this experiment for two reasons. First, this allowed us to use size to determine if each given pulse corresponded to a nanorod or a nanosphere, thus giving us access to ground truth about particle shape. Recall that model D uses only shape data from depth normalised pulses, thus we can compare the performance of this model for shape detection of particles in solution with the ground truth as given by size. Thus our predicted ratio of particles using size and the Model D classification ratio should be comparable. A second reason was to demonstrate that a model fit using nanospheres of one size could be used to determine the shape of nanospheres of a different size, thus further demonstrating that shape can be studied independent of size.
These samples were analysed and the results shown in Fig. 5a. Model D gave a linear response, that is, as the ratio of nanorods increased the model identified an increase, with an impressive linear relationship (R squared value of 0.9264). The gradient here is greater than 1, due the imperfect specificity and sensitivity reported earlier. It is important to note that the presence of both particles in solution does not affect the model's ability to predict particle shape. This data can also be used to estimate the optimal sensitivity and specificity possible with this current experimental setup when the run conditions for all particles in the testing and training sets are identical. Using model D we get a model with sensitivity of 72% and specificity of 91%. I.e., 91% of nanospherical (large) particles are identified by just their shape data as nanospheres, and 72% of nanorod (small) particles are identified by just their shape data as nanorods (blue curve in Fig. 5c).
We expect that future variations of this work will push the RPS-LRM method further towards a rudimentary individual particle imaging system.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c9nr09100a |
| This journal is © The Royal Society of Chemistry 2020 |