Molecular Insights into the pH-Dependent Adsorption and Removal of Ionizable Antibiotic Oxytetracycline by Adsorbent Cyclodextrin Polymers

Effects of pH on adsorption and removal efficiency of ionizable organic compounds (IOCs) by environmental adsorbents are an area of debate, because of its dual mediation towards adsorbents and adsorbate. Here, we probe the pH-dependent adsorption of ionizable antibiotic oxytetracycline (comprising OTCH2 +, OTCH±, OTC−, and OTC2−) onto cyclodextrin polymers (CDPs) with the nature of molecular recognition and pH inertness. OTCH± commonly has high adsorption affinity, OTC− exhibits moderate affinity, and the other two species have negligible affinity. These species are evidenced to selectively interact with structural units (e.g., CD cavity, pore channel, and network) of the polymers and thus immobilized onto the adsorbents to different extents. The differences in adsorption affinity and mechanisms of the species account for the pH-dependent adsorption of OTC. The mathematical equations are derived from the multiple linear regression (MLR) analysis of quantitatively relating adsorption affinity of OTC at varying pH to adsorbent properties. A combination of the MLR analysis for OTC and molecular recognition of adsorption of the species illustrates the nature of the pH-dependent adsorption of OTC. Based on this finding, γ-HP-CDP is chosen to adsorb and remove OTC at pH 5.0 and 7.0, showing high removal efficiency and strong resistance to the interference of coexisting components.


Introduction
Ionizable organic compounds (IOCs) occupy a large fraction of the pre-registered REACH compounds [1]. They have one or more pK a values and are present in the form of a mixture of ionized and unionized species. The coexistence of these species is pH-dependent, which complicates the removal of IOCs in aquatic systems [1,2]. Particularly, veterinary antibiotics (e.g., oxytetracycline (Figure 1)), one emerging IOCs, are ubiquitously detected [3] and have aroused serious concerns about the spread of antibiotic resistance genes [4,5].
Adsorption is a widely used approach in the industry to remove organic pollutants from water [6,7]. This technique is effective for antibiotic removal [7]. Commonly, adsorption efficiency of antibiotics is controlled by solution pH. It is probably due to dual mediation of pH toward antibiotics and adsorbents. Specifically, the solution pH mediates the fraction of the species by altering ionization degree of antibiotics. The species of many antibiotics (e.g., tetracyclines [8,9], sulfathiazole [10] and sulfamethazine [11]) have been proven to differ in adsorption affinity. Meanwhile, conventional natural and synthetic adsorbents themselves (e.g., clays [8,12], humic substances [13,14], activated carbons [15], and CNTs [16]) probably suffer changes in structure and/or chemistry from the variation of pH. In turn, these changes may render alteration of adsorption sites of adsorbents. The dual mediation of pH complicates adsorption behavior of antibiotics, hereby making an impediment to the understanding of their pH-dependent adsorption.
Cyclodextrin polymers (CDPs) are growingly concerned environmental adsorbents, because they are endowed with a unique property of molecular recognition [17,18]. They are widely applied in adsorption and removal of various organic pollutants [19][20][21][22][23][24][25][26]. The polymers are composed of functional unit (i.e., CD monomer) and cross-linking agents (e.g., epichlorohydrin (EPI)). CD monomer with an axial open cavity of hydrophobic character ( Figure 1) is capable of including organics (or a moiety) in terms of geometric compatibility [27][28][29]. This inclusion is driven by various intermolecular interactions (e.g., hydrogen bonding, hydrophobic interactions, van der Waals forces, and electrostatic interactions) [30][31][32], indicating a unique property of molecular recognition. For example, the species of some IOCs with one or two pK a values are included by CD in different patterns [33,34]. This uniqueness is proved to remain in the polymers and considered as the primary adsorption mechanism [17][18][19]. Furthermore, the cross-linking agent EPI containing two reactive functional groups can form bonds with CD monomer and/or itself. The as-prepared polymer is a mixture of different CD units joined by repeating glyceryl linkers, representing a threedimensional polymer network [18,19,21]. As the linkages of CD moieties are relatively polar, the polymer is highly hydrophilic and readily swells in water. This character makes adsorption sites of the polymer (e.g., CD cavity, pore channel, and network) easily accessible [18,19,21]. Additionally, the polymers are stable in both acidic and alkaline solutions [21], showing poor susceptibility to pH variations. It is due to less susceptibility of CD (pK a = ca. 12) and isopropyl alcohol transformed by EPI (pK a = ca. 17.1) to ionization under general conditions.
With these advantages, CDPs are suitable as model adsorbent for illustrating adsorption behavior of IOCs at varying pH. Therefore, this study was aimed to probe the pH-dependent adsorption of oxytetracycline (OTC, as model compound) onto CDPs and thus to select CDPs for removing the antibiotic. To this end, we investigated inclusion complexation (i.e., binding constants and complexation ratios) of CD and OTC through a spectroscopic titration technique and address CD inclusion mechanisms for the species on the basis of structural characterization and molecular simulation of the complexes. We conducted batch adsorption experiments of OTC onto CDPs for measuring adsorption affinity of OTC over pH. Mathematical equations for relating adsorption affinity of OTC at varying pH (or the species) to adsorbent properties were derived to reveal the nature of pHdependent adsorption. The equations were used to select CDPs to remove OTC at different pH from water.

Materials
Oxytetracycline dihydrate (OTC, purity .95%) was obtained from AF Pharma LLC (US). Information on four cyclodextrin monomers (i.e., b-CD, RMCD, HPCD and c-CD) were listed in Table S1. b-Cyclodextrin was recrystallized twice in deionized water and dried under vacuum. Other cyclodextrins were used as received. Other chemicals used were of analytical reagent grade.
Adsorbent CDPs were synthesized in our previous work [19] where the ratio of cyclodextrin to cross-linking agent EPI was set at 1:5. They comprised four onefold polymers (i.e., b-CDP, c-CDP, HP-CDP and RM-CDP) with single CD as complex and three composite polymers (i.e., b-c-CDP, b-HP-CDP and c-HP-CDP) with an equimolar mixture of two CDs as double complexes. Their physiochemical properties (including CD content, crosslinking degree, swelling ratio, particle size, pore size, surface area, and pore volume) were characterized and compiled in Table 1 in reference 19 [19]. The parameter CD content of the as-prepared polymers varied from 38.62% to 71.91%, depending on type and number of the complex used. The corresponding cross-linking degree, representing the mean amounts of EPI per CD unit, was in the range of 19.82 to 4.57.

Inclusion complexation between CD and OTC at varying pH
The solutions of OTC and CD were adjusted to five pH values ranging from 4.17 to 9.30, by adding appropriate amounts of 1mol/L HCl or NaOH. The concentration of OTC was 20 mmol/ L. A series of concentrations of CD were set, e.g., 0 to 2000 mmol/ L for b-CD, 0 to 5200 mmol/L for RMCD and HPCD, and 0 to 6000 mmol/L for c-CD. The solutions were incubated at room temperature (ca. 15uC) and protected from light. Absorption spectra of the solutions were recorded by a UV-2300 spectrophotometer (Techcomp Limited, Shanghai, China). Peak wavelengths varied from 353 nm to 372 nm, depending on pH. Data were used to obtain inclusion constant (K C , L/mol) and inclusion ratio n (namely a ratio of CD to OTC) in terms of the extended Benesi-Hildebrand equation (Equation 1) [35]. The parameter n is commonly 1/2, 1 or 2, which should yield the highest values of fitting degree of the equation. As n is determined, the parameter K C can be calculated as a ratio of the intercept to the slope of the corresponding equation (Equation 1).
where A 0 and A are the absorbance of OTC and the mixture at peak wavelength, respectively; [I] 0 and [M] 0 are the initial concentrations of OTC and CD, respectively; De is the difference of the molar absorptivities of complexed and free OTC; L is optical distance; and n is inclusion ratio.
Moreover, inclusion complexes of CD with OTC were prepared by the grinding method [36]. Briefly, equivalent molar mixes of OTC and CD were grinded for 1 h with addition of three-fold deionized water. The complexes were lyophilized and characterized by FTIR (Shimadzu IRprestige-21, Japan) and NMR (Varian INOVA 400 MHz, USA). The complexes were prepared as KBr disks for FTIR measurement. All NMR spectra of the complexes were recorded in D 2 O with the chemical shift of HOD as reference.

Adsorption of OTC onto cyclodextrin polymers at varying pH
Batch adsorption experiments were conducted in 25-mL glass centrifuge tubes and at 25uC, with a solid-liquid ratio of 1:50. Concentrations of OTC varied from 1 to 30 mg/L. The solutions were adjusted to five pH values (varying from 4.65 to 10.11) by adding 1-mol/L HCl or NaOH. Appropriate amounts of 3-mol/L NaCl were added to all the solutions to obtain 0.02 mol/L Na + or Cl 2 , in order to minimize interference from changes of ionic strength upon the pH adjustment. Then the solutions were kept in the dark and shaken by a rotary shaker (180 rpm) for 30 min. The 30-min contact time was sufficient to ensure adsorption equilibrium on the basis of the adsorption kinetic experiments (Text S1). The samples were centrifuged at 4000 rpm for 15 min. A drop of 6-mol/L HCl was added to the supernatant (acidified to pH ca. 2.5), in order to prevent OTC decomposition (Text S2). The supernatant was filtered through a 0.45-mm Millipore membrane. OTC in the filter was measured on an L-2000 HPLC (Hitachi, Japan). The stationary phase was a Hypersil ODS C-18 column (4.66250 mm, 5 mm). The mobile phase comprised 22% acetonitrile and 78% NaH 2 PO 4 -H 3 PO 4 buffer solution (pH 2.30). The flow rate was set at 1.0 mL/min. The column oven was maintained at 25uC. The injection volume was 10 mL. The detection wavelength was set at 355 nm. OTC was eluted at about 5 min. Data was processed with Langmuir and Freundlich models to obtain adsorption constants (Text S1). And the adsorption distribution coefficients (K d , L/Kg) were also calculated (Text S1).
Additional adsorption experiments were carried out at three pH values (i.e., 5.0, 7.0 and 9.0) for all the polymers, to confirm loading of OTC at varying pH onto the adsorbents. The initial concentration of OTC was set as 20 mg/L. Once approaching adsorption equilibrium, the polymers were centrifuged and rinsed with 10-mL deionized water. Then, the polymers were lyophilized for FTIR characterization.

Computational simulation of inclusion complexation of CD and OTC
The atomic coordinates of b-CD (refcode POBRON) and c-CD (refcode CIWMIE10) were selected from Cambridge Structural Database (CSD). Starting geometries of RMCD and HPCD molecules were built on the basis of the structure of b-CD using Chem3D Ultra software (version 8.0.3). RMCD was an optional derivative of b-CD with methylation of all OH groups at C2, C3 and C6 positions ( Figure 1). HPCD was the hydropropylated derivative at C2 position ( Figure 1) with substitution degree of 6.4 (Table S1). Starting geometries of the species of OTC were also built using the same software. Energy minimization of all the molecules was separately performed using PM3 semi-empirical method of Gaussian 09 software [37]. Based on the optimized geometries, molecular docking between CD and the species of OTC was carried out on CDOCKER module of Discovery Studio (DS, version 2.5) with the CHARMm force field. Ten most stable poses are obtained by the molecular docking technology. Among them, the pose of CD-OTC in which inclusion sites and binding orientation are most consistent with those of the complex evidenced by FTIR and NMR characterization would be chosen for further optimization. The selected complex structures were put into Solvation module of DS with the CHARMM force field to simulate the real water environment. The obtained structure then was optimized using Minimization and Minimization (QM-MM) protocols of DS. For the resulting structures, intermolecular interactions involved in the inclusion were discriminated using the Multiwfn 2.4 package developed by Lu and Chen [38]. In addition, energies of van der Waals forces and electrostatic interactions were recorded from the docking results of DS with the consideration of solvent contributions.

Computational simulation of interactions of EPI oligomers with OTC
The polymer actually is a mixture of the materials comprising various ratios of CD to EPI, since the polymerization is rapid and relatively directionless. Due to the complexity of the cross-linking reaction, the structure of cyclodextrin polymers remains ambiguous. Here, we constructed some oligomers of EPI (e.g., ringopened, tricyclo-, pentacyclo-and heptacyclo-ones) that are essential units of network and pore channel of the polymers. These oligomers may be indicative of structural units of the polymers. And hence, they were used as structure mode to bind the species of OTC. Their geometries were optimized with the aforesaid method for revealing molecular interactions between the polymer and guest compound in a qualitative way.

Data analysis
All experiments were performed in triplicate and the resulting means were used in the following model analysis (Equation 2). Notably, the stability of OTC was ensured under all experimental circumstances (Text S2). The antibiotic comprises four species (i.e., OTCH 2 + , OTCH 6 , OTC 2 and OTC 22 ) with three pK a (i.e., 3.27, 7.32 and 9.11). Apparent interaction constants of OTC (i.e., CD inclusion and adsorption affinity) are considered as overall contribution of the four species (Equation 2). In Equation 2, the fractions of the species can be calculated by substituting pH and pK a into Equations S1-S4 (Text S3). To reduce system error and randomness in the pH designation, we used LINGO 9.0 optimization software (Lingo System Inc., USA) to obtain the interaction constants of the species. The model LINGO1 was chosen and all coefficients were constrained to be $0 (Text S3). The interaction constants were non-weighted. The application of the LINGO software would be tested by comparing the apparent interaction constants calculated on the basis of Equation 2 and obtained from the respective experiments.
where K is the apparent interaction constant of OTC; K + , K 6 , K 2 and K 22 are the interaction constants of OTCH 2 + , OTCH 6 , OTC 2 and OTC 22 , respectively; and a + , a 6 , a 2 and a 22 are corresponding mass fractions of the species.

Results and Discussion
Inclusion of CD with the species of OTC Addition of CD caused obvious changes in the UV-vis spectrum of OTC ( Figure S1), evidencing the formation of inclusion complex. Comparison of regression coefficients (R 2 ) derived from Equation 1 showed the complex had an inclusion ratio of 1:1 (Table S2). The corresponding inclusion constant (i.e., K C ) was compiled in Figure 2 and Table S2, relating to both CD type and solution pH. For instance, the inclusion capacity of b-CD with OTC, initially increased with pH increasing, approached the maximum at pH 4.75, and subsequently reduced. An inverse profile of inclusion capacity over pH, however, was observed for RMCD with the minimum at pH 6.36.
The K C value of the species (Table 1) was obtained through the LINGO optimization of apparent K C values for each CD at varying pH. A good linear relationship (R 2 $0.74 and P,0.05) was observed between the apparent K C derived from spectroscopic titration experiments and calculated on the basis of Equation 2 into which the K C and fraction of the species were substituted (Table S3). This relationship evidenced that the apparent K C of OTC at varying pH (lines in Figure 2) could be obtained by Equation 2. Comparison of the K C values (Table 1) demonstrated that molecular recognition of the species was related to cavity size [39] (e.g., b-CD versus c-CD) and structural modification (b-CD versus HPCD/RMCD) of CD (Figure 1). It is due to the fact that the species differing in geometry ( Figure 1) theoretically match the cavity of CD in different patterns. Overall, the differences in the K C values of the species account for the pH-dependent inclusion of OTC, taking into account the nature of the pH inertness of CD.

Inclusion mechanisms of CD with the species of OTC
Both FTIR and NMR characterization also convinced inclusion complexation of CD (e.g., b-CD) with the antibiotic. In the FTIR spectrum of the complex ( Figure S2), stretching vibration peaks of O-H (3488 cm 21 ) and dimethylamino groups (3372 cm 21 ) of OTC turned to a strong and sharp peak at 3370 cm 21 . In the NMR spectrum of the complex (Table S4), chemical shift of the protons at C4 and C6 of the OTC molecule transported from d2.686 and d1.601 to d2.698 and d1.635, respectively. And the concomitants were both downshield shift of H1-H6 protons of b-CD and reduction of splitting numbers of H3 peak (Table S4, Figure S3). The deshielding effect of b-CD is commonly attributable to the formation of hydrogen bonding. These two results clearly reveal that dimethylamino and its adjacent hydroxyl (i.e., C3, C5, C6 and/or C12) of the molecule are involved in the inclusion by forming hydrogen bonding.
The inclusion complex was visualized using molecular docking technique ( Figure S4). At least four interactions (i.e., hydrogen bonding, hydrophobic interactions, van der Waals forces, and electrostatic interactions) are identified to form inclusion complexes but to different extents. Specifically, hydrogen bonding was observed to form between hydroxyl group/ether bond at the wide rim of b-CD and C5/C12-hydroxyl groups of all the species with an exception of OTC 2 . It supported the deshielding effect of CD in the NMR characterization. The number of hydrogen bonding involved followed an order of OTCH 2 + (2).OTCH 6 (1) = OTC 22 (1).OTC 2 (0)) ( Figure S4), which was in accordance with the decreasing order of K C values (i.e., OTCH 2 + .-OTCH 6 .OTC 22 .OTC 2 ( Table 1)). It can be concluded that hydrogen bonding plays the predominant role in the inclusion complexation of CD with the antibiotic. It was worthy noting that OTCH 6 with less solubility and higher hydrophobicity [40] had stronger inclusion potential than OTC 22 , though both species had the same number of hydrogen bonding. It was ascribed to the presence of hydrophobic interactions that commonly are recognized as driving forces for the inclusion [31]. Furthermore, total interaction energy of the complex, denoted a sum of van der Waals forces and electrostatic interactions, followed an order of OTC 22 .OTC 2 .OTCH 2 + .OTCH 6 that is independent of the K C of the species (Table S5, Figure S5). This finding implied that neither van der Waals forces nor electrostatic interactions played significant roles in the inclusion.
Few attempts to investigate CD inclusion of IOCs with one-or two-pK a also demonstrate that the molecular recognition of the species involves various intermolecular interactions depending on CD type and compound structure [41][42][43][44]. Sebestyen and coworkers summarize that the differences in the K C of the species for 8 amino acids are a consequence of the coordination of hydrogen bonding, steric effects and electrostatic interactions [42]. For 3 aromatic amino acids, phenylalanine, tyrosine and tryptophan, the K C is in the range of 80-120 L/mol for the anion while it is much smaller for the zwitterion and negligible for the protonated species. However, the differences were less pronounced for the other five amino acids.

Adsorption of OTC onto CDPs
The kinetics of adsorption was fast and the plateau reached within 15-30 min ( Figure S6). The data of adsorption kinetics was well fitted to pseudo-second order model ( Figure S7) relative to pseudo-first order model (Table S6). The rate controlling step was recognized as intraparticle diffusion for RMCDP and c-HP-CDP on the basis of the moving boundary model (Text S1) while liquid film diffusion for the other five polymers (Table S7). It is mainly due to rather high solubility [40] and thus hydrophilic character of OTC, limiting mass transfer of the antibiotic to hydrophobic surfaces of materials. Figure S8 represented adsorption isotherms of OTC. The adsorption constants derived from Langmuir (K L ) and Freundlich (K F ) models were presented in Tables S8 and S9, respectively. In some cases, the sign of K L is negative (Table S8), indicating the poor application of Langmuir model. In contrast, Freundlich model provided good fit to all data with R 2 $0.92 (Table S9). The values of heterogeneity factor (n) of Freundlich model ranged from 0.59 to 2.04 (Table S9). Obviously, most of them deviated from unit (n = 1) at which the Freundlich isotherm translates to a linear relationship identical to the linear isotherm governed by hydrophobic partitioning. Such deviation demonstrated that the OTC adsorption process was governed by complicated interactions other than hydrophobic interactions. Even for the same adsorbent, the values of n also varied with solution pH, indicating that adsorption interactions should differ at varying pH. The resulting K F was of dimension, so that it would be a poor indicator for adsorption potential (Table S9). Instead of K F , the dimensionless K d , defined as the mean of K d at all OTC levels (Text S1), was used to represent adsorption affinity of OTC in this study. Comparison of the K d values indicated that this term varied with pH and the maxima commonly approached at pH 5.5-7.5 ( Figure 3). Here, the K d values were comparable to those reported for OTC onto other adsorbents, such as montmorillonite, kaolinite and goethite (ca. 769, 100 and 900 L/Kg at pH 6.5-6.7) [45], and native montmorillonite (ca. 1000 L/Kg at pH 5.0) [46].

Selective adsorption and immobilization of the species
The term K d of the species (Table 2) was obtained through the LINGO optimization of the apparent K d values of OTC at varying pH ( Figure 3). There was good similarity between the apparent K d of OTC calculated on the basis of Equation 2 and obtained from batch adsorption experiments (Table S10) The species-specific adsorption may be attributed to both ionization of OTC at varying pH and different adsorption interactions of the species, similar to the finding of the pHdependent interactions of OTC in clay and organic matter [46].

Recognition of interactions of the species of OTC onto CDPs
As it is well recognized, CD inclusion is the dominant adsorption interaction for CD-based polymers toward most guest compounds [17][18][19]38]. In most cases, high CD content corresponds to high adsorption affinity [17][18][19]38]. In contradiction to this finding, the K d of the species in this study ( Table 2) was not correlated with CD content of the polymers that is compiled in Table 1 in our previous work [19]. Especially, c-HP-CDP that has the second lowest CD content [19] exhibited the highest affinity to the species OTCH 2 + , OTCH 6 , and OTC 2 ( Table 2). Such deviation suggested that the adsorption of the species be governed by a combination of CD inclusion and other mechanisms.
The adsorption interactions were attempted to distinguish using multiple linear regression (MLR) analysis for the K d of OTCH 6 and OTCH 2 with CDPs properties. CD content, cross-linking degree, swelling ratio, particle size, pore size, surface area, pore volume, and species-specific inclusion constant were respectively set as variables x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , and x 8 . The MLR analysis was performed by the stepwise method in the SPSS Statistic program (SPSS 17.0 version) (Equations 3 and 4). Collinearity between variables was excluded in terms of variance inflation factor ,10. The analysis was repeated (Equations 4 and 5) when substituting variables x 1 and x 8 with variable x 9 (denoted x 1 6x 8 ). Variable x 9 , relating to both CD content and species-specific inclusion constant, may be indicative of effective inclusion.
OTC -: K d~0 :819x 2 (R 2~0 :670) ð4Þ OTCH + : K d~1 :075x 7 z0:453x 9 (R 2~0 :930) ð5Þ The adsorption of OTCH 6 was correlated with pore volume and effective inclusion capacity (Equations 3 and 5). It highlights the presence of pore-filling [48,49] and CD inclusion mechanisms for the species OTCH 6 . Comparison of standardized coefficients showed that pore-filling mechanism contributed more than CD inclusion one. It should be noted that the pore of the polymer is composed of pore channel and CD cavity. Analysis of the complex of the species with CD ( Figures S4, S5) or EPI oligomers ( Figure   S10) demonstrated that both units could bind the species via hydrogen bonding and/or s-p interactions. The substitution of variable x 9 rendered a significant increase in the regression coefficient (R 2 ) (Equation 5 versus Equation 3). This improvement suggested that variable x 9 was a good indicator of inclusion potential relative to variables x 1 and x 8 . A concomitant of the regression improvement was an increase in the weight coefficient of pore volume. It is probably due to the dual roles (i.e., inclusion complexation and pore donator) of CD, both which are derived from the cavity of CD.
The adsorption of OTC 2 was mainly related to cross-linking degree (Equation 4). High cross-linking degree commonly causes a dense network structure [19,50]. Hydrogen bonding forces and sp interactions were identified for the species OTC 2 binding to EPI oligomers ( Figure S10). This correlation probably illustrates the presence of network capture mechanism. In contrast, the substitution of variable x 9 didn't affect the MLR analysis of OTC 2 . This indicated an insignificant role of CD inclusion for OTC 2 , which was in agreement with the finding of its poor CD inclusion potential (Table 1).
For the other two species that are present in acidic or basic solutions (i.e., OTCH 2 + and OTC 22 ), the MLR analysis wasn't performed because the K d almost approached zero ( Table 2). The poor adsorption affinity of the two species could be attributed to the fact that the adsorption sites (i.e., CD cavity) were not readily accessible for the two species, because they were difficult to diffuse to hydrophobic surfaces (Table S7). Furthermore, they could interact with hydroxyl groups in the network but in the less magnitude, as compared with the other two species. It is due to the decrease in the number of hydrogen bonding for the species following the order of OTCH 2 (10).OTCH 6 (8).both species (3) ( Figure S10).

The nature of pH-dependent adsorption of OTC
The MLR analysis of apparent K d of OTC with CDP properties was performed, in which the apparent K d values of OTC was calculated on the basis of Equation 2 (showed as lines in Figure 3) when pH varied from 3.0 to 10.0 with an interval of 0.2. The apparent adsorption of OTC at varying pH was correlated with one to five of all eight factors of concern, when CD content and inclusion constant were set as independent variable ( Figure 5I). Surprisingly, the substitution of variable x 9 reduced the number of influencing factors, only including pore volume and/or effective inclusion, while regression coefficients (R 2 ) were comparable ( Figure 5II).
A combination of the MLR analysis for OTC ( Figure 5II) and molecular recognition in adsorption of the species (Equations 3-5) may illustrate pH-dependent adsorption of OTC onto CDPs ( Figure 6). Specifically, the term K d correlates with pore volume at pH 3.0-4.2, where OTCH 6 and OTCH 2 + prevail and are adsorbed via pore-filling mechanism. As the solution pH rises up to 6.8, pore volume and effective inclusion capacity become the influencing factors corresponding to pore-filling and CD inclusion mechanisms, respectively. It is reasonable because OTCH 6 is the predominant species in the range of pH. However, the standardized coefficients of both influencing factors ( Figure 5II) gradually deviate from those in Equation 5 with pH increasing. It is probably due to an increase in the fraction of OTC 2 over pH. Surprisingly, the term pore volume is the sole influencing factor at pH 7.0-10.0, indicating the presence of pore-filling mechanism. It is in contradiction to the fact that the species OTC 2 predominates at pH 8.0-9.0 and its adsorption proceeds via network capture mechanism. This may be due to a combination of pore-filling, CD inclusion and network capture mechanisms, associated with the coexistence of two or three species in this pH range.

Selection of adsorbent to separate OTC from water
The selection of the polymers may be determined on the basis of the understanding of the nature of pH-dependent adsorption of OTC ( Figure 5II). Adsorbent c-HP-CDP with advantages of high pore volume and efficient inclusion potential was selected to separate and remove OTC from simulated wastewater (pH 5.0) and natural water (pH 7.0) (Text S4). The adsorbent removed 80% of OTC at 20 mg/L (pH 5.0) (Figure 7). In contrast, other conventional natural and synthetic adsorbents are observed to exhibit similar or lower adsorption efficiency for OTC, such as activated carbon (ca. 68%) [7,51], multi-wall carbon nanotubes (ca. 70%) [52], and aluminum oxide (44%) [53]. Furthermore, the  polymer exhibits a removal efficiency of above 70% for the antibiotic even at environmentally relevant levels (Figure 7).

Conclusions
The finding of the species-specific adsorption reveals the pHdependent adsorption behavior of antibiotic oxytetracycline onto adsorbent CDPs. The adsorption sites of the polymer such as CD cavity, pore channel, and network are accessible for the species to different extents, associated with the coordination of a variety of weak interactions (e.g., hydrogen bonding, hydrophobic interactions, van der Waals forces and electrostatic interactions). The species OTCH 6 commonly has high adsorption affinity via porefilling and CD inclusion mechanisms, OTC 2 exhibits moderate affinity via network capture mechanism and the other two species are poorly adsorbed due to the limiting accessibility of the inclusion sites as well as other adsorption sites. The resulting immobilization of the species onto CDPs is determined by the K d and fraction of the species in solution, representing a function of pH. This may be the cause for the pH-dependence of adsorption of OTC. Moreover, the mathematical equations are derived to quantitatively relate adsorption affinity of OTC at varying pH to adsorbent properties. These relationships, together with the molecular recognition into adsorption of the species, provide theoretical guidelines for the selection of CDPs to remove the antibiotic. Additionally, this study indicates that CDPs are promising adsorbents for IOCs removal, taking advantage of discernible adsorption interactions, high removal efficiency, and strong resistance to the interference of coexisting natural components.

Supporting Information
Text S1 Adsorption kinetics and isotherms of OTC.