Core−Shell Molecularly Imprinted Polymers on Magnetic Yeast for the Removal of Sulfamethoxazole from Water

In this work, magnetic yeast (MY) was produced through an in situ one-step method. Then, MY was used as the core and the antibiotic sulfamethoxazole (SMX) as the template to produce highly selective magnetic yeast-molecularly imprinted polymers (MY@MIPs). The physicochemical properties of MY@MIPs were assessed by Fourier-transform infrared spectroscopy (FT-IR), a vibrating sample magnetometer (VSM), X-ray diffraction (XRD), thermogravimetric analysis (TGA), specific surface area (SBET) determination, and scanning electron microscopy (SEM). Batch adsorption experiments were carried out to compare MY@MIPs with MY and MY@NIPs (magnetic yeast-molecularly imprinted polymers without template), with MY@MIPs showing a better performance in the removal of SMX from water. Adsorption of SMX onto MY@MIPs was described by the pseudo-second-order kinetic model and the Langmuir isotherm, with maximum adsorption capacities of 77 and 24 mg g−1 from ultrapure and wastewater, respectively. Furthermore, MY@MIPs displayed a highly selective adsorption toward SMX in the presence of other pharmaceuticals, namely diclofenac (DCF) and carbamazepine (CBZ). Finally, regeneration experiments showed that SMX adsorption decreased 21 and 34% after the first and second regeneration cycles, respectively. This work demonstrates that MY@MIPs are promising sorbent materials for the selective removal of SMX from wastewater.


Introduction
Antibiotics are intensively used as human and veterinary medicines for the treatment and prevention of infectious diseases [1]. Among them, sulfamethoxazole (SMX) is a sulfonamide bacteriostatic antibiotic that has been commonly used during the last 80 years to treat urinary tract infections due to its low cost and broad spectrum of activity to treat bacterial diseases [2,3]. However, the widespread and indiscriminate use of SMX, as of other antibiotics, constitutes a huge potential threat to human health and contaminates natural ecosystems by affecting aquatic and soil organisms [4,5].
Yeasts, which belong to the fungus kingdom, are relatively large eukaryotic and single-celled microorganisms (diameters typically measuring 2.0-4.0 µm). Their cell wall includes glucan, mannan, chitin protein, and a small amount of lipids, and it has many surface chemical groups such as carboxyl (-COOH), carbonyl (-C=O), amino (-NH 2 ), hydroxyl (-OH), and phosphoryl (-P=O) groups. Moreover, yeast is very cheap, easy to obtain, and environmentally friendly. These advantages make yeasts appropriate and widely used as supports for bio-nanocomposites [31].
In the above-described context, the objectives of this study were to: (1) Prepare a bio-nanocomposite of yeast-Fe 3 O 4 (magnetic yeast, MY) using an in situ one-step preparation of nano-Fe 3 O 4 ; (2) use MY as the core to synthesize magnetic yeast-molecularly imprinted polymers (MY@MIPs) by a surface-imprinted polymerization method with MIPs as the shell and SMX as the template molecule; (3) characterize the resulting materials by Fourier-transform infrared spectroscopy (FT-IR), a vibrating sample magnetometer (VSM), X-ray diffraction (XRD), thermogravimetric analysis (TGA), specific surface area (S BET ) determination, and scanning electron microscopy (SEM); (4) test the removal performance of MY@MIPs toward SMX and compare it with those of MY and MY@NIPs (magnetic molecularly imprinted polymers without template); and (5) explore the selective sorption capacity of MY@MIP in a real complex matrix (wastewater collected at a STP) and in the presence of other pharmaceuticals (diclofenac and carbamazepine).

Preparation of Magnetic Yeast (MY)
Nano-Fe 3 O 4 was loaded onto the yeast cell surface by a one-step method as described by Tian et al. [32]. Briefly, the yeast cells were cultured in ultrapure water with glucose. After reaching the exponential growth phase (6-10 h), the yeast cells were collected by centrifugation (4000 rpm). Then, collected cells (1.0 g) were suspended in 40 mL of 0.125 M FeCl 3 solution in a three-necked flask and stirred for 1 h at room temperature. After that, 0.6 g of FeCl 2 ·4H 2 O was added under nitrogen atmosphere and stirred for another 1 h. The mixture was then heated in a water bath at 80 • C for 15 min, and the pH was adjusted to approximately 11 with 25% (w/v) ammonium hydroxide. Stirring was kept for 30 min and then stopped to age for 1 h. The resulting magnetic yeast (MY) was then washed, separated by applying a magnetic field, and then dried in an oven (35 • C, 4 h).

Preparation of Magnetic Yeast-Based Molecularly Imprinted Polymer (MY@MIPs)
MY was treated as the core and the MIPs as the shell. The process used for the production of MY@MIPs was as follows: 1 mg of SMX (template molecule) and 4 mmol of 2-vpy (monomer) were dissolved in 60 mL of ACN/toluene (3/1; v/v). This solution was then self-polymerized for 8 h at room temperature (25 • C). Subsequently, 100 mg of MY (polymer supporter), 0.36 mmol of AIBN (initiator), and EGDMA (crosslinker) were added into the polymerized solution (template:monomer:crosslinker, 1:4:20), which was ultrasonicated for 10 min. The mixture was heated and maintained at 60 • C for 24 h under stirring with nitrogen protection. At last, the MY@MIPs were washed with methanol/formic acid (9/1; v/v) for 12 h and purified for 24 h by a Soxhlet extraction method (the extraction solution was methanol). Meanwhile, the MY@NIPs were also produced by following the above-described procedure but in the absence of the template.

Characterization of MY, MY@MIPs, and MY@NIPs
Fourier-transform infrared spectra of the produced materials were obtained in a Shimadzu-IRaffinity-1 equipment, using an ATR module (FTIR-ATR), under a nitrogen purge. The measurements were recorded in the range 500-4000 cm −1 , 4.0 of resolution, 256 scans, and applying atmosphere and background correction.
A vibrating sample magnetometer (VSM EV9) with an oscillatory applied magnetic field (H) to a maximum of 22 kOe was used to determine the saturation magnetization (M S ). The M S was calculated by plotting the magnetic moment versus the applied magnetic field, and it corresponded to the plateau value of the magnetic moment reached divided by the sample mass (10 mg). The sample was encapsulated in an acrylic cylindrical container (5.85 mm of diameter and 2.60 mm of height), which was coupled to the lineal motor of the VSM EV9 instrument, centered between the two polar heads of the electromagnet used to fluctuate the magnetic field. The instrument was calibrated with a disk of pure nickel (8 mm of diameter) using a procedure that establishes the determination of the magnetic field, applied at around 1 Oe, while the dispersion of the magnetic moment is inferior to 0.5%. X-ray diffraction (XRD, 5-90 • ) was measured on a D8-Focus X-ray diffractometer (Bruker Optics) with a test rate of 10 • ·min -1 . The results were analyzed by Jade program (9.0) and Origin (9.0). Thermogravimetric analysis (TGA) was performed in a thermogravimetric balance Setsys Evolution 1750, Setaram, TGA mode (S type sensor). The samples were heated at a heating rate of 10 • C min −1 , under nitrogen atmosphere, from room temperature to 105 • C and from 105 • C to 900 • C, maintaining constant temperature until total stabilization of the sample mass at the end of both stages (approximately 30 min).
The S BET and micropore volume (W 0 ) were determined by nitrogen adsorption isotherms, acquired at 77 K using a Micromeritics Instrument, Gemini VII 2380, after outgassing the materials overnight at 120 • C. S BET was calculated from the Brunauer-Emmett-Teller equation in the relative pressure range 0.01-0.1. Pore volume (V p ) was estimated from the amount of nitrogen adsorbed at a relative pressure of 0.99.
The surface morphology of the materials was analyzed by scanning electron microscopy (SEM) using a Hitachi S4100. The images were obtained at magnifications of 500, 3000, and 10,000×.

Adsorptive Removal of SMX by the Produced Materials
The produced materials (MY, MY@MIPs, MY@NIPs) were used as adsorbents for the removal of SMX under batch operation conditions. Summarizing, the materials were put in contact with a 5 mg L −1 SMX solution in polypropylene tubes, which were shaken in a head-over-head shaker (80 rpm) for a predetermined period of time at controlled temperature (32 • C). The corresponding adsorbent material was separated from the suspension liquid by an external magnetic field. At last, the concentration of SMX in the liquid phase was measured by micellar electrokinetic chromatography (MEKC), using a methodology adapted from Silva et al. (2019) [33]. The experiments were conducted in triplicate, and control experiments without adsorbent were run in parallel. The performance of the materials was evaluated by carrying out kinetic, equilibrium, pH, selectivity, and regeneration/reutilization studies, described in detail in the next subsections.

Kinetic Adsorption Studies in Ultrapure Water
In the kinetic study, tubes containing 250 mg of adsorbent material (MY, MY@NIPs, or MY@MIPs), together with 10 mL of a 5 mg L −1 SMX solution in ultrapure water, were incubated and shaken as described above. After shaking during defined periods of time (t, min), at intervals from 0 to 24 h, the materials were separated from the aqueous phase and the remaining SMX concentration in solution was measured by MEKC. At each time, the corresponding value of the adsorbed concentration (q t , mg·g −1 ) was determined as follows: where C t (mg L −1 ) is the residual SMX concentration at time t, C 0 is the initial SMX concentration (mg L −1 ), and C m is the adsorbent dosage (mg·L −1 ). When adsorption equilibrium was attained, the percentage of adsorption R (%) was determined as: where C e (mg·L −1 ) is the residual SMX concentration at equilibrium.

Equilibrium Adsorption Studies in Ultrapure Water
For the equilibrium studies, the corresponding adsorbent material (MY, MY@NIPs, or MY@MIPs), with doses ranging from 50 to 2000 mg L −1 , was added to 10 mL of a 5 mg L −1 solution of SMX in ultrapure water. Tubes with the mixtures were shaken for 16 h, which allowed equilibrium to be reached. The materials were recovered from the suspension by the application of a magnetic field and the residual concentration of SMX was determined by MEKC. Then, for the different doses of material, the adsorbed concentration at the equilibrium (q e , mg·g −1 ) was determined as follows: where C e (mg L −1 ) is the SMX concentration in the liquid phase at equilibrium.

Adsorptive Performance of MY@MIPs
From the results of the above-mentioned kinetic and equilibrium studies in ultra-pure water, the most efficient material for removal of SMX was MY@MIPs. Thus, in order to assess the practical application of this material, further studies were carried out on the adsorptive performance of MY@MIPs under different experimental conditions.

Kinetic and Equilibrium Adsorption Studies in STP Effluent
The kinetic and equilibrium procedures described in Section 2.4. were carried out using MY@MIPs for the adsorptive removal of SMX from a real matrix, namely the effluent from a STP. In this case, 5 mg L −1 solutions of SMX were prepared using a STP effluent instead of ultrapure water. The effluent was collected from an urban STP in Aveiro (Portugal) that is designed to serve 159,700 population equivalents. This STP consists of primary and biological treatment stages. For this work, water was collected at the outlet of the biological decanter, as this is the final treated effluent that is discharged from the STP into the aquatic environment. Immediately after collection, the effluent was filtered through 0.45 µm, 293 mm Supor ® membrane disk filters (Gelman Sciences) and stored at 4 • C until use, which occurred within a maximum of 15 days. The collected effluent had a pH of 7.99, conductivity of 3.03 mS cm −1 , and total organic carbon content of 21.5 mg L −1 .

pH Study
Adsorption studies on the effect of pH were carried out at 32 • C with the initial conditions of C 0 = 5 mg L −1 in ultrapure water and C m = 300 mg L −1 . Experiments were carried out at three different pHs, namely 4, 7, and 8 (pH was adjusted by adding HCl or NaOH, 1 M). After shaking during 16 h, MY@MIPs were separated from the liquid suspensions, the residual concentration of SMX was analyzed by MEKC, and the corresponding q e (mg g −1 ) at each pH was determined using Equation (3).

Selective Adsorption
To study the selective capacity of MY@MIPs toward SMX, diclofenac (DCF) and carbamazepine (CBZ) were used as competing species in the adsorption experiments. These pharmaceuticals were selected due to their high global frequency of occurrence in wastewater, surface water, and groundwater and their recalcitrant properties, with low removal rates after conventional STP treatments [34]. The concentration of DCF and CBZ in ultrapure water solution was the same as that of SMX (5 mg L −1 ), the C m was 300 mg L −1 , the incubation temperature was 32 • C, the pH was 4, and shaking was maintained during 16 h. Then, the residual concentration of SMX at equilibrium was analyzed and the corresponding q e (mg g −1 ) was determined with Equation (3).

Regeneration and Reutilization
In order to evaluate the adsorptive performance after regeneration, after SMX saturation in ultrapure water, MY@MIPs were regenerated and then tested for the adsorption of SMX in four subsequent cycles. For the regeneration, saturated MY@MIPs were washed by methanol/acetic acid (9/1, v/v) through Soxhlet extraction during 72 h. Then, the regenerated material was used in adsorption experiments as described in previous sections (shaking during 16 h at 32 • C with the initial conditions of C 0 = 5 mg L −1 in ultrapure water and C m = 300 mg L −1 ). The residual concentration of SMX at the equilibrium was analyzed and the corresponding R (%) was determined as for Equation (2).

Preparation of MY
In this study, an in situ one-step method was carried out to load nano-Fe 3 O 4 particles on the surface of yeast, which was used as a biological solid support. Under alkaline conditions, Fe 2+ and Fe 3+ co-precipitated on the surface of yeast and then Fe(OH) 2 or Fe(OH) 3 was converted to nano-Fe 3 O 4 at 80 • C, according to the following chemical reactions: The observation of MY by a high-power optical microscope (Olympus CX22, Japan) clearly showed the loading of magnetic nanoparticles over yeast cells, as shown in Figure 1a. Meanwhile, Figure 1b represents the picture of MY at the actual size. Compared to other methods used to anchor Fe 3 O 4 nanoparticles on the surface of yeast biomass, such as cross-linking or electrostatic-interaction-driven hypercoagulation, the one-step method applied here only took approximately 3.5 h, in opposition to the referred methods, which can take up to 13.5 and 6 h, respectively [35], without considering the time of washing and drying. Hence, the results suggest that the one-step method is an interesting synthesis option.

Preparation of MY@MIPs
MY was used as support of a MIP-based material for the selective adsorption of SMX. Compared to other support materials such as SiO2, carbon nanotubes, or Fe3O4-SiO2 used in the literature [36,37], MY particles can act as support without the need of an intermediate chemical modification step and can distinctly improve grafting efficiency.
The shell of MIPs was co-polymerized on the surface of MY. Hence, to synthetize MIPs with affinity, selectivity, and appreciable removal capacity toward the target compound, the monomer and crosslinker types and the ratio of the reagents should be taken into account. Normally, if the template molecule has an alkaline chemical group, the monomer should be methacrylate (MAA), but if it has an acidic group, the monomer should be vinyl pyrimidine (vpy) [38]. As SMX has an oxazole moiety that displays acidity, 2-vpy was chosen as it has both a hydrogen-bond acceptor (N atom of pyridine) and alkalinity [39]. In this work, the molar ratio of the mixture of template and monomer was 1:4, as the monomer and template were in dynamic equilibrium, and it is not useful to add the monomer indiscriminately. Indeed, an excessive monomer may increase the non-selective sites, resulting in a selectivity decrease. On the other hand, during the synthesis of MIPs, in order to immobilize the template into the polymer without changing the spatial configuration of pores in the polymer, this must have a high rigidity. Therefore, it was necessary to use a crosslinker for increasing rigidity, EDGMA being selected due to its appropriate cost and solubility. However, if the ratio of monomer to crosslinker is too high, it will make the extraction of the template difficult due to the excessive rigidity of the MIP. Considering the referred considerations and conclusions from other studies [40,41], the ratio of monomer and crosslinker was selected to be 1:5.

Characterization of MY, MY@MIPs, and MY@NIPs
FTIR spectra of the produced materials (MY, MY@MIPs, and MY@NIPs), which were obtained in order to shed some light about the chemical groups present on their surface, are depicted in Figure  2. At 548 cm −1 , a characteristic adsorption peak belonging to the Fe-O chemical bond was observed for all materials. Compared to MY, MY@MIPs and MY@NIPs had some new peaks. Among them were the absorption bands at 2363 and 2328 cm −1 (MY@MIPs) and at 2377 and 2337 cm −1 (MY@NIPs), which were attributed to the stretching vibrations of -CN or -NC, respectively. The peak at 1758 cm −1 (MY@MIPs) or at 1727 cm −1 (MY@NIPs) belongs to the stretching vibration of C=O in the EGDMA ester group and the carboxyl group, suggesting that EGDMA worked on the surface. Moreover, MY@MIPs had new peaks at 1118 and 955 cm −1 , which belonged to the symmetrical and asymmetric stretching vibration of C-O in EGDMA, respectively, and reflected that it had a cross-linking polymerization on the surface of MY@MIPs. In the spectra of MY@MIPs and MY@NIPs, adsorption

Preparation of MY@MIPs
MY was used as support of a MIP-based material for the selective adsorption of SMX. Compared to other support materials such as SiO 2, carbon nanotubes, or Fe 3 O 4 -SiO 2 used in the literature [36,37], MY particles can act as support without the need of an intermediate chemical modification step and can distinctly improve grafting efficiency.
The shell of MIPs was co-polymerized on the surface of MY. Hence, to synthetize MIPs with affinity, selectivity, and appreciable removal capacity toward the target compound, the monomer and crosslinker types and the ratio of the reagents should be taken into account. Normally, if the template molecule has an alkaline chemical group, the monomer should be methacrylate (MAA), but if it has an acidic group, the monomer should be vinyl pyrimidine (vpy) [38]. As SMX has an oxazole moiety that displays acidity, 2-vpy was chosen as it has both a hydrogen-bond acceptor (N atom of pyridine) and alkalinity [39]. In this work, the molar ratio of the mixture of template and monomer was 1:4, as the monomer and template were in dynamic equilibrium, and it is not useful to add the monomer indiscriminately. Indeed, an excessive monomer may increase the non-selective sites, resulting in a selectivity decrease. On the other hand, during the synthesis of MIPs, in order to immobilize the template into the polymer without changing the spatial configuration of pores in the polymer, this must have a high rigidity. Therefore, it was necessary to use a crosslinker for increasing rigidity, EDGMA being selected due to its appropriate cost and solubility. However, if the ratio of monomer to crosslinker is too high, it will make the extraction of the template difficult due to the excessive rigidity of the MIP. Considering the referred considerations and conclusions from other studies [40,41], the ratio of monomer and crosslinker was selected to be 1:5.  The magnetic properties of the MY, MY@MIPs, and MY@NIPs were studied by VSM at room temperature, the MS of each material being shown in Table 1. The MS values were determined to be between 26 and 34 emu g −1 , which were compatible with good magnetization. Indeed, Figure S1, within Supplementary Information, shows that MY@MIPs can be easily separated by an external magnetic field, which is beneficial for the after-use separation of the saturated MY@MIPs from treated water, achieving one of the major goals of this study. The XRD spectrum of MY in the 2θ range of 20 to 80° is shown in Figure 3, where the (220), (311), (400), (422), (511), and (440) planes of Fe3O4 may be observed at 2θ = 30.22°, 35.40°, 43.36°, 53.68°, 57.21°, and 62.43°. This pattern is consistent with the standard XRD data of Fe3O4 in the JCPDS-International Centre for Diffraction Data (JCPDS Card: PDF#75-0033). Therefore, XRD results evidenced that Fe3O4 was successfully loaded onto the yeast surface during the production of MY and that, subsequently, surface molecular imprinting did not change the crystalline structure of magnetic nanoparticles. Similar patterns confirming the effective loading of magnetite have been reported in the literature on magnetic MIPs (MMIPs), including MMIPs produced for melamine analysis in milk [42], PEGylated magnetic core−shell structure-molecularly imprinted polymers (PMMIPs) for the specific adsorption of bovine serum albumin (BSA) [43], or core−shell MMIPs for the selective adsorption of tetracycline [44]. The magnetic properties of the MY, MY@MIPs, and MY@NIPs were studied by VSM at room temperature, the M S of each material being shown in Table 1. The M S values were determined to be between 26 and 34 emu g −1 , which were compatible with good magnetization. Indeed, Figure S1, within Supplementary Information, shows that MY@MIPs can be easily separated by an external magnetic field, which is beneficial for the after-use separation of the saturated MY@MIPs from treated water, achieving one of the major goals of this study.   : PDF#75-0033). Therefore, XRD results evidenced that Fe 3 O 4 was successfully loaded onto the yeast surface during the production of MY and that, subsequently, surface molecular imprinting did not change the crystalline structure of magnetic nanoparticles. Similar patterns confirming the effective loading of magnetite have been reported in the literature on magnetic MIPs (MMIPs), including MMIPs produced for melamine analysis in milk [42], PEGylated magnetic core−shell structure-molecularly imprinted polymers (PMMIPs) for the specific adsorption of bovine serum albumin (BSA) [43], or core−shell MMIPs for the selective adsorption of tetracycline [44]. Polymers 2020, 12, x FOR PEER REVIEW 9 of 21 The thermogravimetric (TG) and derivative thermogravimetric (DTG) curves of MY, MY@NIPs, and MY@MIPs are shown in Figure 4. All the materials evidenced three main weight loss peaks: The first at ~100 °C related to moisture; the second at ~300 °C related to the most thermolabile organic fraction; and the third centered at ~700 °C related to less thermolabile organic or inorganic fractions. For MY, a weight loss of approximately 62% was reached at 900 °C ( Figure 4a); the second weight loss peak is particularly accentuated in this material, as it is the one with the highest amount of yeast per unit mass of material and, thus, yeast cells carbonized with increasing temperature. Meanwhile, MY@NIPs ( Figure 4b) and MY@MIPs (Figure 4c) suffered, globally, a lower weight loss than MY, reaching 30% and 50% of weight loss, respectively, at 900 °C. This might be due to the introduction of less thermolabile structures in the composition of these materials (such as the magnetic nanoparticles and polymers). The results of SBET are shown in Table 1. The SBET of each material was as follows: MY-38.8 m 2 g −1 , MY@MIPs-43.2 m 2 g −1 , and MY@NIPs-39.2 m 2 g −1 . Similar SBET (47 m 2 g −1 ) were determined for magnetic sorbents with a metal-organic framework core and MIP shell [45]. Meanwhile, lower SBET, between 6 and 11 m 2 g −1 , have been measured for magnetic sorbents based on the iron oxide (Fe3O4) core and MIP shell [46,47]. Regarding the average pore diameter (D), it was 5.71, 5.41, and 5.06 nm respectively for MY, MY@MIPs, and MY@NIPs. Therefore, the three produced materials are mesoporous with no significant differences between them in terms of porosity. The thermogravimetric (TG) and derivative thermogravimetric (DTG) curves of MY, MY@NIPs, and MY@MIPs are shown in Figure 4. All the materials evidenced three main weight loss peaks: The first at~100 • C related to moisture; the second at~300 • C related to the most thermolabile organic fraction; and the third centered at~700 • C related to less thermolabile organic or inorganic fractions. For MY, a weight loss of approximately 62% was reached at 900 • C (Figure 4a); the second weight loss peak is particularly accentuated in this material, as it is the one with the highest amount of yeast per unit mass of material and, thus, yeast cells carbonized with increasing temperature. Meanwhile, MY@NIPs ( Figure 4b) and MY@MIPs (Figure 4c) suffered, globally, a lower weight loss than MY, reaching 30% and 50% of weight loss, respectively, at 900 • C. This might be due to the introduction of less thermolabile structures in the composition of these materials (such as the magnetic nanoparticles and polymers). The thermogravimetric (TG) and derivative thermogravimetric (DTG) curves of MY, MY@NIPs, and MY@MIPs are shown in Figure 4. All the materials evidenced three main weight loss peaks: The first at ~100 °C related to moisture; the second at ~300 °C related to the most thermolabile organic fraction; and the third centered at ~700 °C related to less thermolabile organic or inorganic fractions. For MY, a weight loss of approximately 62% was reached at 900 °C ( Figure 4a); the second weight loss peak is particularly accentuated in this material, as it is the one with the highest amount of yeast per unit mass of material and, thus, yeast cells carbonized with increasing temperature. Meanwhile, MY@NIPs ( Figure 4b) and MY@MIPs (Figure 4c) suffered, globally, a lower weight loss than MY, reaching 30% and 50% of weight loss, respectively, at 900 °C. This might be due to the introduction of less thermolabile structures in the composition of these materials (such as the magnetic nanoparticles and polymers). The results of SBET are shown in Table 1. The SBET of each material was as follows: MY-38.8 m 2 g −1 , MY@MIPs-43.2 m 2 g −1 , and MY@NIPs-39.2 m 2 g −1 . Similar SBET (47 m 2 g −1 ) were determined for magnetic sorbents with a metal-organic framework core and MIP shell [45]. Meanwhile, lower SBET, between 6 and 11 m 2 g −1 , have been measured for magnetic sorbents based on the iron oxide (Fe3O4) core and MIP shell [46,47]. Regarding the average pore diameter (D), it was 5.71, 5.41, and 5.06 nm respectively for MY, MY@MIPs, and MY@NIPs. Therefore, the three produced materials are mesoporous with no significant differences between them in terms of porosity. The results of S BET are shown in Table 1. The S BET of each material was as follows: MY-38.8 m 2 g −1 , MY@MIPs-43.2 m 2 g −1 , and MY@NIPs-39.2 m 2 g −1 . Similar S BET (47 m 2 g −1 ) were determined for magnetic sorbents with a metal-organic framework core and MIP shell [45]. Meanwhile, lower S BET , between 6 and 11 m 2 g −1 , have been measured for magnetic sorbents based on the iron oxide (Fe 3 O 4 ) core and MIP shell [46,47]. Regarding the average pore diameter (D), it was 5.71, 5.41, and 5.06 nm respectively for MY, MY@MIPs, and MY@NIPs. Therefore, the three produced materials are mesoporous with no significant differences between them in terms of porosity.

Characterization of MY, MY@MIPs, and MY@NIPs
The surface of MY, MY@MIPs, and MY@NIPs was examined by SEM ( Figure 5). All the figures suggested that the particles (either MY, MY@NIPs, or MY@MIPs) were elliptical, which is due to the use of yeast as support, as it has been observed to have an ellipsoid shape with uniform size [48]. As it may be seen, MY@NIPs and MY@MIPs have a dispersed and comparatively smoother appearance than MY, which is rough-faced due to the magnetic nanoparticles coating the smooth-faced yeast [49]. Moreover, under 3000× magnification, results showed that MY@MIPs had a better dispersion compared to MY and MY@NIPs. Under 10,000×, MY@MIPs showed a bigger porosity than the other materials, which may benefit the adsorption of SMX and improve the mass transfer rate from the aqueous phase.
Globally, the characterization results demonstrated the successful loading of Fe 3 O 4 on the yeast surface and the preparation of MIPs on the surface of MY.
Polymers 2020, 12, x FOR PEER REVIEW 10 of 21 The surface of MY, MY@MIPs, and MY@NIPs was examined by SEM ( Figure 5). All the figures suggested that the particles (either MY, MY@NIPs, or MY@MIPs) were elliptical, which is due to the use of yeast as support, as it has been observed to have an ellipsoid shape with uniform size [48]. As it may be seen, MY@NIPs and MY@MIPs have a dispersed and comparatively smoother appearance than MY, which is rough-faced due to the magnetic nanoparticles coating the smooth-faced yeast [49]. Moreover, under 3000× magnification, results showed that MY@MIPs had a better dispersion compared to MY and MY@NIPs. Under 10,000×, MY@MIPs showed a bigger porosity than the other materials, which may benefit the adsorption of SMX and improve the mass transfer rate from the aqueous phase.
Globally, the characterization results demonstrated the successful loading of Fe3O4 on the yeast surface and the preparation of MIPs on the surface of MY.

Adsorption Kinetics
The kinetic results on the adsorption of SMX onto the produced materials are shown in Figure  6, which evidences that, in all cases, the adsorbed concentration qt (mg g −1 ) rapidly increased until 360 min of contact and then slowly increased until becoming stable. Moreover, all the materials performed quite similarly from a kinetic point of view.
Comparing the results obtained here to those reported in the literature, it may be said that a shorter equilibrium time (around 20 min, at room temperature) was determined for the adsorption of SMX onto core−shell MIPs on the surface of magnetic carbon nanotubes (MCNTs@MIP) by Zhao et al. [25]. Meanwhile, using MIPs on the surface of yeast (yeast@MIPs), Wang et al. [16] found that

Adsorption Kinetics
The kinetic results on the adsorption of SMX onto the produced materials are shown in Figure 6, which evidences that, in all cases, the adsorbed concentration q t (mg g −1 ) rapidly increased until 360 min of contact and then slowly increased until becoming stable. Moreover, all the materials performed quite similarly from a kinetic point of view.
Comparing the results obtained here to those reported in the literature, it may be said that a shorter equilibrium time (around 20 min, at room temperature) was determined for the adsorption of SMX onto core−shell MIPs on the surface of magnetic carbon nanotubes (MCNTs@MIP) by Zhao et al. [25]. Meanwhile, using MIPs on the surface of yeast (yeast@MIPs), Wang et al. [16] found that (at 298 to 318 K) 200 min were necessary to attain equilibrium for the adsorption of ciprofloxacin (CIP), and Pan et al. [50] observed an equilibrium time around 375 min for the adsorption (at 303 K) of cephalexin. In any case, it has been noticed that surface-imprinting improves the binding kinetics as compared to traditionally imprinted materials, which take longer (usually around 12-24 h) to attain adsorption equilibrium [51].
Polymers 2020, 12, x FOR PEER REVIEW 11 of 21 (at 298 to 318 K) 200 min were necessary to attain equilibrium for the adsorption of ciprofloxacin (CIP), and Pan et al. [50] observed an equilibrium time around 375 min for the adsorption (at 303 K) of cephalexin. In any case, it has been noticed that surface-imprinting improves the binding kinetics as compared to traditionally imprinted materials, which take longer (usually around 12-24 h) to attain adsorption equilibrium [51]. Pseudo-first-order [52] and pseudo-second-order [53] kinetic models were applied to describe the adsorption kinetics of SMX onto the produced materials. The formulation of the models is as follows: Pseudo-first-order Pseudo-second-order where k1 (min −1 ) and k2 (g mg −1 min −1 ) are the pseudo-first-order and the pseudo-second-order rate constants.
The non-linear fitting kinetic parameters are summarized in Table 2. According to the correlation coefficient (R 2 ) and concordance between experimental and fitted qe values, both models described Pseudo-first-order [52] and pseudo-second-order [53] kinetic models were applied to describe the adsorption kinetics of SMX onto the produced materials. The formulation of the models is as follows: Pseudo-first-order Pseudo-second-order where k 1 (min −1 ) and k 2 (g mg −1 min −1 ) are the pseudo-first-order and the pseudo-second-order rate constants. The non-linear fitting kinetic parameters are summarized in Table 2. According to the correlation coefficient (R 2 ) and concordance between experimental and fitted q e values, both models described the SMX adsorption onto the produced materials, with the pseudo-second-order model describing slightly better the results onto MY and MY@NIPs and the pseudo-first-order model onto MY@MIPs.

Adsorption Isotherm
Equilibrium results on the adsorption of SMX onto the produced materials are shown in Figure 7. With the aim of describing these results, four isotherm models were used: Langmuir [54] and Freundlich [55] isotherm models for the adsorption of SMX onto MY@MIPs; BET isotherm [56] for the adsorption onto MY@NIPs; and Zhu−Gu isotherm [57] for the adsorption onto MY. The equations of these models are as follows: Langmuir isotherm Freundlich isotherm BET isotherm Zhu−Gu isotherm where q m is the maximum adsorption capacity (mg g −1 ); b (L mg −1 ) is the Langmuir equilibrium constant; k f (mg g −1 (mg L −1 ) −1/n ) is the Freundlich constant; n is the degree of non-linearity in the Freundlich isotherm; c is the BET constant, related to the energy of adsorption in the first adsorbed layer; g is the Zhu−Gu constant related to the first adsorption step (the first layer of molecules on the materials); e is the Zhu−Gu constant related to the subsequent layers adsorbed; and r is the aggregation number in the Zhu−Gu isotherm. Experimental results on the equilibrium of SMX adsorption onto the produced materials are depicted in Figure 7 together with fittings to the above-mentioned isotherm models. From Figure 7, it is evident that, contrarily to the adsorption onto MY@MIPs, in the case of MY and MY@NIPs, the q e did not tend to stabilization. Furthermore, in the C e range between 0 and 3 mg L −1 , a lower q e occurred for MY@NIPs than for MY. This may be related to the presence of chemical groups on the surface of MY, which were able to bind SMX groups, but became inaccessible in MY@NIPs due to molecular imprinting. In addition, a first stage with stabilization of q e at C e around 3 mg L −1 may be observed in the MY isotherm, which could be associated with the saturation of the chemical adsorption sites. In the case of MY@MIPs, the isotherm showed an increase in q e with C e with a stabilization trend from C e~3 mg L −1 . Furthermore, it should be noted that, at relatively low C e , the q e values determined for MY@MIPs are higher than those for MY and MY@NIPs, which points to their larger affinity for SMX. Polymers 2020, 12, x FOR PEER REVIEW 13 of 21 The fitted equilibrium parameters are shown in Table 3 together with the correlation coefficients of the fittings (R 2 ). In the case of SMX adsorption onto MY@MIPs, the Langmuir isotherm provided the best fitting of equilibrium results with a higher R 2 than the Freundlich model. From the results of non-linear fittings for the different isotherm models, the values of qm for MY, MY@NIPs, and MY@MIPs were, respectively, 23 ± 1, 3.8 ± 0.3, and 77 ± 3 mg g −1 . These values indicate that molecularly imprinted polymers with the template resulted in a substantial increase in the monolayer adsorption capacity, SMX adsorbing onto the surface of MY@MIPs in a homogeneous distribution by occupying specific sites. Similarly, equilibrium results on the adsorption of SMX onto the MCNTs@MIP produced by Zhao et al. [25] also fitted the Langmuir isotherm, but with a considerably lower qm (0.87 mg g −1 ). Indeed, compared to other materials used for the adsorption of SMX (Table 4), MY@MIPs are competitive in terms of SMX adsorption capacity.  The fitted equilibrium parameters are shown in Table 3 together with the correlation coefficients of the fittings (R 2 ). In the case of SMX adsorption onto MY@MIPs, the Langmuir isotherm provided the best fitting of equilibrium results with a higher R 2 than the Freundlich model. From the results of non-linear fittings for the different isotherm models, the values of q m for MY, MY@NIPs, and MY@MIPs were, respectively, 23 ± 1, 3.8 ± 0.3, and 77 ± 3 mg g −1 . These values indicate that molecularly imprinted polymers with the template resulted in a substantial increase in the monolayer adsorption capacity, SMX adsorbing onto the surface of MY@MIPs in a homogeneous distribution by occupying specific sites. Similarly, equilibrium results on the adsorption of SMX onto the MCNTs@MIP produced by Zhao et al. [25] also fitted the Langmuir isotherm, but with a considerably lower q m (0.87 mg g −1 ). Indeed, compared to other materials used for the adsorption of SMX (Table 4), MY@MIPs are competitive in terms of SMX adsorption capacity.

Kinetic and Equilibrium Adsorption Studies from STP Effluent
In order to assess the practical applicability of MY@MIPs, kinetic and equilibrium experiments were carried out in a real matrix, namely the effluent from a STP. The obtained results together with fittings to the considered kinetic and equilibrium models are in shown in Figure 8, and the fitted parameters are depicted in Table 5. As it may be seen, the pseudo-second-order and the Langmuir isotherm models were those that best described the kinetic and equilibrium experimental results, respectively. On the other hand, it is evident in Figure 8 that, under identical experimental conditions, the adsorption velocity was slower and the q e values were lower for the STP effluent than they were for ultrapure water. This was confirmed by the parameters in Table 5, especially by the comparatively lower q m (24 ± 2 mg g −1 ) than in ultrapure water, which might be related to interferences due to the complex composition of the STP effluent.
isotherm models were those that best described the kinetic and equilibrium experimental results, respectively. On the other hand, it is evident in Figure 8 that, under identical experimental conditions, the adsorption velocity was slower and the qe values were lower for the STP effluent than they were for ultrapure water. This was confirmed by the parameters in Table 5, especially by the comparatively lower qm (24 ± 2 mg g −1 ) than in ultrapure water, which might be related to interferences due to the complex composition of the STP effluent.

pH Study
Results from the study of pH effects on the adsorption of SMX onto MY@MIPs are shown in Figure 9. Under identical experimental conditions, except for the pH, decreasing q e values were obtained at pH 4 > pH 7 > pH 9, thus indicating that SMX adsorption onto MY@MIPs was favored under acidic conditions. This may be related to the pH influence on the status of not only the adsorbate (by protonation/deprotonation) but also the adsorbent (by surface charge). For SMX, the pK a values are 1.97 and 6.16 (Table S1, as Supplementary Information), which means that SMX is mostly positively charged (protonated NH 2 groups, NH 3 + groups) at pH < 1.97 but predominantly negatively charged (deprotonated NH groups, N−) at pH > 6. 16. Therefore, at the experimental pH 4, adsorption of SMX in the non-ionic form was favored, while at pH 7 and 9, SMX was mostly present in the anionic form, which partially hindered its adsorption. Indeed, the decrease in SMX adsorption from wastewater (with pH > 7) has already been related to electrostatic repulsion between the negatively charged SMX and the negatively charged surface of the waste-based adsorbents [59]. Moreover, the monomer used in the synthesis of MIPs was 2-vpa and the pK a of pyridine was 5.21, which is, therefore, negatively charged when pH > 5. 21. Thus, at the experimental pH 7, electrostatic repulsion forces between SMX and MY@MIPs cannot be disregarded, these increasing at pH 9.
mostly positively charged (protonated NH2 groups, NH3 + groups) at pH < 1.97 but predominantly negatively charged (deprotonated NH groups, N−) at pH > 6. 16. Therefore, at the experimental pH 4, adsorption of SMX in the non-ionic form was favored, while at pH 7 and 9, SMX was mostly present in the anionic form, which partially hindered its adsorption. Indeed, the decrease in SMX adsorption from wastewater (with pH > 7) has already been related to electrostatic repulsion between the negatively charged SMX and the negatively charged surface of the waste-based adsorbents [59]. Moreover, the monomer used in the synthesis of MIPs was 2-vpa and the pKa of pyridine was 5.21, which is, therefore, negatively charged when pH > 5. 21. Thus, at the experimental pH 7, electrostatic repulsion forces between SMX and MY@MIPs cannot be disregarded, these increasing at pH 9. Figure 9. Effect of pH on the percentage of SMX adsorption (R (%)) onto MY@MIPs. Figure 9. Effect of pH on the percentage of SMX adsorption (R (%)) onto MY@MIPs.

Selective Adsorption
In order to find out the selectivity of MY@MIPs toward SMX, its adsorption was compared to those of DCF and CBZ from their single solution and then from their ternary solution. The values of percentage of adsorption (R (%)) for the single adsorption of each pharmaceutical are shown in Figure 10a and for adsorption from their ternary solution in Figure 10b.

Selective Adsorption
In order to find out the selectivity of MY@MIPs toward SMX, its adsorption was compared to those of DCF and CBZ from their single solution and then from their ternary solution. The values of percentage of adsorption (R (%)) for the single adsorption of each pharmaceutical are shown in Figure  10a and for adsorption from their ternary solution in Figure 10b. The results in Figure 10a evidence that MY@MIPs have a larger R (%) for SMX than for DCF or CBZ. Furthermore, under identical experimental conditions but from the ternary solution of the considered pharmaceuticals (Figure 10b), selective adsorption of SMX onto MY@MIPs occurred. Indeed, the R (%) determined for SMX from the ternary solution was just slightly lower than from its single solution, which points to the selectivity of MY@MIPs. Moreover, the results reflected that the adsorption of MY@MIPs was SMX > DCF > CBZ.
In this work, 2-vpy was the monomer and it was combined through -NH2. Moreover, considering the structure and properties of SMX, DCF, and CBZ, which are depicted in Table S1, they all have -NH2 and/or -NH groups. However, SMX has two amino groups: -NH and -NH2, which is probably the main reason for its selective adsorption onto MY@MIPs under the presence of DCF and CBZ. For DCF and CBZ, the pKa was 4 and 15.96 (Table S1), respectively. Meanwhile, the pKa value of 2-vpy (monomer) is 5.21, which may explain why MY@MIPs had a better removal ability for DCF The results in Figure 10a evidence that MY@MIPs have a larger R (%) for SMX than for DCF or CBZ. Furthermore, under identical experimental conditions but from the ternary solution of the considered pharmaceuticals (Figure 10b), selective adsorption of SMX onto MY@MIPs occurred. Indeed, the R (%) determined for SMX from the ternary solution was just slightly lower than from its single solution, which points to the selectivity of MY@MIPs. Moreover, the results reflected that the adsorption of MY@MIPs was SMX > DCF > CBZ.
In this work, 2-vpy was the monomer and it was combined through -NH 2 . Moreover, considering the structure and properties of SMX, DCF, and CBZ, which are depicted in Table S1, they all have -NH 2 and/or -NH groups. However, SMX has two amino groups: -NH and -NH 2 , which is probably the main reason for its selective adsorption onto MY@MIPs under the presence of DCF and CBZ. For DCF and CBZ, the pKa was 4 and 15.96 (Table S1), respectively. Meanwhile, the pK a value of 2-vpy (monomer) is 5.21, which may explain why MY@MIPs had a better removal ability for DCF than for CBZ. Selectivity toward SMX was also verified by Zhao et al. [25], who prepared MCNTs@MIP by using SMX as the template molecule and copolymerization of vinyl end groups on the surface of MCNTs [25]. These authors demonstrated the selective adsorption of SMX under the presence of other sulfonamides (SAs), namely sulfamethazine (SMZ), sulfamerazine (SMR), sulfadimethoxine (SDM), and sulfameter (SME). Still, the adsorbed concentration of SMX from the quinary solution was lower than from its single solution, which was ascribed to the close structure of the other SAs, which, therefore, could competitively occupy the imprinted sites.

Regeneration and Reutilization
Saturation of the produced MY@MIPs with SMX was carried out as described in Section 2.4.1. At that moment, R (%) calculated by Equation (2) was 92 ± 4%. Then, saturated MY@MIPs was regenerated as indicated in Section 2.5.4. and reused for the adsorption of SMX until saturation. A total of four regeneration/reutilization cycles were performed and the R (%) calculated for each of them are shown in Table 6. Table 6. SMX adsorption onto MY@MIPs in subsequent cycles after regeneration.

Cycles R (%)
SMX saturation 92 ± 4 1 73 ± 3 2 61 ± 2 3 58 ± 1 4 55 ± 2 As it may be seen, after cycles 1 and 2, the R (%) values decreased to 73 ± 3 and 61 ± 2%, respectively. Such decreases (21 and 34%, respectively, in cycle 1 and 2) indicate that the regeneration procedure affected the adsorption sites on the surface of MY@MIPs, which lost efficiency in the removal of SMX. After cycle 2, just a slight decrease in R (%) occurred, its value being similar in cycles 3 and 4 (58 ± 1 and 55 ± 2%, respectively). Therefore, deterioration of MY@MIPs adsorptive properties was not progressive with successive regenerations but occurred initially, with the performance remaining stable after cycle 2. MIPs sorbents are known to be easily regenerated by washing with organic solvents, with mixtures of methanol and acetic acid having been successfully employed to remove adsorbed pharmaceuticals [60]. Using the same regeneration agent as in this work, namely methanol/acetic acid (9/1, v/v), Dai et al. [61] regenerated MIPs synthesized for the adsorption of diclofenac and carried out thirty cycles with ≥95% recovery. Likewise, Duan et al. [62] also used this mixture for the regeneration of a multitemplate MIP, which was used in twenty regeneration/reutilization cycles, giving ≥95% removal of ibuprofen, naproxen, ketoprofen, diclofenac, and clofibric acid. Wang et al. [48], who used MIPs on the surface of yeast (yeast@MIPs), desorbed ciprofloxacin using the same mixture with losses of only about 8.5% of initial capacity after five cycles. Therefore, the relatively larger deterioration of the adsorptive performance observed in the present work may be related to the fact that magnetic yeast was used here as MIPs support. Thus, further work is to be carried out on the regeneration of the produced MY@MIPs, to maintain a high R (%) upon cyclic operation.

Conclusions
This work developed an efficient strategy to prepare yeast-Fe 3 O 4 (magnetic yeast, MY) and then used molecularly imprinted technology (MIT) to modify MY. The characterization of the produced magnetic yeast-molecularly imprinted polymers (MY@MIPs) showed that elliptical and monosized imprinted polymeric nanospheres with a surface area of about 43.2 m 2 g −1 were successfully produced. Sulfamethoxazole (SMX) adsorption studies using MY@MIPs indicated that the equilibrium was attained in 360 min either in ultrapure water or in a sewage treatment plant (STP) effluent. The Langmuir isotherm model provided the best fitting of equilibrium results and pointed to the monolayer and favorable adsorption of SMX onto MY@MIPs. In addition, the fitted parameters of the Langmuir isotherm model indicated that the maximum SMX adsorption capacity of MY@MIPs was 77 and 24 mg g −1 in ultrapure water and STP effluent, respectively. The pH study pointed out that hydrogen binding was underneath the SMX adsorption onto MY@MIPs. Moreover, MY@MIPs showed successful selective adsorption of SMX from ternary solution under competition by other pharmaceuticals, namely diclofenac (DCF) and carbamazepine (CBZ). Finally, regeneration implied a reduction in SMX removal by MY@MIPs in the first two cycles, then tending to stabilization. Overall, it may be concluded that the MIPs-coated magnetic yeast designed here could be an alternative adsorbent for the selective removal of SMX from complex matrices such as wastewaters.