RSC Advances: Paper
RSC Advances: Paper
RSC Advances: Paper
Nowadays there is growing concern about the adulteration of synthetic anorexics in natural slimming
products (NSP) due to the increasing use of weight-reducing aids. Here, we report surface-enhanced
Raman spectroscopy (SERS) coupled with chemometrics for rapid discrimination and detection of
sibutramine (SIB) and its analogues (mono-desmethylsibutramine, MDS; di-desmethylsibutramine, DDS)
in konjac pressed candy (KPC), which is a NSP used in the treatment of obesity. It has been
demonstrated that the two characteristic peaks at 1060 cm1 (in SERS spectrum of SIB) and 1407 cm1
(in SERS spectrum of MDS) were used to discriminate between SIB, MDS and DDS. Moreover, the two
spectra for KPC with MDS and KPC with DDS only exhibited some slight differences. After the second
derivative transformation, the differences were much more clear, especially for the peak at 1407 cm1 in
the SERS spectrum of KPC with MDS. Meanwhile, partial least-squares regression was used for
Received 2nd September 2014
Accepted 11th December 2014
quantitative analysis of the samples. The limit of detection for SIB, MDS and DDS was 5.00 108 M,
5.00 107 M and 1.00 106 M, respectively. Compared with previously reported studies, the method
DOI: 10.1039/c4ra09636c
is simple (no pretreatment), less time-consuming and inexpensive, and can be applied to rapidly screen
www.rsc.org/advances SIB and its analogues in NSP.
5886 | RSC Adv., 2015, 5, 5886–5894 This journal is © The Royal Society of Chemistry 2015
View Article Online
HPLC,7 ultra-HPLC,8 liquid chromatography-mass spectrometry spectral data and develop qualitative and quantitative models
(LC-MS),9 gas chromatography-mass spectrometry (GC-MS),10 for SIB, MDS and DDS.
ion mobility spectrometer (IMS),11 IMS-MS,12
chemiluminescence-molecularly imprinted sol–gel lms
(CL-MIFs),13 ow injection electrochemiluminescence (FI- 2 Experimental
ECL),14 capillary electrophoresis,15 infrared spectroscopy,16 and 2.1 Chemicals and materials
so on. Due to the structural similarity and the trace amount (0.1
Sodium citrate (Na3C6H5O7$2H2O, 99.8%), silver nitrate
mg per capsule11 or tablet) in pharmaceutical formulations, the
(AgNO3, 99%), KI, KNO3, KBr, K2CO3, KCl, ammonium
identication of SIB and its analogues usually requires
formate (HPLC grade) and formic acid (HPLC grade) were
sophisticated sampling beforehand separation/extraction and
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
This journal is © The Royal Society of Chemistry 2015 RSC Adv., 2015, 5, 5886–5894 | 5887
View Article Online
2.4 Theoretical calculation between three structural analogues in spiked sample, second
derivative transformation and smoothing were also applied
The molecular vibrational spectra of the investigated
when necessary to capture a more signicant variance. The
compounds (including SIB, MDS and DDS) were simulated by
spectral region used for the DA and PLS was from 1800 to 900
density functional theory (DFT) within Gaussian 09 program
cm1.
package. On the basis of the molecular geometries optimized by
B3LYP/6-31+G** method, the frequency calculations were
carried out at the same level to validate the minima and to 2.7 SERS detections of SIB, MDS and DDS in spiked samples
obtain the vibrational information. Meanwhile, the interaction Briey, SIB with different amounts was rst spiked into one
between Ag atom and investigated compounds was simulated gram of KPC powders, and then the mixture was dissolved in 10
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
through the basis sets of 6-31+G** and LanL2DZ, which per- mL ultrapure water by sonication, statically cooled down and
formed to non-metal (C, N and O) and Ag atoms. Based on the ltered. 10 mL of ltrate was added into 90 mL of Ag NPs colloid.
optimized geometries, the Raman spectra calculations were The Raman spectra were recorded with the same condition as
carried out at the same computational level. The assignments of the detection of analytes in aqueous solution. Meanwhile,
calculated vibrational modes were made with the help of Gauss accuracy of the method was determined by standard addition
View 5.0. The calculated wavenumbers of all normal modes technique. Known amounts of reference standard SIB in a range
were scaled by a factor of 0.987 to correct systematic errors of of low, medium and high concentrations were added to the KPC
the theoretical methodology. powders and analyzed under the above measuring conditions.
The experiments were replicated six times for each concentra-
tion and the accuracy was calculated as the % of analyte
2.5 SERS measurements of SIB and its analogues
recovered. MDS and DDS was also tested using the same
In brief, 5.00 mM SIB solution was prepared and diluted with method, respectively.
ultrapure water to the predetermined concentrations. Then,
10 mL of the analyte solution was added into 90 mL of Ag NPs
colloid in a 1.5 mL centrifuge tube, followed by shaking for 1 3 Results and discussion
min. The mixture was rst sucked into a capillary glass tube, 3.1 Characterizations of Ag NPs
and then xed onto a glass slide. The Raman spectra were UV-vis spectroscopy and transmission electron microscopy
recorded using 532 nm laser with 10 mW power and 10 (TEM) were applied to characterize the size and morphology of
objective. The collecting time was 10 s with 2 rounds of accu- the synthesized Ag NPs. The maximum absorption of the
mulations and the pinhole aperture was 25 mm. The experi- measured UV-vis spectrum of Ag NPs colloid provides the
ments were repeated six times for every concentration. MDS and information on the diameter of particles, moreover, its full
DDS were also tested using the same method. width at half-maximum (FWHM) is used to estimate the NPs
dispersion.35 As shown in Fig. S1A,† the maximum UV-vis
2.6 Spectral data analysis absorption for Ag NPs colloid was at 409 nm, and the FWHM
was approximately 90 nm, which indicated that the size
Multivariate statistical tools (DA and PLS) were adopted to
dispersion of Ag NPs was narrow. Meanwhile, TEM image
analyze the spectral data with TQ Analyst soware (7.2.0.161
clearly displayed that the diameter of Ag NPs was 42 9 nm
Release, Thermo Fisher Scientic). Through computing the
(Fig. S1B†). These results revealed that the size of Ag NPs was
unknown's distance from each class center in Mahalanobis
mainly homogeneous and the shape of the most of NPs was
distance units, DA classication technique can be used to
spherical.
determine the class or classes of known materials which are
most similar to an unknown material. The closer each distance
value is to zero, the better is the match. Aer the material was 3.2 Optimization of experimental conditions for SERS
successfully classied by DA, PLS was applied to adulterated The optimum of the experimental conditions for SERS detection
samples (which only contain one synthetic anorexic per sample) is related to various parameters including excitation wave-
to obtain a calibration based on the level of adulteration. For length, aggregating agent, pH value, etc. In order to select the
PLS analysis, spectra of 40 samples with ve different concen- optimal detection condition, a series of tests were performed in
trations of constituent were used. Among these samples, thirty this study.
samples were used for calibration and ten samples were used Excitation wavelength. It is known that the basic mecha-
for validation. The calibration data were used to create the nisms of SERS mainly include charge-transfer (CT) mechanism
method calibration model, while the validation data were used and electromagnetic mechanism. For the CT mechanism, it is
to validate the model. Evaluation of the models was made related to the chemical complexation of the analyte molecules
through correlation coefficient (R) and root-mean-square error to the surface of the metallic NPs. In the case of surface
of prediction (RMSEP). The higher the R value or the lower the complexation, electronic absorption of the complexed analyte
RMSEP value is, the better predictability the model has. Before may shi with respect to the free molecule, and even overlap the
the analysis, spectral data were normalized using a standard excitation wavelength leading to resonant effect.36 Fig. 1A
normal variate (SNV) algorithm and mean centered to correct showed the SERS spectra of SIB absorbed on Ag NPs for three
for instrumental effects on the data. To further discriminate excitations (532, 633 and 780 nm). The intensity of peak at 1691
5888 | RSC Adv., 2015, 5, 5886–5894 This journal is © The Royal Society of Chemistry 2015
View Article Online
cm1 in the spectrum recorded with the 532 nm laser line was indicator of particles-aggregating degree), it has been noted that
stronger than those obtained by 633 and 780 nm laser. The extra addition of anions can enhance the SERS signal of SIB to a
surprising differences on the peak intensities should be caused certain extent. The SERS activation induced by anions can be
by the charge-transfer Ag complex with an electronic absorption ascribed to two effects: one is the aggregation of Ag NPs
in the spectral region of the 532 nm laser. The UV-vis spectrum producing “hot spots”; the other is modication of the Ag
of SIB (Fig. S2a†) showed a strong absorption at 223 nm, very far surface chemistry, resulting from the addition of strongly
from the 633 and 780 nm to have any resonant effect. However, binding anion. As clearly presented in Fig. 1B, the Br, as the
the UV-vis spectrum of the SIB complexed with Ag NPs displayed more strongly binding anion, gave the greatest enhancement
a broad shoulder that extended up to 500 nm region among different kinds of anions, while Cl gave the least
(Fig. S2b†). Thus, the 532 nm laser was in pre-resonance with enhancement for SIB. So, KBr was chosen as the optimal
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
the absorption of the complex, and was used as excitation aggregating agent to conduct the following optimization
source in the following tests. experiments. Additionally, the stability of Ag NPs colloid with
Aggregating agent. The aggregating agent is crucial for the addition of KBr was also examined, as shown in Fig. S5.†
application of SERS using metallic NPs. It is usually an elec- Clearly, the absorbance of Ag NPs with KBr is very stable and the
trolyte, and can induce the NPs aggregation by increasing the Raman intensity of SIB adsorbed on Ag NPs is nearly unchanged
ionic concentration inside the metallic colloid and decreasing within an hour, which ensures the measurement reliability.
the electrostatic barrier. Fig. S3† showed the UV-vis absorption Substrate concentration. In order to maximize the Raman
spectra of Ag NPs colloid blended with different aggregating signal enhancement, Ag NPs colloid at different concentrations
agents (the concentration of each aggregating agent is 10 mM). was obtained by the treatment of centrifugation or dilution.
The absorbance of Ag NPs colloid decreased aer mixed with KI, Fig. S6† revealed that the Raman intensity of peak at 1691 cm1
KNO3, KBr, K2CO3 and KCl, respectively, and maximum continuously increased with the Ag NPs concentration in range
absorption at 409 nm shied to 415 nm, 415 nm, 414 nm, 407 of 0.007–0.2 nM. It could be ascribed to the increasing “hot
nm and 406 nm, respectively. From the previous literature, the spots”, which could maximize the localized electric eld and
evolution of the absorbance of Ag NPs colloid could be related to then greatly amplify Raman signals of SIB in the colloid. The
the aggregation of colloidal NPs.41 When aggregating agents further increase of substrate concentration led to the decrease
containing anions are mixed with Ag NPs colloid, the energy of the Raman peak intensity. It was because that Ag NPs colloid
states on particles' surface are affected by chemical contribu- with high concentration was easily precipitated. Therefore,
tions of the anions and the inter-band transitions between NPs optimal concentration of Ag NPs colloid was 0.2 nM.
through anions. The change of particles' surface energy states pH value. When the Ag NPs colloid at different pH values was
caused by additional anions is regarded as energy hybridization used, it was found that considerable changes occurred in the
effect on particles' surface.42 For this reason, the maximum Raman spectra of SIB. As shown in Fig. S7,† the intensity of
absorption of Ag NPs colloid changed in the presence of addi- band (at 1691 cm1) obtained at pH 7.4 was stronger than that
tional anions compared with the absence of additional anions. at pH 2.5, but this band nearly vanished at pH 11.5. These
The different interactions of Ag NPs with aggregating agents phenomena may be related to variations in molecular proton-
were also conrmed by TEM and dynamic light scattering ated and deprotonated forms as pH values under ambient
(Fig. S4†). conditions37–39 or slight differences on adsorption geometries of
Fig. 1B displayed a series of SERS spectra of SIB absorbed on target molecule on nano-roughened metallic surface.40 To ach-
Ag NPs colloid to which the same amount of KI, KNO3, KBr, ieve maximum SERS intensity, we chose Ag NPs colloid with pH
K2CO3 and KCl was added, respectively. Through comparison of 7.4 as substrate for the detection of SIB. Meanwhile, for the
the six spectra (the feature peak at 1691 cm1 was chosen as the
Fig. 1 (A) SERS spectra of 5.00 104 M SIB in Ag NPs colloid excited by different excitation wavelengths. (B) Effects of aggregating agents on
SERS intensity of 1.00 105 M SIB.
This journal is © The Royal Society of Chemistry 2015 RSC Adv., 2015, 5, 5886–5894 | 5889
View Article Online
other two target analytes (MDS and DDS), similar optimization ring, respectively. In comparison to MDS and DDS, the SIB
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
results were also obtained (Fig. S8–11†). showed the characteristic peak at 742 cm1, corresponding to
the symmetry and asymmetry stretch vibrations of the N1–C2
bond. A close analysis of spectra of MDS and DDS solids
3.3 SERS discrimination of SIB, MDS and DDS
revealed no dramatic differences except the two weak feature
Fig. 2 showed the molecular structures of SIB, MDS and DDS. peaks at 1311 cm1 (deformation vibration of the C2–N1–C7
The structural difference between three molecules is only in the bond) and 1390 cm1 (deformation vibration of the N–H bond
amounts of substituent (methyl) of amino-group. To discrimi- and in-plane deformation vibration of C–H bond from phenyl
nate between three structural analogues, the Raman spectra of ring) for MDS.
the three analytes were determined. The observed Raman bands Fig. 3B presented the normal Raman spectra of three analyte
with their assignment accomplished with the help of theoretical solutions. Compared with the Raman spectra of solids (Fig. 3A),
calculations were listed in Table 1. it was observed that all bands of three analytes were too weak to
Fig. 3A displayed the Raman spectra of three analyte solids. discriminate between three structural analogues.
Their dominant bands appeared at 1595 cm1, 1100 cm1 and Fig. 3C exhibited the SERS spectra of three analytes in Ag NPs
629 cm1, which was assigned to the stretching vibration of C–C colloid. In comparison to the normal Raman spectra (Fig. 3B),
bonds from phenyl ring, deformation vibration of the CH2 we found that the intensities of some bands were obviously
groups and out-of-plane deformation vibration of the phenyl
Table 1 Raman (experimental and calculated) wavenumbers (in cm1) of SIB, MDS and DDS with their vibrational assignmenta
Raman
5890 | RSC Adv., 2015, 5, 5886–5894 This journal is © The Royal Society of Chemistry 2015
View Article Online
Fig. 4 (A) SERS spectra of KPC with SIB, KPC with MDS and KPC with
DDS (2.5 mg g1), (B) second derivative transformations of above
Fig. 3 The Raman spectra of SIB, MDS and DDS: (A) solid, (B) water spectra and (C) Cooman plot for the classification of KPC with SIB, KPC
solution and (C) in Ag NPs colloid. with MDS and KPC with DDS.
enhanced. The enhancement effects can be explained by two bonds in phenyl ring and CH3 scissoring vibration). By further
mechanisms: (1) the electromagnetic mechanism based on the inspecting spectra of MDS and DDS, we observed that there was
coupling of the radiation eld with surface plasmons of Ag NPs; one obvious feature peak at 1407 cm1 (deformation vibrations
(2) the chemical mechanism based on a change in the molec- of the NH bond and CH bonds in phenyl ring) in spectrum of
ular polarizability arising from the interaction with the metal MDS.
substrate. As compared to MDS and DDS, SIB displayed two In order to apply SERS spectroscopy to differentiate the three
characteristic peaks including 1060 cm1 (stretching vibration analytes even in complicated matrix, we have to identify the
of C2–N1–C8 bond) and 1415 cm1 (stretching vibrations of C–C feature bands that allow an easy discrimination of SIB and its
This journal is © The Royal Society of Chemistry 2015 RSC Adv., 2015, 5, 5886–5894 | 5891
View Article Online
Fig. 5 (A) Schematic drawing for the detection of SIB and its analogues from the KPC by the SERS spectroscopy using Ag NPs as Raman sensor,
and the SERS spectra of KPC spiked with (B) SIB, (C) MDS and (D) DDS by the enhancement of Ag NPs colloid.
analogues. In general, a good discriminant peak will be rela- compared the SERS spectra of KPC (blank), KPC with SIB, KPC
tively intense and isolated from nearby peaks which can with MDS and KPC with DDS (Fig. 4A), and discriminant peaks
possibly interfere. Such peaks can serve as markers to identify of each target analyte in KPC were determined. For KPC with
the presence of each compound in mixtures. Herein, we SIB, the peaks at 1060 and 1415 cm1 were relatively uncluttered
Table 2 Real detections and recovery tests of SIB, MDS and DDS spiked in KPC, respectively (n ¼ 6)
Sample Spiked concentration (mg g1) Raman method LC-MS method Raman method RSD (%) LC-MS method RSD (%)
5892 | RSC Adv., 2015, 5, 5886–5894 This journal is © The Royal Society of Chemistry 2015
View Article Online
by nearby peaks, and served as discriminant peaks. For KPC In order to verify the usefulness of the proposed method
with MDS and KPC with DDS, the two spectra exhibited some for the identication and detection of three anorexics in
slight differences. Aer the second derivative transformation practical applications, the selectivity of the SERS method for
(Fig. 4B), the differences looked much more clear, especially for the analytes was rst measured using the structure analogues
the peak at 1407 cm1 which was used to discriminate between (including genistein, daidzein, salidroside, chlorogenic acid,
MDS and DDS. These ndings demonstrated that SERS was able caffeic acid and gallic acid) because these aromatic
to differentiate spectral patterns between the three analogues in compounds (from herbs) might affect the detection of target
KPC samples. To be more intuitive, DA model, as a multivariate analytes in real sample analysis. As shown in Fig. S15,† it is
statistical tool, was also introduced based on the spectra of KPC clearly seen that there is the peak at 999 cm1 in few SERS
spiked with SIB (n ¼ 50), KPC spiked with MDS (n ¼ 50) and KPC spectra of aromatic compounds, but there is not the peak in
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
spiked with DDS (n ¼ 50) without the second derivative trans- the SERS spectra of KPC (blank) (Fig. 4A). The result indicates
formation. DA accurately classied 100% of all samples either that a highly selective SERS strategy is attained with an
as the SIB or its two analogues. Results of DA were shown in a impressive specicity for three target analytes. Subsequently,
Cooman plot (Fig. 4C), which was constructed by plotting the the analysis of three anorexics in real samples was demon-
Mahalanobis distance between each anorexic drug. Eigen strated here by the use of the SERS sensor (Fig. 5). Similar to
analysis, one of the diagnostic checks, denes how much vari- three anorexics in water solution, PLS was also applied for
ation of spectral information is explained by the use of principal quantitative analysis of three anorexics spiked in KPC
components. For this case, eigen analysis revealed that it was (Fig. S16 and S17†). The relative data of recoveries were listed
possible to obtain 99.4% of desired information with only six in Table 2, which indicated that the obtained results were
principal components. The obtained results illustrated that highly consistent with those from the LC-MS method.
SERS combined with DA could be used to differentiate the three However, the SERS method does not need to pretreat samples
analytes intuitively. before assay, and the whole test can be accomplished within
3 min. These conrm that the SERS-based Ag NPs sensor
has a high accuracy and reliability to meet the requirements
3.4 SERS detections of SIB, MDS and DDS in spiked samples of practical application.
Before the detection of trace SIB, MDS and DDS in KPC, the
detection sensitivity to analytes in water using the SERS method
has been tested by the addition of 10 mL analyte aqueous solu- 4 Conclusion
tion into 90 mL of Ag NPs colloid. Fig. S12A† showed SERS
spectra of SIB. With the increase of SIB concentration from 5.00 We have developed a rapid and effective method which can
108 M to 5.00 106 M, the Raman intensity gradually clearly discriminate between SIB and its analogues adulterated
increases. The intensity of the strongest peak at 1691 cm1, in KPC based on SERS data, DA and second derivative trans-
assigned to u(C7–N–C8) + n(C–C) (phenyl ring), was used for the formation. Citrate-reduced Ag NPs were used as the SERS
quantitative evaluation of SIB level and the limit of detection substrate. Optimal excitation wavelength (532 nm), substrate
(LOD) was determined to be 5.00 108 M (equal to 16.7 ng concentration (0.2 nM), pH value (7.4) and aggregating agent
mL1) from three standard deviations above the background. (KBr) for high Raman enhancement were obtained. The SERS
Meanwhile, similar tests were also performed for MDS and DDS. discrimination of SIB and its analogues was achieved by DA
As shown in Fig. S12B and C,† the LOD of MDS and DDS was method, mainly relying on the differences in Raman spectral
5.00 107 M (equal to 151 ng mL1) and 1.00 106 M (equal window between 1000 and 1600 cm1. SIB, MDS and DDS was
to 288 ng mL1), respectively (herein, the intensity of the successfully quantied using PLS model, respectively, and the
stronger peak at 999 cm1, assigned to phenyl ring breathing, LOD for SIB, MDS and DDS was 5.00 108 M (equal to 16.7 ng
was used for the quantitative evaluation of MDS and DDS mL1), 5.00 107 M (equal to 151 ng mL1) and 1.00 106 M
levels). Subsequently, PLS analysis was applied on spectra of (equal to 288 ng mL1), respectively. Moreover, the recoveries of
three target analytes. RMSEP values obtained from the PLS SIB, MDS and DDS in the spiked KPC were in the ranges of 98.8–
model of SIB with different latent variables were shown in 99.7%, 97.2–98.7% and 98.0–99.2%, respectively. Compared
Fig. S13A.† The lowest RMSEP value was achieved when 4 latent with LC-MS method, this new technique was simple, rapid and
variables were used, suggesting that the optimal number of inexpensive (each test only needed a droplet of 10 mL sample
latent variables to construct a PLS model was 4. The optimal and was accomplished within 3 min), and thus can be suitable
number for building PLS model was 3 for MDS and 4 for DDS for high throughput screening analysis, which is performed
(Fig. S13B and C†). PLS prediction result for SIB was shown in with a motorized sample holder stage together with a 96 wells
Fig. S14A† by plotting calculated SIB concentrations against plate to make the method automatic and reduce the total
actual SIB concentrations (R ¼ 0.993; RMSEP ¼ 0.102 mM). For analysis time of batch samples. Although this work is mainly
MDS and DDS, R value was 0.989 and 0.981, and RMSEP value focused on the SERS analysis of SIB, MDS and DDS in KPC, we
was 1.08 mM and 1.22 mM, respectively (Fig. S14B and C†). With have demonstrated that the SERS technique is generally appli-
R values all larger than 0.980, results of PLS analysis suggested cable to detect SIB, MDS and DDS of other NSP such as HST
that PLS could serve as a reliable method to quantify SIB, MDS (Fig. S18–20 and Table S1†). Importantly, the monitoring and
and DDS. tracking of synthetic anorexic drugs illegally adulterated in NSP
This journal is © The Royal Society of Chemistry 2015 RSC Adv., 2015, 5, 5886–5894 | 5893
View Article Online
may greatly reduce the possibility of occurrence of negative 19 E. G. Anastasaki, C. D. Kanakis, C. Pappas, L. Maggi,
effects on public health. A. Zalacain, M. Carmona, G. L. Alonso and M. G. Polissiou,
J. Agric. Food Chem., 2010, 58, 6011–6017.
Acknowledgements 20 F. J. Garcı́a de Abajo, Rev. Mod. Phys., 2007, 79, 1267–1290.
21 D. Kurouski, S. Zaleski, F. Casadio, R. P. Van Duyne and
This work is supported by the National Natural Science Foun- N. C. Shah, J. Am. Chem. Soc., 2014, 136, 8677–8684.
dation of China (no. 21275075) and the Specialized Research 22 H. Zhu, M. L. Du, M. Zhang, P. Wang, S. Y. Bao, M. Zou,
Fund for the Doctoral Program of Higher Education of China Y. Q. Fu and J. M. Yao, Biosens. Bioelectron., 2014, 54, 91–101.
(no. 20113234110001). 23 C. Yuen, W. Zheng and Z. W. Huang, Biosens. Bioelectron.,
2010, 26, 580–584.
Published on 12 December 2014. Downloaded by Carleton University on 12/05/2015 08:23:02.
5894 | RSC Adv., 2015, 5, 5886–5894 This journal is © The Royal Society of Chemistry 2015