Research Article | Open Access
Vesna Antunović, Slavna Tešanović, Danica Perušković, Nikola Stevanović, Rada Baošić, Snežana Mandić, Aleksandar Lolić, "Development of a Flow Injection System for Differential Pulse Amperometry and Its Application for Diazepam Determination", Journal of Analytical Methods in Chemistry, vol. 2018, Article ID 6121489, 6 pages, 2018. https://doi.org/10.1155/2018/6121489
Development of a Flow Injection System for Differential Pulse Amperometry and Its Application for Diazepam Determination
This work presents the development of a flow injection system for differential pulse amperometry (DPA) for diazepam determination in the presence of oxygen. The thin flow cell consisted of the bare glassy carbon electrode, reference silver/silver chloride, and stainless steel as the auxiliary electrode. Electrochemical reduction of diazepam (DZP) was characterised by cyclic voltammetry. Azomethine reduction peak was used for DZP quantification. The detector response was linear in the range 20–250 µmol/dm3 of diazepam, with a calculated detection limit of 3.83 µg/cm3. Intraday and interday precision were 1.53 and 10.8%, respectively. The method was applied on three beverage samples, energetic drink, and two different beer samples, and obtained recoveries were from 93.65 up to 104.96%. The throughoutput of the method was up to 90 analyses per hour.
Amperometric detection can operate in three different modes: constant potential, pulse mode, and differential pulse (or multi pulse mode). Detection at constant potential is usually applied with flow injection analysis systems because simple and cheap potentiostats can be used. When an analyte passes through the detector, it generates a current proportional to the analyte concentration .
Differential pulse amperometry (DPA) is a technique described and first applied by Marcenac and Gonon in 1985 . In DPA, a clean potential is applied for electrode cleaning without current sampling. Then two potential pulses are applied after the cleaning step, and the current at the end of each pulse is recorded as the function of time. During the experiment, only the difference of two current samples is displayed. Obtained dependence of ΔI (the difference of the currents recorded at the end of each potential pulse) with time is named differential pulse amperogram . The main advantage of DPA is that its signal is less affected by interferences than with the other amperometric modes.
Throughout the years, DPA has been used for determination of various substances in a flow injection (FI) system. Quadrupole-pulsed amperometric detection was applied for simultaneous determination of glucose and fructose in a FI system . The first pulse was used for the oxidation of glucose and the second for both analytes. The subtraction of obtained currents was used for quantification of fructose. Three-pulse FI-DPA was used for simultaneous determination of phenolic antioxidants (butylated hydroxyanisole and butylated hydroxytoluene) at a boron-doped diamond electrode . Multiple-pulse amperometric flow injection system was described for simultaneous determination of pharmaceutical active substances, paracetamol and ascorbic acid , caffeine, ibuprofen, and paracetamol , and dypirone and paracetamol in pharmaceutical formulations on a bare glassy carbon electrode (GCE) . A simple, fast and low-cost multiple-pulse amperometric FI system was applied for determination of sildenafil citrate (SC) in Viagra . Three sequential potential pulses were applied as a function of time. SC is detected at two different irreversible oxidation processes (at 1.6 and 1.9 V), and the third pulse (1.0 V) was used for regeneration of the surface of the boron-doped diamond electrode.
Diazepam (DZP) (7-chloro-1,3-dihydro-1-methyl-5-phenyl-2H-1,4-benzodiapin-2-one) belongs to the group of 1,4-benzodiazepines. It is one of the most prescribed 1,4-benzodiazepines for a number of conditions such as anxiety, insomnia, epilepsy, and muscular spasms . Its low price and availability on the black market increase the risk of its abuse. Since it is highly liposolubile, it distributes rapidly inside tissues and has the high absorption rate. The absorption rate is increased when it is mixed with alcohol which enhances its sedative effect . That is why diazepam is well known as “date rape” drug, and it is classified as a drug-facilitated crimes (DFC) drug . A fast and simple system for diazepam determination would be of great importance for forensic science and resolving criminal activities which are related with benzodiazepine drug abuse .
Determination of benzodiazepines is well documented in the literature, and various spectrophotometric  and electrochemical detectors [15–20] were employed. Mass spectrometers were also applied to benzodiazepine determination ; their wide linear ranges and low detection limits are advantageous but are still expensive compared to simple flow injection systems.
Although there are many references describing determination of benzodiazepine by electrochemical detectors, there is no paper of DPA application on diazepam determination.
Electrochemical determination of diazepam is based on reduction of 4,5-azomethine group. This simple reduction yielding dihydro specie is used for quantification of diazepam . Applying more negative potentials increases interferences from oxygen, hydrogen ions, and other metal ions present in the sample . Due to high background currents, sensitive determination is almost impossible. Concentration of hydrogen and metal ions is a matter of optimization process and can be controlled. However, the presence of oxygen is more problematic and can be decreased either by removing with nitrogen or by applying modifications of electrode surfaces. Often purging of nitrogen or argon prior the analyses for a few minutes or even overnight is usually enough to remove dissolved oxygen from working solutions . When flow injection systems are applied, another problem arises: the choice of tubing material. Usually, standard Teflon tubings should be avoided due to absorption of oxygen through the polymer. Other polymer materials are more or less chemically inert and also affect the price of the flow system or even using glass or metallic materials but they are too rigid for handling the simple flow system. Lozano-Chaves  applied modified carbon-paste electrodes as sensors for the determination of diazepam and its metabolites temazepam and oxazepam in biological fluids. As a modifier, they used 5% bentonite and 5% zeolite. Bentonite-modified electrodes showed better sensitivities. Diazepam was determined at pH 10 in Britton–Robinson buffer, with the concentration range of 0.025–3.0 mg/dm3 and the detection limit of 0.021 mg/dm3; all experiments were performed under nitrogen atmosphere. Multiwall carbon nanotube-ionic liquid modified paste electrode was used for diazepam determination in real samples . Analytical parameters for diazepam determination were linearity in the range of 0.02–0.76 mg/dm3 and the detection limit of 4.1 µg/dm3. Experiments were performed under nitrogen atmosphere. As mentioned earlier, reduction of 4,5-azomethine group is often used for quantification of diazepam. Group of authors  noticed that there is an oxidation peak in a reverse scan in cyclic voltammetry and they determined the oxidation potential (+1.0 V versus Ag/AgCl) for diazepam in drinks on a screen-printed electrode by adsorptive stripping voltammetry. The detector response was linear in the range of 7.1–285 mg/dm3 and a detection limit of 1.8 mg/dm3.
This work presents the development of a flow injection for diazepam determination with differential pulse amperometry and its application in spiked beverage samples without removal of dissolved oxygen.
2. Materials and Methods
2.1. Reagents and Chemicals
All reagents used were of analytical grade quality, and all solutions were prepared in degassed and filtered deionised water. Buffer solution of HCl/KCl pH 1 was prepared by mixing appropriate volumes of the diluted solution of the hydrochloric acid (0.2 mol/dm3, Carlo Erba, Val de Reuil, France) and potassium chloride (0.2 mol/dm3, Betahem, Belgrade, Serbia). Britton-Robinson buffer was prepared by mixing 0.04 mol/dm3 CH3COOH (Carlo Erba, Val de Reuil, France), 0.04 mol/dm3 H3PO4 (Carlo Erba, Val de Reuil, France), and 0.04 mol/dm3 H3BO3 (Betahem, Belgrade, Serbia), and the pH value is adjusted with 0.02 mol/dm3 NaOH (Betahem, Belgrade, Serbia). The stock solution of diazepam (Hemofarm, Vršac, Serbia) was prepared by dissolving the required mass in methanol to give a concentration of 10 mmol/dm3. Stock solution was maintained under refrigerated conditions in the absence of light, and it was prepared once a week. Working standards were prepared daily by dilution of this solution with appropriate buffers.
Cyclic voltammetry was used for electrochemical study of diazepam. For the cyclic voltammetry, a CHI 800C potentiostat was used. A three-electrode system consisted of the glassy carbon working electrode (CH Instruments, USA, model CHI104, 3 mm in diameter), Ag/AgCl reference electrode (CH Instruments, USA; CHI111), and platinum wire as auxiliary electrode. pH measurements were performed using WTW720 pH-meter equipped with SenTix 81 pH electrode. Prior the cyclic voltammetry experiments, all solutions were purged with oxygen free nitrogen for 15 minutes.
Cleaning of glassy carbon electrodes was performed mechanically on a polishing pad using aluminium paste of different grain sizes (1, 0.3, and 0.05 µm, Buehler, USA). It was washed with distilled water between each paste, and after the finest paste, it was washed with methanol, distilled water, and then air-dried. This process was repeated each day before the start of the recording and in the case if detector response was not reproductive.
2.3. Flow Injection Analysis System
Peristaltic pump model Mini S 840 (Ismatec, Switzerland) was used. Two-position injection valve model 5020 (Rheodyne, USA) equipped with a sample loop was used. Amperometric flow cell (thin layer, BASi, USA) consists of working, reference, and auxiliary electrodes . The working electrode was glassy carbon electrode (BAS Instruments, model MF-1008, USA) consisting of two circular, dual series electrodes which are embedded in a polymer based on fluorocarbons. The reference Ag/AgCl (BAS Instruments, model MF-2021, USA) electrode is filled with 3 mol/dm3 NaCl solution, and the auxiliary electrode is made of stainless steel. Teflon gasket was placed between the auxiliary and working electrodes; its thickness regulates the volume of the working solution. Electrochemical measurements were performed on model CHI760b. The CHI software was used for data acquisition.
2.4. Electrochemical Procedure
Each DPA cycle applied to the working electrode consisted of three steps: first step, cleaning potential of +0.7 V for 50 ms, second −0.75 V for 50 ms (first potential pulse), and then −0.95 V for 50 ms (second potential pulse). The analytical signal was calculated as the difference between the average signal at last 25 ms of the second and the first potential pulse. DPA measurements were performed on CHI760b potentiostat. For this set of experiments, there was no need to expel the oxygen from solutions.
3. Results and Discussion
3.1. Cyclic Voltammetry of Diazepam and pH Effect
Cyclic voltammograms of a 0.5 mol/dm3 diazepam in HCl/KCl buffer (pH 1) were recorded in the potential range of −0.2 to −0.9 V. Voltammograms were recorded for different scan rates from 20 to 150 mV s−1. In the negative scan, one reduction peak can be observed at around −0.76 V (Figure 1). The value of the peak potential became more cathodic with increasing scan rate and varied linearly with square root of scan rate, indicating irreversibility of the reduction and diffusion-controlled process on the electrode surface.
The effect of pH on reduction of diazepam was investigated when 100 µmol/dm3 diazepam was recorded in HCl/KCl buffer (pH 1) and Britton–Robinson buffer of different pH (3–8) values in the potential range −0.2 to −1.2 V. With increase of pH, the cathodic peak shifts to more negative values indicating the dependence of the reduction potential on the pH value of the medium (Figure 2). As the media becomes less acidic, the FIA signal becomes less stable and reproducible. Hence, the further experiments were performed in pH 1 solutions.
3.2. Optimization of Flow Injection Parameters
The following parameters were investigated for the sensitivity of the thin-layer flow cell: the volume of the sample loop and the gasket thickness. The optimal conditions in FIA system were investigated with 100 µmol/dm3 diazepam solution. Sample loops of 0.050, 0.075, and 0.120 cm3 were tested. The detector response was the most sensitive with 0.050 cm3 of the injected sample. The gasket determines the working volume of the flow cell. All four commercially available gaskets from BASi were tested (13, 51, 127, and 381 µm thickness). For this type of flow cell, the best sensitivity was obtained for 51 µm thick gasket, which was already confirmed in our previous work .
3.3. Optimization of Differential Pulse Amperometry
The potentials used for the two DPA pulses were chosen according to the results obtained by cyclic voltammetry. Hence, the DPA pulses should be a little lower and a little higher than the peak potential of diazepam reduction. We investigated the following values of potential pulses, Eclean 0 to +1 V, E1 −0.6 to −0.75 V, and E2 −0.8 to −0.95 V; each potential was held for 50 ms, and the current was collected at last 25 ms for E2 and E1. All experiments were performed with 0.1 mmol/dm3 diazepam solution in HCl/KCl buffer pH 1. The 0.7 V was applied for cleaning potential due to the signal shape and the system reproducibility. The sensitivity of the detector for different values of E1 and E2 is presented in Table 1, and it shows that the best response was obtained for E1 = −0.75 V and E2 = −0.95 V. Figure 3 presents the DPA waveform applied for diazepam.
Under these conditions, there was no need for purging the solutions prior the experiments or working in the inert atmosphere.
3.4. Analytical Parameters
The calibration curve obtained from the DPA peak current for diazepam reduction at different diazepam concentration shows linear behaviour between 20 and 250 µmol/dm3 with the regression equation ΔI (µA) = 0.063c (µM) + 0.430 (Figure 4). The presented calibration curve is the average obtained from the three curves recorded for the same set of standards; error bars in each point of the curve are standard deviations. The detection limit was calculated to be 13.4 µmol/dm3, but measured as 3 s/m, where s is the standard deviation and m is the slope, and the limit of quantification was calculated to be 44.8 µmol/dm3. Obtained limit of the detection is not comparable to LODs obtained by hyphenated techniques (chromatography/mass spectrometry) which often require complicated sample preparation [13, 27–30]. Electrochemical detection is seldom described for determination of diazepam in beverage samples. Honeychurch et al.  applied adsorptive stripping voltammetry with medium exchange limit, and the limit of detection was 1.8 µg cm−3.
The absolute limit of detection for the sample loop volume of 0.050 cm3 was 0.19 µg of diazepam. Reproducibility of the system was checked by six consecutive injections of 100 µmol/dm3 DZP on the same day for investigation of intraday precision (insert picture on Figure 4), and the relative standard deviation for intraday precision was 1.53%. Interday precision was obtained by recording peak current of 100 µmol/dm3 DZP during three consecutive days; the results showed decrease of precision, presented by the relative standard deviation of 10.8% for interday measurements. Since simple flow injection system was used and the signal frequency was dependent on tubing length, only the throughoutput was up to 90 analyses per hour (Table 2).
3.5. Interference Study
Ascorbic acid, lactose, glucose, and citric acid are the common compounds present in the beverage samples. The effect of interferants on diazepam determination was investigated by injecting solutions containing 100 µmol dm−3 of each interferant and mixtures with diazepam of the same concentration. Results showed that interferants decrease diazepam signal by 4–10% (Table 3).
3.6. Determination of Diazepam in Beverage Samples
The purpose of this method was to investigate concentrations relevant for forensic studies. On Serbian market, diazepam tablets containing 2, 5, or 10 mg of active substance are available; when a tablet is dissolved in 200 cm3 of a beverage, obtained concentrations are from 35–175 µmol/dm3 of diazepam, which fits in the linear range of the detector.
Beverage samples (Guarana and two beer samples) were purchased at a local store. 12.5 cm3 of a sample was diluted with buffer and appropriate volume of diazepam standard solution, ultrasonicated for 30 minutes, and made up with distilled water to 25.00 cm3. The samples were spiked with 50 µmol/dm3 of diazepam. Solutions were injected in triplicate, and the results of recoveries for spiked samples of different concentrations are presented in Table 4.
Flow injection method with differential pulse amperometry was developed and applied for quantification of diazepam in investigated set of beverage samples. Differential pulse amperometry enables the determinations in the presence of oxygen without surface modification. Under the optimal conditions, the linear range of the system was from 20 to 250 µmol/dm3 of diazepam, with the detection limit of 3.83 µg/cm3, the intraday precision of 1.53%, and a throughoutput of 90 analyses per hour. The method was applied on diazepam determination in beverage samples; 0.050 cm3 of sample was enough for fast, precise, and sensitive quantification. With the flow rate of 1 cm3/min, the method is economical, producing small volumes of waste. The simple sample preparation enables injection of samples without extraction or filtration.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was done within the framework of the research Project no. 172051 supported by the Ministry of Education, Science and Technological Development of Serbia.
- N. R. Stradiotto, H. Yamanaka, and M. V. B. Zanoni, “Electrochemical sensors: a powerful tool in analytical chemistry,” Journal of the Brazilian Chemical Society, vol. 14, no. 2, pp. 159–173, 2003.
- F. Marcenac and F. Gonon, “Fast in vivo monitoring of dopamine release in the rat brain with differential pulse amperometry,” Analytical Chemistry, vol. 57, no. 8, pp. 1778-1779, 1985.
- F. J. Heredia-López, J. L. Góngora-Alfaro, F. J. Alvarez-Cervera, and J. L. Bata-Garcı́a, “A PC-controlled voltage pulse generator for electroanalytical applications,” Review of Scientific Instruments, vol. 68, no. 4, pp. 1879–1885, 1997.
- W. Surareungchai, W. Deepunya, and P. Tasakorn, “Quadruple-pulsed amperometric detection for simultaneous flow injection determination of glucose and fructose,” Analytica Chimica Acta, vol. 448, no. 1-2, pp. 215–220, 2001.
- R. A. Medeiros, B. C. Lourenção, R. C. Rocha-Filho, and O. Fatibello-Filho, “Simple flow injection analysis system for simultaneous determination of phenolic antioxidants with multiple pulse amperometric detection at a boron-doped diamond electrode,” Analytical Chemistry, vol. 82, no. 20, pp. 8658–8663, 2010.
- W. T. P. Dos Santos, E. G. N. De Almeida, H. E. A. Ferreira, D. T. Gimenes, and E. M. Richter, “Simultaneous flow injection analysis of paracetamol and ascorbic acid with multiple pulse amperometric detection,” Electroanalysis, vol. 20, no. 17, pp. 1878–1883, 2008.
- S. C. Chaves, P. N. C. Aguiar, L. M. F. C. Torres et al., “Simultaneous determination of caffeine, ibuprofen, and paracetamol by flow-injection analysis with multiple-pulse amperometric detection on boron-doped diamond electrode,” Electroanalysis, vol. 27, no. 12, pp. 2785–2791, 2015.
- W. T. P. dos Santos, D. T. Gimenes, E. G. N. de Almeida, S. P. Eiras, Y. D. T. Albuquerque, and E. M. Richter, “Simple flow injection amperometric system for simultaneous determination of dipyrone and paracetamol in pharmaceutical formulations,” Journal of the Brazilian Chemical Society, vol. 20, no. 7, pp. 1249–1255, 2009.
- A. Carlos, V. Lopes Júnior, R. De Cássia et al., “Determination of sildenafil citrate (Viagra®) in various pharmaceutical formulations by flow injection analysis with multiple pulse amperometric detection,” Journal of the Brazilian Chemical Society, vol. 23, no. 10, pp. 1800–1806, 2012.
- J. Riss, J. Cloyd, J. Gates, and S. Collins, “Benzodiazepines in epilepsy: pharmacology and pharmacokinetics,” Acta Neurologica Scandinavica, vol. 118, no. 2, pp. 69–86, 2008.
- F. Charlson, L. Degenhardt, J. McLaren, W. Hall, and M. Lynskey, “A systematic review of research examining benzodiazepine-related mortality,” Pharmacoepidemiology and Drug Safety, vol. 18, no. 2, pp. 93–103, 2009.
- D. S. M. Ribeiro, J. A. V. Prior, J. L. M. Santos, and J. L. F. C. Lima, “Automated determination of diazepam in spiked alcoholic beverages associated with drug-facilitated crimes,” Analytica Chimica Acta, vol. 668, no. 1, pp. 67–73, 2010.
- M. Acikkol, S. Mercan, and S. Karadayi, “Simultaneous determination of benzodiazepines and ketamine from alcoholic and nonalcoholic beverages by GC-MS in drug facilitated crimes,” Chromatographia, vol. 70, no. 7-8, pp. 1295–1298, 2009.
- A. M. Gil Tejedor, P. Fernández Hernando, and J. S. Durand Alegría, “A rapid fluorimetric screening method for the 1,4-benzodiazepines: determination of their metabolite oxazepam in urine,” Analytica Chimica Acta, vol. 591, no. 1, pp. 112–115, 2007.
- M. E. Lozano-Chaves, J. M. Palacios-Santander, L. M. Cubillana-Aguilera, I. Naranjo-Rodríguez, and J. L. Hidalgo-Hidalgo-de-Cisneros, “Modified carbon-paste electrodes as sensors for the determination of 1,4-benzodiazepines: application to the determination of diazepam and oxazepam in biological fluids,” Sensors and Actuators B: Chemical, vol. 115, no. 2, pp. 575–583, 2006.
- M. M. Correia dos Santos, V. Famila, and M. L. Simões Gonçalves, “Copper–psychoactive drug complexes: a voltammetric approach to complexation by 1,4-benzodiazepines,” Analytical Biochemistry, vol. 303, no. 2, pp. 111–119, 2002.
- G. B. El-Hefnawey, I. S. El-Hallag, E. M. Ghoneim, and M. M. Ghoneim, “Voltammetric behavior and quantification of the sedative-hypnotic drug chlordiazepoxide in bulk form, pharmaceutical formulation and human serum at a mercury electrode,” Journal of Pharmaceutical and Biomedical Analysis, vol. 34, no. 1, pp. 75–86, 2004.
- V. K. Gupta, R. Jain, K. Radhapyari, N. Jadon, and S. Agarwal, “Voltammetric techniques for the assay of pharmaceuticals—a review,” Analytical Biochemistry, vol. 408, no. 2, pp. 179–196, 2011.
- C. G. Amorim, A. N. Araújo, M. C. B. S. M. Montenegro, and V. L. Silva, “Cyclodextrin-based potentiometric sensors for midazolam and diazepam,” Journal of Pharmaceutical and Biomedical Analysis, vol. 48, no. 4, pp. 1064–1069, 2008.
- K. C. Honeychurch, A. Crew, H. Northall et al., “The redox behaviour of diazepam (Valium®) using a disposable screen-printed sensor and its determination in drinks using a novel adsorptive stripping voltammetric assay,” Talanta, vol. 116, pp. 300–307, 2013.
- R. Wang, X. Wang, C. Liang et al., “Direct determination of diazepam and its glucuronide metabolites in human whole blood by μElution solid-phase extraction and liquid chromatography–tandem mass spectrometry,” Forensic Science International, vol. 233, no. 1–3, pp. 304–311, 2013.
- N. Nagappa, T. Mimani, B. Sheshadri, and S. M. G. Mayanna, “Cyclic voltammetric studies of diazepam using glassy carbon electrode-estimation of diazepam in pharmaceutical samples,” Chemical & Pharmaceutical Bulletin, vol. 46, no. 4, pp. 715–717, 1998.
- W. Lund, M. Hannisdal, and T. Greibrokk, “Evaluation of amperometric detectors for high-performance liquid chromatography: analysis of benzodiazepines,” Journal of Chromatography A, vol. 173, no. 2, pp. 249–261, 1979.
- M. A. Zare, M. S. Tehrani, S. W. Husain, and P. A. Azar, “Multiwall carbon nanotube-ionic liquid modified paste electrode as an efficient sensor for the determination of diazepam and oxazepam in real samples,” Electroanalysis, vol. 26, no. 12, pp. 2599–2606, 2014.
- A. Lolić, S. Nikolić, and P. Polić, “Optimization and application of the gas-diffusion flow injection method for the determination of chloride,” Journal of the Serbian Chemical Society, vol. 66, no. 9, pp. 637–646, 2001.
- L. Thomas, J. L. Vilchez, G. Crovetto, and J. Thomas, “Electrochemical reduction of diazepam,” Journal of Chemical Sciences, vol. 98, no. 3, pp. 221–228, 1987.
- M. Gros, S. Rodríguez-Mozaz, and D. Barceló, “Fast and comprehensive multi-residue analysis of a broad range of human and veterinary pharmaceuticals and some of their metabolites in surface and treated waters by ultra-high-performance liquid chromatography coupled to quadrupole-linear ion trap tandem mass spectrometry,” Journal of Chromatography A, vol. 1248, pp. 104–121, 2012.
- K. Soltaninejad, M. Karimi, A. Nateghi, and B. Daraei, “Simultaneous determination of six benzodiazepines in spiked soft drinks by high performance liquid chromatography with ultra violet detection (HPLC-UV),” Iranian Journal of Pharmaceutical Research, vol. 15, no. 2, pp. 457–463, 2016.
- G. Famigliani, V. Termopoli, P. Palma, and A. Cappielo, “Liquid chromatography-electron ionization tandem mass spectrometry with the Direct-EI interface in the fast determination of diazepam and flunitrazepam in alcoholic beverages,” Electrophoresis, vol. 37, no. 7-8, pp. 1048–1054, 2016.
- L. Gautam, S. Sharratt, and M. Cole, “Drug facilitated sexual assault: detection and stability of benzodiazepines in spiked drinks using gas chromatography-mass spectrometry,” PLoS One, vol. 9, no. 2, Article ID e89031, 2014.
Copyright © 2018 Vesna Antunović et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.