Selective Enrichment of Clenbuterol onto Molecularly Imprinted Polymer Microspheres with Tailor-made Structure and Oxygen Functionalities

The noxious clenbuterol misapplied as the feed additive has posed an enormous threat to humans who actively rely on the food chains with high potential of contamination by clenbuterol, such as pork and beef. It is, therefore, highly desirable to develop novel materials and strategies for dealing with the clenbuterol. Herein, functional polymer microspheres prepared by Pickering emulsion polymerization were explored for the selective enrichment of the clenbuterol, and their structure and oxygen functionalities could be tailor-made by a molecular imprinting process. The clenbuterol imprinting was adequately demonstrated to not only increase the particle size (~52 nm vs. ~42 nm) and create cavities for the accommodation of the clenbuterol molecules, but also reduce the oxygen functionalities of the resulting molecularly imprinted polymer microspheres (MIPMs) by approximately 4 at.%, which is believed to correlate with the high specificity of the MIPMs. Various characterization methods were employed to evidence these findings, including scanning electron microscopy, BET measurements, Fourier transform infrared spectroscopy, X-ray photoelectron spectroscopy, and elemental mapping examination. More importantly, the MIPMs showed a markedly superior enrichment capability towards clenbuterol to the counterpart, that is, non-molecularly imprinted polymer microspheres (NIPMs). Compared to the NIPMs without specificity for clenbuterol, the MIPMs exhibited an impressive selectivity to clenbuterol, with the relative selectivity coefficient (k′) values largely exceeding 1, thus corroborating that the useful molecular imprinting led to the generation of the binding sites complementary to the clenbuterol molecule in the size and functionalities. The MIPMs were also employed as the stationary phase to fabricate molecularly imprinting solid-phase extraction column, and the spike recovery was demonstrated to be not significantly decreased even after nine cycles. Furthermore, the reliability of the method was also evidenced through the comparison of the MIPMs prepared from different batches.


Introduction
Clenbuterol, a β 2 -adrenoceptor agonist, has raised worldwide concern due to its widespread applications as remedial drugs for the treatment of bronchial diseases (e.g., asthma, pulmonary emphysema, and chronic bronchitis) [1], and even more extensive attention is paid to its function as "lean meat powder", which can stimulate the transformation of fat tissue to muscular tissue in livestock [2]. Since lean meat is the more favorable food relative to the fat, many illegal producers are powder with the desired range of the particle size. Consequently, this post-treatment procedure is time-consuming and usually leads to irregular particles with a wide distribution of particle size [13]. Precipitation polymerization involves the use of a large amount of the organic solvent, which is undesirable from the environmental viewpoint. Swelling treatment combined with thermal polymerization is multistep, tedious, and time-consuming. By contrast, through the Pickering emulsion polymerization possessing advantages such as emulsion stability and easy operation without further grounding post-treatment, polymer microspheres are readily synthesized with a narrower size distribution. The Pickering emulsion polymerization lies in the utilization of colloidal solid particles with a desired wettability, which results in the formation of oil-in-water or water-in-oil emulsions. They can thus assist in significantly reducing the usage amount of the traditional hazardous organic surfactant for the emulsion preparation, thereby substantially reducing the adverse impact on the environment and humans. Pickering emulsion polymerization has been explored for the preparation of polymer microspheres imprinted with malachite green [32], bisphenols [13], protein [46], λ-cyhalothrin [47] and erythromycin [48]. Nevertheless, no reports can be found on the systematic investigation of the clenbuterol-imprinted polymer microspheres prepared by Pickering emulsion polymerization, e.g., the relationship between their microstructure and surface properties and the adsorption performance toward clenbuterol (including adsorption kinetics, isotherms, thermodynamics, and selectivity, as well as molecularly imprinted solid-phase extraction (MISPE) column applications).

Materials
Fumed silica with a mean diameter of 12 nm (Aerosil 200, 99% purity) was supplied by Evonik Degussa (Evonik Industries AG Inorganic Materials, Hanau, Germany). The initiator, 2,2 -azoisobutyronitrile (AIBN, 99% purity), and the monomer, methacrylic acid (MAA, analytical reagent) were purchased from Tianjin Fucheng Chemical Reagent Factory (Tianjin, China). Clenbuterol, terbutaline, salbutamol, and methyl red were analytical reagents and supplied by Sigma Chemicals Co. Ltd. (Poole, UK). The cross-linking agent, ethylene glycol dimethyl acrylate (EGDMA, 98% purity), was purchased from Aladdin Reagent (Shanghai, China). All of the other chemicals were of the analytical grade and obtained from Tianjin Fucheng Chemical Reagent Factory (Tianjin, China). All of the chemicals were used without further purification unless otherwise stated.

Synthesis of MIPMs
The procedures for the synthesis of MIPMs are illustrated in Figure 1. Typically, the water phase was composed of the monomer, MAA (0.12 mL, 1.41 mmol), Triton X-100 (0.3%, 6 mL, 0.028 mmol), and silica nanoparticles (20 mg, 0.33 mmol), which was suffered to sonication treatment to disperse these silica nanoparticles thoroughly. Separately, the oil phase was prepared by mixing EGDMA (1.88 mL, 9.97 mmol), the template material (i.e., clenbuterol, 10 mg, 0.032 mmol), toluene (0.2 mL, 1.88 mmol) and the initiator (i.e., AIBN, 10 mg, 0.061 mmol) with the aid of sonication treatment for 10 min. The prepared water and oil phases were then mixed by intense agitation for 2 min, leading to the generation of the oil-in-water Pickering emulsion where the silica nanoparticles were dispersed along the interface between oil droplets and continuous water phase. The violent agitation at 70 • C allowed the initiation of the free-radical polymerization of the monomer (MAA) by AIBN in the oil droplets, and in the meanwhile, the cross-linking reaction also proceeded between the polymer chains and EGDMA. The polymerization and cross-linking reactions were performed for 16 h. The finally precipitated product was collected and dipped into a 30% HF solution for 12 h to etch away the silica nanoparticles on the surface of the polymer microspheres, which was then subjected to the Soxhlet solvent extraction with 50% methanol solution for 48 h. Such long extraction duration guaranteed that the template molecules (i.e., clenbuterol) were completely removed from the MIPMs, which could no longer be detected by electrospray ionization mass spectrometry (ESI-MS) and elemental analysis. The resulting polymer microspheres, with the average particle size of around 52 nm, were designated as the MIPMs. The non-molecularly imprinted counterpart (i.e., NIPMs), with the average particle size of about 42 nm, were also synthesized under the same preparation conditions except that no template was added during the Pickering emulsion polymerization. The mole ratio used for the preparation of the MIPMs was the optimized one. The insufficient monomer (MAA) did not give rise to the homogeneous polymer microspheres but bulk polymer, while excess MAA produced the microspheres showing nonspecific adsorption towards clenbuterol. Besides, insufficient crosslinker (EGDMA) could not result in a sufficient number of useful imprinting sites, whereas excess crosslinker caused the template molecules to be deeply embedded in the network and thus reduced the number of active imprinting sites. Mixing the two phases to form oil-in-water droplets, where the oil droplets contained the cross-linking agent (i.e., EGDMA), initiator (i.e., AIBN) and part of clenbuterol molecules (as the molecular imprinting agent) and water phase was mainly comprised of the monomer (i.e., MAA); silica nanoparticles were dispersed along the interface between the oil phase and water phase.  From the SEM images of the MIPM and NIPM samples presented in Figure 2 and 3, respectively, notable differences can be observed. The average particle size of MIPMs was estimated to be approximately 52 nm, larger than that of NIPMs (around 42 nm). The particle size distribution histograms are provided in the insets of Figure 2 and 3. Without the molecular imprinting, many tiny particles were noted on the NIPM surface, in stark contrast to the rather clean surface of the MIPMs. These results imply that the addition of the template molecules (i.e., clenbuterol) could improve the Pickering emulsion polymerization quality, most likely resulting from the Mixing the two phases to form oil-in-water droplets, where the oil droplets contained the cross-linking agent (i.e., EGDMA), initiator (i.e., AIBN) and part of clenbuterol molecules (as the molecular imprinting agent) and water phase was mainly comprised of the monomer (i.e., MAA); silica nanoparticles were dispersed along the interface between the oil phase and water phase. (c) Polymerization of the MAA to form the PMAA cross-linked by EGDMA. (d,e) Precipitation (d) and thorough washing (e) of the cross-linked PMAA microspheres. (g) The wastewater contaminated by clenbuterol. (h,i) Using MIPs to bind clenbuterol selectively (i).

Characterizations
The characterizations of the structural functionalities of the MIPM and NIPM samples were carried out by Fourier transform infrared (FTIR) spectroscopy and X-ray photoelectron spectroscopy (XPS) using an IRAffinity-1S spectrometer (Shimadzu Co., Kyoto, Japan) and XSAM800-XPS equipment (Kratos, UK), respectively. Scanning electron microscopy (SEM) images and elemental mapping images of the MIPMs and NIPMs were obtained using a HITACHI S-4800 Scanning Electron Microscope equipped with an energy dispersive X-ray (EDX) spectrometer (Quantax 70, Bruker Nano GmbH, Berlin, Germany). The Brunauer-Emmett-Teller (BET) measurements of the MIPMs and NIPMs were carried out with a Micromeritics ASAP 2020 nitrogen adsorption apparatus (Norcross, GA, USA). Before the BET test, the MIPMs and NIPMs samples were pre-degassed at 150 • C. The corresponding pore size distribution was determined by analyzing the desorption branch based on the Barrett-Joyner-Halenda (BJH) method. The mass spectrometric analysis was conducted using a 6540 UHD mass spectrometer (Agilent Technologies, Santa Clara, USA).

Adsorption Kinetics
To a glass vial, the MIPMs or the NIPMs (50 mg) were added, followed by the addition of a water solution of clenbuterol (10 mg/L, 10 mL). After the homogenization of the mixture imparted by shaking, the glass vial was equipped with a thermostatic water bath oscillator at an oscillation rate of 120 rpm for the adsorption processing. At different durations, the mixed solution was drawn for the measurements of the clenbuterol concentration in the solution, and the adsorbed amount was calculated on the basis of Equation (1) by LC-MS. The amount of the supernatant injected in the HPLC/MS system was 5 µL, and the test conditions were adopted as follows: electrospray ionization as the ion source, dry gas flow rate of 40 L/min, atomization pressure of 10 psi, and methanol as the mobile phase at the flow rate of 0.5 mL/min.
where C 0 and C t represent the initial clenbuterol concentration and the concentration at time t, respectively. The pseudo-first-order and pseudo-second-order kinetic models were adopted to analyze the adsorption kinetics according to Equations (2) and (3), respectively.
where q e represents the equilibrium adsorption amount (mg/g), q t (mg/g) is adsorption amount at time t, and k 1 and k 2 designate the pseudo-first-order rate constant (1/min) and pseudo-second-order rate constant (g/mg min), respectively.

Adsorption Isotherm
To each of five glass vials, the MIPMs or the NIPMs (50 mg) were added. Then, the clenbuterol solutions at a volume of 10 mL, with different concentrations (10,20,30,40, and 50 mg/L), were respectively charged into the five vials pre-loaded with polymer microspheres. After shaking for homogenization, these mixtures were placed into a thermostatic water bath oscillator, and the adsorption interactions proceeded at the oscillation rate of 120 rpm for 6.5 h which was sufficient for reaching equilibrium as proven by the test on the adsorption kinetics. Three different temperatures were considered, namely 30, 45 and 60 • C corresponding to the temperatures in Kelvin of 303, 318 and 333 K, respectively. After the equilibrium was reached, the resulting mixtures were centrifuged at 4000 rpm for 5 min before the measurement of the supernatants by LC-MS with the same test conditions as those described in the Section 2.4, and equilibrium adsorption amount (q e , mg/g) was calculated based on Equation (4).
where C 0 and C e are the initial and equilibrium concentrations (mg/L) of clenbuterol, respectively, m represents the mass of the MIPM and NIPM adsorbents, and V is the volume (L) of the clenbuterol solution. The commonly-used adsorption isotherm models were considered to analyze the present adsorption test results, in this case, Langmuir and Freundlich isothermal models that were described by Equations (5) and (6), respectively.
where q e and C e are the equilibrium adsorption amount (mg/g) and clenbuterol concentration (mg/L), respectively, k L and k F are the Langmuir constant (L/g) and Freundlich constant (L/g), respectively, q max is the maximum adsorption capacity (mg/g), and 1/n is an empirical parameter related to adsorption intensity.
To assess whether the adsorption is favorable or unfavorable, Equation (7) can be referenced.
where k L is the Langmuir constant (L/g), C 0 presents the initial maximum clenbuterol concentration (mg/L), and R L can be used to indicate the favorability of adsorption process, which are described below:

Adsorption Thermodynamics
To study the inherent energy changes in the adsorption course, three basic thermodynamic parameters were measured, i.e., enthalpy change ∆H 0 , Gibbs free energy change ∆G 0 , and entropy change ∆S 0 . These thermodynamic parameters, ∆H 0 , ∆S 0 and ∆G 0 estimated for the adsorption of clenbuterol were described by Equations (8)-(10): where C e is the equilibrium concentration of clenbuterol in solution (mg/L), q e (mg/g) is the adsorbed amount at equilibrium. T represents the temperature in Kelvin (K), and R is the universal gas constant (8.314 J/(mol K)).

Competitive Adsorption Tests
A glass vial was loaded with the MIPMs or NIPMs (20 mg), and then a 10 mL mixed solution containing clenbuterol, terbutaline, salbutamol, and methyl red (each with a concentration of 1 mg/L) was poured into the glass vial pre-loaded with MIPMs or NIPMs. After homogenization imparted by shaking, the mixture was placed into a thermostatic water bath oscillator, and the adsorption interactions proceeded at the oscillation rate of 120 rpm for 6.5 h which was sufficiently long for making the adsorption system reach equilibrium. Then, the resulting mixture was centrifuged at 4000 rpm for where q e is the equilibrium adsorption amount (mg/g), C e is the equilibrium concentration of the adsorbate in the solution (mg/L), k d represents the distribution coefficient (L/g), k is the selectivity coefficient, x represents competing species against clenbuterol, including terbutaline, salbutamol, and methyl red, and k is the relative selectivity coefficient.

Impact of the Solution as the Medium for the Adsorption of Clenbuterol onto MIPs and NIPMs
Various media were employed to prepare 10 mg/L clenbuterol solutions, including water, a water solution of ammonium acetate (4.6 mmol/L), a water solution of 0.3%TX-100 (4.6 mmol/L), acetonitrile, and an acetonitrile solution of ammonium acetate (4.6 mmol/L). These prepared solutions were separately added into the vial pre-loaded with 50 mg MIPMs or NIPMs. After the mixtures had been homogenized by shaking, they were placed in a thermostatic water bath oscillator with an oscillation rate of 120 rpm/min for the isothermal adsorption for 5 h. The processed clenbuterol solutions were centrifuged for 5 min, and the supernatants were withdrawn for the LC-MS measurement. Based on the equilibrium adsorption quantity calculated, the imprinting factor (IF) was estimated according to the following Equation (14).

Reliability of the Method Used for the Fabrication of MIPMs
A new batch of the clenbuterol-imprinted polymer microspheres was prepared to prove that the present method for the fabrication of MIPMs was reliable, and the newly prepared sample was designated as MIPMs-2. The adsorption isotherm was then measured for MIPMs-2 according to the isothermal adsorption test on MIPMs. Through comparison of the parameters obtained by fitting the adsorption isotherms for the adsorption systems with MIPMs and MIPMs-2, the reliabilities of the present fabrication method and isothermal adsorption test could be assessed.

Test on the Recyclability of the MIPMs
The prepared MIPMs were used as the stationary phase of the molecularly imprinted solid-phase extraction (MISPE) column. Briefly, the deionized water and methanol were firstly used to thoroughly wash the MISPE column walls which were then dried in air. Sieve plates were employed to seal the open ends of the MISPE column, and thus the falling of the MIPMs was prevented. The filled MISPE column was further washed with a mixture of methanol and acetic acid (9/1, v/v) at a volume of 10 mL and then with 2 mL methanol. A water solution of clenbuterol (10 µg/L, 100 mL) subsequently passed through the washed MISPE column at a flow rate of about 5 mL/min under negative pressure. Finally, the MISPE column adsorbed with clenbuterol was successively eluted with a mixture of methanol/acetic acid (9/1, v/v) at volumes of 1 mL, 2 mL, and 2 mL. After the residual solution had been drained out of the MISPE column, the eluted solution was examined using LC-MS based on the test conditions as mentioned earlier. The MISPE test was repeated nine times, and the spike recovery was calculated for each cycle. Consequently, the recycling performance of the MISPE column with the MIPMs was evaluated based on the comparison of the spike recoveries calculated. Figure 1 depicts the main content of this study, that is, investigation of the MIPMs for the enrichment of clenbuterol. More specifically, through the Pickering emulsion polymerization and molecular imprinting technique, the MIPMs were fabricated. The comparison between the prepared MIPMs and NIPMs allowed the clarification of the impact of the molecular imprinting on the polymer microspheres synthesized by Pickering emulsion polymerization. We found that the molecular imprinting could not only change the size of these polymer microspheres but alter their surface functionalities. The molecular imprinting should possess a unique capability to leave recognition cavities as depicted in Figure 1 and evidenced by the BET measurements. The N 2 adsorption-desorption isotherms and BJH pore size distribution plots are provided in Figures S1 and S2 in the Electronic Supporting Information (ESI), respectively. The results obtained via BET measurements are summarized in Table 1. The MIPMs showed a significantly smaller average pore width that was estimated to be around 6 nm relative to that of NIPMs (more than 150 nm), as shown in Figure S2 in the ESI. Notably increased cumulative pore volume and specific surface area can also be measured for MIPMs as compared to NIPMs (Table 1). For NIPMs, the type II isotherm can be assigned ( Figure S1 in the ESI), indicating that NIPMs belong to a non-porous or macroporous material [49], and we hypothesized that the microsphere accumulation could generate gaps among these microspheres, which were tested as the macropores in the N 2 sorption test. By contrast, the N 2 adsorption onto MIPMs results in the formation of the type IV isotherm, which suggests the existence of mesoporous structures due to the generated recognition sites. The non-closure of the hysteresis loop implies incomplete removal of adsorbate from narrow pores [49]. The comparison of the N 2 sorption measurement results between MIPMs and NIPMs revealed the success in the manufacturing of the polymer microspheres via clenbuterol imprinting which left recognition cavities complementary to the size of clenbuterol molecules. Solvent extraction processing was most likely to loosen the structure of MIPMs to a certain extent, resulting in improved specific surface area and pore volume. These beneficial results imparted via clenbuterol imprinting endow the MIPMs with great potential for the clenbuterol enrichment applications. From the SEM images of the MIPM and NIPM samples presented in Figures 2 and 3, respectively, notable differences can be observed. The average particle size of MIPMs was estimated to be approximately 52 nm, larger than that of NIPMs (around 42 nm). The particle size distribution histograms are provided in the insets of Figures 2 and 3. Without the molecular imprinting, many tiny particles were noted on the NIPM surface, in stark contrast to the rather clean surface of the MIPMs. These results imply that the addition of the template molecules (i.e., clenbuterol) could improve the Pickering emulsion polymerization quality, most likely resulting from the clenbuterol-assisted promotion of the emulsion stability. Apart from the molecular imprinting template, clenbuterol might also serve as a surfactant, in that a hydrophobic moiety (-C(CH 3 )) at one end and a hydrophilic moiety (-NH 2 ) at the other exist in the molecular structure of clenbuterol ( Figure S3 in the ESI). Such dual roles of the surfactant and molecular imprinting template played by clenbuterol can impart an advantage of the present strategy for the fabrication of functional polymeric materials for processing clenbuterol. Importantly, small microcracks and micropores could be found on the surface of MIPMs (Figure 2c,d), in contrast to the more seamless surface of NIPMs (Figure 3c), which was an indication of the useful molecular imprinting generated on the polymer microspheres. These findings are also in good consistency with the results obtained via N 2 sorption measurements. These generated microcracks and micropores would help to the selective enrichment of clenbuterol due to the molecular imprinting effect. Elemental mapping images presented in Figure 4a-d confirmed that only C and O elements could be detected in both the MIPMs (Figure 4c) and NIPMs (Figure 4d), indicating that the template molecules were extracted entirely out of the MIPMs since the N element exists in the template ( Figure S3 in the ESI). We also measured the elemental mapping images of the intermediate samples before the procedure of etching treatment with HF and the procedure of the Soxhlet extraction treatment, which are shown in Figure 4a,b, respectively. For both the intermediate products, the N element could be probed, revealing that the successful incorporation of clenbuterol molecules (containing the N element) into the polymer microspheres. Even though the HF is a strong etching agent, it could not exert a significant impact on the clenbuterol template located in the polymer microspheres, as the N element was clearly presented in Figure 4b. Therefore, it is necessary to further extract the clenbuterol template out of the polymer microspheres by an approach following the etching treatment with HF.

Results and Discussion
To this end, we adopted the Soxhlet extraction by considering its ease in operation, albeit with a long extraction duration. After the Soxhlet extraction, the N element could no longer be detected (Figure 4c), confirming that the Soxhlet extraction was powerful for the entire extraction of the clenbuterol template out of the cross-linked polymer microspheres.
Polymers 2019, 11, x FOR PEER REVIEW 9 of 24 end and a hydrophilic moiety (-NH2) at the other exist in the molecular structure of clenbuterol ( Figure S3 in the ESI). Such dual roles of the surfactant and molecular imprinting template played by clenbuterol can impart an advantage of the present strategy for the fabrication of functional polymeric materials for processing clenbuterol. Importantly, small microcracks and micropores could be found on the surface of MIPMs (Figure 2c,d), in contrast to the more seamless surface of NIPMs (Figure 3c), which was an indication of the useful molecular imprinting generated on the polymer microspheres. These findings are also in good consistency with the results obtained via N2 sorption measurements. These generated microcracks and micropores would help to the selective enrichment of clenbuterol due to the molecular imprinting effect. Elemental mapping images presented in Figure 4a-d confirmed that only C and O elements could be detected in both the MIPMs (Figure 4c) and NIPMs (Figure 4d), indicating that the template molecules were extracted entirely out of the MIPMs since the N element exists in the template ( Figure S3 in the ESI). We also measured the elemental mapping images of the intermediate samples before the procedure of etching treatment with HF and the procedure of the Soxhlet extraction treatment, which are shown in Figure 4a,b, respectively. For both the intermediate products, the N element could be probed, revealing that the successful incorporation of clenbuterol molecules (containing the N element) into the polymer microspheres. Even though the HF is a strong etching agent, it could not exert a significant impact on the clenbuterol template located in the polymer microspheres, as the N element was clearly presented in Figure 4b. Therefore, it is necessary to further extract the clenbuterol template out of the polymer microspheres by an approach following the etching treatment with HF. To this end, we adopted the Soxhlet extraction by considering its ease in operation, albeit with a long extraction duration. After the Soxhlet extraction, the N element could no longer be detected (Figure 4c), confirming that the Soxhlet extraction was powerful for the entire extraction of the clenbuterol template out of the cross-linked polymer microspheres.   Functional groups on the MIPMs and NIPMs were then probed by FTIR spectroscopy, and the results are presented in Figure 5a,b. We marked most of the characteristic FTIR absorption peaks and their corresponding functionalities in Figure 5. For the MIPMs, the presence of the carboxyl groups of the poly(methacrylic acid) (PMAA) moiety was most clearly reflected by the FTIR absorption around 3433 and 1728 cm −1 which were indexed to the stretching vibrations of -OH and -COOH, respectively. The FTIR absorptions revealed the methyl, methylene, and methane in the MIPMs at approximately 2991, 2958, 2930, 1456 and 1391 cm −1 which were assigned to the stretching and bending vibrations of the C-H group. We also noted the presence of carboxylate functionalities at 1522 and 1477 cm −1 corresponding to the symmetric and asymmetric stretching vibrations of carboxylate as formed by the partial ionization of carboxyl groups. In comparison between the MIPMs and NIPMs, the apparent difference in the characteristic FTIR absorption was observed due to the stretching vibration of carboxyl groups, as highlighted in Figure 5b. For the MIPMs, this absorption was detected at about 1728 cm −1 , which was shifted to approximately 1719 cm −1 as for the NIPMs. Such a kind of redshift has been elucidated to be the enhancement of hydrogen bonding interactions [50]; in other words, molecular imprinting treatment caused the breakage of some hydrogen bonding between PMAA chains such as COOH . . . HOOC. This weakened hydrogen bonding interactions might also be a result of a decrease in the number of carboxyl groups. The molecular imprinting not only created cavities with the size matched with the clenbuterol molecule but also modulated the surrounding chemical properties, and in this study, the molecular imprinting was demonstrated to change the oxygen functionalities and their hydrogen-bonding interactions. The weakened hydrogen bonding is favorable for the adsorption applications since less diffusion resistance to the adsorbate molecules would be generated by the hydrogen bonding.  Functional groups on the MIPMs and NIPMs were then probed by FTIR spectroscopy, and the results are presented in Figure 5a,b. We marked most of the characteristic FTIR absorption peaks and their corresponding functionalities in Figure 5. For the MIPMs, the presence of the carboxyl groups of the poly(methacrylic acid) (PMAA) moiety was most clearly reflected by the FTIR absorption around 3433 and 1728 cm −1 which were indexed to the stretching vibrations of -OH and -COOH, respectively. The FTIR absorptions revealed the methyl, methylene, and methane in the MIPMs at approximately 2991, 2958, 2930, 1456 and 1391 cm −1 which were assigned to the stretching and bending vibrations of the C-H group. We also noted the presence of carboxylate functionalities at 1522 and 1477 cm −1 corresponding to the symmetric and asymmetric stretching vibrations of carboxylate as formed by the partial ionization of carboxyl groups. In comparison between the MIPMs and NIPMs, the apparent difference in the characteristic FTIR absorption was observed due to the stretching vibration of carboxyl groups, as highlighted in Figure 5b. For the MIPMs, this absorption was detected at about 1728 cm −1 , which was shifted to approximately 1719 cm −1 as for the NIPMs. Such a kind of redshift has been elucidated to be the enhancement of hydrogen bonding XPS technique was further employed to consolidate the above analysis based on the FTIR data, and the results are provided in Figure 6a-d and Table 2. From the XPS survey spectra presented in XPS technique was further employed to consolidate the above analysis based on the FTIR data, and the results are provided in Figure 6a-d and Table 2. From the XPS survey spectra presented in Figure 6a, the atomic percentage of oxygen was estimated, and the results were exhibited by a comparison histogram shown in Figure 6b. The oxygen content was lowered for MIPMs relative to NIPMs, which thereby elucidates that the weakened hydrogen bonding as proven by FTIR spectroscopy was attributed to the decease of oxygen-containing functional groups. High-resolution XPS C1s and O1s core-level spectra are further provided in Figure 6c,d, respectively. The peak around 295 eV assigned to the carboxyl group was significantly weakened for the MIPMs relative to the NIPMs (as highlighted by light blue shading in Figure 6c), which thus suggested that the decrease of the carboxyl groups at least partially caused the lowered oxygen atomic content. The lowered concentration of the oxygen-containing functional groups was also unambiguously ascertained by the XPS O1s spectra (Figure 6d) since the intensity of the primary peak on the O1s spectrum was decreased for the MIPMs as compared to that for the NIPMs. An appropriate decrease in the oxygen group concentration is necessary to impart the MIPMs with high specificity by considering that most of the oxygen-containing functional groups interact with the organic species in a nonspecific way, e.g., through hydrogen bonding interactions [51]. The clenbuterol molecular imprinting led to the generation of the binding sites with an appropriate size and surface properties, facilitating the highly selective enrichment of clenbuterol.  Figure 6a, the atomic percentage of oxygen was estimated, and the results were exhibited by a comparison histogram shown in Figure 6b. The oxygen content was lowered for MIPMs relative to NIPMs, which thereby elucidates that the weakened hydrogen bonding as proven by FTIR spectroscopy was attributed to the decease of oxygen-containing functional groups. High-resolution XPS C1s and O1s core-level spectra are further provided in Figure 6c,d, respectively. The peak around 295 eV assigned to the carboxyl group was significantly weakened for the MIPMs relative to the NIPMs (as highlighted by light blue shading in Figure 6c), which thus suggested that the decrease of the carboxyl groups at least partially caused the lowered oxygen atomic content. The lowered concentration of the oxygen-containing functional groups was also unambiguously ascertained by the XPS O1s spectra (Figure 6d) since the intensity of the primary peak on the O1s spectrum was decreased for the MIPMs as compared to that for the NIPMs. An appropriate decrease in the oxygen group concentration is necessary to impart the MIPMs with high specificity by considering that most of the oxygen-containing functional groups interact with the organic species in a nonspecific way, e.g., through hydrogen bonding interactions [51]. The clenbuterol molecular imprinting led to the generation of the binding sites with an appropriate size and surface properties, facilitating the highly selective enrichment of clenbuterol.   All of the results confirmed the successful preparation of polymer microspheres via Pickering emulsion polymerization and effective molecular imprinting of the as-prepared polymer microspheres, with binding sites created for the enrichment of clenbuterol. The adsorption performance of the MIPMs was then systematically studied and compared with that of the NIPMs, and the results are presented in Figures 7-10, and Tables S1-S3 in the ESI. The adsorption kinetics was firstly investigated (Figure 7 and Table S1 in the ESI). The adsorption of clenbuterol onto both the MIPMs and NIPMs was fast at the beginning 2.5 h and gradually became slow until the equilibrium was reached. There were abundant binding sites unoccupied by the clenbuterol molecules at the initial stage, and when clenbuterol molecules gradually occupied the binding sites, the adsorption showed down. Besides, MIPMs better followed the pseudo-first-order kinetic model relative to the pseudo-second-order kinetic model, which was opposite to the adsorption kinetics in the case of the NIPMs as the adsorbent that better obeyed the pseudo-second-order kinetic model, as judged by the correlation efficient (R 2 ). Such a discrepancy between the adsorption kinetics as for the adsorption systems with the MIPM and NIPMs was assumed to stem from the reduced content of oxygen-containing functional groups in the MIPMs relative to the NIPMs, which weakened the chemical interactions between oxygen groups (especially carboxyl groups, as evidenced by Figures 5 and 6) of the adsorbent and clenbuterol via electron sharing, transferring and exchanging. The better fitting to the pseudo-second-order kinetic model has been presumed to be probably originated from the chemisorption [52]. The reduction of the oxygen functional groups on the MIPMs weakened the chemical interactions between the MIPMs and clenbuterol, which enabled the change of the adsorption kinetics from the pseudo-second-order kinetic model (for NIPMs) to the pseudo-first-order kinetic model (for MIPMs). Nevertheless, the adsorption capacity of the MIPMs was overwhelmingly more significant than that of NIPMs (Figure 7), probably resulting from the generation of many molecular imprinting sites complementary to clenbuterol in the size and surface properties. These created sufficient imprinting sites would endow the MIPMs with high binding specificity to the clenbuterol molecules, together with a higher enrichment capacity as compared to the NIPMs, albeit with reduced oxygen-containing functional groups. The isothermal adsorption results are provided in Figure 8, and Table S2 in the ESI and three different temperatures including 30, 45, 60 °C were investigated, with the corresponding results provided in Figure 8a-c, respectively. For the adsorption system with the MIPMs at all the investigated temperatures, the isothermal adsorptions were better followed by the Langmuir isotherm relative to that by the Freundlich isotherm. The Langmuir model assumes that the adsorption sites on the MIPMs are homogeneous and energetically equivalent and that the adsorption proceeds at specific homogeneous sites within the MIPMs [52]. We can thus speculate  The above adsorption studies clarified that the MIPMs exhibited an improved adsorption capacity relative to the NIPMs, indicative of the useful molecular imprinting effect on the polymer microspheres. The promoter for this improvement was the recognition sites as generated by the effective molecular imprinting treatment of the polymer microspheres. Whether these recognition sites showed selectivity to clenbuterol was the next that we would be demonstrated, and the results are shown in Figure 10, and Figure S6 and Table S6 in the ESI. The MIPMs showed much larger selectivity coefficient (k) values in comparison with those of the NIPMs, together with the relative selectivity coefficient (k') values largely exceeding 1. This result is an indication that, in comparison with the NIPMs, the MIPMs exhibited an enhanced selectivity to clenbuterol among potential competing analogs with high structural similarities ( Figure S1 in the ESI). The specific EIS-MS spectra of the mixed analyte solution with clenbuterol, methylene blue, terbutaline, and Salbutamol before and after processing with NIPMs or with MIPMs are presented in Figure S6 in the ESI.
We also fabricated a MISPE column using the MIPMs as the stationary phase, and the spike recovery was calculated. A total of 9 cycles were examined, and the spike recovery calculated for each cycle was compared, with the result provided in Figure 11. The results show that the MISPE column with the MIPMs can be reused for many cycles without significant degradation of the MISPE column performance. The slight decrease in the spike recovery might be caused by the somewhat contamination of the imprinting sites on the MIPMs during the repeated usage of the MISPE column equipped with the MIPMs.  The isothermal adsorption results are provided in Figure 8, and Table S2 in the ESI and three different temperatures including 30, 45, 60 • C were investigated, with the corresponding results provided in Figure 8a-c, respectively. For the adsorption system with the MIPMs at all the investigated temperatures, the isothermal adsorptions were better followed by the Langmuir isotherm relative to that by the Freundlich isotherm. The Langmuir model assumes that the adsorption sites on the MIPMs are homogeneous and energetically equivalent and that the adsorption proceeds at specific homogeneous sites within the MIPMs [52]. We can thus speculate that the molecularly imprinted sites are homogeneously distributed over the MIPMs, and these sites are energetically equivalent.
Additionally, the Langmuir model also indicates the monolayer coverage of clenbuterol onto the MIPMs [53]. Only are the recognition sites not occupied by the clenbuterol molecules; they are available for the adsorption. The results can also reflect that molecular imprinting treatment enables the MIPMs to be more homogenous in the structure by reducing the heterogeneous sites on the MIPMs such as carboxyl groups, as evidenced by the XPS spectra. To examine the reliability of the present method for the preparation of MIPMs, we also prepared the second batch of MIPMs, namely MIPMs-2, for which the adsorption isotherm was also tested, as shown in Figure S4 and Table S3 in the ESI. As expected, the adsorption of clenbuterol on the MIPMs-2 sample also resulted in the adsorption isotherm, which also exhibited the better fitting to the Langmuir model as compared to the Freundlich model. This result thus elucidates the reliability of the present fabrication method for obtaining similar MIPMs from different batches.
On the contrary, the isothermal adsorption of clenbuterol onto the NIPMs better obeyed the Freundlich isotherm in comparison with the Langmuir isotherm. The Freundlich equation is an empirical equation adopted to describe heterogeneous systems, and it is not restricted to the monolayer adsorption [52]. The more abundant oxygen groups on NIPMs (relative to those on MIPMs) make the oxygen groups more likely to be heterogeneously distributed over the NIPMs. The more significant number of oxygen groups make them more possible to exist in different types, with different binding affinities, on the NIPMs. The heterogeneous distribution of the different types of oxygen groups on the NIPMs enables the adsorption isotherm to follow the Freundlich equation.
The difference in the concentration of oxygen groups and the intensity of hydrogen bonding interactions caused a discrepancy between the MIPMs and NIPMs in both the adsorption kinetics and isotherms. The stronger hydrogen bonding interactions between the surface oxygen groups might deactivate their adsorption activity towards clenbuterol molecules which could be more favorably adsorbed onto the sites free from the hydrogen bonding interactions with other oxygen groups. It was also noted that both the adsorption systems with the MIPMs and NIPMs exhibited promotion of the maximum adsorption capacity with increasing the temperature, implying both endothermic adsorption systems. From the thermodynamic analysis results provided in Figure 9 and Table S4 in the ESI, significant parameters including ∆H • , ∆S • , and ∆G • were calculated for the adsorption of clenbuterol onto the MIPMs. The positive ∆H • and ∆S • values indicate the endothermic nature and randomness of the adsorption of clenbuterol onto the MIPMs [54]. The positive ∆G • value reveals that the adsorption reaction requires energy. Since a decrease of ∆G • could increase the reaction rate, higher temperatures were beneficial to the adsorption in this case, in consistence with the aforementioned endothermic adsorption process.
We also investigated the influence of the type of solution used as the medium for the adsorption experiment, and the results are presented in Figure S5 and Table S5 in the ESI. It can be noted that MIPMs exhibited a higher adsorption capacity toward clenbuterol in comparison with NIPMs, irrespective of the media employed. Additionally, a higher capacity of the adsorption of clenbuterol was achieved in the water-based media (including pure water, a water solution of ammonium acetate and a water solution of Triton X-100), in comparison to that in the acetonitrile-based media (including pure acetonitrile and an acetonitrile solution of ammonium acetate). Nevertheless, the higher selectivity was obtained in the acetonitrile-based media, as evidenced by the much-improved imprinting factor for the adsorption system with the acetonitrile-based media relative to the water-based media. These results thus reveal that, in addition to the structure and properties of the MIPMs, the medium used for the adsorption experiment also plays a significant role in deciding the final adsorption performance, including adsorption capacity and selectivity.
The above adsorption studies clarified that the MIPMs exhibited an improved adsorption capacity relative to the NIPMs, indicative of the useful molecular imprinting effect on the polymer microspheres. The promoter for this improvement was the recognition sites as generated by the effective molecular imprinting treatment of the polymer microspheres. Whether these recognition sites showed selectivity to clenbuterol was the next that we would be demonstrated, and the results are shown in Figure 10, and Figure S6 and Table S6 in the ESI. The MIPMs showed much larger selectivity coefficient (k) values in comparison with those of the NIPMs, together with the relative selectivity coefficient (k ) values largely exceeding 1. This result is an indication that, in comparison with the NIPMs, the MIPMs exhibited an enhanced selectivity to clenbuterol among potential competing analogs with high structural similarities ( Figure S1 in the ESI). The specific EIS-MS spectra of the mixed analyte solution with clenbuterol, methylene blue, terbutaline, and Salbutamol before and after processing with NIPMs or with MIPMs are presented in Figure S6 in the ESI.
We also fabricated a MISPE column using the MIPMs as the stationary phase, and the spike recovery was calculated. A total of 9 cycles were examined, and the spike recovery calculated for each cycle was compared, with the result provided in Figure 11. The results show that the MISPE column with the MIPMs can be reused for many cycles without significant degradation of the MISPE column performance. The slight decrease in the spike recovery might be caused by the somewhat contamination of the imprinting sites on the MIPMs during the repeated usage of the MISPE column equipped with the MIPMs.

Conclusions
We have effectively employed Pickering emulsion polymerization and molecular imprinting manufacturing to produce the MIPMs with recognition sites complementary to the clenbuterol template in the size and functionalities. Numerous cavities are generated on the surface of MIPMs as a result of molecular imprinting, in stark contrast to the NIPM surface. On the other hand, Figure 11. Recycling performance of the MIPMs when used as the stationary phase of the MISPE column for processing of clenbuterol solution through the evaluation of the spike recovery for each cycle.

Conclusions
We have effectively employed Pickering emulsion polymerization and molecular imprinting manufacturing to produce the MIPMs with recognition sites complementary to the clenbuterol template in the size and functionalities. Numerous cavities are generated on the surface of MIPMs as a result of molecular imprinting, in stark contrast to the NIPM surface. On the other hand, relative to the clean MIPM surface, many tiny particles are observed on the NIPM surface, which indicates that the presence of the clenbuterol template during the Pickering emulsion polymerization helps to enhance the polymerization quality attributed to the clenbuterol-aided stabilization of the emulsion. Besides, the existence of the clenbuterol template also enables the average particle size of the resulting MIPMs to be increased as compared to that of the NIPMs; the clenbuterol template can help to stabilize the oil droplets, even in larger sizes in the continuous water phase. Interestingly, molecular imprinting also causes a reduction of the oxygen-containing functional groups (such as carboxyl groups) on the MIPMs, relative to those on the NIPMs, which is beneficial for the selective uptake of the clenbuterol in the subsequent adsorption tests since too many oxygen groups will normally debase the specificity. The MIPMs show much better adsorption performance towards clenbuterol in comparison with the NIPMs, and more importantly, they also exhibit an excellent selectivity to clenbuterol among some potential analogs with high structural similarities. The recyclability of the MIPMs in the clenbuterol enrichment is also manifested through running the MISPE column equipped with the MIPMs nine times. Furthermore, the reliability of the fabrication method is also verified by the comparison of two batches of the MIPMs. Therefore, these MIPMs possess great potential for environmental monitoring and purification applications towards clenbuterol and other threatening species posing a severe threat to human life safety.