# High-Performance ZIF-7@PANI Electrochemical Sensor for Simultaneous Trace Cadmium and Lead Detection in Water Samples: A Box-Behnken Design Optimization Approach

Sondes Guesmi, Nashwan Ali, Nadhem Missaoui, Zouhaier Aloui, Chama Mabrouk, Carlos Castilla-Martinez, Fraj Echouchene, Houcine Barhoumi, Nicole Jaffrezic-Renault, Hamza Kahri

► To cite this version:

Sondes Guesmi, Nashwan Ali, Nadhem Missaoui, Zouhaier Aloui, Chama Mabrouk, et al. High-Performance ZIF-7@PANI Electrochemical Sensor for Simultaneous Trace Cadmium and Lead Detection in Water Samples: A Box-Behnken Design Optimization Approach. Sensors, 2025, 25 (5), pp.1336. 10.3390/s25051336. [hal-05225431]

HAL Id: hal-05225431

https://hal.science/hal-05225431v1

Submitted on 27 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# High-Performance ZIF-7@PANI Electrochemical Sensor for Simultaneous Trace Cadmium and Lead Detection in Water Samples: A Box-Behnken Design Optimization Approach

Sondes Guesmi <fnref n="1" />, Nashwan H. Ali <fnref n="2" />,* ¹D, Nadhem Missaoui ¹D, Zouhaier Aloui <fnref n="3" />, Chama Mabrouk ¹D, Carlos Castilla Martinez <fnref n="4" />, Fraj Echouchene <fnref n="5" />,<fnref n="6" /> ¹D, Houcine Barhoumi ¹D, Nicole Jaffrezic-Renault <fnref n="7" /> ¹D and Hamza Kahri <fnref n="1" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

Academic Editor: Michele Penza

Received: 19 January 2025

Revised: 4 February 2025

Accepted: 10 February 2025

Published: 22 February 2025

**Citation:** Guesmi, S.; Ali, N.H.; Missaoui, N.; Aloui, Z.; Mabrouk, C.; Martinez, C.C.; Echouchene, F.; Barhoumi, H.; Jaffrezic-Renault, N.; Kahri, H. High-Performance ZIF-7@PANI Electrochemical Sensor for Simultaneous Trace Cadmium and Lead Detection in Water Samples: A Box-Behnken Design Optimization Approach. *Sensors* **2025**, *25*, 1336.
https://doi.org/10.3390/s25051336

**Copyright:** © 2025 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license
(https://creativecommons.org/
licenses/by/4.0/).

## Highlights

### What are the main findings?

* Development of an innovative sensor based on ZIF-7 and polyaniline for the simultaneous detection of cadmium (Cd<fnref n="2" />⁺) and lead (Pb<fnref n="2" />⁺) ions in water.

* Optimization of the sensor's performance using the Box-Behnken design (BBD), leading to enhanced sensitivity and selectivity.

### What is the implication of the main finding?

* The ZIF-7@PANI/GCE sensor demonstrated low detection limits (2.96 nM for Pb<fnref n="2" />⁺ and 10.6 nM for Cd<fnref n="2" />⁺), making it highly suitable for trace-level heavy metal detection.

* The successful validation of the sensor using real water samples highlights its potential for practical applications in environmental monitoring and water quality assessment.

**Abstract:** This study presents the development of an innovative electrochemical sensor for the simultaneous detection of cadmium (Cd<fnref n="2" />⁺) and lead (Pb<fnref n="2" />⁺) ions in environmental samples. The sensor is developed based on a composite material of zeolite imidazolate framework ZIF-7 and polyaniline (PANI), referred to as ZIF-7@PANI, where ZIF-7 is rapidly synthesized at room temperature and polyaniline used to enhance the conductivity of the composite. Characterization via X-ray diffraction (DRX), scanning electron microscopy (SEM), and Fourier transform infrared spectroscopy (FTIR) confirmed successful synthesis. The composite was applied to a glassy carbon electrode (GCE) using drop-casting for heavy metal ion detection. Experimental parameters—including pH, incubation time, deposited quantity, and drying time—were optimized using the Box-Behnken design. Under optimal conditions, the ZIF-7@PANI/GCE sensor demonstrated a broad dynamic concentration range (0.002–1 µM for Pb<fnref n="2" />⁺ and 0.02–30 µM for Cd<fnref n="2" />⁺), with low detection limits (2.96 nM for Pb<fnref n="2" />⁺ and 10.6 nM for Cd<fnref n="2" />⁺). It also exhibited strong anti-interference properties and high recovery rates (85–110%), highlighting its potential for real practical applications.

<footnote n="1">¹ Laboratory of Interfaces and Advanced Materials, University of Monastir, Monastir 5000, Tunisia; missaoui.nadhem1@gmail.com (N.M.); houcine.barhoumi@fsm.rnu.tn (H.B.)</footnote>

<footnote n="2">² Applied Chemistry Department, Applied Sciences Collage, University of Samarra, Samarra 34010, Iraq</footnote>

<footnote n="3">³ Chemistry Department, College of Science, King Khalid University (KKU), P.O. Box 9004, Abha 61413, Saudi Arabia</footnote>

<footnote n="4">⁴ Institut Européen des Membranes, IEM—UMR 5635, University of Montpellier, CNRS, 34090 Montpellier, France</footnote>

<footnote n="5">5 University of Sousse, Higher Institute of Applied Sciences and Technology of Sousse, Cité Ettafala, Ibn Khaldoun, Sousse 4003, Tunisia</footnote>

<footnote n="6">6 Laboratory of Electronics and Microelectronics LR99ES30, Faculty of Sciences, University of Monastir, Monastir 5000, Tunisia</footnote>

<footnote n="7">⁷ Institute of Analytical Sciences, University of Lyon, 69622 Villeurbanne, France</footnote>

<footnote>* Correspondence: dr.nashwan@uosamarra.edu.iq (N.H.A.); hamza.kahri@univ-tours.fr (H.K.)</footnote>
**Keywords:** ZIF-7@PANI composite; electrochemical sensor; heavy metal detection; Box-Behnken design; cadmium and lead sensing

# 1. Introduction

Water scarcity, exacerbated by climate change—particularly global warming—poses a significant challenge for emerging countries [1]. Among the many concerns, heavy metal contamination in water reservoirs stands out [2], primarily caused by industrial activity and agriculture. This issue impacts billions of people worldwide who lack access to clean and safe drinking water [3]. Lead (Pb²⁺) and cadmium (Cd²⁺) ions, in particular, pose severe health risks, contributing to lead poisoning, digestive issues, increased cancer susceptibility (notably breast cancer), and, in acute cases, abdominal muscle cramps and toxic hepatitis [4,5]. Traditional analytical methods for detecting heavy metal ions in water, such as atomic emission spectroscopy (AES) [6], atomic absorption spectroscopy (AAS) [7], X-ray fluorescence spectrometry [8], and inductively coupled plasma mass spectrometry (ICP-MS) [9], are often expensive and time-consuming. Therefore, developing sensitive sensors capable of simultaneous, real-time detection of Cd²⁺ and Pb²⁺ ions is essential. Recent advancements in technology have significantly improved detection capabilities, particularly in enhancing both sensitivity and selectivity for environmental monitoring applications [10]. Electrochemical detection offers an ideal solution due to its cost-effectiveness, user-friendly operation, high sensitivity, and compatibility with portable testing equipment [11,12]. Electrode materials play a key role in electrochemical sensor performance [13], leading to extensive research aimed at improving their efficiency. Metal-organic frameworks (MOFs) have attracted considerable attention due to their large specific surface area, abundant adsorption sites, and tunable functionality [14]. They have been widely applied in catalysis [15], gas storage [16], separation and purification processes [17], drug delivery [18], and environmental remediation [19]. More recently, MOFs have shown great promise in developing highly sensitive electrochemical sensors for heavy metal detection, contributing to water quality monitoring [20,21]. Nonetheless, MOF-based electrodes often suffer from limitations such as low conductivity, challenges in large-scale preparation, and poor repeatability. To address these issues, researchers have explored incorporating conductive materials into MOFs, including metals [22], carbon nanomaterials [23,24], biopolymers [25], and polymers [26].

Zeolitic Imidazolate Frameworks (ZIFs), a subclass of MOFs, are known for their tetrahedral structure, formed through strong bonds between metal ions and four imida-zolate ligands [27]. Various methods have been reported for producing ZIFs, including solvothermal [28], sol-gel [29], surfactant-assisted methods [30], and microwave/ultrasound-assisted techniques [31]. However, these methods often demand long reaction times due to slow interactions, highlighting the need for high-performance, eco-friendly, and easily accessible materials with controlled crystal properties. Templates, like polymers or surfac-tants, offer a way to achieve ZIFs with controlled morphology, preventing non-uniform agglomeration. Growing interest in room-temperature synthesis has led to more efficient ap-proaches, such as the polyethylene glycol (PEG)-templated method introduced by Nadhem et al. for producing highly microporous nanoZIF-8 [32,33].

In this study, a PEG-based soft template was used to rapidly and easily synthesize porous nano-ZIF-7. This material was then combined with polyaniline (PANI) to form a ZIF-7@PANI composite, which was applied as a sensitive layer for the simultaneous elec-trochemical detection of Cd²⁺ and Pb²⁺ ions. Key experimental parameters were optimized using the Box-Behnken design to improve the sensor's electroanalytical performance. This
study focused on the investigation of sensitivity, anti-interference properties, and their applicability to real samples for the detection of Cd²⁺ and Pb²⁺ ions in aqueous solutions.

## 2. Experimental

### 2.1. *Chemicals*

Zn(NO₃)₂·6H₂O, benzimidazole (blm; C₇H₆N₂; 98%), polyethylene glycol (PEG, 20,000 g/mol), dimethylformamide (DMF, 99%), methanol (99%), triethylamine (TEA; (C₂H₅)₃N; 99%), and aniline (99%) were purchased from Sigma-Aldrich (Darmstadt, Germany). (NH₄)₂S₂O₈ (APS), Pb(NO₃)₂·6H₂O, and Cd(NO₃)₂·6H₂O were obtained from Fluka Chemika (Buchs, Switzerland). The acetate buffer was used to dissolve Pb(NO₃)₂ and Cd(NO₃)₂·3H₂O, resulting in the preparation of the stock solution of Pb²⁺ and Cd²⁺ ions.

### 2.2. Preparation of ZIF-7@PANI Nanocomposite

#### 2.2.1. Synthesis of Nano-ZIF-7

Nano-ZIF-7 was synthesized following a methodology similar to a previously reported publication for nanocrystalline ZIF-7 [32]. Zn(NO₃)₂·6H₂O and blm served as precursors, while PEG (20,000 g/mol) acted as a soft template. In the initial step, 400 mg of PEG powder was dispersed in 20 mL of DMF, followed by the addition of Zn(NO₃)₂·6H₂O (2.5 g, 8.40 mmol) with mild stirring for 2 min to form solution A. Separately, benzimidazole (3.08 g, 26.07 mmol) was dispersed in DMF (20 mL), and TEA, employed as a base, was gradually added (7.26 mL) under stirring over 2 min to prepare solution B. Upon complete dissolution, solution A was slowly added to solution B and stirred for 5 min at room temperature. The resulting white precipitated was separated by filtration, rinsed with ethanol to remove residual reagents (including PEG), and collected via centrifugation (4000 rpm, 15 min). Finally, the obtained white material was dried in a vacuum oven at 150 °C for 24 h, yielding the final ZIF-7 product.

#### 2.2.2. Synthesis of ZIF-7@PANI Nanocomposite

The ZIF-7@PANI nanocomposite was synthesized according to a prior report [26]. Initially, 0.3 g of ZIF-7 was dispersed in 15 mL of 1 M HCl via sonication. Subsequently, 30 µL of aniline was added to the mixture, which was then sonicated for 20 min. Following this, 0.114 g of APS in 5 mL of 1 mol L⁻¹ HCl was added dropwise to the mixture at 0 °C, and then the mixture was stirred overnight to ensure thorough polymerization. The resulting solid was washed three times alternately with distilled water and ethanol via centrifugation then dried at 60 °C for 24 h. A similar method was used to synthesize pristine PANI by polymerizing aniline (30 µL) in 15 mL of 1 mol L⁻¹ HCl, followed by the addition of APS, without incorporating ZIF-7.

### 2.3. Preparation of Modified Electrodes

Before measurements, the surface of the glassy carbon electrode (GCE) was carefully polished with 0.3 µm alumina powder then immersed in an ethanol-water solution and sonicated for 2 min to ensure thorough cleaning. To modify the electrode surface, a homogeneous suspension of the nanocomposite was prepared by dispersing 2 mg of ZIF-7@PANI in 1 mL of distilled water. From this suspension, a 10 µL was drop-costed onto the electrode surface and left to dry for 3 h. The ZIF-7/GCE and PANI/GCE electrodes were prepared using the same procedure by drop-casting the ZIF-7 and PANI suspensions, respectively, onto the electrode surface for comparative electrochemical analysis.
## 2.4. Characterization Methods

The morphological and structural properties of PANI, ZIF-7, and ZIF-7@PANI were investigated using various techniques. X-ray diffraction (XRD) analysis was carried out using a Bruker D8 Discover diffractometer (Billerica, MA, USA), while Fourier transform infrared spectroscopy (FT-IR) spectra were obtained employing a Perkin Elmer Spectrum 100 spectrometer (Waltham, MA, USA) featuring a universal ATR sampling accessory. Scanning electron microscopy (SEM) imaging was performed with a Supra 40VP Field Emission microscope (Oberkochen, Germany).

## 2.5. Electrochemical Tests

Electrochemical measurements were conducted using a three-electrode configuration within an electrochemical cell. A glassy carbon electrode (diameter: 3 mm) was used as the working electrode, while an Ag/AgCl/KCl electrode and a platinum wire functioned as the reference and counter electrodes, respectively. Electrochemical techniques, including electrochemical impedance spectroscopy (EIS), cyclic voltammetry (CV), and differential pulse voltammetry (DPV), were performed using an Autolab PGSTAT 320N Potentiostat (Metrohm, Milan, Italy), controlled via NOVA 1.5 software for data analysis. Metal ion accumulation was carried out under open circuit conditions by immersing the modified electrode in a stirred buffer solution containing varying metal ion concentrations for an optimized accumulation period. To investigate the electrochemical behavior of the prepared electrodes, cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) were conducted. These analyses were performed in a 5 mM ferrocyanide ($Fe(CN)_6^{3-}/4-)$ solution as probe in 0.1 M KCl.

## 2.6. Water Sample Analysis

Seawater samples were obtained from the SKANES region in the city of Monastir, Tunisia. The seawater supernatant, along with tap and mineral water, was diluted twenty-fold for electrochemical detection without prior treatment. The concentrations of heavy metal ions were assessed by using a standard calibration curve, and the accuracy was verified by adding standard solution to the actual water sample solution.

# 3. Results and Discussions

## 3.1. Structural and Morphological Characterization

The chemical composition of ZIF-7 was confirmed by FTIR analysis (Figure 1). Spectral features between 600 and 1500 cm⁻¹ correspond to the bending and stretching modes of the benzimidazole ring structure [34]. Specifically, the peaks observed at 750 and 1455 cm⁻¹ were attributed to C-H and C=C bonds, respectively [35]. The stretching of =C-H in aromatics appeared at 3032 and 3068 cm⁻¹, while the C-C stretching of the aromatic ring was observed at 1465 and 1611 cm⁻¹ [36]. Additionally, the band at 422 cm⁻¹ was associated with the Zn-N stretching vibration [37]. These findings confirm the successful coordination between benzimidazole and Zn²⁺ in the ZIF-7 nanoparticles.

The FTIR spectrum of polyaniline (PANI) shows absorption peaks at 1560, 1472, 1297, and 1117 cm⁻¹. The peaks observed at 1560 and 1472 cm⁻¹ are associated with the C=C stretching deformations of the quinoid and benzene rings, indicating the presence of emeraldine salt state in PANI. The band at 1297 cm⁻¹ is assigned to the C-N stretching vibration, while the band at 1117 cm⁻¹ is attributed to the quinonoid unit in doped PANI [38]. In the ZIF-7@PANI spectrum, the peak at 1117 cm⁻¹ experiences a blue shift to 1135 cm⁻¹, attributed to the strong electrostatic interaction between PANI and ZIF-7. Additionally, the C=C elongation of the quinoid ring and benzene rings in ZIF-7@PANI shifts from
1560 and 1472 cm⁻¹ in pure PANI to 1570 and 1493 cm⁻¹, suggesting enhanced π-π interactions between the ZIF-7 framework and PANI [39].

<figure><img src="image_3.png" /><figcaption>**Figure 1.** FTIR spectra of PANI, synthesized ZIF-7 and ZIF-7@PANI composite.</figcaption></figure>

The three samples were analyzed by the powder XRD technique. The strong similarity between the simulated and synthesized XRD patterns suggests the production of pure ZIF-7, free from impurities or additional phases (Figure 2). Indexing of the diffractogram reveals significant alignment with the corresponding lattice parameters and space group (R-3). Notably, the observed broadening of the diffraction peaks indicates the formation of smaller ZIF-7 crystals. Figure 2 displays the XRD patterns for the ZIF-7@PANI composite, with peaks at 7.18°, 7.62°, 10.53°, 13.28°, 15.44°, and 19.69° (in 2θ), confirming the presence of ZIF-7 crystals with a cubic structure (R-3) [34,36,40]. Additionally, the peak at 20° confirms the inclusion of PANI in the composite, validating the presence of semi-crystalline PANI planes [41].

Figure 3 shows the SEM images of the different modified screen-printed electrodes (SPCE) constructed using ZIF-7, PANI, and the composite ZIF-7@PANI.

The SEM image of the PANI-modified electrode reveals a non-uniform topography and a spongy morphology (Figure 3a,b). The SEM image of the pure ZIF-7 matrix deposited on the SPCE surface shown in Figure 3c,d displays well-defined octahedral and cubic nanoparticles, ranging in size from 20 nm to 60 nm. After the incorporation of pure ZIF-7 into the polymer, noticeable deformations in the morphology of the composite surface are observed (Figure 3e,f). This change can be attributed to the interlinking between PANI and the pure ZIF-7 matrix. The SEM images in Figure 3e,f strongly indicate the successful formation of the ZIF-7@PANI composites, with small octahedral and cubic crystals intertwined. The resulting composite sample also displays a significant degree of agglomeration.
Following the interpretation of the results from the characterization techniques, the synthesis process of the ZIF-7@PANI composite is illustrated in Scheme 1.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** The powder X-ray diffraction patterns of the ZIF-7 and ZIF-7@PANI.</figcaption></figure>

## 3.2. Electrochemical Characterization of the Modified Electrodes

As shown in Figure 4a, the cyclic voltammogram (CV) of the bare GCE exhibited clear and reversible redox peaks. After functionalizing the working electrode surface with the ZIF-7 matrix, a slight reduction in the electron transport rate was observed at the modified electrode interface, leading to a decrease in peak current compared to the bare electrode. This reduction can be attributed to the relatively low conductivity of ZIF-7, as discussed by Butler et al. (2017) [42]. In contrast, the PANI-modified electrode displayed increased current, due to the conductive nature and fast charge transfer capabilities of polyaniline. Notably, the PANI@ZIF-7-modified electrode displayed the highest current for redox peaks, indicating that combining ZIF-7 and PANI into a mixed matrix enhances electron transfer at the electrode surface. This improvement is due to the high conductivity of PANI and the considerable specific surface area of ZIF-7, making it a suitable choice as an electrochemical sensing platform.

In addition to CV, EIS was employed to obtain additional electrochemical information on the modified electrodes. Under a 0.2 V applied potential, EIS parameters were measured across a frequency range from $1 \times 10^{-3}$ to $1 \times 10^5$ Hz with a 5 mV amplitude. The EIS spectra of different modified GCEs are shown in Figure 4b, featuring both semicircular and linear sections, with the fitting circuit model provided in the inset. The unmodified GCE exhibited a relatively small $R_{ct}$ value of approximately 190 Ω. For the ZIF-7/GCE, the $R_{ct}$ increased to 295 Ω due to its mesoporous structure, which hinders the charge transfer process due to the low conductivity of the ZIF-7 matrix. Conversely, the $R_{ct}$ at the PANI/GCE interface was lower than that of the bare electrode, attributed to the enhanced electrical conductivity of polyaniline. On the ZIF-7@PANI/GCE, the $R_{ct}$ decreased to 140 Ω, suggesting that the ZIF-7 and PANI mixed matrix efficiently coated the surface, creating a porous and highly conductive structure with a significant active surface area.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** SEM images of PANI/SPCE (a,b), ZIF-7/SPCE (c,d), and ZIF-7@PANI/SPCE (e,f). Images (a,c,e) are taken at 25,000 × magnification with a scale bar of 400 nm, while images (b,d,f) are taken at 100,000 × magnification with a scale bar of 100 nm.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Scheme 1.** Schematic illustration of ZIF-7@PANI composite synthesis.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** (a) CV and (b) EIS curves of bare GCE, ZIF-7/GCE, PANI/GCE, and ZIF-7@PANI/GCE in 0.1 M KCl solution containing 5 mM Fe(CN)₆⁻³/⁻⁴.</figcaption></figure>

The feasibility of the ZIF-7@PANI/GCE sensing electrode was evaluated for the simultaneous electrochemical determination of Cd²⁺ and Pb²⁺ ions in the same solution. DPV was employed to investigate the electrochemical response of 10 µM of both Cd²⁺ and Pb²⁺ using various modified electrodes. As depicted in Figure 5, the electrochemical behavior of these ions on different electrodes reveals two distinct peak voltages. For the bare GCE, two relatively weak peaks corresponding to the oxidation of each metal are observed at potentials of -0.92 V and -0.62 V, respectively. In contrast, the ZIF-7/GCE electrode shows a slight decrease in the response currents, and the oxidation potentials shift to -0.86 V and -0.61 V, respectively. When the electrode surface is modified with PANI, the current values increase to 10.27 µA for Cd²⁺ oxidation and 2.06 µA for Pb²⁺ oxidation. Remarkably, the oxidation responses of Cd²⁺ and Pb²⁺ ions increase by 2.5 times on the ZIF-7@PANI/GCE compared to the bare GCE. These results strongly suggest that the combination of PANI and ZIF-7 creates novel pathways for the oxidation of these ions on the electrode surface. The synergy between PANI and ZIF-7 enhances the electrocatalytic activity for Cd²⁺ and Pb²⁺, resulting in higher current responses at lower oxidation potentials for these ions.
<figure><img src="image_8.png" /><figcaption>Figure 5. DPV curves of a 10 µM mixture of Cd²⁺ and Pb²⁺ ions in 0.1 m acetate buffer (pH 5) using different sensors.</figcaption></figure>

### 3.3. Optimization of Experimental Parameters

Electrochemical sensors for Pb²⁺ and Cd²⁺ ions are highly sensitive to operating conditions such as pH, deposited drop volume, incubation time, and drying time. Therefore, optimizing these detection parameters is crucial to achieve a highly efficient and reliable electrochemical sensor. Response Surface Methodology (RSM) is a statistical analysis technique that utilizes experimental data to optimize processes and improve the performance of detection systems [43]. RSM is based on a fractional factorial design, meaning that only certain possible combinations of factor levels are tested, significantly reducing the number of required experiments. The Box-Behnken design is a specific type of experimental design used in RSM [44]. In this design (BBD), each factor is generally studied at three different levels (low, centered, and high) along with three central test points to estimate quadratic effects in the model. In this study, the optimization of four control factors was achieved using the BBD approach to maximize the current response and examine interaction between independent variables. The selected factors included the pH of the medium (X₁, 3–7), deposited drop volume (X₂, 3–15 µL), incubation time (X₃, 1–3 min), and drying time (X₄, 1–3 h). The relationship between the response and independent variables is described by the multiple quadratic regression model (Equation (1)) [45]:

$$Y = a_0 + \sum_{i=1}^{n} a_i X_i + \sum_{i=1}^{n} a_{ii} X_i^2 + \sum_{i=1}^{n} \sum_{j=i+1}^{n} a_{ij} X_i X_j + \epsilon \quad (1)$$

where Y is the predicted response, a₀ is the constant, (Xᵢ and Xⱼ) are the independent variables, aᵢ is the linear effect, aᵢᵢ is the interaction coefficient, aᵢᵢᵢ is the quadratic coefficient, n is the number of optimized factors, and ε is the random error for uncertainties between predicted and determined data. The BBD design required 27 experimental runs, incorporating three central points to estimate quadratic effects. Table 1 outlines the experimental parameters along with their respective levels.
<table><caption>**Table 1.** Levels of the experimental parameters for the Box-Behnken design.</caption><thead><tr><th rowspan="2">Parameters</th><th rowspan="2">Symbol</th><th colspan="3">Parameter Level</th></tr><tr><th>Low (-1)</th><th>Medium (0)</th><th>High (+1)</th></tr></thead><tbody><tr><td>pH</td><td>X<sub>1</sub></td><td>3</td><td>5</td><td>7</td></tr><tr><td>Drop volume (µL)</td><td>X<sub>2</sub></td><td>3</td><td>9</td><td>15</td></tr><tr><td>Incubation time (min)</td><td>X<sub>3</sub></td><td>1</td><td>2</td><td>3</td></tr><tr><td>Drying time (h)</td><td>X<sub>4</sub></td><td>1</td><td>2</td><td>3</td></tr></tbody></table>

The uncoded unit regression equations for the two responses (Cd²⁺ (Equation (2)) and Pb²⁺ (Equation (3))) are given as follows:

$$ \begin{aligned} \hat{I}_{Cd^{2+}} = & -9.84 + 3.81X_1 + 0.24X_2 + 1.73X_3 + 0.56X_4 - 0.44X_1 * X_1 - 0.017X_2 * X_2 \\ & - 0.21X_3 * X_3 - 0.15X_4 * X_4 + 0.0018X_1 * X_2 + 0.078X_1 * X_3 \\ & + 0.122X_1 * X_4 - 0.035X_2 * X_3 + 0.067X_2 * X_4 - 0.60X_3 * X_4 \end{aligned} \quad (2) $$

$$ \begin{aligned} \hat{I}_{Pb^{2+}} = & -24 + 8.68X_1 + 0.91X_2 + 1.73X_3 + 4.31X_4 - 1.02X_1 * X_1 - 0.021X_2 * X_2 \\ & + 0.44X_3 * X_3 - 0.54X_4 * X_4 - 0.0076X_1 * X_2 + 0.064X_1 * X_3 \\ & + 0.36X_1 * X_4 - 0.15X_2 * X_3 - 0.086X_2 * X_4 - 1.32X_3 * X_4 \end{aligned} \quad (3) $$

The statistical analysis results in Table 2 and Figure 6 demonstrate that both models ($\hat{I}_{Cd^{2+}}$ and $\hat{I}_{Pb^{2+}}$) have high R² values and, thus, an excellent overall fit to the data. The values of R²<sub>Adj</sub> suggest that both prediction models are reasonably robust. The performance of two models on new unseen data is reasonably good.

<table><caption>**Table 2.** Statistical analysis of both model ($\hat{I}_{Cd^{2+}}$ and $\hat{I}_{Pb^{2+}}$).</caption><thead><tr><th></th><th>R<sup>2</sup></th><th>R<sup>2</sup><sub>Adj</sub></th><th>R<sup>2</sup><sub>Pred</sub></th></tr></thead><tbody><tr><td>Î<sub>Cd<sup>2+</sup></sub></td><td>78%</td><td>52.3%</td><td>47.3%</td></tr><tr><td>Î<sub>Pb<sup>2+</sup></sub></td><td>88.93%</td><td>76.01%</td><td>36.33%</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Statistical analysis fit of the RSM model (a) for Cd²⁺ current response and (b) for Pb²⁺ current response.</figcaption></figure>
ANOVA analysis was also conducted for each model (Table 3), revealing high F-values and low *p*-values (<0.05) for both, indicating strong statistical significance. The ANOVA table summarizes linear terms, squared terms, and interactions terms. The *p*-values of the linear terms (X₁..₄) show that pH (X₁) is the only significant factor. Regarding the square terms, only X₁ * X₁ is statistically significant for both models while the square term X₂ * X₂ is significant for the 1st model. For interaction terms, the X₃ * X₄ (*p* < 0.05) interaction was the only statistically significant term, indicating an interaction between incubation time (X₃) and drying time (X₄).

<table><caption>**Table 3.** Analysis of variance of both models.</caption><thead><tr><th rowspan="2">Source</th><th rowspan="2">DF</th><th colspan="4">Model ̂I<sub>Cd<sup>2+</sup></sub></th><th colspan="4">Model ̂I<sub>Pb<sup>2+</sup></sub></th></tr><tr><th>SS</th><th>MS</th><th>F-Value</th><th>p-Value</th><th>SS</th><th>MS</th><th>F-Value</th><th>p-Value</th></tr></thead><tbody><tr><td>Model</td><td>14</td><td>22.3247</td><td>1.5946</td><td>8.73</td><td>0.000 **</td><td>155.690</td><td>11.12</td><td>6.89</td><td>0.001 **</td></tr><tr><td>X<sub>1</sub></td><td>1</td><td>0.9391</td><td>0.9391</td><td>5.14</td><td>0.043 *</td><td>28.170</td><td>28.17</td><td>17.44</td><td>0.001 **</td></tr><tr><td>X<sub>2</sub></td><td>1</td><td>0.0095</td><td>0.0095</td><td>0.05</td><td>0.824</td><td>0.143</td><td>0.142</td><td>0.09</td><td>0.771</td></tr><tr><td>X<sub>3</sub></td><td>1</td><td>0.6675</td><td>0.6675</td><td>3.66</td><td>0.080</td><td>0.388</td><td>0.388</td><td>0.24</td><td>0.633</td></tr><tr><td>X<sub>4</sub></td><td>1</td><td>0.0037</td><td>0.0037</td><td>0.02</td><td>0.889</td><td>3.900</td><td>3.899</td><td>2.41</td><td>0.146</td></tr><tr><td>X<sub>1</sub> * X<sub>1</sub></td><td>1</td><td>16.0796</td><td>16.2956</td><td>89.26</td><td>0.000 **</td><td>102.381</td><td>89.29</td><td>55.29</td><td>0.000 **</td></tr><tr><td>X<sub>2</sub> * X<sub>2</sub></td><td>1</td><td>1.7649</td><td>2.0175</td><td>11.05</td><td>0.006 **</td><td>3.493</td><td>3.113</td><td>1.93</td><td>0.190</td></tr><tr><td>X<sub>3</sub> * X<sub>3</sub></td><td>1</td><td>0.1557</td><td>0.2368</td><td>1.30</td><td>0.277</td><td>2.272</td><td>1.015</td><td>0.63</td><td>0.443</td></tr><tr><td>X<sub>4</sub> * X<sub>4</sub></td><td>1</td><td>0.1183</td><td>0.1183</td><td>0.65</td><td>0.437</td><td>1.541</td><td>1.540</td><td>0.95</td><td>0.348</td></tr><tr><td>X<sub>1</sub> * X<sub>2</sub></td><td>1</td><td>0.0019</td><td>0.0019</td><td>0.01</td><td>0.921</td><td>0.034</td><td>0.033</td><td>0.02</td><td>0.888</td></tr><tr><td>X<sub>1</sub> * X<sub>3</sub></td><td>1</td><td>0.0969</td><td>0.0969</td><td>0.53</td><td>0.480</td><td>0.066</td><td>0.066</td><td>0.04</td><td>0.843</td></tr><tr><td>X<sub>1</sub> * X<sub>4</sub></td><td>1</td><td>0.2396</td><td>0.2396</td><td>1.31</td><td>0.274</td><td>2.107</td><td>2.106</td><td>1.30</td><td>0.276</td></tr><tr><td>X<sub>2</sub> * X<sub>3</sub></td><td>1</td><td>0.1793</td><td>0.1793</td><td>0.98</td><td>0.341</td><td>3.204</td><td>3.203</td><td>1.98</td><td>0.184</td></tr><tr><td>X<sub>2</sub> * X<sub>4</sub></td><td>1</td><td>0.6373</td><td>0.6373</td><td>3.49</td><td>0.086</td><td>1.063</td><td>1.062</td><td>0.66</td><td>0.433</td></tr><tr><td>X<sub>3</sub> * X<sub>4</sub></td><td>1</td><td>1.4315</td><td>1.4315</td><td>7.84</td><td>0.016 *</td><td>6.930</td><td>6.929</td><td>4.29</td><td>0.061</td></tr><tr><td>Lack-of-Fit</td><td>10</td><td>2.1053</td><td>0.2105</td><td>4.93</td><td>0.180</td><td>19.336</td><td>1.933</td><td>84.85</td><td>0.012 *</td></tr><tr><td>Pure Error</td><td>2</td><td>0.0854</td><td>0.0427</td><td></td><td></td><td>0.046</td><td>0.022</td><td></td><td></td></tr><tr><td>Total</td><td>26</td><td>24.5154</td><td></td><td></td><td></td><td>175.071</td><td></td><td></td><td></td></tr></tbody></table>

DF: degree of freedom, SS: sum of squares, MS: mean of squares, * significant at 5%, ** significant at 1%.

The desirability function approach was used to optimize the heavy metal sensor by determining the ideal peak current values for Cd²⁺ and Pb²⁺ ions. The response variable (peak current for both ions) was set to “maximum” to attain the best sensor performance or desirability value, while the desired operational parameters (pH, deposited drop volume, incubation time, and drying time) were chosen within the investigated experimental range. Based on the response objective, this methodology was applied to maximize the overall desirability function. Regression analysis identified the optimal conditions as pH = 4.8, DV drop volume = 11.73 µL, incubation time IT = 1 min, and drying time DT = 3 h (Figure 7). The optimal pH of 4.8 for the simultaneous electrochemical detection of Pb²⁺ and Cd²⁺ is supported by electrochemical principles. At this slightly acidic pH, both metal ions remain in their electroactive forms, preventing precipitation as hydroxides, which occurs at higher pH levels [46]. Lower pH values can lead to competition from H⁺ ions, reducing detection sensitivity by interfering with the adsorption of metal ions onto the electrode surface [47]. Additionally, pH 4.8 minimizes interferences from other species in complex matrices, ensuring selective detection of Pb²⁺ and Cd²⁺ [48]. Research indicates that the optimal pH range for detecting these ions typically falls between 4.5 and 5.0, where both ions remain soluble and electroactive [49]. Furthermore, modified electrode surfaces at this pH facilitate efficient adsorption of these ions through favorable interactions with
functional groups on the electrode material, enhancing the electrochemical response and improving sensitivity [50].

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Optimization plot using desirability function.</figcaption></figure>

Under these ideal operating conditions, the RSM model estimated a response of 8.43 µA for Pb²⁺ ions and 3.11 µA for Cd²⁺ ions. According to Roslan et al. [51], a model is considered reliable when the desirability value approaches unity. The model's acceptability and applicability were confirmed by a desirability of 0.9095.

An experiment was carried out under the optimized conditions (pH = 4.8, DV = 11.73 µL, IT = 1 min, and DT = 3 h) to verify the predicted responses for Cd²⁺ and Pb²⁺ ions. To ensure accuracy, experimental errors were minimized, and the RSM model's predictions were validated, the experiment was repeated three times under these optimal conditions. As shown in Table 4, the response for Cd²⁺ ions is 1.9658 µA while the response for Pb²⁺ ions is 8.4699 µA under the optimized conditions.

<table><caption>**Table 4.** Results of experimental check of the predicted response for both ions at optimal conditions.</caption><thead><tr><th></th><th>pH</th><th>VD (µL)</th><th>IT (min)</th><th>DT (h)</th><th>I<sub>Cd</sub><sup>2+</sup> (µA)</th><th>I<sub>Pb</sub><sup>2+</sup> (µA)</th></tr></thead><tbody><tr><td>Predicted</td><td>4.8</td><td>11.73</td><td>1</td><td>3</td><td>3.11 ± 0.78</td><td>8.43 ± 2.32</td></tr><tr><td>Experimental</td><td>4.8</td><td>11.73</td><td>1</td><td>3</td><td>1.9658</td><td>8.4699</td></tr></tbody></table>

### 3.4. Electrochemical Determination of Cd²⁺ and Pb²⁺ Ions

Using the optimal experimental parameters recommended by the model, the distinct detection capability of the ZIF-7@PANI/GCE sensor for Pb²⁺ and Cd²⁺ ions was investigated through DPV measurements. As illustrated in Figure 8, the dissolution peaks current increased linearly with the rising concentration of the target heavy metal ions. The peak current for Pb²⁺ shows a linear correlation with concentrations ranging from 0.002 µM to 30 µM, following the regression equation $I_p (\mu A) = 0.8398 C (\mu M) + 0.6659 (R^2 = 0.9826)$. Similarly, Figure 8b,d depicts a linear calibration curve for Cd²⁺ over the
same concentration range, described by the equation $I_p$ (µA) = 0.5572 C (µM) + 0.7079 ($R^2$ = 0.9838). These results confirm the ZIF-7@PANI/GCE sensor's high sensitivity and effective detection capabilities for $Pb^{2+}$ and $Cd^{2+}$ ions. This enhanced performance is attributed to the significant specific surface area and structure of the ZIF-7 matrix, which provides numerous active sites for ion interaction. Additionally, the ZIF-7@PANI composite benefits from the synergistic effects of PANI and ZIF-7, where PANI facilitates the electron transport process and the abundant $-COO^-$ groups in ZIF-7 promote the pre-enrichment of $Pb^{2+}$ and $Cd^{2+}$ on the electrode surface.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** DPV curves of the ZIF-7@PANI/GCE sensor for the individual detection of (a) $Pb^{2+}$ and (c) $Cd^{2+}$ ions at different concentration ranges from 0.02 µM to 30 µM. Calibration curves of (b) $Pb^{2+}$ and (d) $Cd^{2+}$.</figcaption></figure>

A simultaneous analysis of $Cd^{2+}$ and $Pb^{2+}$ was conducted under the optimized conditions using the ZIF-7@PANI/GCE sensor, with a 1 min accumulation time. The tested ion concentrations ranged from 0 to 30 µM. As shown in Figure 9a, the DPV responses of $Pb^{2+}$ and $Cd^{2+}$ increased progressively with rising ion concentrations. Calibration curves for $Pb^{2+}$ and $Cd^{2+}$ were subsequently plotted, as illustrated in Figures 9b and 9c, respectively. For $Pb^{2+}$ detection, the analysis revealed two distinct sections, likely due to different adsorption behavior at low and high concentrations on the ZIF-7@PANI/GCE sensor surface. Notably, $Pb^{2+}$ ions exhibited higher sensitivity at lower concentrations, suggesting that metal ions initially interact with the most active sites, while at higher concentrations, they progressively occupy less active sites [52]. The regression equations for $Pb^{2+}$ were $I_p$ (µA) = 4.0865 C (µM) + 0.4339 ($R^2$ = 0.9847) within the 0.002–1 µM range and
Ip (µA) = 0.6955C (µM) + 5.1406 (R² = 0.9631) for 1–30µM (Figure 9b). For Cd²⁺, a single linear response was observed over the entire concentration range, following the equation Ip (µA) = 0.3344 C (µM) + 0.5412 (R² = 0.9821) from 0.02 µM to 30 µM (Figure 9c). The limits of detection (LOD) were calculated using the 3S/N method, yielding values of 2.93 nM for Pb²⁺ and 10.6 nM for Cd²⁺, confirming the sensor's high sensitivity and suitability for trace metal ion detection.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** (a) DPV curves of ZIF-7@PANI/GCE sensor for the simultaneous analysis of Pb²⁺ and Cd²⁺ across concentrations from 0.002 µM to 30 µM. The calibration curves of (b) Pb²⁺ and (c) Cd²⁺.</figcaption></figure>

The performance of the sensor ZIF-7@PANI/GCE was compared with other sensors for the simultaneous determination of Pb²⁺ and Cd²⁺ ions. The results, summarized in Table 5, demonstrate the good electroanalytical capabilities of the sensor developed in this study.
<table><caption>**Table 5.** Comparison of the electroanalytical performance of different sensors used for the simultaneous detection of Cd²⁺ and Pb²⁺ ions.</caption><thead><tr><th rowspan="2">Materials</th><th rowspan="2">Method</th><th colspan="2">LOD (nM)</th><th colspan="2">Linear Range (µM)</th><th rowspan="2">References</th></tr><tr><th>Pb<sup>2+</sup></th><th>Cd<sup>2+</sup></th><th>Pb<sup>2+</sup></th><th>Cd<sup>2+</sup></th></tr></thead><tbody><tr><td>CUiO-66/Bi/GCE</td><td>SWASV</td><td>5.98</td><td>41.57</td><td>0.178–1.067</td><td>0.178–1.067</td><td>[53]</td></tr><tr><td>Bi/GDY/GCE</td><td>DPASV</td><td>0.146</td><td>0.171</td><td>0.010–100.0</td><td>0.010–100.0</td><td>[54]</td></tr><tr><td>UiO-66-NH<sub>2</sub>/rGO</td><td>DPASV</td><td>16.39</td><td>-</td><td>0.2–3.2</td><td>-</td><td>[55]</td></tr><tr><td>UiO-66-NH<sub>2</sub>/GaOOH</td><td>DPASV</td><td>16</td><td>28</td><td>0.55–2.50</td><td>0.35–1.60</td><td>[56]</td></tr><tr><td>Pd NPs-porous carbon/GCE</td><td>DPV</td><td>50</td><td>41</td><td>0.5–8.9</td><td>0.5–5.5</td><td>[57]</td></tr><tr><td>ZIF-7@PANI/GCE</td><td>DPV</td><td>2.96</td><td>10.6</td><td>0.002–11–30</td><td>0.02–30</td><td>This work</td></tr></tbody></table>

SWASV: square wave stripping anodic voltammetry, GCE: glassy carbon electrode, DPASV: differential pulse anodic stripping voltammetry, CUiO-66: carbonized Zr-MOF, Bi: bismuth, GDY: graphdiyne, rGO: reduced graphene oxide, DPV: differential pulse voltammetry, Pd NPs: palladium nanoparticles.

### 3.5. Interference Study

In real water samples, multiple ions coexist, potentially interfering with detection results. To assess the anti-interference capacity of the ZIF-7@PANI/GCE sensor, interfering ions such as Ca²⁺, Co²⁺, Mg²⁺, Mn²⁺, Zn²⁺, and Cu²⁺ were introduced at concentrations five times higher than a 5 µM Cd²⁺ and Pb²⁺ solution. As depicted in Figure 10, the current response for both lead and cadmium ions remains largely unaffected in the presence of these ions, except for Cu²⁺. This substantial interference from Cu²⁺ can be attributed to its strong competition for adsorption sites on the electrode surface, which reduces the availability of these sites for Cd²⁺ ions. Notably, the fivefold higher concentration of Cu²⁺ resulted in a significant 19% decline in the response current for Cd²⁺. In contrast, the impact of other ions was relatively minor, with reductions ranging from 0.5% to 9% for Pb²⁺ and 3% to 8% for Cd²⁺. These findings confirm that the ZIF-7@PANI/GCE sensor exhibits a good anti-interference capacity, maintaining reliable detection even in the presence of high concentrations of interfering species.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** The specificity of the ZIF-7@PANI/GCE sensor toward 5 µM Cd²⁺ and Pb²⁺ against 50 µM Ca²⁺, Co²⁺, Mg²⁺, Mn²⁺, Zn²⁺, and Cu²⁺.</figcaption></figure>
### 3.6. Real Sample

To assess the real-world applicability of the proposed sensor and evaluate its electro-analytical performance, the ZIF-7@PANI/GCE electrochemical sensor was used to quantify Pb²⁺ and Cd²⁺ in real water samples, including tap water, mineral water, and seawater from the Monastir region. The water samples were filtered to remove impurities, and initial analyses detected no signals for Pb²⁺ and Cd²⁺, indicating their absence. To validate the sensor's performance, known concentrations of Pb²⁺ and Cd²⁺ ions were added to the samples, followed by electrochemical measurements using the ZIF-7@PANI/GCE sensor. The results, summarized in Table 6, revealed recovery rates ranging from 85% to 110.1%, confirming the sensor's accuracy and reliability for detecting Pb²⁺ and Cd²⁺ ions in environmental water samples.

<table><caption>**Table 6.** Determination of the content of Cd²⁺ and Pb²⁺ ions in tap water, mineral water, and seawater from the Monastir region using the ZIF-7@PANI/GCE sensor.</caption><thead><tr><th>Samples</th><th>Metal Ions</th><th>Added (μM)</th><th>Found (μM)</th><th>Recovery (%)</th></tr></thead><tbody><tr><td rowspan="4">Mineral water</td><td rowspan="2">Pb<sup>2+</sup></td><td>1</td><td>0.91</td><td>91</td></tr><tr><td>10</td><td>9.2</td><td>92</td></tr><tr><td rowspan="2">Cd<sup>2+</sup></td><td>1</td><td>1.01</td><td>101</td></tr><tr><td>10</td><td>10.91</td><td>109.1</td></tr><tr><td rowspan="4">Tap water</td><td rowspan="2">Pb<sup>2+</sup></td><td>1</td><td>0.85</td><td>85</td></tr><tr><td>10</td><td>9.3</td><td>93</td></tr><tr><td rowspan="2">Cd<sup>2+</sup></td><td>1</td><td>0.89</td><td>89</td></tr><tr><td>10</td><td>10.57</td><td>105.7</td></tr><tr><td rowspan="4">Sea water</td><td rowspan="2">Pb<sup>2+</sup></td><td>1</td><td>0.94</td><td>94</td></tr><tr><td>10</td><td>10.09</td><td>100.9</td></tr><tr><td rowspan="2">Cd<sup>2+</sup></td><td>1</td><td>1.088</td><td>108.8</td></tr><tr><td>10</td><td>11.01</td><td>110.1</td></tr></tbody></table>

### 4. Conclusions

This study presents a high-performance modified electrode based on ZIF-7@PANI for the simultaneous detection of Cd²⁺ and Pb²⁺ ions in water samples. The ZIF-7@PANI nanocomposite was designed to leverage the excellent electrical conductivity of polyaniline (PANI) alongside the enhanced ZIF-7 nanomaterials for trace metal ion detection. The sensor's performance was further optimized by tuning different experimental parameters, including accumulation time, deposited quantity, and drying time, using the Box-Behnken experimental design. Although pH influences sensor performance, it was controlled prior to measurements. The synergistic combination of PANI and ZIF-7 significantly improved the sensor's sensitivity, enabling the detection of Pb²⁺ and Cd²⁺ ions at ultra-trace levels, with detection limits of 2.96 nM for Pb²⁺ and 10.6 nM for Cd²⁺. The developed sensor demonstrated outstanding selectivity, achieving strong recovery rates (85–110.1%) in real water samples. In summary, the ZIF-7@PANI/GCE sensor proves to be a promising tool for heavy metal detection, offering a reliable and efficient solution for environmental monitoring and public health protection.

**Author Contributions:** Conceptualization, N.M. and F.E.; Methodology, S.G., C.M. and H.K.; Software, H.B. and H.K.; Validation, N.H.A. and Z.A.; Formal analysis, N.M. and C.M.; Investigation, S.G., H.B. and N.J.-R.; Resources, C.M. and F.E.; Data curation, N.M.; Writing—original draft, S.G. and Z.A.; Writing—review & editing, N.H.A., Z.A. and H.B.; Visualization, C.C.M.; Supervision, N.J.-R. and H.K.; Project administration, C.C.M.; Funding acquisition, C.C.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.
**Data Availability Statement:** Data are contained within the article.

**Acknowledgments:** The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through a Large Research Project under grant number RGP2/209/45. Authors thank the DGRS and Ministry of Higher Education and Scientific Research for partial support. Thanks to the Iraqi Ministry of Higher Education and Scientific Research as well as the Service for Cooperation and Cultural Action of the French Embassy in Iraq and Campus France for their support.

**Conflicts of Interest:** The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

## References

1. Emile, R.; Clammer, J.R.; Jayaswal, P.; Sharma, P. Addressing Water Scarcity in Developing Country Contexts: A Socio-Cultural Approach. *Humanit. Soc. Sci. Commun.* **2022**, *9*, 144. [CrossRef]

2. Kandasamy, G.; Singh, Y.R.; Narayanan, M. Impacts of Heavy Metal Enriched Tailings of Magnesite Mine on Surrounding Water Reservoirs. *J. Environ. Treat. Tech.* **2021**, *9*, 521–527. [CrossRef]

3. Kamba, F.; Sangija, F.; Wei, S. The Poor Sanitation and Access to Clean Water in Rural Areas: Case of Bossangoa (Central African Republic). *Adv. Soc. Sci. Res. J.* **2016**, *3*. [CrossRef]

4. Ali, H.; Khan, E.; Ilahi, I. Environmental Chemistry and Ecotoxicology of Hazardous Heavy Metals: Environmental Persistence, Toxicity, and Bioaccumulation. *J. Chem.* **2019**, *2019*, 6730305. [CrossRef]

5. You, J.; Li, J.; Wang, Z.; Baghayeri, M.; Zhang, H. Application of Co₃O₄ Nanocrystal/rGO for Simultaneous Electrochemical Detection of Cadmium and Lead in Environmental Waters. *Chemosphere* **2023**, *335*, 139133. [CrossRef] [PubMed]

6. Yuan, M.; Peng, X.; Ge, F.; Li, Q.; Wang, K.; Yu, D.-G.; Wang, Z. Simplified Design for Solution Anode Glow Discharge Atomic Emission Spectrometry Device for Highly Sensitive Detection of Ag, Bi, Cd, Hg, Pb, Tl, and Zn. *Microchem. J.* **2020**, *155*, 104785. [CrossRef]

7. Ling, Y.; Luo, F.; Zhu, S. A Simple and Fast Sample Preparation Method Based on Ionic Liquid Treatment for Determination of Cd and Pb in Dried Solid Agricultural Products by Graphite Furnace Atomic Absorption Spectrometry. *LWT* **2021**, *142*, 111077. [CrossRef]

8. Jurowski, K. The Toxicological Assessment of Hazardous Elements (Pb, Cd and Hg) in Low-Cost Jewelry for Adults from Chinese E-Commerce Platforms: In Situ Analysis by Portable X-Ray Fluorescence Measurement. *J. Hazard. Mater.* **2023**, *460*, 132167. [CrossRef]

9. Chen, Y.; He, M.; Chen, B.; Hu, B. Thiol-Grafted Magnetic Polymer for Preconcentration of Cd, Hg, Pb from Environmental Water Followed by Inductively Coupled Plasma Mass Spectrometry Detection. *Spectrochim. Acta Part B At. Spectrosc.* **2021**, *177*, 106071. [CrossRef]

10. Zhang, L.; Zhen, Y.; Tong, L. Optical Micro/Nanofiber Enabled Tactile Sensors and Soft Actuators: A Review. *Opto-Electronic Sci.* **2024**, *3*, 240005–240017. [CrossRef]

11. Dahaghin, Z.; Kilmartin, P.A.; Mousavi, H.Z. Novel Ion Imprinted Polymer Electrochemical Sensor for the Selective Detection of Lead(II). *Food Chem.* **2020**, *303*, 125374. [CrossRef] [PubMed]

12. Sengupta, P.; Pramanik, K.; Sarkar, P. Simultaneous Detection of Trace Pb(II), Cd(II) and Hg(II) by Anodic Stripping Analyses with Glassy Carbon Electrode Modified by Solid Phase Synthesized Iron-Aluminate Nano Particles. *Sens. Actuators B Chem.* **2021**, *329*, 129052. [CrossRef]

13. Pasakon, P.; Mensing, J.P.; Phokaratkul, D.; Karuwan, C.; Lomas, T.; Wisitsoraat, A.; Tuantranont, A. A High-Performance, Disposable Screen-Printed Carbon Electrode Modified with Multi-Walled Carbon Nanotubes/Graphene for Ultratrace Level Electrochemical Sensors. *J. Appl. Electrochem.* **2019**, *49*, 217–227. [CrossRef]

14. Adil, H.I.; Thalji, M.R.; Yasin, S.A.; Saeed, I.A.; Assiri, M.A.; Chong, K.F.; Ali, G.A.M. Metal–Organic Frameworks (MOFs) Based Nanofiber Architectures for the Removal of Heavy Metal Ions. *RSC Adv.* **2022**, *12*, 1433–1450. [CrossRef] [PubMed]

15. Mansour, M.; Kahri, H.; Guergueb, M.; Barhoumi, H.; Gutierrez Puebla, E.; Ayed, B.; Demirci, U.B. Copper-Based MOF, Cu₃(SDBA)₂ (HSDBA), as a Catalyst for Efficient Reduction of 4-Nitrophenol in the Presence of Sodium Borohydride. *React. Chem. Eng.* **2022**, *7*, 908–916. [CrossRef]

16. Missaoui, N.; Chrouda, A.; Kahri, H.; Gross, A.J.; Rezaei Ardani, M.; Pang, A.L.; Ahmadipour, M. PEG-Templated Synthesis of Ultramicroporous n-ZIF-67 Nanoparticles with High Selectivity for the Adsorption and Uptake of CO2 over CH4 and N2. *Sep. Purif. Technol.* **2023**, *316*, 123755. [CrossRef]
17. Liu, S.; Hu, Z.; Yan, X.; Geng, S.; Zhao, X.; Ou, H.; Pan, J. Separation and Purification of Target Flavonoids Using Covalently Connected MOFs@boronic Acid-Functionalized-COFs Magnetic Hybrids: Precise Identification and Enhanced Stability. *Sep. Purif. Technol.* **2023**, *320*, 124061. [CrossRef]

18. Ma, D.; Wang, G.; Lu, J.; Zeng, X.; Cheng, Y.; Zhang, Z.; Lin, N.; Chen, Q. Multifunctional Nano MOF Drug Delivery Platform in Combination Therapy. *Eur. J. Med. Chem.* **2023**, *261*, 115884. [CrossRef] [PubMed]

19. Musarurwa, H.; Tavengwa, N.T. Smart Metal-Organic Framework (MOF) Composites and Their Applications in Environmental Remediation. *Mater. Today Commun.* **2022**, *33*, 104823. [CrossRef]

20. Shafqat, S.S.; Rizwan, M.; Batool, M.; Shafqat, S.R.; Mustafa, G.; Rasheed, T.; Zafar, M.N. Metal Organic Frameworks as Promising Sensing Tools for Electrochemical Detection of Persistent Heavy Metal Ions from Water Matrices: A Concise Review. *Chemosphere* **2023**, *318*, 137920. [CrossRef]

21. Wannassi, J.; Missaoui, N.; Mabrouk, C.; Barhoumi, H.; Crapnell, R.D.; Jaffrezic, N.; Banks, C.E.; Kahri, H. Electrochemical Sensors Based on Metal-Organic Framework and Conductive Polymer HKUST-1@PANI for High-Performance Detection of Lead Ions. *J. Electrochem. Soc.* **2023**, *170*, 117504. [CrossRef]

22. Dong, J.; Zhang, D.; Li, C.; Bai, T.; Jin, H.; Suo, Z. A Sensitive Electrochemical Sensor Based on PtNPs@Cu-MOF Signal Probe and DNA Walker Signal Amplification for Pb2+ Detection. *Bioelectrochemistry* **2022**, *146*, 108134. [CrossRef] [PubMed]

23. Wang, C.; Niu, Q.; Liu, D.; Dong, X.; You, T. Electrochemical Sensor Based on Bi/Bi2O3 Doped Porous Carbon Composite Derived from Bi-MOFs for Pb2+ Sensitive Detection. *Talanta* **2023**, *258*, 124281. [CrossRef]

24. Tran, T.Q.N.; Das, G.; Yoon, H.H. Nickel-Metal Organic Framework/MWCNT Composite Electrode for Non-Enzymatic Urea Detection. *Sens. Actuators B Chem.* **2017**, *243*, 78–83. [CrossRef]

25. Guesmi, S.; Moulaee, K.; Bressi, V.; Kahri, H.; Khaskhoussi, A.; Espro, C.; Barhoumi, H.; Neri, G. Non-Enzymatic Amperometric Glucose Sensing by Novel Cu-MOF Synthesized at Room Temperature. *Mater. Adv.* **2024**, *5*, 1160–1170. [CrossRef]

26. Wang, Y.; Wang, L.; Huang, W.; Zhang, T.; Hu, X.; Perman, J.A.; Ma, S. A Metal–Organic Framework and Conducting Polymer Based Electrochemical Sensor for High Performance Cadmium Ion Detection. *J. Mater. Chem. A* **2017**, *5*, 8385–8393. [CrossRef]

27. Shahsavari, M.; Mohammadzadeh Jahani, P.; Sheikhshoaie, I.; Tajik, S.; Aghaei Afshar, A.; Askari, M.B.; Salarizadeh, P.; Di Bartolomeo, A.; Beitollahi, H. Green Synthesis of Zeolitic Imidazolate Frameworks: A Review of Their Characterization and Industrial and Medical Applications. *Materials* **2022**, *15*, 447. [CrossRef] [PubMed]

28. Ban, Y.; Li, Y.; Liu, X.; Peng, Y.; Yang, W. Solvothermal Synthesis of Mixed-Ligand Metal–Organic Framework ZIF-78 with Controllable Size and Morphology. *Microporous Mesoporous Mater.* **2013**, *173*, 29–36. [CrossRef]

29. Hunter-Sellars, E.; Saenz-Cavazos, P.A.; Houghton, A.R.; McIntyre, S.R.; Parkin, I.P.; Williams, D.R. Sol–Gel Synthesis of High-Density Zeolitic Imidazolate Framework Monoliths via Ligand Assisted Methods: Exceptional Porosity, Hydrophobicity, and Applications in Vapor Adsorption. *Adv. Funct. Mater.* **2021**, *31*, 2008357. [CrossRef]

30. Yan, X.; Yang, Y.; Wang, C.; Hu, X.; Zhou, M.; Komarneni, S. Surfactant-Assisted Synthesis of ZIF-8 Nanocrystals for Phthalic Acid Adsorption. *J. Sol-Gel Sci. Technol.* **2016**, *80*, 523–530. [CrossRef]

31. Lai, L.S.; Yeong, Y.F.; Lau, K.K.; Shariff, A.M. Effect of Synthesis Parameters on the Formation of ZIF-8 Under Microwave-Assisted Solvothermal. *Procedia Eng.* **2016**, *148*, 35–42. [CrossRef]

32. Missaoui, N.; Kahri, H.; Demirci, U.B. Rapid Room-Temperature Synthesis and Characterizations of High-Surface-Area Nanoparticles of Zeolitic Imidazolate Framework-8 (ZIF-8) for CO2 and CH4 Adsorption. *J. Mater. Sci.* **2022**, *57*, 16245–16257. [CrossRef]

33. Missaoui, N.; Mohamed, B.; Chrouda, A.; Kahri, H.; Barhoumi, H.; Pang, A.L.; Ahmadipour, M. Interpreting of the Carbon Dioxide Adsorption on High Surface Area Zeolitic Imidazolate Framework-8 (ZIF-8) Nanoparticles Using a Statistical Physics Model. *Microporous Mesoporous Mater.* **2023**, *360*, 112711. [CrossRef]

34. He, M.; Yao, J.; Li, L.; Wang, K.; Chen, F.; Wang, H. Synthesis of Zeolitic Imidazolate Framework-7 in a Water/Ethanol Mixture and Its Ethanol-Induced Reversible Phase Transition. *ChemPlusChem* **2013**, *78*, 1222–1225. [CrossRef] [PubMed]

35. Kang, C.-H.; Lin, Y.-F.; Huang, Y.-S.; Tung, K.-L.; Chang, K.-S.; Chen, J.-T.; Hung, W.-S.; Lee, K.-R.; Lai, J.-Y. Synthesis of ZIF-7/Chitosan Mixed-Matrix Membranes with Improved Separation Performance of Water/Ethanol Mixtures. *J. Membr. Sci.* **2013**, *438*, 105–111. [CrossRef]

36. He, M.; Yao, J.; Liu, Q.; Zhong, Z.; Wang, H. Toluene-Assisted Synthesis of RHO-Type Zeolitic Imidazolate Frameworks: Synthesis and Formation Mechanism of ZIF-11 and ZIF-12. *Dalton Trans.* **2013**, *42*, 16608–16613. [CrossRef]

37. Hu, Y.; Kazemian, H.; Rohani, S.; Huang, Y.; Song, Y. In Situ High Pressure Study of ZIF-8 by FTIR Spectroscopy. *Chem. Commun.* **2011**, *47*, 12694–12696. [CrossRef] [PubMed]

38. Wang, S.; Ma, L.; Gan, M.; Fu, S.; Dai, W.; Zhou, T.; Sun, X.; Wang, H.; Wang, H. Free-Standing 3D Graphene/Polyaniline Composite Film Electrodes for High-Performance Supercapacitors. *J. Power Sources* **2015**, *299*, 347–355. [CrossRef]

39. Dhibar, S.; Das, C.K. Electrochemical Performances of Silver Nanoparticles Decorated Polyaniline/Graphene Nanocomposite in Different Electrolytes. *J. Alloys Compd.* **2015**, *653*, 486–497. [CrossRef]

40. Park, K.S.; Ni, Z.; Côté, A.P.; Choi, J.Y.; Huang, R.; Uribe-Romo, F.J.; Chae, H.K.; O’Keeffe, M.; Yaghi, O.M. Exceptional Chemical and Thermal Stability of Zeolitic Imidazolate Frameworks. *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 10186–10191. [CrossRef]
41. Yan, J.; Wei, T.; Shao, B.; Fan, Z.; Qian, W.; Zhang, M.; Wei, F. Preparation of a Graphene Nanosheet/Polyaniline Composite with High Specific Capacitance. *Carbon* **2010**, *48*, 487–493. [CrossRef]

42. Butler, K.T.; Worrall, S.D.; Molloy, C.D.; Hendon, C.H.; Attfield, M.P.; Dryfe, R.A.W.; Walsh, A. Electronic Structure Design for Nanoporous, Electrically Conductive Zeolitic Imidazolate Frameworks. *J. Mater. Chem. C* **2017**, *5*, 7726–7731. [CrossRef]

43. Chelladurai, S.J.S.; Murugan, K.; Ray, A.P.; Upadhyaya, M.; Narasimharaj, V.; Gnanasekaran, S. Optimization of Process Parameters Using Response Surface Methodology: A Review. *Mater. Today Proc.* **2021**, *37*, 1301–1304. [CrossRef]

44. Pinheiro, D.R.; Neves, R.d.F.; Paz, S.P.A. A Sequential Box-Behnken Design (BBD) and Response Surface Methodology (RSM) to Optimize SAPO-34 Synthesis from Kaolin Waste. *Microporous Mesoporous Mater.* **2021**, *323*, 111250. [CrossRef]

45. Puthongkham, P.; Wirojsaengthong, S.; Suea-Ngam, A. Machine Learning and Chemometrics for Electrochemical Sensors: Moving Forward to the Future of Analytical Chemistry. *Analyst* **2021**, *146*, 6351–6364. [CrossRef] [PubMed]

46. Bressi, V.; Celesti, C.; Ferlazzo, A.; Len, T.; Moulaee, K.; Neri, G.; Luque, R.; Espro, C. Waste-Derived Carbon Nanodots for Fluorimetric and Simultaneous Electrochemical Detection of Heavy Metals in Water. *Environ. Sci. Nano* **2024**, *11*, 1245–1258. [CrossRef]

47. Okpara, E.C.; Fayemi, O.E.; Wojuola, O.B.; Onwudiwe, D.C.; Ebenso, E.E. Electrochemical Detection of Selected Heavy Metals in Water: A Case Study of African Experiences. *RSC Adv.* **2022**, *12*, 26319–26361. [CrossRef]

48. Meng, R.; Zhu, Q.; Long, T.; He, X.; Luo, Z.; Gu, R.; Wang, W.; Xiang, P. The Innovative and Accurate Detection of Heavy Metals in Foods: A Critical Review on Electrochemical Sensors. *Food Control.* **2023**, *150*, 109743. [CrossRef]

49. Keramari, V.; Karastogianni, S.; Girousi, S. New Prospects in the Electroanalysis of Heavy Metal Ions (Cd, Pb, Zn, Cu): Development and Application of Novel Electrode Surfaces. *Methods Protoc.* **2023**, *6*, 60. [CrossRef] [PubMed]

50. Bhatia, N.; Mehta, D.; Garg, K.; Kaur, S.; Nagaiah, T.C. Employing Mesoporous Nitrogen Containing Carbon for Simultaneous Electrochemical Detection of Heavy Metal Ions. *Chem. Asian J.* **2024**, e202401315. [CrossRef]

51. Roslan, J.; Mustapa Kamal, S.; Yunos, K.; Abdullah, N. Optimization of Enzymatic Hydrolysis of Tilapia Muscle (Oreochromis Niloticus) Using Response Surface Methodology (RSM). *Sains Malays.* **2014**, *43*, 1715–1723.

52. Moulaee, K.; Raza, M.H.; Pinna, N.; Donato, N.; Neri, G. CNT/Al₂O₃ Core–Shell Nanostructures for the Electrochemical Detection of Dihydroxybenzene Isomers. *Phys. Chem. Chem. Phys.* **2021**, *23*, 14064–14074. [CrossRef] [PubMed]

53. Ding, Y.; Wei, F.; Dong, C.; Li, J.; Zhang, C.; Han, X. UiO-66 Based Electrochemical Sensor for Simultaneous Detection of Cd(II) and Pb(II). *Inorg. Chem. Commun.* **2021**, *131*, 108785. [CrossRef]

54. Ai, Y.; Yan, L.; Zhang, S.; Ye, X.; Xuan, Y.; He, S.; Wang, X.; Sun, W. Ultra-Sensitive Simultaneous Electrochemical Detection of Zn(II), Cd(II) and Pb(II) Based on the Bismuth and Graphdiyne Film Modified Electrode. *Microchem. J.* **2023**, *184*, 108186. [CrossRef]

55. Li, Z.; Zhang, X.; Luo, Y.; Li, Q.; Qin, Y.; Wang, G.; Yang, S.; Liu, Z. An Electrochemical Sensor Based on the Composite UiO-66-NH2/rGO for Trace Detection of Pb(II) and Cu(II). *Chem. Phys. Lett.* **2023**, *830*, 140825. [CrossRef]

56. Ru, J.; Wang, X.; Cui, X.; Wang, F.; Ji, H.; Du, X.; Lu, X. GaOOH-Modified Metal-Organic Frameworks UiO-66-NH2: Selective and Sensitive Sensing Four Heavy-Metal Ions in Real Wastewater by Electrochemical Method. *Talanta* **2021**, *234*, 122679. [CrossRef] [PubMed]

57. Veerakumar, P.; Veeramani, V.; Chen, S.-M.; Madhu, R.; Liu, S.-B. Palladium Nanoparticle Incorporated Porous Activated Carbon: Electrochemical Detection of Toxic Metal Ions. *ACS Appl. Mater. Interfaces* **2016**, *8*, 1319–1326. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.