Edward D. H. Mansfielda,
Yash Pandyaa,
Ellina A. Muna,
Sarah E. Rogersb,
Inbal Abutbul-Ionitac,
Dganit Daninoc,
Adrian C. Williamsa and
Vitaliy V. Khutoryanskiy*a
aSchool of Pharmacy, University of Reading, Whiteknights, Reading, Berkshire RG6 6AD, UK. E-mail: v.khutoryanskiy@reading.ac.uk
bISIS Spallation Neutron Source, Science and Technology Facilities Council, Rutherford Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, OX11 0QX, UK
cTechnion – Israel Institute of Technology, Faculty of Biotechnology and Food Engineering, Israel
First published on 8th February 2018
Functionalising nanoparticles with polymers has gained much interest in recent years, as it aids colloidal stability and manipulation of surface properties. Here, polymer-coated thiolated silica nanoparticles were synthesised by self-condensation of 3-mercaptopropyltrimethoxysilane in the presence of hydroxyethylcellulose. These nanoparticles were characterised by dynamic light scattering, small angle neutron scattering, Nanoparticle Tracking Analysis, Raman spectroscopy, FT-IR spectroscopy, thermogravimetric analysis, Ellman's assay, transmission electron microscopy and cryo-transmission electron microscopy. It was found that increasing the amount of hydroxyethylcellulose in the reaction mixture increased the nanoparticle size and reduced the number of thiol groups on their surface. Additionally, by utilising small angle neutron scattering and dynamic light scattering, it was demonstrated that higher concentrations of polymer in the reaction mixture (0.5–2% w/v) resulted in the formation of aggregates, whereby several silica nanoparticles are bridged together with macromolecules of hydroxyethylcellulose. A correlation was identified between the aggregate size and number of particles per aggregate based on size discrepancies observed between DLS and SANS measurements. This information makes it possible to control the size of aggregates during a simple one-pot synthesis; a prospect highly desirable in the design of potential drug delivery systems.
Thiol-functionalised silica nanoparticles have attracted particular attention due to the opportunities afforded by the SH-functional groups present on their surface. Thiolated silica is highly effective in removing mercury from aqueous solutions,5,6 used as a support to prepare palladium-containing catalysts,7 and as intermediate materials to prepare modified stationary phases for chromatography.8 The ease of functionalisation of these particles, for example fluorescent labelling, also provides opportunities for in vivo imaging and cell culture experiments.9
Thiolated silica nanoparticles can be easily synthesised by using 3-mercaptopropyltrimethoxysilane (MPTS) as a precursor, which undergoes either self-condensation or co-condensation with tetraethyl orthosilicate in water–ethanol mixtures, catalysed by the addition of a base.10,11 Recently we reported a modified synthetic approach,12,13 using self-condensation of MPTS in dimethyl sulphoxide which resulted in highly stable, monodisperse, sub-100 nm nanoparticles that adhered to ocular tissues through disulphide bridge formation with cysteine residues present on the mucosal surface. The presence of thiol groups on the surface of these nanoparticles also allows facile functionalisation via reactions with maleimide, iodoacetamide, and alkyne derivatives. Successful functionalisation of MPTS nanoparticles has been demonstrated through their fluorescent labelling, PEGylation,14 and POZylation.15
The synthesis of silica nanoparticles with prolonged stability against aggregation in aqueous media is important for their applications in a number of areas. To this end, one approach is to functionalise the surface of the particle with polymers, such as with PEG.16–18 However, PEG precursors used for covalent attachment to silica surfaces are often relatively expensive.
In the present study, we have successfully synthesised thiolated silica nanoparticles with a protective layer of hydroxyethylcellulose in a one-pot reaction whilst maintaining a high number of free thiol groups for further functionalisation. These particles have been fully characterised using various analytical techniques. Depending on the concentration of hydroxyethylcellulose in the reaction mixture, it was possible to prepare either individual silica nanoparticles coated with polymer chains or aggregates of nanoparticles, bridged together with polymer chains.
Vitrified specimens were prepared at a controlled temperature and at water saturation in the Vitrobot (FEI, Netherlands), and kept in liquid nitrogen until examination. Cryo-TEM analysis was done with a Tecnai T12 G2 TEM (FEI, Netherlands) operating at 120 kV. Images were recorded digitally on a Gatan UltraScan 1000 camera using the DigitalMicrograph software (Gatan, U.K.). Images are recorded in the low-dose imaging mode to minimize beam exposure and electron-beam radiation damage.20,21
SANS experiments used the Sans2d neutron diffractometer at ISIS, UK. The instrument set up and incident wavelength employed a Q range of 0.004–0.4 Å−1. Q is defined in eqn (1);
(1) |
Initially the parent thiolated silica sample (Y1) was modelled to a sphere only, as it was synthesised without the presence of any polymer, and therefore provided baseline information including the scattering length density (SLD) and radius of the particle core, which was used for further analyses. Following this, the parent particle was also fitted to a core–shell model, along with the remaining data for samples decorated with HEC (Y2–Y5). Data obtained from particle characterisation (including radius, polydispersity and density) were used for the fitting.
DLS is a well-established and widely used technique for determining the size and polydispersity of nanoparticles within a sample. The size distribution profiles for the silica nanoparticles (parent and HEC-functionalised) are in Fig. 1a with numerical values in Table 1. As can be seen, increasing the concentration of HEC added to the reaction mixture increases the particles hydrodynamic radius, suggesting the hydroxyethylcellulose is responsible for this effect.
Sample | % mass of HEC (w/v) | Yield (%) | Hydrodynamic radius, nm (DLS) | PDI (DLS) | Hydrodynamic radius, nm (NTA) | Particle radius, nm (TEM)a | Graftingb density (μg nm−2) | Surface thiol content (μmol g−1) |
---|---|---|---|---|---|---|---|---|
a Exemplary TEM images are shown in Fig. SI1.b It should be noted that grafting density calculations used the surface area of an individual nanoparticle obtained by TEM analysis, as this was most accurate to the raw particle size. It should also be noted, that given the likely presence of free HEC in the dispersion, there may be some error in these calculations. | ||||||||
Y1 | 0 | 58 ± 2 | 25 ± 1 | 0.102 | 37 ± 1 | 15 | — | 552 ± 185 |
Y2 | 0.1 | 61 ± 3 | 27 ± 2 | 0.143 | 50 ± 5 | 19 | 0.058 | 348 ± 41 |
Y3 | 0.5 | 60 ± 1 | 54 ± 2 | 0.210 | 59 ± 3 | 19 | 0.064 | 267 ± 83 |
Y4 | 1 | 62 ± 3 | 97 ± 11 | 0.236 | 78 ± 3 | 17 | 0.202 | 189 ± 24 |
Y5 | 2 | 65 ± 4 | 190 ± 30 | 0.241 | 99 ± 4 | 20 | 0.144 | 207 ± 20 |
A complimentary particle sizing technique, NTA, was also used. NTA maps the movement of individual particles over a given time and generates a diffusion coefficient based on the number of pixels the particles cross.23 From this diffusion coefficient, a particle size is obtained using the Stokes–Einstein equation (eqn (2));
(2) |
DLS is a light scattering technique where a laser is shone through a suspension of particles which causes the laser light to scatter (Rayleigh scattering). Fluctuations in the scattering, due to particles moving under Brownian motion, are measured at specific time intervals, which are then fitted to an autocorrelation function. A diffusion coefficient is then determined using Mie theory (based on this autocorrelation function and the refractive index of the sample), followed by a particle size using the Stokes–Einstein equation.
NTA on the other hand is a particle tracking technique, and allows for a diffusion coefficient to be measured by analysing the movement of individual particles in a specific environment. By tracking individual particles, undergoing random Brownian motion from frame to frame, the average spatial displacement of the particles per unit time can be calculated, and this displacement can be related to the hydrodynamic diameter of the particles through the Stokes–Einstein equation (eqn (2)). Although translational Brownian motion is a three-dimensional process, it is possible to use a one, two, or three dimensional diffusion coefficient to determine a particles hydrodynamic diameter as described by eqn (3):
(3) |
It should be noted that there is no assumption of 2-dimensional movement of particles. All particles are assumed to be moving freely in all 3 dimensions while the measurement is sampling the projection of each x, y and z component of that movement onto the xy observation plane.
As can be seen from the size distributions, there are discrepancies between the results generated by each sizing technique. These are quantified by changes in the average sizes observed (Table 1). The decrease in concentration observed by NTA is likely due to the increasing size of particles, which form aggregates. The larger aggregates are formed in solutions with greater concentrations of HEC, where the macromolecules may act as bridges to bind individual thiolated silica nanoparticles. The larger size of these aggregates leads to a smaller number in suspension, which reduces their total concentration measured by NTA.
Table 1 shows additional DLS and NTA data, along with the particle yield produced and concentration of reactive thiol groups on their surface (as determined by Ellman's assay). As can be seen, yields from the synthesis were relatively similar, and no significant difference was noted (p > 0.05). This could be related to the presence of free HEC, which could not be removed by dialysis. In terms of sizing (based on the z-average produced by DLS and the mean particle size from NTA), there is a discrepancy. This is due to the way the data is collected in each case. The z-average from DLS is a mean particle size based on the intensity, and is bias towards a larger particle size (due to increased scattering from larger particles compared to smaller particles); however, the mean particle size recorded by NTA is the mean of all the particles measured, and therefore shows no bias towards polydisperse samples. Fig. 1 shows the size distributions for both DLS and NTA (a and b respectively), showing a range of the particle sizes in the suspension. The peak of each curve respectively corresponds to the hydrodynamic size quoted in the table.
It is also apparent from Table 1 that polydispersity index (PDI) increases as more HEC is added to the reaction mixture (from 0.102 with no HEC to 0.241 with 2% HEC). This suggests either particle aggregation or free HEC was present in the particle dispersion.
Fig. 2 Thermal decomposition of thiolated and HEC-functionalised nanoparticles, and pure HEC measured using TGA. |
By calculating the weight loss in each sample and subtracting it from the mass difference in the unfunctionalised silica nanoparticles, it is possible to determine the grafting density of HEC in the particle dispersion (Table 1). It can be seen that increasing the amount of HEC in the reaction mixture increases the grafting density on the particle surface. The exception to this is Y5, where (when compared to Y4) a decrease in grafting density can be observed, from 0.202 to 0.144 μg nm−2. Although the % mass loss increases in the 2% HEC, the grafting density appears lower. A likely reason for this is that the HEC is more densely incorporated into the internal structure of the particles, and thus more protected from thermal degradation compared to HEC bound loosely on the particle surface.
FT-IR and FT-Raman spectroscopies confirm the presence of HEC in the particle suspension. Spectra are shown for both FT-IR and Raman in Fig. 3 and 4, respectively.
Fig. 4 FT-Raman spectra of nanoparticles prepared from reaction mixtures with different HEC content. Spectrum of pure HEC is also included for comparison. |
As shown in Fig. 3, adding HEC into the reaction mixture introduces a broad band in the region of 3200 cm−1 to 3600 cm−1. This is attributable to the O–H groups present in HEC. As the concentration of HEC in the nanoparticles increases, a second band at 2880 cm−1 and 2900 cm−1 becomes more distinctive and provides further evidence that increasing the amount of HEC in the reaction mixture increases polymer loading onto the nanoparticles. Another feature evident in the spectra for nanoparticles that contain HEC (as well as pure HEC) is the band at around 1640 cm−1, attributable to bound water within the cellulose.24 The large peak observed at 1000–1200 cm−1 arises from the Si–O–Si linkage, as seen in other publications using the same particles.14,15
Fig. 4 shows FT-Raman spectra of HEC-functionalised silica nanoparticles. These data demonstrate the presence of both thiol groups (at 2559 cm−1) and disulphide bonds (at 490 cm−1) in the nanoparticles synthesised from MPTS. It also exhibits many more similarities to the spectra recorded for uncoated organosilica as reported in previous work.14,15 The presence of disulphide bonds results from thiol groups being oxidised during the synthesis procedure.12,13
The Raman spectra were normalised to the maximum intensity of the C–H band. This allows the data from each of the batches to be overlaid to facilitate comparisons between nanoparticles synthesised with varying concentrations of HEC. Fig. 4 illustrates that the spectra differed significantly between the shifts of 1500 cm−1 and 0 cm−1. These differences appear to follow a trend that may be caused by the increasing amount of HEC. The differences are especially pronounced between 1500 cm−1 and 1000 cm−1. Indeed, it was found that as the intensity of –SH groups decreases, the intensity of aliphatic C–C stretches increases. Given that C–C exists mostly in the HEC, and SH in the silica, this confirms this hypothesis.
From the spectroscopic and thermal data, it is clear that the particle size increase resulting from increasing amounts of HEC in the reaction is due to the presence of this polysaccharide on the particle surface, or bound into the particle core. As the volume of the nanoparticle precursor is the same, it appears that the HEC macromolecules are binding into the core of the particle, causing a larger cluster for which the particles can form around. By doing so, the concentration of surface thiol groups is lower (shown by a smaller peak at 2550 cm−1 in the FT-Raman spectra, and molar concentration shown by Ellman's assay).
Fig. 5 Scattering cross-sections for Y1, Y2, Y3, Y4, and Y5. The insert demonstrates differences in sizes observed by DLS, TEM, and SANS. |
Table 2 provides a summary of the sizes for each population of particles used in this study, based on modelling the SANS data to a core–shell sphere (see ESI† for a full breakdown). Data obtained from previous characterisation techniques (TEM for the core particle size, and DLS for the hydrodynamic size, Table 1), was used as a guideline in the fitting process. The core-radius and shell thickness parameters were obtained directly from modelling the SANS data. The total diameter was calculated by adding together the radius and shell thickness (to obtain total radius), and multiplied by 2. The values calculated matched the values obtained by TEM (Table 1), confirming their accuracy. The value for number of particles per aggregate was estimated from the aggregate volume (by DLS) divided by volume of hydrated particles at each HEC concentration (from SANS data plus an estimated 20 nm for the hydration layer around each particle). Further details (and caveats) of the estimation can be found in the ESI.† For the purposes of this discussion, “aggregate size” refers to the hydrodynamic diameter obtained by DLS, and “particle size” refers to the size calculated from SANS modelling.
Sample | Core radius (nm) | Shell thickness (nm) | Total diameter (nm) | Estimated particles per aggregate | |
---|---|---|---|---|---|
Y1 | 15.2 | 0.5 | 31 | 1 | |
Y2 | 17.2 | 0.5 | 35 | 1 | |
Y3 | 17.9 | 3.0 | 41 | 6 | |
Y4 | 14.6 | 3.0 | 35 | 44 | |
Y5 | 17.5 | 3.4 | 42 | 230 |
Comparing the particle sizes obtained for the unfunctionalised silica (i.e. particles with no HEC present, Y1), generated by TEM/SANS to that of DLS, a discrepancy of 19 nm can be observed (31 nm in SANS/TEM, and 50 nm in DLS). Previously we have reported this to be due to the hydration shell surrounding the particle; DLS provides a radius of hydration, whereas SANS and TEM will provide the raw particle size (i.e. size without the hydration shell). A hydration shell thickness of 10 nm, similar to the value obtained here, was observed.14,25
From the DLS data and TEM images for the remaining samples (Y2–Y5), it is clear that the particles synthesised with increasing amounts of HEC have a tendency to form larger aggregates. It is possible that, given the increasing amount of HEC in the solution, that several particles are interacting with a single HEC chain, which acts as a bridge, thus leading to increased aggregation. By correlating the DLS and SANS data, it is possible to determine the number of particles per aggregate (taking into account a 10 nm radius hydration shell around each particle (Table 2, and ESI†)). Interestingly, a relationship was found between aggregate size (and the consequent number of particles per aggregate) and the mass of HEC used in the reaction (Fig. SI3†).
The determined aggregate volume appears to correlate well with the cryo-TEM images shown in Fig. 6. Given that both the particles per aggregate and aggregate size correlate very well with the concentration of HEC (with R2 values of >0.99), it is likely that HEC is interacting with growing particles, and anchoring several particles together. Given that the values appear relatively consistent, with little error within the same system (from the small PDI values from DLS), it could be that an equilibrium exists, whereby HEC induces the formation of a small cluster of silica particles. Based on the studies of these particles using different methods, we can propose a scheme for particle aggregation dependent on HEC content in the reaction mixture, illustrated in (Fig. 7).
Fig. 6 Cryo-TEM images of thiolated silica nanoparticles functionalised with (a) 0.1% w/v HEC, (b) 0.5% w/v HEC, and (c) 1% w/v HEC. All scale bars = 100 nm. |
Fig. 7 Illustration of the increased packing of particles in aggregates as HEC content increases in the reaction mixture. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ra08716k |
This journal is © The Royal Society of Chemistry 2018 |