1. Introduction
In recent years, the growing demand for devices detecting specific substances has led to looking for new technologies with quick response and detection characteristics. In this field, smartphone-based electrochemical biosensors are emerging as powerful tools for the fast detection of markers in healthcare, environmental monitoring, and food safety [
2]. These biosensors can provide rapid, sensitive, and selective analysis of food matrices with minimal sample preparation and low cost. Moreover, smartphone-based electrochemical biosensors can be integrated with wireless communication and cloud computing to enable real-time data processing, storage, and sharing [
3].
Smartphone-based electrochemical biosensors have been mainly applied to point-of-care (POC) diagnostic tools for health monitoring [
4]. At the same time, still, few works in the literature are focused on food quality evaluation [
5]. Analyzing the characteristics of food products results of crucial importance for the determination of their quality (e.g., pH and temperature standards [
6]), nutritional value (e.g., nutrients level [
7]), authenticity (e.g., adulterants [
8] detection) and safety (e.g. pathogens, allergens, contaminants detection [
9]).
Currently, particular attention is given to promoting body wellness by following optimal dietary protocols specific to the individual, providing the proper micronutrient intake and preventing chronic diseases such as diabetes, obesity and heart diseases [
10,
11]. A key role is played by polyphenols (Phs), non-enzymatic exogenous antioxidants present in natural phytochemical compounds (
e.g. coffee, tea, fruits) able to direct inactivate free radicals, which in the human body are responsible for cellular oxidative stress, involving the damage of membranes, proteins, DNA and lipids [
12]. Phs showed the ability to promote immunity towards foreign pathogens through numerous biological activities [
13,
14], reducing risks of arteriosclerosis, cancer, neurodegenerative diseases, and osteoporosis [
15], as well as regulating intestinal mucosal immune responses and allergic diseases [
16,
17,
18].
Coffee beverage is widely consumed worldwide, representing one of the principal sources of Phs in the human diet [
19]. Several scientific studies directly correlate coffee intake with a lower risk of type 2 diabetes mellitus, kidney stones, Parkinson's disease, liver cancer or fibrosis etc. [
20], making the Phs content an important index of the nutritional quality of coffee at industrial and commercial scale.
Phs concentration in coffee primarily depends on its variety. Hundred coffee species are available and two of them are the most cultivated in the world, "
Arabica" (Coffea Arabica L) and "
Robusta" (Coffea Canephora) [
21]. When the coffee is not roasted, it is called "green", specifically
"Robusta" quality, in this case, shows the maximum polyphenol content (7-14.4 % of the dry weight), which is higher than that for Arabica variety (6-7% of the dry weight) [
22,
23]. The roasting degree of the beans is also an important parameter since both time and T directly affect the Phs content and the airflow manipulation in commercial coffee-roasting setups [
24]. Green coffee has the highest range of Phs [
25]. In contrast, high-temperature treatments induce Phs degradation, sugar decomposition, lipid oxidations and pyrolysis, determining the color, aroma and flavor, and its bioactive compounds content.
Typically measuring the polyphenol content in coffee samples is done in laboratories using expensive and complex instruments. Smartphone-based electrochemical biosensors offer a low-cost, portable and user-friendly alternative for polyphenol detection in various settings [
26]. Herein, to the best of our knowledge, we demonstrate for the first time a low-cost smartphone-based amperometric biosensor for the on-site quantification of total Phs content in coffee blends. The device is based on a previous electrochemical platform developed by our laboratory [
26], where a carbon screen-printed electrode (SPE) was modified via drop-casting with multi-walled carbon nanotubes (CNTs) and gold nanoparticles (GNPs). After that, Laccase enzyme from
Trametes versicolor (
TvLac) was immobilized through PVA-SbQ photopolymerization. In the present work, a commercial SPE already modified at the surface by CNTs-GNPs nanostructures was employed, quickly immobilizing the
TvLac via glutaraldehyde (GA) cross-linking. As well known, this enzyme is able to catalyze the oxidation of ortho- and para-diphenols, aminophenols, polyphenols, polyamines, lignins and aryl diamines [
27]. The developed biosensor is optimized for
TvLac and GA concentrations by monitoring the enzyme kinetic efficiency through electrochemical tests. The commercial electrodes modifications were investigated by scanning electron microscopy (SEM) and Fourier-transform infrared (FT-IR) measurements, highlighting improvements of electroactive area (A
el), surface roughness index (ρ) and also homogeneity as
TvLac is immobilized through the formation of GA cross-linked aggregates. The smartphone-based device proposed provides a sensitivity of 0.12 μA/μM, a linear response ranging between 5 and 70 μM and a limit of detection (LOD) of 2.99 μM towards Phs content in coffee samples. The influence of the coffee variety and the roasting degree, both time of roasting and temperature, are further evaluated. The results obtained are validated through the Folin-Ciocâlteu standard method. These biosensors have shown promising results in terms of accuracy, sensitivity, specificity and reproducibility, and have potential applications in food quality control.
2. Results and Discussions
Laccase from
Trametes versicolor (
TvLac) is a blue multicopper oxidase enzyme (BMCO) of the oxidoreductase family, with an active site consisting of three copper centers, namely T1Cu, T2Cu, T3Cu. The T1Cu site is responsible for the oxidation of the substrate, while T2Cu-T3Cu center binds the oxygen, which is converted into water through a four-electron reduction process [
28]. Thus, we based our biosensor on the activity of
TvLac that can mediate the reduction of O
2 to H
2O, which is possible through the catalytic oxidation of specific substrates such as Phs. It is possible to measure the concentration of Phs from the cathodic current produced by the Phs
(ox) reduction during the enzyme regeneration. The
TvLac has been immobilized through GA intermolecular cross-linking and physical adsorption, evaluating the optimal concentration in order not to obstacle the substrate movement through the membrane and improve the sensor stability over time.
Improvement of biosensor sensitivity and enzyme loading are obtained by testing commercial SPEs already modified with different nanomaterials combinations: i) unmodified carbon SPE (DRP-110), ii) multi-walled carbon nanotubes (CNTs) modified carbon SPE (DRP-110CNT), iii) gold nanoparticles (GNPs) modified carbon SPE (DRP-110GNP), iv) CNTs and GNPs modified carbon SPE (DRP-110CNT-GNP).
The optimal amount of
TvLac to immobilize was assessed onto the DRP-110CNT-GNP SPE by testing three different amounts of the enzyme: 1, 1.25 and 1.5 enzyme units (U). The catalytic current,
Figure 1a, increases with increasing the immobilized enzyme units onto the working electrode (WE) surface. However, despite 1.50 U provide for the highest catalytic current, 1.25 U was chosen as the optimal amout, giving a current magnitude almost comparable to that obtained with 1.50 U, but giving a greater reproducibility.
GA concentration was also evaluated, namely 0.8%, 1.3% and 1.6% v/v. Chronoamperometry tests (CA) were performed applying a constant potential of -0.1 V
vs. silver pseudo-reference electrode (AgPRE). In
Figure 1b is shown that 0.8% v/v GA is responsible for a less stability of the current over time, losing the CA response at 6 μM of catechol. This is probably due to a weaker cross-linking, exposing more the enzyme towards the solution and inducing eventually the bioreceptor loss.
The biosensors modified with GA at 1.3% and 1.6% v/v showed comparable responses towards the substrate however, in the case of GA at 1.6% v/v the stronger cross-linking produced a denser GA matrix, hindering the substrate diffusion and causing the saturation of the signal at a catechol concentration of 20 μM. The 1.3% v/v GA solution was selected as optimal, ensuring both a wider concentration range of analysis (0-2000 μM) and good sensitivity.
The surface modification of SPE electrodes was investigated through FT-IR/imaging analysis (
Figure 2) and scanning electron microscopy (SEM) (
Figure S1). Specifically, FT-IR measurements have been performed by collecting 25 different spots in a 5x5 grid array area on the electrode surface (25 pixels).
Figure 2a reports the FT-IR responses of DRP-110, DRP-110CNT and DRP-110CNT-GNP electrodes, at grid central spot. The DRP-110 electrode shows two characteristic broad bands between 879-1103 cm
-1 and 1261-1458 cm
-1, corresponding to the alkoxy C–O stretching and C–C stretching respectively. After introducing CNTs at the electrode surface (DRP-110CNT) there is a slight increase of the aforementioned bands, while the presence of GNPs (DRP-110CNT-GNP) enhances the absorbance response inducing a surface IR enhancement effect (SEIRS) [
29]. GNPs presence further promote intensity of response and also appearance of two broad bands at 2728-2983 cm
-1 and 3056-3446 cm
-1, associated to the C–H asymmetric and symmetric stretching and alcohol/carboxylic O–H stretching respectively.
The FT-IR profile when
TvLac is immobilized onto the WE surface by physical adsorption shows five characteristic bands,
Figure 2b, namely the C – O – C asymmetric and symmetric stretching at 1014 cm
-1 and 1002 cm
-1, the C – N stretching at 1373 cm
-1, the CO – NH stretching of peptide linkage at 1648 cm
-1, the alkylic C – H stretching at 2932 cm
-1 and N – H/O – H stretching at 3310 cm
-1 [
30,
31,
32]. FT-IR tests are also performed after depositing
TvLac onto the DRP-110CNT-GNP electrode. The presence of hydrophobic forces normally induces a difficulty in the identification of the C–O–C and CO–NH FT-IR spectrum [
31]. However, in our case the presence of a
TvLac multilayer onto the electrode, where the outer protein layer has no direct interaction with the carbon-based electrode, it is possible to distinguish vibration of the C–O–C and CO–NH functionalities (1002-1140 cm
-1 and 1648 cm
-1 respectively). By integrating the area of the CO–NH stretching band over the entire electrode surface (i.e., considering the 25 different collected spectra), it is possible to map the surface distribution of
TvLac,
Figure 2c,d. The CO–NH stretching band was selected as reference signal for the FT-IR imaging, being the peptide linkage present in the whole
TvLac structure. In fact, the choose of other functionalities present only in specific regions of the protein could only give information about the
TvLac enzyme orientation respect to the electrode surface. Since the CO–NH stretching vibration is directly influenced by the hydrophobic interaction strength between the protein and the carbon surface, the less the amount of
TvLac onto a specific region, the less the possibility to observe the CO–NH stretching band when the FT-IR spectrum is collected in that spot. Therefore, by observing the map an inhomogeneous distribution of the CO–NH stretching vibration intensity highlights not homogeneous deposition of enzyme on the electrode surface, thus the protein doesn’t cover the electrode homogeneously when immobilized in absence of GA. Differently, when
TvLac is immobilized through the formation of GA cross-linked aggregates, intensity of FT-IR absorbance is higher. In this case, the O–H stretching vibration band is influenced by different contributes: i) the vibration of the OH present in the protein structure; ii) the vibration of the OH related to the incorporated water molecules; iii) the vibration of the OH present in the GA structure [
33]. Also, the alkylic C – H stretching at 2932 cm
-1 is due to both GA and
TvLac. The remaining vibrational modes are only related to the protein; therefore, the increase of the spectrum intensity is due to a higher concentration of the biomolecule in the central spot, which suggests a different distribution over the electrode surface respect to the previous described. This is confirmed by the distribution map in
Figure 2d, where the vibration intensity of the CO–NH stretching mode is strong in the central region of the electrode surface, index of a homogeneous immobilization of
TvLac via GA cross-linking.
CV tests,
Figure S2, were performed to calculate the electroactive area (A
el) of the four SPEs. Specifically, CV at different scan rates have been performed in a solution of 1 mM [Fe (CN)6]
3−/4− with KCl 0.1 M. According to the Randles-Ševčík equation and the obtained plot of I
p,a vs v
1/2 [
34] the A
el values have been calculated and summarized in
Table S1 with the relative roughness factors (ρ) defined by the electroactive/geometric area ratio (A
el/A
geo). Both A
el and ρ increase in the order DRP-110 < DRP-110GNP < DRP-110CNT < DRP-110CNT-GNP. Further CV measurements have been conducted to evaluate the catalytic response of
TvLac depending on the specific surface modification,
Figure 3, in sodium acetate buffer (AcONa) 0.01 mM, pH 5 with catechol 0.02 mM, in presence and absence of enzyme.
When the enzyme is employed the voltammograms showed diffusion controlled anodic and cathodic peaks for the oxidation and reduction of the catechol at the electrode surface. In the presence of
TvLac, if the catalysis regularly takes place, the voltammograms shapes become characteristic sigmoidal curves. Because of the oxidation of the catechol is catalyzed by the enzyme, the anodic peak disappears, while the cathodic peak becomes broader, indicating higher amount of catechol reduced at the electrode surface [
35]. Following the steady state catalytic current and the sigmoid morphology, the performance of each electrode was evaluated. The platform DRP-110 showed a cathodic current of 0.8 μA, while the additional presence of CNTs in DRP-110CNT (
Figure 3a,b) induces improvement of current to 1.7 μA, broadening the cathodic peak. In fact, the higher electronic conductivity and electroactive area provided by the presence of CNTs result in a faster electron transfer (ET) at the electrode surface.
The DRP-110GNP electrode showed an increased catalytic current of 3 μA,
Figure 3c, despite the lower electroactive area in comparison with the DRP-110CNT electrode. This behavior is due to a stronger direct ET communication between the
TvLac and the WE surface promoted by GNPs, known to work as highly efficient electron-conducting tunnels between WE and enzymes [
36]. In addition, GNPs have remarkable affinity with proteins, allowing a suitable surface coverage and preserving the catalytic activity at the same time [
37,
38,
39,
40]. By combining the GNPs and CNTs in the DRP-110CNT-GNP platform, the highest catalytic current of 4.05 μA is reached,
Figure 3d. The catalytic behavior is further investigated through CA tests, carried out measuring the cathodic catalytic current produced through successive additions of equal volumes of a catechol solution at 255 mM over time,
Figure 4a,b.
The DRP-110/
TvLac+GA sensor showed the lowest signal-to-noise ratio (S/N) in the CA profile and the less stability of signal over time,
Figure 4a in grey, as the catalytic current was constant after 20 μM of catechol, impeding the production of an acceptable Michaelis-Menten hyperbole and the kinetic parameters calculation. This behavior can be explained as consequence of three factors: i) immobilization of lower amount of enzyme as the smallest electroactive area and thus the smallest roughness and porosity of the unmodified carbon electrode surface (A
el=2.04 mm
2, ρ=0.16); ii) consequent fast saturation of enzyme active sites; iii) possible enzyme loss in solution as weaker physical adsorption on electrodes material and the poor adhesion of the glutaraldehyde matrix on the same. The DRP-110GNP/
TvLac+GA platform revealed better sensitivity to the substrate as the higher slope of CA profile in respect to that of DRP-110,
Figure 4a in black, as highlighted before. However, there is a low stability over time of the CA signal and the S/N ratio, probably meaning that the roughness (ρ=0.68) and porosity of the WE surface are not suitable enough to guarantee a long permanence of the GA cross-linked aggregate. Also in this case the kinetic parameters could not be determined. Differently, the presence of CNTs in DRP-110CNT/
TvLac+GA significantly favors the stability of the catalytic activity over time till 400 μM of catechol,
Figure 4b in grey. In fact, the well-known high porosity of CNTs plays a dramatic role in the physical adsorption mechanism as well as in the immobilization stability [
41,
42,
43,
44,
45]. However, catalytic current increases slower in comparison to DRP-110GNP/
TvLac+GA electrode, index of lower sensitivity. The co-presence of GNPs and CNTs in DRP-110CNT-GNP/
TvLac+GA,
Figure 4b in black, guarantees an optimal direct ET at the surface, high enzyme loading and good immobilization stability at the same time, providing for a sensitive and stable CA response over time. For the last two platforms the Michaelis-Menten hyperboles were obtained (
Figure 4c), calculating the kinetic parameters and range of linear response (
Figure 4d) for both.
The DRP-110CNT-GNP/
TvLac+GA was chosen for the realization of the smartphone-based biosensor for Phs detection in coffee blends, with a linear range of 5-70 μM, analytical sensitivity (a) equal to 0.124 μA/μM and limit of detection (LOD) of 2.99 μM. LOD calculated from the analytical sensitivity through the following equation: LOD = kσ
B/a, where σ
B is the standard deviation of the blank measurements and k is chosen relatively to the confidence level required [
46], in this work k = 3. Also, the stability over time of the biosensor was evaluated giving a 94% retention of the response after 21 days (
Figure S4). Total Phs content was then evaluated in seven different coffee blends (
Table S2), differing on roasting degrees and origin. More specifically, the roasting degree influence was studied by testing a green coffee sample treated in oven under different temperatures and times. All coffee blends were analysed using the biosensor device and recording a CA with constant applied potential of -0.1 V vs. AgPRE.
The method used for the analysis was the Standard Addition Method: the current intensity value of the blank (I
0) is recorded and, sequentially, the unknown sample followed by the standard catechol solutions at 0, 25, 40, 50, 65 μM were added. As a reference, the CA profile for the analysis of the coffee blend II and the corresponding graph of standard additions are reported in
Figure 5a,b. The I
cat related to each standard sample was plotted
vs the corresponding catechol standard concentration and the one related to the unknown sample was plotted in the graph at [catechol]=0 μM. The Phs concentration was determined by extrapolation of the line passing through the points, taking the dilution factor into account. The reference method employed for the validation of the biosensor is the Folin-Ciocâlteu spectrophotometric assay (calibration curve in
Figure 5c) by diluting each sample 200:1. All the values of Phs content obtained with the two methods are reported in
Table S3. The experimental data were finally converted into mg/g defined as the total concentration of Phs expressed as gallic acid on the total grams (7g) of fresh coffee used to prepare the sample. The analysis of the actual samples thus led to satisfactory results in agreement with work in the literature [
19,
47], in which the concentration of polyphenols was determined by spectrophotometric methods.
The values determined using the optimized biosensor resulted in great agreement with those obtained by Folin-Ciocâlteu, with a mean recovery of 95%. Looking at the results, the dependence of Phs content on both the coffee variety (
C. robusta or
C. arabica) and roasting degree (light, medium, dark) was confirmed. In fact, the roasting temperature causes the chlorogenic acid degradation, reducing the amounts of malic and citric acid thus quinic acid concentration [
48,
49]. Specifically, by considering the blends I, II and III, characterized by the same variety (100%
C. arabica), a higher roasting degree led to greater degradation of Phs resulting in a lower content. The samples IV, V and VI belong to the same manufacturer are characterized by different percentages of
C. arabica and
C. robusta. The highest Phs concentrations were found in coffee blends with higher amount of
C. robusta, as reported in the literature [
50,
51]. The sample VI and II, despite belonging to different manufacturers, showed similar Phs concentrations, as characterized by a medium roasting degree and a composition of 100%
C. arabica. Interesting the case of coffee blend VII, which is composed by 98% of
C. arabica variety and 2% of
C. robusta green, presenting a higher Phs content than blends II and I, but similar to blend III; in fact, the presence of green coffee in a very low percentage did not significantly influence the Phs content. The last blend analyzed is the VIII, a green coffee 100%
C. arabica, roasted under different temperature and time conditions. The results reported in
Table S3 and
Figure 5d besides being in good agreement with those obtained through Folin-Ciocâlteu reference method, are also consistent with the behavior largely treated in literature [
47,
52]. The light roasted samples generally show the highest Phs content [
53], even compared with the relative unroasted samples. In fact, thermal treatments are responsible for both the antioxidants degradation and the production of new antioxidants species such as heterocyclic compounds obtained through the Maillard reaction. However, as breakdown of cellular components during thermal processes is responsible for the release of bound phenolic acids, the higher the roasting time the lower the total Phs content, under all roasting temperatures.