Determination of β-Agonist Residues in Animal-Derived Food by a Liquid Chromatography-Tandem Mass Spectrometric Method Combined with Molecularly Imprinted Stir Bar Sorptive Extraction

A novel clenbuterol molecularly imprinted polymer (MIP)-coated stir bar was prepared and applied to the determination of six β-agonists in animal-derived food. Characterization and various parameters affecting adsorption and desorption behaviours were investigated. The extraction capacities of clenbuterol, salbutamol, ractopamine, mabuterol, brombuterol, and terbutaline for MIP coating were 3.8, 2.9, 3.1, 3.5, 3.2, and 3.3 times higher, respectively, than those of the NIP coating, respectively. The method of MIP-coated SBSE coupled with HPLC-MS/MS was developed. The recoveries in pork and liver samples were 75.8–97.9% with RSD from 2.6 to 5.3%. Limits of detection (LODs) and limits of quantification (LOQs) were 0.05–0.15 μg/kg and 0.10–0.30 μg/kg, respectively. Good linearities were obtained for six β-agonists with correlation coefficients (R 2) higher than 0.994. These results indicated the superiority of the proposed method in the analysis of β-agonists in a complex matrix.


Introduction
β-adrenergic agonists are a kind of chemical synthesized from the benzyl ethanol amines. ese compounds have been extensively used as growth promoters in livestock production, with advantages of increasing the animal's lean meat, improving feed efficiency, and reducing the number of days to market [1,2]. Consumption of animal tissues with high β-adrenergic agonist content by humans leads to symptoms such as muscle tremor, muscle pain, nausea, and dizziness and can pose a serious threat to life [3]. Owing to illegal β-agonist use and its potential hazardous effects on human health, some countries have banned its use in farm animals [4]. e wide variety of these compounds and their extensive use in farm animals give rise to enormous challenges of daily supervision and routine monitoring for the government. It is therefore necessary to develop a validated method with characteristics of high throughput and high sensitivity for the detection and quantification of β-adrenergic agonists.
Recently, several analytical methods have been developed for the determination of β-agonist drugs in foodstuffs and animal urine [5][6][7][8][9][10][11][12][13]. e analytical methods are mainly based on two different immune analysis techniques [6,7,9,11] and on the chromatographic technique [5,8,10,12,13]; both have advantages and drawbacks. Immune analysis techniques such as enzyme-linked immunosorbent assay (ELISA) exhibit the features of rapid and simple analysis but come with drawbacks of giving rise to false positives. Chromatographic techniques such as HPLC, GC-MS, and LC-MS/MS require expensive instrumentation and cumbersome sample pretreatments but have advantages of high accuracy and sensitivity. High-performance liquid chromatography-mass spectrometry/mass spectrometry (HPLC-MS/MS) is a feasible method that is commonly used due to its advantages such as sensitivity and accuracy [14]. For all mass spectrometry analytical procedures, the sample matrix effect is a major problem that influences the accuracy and precision of the testing results [15,16]. erefore, sample preparation is important and crucial, occupying nearly seventy percent of the total analysis time.
Research into sample preparation techniques has become a frontier of modern analytical chemistry, motivated by the demand for preparation methods with advantages of reduced or eliminated solvent, high selectivity, simplicity, speed, and so on. Traditional techniques such as liquid-liquid extraction (LLE) [17,18], immune-affinity chromatography [19], and solid-phase extraction (SPE) [20,21] are time-consuming and labour-intensive. e advanced and nonexhaustive techniques such as SPME (solid-phase microextraction) are solvent-free sample preparation techniques and combine sampling, analyte enrichment, and purification into one step, therefore substantially reducing the cost and total time of analysis [22]. e SPME contains various extraction concepts such as stirrers, vessels, coated fibres, and membranes [23]. Among these techniques, stir bar sorptive extraction (SBSE), first introduced in 1999 by Baltussenetal [24], is an environmentally friendly sample preparation technique using a stir bar coated with polydimethylsiloxane (PDMS). However, SBSE with a PDMS-coated stir bar shows low selectivity and specificity for target analyte and is only suitable for the analysis of nonpolar and weakly polar compounds [25]. Molecularly imprinted polymer (MIP) coatings for SBSE were proposed to solve the problem of polar molecules and selectivity. MIPs are prepared by copolymerization of functional monomers and cross-linkers in the presence of template molecules (target analytes). Extraction of the template molecules leaves the recognition sites in the polymers with specific size, shape, and functional group complementarity to the original print molecule [26]. us, MIPs show selective recognition and high affinity to the template molecule and its structurally related compounds [27]. Moreover, the imprinted polymers are robust, stable, and resistant to a wide range of pH, temperature, and solvents [28,29]. e MIP-SBSE technology that combines MIP with SBSE plays an important role in the aspects of removing or decreasing sample matrix effect. e combination of MIP-SBSE technology with liquid chromatography has also been reported [30][31][32][33][34][35][36][37][38]. However, the combination of MIP-SBSE technology with liquid chromatography-tandem mass spectrometry has been rarely reported. e aim of this present work was to prepare an MIP coating on the glass stir bar and apply it to the selective SBSE of clenbuterol and its analogue compounds. e preparation parameters and extraction conditions, which affected the extraction performance of the stir bar, were studied in detail. e extraction capacity and selectivity of the clenbuterol MIP-coated stir bar were also evaluated. en, a method of determination of structural-related β-agonists by MIP-SBSE followed by HPLC-MS/MS was developed. e proposed method was successfully applied to the analysis of β-agonists in animal-derived food.

2.2.
Apparatus. e coating surface of the stir bar was characterized by a VEGA3 scanning electron microscope (TESCAN, Czech Republic), and the composition of the coating was investigated using a Fourier transform infrared (FT-IR) spectrometer (Nicolet iS5, ermo, USA). Measurements were taken using a TSQ Quantum Ultra triple quadrupole mass spectrometer ( ermo, USA) equipped with an electrospray ion source and coupled with an Accela liquid chromatography ( ermo, USA), equipped with a Hypersil GOLD (100 mm × 2.1 mm i. d., 3 μm) column.

HPLC-MS/MS Conditions.
e mobile phase consisted of mobile phase A (0.1% formic acid in water) and mobile phase B (acetonitrile) at the flow rate of 0.3 mL·min −1 using a gradient elution program. e solvents were degassed by an in-line vacuum degasser. e gradient conditions were as follows: initial time, 10% B; 1 min, 40% B; 3 min, 95% B; 5 min, 100% B, 5.1 min, 10% B; and re-equilibration to 7 min. e column was kept at 35°C. e injection volume was 10 μL. e compound was ionized in an electrospray ionization (ESI) instrument operated in positive mode. Capillary voltage was set at 3.5 kV. e sheath gas and aux gas were nitrogen at the pressures of 45 and 10 arb, respectively. e collision gas was argon at the vacuum pressure of 1.5 mTorr. Desolvation and source temperatures were set at 300°C and 350°C, respectively. To ensure high specificity and sensitivity, two or three of the most abundant precursor → fragment ion transitions (selected reaction monitoring, SRM) were observed for each analyte. SRM transitions, collision energies, and fragmentor for each analyte are given in Table S1.

Stir Bar Preparation.
e glass capillary (2 mm diameter, 100 mm length) was cut into 20 mm long substrate. en, a 16 mm iron rod was inserted in the substrate, and both of its ends were sealed by flame to generate a stir bar (20 mm × 2 mm). e pretreatment of the stir bar was carried out by successively placing it into 1 mol/L sodium hydroxide for 10 h and 0.1 mol/L hydrochloric acid for 1 h. en, the stir bar was dried in an oven at 150°C for 1 h and immersed in a 30% (v/v) KH-570 solution in acetone for 3 h. en, it was pulled out and cleaned with methanol. Approximately 139.0 mg clenbuterol and 0.20 mL MAA were dissolved in 5 mL acetonitrile. e solution was mixed thoroughly and incubated overnight at room temperature. en, 1.60 mL EGDMA and 20.0 mg AIBN were added, and the mixture was degassed by an ultrasonic device for 5 min.
en, a silylated stir bar was inserted into a glass tube, and 0.4 mL of solution was transferred into the tube. e glass tube was purged with nitrogen for 3 min and sealed by flame. e polymerization was accomplished in a water bath at 70°C for 24 h. e coating procedure was repeated two times. en, the stir bar was pulled out by peeling off the outside glass tube and washed with 10% (v/v) acetic solution in methanol to remove the template molecules. e nonimprinted polymer-coated stir bar was prepared under the same conditions but without the addition of template clenbuterol in the synthesis.

Extraction Experiment.
A round-bottom flask was used to conduct an extraction experiment in order to reduce the wear of the stir bar. e schematic diagram for the experiment is shown in Figure S1. e stir bar was immersed in solution and stirred at 500 rpm for 60 min at room temperature. Subsequently, the stir bar was removed from the sample solution, inserted in a 200 μL glass vial with a conical insert, and desorbed with 100 μL 10% (v/v) acid in water solution by an ultrasonic bath for 10 min. en, 10 μL desorption solution was injected into HPLC-MS/MS for analysis.

Sample Preparation.
Samples of pork and liver were purchased from a local supermarket. A total of 2.0 g of each homogenized sample was spiked with β-agonist mixed standard solution.
e spiked sample was extracted with 10 mL 2% (v/v) trichloroacetic acid solution in an ultrasonic bath for 10 min. e extraction procedure was repeated twice. e combined extracted solution was centrifuged at 5000 rpm for 5 min, and the pH of the supernatant was adjusted to 9-10 by 1 mol/L sodium hydroxide. en, the solution was centrifuged again, and the supernatant was transferred to another Teflon centrifugal tube and extracted by 10 mL isopropyl alcohol/ethyl acetate (60 : 40, v/v) two times. e organic layers were combined and dried with nitrogen steam. en, 5 mL toluene was added to dissolve the residual for SBSE.

Preparation of the Stir Bar.
e pretreatment and silanization of the glass bar were carried out to immobilize the MIP coatings on the bar. Soaking of the glass bar in alkaline solution is performed in order to expose the maximum number of silanol groups on the surface to promote silylation with the silane coupling agent. us, double bonds were present in the molecular structure, forming an active site for the MIP coating. After the overnight self-assembled polymerization of the template molecule and functional monomer through intermolecular hydrogen bonding interactions, the cross-linking agent and initiator were added to the prepolymerization solution, and thus, robust coatings were formed on the glass bar through thermal initiation polymerization. en, the cavity structure with high selectivity for a target molecule was formed after eluting the template molecule by breaking the hydrogen bonding interaction between the functional monomer and template. e preparation schematic is illustrated in Figure 1. e experimental parameters of MIP coating preparation such as the polymerization reagents and their proportion, temperature, and time were optimized. e experiment finally confirmed that the optimum mole ratio of template (clenbuterol), functional monomer (MAA), and cross-linker (EGDMA) was 1 : 4 : 16. e coatings of the MIP-coated stir bars prepared under optimum conditions were uniform and compact with a certain thickness (as shown in Figure S2).

Characterization of MIP Coatings.
e morphological structure of the stir bar coating surface was characterized by scanning electron microscopy. Figure 2 shows the surface structure of NIP coatings and MIP coatings under the magnification of 200x and 5000x. It is obvious that both MIP coatings and NIP coatings show homogeneous surfaces (Figures 2(a) and 2(c)), whereas the MIP coating exhibited a more porous structure than the NIP coating (Figures 2(b) and 2(d)). is is because the porous surface of the MIP coatings was more beneficial for analyte adsorption and desorption.
e molecular type and structure of the stir bar coatings were investigated by infrared spectroscopy. As shown in Figure 3, the absorption peaks of the NIP coating (Figure 3(b)), MIP coating before eluting template (Figure 3(c)), and MIP coating after eluting template (Figure 3(d)) show no obvious differences except for the peak at 1520 cm −1 , which was only present in the spectra of template (Figure 3(a)) and MIP coating before eluting template (Figure 3(c)). e absorption peak at 1520 cm −1 was derived from a C � C stretching vibration in the benzene ring of the template (clenbuterol). erefore, it was reasonable to conclude that the template (clenbuterol) did not participate in the polymerization but only interacted with the functional monomer by hydrogen bonding.

Extraction Capability.
e adsorption capacity of the MIP-coated stir bar was studied with various clenbuterol standard solutions in toluene through static adsorption experiments.
e NIP-coated stir bar was used for comparison. e extraction solution volume was 5 mL. e extraction time and desorption time were 60 and 10 min, respectively. As shown in Figure 4, the extraction amount of both the MIP-coated stir bar and NIP-coated stir bar increased with increasing clenbuterol concentration in the 1.0-60.0 μg/L range, and the extraction amounts were balanced over the concentration of 35.0 μg/L. e saturated adsorption amounts of the MIP coating and NIP coating were 67.0 and 17.8 ng, respectively. e adsorption amount of the MIP coating was approximately 3.76 times higher than that of the NIP coating. e excellent adsorption capacity of the MIP coating could be attributed to the molecular imprinting function of the clenbuterol template. e cavities of the MIP coating have a special affinity to the template, while no such cavities were present in the NIP coating that interacts with the template by nonspecific sorption only.

Selectivity Evaluation.
Clenbuterol and its structural analogues including ractopamine, salbutamol, mabuterol, terbutaline, and brombuterol were applied to extraction evaluation. Two reference compounds, namely, benzylalcohol and acrylamide, were used for comparison. All analytes were prepared individually with toluene to avoid competitive adsorption. e concentrations of all analytes were 10 μg/L, and the extraction time was 60 min. As shown in Figure 5, the extraction amounts of the MIP coatings for clenbuterol, salbutamol, ractopamine, mabuterol, terbutaline, and brombuterol were 3.8, 2.9, 3.1, 3.5, 3.3, and 3.2 times higher, respectively, than those of the NIP coatings. By contrast, there was hardly any sorption of both kinds of coatings for benzyl-alcohol and acrylamide, which show different characteristics and structure. e selectivity of the MIP coating was closely related to the shape and size of the cavity, and the strength of interaction between the target molecules and the binding sites [39]. e steric shape complementary between the analytes and the recognition sites of MIPs plays a key role in the molecular recognition process. As a result, the MIP-coated stir bars exhibited higher selectivity toward the template of clenbuterol than its structural analogues. e analogues of salbutamol, mabuterol, terbutaline, and brombuterol with molecular shape and size similar to the template showed better imprinted effect than ractopamine, which had a little different molecular structure and functional groups. e reference compounds of benzyl-alcohol and acrylamide with less similarity with the template could hardly be extracted by both coatings.

Investigation of SBSE Conditions.
To optimize the SBSE conditions, factors that affected the extraction results, such as solvent, time, stirring speed, and temperature, were also investigated in this work. e clenbuterol template was dissolved in the acetone, chloroform, acetonitrile, methanol, and toluene solvents to prepare individual solutions with a concentration of 10 μg/L. e effect of extraction solvents on the extraction capacity was investigated. As shown in Figure S3, the maximum extraction amount was attained in the case of toluene solvent. e above results indicate that the primary driving forces behind the rebinding process, namely, the hydrophobic interaction and hydrogen bonding, were strongly related to the polarity of the extraction solvents.
e hydrogen bonding was liable to be formed in noncovalently imprinted polymers, and higher adsorptive capacity was obtained in the solvents with weaker polarity. us, the matrix effect of water should be removed, which influenced the absorption efficiency of MIP-coated stir bars. e water, methanol, methanol-acetic acid (9 : 1, v/v), and water-acetic acid (9 : 1, v/v) solvents were selected for the investigation of the solvent effect on the desorption. e desorbed amounts were 26.5, 25.6, 27.5, and 30.5 ng, respectively. erefore, water-acetic acid (9 : 1, v/v) was chosen as the best desorption solvent. e extraction and desorption amounts increased with increasing time, and the times at which each reached equilibrium were 60 min and 10 min, respectively (Figures S4 and S5).  erefore, 60 min and 10 min were adopted for extraction and desorption procedures.
e stirring speed varying from 0 to 800 rpm was investigated. e results indicate that the extraction amounts increased with the stirring speed and the maximum extraction amount was obtained at 300 rpm. Beyond 300 rpm, no obvious increase was attained. us, 300 rpm was chosen as the optimal stirring speed. In addition, the change in the extraction temperature exerted no obvious influence on the extraction amounts. For convenience, the adsorption experiments were carried out at room temperature.

Durability Investigation.
To investigate the stability of the coating, the MIP-coated stir bar was immersed into conventional organic solvents, such as toluene, acetonitrile, methanol, n-hexane, ethanol, chloroform, ethyl acetate, dimethylsulfoxide, and tetrahydrofuran. After stirring for 24 h, no decomposition and exfoliation were obtained, indicating that the coatings of MIP-coated stir bar possessed excellent antisolvent fastness properties. e MIP-coated stir bar could be reused after the analyte was removed by methanol-acetic acid (9 : 1, v/v). As the extraction was performed in a round-bottom flask instead of a flat-bottom bottle, the wearing of the stir bar during the stirring has been dramatically reduced. e coating was kept intact, and no decrease in the extraction efficiency was obtained after extraction for 60 times.

Method Validation.
Assay of clenbuterol and its structural analogues by MIP-coated SBSE coupled with HPLC-MS/MS was performed. e blank pork and liver samples were spiked with six β-agonists and taken for analysis. e linear dynamic ranges, correlation coefficients, recoveries, LODs, LOQs, and reproducibility under the optimized experimental conditions are listed in Table 1. e linear range was determined to be 0.50-35 μg/kg for clenbuterol and 1.0-35 μg/kg for ractopamine, mabuterol, brombuterol, or terbutaline. Good linearity was obtained for six β-agonists with correlation coefficients (R 2 ) higher than 0.994. e limits of detection (LODs, signal-to-noise ratio   Journal of Analytical Methods in Chemistry equal to 3) and limits of quantification (LOQs, signal-tonoise ratio equal to 10) were in the ranges of 0.05-0.15 μg/kg and 0.10-0.30 μg/kg, respectively. e recoveries for the spiked pork and liver samples ranged from 75.8% to 97.9%. e RSD for intraday and interday assays of six β-agonists varied from 2.6 to 5.3% (n � 5).

Application of the MIP-Coated Stir Bar.
To validate the selectivity of the method for assay of complex sample matrix, the blank pork samples were spiked with six β-agonists at LOD level and 5.0 μg/kg. Furthermore, the real sample collected from the local market was analyzed by the proposed method. e multiple reaction monitoring (MRM) chromatograms showing the transitions corresponding to the quantifier ions of six β-agonists are depicted in Figure 6. β-Agonists could not be quantitatively assayed in the blank pork sample (Figure 6(a)). However, the six β-agonists were accurately analyzed based on MIP-SBSE (Figure 6(b)). Using the MIP-coated stir bar, the matrix interference has been largely eliminated. As shown in Figure 6(c), the peak at the retention time of 2.33 min was detected in the channel of m/z      277.0/203.1, suggesting that clenbuterol might exist in the pork sample. e extracts were subjected to further analysis using accurate mass AB Sciex Triple Quad 5500 LC/MS/MS in MS/MS/MS scan mode, and the content of clenbuterol in the positive pork sample was determined to be 1.31 μg/kg. e analysis of the pork sample spiked with each β-agonist at 5.0 μg/kg is shown in Figure 6(d). e appearance time for the six β-agonists was within 3 min, and the total analysis time for each specimen was only 5 min. e performance comparison of the proposed method based on MIP-SBSE and other analytical methods for the determination of β-agonists is shown in Table 2. Our work was advantageous in terms of detection throughput, sensitivity, robustness, and reproducibility. e proposed method has been successfully applied to the determination of trace β-agonists in foodstuff of animal origin with improved enrichment ratio of target analytes and reduced sample matrix effect.

Conclusions
A new clenbuterol MIP-coated stir bar for the assay of β-agonists with similar structures was constructed. e new coatings possessed homogeneous and porous morphology and exhibited excellent extraction capability and selectivity for six β-agonists.
e MIP-coated SBSE coupled with HPLC-MS/MS has been successfully applied to the analysis of pork and liver samples spiked with β-agonists, and the recoveries ranged from 75.8 to 97.9%. In comparison with other analytical methods, our work possesses significant advantages of high throughput, high sensitivity, high efficiency, and low detection limit. e approach described herein may be extended to the preparation of a series of selective MIP-coated stir bars with appropriate template molecules for complex sample analysis.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.