Next Article in Journal
In Vitro and Clinical Evaluation of Cannabigerol (CBG) Produced via Yeast Biosynthesis: A Cannabinoid with a Broad Range of Anti-Inflammatory and Skin Health-Boosting Properties
Next Article in Special Issue
Electrochemical Analysis of Sulfisoxazole Using Glassy Carbon Electrode (GCE) and MWCNTs/Rare Earth Oxide (CeO2 and Yb2O3) Modified-GCE Sensors
Previous Article in Journal
On the Nature of the Bonding in Coinage Metal Halides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modelling of Glyphosate Molecularly Imprinted Polymer-Based Microsensor with Multiple Phenomena

by
Fares Zouaoui
1,
Saliha Bourouina-Bacha
2,
Mustapha Bourouina
2,
Nadia Zine
1,
Abdelhamid Errachid
1 and
Nicole Jaffrezic-Renault
1,*
1
Institute of Analytical Sciences, University of Lyon, 69100 Villeurbanne, France
2
Faculty of Technology, University of Bejaia, 06000 Bejaia, Algeria
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(2), 493; https://doi.org/10.3390/molecules27020493
Submission received: 10 December 2021 / Revised: 9 January 2022 / Accepted: 10 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Electrochemical Sensors and Cells for Environmental Applications)

Abstract

:
The massive and careless use of glyphosate (GLY) in agricultural production raises many questions regarding environmental pollution and health risks, it is then important to develop simple methods to detect it. Electrochemical impedance spectroscopy (EIS) is an effective analytical tool for characterizing properties at the electrode/electrolyte interface. It is useful as an analytical procedure, but it can also help in the interpretation of the involved fundamental electrochemical and electronic processes. In this study, the impedance data obtained experimentally for a microsensor based on molecularly imprinted chitosan graft on 4-aminophenylacetic acid for the detection of glyphosate was analyzed using an exact mathematical model based on physical theories. The procedure for modeling experimental responses is well explained. The analysis of the observed impedance response leads to estimations of the microscopic parameters linked to the faradic and capacitive current. The interaction of glyphosate molecules with the imprinted sites of the CS-MIPs film is observed in the high frequency range. The relative variation of the charge transfer resistance is proportional to the log of the concentration of glyphosate. The capacitance decreases as the concentration of glyphosate increases, which is explained by the discharging of the charged imprinted sites when the glyphosate molecule interacts with the imprinted sites through electrostatic interactions. The phenomenon of adsorption of the ions in the CMA film is observed in the low frequency range, this phenomenon being balanced by the electrostatic interaction of glyphosate with the imprinted sites in the CS-MIPs film.

1. Introduction

Glyphosate (GLY) is a non-selective foliar systemic herbicide [1]. GLY has become the most widely used phytosanitary molecule, especially as a weedkiller in agriculture. It is further used in various non-agricultural contexts such as forestry and household maintenance [2]. Glyphosate can cause harmful effects on animals [3,4], insects and birds [5], variety of animals in the aquatic system [6], microbial, fungal and bacterial biomasses [7]. The relationship between the use of GLY and adverse effects on human health has been the subject of numerous studies. In March 2015, the International Agency for Research on Cancer (IARC, Lyon, France), a department of the World Health Organization (WHO), classified glyphosate as “probably carcinogenic to humans”, in l ‘inserting into Category 2A, representing substances with limited evidence of carcinogenicity to humans and sufficient evidence to animals [8]. In May 2016, a joint expert meeting of WHO and the Food and Agriculture Organization of the United Nations (FAO) on pesticide residues in the environment and foods, concluded that “glyphosate is unlikely to present a carcinogenic risk to humans through the diet” [9].
Consequently, there is more interest in the development of selective, rapid, simple and inexpensive GLY detection methods. Electrochemical sensors are devices capable of meeting this demand. The main challenge for the electrochemical detection of glyphosate is its low electroactivity below an accessible potential window [10]. Therefore, in order to achieve low concentrations for the detection of GLY, the technique of immobilizing biological elements (antibody, enzyme, microorganism, DNA), on the surface of the working electrode, are generally used for GLY recognition to increase sensitivity and specificity, which are typically lacking in chemical sensors. However, as they are biomolecules, their use is limited to mild conditions of pH, temperature and nature of the medium in order to avoid denaturation of the protein. In addition, the use of these biological elements is not that simple. Their supply and immobilization at the surface of the electrodes are complex [11,12].
Due to the need to overcome these limitations, an alternative has been developed which consists of using molecularly imprinted polymers (MIPs). MIPs were used as a specific recognition material having binding sites of size and shape complementary to target molecules (template). To generate cavity-selective MIPs, a monomer is simply polymerized in the presence of the target molecule and a crosslinking agent resulting in the formation of a highly crosslinked template/monomer complex. MIPs are characterized by very interesting properties such as physical resistance, robustness, resistance to high pressures and temperatures and high inertia towards various chemicals. In addition, MIPs have high affinity and selectivity towards target molecules; they are comparable to natural receptors [13,14,15,16,17].
Sensors based on molecularly imprinted polymers (MIPs) are being developed for the detection of GLY. The functional monomers most often used in synthesis are MCA (N-methacryloyl-L-cysteine) [18] and PPy (Polypyrrole) [19,20].
In these works, the electrochemical detection of glyphosate was obtained by voltammetric methods. A more recent study, using MIP based on polyacrylamide cryogels, glyphosate was detected through capacitive measurements [21].
The choice of the functional monomer in the synthesis of MIPs is an essential step given its capacity to provide complementary interactions with the target molecules [22]. Chitosan (CS) is a natural polysaccharide obtained by deacetylation of chitin which is the most abundant non-toxic, biodegradable and biocompatible natural amino polysaccharide [23]. Chitosan has three types of reactive functional groups, an amine group and two primary and secondary hydroxyl groups. The presence of reactive functional groups on the polysaccharide chain of CS gives it flexibility, excellent film-forming capacity, good adhesion, biocompatibility, high mechanical resistance, possibility of undergoing structural modifications, which makes it suitable for the preparation of MIPs and for the construction of electrochemical sensors [24,25].
CS-MIPs electrochemical sensors have been reported in the state of the art. They are based on the electrodeposition of CS on the surface of the electrodes [26,27]. The electrodeposition of chitosan on a conductive surface by pH gradient causes release of H2. The blogging of these molecules inside the membrane makes it non-homogeneous [28]. The grafting of chitosan on the diazonium salt by the formation of strong interactions probably leads to more homogeneous membranes and to the resolution of interface problems [29].
Electrochemistry describes the chemical phenomena coupled with reciprocal exchanges of electrical energy when an electrolytic chemical medium interacts with an electrical circuit. It allows atomic scale transformations of matter and generations of different reactive species by transfer of electrons between an electrode and a substrate in an electrolytic solution controlled by an electric current or by a voltage [30]. Consequently, the exchange of electrons at the interface between an electronic conductor and the electrolytic medium containing oxidizer/reducer couples (Ox/Red) corresponds to the field of electrochemical reactions, while the reciprocal electron exchange reaction between different Ox/Red couples correspond to the field of oxidation-reduction reactions [31].
Electroactive species have several “modes of transport” to reach the site of the electron transfer reaction, namely the electrode / electrolyte interface such as convection, diffusion and migration. The phenomena related to the movement of matter can be approached at the microscopic level [30]. The electrochemical process can be affected by the phenomenon of adsorption on the surface of the electrode in contact with the solution. In terms of specific adsorption, electrosorption refers to the formation of a chemical bond between the charged species of the electrolyte and the surface of the electrode [32].
Several electrochemical methods are frequently used for physico-chemical analyzes of the electrode/electrolyte interface, such as conductimetry [33], potentiometry [34], amperometry [35], voltammetry [36], and impedimetric. These electrochemical techniques involve three essential quantities which are the current (I), the potential (E) and the time (t).
Electrochemical impedance spectroscopy (EIS) is one of the non-stationary techniques that allow the analysis of complex electrochemical systems. The EIS has established itself as a very powerful analytical tool given its great ability to study the interfacial properties of electrode/electrolyte [37]. The impedance method consists in imposing a sinusoidal disturbance in low amplitude potential, around a chosen voltage E0, between the working electrode and the reference electrode. The sinusoidal current response is thus measured with a continuous component I0, which corresponds to the potential E0 (I0 having a zero value when E0 is the dropout potential. This response is accompanied by a phase shift (φ) depending on the electrical elements present in the system. This system can be considered as a ≪black box≫, which reacts by emitting a signal I when it is subjected to a disturbance E [38].
This work was focused on the modeling of a glyphosate impedimetric microsensor. It was based on the covalent clawing of molecularly imprinted chitosan on 4-aminophenylacetic acid. Mathematical models based on physical theories are developed to analyze the data obtained experimentally. Analysis of the observed impedimetric response leads to the estimation of the microscopic parameters of the sensors. The validation of the different models is carried out by comparing the experimental and theoretical impedance. The response of ‘miniaturized developed’ sensors is closely related to the presence of specific molecular sites on the film of the biopolymer at the surface of the working macro-electrode, and also related to factors inherent in the structure of the sensor itself and in its mode of operation. Explaining the operation in the depths of this sensor requires in-depth knowledge of all the parameters involved in such a design.
Modeling is a simplified representation of a real physical system or phenomenon, making it possible to reproduce its functioning, to analyze, to explain and to predict certain aspects. Modeling is a tool which helps to understand the intrinsic mechanisms of these analytical instruments. It makes it possible to find relationships between the variables and the parameters that are supposed to influence the metrological characteristics of these sensors, for example the resistance of the solution, the resistance to charge transfer in the different ‘layers’ that make up the sensor, or geometry. Modeling is essential for optimization before a device or process goes to market. A “reliable” model is one that will simulate a sensor in real conditions. Therefore, there arises the problem of determining the parameters, which intervene in the equations of the model [39,40].

2. Materials and Methods

2.1. Reagents and Apparatus

Chitosan (CS, molecular weight 250 kDa with a degree of deacetylation 80–95%), ethanol (99.8%), methanol (purity 99.9%), hydrochloric acid (37%, HCl), Sulfuric acid (purity 95%), 1-Ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC), 4- aminophenylacetic acid (CMA), sodium nitrite (NaNO2), sodium hydroxide (NaOH) and glyphosate (GLY) were obtained from Sigma Aldrich, France. N-Hydroxysuccinimide (NHS) was purchased from Acros Organics, France.
Multichannel potentiostat (Biologic-EC-Lab VMP3) analyser was used for all electrochemical measurements. Transducer with four bare-gold working microelectrodes (WE, S = 0.64 mm²/unit), two Ag/AgCl reference microelectrodes (RE, S = 0.13 mm²/unit) and one counter microelectrode (CE, S = 1.37 mm²). It was fabricated at the National Center for Microelectronics (CNM), CSIC, Barcelona, Spain.

2.2. Preparation of the Microsensor

The Au microelectrodes were rinsed with ethanol and deionized water in ultrasound for 10 min then exposed to UV/ozone for 30 min.
CMA solution (3 mM) was prepared in an aqueous solution of NaNO2 (20 mM) and HCl (20 mM). The mixture was stored at 4 °C for 5 min to slow down the chemical reaction between the different reagents. CMA was electrodeposited into gold microelectrodes by cyclic voltammetry for 5 cycles from −1.2 to 0.0 V at scan rate of 50 mV/s. After electrodeposition, the microelectrodes were rinsed with deionized water and dried using nitrogen. Carboxylic functions of CMA were activated by incubating of electrodes for 1 h in 1-ethyl-3- (3 dimethylaminopropyl) carbodiimide (EDC) (0.4 M) and N-Hydroxysuccinimide (NHS) (0.1 M) prepared in absolute ethanol. Then, the devices are rinsed carefully with pure ethanol and dried under light flow.
0.01 g/mL of CS was prepared in 0.1 M acetic acid and ultrasonicated until complete solubilization of CS. Then, 1 mg/mL of GLY was prepared in CS solution. This mixture was mixed for 2 h to obtain GLY/CS complex, then the pH of the solution was adjusted to 5 using 0.1 M NaOH. Drops of GLY/CS were deposited on the microelctrodes modified with CMA for 24 h. Microelectrodes were rinsed with deionized water then they were incubated in 0.5 M H2SO4 solution for 1 h to cross-link CS. In the end, GLY has been removed from CS by incubating microelectrodes in acetic acid/methanol solution (1:1, v/v) for 30 min. Thus, an electrochemical microsensor based on CS molecularly imprinted film modified with CMA was developed for the specific recognition of GLY.
The analytical results and characteristics of the developed microsensor are well presented and detailed in an earlier publication [29].

2.3. EIS detection Method

Microsensors was characterized using EIS (Initial potential E = 0.2 V. Highest Freq = 100 khz, Lowest Freq = 100 mHz) in phosphate buffer saline solution (PBS) with 5 mM ferro-ferricyanide ([Fe(CN)6]3−/4−) before and after incubation in deionized water containing different concentration of GLY (0.31 pg/mL to 50 ng/mL) for 30 min. EIS was used to investigate changes in charge transfer resistance of the film depending on the GLY concentration.

3. Modeling of the Microsensor

3.1. Mathematical Model

In the models discussed below, the constituent phases have different conductivity values, but the conductivities are of the same type, they are all either ionic or electronic. The reaction at the electrode may include a charge transfer step with formation of ions adsorbed on the surface of the membranes. This is usually visualized by the presence of a three-phase boundary where the electronic conductor, ionic conductor, and pores meet [41].
The impedance spectrum characterizing the microsensor designed by grafting the molecularly imprinted chitosan onto the diazonium salt (CMA) and its equivalent circuit are shown in Figure 1.
The overall impedance spectrum forms two well separated semicircles: From 100 KHz to 17 Hz: This first semicircle is associated with a load transfer resistor (Rct1) in parallel with the constant phase elements (CPE1). At high frequency (100 KHz), the intersection of the impedance curve with the abscissa axis enables the electrolyte resistance Rs to be determined. From 17 Hz to 100 mHz: this second half-circle is modeled by a load transfer resistor (Rct2) in parallel with the constant phase elements (CPE2) [42].
The capacitive current is that associated with the charging or discharging of the electrochemical double layer. The capacitive current does not cause a change in the chemical composition, but only in the distribution of the interfacial electric charges. When a metal is placed in contact with an electrolyte at a given potential, an electrical interface is immediately developed. The area of the metal/solution interface where charged species exist is called the electrical double layer, which then behaves as if an electrical capacitor characterized by a capacitor Cdl [43]. Indeed, it has been established that surface heterogeneities of modified electrodes (roughness, porosity, presence of impurities, absorption on the double layer.) can be the cause of a frequency dispersion of the impedance spectra [44]. In these cases, Constant Phase Element (CPE) are used instead of capacitors to take into account the non-ideality of the system. In the literature, different definitions have been proposed. The CPE formula often used for adjustment is defined by the equation [45]:
Z CPE = 1 Q ( j ω ) n
j = 1 = cos ( π 2 ) + j sin ( π 2 )
Q en (sn−1) ou (F s(n−1)) represents the CPE coefficient, −1 ≤ n ≤ 1 is the correction factor. ω = 2πf, where f represents the frequency measured in Hz.
In particular, ZCPE represents an inductance if n = −1, a resistance if n = 0, a Warburg impedance if n = 0.5 and a capacitance if n = 1 [46,47]. Equation (1) can also be written:
Z CPE = 1 Q ω n [ cos ( n π 2 ) j sin ( n π 2 ) ]
The global real and imaginary impedances equivalent to the circuit of Figure 1 are given by the following equations:
Re ( Z ) = R s + R ct 1 ( 1 + R ct 1 Q 1 ω n 1 cos n 1 π 2 ) 1 + ( R ct 1 Q 1 ω n 1 ) 2 + 2 R ct 1 Q 1 ω n 1 cos n 1 π 2 + R ct 2 ( 1 + R ct 2 Q 2 ω n 2 cos n 2 π 2 ) 1 + ( R ct 2 Q 2 ω n 2 ) 2 + 2 R ct 2 Q 2 ω n 2 cos n 2 π 2
Im ( Z ) = R ct 1 2 Q 1 ω n 1 sin n 1 π 2 1 + ( R ct 1 Q 1 ω n 1 ) 2 + 2 R ct 1 Q 1 ω n 1 cos n 1 π 2 + R ct 2 2 Q 2 ω n 2 sin n 2 π 2 1 + ( R ct 2 Q 2 ω n 2 ) 2 + 2 R ct 2 Q 2 ω n 2 cos n 2 π 2
where Re (Z) is the real part of the impedance, −Im (Z) is the imaginary part of the impedance. Q1 represents the coefficient of CPE1, n1 the correction factor of CPE1, Q2 is the coefficient of CPE2, n2 the correction factor of CPE2.
Changes in potential and current in the electrolyte lead to the concept of ohmic drop described as the electrolyte resistance Rs. This resistance is an important factor contributing to the global impedance of the cell [48]. The impedance of the resistance of the solution is:
Z R S =   R s
R s ( ) =   ρ l A
where ρ resistivity of the solution (Ω cm), l the coating thickness and A is the area of the coating [49].
The resistivity of the ferri/ferrocyanide solution is calculated by the following relationships:
ρ = 1 σ
σ = i q i λ i C i
With, σ: conductivity of the solution (s/m) who can be calculate from redox ion couple, Ci the concentration of the ion (mol/m3), qi the number of charges of the ion, λi the equivalent molar ionic conductivity.
The reactions at the electrodes involve the redox couple ferri/ferrocyanide, in which ferricyanide is the oxidant and the ferrocyanide is the reducing agent. For COx = CRed = C, and for a simple one-electron process (n = 1), Rct is given by Equation (10):
R tc = RT F 2 Ak 0 C
With; R: Ideal gas constant (J·mol−1·K−1), T: Temperature (K), F: Faraday constant (C·mo1−1), A: surface of the working electrode (cm²), k0 standard electron transfer rate constant (cm/s).

3.2. Numerical Simulation

To determine the various parameters and validate the proposed model, the response of the sensor was simulated by using Equations (1) and (2) using Matlab software. The parameters were varied with a very small step (Table 1), then the optimal values were determined by minimizing the error between the experimental values and the predicted values. The simulations should make it possible to find the link between the impedance measurement and the various parameters exposed. The developed program consists of two distinct parts. The first part consists in determining the parameters of the circuit Rs + CPE1/Rct1, after validation of this step the program proceeds to determine the parameters of the circuit CPE2/Rct2.
The error between the experimental and theoretical values is calculated according to the relations (11) and (12). The iterations are stopped when e1,e2 ≤ e0 (10−3).
Error ( e 1 ) = | Re ( Z ) theoretical Re ( Z ) experimental Re ( Z ) theoretical |
Error ( e 2 ) = | ( Im ( Z ) ) theoritical ( Im ( Z ) ) experimental ( Im ( Z ) ) theoretical |

4. Results and Discussions

4.1. Model Validation

In order to validate the physical model, the simulation results are compared with the experimental data. We have obtained the results presented in Figure 2. We can first see that the results obtained theoretically are consistent with the experimental values obtained in ref [28]. Thus, the proposed physical model is in agreement with our simulations and it showed good agreement with the results obtained in the laboratory. Overall, the errors are less than 10−3.

4.2. Analysis of Theoretical Results

The values of each parameter of the model are presented in Table 2. The thickness of the membranes (CS + CMA) deposited on the electrode determined by the model is 36 µm. This value was found to vary of less than 1 µm from one sensor to another. This thickness is constant throughout the experiment indicating the stability of the sensor in the analysis medium. Thus, the resistance of the solution Rs is 151,14 Ohm, this resistance is also constant during the experiment.
For the parameters of the first semi-circle, the value of the standard rate constant k1° decreases with the increase in the concentration of the incubated GLY, going from 5,5 µm/s to 2 µm/s according to a linear law whose equation is (k 1° = −0.284 log [GLY] +4.19) with a correlation coefficient R² = 0.997 (Figure 3a). The decrease in the electron transfer rate caused an increase in the charge transfer resistance Rct1, this linear variation is given by the following regression equation (k 1° = −8.285 Rct1 + 5.35; R² = 0.961) (Figure 3b). The increase in Rct1 is proportional to the increase in GLY concentration. It went from 15132 Ohm for a zero GLY concentration to 41614 Ohm for a concentration of 50 ng/mL of incubated GLY (Rct1 = 3343.5 log [GLY] +14507.2; R² = 0.966) (Figure 3c). Moreover, the value of Q1, associated with the behavior of CPE1, is inversely proportional to the concentration of GLY going from 5.6 × 10−7 Sn·Ω−1 ([GLY] = 0) to 4.56 × 10−7 Sn·Ω−1 ([GLY] = 50 ng/mL) (Q1 = −0.16 × 10−7 log [GLY] + 5.83 × 10−7; R² = 0.977) (Figure 3d). The value of n1 which is 0.79 indicates that the surface is inhomogeneous. n1 being close to 1, it means that CPE1, included in the model, is nearly a pure capacitance.
The value of Q1 decreases with increasing GLY concentration. This decrease makes it possible to conclude that the capacitance of the CS-MIPs film and the concentration of GLY are inversely proportional. The same behavior has been observed in the state of the art where the capacitance of a MIPs sensor based on N, N-methylenebisacrylamide cryogel varies inversely proportional with the concentration of glyphosate [21]. The phenomenon was explained by the discharging of the charged imprinted sites when the glyphosate molecule interacts with the imprinted sites through electrostatic interactions.
Regarding the parameters of the second semi-circle, the rate constant k2° and the coefficient Q2 increase proportionally with the increase in the concentration of the incubated GLY. k2° went from 1.92 µm/s to 5.21 µm/s (k2° = 0.55 log [GLY] +0.78; R² = 0.986) (Figure 4a), while Q went from 1.07 × 10−5 Sn·Ω−1 to 2.3 × 10−5 Sn·Ω−1 (Q2 = 0.1714 × 10−5 log [GLY] + 0.966 × 10−5; R² = 0.983) (Figure 4d). While the resistance to load transfer decreases and it varies linearly as a function of the speed constant k2° (k2° = −1.563 Rct2 + 7.334; R² = 0.993) (Figure 4b) and as a function of the logarithm of the GLY concentration (Rct2 = −3547.1 log [GLY] +41865.1; R² = 0.989) (Figure 4c). These variations are recorded with the change from a zero incubated concentration to a concentration of 50 ng/mL of glyphosate. As for the factor n2, its value is estimated at 0.96 indicating that the surface is more relatively homogeneous.
Note that the constituent phases have different values of conductivity. Rct1, k1°, n1 and Q1 describe the processes that occur in the membrane of chitosan. Consequently, all the variations recorded on these parameters are due to the variation in the concentration of glyphosate in the membrane of MIPs, the greater this concentration; the more the film is resistive, thus causing a decrease in the faradic and capacitive current.
While Rct2, k2°, n2 and Q2 correspond to the film of CMA which separates the membrane MIPs from the surface of the sensor, as it can also represent the polarization of the electrode. The impedance of the polarizable electrode can be represented as a parallel connection of the charge transfer resistor and the double layer capacitor. It was assumed in the literature that such a distribution could result from a microscopic roughness always present on the solid surfaces, as it can be caused by a dispersion of the capacitance of interfacial origin, linked to the adsorption of ions and to in homogeneities surface chemicals. The concentration of the redox couple changes in the pores with an exponential gradient. This problem has been addressed in the literature [50,51]. The accumulation of ions on the electrode surface causes an increase in faradic and capacitive current.
The relative variation of the resistance to charge transfer of the electrode is presented using the following equation |Rct1 − Rct1blank|/Rct1blank (ΔR/R). The calculations of ΔR/R were carried out with Rct1, which is related to the variation in the concentration of GLY. This parameter is found to be linearly proportional to the logarithmic value of GLY concentrations in the range of 0.31 pg/mL to 50 ng/mL, as shown in Figure 5, with a correlation coefficient of 0.999. The linear regression equation is: ΔR/R = 0.221 log[GLY] − 0.02.
The detection limit (LOD) of the considered sensor is estimated at 1 fg/mL. LOD was obtained from the equation 3 S/m [52], where S is the residual standard deviation of the linear regression and m is the slope of the regression line. The results showed a low detection limit and high sensitivity.
The parameters of the CS-MIPs/CMA/Au sensor are compared to another impedance-meter sensor for the detection of GLY which is based on electropolymerized polypyrrole films, as a conductive layer layer between the gold microelectrode surface and the molecularly imprinted chitosan film (CS-MIPs/PPy/Au) [27]. The impedance spectrum of the PPy-based sensor is formed by a single semicircle. This result allows concluding that in the present work, the formation of two half-circles in the impedance spectrum is due to the presence of the CMA. The electron transfer rate in the MIPs films of the CS-MIPs/PPy/Au sensor is much higher (55 µm/s) than that of the CS-MIPs/CMA/Au (5.5 µm/s) sensor due to the presence of the PPy conductive layer which facilitates the electron transfer. The interest of the covalent attachment of the CS-MIPs film on an underlayer is to eliminate the inhomogeneity of the CS-MIPs film caused by the trapping of the H2 molecules released during the electrodeposition of the CS film. Due to the insulating property of the CMA film, accumulation of charged species occurred (second semi-circle in the low frequency range), which is not the case using polypyrrole film [27].
Despite this enormous difference, the sensitivities of both sensors (with PPy and with CMA), calculated with the relative resistance variation (R/R), are close (0.30 for CS-MIPs/PPy/Au and 0.221 for CS-MIPs/CMA/Au), the detection limits are the same: 1 fg/mL and their possible recyclability is the same.

4.3. Comparison between MIPs and NIP

The response of the NIP sensor was modeled with the model developed. Figure 6 shows the theoretical and experimental evolution of Re (Z) as a function of −Im (Z). A good fit is achieved between the experimental data and the theoretical data with a fairly low error (<10−3). The parameters of the NIP sensor are listed in Table 3. These parameters vary insignificantly depending on the concentration of GLY, in particular Rct1. From this point, it comes that the non-specific adsorption of glyphosate is very low.
The coating thickness of the NIP sensor (56.8 µm) is greater than that of the MIPs (36 µm). Therefore, the resistance of the resulting solution is greater on the NIP sensor (238.5 Ohm) compared to MIP (151.14 Ohm).
The value n1 of the MIPs film (n1 = 0.79) is lower than that of the NIP (n1 = 0.85) which makes it possible to say that the MIPs films are more porous due to the presence of free sites, specific to GLY. This point can also explain the compared values of the electron transfer rate constants: the electron transfer rate constant in the MIPs film is greater (k1° = 5.5 µm/s) compared to the NIP film (k1° = 1.4 µm/s). The values of n2 characterizing the CMA films are close indicating that there is not a big difference in the morphology of the two films (n2 (MIPs) = 0.96 − n2 (NIP) = 0.98).

4.4. Effects of Differents Parameters on the Impedance Response of CS-MIPs/CMA/Au

4.4.1. Effect of Coefficient n

The coefficient n1 has an effect on the impedance of the first semicircle which is found at high frequencies. The more n1 decreases the more the maximum of –Im (Z) of the first semicircle and its value at the end towards the middle frequencies decreases, in contrast with the minimum value of Re (Z) which increases and its maximum remains constant. The coefficient n1 has no effect on the second semi-circle, except for the initial value of the imaginary impedance −Im (Z) (Figure 7a).
With the reduction of the coefficient n2 we notice that the maximum of Re.
Experimentally, n1 and n2 can be modified by varying the thickness of the CMA and MIP films, respectively. These parameters can be controlled by varying the number of cycles during the electrodeposition. In addition, n2 can be decreased by the use of a porogene solvent in the synthesis of the CS-MIPs [53], which can cause increase the porosity and tortuosity on the CS film.
(Z) and of –Im (Z) of the first semicircle increases, and that a better definition of the semicircle is obtained for a maximum value n2. For the second semicircle, the maxima of Re (Z) and –Im (Z) in the low frequency domain decrease with the decrease in n2. The larger n2, the more well defined this semicircle (Figure 7b).

4.4.2. Effect of the Standard Rate Constant k°

With the decrease in the value of k1° from 10−2 to 10−3 cm/s we notice that the real and imaginary impedance of the first semi-circle has increased, for the second semi-circle the maximum of –Im (Z) remained constant and the real impedance increased. The passage at a speed k1° of 10−4 cm/s recorded a significant increase in the imaginary impedance of the first semicircle, while the maximum of the real impedance remained the same as for k1° equal to 10−3 cm/s. For the second semicircle, while it is less defined, the maximum value of –Im (Z) is larger and that of Re (Z) is smaller than the values recorded for k1° equal to 10−3 cm/s. With the decrease in k1°, we notice the appearance of a single semicircle with greater resistance (Figure 8), which is observed for the NIP film.
In Figure 9, we have shown the effect of varying the rate constant k2° on the impedance spectrum. For large values of k2° we notice the appearance of a spectrum with a single semi-circle. With the decrease in k2°, two semicircles appear on the spectrum. The maximum of Re (Z) of the first semicircle remains constant, and the maximum of –Im (Z) decreases with the decrease of k2°. Whereas, for the second semicircle, the maximum of –Im (Z) decreases and that of Re (Z) increases with the decrease k2°. For sufficiently low values of k2°, we notice that the semi-circle at low frequencies approaches the shape of a semi-infinite line.
For obtaining a higher value of k1° and k2°, the thickness of the two films should be decreased allowing to reduce the charge transfer resistance. On another side, increasing of the concentration of the redox couple should be increased the values of k1° and k2°. Or, for example, the integration of conductive nanomaterials in the CS and CMA films could also increase the value of k°.

4.4.3. Effect of Coefficient Q

For larger values of Q1, the model gives a spectrum with a single semicircle (Figure 10a). With the decrease in the value of Q1, we obtain spectra with two semi-circles. The latter is more defined for a minimum Q1 value. On the two semicircles, the maximum of −Im (Z) varies little, and the maximum of Re (Z) remains constant.
With the variation of the coefficient of CPE (Q2), we obtain a single semicircle for low values of Q2. With the increase of Q2 the model gives two increasingly better-defined semicircles. The maximum of –Im (Z) varies little, and the maximum of Re (Z) remains constant. For sufficiently large values in Q2 we notice the disappearance of part of the second semi-circle which is located towards the low frequency areas (Figure 10b).
So, the sensor is better if the thickness of the electronic double layer is low. Q1 and Q2 can be decreased by the variation of the movement of the ions by agitation of the electrolyte, by increasing of the temperature, by the variation of the applied potential, or also the changes in the morphology of the CS and CMA films [54,55,56].

5. Conclusions

The desire to understand the intrinsic operating mechanisms of the developed microsensor, or to open up new avenues in the implementation of new design processes, requires in-depth knowledge of all the parameters of influence. In this study, the impedimetric response of the CS-MIPs/CMA/Au microsensors has been modeled using mathematical model based on physical laws. The phenomena at the electrode/electrolyte interface were described. The various parameters were determined by minimizing the error between the experimental and theoretical values. Then, the effect of each parameter and GLY concentration on the response signal was shown. The first half-circle in the high frequency range is due to the interaction of glyphosate molecules with the imprinted sites of the CS-MIPs film. The relative variation of the charge transfer resistance is proportional to the log of the concentration of glyphosate. The capacitance decreases as the concentration of glyphosate increases, which is explained by the discharging of the charged imprinted sites when the glyphosate molecule interacts with the imprinted sites through electrostatic interactions. The second half-circle of impedance in the low frequency range is due to the phenomenon of adsorption of the ions in the CMA film, this phenomenon being balanced by the electrostatic interaction of glyphosate with the imprinted sites in the CS-MIPs film. This phenomenon is due to the insulating properties of the CMA film but is not observed when CS-MIPs film is attached to a conductive film such as polypyrrole. Despite these differences of electrochemical behavior, CS-MIPs/CMA/Au and CS-MIPs/PPy/Au microsensors present the same detection limit: 1 fg/mL.
The developed model here for a MIP-chitosan based sensor is applicable to other MIP based sensors. It appears here that electrostatic interactions between template and imprinted sites induce changes in capacitance. The effect of the electrical properties of an underlayer can change the morphology the Nyquist plot. The morphology of the Nyquist plot of the NIP based sensor is totally different from that of the corresponding MIP based sensor.
In order to increase the sensitivity of detection, some improvements could be brought in the preparation of the MIP-based sensor: the thickness of MIP film and of its underlayer should be obtained thinner, the conductivity of the underlayer should be improved by the addition of conductive nanomaterials, the porosity of the MIP film should be improved by using a porogene, the concentration of the redox probe could be increased.

Author Contributions

Investigation, F.Z.; resources, A.E.; writing—original draft preparation, N.Z.; writing—review and editing, N.J.-R.; supervision, S.B.-B. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the financial support of the EU H2020 WIDESPREAD Program entitled Bionanosens grant agreement # 951887.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Valavanidis, A. Glyphosate, the most widely used herbicide. In Health and Safety Issues. Why Scientists Differ in Their Evaluation of Its Adverse Health Effects; University Campus Zografou: Athens, Greece, 2018. [Google Scholar]
  2. Torretta, V.; Katsoyiannis, I.; Viotti, P.; Rada, E. Critical Review of the Effects of Glyphosate Exposure to the Environment and Humans through the Food Supply Chain. Sustainability 2018, 10, 950. [Google Scholar] [CrossRef] [Green Version]
  3. Krüger, M.; Schledorn, P.; Schrödl, W.; Hoppe, H.-W.; Lutz, W.; Shehata, A.A. Detection of Glyphosate Residues in Animals and Humans. J. Environ. Anal. Toxicol. 2014, 4, 1–5. [Google Scholar]
  4. Séralini, G.-E.; Clair, E.; Mesnage, R.; Gress, S.; Defarge, N.; Malatesta, M.; Hennequin, D.; de Vendômois, J.S. Republished Study: Long-Term Toxicity of a Roundup Herbicide and a Roundup-Tolerantgenetically Modified Maize. Environ. Sci. Eur. 2014, 26, 1–17. [Google Scholar] [CrossRef] [Green Version]
  5. Herbert, L.T.; Vázquez, D.E.; Arenas, A.; Farina, W.M. Effects of Field-Realistic Doses of Glyphosate on Honeybee Appetitive Behaviour. J. Exp. Biol. 2014, 217, 3457–3464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Van Bruggen, A.H.C.; He, M.M.; Shin, K.; Mai, V.; Jeong, K.C.; Finckh, M.R.; Morris, J.G. Environmental and Health Effects of the Herbicide Glyphosate. Sci. Total Environ. 2018, 616–617, 255–268. [Google Scholar] [CrossRef] [PubMed]
  7. Ohtonen, R.; Munson, A.; Brand, D. Soil Microbial Community Response to Silvicultural Intervention in Coniferous Plantation Ecosystems. Ecol. Appl. 1992, 2, 363–375. [Google Scholar] [CrossRef]
  8. IARC Monograph on Glyphosate. Available online: Https://Www.Iarc.Fr/En/Media-Center/Iarcnews/2016/Glyphosate_IARC2016.Ph (accessed on 9 December 2021).
  9. Xu, J.; Smith, S.; Smith, G.; Wang, W.; Li, Y. Glyphosate Contamination in Grains and Foods: An Overview. Food Control 2019, 106, 106710. [Google Scholar] [CrossRef]
  10. Cahuantzi-Muñoz, S.L.; González-Fuentes, M.A.; Ortiz-Frade, L.A.; Torres, E.; Ţălu, Ş.; Trejo, G.; Méndez-Albores, A. Electrochemical Biosensor for Sensitive Quantification of Glyphosate in Maize Kernels. Electroanalysis 2019, 31, 927–935. [Google Scholar] [CrossRef] [Green Version]
  11. Yáñez-Sedeño, P.; Campuzano, S.; Pingarrón, J.M. Electrochemical Sensors Based on Magnetic Molecularly Imprinted Polymers: A Review. Anal. Chim. Acta 2017, 960, 1–17. [Google Scholar] [CrossRef]
  12. Scheller, F.W.; Zhang, X.; Yarman, A.; Wollenberger, U.; Gyurcsányi, R.E. Molecularly Imprinted Polymer-Based Electrochemical Sensors for Biopolymers. Curr. Opin. Electrochem. 2019, 14, 53–59. [Google Scholar] [CrossRef] [Green Version]
  13. BelBruno, J.J. Molecularly Imprinted Polymers. Chem. Rev. 2018, 119, 94–119. [Google Scholar] [CrossRef]
  14. Janczura, M.; Luliński, P.; Sobiech, M. Imprinting Technology for Effective Sorbent Fabrication: Current State-of-Art and Future Prospects. Materials 2021, 14, 1850. [Google Scholar] [CrossRef]
  15. Leibl, N.; Haupt, K.; Gonzato, C.; Duma, L. Molecularly Imprinted Polymers for Chemical Sensing: A Tutorial Review. Chemosensors 2021, 9, 123. [Google Scholar] [CrossRef]
  16. Spychalska, K.; Zając, D.; Baluta, S.; Halicka, K.; Cabaj, J. Functional Polymers Structures for (Bio) Sensing Application—A Review. Polymers 2020, 12, 1154. [Google Scholar] [CrossRef]
  17. Saylan, Y.; Akgönüllü, S.; Yavuz, H.; Ünal, S.; Denizli, A. Molecularly Imprinted Polymer Based Sensors for Medical Applications. Sensors 2019, 19, 1279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Prasad, B.B.; Jauhari, D.; Tiwari, M.P. Doubly Imprinted Polymer Nanofilm-Modified Electrochemical Sensor for Ultra-Trace Simultaneous Analysis of Glyphosate and Glufosinate. Biosens. Bioelectron. 2014, 59, 81–88. [Google Scholar] [CrossRef] [PubMed]
  19. Mazouz, Z.; Rahali, S.; Fourati, N.; Zerrouki, C.; Aloui, N.; Seydou, M.; Yaakoubi, N.; Chehimi, M.; Othmane, A.; Kalfat, R. Highly Selective Polypyrrole MIP-Based Gravimetric and Electrochemical Sensors for Picomolar Detection of Glyphosate. Sensors 2017, 17, 2586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Xu, J.; Zhang, Y.; Wu, K.; Zhang, L.; Ge, S.; Yu, J. A Molecularly Imprinted Polypyrrole for Ultrasensitive Voltammetric Determination of Glyphosate. Microchim. Acta 2017, 184, 1959–1967. [Google Scholar] [CrossRef]
  21. Han, X.-J.; Ji, X.-F.; Zhang, Q.; Sun, J.-W.; Sun, P.-X.; Pan, W.-J.; Wang, J.; Yang, C. Giant “Molecular Capacitor” Arrays—Portable Sensors to Determine Ionizable Compounds. J. Electroanal. Chem. 2020, 865, 114108. [Google Scholar] [CrossRef]
  22. Turiel, E.; Martín-Esteban, A. Molecularly Imprinted Polymers for Sample Preparation: A Review. Anal. Chim. Acta 2010, 668, 87–99. [Google Scholar] [CrossRef]
  23. Dutta, P.K.; Dutta, J.; Tripathi, V. Chitin and Chitosan: Chemistry, Properties and Applications. J. Sci. Ind. Res. 2004, 63, 20–31. [Google Scholar]
  24. Wang, Y.; Wang, E.; Wu, Z.; Li, H.; Zhu, Z.; Zhu, X.; Dong, Y. Synthesis of Chitosan Molecularly Imprinted Polymers for Solid-Phase Extraction of Methandrostenolone. Carbohydr. Polym. 2014, 101, 517–523. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, X.; Wu, M.; Tang, W.; Zhu, Y.; Wang, L.; Wang, Q.; He, P.; Fang, Y. Simultaneous Electrochemical Determination of Ascorbic Acid, Dopamine and Uric Acid Using a Palladium Nanoparticle/Graphene/Chitosan Modified Electrode. J. Electroanal. Chem. 2013, 695, 10–16. [Google Scholar] [CrossRef]
  26. Zouaoui, F.; Bourouina-Bacha, S.; Bourouina, M.; Abroa-Nemeir, I.; Ben Halima, H.; Gallardo-Gonzalez, J.; El Alami El Hassani, N.; Alcacer, A.; Bausells, J.; Jaffrezic-Renault, N.; et al. Electrochemical Impedance Spectroscopy Determination of Glyphosate Using a Molecularly Imprinted Chitosan. Sens. Actuators B Chem. 2020, 309, 127753. [Google Scholar] [CrossRef]
  27. Zouaoui, F.; Bourouina-Bacha, S.; Bourouina, M.; Alcacer, A.; Bausells, J.; Jaffrezic-Renault, N.; Zine, N.; Errachid, A. Experimental Study and Mathematical Modeling of a Glyphosate Impedimetric Microsensor Based on Molecularly Imprinted Chitosan Film. Chemosensors 2020, 8, 104. [Google Scholar] [CrossRef]
  28. Zouaoui, F.; Bourouina-Bacha, S.; Bourouina, M.; Jaffrezic-Renault, N.; Zine, N.; Errachid, A. Electrochemical Sensors Based on Molecularly Imprinted Chitosan: A Review. TrAC Trends Anal. Chem. 2020, 130, 115982. [Google Scholar] [CrossRef]
  29. Zouaoui, F.; Bourouina-Bacha, S.; Bourouina, M.; Alcacer, A.; Bausells, J.; Jaffrezic-Renault, N.; Zine, N.; Errachid, A. Electrochemical Impedance Spectroscopy Microsensor Based on Molecularly Imprinted Chitosan Film Grafted on a 4-Aminophenylacetic Acid (CMA) Modified Gold Electrode, for the Sensitive Detection of Glyphosate. Front. Chem. 2021, 9, 621057. [Google Scholar] [CrossRef] [PubMed]
  30. Lefrou, C.; Fabry, P.; Poignet, J.-C. Électrochimie: Concepts Fondamentaux Illustrés, Nouvelle Edition; Collection Grenoble Sciences; EDP Sciences: Les Ulis, France, 2013. [Google Scholar]
  31. Miomandre, F.; Sadki, S.; Audebert, P.; Méallet-Renault, R.; Amatore, C. Électrochimie: Des Concepts Aux Applications, 4th ed.; Sciences sup. Chimie; Dunod: Paris, France, 2019. [Google Scholar]
  32. Sakar, H.; Celik, I.; Balcik Canbolat, C.; Keskinler, B.; Karagunduz, A. Electro-Sorption of Ammonium by a Modified Membrane Capacitive Deionization Unit. Sep. Sci. Technol. 2017, 52, 2591–2599. [Google Scholar] [CrossRef]
  33. Jaffrezic-Renault, N.; Dzyadevych, S.V. Conductometric Microbiosensors for Environmental Monitoring. Sensors 2008, 8, 2569–2588. [Google Scholar] [CrossRef] [Green Version]
  34. Moreno, T.V.; Malacarne, L.C.; Baesso, M.L.; Qu, W.; Dy, E.; Xie, Z.; Fahlman, J.; Shen, J.; Astrath, N.G. Potentiometric Sensors with Chalcogenide Glasses as Sensitive Membranes: A Short Review. J. Non-Cryst. Solids 2018, 495, 8–18. [Google Scholar] [CrossRef]
  35. Kanyong, P.; Krampa, F.D.; Aniweh, Y.; Awandare, G.A. Enzyme-Based Amperometric Galactose Biosensors: A Review. Microchim. Acta 2017, 184, 3663–3671. [Google Scholar] [CrossRef] [Green Version]
  36. Farghaly, O.A.; Hameed, R.A.; Abu-Nawwas, A.-A.H. Analytical Application Using Modern Electrochemical Techniques. Int. J. Electrochem. Sci. 2014, 9, 3287–3318. [Google Scholar]
  37. Park, S.-M.; Yoo, J.-S. Peer Reviewed: Electrochemical Impedance Spectroscopy for Better Electrochemical Measurements. Anal. Chem. 2003, 75, 455A–461A. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Sacco, A. Electrochemical Impedance Spectroscopy: Fundamentals and Application in Dye-Sensitized Solar Cells. Renew. Sustain. Energy Rev. 2017, 79, 814–829. [Google Scholar] [CrossRef]
  39. Gao, M.; Hazelbaker, M.S.; Kong, R.; Orazem, M.E. Mathematical Model for the Electrochemical Impedance Response of a Continuous Glucose Monitor. Electrochim. Acta 2018, 275, 119–132. [Google Scholar] [CrossRef]
  40. Tkach, V.; de Oliveira, S.C.; de Oliveira, S.C.; Neves, V.S.; Ojani, R.; Yagodynets, P.I. The Electrochemical Oxidation of Procarbazine on Cobalt (III) Oxyhydroxide and Its Mathematical Evaluation. Anal. Bioanal. Electrochem 2016, 8, 432. [Google Scholar]
  41. Los, P.; Laviron, E. The Ac Faradaic Impedance of a Redox System with Langmuirian Adsorption of the Reactants. J. Electroanal. Chem. 1997, 432, 85–99. [Google Scholar] [CrossRef]
  42. Lee, M.; Zine, N.; Baraket, A.; Zabala, M.; Campabadal, F.; Caruso, R.; Trivella, M.G.; Jaffrezic-Renault, N.; Errachid, A. A Novel Biosensor Based on Hafnium Oxide: Application for Early Stage Detection of Human Interleukin-10. Sens. Actuators B Chem. 2012, 175, 201–207. [Google Scholar] [CrossRef]
  43. Sapountzi, E. Développement de Nouveaux Supports Basés Sur Des Nanofibres de Matériaux Hybrides Électrofilées Pour Le Développement de Biocapteurs Électrochimiques; HAL Publisher: Lyon, France, 2015. [Google Scholar]
  44. Brug, G.; Van Den Eeden, A.; Sluyters-Rehbach, M.; Sluyters, J. The Analysis of Electrode Impedances Complicated by the Presence of a Constant Phase Element. J. Electroanal. Chem. 1984, 176, 275–295. [Google Scholar] [CrossRef]
  45. Huang, J.; Li, Z.; Liaw, B.Y.; Zhang, J. Graphical Analysis of Electrochemical Impedance Spectroscopy Data in Bode and Nyquist Representations. J. Power Sources 2016, 309, 82–98. [Google Scholar] [CrossRef]
  46. Stoynov, Z.; Vladikova, D. Differential Impedance Analysis; Marin Drinov Academic Publishing House: Sofia, Bulgaria, 2005. [Google Scholar]
  47. Zoltowski, P. On the Electrical Capacitance of Interfaces Exhibiting Constant Phase Element Behaviour. J. Electroanal. Chem. 1998, 443, 149–154. [Google Scholar] [CrossRef]
  48. Jorcin, J.-B. Spectroscopie d’impédance Electrochimique Locale: Caracterisation de La Delamination Des Peintures et de La Corrosion Des Alliages Al-Cu. Ph.D. Thesis, University of Toulouse, Toulouse, France, 2007. [Google Scholar]
  49. Hinderliter, B.R.; Croll, S.G.; Tallman, D.E.; Su, Q.; Bierwagen, G.P. Interpretation of EIS Data from Accelerated Exposure of Coated Metals Based on Modeling of Coating Physical Properties. Electrochim. Acta 2006, 51, 4505–4515. [Google Scholar] [CrossRef]
  50. Terbouche, A.; Lameche, S.; Ait-Ramdane-Terbouche, C.; Guerniche, D.; Lerari, D.; Bachari, K.; Hauchard, D. A New Electrochemical Sensor Based on Carbon Paste Electrode/Ru(III) Complex for Determination of Nitrite: Electrochemical Impedance and Cyclic Voltammetry Measurements. Measurement 2016, 92, 524–533. [Google Scholar] [CrossRef]
  51. Tanujit, B.; Asokan, S. Electrochemical Impedance Spectroscopy Study of AgI-Ag2O-MoO3 Glasses: Understanding the Diffusion, Relaxation, Fragility and Power Law Behavior. arXiv 2019, arXiv:1906.07438. [Google Scholar] [CrossRef]
  52. Shrivastava, A.; Gupta, V. Methods for the Determination of Limit of Detection and Limit of Quantitation of the Analytical Methods. Chron. Young Sci. 2011, 2, 21. [Google Scholar] [CrossRef]
  53. Sellergren, B.; Shea, K.J. Influence of Polymer Morphology on the Ability of Imprinted Network Polymers to Resolve Enantiomers. J. Chromatogr. A 1993, 635, 31–49. [Google Scholar] [CrossRef]
  54. Fellner, P.; Matiašovský, K. On the Temperature Dependence of the Double-Layer Capacitance. Electrochim. Acta 1972, 17, 233–236. [Google Scholar] [CrossRef]
  55. Pangarov, N.; Christova, I.; Atanasov, M.; Kertov, V. Double-Layer Capacitance Measurements by a Pulse Potentiostatic Method. Electrochim. Acta 1967, 12, 717–721. [Google Scholar] [CrossRef]
  56. De Levie, R. On Porous Electrodes in Electrolyte Solutions: I. Capacitance Effects. Electrochim. Acta 1963, 8, 751–780. [Google Scholar] [CrossRef]
Figure 1. (a) EIS of the microsensor (E = 0.2 V. Freq =100 kHz−100 mHz) in PBS with 5 mM [Fe(CN)6]3−/4−. (b) Equivalent electrical circuit (Rs + CPE1/Rct1+ CPE2/Rct2).
Figure 1. (a) EIS of the microsensor (E = 0.2 V. Freq =100 kHz−100 mHz) in PBS with 5 mM [Fe(CN)6]3−/4−. (b) Equivalent electrical circuit (Rs + CPE1/Rct1+ CPE2/Rct2).
Molecules 27 00493 g001
Figure 2. Experimental and theoretical EIS of CS-MIPs/CMA/Au in 5 mM [Fe(CN)6]3−/4− initial potential E = 0.2 V; f = 100 kHz–100 mHz, for different concentrations of GLY (0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL and 50 ng/mL).
Figure 2. Experimental and theoretical EIS of CS-MIPs/CMA/Au in 5 mM [Fe(CN)6]3−/4− initial potential E = 0.2 V; f = 100 kHz–100 mHz, for different concentrations of GLY (0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL and 50 ng/mL).
Molecules 27 00493 g002
Figure 3. Calibration curves of CS-MIPs/CMA/Au. (a) k1° = f(log[GLY]), (b) k1° = f(Rct1) for different [GLY], (c) Rct1 = f(log[GLY]), (d) Q1= f(log[GLY]). ([GLY] = 0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL–50 ng/mL).
Figure 3. Calibration curves of CS-MIPs/CMA/Au. (a) k1° = f(log[GLY]), (b) k1° = f(Rct1) for different [GLY], (c) Rct1 = f(log[GLY]), (d) Q1= f(log[GLY]). ([GLY] = 0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL–50 ng/mL).
Molecules 27 00493 g003
Figure 4. Calibration curves of CS-MIPs/CMA/Au. (a) k2° = f(log[GLY]), (b) k2° = f(Rct2) for different [GLY], (c) Rct2 = f(log[GLY]), (d) Q2 = f(log[GLY]). ([GLY] = 0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL–50 ng/mL).
Figure 4. Calibration curves of CS-MIPs/CMA/Au. (a) k2° = f(log[GLY]), (b) k2° = f(Rct2) for different [GLY], (c) Rct2 = f(log[GLY]), (d) Q2 = f(log[GLY]). ([GLY] = 0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL–50 ng/mL).
Molecules 27 00493 g004
Figure 5. Calibration curves of CS-MIPs/CMA/Au (ΔR/R = 0.221 log[GLY] − 0.02; R² = 0.975).
Figure 5. Calibration curves of CS-MIPs/CMA/Au (ΔR/R = 0.221 log[GLY] − 0.02; R² = 0.975).
Molecules 27 00493 g005
Figure 6. Experimental and theoretical EIS of CS-MIPs/CMA/Au in 5 mM [Fe(CN)6]3−/4− initial potential E = 0.2 V; f = 100 kHz–100 mHz, for different concentrations of GLY (0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL and 50 ng/mL).
Figure 6. Experimental and theoretical EIS of CS-MIPs/CMA/Au in 5 mM [Fe(CN)6]3−/4− initial potential E = 0.2 V; f = 100 kHz–100 mHz, for different concentrations of GLY (0.31 pg/mL–6.25 pg/mL–125 pg/mL–2.5 ng/mL and 50 ng/mL).
Molecules 27 00493 g006
Figure 7. (a) Effect of coefficient n1 on the variation of the impedance response of CS-MIPs/CMA/Au for n2 = 0.95. (b) Effect of coefficient n2 on the variation of the impedance response of CS-MIPs/CMA/Au for n1 =0.8. In 5 mM [Fe(CN)6]3−/4−. f = 100 kHz–100 mHz. A = 0.0064 cm², T = 298 K, l = 36 µm, k1° = 5 µm/s, Q1 = 5 × 10−7 sn−1, k2° = 2 µm/s, Q2 = 10−5 sn−1.
Figure 7. (a) Effect of coefficient n1 on the variation of the impedance response of CS-MIPs/CMA/Au for n2 = 0.95. (b) Effect of coefficient n2 on the variation of the impedance response of CS-MIPs/CMA/Au for n1 =0.8. In 5 mM [Fe(CN)6]3−/4−. f = 100 kHz–100 mHz. A = 0.0064 cm², T = 298 K, l = 36 µm, k1° = 5 µm/s, Q1 = 5 × 10−7 sn−1, k2° = 2 µm/s, Q2 = 10−5 sn−1.
Molecules 27 00493 g007
Figure 8. Effect of the rate constant k1° on the impedance response of CS-MIPs/CMA/Au. For T = 298 K, A = 0.0064 cm², l = 60µm, Q1 = 5 × 10−7 sn−1, n1 = 0.8, k2° = 2µm/s, Q2 = 10−5 sn−1, n2 = 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Figure 8. Effect of the rate constant k1° on the impedance response of CS-MIPs/CMA/Au. For T = 298 K, A = 0.0064 cm², l = 60µm, Q1 = 5 × 10−7 sn−1, n1 = 0.8, k2° = 2µm/s, Q2 = 10−5 sn−1, n2 = 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Molecules 27 00493 g008
Figure 9. Effect of the rate constant k2° on the impedance response of CS-MIPs/CMA/Au. For T = 298 K, A = 0.0064 cm², l = 60µm, k1° = 5µm/s, Q1 = 5 × 10−7, n1 = 0.8 sn−1, Q2 = 10−5 sn−1, n2 = 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Figure 9. Effect of the rate constant k2° on the impedance response of CS-MIPs/CMA/Au. For T = 298 K, A = 0.0064 cm², l = 60µm, k1° = 5µm/s, Q1 = 5 × 10−7, n1 = 0.8 sn−1, Q2 = 10−5 sn−1, n2 = 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Molecules 27 00493 g009
Figure 10. (a) Effect of the rate constant Q1 on the impedance response of CS-MIPs/CMA/Au for Q2 = 10−5 sn−1. (b) Effect of the rate constant Q2 on the impedance response of CS-MIPs/CMA/Au for Q1 = 5 × 10−7. T = 298 K, A = 0.0064 cm², l = 60µm, k1° = 5µm/s, Q1 = 5 × 10−7 sn−1, Q2 = 10−5 sn−1, n1 = 0.8, n2= 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Figure 10. (a) Effect of the rate constant Q1 on the impedance response of CS-MIPs/CMA/Au for Q2 = 10−5 sn−1. (b) Effect of the rate constant Q2 on the impedance response of CS-MIPs/CMA/Au for Q1 = 5 × 10−7. T = 298 K, A = 0.0064 cm², l = 60µm, k1° = 5µm/s, Q1 = 5 × 10−7 sn−1, Q2 = 10−5 sn−1, n1 = 0.8, n2= 0.95. [Fe(CN)6]3−/4− = 5 mM. f = 100 KHz–100 mHz.
Molecules 27 00493 g010
Table 1. Variation interval of different parameters.
Table 1. Variation interval of different parameters.
ParametersVariation RangeUnit
l[10−3,10−5]cm
k1°[10−3,10−5]cm·s−1
Q1[10−5,10−8]sn·Ω−1
n1[−1,1]/
k2°[10−3,10−5]cm·s−1
Q2[10−5,10−8]sn·Ω−1
n2[−1,1]/
Table 2. Different parameters of CS-MIPs/CMA/Au estimated by the model for different concentrations of GLY.
Table 2. Different parameters of CS-MIPs/CMA/Au estimated by the model for different concentrations of GLY.
[GLY] l
(μm)
Rs
(Ω)
k1° (μm/s)Rct1
(Ω)
Q1·107
(Sn−1)
n1k2°
(μm/s)
Rct2
(Ω)
Q2·105
(Sn−1)
n2e1 × 104e2 × 104
Blank36151.145.515,1325.600.791.9243,3481.070.964.203.59
0.31 pg/mL36151.143.523,7795.450.792.2532,6381.400.967.043.88
6.25 pg/mL36151.143.126,8485.170.792.8728,9991.650.966.526.84
125 pg/mL36151.142.730,8255.020.793.5823,2481.770.967.509.94
2.5 ng/mL36151.142.434,6784.860.794.1719,9592.080.960.980.38
50 ng/mL36151.142.041,6144.560.795.2114,0832.300.962.889.57
Table 3. Different parameters of CS-NIP/CMA/Au estimated by the model for different concentrations of GLY.
Table 3. Different parameters of CS-NIP/CMA/Au estimated by the model for different concentrations of GLY.
[GLY] l
(μm)
Rs
(Ω)
k1°
(μm/s)
Rct1
(Ω)
Q1·107
(Sn−1)
n1k2°
(μm/s)
Rct2
(Ω)
Q2·107
(Sn−1)
n2e1 × 104e2 × 104
Blank56.8238.51.3959,5634.920.8517.0248853.560.982.224.96
0.31 pg/mL56.8238.51.4158,8504.650.8518.6044693.810.984.153.58
6.25 pg/mL56.8238.51.3661,2714.710.8513.0963494.020.983.545.51
125 pg/mL56.8238.51.4557,2344.440.858.9892574.210.985.638.36
2.5 ng/mL56.8238.51.3959,5714.860.859.2889533.850.985.454.48
50 ng/mL56.8238.51.4358,2234.640.859.5087423.660.984.484.22
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zouaoui, F.; Bourouina-Bacha, S.; Bourouina, M.; Zine, N.; Errachid, A.; Jaffrezic-Renault, N. Mathematical Modelling of Glyphosate Molecularly Imprinted Polymer-Based Microsensor with Multiple Phenomena. Molecules 2022, 27, 493. https://doi.org/10.3390/molecules27020493

AMA Style

Zouaoui F, Bourouina-Bacha S, Bourouina M, Zine N, Errachid A, Jaffrezic-Renault N. Mathematical Modelling of Glyphosate Molecularly Imprinted Polymer-Based Microsensor with Multiple Phenomena. Molecules. 2022; 27(2):493. https://doi.org/10.3390/molecules27020493

Chicago/Turabian Style

Zouaoui, Fares, Saliha Bourouina-Bacha, Mustapha Bourouina, Nadia Zine, Abdelhamid Errachid, and Nicole Jaffrezic-Renault. 2022. "Mathematical Modelling of Glyphosate Molecularly Imprinted Polymer-Based Microsensor with Multiple Phenomena" Molecules 27, no. 2: 493. https://doi.org/10.3390/molecules27020493

APA Style

Zouaoui, F., Bourouina-Bacha, S., Bourouina, M., Zine, N., Errachid, A., & Jaffrezic-Renault, N. (2022). Mathematical Modelling of Glyphosate Molecularly Imprinted Polymer-Based Microsensor with Multiple Phenomena. Molecules, 27(2), 493. https://doi.org/10.3390/molecules27020493

Article Metrics

Back to TopTop