Computational and experimental designing of imprinted sorbent for the determination of nitroxidative stress products: an analysis of 4-hydroxyphenylacetic acid conversion

It this paper, the comprehensive design process was carried out to fabricate selective, molecularly imprinted polymer (MIP). The material was used as a sorbent in the optimized analytical method, aimed at verifying the hypothesis that the conditions of an analytical process could convert 4-hydroxyphenylacetic acid to 4-hydro-3-nitrophenylacetic acid, comprising a series of novel studies. The design stage consisted of the analysis of specificity of the MIP, synthesized from the cross-linker and used five various functional monomers independently, in the presence of four different templates. The MIP from 1-vinylimidazole, imprinted by 4-hydroxyphenylacetic acid, revealed the highest specificity in relation to 4-hydro-3-nitrophenylacetic acid, with an affinity factor equal to 3, and the highest selectivity from a group of structurally similar and biologically important biomolecules. The theoretical analysis revealed that electrostatic interaction between the analyte and the polymer matrix enhanced selectivity. The physicochemical characterization showed the specific surface area of the MIP as being equal to 368.6 m2 g−1, and the presence of nitrogen atoms at the level of 6.80% wt., confirming the monomer residue in the material structure. The MIP was applied in the solid phase extraction protocol, allowing for the analysis of 4-hydroxy-3-nitrophenylacetic acid in a human urine sample. Finally, the conversion of 4-hydroxyphenylacetic acid in human urine in nitrate and nitrite salts at low pH conditions revealed an almost twofold increase in 4-hydro-3-nitrophenylacetic acid to 775 ± 81 ng L−1. The results also confirmed the applicability of the new MIP sorbent for the purpose of analysis of low levels of analyte, present in the complex sample.


Introduction
Advanced materials are in high demand in the field of analytical chemistry, as they could improve the characterization of analytical methods, viz. they allow for lower limits of quantification (LLOQ), improve accuracy, enhance selectivity and minimize matrix effects. Among the various materials, the polymeric sorbents, with their predetermined capabilities, are considered the most attractive. Here, molecularly imprinted polymers (MIPs) are recognized as a valuable class of materials, characterized by prearranged properties, due to the three-dimensional cavities formed in the polymer network by the template-tailored synthesis [1]. One of the main attributes of imprinted materials is high specificity, derived from the spatial modification of the polymer surface by the template, an advantage that has paved the way for widespread application of MIPs in different kinds of sensors, as recognition elements [2][3][4] in synthesis and catalysis [5,6] and as potential drug delivery devices [7][8][9]. Nevertheless, MIPs are predominantly known as highly useful materials for separation purposes, due to their capability to clean up complex samples [10][11][12][13]. However, the adsorption on the resulting MIPs is governed by different factors, making the process less predictable and often unsatisfactory. For that reason, the scrupulous designing step and the comprehensive optimization of synthetic protocols are required to obtain the most effective sorbent. Alongside the experimental optimization, a theoretical analysis could be used as versatile tool to rationalize the process and to explain the behaviour of the material [14,15].
The analysis of nitratively-altered macromolecules in biological samples has attracted attention due to their role in the pathogenesis of various diseases, such as asthma, diabetes or neurological disorders [16]. The elevated levels of reactive nitrogen species in cells could lead to the damage of cellular components, provoking pathological actions. In proteins, most of nitration reactions appear on a tyrosine aromatic ring, which is easily nitrated to 3-nitrotyrosine [17][18][19][20][21][22]. This compound is characterized by higher acidity than tyrosine, which affects the conformational stability of protein. As a result, an accelerated degradation of modified proteins is observed and higher levels of free 3-nitrotyrosine are detected, together with its main metabolite, 4-hydroxy-3-nitrophenylacetic acid (NHPA), which is excreted by urine. However, the NHPA is not exclusively formed from 3-nitrotyrosine. Pannala and co-workers [23] found that the levels of NHPA could double by the endogenous nitration of 4-hydroxyphenylacetic acid, present in saliva. The conversion of 4-hydroxyphenylacetic acid to NHPA occurs in the stomach, following the consumption of nitrate-rich products. Thus, it was stated that the urinary NHPA levels may be a good marker for the formation and metabolism of 3-nitrotyrosine, but dietary factors must be taken into account, assessing the nitroxidative stress products in vivo. The conversion of 4-hydroxyphenylacetic acid could also take place during in vitro analysis, affecting analytical results. Thus, to verify the hypothesis that the conditions of the analytical process could lead to an increased level of NHPA, the new MIP was designed. It was applied in the solid phase extraction (SPE) as a selective sorbent for the purpose of in vitro analysis of the conversion of 4-hydroxyphenylacetic acid to NHPA in human urine samples, under nitrite/nitrate-low pH conditions. The process was supported by theoretical analysis, the physicochemical characterization and SPE optimization.
Human urine samples were obtained from a healthy volunteer (a woman, of 20 years of age), who was enrolled in this study after signing a written informed consent. The subject, a non-smoker, did not consume alcohol or take any medication within two weeks of the start of the study. The urine sample was collected in the morning, after overnight fasting.

Polymers
Radical bulk polymerization was performed to obtain the polymers coded as MIP1-MIP5, imprinted by HPA, PHE, NBA and NFE-codes are extended in subscript. The templates acted as structural analogues of NHPA in a so-called pseudo-template approach [11]. The corresponding non-imprinted polymers, NIP1-NIP5, were prepared in the same way but without the template. The experimental quantities of the reagents, used for the preparation of polymers, are listed in Supplementary Information (Table S1).
Briefly, the chosen template, the appropriate functional monomer and the EGDMA were dissolved in a mixture of toluene and methanol (3:2 v/v) in a thickwalled glass tube. The molar ratio of the template to the functional monomer and the cross-linker was 1:4:20. Finally, the polymerization initiator was added. The homogeneous solutions were purged with nitrogen for around five minutes (min.), then the glass tubes were sealed. Subsequently, polymerization was carried out in a nitrogen atmosphere for 24 h at 88°C. The bulk rigid polymers were ground and wet-sieved into particles below 45 lm in diameter. Fine particles were separated by repeated decantation from acetone. The templates were removed from the polymers, using continuous extraction in a Soxhlet apparatus (24-36 h, 80 mL, methanol), followed by a washing sequence with 1% aq. ammonium hydroxide in methanol (5 mL/10 mg of sorbent) and methanol (5 mL/10 mg of sorbent). The template removal was monitored by liquid chromatography, coupled with tandem mass spectrometry (LC-MS/MS). The particles were dried in a vacuum at room temperature.

Instruments
Instrumental analysis was performed using an Agilent 1260 Infinity System (Agilent Technologies, Santa Clara, CA, USA), equipped with a degasser, an autosampler and a binary pump, coupled to a QTRAP 4000 hybrid triple quadrupole/linear ion trap mass spectrometer (AB Sciex, Framingham, MA, United States). The turbo ion spray source was operated in negative mode. The curtain gas, ion source gas 1 and ion source gas 2, were set at 241 kPa, 414 kPa, 276 kPa and ''high'' instrument units (4.6 9 10 -5 Torr), respectively. The ion spray voltage and source temperature were 4,500 V and 600°C, respectively. The target compounds were analysed in multiple reaction monitoring (MRM) mode. The quantitative MRM transitions, declustering potential (DP) and collision energy (CE) are presented in the Supplementary Information (Table S2). Chromatographic separation was achieved with a Bionacom BionaCore C18 UFPLC Column (100 mm 9 4.6 mm, 2.7 lm) from Bionacom (Coventry, United Kingdom). The column was maintained at 40°C at a flow rate of 0.5 mL min -1 . The mobile phases consisted of water with 0.2% formic acid as eluent A, and acetonitrile with 0.2% formic acid as eluent B. The gradient (%B) was as follows: 0 min. 5%, 0.5 min. 5%, 7 min. 90% and 9 min. 90%. The re-equilibration of the column to the initial conditions lasted for 2 min. The injection volume was 10 lL.
The method was characterized by acceptable precision (5.8% within-run; 9.9% between-runs) and accuracy (110% within-run; 114% between-runs). The limit of detection (LOD) and the limit of quantitation (LOQ), determined as the minimal concentration providing a peak height three and five times the baseline noise, was calculated in the urine as 85 ng L -1 and 142 ng L -1 , respectively.
The surface morphology analysis, using scanning electron microscopy (SEM) with a Merlin FE-SEM (Zeiss, Germany) and the X-ray electron dispersive spectroscopy (EDS) analysis, using an EDS X-ray detector (Brucker, Germany) were performed at the Faculty of Chemistry, University of Warsaw, Poland. The samples were Au/Pd sputter-coated, prior to the SEM analysis. The porosity data were determined using the adsorption isotherm of N 2 at 77 K (BET) on an ASAP 2420 system (Micromeritics Inc., USA) at the Faculty of Chemistry, Maria Curie-Skłodowska University, Lublin, Poland. The infra-red (FT-IR) spectra were recorded on a Nicolet iS50 FT-IR (Thermo Fisher Scientific, Waltham, MA, United States) at the Biological and Chemical Research Centre, University of Warsaw, Poland. The thermogravimetry analyses (TGA) were performed at the Faculty of Chemistry, Warsaw University of Technology, Poland on a Q600 thermogravimetric analyser (TA Instruments, USA) in an argon atmosphere, with a heating rate of 5°C min -1 . Fluorescence detection was performed on the spectrofluorophotometer RF-6000 (Shimadzu, Duisburg, Germany), using a 10-mm crystal cuvette. The 3D spectra were recorded using LabSolution software.

Binding studies
Stationary binding experiments were performed to evaluate the binding abilities of the MIPs and the NIPs in relation to the NHPA. These were carried out as previously described [24]. Regarding the NHPA, polypropylene tubes with a volume of 10 mL were filled with 10 mg of MIP1 HPA /MIP1 PHE /MIP1 NBA / MIP1 NFE -MIP5 HPA /MIP5 PHE /MIP5 NBA /MIP5 NFE or NIP1-NIP5 particles. A volume of 5 mL of 100 lg L -1 methanol-water (85:15 v/v) standard solution of NHPA was added to each tube. The tubes were sealed and oscillated, using a shaker at room temperature for 2 h. Then, the tubes were centrifuged, and the aliquots of supernatant were used to analyse the unbound amounts of each compound by LC-MS/ MS, using the calibration lines. The amounts of NHPA bound to each polymer were calculated by subtracting the unbound amounts from the initial amounts.
For the isotherm analysis, the same procedure was applied to the MIP5 HPA and NIP5 with the addition of different concentrations of NHPA in methanolwater (85:15 v/v), ranging between 25 and 1000 lg L -1 .
The aforementioned selectivity studies were performed on the MIP5 HPA in the non-competitive experiments for the analytes: HPA, PHE, NBA, NFE, 3NT, 3CT. The polypropylene tubes, with a volume of 10 mL, were filled with 10 mg of MIP5 HPA, and a volume of 5 mL of 7.5 lmol L -1 methanol-water (85:15 v/v) standard solution of each analyte, was added to each tube.

Solid phase extraction
For the optimization of the SPE process, an empty SPE 1-mL cartridge was filled with a mass of 10 mg of MIP5 HPA and secured by fibreglass frits. Subsequently, the following steps of the optimization of SPE were carried out as follows: conditioning with 1 mL of water, adjusted to pH 4, loading of 1 mL of human urine (adjusted to pH 4), spiked with NHPA standard solution to 3.33 ng mL -1 , washing with 1 mL of water, adjusted to pH 4 prior to adding 1 mL of hexane or 15% methanol in water or 15% methanol in water, adjusted to pH 4 or 15% methanol in water with 1% of ammonium hydroxide (before applying each solvent the sorbent was air-dried for 5 min.) followed by eluting either by t-butyl-methyl ether, ethyl acetate or methanol. Following elution, the solvent was evaporated to dryness and was reconstituted in 100 lL of methanol-water (1:1 v/v) with IS (concentration of 100 ng mL -1 ), prior to analysis on the UPLC-MS/MS system. The recovery was calculated, based on an analysis of the unspiked sample of human urine, which was carried out in the same manner as described above, except for the addition of 70 lL of IS and 33 lL of NHPA standard solution (concentration of 100 ng mL -1 ) to eluate.

Analysis of HPA conversion
For the purpose of analysis of the HPA conversion to a volume of 900 lL of human urine, adjusted to a physiological pH of the gastric environment (pH 2), a volume of 100 lL of sodium nitrite (concentration of 100 lmol L -1 ), sodium nitrate (concentration of 10 mmol L -1 ) and HPA (concentration of 100 mg L -1 ) were added. The mixture was left for 30 min. prior to SPE on the MIP5 HPA , according to the optimized protocol. The eluates were collected, and the solvent was evaporated to dryness and was reconstituted in 100 lL of methanol-water (1:1 v/v) prior to analysis on the UPLC-MS/MS system. The human urine samples were analysed prior to testing, to assess the level of NHPA.

Theoretical simulations
All atom molecular dynamic (MD) simulations were performed using the BIOVIA Discovery Studio 2019 software package [25]. Before the molecular simulations of the prepolymerization, polymerization and adsorption processes, the optimization of the molecular structure and so-called ESP (electro-static potential) charges, a calculation of all compounds was performed using density functional theory (DFT) with a B3LYP/6-311 ? G(d,p) hybrid functional, implemented in the Gaussian 09 program [26]. For all molecular systems, the construction and simulation processions were carried out and parametrized with a CHARMM force field [27]. Simulated systems were initially constructed using the Packmol software [28] to obtain random starting geometries. Firstly, the systems were energy-minimized (molecular mechanics (MM) method) with 100 steepest descent and 10,000 conjugate gradient steps to remove bad contacts. The processes were applied until the RMS (root-mean-square) gradient of the structure fell below 0.01 kcal mol -1 Å -1 . The MD protocols contained a heating step, performed for 100 ps, with time steps of 1 fs. The system was heated from 0 to 300 K. Prior to the production stage, isothermal equilibration was performed for 100 ps, at a temperature of 300 K. The Leapfrog Verlet integration and SHAKE [29] algorithms were used during the simulation process. The production run was conducted for 5 ns in the NVT ensemble (constant-volume/constanttemperature dynamics) at 300 K, and the coordinates were recorded every 10 ps. Trajectory file data, generated from the NVT MD simulation, have been used in all the calculations and analyses presented in this research.
In the first stage of polymers creation, the models of prepolymerization complexes (coded as PC) were constructed between the template and the monomers. To obtain the starting structures, the boxes with the chosen template molecule (HPA, PHE, NBA or NFE), surrounded by four molecules of appropriate monomer molecules (1) -(5) were created. Next, 20 molecules of the cross-linker and an appropriate number of solvent molecules (Table S1) were added, and complete starting structures were created using Packmol. Next, MM and MD procedures were applied to obtain models of prepolymerization systems with intermolecular interactions, formed between the template and monomer or cross-linking agents. The number of reagents and solvent molecules were chosen to mimic their molar ratio, used during the synthetic process.
In the next stage of the modelling procedure, the creation of polymer chains with the binding sites was performed (coded as MIP). The optimized structures of the prepolymerization complexes with minimum potential energy values, were chosen for further investigation. Single bonds between vinyl groups of monomer and cross-linker molecules in the prepolymerization systems were created. The C atoms of the monomer or cross-linker vinyl groups which were closest to each other were connected, but also took into account the fact that all molecules of the monomer and cross-linker in the studied system, should create one cross-linked polymeric chain. Then, the H atoms were added to the chain structure. The operation described, mimicked the polymerization reaction in the synthetic process. So called ESP charges were calculated for the creation of the polymeric chain structures. Next, the MM and MD procedures were repeated for the constructed systems, which consisted of the polymeric chain, the template and the solvent, to form a specific binding site in the polymeric matrix.
In the last step of computational calculations, the adsorption process was simulated. The optimized structures of the polymeric chains with the minimum potential energy values, were selected for further investigation. The template was removed from the models of cavities and the empty spaces were proposed as the computer models of binding sites. Then, the molecule of the NHPA was inserted into the models of the MIPs' cavities (replacing the template molecule). The solvent, consisting of 418 molecules of methanol and 167 molecules of water, was added to form starting structures using Packmol, then the MD simulations were carried out (mimicking the experimental conditions of adsorption). During selectivity studies, the chosen analytes, namely HPA, PHE, NBA, NFE, 3NT and 3CT were inserted into the MIP5 HPA cavity model, using the same solvent as was applied to the NHPA adsorption studies. Then starting structures were obtained and finally, the MD calculations were performed.

Calculation of experimental and theoretical binding parameters
The binding capacities (B, lg g -1 ) of the MIPs and NIPs are calculated according to Eq. 1: The distribution coefficients (K D , L g -1 ) for the MIPs and NIPs are calculated according to Eq. 2: where V represents the volume of the solution (L), C i represents the initial solution concentration (lg L -1 ), C f represents the solution concentration after adsorption (lg L -1 ) and m is the mass of the particles (g). The affinity factors (AF) are obtained according to Eq. 3: The adsorption isotherms were characterized using the Langmuir model, transformed into the Scatchard equation, as shown in Eq. 4: where B max is the total number of binding sites, K d is the dissociation constant, B is the bound amount of the analyte and F is the unbound amount of the analyte. The system, which fits well with a Langmuir model, gives a straight line on the Scatchard plot, with a slope equal to -(1/K d ) and a y-intercept, equal to B max /K d . It also fits with the Freundlich model, represented by Eq. 5: where B is the bound amount of the analyte, F is the unbound amount of the analyte, a is a measure of the capacity (B max ) and m is a heterogeneity index.

Results and discussion
Design of sorbent Different functional monomers and templates were analysed in the material designing step. A set of 20 MIPs were synthesized from five different monomers, viz. (1)-(5) involving four different templates, viz. HPA, PHE, NBA and NFE. These templates acted as structural analogues of the target analyte, NHPA, in the so-called pseudo-template strategy [11]. The NIPs were also obtained.
Effect of functional monomer on binding capacity in relation to the NHPA in relation to the NHPA was carried out to evaluate the effect of the composition of the polymer network on the adsorption behaviour. The stationary experiments were carried out and the B (lg g -1 ) values of the MIPs and the NIPs were calculated, according to Eq. 1. The calculation of K D (L g -1 ), according to Eq. 2, followed by the calculation of AF according to Eq. 3, was also carried out (Table S3). Figure 1 presents the B of MIP1 HPA /MIP1 PHE /MIP1 NBA / As can be seen ( Fig. 1), the lowest B with regard to the NHPA were noted in the case of MIPs prepared from (1) (from 3.30 to 12.2 lg g -1 for MIP1 NFE and MIP1 NBA , respectively). Only slightly higher values of B were observed for polymers prepared from (2) (from 9.81 to 14.0 lg g -1 for MIP2 NFE and MIP2 NBA , respectively) and (3) (from 7.96 to 16.6 lg g -1 for MIP3 PHE and MIP3 HPA , respectively). Significantly higher B were noted in the case of the MIPs synthesized from (4) (from 9.83 to 26.90 lg g -1 for MIP4 NFE and MIP1 HPA , respectively). The highest values of B characterized the MIPs build up from (5), acting as the functional monomer (with the highest B among all tested polymers for MIP5 HPA, equal to 64 lg g -1 ). A similar trend was observed for the NIPs. Those results showed that the adsorption process of the NHPA strongly depended on the composition of the polymer and (5) was the most effective monomer among all tested compounds, in providing MIPs in relation to the NHPA. These results also concurred with the observations of Baggiani and co-workers [30] who postulated that if the NIP indicated binding properties towards a target analyte, the similarly composed MIP would reveal a higher B and a significant imprinting effect.
Thus, in the next step, the impact of different templates and the structural analogues of the target analyte on the specificity of resulting MIPs was investigated.

Effect of template on specificity
As shown (Table S3), the highest specificity (expressed as AF) was noted for MIP5 HPA, prepared from (5) in the presence of HPA with AF = 3.03. The polymers prepared from (4) also revealed satisfactory specificity; MIP4 NBA was prepared in the presence of NBA as the template was characterized by sufficient specificity, AF = 2.51, together with MIP4 PHE, prepared in the presence of PHE, AF = 1.84. The results could be explained by the structure of the template and its interactions with respective monomers. It could be postulated that the interactions between the imidazole nitrogen atoms and the methylene carboxylic group from the HPA (which is characterized by higher acidity when compared to NBA) efficiently stabilized the prepolymerization complexes. The remaining tested systems failed to provide the satisfactory specificity of the resulting MIPs.
Thus, to reveal the molecular bases of the interactions between analysed templates and monomers, as well as to understand the mechanism of adsorption, a theoretical analysis was carried out.

Theoretical proof of specificity
In order to provide an insight into the mechanism that governs the adsorption of NHPA on selected polymeric matrices, the following analyses were carried out: (a) the analysis of interactions between the respective templates, monomers and cross-linker in the prepolymerization complexes of PC5 HPA , PC4 NBA , PC4 PHE , PC1 NBA and PC1 NFE (it imitates the formation process of the complex), (b) the analysis of interaction of the template (HPA, NBA, PHE or NFE) in the cavities of MIP5 HPA , MIP4 NBA , MIP4 PHE , MIP1 NBA and MIP1 NFE (it imitates the interaction of the template in polymerized systems only), (c) the analysis of interactions of NHPA in the cavities of MIP5 HPA , MIP4 NBA , MIP4 PHE , MIP1 NBA and MIP1 NFE (it imitates the process of adsorption of the target analyte inside the cavity, after template removal). Those five systems were selected taking into account the following: the highest B and specificity (PC5 HPA / MIP5 HPA ), the second and third highest specificity (PC4 NBA /MIP4 NBA and PC4 PHE /MIP4 PHE ), the fourth highest specificity but low B (PC1 NBA / MIP1 NBA ), as well as the lowest B and lack of specificity (PC1 NFE /MIP1 NFE ). The impact of the solvent on the interactions is also discussed. In the PC5 HPA model (Fig. 2a), two types of interactions between the template (HPA) and the monomer (5) or cross-linker were observed, viz. hydrogen bonds and hydrophobic interactions. Two non-classical hydrogen bonds were formed between the O atom from the -COOH or -OH groups of the HPA, and the H atoms from the imidazole rings of the two monomer molecules, whereas two hydrophobic interactions (p-alkyl) were observed between the aromatic ring of the HPA and the -CH 2 -group of (5) or the -CH 3 group of the EGDMA. Additionally, the interactions of the HPA with the solvent molecules were present (a p-p T-shaped interaction with a toluene ring, as well as a hydrogen bond, and a plone pair interaction with methanol). After mimicking the polymerization process, the model of the MIP5 HPA, with a template inside, was obtained. The types of interactions between the HPA and the chain model of the MIP5 HPA (just after imitated polymerization) were the same as those in the PC5 HPA (Fig. 2b) but the number of interactions varied. Here, three non-classical hydrogen bonds were formed between the O atoms from the -COOH group of the HPA and the H atoms from the imidazole rings of two monomer residues in the polymeric chain. Two hydrophobic interactions (p-alkyl) were also noted between the aromatic ring of the HPA and the -CH 3 groups from one EGDMA residue. Additionally, the HPA interacted with only one molecule of the solvent, viz. toluene, creating three hydrophobic interactions (p-p T-shaped, p-sigma and p-alkyl type). In the final step, the template molecule of the HPA was removed from the system (mimicking template removal step) and the empty cavity of the MIP5 HPA was obtained. Next, the adsorption process of the NHPA (target analyte) on the MIP5 HPA matrix was simulated (Fig. 2c). Following the MD procedure, three non-classical hydrogen bonds were formed between the NHPA (O atoms of -COOH or -NO 2 groups) and the H atoms from the imidazole rings of the two monomer residues or the -CH 3 group from the EGDMA residue. However, it should be underlined that the experimental process of adsorption was carried out in aqueous solution. Thus, the possibility of the ionization of the imidazole groups and the carboxyl group of the NHPA, was taken into account. Hence, the creation of two electrostatic interactions between the N atoms from the imidazole rings and the O atoms from the -COOor -NO 2 groups of NHPA, as well as the formation of an additional, non-classical hydrogen bond between the O atom from the -COOand the H atom from the -CH 3 group of the EGDMA residue were observed.
In the PC4 NBA system, five interactions between the NBA and the monomers or cross-linker were formed: one hydrogen bond between the H atom from the -COOH group of the NBA and the O atom of the EGDMA, two hydrophobic (p-alkyl) and two psulphur-type interactions between the aromatic ring of the NBA and the vinyl groups or S atoms of three monomers. Two interactions (p-p stacked and palkyl) between the aromatic ring of the NBA and only one solvent molecule (toluene) were also present. After the simulation of the polymerization process, only one p-sulphur-type interaction between the NBA and the S atom of the monomer residue from the polymeric chain was observed, whereas three hydrophobic interactions (two p-alkyl and one p-p stacked type) with two toluene molecules were formed. Finally, the adsorption of the NHPA on the MIP4 NBA cavity after NBA removal, revealed seven interactions between the analyte and the polymer chain. Here, five hydrogen bonds were formed between either the O or H atoms from the -COOH group or the O atom from the -NO 2 group of the NHPA and the H atoms from the -NH-or -NH 2 groups, or the S atoms of the three monomer residues, or the O atom from the EGDMA residue. Additionally, two non-classical hydrogen bonds between the O atoms from the -COOH or -NO 2 group of the NHPA and the -CH 2 -group of the EGDMA or monomer residue, and one p-sulphurtype interaction between the analyte and the monomer residue were observed.
In the analysis of the PC4 PHE system, three interactions between the template and the monomers or the EGDMA were observed, viz. two hydrophobic (palkyl) interactions and one p-sulphur-type interaction, created between the aromatic ring of the PHE and the -CH 3 groups of the EGDMA, or the S atom of (4). Additionally, two non-classical hydrogen bonds between the O atoms from the -COOH group of the PHE and the solvent molecule of methanol were formed. During the analysis of the polymeric system, with the chain of MIP4 PHE and the template, more interactions were found between the template and the polymer immediately after the simulation of the polymerization process, than in the prepolymerization complex. One classical and one non-classical hydrogen bond were created between the O atoms from the -COOH group of the PHE and the H atoms from the -NH 2 or -CH 2 -group of the two monomer residues; four hydrophobic interactions (p-alkyl) were observed between the PHE and the -CH 2 -or -CH 3 groups of the two EGDMA residues, and one psulphur-type interaction was also observed between the PHE and the (4) residue. One classical and one non-classical hydrogen bond were formed between the template and the solvent (methanol) molecules. After simulation of NHPA adsorption in the empty cavity of the MIP4 PHE , one hydrogen bond between the -NO 2 group of analyte and the H atom from the -NH 2 group of the monomer residue was found, together with five p-alkyl-type interactions of NHPA with the -CH 2 -or -CH 3 groups of the four EGDMA residues, as well as one p-sulphur-type interaction. Additionally, eight non-classical hydrogen bonds between the O atoms from the -COOH or -NO 2 and the methanol molecules were also observed.
The template molecule of the NBA in the PC1 NBA prepolymerization system created only two interactions with the EGDMA molecules, however, no interactions with the monomer were observed. One classical hydrogen bond was formed by the H atom from the -COOH group and one non-classical hydrogen bond was formed by the O atom from the -NO 2 group of the NBA. The template interacted with the solvent, forming one non-classical hydrogen bond with methanol, one p-p stacked, two p-alkyl and one p-cation-type interaction with two toluene molecules. Following polymerization, the NBA in the chain of the MIP1 NBA still had only one classical hydrogen bond with the EGDMA residue and three p-alkyl interactions with the -CH 2 -or -CH 3 groups of the EGDMA or (1) residues. Interactions of the NBA with the solvent were also present, viz. one classical hydrogen bond between the O atom from the -NO 2 group of the NBA and methanol, two non-classical hydrogen bonds between the O atoms from the -COOH or -NO 2 groups and methanol, and p-p stacked and p-alkyl interactions with the toluene molecule. After the simulation of the adsorption process of the NHPA on MIP1 NBA , two hydrogen bonds between the H atom from the -COOH group of analyte and the EGDMA residue, and between the O atom from the -NO 2 group of analytes and the monomer residue were observed. Four hydrophobic (p-alkyl) interactions of the NHPA with the -CH 3 groups of the EGDMA or monomer residues were also present, as well as one non-classical hydrogen bond, formed between the O atom from the -COOH group of the NHPA and the -CH 2 -group of the EGDMA residue. In this system, the analyte interacted with the solvent, creating three classical hydrogen bonds between the O atoms from the -COOH or the -NO 2 groups and the methanol or water molecules, and four non-classical hydrogen bonds with methanol.
In the PC1 NFE prepolymerization model, one nonclassical hydrogen bond between the O atom from the -OH group of the PHE and the -CH 2 -group of the EGDMA molecule, and two hydrophobic (palkyl) interactions of the PHE with the -CH 3 groups of the two EGDMA molecules, were observed. No interactions of the template with the monomers were found. NFE also interacted with the toluene molecule, forming p-p T-shaped, p-alkyl and p-sigma interactions. After polymerization, the NFE interacted with the polymer chain of MIP1 NFE , creating two non-classical hydrogen bonds between the O atoms of the -OH, -NO 2 and -CH 2 -groups of one EGDMA residue, and one p-alkyl interaction with the -CH 2 -group of the monomer residue. The PHE formed three interactions: p-p stacked, p-anion and pcation with the two solvent molecules of toluene. After the simulation of the NHPA adsorption process in the MIP1 NFE , three types of interactions between the analyte and polymer chain were observed. Two hydrogen bonds were formed between the O or H atoms from the -COOH group of the NHPA and the O or H atoms from the -COOH groups of the two monomer residues. Two non-classical hydrogen bonds were created between the O atom from the -OH group of the analyte and the -CH 2 -groups of one EGDMA residue. Three p-alkyl interactions were formed between the NHPA and the -CH 2 -or -CH 3 groups of monomer or EGDMA residues. The NHPA created numerous interactions with the solvent. Six classical hydrogen bonds were observed between the O atoms from the -COOH or the -NO 2 groups of the NHPA and methanol, or water molecules; two nonclassical hydrogen bonds and two p-lone pair type interactions between the analyte and the methanol molecules were also present in the analysed system.
In summary, it should be emphasized that numerous interactions between the template and the monomer, supported by an additional interaction with the cross-linker and solvent, stabilized the prepolymerization complex. It could also be noted that the analysis of interactions between the template and the polymer chain revealed slight variations, when compared to interactions in the prepolymerization complex. Finally, considering the interactions of the NHPA analyte in tested cavities, it could be concluded that the B decreased with the increasing number of interactions with solvent molecules.
On the basis of the experimental and theoretical results, the polymer MIP5 HPA was selected as the most appropriate candidate for further analysis.

Morphology and structure
In order to identify the surface morphologies of MIP5 HPA /NIP5, the analysis of micrographs from SEM was carried out. The micrographs are presented in Fig. 3.
The polymers were obtained by the thermal radical bulk polymerization. As can be seen, the micrographs of MIP5 HPA (Fig. 3a) and NIP5 (Fig. 3c) revealed a typical format of particles for bulk, crushed and sieved materials, with a high degree of irregularity, as well as heterogeneous morphology. A higher magnification of micrographs revealed negligible variations on the surface of MIP5 HPA and NIP5 (Fig. 3b, d). MIP5 HPA was characterized by a more densely formed surface, with a lesser amount of detectable macropores.
To provide an insight into the pore structure, the nitrogen sorption analysis was used, and the data obtained, were inserted into the correct mathematical models. Firstly, the nitrogen adsorption isotherms (Brunauer-Emmett-Teller isotherm) were analysed for both MIP5 HPA and NIP5 (Fig. 4) [31].
As shown, the tested materials indicated physisorption type IV isotherms that characterized mesoporous sorbents with unrestricted, mono-multilayered adsorption [32]. The hysteresis loops were also revealed with type H3 loops for MIP5 HPA and NIP5. This could indicate the slit-shaped structure of the pores of both materials obtained (Fig. 4a). Low pressure hystereses were observed for all materials, and the lack of closing points could be explained by the experimental conditions, rather than the porosity of the sorbents [33].
Next, the total specific surface area (Brunauer-Emmett-Teller isotherm) was evaluated, together with the cumulative surface area of the pores (Barrett-Joyner-Halenda model), and the volume and area of the micropores (Harkins-Jura equation) [34,35]. A slightly larger, specific surface area was noted in the case of MIP5 HPA (368.6 m 2 g -1 ) when compared with NIP5 (365.2 m 2 g -1 ). This correlated well with the surface morphology analysis, revealing negligible variations (Fig. 3). The micropore area and the external surface area of MIP5 HPA were also greater than those of NIP5 (94.9, 273.6 m 2 g -1 for MIP5 HPA and 87.5, 267.7 m 2 g -1 for NIP5, respectively). A higher micropore volume was noted for  MIP5 HPA than for NIP5 (0.039 to 0.035 cm 3 g -1 respectively) however, a higher total volume of pores was noted for NIP5, than for MIP5 HPA (0.703 to 0.635 cm 3 g -1 respectively). To provide a greater insight into the characterization of pore structures, the pore size distributions was analysed. Figure 4b presents pore size distributions, calculated using the Barrett-Joyner-Halenda model on both the adsorption and the desorption branches of the isotherms for both MIP5 HPA and NIP5.
The pore volumes, set against pore diameter plots and calculated from the desorption isotherms, revealed peaks at 33 and 23 nm for MIP5 HPA and NIP5, respectively. Thus, both materials were characterized by the mesoporous structure of pores but a higher micropore area and external surface area of MIP5 HPA, could be responsible for the imprinting effect that was observed. The phenomenon was explained by the variations in the internal pore system structures and the role of mesopores. Marć and co-workers [36] concluded that the imprinting process changed the surface morphological properties, increasing the availability of the porous structure for analyte (higher availability of particles, defined as mesopores). A high level of available mesopores was a preferable factor for sorption materials because mesopores are responsible for a particular surface area of the high affinity adsorption sites.
Subsequently, to confirm the composition of the polymer matrix, the EDS analysis was employed. This analytical technique allowed for the identification of the elemental composition of the material. Here, MIP5 HPA was analysed, since the nitrogen atoms derived only from the functional monomer of (5), should be easily detected. The EDS spectrum of MIP5 HPA is presented in Fig. 5. Nitrogen atoms were detected in MIP5 HPA in the region of 0.35 keV. Moreover, the quantitative analysis of MIP5 HPA estimated the percentage weight of nitrogen, carbon and oxygen. The results were as follows: N, 6.80 ± 0.70% wt., C, 83.27 ± 0.69% wt. and O, 9.93 ± 0.58% wt. The EDS analysis is evidence of the presence of the residues of (5) in the polymer matrix of MIP5 HPA .
A critical moment during the synthesis of MIPs relates to effective template removal. This process resulted in the formation of spatial geometries on the surface of MIPs that could be responsible for the specificity of the material. However, taking into account the EDS analysis, it is difficult to prove that the template was efficiently washed out from the polymer matrix, since the HPA molecule does not possess any other heterocyclic atoms, other than the oxygen atom. Here, the fluorescence spectroscopy could be a valuable technique for analysing various acidic metabolites of amino acids, such as 4-hydroxyphenyl lactic acid, HPA and 3,4-hydroxyphenylpropionic acid [37]. The eluates that passed through the MIP5 HPA sorbent after extraction (standard extraction in the Soxhlet apparatus without any further steps) and following the subsequent washing sequence with 1% aq. ammonium hydroxide in methanol and methanol were analysed. Figure 6 presents the three-dimensional fluorescence spectra of fractions.
As can be seen, the template removal step, carried out in the Soxhlet apparatus, had limited effectiveness. The fraction of eluates that passed through MIP5 HPA after extraction revealed the presence of leaking HPA. The detectable fluorescence was observed in the regions of k em = 280 nm and k em-= 570 nm (Fig. 6a, b). The subsequent washing sequence resulted in a significant reduction in the HPA leakage, decreasing fluorescence intensities (Fig. 6c, d). Thus, the fluorescence spectroscopy could be a valuable method in proving the completeness of the template removal step.
Next, the FT-IR analysis was employed to confirm the structure of MIP5 HPA and NIP5. Figure 7 presents the spectra.
The spectra of MIP5 HPA and NIP5 have a very similar pattern. The characteristic vibration peaks, derived from structural fragments of MIP5 HPA, could be assigned as follows: at 3454 cm -1 (broad), the -OH stretching vibrations; at 2991 and at 2959 cm -1 , the stretching vibration of the sp 2 and sp 3 bonds of the -CH from the EGDMA; at 1730 cm -1 , the -C=O stretching vibration; at 1638 cm -1 , the -C=C or -C=N stretching vibrations from the aromatic imidazole ring; at 1460 cm -1 , the asymmetric -CH deformations; at 1393 cm -1 , the -C=N stretching modes from Figure 6 Fluorescence spectra of fraction, passed through MIP5 HPA during the template removal: following Soxhlet extraction (a,b) and the subsequent washing sequence (c, d). the ring; at 1263 cm -1 and 1159 cm -1 , the C-O-C stretching vibrations at 1048 cm -1 , contribute to the C-C-C bending vibrations of the backbone chain; at 757 cm -1 , the stretching vibration of the -C=N bond which connects the imidazole ring to the vinyl moiety and the C-C-C bending vibrations of the backbone of the aliphatic chain, and at 668 cm -1 , the puckering vibration of the imidazole ring [38,39].
Finally, the TGA analyses of MIP5 HPA and NIP5 were carried out to observe the degradation of materials as a function of temperature. The TGA curves are shown in Fig. 8.
As shown, the decomposition patterns of MIP5 HPA and NIP5 are very similar (Fig. 8). The process started at around 200°C and continued up to 470°C. The decomposition step consisted of two stages with the first maximum of weight loss at 310.6°C and 310.0°C for MIP5 HPA and NIP5, respectively, with the loss of nearly 40% of the total mass. The second maximum was observed at 399.6°C and at 398.4°C for MIP5 HPA and NIP5, respectively, with an almost entire loss of the initial mass of material. We could suppose that the initial decomposition was attributed to the short chain degradation, as well as the decarboxylation process, which was also responsible for stable decomposition in the range of 300-470°C. The short maximum of weight loss at 50°C could be explained by the loss of intrinsically bound water.
To sum up, the physicochemical characterization proved the presence of monomer residues in the polymer matrix, and confirmed the different sorption capability of MIP5 HPA and NIP5.

Binding properties and selectivity
In the next step, the binding properties of MIP5 HPA and NIP5 in relation to the NHPA were evaluated. Equilibrium adsorption studies were carried out and the data obtained were inserted into two mathematical models, viz. Langmuir and Freundlich.
Firstly, the Langmuir model, transformed to the Scatchard equation, Eq. 4, was examined. The system, which fits well with the model gives a straight line on the Scatchard plot, with a slope equalling -(1/ K d ) and a y-intercept equalling B max /K d . The system which has more than one population of adsorption sites is characterized by more than one line. The Scatchard plots for MIP5 HPA and NIP5 are shown in Fig. 9a. The analysis revealed two straight lines for MIP5 HPA and only one line for NIP5. This means that two classes of adsorption sites for NHPA were predominant in MIP5 HPA, with two K d and two B max values: K d (1) = 30.2 lg L -1 and B max (1) = 153.7 lg g -1 for the higher affinity adsorption sites, and K d (2) = 357 lg L -1 and B max (2) = 416 lg g -1 for the lower affinity adsorption sites. Thus, MIP5 HPA possessed specific adsorption sites for NHPA. By contrast, NIP5 had only one low affinity adsorption site, characterized by the value of K d = 2 mg L -1 and B max = 0.9 mg g -1 .
Next, the Freundlich model, represented by Eq. 5, was used to confirm the heterogeneity of the adsorption sites. That model fits well with the adsorption data in the low concentration regions. The straight lines of log B versus log F with the regression coefficient of r 2 = 0.945 and 0.968 for MIP5 HPA and NIP5, respectively, provided evidence that adsorption could be described by the Freundlich equation (Fig. 9b). The estimated values of m for MIP5 HPA and NIP5 were 0.39 and 0.91, respectively. The results indicated that MIP5 HPA had a more heterogeneous population of adsorption sites with respect to NIP5 (the heterogeneity increased as the value of m decreased).
Finally, the selectivity of MIP5 HPA and NIP5 in relation to the target analyte, NHPA and the group of structurally related compounds were investigated. Six analogues were analysed, viz. PHE, HPA (acting as the template), NBA, NFE, 3CT and 3NT. The noncompetitive stationary binding experiments were carried out, the B were calculated according to Eq. 1, and the K d and AF according to Eq. 2 and Eq. 3, respectively. The results are shown in Fig. 10.
As shown, the B on MIP5 HPA strongly depend on the physicochemical nature of the analysed compound. Both tested amino-acids, 3-NT and 3-CT, were only slightly adsorbed (B = 0.007 and 0.168 lmol g -1 ). This could be explained by the presence of the electron-rich amine group, as well as the steric hindrance of bulky amino-acid function. A higher binding capacity was observed for NFE (B = 1.896 lmol g -1 ), the compound that possesses two functional groups, electron-withdrawing -NO 2 and electron donating -OH, substituted to the aromatic ring, which can act as a hydrogen bond donoracceptor, capable of interacting with (5) residues in the polymer network. Significantly higher B values were noted for all acidic compounds, viz. PHE, NHPA and HPA (B = 2.316, 2.361 and 2.365 lmol g -1 ). It should be emphasized that HPA was used as the template molecule in the synthesis, and NHPA was the target analyte. All those compounds are aliphatic carboxylic acids that possess the -COOH group, linked by the methylene linker with the aromatic ring. The carboxylic group could play a crucial role in the adsorption process, either in neutral or ionizable forms (the adsorption process was carried out in aqueous solution). Nevertheless, the increase in B between those carboxylic compounds were observed as follows: PHE \ NHPA \ HPA. Thus, it could be deduced that the presence of substituted functions in the aromatic ring also governed the adsorption. The highest B = 3.450 lmol g -1 was noted for NBA. This compound is an aromatic carboxylic acid, which is characterized by stronger acidity. The resonance effect and the presence of an electron withdrawing -NO 2 group in meta position, can stabilize ionizable forms of compound. This could explain the high affinity of NBA to MIP5 HPA . The B of tested compounds on NIP5 did not follow the same trend, emphasizing the specificity of MIP5 HPA . Low B of PHE, NHPA and HPA on NIP5 and high B on MIP5 HPA resulted in high specificity (AF = 2.39, 2.54 and 3.43, respectively).
To sum up, it should be stated that the B of the template HPA and the analyte NHPA were very similar in terms of MIP5 HPA but varied in relation to NIP5, emphasizing the specificity of the imprinted material in relation to NHPA. However, the highest B were noted for NBA, affecting the selectivity of MIP5 HPA . To obtain an insight into the molecular basis of the specificity and selectivity of MIP5 HPA , a theoretical analysis of the interactions between the selected analytes and the polymer, together with an evaluation of binding energies in the in silico generated cavities, were carried out.

Theoretical explanation of selectivity
The theoretical simulations of the adsorption process of different analytes, viz. NHPA, HPA, PHE, NBA, NFE, 3CT and 3NT on MIP5 HPA were performed to assess the selectivity of the MIP5 HPA cavity. The interactions of the NHPA in the cavity of MIP5 HPA, have been previously discussed.
Therefore, firstly the template molecule of the HPA was analysed (Fig. 11a). This created three interactions with the polymeric chain during adsorption: one p-p T-shaped interaction with the imidazole ring, one non-classical hydrogen interaction between the O atom from the -COOH group and the H atom from the imidazole ring of the polymeric chain, and one hydrophobic p-alkyl type interaction with the -CH 3 group of the EGDMA residue. However, considering adsorption conditions and the pH dependent equilibrium, the interactions between the ionizable forms were also considered. Here, one additional electrostatic interaction between the -COOgroup and the imidazolium cation from the polymeric chain, and one additional non-classical hydrogen bond were observed. Different types of interactions between the analyte and the solvent were also present in the analysed system: three classical and three non-classical hydrogen bonds between the H or O atoms from the -COOH or -OH groups of the HPA and the methanol molecules.
Next, the PHE adsorption on the MIP5 HPA chain was evaluated. Here, no interactions between the analyte and the polymer were observed. However, given the possibility of ionization of appropriate functional groups in the analysed system, the electrostatic interaction between the -COOanion from the analyte and the imidazolium cation from the chain, as well as that of the hydrogen bond between the O atom from the -COOgroup of the PHE and the H atom from the aromatic ring of the monomer residue were observed. Two non-classical hydrogen bonds were observed between the COOH group or aromatic ring of the PHE and the methanol molecules.
In the case of the NBA, four non-classical hydrogen bonds between the O atoms from the -COOH or the -NO 2 groups, and the H atoms from the imidazole rings or the -CH 2 -groups of the EGDMA residues in the polymeric chain, were noted (Fig. 11b). Additionally, one p-alkyl type interaction was observed between the analyte and the -CH 3 group of the EGDMA residue. Once again, taking into account the possibility of ion creation in the analysed system, three additional electrostatic interactions between the NBA and the chain were observed. Moreover, numerous interactions of the NBA with the solvent were observed. Three classical hydrogen bonds were formed with the participation of the O atoms from the -COOH and -NO 2 groups of the NBA and water, or methanol molecules. Five non-classical hydrogen bonds were created between the O atoms or the aromatic ring of the analyte and the H atoms of the water or methanol, and one p-lone pair-type interaction was formed between the aromatic ring of the NBA and the methanol.
During the theoretical simulation of the NFE adsorption process on MIP5 HPA , three non-classical hydrogen bonds, formed between the O atoms from the -NO 2 group of the NFE and the H atoms from the imidazole ring or the -CH 2 -group of the EGDMA residue, were noted. Two hydrophobic p-alkyl-type interactions were also formed between the aromatic ring of the NFE and the -CH 3 groups of the EGDMA residue, and one p-lone pair type interaction was created between the aromatic ring of the NFE and the O atom of the EGDMA residue. Considering the possibility of creation of the imidazolium cation in the polymeric chain, additional electrostatic interaction was formed between the cation and the -NO 2 group of the NFE. The analyte interacted with the solvent, creating four classical and four non-classical hydrogen bonds between the O or H atoms of the - OH and -NO 2 groups and the methanol or water molecules, and one p-sigma-type interaction between the aromatic ring of the NFE and the methanol molecule.
The analysis of 3CT adsorption on MIP5 HPA provided information relating to three non-classical hydrogen bonds, which were found between the chlorine atom or the O atom from the -COOgroup of the analyte, and the H atoms from the monomer residue ring or the -CH 2 -groups of the EGDMA residues (Fig. 11c). One interaction between the chlorine atom and the aromatic ring of the imidazole was also observed. Considering the possibility of imidazolium cation creation, no additional interactions of 3CT with the chain were observed. Here, numerous interactions of the analyte with solvent were formed: five classical and three non-classical hydrogen bonds between the H atoms from the -NH 3 ? group or the O atoms from the -OH or -COOgroups of the 3CT and methanol molecules.
Finally, the analyte of 3NT created six non-classical hydrogen bonds with a polymeric matrix between the O atoms from the -COOor -OH groups and the H atoms from the aromatic ring or alkyl chain of monomer residues, or H atoms from the -CH 2groups of the EGDMA residue (Fig. 11d). Similar to the 3CT adsorption simulation on MIP5 HPA , the creation of imidazolium ions in the polymeric chain, did not produce additional interactions. The analyte interacted with the solvent, creating three classical and six non-classical hydrogen bonds, with the participation of the H atoms from the -NH 3 ? group, or the O atoms from the -COOor -NO 2 groups and the methanol molecules.
It could be concluded that the number of electrostatic interactions correlated well with the experimental B values. The highest B of the NBA (3.450 lg g -1 ) in relation to MIP5 HPA are associated with three strong electrostatic interactions with the polymer residues, and the lowest B for 3CT and 3NT correlate with no electrostatic interactions. However, it could be also concluded that the interactions of the analyte with the solvent molecules impacted negatively on the selectivity. The highest selectivity was noted for the NHPA, the system with no interactions with solvent molecules. Additionally, mutual interactions with the monomer and cross-linker residues enhanced selectivity. By contrast, the lower selectivity was noted in the case of the NBA due to the presence of classical and non-classical hydrogen bonds, as well as the p-lone interactions with solvent molecules.

Optimization of solid phase extraction
The optimization of the SPE protocol is crucial for the characterization of recognition ability of MIPs. The appropriate solvents, used in the washing and eluting steps of the SPE procedure, should favour the development of interactions between the analyte and the monomer residues present in the cavities, but exclude the non-specific adsorption on the surface of the polymer outside the cavities. Therefore, the effect of washing and eluting solvents was investigated in this area. The following solvents were tested: water adjusted to pH 4, hexane, 15% methanol in water, 15% methanol in water adjusted to pH 4 and 15% methanol in water with 1% of ammonium hydroxide as washing solvents and t-butyl-methyl ether, ethyl acetate and methanol as eluting solvents. The results are presented in Table 1.
As shown, the most significant recovery was noted after applying water, adjusted to pH 4, or hexane as the washing solvent. The greatest recovery, obtained after washing with water, adjusted to pH 4, could be explained by the similar conditions when comparing the loading step, minimizing non-specific adsorption. The recovery obtained after washing with hexane could be explained by the fact that aprotic and nonpolar hexane could reduce the non-specific adsorption in the most effective way. All tested methanolwater washing systems resulted in lower recoveries. This means that the desorption of specifically adsorbed NHPA also occurred. The highly basic, methanol-water washing system, adjusted by 1% ammonium hydroxide, resulted in an almost total desorption of NHPA and the lowest recovery. The elution solvent was characterized by strong desorption capabilities. Here, methanol was characterized by the highest NHPA desorption percentage, probably due to its polar and protic properties.
To sum up, the impact of the washing solvent was significant for the recovery of the NHPA. Nevertheless, to simplify the separation step and shorten the time of the analysis, the following SPE procedure was selected as the most efficient: loading with 1 mL of human urine (adjusted to pH 4), washing with 1 mL of water, adjusted to pH 4, followed by 1 mL of methanol elution.

Analysis of 4-hydroxyphenylacetic acid conversion
Finally, the in vitro analysis of the conversion of HPA to NHPA in human urine was carried out to ascertain whether the analysis conditions could affect the levels of NHPA in the nitrate-rich acidified human urine sample. This analysis could also prove the applicability of the newly designed imprinted sorbent. For that purpose, the nitrate and nitrite salts were added to the sample of human urine, spiked with HPA, and adjusted to pH 2. The optimized SPE protocol was applied and the elution fraction was analysed by the UPLC-MS/MS system. The chromatogram of the human urine sample before the addition of nitrate and nitrite salts, the 4-hydroxyphenylacetic acid, adjusted to pH 2 and the chromatogram of the elution fraction after the test, are presented in Fig. 12.
As shown, the NHPA was detected in the human urine sample (peak in 6.37 min.) with a concentration of 357 ± 85 ng L -1 (n = 3). The elution fraction after the SPE protocol on the MIP5 HPA of the human urine sample, spiked with nitrate and nitrite salts and the HPA adjusted to pH 2, revealed an almost twofold increase of NHPA to 775 ± 81 ng L -1 (n = 3). The results could be explained by the conversion of HPA to NHPA in the presence of nitrate and nitrite salts, under the acidic conditions of analysis. Thus, the results could prove that the conditions of SPE could affect the concentration of the NHPA in the sample.

Conclusions
Among 20 various imprinted polymers, obtained during the design process, the polymer synthesized from 1-vinylimidazole in the presence of 4-hydroxyphenylacetic acid, revealed the highest specificity in relation to 4-hydroxy-3-nitrophenylacetic acid, a metabolite of 3-nitrotyrosine. The theoretical analysis revealed that numerous interactions between the template and monomer, supported by additional interaction with the cross-linker and the solvent, stabilized the prepolymerization complex. The physicochemical characterization proved the presence of monomer residues in the polymer matrix and confirmed the different sorption capability between the imprinted and non-imprinted polymers. The binding capacities of the template molecule of the Figure 12 Chromatogram of human urine sample before addition of nitrate and nitrite salts, 4-hydroxyphenylacetic acid and adjustment to pH 2 a and chromatogram of elution fraction, after SPE protocol of human urine sample, with added nitrate and nitrite salts, spiked with HPA and adjusted to pH 2 b. 4-hydroxyphenylacetic acid and the analyte of the 4-hydroxyphenylacetic acid, were very similar on the imprinted polymer, but varied significantly when compared to the binding capacities on the non-imprinted polymer, emphasizing the specificity of the imprinted material. The highest binding capacity, noted for 3-nitrobenzoic acid, was explained by the theoretical analysis, revealing that electrostatic interaction, together with the lack of interactions of the analyte with the solvents, supports selectivity. The optimized solid phase extraction protocol allowed for the analysis of the conversion of the 4-hydroxyphenylacetic acid to 4-hydroxy-3-nitrophenylacetic acid, in the presence of nitrate and nitrite salts, under the acidic conditions of analysis. The results proved that the conditions of the solid phase extraction process could affect the concentration of the 4-hydroxy-3-nitrophenylacetic acid in the sample.