Synthesis and Characterization of Molecularly Imprinted Polymers for Phenoxyacetic Acids

Abstract

2-methylphenoxyacetic acid (2-MPA), 2-methyl-4-chlorophenxyacetic acid (MCPA) and 4-chlorophenoxyacetic acid (4-CPA) were imprinted to investigate the cross-selectivities of molecularly imprinted polymers (MIPs). The result indicates that 2-MPA, which is similar in shape, size and functionality with phenoxyacetic herbicides, are suitable to be used as a suitable template to prepare the MIPs for retaining phenoxyacetic herbicides. To study the ion-pair interactions between template molecules and functional monomer 4-vinylpiridine (4-VP), computational molecular modeling was employed. The data indicate that the cross-selectivities of MIPs for phenoxyacetic acid herbicides depend on the binding energies of complexes.

1. Introduction

Phenoxyacetic acids are widely used in agriculture as herbicides or hormone mimics. This series of compounds include 2,4-dichlorophenoxyacetic acid (2,4-D), 2-methyl-4-chlorophenxyacetic acid (MCPA), 4-chlorophenoxyacetic acid (4-CPA), etc. Residues of these compounds in soil, ground water, and plants are of concern due to their polarity. Methods typically used to measure these herbicides and their metabolites in environmental samples rely on extraction and analysis via gas chromatography [1–3]. However, due to the physical and chemical properties of these compounds, they tend to be difficult to separate from complex sample matrices. The synthesis of a material that selectively recognizes these serial compounds would be significant.

Molecularly imprinted polymers (MIPs) were polymerized as shown in figure 1. In solvent, the aim is to develop a complex in which the target molecule is recognized based on whether the functional group is a carboxyl, amido or acyl to form a complex by hydrogen bonding, and electrostatic interaction or hydrophobic effects. A crosslink polymerization is performed to fix the desired complex. After the template molecules are removed, binding cavities that recognize the applicable template are left to rebind selected molecules. Generally, a functional monomer has an appreciative group to match the active site of the template molecule, which makes the complex stable and improve the selectivity of the MIPs. The strength of the interaction between the template and the monomer determines the accuracy of the binding cavity and its selectivity.

In previous work, 2,4-D and 2,4,5-trichlorophenoxyacetic acid (2,4,5-T) were imprinted [4–16]. The result indicates that MIPs showed recognition capability for targeted molecules because of the ion-pair and hydrophobic interactions between the template and the monomer 4-vinylpyridine (4-VP). Baggiani's research demonstrates that 2,4,5-T MIPs have a good potential to serve as phenoxyacetic acid herbicides templates in solid-phase extraction (SPE) [7].

However, MIPs are synthesized by using a target molecule as template. When they are used for solid phase extraction (SPE), the leakage of a trace amount of the imprinted molecules remaining in the MIPs has hindered the accurate and precise assay of the target analyte [17]. To overcome this problem, a structural analogue of the analyte is used as a “dummy template” molecule and separated from the target molecules during chromatographic analysis [18, 19]. The MIPs prepared by dummy template are dummy molecularly imprinted polymers (DMIPs). In this work, herbicides MCPA, growth regulator 4-CPA and expectorant 2-MPA (Figure 2) were imprinted to investigate crossselectivity of these MIPs and choose the suitable dummy template to prepare the MIPs for retaining phenoxyacetic herbicides.

Computer modeling was used previously to design MIPs by calculating the binding energy between the template and the functional monomers [20–22]. In order to investigate the recognition mechanism of ion-pair interactions, computational molecular modeling is employed here to simulate the combination of a 1:1 template: monomer complex.

2. Results and Discussion

2.1 HPLC and Cross-selectivity

MIPs for 2,4-D have been polymerized and characterized by radioligand binding assay [4, 5], SPE [7] or HPLC previously [6, 9, 10]. However, the bleeding template phenomenon could result in inaccurate quantification, especially at low concentration residue levels. In this work, a series of MIPs for MCPA, 4-CPA and 2-MPA are studied by HPLC to investigate the cross-selectivities of these MIPs and choose the suitable dummy template to prepare the MIPs for retaining phenoxyacetic herbicides. The selectivity of MIPs and blank polymers columns were evaluated by eluting 4-CPA, MCPA, 2-MPA, 2,4-D and benzoic acid (BA) with acetonitrile/water/acetic acid (93+5+2, v/v). Their molecular structures are shown in Figure 2. The structure of BA differs from those of the other templates to investigate the imprinting effect by this polymerization. Shown in

Table 1. Capacity factors(k′), separation factors (α) and retention indices (RI) on MCPA MIPs and blank polymer.

	k′ blank	k′ MIPs	αblank	αMIPs	RI
MCPA	0.668	1.411	1	1	1
4-CPA	0.778	1.581	0.859	0.892	0.962
2-MPA	0.468	0.864	1.427	1.633	0.874
BA	0.504	0.535	1.325	2.637	0.503
2,4-D	0.982	2.404	0.680	0.587	1.158

–

Table 3. Capacity factors(k′), separation factors (α) and retention indices (RI) on 2-MPA MIPs and blank polymer.

	k′ blank	k′ MIPs	αblank	αMIPs	RI
2-MPA	0.468	0.943	1	1	1
4-CPA	0.778	1.702	0.602	0.554	1.086
MCPA	0.668	1.410	0.700	0.669	1.048
BA	0.504	0.605	0.928	1.559	0.596
2,4-D	0.982	2.475	0.476	0.381	1.251

are the capacity factors (k′), separation factors (á) and retention indices of various analytes on blank polymers and MIP columns. The k′ of the prepared MIPs are evaluated as:

$$\text{capacity factors}: k\prime =(t-t_0)/t_0,$$

(1)

Where t is the retention time of the analyte and t₀ is the retention time of the void marker [23]. The t₀ is indicated by the methanol signal. The separation factor and retention index are evaluated as:

$$\text{Separation factor: }\alpha =k_{\text{template}}^{\text{'}}/k_{\text{analog}}^{\text{'}}$$

(2)

$$\text{Retention index: RI}={\alpha }_{\text{blank}}/{\alpha }_{\text{MIPs}}$$

(3)

The k′ of 4-CPA, MCPA and 2-MPA on MIPs are higher than those on a blank column (Table 1–3), but the improvement of k′ of BA is subtle. The higher k′ value means a higher retention capability of the MIP. The fact that BA displays higher affinity than 2-MPA on the blank polymers but is opposite on the 2-MPA MIP demonstrates that the oxygen bridge sub-structure of 2-MPA is imprinted successfully. Without the oxygen bridge sub-structure, BA whose phenyl is too big to fit into the narrow part of the cavity matching the oxygen bridge sub-structure, shows poor retention on the MIPs due to the weakened interaction between the carboxyl and the basic nitrogen of the pyridyl. The retention index of 4-CPA on MCPA and that of MCPA on 4-CPA imply that the selective affinity of MIPs for template is improved by shape selectivity.

Although generally DMIPs are inferior to MIPs in terms of selectivity for the target analyte [17], past research indicates that the property of the probe molecule can disturb the shape selectivity of the MIPs when the acidity of molecule is stronger than that of the template [21]. In this case, the retention indices of 2,4-D, 4-CPA and MCPA are higher than those of 2-MPA on 2-MPA MIP columns. Thanks to the absence of a strong active site on the aromatic parts of the templates, exact imprints are not possible. The binding cavities can't discriminate between 2,4-D, 4-CPA, MCPA and 2-MPA when the interaction between the analogs and the binding site is too strong. The k′ of 2,4-D, MCPA or 4-CPA on blank polymers also indicate their natural affinities with 4-VP.

The chromatographic results (Figure 3, 2-MPA MIPs (a), blank polymer (b)) indicate that the MIPs prepared using 2-MPA as a dummy template show selective affinity for phenoxyacetic herbicides. On the DMIPs, BA can be separated from the phenoxyacetic acids and 2-MPA can be separated from phenoxyacetic herbicides such as 4-CPA, MCPA or 2,4-D. So, 2-MPA, which is similar in shape, size and functionality with phenoxyacetic herbicides, is suitable to be used as a dummy template to prepare MIPs for retaining group phenoxyacetic herbicides.

2.2 Computational molecular modeling

Computational molecular modeling is employed to study the selectivity of these MIPs. Previous studies have shown that pbe1pbe [24, 25], a hybrid functional which contains a portion of Hatree-Fock exchange, yields the best performance for predicting hydrogen binding energies [26]. So we chose this method in our study. The basis set 6–311g (d, p) is also good enough to provide the best results.

For each template the most stable template–monomer complexes are located. The optimized geometries for the most stable complexes are show in Figure 4, which are the same as we expected. Their electronic stabilization energies relative to isolated fragments, ΔE, are calculated through the equation:

$$\Delta \text{E}=\text{E}(\text{template}–\text{monomer})-[\text{E}(\text{template})+\text{E}(\text{monomer})]$$

(4)

The interaction energies obtained by docking the template and the monomer structures minimized in vacuum are reported in

Table 4. Computational modeling data: binding energy ΔE(kcal/mol) and bond length r(Å) in vacuum

	ΔE (kcal/mol)	rA(Å)	rB(Å)	rC (Å)
4-CPA	-14.47979	0.967	1.016	1.664
MCPA	-14.56388	0.967	1.016	1.663
2-MPA	-14.15600	0.967	1.013	1.666

^Ar: bond length of O-H on carboxyl of template

^Br: bond length of O-H on carboxyl of complex

^Cr: bond length of H…N of complex

. According to the binding energies and the bond length of H…N, the ion-pair interactions between 4-VP and 4-CPA or MCPA are higher than that of 2-MPA. This evidently leads to 4-CPA and MCPA retaining on 2-MPA MIPs. When the binding energies of the complexes are the same, the retention behavior of the templates on each MIP column depends on the shape selectivity, e.g., 4-CPA and MCPA. The bond length changes shown in Table 4 indicate that the acidic protons leave the oxygen atom and approach the nitrogen atom to form an ion-pair.

3. Experimental Section

3.1 Chemicals and apparatus

Compounds 4-CPA(99%), MCPA(95%), 2,4-D(99%) and 2-MPA(99%) were supplied by the Pesticide Analysis Laboratory, China Agricultural University 4-VP. Ethylene glycol dimethacrylate (EGDMA) and chromatographic solvents were obtained from J&K Chemical Ltd. All other reagents were from Beijing Chemical Reagents Company. 4-VP was distilled at reduced pressure and stored in the dark at –18 °C. EGDMA was purified by basic alumina column immediately before use. Standard samples were prepared by dissolving 10 mg templates and analogs into 10 mL methanol. HPLC analyses were performed using an Agilent Model 1100 Series liquid chromatograph.

3.2 Preparation of Polymers

The following procedure was used for the preparation of all MIPs. In a glass bottle, 1mmol of template was dissolved in water/methanol (1 mL+4 mL). To this solution was added 4mmol 4-VP, 20mmol EGDMA, and 0.51mg 2,2′-azobisisobutyronitrile (AIBN). The blank polymer was formulated in a similar method, without template molecule. The bottle was subjected to sonication, purged by nitrogen gas, and sealed. The sample was placed in a 60 °C water bath for 24 h.

3.3 HPLC experiments

The polymers were ground by a machine mill, the particles were sized using wet sieving by acetone, and the particles between 20 and 80 μm were collected. The particles were slurry packed into stainless steel columns (length, 15.0 cm, i. d. 4.6 mm) for the chromatographic experiments. The templates were removed online by acetonitrile/acetic acid 90+10 (v/v), then by mobile phase acetonitrile/water/acetic acid 93+5+2 (v/v) until baseline equilibrated. The flow rate in all cases was 1.0 mL/min, with a UV detector at 270 nm for 2-MPA and BA, 285 nm for 2,4-D, 4-CPA and MCPA. Sample injections were 5 μL of a 1.0 g/L mixed solution of analytes.

3.4 Computational molecular modeling

The theoretical binding energies of templates 4-CPA, MCPA and 2-MPA with 4-VP in vacuum were calculated using Gaussian03 [27]. The geometry of each structure has been fully optimized by using density functional method (pbe1pbe) with 6–311g (d, p) basis set. All the geometric parameters of the possible stationary points that we found have been located and characterized by the number of imaginary frequencies. All the energies of the structures located have Zero-point Energy (ZPE) correction.