Graphite Studded with Facile-Synthesized Cu 2 O Nanoparticle-Based Cubes as a Novel Electrochemical Sensor for Highly Sensitive Voltametric Determination of Mebeverine Hydrochloride

: Herein, a feasible chemical reduction method followed by intensive mixing was applied for the preparation of an attractive material based on graphite studded with cuprous oxide nanoparticle-based cubes (Cu 2 ONPs–C@G). Transmission electron microscope (TEM), scanning electron microscope (SEM), X-ray diffraction (XRD) and cyclic voltammetry (CV) were utilized for characterization. Cuprous oxide nanoparticles (Cu 2 ONPs), with a diameter range mainly distributed from 4 to 20 nm, aggregate to form microcubes (Cu 2 ONPs–C) with an average diameter of about 367 nm. Paste electrode was prepared using Cu 2 ONPs–C@G (Cu 2 ONPs–C@G/PE) for voltametric quantiﬁcation of the musculotropic antispasmodic drug: mebeverine hydrochloride (MEB). The electrochemical behavior of MEB was studied using CV, and the optimum analytical parameters were investigated using square wave adsorptive anodic stripping voltammetry (SWAdASV). Moreover, density functional theory (DFT) was used to emphasize the ability of MEB to form a complex with Cu 2+ , conﬁrming the suggested electrochemical behavior of MEB at Cu 2 ONPs–C@G/PE. With good stability and high reproducibility, SWAdASV of Cu 2 ONPs–C@G/PE shows successful quantiﬁcation of MEB over the concentration range of 5.00 × 10 − 11 –1.10 × 10 − 9 M with lower limit of detection (LOD) and lower limit of quantiﬁcation (LOQ) values of 2.41 × 10 − 11 M and 8.05 × 10 − 11 M, respectively. Finally, accurate quantiﬁcation of MEB in dosage forms (tablets) and biological ﬂuids (spiked human urine and plasma samples) was achieved using Cu 2 ONPs-C@G/PE.


Introduction
Nanoparticles offer interesting and peculiar properties, such as size-and shapedependent and catalytic activities, crystal shape-dependent thermodynamics, novel mechanical properties, unique magnetic phenomena and quantum confinement phenomena [1]. Particularly copper nanoparticles show attractive unusual properties, such as small dimensions, large specific surface area, low electric resistance and lower cost than other favored metal nanoparticles, e.g., gold and silver nanoparticles [2][3][4]. With special regard to cuprous oxide nanoparticles (Cu 2 ONPs), it is considered as one of the cheapest and abundant p-type semiconductor metal oxides with a narrow band gap and is environmentally friendly, since it is a relatively non-toxic material [5].
Among cuprous oxide transversal forms, cubes may be the most interesting because of their tailor-made properties and their high symmetries, other more complex structures the cheapest and abundant p-type semiconductor metal oxides with a narrow band gap and is environmentally friendly, since it is a relatively non-toxic material [5].
Among cuprous oxide transversal forms, cubes may be the most interesting because of their tailor-made properties and their high symmetries, other more complex structures can be obtained from these simple structural forms [6,7]. It was reported about the self-assembly of Cu2ONPs to form microspheres and microcubes [8]. Different synthesis methods were applied to obtain the microcube morphology, such as hydrothermal [8,9], metal-organic framework calcination [10], surface ion exchange strategy followed by reductive deposition [11] and chemical reduction [12]. However, it is important to introduce a more facile and safe synthesis method where these previously mentioned methods require special equipment [8,9]; use special chemicals [9,[10][11][12], such as surfactants; or use dangerous chemicals [11].
Carbon-based metal nanoparticles have attracted considerable attention due to their mechanical and thermal properties and their ability for electrode modifications [13,14]. Because of these properties, metal and metal oxide nanoparticles have been used successfully for increasing carbon paste electrode sensitivity toward various electroanalytical applications [15,16]. It was found in the literature that pure Cu2ONPs [11] and Cu2ONPs supported on carbon materials were applied successfully in different voltametric determinations [17,18].
Mebeverine hydrochloride (MEB) (Scheme 1) is used in the treatment of various medicinal abdominal troubles. MEB is an effective musculotropic antispasmodic, which is used as a medicine for spasms associated with gastrointestinal disorders, such as irritable bowel syndrome and gastrointestinal spasm secondary [19,20]. MEB relieves spasms without affecting gut motility due to its action not interfering with the autonomic nervous system. The ideal dose of MEB is to be taken in the range of 135-270 mg, and this is to avoid serious health complications in the form of cardiovascular and neurological symptoms. After oral administration, MEB is rapidly and completely metabolized to mebeverine alcohol and veratric acid. Moreover, as a result of MEB hydrolysis and desmethylation steps, mebeverine acid and desmethylmebeverine acid are the major metabolites of MEB found in real human samples [20,21]. Due to its vital importance, MEB quantification in dosage forms and biological fluids at low concentrations is of great importance. Although different analytical methods have been applied for the quantification of MEB, there is an agreement that MEB cannot be detected in real biological samples due to the previously mentioned MEB metabolism. Additionally, most of these methods have been applied for the MEB quantification in dosage forms and spiked biological samples. The literature survey reveals several analytical methods for MEB quantification, such as chromatographic [21,22], GC-MS [23], spectrophotometric [24,25], spectrofluorometric [25,26] and electrochemical methods [27][28][29][30][31][32][33]. Carbon paste electrodes (CPEs) have been applied successfully for MEB quantification either using potentiometric [31] or voltametric techniques [32,33]. This can be attributed to their low cost, low background currents, possibility to be used over a wide Due to its vital importance, MEB quantification in dosage forms and biological fluids at low concentrations is of great importance. Although different analytical methods have been applied for the quantification of MEB, there is an agreement that MEB cannot be detected in real biological samples due to the previously mentioned MEB metabolism. Additionally, most of these methods have been applied for the MEB quantification in dosage forms and spiked biological samples. The literature survey reveals several analytical methods for MEB quantification, such as chromatographic [21,22], GC-MS [23], spectrophotometric [24,25], spectrofluorometric [25,26] and electrochemical methods [27][28][29][30][31][32][33]. Carbon paste electrodes (CPEs) have been applied successfully for MEB quantification either using potentiometric [31] or voltametric techniques [32,33]. This can be attributed to their low cost, low background currents, possibility to be used over a wide studied potential window, ease of surface renewability as well as modification simplicity and utilization variety when compared with other solid electrodes [34][35][36][37]. However, electrodes used for reported voltametric techniques are prepared from high price materials with multistep time-consuming preparation, with a sensitivity that needs improvement [32,33].
Quantum chemical calculations are widely used to predict the optimized structure of a compound and its reactivity to bind the metal ion. Recently, with the massive advancement of computational hardware, density functional theory (DFT) has attracted considerable attention from the scientific community for predicting the structural behavior and reaction mechanisms of several chemical species. Using this advanced computational method, it is easy to predict a systematic method for the analysis and describe the reaction mechanism of different organic compounds and their metal complexes [38,39].
The current work is a continuation of our efforts in the field of voltametric quantification of pharmaceuticals using bare [34,40] and modified [41,42] electrodes. Herein, we propose a simple and straightforward strategy for the synthesis of cuprous oxide microcubes (Cu 2 ONPs-C). The resulting Cu 2 ONPs-C is used for the studding of graphite (Cu 2 ONPs-C@G) via a simple hand mixing method. Both Cu 2 ONPs-C and Cu 2 ONPs-C@G are characterized using different structural and electrochemical techniques. Moreover, Cu 2 ONPs-C@G is used to prepare a working paste electrode, which is denoted as Cu 2 ONPs-C@G/PE. The Cu 2 ONPs-C@G/PE is applied to study the MEB electrochemical behavior using cyclic voltammetry (CV). DFT is used for investigating the behavior of MEB during the voltametric measurement. Thanks to its high sensitivity, square wave adsorptive anodic stripping voltammetry (SWAdASV) is used to investigate the optimum conditions required for the voltametric quantification of MEB. The proposed Cu 2 ONPs-C@G/PE is applied over a wide potential window (−0.6 to +1.0 V vs. Ag/AgCl) in an easier and highly sensitive procedure using SWAdASV for the quantification of MEB in dosage forms (tablets) and biological fluids (spiked human urine and plasma samples).

Apparatus
All voltametric experiments were carried out using the EG&G Princeton Applied Research Potentiostat/Galvanostat Model 273A (USA) combined with a magnetic stirrer (KIKA Labortechinik, Staufen, Germany), controlled by an IBM compatible computer running the electroanalytical software model 270/250 version 4.30. This system was accomplished with a three-electrode system where handmade carbon paste, Ag/AgCl and Pt wire were used as the working, reference, and counter electrodes, respectively. VWR Scientific Products, Model 2000, USA, were used for all pH measurements. The used pH meter was calibrated using standard buffer pH 4.01, 7.00 and 9.00.

Preparation and Characterization of the Cu 2 ONPs-C and Cu 2 ONPs-C@G Composites
Herein, a previous procedure for cupric oxide nanoparticles (CuONPs) preparation was applied with some modifications that resulted in Cu 2 ONPs-C production as we prove later in the discussion of the characterization results [4]. In a typical procedure, 200 mg of Cu(NO 3 ) 2 ·3H 2 O was dissolved in 100 mL distilled water. Then, at room temperature and under vigorous stirring for 10 min, 4 mL of 10% N 2 H 4 ·H 2 O was added dropwise to the previous resulted solution. The addition process was accompanied by a rapid formation of dark blue suspension. The whole solution was heated at 80 • C until complete evaporation. During heating, the dark blue suspension is changed gradually to dark red and finally stabilized at reddish-brown color. The resulted precipitate after drying is characterized as a fine precipitate; it was collected, washed repeatedly with distilled water to ensure the highest purity degrees, then filtered and left till complete drying at room temperature. The as-prepared reddish-brown precipitate was collected and labeled as Cu 2 ONPs-C.
For the preparation of Cu 2 ONPs-C@G, 1 g of graphite was hand-mixed intensely with 80 mg of Cu 2 ONPs-C. To ensure the good distribution and stabilization of Cu 2 ONPs-C on the graphite grain surface, an Eppendorf tube containing this mixture was ultrasonicated for 10 min. A small part of the resulting mixture was kept for the characterization process, and the rest was used in the preparation of a paste working electrode (Cu 2 ONPs-C@G/PE) as previously mentioned [34]. For comparison purposes, another paste electrode based on bare graphite (CPE) was prepared using the same procedure. For the investigation of the prepared materials, different analyses were accomplished. The phase of the prepared nanoparticles was identified by recording X-ray diffraction (XRD) patterns using the X-Ray Diffractometer theta:theta D8 Advance-Burker, with Cu-K α radiation operated at λ = 0.1542 nm, 40 kv and 40 mA, with a secondary monochromator of graphite at an angle of 2θ from 10 • to 80 • with a scan speed of 0.02 • /min. The crystallite sizes of the prepared nanoparticles were calculated by application of Scherrer's equation. Morphology of the prepared nanoparticles was investigated using avscanning electron microscope (SEM) studies carried out on JEOL JSM 5400 LV. To examine the inner structure and size of the obtained Cu 2 ONPs-C@G in the composite, transmission electron micrographs were recorded using a JEOL JEM-100 CXII transmission electron microscope (TEM) equipped with a tungsten-filament 100 kv transmission electron microscope; the images were picked out using a Gatan Erlangshen ES500W camera where Gatan Digital Micrograph software was used for the determination of nanoparticle diameters.

Computational Calculations
The theoretical calculations were performed using the DFT at B3LYP (Becke's threeparameter hybrid functional using the correlation functional of Lee, Yang and Parr, which includes both local and non-local term correlation functionals) methods at 3-21G basis set. To perform this computational work, we used the Gaussian 09 program package [43].

Reagents and Materials
Copper nitrate trihydrate (Cu(NO 3 ) 2 ·3H 2 O) and hydrazine monohydrate (N 2 H 4 ·H 2 O) were purchased from Avondale Laboratories (UK) and Alpha Chemika (India), respectively. MEB was purchased from Merck. An appropriate amount of pure MEB was used to prepare MEB stock solution by dissolving in double-distilled water, and the resulting solution was stored in a brown volumetric flask in a refrigerator. Required concentrations of MEB were prepared daily by serial dilution of its stock solution using the same solvent. For MEB quantification in dosage forms, Duspatalin and Colospasmin ® forte tablets were purchased from Pharco Pharmaceutical Company and Egyptian International Pharmaceutical Industries Company (E.I.P.I.CO), respectively. Each tablet is claimed to contain 135 mg MEB. All other used chemicals were of analytical grade, used without further purification, and their solutions were prepared using double distilled water. Graphite powder was purchased from Measurement Ringsdorff Werke, Germany. For all stripping voltametric measurements, the used blank solution was a mixture of 0.1 M phosphate buffer (pH ≈ 8.0) and 0.1 M KCl. The preconcentration step was performed by immersing the prepared Cu 2 ONPs-C@G/PE in a stirred 15 mL of the blank solution at a potential range from −0.6 to +1.0 V for a given period. Once the stirring stopped, the solution was left for 10 s to settle and allow the background current to decrease. The same procedure was repeated in the presence of 2 × 10 −5 M MEB. Then, cyclic and square wave anodic stripping voltammograms were recorded in the positive potential direction. To restore the initial activity of Cu 2 ONPs-C@G/PE, simply removing the outer surface of paste by smoothing its surface against white paper was performed. All measurements were performed at room temperature. A brief description of the working electrode construction and MEB voltametric quantification is shown in Figure 1.
crease. The same procedure was repeated in the presence of 2 × 10 −5 M MEB. Then, cyclic and square wave anodic stripping voltammograms were recorded in the positive potential direction. To restore the initial activity of Cu2ONPs-C@G/PE, simply removing the outer surface of paste by smoothing its surface against white paper was performed. All measurements were performed at room temperature. A brief description of the working electrode construction and MEB voltametric quantification is shown in Figure 1. For analysis of MEB in dosage forms, ten tablets of MEB dosage forms (Duspatalin or Colospasmin ® forte) were weighed, powdered and mixed well in a mortar. An amount of 0.1294 g of powdered MEB tablets was dissolved in 100 mL of double distilled water, shaken well and filtered using filter paper. An aliquot of the filtrate was then transferred into a 10 mL calibrated flask, and it was completed to the mark with the same solvent. Then, 20 µL of this solution was then added into the voltametric cell, and the voltametric procedure was performed as before.

MEB Analysis from Spiked Human Urine and Plasma Samples
The studied biological fluids were collected from healthy donors who had not received any previous treatments. These fluids were subjected to simple preliminary treatment steps before used in MEB determination.
For the MEB quantification spiked human urine samples, the collected blank urine samples were centrifuged for 20 min at 4000 rpm to get rid of protein residues. Then, 1.0 mL of the resulting supernatant solution was drawn off using a micro-pipette and then injected into an electrochemical cell containing a blank solution mixture containing 0.1 M phosphate buffer (pH ≈ 8.0) and 0.1 M KCl. Thereafter, the described voltametric procedure was conducted, and MEB determination voltammogram was recorded.
The investigated plasma samples were obtained by centrifuging the collected blood samples at 4000 rpm for 30 min using tabletop the high-speed universal harmonic centrifuge 12,000 rpm model PLC-012 (Taiwan). After that, a mixture consisting of MEB aliquots, 1.0 mL of obtained human plasma and 3 mL of acetonitrile was transferred quantitatively to centrifugation tubes and allowed to be centrifuged at 4000 rpm for 20 min. For analysis of MEB in dosage forms, ten tablets of MEB dosage forms (Duspatalin or Colospasmin ® forte) were weighed, powdered and mixed well in a mortar. An amount of 0.1294 g of powdered MEB tablets was dissolved in 100 mL of double distilled water, shaken well and filtered using filter paper. An aliquot of the filtrate was then transferred into a 10 mL calibrated flask, and it was completed to the mark with the same solvent. Then, 20 µL of this solution was then added into the voltametric cell, and the voltametric procedure was performed as before.

MEB Analysis from Spiked Human Urine and Plasma Samples
The studied biological fluids were collected from healthy donors who had not received any previous treatments. These fluids were subjected to simple preliminary treatment steps before used in MEB determination.
For the MEB quantification spiked human urine samples, the collected blank urine samples were centrifuged for 20 min at 4000 rpm to get rid of protein residues. Then, 1.0 mL of the resulting supernatant solution was drawn off using a micro-pipette and then injected into an electrochemical cell containing a blank solution mixture containing 0.1 M phosphate buffer (pH ≈ 8.0) and 0.1 M KCl. Thereafter, the described voltametric procedure was conducted, and MEB determination voltammogram was recorded.
The investigated plasma samples were obtained by centrifuging the collected blood samples at 4000 rpm for 30 min using tabletop the high-speed universal harmonic centrifuge 12,000 rpm model PLC-012 (Taiwan). After that, a mixture consisting of MEB aliquots, 1.0 mL of obtained human plasma and 3 mL of acetonitrile was transferred quantitatively to centrifugation tubes and allowed to be centrifuged at 4000 rpm for 20 min. After finishing the centrifugation, 20 µL of the resulting mixture was injected into the voltametric cell, and the recommended voltametric procedure was applied for the MEB voltametric determination.  [45]. The XRD pattern of Cu 2 ONPs-C proved that there are no characteristic peaks of either Cu or CuO nanoparticles, indicating that the obtained Cu 2 ONPs are a pure crystalline single cubic phase with typical face-centered cubic (fcc) lattice [3,45]. Approximately, all diffraction peaks of bare graphite grains and Cu 2 ONPs-C were recorded in the XRD pattern of Cu 2 ONPs-C@G, as illustrated in Figure 2C, indicating the presence of a combination of nanoparticles and graphite in this sample. The broadness of the peaks can be used to obtain the crystallite size of Cu 2 ONPs by using the Debye-Scherrer formula in Equation (1) as follows [46]:

Characterization of Cu 2 ONPs-C and Cu 2 ONPs-C@G Composite
where D is the average crystallite size of the prepared particles, λ is the X-ray wavelength, typically 0.1541 nm, β is the full width at half maximum (FWHM) in radians, K is the Scherrer's coefficients (0.89) and θ is Bragg's angle. The produced Cu 2 ONPs have an average calculated crystallite size of 14.52 ± 1.08 nm. SEM was used to identify the morphology of bare graphite grains, pure synthesized Cu 2 ONPs-C and graphite grains studded with Cu 2 ONPs-C@G ( Figure 3). As can be seen in Figure 3A, bare graphite grains appeared as overlapped peels of different sizes with rugged and heterogeneous surfaces. At the same magnification, pure synthesized Cu 2 ONPs-C showed nearly monodisperse cubes ( Figure 3B). In the same examined area, it can be seen that some of the prepared Cu 2 ONPs-C can be aggregated to form groups of cubes or to produce a larger cube, as reported before [8]. Figure 3C shows the SEM of Cu 2 ONPs-C@G where graphite peels are apparently studded with brighter cubes of Cu 2 ONPs-C, which are distributed everywhere in the examined area of graphite grains. A few agglomerates in specific locations on the surface of the grains are clearly seen. This is a normal phenomenon, since the graphite grains do not have a smooth surface texture or perfect shape. Additionally, this demonstrates that Cu 2 ONPs-C is in intimate contact with the graphite surface. Cu 2 ONPs-C showed a diameter of more than 200 nm with about 367 nm as the average diameter ( Figure 3D). Chemosensors 2021, 9, x FOR PEER REVIEW 7 of 22      Figure 4 shows the representative TEM image for Cu 2 ONPs-C@G. As can be seen in Figure 4A, two overlapped microcubes (dark black particles; the diameter of each cube is more than 400 nm) with rough surfaces constructed from nanoparticles (light small black particles). Taking into account the resulting SEM analysis ( Figure 3B), and since we used an ordinary mixing method (without using special tools or chemicals), we can assume that there is no significant effect on the morphology or chemical composition of the produced monodisperse cuprous oxide microcubes (Cu 2 ONPs-C). From the clear surface topography of the right TEM image in Figure 4A, we can conclude that the surface nanoparticles are Cu 2 ONPs, which self-assembled in an ordered way to build Cu 2 O microcubes (Cu 2 ONPs-C). Based on measuring the longest diameter for Cu 2 ONPs at the cube surfaces, the diameter range was found to be mainly distributed from 4 to 20 nm with an average diameter of about 9 nm ( Figure 4B). Due to the high magnification and diameter of the cubes, at the downside of the left TEM image ( Figure 4A), we can clearly see part of the graphite layers adherent with the two overlapped cubes.
Chemosensors 2021, 9, x FOR PEER REVIEW 8 of 22 Figure 4 shows the representative TEM image for Cu2ONPs-C@G. As can be seen in Figure 4A, two overlapped microcubes (dark black particles; the diameter of each cube is more than 400 nm) with rough surfaces constructed from nanoparticles (light small black particles). Taking into account the resulting SEM analysis ( Figure 3B), and since we used an ordinary mixing method (without using special tools or chemicals), we can assume that there is no significant effect on the morphology or chemical composition of the produced monodisperse cuprous oxide microcubes (Cu2ONPs-C). From the clear surface topography of the right TEM image in Figure 4A, we can conclude that the surface nanoparticles are Cu2ONPs, which self-assembled in an ordered way to build Cu2O microcubes (Cu2ONPs-C). Based on measuring the longest diameter for Cu2ONPs at the cube surfaces, the diameter range was found to be mainly distributed from 4 to 20 nm with an average diameter of about 9 nm ( Figure 4B). Due to the high magnification and diameter of the cubes, at the downside of the left TEM image ( Figure 4A), we can clearly see part of the graphite layers adherent with the two overlapped cubes.  Since there are different articles reporting the electroanalytical applications of copper nanoparticles, we tried to prove the existence of Cu 2 ONPs localized on graphite grain surfaces using CV. In blank phosphate buffer (pH ≈ 8), a typical cyclic voltammogram was recorded using CPE (dashed line) or Cu 2 ONPs-C@G/PE (solid line) ( Figure 5). CPE showed no peaks in the potential range from −0.6 to +1.0 V. At the same conditions, Cu 2 ONPs-C@G/PE showed a cyclic voltammogram with one well-defined oxidation peak at +0.120 V, which could refer to the oxidation of Cu + to Cu 2+ [47,48]. This is in accordance with XRD measurements which prove the existence of only pure crystalline Cu 2 ONPs and the absence of Cu and CuO nanoparticles. In the reverse direction, a broad and very weak peak was recorded starting from +0.392 V to −0.114 V, which could refer to the reduction in Cu 2+ to Cu + . The appearance of such broad peak can be explained as Cu 2+ was not fully reduced to Cu + [49,50]. The 2nd cathodic peak, which was recorded at E p = −0.275 V, was very distinctive with a larger peak current. This result can be explained as this peak expresses the existence of dual reduction processes: the remaining amount of Cu 2+ and the whole amount of Cu + are reduced to Cu 0 [49,50]. Furthermore, we can assume that the 2nd reduction process (Cu + to Cu 0 ) is the main contributor in the recorded current for this peak thanks to the high concentration of Cu 2 ONPs, which build the cubes (Cu 2 ONPs-C), distributed widely on graphite grain surfaces, which in turn were used to prepare the proposed Cu 2 ONPs-C@G/PE. Chemosensors 2021, 9, x FOR PEER REVIEW Since there are different articles reporting the electroanalytical applications per nanoparticles, we tried to prove the existence of Cu2ONPs localized on graphi surfaces using CV. In blank phosphate buffer (pH ≈ 8), a typical cyclic voltamm was recorded using CPE (dashed line) or Cu2ONPs-C@G/PE (solid line) (Figure showed no peaks in the potential range from −0.6 to +1.0 V. At the same con Cu2ONPs-C@G/PE showed a cyclic voltammogram with one well-defined o peak at +0.120 V, which could refer to the oxidation of Cu + to Cu 2+ [47,48]. This cordance with XRD measurements which prove the existence of only pure cry Cu2ONPs and the absence of Cu and CuO nanoparticles. In the reverse direction, and very weak peak was recorded starting from +0.392 V to −0.114 V, which could the reduction in Cu 2+ to Cu + . The appearance of such broad peak can be explained was not fully reduced to Cu + [49,50]. The 2nd cathodic peak, which was recorded −0.275 V, was very distinctive with a larger peak current. This result can be expla this peak expresses the existence of dual reduction processes: the remaining am Cu 2+ and the whole amount of Cu + are reduced to Cu 0 [49,50]. Furthermore, we sume that the 2nd reduction process (Cu + to Cu 0 ) is the main contributor in the r current for this peak thanks to the high concentration of Cu2ONPs, which build th (Cu2ONPs-C), distributed widely on graphite grain surfaces, which in turn were prepare the proposed Cu2ONPs-C@G/PE.

Electrochemical Behavior of MEB on the Cu2ONPs-C@G/PE
The electrochemical behavior of MEB was investigated at Cu2ONPs-C@G/P CV. For a comparison, Figure 6A shows the recorded cyclic voltammograms in sence and presence of MEB in phosphate buffer solutions (pH ≈ 8) containing 0.1 using CPE and Cu2ONPs-C@G/PE as working electrodes. Using CPE, blank solut not show any peaks in the oxidation and reduction direction, as shown in cycle a the addition of 2 × 10 −5 M MEB, only one clear irreversible anodic peak at Ep = + was recorded, which can be attributed to the oxidation of the amino group (cycl

Electrochemical Behavior of MEB on the Cu 2 ONPs-C@G/PE
The electrochemical behavior of MEB was investigated at Cu 2 ONPs-C@G/PE using CV. For a comparison, Figure 6A shows the recorded cyclic voltammograms in the absence and presence of MEB in phosphate buffer solutions (pH ≈ 8) containing 0.1 M KCl using CPE and Cu 2 ONPs-C@G/PE as working electrodes. Using CPE, blank solutions did not show any peaks in the oxidation and reduction direction, as shown in cycle a1. After the addition of 2 × 10 −5 M MEB, only one clear irreversible anodic peak at E p = +0.800 V was recorded, which can be attributed to the oxidation of the amino group (cycle a2). In the case of Cu 2 ONPs-C@G/PE, the recorded cyclic voltammogram in blank solution, which is shown in cycle b1, was previously explained in Figure 5. After the addition of 2 × 10 −5 M MEB (cycle b2), the Cu + oxidation peak was increased, while no other peaks were observed in the same direction, which can be explained as a formation of an electroactive complex between MEB and the produced Cu 2+ . Additionally, the two recorded reduction peaks shown in cycle b1 completely disappeared (cycle b2). This can be explained as all the produced Cu 2+ was consumed in the formed MEB-Cu complex, which has an irreversible electrochemical behavior. The net current of the added MEB, which is proportional to MEB concentration, can be calculated by Equation (2) as follows: where I net is the calculated current, which is related to the added concentration of MEB; I t is the total recorded current of the formed MEB-Cu complex; and I i is the initial current, which corresponds to oxidation of Cu + to Cu 2+ in the absence of MEB. The proposed MEB-Cu complex structure is elucidated in Scheme 2. Due to Cu + oxidation under the electrochemical conditions, it converts to the preferred Cu 2+ , which may coordinate with MEB in the solution in two possible pathways: the 1st pathway is through the formation of a five-membered ring with the two donor methoxy groups at the edge of MEB structure (complex I), while the 2nd pathway is through the formation a seven-membered ring with the two potential donor atoms (N and O) in the middle of the MEB structure (complex II). Considering the highest stability of the five-membered ring geometrical shape and the steric hindrance at the nitrogen atom, we can confirm that complex I produced by the 1st pathway is the sole product of the suggested reaction.
So, we can assume that our suggested procedure can be applied for the indirect quantification of MEB. Well-dispersed Cu 2 ONPs-C over graphite grain surfaces was beneficial for the electrocatalytic reaction that may occur at the electrode surface [18]. The markable enhancement of Cu + oxidation current indicates the electrocatalytic activity of the Cu 2 ONPs-C@G, which can be attributed to synergistic electrochemical activity between Cu 2 ONPs-C and graphite grain surfaces, which in turn can be illustrated as the existence of the Cu + /Cu 2+ redox couple, making the electrode surface highly sensitive.
A thorough study of the electrochemical behavior of MEB was accomplished by examining the effect of scan rate on the recorded oxidation current of the formed MEB-Cu complex. Cyclic voltammograms of the formed MEB-Cu complex at different scan rate values are shown in Figure 6B (b-h); we can see that with increasing scan rate values, the recorded peak current increased gradually. A linear relationship between peak current (I net ) and square root of scan rate (ν 1/2 ) was obtained, which indicates that the electrochemical oxidation of the formed MEB-Cu complex is a diffusion-controlled process ( Figure 6C). This relationship can be expressed by Equation (3) as follows: Moreover, Figure 6D (a) shows a linear relationship between the logarithm of peak current (log I net ) and logarithm of scan rate (log ν), which can be expressed by Equation (4) as follows: logI net (µA) = 0.4856 logν (mV/s) + 1.5605 R 2 = 0.9995 (4) versible electrochemical behavior. The net current of the added MEB, which is proportional to MEB concentration, can be calculated by Equation (2) as follows: where Inet is the calculated current, which is related to the added concentration of MEB; It is the total recorded current of the formed MEB-Cu complex; and Ii is the initial current, which corresponds to oxidation of Cu + to Cu 2+ in the absence of MEB. The proposed MEB-Cu complex structure is elucidated in Scheme 2. Due to Cu + oxidation under the electrochemical conditions, it converts to the preferred Cu 2+ , which may coordinate with MEB in the solution in two possible pathways: the 1st pathway is showed that the Cu + oxidation peak increased after MEB addition, but with the repeating cyclic sweeps, the recording peak current decreased gradually (cycles b and c). This could be due to the fact that the formed MEB-Cu complex is adsorbed on the Cu2ONPs-C@G/PE surface, causing blocking on its surface and hence a reduction in the exposed effective sites. So, as we indicated before, the Cu2ONPs-C@G/PE surface was renewed after each measurement. Scheme 2. Suggested scheme for the formed MEB-Cu complex.

Computational Analysis Data
To confirm our suggestion that MEB can form the expected complex with produced Cu 2+ , quantum chemical calculations were performed with DFT using the B3LYP function. The resulting quantum chemical parameters are tabulated in Table 1. Additionally, the optimized geometric structures and Mulliken charges of MEB were studied and are presented in Figure 7A,B. The slope of 0.4856 confirms that the occurred electron transfer was a diffusioncontrolled process, since it is very close to the same theoretically expected value of 0.5 for a pure diffusion-controlled process [51]. In the same voltammogram, we can see that a positive shifting occurred in the peak potential with increasing scan rate values, which confirms the irreversibility nature of the recorded oxidation process [51]. Figure 6D (b) shows a linear relationship between peak potential (E p ) and scan rate logarithm (log ν), which can be expressed by Equation (5) as follows: The resulting relationship proves that the oxidation of the MEB-Cu complex at Cu 2 ONPs-C@G/PE was mainly a diffusion-adsorption-controlled process. This behavior was verified by recording successive cyclic voltametric sweeps in the presence of MEB at Cu 2 ONPs-C@G/PE ( Figure 6E). The dashed cycle represents the recorded cyclic voltammogram by Cu 2 ONPs-C@G/PE in blank solution. Cycle a, as explained before, showed that the Cu + oxidation peak increased after MEB addition, but with the repeating cyclic sweeps, the recording peak current decreased gradually (cycles b and c). This could be due to the fact that the formed MEB-Cu complex is adsorbed on the Cu 2 ONPs-C@G/PE surface, causing blocking on its surface and hence a reduction in the exposed effective sites. So, as we indicated before, the Cu 2 ONPs-C@G/PE surface was renewed after each measurement.

Computational Analysis Data
To confirm our suggestion that MEB can form the expected complex with produced Cu 2+ , quantum chemical calculations were performed with DFT using the B3LYP function. The resulting quantum chemical parameters are tabulated in Table 1. Additionally, the optimized geometric structures and Mulliken charges of MEB were studied and are presented in Figure 7A,B.  The ability of MEB molecules to connect with metal orbital may be attributed to their frontier molecular orbital (FMO), i.e., HOMO and LUMO [52]. HOMO is related to the capability of a molecule to donate electrons, whereas LUMO indicates its capability to accept electrons. From Figure 7C,D, it can be noted that the electronic density in HOMO and LUMO on the entire area of each molecule is quite uniform, which is due to its πelectron cloud density. Analysis of HOMO and LUMO is shown in Figure 7C,D, the preferential area of nucleophilic attack in the MEB molecule is the two methoxy groups (O-CH3). Theoretically, the methyl group is a donating electron group that increases the electron density on the oxygen atom, which in turn is considered as a reason to enhance MEB to form a complex with Cu 2+ . The low value of the energy gap (ΔE = ELUMO − EHOMO) indicates the good activity of the studied ligand to react with other species, because the required energy to remove an electron from the last occupied orbital is low [53]. Additionally, the literature reveals that a larger value of the ΔE signifies low reactivity to a chemical species, because the ΔE is related to the softness or hardness of a molecule. A soft molecule is more reactive than a hard molecule because a hard molecule has a larger energy gap [54,55]. According to this, MEB shows a high reactivity to react with Cu, which is evidenced from the results tabulated in Table 1, where ΔE, absolute softness and hardness of MEB were calculated to be 0.1203, 16.6251 and 0.0602 ev, respectively.
Dipole moment is the measure of bond polarity, and it is related to the distribution of electrons in a molecule [56]. As the dipole moment has a high value, the ligand has a high efficiency to react with other species, as predicted in MEB to react with Cu. Furthermore, the total energy (TE) of a molecule includes all forms of potential and kinetic energy, and it is used to determine the molecule's reactivity and stability. The high sta- The ability of MEB molecules to connect with metal orbital may be attributed to their frontier molecular orbital (FMO), i.e., HOMO and LUMO [52]. HOMO is related to the capability of a molecule to donate electrons, whereas LUMO indicates its capability to accept electrons. From Figure 7C,D, it can be noted that the electronic density in HOMO and LUMO on the entire area of each molecule is quite uniform, which is due to its π-electron cloud density. Analysis of HOMO and LUMO is shown in Figure 7C,D, the preferential area of nucleophilic attack in the MEB molecule is the two methoxy groups (O-CH 3 ). Theoretically, the methyl group is a donating electron group that increases the electron density on the oxygen atom, which in turn is considered as a reason to enhance MEB to form a complex with Cu 2+ .
The low value of the energy gap (∆E = E LUMO − E HOMO ) indicates the good activity of the studied ligand to react with other species, because the required energy to remove an electron from the last occupied orbital is low [53]. Additionally, the literature reveals that a larger value of the ∆E signifies low reactivity to a chemical species, because the ∆E is related to the softness or hardness of a molecule. A soft molecule is more reactive than a hard molecule because a hard molecule has a larger energy gap [54,55]. According to this, MEB shows a high reactivity to react with Cu, which is evidenced from the results tabulated in Table 1, where ∆E, absolute softness and hardness of MEB were calculated to be 0.1203, 16.6251 and 0.0602 ev, respectively.
Dipole moment is the measure of bond polarity, and it is related to the distribution of electrons in a molecule [56]. As the dipole moment has a high value, the ligand has a high efficiency to react with other species, as predicted in MEB to react with Cu. Furthermore, the total energy (TE) of a molecule includes all forms of potential and kinetic energy, and it is used to determine the molecule's reactivity and stability. The high stability of a molecule and hence its low ability of electron donation can be confirmed by the highest TE values. From the elaborated results in Table 1, it is evident that MEB has a low TE value (-1396.9 a.u.), which indicates the high tendency of MEB to react with Cu by donating electrons.
Mulliken atomic charge represents the charge of each atom in MEB, and the obtained Mulliken atomic charges are shown in Figure S1A (see Supplementary Material). The calculated atomic charge values of MEB show that the oxygen atom of the methoxy groups (O-CH 3 ) has a higher charge value than the N atom present in MEB. This result indicates that MEB has the ability to form bonds with Cu through O-CH 3 groups, which is emphasized by Molecular Electrostatic Potential (MEP).
Molecular Electrostatic Potential (MEP) was calculated for the suggested complex structures ( Figure S1B; see Supplementary data file). It was found that the negative charge density localized on oxygen atoms (atoms No. 2 and 9; 1st pathway) is higher than the negative charge density localized on the nitrogen atom (atom No. 18; 2nd pathway). High negative charge density was shown as a red color area around oxygen atoms No. 2 and 9, which attached to two methyl groups. This result indicates that a higher ability of MEB to form a metal complex through these two oxygen atoms than only one nitrogen atom exists in MEB. This interpretation confirms the suggested MEB-Cu complex structure expressed as complex (I) in Scheme 2.
Similar quantum calculations were performed for the suggested structure of the MEB-Cu complex to ascertain its stability. The obtained results are tabulated in Table S1 (see  Supplementary data file), and the optimized geometric structures and Mulliken charges of the MEB-Cu complex were studied and are presented in Figure S2 (see Supplementary data file). The illustrated results confirm the high stability of the proposed MEB-Cu complex (I) (Scheme 2), which is in accordance with the fact that the stability of the 5-membered ring (Complex (I)) is higher than the 7-membered ring (Complex (II)).

Optimization of Instrumental Conditions
In the present work, SWAdASV was selected for MEB voltametric investigation because of its high sensitivity [51,57]. The variation of peak current and shape with instrumental parameters, such as accumulation potential, frequency and pulse height, was studied using 2 × 10 −5 M MEB in Britton-Robinson buffer (pH ≈ 8.0) and using Cu 2 ONPs-C@G/PE as the working electrode ( Figure S3A (see Supplementary data file)). In the case of accumulation potential, there was no significant change in peak potential with changing accumulation potential from −1.0 to −0.2 V. It was observed that peak current increased from −1.0 to −0.6 V and then decreased in the range from −0.5 to −0.2 V. The peak current reached its maximum value at accumulation potential −0.6 V, which was used in the subsequent examinations.
The importance of studying frequency and pulse height arises in the current article because we decided to use SWAdASV as the selected technique for the MEB quantification.
It is well known that higher values of frequency and pulse height can increase SWAdASV sensitivity. Additionally, at their high values, both of them have a great effect on the generated capacitive current, which in turn causes distortion in the SWAdASV curves and lowers its sensitivity [57,58]. The effect of frequency on the peak current intensities of 2 × 10 −5 M MEB was studied. The recorded peak current increased with increasing frequency from 10 to 80 Hz. At the highest frequency values, up to 120 Hz, the recorded peak current values became approximately leveled off. So, 100 Hz was selected as the best value for more investigations ( Figure S3B (see Supplementary data file)). Additionally, the effect of pulse height on the same MEB concentration was investigated by varying the pulse height value from 0.01 to 0.30 V. Results illustrated in Figure S3C (see Supplementary data file) reveal that the peak current increased with increasing pulse height from 0.01 to 0.15 V, then became approximately stable with increasing pulse height values up to 0.30 V. So, a pulse height of 0.20 V was chosen for subsequent experiments.

Optimization of Experimental Parameters
Several experimental parameters, such as pH, buffer type, supporting electrolyte type and concentrations, may affect the peak current and its shape. So, the study of the effects of such parameters is very important.
Using Cu 2 ONPs-C@G/PE as the working electrode, the peak current of 2 × 10 −5 M MEB was recorded by SWAdASV in Britton-Robinson buffer at different pH values. This study was carried out over the pH range of 2.25-11.5, as shown in Figure S4A (see Supplementary data file). There was no significant change that occurred in the peak potential of MEB with increasing the pH value. A small peak current was observed at pH ≈ 2.25, which increased gradually up to pH ≈ 8 and then decreased at higher values. Thus, pH ≈ 8 was selected in all consecutive measurements.
Different buffer types of pH ≈ 8 were investigated to select the most suitable buffer type for the SWAdASV quantification of MEB. This study was accomplished in sodium borate-HCl buffer, phosphate buffer, Britton-Robinson and sodium borate-boric acid buffer. Both peak shape and height were taken into account when choosing the buffer type. Figure S4B (see Supplementary data file) shows that the best buffer selected as a medium is phosphate buffer. The effect of phosphate buffer ionic strength on the recorded voltametric peak revealed that minimal background current, the best curve and the highest current were obtained in 0.1 M phosphate buffer.
Among three different supporting electrolytes (KCl, NaNO 3 and NaClO 4 ), the best response of 2 × 10 −5 M MEB in phosphate buffer (pH ≈ 8) was obtained using 0.1 M KCl as the supporting electrolyte ( Figure S4C (see Supplementary data file)). The effect of supporting electrolyte ionic strength was studied at the same conditions, considering minimal background current, best curve shape and the highest sensitivity. In 0.1 M phosphate buffer (pH ≈ 8), the peak current of 2 × 10 −5 M MEB was recorded in the presence of different KCl concentrations. The obtained data show that the best choice, which gives the best shape and highest peak current, was 0.1 M KCl ( Figure S4D, (see Supplementary data file)).

Effect of Accumulation Time and Reproducibility
Accumulation time significantly affects MEB oxidation peak current. This can be obvious through studying the time effect on 6.01 × 10 −9 M and 6.01 × 10 −10 M MEB at the Cu 2 ONPs-C@G/PE working electrode. It was found that with increasing accumulation time up to 240 s, the peak current increased greatly, which induced the existence of MEB-Cu complex adsorption on the Cu 2 ONPs-C@G/PE surface. Upon increasing time values up to 600 s, the recorded peak current increased slowly and approximately leveled off. However, 60 s was generally chosen as the accumulation time to build up the MEB calibration curve using Cu 2 ONPs-C@G/PE. Repeating three experiments on the studied MEB concentrations at 60 s accumulation time to check the reproducibility of the proposed procedure, the relative standard deviation was calculated to be 4%.

Analytical Characterization and Validation of the Proposed Method
In the current article, SWAdASV was selected for the quantification of MEB. The obtained SWAdASV for the successive additions of MEB and its corresponding calibration graph are shown in Figure 8. This voltammogram was conducted using Cu 2 ONPs-C@G/PE as the working electrode in 0.1 M phosphate buffer solutions (pH ≈ 8) containing 0.1 M KCl. As previously explained in the electrochemical behavior of MEB on the Cu 2 ONPs-C@G/PE, the current value corresponding to the added concentration of MEB (I net ) is equal to the difference between the current value resulting from Cu + oxidation (I i ) and the current resulting from the oxidation of the formed MEB-Cu complex (I t ). The obtained results were applied to draw the calibration graph shown in Figure 8 (inset), which shows a linear dependency between the recorded current response and MEB concentration in the range of 5.00 × 10 −11 M-1.10 × 10 −9 M and can be expressed by Equation (6)    Relative standard deviation (% RSD) and other statistical parameters of the proposed method were calculated and are tabulated in Table 2. The low value of % RSD indicates the good precision and accuracy of the method. The lower limit of detection (LOD) and quantification (LOQ) were calculated using the following two equations: LOD = 3σ/b and LOQ = 10σ/b, where σ is the standard deviation of the intercept and b is the slope of the calibration curve [13,14,34,[40][41][42]. The calculated values of LOQ and LOQ for ten repeated quantifications of MEB and related regression data of the elaborated calibration plot are presented in Table 2. Validation of the proposed procedure for SWAdASV quantification of MEB was achieved through investigation of the repeatability of the peak current for three different MEB concentrations (1.00 × 10 −9 , 5.00 × 10 −9 and 6.01 × 10 −8 M). MEB repetitive quantification was recorded on the same day (an intraday) and after five days (an interday). A Cu 2 ONPs-C@G/PE cleaning step was performed for each measurement by wiping the electrode surface against white paper until smooth, and a shiny surface was obtained followed by washing the electrode body with twice distilled water. The standard addition method for five repetitive measurements was applied under the optimized instrumental and experimental conditions; the recovery and % RSD proposed were calculated [59]. Results elaborated in Table 2 confirm that the MEB voltametric quantification using our proposed procedure shows good repeatability and reproducibility within its decided optimized conditions.

Comparison with Other Analytical Methods
The current analytical was more sensitive than other electrochemical methods reported for MEB quantification. Some of these methods are classical electrochemical methods (conductometric [27,28] and potentiometric [29][30][31] methods), which did not offer a sufficient sensitive quantitation limit for pharmacokinetic studies and therapeutic drug monitoring. Additionally, Cu 2 ONPs-C@G/PE was successfully applied for the voltametric quantification of MEB at a lower detection limit than other voltametric sensors previously used for the same purpose [32,33]. Table 3 summarizes a comparison of our results with other voltametric methods previously applied for the quantification of MEB. The tabulated data prove that Cu 2 ONPs-C@G/PE has a sensitivity greater than that of other sensors used for the same purpose. Moreover, we can confirm the great benefits of Cu 2 ONPs-C@G/PE over other analytical methods used for MEB quantification. Our opinion is attributed to an economic point of view due to the low prices of the materials used in Cu 2 ONPs-C@G/PE preparation or its ease of preparation. In contrast, MEB was quantified previously using analytical methods based on expensive nanomaterials or using cumbersome separation methods. In fact, we consider Cu 2 ONPs-C@G as a promising material, since it was successfully applied in the preparation of the current electrochemical sensor (Cu 2 ONPs-C@G/PE) and can be used in the future to prepare various forms of electrochemical sensors with different voltametric applications.

Interference Studies
To demonstrate the selectivity of the proposed method under the optimized conditions, the suggested procedure was performed to record the peak current of MEB in the MEB-Cu complex either in the absence or in presence of different interferences. The studied interferences include some excipients that are used in pharmaceutical industries (e.g., glucose, sucrose, lactose, fructose, starch and talc) and some inorganic salts (e.g., Na + , K + , Zn 2+ , Al 3+ , Ca 2+ , Sr 2+ , Cd 2+ , Co 2+ , Ni 2+ and Mg 2+ ). The experimental results show that up to a thousand-fold excess molar concentration of the studied interferences has no significant effect on the recorded peak current. This is clear from the obtained % recovery with low standard deviation values for five determinations (recovery (%) ± SD = 101.04 ± 1.09 for excipients and 99.85 ± 2.12 for inorganic salts).
To confirm the selectivity of Cu 2 ONPs-C@G/PE toward MEB voltametric determination, a quantum theoretical investigation was achieved aiming to prove the highest ability of MEB to form a complex with Cu 2+ in the presence of the studied excipients. The performed DFT theoretical calculations showed that a ∆E value of MEB is lower than glucose, sucrose and fructose (as an example of studied excipients) (Table S2; see Supplementary  data file). Additionally, Figure S5 (see Supplementary data file) shows the optimized geometric structures of glucose, fructose and sucrose. Another indication for MEB having the highest chemical reactivity in comparison to the studied excipients is that MEB has N and O as heterocyclic atoms that give MEB more negativity than the studied excipients that have only O atoms.

Application of the Cu 2 ONPs-C@G/PE Sensor to Dosage Forms and Biological Fluids
The proposed analytical method was successfully applied for the quantification of MEB in dosage forms and biological fluids (spiked human urine and plasma samples).
Thanks to its sensitivity and selectivity, Cu 2 ONPs-C@G/PE exhibits an excellent analytical behavior towards the voltametric quantification of MEB in two sample dosage forms (Duspatalin and Colospasmin ® forte tablets). The achieved recovery values for studied dosage forms spiked with various MEB concentrations using the standard addition method were in the range of 97.50-103.33%, as shown in Table 4. The achieved results were compared statistically with previous results obtained by SWAdASV [32] and differential pulse voltametric (DPV) [33] reference methods by applying a Student's t-test for accuracy and an F-test for precision at 95% confidence level. Values from the Student's t-test and the variance ratio F-test (1.53 and 1.79, respectively) did not exceed their theoretical values at p = 0.05 (2.75 and 7.84, respectively), indicating no significant differences between the proposed and reference methods. Another application of Cu 2 ONPs-C@G/PE is to quantize MEB in biological fluids (spiked human urine and plasma samples). At defined optimum conditions and using the standard addition method, the recovery percentage values were found in the range of 96.10-104.25%, as shown in Table 4. The t-test and F-test were applied to the obtained dataset, at the 95% confidence level, aiming to prove the reliability between the acquired data using Cu 2 ONPs-C@G/PE and the reference voltametric methods [32,33]. Values obtained from the Student's t-test and the variance ratio F-test (1.15 and 1.59, respectively) did not exceed their theoretical values at p = 0.05 (2.05 and 4.14, respectively). The obtained data elaborated that Cu 2 ONPs-C@G/PE delivered accurate and precise data and showed no significant differences between the proposed and reference methods.

Conclusions
In the current work, we present a facile chemical reduction method followed by intensive mixing to produce a new material of graphite grains studded with cuprous oxide microcubes (Cu 2 ONPs-C@G), which was characterized and verified using various methods. The produced material was used to fabricate a paste electrode (Cu 2 ONPs-C@G/PE), which was successfully applied as a sensitive and selective voltametric quantification of MEB. The electrochemical behavior of MEB was studied using CV, and the optimum analytical parameters were achieved using square wave adsorptive anodic stripping voltammetry (SWAdASV). MEB was successfully quantified over the concentration range of 5.00 × 10 −11 −1.10 × 10 −9 M with LOD and LOQ values of 2.41 × 10 −11 M and 8.05 × 10 −11 M, respectively. Cu 2 ONPs-C@G/PE exhibited a markable analytical sensitivity towards the analytical quantification of MEB in dosage forms (tablets) and biological fluids (spiked human urine and plasma samples). Owing to the simplicity, sensitivity and selectivity of the proposed method, it can be applied easily for the quantification of MEB in various research aspects, such as quality control centers and pharmaceutical research authorities.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.