Development of a MIP-Based QCM Sensor for Selective Detection of Penicillins in Aqueous Media

: Pharmaceuticals wastes have been recognized as emerging pollutants to the environment. Among those, antibiotics in the aquatic environment are one of the major sources of concern, as chronic, low-dose exposure can lead to antibiotic resistance. Herein, we report on molecularly imprinted polymers (MIP) to recognize penicillin V potassium salt (PenV-K), penicillin G potassium salt (PenG-K), and amoxicillin sodium salt (Amo-Na), which belong to the most widespread group of antibiotics worldwide. Characterization and optimization led to two MIPs comprising methacrylic acid as the monomer and roughly 55% ethylene glycol dimethacrylate as the crosslinker. The obtained layers led to sensitive, selective, repeatable, and reusable sensor responses on quartz crystal microbalances (QCM). The LoD for PenV-K, PenG-K, and Amo-Na sensors are 0.25 mM, 0.30 mM, and 0.28 mM, respectively; imprinting factors reach at least around three. Furthermore, the sensors displayed relative selectivity factors of up to 50% among the three penicillins, which is appreciable given their structural similarity.


Introduction
Antibiotics have revolutionized medicine during the last decades and thus changed the pattern of modern life [1,2]. Among all groups of antibiotics, penicillins belong to the most widely used ones, since they are highly effective against microbes and hardly toxic to mammals [3,4].
Penicillins are also known as β-lactam antibiotics and capable to treat the most common bacterial infections in humans and animals [3,5]. Generally, penicillin V (PenV) and penicillin G (PenG) are classified as natural penicillins, which are highly active against gram-positive bacteria. Usually, they are produced directly from fermentation of Penicillium chrysogenum [6,7]. The desire to enhance the efficiency of penicillins against gram-negative pathogens led to the aminopenicillin group of drugs through adding an amino group as a side chain to the penicillin core structure. Amoxicillin (Amo) is a common aminopenicillin. Usually, one can obtain it semi-synthetically from natural penicillin [8].
To date, it is the most prescribed antibiotic for children [9,10]. On the other side, misuse and abuse of antibiotics are becoming an increasing concern. For instance, industrial production, laboratory research, high drug consumption, and health care facilities are sources of water pollution with antibiotics. Among them, hospital sewage plays a special role regarding antibiotic contamination, because it contains various pharmaceuticals and diagnostic agents [11]. This increases both the probability that they inadvertently enter aquatic environments and thus the risk of antibiotic resistance in microbes, which result environments and thus the risk of antibiotic resistance in microbes, which result in destructive effects in a variety of ecosystems and pose a substantial threat to human health [12,13]. In addition to these health effects-affecting not only humans, but also animalsantibiotic resistance also represents a substantial economic threat since it results in extra health care costs and productivity losses, for instance, at least around €1. 5 billion each year in the EU [14].
Therefore, the World Health Organization (WHO) has declared antibiotic resistance a "major threat to public health" and demands strengthening global surveillance of this issue [15]. This includes monitoring the effectiveness of public health interventions and detecting new trends to minimize the sources of antibiotics in the environment [15,16].
A range of analytical methods has been used to detect penicillins, including mass spectrometry, chromatography, or their combination [17,18]. However, they suffer from some limitations mainly relating to sample preparation and overall complexity [19]. Hence, the need for straightforward screening measurements makes it interesting to design sensors.
Molecularly imprinted polymers (MIPs) constitute a possible strategy to obtain the artificial receptors needed for that purpose [20,21]. MIPs rely on-usually-highly crosslinked polymers that are able to bind their target species based on a key-lock mechanism [22,23]. MIPs are synthesized to contain artificial binding sites that mimic the 3D structure of a given template within the polymer bulk or on its surface [24,25]. Generally, preparing imprints in bulk polymers comprises three steps. During the first step, a functional monomer forms an adduct with the template molecule(s) in a suitable solvent through noncovalent interactions. Then, polymerization proceeds in the presence of a cross-linker and an initiator to form the polymeric matrix. Finally, the template is removed by an adequate solvent to create binding sites with selective recognition properties [26,27]. MIPs have several advantages over natural receptors such as antibodies: They combine selectivity, stability, ruggedness, and sensitivity with comparably low cost [27]. Among others, MIPs can be used as receptors on quartz-crystal microbalances (QCMs) [28]. In fact, binding of the target molecule to a binding site on the MIP increases the mass of the coating on the QCM electrode and results in decreasing resonance frequency [29].
Herein, we report a MIP thin film sensor system for the detection of Penicillin V potassium (PenV-K), Penicillin G potassium (PenG-K), and Amoxicillin sodium (Amo-Na) in aqueous media, which are the standard forms to administer these drugs to ensure solubility in body fluids; the free acids are hardly water soluble (see Figure 1 for the respective chemical structures).

Synthesis of MIPs' Thin Films
Synthesis of polymer thin films relied on free radical polymerization. Our syntheses comprised six different imprinted and non-imprinted polymers (NIPs) utilizing different functional monomer-to-cross-linker ratios. For this purpose, we applied MAA as a functional monomer and EGDMA as cross-linker, as well as PenV-K, PenG-K, and Amo-Na as the templates. Table 1 summarizes the corresponding polymer compositions. To synthesize MIP thin films, we first dissolved the functional monomer and 8 mg (0.028 mmol) of the respective template in 800 µL DMSO as the porogenic solvent in an Eppendorf vial and sonicated for about 15 min. In the next step, we added 4 mg (0.024 mmol) AIBN as the initiator and the respective amount of cross-linker according to Table 1. Thereupon, the monomer solution was kept in an ice bath while flushing Argon through the vial for about 20 min to remove oxygen. Subsequently, this mixture was thermally polymerized while stirring with a magnetic stirrer at 60 • C for around 30 min just before reaching the gel point. This oligomer solution then served for coating the respective QCM electrode. The corresponding non-imprinted (NIP) sample was prepared in the same manner but without adding the template. Figure 2 shows the dual-electrode geometry of QCMs with a fundamental resonance frequency of 10 MHz manufactured by screen printing onto commercially available AT-cut quartz plates (168 µm thick, 13.8 mm diameter; purchased from Roditi Inc., London, UK) as described before [30]. Briefly, the first step comprised of preparing the necessary stencils from artificial silk, which we strained and glued onto a metal frame. Next, the textile was covered with photoresist lacquer (Azocol Poly-Plus S, KIWO, Wiesloch, Germany) and dried in the dark for around 1 hour. Last, the stencil was hardened under UV light (UV-Belichtungsgeraet, 220 Volt, isel, Eichenzell, Germany) while exposing it through a mask that shadowed off the further electrode design. After illuminating for around 1 minute and washing with water to remove the non-hardened photoresist, the mask was ready to use for screen printing. Two different sieves were prepared for the air side and sample side (fully metalized side) of QCMs, respectively. During the final step, the quartz was fixed on a disc-shaped Teflon socket by using a vacuum pump. Then, the sieve, comprising the respective pattern, was placed over the QCMs followed by screen printing commercially available brilliant gold paste (Heraeus GmbH, Hanau, Germany) through the stencil. Finally, we heated QCMs in the oven at 400 • C for 4 h to remove organic residues and, thus, to generate the final gold electrodes. Consecutively, the same procedure was repeated to deposit the gold electrode patterns on the second face of the corresponding quartzes. side and sample side (fully metalized side) of QCMs, respectively. During the fi the quartz was fixed on a disc-shaped Teflon socket by using a vacuum pump. T sieve, comprising the respective pattern, was placed over the QCMs followed b printing commercially available brilliant gold paste (Heraeus GmbH, Hanau, G through the stencil. Finally, we heated QCMs in the oven at 400 °C for 4 h to organic residues and, thus, to generate the final gold electrodes. Consecutively, procedure was repeated to deposit the gold electrode patterns on the second fa corresponding quartzes.

Coating Process of Thin Films on QCM
For spin-coating, we used a PI-KEM LTD G3P-8 spin coater. To generate th we dropped 5 µL of oligomer solution containing the template and without a respectively, on one electrode of a dual-electrode QCM by a piston-pipet at a sp of 2000 rpm and kept the QCM at this speed for about 15 s to establish a homo sensor layer. One of the electrodes served as a working electrode to be coated respective MIP; the other one was the reference electrode coated with NIP. In all c first coated one electrode, hardened the polymer, and then coated the other electr terwards, the sensor was incubated at 80 °C overnight to evaporate the DMSO and the polymer. Finally, to remove the template molecules from the polymer ma QCMs were rinsed with distilled water in a stirred beaker overnight at room temp To determine the layer thickness of the polymer films on the QCM, we ut Agilent ENA series network analyzer (model no. E5062A). The change of frequen damping values before/after coating and after washing the polymer was recor allowed us to estimate the respective layers' heights based on the experience of th [31]. In addition, removal of the template molecules from the MIP film was also m by the network analyzer: Successful removal resulted in an increase of the reson quency of QCM after washing.

AFM Analysis
The surface morphologies of imprinted and non-imprinted polymer thin fi washing and template removal, respectively, were examined by atomic force mic (AFM) at ambient temperature. These measurements took place on a Bruker Inst NanoScope VIII with a TESPA-V2 tip in tapping mode (TM) at a scan rate of 1 H scan size of 5 µm, using the AFM operating and evaluation software provided by ducer. Figure 3 shows the schematic of the used QCM measuring setup. To ensu mass-sensitive sensor measurements, it utilizes dual-electrode QCMs to compe unspecific frequency effects (such as temperature changes and nonspecific through the non-imprinted film on the reference electrode. The electrode showin resonant frequency was used for MIP deposition to prevent crosstalk of the two c since frequency decreases on the MIP channel were expected to be larger than the NIP channel.

Coating Process of Thin Films on QCM
For spin-coating, we used a PI-KEM LTD G3P-8 spin coater. To generate thin films, we dropped 5 µL of oligomer solution containing the template and without adding it, respectively, on one electrode of a dual-electrode QCM by a piston-pipet at a spin speed of 2000 rpm and kept the QCM at this speed for about 15 s to establish a homogeneous sensor layer. One of the electrodes served as a working electrode to be coated with the respective MIP; the other one was the reference electrode coated with NIP. In all cases, we first coated one electrode, hardened the polymer, and then coated the other electrode. Afterwards, the sensor was incubated at 80 • C overnight to evaporate the DMSO and harden the polymer. Finally, to remove the template molecules from the polymer matrix, the QCMs were rinsed with distilled water in a stirred beaker overnight at room temperature.
To determine the layer thickness of the polymer films on the QCM, we utilized an Agilent ENA series network analyzer (model no. E5062A). The change of frequencies and damping values before/after coating and after washing the polymer was recorded and allowed us to estimate the respective layers' heights based on the experience of the group [31]. In addition, removal of the template molecules from the MIP film was also monitored by the network analyzer: Successful removal resulted in an increase of the resonance frequency of QCM after washing.

AFM Analysis
The surface morphologies of imprinted and non-imprinted polymer thin films after washing and template removal, respectively, were examined by atomic force microscopy (AFM) at ambient temperature. These measurements took place on a Bruker Instruments NanoScope VIII with a TESPA-V2 tip in tapping mode (TM) at a scan rate of 1 Hz and a scan size of 5 µm, using the AFM operating and evaluation software provided by the producer. Figure 3 shows the schematic of the used QCM measuring setup. To ensure actual mass-sensitive sensor measurements, it utilizes dual-electrode QCMs to compensate for unspecific frequency effects (such as temperature changes and nonspecific binding) through the non-imprinted film on the reference electrode. The electrode showing lower resonant frequency was used for MIP deposition to prevent crosstalk of the two channels, since frequency decreases on the MIP channel were expected to be larger than those on the NIP channel.

QCM Measurement Setup
around 30-60 s to avoid spikes resulting from sample changes). After the frequency reached a stable value or established equilibrium, we flushed the cell three times with 150 µL of 50 mM KCl solution to remove analyte molecules from the sensor layer. Herein, the QCM measurements served to characterize PenV-K, PenG-K, and Amo-Na MIP thin films according to their sensitivity and selectivity of sensor performance. Using KCl solution as the background rather than water helped us to exclude/compensate for any effects caused by changing ion concentrations in the solution.

Optimization of MIP Thin Film
The first step in designing MIP sensors is to optimize polymerization conditions: Figure 4 summarizes the sensor effects of the six polymers mentioned in Table 1 when exposing them to aqueous solutions containing 2.50 mM of PenV-K and PenG-K, respectively. "Sensor response" in this case refers to the signal difference between MIP channel and NIP channel. It is worth mentioning that all sensor responses were normalized to the Briefly, the measurement setup comprised a custom-made polymethylmethacrylate (PMMA) measuring cell with a polydimethylsiloxane (PDMS) cap to fix the QCM sensor inside the cell. The overall sample volume of this cell was around 130 µL. The PDMS contained inlet and outlet tubes to expose the sensor to sample solutions. Furthermore, due to its comparably low thermal conductance, it usually does not influence the temperature in the measuring cell within the time of sample exposure. The cell connects to a custommade oscillator via SMA sockets and coaxial cables. The oscillator uses the QCM as the frequency-determining element. To read out the respective resonance frequencies of the two channels as a function of time and transfer the output data to the PC, we used a LabVIEW routine prepared in house.
Every measurement started with injecting 150 µL of 50 mM KCl solution into the measuring cell using a piston-pipette and waiting until the frequency showed stable baseline at room temperature. Then, we recorded the frequency changes resulting from injecting analyte solutions in stopped flow (during each injection the read-out was paused for around 30-60 s to avoid spikes resulting from sample changes). After the frequency reached a stable value or established equilibrium, we flushed the cell three times with 150 µL of 50 mM KCl solution to remove analyte molecules from the sensor layer. Herein, the QCM measurements served to characterize PenV-K, PenG-K, and Amo-Na MIP thin films according to their sensitivity and selectivity of sensor performance. Using KCl solution as the background rather than water helped us to exclude/compensate for any effects caused by changing ion concentrations in the solution.

Optimization of MIP Thin Film
The first step in designing MIP sensors is to optimize polymerization conditions: Figure 4 summarizes the sensor effects of the six polymers mentioned in Table 1 when exposing them to aqueous solutions containing 2.50 mM of PenV-K and PenG-K, respectively. "Sensor response" in this case refers to the signal difference between MIP channel and NIP channel. It is worth mentioning that all sensor responses were normalized to the respective layer heights because the sensitivity of sensors coated with bulk-imprinted materials varies with thickness. number of available binding sites obviously decreases and, thus, also the respective senso responses. Obviously, both binding capacity and selectivity of MIPs have been optimize by finding the ideal ratio between the functional monomer(s), crosslinker, and templa in a suitable solvent [32]. Overall, polymer 4 turned out optimal for binding PenV-K mo ecules and polymer 5 for PenG-K, even though the differences between them are withi the error margins. Obviously, hydrophilic (hydrogen bond and polar) interactions play a key role i sensor recognition: While the PenV-K molecule provides 1 hydrogen donor and 6 hydro gen acceptors [33], PenG molecules comprises 1 hydrogen donor and 5 hydrogen accep tors [34], which may be the reason that the slightly more cross-linked polymer yield somewhat higher sensor responses for PenG-K. Therefore, further experiments relied o those two polymers. Additionally, polymer 4 was also selected to establish Amo-Na MI and NIP layers. It is worth noting the amoxicillin molecule also offers three hydroge donors and seven hydrogen acceptors [35], which allowed us to use the lower cross-linke polymer recipe. Figure 5 shows the estimated polymer film thicknesses of NIPs and MIPs layer base on the Sauerbrey equation for QCMs (Equation (1)) [36] after coating polymer and afte washing/template removal. To calculate those layer heights from frequency shifts, on needs to use the well-known Sauerbrey equation: Here, Δf is the measured frequency change (Hz), f0 the fundamental frequency o QCM (10 × 10 6 Hz), Δm is the mass change (g), A is the electrode area (0.125 cm 2 ), ρ is th density of quartz (2.643 g/cm 3 ), and µ is the shear module of quartz (2.947 × 10 11 g cm −1 s −2 Inserting the respective values into Equation (1) and assuming a rigid polymer layer wit Clearly, polymer 1 and polymer 2 show only low affinity to the respective template molecule due to the low monomer-to-template ratio. There is almost no difference in sensor signal between the MIP and the NIP sensors, respectively, which indicates only nonspecific binding of the penicillin molecules to the sensor layers. Polymer 3, polymer 4, and polymer 5 exhibit larger imprinting effects against the templates by factors of five to eight, respectively. Polymer 6 containing a higher amount of cross-linker again exhibits lower sensor responses: When increasing the amount of cross-linker above the optimum, the number of available binding sites obviously decreases and, thus, also the respective sensor responses. Obviously, both binding capacity and selectivity of MIPs have been optimized by finding the ideal ratio between the functional monomer(s), crosslinker, and template in a suitable solvent [32]. Overall, polymer 4 turned out optimal for binding PenV-K molecules and polymer 5 for PenG-K, even though the differences between them are within the error margins.
Obviously, hydrophilic (hydrogen bond and polar) interactions play a key role in sensor recognition: While the PenV-K molecule provides 1 hydrogen donor and 6 hydrogen acceptors [33], PenG molecules comprises 1 hydrogen donor and 5 hydrogen acceptors [34], which may be the reason that the slightly more cross-linked polymer yields somewhat higher sensor responses for PenG-K. Therefore, further experiments relied on those two polymers. Additionally, polymer 4 was also selected to establish Amo-Na MIP and NIP layers. It is worth noting the amoxicillin molecule also offers three hydrogen donors and seven hydrogen acceptors [35], which allowed us to use the lower cross-linked polymer recipe. Figure 5 shows the estimated polymer film thicknesses of NIPs and MIPs layer based on the Sauerbrey equation for QCMs (Equation (1)) [36] after coating polymer and after washing/template removal. To calculate those layer heights from frequency shifts, one needs to use the well-known Sauerbrey equation: ymer film will hinder access to the binding sites in the polymer bulk, leading to increasin damping of the electrodes and, thus, making them useless for sensor measurements. Here for example, the MIP PenV-K film was approximately 285 nm thick. After rinsing an template removal, it reached 195 nm (which shows that washing also removed part of th matrix), while the NIP PenV-K film after polymer coating was about 230 nm, and, afte washing, the polymer film revealed a thickness of 210 nm.  Figure 6 shows tapping mode AMF images of the surface topographies of washe PenV-K imprinted and non-imprinted coated thin films on QCM electrodes, respectively It can be seen clearly that both polymers have homogenous surfaces, but with differen textures. The NIP surface in Figure 6a is comparably smooth and flat. The z-axis, in thi case, was 135 nm, which indicates that the wave-like structures are roughly ±65 nm high but most probably are an artifact in the measurement. The immediate roughness of th polymer is much lower. In contrast, the morphology of PenV-K MIP after rinsing an template removal in Figure 6b shows similar structures with lateral dimensions in th micrometer range. However, the surface is evidently much rougher, than that of the NIP Obviously, there are cavities on the MIP surface, which are not present in the NIP surface Of course, these do not represent the imprints as such: This is impossible given the size of the template molecules, since the lateral resolution of AFM depends on tip aperture which was around 5-10 nm for this setup. However, it still evidently demonstrates tha PenV-K molecules strongly influence the morphology of the polymer. Removing them from the polymer matrix affects the overall topography of the corresponding polymer. Here, ∆f is the measured frequency change (Hz), f 0 the fundamental frequency of QCM (10 × 10 6 Hz), ∆m is the mass change (g), A is the electrode area (0.125 cm 2 ), ρ is the density of quartz (2.643 g/cm 3 ), and µ is the shear module of quartz (2.947 × 10 11 g cm −1 s −2 ). Inserting the respective values into Equation (1) and assuming a rigid polymer layer with density ρ = 1 g/cm 3 leads to the conclusion that a frequency shift of ∆f = −1 kHz corresponds to a 40-nm polymer layer height.

Polymer Surface Analysis
This aspect is important in so far as bulk-polymerized MIP thin films need to provide sufficient analyte diffusion into the "inner" parts of the polymer matrix [37]. Hence, increasing the layer height usually enhances the sensor response. Nevertheless, a thick polymer film will hinder access to the binding sites in the polymer bulk, leading to increasing damping of the electrodes and, thus, making them useless for sensor measurements. Here, for example, the MIP PenV-K film was approximately 285 nm thick. After rinsing and template removal, it reached 195 nm (which shows that washing also removed part of the matrix), while the NIP PenV-K film after polymer coating was about 230 nm, and, after washing, the polymer film revealed a thickness of 210 nm. Figure 6 shows tapping mode AMF images of the surface topographies of washed PenV-K imprinted and non-imprinted coated thin films on QCM electrodes, respectively. It can be seen clearly that both polymers have homogenous surfaces, but with different textures. The NIP surface in Figure 6a is comparably smooth and flat. The z-axis, in this case, was 135 nm, which indicates that the wave-like structures are roughly ±65 nm high, but most probably are an artifact in the measurement. The immediate roughness of the polymer is much lower. In contrast, the morphology of PenV-K MIP after rinsing and template removal in Figure 6b shows similar structures with lateral dimensions in the micrometer range. However, the surface is evidently much rougher, than that of the NIP. Obviously, there are cavities on the MIP surface, which are not present in the NIP surface. Of course, these do not represent the imprints as such: This is impossible given the sizes of the template molecules, since the lateral resolution of AFM depends on tip aperture, which was around 5-10 nm for this setup. However, it still evidently demonstrates that PenV-K molecules strongly influence the morphology of the polymer. Removing them from the polymer matrix affects the overall topography of the corresponding polymer.   Table 1. First, the data reveals concentration-dependent sensor signals. Second, MIP yield substantially higher frequency shifts at all concentrations compared to the corresponding NIP: For instance, injection of PenV-K solution at c = 2.50 mM leads to frequency decreasing by −286 Hz on the MIP side and −132 Hz on the NIP side, respectively. This corresponds to a −154-Hz difference signal and, thus, an imprinting factor of 2.2 (the imprinting factor is calculated by dividing the sensor signal of the MIP through the one of the NIP). Additionally, the signal of NIP shows some inherent affinity of the polymer to PenV-K, which also can be desirable to achieve MIP with high affinity [38]. Furthermore, the PenV-K-MIP leads to linear sensor characteristic within the concentration range of 0.30-2.00 mM PenV-K with a correlation coefficient of R 2 = 0.98. Furthermore, the limit of detection (LOD) is calculated based on (3 S/N) (with QCM baseline noise level standing at around −10 Hz). Based on this, the PenV-K MIP sensor offers LOD = 0.25 mM. In addition, repeatability is an important parameter for any chemical sensor to ensure analytical quality. Figure 8 gives the QCM responses obtained for three consecutive exposures with the same sample and washing steps in between them. Here, in terms of quantitative detection, the results are promising: All frequency responses are fully reversible after washing the sensors with 50 mM KCl solution. Imprinting factors here are around   Table 1. First, the data reveals concentration-dependent sensor signals. Second, MIP yield substantially higher frequency shifts at all concentrations compared to the corresponding NIP: For instance, injection of PenV-K solution at c = 2.50 mM leads to frequency decreasing by −286 Hz on the MIP side and −132 Hz on the NIP side, respectively. This corresponds to a −154-Hz difference signal and, thus, an imprinting factor of 2.2 (the imprinting factor is calculated by dividing the sensor signal of the MIP through the one of the NIP). Additionally, the signal of NIP shows some inherent affinity of the polymer to PenV-K, which also can be desirable to achieve MIP with high affinity [38]. Furthermore, the PenV-K-MIP leads to linear sensor characteristic within the concentration range of 0.30-2.00 mM PenV-K with a correlation coefficient of R 2 = 0.98. Furthermore, the limit of detection (LOD) is calculated based on (3 S/N) (with QCM baseline noise level standing at around −10 Hz). Based on this, the PenV-K MIP sensor offers LOD = 0.25 mM.   Table 1. First, the data reveals concentration-dependent sensor signals. Second, MIP yield substantially higher frequency shifts at all concentrations compared to the corresponding NIP: For instance, injection of PenV-K solution at c = 2.50 mM leads to frequency decreasing by −286 Hz on the MIP side and −132 Hz on the NIP side, respectively. This corresponds to a −154-Hz difference signal and, thus, an imprinting factor of 2.2 (the imprinting factor is calculated by dividing the sensor signal of the MIP through the one of the NIP). Additionally, the signal of NIP shows some inherent affinity of the polymer to PenV-K, which also can be desirable to achieve MIP with high affinity [38]. Furthermore, the PenV-K-MIP leads to linear sensor characteristic within the concentration range of 0.30-2.00 mM PenV-K with a correlation coefficient of R 2 = 0.98. Furthermore, the limit of detection (LOD) is calculated based on (3 S/N) (with QCM baseline noise level standing at around −10 Hz). Based on this, the PenV-K MIP sensor offers LOD = 0.25 mM. In addition, repeatability is an important parameter for any chemical sensor to ensure analytical quality. Figure 8 gives the QCM responses obtained for three consecutive exposures with the same sample and washing steps in between them. Here, in terms of quantitative detection, the results are promising: All frequency responses are fully reversible after washing the sensors with 50 mM KCl solution. Imprinting factors here are around In addition, repeatability is an important parameter for any chemical sensor to ensure analytical quality. Figure 8 gives the QCM responses obtained for three consecutive exposures with the same sample and washing steps in between them. Here, in terms of quantitative detection, the results are promising: All frequency responses are fully reversible after washing the sensors with 50 mM KCl solution. Imprinting factors here are [%] = × (2) One can see that the sensing results are in a similar range, and CV values are 7%, 12%, and 9% against PenV-K solutions at c = 2.50 mM, 1.25 mM, and 0.60 mM, respectively.  . As previously mentioned, the PenG-K MIP composition slightly differs from that of PenV-K (see Table 1: PenV-K: polymer 4; PenG-K: polymer 5). Similar to PenV-K sensors, the frequency responses of the sensors decreases linearly as PenG-K concentration increases both on the MIP and NIP channels. One can expect such a decrease: The more concentrated the solution, the more molecules are available to bind to the imprinted cavities present at the polymer matrix. For instance, at 2.50 mM PenG-K, the MIP exhibits a frequency shift of −235 Hz while the NIP signal reaches only around −55 Hz. Thus, frequency shifts on MIPcoated electrodes against NIP-coated ones reveal more than four times higher sensitivity, strongly indicative of successful imprinting and, thus, the incorporation of PenG-K molecules into recognition sites. Additionally, the outcomes indicate that imprinting turned out to be successful with a linear correlation of sensor responses, with LOD = 0.30 mM based on the signal to the noise level of the system.  Figure 10 shows the repeatability evaluation of PenG-K sensor response at the concentrations of 2.50 mM, 2.00 mM, and 1.25 mM, respectively. For that purpose, we injected One can see that the sensing results are in a similar range, and CV values are 7%, 12%, and 9% against PenV-K solutions at c = 2.50 mM, 1.25 mM, and 0.60 mM, respectively. Figure 9 exhibits sensor characteristics results of PenG-K sensor against five different concentrations (2.50 mM, 2.00 mM, 1.25 mM, 0.60 mM, and 0.50 mM). As previously mentioned, the PenG-K MIP composition slightly differs from that of PenV-K (see Table 1: PenV-K: polymer 4; PenG-K: polymer 5). Similar to PenV-K sensors, the frequency responses of the sensors decreases linearly as PenG-K concentration increases both on the MIP and NIP channels. One can expect such a decrease: The more concentrated the solution, the more molecules are available to bind to the imprinted cavities present at the polymer matrix. For instance, at 2.50 mM PenG-K, the MIP exhibits a frequency shift of −235 Hz while the NIP signal reaches only around −55 Hz. Thus, frequency shifts on MIP-coated electrodes against NIP-coated ones reveal more than four times higher sensitivity, strongly indicative of successful imprinting and, thus, the incorporation of PenG-K molecules into recognition sites. Additionally, the outcomes indicate that imprinting turned out to be successful with a linear correlation of sensor responses, with LOD = 0.30 mM based on the signal to the noise level of the system. [%] = × (2) One can see that the sensing results are in a similar range, and CV values are 7%, 12%, and 9% against PenV-K solutions at c = 2.50 mM, 1.25 mM, and 0.60 mM, respectively.  . As previously mentioned, the PenG-K MIP composition slightly differs from that of PenV-K (see Table 1: PenV-K: polymer 4; PenG-K: polymer 5). Similar to PenV-K sensors, the frequency responses of the sensors decreases linearly as PenG-K concentration increases both on the MIP and NIP channels. One can expect such a decrease: The more concentrated the solution, the more molecules are available to bind to the imprinted cavities present at the polymer matrix. For instance, at 2.50 mM PenG-K, the MIP exhibits a frequency shift of −235 Hz while the NIP signal reaches only around −55 Hz. Thus, frequency shifts on MIPcoated electrodes against NIP-coated ones reveal more than four times higher sensitivity, strongly indicative of successful imprinting and, thus, the incorporation of PenG-K molecules into recognition sites. Additionally, the outcomes indicate that imprinting turned out to be successful with a linear correlation of sensor responses, with LOD = 0.30 mM based on the signal to the noise level of the system.  Figure 10 shows the repeatability evaluation of PenG-K sensor response at the concentrations of 2.50 mM, 2.00 mM, and 1.25 mM, respectively. For that purpose, we injected  Figure 10 shows the repeatability evaluation of PenG-K sensor response at the concentrations of 2.50 mM, 2.00 mM, and 1.25 mM, respectively. For that purpose, we injected standards at the same concentration three times before shifting to the next concentration. The data revealed CV values of 9%, 15%, and 11%, respectively. Chemosensors 2021, 9, x FOR PEER REVIEW 10 of 17 standards at the same concentration three times before shifting to the next concentration. The data revealed CV values of 9%, 15%, and 11%, respectively.  Table 2 summarizes the average sensor responses on the MIP and NIP electrodes, respectively, of the Amo-Na sensor, derived from three QCM measurements in the concentration range of 0.30-2.50 mM each. Like the MIP sensors targeting the unmodified penicillins, also in this case, the MIP exhibit significant binding affinity against its target. Nonspecific interaction with the NIP layers is lower by a factor of three. For example, the average frequency shift on the MIP channel when exposing Amo-Na sensor to a solution at c = 2.00 mM analyte reaches −265 Hz, while the NIP channel shows a signal around −63 Hz, equal to around −200-Hz signal difference. Likewise, Figure 11 presents the Amo-Na MIP and NIP sensor characteristics for the concentration series of 0.30-2.50 mM. The correlation coefficients were calculated for MIP to R 2 = 0.98 and R 2 = 0.96 on the NIP side, which indicates the analytical feasibility of the sensor in the concentration range demonstrated. The resulting sensor characteristic shows LOD = 0.28 mM, which is very similar to MIP PenV-K senor. This is even more appreciable since both MIPs rely on the same polymer system, which is rather unusual: Switching between even closely related templates usually makes it necessary to re-optimize the polymer, even if this requires only minor amendments [39].  Table 2 summarizes the average sensor responses on the MIP and NIP electrodes, respectively, of the Amo-Na sensor, derived from three QCM measurements in the concentration range of 0.30-2.50 mM each. Like the MIP sensors targeting the unmodified penicillins, also in this case, the MIP exhibit significant binding affinity against its target. Nonspecific interaction with the NIP layers is lower by a factor of three. For example, the average frequency shift on the MIP channel when exposing Amo-Na sensor to a solution at c = 2.00 mM analyte reaches −265 Hz, while the NIP channel shows a signal around −63 Hz, equal to around −200-Hz signal difference. Likewise, Figure 11 presents the Amo-Na MIP and NIP sensor characteristics for the concentration series of 0.30-2.50 mM. The correlation coefficients were calculated for MIP to R 2 = 0.98 and R 2 = 0.96 on the NIP side, which indicates the analytical feasibility of the sensor in the concentration range demonstrated. The resulting sensor characteristic shows LOD = 0.28 mM, which is very similar to MIP PenV-K senor. This is even more appreciable since both MIPs rely on the same polymer system, which is rather unusual: Switching between even closely related templates usually makes it necessary to re-optimize the polymer, even if this requires only minor amendments [39]. Figure 11. Sensor characteristics were obtained for MIP and NIP sensors against Amo-Na.

Sensing Thermodynamics
Equation (3) shows the feasible recognition reaction of the target analyte (PenV-K, PenG-K, and Amo-Na) interaction within synthesized MIP binding sites.

MIP
Analyte ⇌ complex Hence, in such a system at equilibrium, Equation (4) expresses the association constant (Ka).

= [ ]
Here, CComplex is the concentration of occupied MIP binding sites within the film volume, CMIP is the concentration of unoccupied MIP binding sites after binding, and CAnalyte is the analyte concentration in solution.
By applying the equation for the Langmuir adsorption isotherm and setting the fractional occupancy (θ) as a function of frequency changes for such sensing approaches on QCM (based on the Sauerbrey equation), one obtains Equation (5) [40].
Equation (5) can be rearranged to Equation (6) [41]. Thus, plotting the frequency shifts obtained for the respective penicillin standards against ΔfAnalyte/CAnalyte revealed Ka and the frequency shift corresponding to the maximum number of binding sites in the MIP (Δ ): the slope of the curve corresponded to 1/Ka and the y-intercept to ΔfMIP (Max).

Sensing Thermodynamics
Equation (3) shows the feasible recognition reaction of the target analyte (PenV-K, PenG-K, and Amo-Na) interaction within synthesized MIP binding sites.
Hence, in such a system at equilibrium, Equation (4) expresses the association constant (K a ).
here, C Complex is the concentration of occupied MIP binding sites within the film volume, C MIP is the concentration of unoccupied MIP binding sites after binding, and C Analyte is the analyte concentration in solution.
By applying the equation for the Langmuir adsorption isotherm and setting the fractional occupancy (θ) as a function of frequency changes for such sensing approaches on QCM (based on the Sauerbrey equation), one obtains Equation (5) [40].
Equation (5) can be rearranged to Equation (6) [41]. Thus, plotting the frequency shifts obtained for the respective penicillin standards against ∆f Analyte /C Analyte revealed K a and the frequency shift corresponding to the maximum number of binding sites in the MIP (∆f max MIP ): the slope of the curve corresponded to 1/K a and the y-intercept to ∆f MIP (Max) .  All three MIPs show negative Gibbs energies for binding that is in the same range, namely, around −13.5 kJ/mol. This is well within the range that is expected for reversible sensor signals. As a rule of thumb, ∆G ≥ −25 kJ/mol means reversible signals. The respective equilibrium constant is K = 10 4 .
The values obtained, presented in Table 3, are in line with the literature: Previous research on acrylic-based MIPs binding molecules in roughly similar size range led to ∆G = −11.70 kJ/mol for carbofuran and ∆G a = −5.05 kJ/mol for profenofos, respectively [30]. In polypyrroles, caffeine and theophylline yielded Gibbs energies of ∆G = −10.81 kJ/mol and ∆G = −6.6 kJ/mol, respectively [42]. The energy values indicate that analyte binding resulted from hydrogen bonds or electrostatic interactions.

Selectivity
Of course, selectivity is a key issue for any sensor and preferably each MIP should specifically bind its own analyte. Figure 12 shows the QCM response curves of the selectivity test of Amo-Na sensor against PenG-K and PenV-K at c = 2.50 Mm (50 mM KCl used as background) for both MIP and NIP channels. Firstly, it reveals the higher binding affinity of the sensor to bind its templated analyte (Amo-Na) towards analogous penicillins. Next, the frequency signals on MIP channel decreased by −290 Hz, −175 Hz, and -160 Hz for Amo-Na, PenG-K, and PenV-K, respectively. In addition, the NIP channel showed the same sensor response of around −90 Hz for all injected penicillin (at c = 2.50 Mm). The latter demonstrates that the inherent affinity of the polymer is the same for all three penicillins. This is very noticeable also in terms of proving the imprinted approach: Templating the material indeed leads to substantially different signals for different target molecules. All three MIPs show negative Gibbs energies for binding that is in the same range, namely, around −13.5 kJ/mol. This is well within the range that is expected for reversible sensor signals. As a rule of thumb, ΔG ≥ −25 kJ/mol means reversible signals. The respective equilibrium constant is K = 10 4 .

Type of Sensor
The values obtained, presented in Table 3, are in line with the literature: Previous research on acrylic-based MIPs binding molecules in roughly similar size range led to ΔG = −11.70 kJ/mol for carbofuran and ΔGa = −5.05 kJ/mol for profenofos, respectively [30]. In polypyrroles, caffeine and theophylline yielded Gibbs energies of ΔG = −10.81 kJ/mol and ΔG = −6.6 kJ/mol, respectively [42]. The energy values indicate that analyte binding resulted from hydrogen bonds or electrostatic interactions.

Selectivity
Of course, selectivity is a key issue for any sensor and preferably each MIP should specifically bind its own analyte. Figure 12 shows the QCM response curves of the selectivity test of Amo-Na sensor against PenG-K and PenV-K at c = 2.50 Mm (50 mM KCl used as background) for both MIP and NIP channels. Firstly, it reveals the higher binding affinity of the sensor to bind its templated analyte (Amo-Na) towards analogous penicillins. Next, the frequency signals on MIP channel decreased by −290 Hz, −175 Hz, and -160 Hz for Amo-Na, PenG-K, and PenV-K, respectively. In addition, the NIP channel showed the same sensor response of around −90 Hz for all injected penicillin (at c = 2.50 Mm). The latter demonstrates that the inherent affinity of the polymer is the same for all three penicillins. This is very noticeable also in terms of proving the imprinted approach: Templating the material indeed leads to substantially different signals for different target molecules. Figure 12. Selectivity measurement of Amo-Na Sensor against PenG-K and PenV-K at c = 2.50 mM, respectively. Figure 12. Selectivity measurement of Amo-Na Sensor against PenG-K and PenV-K at c = 2.50 mM, respectively. Furthermore, Figure 13 summarizes the selectivity pattern of PenV-K, PenG-K, and Amo-Na sensors as the relative sensor effect towards each other at c = 2.50 mM. To better interpret the selectivity results, Figure 13 uses the relative selectivity factor (%β) defined by Equation (7), which is the quotient of the net sensor responses (i.e., the differences in frequency shifts between MIP and NIP channels for the two competing species).

β[%] = 100 ×
Sensor e f f ect (competitor) Sensor e f f ect (template) (7) (see Figure 12). This amino-penicillin contains an additional amino group and hydroxyl-group in its structure, which leads to substantial differences in selectivit versus PenV-K and PenG-K. Therefore, Amo-Na sensor demonstrate the highest ity factors in comparison with sensors targeting the natural penicillins by achievi ues of around 52% and 49% against PenV-K and PenG-K, respectively. Again, it mentioning that the core of the respective structures remained constant. Figure 13. Relative sensor responses for PenV-K, PenG-K, and Amo-Na sensors, respec ward the target penicillins at c = 2.50 mM.

Comparison of Data with Previous Studies
Though the analytical literature on detecting penicillins is, of course, extens a limited number of papers has been published so far on designing chemical se such compounds. Table 4 summarizes biomimetic/polymer-based approaches purpose, including three examples of non-penicillin antibiotics in the last three the best of our knowledge, there is no literature on Penicillin V yet. Accordingly, the results reveal that MIPs for PenV-K and PenG-K are selective to their respective compounds, but not very strongly. Therefore, the relative selectivity factors were recorded at 77% and 83% for the two chemicals against each other. However, even these small selectivity numbers are remarkable, given how similar the molecular structures of these compounds are: The two natural penicillins differ by an oxygen atom, which only leads to a slight increase of the side chain length of PenV-K against PenG-K molecules while the core structure is identical, as shown in Figure 1.
In contrast, amoxicillin show higher selectivity, both as a template and as an analyte (see Figure 12). This amino-penicillin contains an additional amino group and an extra hydroxyl-group in its structure, which leads to substantial differences in selectivity factors versus PenV-K and PenG-K. Therefore, Amo-Na sensor demonstrate the highest selectivity factors in comparison with sensors targeting the natural penicillins by achieving β values of around 52% and 49% against PenV-K and PenG-K, respectively. Again, it is worth mentioning that the core of the respective structures remained constant.

Comparison of Data with Previous Studies
Though the analytical literature on detecting penicillins is, of course, extensive, only a limited number of papers has been published so far on designing chemical sensors for such compounds. Table 4 summarizes biomimetic/polymer-based approaches for this purpose, including three examples of non-penicillin antibiotics in the last three lines. To the best of our knowledge, there is no literature on Penicillin V yet.
Obviously, the detection limits reported are both in the range of this work and also much below that. These, in part, originates from the transducer technology utilized: For instance, SPR reaches the ppb range by using gold nanoparticles that amplify the SPR signal. Similarly, it is known that electrochemical techniques achieve very good sensitivity at low concentrations. On the other hand, vancomycin and doxycycline are much larger molecules than the penicillins and, thus, lead to larger measuring effects per binding step. The main difference of this work to all of them is the fact that it offers a platform technology that is useful to not only detect PenV, PenG, and Amo salts, but also differentiate them among each other.

Conclusions
This work has demonstrated successful design of QCM sensors comprising MIP thin films' sensor based on acrylates for the detection of PenV-K, PenG-K, and Amo-Na penicillins. AFM imaging of MIP and NIP surfaces clearly reveal observable differences in the textures and roughness of MIP surfaces compared to NIP surfaces, which clearly demonstrates the influence of penicillin on the polymer morphology after template removal, even though it does not show the imprinted cavities per se. Interestingly, the same polymer leads to optimal results for penicillin V and amoxicillin, whereas it turned out necessary to amend the recipe slightly for obtaining the highest possible sensor responses for penicillin G.
QCM results clearly demonstrate that sensors using these MIPs offer both appreciable sensitivity and selectivity. Especially the latter is worth mentioning given the rather minor differences in the respective molecular structures: The beta-lactam core is, of course, the same for all three molecules. However, PenV and PenG differ only by a single oxygen atom. Nonetheless, each MIP preferably incorporates its respective template. Amoxicillin shows lagers selectivity factors towards both other compounds, respectively. Overall, the synthesized sensor materials offer rapid response and concentration-dependent sensing. Together with the fact that the MIPs rely on very similar polymers, this opens up the way for low-cost sensing arrays to detect different components in mixtures.